example7_1a.m

来自「这是一个MATLAB6。5 辅助控制系统设计与仿真例题的所有源码」· M 代码 · 共 16 行

s1=-6.0001+15.6492i;
ng=1;dg=[1 20 75 0];
ngv=polyval(ng,s1);dgv=polyval(dg,s1);   %多项式求值
g=ngv/dgv;zeta=angle(g);
if zeta>0;phic=pi-zeta;end;
if zeta<0;phic=-zeta;end;
phi=angle(s1);zetaz=(phi+phic)/2;zetap=(phi-phic)/2;
zc=real(s1)-imag(s1)/tan(zetaz);  %计算校正器零、极点
pc=real(s1)-imag(s1)/tan(zetap);
nc=[1 -zc];dc=[1 -pc];
nv=polyval(nc,s1);dv=polyval(dc,s1);kv=nv/dv;
kc=abs(1/(g*kv));    %确定校正器增益
if zeta<0;kc=-kc;end;
kc
Gc=tf(nc,dc)