home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 9.ddi / IDENT.DI$ / TRFSAUX.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  589 b   |  28 lines
  1. function [mag,phas] = trfsaux(a,b,c,d,iu,w,Tsamp)
  2. %TRFSAUX    Auxiliary function to TRFSS
  3. %
  4. %    [mag,phase] = trfsaux(a,b,c,d,ku,w,T)
  5.  
  6. %    L. Ljung 10-2-90
  7. %    Copyright (c) 1990 by the MathWorks, Inc.
  8. %    All Rights Reserved.
  9.  
  10. [no,ns] = size(c);
  11. nw = max(size(w));
  12.  
  13. [t,a] = balance(a);
  14. b = t \ b;
  15. c = c * t;
  16. [p,a] = hess(a);    
  17. b = p' * b(:,iu);
  18. c = c * p;
  19. d = d(:,iu);
  20. if Tsamp>0 w = exp(Tsamp*w * sqrt(-1)); else w = w * sqrt(-1);end
  21. g = ltifr(a,b,w);
  22. g = c * g + diag(d) * ones(no,nw);
  23. mag = abs(g)';
  24. [ny,nx]=size(c);
  25. for ky=1:ny
  26. phas(:,ky) = 180*phase(g(ky,:))'/pi;
  27. end
  28.