www.pudn.com > mrascodes.rar > LypRoots.m, change:2003-01-29,size:502b


ac1=0.2+2*s1.signals.values; 
ac2=3+2*s0.signals.values; 
C=[ones(length(ac1),1) ac1 ac2]; 
 
hold; 
for m=2:length(ac1) 
   rplace=roots(C(m-1,:)); 
   rplace2=roots(C(m,:)); 
   L1X=[real(rplace(1)) real(rplace2(1))]; 
   L1Y=[imag(rplace(1)) imag(rplace2(1))]; 
    
   L2X=[real(rplace(2)) real(rplace2(2))]; 
   L2Y=[imag(rplace(2)) imag(rplace2(2))]; 
    
   line(L1X,L1Y); 
   line(L2X,L2Y); 
    
   %pause; 
   if mod(m,50) == 0 
   %   pause; 
   end 
    
%   plot(rplace,'k.'); 
end