matlab设计各类模拟滤波器.docVIP

1一阶低通有源滤波器 syms x s R=10000; Rf=2000;R1=1000;C=0.0001;%参数设置，方便绘图，使R*C为1 K=(1+Rf/R1); H=K/(1+s*R*C); %转移函数，为一阶负反馈电路。特别注意含集成运放的放大电路都认为满足负反馈条件，（所以计算时才可以用“虚短”续断） figure(1),ezplot(H,[0,100*pi]); %利用roots绘制零极点图 zs为零点，ps为极点 b=[1]; a=[R*C 1]; zs=roots(b); ps=roots(a) figure(2),plot(real(zs),imag(zs),o),hold on plot(real(ps),imag(ps),x); grid axis ([-3 3 -3 3]); %建立坐标系 2，一阶有源高通高通（对比一阶低通：输入端R C位置互换） syms x s R=10000; Rf=2000;R1=1000;C=0.0001; K=(1+Rf/R1); H=s*K/(s+1/(R*C)); %转移函数，注意一阶高通转移函数相比低通多一个零点 figure(1),ezplot(H,[0,100*pi]); %绘制零点极点图 b=[K 0]; a=[1 1/(R*C)]; zs=roots(b) ps=roots(a) figure(2),plot(real(zs),imag(zs),o),hold on plot(real(ps),imag(ps),x); grid axis ([-3 3 -3 3]); 3，二阶有源低通滤波器高通 syms x s R=1000; Rf=20000;R1=10000; C2=0.0001;C3=0.0001;R2=1000;R3=1000; K=(1+Rf/R1); H=(K/R2*R3*C2*C3)/(s.^2+s*(R3*C3+R2*C3-R2*C2+R2*C2-K*R2*C2)/R2*R3*C2*C3+1/R2*R3*C2*C3); %转移函数，（注意各电路K（即Af=1/B）未变，原因是反馈电路相同，则B=R1/R1+Rf始终未变） figure(1),ezplot(H,[0,100*pi]); b=[1]; a=[1 (R3*C3+R2*C3-R2*C2+R2*C2-K*R2*C2)/(R2*R3*C2*C3) 1/(R2*R3*C2*C3)]; zs=roots(b); %绘制零基点图 ps=roots(a) figure(2),plot(real(zs),imag(zs),o);hold on plot(real(ps),imag(ps),x); grid axis ([-10 10 -10 10]); %坐标范围更大 4，全通滤波器（应用为对不同频率信号相频特性不同） syms x s R=10000; H=(1-s*R*C)/(1+s*R*C) figure(1),ezplot(H,[0,100*pi]); %转移函数 b=[-R*C 1]; %绘制零极点图 a=[R*C 1]; zs=roots(b); ps=roots(a) figure(2),plot(real(zs),imag(zs),o),hold on plot(real(ps),imag(ps),x); grid axis ([-3 3 -3 3]); 5， 5，二阶带阻有源滤波器 syms x s R=10000; C=0.0001; H=(s.^2+(1/(R*C)).^2)/(s.^2+(2*s/3*R*C)+1/(3*R.^2*C.^2)); %依电路图求得的转移函数（深度负反馈条件下） % b=[1 0 (1/M).^2]; a=[1 2/(3*R*C) 1/(3*R.^2*C.^2)]; %%利用roots绘制零极点图 zs为零点，ps为极点 zs=roots(b) ps=roots(a) figure(2),plot(real(zs),imag(zs),o);hold on plot(real(ps),i