數(shù)字通信原理課程設(shè)計(jì)---線性相位fir 數(shù)字濾波器的設(shè)計(jì)

1、<p>　　課 程 設(shè) 計(jì) 報(bào) 告</p><p>　　課程名稱： 數(shù)字通信課程設(shè)計(jì) </p><p>　　設(shè)計(jì)名稱： 線性相位FIR 數(shù)字濾波器的設(shè)計(jì) </p><p>　　課 程 設(shè) 計(jì) 任 務(wù) 書(shū)</p><p>　　設(shè)計(jì)名稱： 線性相位FIR 數(shù)字濾波器的設(shè)計(jì)

2、 </p><p>　　課 程 設(shè) 計(jì) 學(xué) 生 日 志</p><p>　　課 程 設(shè) 計(jì) 考 勤 表</p><p>　　課 程 設(shè) 計(jì) 評(píng) 語(yǔ) 表</p><p>　　線性相位FIR數(shù)字濾波器的設(shè)計(jì)</p><p><b>　　設(shè)計(jì)目的和意義</b><

3、/p><p>　　數(shù)字濾波器具有高精度、高可靠性、可程控改變特性或復(fù)用、便于集成等優(yōu)點(diǎn)。數(shù)字濾波器在語(yǔ)言信號(hào)處理、圖像信號(hào)處理、醫(yī)學(xué)生物信號(hào)處理以及其他應(yīng)用領(lǐng)域都得到了廣泛應(yīng)用。 </p><p>　　數(shù)字濾波器有低通、高通、帶通、帶阻和全通等類型。它可以是時(shí)不變的或時(shí)變的、因果的或非因果的、線性的或非線性的。應(yīng)用最廣的是線性、時(shí)不變數(shù)字濾波器，以及FIR濾波器。</p><

4、;p>　　有限長(zhǎng)單位沖激響應(yīng)（FIR）濾波器，是數(shù)字信號(hào)處理系統(tǒng)中最基本的元件，它可以在保證任意幅頻特性的同時(shí)具有嚴(yán)格的線性相頻特性，同時(shí)其單位抽樣響應(yīng)是有限長(zhǎng)的，因而濾波器是穩(wěn)定的系統(tǒng)。因此，F(xiàn)IR濾波器在通信、圖像處理、模式識(shí)別等領(lǐng)域都有著廣泛的應(yīng)用。</p><p>　　有限長(zhǎng)單位沖激響應(yīng)（FIR）濾波器有以下特點(diǎn)： </p><p>　　(1) 系統(tǒng)的單位沖激響應(yīng)h (n)在

5、有限個(gè)n值處不為零； </p><p>　　(2) 系統(tǒng)函數(shù)H(z)在|z|>0處收斂，極點(diǎn)全部在z = 0處（因果系統(tǒng)）； </p><p>　　(3) 結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒(méi)有輸出到輸入的反饋，但有些結(jié)構(gòu)中（例如頻率抽樣結(jié)構(gòu)）也包含有反饋的遞歸部分。</p><p>　　本次課程設(shè)計(jì)的意義在于通過(guò)完成對(duì)數(shù)字濾波器的設(shè)計(jì)而讓我進(jìn)一步的了解和學(xué)習(xí)FIR數(shù)字

6、濾波器的線性相位所具有的特點(diǎn)；濾波器的窗函數(shù)的選取與設(shè)計(jì)性能指標(biāo)的相互關(guān)系；高通和低通濾波器的濾波特點(diǎn)，包括幅頻特性和相頻特性等與數(shù)字濾波器的設(shè)計(jì)相關(guān)的知識(shí)。</p><p><b>　　設(shè)計(jì)原理</b></p><p>　　用窗函數(shù)設(shè)計(jì)法設(shè)計(jì)FIR數(shù)字濾波器的基本原理是用一定寬度窗函數(shù)截取無(wú)限脈沖響應(yīng)序列獲得有限長(zhǎng)的脈沖響應(yīng)序列。 </p><p

7、>　　FIR濾波器的設(shè)計(jì)一般是先給出所要求的理想的濾波器頻率響應(yīng)，要求設(shè)計(jì)一個(gè)FIR濾波器頻率響應(yīng)</p><p><b>　?。?）</b></p><p>　　來(lái)逼近。但是設(shè)計(jì)是在時(shí)域進(jìn)行的，因而先由的傅里葉反變換導(dǎo)出，即</p><p><b>　?。?）</b></p><p>　　由

8、于是矩形頻率特性，故一定是無(wú)限長(zhǎng)的序列，且是非因果的，而我們要設(shè)計(jì)的是FIR濾波器，其h(n)必然是有限長(zhǎng)的，所以要用有限長(zhǎng)的h(n)來(lái)逼近無(wú)限長(zhǎng)的。最有效的方法是截取或者說(shuō)，是用一個(gè)有線長(zhǎng)度的窗函數(shù)序列來(lái)截取，即</p><p><b>　?。?）</b></p><p>　　因而窗函數(shù)序列的形狀及長(zhǎng)度的選擇就很關(guān)鍵。</p><p>　　在

9、選取窗函數(shù)時(shí)，一般希望窗函數(shù)滿足一下兩項(xiàng)要求：窗譜主瓣盡可能地窄，以獲得較陡的過(guò)渡帶；盡量減少窗譜的最大主瓣的相對(duì)幅度，也就是能量盡量集中于主瓣，這樣使肩峰和波紋減小，就可增大阻帶的衰減。但是這兩項(xiàng)要求是不能同時(shí)得到滿足的，往往是增加主瓣寬度以換取對(duì)旁瓣的抑制。因而選用不同形狀的窗函數(shù)都是為了使得到平坦的通帶幅度響應(yīng)和較小的阻帶波紋（也就是加大阻帶衰減）。</p><p>　　常用的窗函數(shù)有矩形窗、Bartlet

10、t窗、Hanning窗、Hamming窗、Blackman窗、Kaiser窗。</p><p>　　下表歸納了以上幾種窗的主要性能，供設(shè)計(jì)FIR濾波器時(shí)參考。</p><p><b>　　表1</b></p><p>　　由上表可知：不同的窗有不同的性質(zhì)：不同的窗函數(shù)，產(chǎn)生泄漏的大小不一樣，頻率分辨能力也不一樣。信號(hào)的截?cái)喈a(chǎn)生了能量泄漏，而用F

11、FT算法計(jì)算頻譜又產(chǎn)生了柵欄效應(yīng)，從原理上講這兩種誤差都是不能消除的，但是我們可以通過(guò)選擇不同的窗函數(shù)對(duì)它們的影響進(jìn)行抑制。(矩形窗主瓣窄，旁瓣大，頻率識(shí)別精度最高，幅值識(shí)別精度最低；布萊克曼窗主瓣寬，旁瓣小，頻率識(shí)別精度最低，但幅值識(shí)別精度最高)。</p><p>　　實(shí)際設(shè)計(jì)過(guò)程中，上述幾種窗函數(shù)可以根據(jù)對(duì)濾波器的過(guò)渡帶寬度和阻帶衰減的要求，適當(dāng)選取窗函數(shù)的類型和長(zhǎng)度N，以得到比較滿意的設(shè)計(jì)效果。如何根據(jù)濾波

12、器長(zhǎng)度N的奇偶性，選擇h(n)的奇偶對(duì)稱性則是另外一個(gè)需要考慮的問(wèn)題。線性相位實(shí)系數(shù)FIR濾波器按其N值奇偶和h(n)的奇偶對(duì)稱性，可以分為四種，它們具有不同的幅頻和相位特性：</p><p>　　1．h(n)為偶對(duì)稱，N為奇數(shù)：</p><p>　　H(e)=[h()+ （4）</p><p>　　它的幅度是關(guān)于點(diǎn)成偶對(duì)稱。</p>&l

13、t;p>　　2．h(n)偶對(duì)稱，N為偶數(shù)：</p><p>　　H(e)=[h()+ （5）</p><p>　　它的幅度是關(guān)于成奇對(duì)稱，處有零點(diǎn)，所以它不適合于作高通濾波器。</p><p>　　3．h(n)為奇對(duì)稱，N為奇數(shù)</p><p>　　H(e)={ （6）</p><p>　　

14、它的幅度是關(guān)于點(diǎn)成奇對(duì)稱, H(e)在處都有零點(diǎn)。因此它不適合于低通和高通。</p><p>　　4．h(n)為奇對(duì)稱，N為偶數(shù)</p><p>　　H(e)={ + （7）</p><p>　　它的幅度是關(guān)于點(diǎn)成奇對(duì)稱。h(e)在處有零點(diǎn)。因此它不適于低通。</p><p><b>　　詳細(xì)設(shè)計(jì)步驟</b><

15、/p><p>　　窗函數(shù)法設(shè)計(jì)線性相位FIR濾波器可以按如下步驟：</p><p>　　1．根據(jù)數(shù)字濾波器的性能指標(biāo)，確定各臨界頻率的歸一化頻率；</p><p>　　2．由所需低通濾波器的過(guò)渡帶球理想低通濾波器的截止頻率，并得出濾波器頻率響應(yīng)的傅里葉反變換導(dǎo)出；</p><p>　　3．求窗函數(shù)。由阻帶衰減確定窗形狀，由過(guò)渡帶寬確定單位脈沖響應(yīng)

16、長(zhǎng)度N，根據(jù)四種線性相位FIR濾波器特性得出，只能采用偶對(duì)稱單位沖激響應(yīng)，這是一般通帶濾波器（包括低通、高通、帶通、帶阻）所需求的單位沖激響應(yīng)，另外，對(duì)高通濾波器來(lái)說(shuō)N只能取奇數(shù)（N=2m+1，m=0，1，…）。</p><p>　　4．求h(n)。由窗函數(shù)表達(dá)式確定FIR濾波器的h(n)；</p><p>　　5．求H(e)。根據(jù)求得的h(n)，用傅里葉變換求得其頻率響應(yīng)H(e)，分析他

17、的頻率特性，若不滿足要求，可適當(dāng)改變窗函數(shù)形式或長(zhǎng)度N,重復(fù)上述過(guò)程，直至得到滿意結(jié)果。</p><p>　　以上是窗函數(shù)法設(shè)計(jì)的一般步驟，當(dāng)采用Matlab函數(shù)設(shè)計(jì)時(shí)就沒(méi)有這么麻煩了，因?yàn)橹苯诱{(diào)用Matlab中的函數(shù)，再輸入相關(guān)參數(shù)就可以得到一個(gè)濾波器傳輸函數(shù)系數(shù)的序列。在程序中主要調(diào)用的函數(shù)是fir1，它是用來(lái)設(shè)計(jì)傳統(tǒng)的低通、高通、帶通、帶阻等類型的FIR濾波器，其調(diào)用的一般格式如下：</p>

18、<p>　　b = fir1(N,Wn)；</p><p>　　b = fir1(N,Wn,‘high’,win)；</p><p>　　參數(shù)說(shuō)明：N：階次，濾波器長(zhǎng)度為N＋1；Wn：通帶截止頻率，為歸一化過(guò)后的頻率，其值在0～1之間；High：濾波器的類型，可以是低通（low）、高通（high）、帶阻（stop）等；win：窗函數(shù)，在不指定窗函數(shù)類型的情況下，fir1自動(dòng)選擇H

19、amming窗； b：濾波器系數(shù)。</p><p>　　特別地，在設(shè)計(jì)高通濾波器時(shí)，由于N的取值必須為奇數(shù)，于是調(diào)用了mod函數(shù)，其作用是求余數(shù)，函數(shù)中的語(yǔ)句為：N=N0+mod(N0+1,2)，其作用是：當(dāng)N0是奇數(shù)時(shí)，余數(shù)為零，此時(shí)N的取值即為N0，當(dāng)N0是偶數(shù)時(shí)，余數(shù)為一，此時(shí)N的取值為N0+1。 </p><p>　　在編寫(xiě)程序時(shí)，有如下五個(gè)步驟：</p><

20、;p>　　1.根據(jù)所要求的濾波器的性能指標(biāo)，求得過(guò)渡帶寬以及濾波器的截止頻率，過(guò)渡帶寬為通帶截止頻率和阻帶截止頻率的差值，即，濾波器的截止頻率為；</p><p>　　2．由過(guò)渡帶寬求取濾波器的窗口長(zhǎng)度，若所要求的過(guò)渡帶寬為：，以海明窗為例，由于海明窗過(guò)渡帶寬滿足，則N的取值為：；</p><p>　　3．根據(jù)濾波器性能指標(biāo)中的阻帶衰減，選擇適當(dāng)?shù)拇昂瘮?shù)，例如要求阻帶衰減不小于-50

21、dB時(shí)，差表1可知，海明窗的最小衰減為-53dB，因此可以選擇海明窗作為窗函數(shù)；</p><p>　　選擇好窗函數(shù)后，調(diào)用fir1函數(shù)，求得濾波系數(shù)，也即是沖激響應(yīng)h(n)；</p><p>　　5．利用求出的h(n)，求取濾波器的頻率響應(yīng)，然后再畫(huà)出圖形。</p><p><b>　　設(shè)計(jì)結(jié)果及分析</b></p><p&

22、gt;　?。ㄒ唬IR高通濾波器</p><p>　　高通濾波器的通帶截止頻率為0.5，阻帶截止頻率為0.3，通帶衰減頻率為1dB，阻帶衰減頻率為40dB。</p><p>　　Matlab程序如下：</p><p>　　wp=0.5*pi; ws=0.3*pi; </p><p>　　width=wp-ws; <

23、/p><p>　　M=2*ceil(3.3*pi/width) </p><p>　　N=N0+mod(N0+1,2)-4; </p><p>　　wc=(ws+wp)/2; </p><p>　　win=hamming(N); </p><p>　　b=fir1(N-1

24、,wc/pi,'high',win);</p><p>　　[hf w]=freqz(b,1,1024); </p><p>　　subplot(2,1,1); </p><p>　　plot(w/pi,20*log10(abs(hf)));</p><p>　　xlabel('頻率(w/pi)'); <

25、;/p><p>　　ylabel('幅度(dB)');</p><p>　　title('幅頻特性曲線');grid</p><p>　　subplot(2,1,2); </p><p>　　plot(w/pi,180/pi*unwrap(angle(hf)));</p><p>　　xla

26、bel('頻率(w/pi)'); </p><p>　　ylabel('相位(degrees)');</p><p>　　title('相頻特性曲線');grid</p><p>　　輸出頻率特性曲線如下圖所示：</p><p><b>　　圖1</b></p>

27、<p><b>　　圖2</b></p><p>　　因?yàn)樽鑾p為-40dB，而由表一可得，hanning窗的阻帶最小衰減為-44dB，于是選取hanning窗作為窗口函數(shù)，計(jì)算得出窗口長(zhǎng)度為32，由于在設(shè)計(jì)高通濾波器時(shí)，N值只能為奇數(shù)，所以取N=33，繪制出的頻率特性與要求的性能指標(biāo)不符，改變N的數(shù)值，直到N=29時(shí)，得到的頻率特性曲線才相近似，但是任然不能完全匹配，繪制的

28、圖形如圖1所示。于是改選衰減為-53dB的hamming窗作為窗口函數(shù)，最初得出的窗口長(zhǎng)度為34，取N=35，但是所得頻率特性不滿足性能指標(biāo)，改變N的值為31后，達(dá)到設(shè)計(jì)要求，繪制的圖形如圖2所示，在阻帶截止頻率處幅度為-41dB，在通帶截止頻率處幅度為-1dB，符合阻帶衰減為-40dB的指標(biāo)要求。</p><p><b>　?。ǘ?、低通濾波器</b></p><p&g

29、t;　　低通濾波器的通帶截止頻率為0.2，阻帶截止頻率為0.4，通帶衰減頻率為1dB，阻帶衰減頻率為50dB。</p><p>　　Matlab程序如下：</p><p>　　wp=0.2*pi; ws=0.4*pi; </p><p>　　width=ws-wp; </p><p>　　N=2*ceil(3.3*pi/wi

30、dth); </p><p>　　wc=(ws+wp)/2; </p><p>　　win=hamming(N); </p><p>　　b=fir1(N-1,wc/pi,'low',win);</p><p>　　[hf w]=freqz(b,1,1024); </p&g

31、t;<p>　　subplot(2,1,1); </p><p>　　plot(w/pi,20*log10(abs(hf)));</p><p>　　xlabel('頻率(w/pi)'); </p><p>　　ylabel('幅度(dB)');</p><p>　　title('幅頻特性

32、曲線');grid</p><p>　　subplot(2,1,2); </p><p>　　plot(w/pi,180/pi*unwrap(angle(hf)));</p><p>　　xlabel('頻率(w/pi)'); </p><p>　　ylabel('相位(degrees)');</

33、p><p>　　title('相頻特性曲線');grid</p><p>　　輸出頻率特性如下圖所示：</p><p><b>　　圖3</b></p><p><b>　　圖4</b></p><p>　　因?yàn)樽鑾p為-50dB，而由表一可得，hanning窗

34、的阻帶最小衰減為-44dB，于是選取hanning窗作為窗口函數(shù)，計(jì)算得出窗口長(zhǎng)度為32，而繪制出的頻率特性與要求的性能指標(biāo)不能完全匹配，繪制的圖形如圖3所示。于是改選衰減為-53dB的hamming窗作為窗口函數(shù)，得出的窗口長(zhǎng)度為34，達(dá)到設(shè)計(jì)要求，繪制的圖形如圖4所示，在阻帶截止頻率處幅度為-51dB，在通帶截止頻率處幅度為-1dB，符合阻帶衰減為-50dB的指標(biāo)要求。</p><p><b>　　

35、體會(huì)</b></p><p>　　本次課程設(shè)計(jì)是用窗函數(shù)分別設(shè)計(jì)一個(gè)高通和低通濾波器，雖然是《數(shù)字通信》課程設(shè)計(jì)，但是涉及到得的知識(shí)卻主要是《數(shù)字信號(hào)處理》中的內(nèi)容，以至于在整個(gè)設(shè)計(jì)過(guò)程中都沒(méi)有怎么參考《現(xiàn)代通信原理》這本書(shū)。</p><p>　　在最開(kāi)始看到這個(gè)題目的時(shí)候以為這個(gè)設(shè)計(jì)題目很簡(jiǎn)單，因?yàn)槲覀冊(cè)?、學(xué)習(xí)過(guò)相關(guān)知識(shí)，并且在實(shí)驗(yàn)課中，也不止一次的進(jìn)行了設(shè)計(jì)，然而，經(jīng)過(guò)認(rèn)真

36、分析過(guò)后發(fā)現(xiàn)，題目并不像想象中的那么簡(jiǎn)單，尤其是對(duì)于我來(lái)說(shuō)，因?yàn)閷?duì)《數(shù)字信號(hào)處理》這本書(shū)中的內(nèi)容并沒(méi)有完全熟悉和掌握，而在實(shí)驗(yàn)課中也沒(méi)有認(rèn)真分析過(guò)實(shí)驗(yàn)程序，同時(shí)，對(duì)設(shè)計(jì)要用到的Matlab的知識(shí)也遺忘了很多，所以在本次設(shè)計(jì)過(guò)程中對(duì)用窗函數(shù)法設(shè)計(jì)濾波器的相關(guān)知識(shí)進(jìn)行了重新學(xué)習(xí)，并且對(duì)在設(shè)計(jì)中遇到的問(wèn)題虛心請(qǐng)教和積極討論。最終，伴隨著課程設(shè)計(jì)的完成，關(guān)于濾波器的知識(shí)和Matlab的使用也逐漸從淡忘到熟悉。因此我覺(jué)得從這方面來(lái)說(shuō)本次課程設(shè)計(jì)給

37、了我很大的收獲，同時(shí)我也認(rèn)識(shí)到自身的學(xué)習(xí)還有很多的不足之處，我重新了解的知識(shí)也僅僅限于設(shè)計(jì)所需要的，而對(duì)其它方面的知識(shí)依然不熟悉，這無(wú)疑會(huì)讓我以后更加注重對(duì)書(shū)本上的知識(shí)的學(xué)習(xí)。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　[1] 程佩青，數(shù)字信號(hào)處理教程（第三版），清華大學(xué)出版社，2007</p><p>　　[2] 張威