吉布斯现象( 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采样已经可以很好 ...
随机推荐
- gcc内嵌汇编详解
[作者:byeyear 首发:cnblogs Email:east3@163.com 转载请注明] 有时候我们希望在C/C++代码中使用嵌入式汇编,因为C中没有对应的函数或语法可用. ...
- 修改mongodb oplog size
转载地址:http://blog.csdn.net/huwei2003/article/details/43307647 修改mongodb oplog size oplog简介: oplog:ope ...
- TMS320C54x系列DSP的CPU与外设——第5章 数据寻址
第5章 数据寻址 C54x DSP提供7种基本寻址方式. ■ Immediate addressing uses the instruction to encode a fixed value. ...
- struts 2 <s:select> 标签
名称 必填 缺省值 类型 描述 emptyOption false false Boolean 是否在题头选项后面添加一个空的(--)选项 headerKey false Object/Strin ...
- 黄聪:让WordPress主题支持语言本地化(使用poedit软件实现中文翻译功能)
如果你的WordPress主题要提交到WordPress官方主题库,使用者来自世界各地的多种语言,那么,你就要让你的WordPress主题支持语言本地化,方便使用者进行语言翻译和制作语言包. 让Wor ...
- Redis 内存使用优化与存储
Redis 常用数据类型 Redis 最为常用的数据类型主要有以下五种: • String • Hash • List • Set • Sorted set 在具体描述这几种数据类型之前,我们先通过一 ...
- de.greenrobot.event.EventBusException: Subscriber class dji.midware.a.e already registered to event class
java.lang.RuntimeException: Unable to create application com.android.tools.fd.runtime.BootstrapAppli ...
- dede的幻灯片除了首页其他页面也显示的方法,
<script language='javascript'> linkarr = new Array(); picarr = new Array(); textarr = new Arra ...
- python urllib模块的urlopen()的使用方法及实例
Python urllib 库提供了一个从指定的 URL 地址获取网页数据,然后对其进行分析处理,获取想要的数据. 一.urllib模块urlopen()函数: urlopen(url, data=N ...
- SVN---脱离SVN控制
创建一个记事本文件,然后吧这句话复制进去for /r . %%a in (.) do @if exist "%%a\.svn" rd /s /q "%%a\.svn&qu ...