nyquist.m

控制系统仿真中的频率特性的MATLAB源码

%绘制根轨迹图

disp('开环传递函数G(s)为：')
printsys(num1,den1,'s');
disp('反馈传递函数H(s)为:')
printsys(num2,den2,'s');
lcwin=figure('Position',[120 49 513 491],'Name','nyquist曲线','NumberTitle','off',...
      'MenuBar','none');

axes('position',[0 0 1 1]);
axis off;




%设置绘图区域
axes('Position',[.15 .33 .8 .56]);

%绘制nyquist曲线
sys1=tf(num1,den1)
sys2=tf(num2,den2)
sys3=series(sys1,sys2)
nyquist(sys3)

grid on
title('nyquist曲线')
xlabel('实轴')
ylabel('虚轴')