home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 7.ddi / ROBUST.DI$ / SFR.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  846 b   |  33 lines
  1. function [am,bm,cm,dm] = sfr(Z1,Z2,Z3,Z4)
  2. % [SS_M] = SFR(SS_) or
  3. % [AM,BM,CM,DM] = SFR(A,B,C,D) produces a right spectral factorization
  4. %      via the duality of SFL such thtat
  5. %
  6. %                 M(s)M'(-s) = I - G(s)G'(-s)
  7. %
  8. %                Input data:  G(s):= SS_ = MKSYS(A,B,C,D);
  9. %                Output data: M(s):= SS_M = MKSYS(AM,BM,CM,DM);
  10. %
  11. %  The regular state-space can be recovered by "branch".
  12.  
  13. % R. Y. Chiang & M. G. Safonov 7/85
  14. % Copyright (c) 1988 by the MathWorks, Inc.
  15. % All Rights Reserved.
  16. % ---------------------------------------------------------------------
  17. %
  18.  
  19. inargs = '(a,b,c,d)';
  20. eval(mkargs(inargs,nargin,'ss'))
  21.  
  22. [amm,bmm,cmm,dmm] = sfl(a',c',b',d');
  23. am = amm';
  24. bm = cmm';
  25. cm = bmm';
  26. dm = dmm';
  27. %
  28. if xsflag
  29.    am = mksys(am,bm,cm,dm);
  30. end
  31. %
  32. % ------ End of SFR.M ---- RYC/MGS 7/85 %
  33.