Jakes模型matlab仿真验证过

% Simulation Of Jakes Model 
clear all; 
f_max = 10; 
M = 9;                                 % # of low frequency oscillators - 1 
N = 4*M+2; 
Ts=1e-03; 
sq = 2/sqrt（N）; 
sigma = 1/sqrt（2）; 
theta = 0;                              % Fixed Phase 
count = 0; 
t0=0.24; %（0.24 0.45）; 
    for t = 0:Ts:3            % Varying time 
               
        count = count + 1 
        g（count） = 0; 
        for n = 1 : M+1 
            if n == M 
               c_q（countn） = 2*sigma*sin（pi*n/M）;  % Gain associated with quadrature component 
               c_i（countn） = 2*sigma*cos（pi*n/M）;  % Gain associated with inphase component 
               f_i（countn） = f_max*cos（2*pi*n/N）;  % Discrete doppler frequencies of inphase component 
               f_q（countn） = f_max*cos（2*pi*n/N）;  % Discrete doppler frequencies of quadrature component 
           else 
               c_i（countn） = sqrt（2）*cos（pi/4）; 
               c_q（countn） = sqrt（2）*sin（pi/4）; 
               f_i（countn） = f_max; 
               f_q（countn） = f_max; 
           end;  % end if 
                
               g_i（countn） =  c_i（countn）*cos（2*pi*f_i（countn）*（t-t0） + theta）; % Inphase component for one oscillator 
              
               g_q（countn） =  c_q（countn）*cos（2*pi*f_q（countn）*（t-t0） + theta）; % Quadrature component for one oscillator 
                
       end;   %end n 
       tp（count） =  sq*sum（g_i（count1:M+1））;% Total Inphase component  
       tp1（count） = sq*sum（g_q（count1:M+1））;% Total quadrature component  
       
   end;  % end count no n again 
   envelope=sqrt（tp.^2+tp1.^2）; 
   rmsenv=sqrt（sum（envelope.^2）/count）; 
    
   [auto_ilag_i] = xcorr（tp‘coeff‘） ;   % Auto-correlation associated with inphase component