diff options
Diffstat (limited to 'Optical_fiber_communication_by_gerd_keiser/chapter11_2.ipynb')
| -rwxr-xr-x | Optical_fiber_communication_by_gerd_keiser/chapter11_2.ipynb | 416 |
1 files changed, 416 insertions, 0 deletions
diff --git a/Optical_fiber_communication_by_gerd_keiser/chapter11_2.ipynb b/Optical_fiber_communication_by_gerd_keiser/chapter11_2.ipynb new file mode 100755 index 00000000..c71b9640 --- /dev/null +++ b/Optical_fiber_communication_by_gerd_keiser/chapter11_2.ipynb @@ -0,0 +1,416 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:4930bebda7035fefc5221e84cdee44dfb2c772873b238e2ca2714c2de4369a09"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 11: Optical amplifires"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.1, Page Number: 397"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declaration\n",
+ "Vg = 2*10**8 #group velocity (m/s)\n",
+ "h = 6.625*10**-34 #planks constant (J*s)\n",
+ "C = 3*10**8 #free space velocity (m/s)\n",
+ "Lam_bda = 1550*10**-9 #operating wave length(nm)\n",
+ "V = C/Lam_bda #frequency (Hz)\n",
+ "w = 5*10**-6 #width of optical amplifier (meters)\n",
+ "d = 0.5*10**-6 #thickness of optical amplifier (meter)\n",
+ "Ps = 10**-6 #optical signal of power\n",
+ "\n",
+ "#calculation\n",
+ "Nph = Ps/(Vg*h*V*w*d) #photon density\n",
+ "\n",
+ "#result\n",
+ "print \"The photon density Nph = \" ,round(Nph*1e-16,2),\"e+6 photons/m3\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The photon density Nph = 1.56 e+6 photons/m3\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.2, Page Number: 397"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 11.2(a)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#varible declaration\n",
+ "I = 100.0*10**-3 #bias current (Amps)\n",
+ "w = 3.0*10**-6 #active area width (meters)\n",
+ "L = 500.0*10**-6 #amplifier lenght (meters)\n",
+ "d = 0.3*10**-6 #active area thick ness(meters)\n",
+ "q = 1.6*10**-19 #charge (coulombs)\n",
+ "Tuo = 0.3 #The confinement factor\n",
+ "a = 2*10**-20 #gain coefficient (square meter)\n",
+ "J = I/(w*L) #3bias current density (Amp/squre meter)\n",
+ "nth = 10**24 #threshold density (per cubic meter)\n",
+ "Tuor = 10**-9; #Time constant (seconds)\n",
+ "\n",
+ "\n",
+ "#calculation\n",
+ "Rp = I/(q*d*w*L) # The pumping rate((electron/m3)/s)\n",
+ "\n",
+ "#result\n",
+ "print \"The pumping rate Rp = \" , round(Rp*1e-33,2)*10**33,\" (electron/m3)/s\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The pumping rate Rp = 1.39e+33 (electron/m3)/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 11.2(b)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#varible declaration\n",
+ "I = 100.0*10**-3 #bias current (Amps)\n",
+ "w = 3.0*10**-6 #active area width (meters)\n",
+ "L = 500.0*10**-6 #amplifier lenght (meters)\n",
+ "d = 0.3*10**-6 #active area thick ness(meters)\n",
+ "q = 1.6*10**-19 #charge (coulombs)\n",
+ "Tuo = 0.3 #The confinement factor\n",
+ "a = 2*10**-20 #gain coefficient (square meter)\n",
+ "J = I/(w*L) #3bias current density (Amp/squre meter)\n",
+ "nth = 10**24 #threshold density (per cubic meter)\n",
+ "Tuor = 10**-9; #Time constant (seconds)\n",
+ "\n",
+ "\n",
+ "#calculation\n",
+ "Rp=I/(q*d*w*L) #The pumping rate((electron/m3)/s)\n",
+ "g0 = Tuo*a*Tuor*(round(Rp*1e-33,2)*10**33-(nth/Tuor)) #The zero singal(1/cm)\n",
+ "\n",
+ "#result\n",
+ "print \"The zero singal g0 = \" ,round(g0),\"1/m =\", round(g0/100,1),\"1/cm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The zero singal g0 = 2340.0 1/m = 23.4 1/cm\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.3, Page Number: 404"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declaration\n",
+ "Lambda_p = 980.0*10**-9 #pump wavelength(nm)\n",
+ "Lambda_s = 1550.0*10**-9 #signal wavelength(nm)\n",
+ "Pp_in = 30.0*10**-3 #input pump power (watts)\n",
+ "G = 1.0*10**2 #gain\n",
+ "\n",
+ "#calculation\n",
+ "Ps_in = (Lambda_p/Lambda_s)*Pp_in/(G-1) #maximum input power(W)\n",
+ "Ps_out = Ps_in+(Lambda_p/Lambda_s)*Pp_in #maximum output power(W)\n",
+ "Ps_out_db = 10*(math.log10(Ps_out*10**3)) #maximum output power(dBm)\n",
+ "\n",
+ "#result\n",
+ "print \"The maximum input power = \" , round(Ps_in*10**6) , \"uW\"\n",
+ "print \"The maximum output power = \" , round(Ps_out*10**3,1),\"mW\"\n",
+ "print \"The maximum output power = \" , round(Ps_out_db,1),\"dBm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The maximum input power = 192.0 uW\n",
+ "The maximum output power = 19.2 mW\n",
+ "The maximum output power = 12.8 dBm\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.6, Page Number: 412"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declarion\n",
+ "Q = 6 #Q factor of 6\n",
+ "\n",
+ "#calculation\n",
+ "OSNR = 0.5*Q*(Q+math.sqrt(2))\n",
+ "OSNR_DB = 10*(math.log10(OSNR)) #The optical signal to noise ratio(dB)\n",
+ "\n",
+ "#result\n",
+ "print \"The optical signal to noise ratio (OSNR) = \" ,round(OSNR,2),\"=\", round(OSNR_DB,1),\"dB\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The optical signal to noise ratio (OSNR) = 22.24 = 13.5 dB\n"
+ ]
+ }
+ ],
+ "prompt_number": 7
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.7, Page Number: 413"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declaration\n",
+ "Lambda_p = 980*10**-9 #pump wavelength (meters)\n",
+ "Lambda_s = 1540*10**-9 #signal wavelength (meters)\n",
+ "Ps_out = 10*10**-3 #output signal power(mW)\n",
+ "Ps_in = 10**-3 #input signal power(mW)\n",
+ "\n",
+ "#calculation\n",
+ "Pp_in = (Lambda_s/Lambda_p)*(Ps_out-Ps_in) #pump power at input(mW)\n",
+ "\n",
+ "#result\n",
+ "print \"Pump power = \" , round(Pp_in*10**3) ,\"mW\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Pump power = 14.0 mW\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.8, Page Number: 413"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declaration\n",
+ "P_ASE1 = -22 #ASE level (dBm)\n",
+ "P_ASE2 = -16 #ASE level (dBm)\n",
+ "Pout = 6 #amplified signal level (dBm)\n",
+ "\n",
+ "#calculation\n",
+ "OSNR1 = Pout-P_ASE1 \n",
+ "OSNR2 = Pout-P_ASE2 #The optical signal to noise ratio(dBm)\n",
+ "\n",
+ "#result\n",
+ "print \"Optical SNR OSNR1 = \" , round(OSNR1) , \"dBm\"\n",
+ "print \"Optical SNR OSNR2 = \" , round(OSNR2) , \"dBm\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Optical SNR OSNR1 = 28.0 dBm\n",
+ "Optical SNR OSNR2 = 22.0 dBm\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.9, Page Number: 414"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declaration\n",
+ "G1 = 10**(30/10) #gain(dB)\n",
+ "G2 = 10**(20/10)\n",
+ "\n",
+ "#calculation\n",
+ "Fpath1 = (((G1-1)/math.log(G1))**2)/G1 #noise penalty factor for G1\n",
+ "fpath_db1=10*math.log10(Fpath1) #noise penalty factor(dB)\n",
+ "Fpath2 = (((G2-1)/math.log(G2))**2)/G2 #noise penalty factor for G2\n",
+ "fpath_db2=10*math.log10(Fpath2) #noise penalty factor(dB)\n",
+ "\n",
+ "#result\n",
+ "print \"Noise penalty factor for G1 = \",round(fpath_db1,1),\"dB\"\n",
+ "print \"Noise penalty factor for G2 = \",round(fpath_db2,1),\"dB\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Noise penalty factor for G1 = 13.2 dB\n",
+ "Noise penalty factor for G2 = 6.6 dB\n"
+ ]
+ }
+ ],
+ "prompt_number": 10
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 11.10, Page Number: 415"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "\n",
+ "#variable declaration\n",
+ "etta = 0.65 #quantum efficiency\n",
+ "nsp = 2 #population inversion\n",
+ "R =50 #load resistance(ohms)\n",
+ "Lambda = 1550*10**-9 #oprating wavelength(meters)\n",
+ "T = 300 #room temperature(kelvin)\n",
+ "h = 6.625*10**-34 #planks constant(J*s)\n",
+ "C = 3*10**8 #free space velocity(m/s)\n",
+ "kB = 1.38*10**-23 #boltzmann's constant \n",
+ "V = C/ Lambda #(Hz)\n",
+ "q = 1.6*10**-19 #Charge (coulombs)\n",
+ "\n",
+ "#calculation\n",
+ "Ps_in = kB*T*h*V/(R*nsp*(etta**2)*(q**2)) #maxiamum input optical power level(Watt)\n",
+ "\n",
+ "#result\n",
+ "print \"Upper bound input otical power level <\",round(Ps_in*10**6,1),\"uW\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Upper bound input otical power level < 490.8 uW\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |