nyquistzp.m

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

ng = input('"zero.eg:[-1 -2 -3]="');
dg = input('"polar.eg:[-1 -2 -3 -4]="');
k = input('"gain k.eg:2="');
G = zpk(ng,dg,k)
nyquist(G);grid;
text(-1,0,'(-1,j0)');
[mun,den]=zp2tf(ng',dg',k);
[mag,pha,w] = bode(G);
disp('"the roots of den="'); roots(den)
[gm,pm,wpc,wgc]=margin(mag,pha,w);
disp('[Wc  gama]');[wgc,pm]
disp('Compare in one graph, input "hold on"');