二阶无源RC滤波的MATLAB仿真

在网上搜了很长时间都没找到二阶无源RC滤波参数的MATLAB仿真，没办法，自己编了。二阶RC滤波回路的分析可采用拉式变换法进行比较简单，令1/CS代替C即得容抗，通过分压电路分析即得二阶RC滤波电路的传递函数，另wj代替S，即得二阶RC滤波电路的幅频特性。上图：（无法上图，可惜）

%本文件模拟二阶无源RC滤波的幅频特性
%figure上面的data cursor工具。在data cursor模式下，右击，再选择section style里的mouse
%position。拖动data cursor时，显示的就是鼠标所在点的位置
%数字滤波器的设计使用FDATOOL工具
%本文件利用abs函数自动计算复数的幅值
%对不同滤波器，只需要修改传递函数H即可
clc
clear
C=0.0022*1e-6; %UF
R=10*1e3;%10k
R1=4.3*1e3;%k
R2=10*1e0%
C1=5.6*1e-9; %UF
C2=3.3*1e-9; %UF

Wn=1/sqrt(R1*C1*R2*C2);
f=0:0.1:5000; %
W=2*pi*f;
H=zeros(size(W));
for i=1:1:length(W)-1
% H(i)=1/(1-(W(i)/Wn)^2+j*3*(W(i)/Wn));
H(i)=1/(1-(R1*R2*C1*C2)*W(i)*W(i)+(R1*C1+R2*C1+R1*C2)*W(i)*j);
end
A=20*log10(abs(H));
% axis([0 35 100 600]);grid;
plot(f,A,'b:o')