经典滤波器的MATLAB仿真源程序汇.docVIP

1、%巴特沃斯低通模拟圆形滤波器 clear all; n=0:0.01:2; for i=1:4 switch i case 1 N=2; case 2 N=5; case 3 N=10; case 4 N=20; end [z,p,k]=buttap(N); %函数buttap--设计巴特沃斯低通滤波器 [b,a]=zp2tf(z,p,k); %函数zp2tf--零极点增益模型转换为传递函数模型 [H,w]=freqs(b,a,n); %函数freqs--求解模拟滤波器频率响应 magH2=(abs(H)).^2; %函数abs--取模值函数 hold on %函数hold--控制是否保持当前图形 plot(w,magH2) %函数plot--画二维线性图 axis([0 2 0 1]); %函数axis--控制坐标轴比例和外观 end xlabel(w/wc); ylabel(|H(jw)|^2); title(巴特沃斯低通模拟滤波器); text(0.72,0.63,N=2) %对不同曲线做标记 text(0.98,0.85,N=20) grid on; 2、%绘制切比雪夫I型低通模拟滤波器的平方幅频响应曲线，滤波器的阶数分别为2,4,6,8. clear all; n=0:0.01:2; for i=1:4 switch i case 1 N=2; case 2 N=4; case 3 N=6; case 4 N=8; end Rs=10; [z,p,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([0 2 0 1]); xlabel(w/wc); ylabel(H(jw)^2); title([N= num2str(N)]); grid on end 3、%切比雪夫II型低通模拟滤波器 clear all; n=0:0.01:2; for i=1:2 switch i case 1 N=7; case 2 N=8; end Rs=10; %阻带文波系数为10dB [z,p,k]=cheb2ap(N,Rs); %函数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([0 2 0 1.1]); xlabel(w/wc); ylabel(|H(jw)|^2); title([N= num2str(N)]); end 4、%运用冲击响应不变法设计一个低通Chebshev1型数字滤波器，其通带上限临界频率是3Hz，阻带临界频率是5H，采样频率是1000Hz，在通带内的最大衰减为0.3dB，阻带内的最小衰减为80dB。 MATLAB程序如下： clc; clear all; %把数字滤波器的频率特征转换成模拟滤波器的频率特征 wp=300*2*pi; ws=400*2*pi; rp=0.3; rs=80; Fs=1000; %选择滤波器的最小阶数。 [N,Wn]=cheb1ord(wp,ws,rp,rs,s); %创建Chebyshev1低通滤波器的原型 [Z,P,K]=cheb1ap(N,rp); [A,B,C,D]=zp2ss(Z,P,K); %实现低通向低通的转换 [AT,BT,CT,DT]=lp2lp(A,B,C,D,Wn); [num1,den1]=ss2tf(AT,BT,