root.m

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

绘制根轨迹图

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

axes('position',[0 0 1 1]);
axis off;
dd=text(0.1,0.18,'稳定性');
set(dd,'fontsize',15,'color',[0 0 0],'fontname','宋体');
numedit=uicontrol(lcwin,'style','edit','units','normalized','position',[0.3 0.14 0.4 0.1],...
        'fontsize',15,'back',[1 1 1]);



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

%绘制根轨迹图
sys1=tf(num1,den1)
sys2=tf(num2,den2)
sys3=series(sys1,sys2)
sys=feedback(sys1,sys2)
rlocus(sys3)
Pole_sys=pole(sys)
if Pole_sys<0
    set(numedit,'string','稳定')
else
    set(numedit,'string','不稳定')
end

ylabel('虚     轴');
xlabel('实     轴');
title('根轨迹图');