raymodel.m

clear;
clc;
ht = input('Enter transmitter antenna height in meters\n');
hr = input('\nEnter reciever antenna height in meters\n');
d = input('\nEnter T-R separation distance in meters\n');
f = input('\nEnter carrier frequency in MHz\n');

%Consider following reference values
Eo = 10^-3;
do = 10^+3;
c = 3*10^8;
n = 120*pi;
Ae = 0.016;

%Received Electric field
l = c/(f*10^6);
pd = (2*ht*hr/d);
fd = (2*pi*pd/l);

Ed = (2*Eo*do/d)*sin(fd/2);
sprintf('\nThe received Electric field intensity is %e V/m',Ed)

%Received Power
Pr = (Ed^2*Ae)/n;
PrdB = 10*log10(Pr);
sprintf('\nThe received Power is %e Watts or %f dB',Pr,PrdB)