《無(wú)限脈沖響應(yīng)濾波器》由會(huì)員分享,可在線閱讀,更多相關(guān)《無(wú)限脈沖響應(yīng)濾波器(37頁(yè)珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、第 10 章 無(wú)限脈沖響應(yīng)濾波器,10.1 無(wú)限脈沖響應(yīng)濾波器基礎(chǔ) 10.2 低通模擬濾波器 10.3 雙線形變換 10.4 巴特沃斯濾波器設(shè)計(jì) 10.5 切比雪夫型濾波器設(shè)計(jì) 10.8 帶通、高通和帶阻 IIR 濾波器,返回,,,infinite impulse response filter(IIR)無(wú)限脈沖響應(yīng)濾波器 bilinear transformation 雙線性變換 prewarping equation 預(yù)扭曲方程 Butterworth filter 巴特沃斯濾波器 Chebyshev Type I filter 切比雪夫I 型濾波器 Chebyshev Type II
2、filter 切比雪夫II 型濾波器 elliptic filter 橢圓濾波器 Impulse invariance method 脈沖響應(yīng)不變法,10.1 無(wú)限脈沖響應(yīng)濾波器基礎(chǔ),濾波器新的輸出,和過(guò)去的輸出及過(guò)去的輸入和現(xiàn)在 輸入有關(guān)。,差分方程 akyn-k= bkxn-k,N k=0,M k=0,若 a0=1,則 yn=-a1yn-1 - - aNyn-N + b0 xn+b1n-1++bMxn-M,傳輸函數(shù),H(z)=,b0+b1z-1+ b2z-2++bMz-M 1+ a1z-1+a2z-2++aNz-N,,這種遞歸濾波器的極點(diǎn)由分母多項(xiàng)式確定,這就可能使濾波器不
3、穩(wěn)定,同時(shí)遞歸濾波器很難實(shí)現(xiàn)線性相位,也就是遞歸濾波器使相位失真,但實(shí)現(xiàn)某種性能要求時(shí)比非遞歸濾波器所需要的系數(shù)少。 設(shè)計(jì)遞歸濾波器的方法是選擇具有待求特性的原型模擬濾波器,然后將其轉(zhuǎn)換為數(shù)字濾波器。,返回,,10.2 低通模擬濾波器,在連續(xù)域中,濾波器用 S 而不用 Z 來(lái)描述。 簡(jiǎn)單低通模擬濾波器的傳輸函數(shù)為 H(s)=,1 s+1,,頻率響應(yīng)為 H()= s用 j代替,1 j +1,,幅度 |H()|=,1 2 +1,,, |H()|0 |H()|1,具有低通特性,,=1弧度/秒 |H()|= 20log( )=-3dB,1 2,,,即 f=/
4、2=1/2 Hz 為它的帶寬 為使濾波器更通用,傳輸函數(shù)變?yōu)?H(s)=,p1 s+p1,,H()= =,p1 j+p1,1 1+j/p1,,,其幅度|H()|=,1 (/p1)2+1,,,=p1 時(shí),出現(xiàn) 3dB,圖10.3,通常模擬濾波器類(lèi)型包括:巴特沃斯、切比雪夫型、 切比雪夫型和橢圓濾波器等。,圖10.4,返回,,10.3 雙線性變換,它為模擬濾波器和數(shù)字濾波器之間的轉(zhuǎn)換提供了一種 方法。 s 2fs,,z 1 z+1,,z用 ej 代替可得到離散時(shí)間傅立葉變換,逆雙線性變換 2tan-1 (/2fs) 圖 10.5 給出了數(shù)字頻率 和模擬頻率 之間的關(guān)系。,,雙線性變
5、換建立了數(shù)學(xué)域和模擬域之間的另一種聯(lián)系: 數(shù)字濾波器的穩(wěn)定區(qū)域是 Z 平面單位圓內(nèi),模擬濾波器 的穩(wěn)定區(qū)域在 S 復(fù)平面的左半部,圖 10.6,,例1:一階模擬低通濾波器的傳輸函數(shù) H(s)=wp1/(s+wp1)。濾波器的-3dB頻率是2000弧度/秒。求出與此模擬濾波器相對(duì)應(yīng)的數(shù)字濾波器的傳輸函數(shù)H(z),采樣頻率為1500Hz.,預(yù)扭曲的數(shù)字濾波器幅度響應(yīng)如圖 10.10 所示,可 看出預(yù)扭曲的濾波器的截止頻率非常接近所要求的值 318.3 Hz。,圖 10.10,返回,,10.4 巴特沃斯濾波器設(shè)計(jì),巴特沃斯濾波器是 IIR 無(wú)限脈沖響應(yīng)濾波器中最簡(jiǎn)單的。,它的傳輸函數(shù) H(s)=,
6、p1 S+p1,,|H()|=,1 (/p1)2+1,,,一階模擬巴特沃斯濾波器,n階模擬巴特沃斯濾波器具有負(fù)責(zé)的傳輸函數(shù) H(s),,其濾波器形狀|H()|= ,與一階相似。,1 (/p1)2n+1,,,階數(shù)越高,滾降越陡,圖10.12 顯示 n=1,2,3,4時(shí)模擬 巴特沃斯濾波器形狀。,圖10.12,描述數(shù)字 IIR 濾波器的術(shù)語(yǔ):(由圖10.13(a)、(b)看),圖10.13,p 通帶波紋:通帶內(nèi)最大和最小增益之差 p 取 3dB (巴特沃斯濾波器) 1 p : 通帶邊緣增益(對(duì)應(yīng)的頻率fp1) (或20log(1 p)dB) s :阻帶波紋,對(duì)應(yīng)的阻帶邊緣增益,
7、 對(duì)應(yīng)頻率 fs1 阻帶邊緣頻率。 阻帶邊緣的增益為20logsdB,設(shè)計(jì)時(shí)阻帶邊緣頻率和阻帶衰減(-20logs)或阻 帶增益(20logs)要給出。圖 10.14 說(shuō)明非線性相 位響應(yīng),表明濾波后的信號(hào)在某種程度上有相位失真。 對(duì)于 n 階巴特沃斯濾波器,滿足通帶和阻帶設(shè)計(jì)要求 所需的階數(shù) n,s1 : 模擬濾波器預(yù)扭曲阻帶邊緣頻率。 p1 : 模擬濾波器預(yù)扭曲通帶邊緣頻率。,設(shè)計(jì)低通巴特沃斯濾波器步驟: 1) 確定待求通帶邊緣頻率 fp1 Hz 、待求阻帶邊緣頻率 fs1 Hz 和待求阻帶衰減 - 20logsdB(或待求阻帶增 益 20logsdB)。通帶邊緣頻率必須對(duì)應(yīng) 3dB
8、增益。,4) 由已給定的阻帶衰減 - 20logs(或增益- 20logs) 確定阻帶邊緣增益 s 。,3) 計(jì)算預(yù)扭曲模擬頻率以避免雙線性變化帶來(lái)的失真。 由 =2fs tan(/2) 求得 p1和 s1,單位是弧度/秒。,2) 用式 =2f/fs 把由 Hz 表示的待求邊緣頻率轉(zhuǎn)成 由弧度表示的數(shù)字頻率,得到 p1 和s1 。,5) 用下式計(jì)算所需濾波器的階數(shù):,n 取整數(shù),例 10.6 設(shè)計(jì)具有巴特沃斯特性的低通 IIR 濾波器, - 3 dB 頻率為 1 200 Hz。在 1 500 Hz 處 增益降到 25 dB。采樣速率為 8 000 Hz。 選擇合適的濾波
9、器階數(shù)并畫(huà)出濾波器形狀。,解: 1) fp1=1200Hz, fs1=1500HZ -20logs=25, fs=8000Hz 2) p1= 2 fp1/fs=0.3 弧度 s1= 2 fs1/fs=0.375 弧度 3) p1=2fs tan(p1/2)=8152.4 弧度/秒 p1= 2fs tan(s1/2)=10690.9弧度/秒,4) 20logs =-25 s=10-25/20=0.0562,5) =10.6 取 n=11,6) 模擬濾波器形狀,1 (/p1)2n+1,,,|H()|= =,1 (/8152.4)22+1,,,|H()|=,1
10、(16000tan(/2)/8152.4)22+1,,,返回,,10.5 切比雪夫型濾波器設(shè)計(jì),切比雪夫?yàn)V波器傳輸函數(shù)非常復(fù)雜,但濾波器形狀 表達(dá)式非常簡(jiǎn)單,n 階切比雪夫型濾波器的形狀定義為,|H()|=,1 1+2Cn2(/p1),,,Cn(x)=,cos(ncos-1(x)) 對(duì)于|x|1 cosh(ncosh-1(x)) 對(duì)于|x|1,,參數(shù) 取決與通帶波紋,濾波器階數(shù)大于 1,切比雪夫?yàn)V波器具有比巴特沃夫 濾波器更陡峭的滾降特性。,圖 10.17,切比雪夫?yàn)V波器的參數(shù)示于圖 10.18,圖 10.18,切比雪夫?yàn)V波器的階數(shù)可由,n,cosh-1(/) cosh-1(s1/p1 ),
11、,= (1/s2) 1,,低通切比雪夫型濾波器的設(shè)計(jì)步驟:,3)對(duì)數(shù)字頻率采用預(yù)扭曲以避免雙線性變換引起的誤差。 由 =2fs tan(/2) 得到p1和 s1,單位是弧度/秒。,2)用公式 =2f/fs 將待求的邊緣頻率轉(zhuǎn)換為數(shù)字頻率 (用弧度表示),得到 p1 和 s1 。,1)確定待求的通帶與阻帶邊緣頻率 fp1 和fs1 、待求的 通帶邊緣增益 20log(1- p) 和待求的阻帶衰減 -20logs(或待求的阻帶增益 20logs )。,5)由指定的衰減-20logs(或增益 20logs),確定 阻帶邊緣增益 s 。計(jì)算式如下:,例 10.9 對(duì)于采樣頻率為 20 KHz
12、 的系統(tǒng),設(shè)計(jì)具有切 比雪夫型特性的 IIR 濾波器。通帶最大增益 為 0 dB,通帶邊緣在 5 kHz 處,其增益為 -1 dB。阻帶邊緣在 7.5 kHz處,其增益為 - 32 dB。,解: fp1=5000 Hz , fs1=7500Hz , 20log(1- p) =-1 fs =20000Hz , 20logs= - 32,p1 = 2fp1/fs =0.5 弧度 s1 = 2fs1/fs =0.75 弧度,p1 = 2fs tan(p1/2) =40000弧度/秒 s1 = 2fs tan(s1/2) =96568.5弧度/秒,7),|H(
13、)|=,1 1+2Cn2(/p1),,,=,1 1+0.2589C42(/40000),,,=2fs tan(/2) 代入上式,|H()|=,1 1+ 0.2589C42(/2),,,C4(x)=,cos(4cos-1(x)) 對(duì)于|x|1 cosh4(4cosh-1(x)) 對(duì)于|x|1,,返回,,10.8 帶通、高通和帶阻 IIR 濾波器,帶通和高通濾波器的階數(shù)要通過(guò)對(duì)低通濾波器的原 型的計(jì)算來(lái)選擇。 圖 10.22 表明了設(shè)計(jì)步驟,圖 10.22,傳輸函數(shù),高通 HH(s)=HL(pp/s) p 為低通濾波器的截止頻率 p為高通濾波器的截止頻率,帶通 HBP(s)=HL(p
14、(s2+lu)/s(u - l)) l:帶通濾波器的低端截止頻率 l u:帶通濾波器的高端截止頻率 u,帶阻 HBS(s)= HL(ps(u- l)/(s2+ul)),從上的p ,p ,l ,u 都要進(jìn)行預(yù)扭曲計(jì)算,避免濾波器失真。 =2fs tan(/2),* 帶通和帶阻濾波器的階數(shù)都為其低通原型的兩倍。,,IIR濾波技術(shù)的一個(gè)應(yīng)用就是產(chǎn)生和恢復(fù)用于按鍵電話機(jī)的雙音多頻信號(hào)(DTMF),FIGURE 10-25 Touch-Tone keypad.,Joyce Van de VegteFundamentals of Digital Signal Processing,Copyri
15、ght 2002 by Pearson Education, Inc.Upper Saddle River, New Jersey 07458All rights reserved.,FIGURE 10-26 Tone generation filter.,Joyce Van de VegteFundamentals of Digital Signal Processing,Copyright 2002 by Pearson Education, Inc.Upper Saddle River, New Jersey 07458All rights reserved.,FIGURE 10-30
16、DTMF generator.,Joyce Van de VegteFundamentals of Digital Signal Processing,Copyright 2002 by Pearson Education, Inc.Upper Saddle River, New Jersey 07458All rights reserved.,FIGURE 10-31 DTMF receiver.,Joyce Van de VegteFundamentals of Digital Signal Processing,Copyright 2002 by Pearson Education, Inc.Upper Saddle River, New Jersey 07458All rights reserved.,