合工大 數(shù)字信號處理 第六章

1、1第一節(jié)第一節(jié) 數(shù)字濾波器的基本概念數(shù)字濾波器的基本概念一、數(shù)字濾波器基本概念一、數(shù)字濾波器基本概念數(shù)字濾波器數(shù)字濾波器: 輸入輸出均為數(shù)字信號，經(jīng)過一定運(yùn)算輸入輸出均為數(shù)字信號，經(jīng)過一定運(yùn)算 關(guān)系改變輸入信號所含頻率成分的相對關(guān)系改變輸入信號所含頻率成分的相對 比例或者濾除某些頻率成分的器件。比例或者濾除某些頻率成分的器件。優(yōu)點(diǎn)：優(yōu)點(diǎn)： 精度高，穩(wěn)定，體積小，重量輕，靈活，不精度高，穩(wěn)定，體積小，重量輕，靈活，不 要求阻抗匹配，能實(shí)現(xiàn)模擬濾波器要求阻抗匹配，能實(shí)現(xiàn)模擬濾波器(AF)無法無法 實(shí)現(xiàn)的特殊濾波功能。實(shí)現(xiàn)的特殊濾波功能。1. 數(shù)字濾波器的分類（回憶）數(shù)字濾波器的分類（回憶）2.數(shù)字

2、濾波器的技術(shù)要求數(shù)字濾波器的技術(shù)要求3. 數(shù)字濾波器設(shè)計(jì)方法概述數(shù)字濾波器設(shè)計(jì)方法概述21、數(shù)字濾波器的分類、數(shù)字濾波器的分類 經(jīng)典濾波器經(jīng)典濾波器：即一般濾波器（輸入信號中有用的頻即一般濾波器（輸入信號中有用的頻 率成分和希望濾除的頻率成分各占有率成分和希望濾除的頻率成分各占有 不同的頻帶，通過以合適的選頻濾波不同的頻帶，通過以合適的選頻濾波 器達(dá)到目的）器達(dá)到目的） 現(xiàn)代濾波器現(xiàn)代濾波器：如維納濾波器，卡爾曼濾波器，自適如維納濾波器，卡爾曼濾波器，自適 應(yīng)濾波器等最佳濾波器（按隨機(jī)信號應(yīng)濾波器等最佳濾波器（按隨機(jī)信號 內(nèi)部的統(tǒng)計(jì)分布規(guī)律，從干擾中最佳內(nèi)部的統(tǒng)計(jì)分布規(guī)律，從干擾中最佳 提取信

3、號）提取信號）3經(jīng)典濾波器的幾種類型：經(jīng)典濾波器的幾種類型：4 數(shù)字濾波器從實(shí)現(xiàn)的網(wǎng)絡(luò)結(jié)構(gòu)或者從單位脈沖響應(yīng)數(shù)字濾波器從實(shí)現(xiàn)的網(wǎng)絡(luò)結(jié)構(gòu)或者從單位脈沖響應(yīng)分類，可以分成無限脈沖響應(yīng)分類，可以分成無限脈沖響應(yīng)(IIR)濾波器和有限脈沖濾波器和有限脈沖響應(yīng)響應(yīng)(FIR)濾波器。它們的系統(tǒng)函數(shù)分別為：濾波器。它們的系統(tǒng)函數(shù)分別為：IIR濾波器濾波器 FIR濾波器濾波器 N1iiiM0iiiza1zb)z(H 1N0nnz )n(h)z(H5用一個(gè)因果穩(wěn)定的離散用一個(gè)因果穩(wěn)定的離散LSI系統(tǒng)的系統(tǒng)的系統(tǒng)函數(shù)系統(tǒng)函數(shù)H(z)逼近此性能指標(biāo)逼近此性能指標(biāo)按設(shè)計(jì)任務(wù)，按設(shè)計(jì)任務(wù)，確定濾波器性能要求確定濾波器性

4、能要求，制定技術(shù)指，制定技術(shù)指標(biāo)標(biāo)利用有限精度算法利用有限精度算法實(shí)現(xiàn)此系統(tǒng)函數(shù)實(shí)現(xiàn)此系統(tǒng)函數(shù)：如運(yùn)算結(jié)：如運(yùn)算結(jié)構(gòu)、字長的選擇等構(gòu)、字長的選擇等實(shí)際技術(shù)實(shí)現(xiàn)：軟件法、硬件法或?qū)嶋H技術(shù)實(shí)現(xiàn)：軟件法、硬件法或DSP芯片法芯片法2、數(shù)字濾波器的設(shè)計(jì)過程、數(shù)字濾波器的設(shè)計(jì)過程63、數(shù)字濾波器的技術(shù)要求數(shù)字濾波器的技術(shù)要求 我們通常用的數(shù)字濾波器一般屬于選頻濾波器。假我們通常用的數(shù)字濾波器一般屬于選頻濾波器。假設(shè)數(shù)字濾波器的傳輸函數(shù)設(shè)數(shù)字濾波器的傳輸函數(shù)H(e j)用下式表示：用下式表示：)( jj j e)H(e)H(e 幅頻特性幅頻特性|H(ej )|： 信號通過濾波器后的各頻率成分衰減情況。信號

5、通過濾波器后的各頻率成分衰減情況。相頻特性相頻特性 ( )：各頻率成分通過濾波器后在時(shí)間上的延時(shí)情況。各頻率成分通過濾波器后在時(shí)間上的延時(shí)情況。7n通帶：通帶： 1| )(|1|1 jpeHan阻帶：阻帶：2| )(|aeHjs n過渡帶：過渡帶：sp| ：阻帶截止頻率：阻帶截止頻率：通帶截止頻率：通帶截止頻率：通帶容限：通帶容限：阻帶容限：阻帶容限p s 1a2a8 通帶內(nèi)和阻帶內(nèi)允許的衰減一般用通帶內(nèi)和阻帶內(nèi)允許的衰減一般用dB數(shù)表示，通帶數(shù)表示，通帶內(nèi)允許的最大衰減用內(nèi)允許的最大衰減用 p表示，阻帶內(nèi)允許的最小衰減表示，阻帶內(nèi)允許的最小衰減用用 s表示，表示， p和和 s分別定義為：分別

6、定義為：dBeHdBeHeHppjjjp)(lg20)()(lg200 dBeHdBeHeHssjjjs)(lg20)()(lg200 通帶截止頻率。通帶截止頻率。為為，為為我們稱此時(shí)的我們稱此時(shí)的，時(shí)，時(shí)，當(dāng)當(dāng)3dBdB 0.707eHccppjp 322)( 上式中，上式中， 歸一化為歸一化為1。| )(|0jeH 9 H(z)H(z-1) 的極點(diǎn)既是共軛的，的極點(diǎn)既是共軛的，又是以單位圓成鏡像對稱的。又是以單位圓成鏡像對稱的。Re zIm jz01aa*a*1/aH(z)的極點(diǎn)的極點(diǎn)：單位圓內(nèi)的極點(diǎn)：單位圓內(nèi)的極點(diǎn)4、表征濾波器頻率響應(yīng)的特征參量、表征濾波器頻率響應(yīng)的特征參量)(e)HH

7、(e)H(ej j 2j jez1-j j )z)H(zH()H(eH(e 幅度平方響應(yīng)幅度平方響應(yīng)10 相位響應(yīng)相位響應(yīng)ImRe)()(eHj)H(ee)H(e)H(ej j jj j ReIm)()H(e)(eHarctanj j 相相位位響響應(yīng)應(yīng)：)( jj j e)H(e)(eH )( j2j j e)(eH)H(e j ez1j -j j j )H(zH(z)j21)H(e)H(ej21)(eH)H(ej21 lnlnln)(11相位對角頻率的導(dǎo)數(shù)的負(fù)值相位對角頻率的導(dǎo)數(shù)的負(fù)值 群延遲響應(yīng)群延遲響應(yīng) jezzH1dzzdHzdd)( )()(Re)(若濾波器通帶內(nèi)若濾波器通帶內(nèi) 為常

8、數(shù)，為常數(shù)， 則為則為線性相位濾波器線性相位濾波器。)( 12 用一因果穩(wěn)定的離散用一因果穩(wěn)定的離散LSI系統(tǒng)逼近給定的性能要求：系統(tǒng)逼近給定的性能要求：s s平面逼近：模擬濾波器的設(shè)計(jì)平面逼近：模擬濾波器的設(shè)計(jì)z z平面逼近：數(shù)字濾波器的設(shè)計(jì)平面逼近：數(shù)字濾波器的設(shè)計(jì)5、IIR數(shù)字濾波器的設(shè)計(jì)方法數(shù)字濾波器的設(shè)計(jì)方法 先設(shè)計(jì)模擬濾波器，再轉(zhuǎn)換為數(shù)字濾波器先設(shè)計(jì)模擬濾波器，再轉(zhuǎn)換為數(shù)字濾波器 N1iiiM0iiiza1zb)z(H 即為求濾波器的各系數(shù)：即為求濾波器的各系數(shù)：iiba 和和13第二節(jié)第二節(jié) 模擬濾波器的設(shè)計(jì)模擬濾波器的設(shè)計(jì) 模擬濾波器的理論和設(shè)計(jì)方法已發(fā)展得相當(dāng)成模擬濾波器的理

9、論和設(shè)計(jì)方法已發(fā)展得相當(dāng)成熟，且有若干典型的模擬濾波器供我們選擇，如：熟，且有若干典型的模擬濾波器供我們選擇，如： 1) 巴特沃斯巴特沃斯(Butterworth)濾波器濾波器 2) 切比雪夫切比雪夫(Chebyshev)濾波器濾波器 3) 橢圓橢圓(Ellipse)濾波器濾波器 4) 貝塞爾貝塞爾(Bessel)濾波器濾波器 這些濾波器都有嚴(yán)格的設(shè)計(jì)公式、現(xiàn)成的曲線這些濾波器都有嚴(yán)格的設(shè)計(jì)公式、現(xiàn)成的曲線和圖表供設(shè)計(jì)人員使用。和圖表供設(shè)計(jì)人員使用。141、模擬低通濾波器的設(shè)計(jì)指標(biāo)及逼近方法、模擬低通濾波器的設(shè)計(jì)指標(biāo)及逼近方法模擬低通濾波器的設(shè)計(jì)指標(biāo)有模擬低通濾波器的設(shè)計(jì)指標(biāo)有 p和和s分別稱

10、為通帶截止頻率和阻帶截止頻率，分別稱為通帶截止頻率和阻帶截止頻率， p是通帶是通帶(0p)中的最大衰減系數(shù)中的最大衰減系數(shù) s是阻帶是阻帶s的最小衰減系數(shù)的最小衰減系數(shù)dBeH20dBeHeH20pp0jajajap)(lg)()(lg dBeH20dBeHeH20ss0jajajas)(lg)()(lg 15c稱為稱為3dB截止頻率：截止頻率：707022jHca.)( dB3jH20ca )(lg 16 濾波器的技術(shù)指標(biāo)給定后，需要設(shè)計(jì)一個(gè)傳輸函濾波器的技術(shù)指標(biāo)給定后，需要設(shè)計(jì)一個(gè)傳輸函數(shù)數(shù)Ha(s)，希望其希望其幅度平方函數(shù)幅度平方函數(shù)滿足給定的指標(biāo)滿足給定的指標(biāo) p和和 s，一般濾波器

11、的單位沖激響應(yīng)為實(shí)數(shù)，有：一般濾波器的單位沖激響應(yīng)為實(shí)數(shù)，有： jsaaaa2asHsHjHjHjH )()()()()(*172、巴特沃斯低通濾波器的設(shè)計(jì)方法、巴特沃斯低通濾波器的設(shè)計(jì)方法(1) 巴特沃斯低通濾波器的幅度平方函數(shù)巴特沃斯低通濾波器的幅度平方函數(shù)|Ha(j)|2 用下式表示：用下式表示：NcajH2211)( N越大，越接近理想濾波器，越大，越接近理想濾波器，N越大，濾波器的實(shí)現(xiàn)也越復(fù)雜。越大，濾波器的實(shí)現(xiàn)也越復(fù)雜。18(2)幅度平方函數(shù)極點(diǎn)分布及幅度平方函數(shù)極點(diǎn)分布及Ha(s)的構(gòu)成的構(gòu)成將幅度平方函數(shù)將幅度平方函數(shù)|Ha(j)|2寫成寫成s的函數(shù)：的函數(shù)：21( )()1(

12、)aaNcHs Hssj上式表明，極點(diǎn)上式表明，極點(diǎn)sk用下式表示：用下式表示：1121()222( 1)()kjNNkccsje 為形成穩(wěn)定的濾波器，為形成穩(wěn)定的濾波器，2N個(gè)極點(diǎn)中個(gè)極點(diǎn)中只取只取s平面左半平面的平面左半平面的N個(gè)極點(diǎn)構(gòu)成個(gè)極點(diǎn)構(gòu)成Ha(s)，而右半平面的而右半平面的N個(gè)極點(diǎn)構(gòu)成個(gè)極點(diǎn)構(gòu)成Ha(-s)。Ha(s)的表示式為的表示式為10( )()NcaNkkHsss19設(shè)設(shè)N=3，極點(diǎn)有極點(diǎn)有6個(gè)，它們分別為個(gè)，它們分別為23012321334135jccjcjccjcsessesesse 取取s平面左半平面的極點(diǎn)平面左半平面的極點(diǎn)s0,s1,s2組成組成Ha(s)：)()

13、()(32323 jcjcccaesesssH 20(3) 頻率歸一化頻率歸一化 由于各濾波器的幅頻特性不同，為使設(shè)計(jì)統(tǒng)一，將所有由于各濾波器的幅頻特性不同，為使設(shè)計(jì)統(tǒng)一，將所有的頻率歸一化。這里采用對的頻率歸一化。這里采用對3dB截止頻率截止頻率c歸一化，歸一歸一化，歸一化后的化后的Ha(s)表示為表示為 101( )()aNkkccHsss式中，式中，s/c=j/c。 令令=/c，稱為歸一化頻率稱為歸一化頻率；令；令 p = j= j/c ， p 稱為歸一化復(fù)變量，這樣歸一化巴特沃斯的傳輸函數(shù)為稱為歸一化復(fù)變量，這樣歸一化巴特沃斯的傳輸函數(shù)為101( )()aNkkHppp21式中，式中，

14、pk為歸一化極點(diǎn)，用下式表示：為歸一化極點(diǎn)，用下式表示：1 21()22,0,1,1kjNkpekN帶入帶入Ha(p)表達(dá)式，得到的表達(dá)式，得到的Ha(p)的分母是的分母是p的的N階多項(xiàng)式，階多項(xiàng)式，用下式表示：用下式表示：NNNappbpbbpH11101)(歸一化的傳輸函數(shù)系數(shù)歸一化的傳輸函數(shù)系數(shù)Ha (p) 的系數(shù)以及極點(diǎn)可以查表得到。的系數(shù)以及極點(diǎn)可以查表得到。22表表6.2.1 巴特沃斯歸一化低通濾波器參數(shù)巴特沃斯歸一化低通濾波器參數(shù) 232425(4)階數(shù)階數(shù)N的確定的確定 階數(shù)階數(shù)N的大小主要影響幅度特性的大小主要影響幅度特性下降的速度下降的速度，它應(yīng)該，它應(yīng)該由技術(shù)指標(biāo)確定。將

15、由技術(shù)指標(biāo)確定。將 = p代入幅度平方函數(shù)中代入幅度平方函數(shù)中:2)(lg10)(lg20ppjapjapeHeH NcppajH2211)( 102101pNcp 26將將 = s代入幅度平方函數(shù)中代入幅度平方函數(shù)中:2)(lg10)(lg20ssjasjaseHeH NcssajH2211)( 102101sNcs 27102101pNcp 102101sNcs 1101101010 spNsp ,pssp 令：令：110110k11010spsp spspkN lglg 用上式求出的用上式求出的N可能有小數(shù)部分，應(yīng)取大于等于可能有小數(shù)部分，應(yīng)取大于等于N的最小整數(shù)。的最小整數(shù)。28 關(guān)于

16、關(guān)于3dB截止頻率截止頻率c，如果技術(shù)指標(biāo)中沒有給出，如果技術(shù)指標(biāo)中沒有給出，可以按照下面兩式求出：可以按照下面兩式求出：102101pNcp 102101sNcs Npcp2110)110( Nscs2110)110( 通常是用一個(gè)算出通常是用一個(gè)算出c，然后用另一個(gè)然后用另一個(gè)(反過來反過來)來檢來檢驗(yàn)。驗(yàn)。29總結(jié)以上，低通巴特沃斯濾波器的設(shè)計(jì)步驟如下：總結(jié)以上，低通巴特沃斯濾波器的設(shè)計(jì)步驟如下： 根據(jù)技術(shù)指標(biāo)根據(jù)技術(shù)指標(biāo)p, p,s, s，求出濾波器的階數(shù)求出濾波器的階數(shù)N。(2) 求出歸一化極點(diǎn)求出歸一化極點(diǎn)pk，得到歸一化傳輸函數(shù)得到歸一化傳輸函數(shù)Ha(p)。 (3) 將將Ha(p

17、)去歸一化。將去歸一化。將p=s/c代入代入Ha(p)，得到實(shí)際的濾波器得到實(shí)際的濾波器 傳輸函數(shù)傳輸函數(shù)Ha(s)。1 21()22,0,1,1kjNkpekNNNNappbpbbpH11101)(spspkN lglg cspaapHsH )()(30例：例： 已知通帶截止頻率已知通帶截止頻率fp=5kHz，通帶最大衰減通帶最大衰減 p=2dB，阻帶阻帶截止頻率截止頻率fs=12kHz，阻帶最小衰減阻帶最小衰減 s=30dB，按照以上技術(shù)指按照以上技術(shù)指標(biāo)設(shè)計(jì)巴特沃斯低通濾波器。標(biāo)設(shè)計(jì)巴特沃斯低通濾波器。解：解： (1) 確定階數(shù)確定階數(shù)N：0.10.11010.024210122.42l

18、g0.02424.25,5lg2.4psaspassppkffNN 31(2) 求極點(diǎn)：求極點(diǎn)：3455016523754,jjjjjsesesesese歸一化傳輸函數(shù)為歸一化傳輸函數(shù)為401( )()akkHppp32 上式分母可以展開成為五階多項(xiàng)式，或者將共軛極點(diǎn)放在上式分母可以展開成為五階多項(xiàng)式，或者將共軛極點(diǎn)放在一起，形成因式分解形式。不如直接查表簡單，由一起，形成因式分解形式。不如直接查表簡單，由N=5，直接直接查表得到：查表得到：極點(diǎn)：極點(diǎn)：-0.3090j0.9511、-0.8090j0.5878、 -1.00005432432101( )aHppb pb pb pb pb其中：

19、其中： b0=1.0000,b1=3.2361,b2=5.2361,b3=5.2361,b4=3.236133(3) 為將為將Ha(p)去歸一化，先求去歸一化，先求3dB截止頻率截止頻率c。)/(2755. 52)110(2110sradkNpcp 將將p=s/c代入代入Ha(p)中得到：中得到：554233245432( )10cacccccHssbsbsbsbsb)/(525.102)110(2110sradkNcss 檢驗(yàn)：檢驗(yàn)： 可以看出，滿足可以看出，滿足 s=30dB 的真實(shí)的真實(shí)fs在在10.525kHz處，與處，與12kHz比，還有富裕量。比，還有富裕量。344、模擬濾波器的頻

20、率變換、模擬濾波器的頻率變換 模擬高通、帶通、帶阻濾波器的設(shè)計(jì)模擬高通、帶通、帶阻濾波器的設(shè)計(jì)(1) 低通到高通的頻率變換低通到高通的頻率變換和和之間的關(guān)系為之間的關(guān)系為 上式即是低通到高通的頻率變換公式，如果已上式即是低通到高通的頻率變換公式，如果已知低通知低通G(j)，高通高通H(j)則用下式轉(zhuǎn)換：則用下式轉(zhuǎn)換： 1 1)()( jGjH35 1 1)()( jGjH36 模擬高通濾波器的模擬高通濾波器的設(shè)計(jì)設(shè)計(jì)步驟如下：步驟如下：A、確定確定高通濾波器的技術(shù)指標(biāo)高通濾波器的技術(shù)指標(biāo)：通帶下限頻率：通帶下限頻率p， 阻帶上限頻率阻帶上限頻率s，通帶最大衰減通帶最大衰減 p，阻帶最小衰阻帶最

21、小衰 減減 s。B、確定相應(yīng)低通濾波器的設(shè)計(jì)指標(biāo)，確定相應(yīng)低通濾波器的設(shè)計(jì)指標(biāo)，將高通濾波器將高通濾波器 的邊界頻率轉(zhuǎn)換成低通濾波器的邊界頻率的邊界頻率轉(zhuǎn)換成低通濾波器的邊界頻率，各項(xiàng)，各項(xiàng) 設(shè)計(jì)指標(biāo)為：設(shè)計(jì)指標(biāo)為：低通濾波器通帶截止頻率低通濾波器通帶截止頻率p=1/p；低通濾波器阻帶截止頻率低通濾波器阻帶截止頻率s=1/s；通帶最大衰減仍為通帶最大衰減仍為 p，阻帶最小衰減仍為阻帶最小衰減仍為 s。37C、設(shè)計(jì)歸一化低通濾波器設(shè)計(jì)歸一化低通濾波器G(p)。D、求模擬高通的求模擬高通的H(s)。將將G(p) 先轉(zhuǎn)換成歸一化先轉(zhuǎn)換成歸一化高通高通H(q)，再去歸一化，得再去歸一化，得H(s)：s

22、pcpGsH )()(cspaapHsH )()(回憶低通濾波器的去歸一化：回憶低通濾波器的去歸一化：38例：例： 設(shè)計(jì)高通濾波器，設(shè)計(jì)高通濾波器，fp=200Hz、fs=100Hz，幅度特性單調(diào)下幅度特性單調(diào)下 降，降，fp處最大衰減為處最大衰減為3dB，阻帶最小衰減阻帶最小衰減 s=15dB。解：解： 高通技術(shù)要求高通技術(shù)要求fp=200Hz p=3dBfs=100Hz, s=15dB歸一化頻率歸一化頻率5 . 0, 1 csscppffff 低通技術(shù)要求：低通技術(shù)要求：dBdBspssp153211 39 設(shè)計(jì)歸一化低通設(shè)計(jì)歸一化低通G(p)。采用巴特沃斯濾波器，故采用巴特沃斯濾波器，故

23、0.10.1321010.181012lg2.47,3lg1( )221psspssppspspkkNNG pppp 求模擬高通求模擬高通H(s)：33223( )( )222cpcccscpsH sG psssf 40(2) 低通到帶通的頻率變換低通到帶通的頻率變換低通與帶通濾波器的幅度特性如圖：低通與帶通濾波器的幅度特性如圖：。限頻率和通帶下限頻率為帶通濾波器的通帶上和l ulu0luB 為通帶中心頻率，有：為通帶中心頻率，有：為通帶帶寬，為通帶帶寬，令：令：20 阻帶下限頻率。為下阻帶上限頻率和上和s21 s41uluulls2ss1sBBBB 2021，令：令：由由與與的對應(yīng)關(guān)系，得到

24、：的對應(yīng)關(guān)系，得到： 202 由上表知由上表知p對應(yīng)對應(yīng)u，代入上式中，有代入上式中，有1202 Bluluuup 42 上式稱為低通到帶通的頻率變換公式。利用該式將帶通的上式稱為低通到帶通的頻率變換公式。利用該式將帶通的邊界頻率轉(zhuǎn)換成低通的邊界頻率。邊界頻率轉(zhuǎn)換成低通的邊界頻率。下面推導(dǎo)由歸一化低通到帶通的轉(zhuǎn)換公式：下面推導(dǎo)由歸一化低通到帶通的轉(zhuǎn)換公式： jp 202 jp 202 將將q=j代入上式，得到：代入上式，得到：qqjjjjp20220222202 43)(22202022luulsssBBsBsBsp 為去歸一化，將為去歸一化，將q=s/B代入代入 ，得到：，得到：qqp202

25、 )(2)()(luulssppGsH 因此有：因此有：44 模擬帶通的設(shè)計(jì)步驟模擬帶通的設(shè)計(jì)步驟A、 確定模擬帶通濾波器的技術(shù)指標(biāo)，即：確定模擬帶通濾波器的技術(shù)指標(biāo)，即：帶通上限頻率帶通上限頻率u，帶通下限頻率帶通下限頻率l下阻帶上限頻率下阻帶上限頻率s1 ,上阻帶下限頻率上阻帶下限頻率s2 通帶中心頻率通帶中心頻率02=lu，通帶寬度通帶寬度B=u-l與以上邊界頻率對應(yīng)的歸一化邊界頻率如下：與以上邊界頻率對應(yīng)的歸一化邊界頻率如下：121220,sslssluuluBBBB45B、 確定歸一化低通技術(shù)要求：確定歸一化低通技術(shù)要求：22222010211,sspssss s與與-s的絕對值可能

26、不相等，一般取絕對值小的的絕對值可能不相等，一般取絕對值小的s，這樣保證在較大的這樣保證在較大的s處更能滿足要求。處更能滿足要求。通帶最大衰減仍為通帶最大衰減仍為 p，阻帶最小衰減亦為阻帶最小衰減亦為 s。C、設(shè)計(jì)歸一化低通設(shè)計(jì)歸一化低通G(p)。D、 將將G(p)轉(zhuǎn)換成帶通轉(zhuǎn)換成帶通H(s)。46例：設(shè)計(jì)模擬帶通濾波器，通帶帶寬例：設(shè)計(jì)模擬帶通濾波器，通帶帶寬B=2200rad/s，中中 心頻率心頻率0=21000rad/s，通帶內(nèi)最大衰減通帶內(nèi)最大衰減 p=3dB， 阻帶阻帶s1=2830rad/s,s2=21200rad/s，阻帶最阻帶最 小衰減小衰減 s=15dB。解解 ：(1) 模擬

27、帶通的技術(shù)要求模擬帶通的技術(shù)要求0=21000rad/s, p=3dBs1 =2830rad/s,s2=21200rad/s, s=15dB B = 2200rad/s;0= 0 /B = 5,s1= s1 /B = 4.15, s2 = s2 /B = 647(2) 模擬歸一化低通技術(shù)要求：模擬歸一化低通技術(shù)要求：874. 1,833. 1, 11202122022 ssssssp 取取s=1.833, p=3dB, s=15dB。(3)設(shè)計(jì)模擬歸一化低通濾波器設(shè)計(jì)模擬歸一化低通濾波器G(p)：0.10.11010.181011.833lg2.83lgpsspssppspspkkN 48取取

28、N=3，查表得：查表得：232()1()221()()luulspsGppppHsGp(4) 求模擬帶通求模擬帶通H(s)：)(2)()(luulssppGsH 160402220403320422056332)23()4()23(2)( BssBsBBsBBssBssH49(3) 低通到帶阻的頻率變換低通到帶阻的頻率變換低通與帶阻濾波器的幅頻特性如圖低通與帶阻濾波器的幅頻特性如圖50 l和和u分別是下通帶截止頻率和上通帶截止頻率，分別是下通帶截止頻率和上通帶截止頻率， s1和和s2分別為阻帶的下限頻率和上限頻率，分別為阻帶的下限頻率和上限頻率， 0為阻帶中心頻率，為阻帶中心頻率，02=ul，

29、 阻帶帶寬阻帶帶寬B=u-l，B作為歸一化參考頻率。作為歸一化參考頻率。 相應(yīng)的歸一化邊界頻率為：相應(yīng)的歸一化邊界頻率為：u=u/B, l=l/B,s1=s1/B, s2=s2/B;02=ul 51根據(jù)根據(jù)與與的對應(yīng)關(guān)系，可得到：的對應(yīng)關(guān)系，可得到： 且且u -l=1，p=1，上式稱為低通到帶阻的頻率變換公式。上式稱為低通到帶阻的頻率變換公式。將上式代入將上式代入p=j，并去歸一化，可得：并去歸一化，可得：上式就是直接由歸一化低通轉(zhuǎn)換成帶阻的頻率變換公式。上式就是直接由歸一化低通轉(zhuǎn)換成帶阻的頻率變換公式。202 ullussssBp 220)(ullussssBppGsH 2202)()()(

30、52 設(shè)計(jì)帶阻濾波器的步驟設(shè)計(jì)帶阻濾波器的步驟A、確定模擬帶阻濾波器的技術(shù)要求：確定模擬帶阻濾波器的技術(shù)要求：下通帶截止頻率下通帶截止頻率l，上通帶截止頻率上通帶截止頻率u阻帶下限頻率阻帶下限頻率s1，阻帶上限頻率阻帶上限頻率s2阻帶中心頻率阻帶中心頻率02=ul，阻帶寬度阻帶寬度B=u-l它們相應(yīng)的歸一化邊界頻率為它們相應(yīng)的歸一化邊界頻率為 l=l/B, u=u/B, s1=s1/B, s2=s2/B, 20=ul以及通帶最大衰減以及通帶最大衰減 p和阻帶最小衰減和阻帶最小衰減 s。53B、 確定歸一化模擬低通技術(shù)要求，確定歸一化模擬低通技術(shù)要求，12222210201,sspssss 取取

31、s和和-s的絕對值較小的的絕對值較小的s；通帶最大衰減為通帶最大衰減為 p，阻帶最阻帶最小衰減為小衰減為 s。C、設(shè)計(jì)歸一化模擬低通設(shè)計(jì)歸一化模擬低通G(p)D、將將G(p)轉(zhuǎn)換成帶阻濾波器轉(zhuǎn)換成帶阻濾波器H(s)ullussssBppGsH 2202)()()(54例：設(shè)計(jì)模擬帶阻濾波器，其技術(shù)要求為：例：設(shè)計(jì)模擬帶阻濾波器，其技術(shù)要求為： l=2905rad/s, s1=2980rad/s, s2= 21020rad/s, u=21105rad/s, p=3dB, s=25dB。試設(shè)計(jì)巴特沃斯帶阻濾波器。試設(shè)計(jì)巴特沃斯帶阻濾波器。解解: (1) 模擬帶阻濾波器的技術(shù)要求：模擬帶阻濾波器的技

32、術(shù)要求： l=2905,u=21105; s1=2980,s2=21020; 02=lu=421000025, B=u-l=2200;55l=l/B=4.525, u=u/B=5.525;s1=s1/B=4.9, s2=5.1;20=lu=25(2) 歸一化低通的技術(shù)要求：歸一化低通的技術(shù)要求：95. 4,05. 5, 12021120222 sssssspdBdBsps25395. 4 ，取取：56(3)設(shè)計(jì)歸一化低通濾波器設(shè)計(jì)歸一化低通濾波器G(p)：0.10.121010.05621014.95lg1.8,2lg1( )21psspssppspspkkNNG ppp (4) 帶阻濾波器的

33、帶阻濾波器的H(s)為：為：40202202244022042)2(22)()(202 sBsBBsssspGsHssBp57第三節(jié)第三節(jié) 用脈沖響應(yīng)不變法用脈沖響應(yīng)不變法 設(shè)計(jì)設(shè)計(jì)IIR數(shù)字低通濾波器數(shù)字低通濾波器 利用模擬濾波器來設(shè)計(jì)數(shù)字濾波器，就是從已利用模擬濾波器來設(shè)計(jì)數(shù)字濾波器，就是從已知的模擬濾波器傳遞函數(shù)知的模擬濾波器傳遞函數(shù)Ha(s)設(shè)計(jì)數(shù)字濾波器的系設(shè)計(jì)數(shù)字濾波器的系統(tǒng)函數(shù)統(tǒng)函數(shù)H(z)。 因此，它歸根結(jié)底是一個(gè)由因此，它歸根結(jié)底是一個(gè)由S平面映射到平面映射到Z平面平面的變換，這個(gè)變換通常是復(fù)變函數(shù)的映射變換，為的變換，這個(gè)變換通常是復(fù)變函數(shù)的映射變換，為了保證了保證轉(zhuǎn)換后的

34、轉(zhuǎn)換后的H(z)穩(wěn)定且滿足技術(shù)要求穩(wěn)定且滿足技術(shù)要求，這個(gè)映，這個(gè)映射變換必須滿足以下兩條基本要求：射變換必須滿足以下兩條基本要求：58(2)數(shù)字濾波器的頻率響應(yīng)模仿模擬濾波器的頻數(shù)字濾波器的頻率響應(yīng)模仿模擬濾波器的頻 響，響，s平面的虛軸映射為平面的虛軸映射為z平面的單位圓，相應(yīng)平面的單位圓，相應(yīng)的頻率之間成線性關(guān)系。的頻率之間成線性關(guān)系。 因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍 是因果穩(wěn)定的。即是因果穩(wěn)定的。即 S平面的左半平面必須映射平面的左半平面必須映射到到Z平面單位圓的內(nèi)部平面單位圓的內(nèi)部。59一、脈沖響應(yīng)不變法的轉(zhuǎn)換原理一、脈沖響應(yīng)不變法的

35、轉(zhuǎn)換原理 通過對連續(xù)函數(shù)通過對連續(xù)函數(shù)ha(t)等間隔采樣得到離散序列等間隔采樣得到離散序列ha(nT)。令令h(n)= ha(nT) ，T為采樣間隔。為采樣間隔。它是一它是一種時(shí)域上的轉(zhuǎn)換方法。種時(shí)域上的轉(zhuǎn)換方法。)()()()()(zHnhnThthsHZaaa 變換變換等間隔采樣等間隔采樣拉氏逆變換拉氏逆變換60設(shè)模擬濾波器的傳輸函數(shù)為設(shè)模擬濾波器的傳輸函數(shù)為Ha(s),相應(yīng)的單位沖激響應(yīng)是相應(yīng)的單位沖激響應(yīng)是ha(t)( )( )aaHsLT h t 設(shè)模擬濾波器設(shè)模擬濾波器Ha(s)只有單階極點(diǎn)，且分母多項(xiàng)式的階只有單階極點(diǎn)，且分母多項(xiàng)式的階次高于分子多項(xiàng)式的階次，將次高于分子多項(xiàng)式

36、的階次，將Ha(s)用部分分式表示：用部分分式表示： 1( )NiaiiAHsss式中式中si為為Ha(s)的單階極點(diǎn)。將的單階極點(diǎn)。將Ha(s)進(jìn)行逆拉氏變換得到進(jìn)行逆拉氏變換得到ha(t)：NitsiatueAthi1)()(61 式中式中u(t)是單位階躍函數(shù)。對是單位階躍函數(shù)。對ha(t)進(jìn)行等間隔采樣，采樣進(jìn)行等間隔采樣，采樣間隔為間隔為T，得到：得到：1( )()()iNs nTaiih nh nTAeu nT對上式進(jìn)行對上式進(jìn)行Z變換，得到數(shù)字濾波器的系統(tǒng)函數(shù)變換，得到數(shù)字濾波器的系統(tǒng)函數(shù)H(z)： 11( )1iNis TiAH zez11 zeAssATsiiii62 一般一

37、般Ha(s)的極點(diǎn)的極點(diǎn)si是一個(gè)復(fù)數(shù)，且以共軛成對的形式出現(xiàn)，是一個(gè)復(fù)數(shù)，且以共軛成對的形式出現(xiàn)，將一對復(fù)數(shù)共軛極點(diǎn)放在一起，可形成一個(gè)二階基本節(jié)。如果將一對復(fù)數(shù)共軛極點(diǎn)放在一起，可形成一個(gè)二階基本節(jié)。如果模擬濾波器的二階基本節(jié)的形式為：模擬濾波器的二階基本節(jié)的形式為： 可以推導(dǎo)出相應(yīng)的數(shù)字濾波器二階基本節(jié)可以推導(dǎo)出相應(yīng)的數(shù)字濾波器二階基本節(jié)(只有實(shí)數(shù)乘法只有實(shí)數(shù)乘法)的形式為：的形式為：1111121211cos12cosTTTz eTz eTz e1112211()sjs 極點(diǎn)為極點(diǎn)為63如果模擬濾波器二階基本節(jié)的形式為：如果模擬濾波器二階基本節(jié)的形式為： 111111221111212

38、1,()s i n12c o sTTTjszeTzeTze 1111112211112121,()s i n12c o sTTTjszeTzeTze極點(diǎn)為極點(diǎn)為 相應(yīng)的數(shù)字濾波器二階基本節(jié)的形式為：相應(yīng)的數(shù)字濾波器二階基本節(jié)的形式為：64( )( ) ()aanhth ttnT設(shè)設(shè)ha(t)的采樣信號用的采樣信號用 表示，表示，)(tha二、存在的問題二、存在的問題拉氏變換：拉氏變換：snTnastnastaaenThdtenTtthdtethsH )()()()()( sTsTezeznnazHznh )()(65 上式表示采樣信號的拉氏變換與相應(yīng)的序列的上式表示采樣信號的拉氏變換與相應(yīng)的序

39、列的Z變換之間變換之間的映射關(guān)系可用下式表示：的映射關(guān)系可用下式表示：sTez 將將s=j代入上式，得代入上式，得 我們知道模擬信號我們知道模擬信號ha(t)的傅里葉變換的傅里葉變換Ha(j)和其采樣信和其采樣信號號 的傅里葉變換的傅里葉變換 之間的關(guān)系滿足：之間的關(guān)系滿足：)(tha)( jHa ksaajkjHTjH)(1)(sTezksaazHjksHTsH )()(1)(66sTez 由由： jrezjsTjTjeere 得得到到： TerT 67 j S平面平面Z平面平面RezjImz =0 (s平面虛軸平面虛軸) r=1 (z平面單位圓平面單位圓) 0 (s左半平面左半平面) r0

40、 (s右半平面右半平面) r1 (z平面單位圓外部平面單位圓外部) r與與 的關(guān)系的關(guān)系68 與與 的關(guān)系：的關(guān)系： = T =0 (s平面實(shí)軸平面實(shí)軸) =0 (z平面正實(shí)軸平面正實(shí)軸) = 0 = 0T (s平面平行于實(shí)軸的直線平面平行于實(shí)軸的直線) (z平面始于原點(diǎn)，輻角為平面始于原點(diǎn)，輻角為 0T的輻射線的輻射線) ：從從 - /T /T ：從：從 - (s平面為平面為2 /T的一個(gè)水平帶的一個(gè)水平帶) (z平面輻角轉(zhuǎn)了一周，覆蓋整個(gè)平面輻角轉(zhuǎn)了一周，覆蓋整個(gè)z平面平面) Z平面平面jImzRez j S平面平面- /T /T3 /T-3 /T多值映射多值映射 69 存在的問題：混疊失

41、真存在的問題：混疊失真可可得得：由由： ksajksHTzH)(1)( TkjHTeHkaj 21)(數(shù)字濾波器的頻率響應(yīng)是模擬濾波器頻率響應(yīng)的周期延拓。數(shù)字濾波器的頻率響應(yīng)是模擬濾波器頻率響應(yīng)的周期延拓。 只有當(dāng)模擬濾波器的頻率響應(yīng)是限帶的，且?guī)抻谡郫B只有當(dāng)模擬濾波器的頻率響應(yīng)是限帶的，且?guī)抻谡郫B頻率以內(nèi)時(shí)，才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)頻率以內(nèi)時(shí)，才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應(yīng)，而不產(chǎn)生混疊失真。即：重現(xiàn)模擬濾波器的頻率響應(yīng)，而不產(chǎn)生混疊失真。即： TjHTeHaj 1)(|70 但是，任何一個(gè)實(shí)際的模擬濾波器頻率響應(yīng)都不但是，任何一個(gè)實(shí)際的

42、模擬濾波器頻率響應(yīng)都不是嚴(yán)格限帶的（非理想），是嚴(yán)格限帶的（非理想）， 變換后就會產(chǎn)生周期延變換后就會產(chǎn)生周期延拓分量的頻譜交疊，即產(chǎn)生頻率響應(yīng)的混疊失真。這拓分量的頻譜交疊，即產(chǎn)生頻率響應(yīng)的混疊失真。這時(shí)數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，時(shí)數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，而是有一定的失真。當(dāng)模擬濾波器的頻率響應(yīng)在折疊而是有一定的失真。當(dāng)模擬濾波器的頻率響應(yīng)在折疊頻率以上處衰減越大、越快時(shí)，變換后頻率響應(yīng)混疊頻率以上處衰減越大、越快時(shí)，變換后頻率響應(yīng)混疊失真就越小。這時(shí)，采用脈沖響應(yīng)不變法設(shè)計(jì)的數(shù)字失真就越小。這時(shí)，采用脈沖響應(yīng)不變法設(shè)計(jì)的數(shù)字濾波器才能得到良好的效果。濾

43、波器才能得到良好的效果。71脈沖響應(yīng)不變法的頻率混疊現(xiàn)象脈沖響應(yīng)不變法的頻率混疊現(xiàn)象72數(shù)字濾波器頻率響應(yīng)幅度還與采樣間隔數(shù)字濾波器頻率響應(yīng)幅度還與采樣間隔T成反比：成反比： TjHTeHaj 1)(| 如果采樣頻率很高，即如果采樣頻率很高，即T很小，數(shù)字濾波器可能具有太很小，數(shù)字濾波器可能具有太高的增益，這是不希望的。為了使數(shù)字濾波器增益不隨采高的增益，這是不希望的。為了使數(shù)字濾波器增益不隨采樣頻率而變化，可以作以下簡單的修正，令：樣頻率而變化，可以作以下簡單的修正，令：h(n)=Tha(nT) 則有則有: NkTskzeTAzHk111)( TjHkTjTjHeHakaj 2)(很小的數(shù)在

44、分子很小的數(shù)在分子比在分母安全！比在分母安全！73例：例： 已知模擬濾波器的傳輸函數(shù)已知模擬濾波器的傳輸函數(shù)Ha(s)為：為：20.5012( )0.64490.7079aHsss用脈沖響應(yīng)不變法將用脈沖響應(yīng)不變法將Ha(s)轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)轉(zhuǎn)換成數(shù)字濾波器的系統(tǒng)函數(shù)H(z)。解：首先將解：首先將Ha(s)寫成部分分式寫成部分分式0.32240.3224( )0.32240.77720.32240.7772ajjHssjsj極點(diǎn)為：極點(diǎn)為：12(0.32240.772),(0.32240.7772)sjsj 74那么那么H(z)的極點(diǎn)為：的極點(diǎn)為：1212,s Ts Tzeze設(shè)設(shè)T

45、=1s時(shí)用時(shí)用H1(z)表示，表示，T=0.1s時(shí)用時(shí)用H2(z)表示，則表示，則111212120.3276( )1 1.03280.2470.0485( )1 1.93070.9375zH zzzzHzzz轉(zhuǎn)換時(shí)，也可以直接按照二階節(jié)的公式來直接寫。轉(zhuǎn)換時(shí)，也可以直接按照二階節(jié)的公式來直接寫。75)()(21 jjeHeH和和)( jH76優(yōu)點(diǎn)：優(yōu)點(diǎn)： 脈沖響應(yīng)不變法使得數(shù)字濾波器的單位脈沖響應(yīng)完脈沖響應(yīng)不變法使得數(shù)字濾波器的單位脈沖響應(yīng)完 全模仿模擬濾波器的單位沖激響應(yīng)，也就是時(shí)域逼全模仿模擬濾波器的單位沖激響應(yīng)，也就是時(shí)域逼 近良好。近良好。 模擬頻率模擬頻率和數(shù)字頻率和數(shù)字頻率之間呈

46、線性關(guān)系之間呈線性關(guān)系=T。 因而，一個(gè)線性相位的模擬濾波器（例如貝塞爾濾因而，一個(gè)線性相位的模擬濾波器（例如貝塞爾濾 波器）通過脈沖響應(yīng)不變法得到的仍然是一個(gè)線性波器）通過脈沖響應(yīng)不變法得到的仍然是一個(gè)線性 相位的數(shù)字濾波器。相位的數(shù)字濾波器。77缺點(diǎn)：缺點(diǎn)： 有頻率響應(yīng)的混疊效應(yīng)。有頻率響應(yīng)的混疊效應(yīng)。 所以，所以， 脈沖響應(yīng)不變脈沖響應(yīng)不變 法只適用于限帶的模擬濾波器法只適用于限帶的模擬濾波器(例如，衰減特性例如，衰減特性 很好的低通或帶通濾波器很好的低通或帶通濾波器)，而且高頻衰減越快，而且高頻衰減越快 ，混疊效應(yīng)越小。，混疊效應(yīng)越小。 對于高通和帶阻濾波器，由于它們在高頻部分對于高通

47、和帶阻濾波器，由于它們在高頻部分 不衰減，不衰減， 因此會產(chǎn)生混疊現(xiàn)象。因此會產(chǎn)生混疊現(xiàn)象。78第四節(jié)第四節(jié) 用雙線性變換法用雙線性變換法 設(shè)計(jì)設(shè)計(jì)IIR數(shù)字低通濾波器數(shù)字低通濾波器一、變換原理一、變換原理 了解決上一節(jié)的混疊問題，我們應(yīng)想辦法讓了解決上一節(jié)的混疊問題，我們應(yīng)想辦法讓S到到Z域的映射域的映射為為單值映射單值映射。方法：方法： 先采用線性壓縮的方法，將模擬角頻率先采用線性壓縮的方法，將模擬角頻率 壓縮至折疊頻率壓縮至折疊頻率以內(nèi)。即：以內(nèi)。即： 若：若： ： - 0 , 壓縮后得到壓縮后得到 1 1：-/T 0 /T )21tan(1T 非線性壓縮非線性壓縮SS1單值映射單值映射

48、S1Z)21tan(1T Tsez1 80找找Z與與S的關(guān)系：的關(guān)系：，得：由：)21tan(1T 2)(2)()21cos()21sin(2222111111TjTjTjTjeejeeTT )()(22221111TjTjTjTjeeeej )1()1()()(1111112222TsTsTsTsTsTseeeeees )1()1(11 zzs)1()1(ssz 81一般來說，為使一般來說，為使AF與與DF的某一頻率有對應(yīng)關(guān)系，可引入常數(shù)的某一頻率有對應(yīng)關(guān)系，可引入常數(shù)C：)21tan(1TC )1()1(11 zzCs)()(sCsCz 及及： 變換常數(shù)變換常數(shù)C的選?。旱倪x?。?常數(shù)常數(shù)

49、C用來調(diào)節(jié)頻率間的對應(yīng)關(guān)系：用來調(diào)節(jié)頻率間的對應(yīng)關(guān)系： 若希望若希望AF與與DF在低頻處有較為確切的對應(yīng)關(guān)系，在低頻處有較為確切的對應(yīng)關(guān)系， 可以選擇：可以選擇：TC2 82 解釋：解釋：在低頻處有較為確切的對應(yīng)關(guān)系，即要求低頻處：在低頻處有較為確切的對應(yīng)關(guān)系，即要求低頻處：1 當(dāng)當(dāng) 1較小時(shí)，處于低頻處，此時(shí)有：較小時(shí)，處于低頻處，此時(shí)有：TTtg1121)21( 而而 和和 1的對應(yīng)公式為：的對應(yīng)公式為：TCTtgC1121)21( TC2 sTsTzzzTs 22)1()1(21183二、優(yōu)點(diǎn)、問題及其解決辦法二、優(yōu)點(diǎn)、問題及其解決辦法 一、優(yōu)點(diǎn)：一、優(yōu)點(diǎn)：)21tan(2)21tan(

50、11TTTC 二、問題：二、問題：為為非非線線性性關(guān)關(guān)系系！的的對對應(yīng)應(yīng)關(guān)關(guān)系系與與 )21tan(2)21tan(21 TTT84 三、解決：三、解決：0 T/0 00)2arctan(2)2arctan(2 TSZ0 指標(biāo)轉(zhuǎn)換指標(biāo)轉(zhuǎn)換預(yù)畸預(yù)畸指標(biāo)轉(zhuǎn)換指標(biāo)轉(zhuǎn)換不預(yù)畸不預(yù)畸2tan20 TSZ002tan22arctan2)2arctan(2 TTT產(chǎn)生畸變產(chǎn)生畸變不產(chǎn)生畸變不產(chǎn)生畸變85 利用模擬濾波器設(shè)計(jì)利用模擬濾波器設(shè)計(jì)IIR數(shù)字低通濾波器的步驟。數(shù)字低通濾波器的步驟。確定數(shù)字低通濾波器的技術(shù)指標(biāo)：通帶截止頻率確定數(shù)字低通濾波器的技術(shù)指標(biāo)：通帶截止頻率p、通通 帶衰減帶衰減 p、阻帶截止

51、頻率阻帶截止頻率s、阻帶衰減阻帶衰減 s將將數(shù)字?jǐn)?shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換換成成模擬模擬低通濾波器的低通濾波器的 技術(shù)指標(biāo)。技術(shù)指標(biāo)。如果采用如果采用脈沖響應(yīng)不變法脈沖響應(yīng)不變法，邊界頻率的轉(zhuǎn)換關(guān)系為，邊界頻率的轉(zhuǎn)換關(guān)系為如果采用如果采用雙線性變換法雙線性變換法，邊界頻率的轉(zhuǎn)換關(guān)系為，邊界頻率的轉(zhuǎn)換關(guān)系為T )21tan(2 T(3) 按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器。按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器。(4) 將模擬濾波器將模擬濾波器Ha(s)，從從s平面轉(zhuǎn)換到平面轉(zhuǎn)換到z平面平面，得到數(shù)字低通，得到數(shù)字低通 濾波器系統(tǒng)函數(shù)濾波器系統(tǒng)函數(shù)H(

52、z)。86例：設(shè)計(jì)低通數(shù)字濾波器，要求在通帶內(nèi)頻率低于例：設(shè)計(jì)低通數(shù)字濾波器，要求在通帶內(nèi)頻率低于0.2rad時(shí)，時(shí)， 容許幅度誤差在容許幅度誤差在1dB以內(nèi)；在頻率以內(nèi)；在頻率0.3到到之間的阻帶衰之間的阻帶衰 減大于減大于15dB。指定模擬濾波器采用巴特沃斯低通濾波器。指定模擬濾波器采用巴特沃斯低通濾波器。 試分別用脈沖響應(yīng)不變法和雙線性變換法設(shè)計(jì)濾波器。試分別用脈沖響應(yīng)不變法和雙線性變換法設(shè)計(jì)濾波器。解：解： (1) 用脈沖響應(yīng)不變法設(shè)計(jì)數(shù)字低通濾波器。用脈沖響應(yīng)不變法設(shè)計(jì)數(shù)字低通濾波器。 數(shù)字低通的技術(shù)指標(biāo)為數(shù)字低通的技術(shù)指標(biāo)為 p=0.2rad, p=1dB; s=0.3rad, s

53、=15dB 模擬低通的技術(shù)指標(biāo)為模擬低通的技術(shù)指標(biāo)為 T=1s,p=0.2rad/s, p=1dB; s=0.3rad/s, s=15dB87設(shè)計(jì)巴特沃斯低通濾波器。先計(jì)算階數(shù)設(shè)計(jì)巴特沃斯低通濾波器。先計(jì)算階數(shù)N及及3dB截止頻率截止頻率c。0.10.1lglg0.31.50.21010.092101lg0.0925.884lg1.5psspspssppspkNkN 88根據(jù)階數(shù)根據(jù)階數(shù)N=6，查表得到歸一化傳輸函數(shù)為：查表得到歸一化傳輸函數(shù)為：234561( )13.86377.46419.14167.46413.8637aHppppppp由式：由式：Npcp2110)110( 得到：得到：

54、c=0.7032rad/s 最后，去歸一化，將最后，去歸一化，將p=s/c代入代入Ha(p)中，得到實(shí)際的傳中，得到實(shí)際的傳輸函數(shù)輸函數(shù)Ha(s)8962652433425665432( )3.86377.46419.14167.46413.86370.12092.7163.6913.1791.8250.1210.1209accccccHsssssssssssss用脈沖響應(yīng)不變法將用脈沖響應(yīng)不變法將Ha(s)轉(zhuǎn)換成轉(zhuǎn)換成H(z)。首先將首先將Ha(s)進(jìn)行部分進(jìn)行部分 分式分解，然后再轉(zhuǎn)換分式分解，然后再轉(zhuǎn)換 z=esT。1112121120.28710.44662.14281.1454( )

55、10.12970.69491 1.06910.36991.85580.630410.99720.2570zzH zzzzzzzz90(2) 用雙線性變換法設(shè)計(jì)數(shù)字低通濾波器。用雙線性變換法設(shè)計(jì)數(shù)字低通濾波器。 數(shù)字低通技術(shù)指標(biāo)仍為數(shù)字低通技術(shù)指標(biāo)仍為 p=0.2rad, p=1dB; s=0.3rad, s=15dB 模擬低通的技術(shù)指標(biāo)為模擬低通的技術(shù)指標(biāo)為21tan,122tan0.10.65/ ,12tan0.151.019/ ,15ppPpssTTrad sdBrad sdB 91 設(shè)計(jì)巴特沃斯低通濾波器。階數(shù)設(shè)計(jì)巴特沃斯低通濾波器。階數(shù)N計(jì)算如下：計(jì)算如下：lglg1.0191.568

56、0.650.092lg0.0925.306lg1.568spspssppspkNkN 取取N=6。并求得：并求得：c=0.7662rad/s92 根據(jù)根據(jù)N=6，查表查表6.2.1得到的歸一化傳輸函數(shù)得到的歸一化傳輸函數(shù)Ha(p)與脈沖與脈沖響應(yīng)不變法得到的相同。為去歸一化，將響應(yīng)不變法得到的相同。為去歸一化，將p=s/c代入代入Ha(p)，得實(shí)際的得實(shí)際的Ha(s)2220.2024( )(0.3960.5871)(1.0830.5871)(1.4800.5871)aHsssssss 用雙線性變換法將用雙線性變換法將Ha(s)轉(zhuǎn)換成數(shù)字濾波器轉(zhuǎn)換成數(shù)字濾波器H(z)111 612121211

57、20.0007378(1)( )( )(1 1.2680.7051)(1 1.0100.358)110.90440.2155azszzH zHszzzzzz93第五節(jié)第五節(jié) 數(shù)字高通、帶通和帶阻濾波器的設(shè)計(jì)數(shù)字高通、帶通和帶阻濾波器的設(shè)計(jì)(1) 確定所需類型數(shù)字濾波器的技術(shù)指標(biāo)確定所需類型數(shù)字濾波器的技術(shù)指標(biāo)(2) 將所需類型數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成所需類型模擬將所需類型數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成所需類型模擬 濾波器的技術(shù)指標(biāo)，轉(zhuǎn)換公式為：濾波器的技術(shù)指標(biāo)，轉(zhuǎn)換公式為： 若用脈沖響應(yīng)不變法：若用脈沖響應(yīng)不變法： 若用雙線性變換法：若用雙線性變換法：T )21tan(2 T (4) 設(shè)計(jì)模擬低

58、通濾波器。設(shè)計(jì)模擬低通濾波器。(5) 將模擬低通通過頻率變換，轉(zhuǎn)換成所需類型的模擬濾波器。將模擬低通通過頻率變換，轉(zhuǎn)換成所需類型的模擬濾波器。(6) 將所需類型的模擬濾波器轉(zhuǎn)換成所需類型的數(shù)字濾波器。將所需類型的模擬濾波器轉(zhuǎn)換成所需類型的數(shù)字濾波器。(3) 將所需類型模擬濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技將所需類型模擬濾波器技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器技 術(shù)指標(biāo)。術(shù)指標(biāo)。94例：設(shè)計(jì)一個(gè)例：設(shè)計(jì)一個(gè)數(shù)字高通數(shù)字高通濾波器，要求通帶截止頻率濾波器，要求通帶截止頻率p=0.8rad，通帶衰減不大于通帶衰減不大于3dB，阻帶截止頻率阻帶截止頻率s=0.44rad,阻帶衰減阻帶衰減不小于不小于15dB。希望采用巴特沃斯型濾波器。希望采用巴特沃斯型濾波器。解：解： (1)數(shù)字高通的技術(shù)指標(biāo)為數(shù)字高通的技術(shù)指標(biāo)為p=0.8rad, p=3dB; s=0.44rad, s=15dB(2) 模擬高通的技術(shù)指標(biāo)計(jì)算如下：模擬高通的技術(shù)指標(biāo)計(jì)算如下： 令令T=1，則有則有12tan6.155/ ,3212