lpnm_opt_gauss.m

《移动衰落信道》Mobiel_Fading_Channels一书后面的相关仿真程序。

%-------------------------------------------------------------------- 
% 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));