matlab设计各类模拟滤波器

《matlab设计各类模拟滤波器》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、1一阶低通有源滤波器symsxsR=10000;Rf=2000;R1=1000;C=0.0001;%参数设置，方便绘图，使R*C为1K=(1+Rf/R1);H=K/(1+s*R*C);%转移函数，为一阶负反馈电路。特别注意含集成运放的放大电路都认为满足负反馈条件，（所以计算时才可以用“虚短”续断）figure(1),ezplot(H,[0,100*pi]);%利用roots绘制零极点图zs为零点，ps为极点b=[1];a=[R*C1];zs=roots(b);ps=roots(a)figure(2),plot(real(zs),imag(zs),'o'),holdon

2、plot(real(ps),imag(ps),'x');gridaxis([-33-33]);%建立坐标系2，一阶有源高通高通（对比一阶低通：输入端RC位置互换）symsxsR=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=[K0];a=[11/(R*C)];zs=roots(b)ps=roots(a)figure(2),plot(real(zs),imag(zs

3、),'o'),holdonplot(real(ps),imag(ps),'x');gridaxis([-33-33]);3，二阶有源低通滤波器高通symsxsR=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始终

4、未变）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');holdonplot(real(ps),imag(ps),'x');gridaxis([-1010-1010]);%坐标范围更大4，全通滤波器（应用为对不同频率信号相频特性不同）symsxsR=10000;H=(1-s*R

5、*C)/(1+s*R*C)figure(1),ezplot(H,[0,100*pi]);%转移函数b=[-R*C1];%绘制零极点图a=[R*C1];zs=roots(b);ps=roots(a)figure(2),plot(real(zs),imag(zs),'o'),holdonplot(real(ps),imag(ps),'x');gridaxis([-33-33]);5，5，二阶带阻有源滤波器symsxsR=10000;C=0.0001;H=(s.^2+(1/(R*C)).^2)/(s.^2+(2*s/3*R*C)+1/(3*R.^2*C.^2));%依电路图

6、求得的转移函数（深度负反馈条件下）%b=[10(1/M).^2];a=[12/(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');holdonplot(real(ps),imag(ps),'x');gridaxis([-66-66]);