第6章IIR數(shù)字濾波器設計(10下)

1、第第 6 章章IIR數(shù)字濾波器設計數(shù)字濾波器設計Design of IIR Digital Filters22021-12-13本章主要內(nèi)容本章主要內(nèi)容n模擬濾波器設計模擬濾波器設計巴特沃思、切比雪夫模擬濾波器設計巴特沃思、切比雪夫模擬濾波器設計模擬濾波器的頻帶變換模擬濾波器的頻帶變換n數(shù)字濾波器（數(shù)字濾波器（Digital Filter，DF）設計）設計用脈沖響應不變法設計用脈沖響應不變法設計IIR數(shù)字濾波器數(shù)字濾波器用雙線性變換法設計用雙線性變換法設計IIR數(shù)字濾波器數(shù)字濾波器IIR數(shù)字濾波器的頻率變換數(shù)字濾波器的頻率變換32021-12-13數(shù)字濾波器的分類數(shù)字濾波器的分類n經(jīng)典濾波器（

2、一般濾波器）：經(jīng)典濾波器（一般濾波器）：信號和干擾的頻帶互不重疊時采用信號和干擾的頻帶互不重疊時采用n現(xiàn)代濾波器：現(xiàn)代濾波器：信號和干擾的頻帶相互重疊時采用（例如：維納濾波信號和干擾的頻帶相互重疊時采用（例如：維納濾波器、卡爾曼濾波器、自適應濾波器等）器、卡爾曼濾波器、自適應濾波器等）n經(jīng)典濾波器的基本功能經(jīng)典濾波器的基本功能高通、低通、帶通、帶阻高通、低通、帶通、帶阻n數(shù)字濾波器分為：數(shù)字濾波器分為：無限長單位脈沖響應濾波器（無限長單位脈沖響應濾波器（IIR DF） 遞歸系統(tǒng)遞歸系統(tǒng) 有限長單位脈沖響應濾波器（有限長單位脈沖響應濾波器（FIR DF）非遞歸系統(tǒng)非遞歸系統(tǒng) 42021-12-1

3、3濾波器的設計方法濾波器的設計方法n直接設計法直接設計法在時域或頻域直接設計數(shù)字濾波器在時域或頻域直接設計數(shù)字濾波器n間接設計法間接設計法先根據(jù)指標要求設計對應的模擬濾波器先根據(jù)指標要求設計對應的模擬濾波器再將模擬濾波器轉換為數(shù)字濾波器再將模擬濾波器轉換為數(shù)字濾波器本章介紹間接設計法本章介紹間接設計法6.1 模擬濾波器設計模擬濾波器設計62021-12-13模擬濾波器設計模擬濾波器設計n模擬濾波器的設計步驟模擬濾波器的設計步驟確定指標確定指標選擇濾波器的類型選擇濾波器的類型計算濾波器的階數(shù)計算濾波器的階數(shù)查表或計算濾波器的參數(shù)，確定其系統(tǒng)函數(shù)查表或計算濾波器的參數(shù)，確定其系統(tǒng)函數(shù)綜合實現(xiàn)及調(diào)試

4、綜合實現(xiàn)及調(diào)試n濾波器的傳輸函數(shù)（頻率響應特性）：濾波器的傳輸函數(shù)（頻率響應特性）：選頻濾波器一般只考慮幅頻特性，對相頻特性不作要求。選頻濾波器一般只考慮幅頻特性，對相頻特性不作要求。幅頻特性體現(xiàn)了各頻率成分在幅度上的衰減幅頻特性體現(xiàn)了各頻率成分在幅度上的衰減,而相頻特性體現(xiàn)的是而相頻特性體現(xiàn)的是不同成分在時間上的延時。不同成分在時間上的延時。對輸出波形有明確要求時，則需考慮對輸出波形有明確要求時，則需考慮線性相位線性相位問題。問題。()() |()|jH jH je 72021-12-136.1.1 模擬濾波器的設計指標模擬濾波器的設計指標n 通帶邊界（截止）頻率，通帶邊界（截止）頻率， 阻

5、帶邊界（截止）頻率，阻帶邊界（截止）頻率， 3db通帶截止頻率通帶截止頻率n系統(tǒng)通帶和阻帶的誤差要求系統(tǒng)通帶和阻帶的誤差要求通帶特性要求：通帶特性要求： 通帶最大衰減通帶最大衰減 ，或通帶峰值波紋，或通帶波紋幅，或通帶峰值波紋，或通帶波紋幅度度 ，用分貝表示：，用分貝表示：p s c 21()1,1apHj 22120lg110lg(1)dBp /10101p p 82021-12-13阻帶特性要求：阻帶特性要求： 阻帶最小衰減阻帶最小衰減 ，或，或 阻帶峰值波紋，或阻帶峰值波紋，或 阻帶波紋幅阻帶波紋幅度度 ，用分貝表示：，用分貝表示：1(),asHjA /2010sA 120lg()20l

6、gdBsAA s A92021-12-13n損耗函數(shù)損耗函數(shù)（衰減函數(shù)衰減函數(shù)）：用來描述濾波器幅頻響應特性）：用來描述濾波器幅頻響應特性n3dB 通帶截止頻率通帶截止頻率 c：n損耗函數(shù)相對放大了小幅度特性：損耗函數(shù)相對放大了小幅度特性：21()20lg20lg()10lg()()aaaHjHjHj ()3dBc 一般幅頻特性一般幅頻特性損耗函數(shù)損耗函數(shù)()1/20.707acHj102021-12-13兩個附加參數(shù)兩個附加參數(shù)n過渡比過渡比或或選擇性選擇性參數(shù)參數(shù)反應過渡帶的性能，過渡帶越窄，反應過渡帶的性能，過渡帶越窄，k 值趨近于值趨近于 1對對“低通濾波器低通濾波器”：n偏離偏離參數(shù)

7、參數(shù) 越小、越小、 越大，越大， 越小，越小， 通帶、阻帶的紋波越小通帶、阻帶的紋波越小一般有一般有 psk 1k 121kA 1k11kA112021-12-13n模擬濾波器的設計模擬濾波器的設計模擬濾波器的理論和設計方法已發(fā)展得相當成熟，且有若干典型模擬濾波器的理論和設計方法已發(fā)展得相當成熟，且有若干典型的模擬濾波器供選擇。的模擬濾波器供選擇。這些濾波器都有嚴格的設計公式、現(xiàn)成的曲線和圖表供設計人員這些濾波器都有嚴格的設計公式、現(xiàn)成的曲線和圖表供設計人員使用。使用。 n典型濾波器典型濾波器巴特沃斯（巴特沃斯（Butterworth）濾波器）濾波器：具有單調(diào)下降的幅頻特性；：具有單調(diào)下降的幅

8、頻特性；切比雪夫（切比雪夫（Chebyshev）濾波器）濾波器：幅頻特性在通帶或阻帶內(nèi)有波：幅頻特性在通帶或阻帶內(nèi)有波動起伏，這可提高選擇性；動起伏，這可提高選擇性；橢圓（橢圓（Ellipse）濾波器）濾波器：幅頻特性在通帶和阻帶都有波動，其選：幅頻特性在通帶和阻帶都有波動，其選擇性最好；擇性最好；貝塞爾（貝塞爾（Bessel）濾波器）濾波器：具有較好的線性相位。：具有較好的線性相位。122021-12-136.1.2 巴特沃斯（巴特沃斯（Butterworth）模擬低通濾）模擬低通濾波器設計波器設計n巴特沃斯巴特沃斯低通低通濾波器濾波器N 階巴特沃斯低通濾波器的階巴特沃斯低通濾波器的幅度平方

9、函數(shù)幅度平方函數(shù)為：為： 其中其中N 為濾波器的階次，為濾波器的階次， 為為 3dB通帶截止頻率通帶截止頻率。221|()|1()aNcHj c 132021-12-13巴特沃斯低通濾波器特點巴特沃斯低通濾波器特點：n在在 點，點， 的的 n ( n p ，隨，隨 ， n當當 =0時，時， N為偶數(shù)，為偶數(shù)， ，min ，N為奇數(shù)，為奇數(shù)， ， max，22202211()1cos arccos(0)1cos ()2aHjNN 2201()1aHj 20()1aHj 2cos ()12N 2cos ()02N 211 1 1,()1,NppC 2a()Hj 2(),()0NapCHj 3920

10、21-12-13n阻帶邊界頻率阻帶邊界頻率nN 大于等于上式的最小整數(shù)大于等于上式的最小整數(shù)2222221()1()111cosharcosh()assNpspHjCAN 2111arcosharcosh()1arcosh()arcoshspkNkA 2arcosh( )ln(1)xxx 其其中中/2010sA /10101p 402021-12-13濾波器的系統(tǒng)函數(shù)濾波器的系統(tǒng)函數(shù)112( ),()NNpakkkNkkHssjss (21)(21)sin,sin22kpkpkkNN 1/2221111,22N412021-12-13實際應用中可以利于查表的方法求得系統(tǒng)函數(shù)：實際應用中可以利于

11、查表的方法求得系統(tǒng)函數(shù)：1. 按按 求出濾波器階次求出濾波器階次N；2. 由由N和和 （或（或 p）查表得到歸一化濾波器分母多項式：查表得到歸一化濾波器分母多項式：3. 歸一化切比雪夫低通濾波器函數(shù)為歸一化切比雪夫低通濾波器函數(shù)為:4. 去歸一化，得到通帶截止頻率為去歸一化，得到通帶截止頻率為 p的低通濾波器函數(shù)：的低通濾波器函數(shù)：21arcosharcoshspAN 1110( )(1)NNNND ppapa paa 11 (2)( )( )NG pD p ( )( )psapHsG p 422021-12-13切比雪夫切比雪夫型型通帶內(nèi)單調(diào)下降，阻帶內(nèi)等波紋通帶內(nèi)單調(diào)下降，阻帶內(nèi)等波紋n幅

12、度平方函數(shù)：幅度平方函數(shù)：n階數(shù)階數(shù) N 的計算的計算 psA 【】2111arcosharcosh()1arcosh()arcosh()spNAkk 2221()()1()asNpsNHjCC 432021-12-13系統(tǒng)函數(shù)系統(tǒng)函數(shù)11()( ),1,2,()NkkakkkNkkszHsCpjkNsp (21)(21)sincos22kpkpkkNN 221/211122NAA(21)cos2skzjkN 2222skskkkkkkk 442021-12-136.1.4橢圓濾波器（考爾濾波器）橢圓濾波器（考爾濾波器）n幅值響應在通帶和阻帶內(nèi)都是等波紋的幅值響應在通帶和阻帶內(nèi)都是等波紋的n對

13、于給定的階數(shù)和給定的波紋要求，橢圓濾波器能獲得較對于給定的階數(shù)和給定的波紋要求，橢圓濾波器能獲得較其它濾波器更窄的過渡帶寬，就這點而言，橢圓濾波器是其它濾波器更窄的過渡帶寬，就這點而言，橢圓濾波器是最優(yōu)的。最優(yōu)的。n通帶和阻帶內(nèi)波紋固定時，階數(shù)越高，過渡帶越窄；通帶和阻帶內(nèi)波紋固定時，階數(shù)越高，過渡帶越窄；n階數(shù)固定，通帶和阻帶紋波越小，過渡帶越寬；階數(shù)固定，通帶和阻帶紋波越小，過渡帶越寬；452021-12-134102040pspspsNdBdBdB =1dB=1dB=0.1dB=0.1dB=0.05dB=0.05dB203, 4, 6psdBN =1dB=1dB462021-12-13n

14、振幅平方函數(shù)為振幅平方函數(shù)為RN( , L) 雅可比橢圓函數(shù)雅可比橢圓函數(shù)L 表示波紋性質(zhì)的參量表示波紋性質(zhì)的參量22221()()1(,)aNAHjRL 472021-12-13 下圖為典型的橢園濾波器振幅平方函數(shù)下圖為典型的橢園濾波器振幅平方函數(shù) 圖中圖中和和A的定義的定義 與與 切比雪夫濾波器切比雪夫濾波器 相同相同482021-12-136.1.7模擬濾波器的比較模擬濾波器的比較n相同階數(shù)的頻率響應特性相同階數(shù)的頻率響應特性相同通帶最大衰減、阻帶最小衰減，巴特沃思濾波器的過渡帶相同通帶最大衰減、阻帶最小衰減，巴特沃思濾波器的過渡帶最寬；橢圓濾波器過渡帶最窄；最寬；橢圓濾波器過渡帶最窄；

15、n相同指標下，橢圓濾波器階次最低，切比雪夫次之，巴特沃思最高相同指標下，橢圓濾波器階次最低，切比雪夫次之，巴特沃思最高，參數(shù)的靈敏度則恰恰相反；，參數(shù)的靈敏度則恰恰相反；n巴特沃思濾波器單調(diào)下降；巴特沃思濾波器單調(diào)下降；n切比雪夫切比雪夫濾波器通帶等波紋幅頻特性，過渡帶、阻帶單調(diào)下降；濾波器通帶等波紋幅頻特性，過渡帶、阻帶單調(diào)下降；n切比雪夫切比雪夫濾波器阻帶等波紋幅頻特性，通帶、過渡帶單調(diào)下降；濾波器阻帶等波紋幅頻特性，通帶、過渡帶單調(diào)下降；n橢圓濾波器通帶、阻帶等波紋幅頻特性，過渡帶單調(diào)下降。橢圓濾波器通帶、阻帶等波紋幅頻特性，過渡帶單調(diào)下降。n貝塞爾（貝塞爾（Bessel）濾波器在整個通

16、帶內(nèi)）濾波器在整個通帶內(nèi)逼近逼近線性相位特性，但過渡線性相位特性，但過渡帶特性較差。帶特性較差。492021-12-136.1.8 頻率變換與高通、帶通及帶阻濾波器頻率變換與高通、帶通及帶阻濾波器的設計的設計高通濾波器性能指標高通濾波器性能指標1()HPHj112shphn通帶邊界頻率：通帶邊界頻率：n阻帶邊界頻率：阻帶邊界頻率：n通帶紋波幅度：通帶紋波幅度：n阻帶紋波幅度：阻帶紋波幅度：2 1 ph sh 502021-12-13n通帶邊界頻率：通帶邊界頻率：上邊界頻率上邊界頻率下邊界頻率下邊界頻率n阻帶邊界頻率：阻帶邊界頻率：上邊界頻率上邊界頻率下邊界頻率下邊界頻率n通帶紋波幅度：通帶紋波

17、幅度：n阻帶紋波幅度：阻帶紋波幅度：帶通濾波器的性能指標帶通濾波器的性能指標2 1 pl pu sl su 1()BPHj112slplsupu512021-12-13帶阻濾波器的性能指標帶阻濾波器的性能指標n通帶邊界頻率：通帶邊界頻率：上邊界頻率上邊界頻率下邊界頻率下邊界頻率n阻帶邊界頻率：阻帶邊界頻率：上邊界頻率上邊界頻率下邊界頻率下邊界頻率n通帶紋波幅度：通帶紋波幅度：n阻帶紋波幅度：阻帶紋波幅度：2 1 1()BSHj112slplsupupl pu sl su 522021-12-13高通、帶通及帶阻濾波器的設計高通、帶通及帶阻濾波器的設計n已有：歸一化低通原型濾波器已有：歸一化低通

18、原型濾波器通帶邊界頻率為通帶邊界頻率為 1 （歸一化頻率）（歸一化頻率）的模擬濾波器的模擬濾波器n高通、帶通及帶阻濾波器的設計過程：高通、帶通及帶阻濾波器的設計過程：采用頻率變換公式，將濾波器的技術指標轉換成歸一采用頻率變換公式，將濾波器的技術指標轉換成歸一化低通濾波器指標；化低通濾波器指標；利用經(jīng)典濾波器的設計公式、圖表，設計歸一化低通利用經(jīng)典濾波器的設計公式、圖表，設計歸一化低通濾波器的系統(tǒng)函數(shù)濾波器的系統(tǒng)函數(shù) ；通過頻率變換公式，由通過頻率變換公式，由 獲得希望濾波器的系統(tǒng)獲得希望濾波器的系統(tǒng)函數(shù)函數(shù) 。( )G p( )G p( )dHs532021-12-13符號規(guī)定：符號規(guī)定：n歸

19、一化低通原型濾波器：歸一化低通原型濾波器：歸一化復變量歸一化復變量 ，歸一化頻率，歸一化頻率通帶邊界頻率（截止頻率）通帶邊界頻率（截止頻率）n期望模擬濾波器的系統(tǒng)函數(shù)：期望模擬濾波器的系統(tǒng)函數(shù)：n頻率變換公式：頻率變換公式：n于是于是( ),G ppj (1)pp 一一般般( ),dHssj 1( ),( )pF ssFp ( )( )( )dp F sHsG p 1()( )( )ds FpG pHs p542021-12-131. 模擬高通濾波器的設計模擬高通濾波器的設計n低通原型到高通濾波器的映射變換：低通原型到高通濾波器的映射變換：在虛軸上的映射關系簡化為頻率變換關系，在虛軸上的映射關

20、系簡化為頻率變換關系，高通濾波器邊界頻率高通濾波器邊界頻率映射關系映射關系： 低通低通 高通高通pphps pph (,)pjsj 0 ,0 ,ppsspp , 00,phphshshph 1()HPHj112shph1()G p112sppph sh 552021-12-13例例 設計巴特沃思模擬高通濾波器設計巴特沃思模擬高通濾波器,通帶邊界頻率為通帶邊界頻率為fp=4kHz，阻帶邊界頻率為阻帶邊界頻率為fs=1kHz，通帶最大衰減為，通帶最大衰減為0.1dB（fp處），處），阻帶最小衰減阻帶最小衰減 s=40dB。解：高通濾波器指標解：高通濾波器指標 1）確定相應低通原型濾波器的指標）確定

21、相應低通原型濾波器的指標( )( )pphHPpsHsG p 4kHz,1kHz,0.1dB,40dBpspsff 24000421000pphss 1p 0.1dB40dBps 562021-12-132) 設計巴特沃斯低通濾波器。設計巴特沃斯低通濾波器。n 確定階次確定階次N：n可得可得Nn按照以上技術指標設計巴特沃斯低通濾波器。按照以上技術指標設計巴特沃斯低通濾波器。 /2040/201010100sA /100.1/101011010.1526204,p 10.254ppssk 120.15262040.0015262899991kA 1lg4.75lgkNNk11/102241.59

22、9999(1)scNA 140.1dB40dBpsps 122(1)scNA 572021-12-13n歸一化的歸一化的 5 階濾波器的系統(tǒng)函數(shù)階濾波器的系統(tǒng)函數(shù)3）高通濾波器的系統(tǒng)函數(shù)）高通濾波器的系統(tǒng)函數(shù)54325( )3.23615.23615.23613.23611Dpppppp 55423324554321( )(/)3.23615.23615.23613.236110.24055.153313.27821.144520.810110.2405NcccccccG pDppppppppppp 54420( )( )5.1073 109.7921 10pphHPpsHsG psss 58

23、2021-12-13n另一設計途徑：另一設計途徑：模擬高通指標模擬高通指標模擬低通指標模擬低通指標(通帶截止頻率與高通一致通帶截止頻率與高通一致)模擬低通濾波器模擬低通濾波器設計歸一化模擬設計歸一化模擬低通濾波器低通濾波器G(p)模擬高通濾波器模擬高通濾波器2( )( )plh cqp 2pqs 592021-12-13MATLABwp=1;ws=4;Rp=0.1;As=40; %設置濾波器指標參數(shù)設置濾波器指標參數(shù)N,wc=buttord(wp,ws,Rp,As,s); %濾波器濾波器G(p)階數(shù)階數(shù)N和和3dB截止截止頻率頻率B,A=butter(N,wc,s); %計算低通濾波器計算低通

24、濾波器G(p)系統(tǒng)函數(shù)分子分母多項系統(tǒng)函數(shù)分子分母多項式系數(shù)式系數(shù)wph=2*pi*4000; %高通模擬濾波器通帶邊界頻率高通模擬濾波器通帶邊界頻率BH,AH=lp2hp(B,A,wph) %低通到高通轉換低通到高通轉換602021-12-132. 低通到帶通的頻率變換低通到帶通的頻率變換:n低通原型到帶通濾波器的頻率變換低通原型到帶通濾波器的頻率變換在虛軸上在虛軸上n通帶邊界頻率：上邊界頻率通帶邊界頻率：上邊界頻率 下邊界頻率下邊界頻率n帶通濾波器的帶寬帶通濾波器的帶寬n帶通濾波器的中心頻率帶通濾波器的中心頻率220pWspBs 220pWB pjsj pl pu WpuplB 0pupl

25、 612021-12-13映射關系映射關系0,ppss 0,plpupuplsuslsusl 1()Gp112spp1()BPHj112slplsupu622021-12-13減少減少 ，或增加，或增加增加增加 ，或減少，或減少帶通濾波器的系統(tǒng)函數(shù)帶通濾波器的系統(tǒng)函數(shù)n對帶通濾波器，可以證明，應該有關系：對帶通濾波器，可以證明，應該有關系：n若給定邊界頻率不滿足條件，則可改變參數(shù)，提高指標若給定邊界頻率不滿足條件，則可改變參數(shù)，提高指標 ， 220( )( )pWsBPpB sHsG p 20plpuslsu pplslusu ppluluss pl sl su pu 632021-12-13

26、例：設計巴特沃思模擬帶通濾波器，通帶上下邊界頻率分別為例：設計巴特沃思模擬帶通濾波器，通帶上下邊界頻率分別為4kHz和和7kHz，阻帶上、下邊界頻率分別為，阻帶上、下邊界頻率分別為2kHz和和9kHz， 通帶內(nèi)最大衰減通帶內(nèi)最大衰減p=1dB，阻帶最小衰減，阻帶最小衰減s=20dB。解：解：1）帶通濾波器指標）帶通濾波器指標 增大增大 2）低通原型濾波器指標）低通原型濾波器指標WpuplBff4,7:12,9:20plpupslsusfkHz fkHzdBfkHz fkHzdB plpuslsufff f ,3.1111plpuslslsuffffkHzf 即即取取222205.53.1111

27、1.96333.1111slspWffffBf 120psdBdB 0puplfff 1p 642021-12-133）設計低通原型濾波器：）設計低通原型濾波器：4）求模擬帶通）求模擬帶通H(s)： 220( )( )pWsBPpB sHsG p 652021-12-133.低通到帶阻的頻率變換低通到帶阻的頻率變換:n低通原型到帶阻濾波器的頻率變換低通原型到帶阻濾波器的頻率變換在虛軸上在虛軸上n阻帶邊界頻率：上邊界頻率阻帶邊界頻率：上邊界頻率 下邊界頻率下邊界頻率n帶阻濾波器的帶寬帶阻濾波器的帶寬n帶阻濾波器的中心頻率帶阻濾波器的中心頻率220WpB sps 220WpB pjsj sl su

28、 WsuslB 0 662021-12-13映射關系映射關系0,ppss 0,slsususlpuplpupl 672021-12-13減少減少 ，或增加，或增加減少減少 ，或增加，或增加帶阻濾波器的系統(tǒng)函數(shù)帶阻濾波器的系統(tǒng)函數(shù)n可以證明：可以證明：n如給定邊界頻率不滿足改條件，改變參數(shù)，提高指標如給定邊界頻率不滿足改條件，改變參數(shù)，提高指標 ， 220( )( )WpB sBSpsHsG p 20plpuslsu plpuslsu plpuslsu pu su sl pl 6.2 IIR數(shù)字濾波器設計數(shù)字濾波器設計692021-12-136.2 IIR數(shù)字濾波器設計數(shù)字濾波器設計n目標：滿足

29、給定頻率響應指標、因果穩(wěn)定的系統(tǒng)函數(shù)目標：滿足給定頻率響應指標、因果穩(wěn)定的系統(tǒng)函數(shù)n間接法設計過程間接法設計過程確定數(shù)字濾波器的指標確定數(shù)字濾波器的指標轉換成轉換成過渡過渡模擬濾波器的指標模擬濾波器的指標設計設計過渡過渡模擬濾波器模擬濾波器將將過渡過渡模擬濾波器轉換為數(shù)字濾波器模擬濾波器轉換為數(shù)字濾波器n指標轉換指標轉換,psps ( )H z( )aHs( )H z22pppsssssssfFFfFF 0,0, 2 sF n實際應用中，常希望用數(shù)字濾波器實際應用中，常希望用數(shù)字濾波器對模擬信號進行處理。如果給定的對模擬信號進行處理。如果給定的指標為模擬指標，則指標轉換為：指標為模擬指標，則指

30、標轉換為：(),(),ppsspsspspsffF (),(),ppsspsff702021-12-13過渡模擬濾波器轉換為數(shù)字濾波器的要求：過渡模擬濾波器轉換為數(shù)字濾波器的要求：n保證因果穩(wěn)定性，保證因果穩(wěn)定性， Ha(s) 的因果穩(wěn)定性映射成的因果穩(wěn)定性映射成 H(z)后保持后保持不變，即不變，即 S 平面的平面的左半平面左半平面 ReS0 應映射到應映射到 Z 平面的平面的單位圓以內(nèi)單位圓以內(nèi)|z|1右半平面右半平面 0單位圓內(nèi)部單位圓內(nèi)部r1左半平面左半平面 0單位圓單位圓r=1虛軸虛軸 =0Z平面平面S平面平面752021-12-13S平面到平面到Z平面的平面的映射是多值映射。映射是

31、多值映射。22/ T2T s s=0 正實軸正實軸 零頻零頻=0T 輻射線輻射線 角度角度: =0 實軸實軸 零頻零頻 =0 平行直線平行直線 頻率頻率:Z平面平面S平面平面/2s /22ssTT結論結論n穩(wěn)定性穩(wěn)定性如果模擬濾波器是因果穩(wěn)定的，則所有極點如果模擬濾波器是因果穩(wěn)定的，則所有極點 都在都在 S 左左半平面，那么變換后半平面，那么變換后 H(z) 的極點也都在單位圓以內(nèi)，的極點也都在單位圓以內(nèi)，因此數(shù)字濾波器保持因果穩(wěn)定。因此數(shù)字濾波器保持因果穩(wěn)定。n頻率響應特性的近似性頻率響應特性的近似性762021-12-131()()jTaskX eXjjkT 12()()jakkX eXj

32、TT 772021-12-13( )()( )at nTh nh nTh t ()jH e ()aHj 12()()jakkH eHjTT 數(shù)字濾波器的頻數(shù)字濾波器的頻響并不是簡單的響并不是簡單的重現(xiàn)模擬濾波器重現(xiàn)模擬濾波器的頻響，而是模的頻響，而是模擬濾波器頻響的擬濾波器頻響的周期延拓！周期延拓！()H j 782021-12-13n無混疊失真的的條件：無混疊失真的的條件：帶限、滿足采樣定理，即帶限、滿足采樣定理，即n混疊失真的影響混疊失真的影響相似性變差、無法滿足阻帶衰減指標相似性變差、無法滿足阻帶衰減指標n脈沖響應不變法的優(yōu)缺點：脈沖響應不變法的優(yōu)缺點：（1）頻率坐標變換是線性的，數(shù)字濾

33、波器在無頻率混疊時能較）頻率坐標變換是線性的，數(shù)字濾波器在無頻率混疊時能較好地重現(xiàn)原模擬濾波器的頻率特性；好地重現(xiàn)原模擬濾波器的頻率特性；（2）由于數(shù)字域的脈沖響應模仿了模擬域的沖激響應，故時域）由于數(shù)字域的脈沖響應模仿了模擬域的沖激響應，故時域特性逼近好；特性逼近好；（3）模擬濾波器頻響在折疊頻率以上衰減越大，失真則越小，）模擬濾波器頻響在折疊頻率以上衰減越大，失真則越小，這時，采用脈沖響應不變法設計的數(shù)字濾波器才能得到良好這時，采用脈沖響應不變法設計的數(shù)字濾波器才能得到良好的效果。的效果。（4）由于頻率混疊現(xiàn)象，故）由于頻率混疊現(xiàn)象，故不適合不適合高通、帶阻濾波器高通、帶阻濾波器的設計。的

34、設計。()0,2saHj 當當792021-12-13增益的補償增益的補償2()1()jakkH eHjTT 11( )1kNks TkTAH zez ( )( )at nTh nTh t 2()()jakkH eHjT 802021-12-13共軛極點的合并共軛極點的合并n避免復數(shù)乘法避免復數(shù)乘法n其中其中11( )1kNks TkTAH zez 101121112111kkkkkks Ts TkkTATAbb zazazezez kkkkkksjsj 01122Re,2Re2cos(),2kkkkTTkkkkTTkkkbTAbTeA eaeTae /210112112( )1Nkkkkkb

35、b zH zazaz 812021-12-13例：例： 將一個具有如下傳遞函數(shù)將一個具有如下傳遞函數(shù) 的模擬濾波器數(shù)字化。的模擬濾波器數(shù)字化。 解：解： 211( )(1)(3)13H sssss 11311( )11TTH zzeze 131342()1()TTTTTzeezeeez 822021-12-13數(shù)字濾波器的頻率響應為數(shù)字濾波器的頻率響應為: 顯然顯然 與采樣間隔與采樣間隔 T （采樣頻率）有關（采樣頻率）有關, 如圖如圖b,T 越小越小,衰減越大衰減越大,混疊越小混疊越小,當當 Fs=24Hz ,混疊可忽略不計混疊可忽略不計,為為什么混什么混疊疊呢呢? 模擬濾波器的頻率響應為模

36、擬濾波器的頻率響應為:222()( )(1)(3)(3)4asjHjH sjjj 3342()()( )1 ()jTTjjTTjTjz eeeeH eH zeeeee ()jH e 832021-12-13842021-12-13例例 6.2.2 （pp. 161162）用脈沖響應不變法設計數(shù)字低通濾器，要求通帶與阻帶具有用脈沖響應不變法設計數(shù)字低通濾器，要求通帶與阻帶具有單調(diào)特性。指標如下：單調(diào)特性。指標如下：設計基本步驟：設計基本步驟：（1）將數(shù)字濾波器指標轉換為）將數(shù)字濾波器指標轉換為過渡過渡模擬濾波器指標：模擬濾波器指標：0.2rad,1dB,0.35rad,10 dBppss/,/,

37、ppsspsTTT 由保持不變（2）設計相應的）設計相應的過渡過渡模擬濾波器：模擬濾波器： （根據(jù)單調(diào)特性，應選用巴特沃思濾波器）（根據(jù)單調(diào)特性，應選用巴特沃思濾波器）（3）由由過渡過渡模擬濾波器模擬濾波器Ha(s)轉換為數(shù)字濾波器轉換為數(shù)字濾波器H(z)： 1( )NkakkAHsss 11( )1kNks TkAH zez 852021-12-13例例 6.2.2（續(xù)）（續(xù)）注意：當注意：當給定數(shù)字濾波器技術指標時，數(shù)字濾波器的設計中，給定數(shù)字濾波器技術指標時，數(shù)字濾波器的設計中，采樣周期采樣周期 T 的取值對頻譜混疊程度影響很?。ㄈ缋械娜≈祵︻l譜混疊程度影響很?。ㄈ缋蠺=1s 和和

38、T=0.1s 的濾波器特性一致）。的濾波器特性一致）。 這是因為數(shù)字濾波器指標中有這是因為數(shù)字濾波器指標中有 ，轉換到過渡模擬，轉換到過渡模擬濾波器指標中，濾波器指標中，T 取不同值時一定有取不同值時一定有 ，模擬濾波器在模擬濾波器在1/2采樣頻率以上響應特性一定都很?。ㄗ璨蓸宇l率以上響應特性一定都很小（阻帶），再變換到數(shù)字濾波器時引起的混疊都較小。帶），再變換到數(shù)字濾波器時引起的混疊都較小。s / 2sssTTF 862021-12-13小結：小結：( )( )()( )( )aaaH sh th nTh nzH z模模擬擬濾濾波波器器采采樣樣變變換換得得數(shù)數(shù)字字濾濾波波器器1( )Nkak

39、kAHsss 若若 1) 脈沖響應不變法的一個重要特點是頻率坐標的變換是線性脈沖響應不變法的一個重要特點是頻率坐標的變換是線性的，的，與與是線性關系。是線性關系。 因此如果模擬濾波器的頻響帶限于折疊頻率以內(nèi)的話，通過變因此如果模擬濾波器的頻響帶限于折疊頻率以內(nèi)的話，通過變換后數(shù)字濾波器的頻響可不失真地反映原響應與頻率的關系。換后數(shù)字濾波器的頻響可不失真地反映原響應與頻率的關系。 例如線性相位的貝塞爾低通濾波器，通過脈沖響應不變法得例如線性相位的貝塞爾低通濾波器，通過脈沖響應不變法得到的仍是線性相位的低通數(shù)字濾波器。到的仍是線性相位的低通數(shù)字濾波器。()()/jaH eHjT 11( )1kNk

40、s TkAH zez 則則872021-12-13 2) 在某些場合，要求數(shù)字濾波器在時域上能模仿模擬濾波器在某些場合，要求數(shù)字濾波器在時域上能模仿模擬濾波器的功能時，如要實現(xiàn)時域沖激響應的模仿，一般使用脈沖響應的功能時，如要實現(xiàn)時域沖激響應的模仿，一般使用脈沖響應不變法。不變法。 3) 如果如果Ha(s) 是穩(wěn)定的，即其極點在是穩(wěn)定的，即其極點在 S 平面的左半平面，映平面的左半平面，映射后得到的射后得到的 H(z) 也是穩(wěn)定的，其極點在也是穩(wěn)定的，其極點在 Z 平面單位圓內(nèi)。平面單位圓內(nèi)。 4) 脈沖響應不變法的最大缺點：有頻響特性的周期延拓效應脈沖響應不變法的最大缺點：有頻響特性的周期延

41、拓效應，因此只能用于帶限的頻響特性。如衰減特性很好的低通或帶，因此只能用于帶限的頻響特性。如衰減特性很好的低通或帶通濾波器，高頻衰減越大，頻響的混疊效應越?。恢劣诟咄ê屯V波器，高頻衰減越大，頻響的混疊效應越?。恢劣诟咄ê蛶ё铻V波器，由于它們在高頻部分不衰減，因此將完全混淆在帶阻濾波器，由于它們在高頻部分不衰減，因此將完全混淆在低頻響應中。因此不能用于設計高通和帶阻濾波器。低頻響應中。因此不能用于設計高通和帶阻濾波器。 882021-12-136.2.2 雙線性變換法設計雙線性變換法設計IIR數(shù)字濾波器數(shù)字濾波器n克服脈沖響應不變法的頻譜混疊現(xiàn)象，采用非線性頻率壓克服脈沖響應不變法的頻譜混疊現(xiàn)

42、象，采用非線性頻率壓縮方法?？s方法。n雙線性變換法的基本設計思想雙線性變換法的基本設計思想脈沖響應不變法：波形逼近脈沖響應不變法：波形逼近雙線性變換法：算法逼近。雙線性變換法：算法逼近。用線性常系數(shù)差分方程用線性常系數(shù)差分方程逼近逼近線性常系數(shù)微分方程線性常系數(shù)微分方程892021-12-13n微分方程到差分方程的近似微分方程到差分方程的近似微分方程微分方程近似近似差分方程差分方程11( )( )NNkaakkkkAHsHsss ( )kakkAHsss ( )( )( )akakadyts ytA xtdt ( ) ( )(1)/,( )()( ) ( )(1)/ 2( ) ( )(1)/

43、2,( )()aaaaadyty ny nTy nynTdtyty ny nxtx nx nx nxnT ( )(1)/ ( )(1)/ 2 ( )(1)/ 2kky ny nTsy ny nA x nx n902021-12-13n系統(tǒng)函數(shù)系統(tǒng)函數(shù)n 有有nS 域到域到 Z 域的映射為雙線性變換域的映射為雙線性變換112 1111( )( )2 11kzkaksTzkAHzHszsTz 112 11( )( )zasTzH zHs 112 11zsTz 912021-12-13 脈沖響應不變法的主要缺點是頻響交疊產(chǎn)生的混淆，這脈沖響應不變法的主要缺點是頻響交疊產(chǎn)生的混淆，這是從是從S平面到平

44、面到Z平面的標準變換平面的標準變換zesT的多值對應關系導致的的多值對應關系導致的,為了克服這一缺點，設想變換分為兩步：為了克服這一缺點，設想變換分為兩步： 第一步：將整個第一步：將整個S平面壓縮到平面壓縮到S1平面的一條橫帶里；平面的一條橫帶里； 第二步：通過標準變換關系將此橫帶變換到整個第二步：通過標準變換關系將此橫帶變換到整個Z平面平面上去。上去。 由此建立由此建立S平面與平面與Z平面一一對應的單值關系，消除多值平面一一對應的單值關系，消除多值性，也就消除了混淆現(xiàn)象。性，也就消除了混淆現(xiàn)象。922021-12-13雙線性變換法的映射關系雙線性變換法的映射關系932021-12-13n 正

45、切變換實現(xiàn)頻率壓縮：正切變換實現(xiàn)頻率壓縮：121tan()2TT 式中式中T仍是采樣間隔，當仍是采樣間隔，當1從從-/T經(jīng)過經(jīng)過0變化到變化到/T時，時，則由則由-經(jīng)過經(jīng)過0變化到變化到+，實現(xiàn)了，實現(xiàn)了S平面上整個虛軸完全壓平面上整個虛軸完全壓縮到縮到S1平面上虛軸的平面上虛軸的/T之間的轉換。之間的轉換?？紤]到：考慮到：令令 s=j ， s1=j 1則有則有111212 1()21s Ts TesthTTTe11111221122sin()22212th()2cos()2TTjjTTjjTeejjTTTTTee 942021-12-13n 再通過再通過 轉換到轉換到Z平面上，得到：平面上，

46、得到：1s Tze112 11zsTz952021-12-13雙線性變換的性能分析雙線性變換的性能分析n穩(wěn)定性穩(wěn)定性映射關系映射關系1122TTzSS1()1()22TTjj 222221/ 2(/ 2)1/ 2(/ 2)TTzTT 010101zzz 因果穩(wěn)定因果穩(wěn)定( )aHs因果穩(wěn)定因果穩(wěn)定( )H z962021-12-13非線性頻率壓縮非線性頻率壓縮n非線性關系消除了頻率混疊，但會帶來幅度和相位失真。非線性關系消除了頻率混疊，但會帶來幅度和相位失真。/2/2/2/2/2/2/2/2/2/22 11222sin(/ 2cos(/ 2)2tan()2jjjjjjjjjjjjejTeeee

47、T eeeeeT eejTjT ) )0112 11zsTz 2tan()2T 即即：972021-12-13頻率響應的畸變頻率響應的畸變n幅頻特性幅頻特性 相頻特性相頻特性0 0| )(|jeH| )(|jHa0 00 00 0)(jeH)( jHa0 00 0982021-12-13雙線性變換法設計數(shù)字低通濾波器的步驟雙線性變換法設計數(shù)字低通濾波器的步驟1）確定數(shù)字低通技術指標：）確定數(shù)字低通技術指標： 通帶邊界頻率通帶邊界頻率 、通帶衰減、通帶衰減 、 阻帶邊界頻率阻帶邊界頻率 、阻帶衰減、阻帶衰減 ；2）將數(shù)字低通指標轉換成）將數(shù)字低通指標轉換成過渡過渡模擬低通指標：模擬低通指標： 不

48、變不變 邊界頻率非線性頻率預畸變校正：邊界頻率非線性頻率預畸變校正：3）設計模擬低通濾波器；）設計模擬低通濾波器；4）轉換成數(shù)字低通濾波器：）轉換成數(shù)字低通濾波器：p s p s 2tan()2T ps22tan()tan()22pspsTT 112 11( )( )zasTzH zHs 992021-12-13T 的選擇的選擇n在給定數(shù)字濾波器在給定數(shù)字濾波器 技術指標情況下技術指標情況下T 可以任意選可以任意選如選如選 T2，則，則n若模擬濾波器已給定，直接變換成數(shù)字濾波器時若模擬濾波器已給定，直接變換成數(shù)字濾波器時T 如何選？如何選？tan()2 tan()tan()22psps 1111zsz 1002021-12-13例：用雙線性變換法設計低通濾波器，通帶阻帶單調(diào)下降，例：用雙線性變換法設計低通濾波器，通帶阻帶單調(diào)下降，指標要求如下：指標要求如下：n解：解：1）轉換成）轉換成過渡過渡模擬低通指標（取模擬低通指標（取 T2） 2）設計過渡模擬濾波器，選用巴特沃思濾波器）設計過渡模擬濾波器，選用巴特沃思濾波器0.2rad1dB0.35rad10dBppss1dB10dB220.2tan()tan()0.32492rad