資源描述:
《IIR數(shù)字濾波器設(shè)計及實現(xiàn)》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、IIR數(shù)字濾波器設(shè)計及實現(xiàn)用雙線性法設(shè)計一個數(shù)字濾波器1編寫雙線性交換法設(shè)計batterworth低通濾波器的程序,要求通帶內(nèi)頻率低于0.2πrad1時,容許幅度誤差在1db之內(nèi),頻率在0.3πrad到πrad之間的阻帶衰減大于10db。2以π/64為采樣間隔,在屏幕上打印出數(shù)字濾波器的頻率區(qū)間[0,π]上的幅頻響應(yīng)特性曲線3在屏幕上打印出H(Z)的分子分母多項式的系數(shù)。源程序代碼為T=64/pi;%設(shè)置采樣周期為π/64fs=1/T;%采樣頻率為周期倒數(shù)wpz=0.2;wsz=0.3;%設(shè)置歸一化通帶
2、和阻帶截止頻率Wp=2*tan(wpz*pi/2);Ws=2*tan(wsz*pi/2);Ap=1;As=10;%設(shè)置通帶最大和最小衰減[N,Wc]=buttord(Wp,Ws,Ap,As,'s');%調(diào)用butter函數(shù)確定巴特沃斯濾波器階數(shù)[B,A]=butter(N,Wc,'s');%調(diào)用butter函數(shù)設(shè)計巴特沃斯濾波器W=linspace(0,pi,400*pi);%指定一段頻率值[D,C]=bilinear(B,A,fs);%調(diào)用雙線性變換法[Hz,W]=freqz(D,C,W);%返回頻率
3、響應(yīng)sys=tf(D,C,T)%給出傳輸函數(shù)H(z)subplot(3,1,1);plot(W,abs(Hz)/abs(Hz(1)));%繪出巴特沃斯數(shù)字低通濾波器的幅頻特性曲線gridon;title('巴特沃斯數(shù)字濾波器1');xlabel('Frequency/Hz');ylabel('Magnitude');subplot(3,1,2);plot(W,20*log10(abs(Hz)));gridon;title('巴特沃斯數(shù)字濾波器1');xlabel('Frequency/Hz');ylab
4、el('dB');[Nd,wdc]=buttord(wpz,wsz,Ap,As);[Bdz,Adz]=butter(Nd,wdc);[Hzz,W]=freqz(Bdz,Adz,W);subplot(3,1,3);plot(W,abs(Hzz)/abs(Hzz(1)));gridon;title('巴特沃斯數(shù)字濾波器2');xlabel('Frequency/Hz');ylabel('Magnitude');其結(jié)果為