資源描述：
《經(jīng)典濾波器的matlab仿真源程序》由會員上傳分享，免費在線閱讀，更多相關內容在行業(yè)資料-天天文庫。
1、1、%巴特沃斯低通模擬圓形濾波器clearall;n=0:0.01:2;fori=1:4switchicase1N=2;case2N=5;case3N=10;case4N=20;end[z,p,k]=buttap(N);%函數(shù)buttap--設計巴特沃斯低通濾波器[b,a]=zp2tf(z,p,k);%函數(shù)zp2tf--零極點增益模型轉換為傳遞函數(shù)模型[H,w]=freqs(b,a,n);%函數(shù)freqs--求解模擬濾波器頻率響應magH2=(abs(H)).^2;%函數(shù)abs--取模值函數(shù)holdon%函數(shù)hold--控制是否保持當前圖形plot(w,magH2)%函數(shù)p
2、lot--畫二維線性圖axis([0201]);%函數(shù)axis--控制坐標軸比例和外觀endxlabel('w/wc');ylabel('
3、H(jw)
4、^2');title('巴特沃斯低通模擬濾波器');text(0.72,0.63,'N=2')%對不同曲線做標記text(0.98,0.85,'N=20')gridon;2、%繪制切比雪夫I型低通模擬濾波器的平方幅頻響應曲線，濾波器的階數(shù)分別為2,4,6,8.clearall;n=0:0.01:2;fori=1:4switchicase1N=2;case2N=4;case3N=6;case4N=8;endRs=10;[z,p
5、,k]=cheb1ap(N,Rs);[b,a]=zp2tf(z,p,k);[H,w]=freqs(b,a,n);magH2=(abs(H)).^2;posplot=['22'num2str(i)];subplot(posplot)plot(w,magH2)axis([0201]);xlabel('w/wc');ylabel('H(jw)^2');title(['N='num2str(N)]);gridonend3、%切比雪夫II型低通模擬濾波器clearall;n=0:0.01:2;fori=1:2switchicase1N=7;case2N=8;endRs=10;%阻帶文
6、波系數(shù)為10dB[z,p,k]=cheb2ap(N,Rs);%函數(shù)cheb2---設計切比雪夫II型低通濾波器[b,a]=zp2tf(z,p,k);[H,w]=freqs(b,a,n);magH2=(abs(H)).^2;%輸出圖形posplot=['12'num2str(i)];subplot(posplot)plot(w,magH2)axis([0201.1]);xlabel('w/wc');ylabel('
7、H(jw)
8、^2');title(['N='num2str(N)]);end4、%運用沖擊響應不變法設計一個低通Chebshev1型數(shù)字濾波器，其通帶上限臨界頻率
9、是3Hz，阻帶臨界頻率是5H，采樣頻率是1000Hz，在通帶內的最大衰減為0.3dB，阻帶內的最小衰減為80dB。MATLAB程序如下：clc;clearall;%把數(shù)字濾波器的頻率特征轉換成模擬濾波器的頻率特征wp=300*2*pi;ws=400*2*pi;rp=0.3;rs=80;Fs=1000;%選擇濾波器的最小階數(shù)。[N,Wn]=cheb1ord(wp,ws,rp,rs,'s');%創(chuàng)建Chebyshev1低通濾波器的原型[Z,P,K]=cheb1ap(N,rp);[A,B,C,D]=zp2ss(Z,P,K);%實現(xiàn)低通向低通的轉換[AT,BT,CT,DT]=lp2
10、lp(A,B,C,D,Wn);[num1,den1]=ss2tf(AT,BT,CT,DT);%運用沖擊響應不變法把模擬濾波器轉換成數(shù)字濾波器[num2,den2]=impinvar(num1,den1,1000);%繪出頻率響應曲線[H,W]=freqz(num2,den2);plot(W*Fs/(2*pi),abs(H));grid;xlabel('幅值');ylabel('頻率');title('沖擊響應不變法低通濾波器');clc;clearall;%把數(shù)字濾波器的頻率特征轉換成模擬濾波器的頻率特征wp=300*2*pi;ws=400*2*pi;rp=0.3;rs=8
11、0;Fs=1000;%選擇濾波器的最小階數(shù)。[N,Wn]=cheb1ord(wp,ws,rp,rs,'s');%創(chuàng)建Chebyshev1低通濾波器的原型[Z,P,K]=cheb1ap(N,rp);[A,B,C,D]=zp2ss(Z,P,K);%實現(xiàn)低通向低通的轉換[AT,BT,CT,DT]=lp2lp(A,B,C,D,Wn);[num1,den1]=ss2tf(AT,BT,CT,DT);%運用沖擊響應不變法把模擬濾波器轉換成數(shù)字濾波器[num2,den2]=impinvar(num1,den1,1000);%繪出頻
