实验四 用窗函数法设计fir滤波器 实验报告.pdfVIP

yu 实验四 用窗函数法设计 FIR 滤波器 （一）实验目的 1. 掌握窗函数法设计 FIR 滤波器的原理和方法，观察用几种常用窗函数设计的 FIR 数字滤波器技术指标； 2. 掌握 FIR 滤波器的线性相位特性； 3. 了解各种窗函数对滤波特性的影响。 （二）实验原理 如果所希望的滤波器的理想频率响应函数为 Hd(ejω) ，则其对应的单位脉冲 1  响应为hd (n)   H (ej  )ej n d ，用窗函数 wN(n)将 hd(n)截断，并进行加权 2  处理，得到实际滤波器的单位脉冲响应 h(n)=hd(n)wN(n) ，其频率响应函数为 N 1 H (ej  )  h(n)ej n 。 如 果 要 求 线 性 相 位 特 性 ， 则 h(n)还 必 须 满 足 n0 h(n)  h(N  1 n) 。可根据具体情况选择h(n) 的长度及对称性。 （三）实验内容 1、生成四种窗函数：矩形窗、三角窗、汉宁窗、海明窗，并观察其频率响应。 实验代码以及运行结果 %矩形窗及其频响 n=15; window1=rectwin(n); [h1,w1]=freqz(window1,1); subplot(2,1,1); stem(window1);title(矩形窗);subplot(2,1,2); plot(w1/pi,20*log(abs(h1))/abs(h1(1)));title(矩形窗频响); %三角窗及其频响 n=15; window2=triang(n);[h2,w2]=freqz(window2,1); subplot(2,1,1);stem(window2);title(三角窗);subplot(2,1,2); plot(w2/pi,20*log(abs(h2))/abs(h2(1)));title(三角窗频响); yu %汉宁窗及其频响 n=15; window3=hann(n); window3=hann(n); [h3,w3]=freqz(window3,1); subplot(2,1,1);stem(window3);title(汉宁窗);subplot(2,1,2); plot(w3/pi,20*log(abs(h3))/abs(h3(1)));title(汉宁窗频响); %海明窗频响 n=15; window4=hamming(n);[h4,w4]=freqz(window4,1); subplot(2,1,1);stem(window4);title(海明窗);subplot(2,1,2); plot(w4/pi,20*log(abs(h4))/abs(h4(1)));title(海明窗频响); 运行结果： yu 2 、根据下列技术指标，设计一个FIR 数字低通滤波器：wp=0.2 π，ws=0.4 π， ap=0.25dB ， as=50dB ，选择一个适当的窗函数，确定单位冲激响应，绘出所设 计的滤波器的幅度响应。 提示：根据窗函数最小阻带衰减的特性表，可采用海明窗可提供大于 50dB 的衰 减，其过渡带为 6.6 π/N ，因此具有较小的阶次。 实验代码 %用窗函数法设计FIR滤波器 clear all; N=33; Wc=0.3*pi; %理想低通滤波器设计(怎样设计理想高通、带通、带阻滤波器？) alpha=(N-1)/2; n=[0:(N-1)]; m=n-alpha+eps; hd=sin(Wc*m)./(pi*m); B=hamming(N);%海明窗 string=[Hamming ,N=,num2str(N)]; h=hd.*(B); %加窗截取 % 以上过程可直接调用FIR1实现 %wc=Wc/pi;%频率归一化 %h=fir1(N-1,wc,hamming(N)); [H,m]=freqz(h,[1],1024,whole); %频率响应 mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(