实验三用窗函数法设计FIR数字滤波器..docVIP

实验三 用窗函数法设计FIR数字滤波器 实验目的及原理方法 参考（301~302）页； 实验内容 1.用MATLAB产生各种窗函数 m=200; a=boxcar(m); %矩形窗 b=hanning(m); %汉宁窗 c=hamming(m); %海明窗 d=blackman(m); %布莱克曼窗 e=kaiser(m,7.865);%凯泽窗 n=1:m; plot(n,a,r*,n,b,g+,n,c,y*,n,d,b.,n,e,mx) %画图矩形窗（红线*点表示）汉宁窗（绿色+点表示） %海明窗（黄色*号表示）布莱克曼窗（蓝色.点表示）凯泽窗（洋红色x点表示） 2.利用窗函数设计FIR滤波器 设计具有下列指标?p=0.25?,Rp=0.25dB，?s=0.3?,Rs=50dB的低通数字滤波器。 要求： 选择合适的窗函数; 画出滤波器的频率特性 几个主要的MATLAB函数 （1）[r,p,k]=residuez(b,a) 将有理分式变换成部分分式表示 （2）[b,a]=residuez(r,p,k) 将部分分式换成有理分式变表示 式中各参数意义如下： （3）freqz(b,a) 画连续系统的频率特性曲线 （4）impz(b,a，101) 画系统的脉冲响应 例如：N=5;Wn=[100 200]/500; [b,a]=butter(N,Wn); impz(b,a,101)%求滤波器系统单位脉冲响应 根据题目要求的FIR滤波器技术指标，由于要求阻带最小衰减为50db，根据书上查表可得，窗函数采用海明窗，具体MATLAB和仿真图形如下： wp=0.25*pi;ws=0.3*pi; tr_width=ws-wp; %求过渡带带宽 M=ceil(6.6*pi/tr_width)+1;%利用海明窗函数求最大阶数 n=[0:1:M-1]; wc=(ws+wp)/2; hd=ideal_lp(wc,M);%求给定理想低通滤波器单位冲激响应 w_ham=(hamming(M));%给定海明窗窗函数 h=hd.*w_ham; [db,mag,pha,grd,w]=freqz_m(h,[1]);%调用自定义幅频函数 delta_w=2*pi/1000; Rp=-(min(db(1:1:wp/delta_w+1))) %Passband Ripple Rs=-round(max(db(ws/delta_w+1:1:501)))%MinStopband attenuation %plots画图 figure(1) subplot(2,1,1);stem(n,hd,.); title(Ideal Impulse response) axis([0 M-1 -0.1 0.3]); ylabel(hd(n)) subplot(2,1,2);stem(n,h,.); title(Actual Impulse Response) axis([0 M-1 -0.1 0.3]); ylabel(h(n)) figure(2) subplot(2,1,1);stem(n,w_ham,.); title(hamming Window) axis([0 M-1 0 1.1]);ylabel(w(n)) subplot(2,1,2);plot(w/pi,db); title(Magnitude Response in dB); grid %增加网格 axis([0 1 -100 10]);ylabel(Decibels) 运行程序验证选择海明窗设计时最终所得通带和阻带衰减如下，Rs=52dB50dB满足题意： 其中程序中用到的两个自定义函数在MATLAB中定义成.m文件如下： function [db,mag,pha,grd,w]=freqz_m(b,a); [H,w]=freqz(b,a,1000,whole);%求连续系统的频率特性曲线 H=(H(1:501));w=(w(1:501)); mag=abs(H); db=20*log10((mag+eps)/max(mag));%幅值用db表示 pha=angle(H);%求相角 grd=grpdelay(b,a,w);%群延迟响应 实验总结 整个实验是运用窗函数法设计FIR数字滤波器并进行相应的MATLAB仿真，从实验内容1中我们可以看到同一阶数中不同窗函数的误差函数曲线有一定差别，而在设计滤波器时我们应该根据所要求的技术指标选择合适的窗函数才能减小设计误差。在此次实验中我收获很多，进一步熟悉了MATLAB软件的应用，同时更加了解了设计FIR数字滤波器的窗函数法。  数字