用雙線性變換法設(shè)計(jì)IIR濾波器

2 IIR濾波器的MATLAB實(shí)現(xiàn)
模擬型濾波器的設(shè)計(jì)過(guò)程，步驟繁瑣，計(jì)算量大。在MATLAB中，我們可以直接根據(jù)所給定的數(shù)字濾波器指標(biāo)調(diào)用特定的函數(shù)來(lái)確定濾波器的系統(tǒng)函數(shù)。具體實(shí)現(xiàn)過(guò)程如下(程序的主要部分，且只返回?cái)?shù)字濾波器的階數(shù)和截止
頻率)：
為調(diào)用階數(shù)選擇函數(shù)，下面的指標(biāo)已經(jīng)把Wp和Ws對(duì)π進(jìn)行了歸一化，且函數(shù)中返回的頻率也是歸一化頻率。

同理可實(shí)現(xiàn)ChebyshevI、ChebyshevII、ellipSe數(shù)字濾波器。Butterworth型和ellipke型系統(tǒng)的幅頻和相頻特性如圖l；ChebyshevI與ChebyshevII系統(tǒng)的幅頻和相頻特性如圖2(頻率為對(duì)π歸一化頻率)。

3 結(jié)論
從圖中我們可以看出，ButterWorth濾波器的頻率特性是平坦的，沒(méi)有出現(xiàn)波動(dòng)，即幅度隨著頻率單調(diào)遞減，它的相頻特性接近于分段特性。而橢圓濾波器幅頻特性存在著波紋(圖中由于坐標(biāo)比例的關(guān)系，只能畫出通帶波紋，阻帶波紋則無(wú)法顯示)；同時(shí)，它不具備線性相位特性。兩種契比雪夫?yàn)V波器的頻率特性如圖所示。從中可知道，這兩種濾波器的幅頻特性都有著波動(dòng)(其中，由于CII型濾波器是在阻帶中出現(xiàn)等波動(dòng)，因此無(wú)法表示出來(lái))，它們不具備線性相位特性。
為了比較這四種濾波器的性能高低，我們比較了實(shí)現(xiàn)相同的指標(biāo)所需的最低濾波器階數(shù)：Nbutt=14；Nchebl=8；Ncheb2=8；Nelli=6。
這樣可以清楚地看出：ButterWorth濾波器需要最高的階數(shù)14；兩種契比雪夫?yàn)V波器所需要的階數(shù)是一樣的8階，說(shuō)明它們具有相同的性能；而橢圓濾波器具有最小的階數(shù)4，說(shuō)明它的性能在這個(gè)意義上是最優(yōu)的。