14 files changed, 298 insertions, 0 deletions
diff --git a/3710/CH8/EX8.1/Ex8_1.sce b/3710/CH8/EX8.1/Ex8_1.sce
new file mode 100644
index 000000000..9e9241ca5
--- /dev/null
+++ b/3710/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,38 @@
+//Example 8.1, Page Number 364
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Field Decay Factor Calculation
+clc;
+clear all;
+
+n1=1.5 //Refraction index of glass/air interface
+n2=1 //Refraction index of air
+
+c1=1*2*3.14*n2
+a=60//given
+d=(a*(%pi/180))
+a1=(sin(d))**2  //By equation 8.8
+a2=(n1/n2)**2
+a3=(a1*a2)-1
+a4=sqrt(a3)
+a5=a4*c1
+
+//assuming wavelength = field distance
+F=exp(-1*a5) //F is the field decay factor
+disp(F,"The Field Decay factor at 60 degrees is:")
+
+//at theta1=42 degrees
+c1=1*2*3.14*n2
+//disp(c1)
+a=42//given
+d=a*(%pi/180)
+b1=(sin(d))**2 //By equation 8.8
+b2=(n1/n2)**2
+b3=(b1*b2)-1
+b4=sqrt(b3)
+c4=b4*c1
+
+//assuming wavelength = field distance
+F1=exp(-1*c4)//F1 is the field decay factor at 42 degrees
+F1=fpround(F1,3)
+
+disp(F1," The Field Decay factor at 42 degrees is:")
diff --git a/3710/CH8/EX8.10/Ex8_10.sce b/3710/CH8/EX8.10/Ex8_10.sce
new file mode 100644
index 000000000..35bd63586
--- /dev/null
+++ b/3710/CH8/EX8.10/Ex8_10.sce
@@ -0,0 +1,30 @@
+//Example 8.10, Page Number 390
+//Material dispersion for a laser and LED source
+clc;
+clf();
+//Fig 8.26 Page No 390
+X=[0.5,0.75,1,1.25,1.4,1.5,1.75,2,2.5]//Values as observed from graph
+V=[0.07,0.04,0.02,0.0,-0.01,-0.02,-0.03,-0.04,-0.06]//Values as observed from graph
+plot(X,V);
+xlabel("Wavelength (um)")
+ylabel("D (Dimensionless Quantity)")
+title("Fig 8.26")
+
+l=10.0**3 //length of the fiber in kilometer
+w1=850.0*(10**-9) //Wavelength in meter
+lw=50.0*(10**-9) //Linewidth in meter
+w2=1550.0*(10**-9) //Wavelength 2 in meter
+lw2=3.0*(10**-9) //Linewidth 2 in meter
+c=3.0*(10**8) //Speed of light in meters per second
+
+
+mprintf(" The Dimensionless Quantity is calculated from Figure 8.26\n")
+d1=2.14*(10**-2) //For w1 after observation
+d2=-1.02*(10**-2) //For w2 after observation
+
+//From equation 8.34
+t1=(l/c)*d1*(lw/w1) //t1 is the material dispersion effects for w1
+t2=(l/c)*d2*(lw2/w2) //t2 is the material dispersion effects for w2
+
+mprintf(" (a) The Material Dispersion Effect Parameter for the LED is %.1e s\n",t1)
+mprintf(" (b) The Material Dispersion Effect Parameter for the Laser is %.1e s",-1*t2)
diff --git a/3710/CH8/EX8.11/Ex8_11.sce b/3710/CH8/EX8.11/Ex8_11.sce
new file mode 100644
index 000000000..e1057542a
--- /dev/null
+++ b/3710/CH8/EX8.11/Ex8_11.sce
@@ -0,0 +1,34 @@
+//Example 8.11, Page Number 396
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Rayleigh scattering loss
+clc;
+
+l=1*(10**-6) //Wavelength in meter
+l1=10**3 //Length of the Fiber in kilometer
+n=1.45 //Refractive Index
+p=0.286
+B=7*(10**-11) //in meter square per Newton
+T=1400 //in kelvin
+k=1.38*(10**-23) //Boltzman Constant in meter square kilogram per second square Kelvin
+
+//From equation 8.38
+ar=((8*(%pi**3))/(3*(l**4)))*((n**8)*B*T*k*(p**2))    //ar is the attenuation due to Rayleigh scattering per meter
+ar=fpround(ar,6)
+mprintf("The Attenuation due to Rayleigh scattering in a silica fiber is %.2e /m\n",ar)
+
+i1=(ar*-1*l1)
+j=exp(i1)
+L=-10*log10(j) //L is the given loss generated from attenuation in dB
+L=fpround(L,5)
+mprintf(" The Total loss due to Attenuation is %.2f db/km\n",L)
+
+//at a wavelength of 1.55 micro meter
+l=1.55*(10**-6) //new Wavelength in meters
+
+ar=((8*(%pi**3))/(3*(l**4)))*((n**8)*B*T*k*(p**2)) //ar is the attenuation due to Rayleigh scattering per meter
+ar=fpround(ar,6)
+mprintf(" The Attenuation due to Rayleigh scattering in a silica fiber is %.1e /m\n",ar)
+
+L=-10*log10(exp(ar*-1*l1))    //L is the given loss generated from attenuation in dB
+L=fpround(L,2)
+mprintf(" The Total loss due to Attenuation is %.2f db/km",L)
diff --git a/3710/CH8/EX8.12/Ex8_12.sce b/3710/CH8/EX8.12/Ex8_12.sce
new file mode 100644
index 000000000..0d706398b
--- /dev/null
+++ b/3710/CH8/EX8.12/Ex8_12.sce
@@ -0,0 +1,20 @@
+//Example 8.12, Page Number 398
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+clc;
+
+n=1.48//Refractive index of fiber
+n0=1//Refractive index between the fibers in air
+
+//From equation 8.39
+Rf=(((n-1)/(n+1))**2) //Rf is the fraction of light
+Rf=fpround(Rf,4)
+
+Tf=((1-Rf)**2)  //Tf is the total transmission for each face due to Fresnel reflection
+Tf=fpround(Tf,4)
+
+L=-10*log10(Tf) //L is the Transmission Loss
+L=fpround(L,2)
+
+mprintf("The Fraction of light reflected back at each end is %.4f\n",Rf)
+mprintf(" The Total Transmission for each face due to Fresnel Reflection is %.3f\n",Tf)
+mprintf(" The Total Transmission loss is %.2f db",L)
diff --git a/3710/CH8/EX8.13/Ex8_13.sce b/3710/CH8/EX8.13/Ex8_13.sce
new file mode 100644
index 000000000..eba00e017
--- /dev/null
+++ b/3710/CH8/EX8.13/Ex8_13.sce
@@ -0,0 +1,16 @@
+//Example 8.13, Page Number 399
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Transmission loss from lateral misalignment
+clc;
+
+l=0.1 //where l=D/2a and occurs due to lateral misalignment where D is the lateral displacement and a is the fiber core radius
+
+//From equation 8.40
+T=(2/%pi)*(acos(l)-l*(sqrt(1-(l**2))))
+L=-10*log10(T) //L is the Transmission loss in dB
+L=fpround(L,2)
+mprintf("The Total Transmission loss is %.2f dB\n",L)
+
+//Including Fresnel loss from the previous example
+LT=L+0.33
+mprintf(" The Total Transmission loss including Fresnel loss is %.2f dB",LT)
diff --git a/3710/CH8/EX8.14/Ex8_14.sce b/3710/CH8/EX8.14/Ex8_14.sce
new file mode 100644
index 000000000..0dd06e2cf
--- /dev/null
+++ b/3710/CH8/EX8.14/Ex8_14.sce
@@ -0,0 +1,19 @@
+//Example 8.14, Page Number 404
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Fiber Couple Losses
+clc;
+
+d1=200*(10**-6) //Core Diameter in meters
+d2=250*(10**-6) //Core + Cladding Diameter in meters
+
+d3=3*d2//Mixing rod diameter in meters
+
+//Power Levels P1=(B*3.14*(3d2**2))/4 & P2=(B*3.14*(d1**2)/)/4 
+
+//From equation 8.41
+L=-10*log10((d1**2)/(d3**2)) //L is the Insertion Loss in dB
+Le=-10*log10((7*(d1**2))/(d3**2)) //Le is the Excess loss in dB
+L=fpround(L,1)
+Le=fpround(Le,2)
+mprintf("The Insertion loss is %.1f dB\n",L)
+mprintf(" The Excess Loss is %d dB",Le)
diff --git a/3710/CH8/EX8.2/Ex8_2.sce b/3710/CH8/EX8.2/Ex8_2.sce
new file mode 100644
index 000000000..cd07c94e8
--- /dev/null
+++ b/3710/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,18 @@
+//Example 8.2, Page Number 369
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Number of modes in a planar dielectic guide
+clc;
+
+w=100*(10**-6) //Thickness in meter
+n1=1.48//refractive index
+n2=1.46//refractive index
+l=1*(10**-6) //Wavelength of light in meters
+
+V=((%pi*w)/l)*sqrt((n1**2)-(n2**2))
+V1=(2*V)/%pi
+mprintf("The value of V is:%.2f\n",V)
+
+//from equation 8.15
+N=2*(1+int(V1)) //N is the total number of possible modes
+N=fpround(N,0)
+mprintf(" The Total number of possible modes is:%d",N)
diff --git a/3710/CH8/EX8.3/Ex8_3.sce b/3710/CH8/EX8.3/Ex8_3.sce
new file mode 100644
index 000000000..c796b705f
--- /dev/null
+++ b/3710/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,13 @@
+//Example 8.3, Page Number 370
+//Single mode guide dimensions
+clc;
+
+n1=1.48//refractive index
+n2=1.46//refractive index
+l=1*(10**-6) //Wavelength in meter
+
+//from equation 8.16
+
+d=l*(1/(2*sqrt((n1**2)-(n2**2))))
+
+disp("meters",d,"The Waveguide core thickness must be less than")
diff --git a/3710/CH8/EX8.4/Ex8_4.sce b/3710/CH8/EX8.4/Ex8_4.sce
new file mode 100644
index 000000000..d7686abca
--- /dev/null
+++ b/3710/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,18 @@
+//Example 8.4, Page Number 375
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Number of modes in a fiber
+clc;
+
+n1=1.48//refractive index 
+n2=1.46//refractive index
+l=0.9*(10**-6) //Wavelength in meters
+r=100*(10**-6) //Core radius in meter
+
+V=(2*%pi*r*sqrt((n1**2)-(n2**2)))/l 
+V=fpround(V,1)
+N=(V**2)/2 //N is the Number of Modes
+N=fpround(N,1)
+mprintf("The Value of V is:%.1f\n",V)
+mprintf(" The Number of Modes able to propagate is:%.1f",N)
+
+//The difference arising is due to approximation
diff --git a/3710/CH8/EX8.5/Ex8_5.sce b/3710/CH8/EX8.5/Ex8_5.sce
new file mode 100644
index 000000000..fe8bfa60d
--- /dev/null
+++ b/3710/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,30 @@
+//Example 8.5, Page Number 377
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+clc;
+
+n1=1.48
+n2=1.46 
+no=1
+
+theta=asin(n2/n1)
+theta=theta*(180/%pi)
+
+//Using equation 8.22
+delta=((n1**2)-(n2**2))/(2*(n1**2))
+delta=fpround(delta,4)
+
+//Using expression 8.22a
+deltaa=(n1-n2)/n1
+deltaa=fpround(deltaa,4)
+
+NA=sqrt((n1**2)-(n2**2))
+NA=fpround(NA,3)
+
+a=asin(NA) //a is a fiber parameter
+a=a*(180/%pi)
+a=fpround(a,2)
+
+mprintf("The Delta Parameter is:%.4f\n",delta)
+mprintf(" The Approximated Delta Parameter is:%.4f\n",deltaa)
+mprintf(" The Numerical aperture is:%.3f\n",NA)
+mprintf(" The Value of Alpha Max is:%d degrees",a)
diff --git a/3710/CH8/EX8.6/Ex8_6.sce b/3710/CH8/EX8.6/Ex8_6.sce
new file mode 100644
index 000000000..608b297f5
--- /dev/null
+++ b/3710/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,12 @@
+//Example 8.6, Page Number 380
+//Intermodal Dispersion in step index fibers
+clc;
+
+n1=1.48//refractive index 
+n2=1.46//refractive index
+L=1*(10**3) //Length of Fiber in kilometer
+c=3*(10**8) //Speed of Light in meters per second
+
+//Using equation 8.24
+td=((L*n1)/(c*n2))*(n1-n2) //t is the time difference due tp intermodal dispersion
+disp(td,"The Time difference due to Intermodal dispersion in s is:")
diff --git a/3710/CH8/EX8.7/Ex8_7.sce b/3710/CH8/EX8.7/Ex8_7.sce
new file mode 100644
index 000000000..ed7fd1882
--- /dev/null
+++ b/3710/CH8/EX8.7/Ex8_7.sce
@@ -0,0 +1,18 @@
+//Example 8.7, Page Number 383
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+//Intermodal Dispersion
+clc;
+
+n1=1.48//refractive index 
+n2=1.46//refractive index
+l=10**3 //Length of the Fiber in kilometer
+c=3*(10**8) //Speed of light meters per second
+
+delta=((n1**2)-(n2**2))/(2*(n1**2)) //delta is one of the fiber parameter
+delta=fpround(delta,4)
+
+//Using Equation 8.26
+t=(l*n1*(delta**2))/(c*8) //t is the minimum pulse broadening
+
+mprintf("The Delta Paramter is:%.4f\n",delta)
+mprintf(" The Minimum Pulse Broadening is %.2e s",t)
diff --git a/3710/CH8/EX8.8/Ex8_8.sce b/3710/CH8/EX8.8/Ex8_8.sce
new file mode 100644
index 000000000..e763b96fa
--- /dev/null
+++ b/3710/CH8/EX8.8/Ex8_8.sce
@@ -0,0 +1,20 @@
+//Example 8.8, Page Number 387
+//Single Mode Fiber in graded index fibers
+clc;
+
+n1=1.48//refractive index
+n2=1.46//refractive index 
+NA=0.242//numerical aperture
+l=1.5*(10**-6) //Wavelength of radiation in meters
+
+//From equation 8.27
+a=(2.405*l)/(2*%pi*NA) // a is the maximum core radius in meters
+
+disp(a,"The Maximum Core Radius is less than(in meters)")
+
+//With NA=0.1
+NA=0.1
+a=(2.405*l)/(2*%pi*NA) //a is the maximum core radius in meters
+
+disp(a," The Maximum Core radius increase in m:")
+mprintf(" With a design criterion of V=2 the core diameter would then be 9.5 um")
diff --git a/3710/CH8/EX8.9/Ex8_9.sce b/3710/CH8/EX8.9/Ex8_9.sce
new file mode 100644
index 000000000..bf84d1960
--- /dev/null
+++ b/3710/CH8/EX8.9/Ex8_9.sce
@@ -0,0 +1,12 @@
+//Example 8.9, Page Number 387
+//The Function fpround(dependency) is used to round a floating point number x to n decimal places
+clc;
+
+V=2
+
+//From equation 8.29a
+w=0.65+(1.619*(V**(-1.5)))+(2.879*(V**-6))
+w=fpround(w,3)
+mprintf("The Mode Field Irradiance Diameter is:%.3f",w)
+
+//The answer provided in the textbook is wrong