rtloci.m

来自「基于ＭＡＴＬＡＢ的辅助设计 简述了矩阵分析的重要性」· M 代码 · 共 66 行

%
% RTLOCI   is the M-file used to plot the root locus of the
%          given system or systems.
%

%    Author:  Ole Barup Sorensen, Rapid Data Ltd 

%    Copyright (c) 1989-94 by Rapid Data Ltd
%    Revision 10:09  07/02/94

okflag=figflag(tit1,0);

if okflag==0,
  Fig14=figure('Numbertitle','off','Name',tit1,'Menubar','none',...
  'Units','normal','Position',[140/Scx 40/Scy 440/Scx 350/Scy]);
  set(0,'currentfig',Fig14);

% design menu
  pt_menu
  uimenu(Han_fil,'Label','Close','Sep','on','Callback',...
  'close(Fig14);clear Fig14');
end
delete(gca)

if (id>5)&(id~=10)&(id~=12)
  Amod=Amod*0;
  Amod(1)=1;
  id=1;
end

config

kinc=(kend-kbeg)/kgain;
k_range = kbeg: kinc: kend;
clear kinc

mmm = rlocus(num,den,k_range);
[p,z]=pzmap(num,den);

if id > 9
   ddd = rlocus(numx,denx,k_range);
   [px,zx]=pzmap(numx,denx);
end

if id <= 9
   plot(real(mmm),imag(mmm),'g*')
   hold on
   plot(real(z),imag(z),'og')
   plot(real(p),imag(p),'+g')
   title(' Root Locus')
   hold off
else
   plot(real(mmm),imag(mmm),'g*',real(ddd),imag(ddd),'r*')
   hold on
   plot(real(z),imag(z),'og')
   plot(real(p),imag(p),'+g')
   plot(real(zx),imag(zx),'or')
   plot(real(px),imag(px),'+r')
   title('Root Locus - Yellow: Compensated')
   hold off
end
xlabel('Real'), ylabel('Imag')
sgrid