nichols.m

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

%绘制根轨迹图

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

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




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

%绘制nichols曲线
sys1=tf(num1,den1)
sys2=tf(num2,den2)
sys3=series(sys1,sys2)
nichols(sys3)
grid on
title('nichols曲线')
xlabel('开环相角')
ylabel('开环增益')