【濾波器】基于脈沖響應(yīng)不變法+雙線性變換法實(shí)現(xiàn)數(shù)字濾波器設(shè)計(jì)含Matlab源碼
1 簡(jiǎn)介
1、設(shè)計(jì)目的
(1)掌握脈沖響應(yīng)不變法和雙線性變換法設(shè)計(jì)IIR數(shù)字濾波器的具體方法和原理,熟悉雙線性變換法和脈沖響應(yīng)不變法設(shè)計(jì)低通、帶通IIR數(shù)字濾波器的計(jì)算機(jī)編程。
(2)觀察雙線性變換法和脈沖響應(yīng)不變法設(shè)計(jì)的數(shù)字濾波器的頻域特性,了解雙線性變換法和脈沖響應(yīng)不變法的特點(diǎn)與區(qū)別。
(3)熟悉Butterworth濾波器、Chebyshev濾波器和橢圓濾波器的頻率特性。

2 部分代碼
%用雙線性變換法設(shè)計(jì)一個(gè)Butterworth數(shù)字低通濾波器
%確定已知的技術(shù)指標(biāo)
Rp=1;Rs=25;T=0.001;fp=200;fs=300;
%求出將要設(shè)計(jì)的數(shù)字濾波器的邊界頻率
wp=2*pi*fp*T;ws=2*pi*fs*T;
%預(yù)畸變
wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2);
%求出最小階數(shù)和截止頻率
[n,wn]=buttord(wp1,ws1,Rp,Rs,'s');
%求出模擬濾波器的分子和分母系數(shù)向量
[b,a]=butter(n,wn,'low','s');
%雙線性變換
[bz,az]=bilinear(b,a,1/T);
[db,mag,pha,grd,w]=freqz_m(bz,az);
%以w/pi為橫軸,以db為縱軸繪圖
plot(w/pi,db,'color','r'); grid;
title('Butterworth數(shù)字低通濾波器(雙線性法)');
%設(shè)定橫軸和縱軸的坐標(biāo)范圍
axis([0,1,-30,2]);
3 仿真結(jié)果





4 參考文獻(xiàn)
[1]王艷文, 崔志娟, 張靜. 基于matlab脈沖響應(yīng)不變法的Butterworth數(shù)字濾波器設(shè)計(jì)[J].? 2022(11).
[2]梁麗娟, and 安蘭珠. "利用MATLAB實(shí)現(xiàn)基于雙線性變換法的IIR數(shù)字濾波器設(shè)計(jì)." 邢臺(tái)學(xué)院學(xué)報(bào) 27.2(2012):2.
博主簡(jiǎn)介:擅長(zhǎng)智能優(yōu)化算法、神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)、信號(hào)處理、元胞自動(dòng)機(jī)、圖像處理、路徑規(guī)劃、無人機(jī)等多種領(lǐng)域的Matlab仿真,相關(guān)matlab代碼問題可私信交流。
部分理論引用網(wǎng)絡(luò)文獻(xiàn),若有侵權(quán)聯(lián)系博主刪除。
