gen_rice_proc.m

一些非常有用的MATLAB源码,对于提高编程水平和学习通信理论知识还有帮助.

%--------------------------------------------------------------------
% gen_Rice_proc.m ---------------------------------------------------
%
% Program for the simulation of deterministic generalized Rice 
% processes (see Fig. 6.29).
%
% Used m-files: parameter_Jakes.m, Mu_i_t.m
%--------------------------------------------------------------------
% xi_t=gen_Rice_proc(N_1,N_2,sigma_1_2,sigma_2_2,kappa_0,...
%                    theta_0,rho,theta_rho,f_max,...
%                    T_s,T_sim,PLOT)
%--------------------------------------------------------------------
% Explanation of the input parameters:
%
% N_1, N_2: number of harmonic functions of the real deterministic  
%           Gaussian processes nu_1(t) and nu_2(t), respectively
% sigma_1_2: average power of the real deterministic Gaussian 
%            process nu_1(t)
% sigma_2_2: average power of the real deterministic Gaussian 
%            process nu_2(t)
% kappa_0: frequency ratio f_min/f_max (0<=kappa_0<=1)
% theta_0: phase shift between mu_1_n(t) and mu_2_n(t)
% rho: amplitude of the LOS component m(t)
% theta_rho: phase of the LOS component m(t)
% f_max: maximum Doppler frequency
% T_s: sampling interval
% T_sim: duration of the simulation
% PLOT: plot of the deterministic generalized Suzuki process xi(t), 
%       if PLOT==1 

function xi_t=gen_Rice_proc(N_1,N_2,sigma_1_2,sigma_2_2,kappa_0,...
                            theta_0,rho,theta_rho,f_max,T_s,...
                            T_sim,PLOT)

if nargin==11,
   PLOT=0;
end

[f1,c1,th1]=parameter_Jakes('es_j',N_1,sigma_1_2,f_max,'rand',0);
c1=c1/sqrt(2);

N_2_s=ceil(N_2/(2/pi*asin(kappa_0)));
[f2,c2,th2]=parameter_Jakes('es_j',N_2_s,sigma_2_2,f_max,'rand',0);
f2 =f2(1:N_2);
c2 =c2(1:N_2)/sqrt(2);
th2=th2(1:N_2);

N=ceil(T_sim/T_s);
t=(0:N-1)*T_s;

xi_t=abs(Mu_i_t(c1,f1,th1,T_s,T_sim)+...
         Mu_i_t(c2,f2,th2,T_s,T_sim)+rho*cos(theta_rho)+...
         j*(Mu_i_t(c1,f1,th1-theta_0,T_s,T_sim)+...
         Mu_i_t(c2,f2,th2+theta_0,T_s,T_sim)+...
         rho*sin(theta_rho)));

if PLOT==1,
   plot(t,20*log10(xi_t),'b-')
   xlabel('t (s)')
   ylabel('20 log xi(t)')
end