1 files changed, 39 insertions, 0 deletions
diff --git a/3638/CH13/EX13.7/Ex13_7.sce b/3638/CH13/EX13.7/Ex13_7.sce
new file mode 100644
index 000000000..cf6258c2a
--- /dev/null
+++ b/3638/CH13/EX13.7/Ex13_7.sce
@@ -0,0 +1,39 @@
+//Introduction to Fiber Optics by A. Ghatak and K. Thyagarajan, Cambridge, New Delhi, 1999
+//Example 13.7
+//OS=Windows XP sp3
+//Scilab version 5.5.2
+clc;
+clear;
+//given
+Id=1e-9;//Dark current of a silicon PIN photodiode in A
+P=500e-9;//Optical power in W
+R=0.65;//Responsivity in A/W
+Rl=1000;//Value of load resistor in Ohms
+e=1.6e-19//Electronic charge in C
+kB=1.38e-23;//Boltzmann constant in SI Units
+Deltaf=100e6;//Detector bandwidth in Hz
+T=300;//Missing data- Temperature in K
+M=50;//Internal gain corresponding to input optical power P
+x=0;//No excess noise
+
+I=M*R*P;//Signal current in A
+mprintf("\n I=%.2f uA",I/1e-6)//Division by 10^(-6) to convert into uA
+//Let the root mean square shot noise current be Ins
+//The rms shot noise current due to signal is:
+Ins=sqrt(2*e*M*I*Deltaf);//As the root mean square shot noise current is the square root of mean square shot noise current in A
+mprintf("\n The rms shot noise current due to signal = %.2f nA",Ins/1e-9);//Division by 10^(-9) to convert into nA
+//The answers vary due to round off error
+
+//The rms shot noise current due to dark current is:
+Ins=sqrt(2*e*(M^2)*Id*Deltaf);//As the root mean square shot noise current is the square root of mean square shot noise current in A
+mprintf("\n The rms shot noise current due to dark current = %.2f nA",Ins/1e-9);//Division by 10^(-9) to convert into nA
+//The answers vary due to round off error
+
+//The rms shot thermal noise current is:
+Ins=sqrt(4*kB*T*Deltaf/Rl);//As the root mean square shot noise current is the square root of mean square shot noise current in A
+mprintf("\n The rms shot thermal noise current = %.2f nA",Ins/1e-9);//Division by 10^(-9) to convert into nA
+//The answers vary due to round off error
+SNR=((M*R*P)^2)/(2*e*(M^(2+x))*(R*P+Id)*Deltaf+4*kB*T*Deltaf/Rl);//Corresponding Signal-to-noise ratio since x=0
+mprintf("\n SNR = %f",SNR);//The answers vary due to round off error
+mprintf("\n SNR in dB = %.2f dB",10*log10(SNR));//For conversion to dB
+//The answers vary due to round off error