数字信号处理实验八.docVIP

数字信号处理实验报告 实验名称：用窗口法设计FIR数字 滤波器 姓名：周云杰 学号 指导老师：杨萌 黄怡 实验时间：2014/12/22 选课时间：周一3—5节 实验目的： 了解一个实际滤波器的设计过程，加深掌握窗口法设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响。 实验原理： 设所希望得到的滤波器的理想频率响应为。那么FIR滤波器的设计就在于寻找一个频率响应为去逼近。在这种逼近中最直接的一种方法是从单位冲激响应序列着手，使逼近理想的单位冲激响应。我们知道可以从理想频率响应通过付利叶反变换得到，即： 但是一般来说，这样得到的单位冲激响应往往都是无限长序列，而且是非因果的，我们以一个截至频率为的线性相位理想低通滤波器为例来说明。设低通滤波器的时延为，即： 则 这是一个以为中心的偶对称的无限长非因果序列。这样一个无限长的序列怎样用一个有限长的序列去接近呢？最简单的办法就是直接截取一段来代替他。例如把到的一段截取来作为，但是为了保证得到的是线性相位滤波器，必须满足的对称性，所以时延应取长度的一半，即。 这种直接截取的办法可以形象地想象为，好比是通过一个窗口看到的一段。为和一个窗口函数的乘积。在这里，窗口函数就是矩形序列，即： 但是一般来说，窗口函数不一定是矩形序列，可以在矩形以内对进行加权处理，因此，一般可表示为 这里的就是窗口函数。这种对理想单位冲激响应加窗处理对频率响应会产生以下三点影响： （1）使理想特性不连续的边沿加宽，形成一过渡带，过渡带的宽度取决于窗口频谱的主瓣宽度。 （2）在过渡带两旁产生的肩峰和余振，他们取决于窗口频谱的旁瓣；旁瓣越多，余振也越多；旁瓣相对值越大，肩峰则越强。 （3）增加截取长度N，只能缩小窗口频谱的主瓣宽度而不能改变旁瓣的相对值；旁瓣与主瓣的相对关系只决定于窗口函数的形状。因此，增加N职能相应减小过渡带宽，而不能改变肩峰值。肩峰值的大小直接决定通带内的平稳和阻带内的衰减，对滤波器性能有很大关系。例如，在矩形窗下，肩峰达8.95%，致使阻带衰减只有21dB，这在工程上往往是不够的。怎样才能改善阻带衰减特性呢？只能从改善窗口形状上找出路，所以希望窗口函数的频谱中应该减少旁瓣，使能量集中在主瓣，这样可以减少肩峰和余振，提高阻带衰减。而且要求主瓣宽度尽量窄，以获得较陡的过渡带，然而这两个要求总不能兼顾，往往需要增加主瓣宽度换取旁瓣的抑制，于是提出了海明窗、凯塞窗、切比雪夫等窗口函数。 实验结果 1 N=21; Wc=0.5*pi; b=fir1(N-1,Wc/pi,hamming(N)); freqz(b,1,512); 2 N=21; wc=0.5*pi; r=(N-1)/2; n=0:N-1; hdn=sin(wc*(n-r))/pi./(n-r); if rem(N,2)~=0; hdn(r+1)=wc/pi; end wn=hamming(N); h=hdn.*wn; H=fft(h,512); w=2*[0:511]/512; subplot(3,1,1); plot(hdn); xlabel(hdn); subplot(3,1,2); plot(w,abs(H)); xlabel(abs(H)); subplot(3,1,3); plot(w,20*log10(abs(H))); xlabel(20*log10(abs(H))); 3 N=21; wc=0.5*pi; r=(N-1)/2; n=0:N-1; hdn=sin(wc*(n-r))/pi./(n-r); if rem(N,2)~=0; hdn(r+1)=wc/pi; end wn=hamming(N); h=hdn.*wn; H=fft(h,512); Wn=fft(wn,100); magWn=abs(Wn); phaWn=angle(Wn); w=2*[0:511]/512; subplot(3,1,1); plot(wn); xlabel(wn); subplot(3,1,2); plot(magWn); xlabel(magWn); subplot(3,1,3); plot(20*log10(magWn)); xlabel(20*log10(magWn)); 实验总结 在处理数据时，选用窗函数一般遵循以下两个原则：一是主瓣应尽量窄，能量尽可