ray_doppler.m

来自「多普勒频移下的Rayleigh衰落,更符合实际信道的带有多普勒频移的Raylei」· M 代码 · 共 16 行

function r=ray_doppler(fm,M,dt,N)
T=N*dt-dt;
t=0:dt:T;
c=sqrt(2/M);w=2*pi*fm;x=zeros(1,N);y=zeros(1,N);
for n=1:M 
alpha=(2*pi*n-pi+(2*pi*rand-pi))/(4*M);
ph1=2*pi*rand-pi;
ph2=2*pi*rand-pi;
x=x+c*cos(w*t*cos(alpha)+ph1);
y=y+c*cos(w*t*sin(alpha)+ph2);
end


r=sqrt(x.^2+ y.^2)/sqrt(2);
r_RdB=20*log10(r);
plot(t,r_RdB)