summaryrefslogtreecommitdiff
path: root/3739/CH3
diff options
context:
space:
mode:
Diffstat (limited to '3739/CH3')
-rw-r--r--3739/CH3/EX3.1/EX3_1.sce30
-rw-r--r--3739/CH3/EX3.10/EX3_10.sce23
-rw-r--r--3739/CH3/EX3.11/EX3_11.sce23
-rw-r--r--3739/CH3/EX3.12/EX3_12.pdfbin0 -> 16025 bytes
-rw-r--r--3739/CH3/EX3.12/EX3_12.sce94
-rw-r--r--3739/CH3/EX3.13/EX3_13.sce20
-rw-r--r--3739/CH3/EX3.14/EX3_14.sce17
-rw-r--r--3739/CH3/EX3.15/EX3_15.sce15
-rw-r--r--3739/CH3/EX3.2/EX3_2.sce19
-rw-r--r--3739/CH3/EX3.3/EX3_3.sce23
-rw-r--r--3739/CH3/EX3.4/EX3_4.sce19
-rw-r--r--3739/CH3/EX3.5/EX3_5.sce13
-rw-r--r--3739/CH3/EX3.6/EX3_6.sce18
-rw-r--r--3739/CH3/EX3.7/EX3_7.sce24
-rw-r--r--3739/CH3/EX3.9/EX3_9.sce27
15 files changed, 365 insertions, 0 deletions
diff --git a/3739/CH3/EX3.1/EX3_1.sce b/3739/CH3/EX3.1/EX3_1.sce
new file mode 100644
index 000000000..c9c7910c8
--- /dev/null
+++ b/3739/CH3/EX3.1/EX3_1.sce
@@ -0,0 +1,30 @@
+//Chapter 3, Example 3.1, page 61
+clc
+//Initialisation
+h=2 //height in Km
+h1=5 //height in Km
+
+
+//Calculation
+t2=290-(6.5*h) //Proposed formula for height h=2Km
+p2=950-117*h
+e2=8-3*h
+t21=294.98-5.22*h-0.007*h**2
+p21=1012.82-111.56*h+3.86*h**2
+p=14.35*2.72**(-0.42*h-0.02*h*h+0.001*h**3)
+
+t5=290-6.5*h1 //proposed formula for height h=5Km
+p5=950-117*h1
+e5=8-3*h1
+t51=294.98-5.22*h1-0.007*h1**2
+p51=1012.82-111.56*h1+3.86*h1**2
+p1=14.35*2.72**(-0.42*h1-0.02*h1**2+0.001*h1**3)
+
+
+//Results
+printf("T(2) = %.1f K",t21)
+printf("\nP(2) = %.2f hpa",p21)
+printf("\np(2) = %.2f hpa",p)
+printf("\nT(5) = %.1f K",t51)
+printf("\nP(5) = %.2f hpa",p51)
+printf("\np(5) = %.2f hpa",p1)
diff --git a/3739/CH3/EX3.10/EX3_10.sce b/3739/CH3/EX3.10/EX3_10.sce
new file mode 100644
index 000000000..0b458aa5b
--- /dev/null
+++ b/3739/CH3/EX3.10/EX3_10.sce
@@ -0,0 +1,23 @@
+//Chapter 3, Example 3.10, page 92
+clc
+//Initialisation
+f=5*10**9 //frequency in Hz
+c=3*10**8 //speed of light
+h1=6 //in metre
+h2=2 //in metre
+pi=3.14
+
+//Calculation
+h=c*f**-1 //wavelength
+w=atan(h1*2250**-1) //grazing angle in radian
+w1=w*180*pi**-1 //grazing angle in degree
+a=((2*pi*h1*h2)*(h*300)**-1)*3.14*180**-1
+e=sin(a)
+F=e*2*180*pi**-1 //PGF value (wrong value calculated in textbook)
+LR=20*log10(F) //Decrease in received signal level
+
+
+//Results
+printf("(1) Grazing angle = %.2f degree",w1)
+printf("\n(2) PGF value = %f",F) //The answer provided in the textbook is wrong
+printf("\n(3) Decrease in received signal level = %.2f dB",LR) //The answer provided in the textbook is wrong
diff --git a/3739/CH3/EX3.11/EX3_11.sce b/3739/CH3/EX3.11/EX3_11.sce
new file mode 100644
index 000000000..32725d46d
--- /dev/null
+++ b/3739/CH3/EX3.11/EX3_11.sce
@@ -0,0 +1,23 @@
+//Chapter 3, Example 3.11, page 98
+clc
+
+//Initialisation
+h=12.5*10**-2 //in meter
+d1=10*10**3 //in meter
+d2=15*10**3 //in meter
+d3=12.5*10**3 //in meter
+d4=12.5*10**3 //in meter
+h=1.25 //in Kilometer
+
+//Calculation
+r1=(((d1*d2)/(d1+d2))*h)**0.5 //radius of first and fourth Fresnel zones
+r4=r1*(4)**0.5
+R1=(((d3*d4)/(d3+d4))*h)**0.5 //radius of first and fourth ellipse zones
+R4=R1*(4)**0.5
+
+//Results
+printf("Radius of first fresnel zones, r1 = %.2f m",r1)
+printf("\nRadius of Second fresnel zones, r4= %.2f m",r4)
+printf("\nh = %.2f x 10^-4 Km",h)
+printf("\nRadius of first ellipse, R1 = %.2f m",R1)
+printf("\nRadius of second ellipse, R4 = %.1f m",R4)
diff --git a/3739/CH3/EX3.12/EX3_12.pdf b/3739/CH3/EX3.12/EX3_12.pdf
new file mode 100644
index 000000000..c2f0ba87d
--- /dev/null
+++ b/3739/CH3/EX3.12/EX3_12.pdf
Binary files differ
diff --git a/3739/CH3/EX3.12/EX3_12.sce b/3739/CH3/EX3.12/EX3_12.sce
new file mode 100644
index 000000000..32464704b
--- /dev/null
+++ b/3739/CH3/EX3.12/EX3_12.sce
@@ -0,0 +1,94 @@
+//Chapter 3, Example 3.12, page 105
+clc
+
+
+//Initialisation
+L=13200 //L parameter in m
+H=10240 //H parameter
+Re=6370000 //actual redius of earth
+ht=30 //height in m
+hr=20 // in m
+re1=8453000 // in metre
+h1=30000 // in metre
+h2=20000 //in metre
+dt1=22.5
+f=10*10**9 //frequency in Hz
+c=3*10**8 //speed of light
+d=30000 //distance in m
+pt=30 //transmitter antenna power
+gt=40 //transmitter antenna gain
+gr=40 //receiver antenna gain
+pi=3.14
+F3=-3
+H=-34
+D=0.75
+
+//Calculation
+dt=sqrt(2*re1*ht)
+X=3*dt*L**-1
+Z1=h1*H**-1
+Z2=h2*H**-1
+vx=10**-3.5 //from fig 3.26
+z1=10**0.95 //from fig 3.27
+z2=10**0.65 //from fig 3.27
+
+//for d=3dt
+lv=20*log10(vx)
+lz1=20*log10(z1)
+lz2=20*log10(z2)
+F=(lv+lz1+lz2)*20**-1
+F1=10**(F)
+F11=20*log10(F1)
+X1=2*dt*L**-1
+d3=3
+f3=-F11
+
+vx1=10**-2.35 //from fig 3.26
+lv1=20*log10(vx1)
+
+//for d=2dt
+F4=1+D
+F5=20*log10(F4)
+d2=2
+f2=-F5
+
+
+//for d=1.1dt
+F6=sqrt(1+D**2)
+F7=20*log10(F6)
+d11=1.1
+f11=-F7
+
+//for d=dt
+d1=1
+f1=0.2
+
+//for plotting graph in terms of points
+
+
+
+for N=0:1:5
+ a=plot(1,0.2,'-o')
+ a1=plot(1.1,-1.9,'-o')
+ a2=plot(2,-4.8,'-o')
+ a3=plot(3,-38,'-o')
+end
+
+title('Path gain F','fontsize',5);
+xlabel("d/dt", "fontsize", 3);
+ylabel("20log(F)(dB)", "fontsize", 3, "color", "blue");
+xstring(1,2,"d/dt",0,0);
+xstring(1.2,0.7,"1.1d/dt",0,0);
+xstring(2,-0.7,"2d/dt",0,0);
+xstring(2.86,-35,"3d/dt",0,0);
+
+
+
+h=c*f**-1 //wavelength
+Pr=pt+gt+gr+H+F3-10*log10(4*pi*d**2) //Received signal power
+
+
+//Results
+printf("(1) Effective receiver path gain F = %.4f",F11)
+printf("\n(2) Path gain F plot is shown")
+printf("\n(3) Received signal power Pr = %.1f dBm",Pr)
diff --git a/3739/CH3/EX3.13/EX3_13.sce b/3739/CH3/EX3.13/EX3_13.sce
new file mode 100644
index 000000000..0f675b1a2
--- /dev/null
+++ b/3739/CH3/EX3.13/EX3_13.sce
@@ -0,0 +1,20 @@
+//Chapter 3, Example 3.13, page 109
+clc
+
+//Initialisation
+eirp=800 //in KW
+d=24 //in Km
+a=0.03 //in radian
+d1=22 //in Km
+d2=2 //in Km
+h=0.4*10**-3 //wavelength in m
+Er=45 //in microvolt
+
+//Calculation
+E=104.8+10*log10(eirp)-20*log10(d) //field intensity
+V=a*sqrt((2*d2*d1)/((d1+d2)*h)) //knife edge obstacle attenuation
+Lke=23 //from table 3.4
+er=10**(Er*20**-1)
+
+//Results
+printf("(1) Electric field intensity = %.3f microV/m",er)
diff --git a/3739/CH3/EX3.14/EX3_14.sce b/3739/CH3/EX3.14/EX3_14.sce
new file mode 100644
index 000000000..9ac627270
--- /dev/null
+++ b/3739/CH3/EX3.14/EX3_14.sce
@@ -0,0 +1,17 @@
+//Chapter 3, Example 3.14, page 115
+clc
+
+//Initialisation
+f1=430 //upper frequency band
+f2=410 //lower frequency band
+d=80 //distance in meter
+
+//Calculation
+Yv=0.1 //Specific attenuation obtained from graph fig 3.34
+Lv=Yv*d //loss of forest trees
+Am=((f1+f2)/2)**0.5 //maximum value for trees excess loss.
+
+//Results
+printf("Specific attenuation index, Yv = %.1f dB/m",Yv)
+printf("\nLoss of forest trees, Lv = %.1f dB",Lv)
+printf("\nMaximum value for trees excess loss = %.1f dB",Am)
diff --git a/3739/CH3/EX3.15/EX3_15.sce b/3739/CH3/EX3.15/EX3_15.sce
new file mode 100644
index 000000000..f567fee49
--- /dev/null
+++ b/3739/CH3/EX3.15/EX3_15.sce
@@ -0,0 +1,15 @@
+//Chapter 3, Example 3.15, page 118
+clc
+
+//Initialisation
+d=40 //length in meter
+Am=2 //area in square meter
+f=10*10**9 //frequency in hertz
+
+//Calculation
+As=40 //using graph fig 3.36, As can be obtained
+As1=30 //using graph fig 3.37, As can be obtained
+
+//Results
+printf("Loss in the summer for trees with leaves, As = %d dB",As)
+printf("\nLoss in winter for trees without leaves, As = %d dB",As1)
diff --git a/3739/CH3/EX3.2/EX3_2.sce b/3739/CH3/EX3.2/EX3_2.sce
new file mode 100644
index 000000000..675268f62
--- /dev/null
+++ b/3739/CH3/EX3.2/EX3_2.sce
@@ -0,0 +1,19 @@
+//Chapter 3, Example 3.2, page 63
+clc
+
+//Initialisation
+h=2 //Height in Km
+T=277 //Tempreture in Kelvin
+p=716
+e=2
+
+
+//Calculation
+er=1+(151.1/T)*(p+(4810*h/T))*10**-6
+n=er**(0.515) //refractive index of the air
+N=(n-1)*10**6 //refractivity number
+
+
+//Results
+printf("n = %.5f",n)
+printf("\nN = %d",N)
diff --git a/3739/CH3/EX3.3/EX3_3.sce b/3739/CH3/EX3.3/EX3_3.sce
new file mode 100644
index 000000000..54c1b2457
--- /dev/null
+++ b/3739/CH3/EX3.3/EX3_3.sce
@@ -0,0 +1,23 @@
+//Chapter 3, Example 3.3, page 67
+clc
+//Initialisation
+er=1.001 //relative permittivity of a medium
+dn=35*10**-6 //vertical gradient of refractive index
+Re=6370 //actual redius of earth
+d=20 //transmitter and receiver distance in Km
+d1=5
+d2=15
+K1=1.3333 //standard atmosphere condition
+
+//Calculation
+R=(er**0.5)/dn
+K=R/(R-Re)
+hm=(125*d**2)/(K*Re) //Earth bulge value in the middle of the path
+h1=(500*d1*d2)/(K*Re) //h1
+h2=(500*d1*d2)/(K1*Re) //h2
+
+
+//Results
+printf("Bulge value = %.1f m",hm)
+printf("\nBulge value, h1 = %.2f m",h1)
+printf("\nBulge value, h2 = %.2f m",h2)
diff --git a/3739/CH3/EX3.4/EX3_4.sce b/3739/CH3/EX3.4/EX3_4.sce
new file mode 100644
index 000000000..42f10c6f3
--- /dev/null
+++ b/3739/CH3/EX3.4/EX3_4.sce
@@ -0,0 +1,19 @@
+//Chapter 3, Example 3.4, page 68
+clc
+//Initialisation
+K=1.33
+d1=24 //heigth in Km
+d2=15 //heigth in Km
+K1=1
+Re=6370 //actual redius of earth
+
+//Calculation
+R=4.12*(d1**0.5+d2**0.5)
+R1=K1*Re
+Rrh=(2*R1*d1)**0.5+(2*R1*d2)**0.5
+
+//Results
+printf("K=1.33")
+printf("\nRrh = %.1f km\n",R)
+printf("K=1")
+printf("\nRrh = %.1f km",Rrh) //The answer provided in the textbook is wrong
diff --git a/3739/CH3/EX3.5/EX3_5.sce b/3739/CH3/EX3.5/EX3_5.sce
new file mode 100644
index 000000000..ace5edaab
--- /dev/null
+++ b/3739/CH3/EX3.5/EX3_5.sce
@@ -0,0 +1,13 @@
+//Chapter 3, Example 3.5, page 74
+clc
+//Initialisation
+No=1 //index of refraction
+N1=1.3*10**-7
+h=20 //height
+
+//Calculation
+wc=asin(((4*No)/((4*No)+((h**2)*N1)))) //critical angle
+
+
+//Results
+printf("Critical angle = %f",wc) //answer is not written in textbook
diff --git a/3739/CH3/EX3.6/EX3_6.sce b/3739/CH3/EX3.6/EX3_6.sce
new file mode 100644
index 000000000..c2a188fc1
--- /dev/null
+++ b/3739/CH3/EX3.6/EX3_6.sce
@@ -0,0 +1,18 @@
+//Chapter 3, Example 3.6, page 76
+clc
+//Initialisation
+dn=-0.2 //air refractivity gradient
+d=20 //height
+b=0.074 //elevation angle from graph 3.10
+f=7 //frequency in Ghz from graph 3.11
+c=2*10**-6
+
+///Calculation
+t=0.156 //1000/6370
+dm=dn+t
+a=(-c*dm*d)**0.5 //elevation angle of waves
+
+//Results
+printf("Elevation angle of waves = %.2f mrad",(a*10**3))
+printf("\nElevation angle = %.3f",b) //from graph 3.10
+printf("\nMinimum frequency of coupling waves into the duct = %d Ghz",f) //from graph 3.11
diff --git a/3739/CH3/EX3.7/EX3_7.sce b/3739/CH3/EX3.7/EX3_7.sce
new file mode 100644
index 000000000..2960a5bb8
--- /dev/null
+++ b/3739/CH3/EX3.7/EX3_7.sce
@@ -0,0 +1,24 @@
+//Chapter 3, Example 3.7, page 80
+clc
+
+//Initialisation
+f=18 //frequency in GHz
+d=30 //in km
+R=25 //rainfall intensity in mm
+
+//Using Table 3.3
+av15=1.128
+av20=1.065
+av18=1.09
+kv15=0.0335
+kv20=0.0691
+kv18=0.0587
+
+//Calculation
+yr=kv18*R**av18 //rain specific attenuation
+de=(90*(90+d)**-1)*d
+A=de*yr //Maximum rain attenuation
+
+//Results
+printf("(1) Rain specific attenuation = %.2f dB/km",yr)
+printf("\n(2) Maximum rain attenuation = %.1f dB",A)
diff --git a/3739/CH3/EX3.9/EX3_9.sce b/3739/CH3/EX3.9/EX3_9.sce
new file mode 100644
index 000000000..c947d6198
--- /dev/null
+++ b/3739/CH3/EX3.9/EX3_9.sce
@@ -0,0 +1,27 @@
+//Chapter 3, Example 3.9, page 89
+clc
+//Initialisation
+rh=-1
+s=4 //sigma in S/m
+f=5*10**9 //frequency in Hz
+eo=8.85*10**-12 //permitivity of free space
+er=75 //permitivity of medium
+w1=30*3.14*180**-1 //in radians
+pi=3.14
+
+
+//Calculation
+w=2*pi*f
+x=s*(w*eo)**-1
+a=sin(w1)-sqrt((er-x)-cos(w1)**2)
+a1=sin(w1)+sqrt((er-x)-cos(w1)**2)
+rh1=a/a1
+b1=(er-x)*sin(w1)-sqrt((er-x)-cos(w1)**2)
+b2=(er-x)*sin(w1)+sqrt((er-x)-cos(w1)**2)
+rv=-b1/b2
+
+
+//Results
+printf("(2) X = %.1f",x)
+printf("\n(3) Rh = %.3f",rh1)
+printf("\n Rv = %.1f",rv)
generated by cgit (git 2.34.1) at 2024-12-22 09:42:58 +0000