O引言濾波器從廣義上來說對特定的頻點或頻點以外的頻率進行有效濾波的電路，這種電路保留輸入信號中的有用信息，濾除不需要的信息，從而達到信號的檢測、提取、識別等不同的目的。如果處理的信號是時域離散信號，那么相應的處理系統(tǒng)就稱為數(shù)字濾波器，由于在實際工作中被處理的信號都是幅度量化的數(shù)字信號，因此，數(shù)字濾波器實際上是用有限精度的算法實現(xiàn)一個線性時不變的時域離散系統(tǒng)。目前，數(shù)字濾波器的應用越來越廣泛，它已深入到很多領域，如圖象處理、醫(yī)學生物信息處理、地質信號處理和模式識別處理等。數(shù)字濾波器的種類很多，分類方法也不同，可以從功能上分類，也可以從實現(xiàn)方法上及設計方法上來分類等等。濾波器在功能上總的可分為四類，即低通(LP).高通(HP)、帶通(BP)、帶阻(BS)濾波器等，從實現(xiàn)方法上，由有限長沖激響應的數(shù)字濾波器被稱為FIR濾波器，具有無限長沖激響應的數(shù)字濾波器稱作IIR濾波器。1IIR濾波器的設計1.1雙線性變換法的基本原理本文中用的是雙線性變換法，因此這里簡單介紹一下。雙線性變換法的S域與Z域間的變換關系為：由式(1)可以看出，z和s之間可以直接代換，由于這是一非線性變換，需考察一下是否能把j。映射成單位圓，以及是否能將s域左半平面映射到單位圓內部。對于5平面上的虛軸j。，用s=j。代入式(1)得：1+(￡)沮 nrZ=——|—=/E2)=必1-(專)沮w=2arctan\z\=lw=2arctan可見:上式表明S平面與Z平面一一單值對應，S平面的虛軸(整個j。)經(jīng)映射后確已成為z平面上的單位圓，但。與為非線性關系，因此，通過雙線性變換后兩個濾波器的頻率特性形狀不能保持相同，雙線性變換不存在混迭效應。對于s域的左半平面，用s=o+j。入式(1)，得到：1+（號）（。+用）N= 1-（§）（°+淚）由此得:2E切￥由此得:知）2此式表明，當。<0,有|z|<1,因而s平面的左半平面被映射在單位圓內部，這意味著穩(wěn)定的模擬濾波器經(jīng)雙線性變換可以映射成穩(wěn)定的數(shù)字濾波器。1.2IIR濾波器的設計IIR數(shù)字濾波器的設計借助模擬濾波器原型，再將模擬濾波器轉換成數(shù)字濾波器，這些過程已經(jīng)成為一整套成熟的設計程序。模擬濾波器的設計已經(jīng)有了一套相當成熟的方法，它不但有完整的公式，而且還有較為完整的圖表查詢，因此，充分利用這些已有的資源將會給數(shù)字濾波器的設計帶來很大的方便。已知數(shù)字低通濾波器的設計要求wp，ws，Rp和As，首先設計一個等效的模擬濾波器，然后再將它映射為所期望的數(shù)字濾波器來確定H（z）。對這個過程所要求的步驟是：1）選取T并確定模擬頻率： 利用設計參數(shù)Op，Os，Rp和As，設計一個模擬濾波器Ha（s）；這可以利用模擬濾波器的原型來完成；3）再將⑴代入Ha（s），求出H（z）。本文給定數(shù)字濾波器特性為：通帶：0.99VIH（必）IVI.010<w<0.4兀阻帶：|H（次）IV0.0010.67<w<7r1） 根據(jù)題意有，凹=0,4兀，飽=0.6江;計算時假定T=lo當w=0.4兀，劣=#tan岑=2tan們尹=1.4530w=0.67t,ns=-^7tan號=2tan^^=2.75282） 設通帶內起伏：1一bi《|H（尸）|W1+?！?則。1=0.01阻帶內起伏：IH（舟），則。2=0?001..?衰減度為（dB）：月,=-20也與斜=0.1737阻帶內衰減度為（dB）:A=-201g三奇=60對于幅度衰減指標，模擬和數(shù)字是一樣的。綜上（1）、（2），模擬濾波器特性如下:

月＞=0.1737,As=60

Cp=l.4530,Q=2.7528現(xiàn)在可以開始設計模擬濾波器，進而確定數(shù)字濾波器。下面給出Butterth濾波器逼近的方法。首先確定階數(shù)N和截止角頻率，根據(jù)指標要求列出：11^00.022)1=z11^00.022)1=z/-L,/0.0221 0.033\2NQcI1/4()0.033)1=-=J1+2N=0.9802=0.0009=14解得:N=「14（10%加-1）/（1（）A加-1）]=14解得:N=2ig（np/a）現(xiàn)要使通帶指標富裕，則Q=——八\2Ny（10Vl0-l）1.6806由此得到巴特沃茲濾波器極點分布為:50=513*=1.6806』28=-0.1120+0.99371$i=si2*=1.6806質=-0.3303+0.94391$2=A”*=1.6806#塘=-0.5320+0.8467i、3=si0*=1?6806必=-0.7071+0.7071154=59*=1.6806播=-0.8467+0.53201$5=$8*=1.6806』繆=—0.9439+0.3303iS6=S7*=1.6806?28=-0.9937+0.11201模擬濾波器的轉移函數(shù)為：H(S)=「4十8.9314/3+39.8850s"+117.7337s"+ 1 256.1214s】°+433.7284s’+589.02067+651.2664s"71 +589.0206s6+433.7284s5+256.1214s1+117.7337^ 1 +39.8850s2+8.1314s+1把、=2(岳M)代入上式得()_0.0001+0.0011zT+0.0071z-2+0.0284/3+HE)--1一i.5395zt+2.9473z一2-2.8363z-3+-0.0782z'4+0.1563g-5+0.2345z-6+2.7482z一4一i.7703/5+1.0616z鄧一0.4647z瑚0.2680/7+0.2345/8+0.1563/9+°to++0.1794z_8-0.0516z~9+0.0123z~100.O284zT+0.OO7"T2+0001*73+0000/14_0.002/11+0.0003/122IIR濾波器的MATLAB實現(xiàn)模擬型濾波器的設計過程，步驟繁瑣，計算量大。在MATLAB中，我們可以直接根據(jù)所給定的數(shù)字濾波器指標調用特定的函數(shù)來確定濾波器的系統(tǒng)函數(shù)。具體實現(xiàn)過程如下（程序的主要部分，且只返回數(shù)字濾波器的階數(shù)和截止頻率）：為調用階數(shù)選擇函數(shù)，下面的指標已經(jīng)把Wp和Ws對n進行了歸一化，且函數(shù)中返回的頻率也是歸一化頻率。殆=0.1737/$=60P %濾波器指標，W?=1.4530,"=2.7528%實現(xiàn)ButterWorth數(shù)字濾波器[Nbutt,Wcbutt]=buttord(Wp,Ws,Rp,As);%確定ButterWorth的階數(shù)和截止頻率[Bbutt,Abutt]=butter(Nbutt,Wcbutt);%確定Butterworth濾波器的系統(tǒng)函數(shù)[Hbutt,Wbutt]=freqz(Bbutt,Abutt,501);%計算頻率響應同理可實現(xiàn)ChebyshevI、ChebyshevII、ellipSe數(shù)字濾波器。Butterworth型和ellipke型系統(tǒng)的幅頻和相頻特性如圖l；ChebyshevI與ChebyshevII系統(tǒng)的幅頻和相頻特性如圖2（頻率為對n歸一化頻率）。頻率/弧度頻率/弧度頻率/弧度頻率/弧度圖1Butterworth和ellipke系統(tǒng)的幅、相頻特性1|—0.8—制0.6-0.4—0.2—°0頻率/弧度 頻率/孤度圖2ChebyshevI與Chebyshev口系統(tǒng)的幅、相頻特性3結論從圖中我們可以看出，ButterWorth濾波器的頻率特性是平坦的，沒有出現(xiàn)波動，即幅度隨著頻率單調遞減，它的相頻特性接近于分段特性。而橢圓濾波器幅頻特性存在著波紋（圖中由于坐標比例的關系，只能畫出通帶波紋，阻帶波紋則無法顯示）；同時，它不具備線性相位特性。兩種契比雪夫濾波器的頻率特性如圖所示。從中可知道，這兩種濾波器的幅頻特性都有著波動（其中，由于CII型濾波器是在阻帶中出現(xiàn)等波動，因此無法表示出來），它們不具備線性相