Scribd
Upload
58 views1 page

Rayleigh Fading Simulation

This document summarizes a Matlab simulation of a wireless signal experiencing Doppler shift and multipath fading. The simulation generates 5000 propagation paths with varying azimuthal angles and Doppler shifts. It calculates the complex envelope of the received signal and analyzes the Rayleigh fading statistics, including the mean power, histogram of amplitudes, and complementary cumulative distribution function of signal attenuation. Plots are generated of the Rayleigh amplitude distribution, time-varying received signal intensity, and exceedance probability versus attenuation.

Uploaded by

tqminh25
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
58 views1 page

Rayleigh Fading Simulation

This document summarizes a Matlab simulation of a wireless signal experiencing Doppler shift and multipath fading. The simulation generates 5000 propagation paths with varying azimuthal angles and Doppler shifts. It calculates the complex envelope of the received signal and analyzes the Rayleigh fading statistics, including the mean power, histogram of amplitudes, and complementary cumulative distribution function of signal attenuation. Plots are generated of the Rayleigh amplitude distribution, time-varying received signal intensity, and exceedance probability versus attenuation.

Uploaded by

tqminh25
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 1

Matlab simulation: the program [3]

clear;
clc;
N=5000;
%Number of paths
t=0.0001:10/N:10;
%Time range
fc=900*10.^6;
%Carrier frequency
wc=2*pi*fc;
v=30/3600;
%Receiver speed[km/h]
c=300*10^3;
%Light speed
wm=wc*(v/c);
%Maximum shift
fm=wm/(2*pi);
%Doppler shift
for i=1:N
A(i)=(2*pi/N)*i;
%Azimuthal angles
wn(i)=wm*cos(A(i));
O(i)=(pi*i)/(N+1);
xc(i)=2*cos(wn(i)*t(i)).*cos(O(i))+cos(wm*t(i));
xs(i)=2*cos(wn(i)*t(i)).*sin(O(i));
T(i)=(1/(2*N+1)^0.5).*(xc(i)+j*xs(i));%Complex envelope
end
M=mean(abs(T));
%Mean
MdB=20*log10(M);
TdB=floor(20*log10(abs(T)));
%Field [dB]
z1=hist(abs(T));
z=hist(TdB,9);
n=0;
for k=1:9
n=n+z(k);
end
for j=1:9
P(j)=z(j)/n;
end
f(1)=P(1);
for x=2:9
f(x)=f(x-1)+P(x);
F(10-x)=f(x);
end
plot(z1)
%Distribution chart
title('Rayleighs distribution');
pause
semilogy(t,abs(T)/max(abs(T)),'r') %Fading graphic
title('Received field');
ylabel('Received field intensity');
xlabel('time');
grid on
pause
loglog(100*F,'m');%Indicates the probability of exceed a given
attenuation
title('Probability of exceed a given attenuation (%)');

Bosch, F 220976 Rev.:0002 Lang=Eng -2001

You might also like