第8章 有限脉冲响应滤波器的的设计 《数字信号处理》课件.pptVIP

程序接上页 程序运行结果如图8.8所示， 相对于阻带来说，其幅度函数A(ω)的数值在通带内近 图8.8 subplot(2,2,1);plot(w/pi,A,k);grid;axis([-1,1,-0.1,1.2]); xlabel(\omega/\pi);ylabel(A(\omega)); subplot(2,2,2);plot(w/pi,p,k);grid;axis([-1,1,-35,35]); ylabel(\theta(\omega)/radian) 似等于1，有些像理想的幅度函数Ad(ω)；滤波器h(n)的相位函数θ(ω)是过原点的直线。 造成实际频谱H(ω)和理想频谱Hd(ω)之间差别的原因是：实际的h(n)只是理想的hd(n)的一段。这一段的长度越长，H(ω)和Hd(ω)的差别就越小。 本题得到的有限脉冲响应低通滤波器h(n)，是用理想滤波器的相位不为零获得的。 44节 （2）理想滤波器的相位为零 在角频率的主值区间[-π, π)的理想频谱 根据离散时间的傅里叶反变换式(3.81)，零相位的理想低通模型式(8.53)的单位脉冲响应 它的波形如图8.9左图所示，对称中心在n=0。为了得到有限长的因果线性相位滤波器h(n)，可以直接剪下 (8.53) (8.54) hd(n)的n=-10~10的最大数值部分，它的长度N=21点；并且，将减下的部分向右移动a=(N-1)/2=10点，才能成为因果序列。这个结果的数学表达式是 它就是按题目要求设计的低通滤波器。 图8.9 (8.55) 8.3.2 截断序列的后果 实际剪裁或截断序列时，为了保证h(n)满足线性相位条件，是用一个对称的有限长序列w(n)乘以被剪裁的原始序列hd(n)，这个有限长序列w(n)叫做窗序列。不同的窗序列相当于不同的观点。 常用的窗序列有矩形窗，写为 其波形如图8.10所示。 常用的窗序列还有汉宁窗，写为 其波形如图8.10所示。 (8.56) (8.57) 常用的窗序列还有汉明窗，写为 常用的窗序列还有布莱克曼窗，写为 图8.10 (8.58) (8.59) 从时域来看，用一个短序列代替一个长序列总是有失真的。这个问题也可以从频域来看。 观察加窗序列h(n)=hd(n)w(n)的频谱H(ω)时，可利用离散时间傅里叶变换的频域卷积定理，即 该积分显示，窗序列的频谱W(ω)对截取序列的频谱H(ω)质量是有影响的。 例题8.4 已知截止频率ωc=0.2π的理想低通滤波器的单位脉冲响应是 (8.60) (8.61) 如果分别用矩形窗和汉明窗对hd(n)加窗，获取长度N=21的因果线性相位低通滤波器，请分析这两种窗口对我们设计的滤波器将带来什么频谱影响？ 解 已知窗口法设计的系统频谱H(ω)正比于理想频谱Hd(ω)和窗序列频谱W(ω)的卷积，而且根据公式(8.51)和(8.50)知道，理想低通滤波器的频谱是 那么，我们就可以分析这两种窗口对设计的滤波器带来的频谱影响。 (8.62) （1）矩形窗的影响 根据离散时间傅里叶变换，矩形窗序列wrectangle(n)的频谱是 若用矩形窗对模型序列hd(n)加窗，将得到有限脉冲响应滤波器的序列 其频谱可用离散时间傅里傅里叶变换或频域卷积定理计算。 (8.63) 根据频域卷积定理，滤波器h(n)的频谱 它的相位函数是线性的，它的幅度函数A(ω)正比于矩形窗的幅度函数Arectangle(ω-φ)在通带[-ωc, ωc]上的积分，也就是矩形窗的幅度函数曲线在通带上的面积。 为了避免计算面积，可观察其幅度函数图8.11， (8.64) 绘制图8.11的MATLAB程序是： 程序转下页 N=21;wc=0.2*pi;a=(N-1)/2;w=0:0.01:2*pi;%窗口长度、截止频率、群延时、角频率变量 Ad=[w=wc|w=(2*pi-wc)];%理想滤波器的幅度函数 subplot(3,1,1);plot([-(2-w/pi),w/pi],[Ad,Ad],k);grid;axis([-2,2,0,1.1]);%画-2π~2π的幅度函数 xlabel(\phi/\pi);ylabel(A_{d}(\phi));box;%理想滤波器的幅度函数的坐标符号 Ar=sin((w+eps)*N/2)./sin((w+eps)/2);%矩形窗序列的幅度函数 subplot(3,1,2);plot([-(2-w/pi),w/pi],[Ar,Ar],k);