基于Matlab的IIR數(shù)字濾波器設(shè)計(jì)方法比較及應(yīng)用
作者：
時(shí)間：2010-08-05
來(lái)源：網(wǎng)絡(luò)
收藏
加入技術(shù)交流群
掃碼加入
和技術(shù)大咖面對(duì)面交流
海量資料庫(kù)查詢(xún)
數(shù)字濾波器(Digital Filter)是指輸入、輸出都是離散時(shí)間信號(hào)，通過(guò)一定運(yùn)算關(guān)系改變輸入信號(hào)所含頻率成分的相對(duì)比例或者濾除某些頻率成分的器件。數(shù)字濾波器在數(shù)字信號(hào)處理中起著非常重要的作用，在信號(hào)的過(guò)濾、檢測(cè)與參數(shù)的估計(jì)等方面，是使用最為廣泛的一種線(xiàn)性系統(tǒng)。
實(shí)現(xiàn)數(shù)字濾波器的方法有兩種，一是采用計(jì)算機(jī)軟件進(jìn)行，就是把所要完成的工作通過(guò)程序讓計(jì)算機(jī)來(lái)實(shí)現(xiàn)；二是設(shè)計(jì)專(zhuān)用的數(shù)字處理硬件。這個(gè)地方主要用到的就是第一種方法。即是用Mafiab提供的信號(hào)處理工具箱來(lái)實(shí)現(xiàn)數(shù)字濾波器。
Matlab信號(hào)處理工具箱提供了豐富的設(shè)計(jì)方法，可以使得繁瑣的程序設(shè)計(jì)簡(jiǎn)化成函數(shù)的調(diào)用，只要以正確的指標(biāo)參數(shù)調(diào)用函數(shù)，就可以正確快捷地得到設(shè)計(jì)結(jié)果。
1 IIR數(shù)字濾波器設(shè)計(jì)方法比較
1．1 概述
IIR數(shù)字濾波器的最通用的方法是借助于模擬濾波器的設(shè)計(jì)方法。模擬濾波器設(shè)計(jì)已經(jīng)有了相當(dāng)成熟的技術(shù)和方法，有完整的設(shè)計(jì)公式，還有比較完整的圖表可以查詢(xún)，因此設(shè)計(jì)數(shù)字濾波器可以充分利用這些豐富的資源來(lái)進(jìn)行。
對(duì)于IIR數(shù)字濾波器的設(shè)計(jì)具體步驟如下：
(1)按照一定的規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬低通濾波器的技術(shù)指標(biāo)。
(2)根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器G(s)(G(s)是低通濾波器的傳遞函數(shù))。
(3)再按照一定的規(guī)則將G(s)轉(zhuǎn)換成H(z)(H(z)是數(shù)字濾波器的傳遞函數(shù))。若設(shè)計(jì)的數(shù)字濾波器是低通的，上述的過(guò)程可以結(jié)束，若設(shè)計(jì)的是高通、帶通或者是帶阻濾波器，那么還需要下面的步驟：
將高通、帶通或帶阻數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為低通模擬濾波器的技術(shù)指標(biāo)，然后設(shè)計(jì)出低通G(s)，再將G(s)轉(zhuǎn)換為H(z)。
Matlab信號(hào)工具箱提供了幾個(gè)直接設(shè)計(jì)IIR數(shù)字濾波器的函數(shù)，直接調(diào)用這些函數(shù)就可以很方便地對(duì)濾波器進(jìn)行設(shè)計(jì)。這里選取巴特沃斯法、切比雪夫I、切比雪夫Ⅱ、橢圓法四種方法進(jìn)行比較。給出用上述方法設(shè)計(jì)數(shù)字濾波器的函數(shù)如下：
其中：Wp表示通帶截止頻率；Ws表示阻帶截止頻率；Rp表示通帶紋波系數(shù)；Rs表示阻帶紋波系數(shù)；N表示濾波器最小階數(shù)；Wn表示截止頻率。b，a分別表示階次為N+1的數(shù)字濾波器系統(tǒng)傳遞函數(shù)的分子和分母多項(xiàng)式系數(shù)向量；Fs為采樣頻率；n為在區(qū)間[O Fs]頻率范圍內(nèi)選取的頻率點(diǎn)數(shù)；f記錄頻率點(diǎn)數(shù)。n取2的冪次方，可以提高運(yùn)算的速度，因?yàn)閒reqz函數(shù)采用基2的FFT算法。ftype=high時(shí)，為高通濾波器；ftype=bandpass時(shí)，為帶通濾波器；ftype=stop時(shí)，為帶阻濾波器。
1．2 四種設(shè)計(jì)方法比較
(1)在低通濾波器中的比較
假如：Wp=20 Hz，Ws=50 Hz，F(xiàn)s=200，Rp=1 dB，Rs=30 dB，分別用Butterworth低通濾波器、Chebyshev I型低通濾波器、ChebyshevⅡ型低通濾波器、橢圓低通濾波器四種方法進(jìn)行設(shè)計(jì)，如圖1所示。
評(píng)論