ex4_1.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 12 行

% ex4_1
% 绘制奈魁斯特图
figure('pos',[30 50 260 200],'color','w');
axes('pos',[0.15 0.2 0.7 0.7]);
G1=tf([1],[10 1]);
w=[0 logspace(-2,3,200)];
G1v=polyval(G1.num{1},w*j)./polyval(G1.den{1},w*j);
G2=zpk([],[-1+0.2j -1-0.2j],1);
G2t=tf(G2);
G2v=polyval(G2t.num{1},w*j)./polyval(G2t.den{1},w*j);
plot(real(G1v),imag(G1v),'k',real(G2v),imag(G2v),'k');