李凈 楊振平

【摘要】濾波器在實(shí)際信號(hào)處理中是非常重要的。任何被檢測(cè)的信號(hào)都含有噪聲，而濾波器是去除這種噪聲的。數(shù)字濾波器是具有一定傳輸選擇特性的數(shù)字信號(hào)處理裝置。數(shù)字信號(hào)的輸入與輸出是一個(gè)線性時(shí)不變的離散系統(tǒng)。它是利用離散系統(tǒng)特性對(duì)輸入信號(hào)進(jìn)行加工處理，讓有用信號(hào)通過(guò)。數(shù)字濾波器和模擬濾波器都是用來(lái)濾波的，根據(jù)頻率響應(yīng)特性可分為低通、高通、帶通、帶阻等類型；數(shù)字濾波器按單位脈沖響應(yīng)的性質(zhì)可分為無(wú)限長(zhǎng)單位脈沖響應(yīng)濾波器IIR和有限長(zhǎng)單位脈沖響應(yīng)濾波器FIR兩種。本文主要介紹IIR數(shù)字濾波器的設(shè)計(jì)。

【關(guān)鍵詞】數(shù)字濾波器 IIR 沖激響應(yīng)不變法 雙線性變換法

一、引言：

隨著信息學(xué)科的快速發(fā)展，數(shù)字信號(hào)處理技術(shù)正在飛速發(fā)展，并且滲透到其它技術(shù)中。數(shù)字濾波技術(shù)是數(shù)字信號(hào)分析、處理技術(shù)的重要分支。它對(duì)信號(hào)安全可靠和有效靈活的傳輸是至關(guān)重要的。

二、IIR數(shù)字濾波器的設(shè)計(jì)方法

一個(gè)數(shù)字濾波器可以用系統(tǒng)函數(shù)表示為：

由這樣的系統(tǒng)函數(shù)可以得到表示系統(tǒng)輸入與輸出關(guān)系的常系數(shù)線性差分方程為：

由此可知數(shù)字濾波器的功能就是把輸入序列x（n）通過(guò)一定的運(yùn)算變換成輸出序列y（n）。不同的處理方法決定了濾波器實(shí)現(xiàn)結(jié)構(gòu)的不同。無(wú)限沖激響應(yīng)濾波器的單位抽樣響應(yīng)h（n）是無(wú)限長(zhǎng)的，其差分方程是遞歸型的，在結(jié)構(gòu)上存在著輸出信號(hào)到輸入信號(hào)的反饋，所以在z平面的有限區(qū)間（0<︱z︱<∞）有極點(diǎn)存在。

目前，IIR數(shù)字濾波器的設(shè)計(jì)最通用的方法是借助于模擬濾波器的設(shè)計(jì)方法。模擬濾波器設(shè)計(jì)已經(jīng)相當(dāng)成熟了，給數(shù)字濾波器的設(shè)計(jì)帶來(lái)了很大的方便。

IIR數(shù)字濾波器的設(shè)計(jì)步驟如下（間接法）：

（1） 按任務(wù)需要將給出的數(shù)字濾波器的性能指標(biāo)轉(zhuǎn)換為模擬濾波器的技術(shù)指標(biāo)；

（2） 根據(jù)轉(zhuǎn)換后的性能指標(biāo)設(shè)計(jì)模擬低通濾波器H（s）；

（3） 再按照一定的規(guī)則將H（s）轉(zhuǎn)換為H（z）。

（4） 若所設(shè)計(jì)的數(shù)字濾波器是低通的，那么上述設(shè)計(jì)工作就可以結(jié)束了，若所設(shè)計(jì)的是高通，帶通或帶阻濾波器，那么需要把它們的性能指標(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)不變法設(shè)計(jì)數(shù)字濾波器就是使單位樣值響應(yīng)h（n）與相應(yīng)的模擬濾波器的沖激響應(yīng) 在抽樣點(diǎn)的量值相等，即 其中T為抽樣周期。

因此用沖激響應(yīng)不變法設(shè)計(jì)IIR數(shù)字濾波器就得首先根據(jù)設(shè)計(jì)要求來(lái)確定相應(yīng)的模擬濾波器的傳遞函數(shù) ，經(jīng)過(guò)拉普拉斯的反變換求出沖激響應(yīng) 后，對(duì)它進(jìn)行抽樣得到的 等于數(shù)字濾波器的單位樣值響應(yīng)h（n），在經(jīng)過(guò)z變換得到的H（z）就是數(shù)字濾波器的傳遞函數(shù)。如果模擬濾波器的 的分母階次大于分子的階次，則可將 展成部分分式形式

那么經(jīng)過(guò)拉普拉斯反變換求出的模擬濾波器的沖激響應(yīng) 則為

相應(yīng)的數(shù)字濾波器的單位樣值響應(yīng)為

對(duì)h（n）進(jìn)行z變換，得到數(shù)字濾波器的系統(tǒng)函數(shù)

由上述推導(dǎo)公式可知，只要模擬濾波器的傳遞函數(shù) 的N個(gè)極點(diǎn) 都是單極點(diǎn)，當(dāng)已經(jīng)求出各個(gè)極點(diǎn)的值 和部分分式的系數(shù) 后，則可以從模擬濾波器的傳遞函數(shù)的表達(dá)式（ ）直接得到數(shù)字濾波器的傳遞函數(shù)H（z）的表達(dá)式（ ）。

雙線性變換法是使數(shù)字濾波器的頻率響應(yīng)與模擬濾波器的頻率響應(yīng)相似的一種變換，它將s平面整個(gè)變換到一個(gè)中介的s1平面的一個(gè)水平窄帶 Ω：-π/T到

π/T之中，然后經(jīng)過(guò) 的變換，將s1平面的映射到z平面。這樣s平面就與s平面一 一對(duì)應(yīng)，消除了多值變換性，避免出現(xiàn)混疊現(xiàn)象。

S平面到z平面的變換可利用

令 jΩ=s， jΩ1=s1得：

從s1平面到z平面的變換即 ，代入上式，得到：

為使模擬濾波器的某一頻率與數(shù)字濾波器的任一頻率有對(duì)應(yīng)關(guān)系，可引入待定系數(shù)c，則Ω=c*tan（Ω1T/2）由此得到

這種s平面與z平面之間的單值映射關(guān)系就是雙線性變換，模擬濾波器的數(shù)字化只須用上述方法置換。

沖激響應(yīng)不變法與雙線性變換法的比較分析：

沖激響應(yīng)不變法的模擬頻率到數(shù)字頻率的轉(zhuǎn)換是線性的但是它在一定條件下會(huì)出現(xiàn)混疊失真，雙線性變換法雖然可以克服多值映射的關(guān)系和消除頻率混疊的現(xiàn)象，但是它是非線性的，在高頻處有較大的失真。所以沖激響應(yīng)不變法的數(shù)字化和雙線性變換法的數(shù)字化要根據(jù)實(shí)際情況運(yùn)用。

三、結(jié)束語(yǔ)：

本文主要介紹了IIR濾波器的性能指標(biāo)，設(shè)計(jì)方法和特點(diǎn)。通常IIR濾波器的設(shè)計(jì)過(guò)程如下：首先根據(jù)濾波器的參數(shù)要求設(shè)計(jì)對(duì)應(yīng)的濾波器（如巴特沃斯濾波器、切比雪夫?yàn)V波器等等），然后通過(guò)映射（如沖激響應(yīng)不變法、雙線性變換法等等）將模擬濾波器變換為數(shù)字濾波器，從而確定IIR濾波器的參數(shù)。IIR數(shù)字濾波器的最大缺點(diǎn)在于反饋的穩(wěn)定性不能得到保證。

【參考文獻(xiàn)】

[1]吳大正，信號(hào)與線性系統(tǒng)分析（第四版）【M】.高等教育出版社

[2]高西全，丁玉美，數(shù)字信號(hào)處理教程（第三版）?！綧】西安電子科技大學(xué)出版社

[3]程佩青，數(shù)字信號(hào)處理教程（第二版）。【M】清華大學(xué)出版社，2001（8）