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

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是通带或阻带，2和4

3、是过渡带，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,要求阻带衰减不小于 40dB

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

7、.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出时域波形和幅频特性。理想滤波器的幅频特性四、实验结果与分析0.51解：求通带和阻带的数字边界频率0.7，求理想高通滤波器的边界频率。理想高通滤波器的边界频率是实际滤波器幅度函数正负 肩峰频率的中心，但这里不知道正负肩峰的频率，因此用通带和阻带边界频率的中心近似0.6选择窗函数和窗口长度。阻带衰减不小于50dB,因此选择海明窗（其阻带衰减为53dB,满足衰减要求）。滤波器的过渡带宽为0.2 ，因此窗口长度为660.2N 33线性相位延迟常数为16根据理想边界频率.和线性相位延迟常数，求理想单位脉冲响应m（n）:hd(n)円 n

9、ej(na)dej (n a)d sin( n a) sin( n a)Wn (n a)窗函数与理想单位脉冲响应相乘，得到线性相位FIR低通滤波器的单位脉冲响应h(n) w( n)hd( n)re c" / n、“sin (n a) sin 0.6 (n 16)1 / 、“0.54 0.46cos( )*Rn (n) n 1616(n 16)0.4n 16时域波形:线性相位高通滤波器FIR的单位脉冲响应h(n)0.40.30.20.1-0.1-0.2-0.3-0.4101520253035.17,1T rr r0幅频特性：0-10-20-30-40B 度-50 幅-60-70-80-

10、90-100线性相位高通滤波器FIR ,海明窗，N=3300.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)-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)');

11、H,w=freqz(h,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB');title('线性相位高通滤波器FIR，海明窗，N=33');2.解：求通带和阻带的数字边界频率：c1 0.3 , C2 0.7 , r1 0-2 , r2 0.9求理想带滤波器的边界频率。理想带通滤波器的边界频率是实际滤波器幅度函数正负肩峰频率的中心，但这里不知道正负肩峰的频率，因此用通带和阻带边界频率的中心近似c1120.25 , 2c

12、2 r20.8选择窗函数和窗口长度。阻带衰减不小于40dB,且阶数最小，因此选择汉宁窗（其阻带衰减为44dB,满足衰减要求）。滤波器的过渡带宽为c10.1r2c20.2，选择较小过渡带宽，因此窗口长度为线性相位延迟常数为6.20.2N 3115根据理想边界频率n和线性相位延迟常数求理想单位脉冲响应m(n):hd( n) 丁 £sin( n ) 2 sin( n (2 1)d)l窗函数与理想单位脉冲响应相乘，得到线性相位h(n) w(n )hd( n)时域波形:0.6幅频特性:2 j (n )V d FIR带通滤波器的单位脉冲响应0.5 0.5cos(n)sin°8 (n 型

13、 sin°.25 (n 理 “ 1515(n 15)0.55n 150.50.40.30.20.10-0.1-0.2-0.320253015n-叫一 12L11T线性相位带通滤波器FIR的单位脉冲响应h(n)线性相位带通滤波器FIR ,汉宁窗，N=310 C2:2" 7:CC、-10 -20f-30 -40 “B /度-5° -'|.r0.91幅-60/-70-80 -90-100 :-00.10.20.30.40.50.60.70.8频率/Hz程序：N=31;n=0:1:N-1;w1=0.25*pi;w2=0.8*pi;a=(N-1)/2;hd(a+1)

14、=0.55;win=ha nnin g(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（O 1 -100 0）;xlabel（'频率 /Hz'）;ylabel（' 幅度 /dB'）;title（'线性相位带通滤波器FIR，汉宁窗，N=3

15、1'）;3.解：求通带和阻带的数字边界频率：c10.3 , C2 0.7 , r1 0.35 , r2 0-6求理想带滤波器的边界频率。理想带阻滤波器的边界频率是实际滤波器幅度函数正负肩峰频率的中心，但这里不知道正负肩峰的频率，因此用通带和阻带边界频率的中心近似c1r120.325 , 2c2r220.65选择窗函数和窗口长度。阻带衰减不小于50dB,且阶数最小，因此选择海明窗（其阻带衰减为53dB,满足衰减要求）。滤波器的过渡带宽为1r10.05r2c2 0.1,选择较大过渡带宽，由于选择的是带阻滤波器,N必须为奇数，因此窗口长度为6.60.1N 67线性相位延迟常数为33根据理想边

16、界频率n和线性相位延迟常数，求理想单位脉冲响应h（n）:ej (n )d 12 j (n )1 j (n )hd(n) e d 一 e d2 1sin(n ) sin(n ) J-sin(n) ? n(n)窗函数与理想单位脉冲响应相乘，得到线性相位FIR带通滤波器的单位脉冲响应h(n) w(n 血(n)sin(n 33) sin( n 33)*0.35 卜sin(n 33)*0.65 “0.54 0.46cos(n)n 3333(n )0.7n 33时域波形:线性相位带阻滤波器FIR的单位脉冲响应h(n)0.70.60.50.40.30.20.1-0.1-0.210203040506070.

17、-Lt 1/1 J. II 1 ：宀 3 0.1L .燧也1；'r r ' r r'7 _频域特性:程序：N=67;n=0:1:N-1;w1=0.35*pi;w2=0.65*pi;a=(N-1)/2;(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;线性相位带阻滤波器FIR的单位脉冲响应h(n)');H,w=freqz(h,1);H=20*log10(abs(H);f

18、igure;plot(w/pi,H);axis(0 1 -90 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB');title('线性相位带阻滤波器FIR，海明窗，N=33');FIR低通滤波器，矩形窗4实验结果：0-10-J-1-1-20-30-40B度-50幅-60-70-80I1Iff-90-10000.10.20.30.40.50.60.70.80.91频率/HzFIR低通滤波器的单位脉冲响应h(n)，矩形窗1.20.80.60.40.2-0.2101520253035404550nJ. R1 1.1 J T I.

19、4 J11 11I '1FIR低通滤波器的单位脉冲响应h(n)，汉宁窗|FIR低通滤波器，汉宁窗-10-20-30-40-50-60-70-80-90-100 匚:00.10.20.30.40.50.60.70.80.91频率/HzFIR低通滤波器的单位脉冲响应h(n)，海明窗FIR低通滤波器，海明窗nilI, 11CCB度幅-390-O00O-260-8FIR低通滤波器的单位脉冲响应h(n),布莱克曼窗FIR低通滤波器，布莱克曼窗-10-20-30-40-50-60-70-80-9000.10.20.30.40.50.60.70.80.9频率/Hz1-100程序：N=51;M=N-1

20、;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',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 低通滤波器的

21、单位脉冲响应 h(n) H,w=freqz(b1,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel(' 频率 /Hz');ylabel(' 幅度 /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);

22、figure;plot(w/pi,H);axis(0 1 -100 0);xlabel(' 频率 /Hz');ylabel(' 幅度 /dB');title('FIR 低通滤波器，汉宁窗 ');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/

23、pi,H);axis(O 1 -100 0);xlabel('频率 /Hz');ylabel(' 幅度 /dB');title('FIR 低通滤波器，海明窗)；figure;stem( 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('频率 /H

24、z');ylabel(' 幅度 /dB');title('FIR低通滤波器，布莱克曼窗');5. 解：由于(2 tt/40) X 5<0.25 n因此下边界在4和5之间，上边界在6和35之间。又因为N为偶数，属于第二类线性相位滤波器，幅度函数在02 n上以n为中心呈奇对称分布。因此频率采样值的幅度为：0.5Hkk 34350.51 36k 39频响采样的相位为频率采样值H(k)为:单位脉冲响应h(n)为:实验结果:39-j40k e-j39 k0.5e0kk45Hk06k34-j即0.5e"°k3539 -j-40k e36k3

25、91 N 1j2 kn频率采样法，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,20*log10(abs(H2),'-b');axis(O 1-70 10); grid

26、;xlabel('omega/pi'); ylabel('幅度 /dB');legend('矩形窗','汉宁窗');title('频率采样法，N=40,过渡带设一个采样点');6. 解：(1)由于(2 n31) X 10<2,3，因此下边界在 9和10之间，上边界在 11和20之间。 又因为N为奇数，高通，属于第一类线性相位滤波器，幅度函数在02 n上以n为中心呈偶对称分布。因此频率采样值的幅度为：0 0k 9,23 k 31Hk0.5k 10,21111 k 20频响采样的相位为：,2N 130 ,kkkN

27、231频率采样值H(k)为：00 k 9,23 k 31Hk30-j3Tk0.5ek 10,21孔k e 3111 k 20单位脉冲响应h(n)为：1N 1j2 knh(n)H (k)e NNk 0(2)由于(2 n32) X 10<2*3，因此下边界在 9和10之间，上边界在 11和20之间。又因为N为偶数，高通，属于第四类线性相位滤波器，幅度函数在02n上以n为中心呈偶对称分布。因此频率采样值的幅度为 ：0 0k 9,23 k 31Hk0.5k 10,21111 k 20频响采样的相位为k频率采样值H(k)为：31200 k 9,23 k 310.5ek 10,21-j30 k-11

28、 k 20单位脉冲响应h(n)为:h(n)j2 knH (k)e N0实验结果:线性相位高通滤波器的单位脉冲响应h(n),N=310.40.30.20.10-0.1-0.2-0.310253015n)0 3¥1Ir ¥ GAAl:.20频率采样法，N=31 ,过渡带设一个采样点1G1 C爭$1 j2 一曲nr11 (眄十1.510.50#X 10-15 线性相位高通滤波器的单位脉冲响应h(n), N=32-0.5-1-1.51025301520n35频率采样法，N=32 ,过渡带设一个采样点程序：N=31;n=0:N-1;k=0:N-1;Hk=zeros(1,10) 0.5

29、 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),':r',w/pi,20*log10(a

30、bs(H2),'-b');axis(0 1 -70 10); grid;xlabel('omega/pi'); ylabel('幅度/dB');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)'.*

31、h1;figure;stem(n,h1);xlabel('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('

32、幅度 /dB');legend(' 矩形窗 ',' 海明窗 ');title('频率采样法，N=32,过渡带设一个采样点');7. 实验结果：第四类滤波器，多带线性相位滤波器的单位脉冲响应h(n)0.30.20.1-0.1-0.2-0.3-0.41015202530354550nOo1,j=a.E1I&冷Jt6%-：-1:4010第四类滤波器，N=50 ,过渡带设一个采样点Ir'l1IAii/1i !ThG I: :>1 II1FT11L Ar0-10-20-30-40-50-60-700.10.20.30.40.50.60.70.80.9/第二类滤波器，多带线性相位滤波器的单位脉冲响应h(n)度 幅第二类滤波器，N=50 ,过渡带设一个采样点程序：clear;N