home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 2.ddi / MUTOOLS1.DI$ / STATECC.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  835 b   |  31 lines
  1. % function sysout = statecc(sysin,t)
  2. %
  3. %   Applies a state coordinate transformation to the
  4. %   system matrix, yielding a new system matrix with:
  5. % 
  6. %       sysout = pck(inv(t)*a*t,inv(t)*b,c*t,d)
  7. %
  8. %   See also: EIG, INV, MINV, MMULT, REORDSYS, SRESID, 
  9. %             STRANS, and TRUNC.
  10.  
  11. function sysout = statecc(sysin,t)
  12. if nargin ~= 2
  13.   disp(['usage: sysout = statecc(sysin,t)']);
  14. else
  15.   [type,mrow,mcol,mnum] = minfo(sysin);
  16.   [nrt,nct] = size(t);
  17.   if nrt ~= mnum | nct ~= mnum
  18.     error(['incompatible dimensions'])
  19.     return
  20.   end
  21.   if strcmp(type,'syst')
  22.     [a,b,c,d] = unpck(sysin);
  23.     nab = t\[a b];
  24.     sysout = pck(nab(1:mnum,1:mnum)*t,nab(1:mnum,mnum+1:mnum+mcol),c*t,d);
  25.   else
  26.     disp(['sysin needs to be a SYSTEM matrix'])
  27.   end
  28. end
  29. %
  30. % Copyright MUSYN INC 1991,  All Rights Reserved
  31.