matlab coding for All propagation model using rayleigh
Download as DOCX, PDF0 likes287 views
This document contains Matlab code to simulate and compare the path loss of radio signals between a transmitter and receiver using different propagation models over a distance of 5 km. It generates a Rayleigh faded signal, calculates the predicted received power using four path loss models, and plots the actual received power versus distance to compare the models.
![5.All models
clc;
clear all;
close all;
f=900 %FREQUENCY in Mega Hertz
%DISTANCE BETWEEN TRANSMITTER AND
RECEIVER in Kilometers
%v=input('MU speed in m/s....5, 10, 15, 20, 25,..>');
D=5000.0;
v=45;
numpaths1 =5;
Fs = 4*f; %sampling frequency
Ts = 1/Fs; %sampling period
k=0.48/Ts;
tmax=D/v;
t=[0:k*Ts:tmax];
d=v*t
wc = 2*pi*f; %radian frequency
ray0 = zeros(1,length(t)); %received signal
for i = 1:numpaths1
wd = 2*pi*v*f*cos(unifrnd(0,2*pi))/3e8;
a = weibrnd(1,3,1,length(t));
ray0 = ray0 + a.*cos((wc+wd)*t+unifrnd(0,2*pi,1,length(t)))
end;
[ray0i ray0q] = demod(ray0,f,Fs,'qam') %demodulated signal
env_ray0 = sqrt(ray0i.^2+ray0q.^2) %envelope of received signal
power_ray0=env_ray0.^2
powerdB0=10*log10(power_ray0)
hb=30 ; %BASE STATION ANTENNA HEIGHT in meter
hr=3; %MOBILE STATION ANTENNA HEIGHT in meter
ht=30; %TRANMITTER ANTENNA
HEIGHT in meter
cm=3;
ahm= (3.2*(log10(11.75*hr))^2-4.97);
PL1=46.3+33.9*log10(f)-13.82*log10(hb)-ahm+(44.9-
6.55*log10(hb))*(log10(d))+cm; %PATH LOSS in dB
pr1=43-PL1;
display(powerdB0)
ahr=3.2*(log10(11.75*hr).^2-4.97); % GAIN FACTOR IN
URBAN AREA in dB
PL2=69.55+(26.16*(log10(f)))-(13.82*log10(ht))-(ahr)+(44.9-
6.55*(log10(ht)))*(log10(d));
pr2=43-PL2;
display(powerdB0)
lambda=3*10^8/(f*10^6); % WAVELENGTH in Meter
d0=0.1 % DISTANCE in KiloMetre](https://image.slidesharecdn.com/allmodelusingrayleigh-150419084728-conversion-gate01/85/matlab-coding-for-All-propagation-model-using-rayleigh-1-320.jpg)
matlab coding for All propagation model using rayleigh
- 1. 5.All models clc; clear all; close all; f=900 %FREQUENCY in Mega Hertz %DISTANCE BETWEEN TRANSMITTER AND RECEIVER in Kilometers %v=input('MU speed in m/s....5, 10, 15, 20, 25,..>'); D=5000.0; v=45; numpaths1 =5; Fs = 4*f; %sampling frequency Ts = 1/Fs; %sampling period k=0.48/Ts; tmax=D/v; t=[0:k*Ts:tmax]; d=v*t wc = 2*pi*f; %radian frequency ray0 = zeros(1,length(t)); %received signal for i = 1:numpaths1 wd = 2*pi*v*f*cos(unifrnd(0,2*pi))/3e8; a = weibrnd(1,3,1,length(t)); ray0 = ray0 + a.*cos((wc+wd)*t+unifrnd(0,2*pi,1,length(t))) end; [ray0i ray0q] = demod(ray0,f,Fs,'qam') %demodulated signal env_ray0 = sqrt(ray0i.^2+ray0q.^2) %envelope of received signal power_ray0=env_ray0.^2 powerdB0=10*log10(power_ray0) hb=30 ; %BASE STATION ANTENNA HEIGHT in meter hr=3; %MOBILE STATION ANTENNA HEIGHT in meter ht=30; %TRANMITTER ANTENNA HEIGHT in meter cm=3; ahm= (3.2*(log10(11.75*hr))^2-4.97); PL1=46.3+33.9*log10(f)-13.82*log10(hb)-ahm+(44.9- 6.55*log10(hb))*(log10(d))+cm; %PATH LOSS in dB pr1=43-PL1; display(powerdB0) ahr=3.2*(log10(11.75*hr).^2-4.97); % GAIN FACTOR IN URBAN AREA in dB PL2=69.55+(26.16*(log10(f)))-(13.82*log10(ht))-(ahr)+(44.9- 6.55*(log10(ht)))*(log10(d)); pr2=43-PL2; display(powerdB0) lambda=3*10^8/(f*10^6); % WAVELENGTH in Meter d0=0.1 % DISTANCE in KiloMetre
- 2. A=20*log10(4*pi*d0/lambda) a=4.6; b=0.0075; c=12.6; hb=30 % BASE STATION ANTENNA HEIGHT in Meter gamma=a-(b*(hb))+(c/hb) Xf=(6.0)*(log10(f/2000)) hr=3 %MOBILE STATION ANTENNA HEIGHT in Meter Xh=-20.80*log10(hr/2000) s=10.4 % LOG NORMALLY DISTRIBUTED FACTOR Amu=18 PL3=A+(10*(gamma)*(log10(d/d0)))+Xf+Xh+s % PATH LOSS in dB pr3=43-PL3 display(powerdB0) Gr=10*log10(hr/3) %RECEIVER ANTENNA HEIGHT GAIN FACTOR in dB Gt=20*log10(ht/200) %TRANSMITTER ANTENNA HEIGHT GAIN FACTOR in dB Garea=9 Lf=-10*log10((lambda.^2)*(((4*pi*d*1000).^2).^-1)) PL4=((Lf)+(Amu)-(Gt)-(Gr)-(Garea)) % PATH LOSS in dB pr4=43-PL4 display(powerdB0) %fadedSig = pr+powerdB0 actual1=pr1-powerdB0 actual2=pr2-powerdB0 actual3=pr3-powerdB0 actual4=pr4-powerdB0 plot(d,actual1,'--',d,actual2,'--',d,actual3,'-.',d,actual4,'-*'); legend('actual1= cost','actual2= hata','actual3= sui','actual4= okumura'); title('comparison of Signal Strength incorporating Rayleigh Fading '); xlabel('Distance in meters'); ylabel('Signal Strength in dBm');