Index: mttroot/mtt/bin/trans/dm2fr_m ================================================================== --- mttroot/mtt/bin/trans/dm2fr_m +++ mttroot/mtt/bin/trans/dm2fr_m @@ -13,10 +13,13 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.12 1998/09/02 12:02:14 peterg +## Added mttinput paramater +## ## Revision 1.11 1998/08/11 13:32:42 peterg ## Lowercase mttLAST etc ## ## Revision 1.10 1998/07/27 20:27:25 peterg ## Now get parameters for simpar file. @@ -94,17 +97,15 @@ if exist('mttwsteps')==0 mttwsteps = 100; end; W = logspace(mttwmin,mttwmax,mttwsteps)'; - if exist('mttinput')==0 mttinput = 1; end; mttu0 = zeros(nu,1); mttu0(mttinput) = 1; - [n,m]=size(W); if m>n W=W'; end; @@ -112,11 +113,11 @@ [n,m]=size(mttu0); if m>n mttu0=mttu0'; end; - [A,B,C,D,E] = $1_dm + [A,B,C,D,E] = $1_dm; fr = dm2fr(A,B,C,D,E,W,mttu0); lw = log10(W); lmfr = log10(abs(fr)); pfr = angle(fr)*180/pi; @@ -132,20 +133,20 @@ % Nyquist style re = real(fr); im = imag(fr); nyq = [re(:,1) im(:,1)]; for i = 2:ny - nyq = [nyq re(:,i) im(:,i)] + nyq = [nyq re(:,i) im(:,i)]; end; write_matrix(nyq, '$1_nyfr'); % Nichols style re = lmfr; im = pfr; nic = [im(:,1) re(:,1)]; for i = 2:ny - nic = [nic im(:,i) re(:,i)] + nic = [nic im(:,i) re(:,i)]; end; write_matrix(nic, '$1_nifr'); EOF