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