rice_proc.m

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

%-------------------------------------------------------------------- 
% Rice_proc.m ------------------------------------------------------- 
% 
% Program for the simulation of deterministic Rice processes xi(t)  
% (see Fig. 4.3). 
% 
% Used m-file: Mu_i_t.m 
%-------------------------------------------------------------------- 
% xi_t=Rice_proc(f1,c1,th1,f2,c2,th2,rho,f_rho,theta_rho,T_s,T_sim,PLOT) 
%-------------------------------------------------------------------- 
% Explanation of the input parameters: 
% 
% f1, c1, th1: discrete Doppler frequencies, Doppler coefficients, 
%              and Doppler phases of mu_1(t) 
% f2, c2, th2: discrete Doppler frequencies, Doppler coefficients,  
%              and Doppler phases of mu_2(t) 
% rho: amplitude of the LOS component m(t) 
% f_rho: Doppler frequency of the LOS component m(t) 
% theta_rho: phase of the LOS component m(t) 
% T_s: sampling interval 
% T_sim: duration of the simulation 
% PLOT: plot of the deterministic Rice process xi(t), if PLOT==1  
 
function xi_t=Rice_proc(f1,c1,th1,f2,c2,th2,rho,f_rho,theta_rho,T_s,T_sim,PLOT) 
 
if nargin==10, 
   PLOT=0; 
end 
 
N=ceil(T_sim/T_s); 
t=(0:N-1)*T_s; 
arg=2*pi*f_rho*t+theta_rho; 
 
xi_t=abs(Mu_i_t(c1,f1,th1,T_s,T_sim)+rho*cos(arg)+j*(Mu_i_t(c2,f2,th2,T_s,T_sim)+rho*sin(arg)) ); 
 
if PLOT==1, 
   plot(t,20*log10(xi_t)) 
   xlabel('t (s)') 
   ylabel('20 log xi(t)') 
end