畢業(yè)設計論文基于MATLAB的模擬濾波器的設計

1、基于matlab的模擬濾波器的設計摘 要幾乎在所有的工程技術領域中都會涉及到信號處理問題 ，濾波器作為信號處理的重要組成部分，已發(fā)展的相當成熟。本論文首先介紹了濾波器的濾波原理以及模擬濾波器、數(shù)字濾波器的設計方法。重點介紹了模擬濾波器的設計和仿真。系統(tǒng)研究了模擬濾波器（包括巴特沃斯濾波器和切比雪夫濾波器）的設計原理和方法，并在此基礎上論述了模擬濾波器（包括低通、高通、帶通、帶阻）的設計。在此基礎上，用matlab虛擬實現(xiàn)模擬濾波器。此設計擴展性好，便于調(diào)節(jié)濾波器的性能，可以根據(jù)不同的要求在matlab上加以實現(xiàn)。通過matlab的仿真與實現(xiàn)，可以看出傳統(tǒng)的模擬濾波器設計方法繁瑣且不直觀，而ma

2、tlab具有較嚴謹?shù)目茖W計算和圖形顯示這一優(yōu)點，使設計結果顯示的更加直觀，而且對濾波器的精度也有了很大的提高，能更好的達到預期效果。同時，又對模擬濾波器低通至高通、帶通、帶阻的轉換進行了理論上的闡述。 關鍵字： 濾波/模擬濾波器/matlab matlab based design of analog filtersabstractin almost all areas of engineering and technology will be related to signal processing, signal processing filter as an important comp

3、onent of the development has reached a mature state. this paper introduces the principle of filter and filter analog filters, digital filter design method. focuses on the design of analog filters and simulation.analog filter system (including the butterworth filter and chebyshev filter) design princ

4、iples and methods, and on this basis on the analog filters (including low-pass, high pass, band-pass, band-stop) design. on this basis, the use of virtual realization of analog filters matlab. good scalability of this design, easy to adjust the performance of filters can be based on different requir

5、ements to be in the realization of matlab.through the matlab simulation and realization of, we can see that the traditional design method of analog filters and do not intuitive red, and matlab with more stringent scientific computing and graphical display of the advantages of the design showed that

6、the more intuitive, but also to filter accuracy has been greatly improved to better achieve the desired results. at the same time, low-pass analog filters of high pass, band-pass, band-stop conversion of a theoreticalkeyword: filtering, analog filters, matlab目 錄中文摘要i英文摘要ii1 濾波器簡介11.1 濾波器的工作原理11.1.1

7、模擬濾波器的工作原理11.1.2 數(shù)字濾波器的工作原理21.2 濾波器的基本特性31.2.1 模擬濾波器與數(shù)字濾波器的基本特性31.2.2 無限沖擊響應iir和有限沖擊響應fir濾波器51.3 濾波器的主要技術指標52 模擬濾波器的設計62.1 模擬濾波器的設計方法72.2 模擬原型濾波器及最小階數(shù)的選擇82.2.1 巴特沃斯濾波器及最小階數(shù)的選擇82.2.2 切比雪夫濾波器及最小階數(shù)的選擇122.2.3 橢圓濾波器及最小階數(shù)的選擇172.2.4貝塞爾濾波器183 matlab仿真183.1 matlab簡介183.2 對低通模擬濾波器的仿真203.3 模擬高通濾波器的仿真213.4 模擬帶通

8、濾波器的仿真223.5 對帶阻模擬濾波器的仿真244 頻率轉換254.1低通至高通的轉換254.2低通至帶通的變換264.3低通至帶阻的變換285 總結與展望30致 謝31參考文獻321 濾波器簡介從廣義上講，任何對某些頻率（相對于其他頻率來說）進行修正的系統(tǒng)稱為濾波器。嚴格地講，對輸入信號通過一定的處理得到輸出信號，這個處理通常是提取信號中某頻率范圍內(nèi)的信號成分，把這種處理的過程稱為濾波。實現(xiàn)濾波處理的運算電路或設備稱為濾波器。在許多科學技術領域中，廣泛應用線性濾波和頻譜分析對信號進行加工處理，模擬濾波是處理連續(xù)信號，數(shù)字濾波則是處理離散信號，而后者是在前者的基礎上發(fā)展起來的。我們知道，無源

9、或有源模擬濾波器是分立元件構成的線性網(wǎng)絡，他們的性能可以用線性微分方程來描述，而數(shù)字濾波器是個離散線性系統(tǒng)，要用差分方程來描述，并以離散變換方法來分析。這些方程組可以用專用的或通用的數(shù)字計算機進行數(shù)字運算來實現(xiàn)。因此，數(shù)字濾波器的濾波過程是一個計算過程，它將輸入信號的序列數(shù)字按照預定的要求轉換成輸出數(shù)列。1.1 濾波器的工作原理1.1.1 模擬濾波器的工作原理我們知道，模擬濾波器是對模擬信號實行線性濾波的一種線性時不變系統(tǒng)，如圖1.1所示。在時域內(nèi)，它的動態(tài)特性可以用系統(tǒng)的單位沖激函數(shù)的響應來描述，也就是該濾波系統(tǒng)在任何時刻對輸入單位沖激信號=（t）的輸出響應。這個函數(shù)從時域上反映了該濾波系統(tǒng)

10、的傳輸特性。對于任意輸入信號，系統(tǒng)的輸出可以卷積表示： = (1.1)上式表明在對線性濾波器系統(tǒng)進行時域分析時，采用了疊加原理，先將任意輸入信號波形分成不同時間的窄脈沖之和，再分別求出各個脈沖通過濾波器之后的響應，并進行線性疊加從而得到總的輸出信號。模擬濾波器,h(s)圖1.1模擬濾波器原理在頻域分析時，線性濾波器的轉移函數(shù)等于系統(tǒng)的單位沖激函數(shù)的響應的拉普拉斯變換： （1.2）很明顯，當s=j，上式就是傅立葉變換的表達式，它反映了濾波器的傳輸特性對各種頻率的響應，也就是濾波器的頻率響應函數(shù)，它決定著濾波特性。當濾波器輸入信號與輸出信號的拉普拉斯變換，得 （1.3）這表明兩信號卷積的變換等于各

11、自變換的乘積。在頻譜關系上，一個輸入信號的頻譜，經(jīng)過濾波器的作用后，被變換成的頻譜。因此，根據(jù)不同的濾波要求來選定，就可以得到不同類型的模擬濾波器。還可以看出，濾波器的濾波過程就是完成信號與它的單位沖激函數(shù)響應之間的數(shù)學卷積運算過程。1.1.2 數(shù)字濾波器的工作原理在數(shù)字濾波中，我們主要討論離散時間序列。如圖1.2所示。設輸入序列為，離散或數(shù)字濾波器對單位抽樣序列的響應為。因在時域離散信號和系統(tǒng)中所起的作用相當于單位沖激函數(shù)在時域連續(xù)信號和系統(tǒng)中所起的作用。數(shù)字濾波器,h(z)圖1.2 數(shù)字濾波器原理數(shù)字濾波器的序列將是這兩個序列的離散卷積，即 （1.4） 同樣，兩個序列卷積的z變換等于個自z

12、變換的乘積，即 （1.5）用代入上式，其中t為抽樣周期，則得到 （1.6）式中和 分別為數(shù)字濾波器輸入序列和輸出序列的頻譜，而為單位抽樣序列響應的頻譜。由此可見，輸入序列的頻譜經(jīng)過濾波后，變?yōu)?，按照的特點和我們處理信號的目的，選取適當?shù)氖沟臑V波后的符合我們的要求。1.2 濾波器的基本特性1.2.1 模擬濾波器與數(shù)字濾波器的基本特性如利用模擬電路直接對模擬信號進行處理則構成模擬濾波器，它是一個連續(xù)時間系統(tǒng)。如果利用離散時間系統(tǒng)對數(shù)字信號（時間離散、幅度量化的信號）進行濾波則構成數(shù)字濾波器。數(shù)字濾波器的差分方程表示為： 系統(tǒng)函數(shù)表示:數(shù)字濾波器的特性通常用其頻率響應函數(shù)來描述，包括幅度特性和相位

13、特性。按信號通過系統(tǒng)時的特性（主要是幅頻特性）來分類：可以有低通、高通、帶通和帶阻四種基本類型。（1） 低通數(shù)字濾波器：圖1.3所示|h(ej)|c-c-22-fs/2-fsfs/2fs-fcfcf圖1.3 低通數(shù)字濾波器的頻譜（2） 高通數(shù)字濾波器：圖1.4所示|h(ej)|c-c-22圖1.4 高通數(shù)字濾波器的頻譜（3） 帶通數(shù)字濾波器：圖1.5所示|h(ej)|2-2-221-1圖1.5 帶通數(shù)字濾波器的頻譜（4）帶阻數(shù)字濾波器：圖1.6所示1|h(ej)|2-2-22-1 圖1.6 帶阻數(shù)字濾波器的頻譜 其他較復雜的特性可以由基本濾波器組合。1.2.2 無限沖擊響應iir和有限沖擊響應

14、fir濾波器按系統(tǒng)沖擊響應（或差分方程）可以分成無限沖擊響應 iir和有限沖擊響應fir濾波器兩類。這兩種濾波器都可以現(xiàn)實各種頻率特性要求，但它們在計算流程、具體特性逼近等方面是有差別的。(1) fir濾波器(非遞歸型):(2) iir濾波器（遞歸型） 還有一些其他的分類方法，例如在特定場合使用的濾波器。1.3 濾波器的主要技術指標濾波器的主要技術指標取決于具體的應用或相互間的相互關系。具體的有最大通帶增益（即通帶允許起伏）；最大阻帶增益；通帶截止頻率；阻帶截止頻率。如圖1.7所示db0dbdb圖1.7 濾波器的主要技術指標2 模擬濾波器的設計模擬濾波器的理論和設計方法已經(jīng)發(fā)展的相當成熟，且有

15、若干典型的模擬濾波器供我們選擇，如巴特沃斯（butterworth濾波器.切比雪夫（chebyshev）濾波器等。這些工作的理論分析和設計方法在20世紀30年代就完成，然而煩瑣.冗長的數(shù)字計算使它難以付諸實用。直到50年代，由于計算機技術的逐步成熟，求出大量設計參數(shù)和圖表，這種方法才得到廣泛應用。這些典型的濾波器各有特點：巴特沃斯濾波器具有單調(diào)下降的幅頻特性；切比雪夫濾波器的幅頻特性在通帶或者阻帶有波動發(fā)，可以提高選擇性。這樣根據(jù)具體要求可以選擇不同類型的濾波器。模擬濾波器按幅度特征可以分成低通、高通、帶通和帶阻濾波器。它們的理想幅度特性如圖2.1所示，但我們設計濾波器時，總是先設計低通濾波器

16、，再通過頻率變換將低通濾波器轉換成希望類型的濾波器帶通帶阻圖2.1 模擬濾波器理想幅度特性低通高通2.1 模擬濾波器的設計方法利用頻率變換設計模擬濾波器的步驟為：（1）給定模擬濾波器的性能指標，如截止頻率或上、下邊界頻率等。（2）確定濾波器階數(shù)（3）設計模擬低通原型濾波器。（4）按頻率變換設計模擬濾波器（低通、高通、帶通、帶阻）。模擬低通濾波器的設計指標有，和，其中和分別稱為通帶截止頻率和阻帶截止頻率。 是通帶（=0）中的最大衰減系數(shù)，是阻帶的最小衰減系數(shù)，和一般用db表示。對于單調(diào)下降的幅度特性，可表示成： （2.1） （2.2）如果=0處幅度已歸一化為一，即，和表示為 （2.3） （2.4

17、）以上技術指標用圖2.2表示，圖中稱為3db 截止頻率，因,-20 10.7070圖2.2 低通濾波器的幅度特性濾波器的技術指標給定以后，需要設計一個傳輸函數(shù)，希望其幅度平方函數(shù)滿足給定的指標和，一般濾波器的單位沖激響應為實數(shù)，因此 = (2.5)如果能由，求出，那么就可以求出所需的，對于上面介紹的典型濾波器，其幅度平方函數(shù)有自己的表達式，可以直接引用。這里要說明的是必須是穩(wěn)定的。因此極點必須落在s平面的左半平面，相應的的極點落在右半平面。2.2 模擬原型濾波器及最小階數(shù)的選擇2.2.1 巴特沃斯濾波器及最小階數(shù)的選擇巴特沃斯濾波器是最基本的逼近方法形式之一。它的幅頻特性模平方為 (2.6)式

18、中n是濾波器的階數(shù)。當=0時，；當=時，是3db截止頻率。不同階數(shù)n的巴特沃斯濾波器特性如圖2.3所示，這一幅頻特性具有下列特點：（1）最大平坦性：可以證明：在=0點，它的前（2n-1）階導數(shù)都等于0，這表明巴特沃斯濾波器在=0附近一段范圍內(nèi)是非常平直的，它以原點的最大平坦性來逼近理想低通濾波器?！白钇巾憫奔从纱硕鴣?。（2）通帶，阻帶下降的單調(diào)性。這種濾波器具有良好的相頻特性。（3）3db的不變性：隨著n的增加，頻帶邊緣下降越陡峭，越接近理想特性，但不管n是多少，幅頻特性都通過-3db點。當時，特性以20ndb/dec速度下降。圖2.3 不同階數(shù)n的巴特沃斯濾波器特性現(xiàn)根據(jù)式（2.6）求巴特

19、沃斯濾波器的系統(tǒng)函數(shù)ha（s）。令=s/j，帶入式（2.6）對應的極點： (2.7)即為的極點，此極點分布有下列特點：(1)的2n個極點以/n為間隔均勻分布在半徑為的圓周上，這個圓稱為巴特沃斯圓。(2)所有極點以j軸為對稱軸成對稱分布，j軸上沒有極點。(3)當n為奇數(shù)時，有兩個極點分布在的實軸上；n為偶函數(shù)時，實軸上沒有極點。所有復數(shù)極點兩兩呈共軛對稱分布。圖2.4畫出了n=3時的極點分布。全部零點位于s=處。j圖2.4 n=3時ha（s）ha（-s）極點分布為得到穩(wěn)定的，取全部左半平面的極點。 （2.8）當n為偶數(shù)時 （2.9）當n為奇數(shù)時 （2.10）為使用方便把式（2.9）和式（2.10

20、）對進行歸一化處理，為此，分子分母各除以，并令，稱為歸一化復頻率：（n為偶數(shù)） （2.11） （n為奇數(shù)）（2.12）用歸一化頻率表示的頻率特性稱為原型濾波特性（即歸一化復頻率s 的虛部）。對式（2.6）所示的低通巴特沃斯特性用表示得到： （2.13）稱為巴特沃斯低通原型濾波器幅頻特性。在低通原型濾波頻率特性上，截止頻率=1。若給出模擬低通濾波器的設計性能指標要求：通帶邊界頻率，阻帶邊界頻率，通帶波紋，阻帶衰減，要確定butterworth ,，低通濾波器最小階數(shù)n及截止頻率。，的意義如圖所示。當=時， 即，以截至頻率（幅值下降3db）為1，化為相對為相對的相對頻率由上式可寫為。同理，當=時，

21、 。由此可見 n應向上取整，再用matlab 編程計算濾波器最小階數(shù)n和截止頻率。2.2.2 切比雪夫濾波器及最小階數(shù)的選擇巴特沃斯濾波器的頻率特性曲線，無論在通帶和阻帶都是頻率的單調(diào)函數(shù)。當通帶邊界處滿足指標要求時，通帶內(nèi)肯定會有余量。因此，更有效的設計方法應該是將精確度均勻地分布在整個通帶內(nèi)。這可通過選擇具有等波紋特性的逼近函數(shù)來達到。切比雪夫濾波器的振幅特性就是具有這種等波紋特性。它有兩種型式：振幅特性在通帶內(nèi)是等波紋的，在阻帶內(nèi)是單調(diào)的切比雪夫i型濾波器；振幅特性在通帶內(nèi)是單調(diào)的，在阻帶內(nèi)是等波紋的切比雪夫ii型濾波器。采用何種型式切比雪夫濾波器取決于實際用途。這種濾波器的幅頻特性模平

22、方為： （2.14）式中是決定通帶內(nèi)起伏的等波紋參數(shù)，是第一類切比雪夫多項式，定義為：= （2.15）表2.1列出了對應不同階數(shù)n時的切比雪夫多項式。圖（2.5）畫出了多項式特性曲線，從這組特性曲線可以看出：x1時，在±1之間波動；n不論為何值都有=1；當x>1，單調(diào)上升。此外，切比雪夫多項式滿足下列遞推公式 n=1，2 （2.16）圖2.6（a）是按式（2.14）畫出的切比雪夫等波紋濾波器的幅頻特性，圖2.6（b）是通帶內(nèi)起伏與的關系。切比雪夫濾波器的濾波特性具有下列特點：（1） 所有曲線在=時通過點，因而把定義為切比雪夫濾波器的截止角頻率。（2） 在通帶內(nèi)/1，在1和之間變

23、化；在通帶外，/>1，特性呈單調(diào)下降，下降速度為20ndb/dec。（3） n為奇數(shù)，=1；n為偶數(shù)，=。通帶內(nèi)誤差分布是均勻的，實際上這種逼近稱為最佳一致逼近。（4） 由于濾波器通帶內(nèi)有起伏，因而使通帶內(nèi)的相頻特性也有相應的起伏波動。即相位是非線性的，這給信號傳輸時帶來線性畸變，所以在要求群時延為常數(shù)時不宜采用這種濾波器?，F(xiàn)根據(jù)式（2.14）求切比雪夫濾波器的系統(tǒng)函數(shù)。將=帶入式（2.14） （2.17）為求極點分布需求解方程： （2.18）表2.1 n=07時切比雪夫多項式tn（x）ntn(x)ntn(x)01231x456711-1-1xt1(x)11-1-1xt2(x)11-1-

24、1xt3(x)11-1-1xt4(x)圖2.5 t1t4切比雪夫特性曲線考慮到是復變量，為解出切比雪夫多項式，設：= （2.19）1n=4n=5(a)11o-111-1o0.590.310.950.75(b)圖2.6 切比雪夫濾波特性及內(nèi)波紋關系另把=cos代入式（2.15），并且令此式等于，求解，： （2. 20）解的滿足上式的，為 （2.21）把，值代回式（2.19），求的極點值：+ ，k=1,2,，2n (2.22)就是切比雪夫濾波器的極點，給定n，即可求的2n個極點分布。由式（2.22）實部與虛部的正弦和余弦函數(shù)平方約束關系可以看出，此極點分布滿足橢圓方程，其短軸和長軸分別為 （2.2

25、3）圖2.7畫出了n=3時切比雪夫濾波器的極點分布。abj 圖2.7極點所在的橢圓可以和半徑為a的圓和半徑為b的圓聯(lián)系起來，這兩個圓分別稱為巴特沃斯小圓和巴特沃斯大圓。n階切比雪夫濾波器極點的縱坐標，而橫坐標等于n階巴特沃斯小圓極點的橫坐標取左半平面的極點： k=1,2，n （2.24）則切比雪夫濾波器的系統(tǒng)函數(shù)： （2.25）其中，常數(shù)a=。因而切比雪夫濾波器的系統(tǒng)函數(shù)表示為： （2.27）切比雪夫濾波器的截止角頻率不是像巴特沃斯濾波器中所規(guī)定的（-3db）處角頻率，而是通帶邊緣的頻率。若波紋參數(shù)滿足，可以求的-3db處的角頻率為 （2.28）將式（2.27）表示的對歸一化，得到切比雪夫i型

26、低通原型濾波器的系統(tǒng)函數(shù) （2.29）對不同的n，式（2.29）的分母多項式已制成表格，供設計參考。和butterworth低通模擬濾波器設計一樣，若給定性能指標要求：，確定chebyshev低通模擬濾波器最小階數(shù)n和截止頻率（-3db頻率）。2.2.2.1 chbbyshevi型由式可得 故階數(shù)n可由下式求得 式中，截至頻率由上面兩式用matlab 編程計算濾波器最小階數(shù)n和截止頻率 2.2.2.2 chbbyshevii型chbbyshevii型通帶內(nèi)是平滑的，而阻帶具有等波紋起伏特性。因此，在階數(shù)n的計算公式上是相同的，而-3db截止頻率則不同。2.2.3 橢圓濾波器及最小階數(shù)的選擇橢圓

27、的模擬低通濾波器圓形的平方幅值響應函數(shù)為 式中，為小于1的正書，表示波紋情況；為截止頻率；）為橢圓函數(shù)，定義為當n為偶數(shù)（n=2m）時,當n為奇數(shù)（n=2m+1）時, 其中橢圓模擬濾波器特點是：在通帶和阻帶內(nèi)均具有等波紋起伏特性。何以上濾波器相比，相同的性能指標所需要的階數(shù)最小。但頻率響應應具有明顯的非線性。由式濾波器的階數(shù)可由下式確定， ， 式中 由上式計算濾波器的最小階數(shù)n和截止頻率。2.2.4貝塞爾濾波器貝塞爾模擬低通濾波器原型的特點是在零頻時具有最平坦的群延遲，并在整個通帶內(nèi)延遲幾乎不變。在零頻時的群延遲為。由于這一特點，貝塞爾模擬濾波器通帶內(nèi)保持信號形狀不變。濾波器傳遞漢書具有下面形

28、式3 matlab仿真3.1 matlab簡介整個matlab系統(tǒng)有五個主要部分: matlab語言。它是基于矩陣/數(shù)組的高級語言，它包括流程控制語句、函數(shù)、數(shù)據(jù)結構和輸入/輸出等，它還具有面向對象編程的特點。它既適合編寫小巧玲瓏的程序，也適合于開發(fā)復雜的大型應用程序。 matlab工作環(huán)境。它集成了一系列的工具和應用，方便用戶管理環(huán)境變量，輸入/輸出數(shù)據(jù)，開發(fā)、管理、調(diào)試用戶自己的m-文件以及matlab的應用程序。 圖形處理。它既包括二維和三維的數(shù)據(jù)可視化、圖像處理、動畫等高層指令，也包括低層的繪圖指令，允許用戶為應用程序設計自己的用戶圖形界面。 matlab數(shù)學函數(shù)庫。它包括數(shù)量龐大的計

29、算函數(shù)，從簡單的基本函數(shù)到復雜的矩陣求逆，矩陣的特征值，貝塞爾函數(shù)和快速傅里葉變換等。 matlab應用程序界面(api)。它是一組動態(tài)的庫函數(shù)，使得用戶在自己的c和fortran程序中可以和matlab交互，調(diào)用matlab的動態(tài)鏈接庫作計算。matlab語言的特點是:起點高 每個變量代表一個矩陣，它可以有n*m個元素;(2)每個元素都看作復數(shù)，這個特點在其它語言中不多見;(3)所有的運算，包括加、減、乘、除、函數(shù)運算都對矩陣和復數(shù)有效。人機界面適合科技人員 語言規(guī)則與科技人員的書寫習慣相近，因此易讀易寫，易于在科技人員之間交流; 矩陣的行數(shù)、列數(shù)無需定義:若要輸入一個矩陣，在用其它語言編程

30、時必須先定義矩陣的階數(shù)，而用matlab語言則不必有階數(shù)定義語句，輸入數(shù)據(jù)的行列數(shù)就決定了它的階數(shù); 鍵入算式立即得到結果，無需編譯:matlab是以解釋方式工作的，即它對每條語句解釋后立即執(zhí)行。若有錯誤也立即作出反應，便于編程者馬上改正。這些都大大減輕了編程和調(diào)試的工作量。強大面簡易的作圖功能能根據(jù)輸入數(shù)據(jù)自動確定繪圖坐標;能繪制三維坐標中的曲線和曲面;可設置不同顏色、線型、視角等;如果數(shù)據(jù)齊全，一條命令即可畫出圖來。智能化程度高、功能豐富、可擴展性強繪圖時自動選擇最佳坐標以及按輸入或輸出變元數(shù)自動選擇算法等;做數(shù)值積分時自動按精度選擇步長;自動檢測和顯示程序錯誤的能力強，易于調(diào)試。matl

31、ab軟件包括基本部分和專業(yè)擴展部分?；静糠职?矩陣的運算和各種變換、代數(shù)和超越方程的求解、數(shù)據(jù)處理和傅里葉變換、數(shù)值積分等等.專業(yè)擴展部分稱為工具箱.它實際上是用matlab的基本語句編成的各種子程序集，用于解決某一方面的專門問題，或實現(xiàn)某一類的新算法。易擴展性是matlab最重要的特點，每一個matlab用戶都可以成為對其有貢獻的人。在matlab的發(fā)展過程中，許多科學家、數(shù)學家、工程技術人員用它開發(fā)出了一些新的、有價值的應用程序，所有的程序完全不需要使用低層代碼來編寫。通過這些工作，已經(jīng)發(fā)展起來的工具箱有控制系統(tǒng)、信號處理、圖像處理、系統(tǒng)辨識、模糊集合、神經(jīng)元網(wǎng)絡、小波分析等20余個。

32、如果使用matlab來開發(fā)光學方面的應用程序，在不久的將來，也可能出現(xiàn)專門用來解決光學問題的工具箱。3.2 對低通模擬濾波器的仿真第三章所討論的設計思想，在matlab對模擬濾波器進行仿真的過程中依然適用，其具體步驟總結如下：（1） 確定模擬濾波器的性能指標，如截止頻率（對于低通和高通）或上、下邊界頻率,；波紋特性；帶阻衰減等。（2） 確定濾波器階數(shù)。（3） 計模擬低通濾波原型濾波器。matlab信號處理工具箱的濾波器原型函數(shù)buttap,cheb1ap。（4） 按頻率變換設計模擬濾波器（低通、高通、帶通、帶阻）。matlab信號處理工具箱的頻率變換函數(shù)lp2lp,lp2hp,lp2bp,lp

33、2bs。但是，按照這種設計思想的編程較為麻煩。matlab信號處理工具箱還提供模擬濾波器的完全設計函數(shù)：butter,cheby1等。用戶只需調(diào)用一次設計函數(shù)就可自動完成全部設計過程，編程十分簡單。下面將以模擬低通butterworth濾波器的設計為例，解釋這種設計方法。設計指標：通帶截止頻率=200，阻帶截止頻率=300，通帶衰減=1db，阻帶衰減=16db。仿真結果見圖3.1%matlab program2.1%design a buttworth analog lowpass filterws=300*pi;wp=200*pi;rp=1;rs=16;%compute oder and c

34、uttoff frequency n,wn=buttord(wp,ws,rp,rs,'s') fc=wn/(2*pi) b,a=butter(n,wn,'s'); %outputw=linspace(1,3000,1000)*2*pi; h=freqs(b,a,w); magh=abs(h); phah=unwrap(angle(h); plot(w/(2*pi),20*log10(magh); xlabel('frequency(hz)'); ylabel('magnidute(db)'); grid onn = 7wn =72

35、5.7292fc =115.5034圖3.1 模擬低通濾波器仿真圖3.1 模擬低通濾波器仿真3.3 模擬高通濾波器的仿真函數(shù)cheby1用于chebyshev i 型模擬濾波器的設計。調(diào)用格式為：b,a=cheby1(n,) b,a=cheby1(n,) 其中，為通帶波紋（db）, 為濾波器截止頻率，s為模擬濾波器，確省時為數(shù)字濾波器。本例給出利用chebyshev函數(shù)設計模擬高通的設計方法：技術指標：通帶截止頻率=1500hz，阻帶截止頻率=1000hz，通帶衰減=1db，阻帶衰減=20db。仿真結果見圖3.2ws=1000*pi; wp=150000*pi; rp=1; rs=20; n,

36、wn=cheb1ord(wp,ws,rp,rs,'s')b,a=cheby1(n,rp,wn,'high','s');w=linspace(1,3000,1000)*2*pi;h=freqs(b,a,w);magh=abs(h);phah=unwrap(angle(h);plot(w/(2*pi),20*log10(magh);xlabel('frequency(hz)');ylabel('magnidute(db)');grid onn =1wn = 4.7124e+005圖3.2 模擬高通濾波器仿真3.4 模擬

37、帶通濾波器的仿真函數(shù)butter用于butterworth濾波器設計，調(diào)用格式： b,a=butter(n,s) b,a=butter(n, ,ftype,s) 其中，n為濾波器階數(shù)；為濾波器截止頻率，s為模擬濾波器，確省時為數(shù)字濾波器。 ftype濾波器類型： high為高通濾波器，截止頻率； stop為帶阻濾波器，=（）； ftype缺省時為低通或帶通濾波器。以下設計一個butterworth模擬帶通濾波器，設計指標為：逼近頻率10002000hz,兩側過渡帶寬500hz，通帶衰減1db，阻帶衰減大于100db。結果見圖3.3。%matlab program4.1%design a but

38、terworth analog bandpass filterwp=1000 2000*2*pi; ws=500 2500*2*pi; rp=1; rs=100; n,wn=buttord(wp,ws,rp,rs,'s') b,a=butter(n,wn,'s'); n =23wn =1.0e+004 *0.6220 1.2695圖3.3 模擬帶通濾波器仿真3.5 對帶阻模擬濾波器的仿真這里設計一個chebyshev i 型模擬帶阻濾波器，設計指標為：阻帶頻率1000hz2000hz，兩側過渡帶寬500hz，通帶衰減1db，阻帶衰減大于50db。仿真結果見圖3.

39、4。%matlab program5.1%design a chebyshev i analog bandstop filterws=1000 2000*2*pi; wp=500 2500*2*pi; rp=1; rs=50; %compute oder and cuttoff frequencyn,wn=cheb1ord(wp,ws,rp,rs,'s') b,a=cheby1(n,rp,wn,'stop','s'); %outputw=linspace(1,3000,1000)*2*pi; h=freqs(b,a,w); magh=abs(h)

40、; phah=unwrap(angle(h); plot(w/(2*pi),20*log10(magh); xlabel('frequency(hz)'); ylabel('magnidute(db)'); grid on圖3.4 模擬帶阻濾波器仿真4 頻率轉換上述模擬濾波器的設計，只是討論了低通濾波器的設計問題，高通、帶通、帶阻濾波器可以通過對濾波器特性的頻率變換，轉換成低通濾波器的設計。這種頻率變換的方法又稱原型變換，變換得到的低通濾波器稱為低通原型濾波器。頻率變換是指低通原型傳遞函數(shù)與其他類型（高通、帶通、帶阻）濾波器傳遞函數(shù)中頻率之間的轉換關系。具體做法

41、是：先根據(jù)對高通、帶通、帶阻等濾波器特性指標要求，導出相應的低通原型的指標來，確定低通原型的，再根據(jù)一定變換關系得出高通、帶通、帶阻濾波器的。4.1低通至高通的轉換設低通濾波器傳遞函數(shù)為，角頻率為，截止頻率為；高通濾波器傳遞函數(shù)為，角頻率，通帶始點角頻率為。設有如下的變換關系 （4.1）而，有 （4.2）上市表明：s平面中的虛軸正好映射到p平面的虛軸上，其頻率變換關系為 （4.3）與之間存在 ， ， ， 相應的關系可表示成如圖4.1所示的曲線0圖4.1 通至高通的頻率變換系 根據(jù)上述的頻率變換，將高通濾波器的特性指標：高通通帶始點頻率，阻帶始點頻率，分別代入式（4.3）中，求出低通原型的通帶截

42、止頻率，阻帶始點頻率；而高通的通帶衰減及阻帶衰減即為對低通原型通帶與阻帶的要求。根據(jù)、和確定低通原型傳遞函數(shù)，即可求出高通傳遞函數(shù)為 4.2低通至帶通的變換設低通原型的傳遞函數(shù)為，角頻率為，截止角頻率為；帶通通帶中心頻率為，帶通濾波器的傳遞函數(shù)為，角頻率與為通帶的上、下邊界頻率，通帶帶寬b=。低通原型與帶通傳遞函數(shù)的變換關系為 （4.4）而，有 （4.5） 說明s平面中的虛軸證號映射到p平面的虛軸上，并有下列頻率變換關系存在 （4.6）其頻率變換關系曲線如圖4.2所示，由上式曲線應對稱于原點，圖中只給出了的部分。圖4.2 低通至帶通的頻率變換關系由式（4.6）和圖4.2，有 顯然，變換的結果是

43、：把低通原型由直流至截止頻率通帶范圍內(nèi)的頻率特性平移至帶通濾波器的之間；而由直流至截止頻率通帶范圍內(nèi)的頻率特性平移至帶通濾波器的之間。由于與變換是非線性的，因此。由圖4.2和式（4.6）可見，低通原型與帶通之間的頻率變換有如下關系： 解上述方程，可得 （4.7） （4.8）上兩式表明：帶通濾波器的中心頻率是其上、下邊界頻率的幾何平均值，寬帶與低通原型的通帶寬度相等。利用頻率變換方法求解帶通傳遞函數(shù)的具體思路是：由給定帶通指標：通帶寬度b、帶通通帶中心頻率、通帶衰減、阻帶始點頻率與、阻帶衰減等，求出低通原型指標及阻帶始點頻率；然后按低通原型指標確定低通原型傳遞函數(shù)；再代入（4.4）可以得到帶通傳遞函數(shù)為 實際上，也可以講一個低通濾波器與一個高通濾波器相聯(lián)接，只要低通濾波器的截止頻率大于高通濾波器的截止頻率即可，如圖4.3所示。0 0 0 低通高通圖4.3 低通與高通相聯(lián)接組成帶通濾波器原理4.3低通至帶阻的變換變換的方法與上述兩種濾波器相似，尤其與帶通濾波器的變換類似。若低通原型的傳遞函數(shù)是，帶阻濾波器的傳遞函數(shù)是，兩者之間有以下的變換關系： （4.9）式中的為阻帶的中心頻率，代入，有 可得到低通原型與帶阻濾波器之間的頻率轉換關系為 （4.10）根據(jù)上述的變換關系，