實驗6FIR濾波器設(shè)計

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

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

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分別求出濾波器的幅度頻率響應(yīng) H()，系統(tǒng)函數(shù) H(z)以及零極點分布，并繪制相應(yīng)的波形和分布圖。 Type2： Type3： Type4： b. 設(shè)計 FIR 數(shù)字低通濾波器，技術(shù)指標(biāo)為：p=0.2，st=0.3，1=0.25dB，2=50dB。 (1) 通過技術(shù)指標(biāo)，選擇一種窗函數(shù)

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

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

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

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

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

9、對該濾波器輸入一個寬度為 10 的矩形序列，求濾波器的輸出信號，并繪制相應(yīng)的波形圖。測試記錄分 析結(jié) 論  理論與實驗值相符。（后附代碼）小 結(jié)認(rèn)真復(fù)習(xí) FIR 數(shù)字濾波器的基本概念，線性相位 FIR 濾波器的條件和特點、幅度函數(shù)特點、零點位置的基本特點與性質(zhì)；窗函數(shù)設(shè)計法的基本概念與方法，各種窗函數(shù)的性能和設(shè)計步驟，線性相位 FIR 低通、高通、帶通和帶阻濾波器的設(shè)計方法，頻率采樣設(shè)計法的基本概念和線性相位的實現(xiàn)方法。 以下由實驗教師填寫記 事評 議  成績評定  平時成績_ 實驗報告成績_ 綜合成績 _ 指導(dǎo)教師簽名：代碼：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