path: root/Optical_fiber_communication_by_gerd_keiser/chapter11_1.ipynb

{
 "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": {}
  }
 ]
}