fun_jakes.m

jakes模型

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

load data

f_i_n=x;
r=acf_mue(f_i_n,c_i_n,tau);
F=norm(abs(Jo-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,Jo,tau,r)
   xlabel('tau (s)')
   ylabel('ACF')
   title(['Error-norm=',num2str(F)])
   pause(0)
end
save x x