無限長(zhǎng)單位脈沖響應(yīng)濾波器的設(shè)計(jì)方法PPT教案

1、無限長(zhǎng)單位脈沖響應(yīng)濾波器的設(shè)計(jì)方法無限長(zhǎng)單位脈沖響應(yīng)濾波器的設(shè)計(jì)方法 下面舉例討論應(yīng)用模擬濾波器低通原型，設(shè)計(jì)各種數(shù)字濾波器的基本原理，著重討論雙線性變換法。 一低通變換 通過模擬原型設(shè)計(jì)數(shù)字濾波器的四個(gè)步驟： 1）確定數(shù)字濾波器的性能要求，確定各臨界頻率k。2）由變換關(guān)系將k映射到模擬域，得出模擬濾波器的臨界頻率值k。3）根據(jù)k設(shè)計(jì)模擬濾波器的Ha(s)4） 把Ha(s) 變換成 H(z)（數(shù)字濾波器傳遞函數(shù)）第1頁(yè)/共62頁(yè) 設(shè)采樣周期 ,設(shè)計(jì)一個(gè)三階巴特沃茲LP濾波器,其3dB截止頻率fc=1khz。分別用脈沖響應(yīng)不變法和雙線性變換法求解。 解：a. 脈沖響應(yīng)不變法 由于脈沖響不變法的頻

2、率關(guān)系是線性的，所以可直接按c =2fc設(shè)計(jì)Ha(s)。根據(jù)上節(jié)的討論，以截止頻率c 歸一化的三階巴特沃茲 濾波器的傳遞函數(shù)為： 以 代替其歸一化頻率，得：322211)(ssssHa32)/()/( 2)/( 211)(cccassssHcs /)4(250khzfsTs第2頁(yè)/共62頁(yè))(sHaccf2cs /第3頁(yè)/共62頁(yè) 為進(jìn)行脈沖響應(yīng)不變法變換，計(jì)算Ha(S)分母多項(xiàng)式的根，將上式寫成部分分式結(jié)構(gòu)： 對(duì)照前面學(xué)過的脈沖響應(yīng)不變法中的部分分式形式 有 將上式部分系數(shù)代入數(shù)字濾波器的傳遞函數(shù)： , -極點(diǎn) 2/ )31 (3/2/ )31 (3/)(6/6/jcsecjcseccscs

3、Hajj6/2113/;,jccecAsA2/)31 (,3/;2/)31 (36/32jsecAjscjcNiTSiZeAZHi111)(iS第4頁(yè)/共62頁(yè)并將 代入，得： 合并上式后兩項(xiàng)，并將 代入，計(jì)算得： 12/ )31(6/12/ )31(6/11)3/(1)3/(1/)(ZeeTZeeTZeTZHjjcjjcCccc5 . 02Tfcc21112079.01905.015541.0571.12079.01571.11)(ZZZZTZHTcc/第5頁(yè)/共62頁(yè)與的數(shù)字濾波器具有相同的傳遞函數(shù)，這一結(jié)論適合于所有的數(shù)字濾波器設(shè)計(jì)。最后得：scff /21112079. 01905.

4、015541. 0571. 12079. 01571. 1)(zzzzZHCKHzfKHzfcs10,40KHzfKHzfcs1,4第6頁(yè)/共62頁(yè) b. 雙線性變換法 （一）首先確定數(shù)字域臨界頻率 （二）根據(jù)頻率的非線性關(guān)系，確定預(yù)畸的模擬濾波器臨界頻率 (三 ) 以 代入歸一化的三階巴特沃模擬器傳遞函數(shù) 并將 代入上式。 （四）將雙線性變換關(guān)系代入，求H(Z)。5 . 02TfccTtgTcc222cs/32)/()/(2)/(211)(cccassssHTc/2第7頁(yè)/共62頁(yè)311211111121111211211)()(11zzzzzzsHZHzzTsa 2311123121111

5、21131311131313111113131311212113131312111312122122211111141111111211111211211zzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzz第8頁(yè)/共62頁(yè) 圖1 三階Butterworth 數(shù)字濾波器的頻響脈沖響應(yīng)不變法雙線性變換法fs/2第9頁(yè)/共62頁(yè)我們也可以用 MATLAB 完成設(shè)計(jì)，在 MATLAB 中相關(guān)的語(yǔ)句有butter（巴特沃茲濾波器） 、impinvar（脈沖響應(yīng)不變法） 、bilinear（雙線性變換） ，具體的程序如下：B,A=butter(3,2*pi*1000,s); num1,de

6、n1=impinvar(B,A,4000); h1,w=freqz(num1,den1); B,A=butter(3,2/0.00025,s); num2,den2=bilinear(B,A,4000); h2,w=freqz(num2,den2); f=w/pi*2000; plot(f,abs(h1),-.,f,abs(h2),-); grid; xlabel(頻率/Hz ) ylabel(幅值/dB)第10頁(yè)/共62頁(yè)020040060080010001200140016001800200000.10.20.30.40.50.60.70.80.91頻率/Hz圖3.14 三階巴特沃茲濾波器

7、的頻率響應(yīng)幅值第11頁(yè)/共62頁(yè) 圖1為兩種設(shè)計(jì)方法所得到的頻響，對(duì)于雙線性變換法，由于頻率的非線性變換，使截止區(qū)的衰減越來越快，最后在折 疊頻率處 形成一個(gè)三階傳輸零點(diǎn),這個(gè)三階零點(diǎn)正是模擬濾波器在 處的三階傳輸零點(diǎn)通過映射形成的。因此,雙線性變換法使過渡帶變窄,對(duì)頻率的選擇性改善,而脈沖響應(yīng)不變法存在混淆,且沒有傳輸零點(diǎn)。, 1Z第12頁(yè)/共62頁(yè) 二.高通變換 設(shè)計(jì)高通、帶通、帶阻等數(shù)字濾波器時(shí)，有兩種方法： 先設(shè)計(jì)一個(gè)相應(yīng)的高通、帶通或帶阻模擬濾波器，然后通過脈沖響應(yīng)不變法或雙線性變換法轉(zhuǎn)換為數(shù)字濾波器。 模擬原型 模擬高通、帶通、帶阻 數(shù)字高通、帶通、帶阻 設(shè)計(jì)方法同上面討論的低通濾

8、波器的設(shè)計(jì)。 即確定 轉(zhuǎn)換為相應(yīng)的 高通、帶通、帶阻 模擬濾波器的設(shè)計(jì) Ha(s) H(Z) 直接利用模擬濾波器的低通原型，通過一定的頻率變換關(guān)系，一步完成各種數(shù)字濾波器的設(shè)計(jì)。 頻率變換 模擬原型 數(shù)字低通、高通、帶通、帶阻kk第13頁(yè)/共62頁(yè)這里只討論第二種方法。因其簡(jiǎn)捷便利，所以得到普遍采用。 變換方法的選用： 脈沖響應(yīng)不變法：對(duì)于高通、帶阻等都不能直接采用，或只 能在加了保護(hù)濾波器后才可使用。因此，使 用直接頻率變換（第二種方法），對(duì)脈沖響 應(yīng)不變法要有許多特殊的考慮，它一般應(yīng)用 于第一種方法中。 雙線性變換法：下面的討論均用此方法，實(shí)際使用中多數(shù)情況 也是如此。基于雙線性變換法的高

9、通濾波器設(shè)計(jì)： 在模擬濾波器的高通設(shè)計(jì)中，低通至高通的變換就是S變量的倒置，這一關(guān)系同樣可應(yīng)用于雙線性變換，只要將變換式中的S代之以1/S，就可得到數(shù)字高通濾波器. 即11112zzTs第14頁(yè)/共62頁(yè)jjjctgTeeTseZjjj22112,時(shí)22ctgT如圖第15頁(yè)/共62頁(yè) 映射到 即 映射到 即 圖1 高通變換頻率關(guān)系 這一曲線的形狀與雙線性變換時(shí)的頻率非線性關(guān)系曲線相對(duì)應(yīng)，只是將 坐標(biāo)倒置，因而通過這一變換后可直接將模擬低通變?yōu)閿?shù)字高通,如圖2。22ctgT01z1z01.01.00第16頁(yè)/共62頁(yè) 圖2 高通原型變換 第17頁(yè)/共62頁(yè)應(yīng)當(dāng)明確： 所謂高通DF，并不是高到 ，

10、由于數(shù)字頻域存在 折疊頻 率 ，對(duì)于實(shí)數(shù)響應(yīng)的數(shù)字濾波器， 部分只是 的鏡象部分，因此有效的數(shù)字域僅是 ，高通也僅指這一段的高端，即到 為止的部分。 高通變換的計(jì)算步驟和低通變換一樣。但在確定模擬原型預(yù)畸的臨界頻率時(shí)，應(yīng)采用 ,不必加負(fù) 號(hào)，因臨界頻率只有大小的意義而無正負(fù)的意義。2由0由022kkctgT第18頁(yè)/共62頁(yè) ： 采樣 設(shè)計(jì)一個(gè)三階切比雪夫高通DF，其通過頻率 （但不必考慮 以上的頻率分量），通帶內(nèi)損耗不大于1dB。 解：首先確定數(shù)字域截止頻率 , 則 切比雪夫低通原型的模函數(shù)為： 為N階切比雪夫多項(xiàng)式,100,10usTkHzfskHzf5 .2kHzfs525.0211Tf

11、22211TctgT1222/11)(NaVjH)(NV第19頁(yè)/共62頁(yè)5089. 01101 . 03212131319883. 0238. 14913. 04913. 0)(ssssHadB1為方便，將 和 S 用T/2歸一化， 則12, 12/11TssT329883. 0238. 14913. 04913. 0)(ssssHa第20頁(yè)/共62頁(yè)于是 321321112041. 06043. 03432. 013311321. 0)()(11zzzzzzsHZHzzsa圖3 三階切比雪夫高通頻響第21頁(yè)/共62頁(yè)例5（書上 ） 設(shè)計(jì)一數(shù)字高通濾波器，它的通帶為400500Hz，通帶內(nèi)容

12、許有0.5dB的波動(dòng)，阻帶內(nèi)衰減在小于317Hz的頻帶內(nèi)至少為19dB，采樣頻率為1,000Hz。確定最小階數(shù) N。 模擬切比雪夫?yàn)V波器設(shè)計(jì)中階數(shù)的確定公式為 求得最小的N： 1220184. 0110,5 . 01 . 02今)(1lg102dBeA9 . 1210A及、2AcScsAN/cosh/1cosh121第22頁(yè)/共62頁(yè)h,w=freqz(num,den);f=w/pi*500;plot(f,20*log10(abs(h);axis(0,500,-80,10);grid;xlabel()ylabel(幅度/dB)第23頁(yè)/共62頁(yè)0501001502002503003504004

13、50500-80-70-60-50-40-30-20-10010頻率/Hz幅度/dB頻率/Hz 切比雪夫高通濾波器幅度/dB第24頁(yè)/共62頁(yè)三帶通變換 如圖1 ,如果數(shù)字頻域上帶通的中心頻率為 ，則帶通變換的目的是將: （頻率映射關(guān)系具有周期性, 幅頻響應(yīng)具有原點(diǎn)對(duì)稱性）。 即將S的原點(diǎn)映射到 ，而將 點(diǎn)映射到 ，滿足這一要求的雙線性變換為：000000映射0jezjS1z11cos2)1(122zzzzzezezsojjoo模擬低通000映射第25頁(yè)/共62頁(yè) 圖1 帶通原型變換 第26頁(yè)/共62頁(yè)jez 當(dāng) 時(shí) 因此 （帶通變換關(guān)系 ）jjojjjojjeeeeeeescos211cos

14、222jsjso又,sincoscossincoscoso第27頁(yè)/共62頁(yè) 圖中 點(diǎn)正好映射在 上，而 映射在 ， 兩端，因此滿足帶通變換的要求。000帶通變換的頻率關(guān)系第28頁(yè)/共62頁(yè)穩(wěn)定性證明： 同時(shí)，這一變換也滿足穩(wěn)定性要求，設(shè) 由于上式完全是實(shí)數(shù)，所以是映射在S平面 軸上。其中分子永遠(yuǎn)非負(fù)的 ， 因此 的正負(fù)決定于分母由此證明了，S左半平面映射在單位圓內(nèi)，而右半平面映射在單位圓外，這種變換關(guān)系是穩(wěn)定的變換關(guān)系，可用它來完成帶通的變換，如圖1。0 rz11cos222rrrso1cos1211cos212222rrrrrroo0)cos1(212orr12r010,1時(shí)，時(shí)，rr第2

15、9頁(yè)/共62頁(yè)設(shè)計(jì)： 設(shè)計(jì)帶通時(shí)，一般只給出上、下邊帶的截止頻率 作為設(shè)計(jì)要求。 為了應(yīng)用以上變換，首先要將上下邊帶參數(shù) 換算成中心頻率 及模擬低通截止頻率 。為此將 代入變換關(guān)系式：由于 在模擬低通中是一對(duì)鏡象頻率，代入上面兩等式，求出21,21,0cco,21求21,111sincoscoso21,210cos2cos2cossinsinsincos212121210222sincoscoso第30頁(yè)/共62頁(yè)又 同時(shí)也就是模擬低通的截止頻率 ，有了這兩個(gè)參數(shù)就可完成全部計(jì)算。 ：采樣 fs=400kHz，設(shè)計(jì)一巴特沃茲帶通濾波器，其3dB邊界頻率分別為f2=90kHz,f1=110kHz

16、,在阻帶f3=120kHz處最小衰減大于10dB。解：確定數(shù)字頻域的上下邊帶的角頻率求中心頻率： 1c11sincoscosoc45.0/222sff55.0/211sff6.0/233sff55. 0sin45. 0sin55. 045. 0sincos05 . 00第31頁(yè)/共62頁(yè)求模擬低通的通帶截止頻率 與阻帶邊界頻率 ： 從 頻率增加了約1.05倍，衰減增加了（10-3）dB，故選用二階巴特沃茲濾波器可滿足指標(biāo)（查表） 歸一化的系統(tǒng)函數(shù)： 代入 ， 代入變換公式 c1584. 055. 0sin55. 0cos5 . 0cosc3249.06 .0sin6 .0cos5 .0cosr

17、rc 到121)(2sssHa1/2/1)(2ccasssH1111cos222202zzzzzs第32頁(yè)/共62頁(yè)f1f2f3fdB10dB30)(jeH例6 帶通濾波器設(shè)計(jì)第33頁(yè)/共62頁(yè)111317. 6211317. 61)()(2222221122zzzzsHazHzzs66.3725.8466.37) 1(2422zzz第34頁(yè)/共62頁(yè) 四帶阻變換 把帶通的頻率關(guān)系倒置就得到帶阻變換。 ozzzscoscossin,1cos21022011coscossincco,21求給定21210sinsin)sin(cos第35頁(yè)/共62頁(yè)grid;xlabel(頻率/Hz)ylabel

18、(幅度/dB)第36頁(yè)/共62頁(yè)5060708090100110120130140150-30-25-20-15-10-50510頻 率 / Hz幅 度 / dB頻率/Hz巴特沃茲帶阻濾波器幅度/dB第37頁(yè)/共62頁(yè) 上一節(jié)討論了由模擬網(wǎng)絡(luò)的低通原型來設(shè)計(jì)各種DF的方法，這種原型變換的設(shè)計(jì)方法同樣也可直接在數(shù)字域上進(jìn)行。 DF低通原型函數(shù) 這種變換是由 所在的Z平面到H（z）所在的Z平面的一個(gè)映射變換。 為便于區(qū)分變換前后兩個(gè)不同的Z平面，我們把變換前的 Z平面定義為u平面，并將這一映射關(guān)系用一個(gè)函數(shù)g表示： 變換)(zHp)(ZHp)(11zgu各種DF的 H（z）第38頁(yè)/共62頁(yè)平面平

19、面zuZHuHzgup)()()(11于是，DF的原型變換可表為： )(11)()(zgupuHzH第39頁(yè)/共62頁(yè) 函數(shù) 的特性： 1) 是 的有理函數(shù)。 2）希望變換以后的傳遞函數(shù)保持穩(wěn)定性不變，因此要求 u的單位圓內(nèi)部必須對(duì)應(yīng)于z的單位圓內(nèi)部。 3） 必須是全通函數(shù)。 為使兩個(gè)函數(shù)的頻響滿足一定的變換要求，Z的單位圓應(yīng)映射到u的單位圓上，若以 分別表示u平面和Z平面的單位圓，則式為 且必有 ，其中 是 的相位函數(shù)， 即函數(shù)在單位圓上的幅度必須恒為1，稱為全通函數(shù)。 )(1zg)(1zg1z)(1zgjjee 和)(jjjjeegege )(jeg1jeg第40頁(yè)/共62頁(yè) 任何全通函數(shù)

20、都可以表示為： 其中 為極點(diǎn)，可為實(shí)數(shù)，也可為共軛復(fù)數(shù)，但必須在單位圓以內(nèi)，即 ，以保證變換的穩(wěn)定性不變，*為取共軛。 的所有零點(diǎn) 都是其極點(diǎn)的共軛倒數(shù) N：全通函數(shù)的階數(shù)。 變化時(shí)，相位函數(shù) 的變化量為 。 不同的N和 對(duì)應(yīng) 各類不同的變換。 Niiizzzg11*111)(i1i)(1zg*/1i0i N第41頁(yè)/共62頁(yè)g(-1) = -1滿足以上要求的映射函數(shù)應(yīng)為：其中是實(shí)數(shù)，且)(jpeHjeH0 N1111)(zzzg10第42頁(yè)/共62頁(yè)圖1(a) LP-LP變換（有對(duì)稱性）c00第43頁(yè)/共62頁(yè) 代入（1）式，可得到上述變換所反映的頻率變換關(guān)系： 由此得 上式把 ， 。頻率特

21、性: 呈線性關(guān)系，其余為非線性。 當(dāng) 時(shí)， ， 帶寬變窄， 當(dāng) 時(shí)， ， 帶寬變寬， 適當(dāng)選擇 ，可使 變換為 ，如上圖所示 。 :低通原型截止頻率， : 變換后截止 頻率 jjeuez及將）（21jjjeeecos12sin122arctg,0zucc時(shí)00cccccccc00第44頁(yè)/共62頁(yè)圖 LP-LP頻率變換特性 cc第45頁(yè)/共62頁(yè) 確定 ： 把變換關(guān)系 帶入（2）式 ，有： 得 （2） 式的 頻率關(guān)系，如圖 cccjjjeee12sin2sincccccc第46頁(yè)/共62頁(yè) LP-HP a .基本思想：上述 LP 變換中的Z代以Z , 則 LP = HP 。)( zH 第47頁(yè)

22、/共62頁(yè)）（21jjjeee或, 00LP-HP變換把cc如圖2（a）， 在上述LP-LP 變換中，將 Z代以Z , 得 LP - HP變換關(guān)系：）（11111111zzzzu第48頁(yè)/共62頁(yè) 原型低通的截止頻率 對(duì)應(yīng)于高通的邊界頻率 ，欲將 變換到 ，由（2）式， 有： cccccccccccjjjjjjeeeeee1,12cos2coscccc:確定第49頁(yè)/共62頁(yè)圖2 (a) LP Hp變換cc00第50頁(yè)/共62頁(yè) LP-BP LP-BP變換把帶通的中心頻率 故 N=2。 由以上分析得變換關(guān)系： 或 00c2c1) 1 (1)(1122211211zrzrrzrzzgu) 2(1122212jjjjjererreree00000,0時(shí)，如圖3(