home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / CGLOCI2.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  1.9 KB  |  68 lines
  1. function [cg,ph] = cgloci2(Z1,Z2,Z3,Z4,Z5,Z6)
  2. %
  3. % Characteristic Gain/Phase Frequency Response.
  4. % 
  5. % [CG,PH] = CGLOCI2(SS_,Type,W) or
  6. % [CG,PH] = CGLOCI2(A, B, C, D, Type, W) produces the matrix CG & PH 
  7. % containing the characteristic Gain/Phase values for the system
  8. %                .
  9. %                x = Ax + Bu
  10. %                y = Cx + Du
  11. %                                                   -1
  12. %     with the frequency response G(jw) = C(jwI - A)  B + D.
  13. % CGLOCI calculates eigenvalues of one of the following depending on the 
  14. % value of Type: 
  15. %
  16. %     Type = 1   ----   G(jw)
  17. %     Type = 2   ----   inv(G(jw)) 
  18. %     Type = 3   ----   I + G(jw)
  19. %     Type = 4   ----   I + inv(G(jw)) 
  20. %
  21. % Vector W contains the frequencies at which the frequency response
  22. % is to be evaluated.  The CG and PH matrices have rows which 
  23. % correspond to the characteristic gain/phase in decending order.
  24.  
  25. % R. Y. Chiang & M. G. Safonov 6/29/87
  26. % Copyright (c) 1988 by the MathWorks, Inc.
  27. % All Rights Reserved.
  28.  
  29. inargs='(a,b,c,d,Type,w)';
  30. eval(mkargs(inargs,nargin,'ss'))
  31.  
  32. [mg] = freqrc(a,b,c,d,w);
  33. [rmg,cmg] = size(mg);
  34. [rb,cb] = size(b);
  35. [rc,cc] = size(c);
  36. gg = ones(rc,cb);
  37. for is = 1 : cmg
  38.   gg(:) = mg(:,is);
  39.   if (Type == 1)
  40.     char = eig(gg);
  41.     cg(:,is) = sort(abs(char)); 
  42.     ph(:,is) = sort(180./pi*imag(log(char))); 
  43.   end
  44.   if (Type == 2)
  45.     char = eig(inv(gg));
  46.     cg(:,is) = sort(abs(char));
  47.     ph(:,is) = sort(180./pi*imag(log(char)));
  48.   end
  49.   if (Type == 3)
  50.     [tmp1,tmp2]=size(gg);
  51.     char = eig(eye(tmp1,tmp2)+gg);
  52.     cg(:,is) = sort(abs(char)); 
  53.     ph(:,is) = sort(180./pi*imag(log(char))); 
  54.   end
  55.   if (Type == 4)
  56.     [tmp1,tmp2]=size(gg);
  57.     char = eig(eye(tmp1,tmp2)+inv(gg));
  58.     cg(:,is) = sort(abs(char)); 
  59.     ph(:,is) = sort(180./pi*imag(log(char))); 
  60.   end
  61. end
  62. cg = cg(cb:-1:1,:);
  63. ph = ph(cb:-1:1,:);
  64. %
  65. % ----- End of CGLOCI2.M ---- RYC/MGS 6/29/87 %
  66.  
  67.  
  68.