实验6FIR滤波器设计

《实验6FIR滤波器设计》由会员分享，可在线阅读，更多相关《实验6FIR滤波器设计（17页珍藏版）》请在装配图网上搜索。

1、课程编号实验项目序号本科学生实验卡和实验报告信息科学与工程学院通信工程专业 2013级 1301 班课程名称：数字信号处理实验项目：FIR滤波器设计20152016学年 第 二 学期学号：0104 姓名：王少丹专业年级班级： 通信1301四合院实验室 组别 实验日期2016 年 6月_12 日课 程 名 称数字信号处理实验课时4实 验 项 目 名 称 和 编 号FIR滤波器设计同组者姓 名实 验 目 的学会运用MATLA产生常用离散时间信号实 验 环 境MATLAB(1)认真复习FIR数字滤波器的基本概念，线性相位FIR滤波器的条件和特实点、幅度函数特点、零点位置的基本特点与性质；窗函数设计法

2、的基本概念与方验法，各种窗函数的性能和设计步骤，线性相位FIR低通、高通、带通和带阻滤内波器的设计方法，频率米样设计法的基本概念和线性相位的实现方法。容(2)掌握几种线性相位的特点，熟悉和掌握矩形窗、三角形窗、汉宁窗、海明和窗、布莱克曼窗、凯塞窗设计IIR 数字滤波器的方法，熟悉和掌握频率抽样原设计法的线性相位的设计方法，并对各种线性相位的频率抽样法的设计给出调整理和改进。(3)熟悉利用MATLAB进行各类FIR数字滤波器的设计方法。a.设线性相位FIR滤波器单位抽样响应分别为h(n )= -4,1, -1，- 2,5,6,5, -2, -1,1, -4h(n )= -4,1, -1，- 2,

3、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)以及零极点分布，并绘制相应的波形和分布图。实Type2:验Type3:步Type4:骤b.设计FIR数字低通滤波器，技术指标为：3 p=0.2n，3 st=0.3n, S 1=0.25dB，方S 2=50dB。法(1)通过技术指标，选择一种窗函数进行设计；关求滤波器的单位抽样响应、频率响应，并绘制波形。键(3)选择凯塞窗函数设计该滤波器，并绘制相应

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

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

6、B,阻带最小衰减为S 2=55dB。 (1)通过技术指标，选择一种窗函数进行设计；求滤波器的单位抽样响应、频率响应，并绘制波形。系统一：g. 设计FIR数字高通滤波器，技术指标为：通带截止频率为3p=0.6n，阻带截止频率为3 st=0.4n，通带最大衰减为S 1=0.25dB，阻带最小衰减为S 2=40dB。(1)通过技术指标，选择一种窗函数进行设计；求滤波器的单位抽样响应、频率响应，并绘制波形。h. 滤波器的技术指标为：通带截止频率为3 p=0.6n，阻带截止频率为3 st=0.4n, 通带最大衰减为S 1=0.25dB，阻带最小衰减为S 2=40dB。(1)通过技术指标，选择一种窗函数设

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

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

9、为 10的矩形序列，求滤波器的输出信号，并绘制 相应的波形图。测 试 记 录 分 析 结 论理论与实验值相符。(后附代码)小 结认真复习FIR数字滤波器的基本概念，线性相位FIR滤波器的条件和特点、幅度函数特点、零点位置的基本特点与性质；窗函数设计法的基本概念与方 法，各种窗函数的性能和设计步骤，线性相位FIR低通、高通、带通和带阻滤波器的设计方法，频率采样设计法的基本概念和线性相位的实现方法。以下由实验教师填写记 事 评 议成 绩 评 疋平时成绩实验报告成绩综合成绩指导教师签名：代码:wp = 0.2*pi;ws = 0.4*pi; tr_width = ws-wp;M = ceil(6.6

10、*pi/tr_width)+1;n = 0:1:M-1;wc = (ws+wp)/2;hd = ideal_lp(wc,M); w_ham = (hamming(M); h = hd.*w_ham;db,mag,pha,H,w = freqz_m3(h,1);delta_w = 2*pi/1000;Rp = -(min(db(1:1:wp/delta_w+1);As = -round(max(db(ws/delta_w+1:1:501); figure(1)subplot(221)stem(n,hd);title( Ideal Impulse Rresponse ) axis(0 M-1 -0

11、.1 0.3);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(n,h);title( Actual Impuse Response ) axis(0 M-1 -0.1 0.3);xlabel( n );ylabel( h(n) ) subplot(224)decibels )plot(w/pi,db);title( Magnitude Response in db

12、 );grid axis(0 1 -100 10);xlabel( frequence in pi unit );ylabel( %solution of problem (b)Rn = stepseq(0,0,9)n = 0:1:9dw = w(2)-(1)w = -fliplr(w),w(2:501)H = fliplr(H),H(2:501)xw = Rn *(exp(-j).A( n*w)answ = H.*xwRnn = an sw*exp(-j).A(w* n)*dw/(2*pi)figure(2) subplot(221)stem(n,Rn)xlabel( n );ylabel(

13、Rn);title(input rectangle signal in T-domainsubplot(222)stem(w/pi,xw)xlabel( w/pi);ylabel(Rw );title(input rectangle signal in F-domainsubplot(223)plot(w/pi,abs(answ)xlabel( w/piafter filteration);ylabel()Magnitude of H(w);title(output signal in F-domainsubplot(224)stem(n,abs(Rnn)xlabel( n );ylabel(

14、Rnn);title(output signal in T-domain after filterationM = 25; alpha = (M-1)/2; n = 0:M-1;hd = (2/pi)*(s in (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);Ideal Impulse Responsesubplot(2,2,1); stem(n,hd); title(axis(-1 M -1.2 1.2); xlabel( su

15、bplot(2,2,2); stem(n,w_han);title(n ); ylabel( hd(n);Hann Window )axis(-1 M 0 1.2); xlabel(n ); ylabel(w(n)subplot(2,2,3); stem(n,h);title(Actual Impulse Responseaxis(-1 M -1.2 1.2); xlabel(n ); ylabel(h(n)w = w; Hr = Hr;w = -fliplr(w), w(2:501); Hr = -fliplr(Hr), Hr(2:501);subplot(2,2,4);plot(w/pi,

16、Hr); title(Amplitude Response);grid;xlabel( frequency in pi units); ylabel(Hr ); axis(-1 1 -1.11.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; wdl = 0,0.25,0.25,1;%Ideal Amp Res for plottingk1 = 0:floor(M-1)/2); k2 = floor(M-1

17、)/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,Hdr);Frequency Samples: M=20axis(0,1,-0.1,1.1); title(xlabel( frequency in pi units); ylabel(

18、subplot(2,2,2); stem(l,h); axis(-1,M,-0.1,0.3)Hr(k);title( Impulse Response); xlabel(n ); ylabel(h(n);o );)subplot(2,2,3); plot(ww/pi,Hr,wl(1:11)/pi,Hrs(1:11),axis(0,1,-0.2,1.2); title(Amplitude Response6文档来源为 :从网络收集整理xlabel( frequency in pi units ); ylabel( Hr(w) ) subplot(2,2,4);plot(w/pi,db); axi

19、s(0,1,-60,10); grid);title( Magnitude Response ); xlabel(frequency in pi unitsylabel( 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,m

20、ag,pha,w=freqz_m2(h,1);delta_w=2*pi/1000; Ap=-(min(db(1:1:Wpl/delta_w+1)As=-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 );ylab

21、el( w(n) )subplot(223)stem(n,h);title(Actual Impuse Response)axis(0 M-1 -0.1 0.7);xlabel(n );ylabel( h(n) )subplot(224)decibels );plot(w/pi,db);title( Magnitude Response in db );grid axis(0 1 -100 10);xlabel( frequence in pi unit );ylabel( clear allfigureWpl=0.2*pi;Wph=0.8*pi; Wsl=0.4*pi;Wsh=0.6*pi;

22、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(Wsl/delta_w+1:1:Wsh/delta_w+1) subplot(221)stem(n,hd);t

23、itle(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(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( Magn

24、itude Response in db );griddecibels );axis(0 1 -100 10);xlabel( frequence in pi unit );ylabel( figurews1=0.2*pi;wp1=0.35*pi;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=(bl

25、ackman(M);h=hd.*w_bla;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( 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*p

26、i/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=( hamming(M);h=hd.*w_bla;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( subplot(2,2,4);plot(W/pi,20*log10(abs(H);title( clear allfigur

27、ews1=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);Qi ?3?3 e ?u)；2? a 3? u， )；?e ?3?3 e ? v)；u ?e ? io | db； ?ai ?3?3 e ? U );o ?e ?3?3 e ? U；u ?e ? io | db； ?n=0:1:M-1;wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);w_bl

28、a=( boxcar(M);h=hd.*w_bla;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(Qi ?3?3 e ?u)；?D? )；?e ?3?3 e ? v)；subplot(2,2,4);plot(W/pi,20*log10(abs(H);title(clear allfigureAs=55;ws=11*pi/27;wp=15*pi/27;tr_width=wp-ws;%?1y?e?M=ceil(As-7.

29、95)*2*pi/(14.36*tr_width)+1)+1;%disp( ?2 - ?+卩?3O ?num2str(M);beta=0.1102*(As-8.7);%?-? o ?beta?卩n=0:1:M-1;disp( ?D? a ?2 - ?- ,num2str(beta); w_kai=(kaiser(M,beta);%?o ?-? o o _ey-u ?e ? 10 | db; ?-? 0?2 - ?+卩?3O ? ewc=(ws+wp)/2;hd=ideal_lp(pi,M)-ideal_lp(wc,M);%?0ai ?3? 1o|h=hd.*w_kai;db,mag,pha,g

30、rd,w=freqz_m(h,1);delta_w=2*pi/1000;Rp=-(min(db(wp/delta_w+1:1:501);disp( ?e？2 ? :nU,iari2str(Rp);As=-round(max(db(1:1:ws/delta_w+1);disp( x ?D?xe ? ?a?,num2str(As);subplot(2,2,1);stei(n,hd);title( ai ?3? 10 I); axis(0 M-1 -0.4 0.8);ylabel( hd(n);subplot(2,2,2);stei(n,w_kai);title( ?-? 6;axis(0 M-1

31、0 1.1);ylabel( wd(n);subplot(2,2,3);stem(n,h);title( ?e ?3? 16 |;axis(0 M-1 -0.4 0.8);xlabel( n );ylabel( h(n);subplot(2,2,4);plot(w/pi,db);title( ? u ?e ? 16 |/dB);axis(0 1 -100 10);grid;xlabel( o ?pi?a 卩丫 ? e；ylabel( ? ?e y/dB);clear allfigure%?u matlab?DD? - o _e yg ?u ? e ?0Wp=0.6*pi;Ws=0.4*pi;t

32、r_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.*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)hd(n)w(n) )stem(n,hd);title(Ideal Impulse Rrespons

33、eaxis(0 M-1 -0.4 0.5);xlabel(n );ylabel(subplot(222)stem(n,w_ham);title(Hamming Window )axis(0 M-1 0 1.1);xlabel(n );ylabel(subplot(223)decibels );decibels );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

34、in db );grid axis(0 1 -100 10);xlabel( frequence in pi unit );ylabel( clear allfigureWp=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_w=2*pi/1000;Ap=-(min(db(Wp

35、/delta_w+1:1:501) As=-round(max(db(1:1:Ws/delta_w+1)figure 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_ham);title( Hamming Window ) axis(0 M-1 0 1.1);xlabel( n );ylabel( w(n) ) subplot(223)stem(n,h);title( Actual Impu

36、se 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 );grid axis(0 1 -100 10);xlabel( frequence in pi unit );ylabel(%o ?Hamming o y ee ?FIR e y X ?2 ？*figure wp=0.2*pi;ws=0.3*piN=61 n=0:1:N-1wc=(ws+wp)/2;?？2 ？*hd=idealp(wc,N);?卩丨丨卩?3

37、? 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_w+1)%e ?e ?厂? ?As=-round(max(db(ws/delta_w+1:1:501)%e ?e ?X ?D?Xe ? ?subplot(221);stem(n,hd);title(I ?3? O ? io )axis(0 N-1 -0.1 0.3);xlabel(n );ylabel(hd(n)subplot(222);s

38、tem(n,w_ham);title(hamming )axis(0 N-1 0 1.1);xlabel(n);ylabel(w(n);subplot(223);stem(n,h);title(?e 3? O ? io) |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(o ?PI?a 卩 ? ii ? ii ? e;ylabel(? e y uclear allfigure%o ?Kaiser o_e y ee ?FIR e yX ?2 - ?-wp

39、=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?卩丨厂卩?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)%?ei ?ei ? ?As=-round(max(db(ws/delta_w+1:1

40、:501)%ei ?ei ?X ?D?Xe ? ?subplot(211);plot(w/pi,db);title(?- e - u ?e ? io |(dB0ridaxis(0 0.5 -100 0)ylabel( ? e y u ?e /db);xlabel( o ?/pi?a i ?i ?i ?e)subplot(212);plot(w/pi,pha);title(? a ?io |);gridaxis(0 0.5 -4 4)ylabel( ? a ? );xlabel( o ?/pi?ai ?i ?i ?e) 13文档来源为:从网络收集整理.word版本可编辑.? e /db);n )

41、;);n ););figure h1=-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 |h1); ylabel( h(n) );xlabel( subplot(2,1,2),plot(w/pi,A); ylabel( A );xlabel( pi title( ? u ?卩? io)j figurerz=roots(h1)for i=1:8 r(i)=1/rz(i);endrzplane(h1,1

42、);title(hi a ?卩? ?);clearallfigureh2=-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( subplot(2,1,2),plot(w/pi,A); ylabel( A );xlabel( pi title( u ? io);|figurerz=roots(h2)for i=1:8 r(i)=1/rz(i)

43、;endrzplane(h2,1);title(?h2 a ?卩? ?);clearallfigure 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( subplot(2,1,2),plot(w/pi,A); ylabel( A );xlabel( pi title( u ?卩? io I? figurerz=roots(h3)for

44、i=1:8r(i)=1/rz(i);endrzplane(h3,1);title(h3 a ?卩?);clearallfigure h4=-4,1,-1,-2,5,6,-6,-5,2,1,-1,4;M=length(h4);n=0:M-1; A,w,type,tao=amplres(h4);type subplot(2,1,1),stem(n,h4);title( 3? O ?io Ih4 ); ylabel( h(n) );xlabel( subplot(2,1,2),plot(w/pi,A); ylabel( A );xlabel( pi title( ? u ? io |?figurer

45、z=roots(h4)for i=1:8 r(i)=1/rz(i);endrzplane(h4,1);title( h4 a ?卩?);n ););n ););h = 1,2,3,2,1;M = 5;n = (M-3)/2;a = 2*h(3-n)k = -500:500;K = 500;w = k*pi/K;p = a*cos(w*n);plot(w,p);ylabel( p(w) ); xlabel( w );function db,mag,pha,H,w=freqz_m3(b,a) H,w=freqz(b,a,1000, whole );H=(H(1:1:501); w=(w(1:1:5

46、01);mag=abs(H); db=20*log10(mag+eps)/max(mag);pha=angle(H);endfunction db,mag,pha,grd,w = freqz_m(b,a);H,w = freqz(b,a,1000,whole );H = (H(1:1:501); w = (w(1:1:501);mag = abs(H);db = 20*log10(mag+eps)/max(mag);pha = angle(H);grd = grpdelay(b,a,w);function db, mag, pha, w=freqz_m2(b,a);%Modified vers

47、ion of freqz subroutine H,w=freqz(b,a,1000, whole );H=(H(1:501); w=(w(1:501);mag=abs(H); db=20*log10(mag+eps)/max(mag);pha=angle(H);function Hr,w,b,L = Hr_Type2(h);M = length(h); L = M/2;b = 2*h(L:-1:1); n = 1:1:L; n = n-0.5;w = 0:1:500*pi/500; Hr = cos(w*n)*b;endfunctionHr,w,c,L = Hr_Type3(h);% Com

48、putes Amplitude response Hr(w) of a Type-3 LP FIR filter % % Hr,w,c,L = Hr_Type3(h)% Hr = Amplitude Response% w = frequencies between 0 pi over which Hr is computed% c = Type-3 LP filter coefficients% L = Order of Hr% h = Type-3 LP impulse response%M = length(h); L = (M-1)/2;c = 2*h(L+1:-1:1); n = 0:1:L;w = 0:1:500*pi/500; Hr = sin(w*n)*c; end