MATLAB-IIR數(shù)字低通濾波器的設計方法

1、MATLAB IIR數(shù)字低通濾波器的設計方法1、IIR數(shù)字低通濾波器濾波器分兩大類：經(jīng)典濾波器和現(xiàn)代濾波器。經(jīng)典濾波的概念，是根據(jù)傅里葉分析和變換提出的一個工程概念。根據(jù)高等數(shù)學理論，任何一個滿足一定條件的信號，都可以被看成是由無限個正弦波疊加而成。換句話說，就是工程信號是不同頻率的正弦波線性疊加而成的，組成信號的不同頻率的正弦波叫做信號的頻率成分或叫做諧波成分?，F(xiàn)代濾波器理論研究的主要內(nèi)容是從含有噪聲的數(shù)據(jù)記錄又稱時間序列中估計出信號的某些特征或信號本身。經(jīng)典濾波器分為低通、高通、帶通、帶阻濾波器。每一種又有模擬濾波器AF和數(shù)字濾波器DF。對數(shù)字濾波器，又有IIR濾波器和FIR濾波器。我們要

2、實現(xiàn)的是IIR 數(shù)字濾波器。IIR數(shù)字濾波器采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。IIR濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成，可以組合成直接型、正準型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。由于運算中的舍入處理，使誤差不斷累積，有時會產(chǎn)生微弱的寄生振蕩。單位響應：IIR濾波器的單位脈沖響應為無限長，網(wǎng)絡中有反饋回路。幅頻特性：IIR數(shù)字濾波器幅頻特性精度很高，不是線性相位的，可以應用于對相位信息不敏感的音頻信號上IIR DF 的轉(zhuǎn)移函數(shù)是：IIR濾波器目前最通用的方法是利用已經(jīng)很成熟的模擬濾波器的設計方法進行設計。2、濾波器的技術(shù)要求低通濾波器：通帶截止頻率又稱通帶上限

3、頻率：阻帶下限截止頻率：通帶允許的最大衰減 ：阻帶允許的最小衰減，的單位dB：通帶上限角頻率：阻帶下限角頻率，即 3、IIR 數(shù)字濾波器的設計方法IIR數(shù)字濾波器的設計方法有兩類：間接設計法和直接設計法。間接設計法是借助模擬濾波器設計方法進行設計的，先根據(jù)數(shù)字濾波器設計指標設計相應的過渡模擬濾波器，再將過渡模擬濾波器轉(zhuǎn)換為數(shù)字濾波器。直接設計法師在時域或頻域直接設計數(shù)字濾波器。IIR數(shù)字濾波器在設計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設計數(shù)據(jù)或圖表可查，其設計工作量比較小，對計算工具的要求不高。在設計一個IIR數(shù)字濾波器時，我們根據(jù)指標先寫出模擬濾波器

4、的公式，然后通過一定的變換，將模擬濾波器的公式轉(zhuǎn)換成數(shù)字濾波器的公式。設計步驟如下：1按一定規(guī)則將給出的數(shù)字濾波器的技術(shù)指標轉(zhuǎn)換為模擬低通濾波器的技術(shù)指標。2根據(jù)轉(zhuǎn)換后的技術(shù)指標設計模擬低通濾波器；3再按一定的規(guī)則將轉(zhuǎn)換成。4假設是高通、帶通或帶阻數(shù)字濾波器則將它們的技術(shù)指標先轉(zhuǎn)化為低通模擬濾波器的技術(shù)指標，然后按上述步驟2設計出低通，再將轉(zhuǎn)換為所需的。4幾種不同類型的濾波器的介紹：因為我們設計的濾波器的沖擊響應一般都為實數(shù)，所以有這樣，如果我們能由，求出，那么就容易得到所需要的。不同類型的的表達式，代表了幾種不同類型的濾波器。1巴特沃思(Butterworth)濾波器：C為待定常數(shù)，N為待定

5、的濾波器階次。2切比雪夫I型(Chebyshev I )濾波器：5巴特沃思模擬低通濾波器的設計由于每一個濾波器的頻率范圍將直接取決于設計者的應用目的，因此必然是千差萬別。為了使設計標準化，我們需要將濾波器的頻率參數(shù)作歸一化處理。設所給的實際頻率為或，歸一化后的頻率為，對低通模擬濾波器，令顯然，。又令歸一化復數(shù)變量為，顯然所以巴特沃思模擬低通濾波器的設計可按以下三個步驟來進行。(1)將實際頻率規(guī)一化(2)求C和N這樣C和N可求。假設令3dB，則C1，這樣巴特沃思濾波器的設計就只剩一個參數(shù)N，這時(3)確定因為，根據(jù)上面公式有由解得，1，2，2N這樣可得求得后，用代替變量，即得實際需要得。6用雙線

6、性Z變換法設計IIR數(shù)字低通濾波器s平面到z平面的映射關(guān)系稱為雙線性Z變換，由此關(guān)系求出及即因為設計濾波器時系數(shù)會被約掉，所以又有7、MATLAB設計低通濾波器相關(guān)函數(shù) 1、Buttord.m本文件用來確定數(shù)字低通或模擬低通濾波器的階次，其調(diào)用格式分別是1N，Wn=buttord(Wp，Ws，Rp，Rs)2N，Wn=buttord(Wp，Ws，Rp，Rs，s)格式1)對應數(shù)字濾波器，式中Wp,Ws分別是通帶和阻帶的截止頻率，實際上它們是歸一化頻率，其值在01之間，1對應抽樣頻率的一半。對低通和高通濾波器，Wp,Ws都是標量，對帶通和帶阻濾波器，Wp,Ws都是12的向量。Rp,Rs分別是通帶和阻

7、帶的衰減，單位為dB。N是求出的相應低通濾波器的階次，Wn是求出的3dB頻率，它和Wp稍有不同。格式2對應模擬濾波器，式中各個變量的含義和格式1)相同，但Wp,Ws及Wn的單位為rad/s ，因此，它們實際上式頻率 。2、Buttap.m本文件用來設計模擬低通原型濾波器 ，其調(diào)用格式是z，p，k=buttap(N)N是欲設計的低通原型濾波器的階次，z，p，k分別是設計出的 的極點、零點及增益。3、Lp2lp.m4、Lp2hp.m5、Lp2bp.m6、Lp2bs.m從文件名可以看出，上述4個文件的功能分別是將模擬低通原型濾波器 轉(zhuǎn)換為實際的低通、高通、帶通及帶阻濾波器，其調(diào)用格式分別為：1B，A

8、lp2lpb，a，Wo或B，Alp2hpb，a，Wo2B，Alp2bpb，a，Wo或B，Alp2bsb，a，Wo式中b，a分別是模擬低通原型濾波器 的分子、分母多項式的系數(shù)向量，B，A分別是轉(zhuǎn)換后的 的分子、分母多項式的系數(shù)向量；在1中，Wo是低通或高通濾波器的截止頻率；在2中，Wo是帶通或帶阻濾波器的中心頻率，Bw是其帶寬。7、Bilinear.m本文件實現(xiàn)雙線性變換，即由模擬濾波器 得到數(shù)字濾波器 ，而s和z的關(guān)系由 給出。其調(diào)用格式是：Bz，Az bilinearB，A，F(xiàn)s 式中B，A分別是 的分子、分母多項式的系數(shù)向量，Bz，Az分別是 的分子、分母多項式的系數(shù)向量， Fs是抽樣頻率

9、。8、Freqz.m本文件用來在已知 , 的情況下求出系統(tǒng)的頻率響應 ，格式是H，wfreqzb，a，N，whole，F(xiàn)s其中N是頻率軸的分點數(shù)，建議N為2的整次冪；w 是返回頻率軸坐標向量，供繪圖用；Fs是抽樣頻率，假設Fs1，頻率軸給出歸一化頻率；whole 指定計算的頻率范圍是從0Fs，缺省時是從0Fs/29、b,a=zp2tf(z,p,k) 零極點型傳遞函數(shù)轉(zhuǎn)換到一般傳遞函數(shù)，它用于在極零點已知時求出 和 的系數(shù)。8、設計一個IIR 數(shù)字濾波器 技術(shù)指標：Hz，Hz，dB， dB，抽樣頻率Hz。 設計步驟：（1） 角頻率 （2） 對角頻率做預畸變 3模擬低通濾波器的階次，利用函數(shù)N，W

10、n=buttord(Wp，Ws，Rp，Rs，s) 注意：Wp，Ws應該為2中的，。4設計模擬低通原型濾波器，其調(diào)用格式是 z，p，k=buttap(N)。 N是欲設計的低通原型濾波器的階次，z，p，k分別是設計出的的極點、零點及增益。5求模擬低通原型濾波器的分子分母系數(shù)，b,a=zp2tf(z,p,k)。6求出的分子、分數(shù)系數(shù)。B，Alp2lpb，a，Wo7求出的分子、分母系數(shù)，利用bilinear函數(shù)。8求頻率響應，利用Freqz函數(shù)。9畫的圖 。程序設計：fp=100;fs=300;ap=3;%通帶最大衰減as=20;%阻帶最小衰減Fs1=1000;%抽樣頻率wp=2*pi*fp/Fs1ws=2*pi*fs/Fs1 %數(shù)字角頻率Fs=Fs1/Fs1;T=1/Fs;Op=2/T*tan(wp/2)Os=2/T*tan(ws/2)%角頻率畸變N,Wn=buttord(Op,Os,ap,as,s)z,p,k=buttap(N)%G(p)的零點,極點,增益b,a=zp2tf(z,p,k);%G(p)的分子、分數(shù)系數(shù)B,A=lp2lp(b,a,Op);%H(s)的分子、分數(shù)