数字信号实验报告.FIR数字滤波器设计说明书.doc

数字信号处理课程报告 题 目： FIR数字滤波器设计 一、实验目的 1．熟悉线性相位FIR数字滤波器特性。掌握用窗函数法设计FIR数字滤波器的原理和方法各种窗函数对滤波性能的影响。 掌握用法设计FIR数字滤波器的原理和方法滤波性能的影响。 仪器及材料计算机，MATLAB软件 实验及用窗函数法设计FIR数字滤波器FIR数字滤波器 选取合适的窗函数编程设计，求滤波器的单位脉冲响应h(n)，出h(n)的幅频特性和相频特性曲线。用设计FIR数字滤波器求滤波器的单位脉冲响应h(n)，出h(n)的幅频特性和相频特性曲线。 clc,clear % wp=0.2*pi;Rp=0.25dB; ws=0.3*pi;As=50dB; % 查第237页表7.1可知用Hamming窗、Blackman窗均可(最小阻带衰减=50dB)， % 但Hamming窗具有较小的过渡带(6.6pi/M)，故选择Hamming窗。 wp=0.2*pi;ws=0.3*pi; tr_width=ws-wp; M=ceil(6.6*pi/tr_width)+1 %利用Hamming窗的过渡带=(6.6pi/M),求窗宽M n=0:1:M-1; wc=(ws+wp)/2; %求通带截止频率wc hd=ideal_lp(wc,M); %求理想低通冲激响应hd(n) w_ham=(hamming(M)); h=hd.*w_ham; %求FIR滤波器低通冲激响应h(n) [H,w]=freqz(h,1,1000,whole); db=20*log10((abs(H)+eps)/max(abs(H))); %求FIR滤波器频响的dB值 delta_w=2*pi/1000; %将2pi等分1000份 Rp1=-min(db(1:1:(wp/delta_w+1))) %求Passband Ripple As1=-max(db(ws/delta_w+1:1:501)) %求Min Stop attenuation %plots figure(3);clf; subplot(131); stem(n,hd); title(Ideal Impulse Reponse hd(n)); axis([0,M-1,-0.1,0.3]); ylabel(hd(n)); set(findobj(gcf,Type,line,Color,[0 0 1]), Color,m, LineWidth,2); subplot(132); stem(n,w_ham); title(Hamming Window w(n)); axis([0,M-1,0,1.3]); ylabel(w(n)); set(findobj(gcf,Type,line,Color,[0 0 1]), Color,r, LineWidth,2); subplot(133); stem(n,h); axis([0,M-1,-0.1,0.3]); ylabel(h(n));title(Actual Impulse Response h(n)); set(findobj(gcf,Type,line,Color,[0 0 1]), Color,b, LineWidth,2); figure(1);clf; subplot(221); stem(n,hd); title(Ideal Impulse Reponse); axis([0,M-1,-0.12,0.3]); ylabel(hd(n)); set(findobj(gcf,Type,line,Color,[0 0 1]), Color,m, LineWidth,1); subplot(222); stem(n,w_ham); title(Hamming Window); axis([0,M-1,0,1.3]); ylabel(w(n)); set(findobj(gcf,Type,line,Color,[0 0 1]), Color,g, LineWidth,1) subplot(223); stem(n,h); axis([0,M-1,-0.1,0.3]); ylabel(h(n));title(Actual Impulse Response); set(findobj(gcf,Type,line,Color,[0 0 1]), Color,r, LineWidth,1); subplot(224); plot(w/pi,db); axis([0,1,-100,10]); yla