吉布斯现象( 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采样已经可以很好 ...
随机推荐
- 【转】libvirt kvm 虚拟机上网 – Bridge桥接
libvirt kvm 虚拟机上网 – Bridge桥接 2013 年 7 月 3 日 / 东东东 / 暂无评论 目录 [hide] 1 Bridge桥接原理 2 在host机器配置桥接网络 2.1 ...
- 【转】String.format详解
一.前言 String.format 作为文本处理工具,为我们提供强大而丰富的字符串格式化功能,为了不止步于简单调用 String.format("Hello %s", " ...
- showdialog窗体不在任务栏显示的问题处理
场景: c#开发的windows窗体用showdialog弹出时,在任务栏中 win7系统显示,win xp和win 2003却不显示. 窗体的ShowInTaskbar已设置为True. 解决: 在 ...
- Hadoop学习1--解决启动过程中的问题
方法:http://www.aboutyun.com/thread-12694-1-1.html http://www.linuxidc.com/topicnews.aspx?tid=13 http: ...
- DG_Oracle DataGuard Switchover主备节点切换(案例)
2014-06-09 Created By BaoXinjian Thanks and Regards http://wenku.baidu.com/view/dc9f00d349649b6648d7 ...
- 前端实现 SVG 转 PNG
http://fex.baidu.com/blog/2015/11/convert-svg-to-png-at-frontend/ 前言 svg 是一种矢量图形,在 web 上应用很广泛,但是很多时候 ...
- powerdesigner中将表的name在生成建表sql时生成注释
1.为powerdesigner的表设置注释方法: powerdesigner默认没有注释: 设置方法: 选择那个表 右键- >Properties- >Columns- >Cust ...
- 通过位运算生成ID
public static void main(String[] args) { long serverId = 65535; System.out.println("ServerId:&q ...
- R %operator% 含义
%foo% is the syntax for a binary operator. In base R: %in%: '"%in%" <- function(x, tabl ...
- windows service的继承类ServiceBase
https://msdn.microsoft.com/zh-cn/library/system.serviceprocess.servicebase.exitcode(v=vs.80).aspx 在停 ...