Documentation of trfsaux

Global Index (short | long) | Local contents | Local Index (short | long)

Function Synopsis

[mag,phas] = trfsaux(a,b,c,d,iu,w,Tsamp)

Help text

TRFSAUX        Auxiliary function to TRFSS

   [mag,phase] = trfsaux(a,b,c,d,ku,w,T)

Cross-Reference Information

Listing of function trfsaux

function [mag,phas] = trfsaux(a,b,c,d,iu,w,Tsamp)

%   L. Ljung 10-2-90
%   Copyright (c) 1986-98 by The MathWorks, Inc.
%   $Revision: 2.3 $  $Date: 1997/12/02 03:40:46 $

[no,ns] = size(c);
nw = length(w);

[t,a] = balance(a);
b = t \ b;
c = c * t;
[p,a] = hess(a);
b = p' * b(:,iu);
c = c * p;
d = d(:,iu);
if Tsamp>0 w = exp(Tsamp*w * sqrt(-1)); else w = w * sqrt(-1);end
g = ltifr(a,b,w);
g = c * g + diag(d) * ones(no,nw);
mag = abs(g)';
[ny,nx]=size(c);
for ky=1:ny
phas(:,ky) = 180*phase(g(ky,:))'/pi;
end