電子信息工程專業課程設計題目：Iir濾波器設計班級：電信071姓名：袁晶學號：200709802摘 要摘 要：本論文首先介紹了濾波器的濾波原理以及模擬濾波器、數字濾波器的設計方法。重點介紹了IIR數字濾波器的設計方法。即脈沖響應不變法和雙線性變換法。在此基礎上，用DSP虛擬實現任意階IIR數字濾波器。此設計擴展性好，便于調節濾波器的性能，可以根據不同的要求在DSP上加以實現。關鍵字： DSP 濾波、IIR數字低通濾波器第1章 濾波器簡介從廣義上講，任何對某些頻率（相對于其他頻率來說）進行修正的系統稱為濾波器。嚴格地講，對輸入信號通過一定的處理得到輸出信號，這個處理通常是提取信號中某頻率范圍內的信號成分，把這種處理的過程稱為濾波。實現濾波處理的運算電路或設備稱為濾波器。在許多科學技術領域中，廣泛應用線性濾波和頻譜分析對信號進行加工處理，模擬濾波是處理連續信號，數字濾波則是處理離散信號，而后者是在前者的基礎上發展起來的。我們知道，無源或有源模擬濾波器是分立元件構成的線性網絡，他們的性能可以用線性微分方程來描述，而數字濾波器是個離散線性系統，要用差分方程來描述，并以離散變換方法來分析。這些方程組可以用專用的或通用的數字計算機進行數字運算來實現。因此，數字濾波器的濾波過程是一個計算過程，它將輸入信號的序列數字按照預定的要求轉換成輸出數列。1.1 濾波器的工作原理1.1.1 模擬濾波器的工作原理我們知道，模擬濾波器是對模擬信號實行線性濾波的一種線性時不變系統，如圖1.1所示。在時域內，它的動態特性可以用系統的單位沖激函數的響應來描述，也就是該濾波系統在任何時刻對輸入單位沖激信號=（t）的輸出響應。這個函數從時域上反映了該濾波系統的傳輸特性。對于任意輸入信號，系統的輸出可以卷積表示： = (1.1)上式表明在對線性濾波器系統進行時域分析時，采用了疊加原理，先將任意輸入信號波形分成不同時間的窄脈沖之和，再分別求出各個脈沖通過濾波器之后的響應，并進行線性疊加從而得到總的輸出信號。模擬濾波器,H(s)圖1.1模擬濾波器原理在頻域分析時，線性濾波器的轉移函數等于系統的單位沖激函數的響應的拉普拉斯變換： （1.2）很明顯，當s=j，上式就是傅立葉變換的表達式，它反映了濾波器的傳輸特性對各種頻率的響應，也就是濾波器的頻率響應函數，它決定著濾波特性。當濾波器輸入信號與輸出信號的拉普拉斯變換，得 （1.3）這表明兩信號卷積的變換等于各自變換的乘積。在頻譜關系上，一個輸入信號的頻譜，經過濾波器的作用后，被變換成的頻譜。因此，根據不同的濾波要求來選定，就可以得到不同類型的模擬濾波器。還可以看出，濾波器的濾波過程就是完成信號與它的單位沖激函數響應之間的數學卷積運算過程。1.1.2 數字濾波器的工作原理在數字濾波中，我們主要討論離散時間序列。如圖1.2所示。設輸入序列為，離散或數字濾波器對單位抽樣序列的響應為。因在時域離散信號和系統中所起的作用相當于單位沖激函數在時域連續信號和系統中所起的作用。數字濾波器,H(z)圖1.2 數字濾波器原理數字濾波器的序列將是這兩個序列的離散卷積，即 （1.4） 同樣，兩個序列卷積的z變換等于個自z變換的乘積，即 （1.5）用代入上式，其中T為抽樣周期，則得到 （1.6）式中和 分別為數字濾波器輸入序列和輸出序列的頻譜，而為單位抽樣序列響應的頻譜。由此可見，輸入序列的頻譜經過濾波后，變為 ，按照的特點和我們處理信號的目的，選取適當的使的濾波后的符合我們的要求。1.2 濾波器的基本特性1.2.1 模擬濾波器與數字濾波器的基本特性如利用模擬電路直接對模擬信號進行處理則構成模擬濾波器，它是一個連續時間系統。如果利用離散時間系統對數字信號（時間離散、幅度量化的信號）進行濾波則構成數字濾波器。數字濾波器的差分方程表示為： 系統函數表示:數字濾波器的特性通常用其頻率響應函數來描述，包括幅度特性和相位特性。按信號通過系統時的特性（主要是幅頻特性）來分類：可以有低通、高通、帶通和帶阻四種基本類型。（1） 低通數字濾波器：圖1.3所示|H(ej)|c-c-22-fs/2-fsfs/2fs-fcfcf圖1.3 低通數字濾波器的頻譜（2） 高通數字濾波器：圖1.4所示|H(ej)|c-c-22圖1.4 高通數字濾波器的頻譜（3） 帶通數字濾波器：圖1.5所示|H(ej)|2-2-221-1圖1.5 帶通數字濾波器的頻譜（4）帶阻數字濾波器：圖1.6所示1|H(ej)|2-2-22-1 圖1.6 帶阻數字濾波器的頻譜 其他較復雜的特性可以由基本濾波器組合。1.2.2 無限沖擊響應IIR和有限沖擊響應FIR濾波器按系統沖擊響應（或差分方程）可以分成無限沖擊響應 IIR和有限沖擊響應FIR濾波器兩類。這兩種濾波器都可以現實各種頻率特性要求，但它們在計算流程、具體特性逼近等方面是有差別的。(1) FIR濾波器(非遞歸型):(2) IIR濾波器（遞歸型） 還有一些其他的分類方法，例如在特定場合使用的濾波器。1.3 濾波器的主要技術指標濾波器的主要技術指標取決于具體的應用或相互間的相互關系。具體的有最大通帶增益（即通帶允許起伏）；最大阻帶增益；通帶截止頻率；阻帶截止頻率。如圖1.7所示dB0dBdB圖1.7 濾波器的主要技術指標第2章 模擬濾波器的設計模擬濾波器的理論和設計方法已經發展的相當成熟，且有若干典型的模擬濾波器供我們選擇，如巴特沃斯（Butterworth濾波器.切比雪夫（Chebyshev）濾波器等。這些工作的理論分析和設計方法在20世紀30年代就完成，然而煩瑣.冗長的數字計算使它難以付諸實用。直到50年代，由于計算機技術的逐步成熟，求出大量設計參數和圖表，這種方法才得到廣泛應用。這些典型的濾波器各有特點：巴特沃斯濾波器具有單調下降的幅頻特性；切比雪夫濾波器的幅頻特性在通帶或者阻帶有波動發，可以提高選擇性。這樣根據具體要求可以選擇不同類型的濾波器。模擬濾波器按幅度特征可以分成低通、高通、帶通和帶阻濾波器。它們的理想幅度特性如圖2.1所示，但我們設計濾波器時，總是先設計低通濾波器，再通過頻率變換將低通濾波器轉換成希望類型的濾波器帶通帶阻圖2.1 模擬濾波器理想幅度特性低通高通2.1 模擬濾波器的設計方法利用頻率變換設計模擬濾波器的步驟為：（1）給定模擬濾波器的性能指標，如截止頻率或上、下邊界頻率等。（2）確定濾波器階數（3）設計模擬低通原型濾波器。（4）按頻率變換設計模擬濾波器（低通、高通、帶通、帶阻）。模擬低通濾波器的設計指標有，和，其中和分別稱為通帶截止頻率和阻帶截止頻率。 是通帶（=0）中的最大衰減系數，是阻帶的最小衰減系數，和一般用dB表示。對于單調下降的幅度特性，可表示成： （2.1） （2.2）如果=0處幅度已歸一化為一，即，和表示為 （2.3） （2.4）以上技術指標用圖2.2表示，圖中稱為3dB 截止頻率，因,-20 10.7070圖2.2 低通濾波器的幅度特性濾波器的技術指標給定以后，需要設計一個傳輸函數，希望其幅度平方函數滿足給定的指標和，一般濾波器的單位沖激響應為實數，因此 = (2.5)如果能由，求出，那么就可以求出所需的，對于上面介紹的典型濾波器，其幅度平方函數有自己的表達式，可以直接引用。這里要說明的是必須是穩定的。因此極點必須落在s平面的左半平面，相應的的極點落在右半平面。2.2 模擬原型濾波器及最小階數的選擇2.2.1 巴特沃斯濾波器及最小階數的選擇巴特沃斯濾波器是最基本的逼近方法形式之一。它的幅頻特性模平方為 (2.6)式中N是濾波器的階數。當=0時，；當=時，是3dB截止頻率。不同階數N的巴特沃斯濾波器特性如圖2.3所示，這一幅頻特性具有下列特點：（1）最大平坦性：可以證明：在=0點，它的前（2N-1）階導數都等于0，這表明巴特沃斯濾波器在=0附近一段范圍內是非常平直的，它以原點的最大平坦性來逼近理想低通濾波器。“最平響應”即由此而來。（2）通帶，阻帶下降的單調性。這種濾波器具有良好的相頻特性。（3）3dB的不變性：隨著N的增加，頻帶邊緣下降越陡峭，越接近理想特性，但不管N是多少，幅頻特性都通過-3dB點。當時，特性以20NdB/dec速度下降。圖2.3 不同階數N的巴特沃斯濾波器特性現根據式（2.6）求巴特沃斯濾波器的系統函數Ha（s）。令=s/j，帶入式（2.6）對應的極點： (2.7)即為的極點，此極點分布有下列特點：(1)的2N個極點以/N為間隔均勻分布在半徑為的圓周上，這個圓稱為巴特沃斯圓。(2)所有極點以j軸為對稱軸成對稱分布，j軸上沒有極點。(3)當N為奇數時，有兩個極點分布在的實軸上；N為偶函數時，實軸上沒有極點。所有復數極點兩兩呈共軛對稱分布。圖2.4畫出了N=3時的極點分布。全部零點位于s=處。j圖2.4 N=3時Ha（s）Ha（-s）極點分布為得到穩定的，取全部左半平面的極點。 （2.8）當N為偶數時 （2.9）當N為奇數時 （2.10）為使用方便把式（2.9）和式（2.10）對進行歸一化處理，為此，分子分母各除以，并令，稱為歸一化復頻率：（N為偶數） （2.11） （N為奇數）（2.12）用歸一化頻率表示的頻率特性稱為原型濾波特性（即歸一化復頻率s 的虛部）。對式（2.6）所示的低通巴特沃斯特性用表示得到： （2.13）稱為巴特沃斯低通原型濾波器幅頻特性。在低通原型濾波頻率特性上，截止頻率=1。若給出模擬低通濾波器的設計性能指標要求：通帶邊界頻率，阻帶邊界頻率，通帶波紋，阻帶衰減，要確定butterworth ,，低通濾波器最小階數N及截止頻率。，的意義如圖所示。當=時， 即，以截至頻率（幅值下降3dB）為1，化為相對為相對的相對頻率由上式可寫為。同理，當=時， 。由此可見 N應向上取整，再用MATLAB 編程計算濾波器最小階數N和截止頻率。就是切比雪夫濾波器的極點，給定N，即可求的2N個極點分布。由式（2.22）實部與虛部的正弦和余弦函數平方約束關系可以看出，此極點分布滿足橢圓方程，其短軸和長軸分別為 （2.23）圖2.7畫出了N=3時切比雪夫濾波器的極點分布。abj 圖2.7極點所在的橢圓可以和半徑為a的圓和半徑為b的圓聯系起來，這兩個圓分別稱為巴特沃斯小圓和巴特沃斯大圓。N階切比雪夫濾波器極點的縱坐標，而橫坐標等于N階巴特沃斯小圓極點的橫坐標取左半平面的極點： k=1,2，N （2.24）則切比雪夫濾波器的系統函數： （2.25）其中，常數A=。因而切比雪夫濾波器的系統函數表示為： （2.27）切比雪夫濾波器的截止角頻率不是像巴特沃斯濾波器中所規定的（-3dB）處角頻率，而是通帶邊緣的頻率。若波紋參數滿足，可以求的-3dB處的角頻率為 （2.28）將式（2.27）表示的對歸一化，得到切比雪夫I型2.2.2低通原型濾波器的系統函數 （2.29）對不同的N，式（2.29）的分母多項式已制成表格，供設計參考。和butterworth低通模擬濾波器設計一樣，若給定性能指標要求：，確定Chebyshev低通模擬濾波器最小階數N和截止頻率（-3dB頻率）。2.2.2.1 Chbbyshev I型由式可得 故階數N可由下式求得 式中，截至頻率由上面兩式用Matlab 編程計算濾波器最小階數N和截止頻率 2.2.2.2 Chbbyshev II型Chbbyshev II型通帶內是平滑的，而阻帶具有等波紋起伏特性。因此，在階數N的計算公式上是相同的，而-3dB截止頻率則不同。2.2.3 橢圓濾波器及最小階數的選擇橢圓的模擬低通濾波器圓形的平方幅值響應函數為 式中，為小于1的正書，表示波紋情況；為截止頻率；）為橢圓函數，定義為當N為偶數（N=2m）時,當N為奇數（N=2m+1）時, 其中橢圓模擬濾波器特點是：在通帶和阻帶內均具有等波紋起伏特性。何以上濾波器相比，相同的性能指標所需要的階數最小。但頻率響應應具有明顯的非線性。由式濾波器的階數可由下式確定， ， 式中 由上式計算濾波器的最小階數N和截止頻率。2.2.4貝塞爾濾波器貝塞爾模擬低通濾波器原型的特點是在零頻時具有最平坦的群延遲，并在整個通帶內延遲幾乎不變。在零頻時的群延遲為。由于這一特點，貝塞爾模擬濾波器通帶內保持信號形狀不變。濾波器傳遞漢書具有下面形式第3章 IIR數字濾波器的設計3.1 IIR數字濾波器的設計方法 IIR濾波器是一種數字濾波器，濾波器的系統函數如式 （3.1）由于它的脈沖響應序列是無限長的，故稱無限沖激響應濾波器。IIR濾波器的設計就是根據濾波器某些性能指標要求，設計濾波器的分子和分母多項式。它和FIR濾波器相比優點是在滿足相同性能指標要求條件下，IIR濾波器的階數要明顯低于FIR濾波器。但IIR濾波器的相位是非線性的。IIR濾波器設計方法可分為三種：模擬濾波器變換（經典設計法）、直接設計法、參數模型設計法、最大平滑濾波器設計。3.2 IIR濾波器經典設計 基于模擬濾波器變換原理，首先是根據濾波器的技術指標設計出相應的模擬濾波器，然后再將設計好的模擬濾波器變換成滿足給定技術指標的數字濾波器。這就是IIR濾波器設計的經典法。經典法可設計出低通、高通、帶通、帶阻等各種IIR濾波器。在MATLAB中，經典法設計IIR數字濾波器采用下面的主要步驟：模擬低通濾波原型頻率變換模擬離散化IIR濾波器由上可見，經典設計法是利用模擬濾波器的設計成果。第二步完成后，一個達到期望性能指標的模擬濾波器（低通、高通、帶通和帶阻）已經設計出來。第三步離散化主要任務就是把模擬濾波器變換成數字濾波器，即把模擬濾波器的系數映射成數字濾波器的系統函數。數字濾波器的設計工作就全部完成。實現系統傳遞函數s域至z域映射有脈沖響應不變法和雙線性映射兩種方法。脈沖響應不變法利用模擬濾波器成熟的理論和設計方法來設計IIR數字低通濾波器是經常用的方法。設計過程是：按照技術要求設計一個模擬低通濾波器，得到模擬低通濾波器的傳輸函數轉換成數字低通濾波器的系統函數H(z)。這樣設計的關鍵問題就是找這種轉換關系，將s平面上的轉換成z平面上的H(z)。為了保證轉換后的H(z)穩定且滿足技術要求，對轉換關系提出兩點要求：（1）因果穩定的模擬濾波器轉換成數字濾波器，仍是因果穩定的。我們知道，模擬濾波器因果穩定要求其傳輸函數的極點全部位于s平面的左半平面；數字濾波器因果則要求H(z)的極點全部在單位圓內。因此，轉換關系應是s平面的左半平面映射z平面的單位圓內部。 （2）數字濾波器的頻率響應模仿模擬濾波器的頻響，s平面的虛軸映射z平面的單位圓，相應的頻率至間成線性關系。 將傳輸函數從s平面傳換到z平面的方法有多種，但工程上常用的是脈沖響應不變法和雙線性變換法。我們先研究脈沖響應不變法。設模擬濾波器的傳輸函數為，相應的單位沖激響應是，=LTLT.代表拉氏變換，對進行等間隔采樣，采樣間隔為T，得到，將h(n)= 作為數字濾波器的單位取樣響應，那么數字濾波器的系統函數H(z)便是h(n)的Z變換。因此脈沖響應不變法是一種時域上的轉換方法，它是h(n)在采樣點上等于。 設模擬濾波器只有單階極點，且分母多項式的階次高于分子多項式的階次，將用部分分式表示： （3.2）式中為的單階極點。將逆拉氏變換得到： （3.3）式中u(t)是單位階躍函數。對進行等間隔采樣，采樣間隔為T,得到： （3.4） 對上式進行Z變換，得到數字濾波器的系統函數H(z)： (3.5) 對比3.2、3.5式，的極點映射到z平面，其極點變成，系數不變化。下面我們分析從模擬濾波器轉換到數字濾波器，s平面和z平面之間的映射關系，從而找到這種轉換方法的優缺點。這里我們以采樣信號作為橋梁，推導其映射關系。設的采樣信號用表示， （3.6）對進行拉氏變換，得到： （3.7）式中是在采樣點t=n T時的幅度值，它與序列h(n)的幅度值相等，它與序列h(n)的幅度值相等，即h(n)= ，因此得到：上式表示采樣信號的拉氏變換與相應的序列的Z變換之間的映射關系可用下式表示： （3.8） 我們知道模擬信號的傅立葉變換和其采樣信號的傅立葉變換和其采樣信號的傅立葉變換之間的關系滿足式： （3.9）將代入上式，得 （3.10）由 （3.8）、（3.10）式得到： （3.11）上式表明將模擬信號的拉氏變換在s平面上沿虛軸按照周期延遲后，再按照（3.8）式映射關系，映射到z平面上，就得到H（z）。（3.8）式可稱為標準映射關系。下面進一步分析這種映射關系。設 按照（3.8）式，得到：因此得到： （3.12）那么 上式關系式說明，s平面的虛軸（）映射z平面的單位圓（），s平面左平面（）映射z平面單位內（），s平面右半平面映射平面單位圓外（）。這說明如果因果穩定，轉換后得到（）仍是因果穩定的。另外，注意到z=是一個周期函數，可寫成M為任意整數當不變，模擬頻率變化2/T的整數倍時，映射值不變。或者說，將s平面沿著j軸分割成一條條寬為2/T的水平帶，每條水平面都按照前面分析的映射關系對應著整個z平面。此時所在的s平面與所在的z平面的映射關系如圖所示。當模擬頻率從變化到時，數字頻率則從變化到，且按照（3.12）式，即與之間成線性關系。但是，從模擬信號到采樣信號，其拉氏變換要按照式。其為周期，沿虛軸方向進行周期化。如果原模擬信號得頻帶不是限于之間，則會在得奇數倍附近產生頻率混疊，從而映射到平面上，在附近產生頻率混疊。脈沖相應不變法的頻率混疊現象如圖所示。這種會使設計出的濾波器在附近的頻率特性，程度不同地偏離模擬濾波在附近的頻率特性，嚴重時使濾波器不滿足給定的技術指標。因此，希望設計的濾波器使帶限濾波器，如果不是帶限的，例如高通濾波器，帶阻濾波器，需要在高通帶阻濾波器之前加保護濾波器，濾除高于折疊頻率以上的頻帶，以避免產生頻率混疊現象。但這樣會增加系統的成本和復雜性，因此，高通與帶阻濾波器不適合用這種方法設計。假設沒有頻率混疊現象，即滿足將關系式代入，得到：說明用脈沖響應不變法設計的數字濾波器可以很好地重現原模擬濾波器的頻響。上式中，的幅度特性與采樣間隔成反比，這樣當T較小時，就會有太高的增益。為避免這一現象，令那么,此時一般的極點是一個復數，且以共軛成對的形式出現，在式中講一對復數共軛極點放在一起，形成一個二階基本節。如果模擬濾波器的二階基本節的形式為，極點為可以推導出相應的數字濾波器二階基本節（只有實數乘法）的形式為如果模擬濾波器二階基本節的形式為，極點為，則對應的數字濾波器二階基本節的具體形式為利用以上這些變換關系，可以簡化設計。雙線性變換法這種變換方法，采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到/T之間，再用轉換到z平面上。設Ha（s），s=j，經過非線性頻率壓縮后用Ha（s1），=j1表示，這里用正切變換實現頻率壓縮： （3.13）式中T仍是采樣間隔，當1從-/T經過0變化到/T時，則由-經過0變化到+，實現了s平面上整個虛軸完全壓縮到平面上虛軸的/T之間的轉換。這樣便有 （3.14）再通過 轉換到z平面上，得到： （3.15） （3.16）式（3.15）或式（3.16）稱為雙線性變換。從s平面映射到平面，再從平面映射到z平面，其映射情況如圖3.1所示。由于從s平面到平面具有非線性頻率壓縮的功能，因此不可能產生頻率混疊現象。另外，從平面轉換到z平面仍然采用標準轉換關系 ，平面的 /T之間水平帶的左半部分映射z平面單位圓內部，虛軸映射單位圓。這樣，Ha（s）因果穩定，轉換成的H（z）也是因果穩定的。下面分析模擬頻率和數字頻率之間的關系。令 ，并代入（3.15）式中，有 （3.17）ImRez平面js平面j11s1平面/T-/T圖3.1 雙線性變換法的映射關系上式說明，s平面上與平面上的成非線性正切關系，如圖3.2所示。在=0附近接近線性關系；當增加是，增加得愈來愈快；當趨近時，趨近于。正是因為這種非線性關系，消除了頻率混疊現象。-0圖3.2 雙線性變換法的頻率變換與之間的非線性關系是雙線性變換法的缺點，直接影響數字濾波器頻響逼真的模仿模擬濾波器的頻響，幅度特性和相位特性失真的情況如圖3.3所示。這種非線性影響的實質問題是：如果的刻度是均勻的，則影射到z平面的刻度不是均勻的，而是隨增加愈來愈密。000000圖3.3 雙線性變換法幅度和相位特性的非線性映射雙線性變換法可由簡單的代數公式（3.15）將Ha（s）直接轉換成H（z），這是該變換法的優點。但當階數稍高時，將H（z）整理成需要的形式，也不是一件簡單的工作。3.3 IIR濾波器直接設計 IIR數字濾波器的經典設計法只限于幾種標準的低通、高通、帶通、帶阻濾波器，而對于具有形狀或多頻帶濾波器的設計是無能為力的。 如果所設計IIR濾波器幅頻特性比較復雜，可采用最小二乘法擬合給定的幅頻響應，使設計的濾波器幅頻特性逼近期望的頻率特性，這種方法稱為IIR濾波器的直接設計法。 MATLAB信號處理工具箱函數YULEWALK采用直接設計IIR數字濾波器。函數YULEWALK采取下面步驟計算分子多項式：（1） 計算與分子多項式相應的幅值平方響應的輔助式；（2） 由輔助分子式和分母多項式計算完全的頻率響應；（3） 計算濾波器的脈沖響應；（4） 采用最小二乘法擬合脈沖響應最終求得濾波器的分子多項式系數。 函數YULEWALK允許我們自由定義濾波器的頻率向量f和幅值向量m，因此該函數可設計具有任意形狀的幅頻響應的濾波器，包括多頻帶濾器。3.4 IIR濾波器的MATLAB設計數字高通、帶通及帶阻IIR濾波器設計除了低通數字濾波器之外，實際中還常常需要高通、帶通及帶阻數字濾波器，這三種數字濾波器的設計步驟如下：步驟1：將數字濾波器H(z)的技術指標p和s，通過tan(/2)轉變為模擬濾波器G(s)的技術指標p和s，作歸一化處理后，得到p1，ss/p；步驟2：化解為模擬原型濾波器G(s)的技術指標；步驟3：設計模擬原型濾波器G(p)；步驟4：將G(p)轉換為模擬濾波器的轉移函數G(s)；步驟5：將G(s)轉換成數字濾波器的轉移函數H(z)，s=(z1)(z+1)。所謂原型濾波器是指歸一化的低通濾波器。本節主要討論通過IIR數字濾波器的原型轉換設計法和IIR數字濾波器的直接設計方法來設計數字高通、帶通及帶阻濾波器其轉換方法主要有3種：一是直接由模擬低通濾波器轉換成數字高通、帶通或帶阻濾波器；二是先由模擬低通濾波器轉換成模擬高通、帶通或帶阻濾波器，然后再把它轉換成相應的數字濾波器；三是將模擬低通濾波器先轉換成數字低通濾波器，再通過變量代換變換成高通、帶通或帶阻濾波器。以先由模擬低通濾波器轉換成模擬高通、帶通或帶阻濾波器，然后再把它轉換成相應的數字濾波器為例首先由頻率變換將模擬低通原型濾波器變換成模擬高通、帶通或帶阻濾波器，然后根據雙線性Z變換將模擬高通濾波器變換成相應的高通、帶通或帶阻數字濾波器。其模擬低通原型濾波器變換成模擬高通、帶通或帶阻濾波器的關系式如表4-1所示。模擬低通原型濾波器變換成模擬高通、帶通或帶阻濾波器頻率變換模擬濾波器類型頻率變換公式高通帶通帶阻將模擬低通濾波器先轉換成數字低通濾波器，再通過變量代換變換成高通、帶通 或帶阻濾波器該方法中的模擬低通濾波器先轉換成數字低通濾波器，可以利用沖激響應不變法和雙線性變換法來實現，再通過變量代換變換成數字高通、帶通或帶阻濾波器。用雙線性變換法設計一個帶通橢圓濾波器使其幅頻特性逼近于一個具有如下技術指標的模擬帶通濾波器：Wp1=10Hz,Wp2=20Hz,在通帶內的最大衰減為0.5dB,在阻帶內的最大衰減為50dB，抽樣頻率為100Hz。MATLAB實現，該數字濾波器的幅頻特性曲線如圖所示。wp1=10;wp2=20;Fs=100;rp=0.5;rs=50;wp1=2*pi*wp1;wp2=2*pi*wp2;Bw=wp2-wp1;Wo=sqrt(wp2*wp1);z,p,k=ellipap(7,rp,rs);A,B,C,D=zp2ss(z,p,k);At,Bt,Ct,Dt=lp2bp(A,B,C,D,Wo,Bw);At1,Bt1,Ct1,Dt1=bilinear(At,Bt,Ct,Dt,Fs);num,den=ss2tf(At1,Bt1,Ct1,Dt1);H,W=freqz(num,den);plot(W*Fs/(2*pi),abs(H);grid;xlabel(頻率/Hz);ylabel(幅值);圖4-6 橢圓濾波器的頻率響應窗體底端使用雙線性Z變換法設計一個低通數字濾波器，給定的數字濾波器的技術指標為fp100Hz，fs=300Hz，p=3dB，p=20dB，抽樣頻率為Fs=1000Hz。實現例4-2的MATLAB程序如例程4-2所示，程序的運行結果如圖4-2所示基于雙線性Z變換法設計數字濾波器：fp=100;fs=300;Fs=1000;rp=3;rs=20;%wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Fs=Fs/Fs;%Firstly to finish frequency prewarping;wap=tan(wp/2);was=tan(ws/2);n,wn=buttord(wap,was,rp,rs,s)%Note:s!z,p,k=buttap(n);bp,ap=zp2tf(z,p,k)bs,as=lp2lp(bp,ap,wap)%Note:s=(2/Ts)(z-1)(z+1);Ts=1,that is 2fs=1,fs=0.5;bz,az=bilinear(bs,as,Fs/2)h,w=freqz(bz,az,256,Fs*1000);plot(w,abs(h);grid on。IIR濾波器的DSP實現.在用定點DSP器件設計數字濾波器時，一個重要的問題就是由于硬件字長精度有限，運算會出現溢出。IIR濾波器可以用較少的階數獲得很高的選擇特性，所用的存儲單元少、運算次數少，具經濟、高效的特點。在相位要求不敏感的場合，如語音通信等，很適合用IIR濾波器。 簡述CCS環境 CCS，即Code Composer Studio，是TI公司在1999年推出的一個開放、具有強大集成開發環境。它最初是由GO DSP公司為TI的C6000系列DSP開發的。在TI收購了GO DSP后，將CCS擴展到了其它系列。現在所有TI的DSP都可以使用CCS進行開發，但是其中的DSPBIOS功能只有C5000和C6000的CCS中才提供。 以前的DSP軟件開發都是在一個分散的開發環境下進行，程序的編寫、代碼的生成以及調試等都是要通過命令來完成，類似于以前的DOS，十分煩雜。而CCS的出現是DSP開發軟件的一次革命性的變化。CCS主要由代碼生成工具、CCS集成開發環境、DSPBIOS和API函數以及RTDX組成。 CCS主要特點集成可視化代碼編輯界面，可以方便地直接編寫C、匯編、.h文件、.cmd文件等。集成代碼生成工具，包括匯編器、優化的C編譯器和連接器等。具有完整的基本調試工具，可以載入執行文件(.out)，查看寄存器窗口、存儲器窗口和變量窗口、反匯編窗口等，支持在C源代碼級進行調試。支持多片DSP聯合調試。斷點工具，支持硬件斷點、數據空間讀/寫斷點、條件斷點等。探針工具，用于進行算法仿真，數據監視等。剖析工具，用于評估代碼執行的時間。數據圖形顯示工具，可繪制時域/頻域波形、眼圖、星座圖等，并可以自動刷新。提供GEI工具，用戶可以根據需要編寫自己的控制面板/菜單，從而方便直觀地修改變量，配置參數。DSP/BIOS和API函數以及RTDX插件DSP/BIOS(Basic Input Output System)和API(Application Program Interface)函數為CCS的主要插件之一。DSPBIOS可以看作是一個準實時操作系統，支持TI DSP芯片的各種實時操作系統都是以DSP/BIOS作為底層軟件，為嵌入式應用提供基本的運行服務。并且，它還能實時獲取目標機的信息，并將其傳遞給主機上的BIOSCOPE工具，對應用程序進行實時分析RTDX (Real Time Data Exchange)插件是CCS中另一個十分重要的插件。實時數據交換技術為CCS提供了一個實時、連續的可視環境，開發人員可以看到DSP應用程序工作的真實過程。RTDX允許系統開發者在不停止運行目標應用程序的情況下在計算機和DSP芯片之間傳輸數據，同時還可以在主機上利用對象鏈接嵌入(OLE)技術分析和觀察數據。RTDX可以在DSPBIOS中使用，也可以脫離DSP/BIOS使用。由于CCS中的Simulator不支持RTDX，所以必須在連接有硬件仿真器或目標板的Emulator下使用。 CCS配置 擊桌面圖標的“setup CCStudio v3.1”圖標，運行CCS設置程序，如圖所示：點擊AddSave&quit完成設置。用標準配置文件設置系統配置 CCS環境中工程文件的使用 建立工程文件在CCS集成環境下開發匯編程序或者C/C+程序，首先要建立一個工程項目文件（*.pjt），再向工程項目文件中添加匯編程序源文件（*.asm），C/C+源文件（*.c）和鏈接命令文件（*.cmd），并設置工程項目選項。使用CCS開發應用程序的一般步驟如下：創建或打開一個工程項目文件（*.pjt），編輯各類文件，可以使用CCS提供的集成編輯環境，對鏈接命令文件和源程序進行編輯。對工程項目進行編譯。在編譯過程中如果出現語法錯誤，將在編譯鏈接信息視窗（build）窗口中顯示錯誤信息，用戶可以根據顯示的信息找到錯誤的位置，更改錯誤。對結果和數據進行分析和算法評估。用戶可以利用CCS提供的探測點.圖形顯示和性能評價等工具，對運行結果及輸出數據進行分析，評估算法的可能性。 下面分別進行介紹建立工程文件，單擊ProjectNew命令，系統將彈出如圖所示的對話框，在該對話框中輸入項目文件名，如iir單擊“完成”系統就會創建一個名為iir.pjt的工程項目文件。創建新文件選擇FileNew可打開一個新的編輯窗口。在新窗口中輸入源代碼。選擇FileSave，在出現的對話框中輸入一個文件名，并選擇一個擴展名（C源代碼選擇*.c，匯編源代碼選擇*.asm），然后單擊保存將源程序保存。 向工程項目中添加文件添加文件：單擊ProjectAdd File to Project命令然后會彈出如圖所示的對話框，單擊打開完成對文件的添加。向工程項目中添加文件對話框 編譯鏈接和運行目標文件 對程序進行編譯鏈接對程序進行編譯：執行ProjectCompile命令就可以對當前的匯編程序進行編譯生成.obj文件。如果程序存在語法錯誤，那么就會在下面的編譯鏈接信息框中顯示錯誤信息。根據錯誤提示，讀者可對程序進行修改2。對程序進行編譯鏈接執行ProjectBuild命令，就可以對當前的項目文件同時進行編譯，匯編和鏈接操作，并生成與工程項目名稱相同的可執行的.out文件。如果有錯誤信息，則會在“編譯鏈接”信息框中顯示。也可以執行Project下的build All命令，所有項目中的文件重新編譯，匯編和鏈接，生成.out文件。 裝載.out文件單擊主菜單“File”中“Load Program”選項，在對話框中，在CCS安裝目錄下，找到構建該工程的Debug目錄，選擇構建生成后的.out文件，并打開CCS裝載完畢后，該“*.out”文件到目標dsp之后，會自動彈出“Disassembly”窗口如圖所示。CCS調試程序時出現的視窗加工程界面代碼如下：extern void InitC5402(void);extern void OpenMcBSP(void);extern void CloseMcBSP(void);extern void READAD50(void);extern void WRITEAD50(void);/* Main Function Program*/ #include stdio.h #include math.h #def