rlocuszp.m

来自「这是一个用MATLAB工具设计自动控制理论中系统的仿真」· M 代码 · 共 18 行

ng=input('Zero,eg:[-1 -2 -3]:"');
dg=input('"Polar.eg:[-1 -2 -3 -4];"');
k=input('"Gain k.eg:1"');
G=zpk(ng,dg,k);
rlocus(G);
grid on ;
hold on;
[num,den]=zp2tf(ng',dg',k)
disp('"the roots of den "'); roots(den)
[mag,pha,w]=bode(G);
[gm,pm,wgc,wpc]=margin(mag,pha,w);
disp('"the value of [Wc,Gama]"');[wpc,pm]
rlocfind(G);
disp('the value of K');k
disp('Compare in one graph.input"hold on"');
disp('Get the value of K, input "zoom on"and "rlocfind(G)"');
disp('Get asymptote;input "rlocasymp"');