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

试验四 用窗函数法设计FIR滤波器 （一）试验目标 1. 掌握窗函数法设计FIR滤波器原理和方法，观察用多个常见窗函数设计FIR数字滤波器技术指标； 2. 掌握FIR滤波器线性相位特征； 3. 了解多种窗函数对滤波特征影响。 （二）试验原理 假如所期望滤波器理想频率响应函数为Hd(ejω)，则其对应单位脉冲响应为，用窗函数wN(n)将hd(n)截断，并进行加权处理，得到实际滤波器单位脉冲响应h(n)=hd(n)wN(n)，其频率响应函数为。假如要求线性相位特征，则h(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(三角窗频响); %汉宁窗及其频响 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(海明窗频响); 运行结果： 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(H); subplot(2,2,1) n=0:N-1; stem(n,h,.) axis([0 N-1 -0.1 0.3]) hold on n=0:N-1; x=zeros(N); plot(n,x,-) hold off xlabel(n) ylabel(h(n)) title(实际低通滤波器h(n)) text(0.3*N,0.27,string) subplot(2,2,2) plot(m/pi,db) axis([0 1 -100 0]) xlabel(w/pi) ylabel(dB) title(副频衰减特征) grid on subplot(2,2,3) plot(m,pha) hold on n=0:7; x=zeros(8); plot(n,x,-) hold off axis([0 3.15 -4 4]) xlabel(频率(rad)) ylabel(相位(rad)) title(相频特征) subplot(2,2,4) plot(m,mag) axis([0 3.15 0 1