nyquistpoly.m

这是一个用MATLAB工具设计自动控制理论中系统的仿真

num = input('"input num.eg:[0 0 1]="');
den = input('"input den.eg:[1 2 3]="');
G = tf(num,den)
nyquist(G);grid;
text(-1,0,'(-1,j0)');
[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"');