home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 2.ddi / MUTOOLS2.DI$ / HINF_SP.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  738 b   |  32 lines
  1. %
  2. %   split up the H-infinity matrices
  3. %
  4. function [ap,bp,cp,dp,b1,b2,c1,c2,d11,d12,d21,d22,ndata]...
  5.           = hinf_sp(p,nmeas,ncon)
  6. %
  7. [pptype,nr,nc,np] = minfo(p);
  8. if pptype ~= 'syst'
  9.   error('matrix is not a system matrix')
  10. end
  11. npout = nr;
  12. npin  = nc;
  13. np1 = npout-nmeas;
  14. np2 = nmeas;
  15. nm1 = npin-ncon;
  16. nm2 = ncon;
  17. ndata = [np1 np2 nm1 nm2];
  18. [ap,bp,cp,dp] = unpck(p);
  19. dum = sqrt(norm(bp,'inf')/norm(cp,'inf'));
  20. bp = bp/dum;
  21. cp = cp*dum;
  22. b1 = bp(:,[1:nm1]);
  23. b2 = bp(:,[(nm1+1):npin]);
  24. c1 = cp([1:np1],:);
  25. c2 = cp([(np1+1):npout],:);
  26. d11 = dp(1:np1,1:nm1);
  27. d12 = dp(1:np1,(nm1+1):npin);
  28. d21 = dp((np1+1):npout,1:nm1);
  29. d22 = dp((np1+1):npout,(nm1+1):npin);
  30. %
  31. % Copyright MUSYN INC 1991,  All Rights Reserved
  32.