1、實驗卡和實驗報告 信息科學與工程學院 課程編號 實驗項目序號本科學生實驗卡和實驗報告信息科學與工程學院通信工程專業 2013級 1301 班課程名稱：數字信號處理實驗項目：FIR濾波器設計20152016學年 第 二 學期學號：201308030104_ 姓名：_ _王少丹_ _ 專業年級班級： _通信1301_四合院_ 實驗室 組別_ 實驗日期 _2016 年_ 6_ 月_12_ 日課程名稱數字信號處理實驗課時4實驗項目名稱和編號FIR濾波器設計同組者 姓 名實驗目的學會運用MATLAB產生常用離散時間信號實驗環境MATLAB實驗內容和原理（1）認真復習 FIR 數字濾波器的基本概念，線性相

2、位 FIR 濾波器的條件和特點、幅度函數特點、零點位置的基本特點與性質；窗函數設計法的基本概念與方法，各種窗函數的性能和設計步驟，線性相位 FIR 低通、高通、帶通和帶阻濾波器的設計方法，頻率采樣設計法的基本概念和線性相位的實現方法。 （2）掌握幾種線性相位的特點，熟悉和掌握矩形窗、三角形窗、漢寧窗、海明窗、布萊克曼窗、 凱塞窗設計 IIR 數字濾波器的方法， 熟悉和掌握頻率抽樣設計法的線性相位的設計方法，并對各種線性相位的頻率抽樣法的設計給出調整和改進。（3）熟悉利用 MATLAB 進行各類 FIR 數字濾波器的設計方法。實驗步驟方 法關鍵代碼a. 設線性相位 FIR 濾波器單位抽樣響應分別

3、為 h(n )= -4,1, -1,- 2,5,6,5, -2, -1,1, -4h(n)= -4,1, -1,- 2,5,6,6,5,- 2, -1,1,- 4h(n)= -4,1,- 1, -2,5,0,- 5,2,1,- 1,4h(n)= -4,1, -1, -2,5,6,- 6, -5,2,1, -1,4分別求出濾波器的幅度頻率響應 H()，系統函數 H(z)以及零極點分布，并繪制相應的波形和分布圖。 Type2： Type3： Type4： b. 設計 FIR 數字低通濾波器，技術指標為：p=0.2，st=0.3，1=0.25dB，2=50dB。 (1) 通過技術指標，選擇一種窗函數

4、進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。 (3) 選擇凱塞窗函數設計該濾波器，并繪制相應的波形圖。c. 設計 FIR 數字帶通濾波器，技術指標為： 下阻帶邊緣：st1=0.2，s1=60dB，下通帶邊緣：p1=0.35，p1=1dB； 上通帶邊緣：p2=0.65，p1=1dB，上阻帶邊緣：st2=0.8，s2=60dB； (1) 通過技術指標，選擇一種窗函數進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。d. 設計 FIR 數字帶通濾波器，技術指標為： 下阻帶邊緣：st1=0.2，s1=60dB，下通帶邊緣：p1=0.4，p1=1dB； 上通帶邊緣：

5、p2=0.6，p1=1dB，上阻帶邊緣：st2=0.8，s2=60dB； (1) 通過技術指標，選擇一種窗函數進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。e. 設計 FIR 數字帶通濾波器，技術指標為： 下阻帶邊緣：st1=0.2，s1=20dB，下通帶邊緣：p1=0.4，p1=1dB； 上通帶邊緣：p2=0.6，p1=1dB，上阻帶邊緣：st2=0.8，s2=20dB； (1) 通過技術指標，選擇一種窗函數進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。f. 設計 FIR 數字高通濾波器，技術指標為：通帶截止頻率為p=15/27，阻帶截止頻率為st=1

6、1/27，通帶最大衰減為1=2.5dB，阻帶最小衰減為2=55dB。 (1) 通過技術指標，選擇一種窗函數進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。系統一：g. 設計 FIR 數字高通濾波器，技術指標為：通帶截止頻率為p=0.6，阻帶截止頻率為st=0.4，通帶最大衰減為1=0.25dB，阻帶最小衰減為2=40dB。 (1) 通過技術指標，選擇一種窗函數進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。 h. 濾波器的技術指標為：通帶截止頻率為p=0.6，阻帶截止頻率為st=0.4，通帶最大衰減為1=0.25dB，阻帶最小衰減為2=40dB。 (1) 通

7、過技術指標，選擇一種窗函數設計一個具有/2 相移的 FIR 高通濾波器； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。 i. 設計 FIR 數字帶阻濾波器，其技術指標為： 低端阻帶邊緣：st1=0.4，s1=40dB，低端通帶邊緣：p1=0.2，p1=1dB； 高端通帶邊緣：p2=0.8，p1=1dB，高端阻帶邊緣：st2=0.6，s2=40dB； (1) 通過技術指標，選擇一種窗函數進行設計； (2) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。 FIR 濾波器的單位抽樣響應為h(n)=1/9，編制 MATLAB 程序求系統的頻率采樣型結構的系數，并畫出頻率抽樣型結構。m. 一個

8、理想差分器的頻率響應為： 用長度為 21 的漢寧窗設計一個數字 FIR 差分器，并繪制其時域和頻率的響應波形。 n. 利用漢寧窗設計一個長度為 25 的數字希爾伯特變換器，并繪制它的時域和頻域的響應波形。 p. FIR 數字低通濾波器的技術指標為：p=0.2，st=0.3，1=0.25dB，2=50dB。利用頻率采樣方法設計 FIR 數字濾波器， 并繪制濾波器的單位沖激響應、 幅度頻率響應的波形。 q. 用窗函數法設計一個線性相位的 FIR 數字低通濾波器，其技術指標為：p=0.2，st=0.4，1=0.25dB，2=50dB。 (1) 求濾波器的單位抽樣響應、頻率響應，并繪制波形。 (2)

9、對該濾波器輸入一個寬度為 10 的矩形序列，求濾波器的輸出信號，并繪制相應的波形圖。測試記錄分 析結 論  理論與實驗值相符。（后附代碼）小 結認真復習 FIR 數字濾波器的基本概念，線性相位 FIR 濾波器的條件和特點、幅度函數特點、零點位置的基本特點與性質；窗函數設計法的基本概念與方法，各種窗函數的性能和設計步驟，線性相位 FIR 低通、高通、帶通和帶阻濾波器的設計方法，頻率采樣設計法的基本概念和線性相位的實現方法。 以下由實驗教師填寫記 事評 議  成績評定  平時成績_ 實驗報告成績_ 綜合成績 _ 指導教師簽名：代碼：wp = 0.2*pi;ws = 0

10、.4*pi;tr_width = ws-wp;M = ceil(6.6*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('

11、;Ideal Impulse Rresponse')axis(0 M-1 -0.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(

12、'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')%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)

13、xw = Rn*(exp(-j).(n'*w)answ = H.*xwRnn = answ*exp(-j).(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 sig

14、nal in F-domain')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')

15、M = 25; alpha = (M-1)/2; n = 0:M-1;hd = (2/pi)*(sin(pi/2)*(n-alpha).2)./(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

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

17、plr(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; wdl = 0,0

18、.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,Hdr);axis

19、(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,wl(1:11)/pi,Hrs(1:11),&

20、#39;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('Decibels'); Wpl=0.2*pi;W

21、ph=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(Wsl/delta_w+1:1:Wsh/de

22、lta_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(n,h);title('Actual Impuse R

23、esponse')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 allfigureWpl=0.2*pi;Wph=0.8*pi;Wsl=0.4*pi;Wsh=0.6*pi;tr_width

24、=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);title(

25、9;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

26、('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');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

27、/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;H,W=freqz(h,1);subplot(2,2,1);stem(n,hd);title('ÀíÏëÂö³å³éÑù');subplot(2,2,2);stem(n,w_bla);title('²

28、¼À³¿ËÂü´°');subplot(2,2,3);stem(n,h);title('ʵ¼ÊÂö³å³éÑù');subplot(2,2,4);plot(W/pi,20*log10(abs(H);title('·ù¶ÈÏìÓ¦£¨db£

29、9;'); 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=(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('

30、ÀíÏëÂö³å³éÑù');subplot(2,2,2);stem(n,w_bla);title('º£Ã÷´°');subplot(2,2,3);stem(n,h);title('ʵ¼ÊÂö³å³éÑù');subplot(2,2,4);plot(W/p

31、i,20*log10(abs(H);title('·ù¶ÈÏìÓ¦£¨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)-ide

32、al_lp(wc1,M);w_bla=( 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('

33、02;µ¼ÊÂö³å³éÑù');subplot(2,2,4);plot(W/pi,20*log10(abs(H);title('·ù¶ÈÏìÓ¦£¨db£©'); clear allfigureAs=55;ws=11*pi/27;wp=15*pi/27; tr_width=wp-ws;%¼ÆËã

34、5;ý¶É´ø M=ceil(As-7.95)*2*pi/(14.36*tr_width)+1)+1; %°´¿­Ôó´°¼ÆËãÂ˲¨Æ÷µÄ³¤¶Èdisp('Â˲¨Æ÷µÄ³¤¶&#

35、200;',num2str(M);beta=0.1102*(As-8.7); %¼ÆËã¿­Ôó´°µÄbetaÖµn=0:1:M-1;disp('ÏßÐÔÏàλÂ˲¨Æ÷',num2str(beta);w_kai=(kaiser(M,beta)'%Çó¿

36、;­Ôó´°º¯Êýwc=(ws+wp)/2;hd=ideal_lp(pi,M)-ideal_lp(wc,M); %ÇóÀíÏëÂö³åÏìÓ¦h=hd.*w_kai; db,mag,pha,grd,w=freqz_m(h,1);delta_w=2*pi/1000;Rp=-(min(db(wp/delta_w+1:1:501);disp('Ê

37、1;¼Êͨ´ø²¨¶¯Îª',num2str(Rp);As=-round(max(db(1:1:ws/delta_w+1);disp('×îС×è´øË¥¼õΪ?',num2str(As);subplot(2,2,1);stem(n,hd);title('ÀíÏë

38、;Âö³åÏìÓ¦');axis(0 M-1 -0.4 0.8);ylabel('hd(n)');subplot(2,2,2);stem(n,w_kai);title('¿­Ôó´°'); axis(0 M-1 0 1.1);ylabel('wd(n)');subplot(2,2,3);stem(n,h);title('ʵ¼ÊÂö&

39、#179;åÏìÓ¦');axis(0 M-1 -0.4 0.8);xlabel('n');ylabel('h(n)');subplot(2,2,4);plot(w/pi,db);title('?·ù¶ÈÏìÓ¦/dB');axis(0 1 -100 10);grid;xlabel('ÒÔpiΪµ¥Î»µ

40、6;ƵÂÊ');ylabel('?·Ö±´Êý/dB'); clear allfigure%ÔÚmatlabÖÐн¨º¯Êýg£¬´úÂëÈçÏ£ºWp=0.6*pi;Ws=0.4*pi;tr_width=Wp-Ws;M=cei

41、l(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)stem(n,hd);title('Ideal Impulse Rresponse')axis(0 M-1 -0.4 0.

42、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 Impuse Response')axis(0 M-1 -0.4 0.5);xlabel('n');ylabel('h(n)')subplot(224)p

43、lot(w/pi,db);title('Magnitude Response in db');gridaxis(0 1 -100 10);xlabel('frequence in pi unit');ylabel('decibels'); 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.*rot

44、90(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)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 Wi

45、ndow')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');gridaxis(0 1 -100 10);xlabel('freq

46、uence in pi unit');ylabel('decibels');%ÓÃHamming´°º¯ÊýÉè¼ÆFIRÊý×ÖÂ˲¨Æ÷figurewp=0.2*pi;ws=0.3*piN=61n=0:1:N-1wc=(ws+wp)/2;%ÀíÏëµÍͨÂ&

47、#203;²¨Æ÷hd=ideal_lp(wc,N);%ÀíÏëµÍͨµÄ³å¼¤ÏìÓ¦w_ham=(hamming(N)'h=hd.*w_ham;%FIRÂ˲¨Æ÷³å¼¤ÏìÓ¦db,mag,pha,grd,w=freqz

48、_m(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)%ʵ¼ÊµÄ×îС×è´øË¥¼õsubplot(221);st

49、em(n,hd);title('ÀíÏë³å¼¤ÏìÓ¦')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');ylabel('w(n)');subplot(223);s

50、tem(n,h);title('ʵ¼Ê³å¼¤ÏìÓ¦')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('ÒÔPIΪµ¥Î»µÄƵÂ&

51、#202;');ylabel('¶ÔÊý·ù¶È/db'); clear allfigure%ÓÃKaiser´°º¯ÊýÉè¼ÆFIRÊý×ÖÂ˲¨Æ÷wp=0.2*pi;ws=0.3*pi;As=50tr_width=ws-wpN=ceil(As-7.95)/(14.3

52、6*tr_width/(2*pi)+1)+1n=0:1:N-1beta=0.1102*(As-8.7)wc=(wp+ws)/2%ÀíÏëµÍͨµÄ½ØֹƵÂÊhd=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

53、)%?ʵ¼ÊµÄͨ´øË¥¼õAs=-round(max(db(ws/delta_w+1:1:501)%ʵ¼ÊµÄ×îС×è´øË¥¼õsubplot(211);plot(w/pi,db);title('¿­É­´

54、°·ù¶ÈÏìÓ¦(dB)');gridaxis(0 0.5 -100 0)ylabel('¶ÔÊý·ù¶È/db');xlabel('ÒÔ/piΪµ¥Î»µÄƵÂÊ')subplot(212);plot(w/pi,pha);title('

55、;ÏàλÏìÓ¦');gridaxis(0 0.5 -4 4)ylabel('Ïàλ');xlabel('ÒÔ/piΪµ¥Î»µÄƵÂÊ')figureh1=-4,1,-1,-2,5,6,5,-2,-1,1,-4;M=length(h1);n=0:M-1;A,w,type,tao=amplr

56、es(h1);typesubplot(2,1,1),stem(n,h1);title('³å¼¤ÏìÓ¦h1');ylabel('h(n)');xlabel('n');subplot(2,1,2),plot(w/pi,A);ylabel('A');xlabel('pi');title('?·ùƵÏìÓ¦');figurerz=roots(h

57、1)for i=1:8r(i)=1/rz(i);endr'zplane(h1,1);title('h1Á㼫µãͼ'); 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);typesubplot(2,1,1),stem(n,h2);title('³å¼¤ÏìÓ¦h2

58、');ylabel('h(n)');xlabel('n');subplot(2,1,2),plot(w/pi,A);ylabel('A');xlabel('pi');title('·ùƵÏìÓ¦');figurerz=roots(h2)for i=1:8r(i)=1/rz(i);endr'zplane(h2,1);title('?h2Á㼫µã

59、05;¼'); clear allfigureh3=-4,1,-1,-2,5,0,-5,2,1,-1,4;M=length(h3);n=0:M-1;A,w,type,tao=amplres(h3);typesubplot(2,1,1),stem(n,h3);title('³å¼¤ÏìÓ¦h3');ylabel('h(n)');xlabel('n');subplot(2,1,2),plot(w/pi,A);ylabel('A');xla