我常用这个程序算环路滤波器,感觉差不多,感兴趣就是试试吧。
只是这程序是在scilab中运行的。用google艘一下scilab,这是个免费的软件。
程序:
s = poly(0,'s');
r1 = 100e3;
r2 = 270e3;
c = 0.1e-6;
kvco = 2*%pi*32;
kn = 1/2;
krssi = 0.6;
f = (s*c*r2+1)/s/c/r1;
// f = r2/r1/(s*c*r2+1);
g = kvco/s * kn * krssi * f;
xset("window",0);xbasc(0);xselect();
set("figure_style","old");
bode(syslin('c', g), 0.1, 9e3, .01);
180+p_margin(g)
h = g/.(1);
h = tf2ss(h, 1e-40);
t = 0:.001:0.2;
xset("window",1);xbasc(1);xselect();
plot2d([t',t'],[(csim('step',t,h))',ones(t')]);