LPNM_opt_Jakes.m
上传用户:hbrsgg3
上传日期:2022-08-01
资源大小:5k
文件大小:1k
源码类别:

通讯/手机编程

开发平台:

Windows_Unix

  1. %--------------------------------------------------------------------
  2. % LPNM_opt_Jakes.m --------------------------------------------------
  3. %
  4. % Program for the computation of the discrete Doppler frequencies
  5. % employing the Jakes PSD by using a numerical optimization method.
  6. %
  7. % Used m-files: parameter_Jakes.m, fun_Jakes.m,
  8. % grad_Jakes.m, acf_mue.m
  9. %--------------------------------------------------------------------
  10. % [f_i_n,c_i_n]=LPNM_opt_Jakes(N,f_max,sigma_0_2,p,N_i,PLOT)
  11. %--------------------------------------------------------------------
  12. % Explanation of the input parameters:
  13. %
  14. % N: length of vector tau
  15. % f_max: maximum Doppler frequency
  16. % sigma_0_2: average power of the real Gaussian process mu_i(t)
  17. % p: parameter of the Lp-norm (here: p=2,4,6,...)
  18. % N_i: number of harmonic functions
  19. % PLOT: display of the intermediate optimization results, if PLOT==1
  20. function [f_i_n,c_i_n]=LPNM_opt_Jakes(N,f_max,sigma_0_2,p,N_i,PLOT)
  21. tau=linspace(0,N_i/(2*f_max),N);
  22. Jo=sigma_0_2*besselj(0,2*pi*f_max*tau);
  23. %c_i_n=sqrt(sigma_0_2)*sqrt(2/N_i)*ones(N_i,1);
  24. % save data Jo tau N_i c_i_n p PLOT
  25. % Initial values:
  26. [f_i_n,dummy1,dummy2]=parameter_Jakes('es_j',N_i,...
  27. sqrt(sigma_0_2),f_max,'none',0);
  28. o=foptions;
  29. o(1)=1;
  30. o(1)=0;
  31. o(2)=1e-9;
  32. o(14)=N_i/10*200;
  33. o(9)=0;
  34. xo=f_i_n;
  35. x=fminu('fun_Jakes',xo,o,'grad_Jakes');
  36. f_i_n=x;