Index: mttroot/mtt/bin/trans/m/dm2fr.m ================================================================== --- mttroot/mtt/bin/trans/m/dm2fr.m +++ mttroot/mtt/bin/trans/m/dm2fr.m @@ -8,10 +8,13 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% $Id$ % %% $Log$ +% %% Revision 1.5 1996/08/16 14:26:37 peter +% %% Check and fix size of u0. +% %% % %% Revision 1.4 1996/08/15 12:50:51 peter % %% Put in a conj to undo effect of transpose. % %% % %% Revision 1.3 1996/08/15 11:53:44 peter % %% Now has u0 input vector @@ -52,11 +55,11 @@ fr = zeros(N,Ny); i = 0; for w = W' i = i+1; - FR = C*( (E*j*w - A) \ B*u0 ) + D*u0 + FR = C*( (E*j*w - A) \ B*u0 ) + D*u0; fr(i,:) = conj(FR'); end;