fun_gauss.m

来自「jakes模型」· M 代码 · 共 37 行

%--------------------------------------------------------------------
% fun_Gauss.m -------------------------------------------------------
%
% Computation of the error function according to Eq.(5.61) for the
% optimization of the discrete Doppler frequencies (Gaussian PSD).
%
% Used m-file: acf_mue.m
%--------------------------------------------------------------------
% F=fun_Gauss(x)
%--------------------------------------------------------------------
% Explanation of the input parameters:
%
% x: parameter vector to be optimized
function F=fun_Gauss(x)

load data

f_i_n=x;

r=acf_mue(f_i_n,c_i_n,tau);
F=norm(abs(r_mm-r),p);
if PLOT==1,
   subplot(1,2,1)
   stem(f_i_n,c_i_n)
   xlabel('f_i_n')
   ylabel('c_i_n')
   
   title(['N_i = ',num2str(N_i)])
   subplot(1,2,2)
   plot(tau,r_mm,tau,r)
   xlabel('tau (s)')
   ylabel('ACF')
   title(['Error-norm=',num2str(F)])
   pause(0)
end

save x x