1、<p>　　數(shù) 字 信 號(hào) 處 理</p><p><b>　　課 程 設(shè) 計(jì)</b></p><p>　　名 稱 基于窗函數(shù)設(shè)計(jì)法線性相位型</p><p>　　FIR數(shù)字低通濾波器設(shè)計(jì) </p><p>　　姓 名

2、 </p><p>　　班 級(jí) </p><p>　　學(xué) 號(hào) </p><p>　　指 導(dǎo) 教 師

3、 </p><p><b>　　目錄 </b></p><p>　　一 內(nèi)容概述 ……………………………………… 3</p><p>　　二 題目要求細(xì)則 ……………………………………… 4</p><p>　　三 設(shè)計(jì)原理 ……………………………………… 5</p>

4、;<p>　　四 設(shè)計(jì)過程 ……………………………………… 6</p><p>　　五 實(shí)驗(yàn)結(jié)果及分析 ……………………………………… 7</p><p>　　六 實(shí)驗(yàn)心得 ……………………………………… 14</p><p>　　七 參考文獻(xiàn) ……………………………………… 15</p><p&

5、gt;　　八 源代碼 ……………………………………… 16</p><p>　　《數(shù)字信號(hào)處理課程設(shè)計(jì)》任務(wù)書</p><p>　　一【內(nèi)容概述】 </p><p>　　FIR數(shù)字濾波器設(shè)計(jì)最簡(jiǎn)單的方法是窗函數(shù)法，通常也稱之為傅里葉級(jí)數(shù)法。它是在屬于進(jìn)行的，因而必須由理想濾波器的頻率響應(yīng)推導(dǎo)出其單位沖擊響應(yīng)，再設(shè)計(jì)一個(gè)FIR數(shù)字濾波器的單位沖擊

6、響應(yīng)去逼近。</p><p>　　窗函數(shù)設(shè)計(jì)法的基本原理是用一定寬度窗函數(shù)截取無限脈沖響應(yīng)序列獲得有限長(zhǎng)的脈沖響應(yīng)序列,先通過傅里葉逆變換獲得理想濾波器的單位脈沖響應(yīng)hd(n)。再由性能指標(biāo)確定窗函數(shù)W(n)和窗口長(zhǎng)度N。最后求得實(shí)際濾波器的單位脈沖響應(yīng)h(n)。</p><p>　　窗函數(shù)在設(shè)計(jì)FIR數(shù)字濾波器中有很重要的作用，正確的選擇窗函數(shù)可以提高所涉及的數(shù)字濾波器的性能，或者在滿足

7、設(shè)計(jì)要求的情況下，減小FIR數(shù)字濾波器的階數(shù)。</p><p><b>　　二【題目要求細(xì)則】</b></p><p>　　基于窗函數(shù)設(shè)計(jì)法線性相位型FIR數(shù)字低通濾波器設(shè)計(jì)</p><p>　　目的：① 熟練掌握MATLAB工具軟件在工程設(shè)計(jì)中的使用；</p><p>　?、?熟練掌握窗函數(shù)法設(shè)計(jì)線性相位型FIR濾波器

8、的方法。</p><p>　　要求：① 根據(jù)給定ALPF幅頻特性指標(biāo)，生成與之相對(duì)應(yīng)的模擬濾波器h(t)； </p><p>　?、?根據(jù)h(t)取得與之想對(duì)應(yīng)的數(shù)字濾波器沖擊響應(yīng)函數(shù)h(n)；</p><p>　?、?選擇適當(dāng)?shù)拇昂瘮?shù)，構(gòu)建線性相位型FIR的沖擊響應(yīng)函數(shù)hd(n)；</p><p>　?、?設(shè)計(jì)與之相對(duì)應(yīng)的DLPF，對(duì)比分析

9、DLPF幅頻特性是否符合要求；</p><p>　?、?產(chǎn)生一有干擾頻率的時(shí)域序列（借助FFT分析說明其有干擾），使之通過所設(shè)計(jì)的DLPF，</p><p>　　對(duì)濾波輸出結(jié)果作出分析，說明輸出結(jié)果。</p><p>　　步驟：① 根據(jù)給定設(shè)計(jì)指標(biāo)（通帶截止頻率fb＝1 000 Hz，通帶最大衰減αp=1dB，阻帶截止頻率fc</p><p>

10、;　?。? 200 Hz，阻帶最小衰減αs=100 dB，令T=2）生成與之相對(duì)應(yīng)的模擬濾波器h(t)；</p><p>　?、?根據(jù)h(t)取得與之想對(duì)應(yīng)的數(shù)字濾波器沖擊響應(yīng)函數(shù)h(n)；</p><p>　?、?選擇適當(dāng)?shù)拇昂瘮?shù)，構(gòu)建線性相位型FIR的沖擊響應(yīng)函數(shù)hd(n)；</p><p>　　④ 設(shè)計(jì)與之相對(duì)應(yīng)的DLPF，給出窗函數(shù)及所設(shè)計(jì)濾波器的

11、幅度特性，并對(duì)比理想沖激響應(yīng)</p><p>　　和實(shí)際沖擊響應(yīng)，對(duì)比分析DLPF幅頻特性是否符合要求；</p><p>　?、?試說明不同窗函數(shù)對(duì)所設(shè)計(jì)數(shù)字濾波器性能的影響；</p><p>　?、?產(chǎn)生一有干擾頻率的時(shí)域序列（借助FFT分析說明其有干擾），使之通過所設(shè)計(jì)的DLPF，</p><p>　　對(duì)濾波輸出結(jié)果作出分析，說明輸出結(jié)果

12、。</p><p><b>　　三【設(shè)計(jì)原理】</b></p><p>　　窗函數(shù)設(shè)計(jì)法的基本原理是用一定寬度窗函數(shù)截取無限脈沖響應(yīng)序列獲得有限長(zhǎng)的脈沖響應(yīng)序列,先通過傅里葉逆變換獲得理想濾波器的單位脈沖響應(yīng)hd(n)。再由性能指標(biāo)確定窗函數(shù)W(n)和窗口長(zhǎng)度N。最后求得實(shí)際濾波器的單位脈沖響應(yīng)h(n)。 </p><p>　　FIR數(shù)字濾波器

13、設(shè)計(jì)最簡(jiǎn)單的方法是窗函數(shù)法，通常也稱之為傅里葉級(jí)數(shù)法。它是在屬于進(jìn)行的，因而必須由理想濾波器的頻率響應(yīng)推導(dǎo)出其單位沖擊響應(yīng)，再設(shè)計(jì)一個(gè)FIR數(shù)字濾波器的單位沖擊響應(yīng)去逼近。</p><p>　　假設(shè)理想濾波器的頻率響應(yīng)的表示為：</p><p>　　由于濾波器的頻率響應(yīng)與單位沖擊響應(yīng)序列是傅里葉變換對(duì)，據(jù)此：</p><p>　　求得理想濾波器的單位沖擊響應(yīng)后，可得

14、到該離散濾波器的系統(tǒng)傳遞函數(shù)：</p><p>　　注意，該為無限長(zhǎng)序列，因此是物理不可實(shí)現(xiàn)的。</p><p>　　為了使系統(tǒng)變?yōu)槲锢砜蓪?shí)現(xiàn)的，且使實(shí)際的FIR濾波器頻率響應(yīng)盡可能逼急理想濾波器的頻率響應(yīng)，采用窗函數(shù)將無線脈沖響應(yīng)截取一段好來近似表示，可得：，從而有。式中N表示窗口長(zhǎng)度，這樣就是物理可實(shí)現(xiàn)的系統(tǒng)。</p><p><b>　　四【設(shè)計(jì)過程

15、】</b></p><p>　　先產(chǎn)生一個(gè)理想低通，分別用三種不同的窗函數(shù)與之時(shí)域內(nèi)點(diǎn)乘后，得到可物理實(shí)現(xiàn)的低通濾波器h(n)，觀察它的幅度相頻特性。</p><p>　　產(chǎn)生一個(gè)干擾序列，加上噪聲后得到一個(gè)有毛刺的正弦波，通過三種加不同窗的濾波器，得出三個(gè)慮完波的時(shí)域序列和頻譜。比較三個(gè)慮后的時(shí)域頻域波形，比較出三種不同濾波器的優(yōu)劣。</p><p>

16、　　五【實(shí)驗(yàn)結(jié)果及分析】</p><p>　　1.設(shè)計(jì)參數(shù)：通帶截止頻率fb＝1 000 Hz，通帶最大衰減αp=1dB，阻帶截止頻率fc</p><p>　?。? 200 Hz，阻帶最小衰減αs=100 dB，令T=2</p><p>　　2.不同窗形狀對(duì)比：</p><p>　　Blackman窗、漢寧窗、海明窗、凱塞窗： &l

17、t;/p><p>　　分析：若要根據(jù)指標(biāo)設(shè)計(jì)出理想的數(shù)字濾波器，只有布萊克曼窗和凱塞窗可以接近設(shè)計(jì)要求由于除凱塞窗之外的窗函數(shù)全部是參數(shù)固定窗函數(shù)，所以若要想達(dá)到給定設(shè)計(jì)指標(biāo)阻帶最小衰減αs=100dB，只有選擇凱塞窗，估算α值，才有可能達(dá)到要求。由圖可知，凱塞窗和布萊克曼窗比較接近實(shí)驗(yàn)要求，漢寧窗和海明窗較差，所以matlab默認(rèn)窗函數(shù)海明窗在此設(shè)計(jì)指標(biāo)不太合適。</p><p>　　3.理

18、想單位脈沖響應(yīng)對(duì)比：</p><p>　　布萊克曼： 漢寧：</p><p>　　海明： 凱塞：</p><p>　　4.加窗后截取實(shí)際單位脈沖響應(yīng)對(duì)比：</p><p>　　布萊克曼： 漢寧：</p><p&g

19、t;　　海明： 凱塞：</p><p>　　分析：用不同的窗函數(shù)去截取序列會(huì)產(chǎn)生不同的誤差，也就是通常說的吉布斯效應(yīng)，會(huì)引起過渡帶加寬以及通帶和阻帶內(nèi)的波動(dòng)，尤其使阻帶衰減減小。由圖可知，布萊克曼窗和凱塞窗圖中最大正、負(fù)峰值對(duì)應(yīng)頻率間距較小，而海明窗和漢寧窗較大，所以選擇后兩者會(huì)增大截?cái)嘈?yīng)</p><p>　　5.不同窗函數(shù)幅度響應(yīng)對(duì)

20、比：</p><p><b>　　布萊克曼：</b></p><p><b>　　漢寧：</b></p><p><b>　　海寧：</b></p><p><b>　　凱塞窗：</b></p><p>　　分析：要想獲得理想的阻帶

21、衰減，需要使用窗譜最大旁瓣相對(duì)幅度較小的窗函數(shù)。由圖可知，漢寧窗的阻帶波紋最大，所以其阻帶衰減就最小。而布萊克曼窗和凱塞窗阻帶波紋相對(duì)較小，所以其阻帶衰減較明顯。而凱塞窗可以通過調(diào)整窗長(zhǎng)度以及α值可以有效控制過渡帶寬度，獲得良好阻帶衰減。所以選擇不同的窗函數(shù)可以減少帶內(nèi)波動(dòng)以增大阻帶衰減。</p><p>　　6.濾波器性能分析：</p><p>　?。╝p:通帶最大衰減；as:阻帶最小衰

22、減）</p><p>　　Blackman窗：Ap1 = 0.0027</p><p><b>　　As1 =73</b></p><p>　　漢寧窗： Ap2 =0.0711</p><p><b>　　As2 =44</b></p><p>　　海明窗： Ap

23、3 =0.0282</p><p>　　As3 =51 </p><p>　　凱塞窗： Ap4 =0.0056</p><p><b>　　As4 =108</b></p><p>　　例如：布萊克曼窗指數(shù)：</p><p>　　通過分別對(duì)系統(tǒng)加Blackman窗、漢寧窗、海明窗、凱

24、塞窗得到以上結(jié)果。從加窗后系統(tǒng)的幅頻響應(yīng)曲線可以看出，Blackman窗凱塞窗引起的過渡帶加寬現(xiàn)象和阻帶波動(dòng)較小，即截?cái)嘈?yīng)較小，而且通過比較可以看出，凱塞窗還具有最小的通帶最大衰減和阻帶最小衰減，綜合比較發(fā)現(xiàn)凱塞窗blackman窗的各項(xiàng)性能優(yōu)于海明窗和漢寧窗，其中凱塞窗性能最好。所以，選擇凱塞窗選擇不同的α值，可以設(shè)計(jì)出性能最優(yōu)的FIR濾波器。</p><p>　　7．產(chǎn)生一個(gè)有干擾頻率的時(shí)域序列，使之通過所

25、設(shè)計(jì)的DLPF，</p><p>　　對(duì)濾波輸出結(jié)果作出分析，說明輸出結(jié)果。</p><p><b>　　布萊克曼：</b></p><p><b>　　漢寧窗：</b></p><p><b>　　海明窗：</b></p><p><b>　

26、　凱塞窗：</b></p><p>　　分析：由圖明顯看出時(shí)域和頻域部分分別經(jīng)過低通濾波器后高頻部分被濾掉，時(shí)域部分毛刺被濾掉。</p><p><b>　　8．?dāng)U展部分：</b></p><p><b>　　第一種：</b></p><p><b>　　第二種：</b&

27、gt;</p><p>　　分析：由圖可知，通過帶通濾波器后低頻和高頻被濾掉</p><p><b>　　七．【參考文獻(xiàn)】</b></p><p><b>　　參考文獻(xiàn) </b></p><p>　　高西全，數(shù)字信號(hào)處理（第三版）西安電子科技大學(xué)出版社</p><p>　　程

28、佩青，數(shù)字信號(hào)處理教程（第三版），清華大學(xué)出版社</p><p>　　MATLAB7.0在數(shù)字信號(hào)處理中的應(yīng)用  羅軍輝  機(jī)械工業(yè)出版社</p><p>　　飛思科技產(chǎn)品研發(fā)中心編著.MATLAB7輔助信號(hào)處理技術(shù)與應(yīng)用[M].北京。 電子工業(yè)出版社，2005. </p><p>　　5、數(shù)字信號(hào)處理原理及其MATLA

29、B實(shí)現(xiàn) 從玉良  電子工業(yè)出版社 </p><p>　　八．【實(shí)驗(yàn)部分代碼】</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　% </p><p>

30、　　% 基于窗函數(shù)設(shè)計(jì)法線性相位型FIR數(shù)字低通濾波器設(shè)計(jì)</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p><b>　　%%</b></p><p>　　%%%%%%%%%%%%%%%%%%%%%% 布萊克曼 %%%%%%%%%%%%%%

31、%%%%%%%%%%</p><p><b>　　%%</b></p><p><b>　　clc,clear</b></p><p>　　fp=1000; </p><p><b>　　fc=1200;</b></p><p>　　fs=4000

32、; %對(duì)信號(hào)的采樣頻率為4000</p><p>　　ws=2*pi*fc/fs %t通帶截至頻率</p><p>　　wp=2*pi*fp/fs %阻帶截至頻率</p><p>　　tr_width=ws-wp; %求過渡帶</p><

33、p>　　M=ceil(11*pi/tr_width)+1 %利用blackman窗的過渡帶=(11pi/M),求窗寬M</p><p>　　n=0:1:M-1;</p><p>　　wc=(ws+wp)/2; %求通帶截止頻率wc</p><p>　　alpha=(M-1)/2;</p><p>　　n=[

34、0:(M-1)];</p><p>　　m=n-alpha;</p><p>　　h1=sin(wc*m)./(pi*m); %求線性低通濾波器的單位脈沖響應(yīng)</p><p>　　hd=ideal_lp(wc,M); %求系統(tǒng)單位脈沖響應(yīng)hd(n)</p><p>　　w_black=(blackman(M))

35、9;; %根據(jù)階數(shù)M求blackman窗</p><p>　　h=hd.*w_black; %對(duì)FIR濾波器低通沖激響應(yīng)加窗，h為低通濾波器的沖擊相應(yīng)</p><p>　　[db,mag,pha,w]=freqz_m2(h,[1]) % 計(jì)算實(shí)際濾波器幅度響應(yīng)</p><p>　　delta_w=2*pi/1000;</p&

36、gt;<p>　　Ap1=-(min(db(1:1:wp/delta_w+1))) %實(shí)際通帶紋波</p><p>　　As1=-round(max(db(ws/delta_w+1:1:501))),</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>

37、;　　%%%%%%%%%%%%%%%% 漢寧窗 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p><b>　　clc,clear</b></p><p><b>　　fb=1000;</b></p><p><b>　　fc=1200;</b></p><

38、;p>　　fs=4000; %對(duì)信號(hào)的采樣頻率為4000</p><p>　　ws=2*pi*fc/fs; %t通帶截至頻率</p><p>　　wp=2*pi*fb/fs %阻帶截至頻率</p><p>　　tr_width=ws-

39、wp; %求過渡帶</p><p>　　M=ceil(6.2*pi/tr_width)+1 %利用Hanning窗的過渡帶=(6.2pi/M),求窗寬M</p><p>　　n=0:1:M-1;</p><p>　　wc=(ws+wp)/2; %求通帶截止頻率wc</p

40、><p>　　hd=ideal_lp(wc,M); %求單位脈沖響應(yīng)hd(n)</p><p>　　w_han=(hanning(M))';</p><p>　　h=hd.*w_han; %求FIR濾波器低通沖激響應(yīng)h(n)</p><p>　　[db,mag,pha

41、,w]=freqz_m2(h,[1]); % 計(jì)算實(shí)際濾波器幅度響應(yīng)</p><p>　　delta_w=2*pi/1000;</p><p>　　Ap2=-(min(db(1:1:wp/delta_w+1))) %實(shí)際通帶紋波</p><p>　　As2=-round(max(db(ws/delta_w+1:1:501))) %%</p>&l

42、t;p>　　%%%%%%%%%%%%%%%%% 海明窗 %%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p><b>　　clc,clear</b></p><p><b>　　fp=1000;</b></p><p><b>　　fc=1200;</b></p>

43、<p>　　fs=4000; %對(duì)信號(hào)的采樣頻率為4000</p><p>　　ws=2*pi*fc/fs; %t通帶截至頻率</p><p>　　wp=2*pi*fp/fs %阻帶截至頻率</p><p>　　tr_w

44、idth=ws-wp; %求過渡帶</p><p>　　M=ceil(6.6*pi/tr_width)+1 %利用Hamming窗的過渡帶=(6.6pi/M),求窗寬M</p><p>　　n=0:1:M-1;</p><p>　　wc=(ws+wp)/2; %求通帶截

45、止頻率wc</p><p>　　hd=ideal_lp(wc,M); %求單位脈沖響應(yīng)hd(n)</p><p>　　w_ham=(hamming(M))'; %海明窗</p><p>　　h=hd.*w_ham; %求FIR濾波器低通沖激響應(yīng)h(n) 實(shí)

46、際單位脈沖響應(yīng)</p><p>　　[db,mag,pha,w]=freqz_m2(h,[1]); % 計(jì)算實(shí)際濾波器幅度響應(yīng)</p><p>　　delta_w=2*pi/1000;</p><p>　　Ap3=-(min(db(1:1:wp/delta_w+1))) %實(shí)際通帶紋波</p><p>　　As3=-round(m

47、ax(db(ws/delta_w+1:1:501)))</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%% 凱塞窗 %%%%%%%%%%%%%%%%%%</p><p><b>　　clc,clear</b></p><p>　　fp=1000; </p><p><b>　　fc

48、=1200;</b></p><p>　　fs=4000; %對(duì)信號(hào)的采樣頻率為4000</p><p>　　ws=2*pi*fc/fs %t通帶截至頻率</p><p>　　wp=2*pi*fp/fs %阻帶截至頻率&

49、lt;/p><p>　　tr_width=ws-wp %求過渡帶</p><p>　　M=ceil(11*pi/tr_width)+1 %利用blackman窗的過渡帶=(11pi/M),求窗寬M</p><p>　　n=0:1:M-1;</p><p>　　wc=(ws+wp)/2;

50、 %求通帶截止頻率wc</p><p>　　hd=ideal_lp(wc,M); %求系統(tǒng)單位脈沖響應(yīng)hd(n)</p><p>　　w_kai=(kaiser(M,10.2256))'; %根據(jù)階數(shù)M求blackman窗 理想的(100-8.7)*0.112 </p><

51、p>　　h=hd.*w_kai; %對(duì)FIR濾波器低通沖激響應(yīng)加窗，h為低通濾波器的沖擊相應(yīng)</p><p>　　[db,mag,pha,w]=freqz_m2(h,[1]); % 計(jì)算實(shí)際濾波器幅度響應(yīng)</p><p>　　delta_w=2*pi/1000;</p><p>　　Ap4=-(min(db(1

52、:1:wp/delta_w+1))) %實(shí)際通帶紋波</p><p>　　As4=-round(max(db(ws/delta_w+1:1:501)))</p><p>　　a=7.9856;L=(100-8)/(2.285*tr_width) </p><p><b>　　%產(chǎn)生干擾序列</b></p><p>

53、<b>　　%產(chǎn)生隨機(jī)白噪聲</b></p><p><b>　　f1=600;</b></p><p><b>　　f2=1800;</b></p><p>　　W1=f1*(2*pi);</p><p>　　W2=f2*(2*pi);</p><p>

54、;　　t=0:(1/4000):0.05;</p><p>　　m=sin(W1.*t)+cos(W2.*t);</p><p><b>　　N=10000;</b></p><p>　　M=abs(fft(m,N));</p><p>　　k=filter(h,1,m);</p><p>　　%

55、k=conv(m,h2);</p><p>　　K=abs(fft(k,N));</p><p><b>　　figure();</b></p><p>　　subplot(2,2,1),plot(m);ylabel('m');title('輸入信號(hào)');</p><p>　　subplo

56、t(2,2,2),stem(M(1:ceil(N/2)));ylabel('fft(m)');</p><p>　　title('輸入信號(hào)幅頻響應(yīng)');</p><p>　　subplot(2,2,3),plot(k);ylabel('k');title('凱塞窗 通過低通濾波器');;</p><p>

57、;　　subplot(2,2,4),stem(K(1:ceil(N/2)));ylabel('fft(k)');</p><p>　　title('通過低通濾波器后幅頻響應(yīng)');</p><p><b>　　擴(kuò)展：</b></p><p>　　1. clf,clc,clear</p><p&g

58、t;　　wph=0.65*pi;</p><p>　　wpl=0.35*pi;</p><p>　　wsth=0.83*pi;</p><p>　　wstl=0.17*pi;</p><p>　　tr_width=wsth-wph;</p><p>　　Mham=ceil(6.6*pi/tr_width)+1 %利

59、用Hamming窗的過渡帶=(6.6pi/M),求窗寬M</p><p>　　n=0:1:Mham-1;</p><p>　　[hd]=ideal_bp(wpl,wstl,wph,wsth,Mham); %求單位脈沖響應(yīng)hd(n)</p><p>　　w_ham=(hamming(Mham))';</p><p>　　h=hd.*w

60、_ham; %求FIR濾波器低通沖激響應(yīng)h(n)</p><p>　　[H,w]=freqz(h,1,1000,'whole');</p><p>　　db=20*log10((abs(H))/max(abs(H))); %求FIR濾波器頻響的dB值</p><p>　　delta_w=2*pi/1000;</p>&l

61、t;p>　　rp=-min(db(1:1:(wpl/delta_w+1))) %求通帶最大衰減</p><p>　　as=-max(db(wstl/delta_w+1:1:501)) %求阻帶最小衰減</p><p>　　Horg=fft(hd,1000);</p><p>　　db_org=20*log10((abs(Horg))

62、/max(abs(Horg))); %求FIR濾波器頻響的dB值</p><p><b>　　2. </b></p><p>　　clear all;wpl=0.4*pi;wph=0.6*pi;wsl=0.2*pi;wsh=0.8*pi;</p><p>　　tr_width=min((wpl-wsl)),(wsh-wph);</p&g

63、t;<p>　　N=ceil(11*pi/tr_width)+1</p><p>　　n=0:1:N-1;</p><p>　　wcl=(wsl+wpl)/2;</p><p>　　wch=(wsh+wph)/2;</p><p>　　hd=ideal_bp1(wcl,wch,N);</p><p>　　

64、w_bman=(blackman(N))';</p><p>　　h=hd.*(w_bman);</p><p>　　[db,mag,pha,w]=freqz_m2(h,[1]); % 計(jì)算實(shí)際濾波器幅度響應(yīng)</p><p>　　delta_w=2*pi/1000;</p><p>　　Ap=-(min(db(wpl/delta_

65、w+1:1:wph/delta_w+1))) %實(shí)際通帶紋波</p><p>　　As=-round(max(db(wsh/delta_w+1:1:501)))</p><p>　　函數(shù)：1 function [hd]=ideal_lp(wc,M)</p><p>　　alpha=(M-1)/2;</p><p>　　n=0:1:(M-1

66、);</p><p>　　m=n-alpha+eps;</p><p>　　hd=sin(wc*m)./(pi*m);%求線性低通濾波器的單位脈沖響應(yīng)</p><p>　　2 function [hd]=ideal_bp(wpl,wstl,wph,wsth,M)</p><p>　　wcl=(wpl+wstl)/2;</p>&

67、lt;p>　　wch=(wph+wsth)/2;</p><p>　　mid=(M-1)/2;</p><p>　　n=0:(M-1);</p><p><b>　　m=n-mid;</b></p><p>　　for i=1:length(m)</p><p>　　if m(i)==0&l

68、t;/p><p><b>　　hd(i)=1;</b></p><p><b>　　else</b></p><p>　　hd(i)=(sin(wch*m(i))-sin(wcl*m(i)))/(pi*m(i));</p><p><b>　　end</b></p>

69、<p><b>　　end</b></p><p>　　3：function [db,mag,pha,w]=freqz_m2(b,a)</p><p>　　[H,w]=freqz(b,a,1000,'whole');</p><p>　　H=(H(1:1:501))';</p><p>