吉布斯现象( Gibbs)
在连续傅里叶级数(或积分)变换中,信号所对应的离散频谱(或连续频谱)为(或),其频率是无限离散分布的(或频谱的分布范围是无限区间的)。很显然,单位时间内,频率较低(简称低频,即较小)的简谐波相对频率较高(简称高频,即较大)的简谐波在空间的变化要平稳得多。例如,时所对应的直流分量在空间是不变化的(信号在整个区间的平均值),其它成分的信号则随频率的增大而更加快速变化。
对于一个在有限区间分布的信号,其连续频谱在频率域的分布往往是无限区间的。实际信号处理时,我们通常只能在有限区间内做傅里叶分析(除非理论分析),也就是说,我们只能取有限区间来替代理论分析中的无限区间,多数情况下,我们总是选择信号的低频部分,而舍弃高频部分。信号的高频部分往往反映的是信号的快速变化特征,如果信号本身是连续的,这样做一般不会引起信号的显著变化;可是,如果信号的高频成分比较丰富、比较重要,特别是在信号本身存在较为明显的剧烈突变时,这样做自然就会引起一定的误差。让我们举例对它进行分析。
设有一方波信号,其表达式为
在有限区间上该信号的连续傅里叶级数变换所对应的离散频谱为
并且有
实际工作中,我们不可能进行无限相的计算问题。通常我们只能采取一定程度的近似逼近,即用一个正整数M来代替上式中的无穷大
这样得到的有限频谱的逼近信号与原始信号见图。
从图中可以看到,在原始信号的突变点处,逼近信号出现了明显的振荡现象,随着M的增大,这些振荡并没有消失,而是更加集中于突变点附近。这种在突变点处出现的振荡现象被称为吉布斯(Gibbs)现象,它是由于在反变换的计算过程中用有限项近似无限项从而丢失原始信号中的高频成分所致。吉布斯现象在信号的变换及滤波器的设计和应用中极为普遍。
(生成此图形的软件名:Gibbs_Phenomena_CFST.m)
在连续傅里叶变换积分变换中,是否存在Gibbs现象?回答是肯定的,通过在CFT中采用有限区间的频率范围计算
可以得到不同的结果(见下图)。
(生成此图形的软件名:Gibbs_Phenomena_CFT.m)
(读者可认真比较以上两类信号的傅里叶级数与傅里叶积分之区别。)
吉布斯现象( Gibbs)的更多相关文章
- Mel倒谱系数
Mel倒谱系数:MFCC Mel频率倒谱系数(Mel Frequency Cepstrum Coefficient)的缩写是MFCC,Mel频率是基于人耳听觉特性提出来的,它与Hz频率成非线性对应关系 ...
- PRML读书会第十一章 Sampling Methods(MCMC, Markov Chain Monte Carlo,细致平稳条件,Metropolis-Hastings,Gibbs Sampling,Slice Sampling,Hamiltonian MCMC)
主讲人 网络上的尼采 (新浪微博: @Nietzsche_复杂网络机器学习) 网络上的尼采(813394698) 9:05:00 今天的主要内容:Markov Chain Monte Carlo,M ...
- MC, MCMC, Gibbs采样 原理&实现(in R)
本文用讲一下指定分布的随机抽样方法:MC(Monte Carlo), MC(Markov Chain), MCMC(Markov Chain Monte Carlo)的基本原理,并用R语言实现了几个例 ...
- 随机采样方法整理与讲解(MCMC、Gibbs Sampling等)
本文是对参考资料中多篇关于sampling的内容进行总结+搬运,方便以后自己翻阅.其实参考资料中的资料写的比我好,大家可以看一下!好东西多分享!PRML的第11章也是sampling,有时间后面写到P ...
- LDA-math-MCMC 和 Gibbs Sampling
http://cos.name/2013/01/lda-math-mcmc-and-gibbs-sampling/ 3.1 随机模拟 随机模拟(或者统计模拟)方法有一个很酷的别名是蒙特卡罗方法(Mon ...
- MC, MCMC, Gibbs採样 原理&实现(in R)
本文用讲一下指定分布的随机抽样方法:MC(Monte Carlo), MC(Markov Chain), MCMC(Markov Chain Monte Carlo)的基本原理,并用R语言实现了几个样 ...
- 如何做Gibbs采样(how to do gibbs-sampling)
原文地址:<如何做Gibbs采样(how to do gibbs-sampling)> 随机模拟 随机模拟(或者统计模拟)方法最早有数学家乌拉姆提出,又称做蒙特卡洛方法.蒙特卡洛是一个著名 ...
- LDA Gibbs Sampling
注意:$\alpha$和$\beta$已知,常用为(和LDA EM算法不同) 1. 为什么可用 LDA模型求解的目标为得到$\phi$和$\theta$ 假设现在已知每个单词对应的主题$z$,则可 ...
- MCMC(四)Gibbs采样
MCMC(一)蒙特卡罗方法 MCMC(二)马尔科夫链 MCMC(三)MCMC采样和M-H采样 MCMC(四)Gibbs采样 在MCMC(三)MCMC采样和M-H采样中,我们讲到了M-H采样已经可以很好 ...
随机推荐
- Mongodb集群搭建及spring和java连接配置记录
一.基本环境: mongdb3.0.5数据库 spring-data-mongodb-1.7.2.jar mongo-java-driver-3.0.2.jar linux-redhat6.3 tom ...
- CentOs of Tomcat commands
1.启动tomcat 进入tomcat目录的bin目录下 执行 [root@iZ253lxv4i0Z bin]# ./startup.sh Using CATALINA_BASE: /usr/loca ...
- PLSQL_闪回操作2_Fashback Version Query
2014-12-09 Created By BaoXinjian
- NeHe OpenGL教程 第二十七课:影子
转自[翻译]NeHe OpenGL 教程 前言 声明,此 NeHe OpenGL教程系列文章由51博客yarin翻译(2010-08-19),本博客为转载并稍加整理与修改.对NeHe的OpenGL管线 ...
- apache/php 开启 gzip压缩
1.php方式开启 原理: header("Content-Encoding: gzip"); echo gzencode('songjiankang'); 示例1: functi ...
- js实现加减乘除
/** ** 除法函数,用来得到精确的除法结果 ** 说明:javascript的除法结果会有误差,在两个浮点数相除的时候会比较明显.这个函数返回较为精确的除法结果. ** 调用:accDiv(arg ...
- C++学习24 虚析构函数
在C++中,构造函数用于在创建对象时进行初始化工作,不能声明为虚函数.因为在执行构造函数前对象尚未创建完成,虚函数表尚不存在,也没有指向虚函数表的指针,所以此时无法查询虚函数表,也就不知道要调用哪一个 ...
- [ActionScript3.0] 逻辑或"||=" ,等于"=="和全等于"==="
function a(o:Object):void { o||=new Object(); trace(o); } //此上下两个方法作用是一样的 function b(o:Object):void ...
- editplus格式化xml文档
使用editplus打开xml文档,发现文件内容没有格式化,非常难看,使用IE9打开这个xml文档是格式化的,但是IE9打开不能编辑. 在网上找了一些资料发现能够通过一个小插件来是editplus格式 ...
- List<IPoint> to IPointCollection to IPolygon
IPointCollection 到 IPolygon的转换 IPoint pPoint = new PointClass(); //IPolygon pPolygon1 = n ...