【正文】 節(jié)系統(tǒng)通過改變輸入系數(shù)即可完成； ? 極點(diǎn)位置可單獨(dú)調(diào)整； ? 運(yùn)算速度快 (可并行進(jìn)行 )； ? 各二階網(wǎng)絡(luò)的 誤差互不影響，總的誤差小，對字長要求低。 缺點(diǎn)： 不能直接調(diào)整零點(diǎn)，因多個(gè)二階節(jié)的零點(diǎn)并不是整個(gè)系統(tǒng)函數(shù)的零點(diǎn)，當(dāng)需要準(zhǔn)確的傳輸零點(diǎn)時(shí)，級聯(lián)型最合適。 (3)IIR 濾波器的并聯(lián)型結(jié)構(gòu)。 優(yōu)點(diǎn)： ? 簡化實(shí)現(xiàn)，用一個(gè)二階節(jié)，通過變換系數(shù)就可實(shí)現(xiàn)整個(gè)系統(tǒng)； ? 極、零點(diǎn)可單獨(dú)控制、調(diào)整，調(diào)整 α1i、 α2i只單獨(dú)調(diào)整了第 i 對零點(diǎn)，調(diào)整 β1i、 β2i則單獨(dú)調(diào)整了第 i 對極點(diǎn)； ? 各二階節(jié)零、極點(diǎn)的搭配可互換位置，優(yōu)化組合以減小運(yùn)算誤差； ? 可流水線操作。 缺點(diǎn)： ? 二階階電平難控制，電平大易導(dǎo)致溢出，電平小則使信噪比減小 。 課程設(shè)計(jì) 第 8 頁 共 44 頁 a、直接型 b、并聯(lián)型 課程設(shè)計(jì) 第 9 頁 共 44 頁 c、串聯(lián)型 圖 1 濾波器的性能指標(biāo) 我們在進(jìn)行濾波器設(shè)計(jì)時(shí)，需要確定其性能指標(biāo)。一般來說，濾波器的性能要求往往以頻率響應(yīng)的幅度特 性的允許誤差來表征。以低通濾波器特性為例，頻率響應(yīng)有通帶、過渡帶及阻帶三個(gè)范圍。 在通帶內(nèi)： 1 AP≤| H(ejω)| ≤1 |ω|≤ωc 在阻帶中： |H(ejω)| ≤ Ast ωst ≤|ω|≤ωc 其中 ωc 為通帶截止頻率 , ωst為阻帶截止頻率， Ap為通帶誤差 , Ast為阻帶誤差。 圖 低通濾波器頻率響應(yīng) 幅度 特性的容限圖 課程設(shè)計(jì) 第 10 頁 共 44 頁 與模擬濾波器類似，數(shù)字濾波器按頻率特性劃分為低通、高通、帶通、帶阻、全通等類型，由于數(shù)字濾波器的頻率響應(yīng)是周期性的，周期為 2π。各種理想數(shù) 字濾波器的幅度頻率響應(yīng)如圖所示： 圖 各種理想數(shù)字濾波器的幅度頻率響應(yīng) IIR 數(shù)字濾波器的設(shè)計(jì)方法 目前， IIR 數(shù)字濾波器設(shè)計(jì)最通用的方法是借助于模擬濾波器的設(shè)計(jì)方法。模擬濾波器設(shè)計(jì)已經(jīng)有了一套相當(dāng)成熟的方法，它不但有完整的設(shè)計(jì)公式，而且還有較為完整的圖表供查詢，因此，充分利用這些已有的資源將會(huì)給數(shù)字濾波器的設(shè)計(jì)帶來很大方便， IIR 數(shù)字濾波器的設(shè)計(jì)步驟是： 課程設(shè)計(jì) 第 11 頁 共 44 頁 (1)按一定規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬濾波器的技術(shù)指標(biāo)； (2)根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器 H(s)； (3)在按一定規(guī)則將 H(s)轉(zhuǎn)換為 H(z)。 若所設(shè)計(jì)的數(shù)字濾波器是低通的，那么上述設(shè)計(jì)工作可以結(jié)束，若所設(shè)計(jì)的是高通、帶通或者帶阻濾波器，那么還有步驟： (4)將高通、帶通或者帶阻數(shù)字濾波器的技術(shù)指標(biāo)先轉(zhuǎn)化為低通濾波器的技術(shù)指標(biāo)，然后按上述步驟 (2)設(shè)計(jì)出模擬低通濾波器 H(s)，再由沖擊響應(yīng)不變法或雙線性變換將 H(s)轉(zhuǎn)換為所需的 H(z)。 s z 映射的方法有 ： 沖激響應(yīng)不變法、階躍響應(yīng)不變法、雙線性變換法等。下面討論 脈沖 響應(yīng)不變法 。 脈沖響應(yīng)不變法是從濾波器的脈沖響應(yīng)出發(fā)，使數(shù)字濾波器的單位脈沖響應(yīng)序 h(n)模仿模擬濾波去的沖擊響應(yīng) ha(t)，使 h(n)正好等于 ha(t)的采樣值，即 h(n)=ha(nT) （ 23） T 為采樣周期。如以 Ha(s)及 H(z)分別表示 ha(t)的拉式變換及 h(n)的 z 變 換，即 Ha(s)=L[ha(t)] （ 24） H(z)=Z[h(n)] （ 25） 則根據(jù)采樣序列 z 變換與模擬信號拉式變換的關(guān)系，得： （ 26） 上式表明，采樣脈沖響應(yīng)不變法將模擬濾波器變換為數(shù)字濾波器時(shí)，它所完成的 S平面到 Z 平面的變換，正是以前討論的拉式變換到 Z 變換的標(biāo)準(zhǔn)變換，即首先對 Ha(s)作周期嚴(yán)拓，然后再經(jīng)過 z=est 的映射關(guān)系映射到 Z 平面上。 應(yīng)當(dāng)指出， Z=est 的映射關(guān)系表明， S 平面上每一條寬為 2pi/T 的橫帶部分，都將重疊課程設(shè)計(jì) 第 12 頁 共 44 頁 地映射到 Z 平面的整個(gè)全部平面上。每一橫帶的左半部分映射到 Z 平面單位圓以內(nèi)，每一橫帶的右半部分映射到 Z 平面單位圓以外 ， jΩ 軸映射在單位圓上，但 jΩ 軸上的每一段2pi/T 都應(yīng)于繞單 位圓一周 。 Z=est 的映射關(guān)系反映的是 Ha(s)的周期嚴(yán)拓與 H(z)的關(guān)系，而不是 Ha(s)b 本身與 H(z)的關(guān)系，因此，使用脈沖響應(yīng)不變法時(shí)，從 Ha(s)到 H(z)并沒有一個(gè)由 S 平面到 Z 平面的簡單代數(shù)映射關(guān)系，即沒有一個(gè) s=f(z)的代數(shù)關(guān)系式。另外，數(shù)字濾波器的頻響也不是簡單的重現(xiàn)模擬濾波器的頻響，而是模擬濾波器頻響的周期嚴(yán)拓，周期為 Ω S=2π /T=2π fs，即 （ 27） 巴特沃斯濾波器 巴特沃斯濾波器是電子濾波器的一種。巴特沃斯濾波器的特點(diǎn)是通頻帶的頻率響應(yīng)曲線最平滑。 巴特沃斯濾波器的特性是通頻帶內(nèi)的頻率響應(yīng)曲線最大限度平坦，沒有起伏，而在組頻帶則逐漸下降為零。在振幅的對數(shù)對角頻率的波得圖上，從某一邊界角頻率開始，振幅隨著角頻率的增加而逐步減少，趨向負(fù)無窮大。 其振幅平方函數(shù)具有如 21 式： （ 28） 式中， N 為整數(shù)，稱為濾波器的階數(shù)， N 越大，通帶和阻帶的近似性越好，過渡帶也越陡。課程設(shè)計(jì) 第 13 頁 共 44 頁 如下圖 所示： 圖 巴特沃茲振幅平方函數(shù) 通帶 ： 使信號通過的頻帶 ；阻帶：抑制噪聲通過的頻帶； 過渡帶：通帶→阻帶間過渡的頻率范圍 ； Ω c：截止頻率。 理想濾波器的過渡帶為Ω，阻帶 |H(jΩ )|=0，通帶內(nèi)幅度 |H(jΩ )|=常數(shù)， H(jΩ )線性相 位。通帶內(nèi)，分母 Ω /Ω c1,相應(yīng) (Ω /Ω c)2N 隨 N 的增加而趨于 0， A(Ω 2)→1 ， 在過渡帶和阻帶， Ω /Ω c1， 隨 N 的增加， Ω e/Ω c1，所以 A(Ω2)快速下降。 Ω =Ω c 時(shí)， ，幅度衰減 ，相當(dāng)于 3bd 衰減點(diǎn)。 振幅平方函數(shù)的極點(diǎn)可寫成如式 22： Ha(s).Ha(s)= （ 29） 課程設(shè)計(jì) 第 14 頁 共 44 頁 可分解為 2N 個(gè)一次因式令分母為零，→ 可見， Butterworth 濾波器的振幅平方函數(shù)有 2N 個(gè)極點(diǎn)，它們均勻?qū)ΨQ地分布在 |s|=Ω c 的圓周上。 第三章 IIR 濾波器的頻率響應(yīng) 實(shí)例 用脈沖響應(yīng)不變法設(shè)計(jì) IIR 低通數(shù)字濾波器實(shí)例。 數(shù)字低通的技術(shù)指標(biāo)為： Wp= rad Ap=1dB T=2s Ws= rad As=15dB 程序?yàn)椋? T=2。 fs=1/T。 %采樣頻率為 采樣周期倒數(shù) Wp=.*pi。 課程設(shè)計(jì) 第 15 頁 共 44 頁 Ws=.*pi。 %設(shè)計(jì)歸一化通帶阻帶截止頻率 Ap=1。 As=15。 %設(shè)置通帶最大最小衰減 [N,Wc]=buttord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用 butter函數(shù)確定巴特沃斯濾波器階數(shù) [B,A]=butter(N,Wc,39。s39。)。 %調(diào)用 butter函數(shù)設(shè)計(jì)巴特沃斯濾波器 W=linspace(0,pi,400*pi)。 %指定一段頻率值 [D,C]=impinvar(B,A,fs)。 %調(diào)用脈沖不變法 168。 Hz=freqz(D,C,W)。 %返回頻率值 166。 plot(W/pi,abs(Hz)/abs(Hz(1)))。 %繪出巴特沃斯數(shù)字低通濾波器的幅頻特性曲線 grid on。 title(39。巴特沃斯低通數(shù)字濾波器 39。)。 xlabel(39。Frequency/Hz39。)。 ylabel(39。Magnitude39。)。 得出幅頻特性如下圖： 課程設(shè)計(jì) 第 16 頁 共 44 頁 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 100 . 20 . 40 . 60 . 811 . 21 . 4巴特沃斯低通數(shù)字濾波器F r e q u e n c y / H zMagnitude 用脈沖響應(yīng)不變法設(shè)計(jì) IIR 高通數(shù)字濾波器實(shí)例。 數(shù)字 高 通的技術(shù)指標(biāo)為： Wp= rad Ap=2dB T=2s Ws= rad As=15dB 程序?yàn)椋? T=2。 %采樣周期 fs=1/T。 %采樣頻率 課程設(shè)計(jì) 第 17 頁 共 44 頁 Wp= .*pi。 Ws= .*pi。 %設(shè)置歸一化通帶和阻帶截止平率 Ap= 2。 As=15。 %設(shè)置通帶最大最小衰減 [N,Wc]=buttord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用 butter函數(shù)確定巴特沃斯濾波器階數(shù) [B,A]=butter(N,Wc, 39。high39。,39。s39。)。 %調(diào)用 butter函數(shù)設(shè)計(jì)巴特沃斯濾波器 W=linspace(0,pi,400*pi)。 %指定一段頻率值 [D,C]=impinvar(B,A,fs)。 %調(diào)用脈沖響應(yīng)不變法 168。 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) 166。 plot(W/pi,abs(Hz)/abs(Hz(1)))。 %繪制巴特沃斯數(shù)字高通濾波器的幅頻響應(yīng)曲線 grid on。 title(39。巴特沃斯高通數(shù)字濾波器 39。)。 xlabel(39。Frequency/Hz39。)。 ylabel(39。Magnitude39。)。 得出幅頻特性如下圖： 課程設(shè)計(jì) 第 18 頁 共 44 頁 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 111 . 0 51 . 11 . 1 51 . 21 . 2 51 . 31 . 3 51 . 4巴特沃斯高通數(shù)字濾波器F r e q u e n c y / H zMagnitude 用脈沖響應(yīng)不變法設(shè)計(jì) IIR 帶通數(shù)字濾波器實(shí)例。 帶通濾波器技術(shù)指標(biāo)為： Wp=[ ]rad Ap=1