diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2409/CH16 | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '2409/CH16')
-rwxr-xr-x | 2409/CH16/EX16.1/Ex16_1.sce | 46 | ||||
-rwxr-xr-x | 2409/CH16/EX16.2/Ex16_2.sce | 55 |
2 files changed, 101 insertions, 0 deletions
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) |