diff options
Diffstat (limited to '3788')
54 files changed, 604 insertions, 0 deletions
diff --git a/3788/CH10/EX10.6.1/Ex10_6_1.PNG b/3788/CH10/EX10.6.1/Ex10_6_1.PNG Binary files differnew file mode 100644 index 000000000..a95137ba1 --- /dev/null +++ b/3788/CH10/EX10.6.1/Ex10_6_1.PNG diff --git a/3788/CH10/EX10.6.1/Ex10_6_1.sce b/3788/CH10/EX10.6.1/Ex10_6_1.sce new file mode 100644 index 000000000..ec1e5f619 --- /dev/null +++ b/3788/CH10/EX10.6.1/Ex10_6_1.sce @@ -0,0 +1,21 @@ +//Example 10.6.1
+//Length and gain of Satellite
+clc
+clear
+height=750
+theta=10
+SEC=theta + 90
+re=6378
+rs=re + 750
+del = asind(re*(sind(SEC)/rs))
+Y=180-100-del
+Yradian=Y*(%pi/180)
+
+ArcEZ=re*Yradian
+Diameter=2*ArcEZ
+printf("Length of coverage region is %f km \n",Diameter)
+
+Beamwidth=2*del
+Gain=33000/Beamwidth^2
+G=10*log10(Gain)
+printf("Gain of Satellite Antenna is %f dB",G)
diff --git a/3788/CH2/EX2.1.1/Ex2_1_1.sce b/3788/CH2/EX2.1.1/Ex2_1_1.sce new file mode 100644 index 000000000..3d3283b7f --- /dev/null +++ b/3788/CH2/EX2.1.1/Ex2_1_1.sce @@ -0,0 +1,13 @@ +//Example 2.1.1
+// Calculate the Orbital radius for a Geostationary Satellite
+
+//Variable Declaration
+T=86164.09 // Time Period for 1 sidereal day in Sec
+
+//Calculation
+u=(3.986004418*10^5)
+a=((T^2*u)/(4*%pi^2))^(1/3)
+
+//Result
+printf("The Radius of the circular orbit with 1 day period is : %d km",a)
+
diff --git a/3788/CH2/EX2.1.2/Ex2_1_2.sce b/3788/CH2/EX2.1.2/Ex2_1_2.sce new file mode 100644 index 000000000..3378dea56 --- /dev/null +++ b/3788/CH2/EX2.1.2/Ex2_1_2.sce @@ -0,0 +1,20 @@ +//Example 2.1.2
+// Calculate the Orbital radius and Linear Velocity of Shuttle along its orbit
+
+//Variable Declaration
+re=6378.14 //radius of earth in km
+h=250 //altitude in km
+a=re+h
+u=(3.986004418*10^5)
+
+//Calculation
+T=sqrt((4*(%pi^2)*(a^3)/u)) // Period of orbit in Sec
+circum=2*%pi*a //circumference of orbut=2*pi*a in km
+Vs=(2*%pi*a)/T //velocity in km/s
+v=sqrt(u/a) //velocity by alternate method
+
+//Result
+disp(T,'Period of orbit in Sec')
+disp(circum,'Circumference of Orbit in km')
+disp(Vs,'Velocity of Satellite in Orbit in km/s')
+disp(v,' Velocity of Satellite in Orbit in km/s(By Alternate Method)')
diff --git a/3788/CH2/EX2.1.3/Ex2_1_3.sce b/3788/CH2/EX2.1.3/Ex2_1_3.sce new file mode 100644 index 000000000..503bb20b7 --- /dev/null +++ b/3788/CH2/EX2.1.3/Ex2_1_3.sce @@ -0,0 +1,19 @@ +//Example 2.1.3
+//Elliptical Orbit
+
+//Variable Declaration
+re=6378.14
+hp=1000
+ha=4000
+u=3.98600*((10)^5)
+
+//Calculation
+a=(2*re+hp+ha)/2
+T=sqrt((4*(%pi^2)*(a^3)/u))
+
+e=1-((re+hp)/a) //Eccentricity of orbit
+
+//Result
+printf("The period of orbit is : %f seconds",T)
+disp(e,'Eccentricity of orbit')
+
diff --git a/3788/CH2/EX2.2.1/Ex2_2_1.sce b/3788/CH2/EX2.2.1/Ex2_2_1.sce new file mode 100644 index 000000000..2629612f8 --- /dev/null +++ b/3788/CH2/EX2.2.1/Ex2_2_1.sce @@ -0,0 +1,20 @@ +//Example 2.2.1
+//Geostationary Satellite Look Angles
+
+//Variable Declaration
+Le=52.0
+le=0
+longs=66.0
+
+//Calculation
+y=acosd(cosd(Le)*cosd(longs-le))
+printf("Central Angle is %f degrees\n",y)
+
+El=atand((6.6107345-cosd(y))/sind(y))-y
+printf("Elevation Angle is %f degrees\n",(El))
+
+alpha=atand(tand(longs-le)/sind(Le))
+printf("Intermediate Angle is %f degrees\n",alpha)
+
+Az=180-alpha
+printf("Azimuth Angle is %f degrees (clockwise from true north )",Az)
diff --git a/3788/CH2/EX2.3.1/Ex2_3_1.sce b/3788/CH2/EX2.3.1/Ex2_3_1.sce new file mode 100644 index 000000000..08a52f005 --- /dev/null +++ b/3788/CH2/EX2.3.1/Ex2_3_1.sce @@ -0,0 +1,16 @@ +//Example 2.3.1
+//Drift with a Geostationary Satellite
+
+//Variables
+T=86400 //T corresponds to 1 solar day
+u=(3.986004418*10^5)
+
+//Calculation
+//Part1 Orbital Radius
+a=((T^2*u)/(4*%pi^2))^(1/3) //4 pi^2=39.4784716
+printf("The radius of the circular orbit is : %f km\n",a)
+
+//Part2 Rate of Drift
+drift=360*(235.9/T)
+printf("The drift is %f degrees/day",drift)
+
diff --git a/3788/CH2/EX2.6.1/Ex2_6_1.sce b/3788/CH2/EX2.6.1/Ex2_6_1.sce new file mode 100644 index 000000000..f90ae3efc --- /dev/null +++ b/3788/CH2/EX2.6.1/Ex2_6_1.sce @@ -0,0 +1,28 @@ +//Example 2.6.1
+//Doppler Shift for A leo Satellite
+
+//Part1 Velocity of Satellite in Orbit
+//Variables
+re=6378 //radius of earth in km
+h=1000 //altitude in km
+a=re+h
+u=(3.986004418*10^5)
+
+//Calculation
+//Part1 Velocity of Satellite in Orbit
+T=sqrt((4*(%pi^2)*(a^3)/u))
+circum=2*%pi*a
+vs=circum/T
+printf("Velocity of Satellite is %f km/s\n",vs)
+
+//Part2 Component of velocity towards the observer
+D=re/a
+vr=(vs*(re/a))
+printf("Velocity of Satellite towards observer is %f km/s\n",vr)
+
+//Part3 Doppler shift in received signal
+VT=6354
+lam=0.1132
+Df=(VT/lam)/1000
+disp(Df,'Doppler shift in kHz')
+
diff --git a/3788/CH4/EX4.2.1/Ex4_2_1.PNG b/3788/CH4/EX4.2.1/Ex4_2_1.PNG Binary files differnew file mode 100644 index 000000000..8b4b221e1 --- /dev/null +++ b/3788/CH4/EX4.2.1/Ex4_2_1.PNG diff --git a/3788/CH4/EX4.2.1/Ex4_2_1.sce b/3788/CH4/EX4.2.1/Ex4_2_1.sce new file mode 100644 index 000000000..745fd4be2 --- /dev/null +++ b/3788/CH4/EX4.2.1/Ex4_2_1.sce @@ -0,0 +1,21 @@ +
+//example 4.2.1
+//Calculate the Flux Density and Power received
+
+//Variables
+clc
+clear
+D = 40000
+gain = 17 //gain is in dB
+Gt = 50
+A = 10 //effective area of antenna
+Pt = 10 //transmitted power
+R = 4*(10)^7
+
+//Calculation
+F=(Pt*Gt)/(4*%pi*R^2) //flux density equation
+Pr = Pt/A //Received Power
+
+//Result
+printf("The flux density is %f Watts per sqm \n",F)
+printf("The power received by antenna is %f Watts",Pr)
diff --git a/3788/CH4/EX4.2.2/Ex4_2_2.PNG b/3788/CH4/EX4.2.2/Ex4_2_2.PNG Binary files differnew file mode 100644 index 000000000..1cc54dd4f --- /dev/null +++ b/3788/CH4/EX4.2.2/Ex4_2_2.PNG diff --git a/3788/CH4/EX4.2.2/Ex4_2_2.sce b/3788/CH4/EX4.2.2/Ex4_2_2.sce new file mode 100644 index 000000000..0a8d1be7e --- /dev/null +++ b/3788/CH4/EX4.2.2/Ex4_2_2.sce @@ -0,0 +1,18 @@ +// Example 4.2.2
+// Calculate the power received by the antenna while the gain of
+// receiving antenna is 52.3 dB
+
+//variables
+clc
+clear
+EIRP=27.0
+Gr=52.3
+R=4*10^7
+lam=2.727*10^(-2)
+
+//calculation
+pathloss=20*log10((4*%pi*R)/lam) //Finding path loss
+Pr=EIRP+Gr-pathloss //Finding Power received
+
+//Result
+printf("Power received at the antenna is %f dBW",Pr)
diff --git a/3788/CH4/EX4.3.1/Ex4_3_2.PNG b/3788/CH4/EX4.3.1/Ex4_3_2.PNG Binary files differnew file mode 100644 index 000000000..29698fd8f --- /dev/null +++ b/3788/CH4/EX4.3.1/Ex4_3_2.PNG diff --git a/3788/CH4/EX4.3.1/Ex4_3_2.sce b/3788/CH4/EX4.3.1/Ex4_3_2.sce new file mode 100644 index 000000000..cc5b4dc99 --- /dev/null +++ b/3788/CH4/EX4.3.1/Ex4_3_2.sce @@ -0,0 +1,22 @@ +//Example 4.3.2
+// Calculate the New System noise temperature
+
+//Variables
+clc
+clear
+Gain = 50
+G1 = 0.631
+attenuation = 2
+Tm = 500
+TIF = 1000
+Tp = 300
+
+//Calculation
+Twg = Tp*(1-G1)
+Tin = G1*25
+T = (Tin + Twg +Gain + (Tm/(10^5)) + (TIF/(10^4)))
+Ts = (T/G1)
+
+//Result
+printf("Waveguide Noise temperature is %fK \n",Twg)
+printf("The new System noise temperature is %fK",Ts)
diff --git a/3788/CH4/EX4.3.3/Ex4_3_3.PNG b/3788/CH4/EX4.3.3/Ex4_3_3.PNG Binary files differnew file mode 100644 index 000000000..01a6a72f1 --- /dev/null +++ b/3788/CH4/EX4.3.3/Ex4_3_3.PNG diff --git a/3788/CH4/EX4.3.3/Ex4_3_3.sce b/3788/CH4/EX4.3.3/Ex4_3_3.sce new file mode 100644 index 000000000..4925a9693 --- /dev/null +++ b/3788/CH4/EX4.3.3/Ex4_3_3.sce @@ -0,0 +1,11 @@ +//Example 4.3.3
+// Calculate Noise Temperature
+
+clc
+clear
+T0 = 290 //reference temp usually 290K
+NF = 1.78 //Noise Figure
+Td = T0*(NF - 1)
+
+//Result
+printf("Value of noise temperature is %dK",Td)
diff --git a/3788/CH4/EX4.3.4/Ex4_3_4.PNG b/3788/CH4/EX4.3.4/Ex4_3_4.PNG Binary files differnew file mode 100644 index 000000000..110fda409 --- /dev/null +++ b/3788/CH4/EX4.3.4/Ex4_3_4.PNG diff --git a/3788/CH4/EX4.3.4/Ex4_3_4.sce b/3788/CH4/EX4.3.4/Ex4_3_4.sce new file mode 100644 index 000000000..f40655044 --- /dev/null +++ b/3788/CH4/EX4.3.4/Ex4_3_4.sce @@ -0,0 +1,25 @@ +//Example 4.3.4
+// Calculate the GT ratio and also if the Noise temperature is risen
+
+//Variables
+clc
+clear
+D = 30
+na = 0.68
+lam = 0.0723
+temp0 = 79
+temp1 = 88
+
+//Calculation
+Gr = na*((%pi*D)/lam)^2
+GrdB = 10*log10(Gr)
+Ts0 = 10*log10(temp0)
+Ts1 = 10*log10(temp1)
+GT = GrdB - Ts0
+GTnew = GrdB - Ts1
+
+//Result
+printf("The Gain of antenna is %f dB \n",GrdB)
+printf("System noise temperature is %f dBK \n",Ts1)
+printf("Earth station GT ratio is %f dB/K \n",GT)
+printf("If the temperature rises to 88K ,then new GT ratio is %f dB/K",GTnew)
diff --git a/3788/CH4/EX4.6.1/Ex4_6_1.PNG b/3788/CH4/EX4.6.1/Ex4_6_1.PNG Binary files differnew file mode 100644 index 000000000..04d535a1e --- /dev/null +++ b/3788/CH4/EX4.6.1/Ex4_6_1.PNG diff --git a/3788/CH4/EX4.6.1/Ex4_6_1.sce b/3788/CH4/EX4.6.1/Ex4_6_1.sce new file mode 100644 index 000000000..d5ea32ff0 --- /dev/null +++ b/3788/CH4/EX4.6.1/Ex4_6_1.sce @@ -0,0 +1,20 @@ +//Example 4.6.1
+//Calculate the power op of an uplink transmitter
+//Variables
+clc
+clear
+Pr = -127
+Gt = 50
+Gr = 26
+Lp = 207.2
+Lta = 1.5
+Lat = 0.5
+Lpt = -2
+Pin = 0-127
+
+//Result
+Pt = Pr - Gt - Gr + Lp + Lat + Lta - Lpt
+printf("The transmitter output power is %f dBW \n",Pt)
+Rainattenuation = 7
+Ptrain = Pt + Rainattenuation
+printf("The transmitter output power due to fading of rain is %f dBW",Ptrain)
diff --git a/3788/CH4/EX4.7.1/Ex4_7_1.PNG b/3788/CH4/EX4.7.1/Ex4_7_1.PNG Binary files differnew file mode 100644 index 000000000..8ec9c5468 --- /dev/null +++ b/3788/CH4/EX4.7.1/Ex4_7_1.PNG diff --git a/3788/CH4/EX4.7.1/Ex4_7_1.sce b/3788/CH4/EX4.7.1/Ex4_7_1.sce new file mode 100644 index 000000000..6f291cbff --- /dev/null +++ b/3788/CH4/EX4.7.1/Ex4_7_1.sce @@ -0,0 +1,17 @@ +//Example 4.7.1
+//Calculate the overall CN ratio
+
+//Variables
+clc
+clear
+CNdnratio = 100
+CNupratio = 100
+
+CIratio = 24
+CIratioindB = 0.004
+CN0 = (1/((1/CNupratio) + (1/CNdnratio)))
+CN1 = (1/((1/CNupratio) + (1/CNdnratio) + CIratioindB))
+
+//result
+printf("CN0 ratio is %f\n",CN0)
+printf("The overall CN0 ratio at earth receiveing station is %f ",CN1)
diff --git a/3788/CH5/EX5.2.1/Ex5_2_1.PNG b/3788/CH5/EX5.2.1/Ex5_2_1.PNG Binary files differnew file mode 100644 index 000000000..88575cf23 --- /dev/null +++ b/3788/CH5/EX5.2.1/Ex5_2_1.PNG diff --git a/3788/CH5/EX5.2.1/Ex5_2_1.sce b/3788/CH5/EX5.2.1/Ex5_2_1.sce new file mode 100644 index 000000000..430068de7 --- /dev/null +++ b/3788/CH5/EX5.2.1/Ex5_2_1.sce @@ -0,0 +1,16 @@ +//Example
+//Calculate the baseband SN ratio for the video signal
+//Variables
+clc
+clear
+Fmax = 4.2
+RFbw = 30
+CNratio = 15
+P = 9
+Q = 8
+
+//Result
+delFpk = (RFbw/2) - Fmax
+Brf = 2*(delFpk + Fmax)
+SN = CNratio + 10*log10(RFbw/Fmax) + 20*log10(delFpk/Fmax)+ 1.5 +P + Q
+printf("The baseband SN ratio for the video signal is %f",SN)
diff --git a/3788/CH5/EX5.2.2/Ex5_2_2.PNG b/3788/CH5/EX5.2.2/Ex5_2_2.PNG Binary files differnew file mode 100644 index 000000000..f5490c622 --- /dev/null +++ b/3788/CH5/EX5.2.2/Ex5_2_2.PNG diff --git a/3788/CH5/EX5.2.2/Ex5_2_2.sce b/3788/CH5/EX5.2.2/Ex5_2_2.sce new file mode 100644 index 000000000..1ce9fcd18 --- /dev/null +++ b/3788/CH5/EX5.2.2/Ex5_2_2.sce @@ -0,0 +1,17 @@ +//Example 5.2.2
+//Calculate the baseband SN ratio for the Voice channel
+
+//Variables
+clc
+clear
+Fmax = 3.4
+Brf = 45
+Rs = 9.6
+P = 8.8
+SNratio = 7
+CNratio = 13
+
+//result
+delFpk = Brf/2 - Fmax
+SNfm = CNratio + 10*log10(Brf/Fmax) + 20*log10(delFpk/Fmax) + 1.8 + P
+printf("The baseband SN ratio for the voice channel is %f dB ",SNfm)
diff --git a/3788/CH5/EX5.2.3/Ex5_2_3.PNG b/3788/CH5/EX5.2.3/Ex5_2_3.PNG Binary files differnew file mode 100644 index 000000000..72270fbfc --- /dev/null +++ b/3788/CH5/EX5.2.3/Ex5_2_3.PNG diff --git a/3788/CH5/EX5.2.3/Ex5_2_3.sce b/3788/CH5/EX5.2.3/Ex5_2_3.sce new file mode 100644 index 000000000..76100e6c6 --- /dev/null +++ b/3788/CH5/EX5.2.3/Ex5_2_3.sce @@ -0,0 +1,16 @@ +//Example 5.2.3
+//Calculate the SN ratio if CN = 10dB
+
+//Variables
+clc
+clear
+delFpk = 3.6
+Fmax = 4.8
+CN = 10
+delFpeak = 3.6
+
+//result
+Brf = 2*(delFpk + Fmax)
+SNout = CN +10*log10(Brf/Fmax) + 20*log10(delFpeak/Fmax) + 1.8
+printf("The SN ratio is %f dB,if the CN ratio for \n the signal from the satellite is 10dB",SNout)
+
diff --git a/3788/CH5/EX5.3.1/Ex5_3_1.PNG b/3788/CH5/EX5.3.1/Ex5_3_1.PNG Binary files differnew file mode 100644 index 000000000..c116b57a2 --- /dev/null +++ b/3788/CH5/EX5.3.1/Ex5_3_1.PNG diff --git a/3788/CH5/EX5.3.1/Ex5_3_1.sce b/3788/CH5/EX5.3.1/Ex5_3_1.sce new file mode 100644 index 000000000..8d24c65fc --- /dev/null +++ b/3788/CH5/EX5.3.1/Ex5_3_1.sce @@ -0,0 +1,11 @@ +//Example 5.3.1
+//Calculate the pulse rate for the link
+//Variables
+clc
+clear
+Bocc = 100000.00
+alpha = 0.5
+
+//Result
+Rs = Bocc/(1 + alpha)
+printf("The correct symbol rate is %f symbols/sec",Rs)
diff --git a/3788/CH5/EX5.3.2/Ex5_3_2.PNG b/3788/CH5/EX5.3.2/Ex5_3_2.PNG Binary files differnew file mode 100644 index 000000000..058e0efce --- /dev/null +++ b/3788/CH5/EX5.3.2/Ex5_3_2.PNG diff --git a/3788/CH5/EX5.3.2/Ex5_3_2.sce b/3788/CH5/EX5.3.2/Ex5_3_2.sce new file mode 100644 index 000000000..1fbf66a3e --- /dev/null +++ b/3788/CH5/EX5.3.2/Ex5_3_2.sce @@ -0,0 +1,16 @@ +//Example 5.3.2
+//Claculate the bandwidth and frequency range
+//Variables
+clc
+clear
+Rs = 16
+fc = 14.125
+alpha = 0.25
+
+//result
+Bocc = Rs*(1 + alpha)
+fl = fc - (Rs/2)*(1+alpha)
+fh = fc + (Rs/2)*(1+alpha)
+
+printf("The bandwidth occupied by RF signal is %f Mhz\n",Bocc)
+printf("The frequecny range is from %f Ghz to %f Ghz",fl,fh)
diff --git a/3788/CH5/EX5.3.3/Ex5_3_3.PNG b/3788/CH5/EX5.3.3/Ex5_3_3.PNG Binary files differnew file mode 100644 index 000000000..c8e664968 --- /dev/null +++ b/3788/CH5/EX5.3.3/Ex5_3_3.PNG diff --git a/3788/CH5/EX5.3.3/Ex5_3_3.sce b/3788/CH5/EX5.3.3/Ex5_3_3.sce new file mode 100644 index 000000000..3f59009a3 --- /dev/null +++ b/3788/CH5/EX5.3.3/Ex5_3_3.sce @@ -0,0 +1,13 @@ +//example 5.3.3
+//calculate bit rate for BPSK and QPSK
+//variables
+clc
+clear
+BW = 36
+alpha = 0.4
+
+//result
+RsBPSK = BW/(1 + alpha)
+RsQPSK = 2*RsBPSK
+printf("The maximum symbol rate for BPSK RF link is %f Msps \n",RsBPSK)
+printf("The maximum symbol rate for QPSK RF link is %f Msps ",RsQPSK)
diff --git a/3788/CH5/EX5.4.1/Ex5_4_1.PNG b/3788/CH5/EX5.4.1/Ex5_4_1.PNG Binary files differnew file mode 100644 index 000000000..a68228fd1 --- /dev/null +++ b/3788/CH5/EX5.4.1/Ex5_4_1.PNG diff --git a/3788/CH5/EX5.4.1/Ex5_4_1.sce b/3788/CH5/EX5.4.1/Ex5_4_1.sce new file mode 100644 index 000000000..ff73ca724 --- /dev/null +++ b/3788/CH5/EX5.4.1/Ex5_4_1.sce @@ -0,0 +1,34 @@ +//example 5.4.1
+//Calculate the bitrate ,Symbol rate BW and BER values
+
+//Variables
+clc
+clear
+CN = 25
+NoiseBw = 1.0
+r=0.3
+Rs = 1
+Bocc = Rs*(1+r)
+rainattenuation = 3
+printf("The occupied bandwidth of the RF signal is %f Mhz \n",Bocc)
+//BPSK
+Rb = 1
+printf("The bit rate is %f Mbps \n",Rb)
+BERclearair = erfc((2*CN)^(1/2))
+printf("BER in clear air for BPSK is %f \n Since the all BER values are -ve high orders \n the BER values are shown zero\n",BERclearair)
+
+//QPSK
+Rb1 = 2*Rs
+printf("The bit rate for QPSK is %f Mbps \n",Rb1)
+BERclearair1 = erfc((CN)^(1/2))
+printf("BER in clear air for QPSK is %f \n",BERclearair1)
+
+CN1 = CN - rainattenuation
+
+//BPSK
+BERrain = erfc((2*CN1)^(1/2))
+printf("BER in rain for BPSK is %f \n",BERrain)
+
+//QPSK
+BERrain1 = erfc((CN1)^(1/2))
+printf("BER in rain for BPSK is %f \n",BERrain1)
diff --git a/3788/CH5/EX5.4.2/Ex5_4_2.PNG b/3788/CH5/EX5.4.2/Ex5_4_2.PNG Binary files differnew file mode 100644 index 000000000..0825593d3 --- /dev/null +++ b/3788/CH5/EX5.4.2/Ex5_4_2.PNG diff --git a/3788/CH5/EX5.4.2/Ex5_4_2.sce b/3788/CH5/EX5.4.2/Ex5_4_2.sce new file mode 100644 index 000000000..88bd4968b --- /dev/null +++ b/3788/CH5/EX5.4.2/Ex5_4_2.sce @@ -0,0 +1,28 @@ +//Example 5.4.2
+//Calculate Bit rate,symbol rate ,BER values for BPSK and QPSK
+
+//Variables
+clc
+clear
+BW = 10
+alpha = 0.25
+CN = 16
+marginbpsk = 0.8
+marginqpsk = 1.2
+Rs = BW/(1+alpha)
+RsQPSk = 2*Rs
+
+printf("The Symbol rate for satelite link is %f Msps \n",Rs)
+printf("The bit rate for BPSK is %f MBps \n",Rs)
+printf("The bit rate for QPSK is %f MBps \n",RsQPSk)
+
+CNeff = 10^((CN - marginbpsk)/10)
+BERBPSk = erfc((2*CNeff)^(1/2))
+printf("C/N effective for BPSk is %f \n ",CNeff)
+printf("The BER value for BPSK is %f \n",BERBPSk)
+
+
+CNeff1 = 10^((CN - marginqpsk)/10)
+BERQPSk = erfc((CNeff)^(1/2))
+printf("C/N effective for QPSk is %f \n ",CNeff1)
+printf("The BER value for QPSK is %f ",BERQPSk)
diff --git a/3788/CH6/EX6.2.1/Ex6_2_1.PNG b/3788/CH6/EX6.2.1/Ex6_2_1.PNG Binary files differnew file mode 100644 index 000000000..5ebe52ea0 --- /dev/null +++ b/3788/CH6/EX6.2.1/Ex6_2_1.PNG diff --git a/3788/CH6/EX6.2.1/Ex6_2_1.sce b/3788/CH6/EX6.2.1/Ex6_2_1.sce new file mode 100644 index 000000000..eee07265c --- /dev/null +++ b/3788/CH6/EX6.2.1/Ex6_2_1.sce @@ -0,0 +1,42 @@ +//example 6.2.1
+//calculate power level at op of transponder
+//variables
+clc
+clear
+pearth = 500
+gain = 105
+backoff = 3
+outputpower = 40
+BWStA = 15
+BWStB = 10
+BWStC = 5
+Pt = 20
+EIRPa = 3.0
+EIRPb = 4.8
+EIRPc = 7.8
+
+PtindB=10*log10(outputpower) - backoff
+printf("Power of tansponder is %fdBW \n",PtindB)
+BWt = BWStA + BWStB + BWStC
+PtA = 10*log10((BWStA/BWt)*Pt)
+PtB = 10*log10((BWStB/BWt)*Pt)
+PtC = 10*log10((BWStC/BWt)*Pt)
+printf("Transponder power output allocated to StA is %f dBW \n",PtA)
+printf("Transponder power output allocated to StB is %f dBW\n",PtB)
+printf("Transponder power output allocated to StC is %f dBW\n",PtC)
+
+PinA = PtA - gain
+PinB = PtB - gain
+PinC = PtC - gain
+printf("Transponder input power for StA signal is %f dBW\n",PinA)
+printf("Transponder input power for StB signal is %f dBW\n",PinB)
+printf("Transponder input power for StC signal is %f dBW \n",PinC)
+
+Pte = 10*log10(250)
+PStA = Pte - EIRPa
+PStB = Pte - EIRPb
+PStC = Pte - EIRPc
+printf("The Earth Station A transmit power is %f dBW \n",PStA)
+printf("The Earth Station B transmit power is %f dBW \n",PStB)
+printf("The Earth Station C transmit power is %f dBW\n",PStC)
+
diff --git a/3788/CH6/EX6.3.1/Ex6_3_1.PNG b/3788/CH6/EX6.3.1/Ex6_3_1.PNG Binary files differnew file mode 100644 index 000000000..780c395f0 --- /dev/null +++ b/3788/CH6/EX6.3.1/Ex6_3_1.PNG diff --git a/3788/CH6/EX6.3.1/Ex6_3_1.sce b/3788/CH6/EX6.3.1/Ex6_3_1.sce new file mode 100644 index 000000000..917bc811a --- /dev/null +++ b/3788/CH6/EX6.3.1/Ex6_3_1.sce @@ -0,0 +1,18 @@ +//Example 6.3.1
+//Calculate the no of digital channels
+//Variables
+clc
+clear
+Tframe=2000 //20ms
+N = 5 //5us
+tg=5
+tpre = 20
+datarate = 64000
+QPSKsymbol = 2
+QPSKTxburst = 30 //30mbaud
+
+Td = (Tframe-N*(tg+tpre))/N
+Rb = QPSKsymbol*QPSKTxburst
+Cb = Td*(((Rb*(10)^6))/Tframe)
+n = (Cb)/64000
+printf("The number of 64kbps digital sppech channels that \n can be carried by one earth station is %d",n)
diff --git a/3788/CH8/EX8.2.1/Ex8_2_1.PNG b/3788/CH8/EX8.2.1/Ex8_2_1.PNG Binary files differnew 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 Binary files differnew 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 Binary files differnew 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 Binary files differnew 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 Binary files differnew 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 Binary files differnew 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)
|
