FIR滤波器的设计.pptxVIP

实验三：FIR 数字滤波器的设计 实验目的 掌握用窗函数法，频率采样法及优化设计法设计FIR滤波器的原理及方法。 熟悉线性相位FIR滤波器的幅频特性和相频特性。 了解各种不同窗函数对滤波器性能的影响。 一、 实验内容 1. N=45，计算并画出矩形窗、汉明窗、布莱克曼窗的归一化的幅度谱，并比较各自的主要 特点。 clearall; N=45; wn1=kaiser(N,0); wn2=hamming(N); wn3=blackman(N); [h1,w1]=freqz(wn1,N); [h2,w2]=freqz(wn2,N); [h3,w3]=freqz(wn3,N); plot(w1/pi,20*log10(abs(h1)),r-,w2/pi,20*log10(abs(h2)),b-,w3/pi,20*log10(abs(h3)),g-); axis([0,1,-120,10]);grid; xlabel(归一化频率/\pi);ylabel(幅度/dB); title(三种窗口函数); legend(矩形窗,汉明窗,布莱克曼窗,3);; ; ; ;些差别。汉宁窗在这方面具有较好的特性。 4. 用Kaiser窗设计一专用线性相位滤波器，N=40，当β=4、6、10时，分别设计、比较它 们的幅频和相频特性，注意β取不同值时的影响。 ;; 总结：Kaiser窗的beta值越大，过渡带的带宽越宽，阻带的衰减越是厉害；beta值越小， 过渡带的带宽越窄，阻带的衰减越是却有所降低。所以在实际应用的时候，要权衡过渡 带带宽和阻带衰减，以获得最优性能！ 5. 用频率采样法设计(4)中的滤波器，过渡带分别设一个过渡点，令H(k)=0.5。比较两种 不同方法的结果。 clearall; N=40; Hk=[zeros(1,3)0.5ones(1,5)0.5zeros(1,1)0.5ones(1,5)0.5zeros(1,5)-0.5-ones(1,5)-0.5zero s(1,1)-ones(1,5)-0.5zeros(1,3)]; k=0:N-1; hn=real(ifft(Hk.*exp(-j*pi*(N-1)*k/N))); [Hw]=freqz(hn,1); plot(w/pi,20*log10(abs(H))); axis([01-8010]);grid; xlabel(归一化频率/\pi) ylabel(幅度/dB) title(频率采样法设计专用线性相位滤波器);;总结：利用频率采样，获得的滤波器通带内波动较好，但衰减有所降低，过渡带也比较 宽！ 6. 用雷米兹(Remez)交替算法设计(4)中的滤波器，并比较(4)、(5)、(6)三种不同方法的 结果。 clearall; N=40; f=[00.150.20.40.450.550.60.80.851]; a=[0011001100]; wt=[21212]; b=remez(N-1,f,a,wt); [h,w]=freqz(b,1); plot(w/pi,20*log10(abs(h))); axis([01-7010]);grid; xlabel( 归 一 化 频 率 /\pi) ylabel(幅度/dB) title(雷米兹交替算法设计专用线性相位滤波器);; ;