有限长单位脉冲响应滤波器设计

1、实验五 有限长单位脉冲响应滤波器设计一、实验目的1、掌握用窗函数法、频率采样法以及优化设计法设计FIR 滤波器的原理及方法，熟悉相应的 MA TLAB 编程。2、熟悉线性相位 FIR 滤波器的幅频特性和相频特性。3、了解各种不同窗函数对滤波器性能的影响。二、实验原理window=ones(1, N) : 产生 N 点矩形窗，行向量。window=hann(N) : 产生 N 点汉宁窗，列向量。window=hanning(N):产生N点非零汉宁窗，列向量。等价于去除hann(N+2)的第一个零元素和最后一个零元素，得到的 N 点非零窗函数。window=hamming(N):产生N点海明窗，列

2、向量。window=blackman(N):产生N点布莱克曼窗，列向量。window=kaiser(N, beta) :产生参数为 beta的N点凯塞窗，列向量。M, Wd, beta, ftype=kaiserord(f, a, dev, fs): 凯塞窗参数估计。 f 为一组边界频率，最高频率为 fs/2 。 a 为 f 中各个频带的幅度值，通带取 1，阻带取 0。如果 f 中有 2 个元素，则形成3个频带，其中第1个和第3个是通带或阻带，第 2个是过渡带，a中也有 2个元素，指明第 1个和第 3个频带是通带还是阻带；如果 f 中有4个元素，则形成 5个频 带，其中1，3和5是通带或阻带，

3、2和4是过渡带，a中有3个元素，指明1，3和5是通 带还是阻带。dev的维数与a相同，指明每个频带上的波动值。fs为采样频率。M为FIR滤波器的阶数，M=N-1。Wd为归一化边界频率，等于数字边界角频率除以n，或者边界频率除以fs/2 。 beta就是凯塞窗的参数 3。ftype为滤波器的类型。b = fir1(M, Wd, 'ftype', window):用窗函数法求 FIR滤波器的系数 b(单位脉冲响应)。M为滤波器的阶数，M=N-1。Wd为一组归一化边界频率，通带和阻带间隔分布，无过渡 带；只有一个元素， 表示低通或高通滤波器； 有两个元素表示带通和带阻滤波器； 有三个

4、及 以上元素，表示多带滤波器。 'ftype' 表示滤波器类型， 'high' 表示高通滤波器， 'stop' 表 示带阻滤波器， 'DC-0' 表示多带滤波器的第一个频带为阻带， 'DC-1' 表示多带滤波器的第一 个频带为通带。 window 为窗口类型，缺省为海明窗。b = fir2(M, f, m, window):用频率采样法求 FIR滤波器的系数b。M为滤波器的阶数， M=N-1。 f 为一组归一化频率，第一个元素必须为 0，最后一个元素必须为 1 (对应奈奎斯特 频率，即采样频率的一半)，中间的元素按

5、升序排列。m的维数与f相同，指明f中每个频率上的理想幅度。window为窗口类型，缺省为海明窗。Fir2可以实现任意幅度特性的滤波器。三、实验内容1、 用窗函数法设计一个线性相位FIR高通滤波器，通带边界频率为0.7 n,阻带边界频率为0.5 n要求阻带衰减不小于 50dB,窗函数从矩形窗、汉宁窗、海明窗和布莱克曼窗中选取，且要求滤波器的阶数最小。给出该滤波器的单位脉冲响应h(n)的解析式，并用MATLAB绘出时域波形和幅频特性。2、 用窗函数法设计一个线性相位FIR带通滤波器，通带边界频率为3kHz和7kHz,阻带边界频率为2kHz和9kHz，采样频率为20kHz,要求阻带衰减不小于 40d

6、B,窗函数从矩形 窗、汉宁窗、海明窗和布莱克曼窗中选取，且要求滤波器的阶数最小。给出该滤波器的单位脉冲响应h(n)的解析式，并用 MATLAB出时域波形和幅频特性。3、 用窗函数法设计一线性相位FIR带阻滤波器，通带边界频率为600Hz和1400Hz，阻带边界频率为700Hz和1200Hz，采样频率为 4000Hz,要求阻带衰减不小于 50dB,窗函数从 矩形窗、汉宁窗、海明窗和布莱克曼窗中选取，且要求滤波器的阶数最小。给出该滤波器的单位脉冲响应h(n)的解析式，并用 MATLAB出时域波形和幅频特性。4、 分别用矩形窗、汉宁窗、海明窗和布莱克曼窗设计FIR低通滤波器，理想滤波器的 边界频率为

7、0.8 n, N=51，用Fir1函数设计。根据所设计滤波器的阻带最小衰减值确定实际阻带边界频率。5、 用频率采样法设计一个线性相位低通滤波器，N=40,通带边界频率为M4，过渡带设置一个采样点| H(k)|=0.5，给出单位脉冲响应 h(n)的序列值，并作出幅频响应。6、 用频率采样法设计一个线性相位高通滤波器，通带边界频率为2 M3，过渡带设置一 个采样点|H(k)|=0.5，分别求N=31和N=32时的单位脉冲响应 h(n)的序列值，并分别作出 幅频响应。7、 用频率采样法设计一个多带线性相位滤波器，理想幅频特性如下图所示，N=50,过 渡带设置一个采样点，过渡点的值设为0.5，窗函数选

8、矩形窗。给出该滤波器的单位脉冲响应h(n)的解析式，并用 MATLAB出时域波形和幅频特性。理想滤波器的幅频特性四、实验结果与分析时域波形:线性相位高通滤波器FIR的单位脉冲响应h(n)0.40.30.20.10-0.1-0.2-0.3-0.41015202535爭.11 '30线性相位高通滤波器FIR ,海明窗，N=33幅频特性：0-10-20-30-40B 度-50 幅-60-70-80-90-10000.10.20.30.40.50.60.70.80.91频率/Hz程序：N=33;n=0:1:N-1;wn=0.6*pi;a=(N-1)/2;hd=(s in (pi*( n-a)-

9、si n(n-a)*w n)./(pi*( n-a);hd(a+1)=0.4;win=hammi ng(N);h=wi n'.*hd;figure;stem( n,h);xlabel(' n');ylabel('h( n)');grid;title('线性相位高通滤波器FIR的单位脉冲响应h(n)');H,w=freqz(h,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB')

10、;title('线性相位高通滤波器FIR，海明窗，N=33');2.解：求通带和阻带的数字边界频率：c1 =03込；：92 = 0.7：,；：；r1 =求理想带滤波器的边界频率。理想带通滤波器的边界频率是实际滤波器幅度函数正负肩 峰频率的中心，但这里不知道正负肩峰的频率，因此用通带和阻带边界频率的中心近似：C0c1 +r1Wcor2 X = = 0.25二，'2 = c2= 0.8二2 2选择窗函数和窗口长度。阻带衰减不小于40dB,且阶数最小，因此选择汉宁窗(其阻带衰减为44dB,满足衰减要求)。滤波器的过渡带宽为，1二0.1 ,A o = o -o =0 2兀2 r

11、2 c2，选择较小过渡带宽，因此窗口长度为：62：N线性相位延迟常数为：N -1 心152根据理想边界频率-.n和线性相位延迟常数：，求理想单位脉冲响应m（n）：hd(= ：二ej:ej(7d 2sin( n 七強_sin( n七时皿n©L tT窗函数与理想单位脉冲响应相乘，得到线性相位FIR带通滤波器的单位脉冲响应时域波形:幅频特性:h(n) = w(n)hd (n) = *0.5.5cos(二 n) sin 08 二(n5)<n °25 二(n5) 15-：（n5）n=150.55nT50.60.50.40.30.20.10-0.1-0.2-0.310202530

12、15n* 1亠11i：|.i11Ay线性相位带通滤波器FIR的单位脉冲响应h（n）线性相位带通滤波器FIR ,汉宁窗，N=3100.10.20.30.40.50.60.70.80.91频率/Hz0-10-20-30-40B度-50幅-60-70-80-90-100程序：N=31;n=0:1:N-1;w1=0.25*pi;w2=0.8*pi;a=(N-1)/2;hd=(0.5-0.5.*cos(pi/15* n).*(si n(w2*( n-a)-si n(w1*( n-a)./(pi* (n-a); hd(a+1)=0.55;win=ha nnin g(N);h=wi n'.*hd;f

13、igure;stem( n,h);xlabel(' n');ylabel('h( n)');grid;title('线性相位带通滤波器FIR的单位脉冲响应h(n)');H,w=freqz(h,1);H=20*log10(abs(H);figure;plot（w/pi,H）;axis（O 1 -100 0）;xlabel（'频率 /Hz'）;ylabel（' 幅度 /dB'）;title（'线性相位带通滤波器FIR，汉宁窗，N=31'）;3.解：求通带和阻带的数字边界频率：c2=0.7二，r1 = 0

14、.35 , r2 = 0.6二求理想带滤波器的边界频率。理想带阻滤波器的边界频率是实际滤波器幅度函数正负肩峰频率的中心，但这里不知道正负肩峰的频率，因此用通带和阻带边界频率的中心近似CO + 0 dco o + o o1 = u = 0.325, 2 =亠 空二 0.652 2选择窗函数和窗口长度。阻带衰减不小于50dB,且阶数最小，因此选择海明窗（其阻带衰减为53dB,满足衰减要求）。滤波器的过渡带宽为.“ =川和一 . c1 = 0.05 ,八1 = r2 - ' C2 = 0.1：,选择较大过渡带宽，由于选择的是带阻滤波器，N必须为奇数，因此窗口长度为：6.6-N=0.1 二N

15、= 677171线性相位延迟常数为N -1根据理想边界频率-，n和线性相位延迟常数：，求理想单位脉冲响应hd（n）：二（n-：）-sin(n-Q!叩 七in( n-QiCq-sin(n ©傀n 恕71窗函数与理想单位脉冲响应相乘，得到线性相位FIR带通滤波器的单位脉冲响应h(n) = w(n)hd (n) 口0.7n 33n ；3时域波形:0.7J.1 -亠-1.2 ntrl 'iA Ge 打二、二.1L .燧也i；1 r F, r rf 上r r - I' ir r'7 _线性相位带阻滤波器FIR的单位脉冲响应h(n)0.60.50.40.30.20.10-

16、0.110205060-0.23040n70频域特性:程序：N=67;n=0:1:N-1;w1=0.35*pi;w2=0.65*pi;a=(N-1)/2;hd=(0.54-0.46.*cos(pi/33* n).*(si n( pi*( n-a)+si n(w1*( n-a)-si n(w2*( n-a)./(Pi* (n-a);hd(a+1)=0.7;win=ha nnin g(N);h=wi n'.*hd;figure;stem( n,h);xlabel(' n');ylabel('h( n)');grid;title('线性相位带阻滤波器F

17、IR的单位脉冲响应h(n)');H,w=freqz(h,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -90 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB');title('线性相位带阻滤波器FIR，海明窗，N=33');FIR低通滤波器，矩形窗4实验结果：0-10fi-1 r1-20-30-40度-50幅-60-70-80-90I1Iff00.10.20.30.40.50.60.70.80.91频率/Hz-1001.2Lj.-'J. T "

18、T 1 -| f.1 1 T-kL-r LM -.J)'-1 T1 1- '-1 11 IFIR低通滤波器的单位脉冲响应h(n)，矩形窗0.80.60.40.2150-0.21020253035404550FIR低通滤波器的单位脉冲响应h(n)，汉宁窗nFIR低通滤波器，汉宁窗-10-20-30-40-50-60-70-80-9000.10.20.30.40.50.60.70.80.9频率/Hz-1001.2-0.20FIR低通滤波器的单位脉冲响应h(n)，海明窗10.80.60.40.205101520253035404550nFIR低通滤波器，海明窗n频率/HzFIR低通滤

19、波器的单位脉冲响应h(n),布莱克曼窗0.80.70.60.50.40.30.10-0.1-0.2101520253035404550 m 二''i J亍1石A1 1*1-0.2nFIR低通滤波器，布莱克曼窗0-10-20-30-40B度-50幅-60-70-80-90-10000.10.20.30.40.50.60.70.80.91频率/Hz程序：N=51;M=N-1;n=0:1:M;Wd=0.8;w1= o nes(N,1);w2=ha nnin g(N);w3=ham min g(N);w4=blackma n( N);b1= fir1(M,Wd,'low'

20、;,w1);b2=fir1(M,Wd,'low',w2);b3=fir1(M,Wd,'low',w3);b4=fir1(M,Wd,'low',w4);figure;stem( n,b1);xlabel(' n');ylabel('h( n)');grid;title('FIR 低通滤波器的单位脉冲响应 h(n) H,w=freqz(b1,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel(' 频率 /Hz');yl

21、abel(' 幅度 /dB');title('FIR 低通滤波器，矩形窗 ');figure;stem(n,b2);xlabel('n');ylabel('h(n)');grid;title('FIR 低通滤波器的单位脉冲响应 h(n) H,w=freqz(b2,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel(' 频率 /Hz');ylabel(' 幅度 /dB');title('FIR 低通滤波器，汉

22、宁窗 ');figure;stem(n,b3);xlabel('n');ylabel('h(n)');grid;title('FIR 低通滤波器的单位脉冲响应 h(n) H,w=freqz(b3,1);矩形窗 ');汉宁窗 ');海明窗 ');H=20*log10(abs(H);figure;plot(w/pi,H);axis(O 1 -100 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB');title('FIR 低通滤波器，海明窗)；figure;ste

23、m( n,b4);xlabel(' n');ylabel('h( n)');grid;title('FIR 低通滤波器的单位脉冲响应 h(n)，布莱克曼窗');H,w=freqz(b4,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB');title('FIR低通滤波器，布莱克曼窗');5. 解：由于(2 tt/40) X 5<0.25 n因此下边界在 4和5之间

24、，上边界在 6和35之间。又因为N为偶数，属于第二类线性相位滤波器，幅度函数在02 n上以n为中心呈奇对称分布。因此频率采样值的幅度为：10kW40.5 k = 5Hk= 06兰 k 兰34-0.5 k =35-1 36 乞 k 39频响采样的相位为：, 2兀 N 139兀 _ = -kkN240频率采样值H(k)为:39 - fk e.39 -0.5e-f单位脉冲响应h(n)为:实验结果:10Hk-j即-0.5e"°-e1 Nj2 knh(nr：严以0-10-20-30-40-50-60-700.10.20.30.40.50.60.70.80.9-/ 二* -矩形窗汉宁窗

25、''A1 _ 1 1I Il严卜1LIIA乔i '| II1 1 1 11 J If 1 / 1!1频率采样法，N=40 ,过渡带设一个采样点程序：clear;N=40;k=0:N-1;Hk=ones(1,5) 0.5 zeros(1,29) -0.5 -ones(1,4);h仁 real(ifft(Hk.*exp(-j*pi*(N-1)*k/N);h2=ha nnin g(N)'.*h1;H1,w=freqz(h1,1);H2,w=freqz(h2,1);figure;plot(w/pi,20*log10(abs(H1),':r',w/pi,2

26、0*log10(abs(H2),'-b');axis(O 1-70 10); grid;xlabel('omega/pi'); ylabel('幅度 /dB');lege nd('矩形窗','汉宁窗');title('频率采样法，N=40,过渡带设一个采样点');6. 解：(1)由于(2 n31) X 10<2,3，因此下边界在 9和10之间，上边界在 11和20之间。 又因为N为奇数，高通，属于第一类线性相位滤波器，幅度函数在02 n上以n为中心呈偶对称分布。因此频率采样值的幅度为 ：00兰

27、k兰9,23兰k兰31叫=0.510,21111兰k兰20频响采样的相位为30 二31频率采样值H(k)为:=0.5e 310乞k空9,23空k乞31k = 10,2111乞k乞20单位脉冲响应h(n)为:1 NVj? knh(U(k)eN(2)由于(2 n32) X 10<2n3，因此下边界在 9和10之间，上边界在 11和20之间。又因为N为偶数，高通，属于第四类线性相位滤波器，幅度函数在02n上以n为中心呈偶对称分布。因此频率采样值的幅度为：00兰k兰9,23兰k兰31出二 0.510,21111兰k兰20频响采样的相位为频率采样值H(k)为:单位脉冲响应h(n)为:实验结果:0.

28、4300.5e3130 -:3120冬k乞9,23乞k乞31k 二 10,211仁k空201 Nj2 knh(nr：严以Ij)+ Q ,;r191: ?：l：0A线性相位高通滤波器的单位脉冲响应h(n), N=310.30.20.10-0.1-0.2-0.31015202530n频率采样法，N=31 ,过渡带设一个采样点-/ 二2:Q1i; 1匕J :i丁11；| O 7 占1.510.5x 10-15 线性相位高通滤波器的单位脉冲响应h(n), N=320-0.5-1-1.51025301520n35频率采样法，N=32 ,过渡带设一个采样点度 幅-/ 二程序：N=31;n=0:N-1;k=

29、0:N-1;Hk=zeros(1,10) 0.5 ones(1,10) 0.5 zeros(1,9);h1= real(ifft(Hk.*exp(-j*pi*(N-1)*k/N);h2=hammi ng(N)'.*h1;figure;stem( n,h1);xlabel(' n');ylabel('h( n)');grid;title('线性相位高通滤波器的单位脉冲响应h(n) , N=31');H1,w=freqz(h1,1);H2,w=freqz(h2,1);figure;plot(w/pi,20*log10(abs(H1),'

30、;:r',w/pi,20*log10(abs(H2),'-b');axis(0 1 -70 10); grid;legend（' 矩形窗 ',' 海明窗 '）;title（'频率采样法，N=31,过渡带设一个采样点'）;clear;N=32;n=0:N-1;k=0:N-1;Hk=zeros(1,10) 0.5 ones(1,11) 0.5 zeros(1,9); h1=ifft(Hk.*exp(-j*pi*(N-1)*k/N);h2=hamming(N)'.*h1;figure;stem(n,h1);xlabel(

31、'n');ylabel('h(n)');grid;title(' 线性相位高通滤波器的单位脉冲响应 h(n) ， N=32'); H1,w=freqz(h1,1);H2,w=freqz(h2,1); figure;plot(w/pi,20*log10(abs(H1),':r',w/pi,20*log10(abs(H2),'-b'); axis(0 1 -120 10); grid;xlabel('omega/pi'); ylabel('幅度 /dB');legend(' 矩形

32、窗 ',' 海明窗 ');title('频率采样法，N=32,过渡带设一个采样点');实验结果：第四类滤波器，多带线性相位滤波器的单位脉冲响应h(n)0.30.20.1-0.1-0.2-0.3-0.4L T1 r1 1'1|1101520253035404550n100-10-20-40-50-600.10.20.30.40.50.60.70.80.9-/二第四类滤波器，N=50 ,过渡带设一个采样点J1 I11A 1i ijnr1 11I Kn1A111ry-30-70第二类滤波器，多带线性相位滤波器的单位脉冲响应h(n)第二类滤波器，N=50 ,过渡带设一个采样点程序：clear;N=5