吉布斯效应研究报告.docVIP

吉布斯效应研究报告 学院：通信工程学院 年级专业： 学号： 学生： 指导教师： 时间：2015.4.13 吉布斯效应研究报告 摘 要： 关键词：将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯现象。 图一k=0.01 N=20 w=8.99% 图二 k=0.01 N=200 W=8.77% 图三 k=0.01 N=2000 w=8.94% 图四 k=0.001 N=20 w=8.99% 做多次重复实验，可得以下表格： 实验次数 k N w 1 0.01 20 8.99% 2 0.01 200 8.77% 3 0.01 2000 8.94% 4 0.001 20 8.99% 5 0.001 200 8.95% 由上可得，随着矩形窗的宽度和采样点数的增加，峰起值占总跳变值的比值越来越接近9%，在x轴方向跳变点越来越接近平衡点，吉布斯效应越来越明显。 b.锯齿波分析 其傅里叶展开形式如下： MATLAB源程序如下： k=input(取点间隔k=); N=input(选取项数N=); t=-2*pi:k:2*pi; s=sawtooth(t-pi); y=0; for n=1:1:N, y=y+2/pi*sin(n*t+(n+1)*pi)/n; %前n项和 end plot(t,s,t,y); w=(max(y)-1)/2 %峰起值占总跳变值的比例 锯齿波产生的吉布斯效应与方波的类似，可以得出相同的结论，此处不再赘述。 频域吉布斯效应 和上面相反，如果一个信号在频域是一个“方波”，由于数字域频率以2*pi为周期，这个“方波”可以在-pi到pi范围讨论。课本上都是以理想低通滤波器作为例子，其实我们可以拿理想带通、理想高通滤波器作“频域方波”的形式。当频域波形确定后，我们需要在物理上实现它，那就得得到它的系统函数。先做逆变换，变成时域信号h(n)，但长度为无限长，不能物理实现，于是采用截断。一开始人们就想“咔嚓”截一大段来逼近理想情况。很快发现，这种做法是比较武断的，逼近的波形出现了波动，而且有了过渡带。 MATAB源程序如下 wp=pi/2;ws=pi/4; Bt=wp-ws; N0=ceil(6.2*pi/Bt); N=N0+mod(N0+1,2); wc=(wp+ws)/2/pi; hn=fir1(N-1,wc,low,boxcar(N)); yn=fft(hn,512); plot(abs(yn)); MATLAB图像如下： 4 吉布斯效应实例——数字滤波中的吉布斯效应 （1）吉布斯现象产生的原因 理想滤波器的时间函数h(n),它的长度是无限的。在实际滤波中，只能取h(n)的有限部分，不可能在无限长的时间内不停进行，应该是限定在某个起止时刻，滤波因子只能取有限长度。用有限长度的滤波因子却达不到理想的频率响应所给出的滤波效果。由于h(n)是偶函数，取-N到N之间的部分，把|h|N的部分截掉，其对应的频谱如图2.从图2中可以看出在f1、f2及-f1、-f2左右处，曲线产生了较为严重的振荡现象，具有许多小的突起，最大突起在截止频率附近。 产生这种现象的原因有两个：（1）h(n)的频谱在f1、f2、-f1、-f2处产生突变（2）曲线截断引起，把无限长的时间函数h(n)截为有限长的时间函数。在第一个原因的内在条件下，有限和无限的矛盾导致了吉布斯现象。 图3为在截止频率f1、f2、-f1、-f2不变的条件下，改变截止长度N对频率响应的影响。从中可以得到，当滤波因子截断长度N改变时，滤波器的频率响应的截止限度相应的改变，突起数目也随之改变，但最大突起的相对幅度与截断长度的大小无关，皆为正负9%左右。 吉布斯效应对滤波效果的影响 滤波器产生吉布斯现象，引起滤波器的泄露效应，造成滤波效果不好，不能有效压制干扰，还可能使通频带信号的频率产生畸变。 吉布斯现象的消除 消除吉布斯的有效方法有： 镶边法：在原来不连续的理想频率响应的每边镶上两个连续变化的边。 乘因子法：在时间域上将原脉冲响应乘上一个衰减因子。 5 心得体会 通过对吉布斯效应的研究，对傅里叶变换有了进一步了解，同时更加熟悉MATLAB软件的操作。在仿真过程中，遇到了一些问题，才发现自己有很多不足，对很多知识都是一知半解，需要加强。所以，边做边温习课本，同时查阅网上的资料，自己