jakes.m

来自「采用matlab的m语言编写」· M 代码 · 共 36 行

%Initial  paraments
Num_tath=300; %Numberofpaths
%Time range
t=0.0000001:0.0000001:0.05;
Wd=140*2*pi %Doppler shift
xc=ones(1,500000);
xs=ones(1,500000);
%Rayleigh distribution
for i=1:500000
      % xc:Re      xs:IM
      xc(i)=cos(Wd*t(i));
      xs(i)=cos(Wd*t(i));
      for k=1:75
            xc(i)=xc(i)+2*cos(p i*k/75)*.
cos((Wd*(cos(2*pi*k/75)))*t(i));
      xs( i)= xs( i)+ 2* sin (pi *k/75)*.
cos((Wd*(cos(2*pi*k/75)))*t(i));
      end
end
Rayleigh=xc+j*xs;
Amp=abs(Rayleigh);
theta=angle(Rayleigh)/pi;

i=1:500000;
title('Rayleigh Distribution');
subplot(2,1,1);
xlabel('time');
ylabel('Amp');

semilogy(i,Amp);
subplot(2,1,2 );
plot(i,theta),
xlabel('time');
ylabel('Frequence Shift):