rayleighgen.m

这是一个wcdma下行功率控制仿真的源代码

%Jakes model
function [JakesFader] = RayleighGen(startTime,dopller,totaltime,sampletime)
N0 = 8;
N = 4 * N0 + 2;
fm = (dopller * 1000.0 / 3600) * 2 * 10^9 / 300000000;
Wm = 2*pi*fm;
colnum = totaltime / sampletime;
%maxphase is phase according to max(wm)
maxphase = pi / 4;       
j = 1:colnum;
%Xc is real of fading ,Xs is imag of fading
Xc = 0;
Xs = 0;
n = 1:N0;
t = startTime + j * sampletime;
Beta = Wm * cos(2 * pi * n' / N) * t;
Alpha = pi * n'/ N0 * ones(1,length(t));
Xcc = (cos(Alpha+Beta)+cos(Alpha-Beta))/2;
Xss = (sin(Alpha+Beta)+sin(Alpha-Beta))/2;
Xc = sum(Xcc,1);
Xs = sum(Xss,1);
Xc = (Xc*2 + sqrt(2.0) * cos(maxphase) * cos(Wm * t))/sqrt(N0 + 1/2);
Xs = (Xs*2 + sqrt(2.0) * sin(maxphase) * cos(Wm * t))/sqrt(N0 + 1/2);
JakesFader = complex(Xc,Xs);