Index: mttroot/mtt/bin/trans/dm2fr_m ================================================================== --- mttroot/mtt/bin/trans/dm2fr_m +++ mttroot/mtt/bin/trans/dm2fr_m @@ -13,10 +13,14 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.6 1996/08/16 14:26:07 peter +## Fixed multi-output bug. +## Made sure that u0 is ok. +## ## Revision 1.5 1996/08/15 16:23:02 peter ## Now uses W in place of w for consistancy with time responses. ## ## Revision 1.4 1996/08/15 11:52:42 peter ## Now creats a number of versions: @@ -49,11 +53,13 @@ echo Using default parameter $PARAMS else PARAMS=$2; fi -$MATRIX << EOF > dm2fr_m.log +$MATRIX << EOF > dm2fr_m.log 2>mtt_error + + $1_numpar; [nx,ny,nu,nz,nyz] = $1_def; $PARAMS %Defaults @@ -69,12 +75,12 @@ [n,m]=size(W); if m>n W=W'; end; - [A,B,C,D,E] = $1_dm($1_numpar) - fr = dm2fr(A,B,C,D,E,W,u0) + [A,B,C,D,E] = $1_dm; + fr = dm2fr(A,B,C,D,E,W,u0); lw = log10(W); lmfr = log10(abs(fr)); pfr = angle(fr)*180/pi; % Complex frequency response @@ -103,6 +109,22 @@ nic = [nic im(:,i) re(:,i)] end; write_matrix(nic, '$1_nifr'); EOF + +# Test for errors and print if any +err_length=$(wc -c