【正文】 ing application. This paper introduces the concept of digital filter, classification and design requirements. Then using MATLAB language programming to design IIR digital filter. Key design Chebyshev type I and II digital Chebyshev lowpass filter, and introduces optimization design.Keywords: Digital Filter 。采用Matlab設(shè)計濾波器，使原來非常繁瑣復(fù)雜的程序設(shè)計變成了簡單的函數(shù)調(diào)用，為濾波器的設(shè)和實現(xiàn)開辟了廣闊的天地，尤其是Matlab工具箱使各個領(lǐng)域的研究人員可以直觀方便地進行科學(xué)研究與工程應(yīng)用。 在設(shè)計中，脈沖響應(yīng)不變法中相位響應(yīng)有嚴格的線性，不存在穩(wěn)定性的問題，設(shè)計簡單。橢圓濾波器的選擇性相對前三種是最好的，下降斜度比較大，通帶和阻帶內(nèi)均為等波紋，同樣的性能指標(biāo)，橢圓濾波器可以最低的階數(shù)來實現(xiàn)。設(shè)計切比雪夫I型和切比雪夫II型數(shù)字低通濾波器，要求通帶邊界頻率fp=，通帶最大衰減Rp=；阻帶邊界頻率fs=8kHZ，阻帶最小衰減Rs=30dB,采樣頻率為Fs=20kHZ。其中固有頻率Wn必須是歸一化頻率。歸一化頻率的范圍在[0, 1]之間。在使用各類濾波器函數(shù)時應(yīng)當(dāng)注意以下重點:（1）階數(shù)和固有頻率的選擇:[N,Wn]=buttord(Wp,Ws,Rp,Rs)可得到符合要求性質(zhì)的濾波器的最小階數(shù)N以及數(shù)字Butterworth濾波器的固有頻率Wn(即3dB )。切比雪夫 II型濾波器的下降斜度比I型小，但其阻帶內(nèi)波紋較大?？梢栽O(shè)計低通、高通、帶通和帶阻的數(shù)字和模擬切比雪夫I型濾波器，其通帶內(nèi)為等波紋，阻帶內(nèi)為單調(diào)。主要有以下幾種方法： 這種方法也稱為零極點累試法。低通巴特沃斯濾波器設(shè)計步驟如下：（1）確定階數(shù)N。在許多國際一流學(xué)術(shù)刊物上，（尤其是信息科學(xué)刊物），都可以看到Matlab的應(yīng)用。數(shù)字信號處理常用的函數(shù)有波形的產(chǎn)生、濾波器的分析和設(shè)計、傅里葉變換、Z變換等，如：波形產(chǎn)生：sawtooth（鋸齒波或三角波）Diric（Dirichlet或周期sinc函數(shù)）rand（白噪聲信號波形）square（方波）sinc（sinc或 函數(shù)）chirp（chirp信號波形）濾波器的分析：abs（求幅值）angle（求相角）conv（求卷積）freqz（數(shù)字濾波器頻率響應(yīng)）impz（數(shù)字濾波器的沖擊響應(yīng)）FIR濾波器設(shè)計：triang（三角窗）blackman（布萊克曼窗）boxcar（矩形窗）hamming（海明窗）hanning（漢寧窗）kaiser（凱塞窗）fir1（基于窗函數(shù)法）fir2（基于頻率抽樣法）firrcos（上升余弦FIR濾波器設(shè)計法）intfilt（內(nèi)插FIR濾波器設(shè)計法）kaiserord（用Kaiser窗設(shè)計FIR濾波器的參數(shù)估計）各種變換：czt（線性調(diào)頻Z變換）dct（離散余弦變換）fft（一維快速傅里葉變換）fft2（二維快速傅里葉變換）idct（逆離散余弦變換）ifft（一維逆快速傅里葉變換）ifft2（二維逆快速傅里葉變換）hilbert（Hilbert變換）可以看出Matlab是一個功能十分強大的系統(tǒng)，是集數(shù)值計算、圖形管理、程序開發(fā)為一體的環(huán)境。功能性工具箱主要用來擴充Matlab的符號計算、圖形可視化、建模仿真、文字處理等功能以及與硬件實時交互的功能。f .有豐富的網(wǎng)絡(luò)資源，從相關(guān)的Web網(wǎng)站可以直接獲得全套的Matlab聯(lián)機幫助文件和說明書的電子文檔，還可以獲得各類技術(shù)支持與幫助。、統(tǒng)計、科學(xué)及工程方面的函數(shù)，使用簡單快捷，并且有很強的用戶自定義函數(shù)的能力。g. Matlab的缺點是，它和其他高級程序相比，程序的執(zhí)行速度較慢。c. MatlabB既具有結(jié)構(gòu)化的控制語句（如for循環(huán)，while循環(huán)，break語句和if語句），又有面向?qū)ο缶幊痰奶匦?。Matlab用更直觀的，符合人們思維習(xí)慣的代碼，代替了C和FORTRAN語言的冗長代碼。3 濾波器的Matlab設(shè)計 Matlab概述 Matlab語言介紹Matlab是矩陣實驗室（Matrix　Laboratory）之意。這是因為S平面與Z平面是單值的一一對應(yīng)關(guān)系。其次，將s=σ+jΩ代入式（212），得 （213）因此 （214）由此看出，當(dāng)σ0時，|z|1；當(dāng)σ0時，|z|1。圖4雙線性變換的映射關(guān)系為了將S平面的整個虛軸jΩ壓縮到S1平面jΩ1軸上的π/T到π/T段上，可以通過以下的正切變換實現(xiàn) （26） 式中,T仍是采樣間隔。這是因為從S平面到Ｚ平面是多值的映射關(guān)系所造成的。至于高通和帶阻濾波器，由于它們在高頻部分不衰減，因此將完全混淆在低頻響應(yīng)中。 脈沖響應(yīng)不變法的優(yōu)缺點從以上討論可以看出，脈沖響應(yīng)不變法使得數(shù)字濾波器的單位脈沖響應(yīng)完全模仿模擬濾波器的單位沖激響應(yīng)，也就是時域逼近良好，而且模擬頻率Ω和數(shù)字頻率ω之間呈線性關(guān)系ω=ΩT。這時數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，而帶有一定的失真。應(yīng)序列h(n)模仿模擬濾波器的沖激響應(yīng)ha(t)，即將ha(t)進行等間隔采樣，使h(n)正好等于ha(t)的采樣值，滿足h(n)=ha(nT) 式,T是采樣周期。（4）用所選的轉(zhuǎn)換方法，將模擬濾波器轉(zhuǎn)換成數(shù)字低通濾波器系統(tǒng)函數(shù)。其設(shè)計步驟是：先設(shè)計過渡模擬濾波器得到系統(tǒng)函數(shù)，然后將按某種方法轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)。其設(shè)計步驟是，先設(shè)計模擬濾波器，再按照某種方法轉(zhuǎn)換成數(shù)字濾波器。2 IIR濾波器的設(shè)計 IIR濾波器介紹IIR數(shù)字濾波器是一種離散時間系統(tǒng)，其系統(tǒng)函數(shù)為 （21）假設(shè)M≤N，當(dāng)M＞N時,系統(tǒng)函數(shù)可以看作一個IIR的子系統(tǒng)和一個(MN)的FIR子系統(tǒng)的級聯(lián)。實現(xiàn)IIR濾波器的階次較低，所用的存儲單元較少，效率高，精度高，而且能夠保留一些模擬濾波器的優(yōu)良特性，因此應(yīng)用很廣。數(shù)字濾波器是數(shù)字信號處理的關(guān)鍵部分之一，它是語音處理、圖像處理、軟件無線電等應(yīng)用中最基本的組成部分。尤其是Matlab中的信號處理工具箱使各個領(lǐng)域的研究人員可以直觀方便地進行科學(xué)研究與工程應(yīng)用。長治學(xué)院學(xué)士學(xué)位論文（設(shè)計）長 治 學(xué) 院2013屆學(xué)士學(xué)位畢業(yè)論文數(shù)字濾波器的設(shè)計學(xué) 號： 09407232 姓 名： 楊 婕 指導(dǎo)教師： 上官晉太 專 業(yè)： 計算機科學(xué)與技術(shù) 系 別： 計算機系 完成時間：2013年5月數(shù)字濾波器的設(shè)計專業(yè)：計算機科學(xué)與技術(shù) 姓名：楊 婕 學(xué)號：09407232指導(dǎo)教師：上官晉太摘 要：現(xiàn)代通信系統(tǒng)中，由于信號中經(jīng)?；煊懈鞣N復(fù)雜成分，所以很多信號分析都是用濾波器進行，而數(shù)字濾波器是通過數(shù)值運算實現(xiàn)濾波，具有處理精度高、穩(wěn)定、靈活、不存在阻抗匹配問題，可以實現(xiàn)模擬濾波器無法實現(xiàn)的特殊濾波功能。Matlab軟件以矩陣運算為基礎(chǔ)，把計算、可視化及程序設(shè)計有機融合到交互式工作環(huán)境中，并且為數(shù)字濾波的研究和應(yīng)用提供了一個直觀、高效、便捷的利器?？梢酝ㄟ^基本的濾波器二階通用濾波器傳遞函數(shù)，推導(dǎo)出最通用的濾波器類型：低通、帶通、高通、帯阻和橢圓型濾波器。數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時域特性，可分為兩種，即無限長沖激響應(yīng)(IIR)數(shù)字濾波器和有限長沖激響應(yīng)(FIR)數(shù)字濾波器。首先介紹了數(shù)字濾波器的概念，分類以及設(shè)計要求，接著利用MATLAB函數(shù)語言編程來設(shè)計濾波器，重點設(shè)計Chebyshev I型和Chebyshev II型數(shù)字低通濾波器,并介紹最優(yōu)化設(shè)計。 IIR數(shù)字濾波器設(shè)計步驟IIR數(shù)字濾波器的設(shè)計一般有兩種方法：一個是利用模擬濾波器的設(shè)計方法進行。另外一種直接在頻率或者時域內(nèi)進行，由于需要解聯(lián)立方程，設(shè)計時需要計算機做輔助設(shè)計。（3）按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計過渡模擬低通濾波器。使數(shù)字濾波器的單位脈沖響應(yīng)完全模仿模擬濾波器的單位沖激響應(yīng)，也就是時域逼近良好，一個線性相位的模擬濾波器通過脈沖響應(yīng)不變法得到的仍是一個線性相位的數(shù)字濾波器。正如采樣定理所討論的，只有當(dāng)模擬濾波器的頻率響應(yīng)是限帶的，且?guī)抻谡郫B頻率以內(nèi)時，即 (24)才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應(yīng)，而不產(chǎn)生混疊失真，即 ， (25)但任何一個實際的模擬濾波器頻率響應(yīng)都不是嚴格限帶的，變換后就會產(chǎn)生周期延拓分量的頻譜交疊，即產(chǎn)生頻率響應(yīng)的混疊失真。圖3脈沖響應(yīng)不變法中的頻響混疊現(xiàn)象對某一模擬濾波器的單位沖激響應(yīng)ha(t)進行采樣，采樣頻率為fs，若使fs增加，即令采樣時間間隔（T=1/fs）減小，則系統(tǒng)頻率響應(yīng)各周期延拓分量之間相距更遠，因而可減小頻率響應(yīng)的混疊效應(yīng)。所以，脈沖響應(yīng)不變法只適用于限帶的模擬濾波器(例如，衰減特性很好的低通或帶通濾波器)，而且高頻衰減越快，混疊效應(yīng)越小。 雙線性變換法設(shè)計IIR數(shù)字濾波器 設(shè)計原理脈沖響應(yīng)不變法的主要缺點是產(chǎn)生頻率響應(yīng)的混疊失真。這樣就使S平面與Z平面建立了一一對應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖4 所示。首先,把z=ejω，可得 （212）即S平面的虛軸映射到Z平面的單位圓。雙線性變換法與脈沖響應(yīng)不變法相比，其主要的優(yōu)點是避免了頻率響應(yīng)的混疊現(xiàn)象。由圖4看出，在零頻率附近，模擬角頻率Ω與數(shù)字頻率ω之間的變換關(guān)系接近于線性關(guān)系；但當(dāng)Ω進一步增加時，ω增長得越來越慢，最后當(dāng)Ω→∞時，ω終止在折疊頻率ω=π處，因而雙線性變換就不會出現(xiàn)由于高頻部分超過折疊頻率而混淆到低頻部分去的現(xiàn)象，從而消除了頻率混疊現(xiàn)象。Matlab最突出的特點就是簡潔。f. Matlab的圖形功能強大。 Matlab的功能Matlab包含的內(nèi)容非常豐富，功能強大，可以概括為以下幾個方面：，如DOS、Windows 95/98/2000/2000/NT、Compaq Alpha、LinuxSun Solaris等。e. 有從外部文件及外部硬件設(shè)備讀入數(shù)據(jù)的能力。其工具箱可以分為兩大類，即功能性工具箱和科學(xué)性工具箱。在數(shù)字信號處理中常用的算法，如FFT，卷積，相關(guān)，濾波器設(shè)計，參數(shù)模型等，幾乎都只用一條語句即可調(diào)用。在國際學(xué)術(shù)界，Matlab已經(jīng)被確認為準確、可靠的科學(xué)計算標(biāo)準軟件。 IIR數(shù)字濾波器的典型設(shè)計法設(shè)計巴特沃斯數(shù)字低通濾波器和橢圓數(shù)字低通濾波器，要求通帶邊界頻fp=，通帶最大衰減Rp=；阻帶邊界頻率fs=8kHZ，阻帶最小衰減Rs=30dB,采樣頻率為Fs=20kHZ。雙線性變換法設(shè)計步驟如下：（1）首先寫出該濾波器的系統(tǒng)函數(shù)（2）利用雙線性變換法轉(zhuǎn)換，數(shù)字濾波器的系統(tǒng)函數(shù)為 用雙線性變換法設(shè)計橢圓數(shù)字低通濾波器的M程序（見附錄代碼2）： 圖6 雙線性變換法設(shè)計的橢圓數(shù)字低通濾波器運行波形運行結(jié)果：N=2b