blob: d3378503c5fdbeb5155991d6e257fb20c5071778 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
//Example 8.5.2
//Calculate the rain attenuation
clc
clear
Upfreq=17.80 //17.80Ghz
//Polarization=vertical
kv=0.0510
av=1.0927
elevationangle=45
hs=0.05
hr=4.00
R=63
latitude=25
Ls=((hr - hs)/(sind(elevationangle)))
printf("The slant path length is %f km",Ls)
Lg=Ls*cosd(elevationangle)
printf("The horizontal projection of slant path length is %f km \n",Lg)
Yr=kv*(R)^av
printf("The specific attenuation is %f dB/km \n",Yr)
r=1/(1+0.78*(sqrt((Lg*Yr)/Upfreq))-0.38*(1-exp(-2*Lg)))
printf("The horizontal reduction factor for Miami is %f \n",r)
eta=atand((hr - hs)/(Lg*r))
Lr=(Lg*r)/(cosd(elevationangle))
X=36-abs(latitude)
v=1/(1+sqrt(sin(elevationangle))*(31*(1-exp(-(elevationangle/(1+X)))))*(sqrt(Lr*Yr)/Upfreq^2)-0.45)
printf("Vertical adjustment factor for Miami is %f \n",v)
Le=Lr*v
printf("Effective path length for Miami %f \n",Le)
A=(Yr*Le)
printf("Rain attenuation for Miami uplink of the average year is %f dB ",A)
|
