home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / HKSV.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  517 b   |  18 lines
  1. function [hsv,p,q] = hksv(a,b,c)
  2. %
  3. % [HSV,P,Q] = HKSV(A,B,C) computes reachability and observability
  4. %    grammians P, Q, and the Hankel singular values "hsv".
  5. %
  6.  
  7. % R. Y. Chiang & M. G. Safonov 3/86
  8. % Copyright (c) 1988 by the MathWorks, Inc.
  9. % All Rights Reserved.
  10. % ------------------------------------------------------------------
  11. p = gram(a,b);
  12. q = gram(a',c');
  13. hsv = sqrt(eig(p*q));
  14. [hsv,index] = sort(real(hsv));
  15. [ma,na] = size(a);
  16. hsv = hsv(ma:-1:1,:);
  17. %
  18. % -------- End of HKSV.M --- RYC/MGS