From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books
---
2409/CH16/EX16.1/Ex16_1.sce | 46 +++++++++++++++++++++++++++++++++++++
2409/CH16/EX16.2/Ex16_2.sce | 55 +++++++++++++++++++++++++++++++++++++++++++++
2 files changed, 101 insertions(+)
create mode 100755 2409/CH16/EX16.1/Ex16_1.sce
create mode 100755 2409/CH16/EX16.2/Ex16_2.sce
(limited to '2409/CH16')
diff --git a/2409/CH16/EX16.1/Ex16_1.sce b/2409/CH16/EX16.1/Ex16_1.sce
new file mode 100755
index 000000000..b3dc9b1ab
--- /dev/null
+++ b/2409/CH16/EX16.1/Ex16_1.sce
@@ -0,0 +1,46 @@
+
+//Varaible Declaration
+
+EIRP=55 //EIRP for satellite(dBW)
+fD=12.5 //Downlink frequency(GHz)
+Pss=-101 //Receiving at ground station direction(degrees west)
+Rb=40*10**6 //Transmission Rate(Hz)
+D=18 //Diameter of antenna(inches)
+n=0.55 //Efficiency of antenna
+Tant=70 //Antenna noise(Kelvin)
+Teq=100 //Equivalent noise temperature at LNA(Kelvin)
+R=6371 //Radius of earth(Km)
+L=2 //Transmission losses(dB)
+aGSO=42164 //Circumference of earth(km)
+k=-228.6 //Boltzmann's constant (dB)
+PE=-90 //Longitude of Earth station(degrees west)
+LE=45 //Latitude of Earth station(degrees north)
+f=14 //Frequency(GHz)
+//Calculation
+B=PE-Pss
+b=acos(cos(B*3.142/180)*cos(LE*3.142/180))
+b=b*180/3.142
+A=asin(sin(abs(B)*3.142/180)/sin(b*3.142/180))
+A=A*180/3.142
+Az=180+A //Azimuth angle of antenna(degrees)
+d=(R**2+aGSO**2-2*R*aGSO*cos(b*3.142/180))**0.5 //Range of antenna(km)
+El=acos(aGSO*sin(b*3.142/180)/d) //Elevation angle of antenna(radians)
+El=El*180/3.142 //Elevation angle of antenna(degrees)
+El=round(El)
+d=round(d)
+FSL=32.4+20*log10(d)+20*log10(f*10**3) //Free space loss(dB)
+LOSSES=FSL+L //Total Transmission Losses
+Ts=Teq+Tant //Total system noise temperature(Kelvin)
+T=10*log10(Ts) //Total system noise temperature(dBK)
+G=n*(3.192*f*(D/(12)))**2
+G=10*log10(G) //Antenna Gain(dB)
+GTR=G-T //G/T ratio(dB)
+CNR=EIRP+GTR-LOSSES-k //Carrier to noise ratio(dB)
+Rb=10*log10(Rb) //Transmission Rate(dBHz)
+EbN0R=CNR-Rb //Eb/N0 ratio at IRD(dB)
+
+//Results
+printf("The Azimuth angle of antenna is %.1f degrees" ,Az)
+printf("The Elevaation Angle of Antenna is %.f degrees",El)
+printf("The Range of Antenna is %.f km",d)
+printf("The Eb/N0 ratio at IRD is %.1f dB",EbN0R)
diff --git a/2409/CH16/EX16.2/Ex16_2.sce b/2409/CH16/EX16.2/Ex16_2.sce
new file mode 100755
index 000000000..600107c53
--- /dev/null
+++ b/2409/CH16/EX16.2/Ex16_2.sce
@@ -0,0 +1,55 @@
+
+//Varaible Declaration
+
+R01=42 //Rainfall at earth station(mm/hr)
+p=0.01 //Percentage of time for which rain exceeds
+LE=45 //Latitue of earth station(degrees)
+hR=3.5 //Rain Height(km)
+h0=0 //Mean Sea level(km)
+Ta=272 //
+El=37 //Elevation angle of the antenna(degrees)
+Ts=170 //Total system noise temperature(Kelvin)
+NCR=2.3*10**-9 //Carrier to noise ratio
+fD=12.5 //Frequency of operation(GHz)
+f12=12 //Frequency 12GHz(GHz)
+f15=15 //Frequency 15GHz(GHz)
+//Coefficients for horizontal and vertical polarizations at 12GHz and 15GHz as given in Table 4.2
+
+ah12=0.0188
+av12=0.0168
+bh12=1.217
+bv12=1.2
+
+ah15=0.0367
+av15=0.0335
+bh15=1.154
+bv15=1.128
+
+//Calculation
+
+//Using Interpolation to find coefficients at 12.5 GHz
+ah= ah12+(ah15-ah12)*(fD-f12)/(f15-f12)
+bh= bh12+(bh15-bh12)*(fD-f12)/(f15-f12)
+av=av12+(av15-av12)*(fD-f12)/(f15-f12)
+bv= bv12+(bv15-bv12)*(fD-f12)/(f15-f12)
+
+//Coefficients for circular polarization
+ac=(ah+av)/2
+bc=(ah*bh+av*bv)/(2*ac)
+Ls1=(hR-h0)/sin(El*3.142/180) //Slant Path Length(km)
+Ls= Ls1 //Slant Path Length(km)
+LG= Ls*cos(El*3.142/180) //Horizontal projection of slant path length(km)
+r011=90/(90+4*LG) //Reduction Factor
+r01= r011 //Reduction Factor
+L= Ls1*r01 //Effective path length(km)
+alpha= ac*R01**bc //Specific attenuation(dB/km)
+A= 10**(alpha*L/(10)) //Total Attenuation(dB)
+Trn=Ta*(1-1/A) //noise temperature with effect of rain
+Tscs=Ts
+NCrain=NCR*(A+(A-1)*Ta/Tscs) //Noise to carrier ratio due to rain
+CNrain=-10*log10(NCrain)//Noise to carrier ratio due to rain(dB)
+Rb=10*log10(40*10**6) //Transmission rate(dB)
+EbN0rain= CNrain-Rb //Upper limit of Eb/N0 ratio in prescence of rain(dB)
+
+//Result
+printf("Hence the upper limit for Eb/N0 for given conditions is %.1f dB",EbN0rain)
--
cgit