【正文】 ??zzzzzzsHZHzzTsa? ?? ? ? ?? ? ? ? ? ? ? ?? ?? ? ? ?? ?? ? ? ?? ?? ? ? ?? ? ? ?? ?? ?? ? ? ?? ?? ?? ?? ?? ?231112312111121131311131313111113131311212113131312111312122122211111141111111211111211211??????????????????????????????????????????????????????????????????????????????????????zzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzz 圖 1 三階 Butterworth 數(shù)字濾波器的頻響 脈沖響應不變法 雙線性變換法 fs/2 我們也可以用 MA T L A B 完成設計，在 MA T L A B 中相關的語句有b u t t e r （ 巴特沃茲濾波器）、 i m p i n v a r （ 脈沖響應不變法 ）、b i l i n e a r （ 雙線性變換），具體的程序如下：[B , A ] = b u t t e r ( 3 , 2 * p i * 1 0 0 0 , 39。為此，實際應用脈沖響應不變法時稍作一點修改，即求出 H（ Z）后，再乘以因子 T,使 H（ Z）只與 有關，即只與 fc和 fs的相對值 有關，而與采樣頻率 fs無直接關系。 解： a. 脈沖響應不變法 由于脈沖響不變法的頻率關系是線性的 ， 所以可直接按 Ωc =2πfc設計 Ha(s)。 2） 由變換關系將 {ωk}映射到模擬域 ， 得出模擬濾波器的臨界頻率值 {Ωk}。 ３ )計算 H(Z) 雙線性變換比脈沖響應法的設計計算更直接和簡單 。1|| ??? zz 時?時 小結 1) 與脈沖響應不變法相比，雙線性變換的主要優(yōu)點： S平 面與 Z平面是單值的一一對應關系（靠頻率的嚴重 非線性關系得到的），即整個 jΩ軸單值的對應于單位 圓一周，關系式為： ????????22 ?tgT可見， ω和 Ω為非線性關系，如圖所示。 實現(xiàn)方法： 通過標準變換關系映射到 Z平面 ， 將 Ω 1映射到 Z平面的單位圓上 Tsez 1?通常取 C=2/T, TsTseecTscs1111)2(th 1 ??????????????????zzTsTsez 1?再將 S1 平面通過標準變換關系映射到 Z平面，即令 將這一關系解析擴展至整個 S平面，則得到 S平面到 S1 最后得 S平面與 Z 雙線性換法的主要優(yōu)點是 S平面與 Z平面一一單值對應 ，S平面的虛軸 (整個 jΩ)對應于 Z平面 單位圓的一周 ， S平面的Ω=0處對應于 Z平面的 ω=0處 ， 對應即數(shù)字濾波器的頻率響應終 止于折疊頻率處 ， 所以雙線性變換不存在混迭效應 。 例如線性相位的貝塞爾低通濾波器，通過脈沖響應不變法 2)在某些場合，要求數(shù)字濾波器在時域上能模仿模擬濾波器的功能時，如要實現(xiàn)時域沖激響應的模仿，一般使用脈沖響應 TjHeHaj /)()( ?? ???? 3)如果 Ha(s)是穩(wěn)定的 ， 即其極點在 S左半平面 ， 映射后得到的 H(Z)也是穩(wěn)定的 。我們只要知道濾波器的階數(shù)，就可直接查出低通原型的系統(tǒng)函數(shù)。模擬濾波器頻響 在折疊頻率以上衰減越大，失真則越小， 這時，采用脈沖響應不變法設計的數(shù)字濾 波器才能得到良好的效果。 STez? 如果模擬濾波器的系統(tǒng)函數(shù)是穩(wěn)定的，其極點應位于左半平面 11 1)( ?? ?? ?zeAzHTskNkk0][ ?ke sR 對應 Z平面的極點有 1?? Tsk kez 位于單位園內。 1．脈沖響應不變法 原理： 從時域響應出發(fā)，使求得的數(shù)字濾波器的單位脈沖響應 h(n)等于模擬濾波器的單位沖激響應 h(t)的抽樣值。由指定的技術指標 利用上述公式和表格進行設計時，最關鍵的 2個參數(shù)是濾波器的節(jié)數(shù) N和 3dB 截止頻率 。 模擬濾波器設計 IIR濾波器的設計是基于模擬濾波器的成熟技術而完成的 簡單介紹模擬濾波器設計的一些基本概念，并介紹兩種常用的濾波器的設計方法 ： ? 巴特沃思 (Butterworth)濾波器 ? 切比雪夫 (Chebyshev)濾波器 6． 2． 1模擬濾波器設計的基本概念 1 . 模擬濾波器的頻率特性與衰減特性 濾波器的頻率特性主要取決于構成濾波器系統(tǒng)的系統(tǒng)函數(shù) ???? jssHjH |)()( 工程設計中給定的指標往往是通帶和阻帶的衰減 ， 它一般用反映功率增益的幅度平方函數(shù)或稱模方函數(shù)來定義 dBjHjHA |)(|lg20|)(|lg10)( 2 ??????? 當要求濾波器具有線性相位特性 （ 延時 τ 為常數(shù) ） 時濾波器的頻率特性為 H j H j e j( ) | ( )| ( )? ? ?? ? ???? ?? )(，2. 歸一化與頻率變換 采用歸一化參數(shù) ? 設計結果具有普遍性 ? 計算方便 歸一化包含： ? 電路參數(shù)歸一化：將系統(tǒng)中無源元件的阻抗或運算阻抗分別除以基準電阻 (系統(tǒng)的負載電阻值 )。 3．設計 IIR濾波器的幾種方法 IIR數(shù)字濾波器的系統(tǒng)函數(shù)可表示為的有理分式 11101)(????????zazbzHkNkkNk? 設計 IIR濾波器的系統(tǒng)函數(shù)，就是要確定 H(z)的階數(shù) N（通常稱 N為濾波器的階數(shù)）以及分子分母多項式的系數(shù) ?? jezjkk zHeHba ?? )()(，使其、 滿足指定的頻率特性 (1)利用模擬濾波器的理論來設計 模擬濾波器研究較早，理論已經十分成熟，有許多簡單而嚴謹?shù)脑O計公式和大量的圖表可以利用，利用這些現(xiàn)有技術來解決數(shù)字濾波器的設計問題 采用這種方法時，要先要設計一個合適的模擬濾波器，然后將它轉換成滿足給定指標的數(shù)字濾波器 這種方法適合于設計幅頻特性比較規(guī)則的濾波器，例如低通、高通、帶通、帶阻等 當把模擬濾波器的 H(s)轉換成數(shù)字濾波器的 H(z) 時，要實現(xiàn) S平面向 Z平面的映射，必須滿足兩個條件 必須保證模擬頻率映射為數(shù)字頻率，且保證兩者的頻率特性基本一致 ? 要求變換后代表 S平面的虛軸 jΩ應映射到 Z平面的單位圓 ? 且數(shù)字濾波器的頻率響應和模擬濾波器頻率響應的形狀應基本保持不變； 因果穩(wěn)定的模擬濾波器系統(tǒng)函數(shù) H(s)轉換成數(shù)字濾波器傳輸函數(shù) H(z)后，仍然是因果穩(wěn)定的 ? 要求 S平面左半平面的極點必須映射到 Z平面的單位圓內 兩種常用的方法 ? 脈沖響應不變法：從時域的角度出發(fā)進行映射 ? 雙線性不變法：從頻域角度出發(fā)進行映射 (2)利用最優(yōu)化技術進行 CAD設計 若需設計濾波器的幅頻特性是任意的或者形狀比較復雜，可采用計算機輔助設計 (CAD)方法進行優(yōu)化設計 設計思想 )( ?jd eH 希望濾波器的幅頻響應 ： )( ?jeH 設計 濾波器的幅頻響應 ： ? 選擇一種最優(yōu)化的準則，例如采用最小均方誤差準則 )()( ?? jjd eHeH 、設在指定的一組離散的頻率點 Mii ,2,1},{ ???的均方誤差 221[ ( ) ( ) ]iiM jjdiH e H e???????? 求解 H(z)的系數(shù) ，、 kk ba 使 均方誤差 最小 ? 當 濾波器階數(shù) N 較高時， 轉換為 一個多變量最優(yōu)化問題，需要大量的迭代運算，因此必須采用 CAD的方法。 若給定的指標 =3dB, 即通帶邊頻 時 , ε =1,可求得 cp ???pA(b)從模方函數(shù)求系統(tǒng)函數(shù) H(s) ①求得極點 ,2,1)/(1 1|)(| 22 ??????? NjjjH Nc帶入上式，得：?? js0)(1 2 ??? Ncjs?)12()1( ??? kje由于Nkes NNkjck 2,2,12/)12( ???? ?? ?分析討論 ?c在歸一化頻率的情況 =1， 極點均勻分布在單位圓上 Nkes NNkjk 2,2,12/)12( ??? ?? ?對于物理可實現(xiàn)系統(tǒng) ,它的所有極點均應在 s的左半平面上 ②系統(tǒng)函數(shù)的構成 濾波器的極點求出后，可取左平面上的所有極點構成系統(tǒng)函數(shù) )(1)(1iNissAsH????對于低通濾波器 ， 為了保證在頻率零點 =0處 ， =1， 可取 ?|)(| ?jH????NiiN sA1)1()()1()(1 iiNiNssssH??? ??因此得 例 621舉例說明系統(tǒng)函數(shù)的構成 設計一巴特沃思濾波器 ， 使其滿足以下指標 : 通帶邊頻 =100k rad/s, 通帶的最大衰減為 =3dB， 阻帶邊頻為 =400k rad/s,阻帶的最小衰減為 =35 dB p?Ap s?As解 :由于通帶邊頻就是 3dB 截止頻率，即 cp ???1110 ??? pA? ??? sA?確定階次 N ， )/l g ( )/l g ( ?????csN ?? 3?N取求左半平面的極點 : s e kk c j k N N? ?? ?? ( ) / , ,2 1 2 1 2 3?,3/21 ?jc es ?? ,2 ?jc es ?? s ec j3 2 3? ?? ? /得極點： 構成巴特沃思濾波器傳輸函數(shù) H(s)為 H s s s ss s s s s s s s scc c c( ) ( )( )( )? ?? ? ? ? ? ? ?1 2 31 2 333 2 2 32 2?? ? ?相對截止頻率 歸一化，得歸一化巴特沃思濾波器傳輸函數(shù) ?c H sa ( )H s s s sa ( ) ? ? ? ?12 2 13 2③一般 N階歸一化巴特沃思濾波器傳輸函數(shù) 表示 H ss s a s a s a s saiNiNN N( )( )???? ? ? ? ?????1 1111 2211?是 =1時的極點，分布在單位圓上 si ?c分母一般稱為巴特沃思多項式，其系數(shù)可通過查表求得 ,見表 521 表 621 巴特沃思多項式系數(shù) N a1 a2 a3 a4 a5 a6 a7 a8 a9 2 3 4 5 6 7 8 9 10 表 622 巴特沃思多項式因式分解 N 巴 特 沃 思 多 項 式 1 s+1 2 s2++1 3 (s+1)(s2+s+1) 4 (s2++1)(s2++1) 5 (s+1)(s2++1)(s2++1) 6 (s2++1)(s2 ++1)(s2 ++1) 7 (s+1)(s2++1)(s2++1)(s2++1) 8 (s2++1)(s2++1)(s2++1)(s2++1) 9 (s+1)(s2++1)(s2+s+1)(s2++1)(s2++1) 上述歸一化公式和表格是相對 3dB 截止頻率 給出的。 2） Ha(s) 的因果穩(wěn)定性映射成 H（ z）后保持不變，即 S平面的左半平面 Re{S}＜ 0 應映射到 Z平面的單位圓以內 |Z|1。 ?? ?? ???????????????????????????nn s TnststnenThdtenTtthdtenTtthsH)()()(])()