mu_i_t.m

jakes模型

%--------------------------------------------------------------------
% Mu_i_t.m ----------------------------------------------------------
%
% Program for the simulation of real deterministic Gaussian processes
% mu_i(t) [see Fig. 4.2(b)].
%--------------------------------------------------------------------
% mu_i_t=Mu_i_t(c,f,th,T_s,T_sim,PLOT)
%--------------------------------------------------------------------
% Explanation of the input parameters:
%
% f: discrete Doppler frequencies
% c: Doppler coefficients
% th: Doppler phases
% T_s: sampling interval
% T_sim: duration of the simulation
% PLOT: plot of the deterministic Gaussian process mu_i(t),
% if PLOT==1

function mu_i_t=Mu_i_t(c,f,th,T_s,T_sim,PLOT)

if nargin==5,
   PLOT=0;
end

N=ceil(T_sim/T_s);
t=(0:N-1)*T_s;
mu_i_t=0;
for k=1:length(f),
    mu_i_t=mu_i_t+c(k)*cos(2*pi*f(k)*t+th(k));
end

if PLOT==1,
   plot(t,mu_i_t)
   xlabel('t (s)')
   ylabel('mu_i(t)')
end