畢業(yè)論文---基于matlab的fir數(shù)字濾波器設(shè)計(jì)及其軟件實(shí)現(xiàn)

1、<p><b>　　摘 要</b></p><p>　　在這個(gè)具有挑戰(zhàn)性的時(shí)代，把計(jì)算機(jī)充分運(yùn)用到教學(xué)及工程運(yùn)算中，雖然具有重要的意義，而隨著計(jì)算機(jī)技術(shù)的發(fā)展，計(jì)算機(jī)軟件在工程設(shè)計(jì)領(lǐng)域應(yīng)用越來越廣。</p><p>　　本設(shè)計(jì)研究是基于MATLAB的FIR數(shù)字濾波器的設(shè)計(jì)，是基于應(yīng)用計(jì)算機(jī)軟件編程。首先了解濾波器的基本工作原理，針對(duì)FIR濾波器，有各種不同的

2、方法。FIR濾波器一般要求信號(hào)具有線性相位，同時(shí)由于脈沖響應(yīng)無限長(zhǎng)，因此具有永遠(yuǎn)穩(wěn)定的特性。根據(jù)這個(gè)特點(diǎn)，。本文主要對(duì)FIR濾波器采用了窗函數(shù)法，最優(yōu)設(shè)計(jì)法，最小二乘設(shè)計(jì)法，升余弦函數(shù)設(shè)計(jì)法和任意響應(yīng)設(shè)計(jì)法來設(shè)計(jì)FIR濾波器，同時(shí)，還對(duì)濾波器進(jìn)行階數(shù)評(píng)估。濾波是信號(hào)處理中最基本有極為重要的技術(shù)，利用濾波器技術(shù)可以從復(fù)雜的信號(hào)中提取出所需要的信號(hào)，抑制不需要的信號(hào)。絕大多數(shù)傳感器輸出的信號(hào)，在使用過程中，都必須進(jìn)行濾波，所以濾波器是具有一

3、定傳輸選擇特性的，對(duì)信號(hào)進(jìn)行加工處理的裝置，它允許輸入信號(hào)中的一些成分通過，抑制或衰減另一些成分。其功能是將輸入信號(hào)變換為人們所需要的輸入信號(hào)。</p><p>　　本設(shè)計(jì)主要給出了FIR數(shù)字濾波器的基本特性和設(shè)計(jì)方法。</p><p>　　關(guān)鍵詞：MATLAB 特性 FIR濾波器 設(shè)計(jì)</p><p><b>　　Abstract</

4、b></p><p>　　In this challenging era, the full use of computers to teaching and engineering operations, although of great significance, and with the development of computer technology, computer software in

5、 engineering design applications more widely. The design study is based on MATLAB for FIR digital filter design, is based on the application of computer software programming. First understand the basic working principle

6、 of the filter for FIR filters, a variety of different ways. General requirements for FI</p><p><b>　　朗讀字典</b></p><p><b>　　目錄</b></p><p><b>　　摘 要I</

7、b></p><p>　　AbstractII</p><p><b>　　第一章緒論1</b></p><p>　　1.1 數(shù)字濾波器定義1</p><p>　　1.2 數(shù)字濾波器的分類1</p><p>　　1.3 實(shí)際濾波器的設(shè)計(jì)指標(biāo)2</p><p&

8、gt;　　1.4 幾種常見的特殊濾波器2</p><p>　　第二章FIR濾波器特性3</p><p>　　2.1 FIR濾波器簡(jiǎn)介3</p><p>　　2.2 FIR濾波器的線性相位特性3</p><p>　　2.2.1 FIR濾波器的第一類線性相位4</p><p>　　2.2.2 FIR濾波

9、器的第二類線性相位6</p><p>　　2.2.3 線性相位FIR濾波器的零點(diǎn)特性8</p><p>　　第三章FIR濾波器的ＭＡＴＬＡＢ設(shè)計(jì)9</p><p>　　3.1 用窗函數(shù)法設(shè)計(jì)FIR數(shù)字濾波器9</p><p>　　3.1.1 FIR數(shù)字濾波器窗函數(shù)設(shè)計(jì)法9</p><p>　　3.1.

10、2 加窗的線性相位FIR數(shù)字濾波器設(shè)計(jì)函數(shù)fir1和fir212</p><p>　　3.2 FIR濾波器的最優(yōu)設(shè)計(jì)法15</p><p>　　3.2.1 firs函數(shù)15</p><p>　　3.2.2 remez函數(shù)16</p><p>　　3.3 FIR濾波器最小二乘設(shè)計(jì)法17</p><p>

11、;　　3.3.1 fircls函數(shù)17</p><p>　　3.3.2 fircls函數(shù)19</p><p>　　3.4 FIR濾波器升余弦函數(shù)設(shè)計(jì)法20</p><p>　　3.5 FIR濾波器的任意響應(yīng)設(shè)計(jì)法21</p><p>　　第四章　FIR 濾波器階數(shù)估計(jì)24</p><p>　　4.1

12、 kaiserord函數(shù)24</p><p>　　4.2 remezord函數(shù)25</p><p>　　第五章 小結(jié)27</p><p><b>　　參考文獻(xiàn)28</b></p><p><b>　　致 謝29</b></p><p><b>　　緒

13、論</b></p><p>　　當(dāng)前，通信技術(shù)高速發(fā)展，業(yè)務(wù)范圍不斷擴(kuò)大，人們對(duì)產(chǎn)品的需求迅速增長(zhǎng)，濾波器在這些產(chǎn)品電路中就扮演著重要的角色，當(dāng)然數(shù)字濾波器更影響著人們生活的方方面面。數(shù)字濾波器既可以是用數(shù)字硬件裝配成德一臺(tái)完成給定運(yùn)算的專用數(shù)字計(jì)算機(jī)，也可將所需的運(yùn)算編成程序，讓通用計(jì)算機(jī)來執(zhí)行。</p><p>　　MATLAB是“矩陣實(shí)驗(yàn)室”（MATrix LABorat

14、oy）的縮寫，它是由美國(guó)Mathworks公司于1984年正式推出的，是一種以矩陣運(yùn)算為基礎(chǔ)的交互式程序語(yǔ)言，專門針對(duì)科學(xué)、工程計(jì)算及繪圖的需求。MATLAB是功能強(qiáng)大的科學(xué)及工程計(jì)算軟件，它不但具有以矩陣計(jì)算為基礎(chǔ)的強(qiáng)大數(shù)學(xué)計(jì)算和分析功能，而且還具有豐富的可視化圖形表現(xiàn)功能和方便的程序計(jì)算能力。MATLAB的應(yīng)用領(lǐng)域極為廣泛，除數(shù)學(xué)計(jì)算和分析外，還被廣泛地應(yīng)用于自動(dòng)控制、系統(tǒng)仿真、數(shù)字信號(hào)處理、圖形圖像分析、數(shù)理統(tǒng)計(jì)、人工智能、虛擬現(xiàn)

15、實(shí)技術(shù)、通信工程、金融系統(tǒng)等領(lǐng)域，因此，MATLAB是面向21世紀(jì)的計(jì)算機(jī)程序設(shè)計(jì)及科學(xué)計(jì)算語(yǔ)言。</p><p>　　1.1 數(shù)字濾波器定義</p><p>　　數(shù)字濾波器通常是指一個(gè)有限精度算法實(shí)現(xiàn)的離散線性是不變系統(tǒng)。通常用的數(shù)字濾波器一般屬于選頻濾波器，可以用下式表示：</p><p>　　稱為濾波器幅頻響應(yīng)， 稱為濾波器的相頻響應(yīng)。幅頻響應(yīng)表示信號(hào)通過該濾

16、波器后各頻率成分的衰減情況，而相頻響應(yīng)反映各頻率成分通過濾波器后在時(shí)間上的延時(shí)情況。</p><p>　　濾波器性能最容易通過它的幅頻響應(yīng)的形狀來描述 。濾波器在某個(gè)頻率的幅度增益決定了濾波器對(duì)此頻率輸入的放大因子，增益可取任意值。增益高的頻率范圍，信號(hào)可以通過，稱之為濾波器的通帶，增益低的頻率范圍，濾波器對(duì)信號(hào)有衰減和阻塞作用，稱這位濾波器的阻帶。</p><p>　　1.2 數(shù)字濾波

17、器的分類</p><p>　　數(shù)字濾波器按照不同的分類方法，有許多種類，總體上可分為兩大類。一類是經(jīng)典濾波器，即一般的線性系統(tǒng)濾波器。另一類即所謂的現(xiàn)代濾波器，它的理論建立在隨機(jī)信號(hào)處理的理論基礎(chǔ)上，利用了隨機(jī)信號(hào)內(nèi)部的統(tǒng)計(jì)特性對(duì)信號(hào)進(jìn)行濾波。</p><p>　　從功能上可以分為低通、高通、帶通和帶阻濾波器，根據(jù)數(shù)字濾波器沖激響應(yīng)的時(shí)域特性，可將數(shù)字濾波器分為兩種，即無限長(zhǎng)沖激響應(yīng)（II

18、R）濾波器和有限長(zhǎng)沖激響應(yīng)（FIR）濾波器[3]。</p><p>　　1.3 實(shí)際濾波器的設(shè)計(jì)指標(biāo)</p><p>　　當(dāng)濾波器形狀為非理想時(shí)，要用一些參數(shù)指標(biāo)來描述其關(guān)鍵特性。濾波器的通帶定義了濾波器允許通過的頻率范圍。在阻帶內(nèi)，濾波器對(duì)信號(hào)嚴(yán)重衰減。 和 分別稱為通帶截止頻率（或通帶上限頻率）和阻帶截止頻率（或阻帶下限頻率）。參數(shù) 定義了通帶波紋，及濾波器通帶內(nèi)偏移單位增益的最大值。

19、參數(shù) 定義了阻帶波紋，及濾波器阻帶內(nèi)偏離零增益的最大值。參數(shù) 定義了過渡帶寬度，即阻帶下限和通帶上限之間的距離， = 。過渡帶一般是單調(diào)下降的，通帶內(nèi)和阻帶內(nèi)允許的衰減一般用單位dB表示，通帶內(nèi)允許的最大衰減用 表示，阻帶內(nèi)允許的最小衰減用 表示，它們分別定義為</p><p><b>　　dB</b></p><p><b>　　dB</b>&

20、lt;/p><p>　　式中， 是通帶內(nèi)的幅度最大值； 是通帶內(nèi)的幅度最小值， 是阻帶內(nèi)最大值。幅度下降到0.707即 時(shí)，w= ，此時(shí) =3dB，稱 為3dB通帶截止頻率。</p><p>　　1.4 幾種常見的特殊濾波器</p><p>　　全通濾波器、數(shù)字陷波器、梳狀濾波器、最小相位濾波器</p><p>　　第二章FIR濾波器特性<

21、;/p><p>　　2.1 FIR濾波器簡(jiǎn)介</p><p>　　FIR系統(tǒng)只有零點(diǎn)，因此這類系統(tǒng)不像IIR系統(tǒng)那樣易取得比較好的通帶和阻帶衰減特性。但FIR系統(tǒng)有自己突出的優(yōu)點(diǎn)，其一是系統(tǒng)總是穩(wěn)定的，其二是易實(shí)現(xiàn)線性相位，其三是允許設(shè)計(jì)多通帶（或多阻帶）濾波器。</p><p>　　FIR濾波器是指系統(tǒng)的單位沖擊響應(yīng)h（n）僅在有限的范圍內(nèi)有非零值的濾波器。N-1階

22、濾波器的系統(tǒng)函數(shù)H（z）可表示為</p><p>　　H（n）是 的N-1次多項(xiàng)式，它在z平面上有N-1個(gè)零點(diǎn)，原點(diǎn)z=0是N-1階重極點(diǎn)。因此，F(xiàn)IR濾波器永遠(yuǎn)穩(wěn)定。</p><p>　　FIR濾波器的頻率響應(yīng)為</p><p>　　由于H（ ）一般為復(fù)數(shù)，因此，可將其表示成</p><p>　　2.2 FIR濾波器的線性相位特性<

23、/p><p>　　FIR濾波器的頻率響應(yīng)為</p><p>　　由于H（ ）一般為復(fù)數(shù)，因此，可將其表示成</p><p>　　線性相位特性是指 是w的線性函數(shù)，即</p><p>　　式中， 是常數(shù)。此時(shí)通過這一系統(tǒng)的各頻率分量的時(shí)延為相同的常數(shù)，系統(tǒng)的群時(shí)延為</p><p>　　即系統(tǒng)的時(shí)延是一個(gè)與w無關(guān)的常數(shù) ，稱

24、系統(tǒng)H（z）具有嚴(yán)格的線性相位。</p><p>　　由于嚴(yán)格線性相位條件在數(shù)學(xué)上處理較為困難，因此在FIR濾波器設(shè)計(jì)中一般使用廣義線性相位。若一個(gè)離散系統(tǒng)的頻率響應(yīng)H（ ）可以寫為</p><p><b>　　H（ ）= </b></p><p>　　其中， 和 是與w無關(guān)的常數(shù)， 是w的是函數(shù)（可任意取值），則稱為廣義幅頻響應(yīng)， 稱為廣義相

25、頻響應(yīng)。</p><p><b>　　上式中， 為</b></p><p>　　當(dāng) =0時(shí)，稱是第一類線性相位，當(dāng) （一般 ）時(shí)，稱是第二類線性相位。</p><p>　　2.2.1 FIR濾波器的第一類線性相位</p><p>　　FIR濾波器滿足第一類線性相位的條件是</p><p>　　由

26、于h（n）的點(diǎn)數(shù)N可分為奇數(shù)和偶數(shù)兩種情況，因而第一類線性相位又可分為兩類。</p><p>　?。?） h（n）為偶對(duì)稱，N為奇數(shù)</p><p><b>　　顯然</b></p><p>　　滿足第一類線性相位的條件，而</p><p>　　由于N為奇數(shù)，因此n-- 為整數(shù)。由于 關(guān)于w=0， ，2 這些點(diǎn)都是偶對(duì)

27、稱的，因此 關(guān)于w=0， ，2 也是偶對(duì)稱的，這種情況適合設(shè)計(jì)低通、高通、帶通、帶阻濾波器。</p><p>　?。?）h（n）為偶對(duì)稱，N為偶函數(shù)</p><p>　　采用類似的方法，可以得到</p><p>　　其相頻響應(yīng)與N為奇數(shù)時(shí)相同，為</p><p><b>　　而幅頻響應(yīng)為</b></p>&

28、lt;p>　　由于N為偶數(shù)，因此 為0.5的奇數(shù)倍。 關(guān)于w= 奇對(duì)稱，當(dāng)w= 時(shí)， =0，因此， =0，即H（z）在z=-1處必然有一個(gè)零點(diǎn)，而且， 關(guān)于w= 奇對(duì)稱這種情況不適合設(shè)計(jì)高通和帶阻濾波器。</p><p>　　2.2.2 FIR濾波器的第二類線性相位</p><p>　　FIR濾波器滿足第二類相位的條件</p><p>　　h（n）為奇對(duì)稱，

29、N為奇數(shù)</p><p>　　h（n）的中間項(xiàng)h（ ）必須為零，因此</p><p><b>　　此時(shí)的相位為</b></p><p>　　滿足第二類線性相位的條件，而幅頻響應(yīng)為</p><p>　　因?yàn)?在w=0， 處都為0，因此 在w=0， 處為0，也即H（z）在z= 1處都有零點(diǎn)，并且 關(guān)于w=0， 奇對(duì)稱。因此，

30、該類型濾波器不適合設(shè)計(jì)低通、高通和帶阻濾波器，只適合設(shè)計(jì)帶通濾波器。</p><p>　?。?）h（n）為奇對(duì)稱，N為偶數(shù)</p><p><b>　　相位為</b></p><p><b>　　幅頻響應(yīng)為</b></p><p>　　由于 在w=0,2 處為0，因此 在w=0,2 處為0，即H（z

31、）在z=1上有零點(diǎn)，且關(guān)于w=0,2 奇對(duì)稱。</p><p>　　h（n）為奇對(duì)稱有兩種情況，對(duì)于任何頻率都有固定的 相移，因此第二類線性相位FIR濾波器一般用于正交移相或微分器。</p><p>　　2.2.3 線性相位FIR濾波器的零點(diǎn)特性</p><p>　　線性相位FIR濾波器的單位脈沖響應(yīng)具有對(duì)稱特性，即</p><p><

32、;b>　　則</b></p><p>　　將m=N-1-n代入上式，得</p><p>　　因此，其系統(tǒng)函數(shù)具有以下特點(diǎn)：</p><p>　　可以看出，若z= 及z= 也必定是零點(diǎn)。</p><p>　　第三章FIR濾波器的ＭＡＴＬＡＢ設(shè)計(jì)</p><p>　　IIR濾波器設(shè)計(jì)利用了模擬濾波器成熟理

33、論以及設(shè)計(jì)圖表，保持了典型模擬濾波器優(yōu)良的幅頻特性，但設(shè)計(jì)中只考慮了幅頻特性，沒有考慮相位特性，因此濾波器的相位特性通常是非線性的，要得到線性相位必須外加相位矯正網(wǎng)絡(luò)，使濾波器變得復(fù)雜。而FIR濾波器在保證幅頻滿足技術(shù)要求的同時(shí)，更榮易做到嚴(yán)格的線性相位。</p><p>　　FIR濾波器的系統(tǒng)函數(shù)為</p><p>　?。龋ǎ┰趜平面上有N-1個(gè)零點(diǎn)，并且z=0是N-1階重極點(diǎn)，位于r

34、=1的單位圓內(nèi)，系統(tǒng)永遠(yuǎn)穩(wěn)定。</p><p>　　FIR濾波器的設(shè)計(jì)方法有窗函數(shù)法、頻率采樣法和契比雪夫等波紋逼近法</p><p>　　3.1 用窗函數(shù)法設(shè)計(jì)FIR數(shù)字濾波器</p><p>　　3.1.1 FIR數(shù)字濾波器窗函數(shù)設(shè)計(jì)法</p><p>　　通常，理想數(shù)字濾波器的單位取樣序列 是無限長(zhǎng)的，是非因果序列，在物理上是不可

35、實(shí)現(xiàn)的或者難以實(shí)現(xiàn)的。為了可物理實(shí)現(xiàn)，需要對(duì) 截取一段h（n）來近似代表，即 h（n）= </p><p>　　窗函數(shù)主瓣寬度要盡可能窄，以獲取最小的過渡帶；旁瓣相對(duì)值盡可能小，以使得通帶波紋小，并且阻帶衰減大。</p><p>　　（1）常用的窗函數(shù)為矩形窗、巴特利窗 漢寧窗、海明窗、布萊克曼窗、凱澤窗和契比雪夫窗。</p><p><b>　　1）矩形

36、窗</b></p><p><b>　　矩形窗的頻域響應(yīng)為</b></p><p>　　其主瓣寬度為 。調(diào)用w=boxcar（n）函數(shù)</p><p><b>　　2）巴特利窗</b></p><p><b>　　其頻率響應(yīng)為</b></p><

37、p>　　其主瓣寬度為 ，調(diào)用w=bartlett（n）函數(shù)</p><p><b>　　3）漢寧窗</b></p><p><b>　　其頻域表達(dá)式為</b></p><p>　　主瓣寬度為 ，調(diào)用w=hann（n）函數(shù)</p><p><b>　　4）海明窗</b>&l

38、t;/p><p><b>　　其頻域表達(dá)式為</b></p><p>　　主瓣寬度為 ，調(diào)用w=hamming（n）函數(shù)</p><p><b>　　5）布萊克曼窗</b></p><p><b>　　其頻域表示式為</b></p><p>　　其主瓣寬度為

39、 ，調(diào)用w=blackman（n）函數(shù)</p><p><b>　　6）凱澤窗</b></p><p>　　式中， ， 為第一類修正貝塞爾函數(shù)，窗函數(shù)的幅度函數(shù)為</p><p>　　其主瓣寬度為 ，調(diào)用w=Kaiser（n，beta）函數(shù)</p><p><b>　　7）契比雪夫窗</b><

40、/p><p>　　調(diào)用w=chebwin（n，r）</p><p>　　表3.1 窗函數(shù)基本參數(shù)</p><p>　?。?）FIR濾波器的設(shè)計(jì)步驟</p><p>　　1）根據(jù)技術(shù)要求確定待求濾波器的單位取樣響應(yīng)或者頻率響應(yīng) 。</p><p>　　2）根據(jù)對(duì)過渡帶和阻帶衰減的要求，選擇窗函數(shù)形式，并估計(jì)窗口長(zhǎng)度N。<

41、;/p><p>　　3）計(jì)算濾波器的單位取樣響應(yīng)h（n）= ，h（n）即為所求。</p><p>　　4）檢測(cè)所設(shè)計(jì)的濾波器性能，驗(yàn)算技術(shù)指標(biāo)是否滿足條件。</p><p>　　3.1.2 加窗的線性相位FIR數(shù)字濾波器設(shè)計(jì)函數(shù)fir1和fir2</p><p>　　1）fir1。其調(diào)用格式為：</p><p>　　b=

42、fir1（n，Wn）</p><p>　　b=fir1（n，Wn，‘ftype’）</p><p>　　b=fir1（n，Wn，window）</p><p>　　b=fir1（n，Wn，‘ftype’，window）</p><p>　　b=fir1（…,‘normalization’）</p><p>　　該函數(shù)實(shí)

43、現(xiàn)加窗的線性相位FIR數(shù)字濾波器設(shè)計(jì)，可以設(shè)計(jì)標(biāo)準(zhǔn)低通、帶通、高通和帶阻濾波器。</p><p>　　b=fir1（n，Wn）得到n階低通FIR濾波器，調(diào)用后返回維數(shù)為n+1的行向量b，它是濾波器的系數(shù)。b與FIR濾波器的系統(tǒng)函數(shù)的關(guān)系為</p><p><b>　　+1 </b></p><p>　　對(duì)于高通、帶阻濾波器，n取偶數(shù)， 為濾波器

44、的截止頻率，范圍為（0,1）對(duì)于帶通和帶阻濾波器， ，且 ；對(duì)于多通帶濾波器，w= ，頻段為0 </p><p>　　用窗函數(shù)法設(shè)計(jì)線性相位FIR低通濾波器，通帶截止頻率 =0.5 ， ，阻帶衰減不小于40dB，通帶衰減不大于3dB。</p><p>　　下面采用漢寧窗設(shè)計(jì)濾波器：</p><p>　　design FIR filter with windows&l

45、t;/p><p>　　wp=0.5*pi;</p><p>　　ws=0.6*pi;</p><p>　　wdel=ws-wp;</p><p>　　N=ceil(8*pi*wdel);</p><p>　　Wn=(0.5+0.6)*pi/2;</p><p>　　Window=hanning(N+

46、1);</p><p>　　b=firl(N,Wn/pi,window);</p><p>　　freqz(b,1,512)</p><p>　　圖3.1 采用漢寧窗設(shè)計(jì)低通濾波器</p><p>　　用窗函數(shù)法設(shè)計(jì)多通帶濾波器，歸一化通帶為 。</p><p>　　由于高頻段為通帶，濾波器的階數(shù)應(yīng)為偶數(shù)，這里定位36

47、，實(shí)現(xiàn)程序代碼如下，濾波器幅頻-相頻特性如圖</p><p>　　design multiband FIR filter with windows</p><p>　　w=[0.2 0.4 0.6 0.8];</p><p>　　b=fir1(36,w,'dc-1');</p><p>　　freqz(b,1,512);<

48、;/p><p>　　圖3.2 多通帶FIR數(shù)字濾波器</p><p>　　2）fir2。其調(diào)用格式為</p><p>　　b=fir2（n，f，m）</p><p>　　b=fir2（n，f，m，window）</p><p>　　b=fir2（n，f，m，npt）</p><p>　　b=fir2

49、（n，f，m，npf，window）</p><p>　　b=fir2（n，f，m，npf，lap）</p><p>　　b=fir2（n，f，m，npf，lap，window）</p><p>　　該函數(shù)用于設(shè)計(jì)基于窗函數(shù)的任意相應(yīng)的FIR濾波器，其頻率響應(yīng)由向量f和m決定，取值在 之間；n為濾波器階數(shù)；b向量為返回濾波器的系數(shù)；window為窗類型，長(zhǎng)度必須為n

50、+1，默認(rèn)時(shí)為hanning窗；npt為對(duì)頻率響應(yīng)進(jìn)行內(nèi)插點(diǎn)數(shù)，默認(rèn)時(shí)為512；lap參數(shù)用于指定fir2在重復(fù)頻率點(diǎn)附近插入的區(qū)域大小。</p><p>　　設(shè)計(jì)多通帶ＦＩＲ濾波器，濾波器階數(shù)為４０，比較理想濾波器和實(shí)際濾波器的頻率響應(yīng)。</p><p>　　濾波器設(shè)計(jì)程序代碼為</p><p>　　design multiband FIR filter with

51、 fir2</p><p>　　m=[0 0 1 1 0 0 0 1 1 0 0 0 1 1 0 0];</p><p>　　f=[0 0.1 0.15 0.2 0.25 0.3 0.4 0.45 0.5 0.55 0.6 0.7 0.75 0.8 0.85 1];</p><p><b>　　N=40;</b></p><

52、p>　　b=fir2(N,f,m,hamming(N+1));</p><p>　　[h,w]=freqz(b,1,128);</p><p>　　plot(f,m,w/pi,abs(h))</p><p>　　xlabel('frequency')</p><p>　　ylabel('magnitude'

53、;)</p><p><b>　　grid;</b></p><p>　　圖3.3 用fir2設(shè)計(jì)的多通帶FIR數(shù)字濾波器</p><p>　　3.2 FIR濾波器的最優(yōu)設(shè)計(jì)法</p><p>　　firs是fir1和fir2函數(shù)的擴(kuò)展，它采用最小二乘法，使指定頻段內(nèi)的理想分段線性函數(shù)與濾波器幅頻響應(yīng)之間的誤差平方和最

54、小。remez函數(shù)則利用Parks-McClellan算法設(shè)計(jì)線性相位FIR濾波器。該算法能使期望頻率響應(yīng)和實(shí)際頻率響應(yīng)之間的最大誤差最小。firs函數(shù)設(shè)計(jì)的濾波器在頻率響應(yīng)上出現(xiàn)等波紋濾波器。函數(shù)firs、remez調(diào)用的語(yǔ)法規(guī)則相同。</p><p>　　3.2.1 firs函數(shù)</p><p><b>　　其調(diào)用格式是</b></p><p

55、>　　b=firs（n，f，a）</p><p>　　b=firs（n，f，a，w）</p><p>　　b=firs（n，f，a，‘ftype’）</p><p>　　b=firs（n，f，a，w，‘ftype’）</p><p>　　該函數(shù)用于設(shè)計(jì)n階FIR濾波器，其幅頻特性由f和a向量確定，調(diào)用后返回長(zhǎng)度為n+1的濾波器系數(shù)向量b

56、，且這些系數(shù)遵循以下偶對(duì)稱關(guān)系：</p><p>　　f是頻率點(diǎn)向量，范圍為 ，頻率點(diǎn)是逐漸增大，允許向量中有重復(fù)的頻率點(diǎn)；a是指定頻率點(diǎn)的幅頻響應(yīng)，期望的頻率響應(yīng)由（F（k），A（k+1））的連線組成，firls（）則把F（k+1）和F（k+2）之間的頻帶視為過渡帶。所以，所需要的頻率響應(yīng)是分段線性的，其總體平方誤差最小。</p><p>　　b=firls（n，f，a，w）則使用權(quán)系數(shù)

57、w給誤差加權(quán)。W的長(zhǎng)度為F和A的一半。</p><p>　　b=firls（n，f，a，w，‘ftype’）調(diào)用形式中，參數(shù)‘ftype’用于指定所設(shè)計(jì)的濾波器類型，ftype=hilbert，為奇對(duì)稱的線性相位濾波器，返回的濾波器系</p><p>　　數(shù)滿足b（k）=-b（n+2-k），k=1，2.，…，n+1；ftype=differentiator，則采用特殊加權(quán)技術(shù)，生成奇對(duì)稱的

58、線性相位濾波器，使低頻段誤差大大小于高頻段誤差。</p><p>　　設(shè)計(jì)24階反對(duì)陳濾波器，使其具有分段線性的通帶，并畫出理想濾波器和實(shí)際濾波器的頻率響應(yīng)。</p><p><b>　　實(shí)現(xiàn)代碼如下</b></p><p>　　Performance of filter</p><p>　　f=[0 .3 .4 .6

59、.7 .9];</p><p>　　m=[0 1 0 0 0.5 0.5];</p><p>　　b=firls(24,f,m,'hilbert')</p><p>　　for i=1:2:6</p><p>　　plot([f(i) f(i+1)],[m(i),m(i+1)],'--'),hold on<

60、;/p><p><b>　　end</b></p><p>　　[h,f]=freqz(b,1,512,2);</p><p>　　plot(f,abs(h));</p><p><b>　　grid</b></p><p>　　圖3.4 理想和實(shí)際分段帶通濾波器特性</

61、p><p>　　注：虛線是理想濾波器</p><p>　　3.2.2 remez函數(shù)</p><p>　　該函數(shù)的調(diào)用格式與參數(shù)意義都與firls函數(shù)一樣，只是采用的算法不同，下面通過實(shí)力比較它們的區(qū)別</p><p>　　分別使用remez和firls函數(shù)設(shè)計(jì)20階帶通濾波器，通常頻率為Wp1=0.4，Wp2=0.6，實(shí)現(xiàn)代碼如下</p

62、><p>　　Design FIR filter with firls and remez</p><p>　　f=[0 0.3 0.4 0.6 0.7 1];</p><p>　　m=[0 0 1 1 0 0 ];</p><p><b>　　n=20;</b></p><p>　　b=firls(

63、n,f,m);</p><p>　　[h,w]=freqz(b);</p><p>　　bb=remez(n,f,m);</p><p>　　[hh,w]=freqz(bb);</p><p>　　axes('position',[0.2 0.2 0.5 0.5]);</p><p>　　plot(w/

64、pi,abs(h),'b:',w/pi,abs(hh),'r-');</p><p>　　xlabel('frequency');</p><p>　　ylabel('magnitude');</p><p><b>　　grid on;</b></p><p&

65、gt;　　圖3.5 兩種最優(yōu)FIR濾波器的比較</p><p>　　3.3 FIR濾波器最小二乘設(shè)計(jì)法</p><p>　　約束最小二乘法就是在給定濾波器幅頻響應(yīng)最大允許波紋閥值約束條件下，使濾波器在整個(gè)幅頻范圍內(nèi)誤差平方最小化。</p><p>　　Fircls和fircls1兩個(gè)函數(shù)，可以實(shí)現(xiàn)此濾波器的設(shè)計(jì)。</p><p>　　3.3

66、.1 fircls函數(shù)</p><p><b>　　其調(diào)用格式為</b></p><p>　　b=fircls（n，f，a，up，lo）</p><p>　　該函數(shù)返回長(zhǎng)度N+1的線性相位濾波器，期望逼近的頻率分段恒定，由向量f和a確定，頻率的上下限由參數(shù)up和lo確定，長(zhǎng)度與a相同。f中元素為臨界頻率，取值范圍 ，且按遞增順序排列。</

67、p><p>　　使用fircls函數(shù)設(shè)計(jì)一帶通濾波器，通帶為 ，并用設(shè)計(jì)的濾波器對(duì)信號(hào) 濾波（信號(hào)采樣頻率為600Hz）。</p><p><b>　　實(shí)現(xiàn)代碼如下</b></p><p>　　Design FIR filter with fircls</p><p>　　t=0:1/600:1;</p>&l

68、t;p>　　sig=sin(2*pi*15*t)+0.5*sin(2*pi*90*t)+0.2*sin(2*pi*300*t);</p><p>　　plot(t,sig);</p><p><b>　　n=50;</b></p><p>　　f=[0 0.2 0.4 1];</p><p>　　a=[0 1 0]

69、;</p><p>　　up=[0.02 1.02 0.01];</p><p>　　lo=[-0.02 0.98 -0.01];</p><p>　　b=fircls(n,f,a,up,lo);</p><p>　　newsig=fftfilt(b,sig)；</p><p><b>　　figure<

70、;/b></p><p>　　ft=t(301:350)；</p><p>　　ns=newsig(301:350);</p><p>　　zns=interp(ns,12);</p><p>　　znt=interp(ft,12);</p><p>　　plot(znt,zns);</p><

71、;p><b>　　時(shí)間/s</b></p><p><b>　　時(shí)間/s</b></p><p>　　圖3.6 濾波前后的波形</p><p>　　3.3.2 fircls函數(shù)</p><p>　　該函數(shù)采用約束最小二乘法設(shè)計(jì)基本的線性相位高通和低通濾波器，其調(diào)用格式為</p>

72、;<p>　　b=fircls（n，wo，dp，ds）</p><p>　　b=fircls1（n,wo,dp,ds,‘high’）</p><p>　　b=fircls1(n,wo,dp,ds,wt)</p><p>　　b=fircls1(n,wo,dp,ds,wt,‘high’)</p><p>　　b=fircls1(n

73、,wo,dp,ds,wp,ws,k)</p><p>　　b=fircls1(n,wo,dp,ds,wp,ws,k,‘high’)</p><p>　　b=fircls1(n,wo,dp,ds,…‘design-flag’)</p><p>　　b=fircls（n，wo，dp，ds）函數(shù)返回長(zhǎng)度為n+1的線性相位低通FIR濾波器，截止頻率為wo，取值0～1之間。通

74、常幅度偏離1.0的最大值為dp，阻帶偏離0的最大值為ds。</p><p>　　b=fircls1（n,wo,dp,ds,‘high’）則返回高通濾波器，n必為偶數(shù)。</p><p>　　b=fircls1(n,wo,dp,ds,wp,ws,k)采用平方誤差加權(quán)，通常的權(quán)值比阻帶的大k倍。Wp為通帶邊緣頻率，ws為阻帶邊緣頻率，其中 ；如果要設(shè)計(jì)高通濾波器，則必須使 。</p>

75、<p>　　設(shè)計(jì)55階低通濾波器，截止頻率在0.3。</p><p><b>　　設(shè)計(jì)代碼如下</b></p><p>　　Design FIR filter with fircls1</p><p><b>　　n=55;</b></p><p><b>　　wo=0.3;

76、</b></p><p>　　dp=0.02;ds=0.008;</p><p>　　b=fircls1(n,wo,dp,ds);</p><p>　　[H,W]=freqz(b);</p><p>　　plot(W/pi,abs(H));grid</p><p>　　圖3.7 fircls1 設(shè)計(jì)的低通

77、濾波器</p><p>　　3.4 FIR濾波器升余弦函數(shù)設(shè)計(jì)法</p><p>　　用fircos函數(shù)，用于設(shè)計(jì)有光滑、升余弦過渡帶的低通線性相位濾波器，其調(diào)用格式為：</p><p>　　b=fircos（n，F(xiàn)o，df，fs）</p><p>　　b=fircos（n，F(xiàn)o，df，fs，‘bandwidth’）</p>

78、<p>　　b=fircos（n，F(xiàn)o，df）</p><p>　　b=fircos（n，F(xiàn)o，r，fs，‘rolloff’）</p><p>　　b=fircos（…,‘type’）</p><p>　　b=fircos（…，‘type’，delay，window）</p><p>　　=fircos（…）</p>

79、<p>　　其中，n為濾波器階數(shù)；Fo為低通濾波器的截止頻率，df為過渡帶頻寬，fs是采樣頻率，單位都是Hz；Fo和Fo+ 都介于0～ 之間；b為返回濾波器系數(shù)向量。</p><p>　　設(shè)計(jì)20階升余弦濾波器，截止頻率為250Hz，過渡帶寬為100Hz，采樣頻率為1000Hz。</p><p><b>　　其程序代碼如下</b></p>&l

80、t;p>　　Design FIR filter with firrcos</p><p><b>　　n=20;</b></p><p><b>　　Fo=250;</b></p><p><b>　　df=100;</b></p><p><b>　　Fs=1

81、000;</b></p><p>　　b=firrcos(n,Fo,df,Fs);</p><p>　　[h,f]=freqz(b,1,512,Fs);</p><p>　　axes('position',[0.2 0.2 0.5 0.5]);</p><p>　　plot(f,abs(h));</p>

82、<p>　　xlabel('frequency(Hz)');</p><p>　　ylabel('magnitude');</p><p><b>　　grid on；</b></p><p>　　圖3.8 升余弦函數(shù)設(shè)計(jì)低通濾波器</p><p>　　3.5 FIR濾波器

83、的任意響應(yīng)設(shè)計(jì)法</p><p>　　Cremez函數(shù)是設(shè)計(jì)FIR濾波器的另一函數(shù)，該函數(shù)可用來設(shè)計(jì)任意響應(yīng)的復(fù)數(shù)和非線性相位等波紋FIR濾波器。它在指定的濾波器頻率響應(yīng)的方式上和其他濾波器設(shè)計(jì)函數(shù)有所不同，接收的是函數(shù)。</p><p>　　Cremez函數(shù)采用擴(kuò)展remez交換算法，能使契比雪夫誤差最小，其調(diào)用格式為</p><p>　　b=cremez（n，f

84、，‘fresp’）</p><p>　　b=cremez（n，f，‘fresp’，w）</p><p>　　b=cremez（n，f，{‘fresp’，p1，p2，…}，w）</p><p>　　b=cremez（n，f，a，w）</p><p>　　b=cremez（…，‘sym’）</p><p>　　b=crem

85、ez（…，‘skip-stage2’）</p><p>　　b=cremez（…，‘debug’）</p><p>　　b=cremez（…，{lg rid}）</p><p>　　[b，delta，opt]=crème（…）</p><p>　　其中，參數(shù)f是頻率邊界向量，范圍是-1到1。W是每個(gè)頻帶的加權(quán)系數(shù)，長(zhǎng)度為f的一半。p

86、1/p2是可選參數(shù)，fresp則指定濾波器類型，可能取值為：</p><p>　　lowpass、highpass、bandpass、bandstop。比如以下引用形式</p><p>　　b=crème（n，f，‘lowpass’，…）</p><p>　　b=cremez（n，f，{‘lowpass’，d}，…）</p><p>

87、;　　multiband。該函數(shù)指定設(shè)計(jì)多通帶濾波器，例如：</p><p>　　b=cremez（n，f，{‘multiband’，a}，…）</p><p>　　b=cremez（n，f，{‘multiband’，a，d}，…）</p><p>　　其中，a為在邊界頻帶向量f各點(diǎn)上的期望幅值。</p><p>　　differentiat

88、or。該函數(shù)指定設(shè)計(jì)線性相位微分器，其零頻必須在過渡帶，而頻帶加權(quán)設(shè)置與頻率成反比，引用方式如下，其中fs用于確定微分器的響應(yīng)斜率，默認(rèn)值為1。</p><p>　　b=cremez（n，f，{‘differentiator’，fs}，…）</p><p>　　b=cremez（n，f，{‘differentiator’，fs，d}，…）</p><p>　　hil

89、bfilt。該函數(shù)指定設(shè)計(jì)線性相位希爾伯特變換濾波器，其零頻必須在過渡帶，比如：</p><p>　　b=crème（n，f，‘hilbfilt’，…）</p><p>　　b=cremez（N，F(xiàn)，{‘hilbfilt’，d}，…）</p><p>　　b=cremez（…，‘sym’）則使用所設(shè)計(jì)濾波器的沖擊響應(yīng)有對(duì)稱性，‘sym’取值如下：</

90、p><p>　　‘none’，默認(rèn)值，若f中有負(fù)值或‘fresp’未賦值則為默認(rèn)值。</p><p>　　‘even’，使單位沖擊響應(yīng)為實(shí)偶對(duì)稱序列。</p><p>　　‘odd’，使單位沖擊響應(yīng)為實(shí)奇對(duì)稱序列。設(shè)計(jì)希爾伯特和差分濾波器時(shí)，取默認(rèn)值。</p><p>　　‘real’，使頻率響應(yīng)為共軛對(duì)稱。</p><p&g

91、t;　　b=crème（…，‘debug’）則將顯示中間計(jì)算結(jié)果，debug可為‘trace’、‘plot’、‘both’或者‘off’之一，默認(rèn)值為‘off’。</p><p>　　設(shè)計(jì)30階任意響應(yīng)多頻帶濾波器。濾波器的頻率邊界頻率向量為</p><p>　　f=[-1，-0.5，-0.4，0.3，0.4，0.8]；各頻段的幅值響應(yīng)為a=[3，1，2，2，2，1]；各頻段最優(yōu)

92、化權(quán)向量為w=[1，10，5]</p><p><b>　　設(shè)計(jì)代碼如下</b></p><p>　　Complex and nonlinear phase FIR filter design</p><p><b>　　n=30;</b></p><p>　　f=[-1 -0.5 -0.4 0.3

93、 0.4 0.8];</p><p>　　a=[3 1 2 2 2 1];</p><p>　　w=[1 10 5];</p><p>　　b=cremez(n,f,{'multiband',a},w);</p><p>　　[h,wf]=freqz(b,1,512,'whole');</p>&

94、lt;p>　　plot(wf/pi-1,fftshift(abs(h)));</p><p>　　xlabel('frequency');</p><p>　　ylabel('maginitude')</p><p><b>　　grid on;</b></p><p>　　圖3.9

95、 任意復(fù)響應(yīng)濾波器頻譜圖</p><p>　　第四章　FIR 濾波器階數(shù)估計(jì)</p><p>　　4.1 kaiserord函數(shù)</p><p><b>　　其調(diào)用格式為</b></p><p>　　[n,Wn,beta,ftype]=kaiserord(f,a,dev)</p><p>　　

96、[n,Wn,beta,ftype]=kaiserord(f,a,dev,fs)</p><p>　　C=kaiserord(f,a,dev,fs,’cell’)</p><p>　　該函數(shù)返回值在函數(shù)b=fir1（n，Wn，kaiser（n+1，beta），ftype，‘noscale’）中需要使用的階數(shù)n、歸一化頻帶邊緣Wn、kaiser窗函數(shù)beta、濾波器類型ftype的估計(jì)值。所設(shè)

97、計(jì)的濾波器近似滿足由輸入?yún)?shù)f、a和dev給定的性能。F為頻率向量，a為f定義的各頻帶的理想頻率響應(yīng)幅值，f的長(zhǎng)度是a的2倍減去2（是偶數(shù)），第一頻帶的起始頻率為0，最后一個(gè)頻帶截止頻率為fs。向量dev的各元素分別是各頻帶允許的最大偏移量。</p><p>　　而c=kaiserord（f，a，dev，fs，‘cell’）返回的是向量組，用于fir1函數(shù)的參數(shù)。</p><p>　　設(shè)計(jì)

98、低通濾波器，通帶范圍為0～1Hz，阻帶范圍為4～1500Hz，通帶波紋為5％，阻帶衰減為40dB，系數(shù)采樣頻率fs=8000Hz。</p><p><b>　　設(shè)計(jì)程序如下</b></p><p>　　Design lowpass filter with kaiserord</p><p>　　fsamp=8000;</p>&l

99、t;p>　　fcuts=[1000 1500];</p><p>　　mags=[1 0];</p><p>　　devs=[0.05 0.01];</p><p>　　[n,wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsamp);</p><p>　　hh=fir1(n,wn,ftype,kai

100、ser(n+1,beta),'noscale');</p><p><b>　　freqz(hh)</b></p><p>　　圖4.1 利用kaiserord設(shè)計(jì)的低通濾波器幅頻特性</p><p>　　4.2 remezord函數(shù)</p><p>　　remezord函數(shù)為remez函數(shù)選擇濾波

101、器階數(shù)，在給定頻域中的性能指標(biāo)后，remezord可產(chǎn)生近似的滿足指標(biāo)的最小階數(shù)。函數(shù)調(diào)用格式如下：</p><p>　　[n，fo，ao，w]=remezord（f，a，dev）</p><p>　　[n，fo，ao，w]=remezord（f，a，dev，fs）</p><p>　　c=remezord（f，a，dev，fs，‘cell’）</p>

102、<p>　　該函數(shù)能找出近似階數(shù)n、歸一化頻率邊界fo、頻帶內(nèi)幅值ao、以及權(quán)向量w，使由remez函數(shù)構(gòu)成的濾波器滿足參數(shù)f、a、dev指定的性能要求。f為頻率向量，長(zhǎng)度是參數(shù)a的2倍減去2（是偶數(shù)），而a為f各頻帶理想頻率響應(yīng)的幅值，向量dev的元素分別是各頻帶允許的最大偏差。</p><p>　　設(shè)計(jì)最小階數(shù)的低通濾波器，通常截止頻率為500Hz，阻帶截止頻率為600Hz采樣頻率為2000Hz，

103、阻帶最小衰減為40dB，通帶最大衰減為3dB。</p><p><b>　　程序代碼如下</b></p><p>　　Design lowpass filter with kaiserord</p><p><b>　　rp=3;</b></p><p><b>　　rs=4.0;<

104、/b></p><p><b>　　fs=2000;</b></p><p>　　f=[500 600];</p><p><b>　　a=[1 0];</b></p><p>　　dev=[(10^(rp/20)-1)/(10^(rp/20)+1) 10^(-rs/20)];</p&g

105、t;<p>　　[n,fo,ao,w]=remezord(f,a,dev,fs);</p><p>　　b=remez(n,fo,ao,w);</p><p>　　freqz(b,1,1024,fs);</p><p>　　title('lowpass filter designed to specifications');</p

106、><p>　　圖4.2 利用remezord設(shè)計(jì)的低通濾波器幅頻特性</p><p>　　估計(jì)函數(shù)還包括buttord、cheblord、cheb2ord、ellipord、kaiserord、firpm，它們分別用于各自不同類型的濾波器的階數(shù)估計(jì)。</p><p><b>　　第五章 小結(jié)</b></p><p>　　

107、濾波是信號(hào)處理中最基本有極為重要的技術(shù)，利用濾波器技術(shù)可以從復(fù)雜的信號(hào)中提取出所需要的信號(hào)，抑制不需要的信號(hào)。</p><p><b>　　FIR濾波器的優(yōu)點(diǎn)</b></p><p>　　FIR濾波器能嚴(yán)格做到線性相位或群延時(shí)為常數(shù)，而IIR濾波器只能逼近線性相位。</p><p>　　FIR濾波器是全零點(diǎn)型濾波器，總是穩(wěn)定的，不會(huì)因?yàn)V波運(yùn)算的

108、舍入誤差而產(chǎn)生極限環(huán)振蕩現(xiàn)象</p><p>　　對(duì)同樣幅度相應(yīng)的濾波器，用FIR濾波器實(shí)現(xiàn)比用IIR濾波器實(shí)現(xiàn)需要較高的階數(shù)，多達(dá)5~10倍。當(dāng)濾波器的特性要求較高時(shí)，用FIR濾波器來實(shí)現(xiàn)，濾波過程需要較多的計(jì)算時(shí)間</p><p><b>　　MATLAB軟件</b></p><p>　　作為強(qiáng)大的計(jì)算軟件， MATLAB在濾波器設(shè)計(jì)方面還

109、提供了編寫圖形用戶界面的功能。所謂圖形用戶界面, 簡(jiǎn)稱為GUI, 是由各種圖形對(duì)象, 如圖形窗口菜單按鈕、文本框等構(gòu)建的用戶界面，它是用戶和計(jì)算機(jī)之間交流的工具。MATLAB 將所有GUl 支持的用戶控件都集成在這個(gè)環(huán)境中并提供界面外觀、屬性和行為響應(yīng)方式的設(shè)置方法,隨著版本的提高,這種能力還會(huì)不斷加強(qiáng)。而且具有強(qiáng)大的繪圖功能,可以輕松的獲得更高質(zhì)量的曲線圖。MATLAB語(yǔ)言進(jìn)行程序設(shè)計(jì)和利用信號(hào)處理工具箱的FDATool工具進(jìn)行界面設(shè)

110、計(jì)可以隨時(shí)對(duì)比設(shè)計(jì)要求和濾波器特性調(diào)整參數(shù)，直觀簡(jiǎn)便，極大的減輕了工作量，有利于濾波器設(shè)計(jì)的最優(yōu)化。還可在MTLAB環(huán)境下的仿真軟件Simulink對(duì)所設(shè)計(jì)的濾波器進(jìn)行模擬仿真。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　陳懷琛等.MATLAB及在電子信號(hào)課程中的應(yīng)用.北京：電子工業(yè)出版社，2002</p><p>　　

111、張賢達(dá)著.現(xiàn)代信號(hào)處理（第二版）.北京：清華大學(xué)出版社，2002</p><p>　　A.V.奧本海姆，R.W.謝弗.數(shù)字信號(hào)處理.北京：科學(xué)出版社，1982</p><p>　　張志涌等.精通MATLAB6.5版.北京：北京航空航天大學(xué)出版社，2003</p><p>　　羅俊輝等.MATLAB7.0在數(shù)字信號(hào)處理中的應(yīng)用.機(jī)械工業(yè)出版社，2005</p&g

112、t;<p>　　顧啟泰.離散事件系統(tǒng)建模與仿真.北京：清華大學(xué)出版社，1999</p><p>　　李人厚等譯校.精通MATLAB綜合輔導(dǎo)與指南.西安：西安交通大學(xué)出版社，1997</p><p>　　林茂庸，柯有安.雷達(dá)信號(hào)理論[M].北京：國(guó)防工業(yè)出版社，1981</p><p>　　張直中.雷達(dá)信號(hào)的選擇與處理[M].北京：國(guó)防工業(yè)出版社，197

113、9</p><p>　　[美]R.L.米切爾.雷達(dá)系統(tǒng)模擬，1982</p><p>　　何友，關(guān)鍵，彭應(yīng)寧等.雷達(dá)自動(dòng)檢測(cè)與恒虛警處理.北京：清華大學(xué)出版社，1999</p><p>　　MATLAB User’s Guide，The Mathworks，inc.2004</p><p>　　Signal Processing Toolbo

114、x User’Guide，The Mathworks，inc.2004</p><p><b>　　致 謝</b></p><p>　　這次畢業(yè)設(shè)計(jì)，是在劉老師的悉心指導(dǎo)下完成的，他以其嚴(yán)謹(jǐn)?shù)膶W(xué)術(shù)態(tài)度、廣博的專業(yè)知識(shí)、獨(dú)到的專業(yè)見解，讓我學(xué)到了好多專業(yè)和非專業(yè)知識(shí)和技能，同時(shí)，鍛煉了自己的能力，為以后的工作和學(xué)習(xí)打下了良好的基礎(chǔ)，在此我要衷心感謝我的指導(dǎo)老師劉老師，向