第4章無限長單位脈沖響應iir濾波器的設計方法

1、第第4 4章章 無限長單位脈沖無限長單位脈沖響應（響應（IIRIIR）濾波器設計）濾波器設計 4.1 濾波器的基本原理濾波器的基本原理4.2 模擬濾波器設計方法模擬濾波器設計方法4.3 根據模擬濾波器設計根據模擬濾波器設計IIR濾波器濾波器4.4 從模擬濾波器低通原型到各種數字從模擬濾波器低通原型到各種數字濾波器的頻率變換濾波器的頻率變換4.5 從低通數字濾波器到各種數字濾波從低通數字濾波器到各種數字濾波器的頻率變換器的頻率變換4.6 IIR數字濾波器的最優(yōu)化設計方法數字濾波器的最優(yōu)化設計方法概述概述: 許多信息處理過程，如信號的過濾，檢測、預許多信息處理過程，如信號的過濾，檢測、預測等都要用

2、到濾波器，數字濾波器是數字信號處理測等都要用到濾波器，數字濾波器是數字信號處理中使用得最廣泛的一種線性系統(tǒng)，是數字信號處理中使用得最廣泛的一種線性系統(tǒng)，是數字信號處理的重要基礎。的重要基礎。 數字濾波器的功能（本質）是將一組輸入的數數字濾波器的功能（本質）是將一組輸入的數字序列通過一定的運算后轉變?yōu)榱硪唤M輸出的數字字序列通過一定的運算后轉變?yōu)榱硪唤M輸出的數字序列。實現方法主要有兩種：數字信號處理硬件和序列。實現方法主要有兩種：數字信號處理硬件和計算機軟件。計算機軟件。 數字濾波器數字濾波器線性時不變系統(tǒng)。線性時不變系統(tǒng)。iaibicid數字濾波器的設計步驟：數字濾波器的設計步驟：1）按照任務要

3、求確定濾波器的性能要求。）按照任務要求確定濾波器的性能要求。2）用一個因果穩(wěn)定的離散線性時不變系統(tǒng)的）用一個因果穩(wěn)定的離散線性時不變系統(tǒng)的 H(z) 或或 h(n) 去逼近這個性能要求，即求去逼近這個性能要求，即求H(z)或或 h(n) 的表的表達式。確定系數達式。確定系數 、 或零極點或零極點 、 ，以，以使濾波器滿足給定的性能要求使濾波器滿足給定的性能要求第四章、第五章第四章、第五章討論討論 3）數字濾波器的實現。包括選擇運算結構，確定）數字濾波器的實現。包括選擇運算結構，確定運算和系數存儲的字長，選用通用計算機及相應的運算和系數存儲的字長，選用通用計算機及相應的軟件或專用數字濾波器硬件實

4、現這一系統(tǒng)。軟件或專用數字濾波器硬件實現這一系統(tǒng)。 選擇合適的字長和有效數字的處理方法等選擇合適的字長和有效數字的處理方法等(第六第六章章)。 4.1 濾波器的基本原理濾波器的基本原理濾波：濾除信號中不需要的分量，保留有用的分量。濾波：濾除信號中不需要的分量，保留有用的分量。離散系統(tǒng)的數字模型：離散系統(tǒng)的數字模型： MiiNiiinxainybny01離散系統(tǒng)的系統(tǒng)函數：離散系統(tǒng)的系統(tǒng)函數： NiiiMiiizbzazH101離散系統(tǒng)的單位脈沖響應：離散系統(tǒng)的單位脈沖響應： zHZnh1 nxnhnhnxny*4.1.1 濾波器的分類濾波器的分類 模擬模擬 遞歸系統(tǒng)遞歸系統(tǒng) IIR 數字數字

5、非遞歸系統(tǒng)非遞歸系統(tǒng) FIR從功能上：從功能上： 低通低通 LP low pass 高通高通 HP high pass 帶通帶通 BP band pass 帶阻帶阻 BS band stopjHccOjHccOjH22O11jH22O11模擬低通濾波器模擬低通濾波器AF LP模擬高通濾波器模擬高通濾波器AF HP模擬帶通濾波器模擬帶通濾波器AF BP模擬帶阻濾波器模擬帶阻濾波器AF BSjeHccOccO數字低通濾波器數字低通濾波器DF LP數字高通濾波器數字高通濾波器DF HP數字帶通濾波器數字帶通濾波器DF BP數字帶阻濾波器數字帶阻濾波器DF BSjeH22O11jeH22O11jeH4

6、.1.2 可實現濾波器的特性可實現濾波器的特性11-12c0r通帶過渡帶阻帶時c11min1lg2011lg201lg20jeHjjjeHeHeHAlg201lg201lg102時r22maxlg201lg201lg20jeHAt通帶波動通帶波動最小阻帶衰耗最小阻帶衰耗設計方法：設計方法： 1）先設計一個合適的模擬濾波器，然后變換）先設計一個合適的模擬濾波器，然后變換成滿足預定指標的數字濾波器。成滿足預定指標的數字濾波器。 由于模擬的網絡綜合理論已經發(fā)展得很成熟由于模擬的網絡綜合理論已經發(fā)展得很成熟模擬濾波器有簡單而嚴格的設計公式，設計起來模擬濾波器有簡單而嚴格的設計公式，設計起來方便、準確、

7、可將這些理論推廣到數字域，作為方便、準確、可將這些理論推廣到數字域，作為設計數字濾波器的工具。設計數字濾波器的工具。2）最優(yōu)化設計方法）最優(yōu)化設計方法分兩步：分兩步： a) 確定一種最優(yōu)準則，如最小均方誤差準則確定一種最優(yōu)準則，如最小均方誤差準則，即使設計出的實際頻率響應的幅度特性，即使設計出的實際頻率響應的幅度特性 （與所要求的理想頻率響應（與所要求的理想頻率響應 的均方誤的均方誤差最小，差最小，| )(|jeH| )(|jdeHmin)()(21MijdjiieHeHiaibiaib此外還有其他多種誤差最小準則，此外還有其他多種誤差最小準則， b) 在此最佳準則下，求濾波的系數在此最佳準則

8、下，求濾波的系數 和和 通過不斷地迭代運算，改變通過不斷地迭代運算，改變 、 ，直到，直到 滿足要求為止。滿足要求為止。 以上兩種設計方法中，著重講第一種，以上兩種設計方法中，著重講第一種，因為數字濾波器在很多場合所要完成的任務因為數字濾波器在很多場合所要完成的任務與模擬濾波器相同，如作低通、高通、帶通與模擬濾波器相同，如作低通、高通、帶通及帶阻網絡等，這時數字濾波也可看作是及帶阻網絡等，這時數字濾波也可看作是“模仿模仿”模擬濾波器。在模擬濾波器。在IIR濾波器設計中，濾波器設計中，采用這種設計方法目前最普遍。由于計算機采用這種設計方法目前最普遍。由于計算機技術的發(fā)展，最優(yōu)化設計方法的使用也逐

9、漸技術的發(fā)展，最優(yōu)化設計方法的使用也逐漸增多。增多。 0)()(dtethjHtjaa0sincos)()(dttjtthjHaa)()(jHjHaa 模擬濾波器的設計就是根據一組設計規(guī)范設計模擬模擬濾波器的設計就是根據一組設計規(guī)范設計模擬系統(tǒng)函數系統(tǒng)函數Ha(s)，使其逼近某個理想濾波器特性。，使其逼近某個理想濾波器特性。 因果系統(tǒng)中因果系統(tǒng)中式中式中ha(t)為系統(tǒng)的沖激響應，是實函數。為系統(tǒng)的沖激響應，是實函數。 不難看出不難看出4.2 模擬濾波器設計方法模擬濾波器設計方法定義幅度平方函數定義幅度平方函數 式中式中 Ha(s)模擬濾波器模擬濾波器 系統(tǒng)函數系統(tǒng)函數 Ha(j)濾波器的頻率

10、響應濾波器的頻率響應 |Ha(j)|濾波器的幅頻響應濾波器的幅頻響應又又 S=j，2=-S2 A(2)=A(-S2)|S=j) 1 ()()()()()()()()()(222jsaaaaaaasHsHjHjHAjHjHjHA問題：由問題：由A(-S2)Ha(S) 對于給定的對于給定的A(-S2)，先在，先在S復平面上標出復平面上標出A(-S2)的極點和零點，由的極點和零點，由(1)式知，式知，A(-S2)的極的極點和零點總是點和零點總是“成對出現成對出現”，且對稱于，且對稱于S平面平面的實軸和虛軸，選用的實軸和虛軸，選用A(-S2)的對稱極、零點的的對稱極、零點的任一半作為任一半作為Ha(s

11、)的極、零點，則可得到的極、零點，則可得到Ha(s)。 為了保證為了保證Ha(s)的穩(wěn)定性，應選用的穩(wěn)定性，應選用A(-S2)在在S左半平面的極點作為左半平面的極點作為Ha(s)的極點，零點可選的極點，零點可選用任一半。用任一半。例例4.1 設已知設已知 ，求對應的，求對應的42212A sHa解：解：42221222ssAsAs sHsHaa sjIm sReO2221j21j21j21j 21212jsjsssHa1222sssNcajjjHA22211)()()(2AN為濾波器階數為濾波器階數其幅度平方函數：其幅度平方函數：特點：具有通帶內最大平坦的振幅特性，且隨特點：具有通帶內最大平坦

12、的振幅特性，且隨f ，幅頻特，幅頻特 性性 單調單調 。4.2.1 巴特沃思巴特沃思(Butterworth)濾波器濾波器 巴特沃思濾波器巴特沃思濾波器 振幅平方函數振幅平方函數圖圖1中，中，N增加，通帶和阻帶的近似性越好，過渡帶增加，通帶和阻帶的近似性越好，過渡帶越陡。越陡。/c1時，時， ( /c)2N 1, 增加，增加， A(2)快速減小?？焖贉p小。=c, ， ,幅度衰減幅度衰減 ，相當于，相當于3dB衰減點。衰減點。 21jH21)(2A振幅平方函數的極點：振幅平方函數的極點：NcaajSSHSH2)(11)()()() 1(21cNPjSNkeeeScNkjcjNkjP21)()(2

13、21222112 令分母為零，得令分母為零，得Butterworth濾波器濾波器 的振幅平方函數有的振幅平方函數有2N個極個極點，它們均勻對稱地分布在點，它們均勻對稱地分布在|S|=c的圓周上。的圓周上。三階三階A(-S2)的極點分布的極點分布 考慮到系統(tǒng)的穩(wěn)定性，知考慮到系統(tǒng)的穩(wěn)定性，知AF的系統(tǒng)函數是由的系統(tǒng)函數是由S平面平面左半部分的極點（左半部分的極點（SP1，SP2，SP3）組成的，它們分別）組成的，它們分別為：為：3232321,jcpcpjcpeSSeS)()()(3213pppcaSSSSSSsH1221)(231SSSsHa1)/(2)/(2)/(1)(23cccassssH

14、系統(tǒng)函數為：系統(tǒng)函數為：1c令令 ，得歸一化的三階，得歸一化的三階BW AF：如果要還原的話，則有如果要還原的話，則有1221)(231SSSsHa歸一化的三階歸一化的三階BW AF：121)(21SSsHa歸一化的二階歸一化的二階BW AF：11)(1SsHa歸一化的一階歸一化的一階BW AF：補充：確定除數補充：確定除數N及及c1. 確定確定N(已知通帶邊頻已知通帶邊頻1、通帶波動、通帶波動、阻帶邊頻阻帶邊頻r、阻帶最小衰耗、阻帶最小衰耗At)NcjH21211lg101lg10NcrrjHAt221lg101lg10NcjHA221lg101lg101 . 021101NcAtNcr1

15、. 021011101 . 021Nc1101 . 02AtNcr1101101 . 01 . 021AtNr110110lglg21 . 01 . 01AtrN11 . 01 . 0lg110110lg21rAtN11 . 01 . 0lg110110lg21rAtN2. 確定確定c1101 . 021Nc1101 . 02AtNcrNc211 . 01110NAtrc211 . 0110若已知衰減為若已知衰減為3dB的頻率點，此頻率點即為的頻率點，此頻率點即為c4.2.2 切比雪夫（切比雪夫（chebyshev）濾波器）濾波器 特點：誤差值在規(guī)定的頻段上等幅變化。特點：誤差值在規(guī)定的頻段上

16、等幅變化。 巴特沃茲濾波器在通帶內幅度特性是單調下降的，巴特沃茲濾波器在通帶內幅度特性是單調下降的，如果階次一定，則在靠近截止頻率如果階次一定，則在靠近截止頻率 處，幅度下降處，幅度下降很多，或者說，為了使通常內的衰減足夠小，需要的很多，或者說，為了使通常內的衰減足夠小，需要的階次（階次（N）很高，為了克服這一缺點，采用切比雪夫）很高，為了克服這一缺點，采用切比雪夫多項式逼近所希望的多項式逼近所希望的 。 切比雪夫濾波器的切比雪夫濾波器的 在通帶范圍內是等幅起伏的，在通帶范圍內是等幅起伏的，所以同樣的通帶衰減，其階數較巴特沃茲濾波器要小。所以同樣的通帶衰減，其階數較巴特沃茲濾波器要小?？筛鶕?/p>

17、要對通帶內允許的衰減量（波動范圍）提出可根據需要對通帶內允許的衰減量（波動范圍）提出要求，如要求波動范圍小于要求，如要求波動范圍小于1db。 c2)( jH2)( jH振幅平方函數為振幅平方函數為)(11)()(2222cNaVjHAc1)coshcosh(1)coscos()(11xxNxxNxVN)(, 11)(,1xVxxxVxNN時有效通帶截止頻率有效通帶截止頻率與通帶波紋有關的參量，與通帶波紋有關的參量， 大大 ，波紋大。，波紋大。 0 1 VN（x）N階切比雪夫多項式，定義為階切比雪夫多項式，定義為 如圖如圖1,通帶內通帶內 變化范圍變化范圍1 c，隨，隨/c ， 0 (迅速趨于零

18、迅速趨于零)當當 =0時，時， N為偶數，為偶數， ，min ， N為奇數，為奇數， ， max，)2(cos11)0arccos(cos11)(22220NNjHa22011)(jHa1)(02jHa1)2(cos2N0)2(cos2N2111c2a)( jH2a)( jH 切比雪夫濾波器的振幅平方特性切比雪夫濾波器的振幅平方特性 2min111lg20)(1lg20jHa)1lg(1021101 . 02)(dB221)(,AjHar時 給定通帶波紋值分貝數給定通帶波紋值分貝數 后，可求后，可求 。有關參數的確定有關參數的確定: a、通帶截止頻率、通帶截止頻率c ，預先給定，預先給定 b、

19、通帶波動為、通帶波動為 c、階數N由阻帶的邊界條件確定。（ 、A事先給定） 222111AVcrNr221)(,AjHar時crNarAarNxNarxVx/cosh/1cosh)coshcosh()(,12得時22maxlg101lg10AjHAtaAtA1 . 0210crAtararNcosh110110cosh1 . 01 . 0 4.2.3 橢圓（Elliptic）濾波器 特點：幅值響應在通帶和阻帶內都是等波紋的，對于給定的階數和給定的波紋要求，橢圓濾波器能獲得較其它濾波器更窄的過渡帶寬，就這點而言，橢圓濾波器是最優(yōu)的。 其振幅平方函數為RN（，L）雅可比橢圓函數L表示波紋性質的參量

20、),(11)()(2222LRjHANa N=5， 的特性曲線 可見，在歸一化通帶內（-11）， 在（0，1）間振蕩，而超過L后， 在 間振蕩。這一特點使濾波器同時在通帶和阻帶具有任意衰減量。 ),(25LR2L),(25LR),(25LR 下圖為典型的橢園濾波器振幅平方函數 橢圓濾波器的振幅平方函數 圖中和A的定義 同切比雪夫濾波器rr當c、r、和A確定后，階次N的確定方法為：)1()()1()(1/212121kKkKkKkKNAkkrc確定參量確定參數2/1222/1210)1 ()1 ()(tktdtkK式中 為第一類完全橢圓積分 上面討論了三種最常用的模擬低通濾波器的特性和設計方法，

21、設計時按照指標要求，合理選用。 一般，相同指標下，橢圓濾波器階次最低，切比雪夫次之，巴特沃茲最高，參數的靈敏度則恰恰相反。 以上討論了由A（2 ）Ha (s),下面討論由Ha(s)H(Z)的變換設計法。4.2.4 模擬高通、帶通及帶阻濾波模擬高通、帶通及帶阻濾波器的設計器的設計設計過程：設計過程：1）將模擬濾波器的參數指標通過頻率變換轉化成）將模擬濾波器的參數指標通過頻率變換轉化成原型模擬低通濾波器的參數指標。原型模擬低通濾波器的參數指標。2）設計滿足指標要求的原型模擬低通濾波器。）設計滿足指標要求的原型模擬低通濾波器。3）通過頻率變換將原型模擬低通濾波器的系統(tǒng)函）通過頻率變換將原型模擬低通濾

22、波器的系統(tǒng)函數變換為其他類型（高通、帶通和帶阻）的模擬濾數變換為其他類型（高通、帶通和帶阻）的模擬濾波器。波器。4.2.5 用用MATLAB設計模擬濾波器設計模擬濾波器1. 巴特沃思模擬濾波器設計巴特沃思模擬濾波器設計1）求巴特沃思濾波器的階數）求巴特沃思濾波器的階數N,Wn=buttord(Wp,Ws,Rp,Rs,s)Wp,Ws：通帶、阻帶邊界頻率：通帶、阻帶邊界頻率(rad/s)Rp,Rs：通帶波動、最小阻帶衰減：通帶波動、最小阻帶衰減(dB)N：階數：階數Wn：3dB點頻率點頻率s：表示模擬域：表示模擬域2）給定階數計算濾波器的零極點和增益因子）給定階數計算濾波器的零極點和增益因子z,p

23、,k=buttap(N)z：零點，實際是一個空向量：零點，實際是一個空向量p：歸一化：歸一化(3dB的邊界頻率為的邊界頻率為1)的的N個極點個極點N：階數：階數k：增益因子：增益因子 NapspspsksH213）給定階數和邊界頻率設計濾波器）給定階數和邊界頻率設計濾波器B,A=butter(N,Wn,type,s)Wn：3dB點邊界頻率；若點邊界頻率；若Wn=W1,W2，則設計則設計2N階帶通或帶阻階帶通或帶阻B：系統(tǒng)函數分子的系數：系統(tǒng)函數分子的系數設計低通和帶通，設計低通和帶通，type可缺省，高可缺省，高通通high，帶阻，帶阻stopA：系統(tǒng)函數分母的系數：系統(tǒng)函數分母的系數2. 切

24、比雪夫模擬濾波器設計切比雪夫模擬濾波器設計1）計算濾波器的階數）計算濾波器的階數N,Wn=cheb1ord(Wp,Ws,Rp,Rs,s)Wp,Ws：通帶、阻帶邊界頻率：通帶、阻帶邊界頻率(rad/s)Rp,Rs：通帶波動、最小阻帶衰減：通帶波動、最小阻帶衰減(dB)N：階數：階數Wn：通帶邊界頻率：通帶邊界頻率s：表示模擬域：表示模擬域2）計算濾波器的零極點和增益因子）計算濾波器的零極點和增益因子z,p,k=cheb1ap(N,Rp)3）給定階數和指標參數設計濾波器）給定階數和指標參數設計濾波器B,A=cheby1(N,Rp,Wn,type,s)3. 橢圓模擬濾波器設計橢圓模擬濾波器設計1）計

25、算濾波器的階數）計算濾波器的階數N,Wn=ellipord(Wp,Ws,Rp,Rs,s)2）計算濾波器的零極點和增益因子）計算濾波器的零極點和增益因子z,p,k=ellipap(N,Rp)3）給定階數和指標參數設計濾波器）給定階數和指標參數設計濾波器B,A=ellip(N,Rp,Rs,Wn,type,s)4. 頻率變換頻率變換1）低通到高通變換）低通到高通變換B,A=lp2hp(num,den,w0)num,den：低通濾波器系統(tǒng)函數分子分母的系數：低通濾波器系統(tǒng)函數分子分母的系數w0：高通濾波器的通帶邊界角頻率：高通濾波器的通帶邊界角頻率B,A：高通濾波器系統(tǒng)函數分子分母的系數：高通濾波器系

26、統(tǒng)函數分子分母的系數2）低通到帶通變換）低通到帶通變換B,A=lp2bp(num,den,w0,Bw)w0：帶通濾波器的中心角頻率：帶通濾波器的中心角頻率Bw：帶通濾波器的通帶帶寬：帶通濾波器的通帶帶寬(rad/s)3）低通到帶阻變換）低通到帶阻變換B,A=lp2bs(num,den,w0,Bw)w0：帶阻濾波器的中心角頻率：帶阻濾波器的中心角頻率Bw：帶阻濾波器的阻帶帶寬：帶阻濾波器的阻帶帶寬(rad/s)例例4.2 設計滿足下列指標的模擬低通濾波器，分設計滿足下列指標的模擬低通濾波器，分別采用巴特沃思濾波器、切比雪夫濾波器和橢圓別采用巴特沃思濾波器、切比雪夫濾波器和橢圓濾波器。濾波器。kH

27、zfc1kHzfr2dBRp1dBRs40clear;fc=1000;fr=2000;rp=1;rs=40;wc=2*pi*fc;wr=2*pi*fr;N,wn=buttord(wc,wr,rp,rs,s);B,A=butter(N,wn,s);omega=0:200:10000*pi;h=freqs(B,A,omega);gain=20*log10(abs(h);plot(omega/(2*pi),gain);xlabel(頻率頻率/Hz);ylabel(幅度幅度/dB)title(巴特沃思濾波器巴特沃思濾波器)grid巴特沃思濾波器：巴特沃思濾波器：N=8 B=0, 0, 0, 0, 0,

28、 0, 0, 0, 6.21871030A=1, 3.6222104, 6.5603108, 7.70931012, 6.40601016, 3.84981020, 1.63601024, 4.51081027, 6.21871030 3074830102187. 6106222. 3102187. 6sssHbutter0500100015002000250030003500400045005000-120-100-80-60-40-200頻 率 /Hz幅度/dB巴 特 沃 思 濾 波 器切比雪夫濾波器：切比雪夫濾波器：N=5 B=0, 0, 0, 0, 0, 1.20281018A=1,

29、5.8862103, 6.6672107, 2.41701011, 9.04791014, 1.20281018 1845181102028. 12 .5886102028. 1sssHcheby0500100015002000250030003500400045005000-90-80-70-60-50-40-30-20-100頻 率 /Hz幅度/dB切 比 雪 夫 濾 波 器橢圓濾波器：橢圓濾波器：N=4 B=0, 0, 0, 0, 0.50171014A=1, 5.9008103, 5.9759107, 1.99361011, 5.63011014 183414106301. 58 .5

30、900105017. 0sssHellip0500100015002000250030003500400045005000-90-80-70-60-50-40-30-20-100頻 率 /Hz幅度/dB橢 圓 濾 波 器4.3 根據模擬濾波器設計根據模擬濾波器設計IIR濾波器濾波器利用模擬濾波器設計數字濾波器，就是從已知的利用模擬濾波器設計數字濾波器，就是從已知的模擬濾波器系統(tǒng)函數模擬濾波器系統(tǒng)函數Ha(s)設計數字濾波器系統(tǒng)函數設計數字濾波器系統(tǒng)函數H(z)，這歸根到底是一個由，這歸根到底是一個由S平面到平面到Z平面的變換，平面的變換，這種映射變換應遵循兩個基本原則：這種映射變換應遵循兩個基

31、本原則：1）H(z)的頻響與的頻響與Ha(s)的頻響保持一致，即的頻響保持一致，即s平面平面的虛軸必須映射到的虛軸必須映射到z平面的單位圓平面的單位圓 上。上。2）因果穩(wěn)定的）因果穩(wěn)定的Ha(s) 應能映射成應能映射成 因果穩(wěn)定因果穩(wěn)定H(z)。也就是也就是s平面的左半平面平面的左半平面 Res0 應該映射到應該映射到z平面的平面的單位圓以內單位圓以內|z| HP 。b. 高通變換高通變換jjjeee1,00LP-HP變換把變換把在上述在上述LP-LP 變換中，將變換中，將 Z代以代以Z , 得得 LP - HP變換關系：變換關系：11111111zzzzzGu 原型低通的截止頻率原型低通的截

32、止頻率 對應于高通的邊界頻率對應于高通的邊界頻率 ，欲將欲將 變換到變換到 ：cccccccjjjeee12cos2coscccc:確定ccccjjjeee1ccLP Hp變換00 LP - HP變換變換例例4.13 低通數字濾波器的通帶邊界頻率低通數字濾波器的通帶邊界頻率 的切比雪夫濾波器的切比雪夫濾波器求邊界頻率求邊界頻率 的數字高通濾波器。的數字高通濾波器。2 . 0c6 . 0c)8482. 04996. 11)(6493. 05548. 11 ()1 (001836. 0)(2121411uuuuuuHp38197. 02/ )2 . 06 . 0cos(2/ )2 . 06 . 0

33、cos(11138197. 0138197. 0zzu)7647. 0556. 01)(4019. 00416. 11 ()1 (02426. 0)()(21214138197. 0138197. 01111zzzzzuHzHzzup4.5.3 數字低通數字低通-數字帶通數字帶通0000000 2,0N即全通函數的階數為時， 所以全通函數取負號。，時, 1) 1 (,0G帶通的中心頻率應對應于低通原型的通帶中心，帶通的中心頻率應對應于低通原型的通帶中心，即即11zGu1*1*111zzzz由以上分析得變換關系：由以上分析得變換關系：1)(1122211211zrzrrzrzzGu1122212

34、jjjjjererreree或：或：LP-BP變換變換102cc004.5.4 數字低通數字低通-數字帶阻數字帶阻0000002N 0故，又又 G(1)=1, 所以，全通函數取正號。所以，全通函數取正號。帶阻的中心頻率帶阻的中心頻率由以上分析得變換關系：由以上分析得變換關系： 或或 1)(1122211211zrzrrzrzzGu1122212jjjjjererrereeLP-BS變換變換201cc)(jeH004.6 IIR數字濾波器的最優(yōu)化設計方法數字濾波器的最優(yōu)化設計方法4.6.1 帕德（帕德（Pade）逼近法）逼近法Hd(n)是一個理想濾波器的單位脈沖響應。是一個理想濾波器的單位脈沖響應。用一個因果濾波器逼近該理想濾波器，其單位脈沖用一個因果濾波器逼近該理想濾波器，其單位脈沖響應為響應為h(n)，則系統(tǒng)函數為：，則系統(tǒng)函數為：0101)()(nNiiiMiiinzbzaznhzH求得求得H(z)的的N+M+1個待定的系數，使得個待定的系數，