12 files changed, 122 insertions, 0 deletions
diff --git a/3788/CH8/EX8.2.1/Ex8_2_1.PNG b/3788/CH8/EX8.2.1/Ex8_2_1.PNG
new file mode 100644
index 000000000..e3006b892
--- /dev/null
+++ b/3788/CH8/EX8.2.1/Ex8_2_1.PNG
diff --git a/3788/CH8/EX8.2.1/Ex8_2_1.sce b/3788/CH8/EX8.2.1/Ex8_2_1.sce
new file mode 100644
index 000000000..68f1b8d6a
--- /dev/null
+++ b/3788/CH8/EX8.2.1/Ex8_2_1.sce
@@ -0,0 +1,12 @@
+//Example 8.2.1
+//find the physical pathlength and path attenuation
+clc
+clear
+elevationangle=35
+height=3
+specificattenation=2
+L=height/sind(35)
+printf("The physical pathlength is %f km \n",L)
+//error for L part
+A=specificattenation*5.23
+printf("The path attenuation is %f dB " ,A)
diff --git a/3788/CH8/EX8.5.1/Ex8_5_1.PNG b/3788/CH8/EX8.5.1/Ex8_5_1.PNG
new file mode 100644
index 000000000..2b5d3361b
--- /dev/null
+++ b/3788/CH8/EX8.5.1/Ex8_5_1.PNG
diff --git a/3788/CH8/EX8.5.1/Ex8_5_1.sce b/3788/CH8/EX8.5.1/Ex8_5_1.sce
new file mode 100644
index 000000000..5c9144ded
--- /dev/null
+++ b/3788/CH8/EX8.5.1/Ex8_5_1.sce
@@ -0,0 +1,11 @@
+//Example 8.5.1
+//Calculate the specific attenuation
+clc
+clear
+rainfallrate=40
+f=10            //10 Ghz
+kv=0.00887
+av=1.264
+yr=kv*(rainfallrate)^av
+printf("The specific attenuation at 10 Ghz is %f dB/km",yr)
+
diff --git a/3788/CH8/EX8.5.2/Ex8_5_2.PNG b/3788/CH8/EX8.5.2/Ex8_5_2.PNG
new file mode 100644
index 000000000..49d389bb0
--- /dev/null
+++ b/3788/CH8/EX8.5.2/Ex8_5_2.PNG
diff --git a/3788/CH8/EX8.5.2/Ex8_5_2.sce b/3788/CH8/EX8.5.2/Ex8_5_2.sce
new file mode 100644
index 000000000..d3378503c
--- /dev/null
+++ b/3788/CH8/EX8.5.2/Ex8_5_2.sce
@@ -0,0 +1,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)
diff --git a/3788/CH8/EX8.5.3/Ex8_5_3.PNG b/3788/CH8/EX8.5.3/Ex8_5_3.PNG
new file mode 100644
index 000000000..63c1eb0c4
--- /dev/null
+++ b/3788/CH8/EX8.5.3/Ex8_5_3.PNG
diff --git a/3788/CH8/EX8.5.3/Ex8_5_3.sce b/3788/CH8/EX8.5.3/Ex8_5_3.sce
new file mode 100644
index 000000000..9698312cf
--- /dev/null
+++ b/3788/CH8/EX8.5.3/Ex8_5_3.sce
@@ -0,0 +1,9 @@
+//Example 8.5.3
+//Rain attenuation at 10 degrees
+clc
+clear
+A45=4                 //Rain atttenuation of 4 dB at elevation angle of 45
+angle=45
+newangle=10
+A=(cscd(10)/cscd(45))*A45
+printf("Rain attenuation of %f dB at elevation angle of 10 degrees",A)
diff --git a/3788/CH8/EX8.5.4/Ex8_5_4.PNG b/3788/CH8/EX8.5.4/Ex8_5_4.PNG
new file mode 100644
index 000000000..ebcab0210
--- /dev/null
+++ b/3788/CH8/EX8.5.4/Ex8_5_4.PNG
diff --git a/3788/CH8/EX8.5.4/Ex8_5_4.sce b/3788/CH8/EX8.5.4/Ex8_5_4.sce
new file mode 100644
index 000000000..c704f8aff
--- /dev/null
+++ b/3788/CH8/EX8.5.4/Ex8_5_4.sce
@@ -0,0 +1,9 @@
+//Example 8.5.4
+//RAin attenuation at 11.4GHZ
+clc
+clear
+carrierfreq=10.7
+newcarrierfreq=11.4
+A=6
+Anew=((newcarrierfreq)^2/(carrierfreq)^2)*A
+printf("Rain attenuation at Carrier frequency \n of 11.4 Ghz is %f dB \n",Anew)
diff --git a/3788/CH8/EX8.6.2/Ex8_6_2.PNG b/3788/CH8/EX8.6.2/Ex8_6_2.PNG
new file mode 100644
index 000000000..45918cb66
--- /dev/null
+++ b/3788/CH8/EX8.6.2/Ex8_6_2.PNG
diff --git a/3788/CH8/EX8.6.2/Ex8_6_2.sce b/3788/CH8/EX8.6.2/Ex8_6_2.sce
new file mode 100644
index 000000000..86f8d1e64
--- /dev/null
+++ b/3788/CH8/EX8.6.2/Ex8_6_2.sce
@@ -0,0 +1,45 @@
+//Example 8.6.2
+//Calculate the XPD 
+clc
+clear
+elevationangle=30
+f=12
+attenuation=7
+t1=20
+t2=0
+sigma=10
+p=0.01
+Cf=30*log10(f)
+printf("Cf = %f \n",Cf)
+
+Vf=12.8*f^0.19
+printf("Vf = %f \n",Vf)
+
+Ca=Vf*log10(attenuation)
+printf("Ca = %f \n",Ca)
+
+Ct1=(-10)*log10(1-0.484*(1+cosd(4*t1)))
+Ct2=(-10)*log10(1-0.484*(1+cosd(4*t2)))
+printf("Ct1 for tilt angle of 20 degrees = %f \n",Ct1)
+printf("Ct2 for tilt angle of 0 degree = %f \n",Ct2)
+
+Ctheta=-40*log10(cosd(elevationangle))
+printf("Ctheta = %f \n",Ctheta)
+
+Csigma=0.0052*sigma^2
+printf("Csigma = %f \n",Csigma)
+
+XPD1=Cf-Ca+Ct1+Ctheta+Csigma
+XPD2=Cf-Ca+Ct2+Ctheta+Csigma
+printf("XPD1 for tilt angle of 20 degrees = %f dB \n",XPD1)
+printf("XPD2 for tilt angle of 0 degree = %f dB \n",XPD2)
+
+Cice1=XPD1*(0.3+0.1*log10(p))/2
+Cice2=XPD2*(0.3+0.1*log10(p))/2
+printf("Cice1 for tilt angle of 20 degrees = %f \n",Cice1)
+printf("Cice2 for tilt angle of 0 degree = %f \n",Cice2)
+
+XPDp1=XPD1-Cice1
+XPDp2=XPD2-Cice2
+printf("XPDp1 for tilt angle of 20 degrees = %f dB \n",XPDp1)
+printf("XPDp2 for tilt angle of 0 degree = %f dB \n",XPDp2)