nt = 50000 positionDetecteur = 30100 ### Creation d'un fbg ### fbgStart = 20000 fbgStop = 39660 # Pas du reseaux Longueur_onde_bragg = 1.55005 Longueur_onde_bragg1 = 1.02 Longueur_onde_bragg2 = 1.03 LBRG = linspace(Longueur_onde_bragg1, Longueur_onde_bragg2, fbgStop - fbgStart) GrandLambda = Longueur_onde_bragg / (2 * 1.45) # Difference d'indice deltaIndex = 3.49e-5 fbg1 = PyOFTK.apodizedFBG(2.5, 62.5, 0.05, 0.0, 50e-2, deltaIndex, GrandLambda) gamma = fbg1.nlGamma(1.03) detZero = fbg1.detuning(LMBD) detVec = zeros(nbrPoints_z, double) detVec[fbgStart:fbgStop] += detZero #for i in arange(fbgStop-fbgStart): # detVec[fbgStart+i] = detZero = fbg1.detuning(LBRG[i]) kapZero = fbg1.kappa(LMBD) kapVec = zeros(nbrPoints_z, double) alphaVec = zeros(nbrPoints_z, double) alphaVec[fbgStart:fbgStop] += 0.0 # Initialisation du champ incident d = linspace(-nbrPoints_z / 2, nbrPoints_z / 2, nbrPoints_z) z = linspace(0, fbg1.length, nbrPoints_z)
# # # # Author: Martin Laprise # Universite Laval # [email protected] # # ********************************************* import pylab as pl import PyOFTK LMBD=(1.0526+0.00047)/(2*1.45) LMBD2=0.34485 LNGTH = 100 fbg1 = PyOFTK.apodizedFBG(3.0, 62.5, 0.04, 0.0, 1e-1, 1e-2, LMBD2) fbg2 = PyOFTK.apodizedFBG(3.0, 62.5, 0.04, 0.0, 1e-1, 5e-3, LMBD2) fbg3 = PyOFTK.apodizedFBG(3.0, 62.5, 0.04, 0.0, 1e-1, 1e-3, LMBD2) print "Longueur d'onde de Bragg: " + str(fbg1.braggWavelength) + " um" wvl1 = pl.linspace(fbg1.braggWavelength-0.00047, fbg1.braggWavelength-0.00050, LNGTH) wvl2 = pl.linspace(fbg1.braggWavelength+0.010, fbg1.braggWavelength+0.050, LNGTH) beta2Grating1 = pl.zeros(LNGTH, float) beta2Grating2 = pl.zeros(LNGTH, float) beta2Grating3 = pl.zeros(LNGTH, float) beta3Grating1 = pl.zeros(LNGTH, float) beta3Grating2 = pl.zeros(LNGTH, float) beta3Grating3 = pl.zeros(LNGTH, float) for i in range(LNGTH):
puissanceCrete = 478.8 LMBD = 1.03 nt = 24000 positionDetecteur = 4999 ### Creation d'un fbg ### fbgStart = 2000 fbgStop = 8000 # Pas du reseaux #Longueur_onde_bragg = 1.03238 Longueur_onde_bragg = 1.03105615 #Longueur_onde_bragg = 1.0289438500000001 GrandLambda = Longueur_onde_bragg/(2*1.45) # Difference d'indice deltaIndex = 3.695e-3 fbg1 = PyOFTK.apodizedFBG(2.5, 62.5, 0.05, 0.0, 40e-2, deltaIndex, GrandLambda) #gamma = fbg1.nlGamma(1.03) gamma = 0.0 detZero = fbg1.detuning(LMBD) detVec = zeros(nbrPoints_z,double) detVec[fbgStart:fbgStop] += detZero kapZero = fbg1.kappa(LMBD) kapVec = zeros(nbrPoints_z,double) #iz = arange(1.0,nbrPoints_z+1)/nbrPoints_z #kapVec[3000:nbrPoints_z-1] = kapZero*(iz[3000:nbrPoints_z-1]) kapVec[fbgStart:fbgStop] = kapZero # Initialisation du champ incident d = linspace(-nbrPoints_z/2, nbrPoints_z/2, nbrPoints_z) z = linspace(0,fbg1.length, nbrPoints_z)
# # # # Author: Martin Laprise # Universite Laval # [email protected] # # ********************************************* import pylab as pl import PyOFTK LMBD = (1.0526 + 0.00047) / (2 * 1.45) LMBD2 = 0.34485 LNGTH = 100 fbg1 = PyOFTK.apodizedFBG(3.0, 62.5, 0.04, 0.0, 1e-1, 1e-2, LMBD2) fbg2 = PyOFTK.apodizedFBG(3.0, 62.5, 0.04, 0.0, 1e-1, 5e-3, LMBD2) fbg3 = PyOFTK.apodizedFBG(3.0, 62.5, 0.04, 0.0, 1e-1, 1e-3, LMBD2) print "Longueur d'onde de Bragg: " + str(fbg1.braggWavelength) + " um" wvl1 = pl.linspace(fbg1.braggWavelength - 0.00047, fbg1.braggWavelength - 0.00050, LNGTH) wvl2 = pl.linspace(fbg1.braggWavelength + 0.010, fbg1.braggWavelength + 0.050, LNGTH) beta2Grating1 = pl.zeros(LNGTH, float) beta2Grating2 = pl.zeros(LNGTH, float) beta2Grating3 = pl.zeros(LNGTH, float) beta3Grating1 = pl.zeros(LNGTH, float) beta3Grating2 = pl.zeros(LNGTH, float) beta3Grating3 = pl.zeros(LNGTH, float)