rice_proc.m

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

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