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

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

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

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

4、綜合實(shí)現(xiàn)及調(diào)試n濾波器的傳輸函數(shù)（頻率響應(yīng)特性）：濾波器的傳輸函數(shù)（頻率響應(yīng)特性）：選頻濾波器一般只考慮幅頻特性，對(duì)相頻特性不作要求。選頻濾波器一般只考慮幅頻特性，對(duì)相頻特性不作要求。幅頻特性體現(xiàn)了各頻率成分在幅度上的衰減幅頻特性體現(xiàn)了各頻率成分在幅度上的衰減,而相頻特性體現(xiàn)的是而相頻特性體現(xiàn)的是不同成分在時(shí)間上的延時(shí)。不同成分在時(shí)間上的延時(shí)。對(duì)輸出波形有明確要求時(shí)，則需考慮對(duì)輸出波形有明確要求時(shí)，則需考慮線性相位線性相位問題。問題。()() |()|jH jH je 72021-12-136.1.1 模擬濾波器的設(shè)計(jì)指標(biāo)模擬濾波器的設(shè)計(jì)指標(biāo)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ù)）：用來描述濾波器幅頻響應(yīng)特性）：用來描述濾波器幅頻響應(yīng)特性n3dB 通帶截止頻率通帶截止頻率 c：n損耗函數(shù)相對(duì)放大了小幅度特性：損耗函數(shù)相對(duì)放大了小幅度特性：21()20lg20lg()10lg()()aaaHjHjHj ()3dBc 一般幅頻特性一般幅頻特性損耗函數(shù)損耗函數(shù)()1/20.707acHj102021-12-13兩個(gè)附加參數(shù)兩個(gè)附加參數(shù)n過渡比過渡比或或選擇性選擇性參數(shù)參數(shù)反應(yīng)過渡帶的性能，過渡帶越窄，反應(yīng)過渡帶的性能，過渡帶越窄，k 值趨近于值趨近于 1對(duì)對(duì)“低通濾波器低通濾波器”：n偏離偏離參數(shù)

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

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

9、函數(shù)幅度平方函數(shù)為：為： 其中其中N 為濾波器的階次，為濾波器的階次， 為為 3dB通帶截止頻率通帶截止頻率。221|()|1()aNcHj c 132021-12-13巴特沃斯低通濾波器特點(diǎn)巴特沃斯低通濾波器特點(diǎn)：n在在 點(diǎn)，點(diǎn)， 的的 n ( n p ，隨，隨 ， n當(dāng)當(dāng) =0時(shí)，時(shí)， 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實(shí)際應(yīng)用中可以利于查表的方法求得系統(tǒng)函數(shù)：實(shí)際應(yīng)用中可以利于

11、查表的方法求得系統(tǒng)函數(shù)：1. 按按 求出濾波器階次求出濾波器階次N；2. 由由N和和 （或（或 p）查表得到歸一化濾波器分母多項(xiàng)式：查表得到歸一化濾波器分母多項(xiàng)式：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 的計(jì)算的計(jì)算 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幅值響應(yīng)在通帶和阻帶內(nèi)都是等波紋的幅值響應(yīng)在通帶和阻帶內(nèi)都是等波紋的n對(duì)

13、于給定的階數(shù)和給定的波紋要求，橢圓濾波器能獲得較對(duì)于給定的階數(shù)和給定的波紋要求，橢圓濾波器能獲得較其它濾波器更窄的過渡帶寬，就這點(diǎn)而言，橢圓濾波器是其它濾波器更窄的過渡帶寬，就這點(diǎn)而言，橢圓濾波器是最優(yōu)的。最優(yōu)的。n通帶和阻帶內(nèi)波紋固定時(shí)，階數(shù)越高，過渡帶越窄；通帶和阻帶內(nèi)波紋固定時(shí)，階數(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的定義的定義 與與 切比雪夫?yàn)V波器切比雪夫?yàn)V波器 相同相同482021-12-136.1.7模擬濾波器的比較模擬濾波器的比較n相同階數(shù)的頻率響應(yīng)特性相同階數(shù)的頻率響應(yīng)特性相同通帶最大衰減、阻帶最小衰減，巴特沃思濾波器的過渡帶相同通帶最大衰減、阻帶最小衰減，巴特沃思濾波器的過渡帶最寬；橢圓濾波器過渡帶最窄；最寬；橢圓濾波器過渡帶最窄；

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

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

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

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

19、一化低通原型濾波器：歸一化低通原型濾波器：歸一化復(fù)變量歸一化復(fù)變量 ，歸一化頻率，歸一化頻率通帶邊界頻率（截止頻率）通帶邊界頻率（截止頻率）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. 模擬高通濾波器的設(shè)計(jì)模擬高通濾波器的設(shè)計(jì)n低通原型到高通濾波器的映射變換：低通原型到高通濾波器的映射變換：在虛軸上的映射關(guān)系簡化為頻率變換關(guān)系，在虛軸上的映射關(guān)

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

21、相應(yīng)低通原型濾波器的指標(biāo)( )( )pphHPpsHsG p 4kHz,1kHz,0.1dB,40dBpspsff 24000421000pphss 1p 0.1dB40dBps 562021-12-132) 設(shè)計(jì)巴特沃斯低通濾波器。設(shè)計(jì)巴特沃斯低通濾波器。n 確定階次確定階次N：n可得可得Nn按照以上技術(shù)指標(biāo)設(shè)計(jì)巴特沃斯低通濾波器。按照以上技術(shù)指標(biāo)設(shè)計(jì)巴特沃斯低通濾波器。 /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另一設(shè)計(jì)途徑：另一設(shè)計(jì)途徑：模擬高通指標(biāo)模擬高通指標(biāo)模擬低通指標(biāo)模擬低通指標(biāo)(通帶截止頻率與高通一致通帶截止頻率與高通一致)模擬低通濾波器模擬低通濾波器設(shè)計(jì)歸一化模擬設(shè)計(jì)歸一化模擬低通濾波器低通濾波器G(p)模擬高通濾波器模擬高通濾波器2( )( )plh cqp 2pqs 592021-12-13MATLABwp=1;ws=4;Rp=0.1;As=40; %設(shè)置濾波器指標(biāo)參數(shù)設(shè)置濾波器指標(biāo)參數(shù)N,wc=buttord(wp,ws,Rp,As,s); %濾波器濾波器G(p)階數(shù)階數(shù)N和和3dB截止截止頻率頻率B,A=butter(N,wc,s); %計(jì)算低通濾波器計(jì)算低通

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

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

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

27、1.96333.1111slspWffffBf 120psdBdB 0puplfff 1p 642021-12-133）設(shè)計(jì)低通原型濾波器：）設(shè)計(jì)低通原型濾波器：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映射關(guān)系映射關(guān)系0,ppss 0,slsususlpuplpupl 672021-12-13減少減少 ，或增加，或增加減少減少 ，或增加，或增加帶阻濾波器的系統(tǒng)函數(shù)帶阻濾波器的系統(tǒng)函數(shù)n可以證明：可以證明：n如給定邊界頻率不滿足改條件，改變參數(shù)，提高指標(biāo)如給定邊界頻率不滿足改條件，改變參數(shù)，提高指標(biāo) ， 220( )( )WpB sBSpsHsG p 20plpuslsu plpuslsu plpuslsu pu su sl pl 6.2 IIR數(shù)字濾波器設(shè)計(jì)數(shù)字濾波器設(shè)計(jì)692021-12-136.2 IIR數(shù)字濾波器設(shè)計(jì)數(shù)字濾波器設(shè)計(jì)n目標(biāo)：滿足

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

30、標(biāo)轉(zhuǎn)換為：(),(),ppsspsspspsffF (),(),ppsspsff702021-12-13過渡模擬濾波器轉(zhuǎn)換為數(shù)字濾波器的要求：過渡模擬濾波器轉(zhuǎn)換為數(shù)字濾波器的要求：n保證因果穩(wěn)定性，保證因果穩(wěn)定性， Ha(s) 的因果穩(wěn)定性映射成的因果穩(wěn)定性映射成 H(z)后保持后保持不變，即不變，即 S 平面的平面的左半平面左半平面 ReS0 應(yīng)映射到應(yīng)映射到 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 正實(shí)軸正實(shí)軸 零頻零頻=0T 輻射線輻射線 角度角度: =0 實(shí)軸實(shí)軸 零頻零頻 =0 平行直線平行直線 頻率頻率:Z平面平面S平面平面/2s /22ssTT結(jié)論結(jié)論n穩(wěn)定性穩(wěn)定性如果模擬濾波器是因果穩(wěn)定的，則所有極點(diǎn)如果模擬濾波器是因果穩(wěn)定的，則所有極點(diǎn) 都在都在 S 左左半平面，那么變換后半平面，那么變換后 H(z) 的極點(diǎn)也都在單位圓以內(nèi)，的極點(diǎn)也都在單位圓以內(nèi)，因此數(shù)字濾波器保持因果穩(wěn)定。因此數(shù)字濾波器保持因果穩(wěn)定。n頻率響應(yīng)特性的近似性頻率響應(yīng)特性的近似性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混疊失真的影響混疊失真的影響相似性變差、無法滿足阻帶衰減指標(biāo)相似性變差、無法滿足阻帶衰減指標(biāo)n脈沖響應(yīng)不變法的優(yōu)缺點(diǎn)：脈沖響應(yīng)不變法的優(yōu)缺點(diǎn)：（1）頻率坐標(biāo)變換是線性的，數(shù)字濾

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

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

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

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

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

38、T=0.1s 的濾波器特性一致）。的濾波器特性一致）。 這是因?yàn)閿?shù)字濾波器指標(biāo)中有這是因?yàn)閿?shù)字濾波器指標(biāo)中有 ，轉(zhuǎn)換到過渡模擬，轉(zhuǎn)換到過渡模擬濾波器指標(biāo)中，濾波器指標(biāo)中，T 取不同值時(shí)一定有取不同值時(shí)一定有 ，模擬濾波器在模擬濾波器在1/2采樣頻率以上響應(yīng)特性一定都很小（阻采樣頻率以上響應(yīng)特性一定都很?。ㄗ鑾В僮儞Q到數(shù)字濾波器時(shí)引起的混疊都較小。帶），再變換到數(shù)字濾波器時(shí)引起的混疊都較小。s / 2sssTTF 862021-12-13小結(jié)：小結(jié)：( )( )()( )( )aaaH sh th nTh nzH z模模擬擬濾濾波波器器采采樣樣變變換換得得數(shù)數(shù)字字濾濾波波器器1( )Nkak

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

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

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

42、象，采用非線性頻率壓縮方法?？s方法。n雙線性變換法的基本設(shè)計(jì)思想雙線性變換法的基本設(shè)計(jì)思想脈沖響應(yīng)不變法：波形逼近脈沖響應(yīng)不變法：波形逼近雙線性變換法：算法逼近。雙線性變換法：算法逼近。用線性常系數(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 脈沖響應(yīng)不變法的主要缺點(diǎn)是頻響交疊產(chǎn)生的混淆，這脈沖響應(yīng)不變法的主要缺點(diǎn)是頻響交疊產(chǎn)生的混淆，這是從是從S平面到平

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

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

46、得到：1s Tze112 11zsTz952021-12-13雙線性變換的性能分析雙線性變換的性能分析n穩(wěn)定性穩(wěn)定性映射關(guān)系映射關(guā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非線性關(guān)系消除了頻率混疊，但會(huì)帶來幅度和相位失真。非線性關(guān)系消除了頻率混疊，但會(huì)帶來幅度和相位失真。/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頻率響應(yīng)的畸變頻率響應(yīng)的畸變n幅頻特性幅頻特性 相頻特性相頻特性0 0| )(|jeH| )(|jHa0 00 00 0)(jeH)( jHa0 00 0982021-12-13雙線性變換法設(shè)計(jì)數(shù)字低通濾波器的步驟雙線性變換法設(shè)計(jì)數(shù)字低通濾波器的步驟1）確定數(shù)字低通技術(shù)指標(biāo)：）確定數(shù)字低通技術(shù)指標(biāo)： 通帶邊界頻率通帶邊界頻率 、通帶衰減、通帶衰減 、 阻帶邊界頻率阻帶邊界頻率 、阻帶衰減、阻帶衰減 ；2）將數(shù)字低通指標(biāo)轉(zhuǎn)換成）將數(shù)字低通指標(biāo)轉(zhuǎn)換成過渡過渡模擬低通指標(biāo)：模擬低通指標(biāo)： 不

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