qbxf8.m

切比雪肤八阶滤波器在MATLAB中的实现

dt=0.0000001;
t=-0.0004*pi:dt:0.0004*pi
R1=10000;
R2=10000;
R3=10000;
R4=10000;
R5=10000;
R6=10000;
R7=10000;
R8=10000;
C1=4.7*10^(-9);
C2=4.7*10^(-9);
C3=4.7*10^(-9);
C4=4.7*10^(-9);
C5=10^(-9);
C6=2.7*10^(-9);
C7=2.2*10^(-9);
C8=10^(-9);
w=8150;
sys1=tf(1,[R1*R2*C2*C1,(R1*C2+R2*C2),1]);
sys2=tf(1,[R3*R4*C4*C3,(R3*C4+R4*C3),1]);
sys3=tf(1,[R5*R6*C6*C5,(R5*C6+R6*C6),1]);
sys4=tf(1,[R7*R8*C8*C7,(R7*C8+R8*C8),1]);
x=sin(w*t);
y=lsim(sys1,x,t);
y1=lsim(sys2,y,t);
y2=lsim(sys3,y1,t);
lsim(sys4,y2,t)
figure(2)
plot(t,x)