lpnm_opt_gauss.m

jakes模型

%--------------------------------------------------------------------
% LPNM_opt_Gauss.m --------------------------------------------------
%
% Program for the computation of the discrete Doppler frequencies
% employing the Gaussian PSD by using a numerical optimization
% method.
%
% Used m-files: parameter_Gauss.m, fun_Gauss.m,
% grad_Gauss.m, acf_mue.m
%--------------------------------------------------------------------
% [f_i_n,c_i_n]=LPNM_opt_Gauss(N,f_max,f_c,sigma_0_2,p,N_i,PLOT)
%--------------------------------------------------------------------
% Explanation of the input parameters:
%
% N: length of vector tau
% f_max: maximum Doppler frequency
% f_c: 3-dB-cutoff frequency
% sigma_0_2: average power of the real Gaussian process mu_i(t)
% p: parameter of the Lp-norm (here: p=2,4,6,...)
% N_i: number of harmonic functions
% PLOT: display of the intermediate optimization results, if PLOT==1
function [f_i_n,c_i_n]=LPNM_opt_Gauss(N,f_max,f_c,sigma_0_2,p,N_i,PLOT)

kappa_c=f_max/f_c;

F_list=[];
save F_list F_list

tau_max=N_i/(2*kappa_c*f_c);
tau=linspace(0,tau_max,N);
r_mm=sigma_0_2*exp(-(pi*f_c/sqrt(log(2))*tau).^2);

[f_i_n,c_i_n,dummy]=parameter_Gauss('es_g',N_i,sigma_0_2,f_max,f_c,'none',PLOT);

save data r_mm tau N_i c_i_n p PLOT
o=foptions;
o(1)=1;
o(1)=0;
o(2)=1e-9;
o(14)=N_i/10*200;
o(9)=0;

xo=f_i_n;

x=fminu('fun_Gauss',xo,o,'grad_Gauss');

load x

f_i_n=sort(abs(x));