实验6FIR滤波器设计

1、课程编号实验项目序号本科学生实验卡和实验报告信息科学与工程学院通信工程专业 2013级 1301 班课程名称：数字信号处理实验项目：FIR滤波器设计20152016学年 第 二 学期教育资料课 程 名 称数字信号处理实验课时4实 验 项 目 名 称 和 编 号FIR滤波器设计同组者姓 名实 验 目 的学会运用MATLA产生常用离散时间信号实 验 环 境MATLAB实 验 内 容 和 原 理（1） 认真复习FIR数字滤波器的基本概念，线性相位 FIR滤波器的条件和特点、幅 度函数特点、零点位置的基本特点与性质；窗函数设计法的基本概念与方法，各种窗函数的性能和设计步骤，线性相位 FIR低通、高通、

2、带通和带阻滤波器的设计方法， 频率米样设计法的基本概念和线性相位的实现方法。（2）掌握几种线性相位的特点，熟悉和掌握矩形窗、三角形窗、汉宁窗、海明窗、布 莱克曼窗、凯塞窗设计IIR数字滤波器的方法，熟悉和掌握频率抽样设计法的线性 相位的设计方法，并对各种线性相位的频率抽样法的设计给出调整和改进。（3）熟悉利用MATLAB进行各类FIR数字滤波器的设计方法。学号：201308030104_姓名：王少丹.专业年级班级： 通信1301四合院实验室 组别实验日期_2016年_6_ 月_12_ 日a. 设线性相位FIR滤波器单位抽样响应分别为h(n )= -4,1, -1，- 2,5,6,5, -2,

3、-1,1, -4h(n )= -4,1, -1，- 2,5,6,6,5,- 2, -1,1，- 4h(n )= -4,1，-，-2,5,0,- 5,2,1，- 1,4h(n )= -4,1, -1, -2,5,6,- 6, -5,2,1, -1,4分别求出滤波器的幅度频率响应H( 3),系统函数H(z)以及零极点分布，并绘制相应的波形和分布图。-0, 98070.19561axis实验步骤方法关 键 代 码D. 4052(1 40521.2169(k 82180.23500. 83OOi厲 830011.237411* 2374i0. 00001厲 OOOOi0.7299i-H 9807-0.

4、55780.23900. 821S1” 21690.23900.7299i0.196610. 830Oi0.830010. 72&9iCL 72&9iCk 000010. OOQOi宅 j j 口SWE)KAffl ZJKTj 舸 *： *KT hiType2:=1,3120 -1.C000 -0.3363 -0.SS6S0.47730.4778 -0-2957 -0.29570.76220.29260.2926十 0,00001+ 0.00001十 0.46221-0.46221+ 1.18511-1.1851L+ 0.96631-0. 95531+ 0.00001+ 0.

5、72581-0.72531ans 0. 7622 + 0. OOOOi -1. 0000 十 0. OOdOi -0. 336S + 0. 422i -CL X868 - 0. 062210. 2926 + 0. 7258i0. 2926 - 0. 725Si -0. 2957 + 0. 95531 -0. 2&57 - 0. 95531Type3:Type4:-E 0000+0+ 00001-D. 8732+0. 4E741-h 87320. 487410. 101C+1.204110. 1010L+ 2041 i1. 0000+0. 000010. 8280+0+ 560710.

6、82300. 5607lQ.Q6&2+0*S247i0. 0695一0. 82471ans-1,0000 + 0.00001-0-3732 + 0_ 48741-0.8732 - 0. 467410.0692 + 0_ 824710. 0692 - 0.S24711.0000 + 0.OOOOi0+8230 + 0.560710- 3280 - 0. 56071figure 6加n JB=i.：E 共伯 fix 7）二flOil： 奇二冋 iEBiHn.l* 很 S.llR11:oJ-1- p*Is d10-til-hici iD 尺抽| F?nXZ =0.2631 十 1. 3394

7、10.2631 - 1. 3394i-0. 9505 十 0, 31081-0. 9506 - 0.310Bi-a 7309 + 0- G825i-0. 7309 - 0,632510. 9021 + 0.431510.9021 - 0_ 43161t. aoco + OOOOQi0.1412 + 0.7LB9i0:1412 - 0.7159ians =0.1412 十 0,71S9i0.1412 - 0.71S91-0. 9505 + 0. 31081-0. 9505 - 0. 31091-0.7309 + 0+ 6S25i-0. 7305 - 0. 682510.9021 + 0,4315

8、i0. 9021 - 0. 43151b. 设计FIR 数字低通滤波器，技术指标为：p=0.2 n,3 st=0.3 n,S 1=0.25dB, S 2=50dB(1)通过技术指标，选择一种窗函数进行设计；求滤波器的单位抽样响应、频率响应，并绘制波形。(3)选择凯塞窗函数设计该滤波器，并绘制相应的波形图。女綽何询辑(曰人jiscn 卓面socw)律助g士冈 口?11蘇茴幅鹿口向应(dE3)C5U2以伸为单位的頻率0 20 250 3決伽为卑位的频率 相位响应Q Tc. 设计FIR数字带通滤波器，技术指标为：下阻带边缘：3 st1=0.2 n,S s1=60dB,下通带边缘：3 p1=0.35

9、n,S p仁 1dB; 上 通带边缘：p2=0.65 n,S p仁1dB,上阻带边缘：st2=0.8 n,S s2=60dB;(1)通过技术指标，选择一种窗函数进行设计；求滤波器的单位抽样响应、频率响应，并绘制波形。d. 设计FIR数字带通滤波器，技术指标为：下阻带边缘：3 st1=0.2 n,S s1=60dB,下通带边缘：3 p1=0.4 n,S p仁1dB; 上通 带边缘： p2=0.6 n,S p仁 1dB,上阻带边缘： st2=0.8 n,S s2=60dB;(1) 通过技术指标，选择一种窗函数进行设计；(2) 求滤波器的单位抽样响应、频率响应，并绘制波形。e. 设计FIR数字带通滤

10、波器，技术指标为：下阻带边缘：3 st1=0.2 n,S s1=20dB,下通带边缘：3 p1=0.4 n,S p仁1dB; 上通 带边缘： p2=0.6 n,S p仁 1dB,上阻带边缘： st2=0.8 n,S s2=20dB;(1)通过技术指标，选择一种窗函数进行设计；求滤波器的单位抽样响应、频率响应，并绘制波形。(1)通过技术指标，选择一种窗函数进行设计； 求滤波器的单位抽样响应、频率响应，并绘制波形h.滤波器的技术指标为：通带截止频率为3p=0.6 n,阻带截止频率为3 st=0.4 n,通带最大衰减为S 1=0.25dB，阻带最小衰减为S 2=40dB(1)通过技术指标，选择一种窗

11、函数设计一个具有n /2相移的FIR高通滤波器；求滤波器的单位抽样响应、频率响应，并绘制波形。i.设计FIR数字带阻滤波器，其技术指标为：低端阻带边缘：3 st1=0.4 n, S s仁40dB,低端通带边缘：3 p1=0.2 n, S p仁1dB;高 端通带边缘：3 p2=0.8 n,S p仁1dB,高端阻带边缘：3 st2=0.6 n,S s2=40dB;(1)通过技术指标，选择一种窗函数进行设计；求滤波器的单位抽样响应、频率响应，并绘制波形。V23 24FIR滤波器的单位抽样响应为 h(n)=1/9亠，编制MATLAB程序求系统的频率采样型结构的系数，并画出频率抽样型结构( < w

12、 n) (n <«/<1一伽m. 一个理想差分器的频率响应为:用长度为21的汉宁窗设计一个数字FIR差分器，并绘制其时域和频率的响应波形Figure 1交性旧扁事氓去蚩M 摘人曲二冥(T) 卓面(D) 奇口 fW)Idaal Impulse Response-u.t)-1-p=% 甥®辺諾皀叮 口Hamming Window0 S010Actual Impulse0 6irequency in pi units10n亘1520S-Fn 2亠盟口肓刍 O-6G 0-一二CI O Q ? f>? ? Q Ocbe-to5G立舍2Ci °115Amph

13、tLida Responn.利用汉宁窗设计一个长度为25的数字希尔伯特变换器，并绘制它的时域和频域的 响应波形。p. FIR 数字低通滤波器的技术指标为：p=0.2 n,3 st=0.3 n,S 1=0.25dB,S2=50dB利用频率采样方法设计FIR数字滤波器，并绘制滤波器的单位冲激响应、幅度频率响应的波形。q.用窗函数法设计一个线性相位的FIR数字低通滤波器，其技术指标为：p=0.2n,3 st=0.4 n,S 1=0.25dB,S 2=50dB求滤波器的单位抽样响应、频率响应，并绘制波形。(2)对该滤波器输入一个宽度为10的矩形序列，求滤波器的输出信号，并绘制相应 的波形图。测 试 记

14、 录 分 析 结 论理论与实验值相符。(后附代码)小 结认真复习FIR数字滤波器的基本概念，线性相位 FIR滤波器的条件和特点、幅 度函数特点、零点位置的基本特点与性质；窗函数设计法的基本概念与方法，各种窗 函数的性能和设计步骤，线性相位 FIR低通、高通、带通和带阻滤波器的设计方法， 频率米样设计法的基本概念和线性相位的实现方法。以下由实验教师填写记 事 评 议成 绩 评 疋平时成绩实验报告成绩综合成绩指导教师签名：代码:wp = 0.2*pi;ws = 0.4*pi;tr_width = ws-wp;M = ceil(6.6*pi/tr_width)+1;n = 0:1:M-1;wc =

15、(ws+wp)/2;hd = idealp(wc,M);w_ham = (hammi ng(M)'h = hd.*w_ham;db,mag,pha,H,w = freqz_m3(h,1);delta_w = 2*pi/1000;Rp = -(mi n( db(1:1:wp/delta_w+1);As = -rou nd(max(db(ws/delta_w+1:1:501);figure(1)subplot(221)stem (n ,hd);title('Ideal Impulse Rresp on se')axis(0 M-1 -0.1 0.3);xlabel('

16、;n' );ylabel('hd(n)')subplot(222)stem(n,w_ham);title( 'Hamming Window' )axis(0 M-1 0 1.1);xlabel( 'n' );ylabel( 'w(n)' ) subplot(223)stem(n,h);title( 'Actual Impuse Response' ) axis(0 M-1 -0.1 0.3);xlabel( 'n' );ylabel( 'h(n)' ) subplot(224

17、)plot(w/pi,db);title( 'Magnitude Response in db' );gridaxis(0 1 -100 10);xlabel( 'frequence in pi unit' );ylabel( 'decibels' ) %solution of problem (b)Rn = stepseq(0,0,9)n = 0:1:9 dw = w(2)-(1) w = -fliplr(w),w(2:501) H = fliplr(H),H(2:501)xw = Rn *(exp(-j).A( n'*w)answ =

18、 H.*xwRnn = an sw*exp(-j).A(w'* n)*dw/(2*pi)figure(2)subplot(221)stem(n,Rn)xlabel( 'n' );ylabel( 'Rn' );title( 'input rectangle signal in T-domain')subplot(222)stem(w/pi,xw)xlabel( 'w/pi' );ylabel( 'Rw' );title( 'input rectangle signal in F-domain'

19、)subplot(223)plot(w/pi,abs(answ)xlabel( 'w/pi' );ylabel( 'Magnitude of H(w)' );title( 'output signal in F-domain after filteration' ) subplot(224) stem(n,abs(Rnn) xlabel( 'n' );ylabel( 'Rnn' );title( 'output signal in T-domain after filteration'M = 25;

20、 alpha = (M-1)/2; n = 0:M-1;hd = (2/pi)*(sin(pi/2)*(n-alpha).A2)./(n-alpha); hd(alpha+1)=0;w_han = (hann(M)' h = hd .* w_han; Hr,w,P,L = Hr_Type3(h);subplot(2,2,1); stem(n,hd); title('Ideal Impulse Response' )axis(-1 M -1.2 1.2); xlabel('n' ); ylabel( 'hd(n)');subplot(2,2

21、,2); stem(n,w_han);title( axis(-1 M 0 1.2); xlabel('n'subplot(2,2,3); stem(n,h);title( axis(-1 M -1.2 1.2); xlabel( w = w' Hr = Hr''Hann Window' ); ylabel( 'w(n)' )'Actual Impulse Response')'n' ); ylabel( 'h(n)' )w = -fliplr(w), w(2:501); Hr =

22、-fliplr(Hr), Hr(2:501); subplot(2,2,4);plot(w/pi,Hr); title('Amplitude Response' );grid;xlabel( 'frequency in pi units'); ylabel( 'Hr' ); axis(-1 1 -1.1 1.1);M = 20; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l;Hrs = 1,1,1,zeros(1,15),1,1;%Ideal Amp Res sampledHdr = 1,1,0,0; w

23、dl = 0,0.25,0.25,1;%Ideal Amp Res for plottingk1 = 0:floor(M-1)/2); k2 = floor(M-1)/2)+1:M-1;angH = -alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2);H = Hrs.*exp(j*angH); h = real(ifft(H,M);db,mag,pha,grd,w = freqz_m(h,1); Hr,ww,a,L = Hr_Type2(h); subplot(2,2,1);plot(wl(1:11)/pi,Hrs(1:11), 'o' ,wdl

24、,Hdr);axis(0,1,-0.1,1.1); title( 'Frequency Samples: M=20' ) xlabel( 'frequency in pi units'); ylabel( 'Hr(k)' );subplot(2,2,2); stem(l,h); axis(-1,M,-0.1,0.3)title( 'Impulse Response' ); xlabel( 'n' ); ylabel( 'h(n)' ); subplot(2,2,3); plot(ww/pi,Hr,w

25、l(1:11)/pi,Hrs(1:11), 'o' ); axis(0,1,-0.2,1.2); title('Amplitude Response' )xlabel( 'frequency in pi units'); ylabel( 'Hr(w)' )subplot(2,2,4);plot(w/pi,db); axis(0,1,-60,10); gridtitle( 'Magnitude Response' ); xlabel( 'frequency in pi units');ylabel(

26、'Decibels' );Wpl=0.2*pi;Wph=0.8*pi;Wsl=0.4*pi;Wsh=0.6*pi;tr_width=min(Wsl-Wpl),(Wph-Wsh);M=ceil(6.2*pi/tr_width) n=0:1:M-1;Wcl=(Wsl+Wpl)/2;Wch=(Wsh+Wph)/2; hd=ideal_bs(Wcl,Wch,M); w_ham=(hanning(M)' h=hd.*w_ham;db,mag,pha,w=freqz_m2(h,1);delta_w=2*pi/1000;Ap=-(min(db(1:1:Wpl/delta_w+1)As

27、=-round(max(db(Wsl/delta_w+1:1:Wsh/delta_w+1) subplot(221)stem(n,hd);title( 'Ideal Impulse Rresponse' ) axis(0 M-1 -0.1 0.7);xlabel( 'n' );ylabel( 'hd(n)' ) subplot(222)stem(n,w_ham);title( 'Hamming Window' )axis(0 M-1 0 1.1);xlabel( 'n' );ylabel( 'w(n)

28、9; ) subplot(223)stem(n,h);title( 'Actual Impuse Response' )axis(0 M-1 -0.1 0.7);xlabel( 'n' );ylabel( 'h(n)' ) subplot(224)plot(w/pi,db);title('Magnitude Response in db');gridaxis(0 1 -100 10);xlabel('frequence in pi unit');ylabel( 'decibels' )clear a

29、llfigureWpl=0.2*pi;Wph=0.8*pi;Wsl=0.4*pi;Wsh=0.6*pi;tr_width=min(Wsl-Wpl),(Wph-Wsh);M=ceil(6.2*pi/tr_width)n=0:1:M-1;Wcl=(Wsl+Wpl)/2;Wch=(Wsh+Wph)/2;hd=ideal_bs(Wcl,Wch,M); w_ham=(hanning(M)' h=hd.*w_ham;db,mag,pha,w=freqz_m2(h,1);delta_w=2*pi/1000;Ap=-(min(db(1:1:Wpl/delta_w+1)As=-round(max(db(

30、Wsl/delta_w+1:1:Wsh/delta_w+1) subplot(221)stem(n,hd);title('Ideal Impulse Rresponse')axis(0 M-1 -0.1 0.7);xlabel('n' );ylabel( 'hd(n)' )subplot(222)stem(n,w_ham);title( 'Hamming Window' )axis(0 M-1 0 1.1);xlabel('n' );ylabel( 'w(n)' )subplot(223)stem(

31、n,h);title('Actual Impuse Response')axis(0 M-1 -0.1 0.7);xlabel('n' );ylabel( 'h(n)' )subplot(224)plot(w/pi,db);title( 'Magnitude Response in db' );gridaxis(0 1 -100 10);xlabel('frequence in pi unit');ylabel( 'decibels' );figure ws1=0.2*pi; wp1=0.35*pi

32、;ws2=0.8*pi; wp2=0.65*pi; Ap=60;Rp=1;tr_width=min(wp1-ws1),(ws2-wp2);M=ceil(11*pi/tr_width);n=0:1:M-1;wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);w_bla=(blackman(M)'h=hd.*w_bla;Qi ?3?3e ?R ''2? a 3?u');'?e ?3?3e ?d'H,W=freqz(h,1); subplot(2,2,1);stem(n,

33、hd);title( subplot(2,2,2);stem(n,w_bla);title( subplot(2,2,3);stem(n,h);title(-u ?e ?io | £Jdb£?'subplot(2,2,4);plot(W/pi,20*log10(abs(H);title( clear allfigurews1=0.2*pi; wp1=0.4*pi; ws2=0.8*pi; wp2=0.6*pi;Ap=60;Rp=1;tr_width=min(wp1-ws1),(ws2-wp2);M=ceil(11*pi/tr_width);n=0:1:M-1;wc1

34、=(ws1+wp1)/2;wc2=(wp2+ws2)/2;hd=ideal_lp(wc2,M)-ideal_lp(wc1,M); w_bla=( hamming(M)' h=hd.*w_bla;Qi ?3?3e ?R ''o £?+，° 'H,W=freqz(h,1); subplot(2,2,1);stem(n,hd);title( subplot(2,2,2);stem(n,w_bla);title(subplot(2,2,3);stem(n,h);title('?e ?3?3e ?d'subplot(2,2,4);plo

35、t(W/pi,20*log10(abs(H);title(' u ?e ?io | £")db£?'clear allfigurews1=0.2*pi;wp1=0.4*pi;ws2=0.8*pi;wp2=0.6*pi;Ap=20;Rp=1;tr_width=min(wp1-ws1),(ws2-wp2);M=ceil(11*pi/tr_width);n=0:1:M-1;wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;hd=ideal_lp(wc2,M)-ideal_lp(wc1,M); w_bla=( boxcar(M)'h=

36、hd.*w_bla;Qi ?3?3e'?D? '° '?e ?3?3eH,W=freqz(h,1);subplot(2,2,1);stem(n,hd);title( subplot(2,2,2);stem(n,w_bla);title(-u ?e ?io | £Jdb£?'subplot(2,2,3);stem(n,h);title( subplot(2,2,4);plot(W/pi,20*log10(abs(H);title( clear allfigureAs=55;ws=11*pi/27;wp=15*pi/27;tr_widt

37、h=wp-ws;%?1ye'?M=ceil(As-7.95)*2*pi/(14.36*tr_width)+1)+1;% ' ?-? o'° ?2？+卩?3O ?edisp( '?2 - ?+卩?3O?i<um2str(M);beta=0.1102*(As-8.7);%?-?o'° ?beta?卩i=0:1:M-1;disp( '?D? a ?2?-,hum2str(beta);w_kai=(kaiser(M,beta)'%?o ?-? o'° o_eywc=(ws+wp)/2;hd=ideal_

38、lp(pi,M)-ideal_lp(wc,M);%?6ai ?3?io|h=hd.*w_kai;db,mag,pha,grd,w=freqz_m(h,1);delta_w=2*pi/1000;Rp=-(mii(db(wp/delta_w+1:1:501);disp( '?e？2?Un?astr(Rp);As=-rouid(max(db(1:1:ws/delta_w+1);disp( ' x ?D?xe' ?¥ ?a?m2str(As);subplot(2,2,1);stem(i,hd);title( ' ai ?3?io);|'axis(0 M

39、-1 -0.4 0.8);ylabel( 'hd(i)');subplot(2,2,2);stem(i,w_kai);title( '?-? 6'°);'axis(0 M-1 0 1.1);ylabel( 'wd(i)');subplot(2,2,3);stem(n,h);title( '?e ?3?i 笔 |'axis(0 M-1 -0.4 0.8);xlabel( 'n' );ylabel( 'h(n)');subplot(2,2,4);plot(w/pi,db);title(

40、 '? u ?e ?io |/dB'axis(0 1 -100 10);grid;xlabel( ' o ?pi?a 卩丫 ? ? )? e'ylabel( '? ?±'e y/QB'clear allfigure%?u matlab?DD?o_e yg£?'u ?e ?£oWp=0.6*pi;Ws=0.4*pi;tr_width=Wp-Ws;M=ceil(6.2*pi/tr_width);n=0:1:M-1;Wc=(Ws+Wp)/2;hd=ideal_lp(pi,M)-ideal_lp(Wc,M);

41、 w_ham=(hanning(M)'h=hd.*w_ham;db,mag,pha,w=freqz_m2(h,1);delta_w=2*pi/1000;Ap=-(min(db(Wp/delta_w+1:1:501) As=-round(max(db(1:1:Ws/delta_w+1) subplot(221)stem(n,hd);title( 'Ideal Impulse Rresponse' ) axis(0 M-1 -0.4 0.5);xlabel('n' );ylabel( 'hd(n)')subplot(222)stem(n,w_

42、ham);title( 'Hamming Window' )axis(0 M-1 0 1.1);xlabel('n' );ylabel( 'w(n)')subplot(223)stem(n,h);title( 'Actual Impuse Response' ) axis(0 M-1 -0.4 0.5);xlabel( 'n' );ylabel( 'h(n)' ) subplot(224)plot(w/pi,db);title( 'Magnitude Response in db'

43、);gridaxis(0 1 -100 10);xlabel( 'frequence in pi unit' );ylabel( 'decibels' ) clear allfigure Wp=0.6*pi; Ws=0.4*pi; tr_width=Wp-Ws;M=ceil(6.2*pi/tr_width) n=0:1:M-1;Wc=(Ws+Wp)/2; hd=ideal_lp(pi,M)-ideal_lp(Wc,M);w_ham=(hanning(M); h=hd.*rot90(w_ham);db,mag,pha,w=freqz_m2(h,1); delta_

44、w=2*pi/1000;Ap=-(min(db(Wp/delta_w+1:1:501) As=-round(max(db(1:1:Ws/delta_w+1) figuresubplot(221)stem(n,hd);title( 'Ideal Impulse Rresponse' ) axis(0 M-1 -0.4 0.5);xlabel( 'n' );ylabel( 'hd(n)' ) subplot(222)stem(n,w_ham);title( 'Hamming Window' ) axis(0 M-1 0 1.1);xl

45、abel( 'n' );ylabel( 'w(n)' ) subplot(223)stem(n,h);title( 'Actual Impuse Response' ) axis(0 M-1 -0.4 0.5);xlabel( 'n' );ylabel( 'h(n)' ) subplot(224)plot(w/pi,db);title( 'Magnitude Response in db' );gridaxis(0 1 -100 10);xlabel( 'frequence in pi un

46、it' );ylabel( 'decibels' )%o ?Hammi ng ° o_e yee ?FIRe yX ?2？*figure wp=0.2*pi;ws=0.3*piN=61n=0:1:N-1wc=(ws+wp)/2;?小？2？*hd=idealp(wc,N);?卩丨厂卩？3?O ? io|w_ham=(hamming(N)'h=hd.*w_ham; %FIR?2' ?-3?O? io| db,mag,pha,grd,w=freqz_m(h,1); delta_w=2*pi/1000;Rp=-(min(db(1:1:wp/delta_

47、w+1)%e ?¥ ?As=-round(max(db(ws/delta_w+1:1:501)%e ?x ?D?xe' ?¥ ?subplot(221);stem(n,hd);title('ai ?3?O ? i 0 |'axis(0 N-1 -0.1 0.3);xlabel('n');ylabel( 'hd(n)')subplot(222);stem(n,w_ham);title('hamming'° 'axis(0 N-1 0 1.1);xlabel('n' );y

48、label( 'w(n)' );subplot(223);stem(n,h);title(.?e 3?O? 1)0 |'axis(0 N-1 -0.1 0.3);xlabel('n');ylabel( 'h(n)')subplot(224);plot(w/pi,db);axis(0 0.8 -100 0);xlabel('0 ?PI?a 卩¥ ?卩? y)?eabel('? e y u ?e /Qb'clear allfigure%o ?Kaiser '° o_e yee ?FIRe y

49、X ?2?- wp=0.2*pi;ws=0.3*pi;As=50tr_width=ws-wpN=ceil(As-7.95)/(14.36*tr_width/(2*pi)+1)+1n=0:1:N-1beta=0.1102*(As-8.7)wc=(wp+ws)/2 %di ?卩丨厂卩？1?卩?ehd=ideal_lp(wc,N)w_kai=(kaiser(N,beta)'h=hd.*w_kaidb,mag,pha,grd,w=freqz_m(h,1) delta_w=2*pi/1000Rp=-(min(db(1:1:wp/delta_w+1)%?e ?e ?i ?¥?As=-ro

50、und(max(db(ws/delta_w+1:1:501)%e ?e ?X ?D?xe' ?¥?subplot(211);plot(w/pi,db);title('?- e -' °u ?e ?16】叫axis(0 0.5 -100 0)ylabel( '? e y u ?e /dfcXlabel( ' o ?/pi?a 严？卩? )?e' subplot(212);plot(w/pi,pha);title('? a ? io )：gridaxis(0 0.5 -4 4)ylabel( '? a ?'

51、);xlabel( ' o ?/pi?a 卩丫 ? ? )? e'figureh1=-4,1,-1,-2,5,6,5,-2,-1,1,-4;M=length(h1);n=0:M-1;A,w,type,tao=amplres(h1);type subplot(2,1,1),stem(n,h1);title( '3? O?io |h)ylabel( 'h(n)' );xlabel( 'n' );subplot(2,1,2),plot(w/pi,A);ylabel( 'A' );xlabel( 'pi');tit

52、le( '? u ?卩? i 0 |'figurerz=roots(h1)for i=1:8r(i)=1/rz(i);endr'zplane(h1,1);title( 'h1 a ?卩? ?' clear allfigureh2=-4,1,-1,-2,5,6,6,5,-2,-1,1,-4; M=length(h2);n=0:M-1;A,w,type,tao=amplres(h2);type subplot(2,1,1),stem(n,h2);title( '3? O?io |h2； ylabel( 'h(n)' );xlabel(

53、 'n' ); subplot(2,1,2),plot(w/pi,A); ylabel( 'A' );xlabel( 'pi');title( ' u ?卩？ i )0 I' figure rz=roots(h2) for i=1:8 r(i)=1/rz(i);end r' zplane(h2,1);title( '?h2 a ?卩？ i ?' clear allfigure h3=-4,1,-1,-2,5,0,-5,2,1,-1,4;M=length(h3);n=0:M-1; A,w,type,tao=amplres(h3);type subplot(2,1,1),stem(n,h3);title( '3? O?io |h3； ylabel( 'h(n)' );xlabel( 'n' ); subplot(2,1,2),plot(w/pi,A); ylabel( 'A' );xlabel(