Optimisation des Réseaux de Bragg en WDM
Optimisation des Réseaux de Bragg en WDM
Présentée par :
Mr KARIM Fethallah
Soutenue le 06/12/2011 devant le jury :
Mes premières pensées vont tout particulièrement à mes parents à qui je dois ce succès
et pour qui le sens du sacrifice au cours de ces longues années d’étude ne s’est jamais tari et a
été le plus précieux de leurs cadeaux. Grâce à leur soutien et à leurs encouragements
sublimement dosés, ils m’ont permis de vivre mes années d’études avec sérénité et efficacité.
Grand merci.
J’exprime ensuite une profonde gratitude, et ces mots sont pesés, à mon directeur de
M. KANDOUCI, Professeur à l’université de Bel Abbes, membres de mon jury pour m’avoir
Télécommunications, ainsi que tous mes amis pour leurs encouragements qui m’ont permis
décennies dans le secteur des télécommunications optiques, ils représentent une alternative
attirante et moins coûteuse grâce à leur simple structure. Nous retrouvons ces composants
dans plusieurs applications comme le filtrage multi canal adaptatif, les multiplexeurs à
chromatique, etc.
de réseaux de Bragg par des métaheuristiques hybrides. Une comparaison entre les
simplexe (AGCS) et le recuit simulé avec recherche heuristique de motif (RSRHM) a été
faite. La RTD a été choisie pour reconstruire les paramètres physiques d’un réseau mono
canal et un réseau multi canal. Nous avons appliqué cet algorithme hybride pour estimer les
coefficients thermo-optique et expansion thermique d’un réseau mono canal soumis à des
gradients de température. La RTD a été aussi utilisée pour reconstruire des profils de tractions
d’un réseau mono canal et un réseau échantillonné soumis à des compressions et des
étirements.
communications area, they represent an attractive and less expensive alternative thanks to
their simple structure. We find these components in several applications like adaptive multi-
compensation, etc.
Search (DTS), Simulated Annealing Heuristic Pattern Search (SAHPS) and Simplex Coding
Genetic Algorithm (SCGA) performances has been made. We have chosen DTS method for
We have applied this hybrid method for estimation of thermal-expansion and thermo-optic
DTS is also used for reconstruction of a single channel grating and a sampled grating strain
profiles, whiles these gratings are subjected to positive and negative strains.
I Introduction ........................................................................................................................ 5
II Principe de base du MDL................................................................................................... 5
III Présentation des différentes technologies des multiplexeurs optiques .............................. 9
III.1 Introduction ................................................................................................................ 9
III.2 Caractéristiques des composants optiques ................................................................ 10
III.2.1 Pertes d’insertion............................................................................................ 10
III.2.2 Dépendance de la polarisation ....................................................................... 10
III.2.3 Diaphonie et isolation .................................................................................... 10
III.2.4 Dispersion ...................................................................................................... 11
IV Présentation des différents types de dé/multiplexeurs appliqués dans un système de
transmission MDL .......................................................................................................... 12
IV.1 Composants MDL passifs de base ............................................................................ 13
IV.2 Filtres diélectriques................................................................................................... 14
IV.2.1 Introduction................................................................................................... 14
IV.2.2 Théorie des filtres interférentiels .................................................................. 14
IV.2.2.1 Introduction.................................................................................... 14
IV.2.2.2 Filtres multi couches anti réflexion (AR) ...................................... 14
IV.2.2.3 Filtres multi couches avec haute réflexion (HR) ........................... 15
IV.2.2.4 Filtre passe bande........................................................................... 16
IV.2.3 Coupleurs à trois ports basés sur des lentilles................................................ 17
IV.2.4 Composants en zigzag.................................................................................... 18
IV.2.5 Caractéristiques de quelques dé/multiplexeurs à base de filtres à couches
minces............................................................................................................. 19
IV.2.5.1 Coupleurs mono canal utilisés dans le MDL ................................. 19
IV.2.5.2 Coupleurs à division de bande ....................................................... 19
IV.3 Composants MDL basés sur des réseaux de diffraction ........................................... 22
IV.3.1 Introduction.................................................................................................... 22
IV.3.2 Principe de fonctionnement des réseaux basés sur la diffraction .................. 23
IV.3.2.1 Ordres de diffraction ....................................................................... 23
IV.3.2.2 Intervalle spectral libre ................................................................... 23
IV.3.3 Exemple d’un démultiplexeur basé sur un réseau de diffraction ................... 24
IV.4 Composants MDL basés sur des réseaux de guides d’onde assemblés .................... 24
IV.5 Composants basés sur des interféromètres de Mach-Zehnder ................................. 26
IV.6 Filtres optiques accordables ...................................................................................... 27
IV.6.1 Introduction.................................................................................................... 27
IV.6.2 Coupleurs directionnels 2×2 accordables ...................................................... 28
IV.6.3 Filtres multi réseaux accordables................................................................... 29
V Conclusion ....................................................................................................................... 30
I Introduction ...................................................................................................................... 33
II Principes fondamentaux des réseaux de Bragg................................................................ 33
II.1 Introduction .............................................................................................................. 33
II.2 Théorie des modes couplés......................................................................................... 36
II.3 Applications des réseaux de Bragg............................................................................. 38
II.4 Modélisation des réseaux de Bragg............................................................................ 38
II.5 Réseau de Bragg uniforme ......................................................................................... 40
II.6 Méthode des matrices de transfert .............................................................................. 41
II.7 Réseaux de Bragg à pas variable ................................................................................ 42
II.7.1 Principe de base d’un réseau à pas variable ..................................................... 42
II.7.2 Résultats de simulation de la réponse spectrale des réseaux de Bragg à pas
variable linéaires............................................................................................... 45
II.7.2.1 Réseaux chirpés linéaires avec différentes variables chirp ................. 45
II.7.2.2 Réseaux chirpés linéaires avec différentes longueurs ......................... 46
II.7.2.3 Réseaux chirpés linéaires avec différents changements d’indice de
réfraction ............................................................................................... 46
II.7.2.4 Relation entre la variable chirp, le paramètre chirp et la longueur
d’onde centrale ................................................................................... 48
II.8 Réseaux de Bragg apodisés ........................................................................................ 48
II.8.1 Principe de base des réseaux apodisés............................................................... 48
II.8.2 Fonctions d’apodisation..................................................................................... 50
II.9 Réseaux de Bragg échantillonnés............................................................................... 53
II.9.1 Principe de base................................................................................................. 53
II.9.2 Techniques de densification des canaux............................................................ 55
II.9.2.1 Changement de la période d’échantillonnage....................................... 55
II.9.2.2 Technique de décalage périodique de phase......................................... 59
II.9.2.3 Technique du changement linéaire de phase ........................................ 62
II.10 Techniques de décalage des canaux ......................................................................... 63
II.10.1 Ajout d’une extra phase.................................................................................. 63
II.10.2 Déplacement des canaux par application de gradients de température et des
tractions mécaniques sur un réseau de Bragg................................................. 65
II.10.2.1 Application du gradient de température ............................................ 65
II.10.2.2 Application d’une traction mécanique............................................... 69
II.11 Multiplexeurs et démultiplexeurs à base de réseaux de Bragg................................. 74
II.11.1 Introduction .................................................................................................... 74
II.11.2 Configuration d’un démultiplexeur optique ................................................... 75
II.11.3 Configuration d’un multiplexeur à insertion/extraction en longueurs d’onde 78
III Conclusion ....................................................................................................................... 81
I Introduction ........................................................................................................................ 84
II Algorithmes métaheuristiques........................................................................................... 85
II.1 Algorithme génétique ................................................................................................. 86
II.1.1 Introduction ...................................................................................................... 86
II.1.2 Codage des individus d’une population ........................................................... 87
II.1.2.1 Codage binaire..................................................................................... 87
II.1.2.2 Codage à caractères multiples ............................................................. 87
II.1.3 Evaluation et sélection...................................................................................... 87
II.1.3.1 Sélection par roulette........................................................................... 88
II.1.3.2 Sélection par rang ................................................................................ 89
II.1.3.3 Sélection par tournoi............................................................................ 90
II.1.3.4 Elitisme................................................................................................ 90
II.1.4 Croisement........................................................................................................ 90
II.1.4.1 Croisement en un point........................................................................ 90
II.1.4.2 Croisement en deux points .................................................................. 91
II.1.5 Mutation ........................................................................................................... 91
II.2 Recuit simulé .............................................................................................................. 92
II.2.1 Introduction ...................................................................................................... 92
II.2.2 Notions ............................................................................................................. 93
II.2.2.1 Probabilité de Boltzmann .................................................................... 93
II.2.2.2 Critère de Metropolis........................................................................... 93
II.2.2.3 Chaine de Markov et équilibre thermique ........................................... 94
II.2.2.4 Algorithme........................................................................................... 94
II.2.3 Paramètres ........................................................................................................ 95
II.2.3.1 Température initiale............................................................................. 95
II.2.3.2 Décroissement de température ............................................................ 95
II.2.3.3 Nombre d’itérations à température constante...................................... 95
II.2.3.4 Critères d’arrêt..................................................................................... 95
II.2.4 Recuit simulé à pas tabulé ................................................................................ 96
II.2.4.1 Prédétermination du vecteur pas ......................................................... 96
II.2.4.2 Etapes de l’algorithme RSPT .............................................................. 96
II.3 Recherche tabou ......................................................................................................... 97
II.3.1 Introduction ...................................................................................................... 97
II.3.2 Principe de l’algorithme ................................................................................... 97
II.3.3 La recherche tabou de Hu................................................................................. 98
III Métaheuristiques hybrides ............................................................................................... 99
III.1 Introduction .............................................................................................................. 99
III.2 Classification hiérarchique........................................................................................ 99
III.3 Classification générale ............................................................................................ 100
III.4 Algorithme génétique avec codage de simplexe (AGCS)....................................... 101
III.4.1 Introduction ................................................................................................... 101
III.4.2 Méthode de Nelder-Mead ............................................................................. 101
III.4.3 Description générale des étapes de l’AGCS ................................................. 103
III.4.3.1 Initialisation...................................................................................... 103
III.4.3.2 Boucle de l’algorithme génétique .................................................... 103
III.[Link] Sélection.......................................................................... 103
III.[Link] Croisement et mutation ................................................... 104
III.[Link] Réduction de la population.............................................. 105
III.[Link] Intensification des membres par l’algorithme de Nelder-
Mead ................................................................................ 105
III.5 Recuit simulé avec recherche heuristique de motif (RSRHM) .............................. 105
III.5.1 Introduction ................................................................................................... 105
III.5.2 Direction de descente approximative (DDA)................................................ 105
III.5.3 Recherche heuristique de motif (RHM) ........................................................ 107
III.5.4 Description générale des étapes du RSRHM ................................................ 108
III.6 Recherche tabou dirigée (RTD) .............................................................................. 109
III.6.1 Les éléments mémoire de la recherche tabou (RTa) ..................................... 110
III.6.1.1 introduction ...................................................................................... 110
III.6.1.2 Liste tabou multi-classée (LT) ......................................................... 110
III.6.1.3 Liste des régions visitées (LRV) ...................................................... 112
III.6.2 Stratégies de recherche locale et de voisinage .............................................. 113
III.6.3 Description générale des étapes de la RTD................................................... 113
III.6.3.1 Introduction ...................................................................................... 113
III.6.3.2 La boucle exploration-diversification .............................................. 114
III.6.3.3 Recherche d’intensification.............................................................. 115
IV Synthèse et optimisation des réseaux de Bragg par la RTD.......................................... 116
IV.1 Introduction ............................................................................................................ 116
IV.2 Comparaison entre les performances des algorithmes RTD, RSRHM et AGCS ... 117
IV.3 Synthèse des paramètres physiques d’un réseau de Bragg mono canal par la RTD118
IV.4 Synthèse des paramètres physiques d’un réseau de Bragg multi canal par la RTD 121
IV.5 Estimation des coefficients thermo-optique et expansion thermique d’un réseau de
Bragg par la RTD................................................................................................... 124
IV.6 Optimisation des différents profils des contraintes mécaniques appliquées sur un
réseau de Bragg par la RTD................................................................................... 129
IV.6.1 Introduction................................................................................................... 129
IV.6.2 Reconstruction d’un profil d’une traction appliquée sur un réseau de Bragg à
partir d’un spectre original déformé............................................................. 132
V Conclusion ...................................................................................................................... 134
Conclusion générale............................................................................................................ 137
Annexe 1............................................................................................................................. 140
Liste des figures et tableaux ............................................................................................... 145
Liste des symboles.............................................................................................................. 150
Bibliographie ...................................................................................................................... 155
INTRODUCTION GENERALE
Introduction générale
2
polarisation et leur compatibilité est parfaite avec les différentes fibres mono mode de
transmission. Ces réseaux sont utilisés dans plusieurs applications comme le filtrage mono
canal et multi canal, ce sont des composants fondamentaux de quelques types de
multiplexeurs à insertion/extraction en longueur d’onde, ils sont utilisés comme des
compensateurs accordables de la dispersion chromatique dans un système MDL, ces
dispositifs sont appliqués aussi dans le routage optique.
Le réseau de Bragg est une variation périodique de son indice de réfraction à travers
une direction de propagation dans le cœur de la fibre. Pour la caractérisation de ce réseau, la
théorie des modes couplés est utilisée pour analyser la nature de propagation de la lumière
dans la fibre perturbée (modes en co propagation et modes en contre propagation). La
méthode des matrices de transfert est appliquée pour résoudre les équations des modes
couplés. Nous avons conçu alors un programme sous l’environnement MATLAB basé sur
cette méthode pour représenter les caractéristiques spectrales d’un réseau de Bragg uniforme,
à pas variable, apodisé et échantillonné.
Le but de notre travail est de reconstruire les paramètres physiques d’un réseau de
Bragg mono canal et multi canal en utilisant des algorithmes hybrides métaheuristiques, et de
comparer entre les performances de ces algorithmes en termes de fonction objectif, nombre
d’évaluations de cette fonction et le temps de calcul. La Recherche Tabou Dirigée (RTD),
l’Algorithme Génétique avec Codage de Simplexe (AGCS) et le Recuit Simulé avec
Recherche Heuristique de Motif (RSRHM) seront utilisés dans cette thèse pour synthétiser un
réseau de Bragg.
Notre thèse englobe trois grandes parties. Dans le chapitre I nous évoquerons les
différents filtres optiques largement utilisés pour concevoir les multiplexeurs et les
démultiplexeurs optiques. Dans le chapitre II, une analyse des différents types de réseaux de
Bragg sera faite en représentant les caractéristiques spectrales de ce réseau en utilisant la
méthode des matrices de transfert. Quelques techniques de densification des canaux dans un
réseau de Bragg échantillonné seront présentées. Ces canaux peuvent être déplacés vers
d’autres longueurs d’onde choisies en appliquant des gradients de température ou des
tractions mécaniques sur la fibre de Bragg. Enfin, le chapitre III sera consacré à la synthèse
d’un réseau de Bragg mono canal et échantillonné. La Recherche Tabou Dirigée sera utilisée
pour reconstruire les paramètres physiques de ces réseaux, cette méthode sera appliquée pour
estimer le coefficient thermo-optique et le coefficient d’expansion thermique d’un réseau de
Bragg mono canal dans le cas où nous choisissons de régler ce réseau sur une certaine
longueur d’onde en lui appliquant un changement de température. La même métaheuristique
hybride sera appliquée pour reconstruire des profils de tractions positives et négatives d’un
réseau mono canal et échantillonné à partir de spectres déformés originaux, dans le but de
déplacer les spectres réfléchis vers d’autres longueurs d’onde prédéfinies.
3
Chapitre I
5
Tableau I.1 Valeurs couramment utilisées pour OC-N et STM-M [I.7]
Fig (1.1). Evolution de la capacité de transmission d’un lien optique: (a) lien optique point à
point simple; (b) lien multiplexé spatialement en quatre canaux; (c) Lien multiplexé en quatre
longueurs d’onde [I.7]
6
Ceci montre que pour maintenir une compatibilité entre SDH et SONET, en pratique N est un
multiple de 3. Au contraire du SONET, la SDH ne fait pas de distinction entre un signal
électrique logique (STS-N dans SONET) et un signal optique physique (OC-N), alors ces
signaux sont tous les deux reconnus par la STM-M comme un seul type de signal [I.7]. Le
tableau (I.1) récapitule les valeurs couramment utilisées pour OC-N et STM-M.
Dans les liens optiques point à point standards, une seule fibre possède une source de
lumière sur son extrémité d’émission et un photo détecteur sur l’autre extrémité de réception,
comme il est montré sur la fig (1.1)(a) [I.8]. Si nous considérons que ce lien optique a une
capacité de transmission de 2.5 Gbit/s, ça veut dire qu’on a multiplexé électriquement 16
STM-1 SDH ou 16 OC-3 SONET avec un débit élémentaire de 155 Mbit/s. Ce signal
électrique obtenu est utilisé pour moduler une diode laser afin de convertir le signal électrique
en un format optique et il sera par la suite émis sur la fibre optique. Dans ce type de liens
optiques, les signaux qui arrivent des différentes sources de lumière utilisent séparément et
uniquement des fibres optiques dédiées. Pour augmenter le nombre de canaux OC-3 à 64 sans
passer par un Multiplexage de Division en Longueurs d’onde (MDL), nous aurons besoin de
quatre fibres monomode séparées, comme il est illustré dans la fig (1.1) (b). Ce schéma de
transmission est traditionnellement appelé Multiplexage de Division en Espace (MDE).
Cependant, dans le MDL nous pouvons utiliser des longueurs d’onde espacées et
indépendantes pour émettre les 64 canaux OC-3/STM-1 sur la même fibre optique, comme il
est présenté sur la fig (1.1) (c) [I.7].
Pour voir le potentiel d’un MDL, examinons les caractéristiques d’une source optique
de haute qualité. Par exemple, un laser à répercussion distribué possède un spectre fréquentiel
très étroit sur un ordre de 1 MHz, qui est équivalent à une largeur spectrale de 10-5 nm [I.7].
Lorsqu’on utilise ce type de sources, un espacement en longueur d’onde variant de 0.4 nm à
1.6 nm est typiquement employé. Cet intervalle en longueur d’onde est utilisé pour prendre en
considération un éventuel déplacement des longueurs d’onde centrales à cause d’une dérive
de température ou d’une traction appliquée. La figure (1.2) présente l’atténuation de la
lumière dans une fibre à base de silice, en fonction de la longueur d’onde. Ce graphe montre
que les deux régions à faibles pertes d’une fibre monomode appartiennent aux deux
intervalles : 1270-1350 nm (la fenêtre 1310 nm) et 1480-1600 nm (la fenêtre 1550 nm) [I.7].
Pour trouver la largeur de bande optique correspondante à une largeur spectrale
particulière concernant ces deux régions, nous utilisons la relation fondamentale c=λυ qui
relie une longueur d’onde λ à une fréquence porteuse υ, où c est la vitesse de la lumière. En
faisant une dérivation de la relation fondamentale, nous obtenons [I.7]
c
(1.1)
2
où la déviation en fréquence correspond à une déviation en longueur d’onde autour
de . A partir de l’équation (1.1), la largeur de bande optique = 14 THz correspond à une
bande spectrale utilisable =80 nm dans la fenêtre 1310 nm. De la même manière, la
7
Fig (1.2). Atténuation de la lumière dans une fibre de silice en fonction de la longueur d’onde
[I.7]
largeur de bande optique = 15 THz correspond à une bande spectrale utilisable =120
nm dans la fenêtre 1550 nm. Ceci implique que la largeur de bande totale disponible d’une
fibre est approximativement égale à 30 THz dans les deux fenêtres de régions à faibles pertes.
Malgré que la largeur spectrale d’une source de haute qualité occupe seulement une
largeur de bande optique étroite, les deux fenêtres à faible pertes fournissent plusieurs régions
additionnelles. En utilisant un certain nombre de sources de lumière, chacune émettant sur une
longueur d’onde différente qui est suffisamment espacée de sa voisine de telle sorte pour
éviter l’interférence. Par exemple, si nous prenons une bande spectrale de 0.8 nm (ou une
bande fréquentielle de 100 GHz), nous pouvons alors émettre 50 signaux indépendants dans
un intervalle de bande allant de 1530 nm jusqu’à 1560 nm sur une seule fibre. Ceci est le
principe de base d’un multiplexage de division en longueurs d’onde (MDL) [I.7].
Le MDL est essentiellement un multiplexage de division en fréquences porteuses
optiques, les standards MDL développés par l’ITU-T spécifient les espacements entre canaux
par les termes de fréquences [I.9, I.10]. La recommandation G.692 spécifie la sélection des
canaux à partir d’une grille de fréquences référencée à 193.100 THz (1552.524 nm dans le
verre) et un espacement de 100 GHz (0.8 nm à 1552 nm) entre ces canaux. La
recommandation inclut aussi d’autres espacements comme 50 GHz (0.4 nm) et 200 GHz (1.6
nm)[I.7].
Malgré que l’installation d’un câble de fibre optique soit très coûteuse, l’extension de
la capacité d’un réseau optique rend cette installation plus économique. Traditionnellement, la
8
capacité d’un lien optique est améliorée en augmentant les débits de transmission.
Initialement, les travaux expérimentaux ont permis d’atteindre des débits de transmission de
2.5 Gbit/s (SONET OC-48 ou SDH STM-16). Les systèmes de transmission qui fonctionnent
à ces débits ont été commercialisés dés l’année 1991. Cependant, en allant au niveau de
multiplexage de 10 Gbit/s (OC-192), beaucoup d’obstacles apparaissent et peuvent
sérieusement dégrader la performance d’un réseau MDL, nous pouvons citer quelques uns
[I.7]
- La dispersion chromatique d’une fibre de transmission, qui limite le débit de transmission en
créant un étalement spectral des impulsions optiques transmises.
- Le gain non uniforme correspondant à la gamme des longueurs d’onde désirées dans des
amplificateurs à fibres dopées erbium;
- Le phénomène d’éparpillement de Raman stimulé (ERS), qui est défini par des intéractions
entre les signaux optiques et les vibrations moléculaires ou acoustiques dans la fibre ;
- Les effets non linéaires de la fibre qui apparaissent à partir de la modulation de l’indice de
réfraction de la silice par l’intermédiaire des changements d’intensité dans le signal, ce qui
induit une production des effets comme : le mixage à quatre ondes, la modulation auto phase
et la modulation de phase croisée;
- La dispersion de mode de polarisation qui apparait à partir d’un déplacement des modes de
polarisation orthogonale à des différentes vitesses, ceci est dû à la biréfringence de la fibre;
- Des réflexions à partir d’épissures et de connecteurs qui peuvent causer une instabilité des
sources laser.
De nouvelles fibres conçues, des techniques de compensation de dispersion et des
isolateurs optiques ont été introduits pour remédier à ces limitations citées au dessus.
9
Fig (1.3). Définition de quelques paramètres caractérisant un composant optique
assemblés, ont été introduites et appliquées [I.14]. Dans la section suivante, quelques
caractéristiques des filtres optiques seront présentées.
10
Fig (1.4). Réponses spectrales de deux canaux correspondants à un filtre diélectrique à
couches minces (Thin Film Filter) et un réseau de Bragg (Fiber Grating) [I.14]
III.2.4 Dispersion
Il existe deux types de dispersion qui sont concernées dans les systèmes MDL ; la
dispersion de mode de polarisation (DMP) et la dispersion chromatique (DCR). La DMP est
la variation du délai de groupe en fonction de l’état de polarisation. La dispersion chromatique
est la variation du délai de groupe en fonction de la longueur d’onde. Dans un démultiplexeur,
le délai de groupe est constant dans la bande passante, mais il varie significativement aux
deux extrémités de la largeur de bande du spectre (voir fig (1.5)). Cependant, un réseau de
Bragg peut avoir quelques rides dans sa bande passante comme il est montré sur la fig (1.6).
Fig (1.5). Délai de groupe d’un coupleur interférométrique fonctionnant à 100 GHz [I.14]
11
Fig (1.6). Dispersion d’un réseau de Bragg (FBG) et d’un filtre diélectrique à couches minces
(TFF) [I.14]
12
Fig (1.7). Implémentation d’un lien optique multiplexé en longueurs d’onde [I.16]
Fig (1.9). Spectres de réflectivité d’un multi couches conçu avec 5, 10 et 15 périodes, centrés
sur la longueur d’onde 1550 nm, avec nB=1.45 et nE=2.1 [I.17]
15
La région de haute réflexion est appelée bande d’arrêt du multi couches HR. La
largeur de bande d’arrêt, 2Δg, ne dépend pas du nombre de périodes, mais de rapport des
indices hauts et bas de réfraction.
Les deux extrémités de bandes de stop sont à gedge=1±Δg avec Δg qui est décrite comme suit
[I.17]
2 1 n E n B
g sin (1.3)
nE nB
La bande passante devient large pour un rapport d’indice de réfraction élevé. Dans
l’exemple présenté sur la fig (1.9), nB est égal à 1.45 alors que nE est égal à 2.1. A partir de
cette figure (spectre conçu à 15 périodes), nous remarquons que l’intervalle de bande d’arrêt
est 1387-1756 nm pour un filtre centré sur 1550 nm.
Le niveau de pique de réflectivité dépend du rapport des indices de réfraction et du
nombre de périodes. Un rapport élevé et un grand nombre de périodes rendent le niveau de
pique de réflectivité plus haut dans la bande d’arrêt [I.18, p. 186].
2p
n ns
R 1 4 B (1.4)
nE nA
16
Fig (1.10). Réponse spectrale d’un filtre passe bande avec différents nombres de cavités [I.17]
c
ISL (1.6)
2nL cos
Fig (1.11). Schéma d’un filtre à couches minces avec trois ports [I.17]
17
La lumière qui émerge de la fibre commune est collimatée par la première lentille et
interagit après avec le filtre interférentiel. La lumière transmise par le filtre est couplée dans la
fibre T par une autre lentille, alors que la lumière réfléchie par le filtre est couplée dans la
fibre R. Comme il est montré sur la fig (1.11), le dispositif permet de séparer deux bandes
spectrales entrantes, il fonctionne alors comme un démultiplexeur. Nous pouvons aussi
utiliser ce dispositif avec une configuration inversée de telle sorte qu’il peut combiner deux
bandes spectrales. Ce composant peut fonctionner aussi comme un multiplexeur [I.17].
Un collimateur réduit la divergence des faisceaux lumineux et rend ces faisceaux sous
forme d’ondes planes. Ce dispositif réduit l’angle d’incidence nécessaire pour séparer les
faisceaux dans un composant multi ports [I.28]. De petits angles d’incidence permettent de
réduire la sensibilité de polarisation d’un filtre. Dans des configurations typiques on utilise
seulement des angles de petits degrés [I.29, I.30].
18
IV.2.5 Caractéristiques de quelques dé/multiplexeurs à base de filtres à couches
minces
IV.2.5.1 Coupleurs mono canal utilisés dans le MDL
Le coupleur mono canal est considéré comme un bloc fondamental dans plusieurs
systèmes MDL. La configuration courante utilise un coupleur à trois ports comme il est décrit
dans la section IV.2.3. La figure (1.13) montre un coupleur à trois ports. Ce composant
possède un diamètre de 5 à 6 mm et une longueur de 30 à 40 mm [I.17].
La façon la plus simple pour concevoir un démultiplexeur à plusieurs canaux est de
rassembler des séries de coupleurs comme il est présenté dans la fig (1.14). Cette architecture
peut être utilisée d’une manière inversée pour implémenter un module de multiplexeur.
Les coupleurs mono canal basés sur les filtres à couches minces, qui sont largement
utilisés et commercialisés, fonctionnent pour des espacements inter-canal de 100 GHz et 200
GHz. Les constructeurs offrent généralement un produit standard, mais ils peuvent adapter
quelques caractéristiques suivant les besoins des consommateurs ; par exemple on peut
concevoir un coupleur caractérisé soit par une haute isolation, soit par une basse dispersion
chromatique. Le tableau (I.2) récapitule les spécifications d’un coupleur 100 GHz standard et
d’un autre coupleur fonctionnant au même espacement, mais caractérisé avec une dispersion
chromatique basse [I.17].
Fig (1.13). Coupleur mono canal à trois ports (conçu par JDSU) [I.17]
Fig (1.14). Démultiplexeur à quatre canaux basé sur quatre coupleurs cascadés à base de
filtres de couches minces [I.17]
19
Tableau I.2. Caractéristiques d’un coupleur à trois ports standard et autre à basse dispersion
chromatique (conçus par JDSU) [I.17]
Ils peuvent être utilisés de la même manière que les coupleurs mono canal, mais toutes
les opérations sont faites pour une certaine bande de canaux. Deux exemples de filtres à
division de bande sont illustrés dans la fig (1.15). Un filtre 4×1 permet de transmettre quatre
canaux, en faisant une isolation entre les canaux réfléchis et les canaux transmis. Cette
isolation est achevée en supprimant un canal de chaque extrémité du gabarit du filtre. Le filtre
4×0 permet aussi de transmettre quatre canaux mais sans rejeter les deux canaux à l’extrémité
du gabarit. Le deuxième type de filtre est plus performant, cependant il présente beaucoup de
contraintes lors de son fabrication surtout au niveau des extrémités du gabarit [I.17].
Fig (1.15). Performance de deux différents coupleurs à division de bande (a) Coupleur 4×1
(b) Coupleur 4×0 [I.17]
20
Fig (1.16). Structure d’un démultiplexeur à 40 canaux basé sur cinq filtres 8×0 à division de
bande combinés avec des filtres mono canal [I.17]
Les coupleurs à division de bande sont souvent combinés avec des coupleurs mono
canal dans une architecture de modules MUX/DEMUX pour des systèmes caractérisés avec
un grand nombre de canaux. La fig (1.16) présente l’architecture d’un démultiplexeur à 40
canaux.
Ce dispositif est composé de cinq filtres 8×0 à division de bande, ce qui permet de
réduire au maximum le nombre de réflexions. La suppression des canaux reste une opération
indésirable parce qu’elle a l’effet de réduire la bande passante et par conséquent le débit d’une
communication disponible. Ceci a conduit à développer d’autres types de filtres avec des
pentes raides pour rejeter seulement un seul canal au pire des cas [I.17]. Le tableau I.3
présente les caractéristiques de trois différents coupleurs à division de bande commercialisés
fonctionnant à 100 GHz.
Pour un espacement inter-canal de 100 GHz, les modules des démultiplexeurs sont
offerts à 4, 8, 16 et 40 canaux. La figure (1.17) présente un exemple d’une réponse spectrale
d’un démultiplexeur à 40 canaux. Les caractéristiques des dé/multiplexeurs sont relativement
similaires à celles d’un seul coupleur à trois ports présentées dans le tableau I.2.
21
Fig (1.17). Réponse spectrale d’un démultiplexeur à 40 canaux fonctionnant à 100 GHz
(conçu par CIRCA) [I.17]
Fig (1.18). Principe d’un démultiplexage par diffraction dans un réseau optique : les longueurs
d’onde λ1, λ2, λ3 venants d’une fibre de transmission sont diffractées dans différentes
directions. Le réseau travaille en (a) réflexion et en (b) transmission [I.34]
22
IV.3.2 Principe de fonctionnement des réseaux basés sur la diffraction
IV.3.2.1 Ordres de diffraction
La figure (1.19) représente un réseau de fentes équidistantes, nous considérons que ce
réseau est déposé dans le vide avec un indice de réfraction égal à 1. Supposant une onde plane
incidente avec un angle θ mesuré par rapport à la normale au réseau. Chaque fente diffracte la
lumière en transmission. Dans la direction θ’, qui est un angle mesuré par rapport à la
perpendiculaire au réseau, les ondes venantes à partir des différentes fentes sont en phase si la
différence de marche optique Δ0 entre des trajets optiques successifs (AB)+(CD) est [I.34]
où m est un entier, λ est la longueur d’onde et Λ est la distance entre deux fentes successives.
Une transmission directe correspond à m= 0, alors que m=±1 correspond aux premiers ordres
de diffraction sur chaque côté de la transmission directe [I.34].
23
Fig (1.20). Configuration d’un démultiplexeur basé sur un réseau de diffraction avec une
correction d’une aberration sphérique causée par un miroir concave sphérique. La correction
est faite avec une combinaison de deux matériaux d’indices n et N présentant la même valeur
de dispersion (Configuration faite par Stimax®) [I.42]
Fig (1.21). Réponse spectrale de trois canaux adjacents avec un espacement inter-canal de 100
GHz (documentation Stimax®) [I.42]
IV.4 Composants MDL basés sur des réseaux de guides d’onde assemblés
Dans cette section, le principe de fonctionnement d’un Réseau de Guides d’onde
Assemblés (RGA) agissant comme un dé/multiplexeur est décrit. La figure (1.22) illustre la
structure de ce type de réseau.
24
Fig (1.22). Structure d’un démultiplexeur à base de réseau de guides d’onde assemblés
[I.45]
L’entrée (a) est constituée de plusieurs canaux, typiquement entre 8 et 40 canaux dans
des composants commercialisés. Les espacements inter-canal couramment utilisés dans ces
composants sont 50 GHz et 100 GHz, alors que des espacements de 25 GHz et 10 GHz ont été
déjà testés sous des conditions de laboratoire [I.43-I.44]. La longueur d’onde de
fonctionnement est 1550 nm où l’atténuation est faible des les fibres optiques. La lumière
pénètre à partir de l’entrée du guide d’onde (a) dans la Région de Propagation Libre (RPL) (b)
et se disperse pour illuminer les guides d’onde assemblés à partir d’un plan focal (voir fig
(1.23)) [I.45].
La longueur de chaque guide d’onde assemblé (GA) augmente par ΔL, où ΔL = m
λw/neff, m est un nombre entier, λw est la longueur centrale de fonctionnement et neff est
25
l’indice effectif de réfraction d’un seul mode supporté par chaque guide d’onde formant le
RGA [I.46]. Par conséquent, à la longueur d’onde centrale, un profil de phase constant expose
à l’extrémité des guides d’onde assemblés un nombre entier de cycles sur le même plan focal,
comme il est montré sur la fig (1.23).
Alors, à la fréquence centrale, la lumière est focalisée au centre du plan (e), où un
guide d’onde de sortie est positionné pour la capturer. Différentes longueurs d’onde
fournissent différentes quantités de changement de phase, et c’est à cause du changement de
la longueur du guide d’onde que les phases vont changer à travers le plan de sortie du guide
d’onde assemblé, ce qui va causer par la suite un déplacement du point focal à travers le plan
focal (e) à la fin de la région de propagation libre (RPL). Un guide d’onde de sortie est
positionné sur le plan de sortie pour capter chaque fréquence (canal) entrante [I.46].
nL
I 0 I i cos 2 (1.9)
c
L+ΔL
26
Un MUX/DEMUX basé sur un interféromètre de Mach-Zehnder présente quelques
inconvénients. Les composants cascadés ne fournissent pas une isolation suffisante par eux
même, on a besoin donc de filtres additionnels. Ces filtres peuvent être des diélectriques ou
des réseaux de Bragg. Un autre désavantage pour quelques fabricants, c’est qu’ils ont besoin
d’une stabilisation de température active [I.14].
27
La figure (1.25) montre le concept de base d’un filtre optique accordable. Ce filtre
opère sur une gamme de fréquences Δf et il est réglé électriquement pour permettre de faire
passer une seule longueur d’onde. Les paramètres principaux de ce système accordable sont
les suivants [I.50]
1. L’intervalle de réglage Δf dans lequel le filtre est accordé. Supposant que ce filtre est
réglé sur une de ces deux longueurs d’onde 1300 nm ou 1500 nm, alors une bande de
25 THz (Δλ=200 nm) sera un intervalle de réglage raisonnable pour cette opération.
Dans les réseaux où on utilise des amplificateurs optiques basés sur des fibres, un
intervalle maximale de réglage de Δλ=35 nm centré sur 1550 nm (ou Δf=4.4 THz
centré sur 193.1 THz) est utilisé.
2. L’espacement inter-canal δf, qui est la séparation fréquentielle minimale nécessaire
entre les canaux pour garantir une dégradation de diaphonie minimale. Le niveau de
signal de la diaphonie par rapport à un canal adjacent est généralement à environ 30
dB au dessous du signal désiré.
3. Le nombre maximal de canaux Nca est défini par le nombre maximal de canaux
équitablement espacés qui peuvent être utilisés dans un intervalle de réglage, à
condition de maintenir un bas niveau de diaphonie des canaux adjacents. Le nombre
de canaux Nca est défini comme le rapport entre l’intervalle de réglage total Δf et
l’espacement inter-canal δf, il est décrit comme suit
f
N ca (1.10)
f
4. La vitesse de réglage désigne le temps que va mettre un filtre pour changer une
fréquence par une autre. Pour des applications où la durée de réglage d’un canal est
relativement longue (quelques minutes à quelques heures), une vitesse de réglage à
milliseconde est suffisante. Cependant, si on décide de commuter un paquet
d’informations rapidement alors une vitesse de réglage à microseconde est nécessaire.
Plusieurs technologies ont été déjà étudiées pour concevoir des filtres optiques
accordables. Durant l’évolution des méthodologies des systèmes MDL, tout l’intérêt a été
focalisé sur des systèmes qui possèdent des espacements inter-canal supérieurs ou égaux à
100 GHz dans la fenêtre de transmission 1550 nm. Les filtres optiques accordables qui
possèdent des espacements inter-canal inférieurs ou égaux à 100 GHz incluent les coupleurs
directionnels 2×2 [I.51], les interféromètres de Mach-Zehnder [I.52, I.53], les filtres Fabry-
Pérot basés sur les fibres [I.54, I.55], les réseaux de guides d’onde assemblés [I.56, I.57], les
filtres Fabry-Pérot à cristaux liquide [I.58, I.59], les filtres multi réseaux [I.60, I.61] et les
filtres accordables acousto-optiques [I.62, I.63]. Dans cette section, nous décrivons trois types
de filtres optiques accordables.
28
Fig (1.26). Coupleur directionnel multi-électrodes fabriqué sur un crystal LiNbO3 [I.7]
électrodes déposé sur un crystal LiNbO3. Pour une opération d’extraction en longueur d’onde,
on a M longueurs d’onde qui vont entrer par le port 1. Une application d’un certain voltage
sur les électrodes va changer l’indice de réfraction des guides d’onde. C’est ainsi qu’on va
sélectionner une longueur d’onde (λi), qui va être couplée à un second guide d’onde, et par la
suite elle va sortir par le port 4. Les longueurs d’onde M-1 restantes passent à travers le
composant et sortent par le port 3 [I.7]. Pour insérer une longueur d’onde et la combiner avec
un faisceau entrant par le port 1, nous insérons une longueur d’onde λi dans le port 2, de telle
façon qu’on va créer un couplage au sommet du guide d’onde. Alors, le faisceau multiplexé
sort par le port 3 combiné avec les autres longueurs d’onde λ1, . . . , λi-1, λi+1, . . . , λM déjà
injectées dans le port 1 [I.7].
Les intervalles de réglage de ce type de filtres sont à environ 60 nm avec des largeurs
de bande aux allants tours de 1 nm pour chaque canal [I.7].
29
Fig (1.27). Principe de fonctionnement des filtres multi réseaux accordables [I.7]
le port 3 du circulateur de droite. Les longueurs d’onde ajoutées sortent par le port 1 de ce
circulateur et traversent les séries de réseaux accordés. Ces réseaux réglables réfléchissent à
leur tour chaque longueur d’onde pour la faire retourner au circulateur de droite dans le but de
la combiner avec d’autres longueurs d’onde [I.7].
V. Conclusion
La demande des abonnés pour de grandes largeurs de bande s’est très vite agrandi
durant ces dernières années. Le taux d’augmentation des débits de transmission a atteint les 70
% en 2007. Le facteur le plus important qui a fait augmenter cette largeur de bande est les
applications vidéo.
La méthode la plus efficace pour augmenter le débit de transmission des données dans
une liaison de transmission est d’envoyer simultanément plusieurs longueurs d’onde sur une
seule fibre monomode. Le multiplexage de division en longueur d’onde (MDL, ou WDM en
anglais) est une technique de transmission qui possède beaucoup d’avantages, parmi lesquels
nous citons la flexibilité. Ceci a permis aux fournisseurs de télécommunications d’installer
juste les équipements qu’ils ont besoin actuellement, et au fur et à mesure leurs réseaux
optiques seront améliorés d’une façon modulaire en fonction des futurs projets.
Les réseaux MDL nécessitent plusieurs composants optiques. Les dispositifs clé d’un
système MDL sont les multiplexeurs et les démultiplexeurs optiques. Ils permettent de
combiner plusieurs longueurs d’onde en une seule fibre et de les séparer à la fin de cette fibre
en réception. Différents types de multiplexeurs ont été présentés dans ce chapitre. Chaque
approche a ses propres caractéristiques et elle a bien sur des avantages et des inconvénients, il
n’existe plus de solution préférée ou une technologie qui présente un net avantage par rapport
aux autres. Le choix d’une technologie d’un multiplexeur dépend principalement des
spécifications du système de transmission utilisé. La technique de division de multiplexage en
longueurs d’onde est considérée comme une étape très importante pour incrémenter la largeur
de bande des systèmes de télécommunications existants, le taux d’augmentation possible peut
30
atteindre un facteur de 160. Les obstacles majeurs qui peuvent empêcher d’augmenter le
nombre de canaux dans un multiplexeur sont les niveaux de diaphonie et d’isolation obtenus,
et les pertes optiques dûes à l’absorption intrinsèque et au phénomène d’éparpillement de
Rayleigh. D’un autre coté, d’autres techniques de modulation et de démodulation plus
sophistiquées peuvent être utilisées, alors que la modulation à saut d’amplitude (on-off
keying, en anglais) a échoué d’apporter des améliorations en termes de gestion de la largeur
de bande de transmission.
31
Chapitre II
33
Fig (2.1).Le changement de l’indice de réfraction dans un réseau de Bragg [II.1]
Bragg ; appelée technologie d’écriture par radiations ultraviolettes (UV), a été démontrée par
Meltz et al. [II.5]. Plusieurs recherches ont été mis en œuvre pour améliorer la qualité et la
durabilité des réseaux de Bragg. Ces réseaux sont considérés comme des composants clé pour
les communications optiques modernes. Les produits commerciaux des fibres de Bragg ont
été disponibles dés l’année 1995 [II.1].
Un réseau de Bragg est une structure avec une perturbation périodique de l’indice de
réfraction dans un guide d’ondes. Ces réseaux sont fabriqués en exposant le cœur d’une fibre
mono mode à un rayonnement périodique de la lumière UV. Cette exposition induit un
changement permanent de l’indice de réfraction dans le cœur de la fibre. L’indice de
modulation dépend du rayonnement d’exposition [II.1]. La figure (2.1) montre le changement
périodique dans l’indice du cœur. Cette fibre optique de courte longueur avec une modulation
d’indice de réfraction est appelée réseau de Bragg.
La modulation de l’indice de réfraction est représentée par l’équation suivante [II.6]
2
n x , y, z n x, y, z n x, y, z cos z (2.1)
où n x , y, z est l’indice de réfraction moyen du cœur, n x, y, z est la modulation de
l’indice de réfraction, et Λ est la période de Bragg.
Une petite quantité de lumière incidente est réfléchie à chaque changement périodique
de l’indice de réfraction. La totalité des ondes lumineuses réfléchies sont combinées en une
seule gamme spectrale large de réflexion à une longueur d’onde particulière, lorsqu’un très
fort mode de couplage ait lieu. Ceci est appelée la condition de Bragg, et la longueur d’onde
sur laquelle la réflexion se produit est appelée la longueur d’onde de Bragg. Seulement les
longueurs d’onde qui vont satisfaire la condition de Bragg sont fortement réfléchies. La
34
Fig (2.2). Diagramme illustrant les propriétés d’un réseau de Bragg [II.1]
D 2neff (2.2)
où n eff est l’indice effectif de réfraction et Λ est la période de réseau. C’est la condition de la
résonance de Bragg. A partir de l’équation (2.2), nous remarquons que la longueur d’onde de
Bragg dépend de l’indice de réfraction et de la période de réseau.
Des réseaux à pas longs avec une petite excursion de l’indice de réfraction présentent
un haut pique de réflexion et une petite largeur de bande, comme nous pouvons voir sur la Fig
(2.2).
Un réseau de Bragg présente beaucoup d’avantages, parmi lesquels nous citons:
simple structure, pertes d’insertion faibles, sélectivité en longueur d’ondes large, insensibilité
à la polarisation et une compatibilité complète avec les fibres optiques monomodes de
transmission. Les réseaux de Bragg uniformes sont des filtres travaillant en réflexion. Suivant
chaque application, ces réseaux peuvent avoir des largeurs de bande variant de 0.1 nm à 10
nm. Une réflectivité centrée sur la longueur d’onde Bragg peut être conçue avec un maximum
de transmittance inférieur à 1% comme elle peut avoir un maximum de transmittance
avoisinant 99%. Les caractéristiques des fibres de Bragg comme la photosensibilité,
l’apodisation, la dispersion, le contrôle de la largeur de bande, la température et les réponses
35
de tractions, la compensation thermique ont tous été utilisés dans les systèmes de transmission
optiques et aussi dans les systèmes de détection [II.7].
En substituant l’équation (2.3) du champ modal dans l’équation (2.4), le champ électrique
E t x, y, z, t est écrit sous la forme [II.1]
E t x, y, z, t A j z expi j z A j z exp i j z e jt x , y exp iwt
(2.5)
j
où Aj+(z) et Aj-(z) sont les amplitudes des ondes incidentes et réfléchies respectivement ; βj est
la constante de propagation ; et e jt x, y est le champ du mode transverse. Et x, y, z , t est
l’une des solutions des équations de Maxwell.
L’indice de réseau est dépendant de la direction z le long de la fibre. L’indice de
réfraction n(x,y,z) dans l’équation (2.1) est écrit comme suit [P1]
2
n z n0 n0 n z cos z z (2.6)
36
où l’indice de réfraction moyen n est représenté par n0 n0 , et n0 n0 . n0 est l’indice de
réfraction du cœur sans perturbation ; n0 est la modulation d’indice moyenne (inDex
Change, DC, en anglais) ; n z est l’amplitude de la modulation d’indice (Amplitude
Change, AC, en anglais) ; z est la phase de réseau et est la période de réseau.
La distribution du champ électrique dans un réseau, Et x, y, z , t , satisfait l’équation de
propagation d’une onde scalaire. Celle-ci est déduite à partir d’une simplification des
équations de Maxwell sous une faible approximation, cette équation est donnée par [II.1]
2
t k 2 n 2 x, y, z 2 Et x, y, z , t 0 (2.7)
2
où k est la constante de propagation dans le vide.
Le champ électrique Et x, y, z , t et l’indice de réfraction n x, y, z sont substitués dans
l’équation de propagation (2.7) pour avoir les équations des modes couplés suivantes [II.1]
dA n
dz
i A m K mnt
K zmn expi m n z i A m K tmn K zmn exp i m n z (2.8)
m m
dAn
dz
i Am K mnt
K mn
z
expi m n z i Am K mn
t
K mn
z
exp i m n z (2.9)
m m
t
où K mn z est le coefficient de couplage transverse entre les modes n et m, K mnt z est donné
par [II.10]
w
t
K mn z dxdy x, y, z emt x, y e *mt x, y (2.10)
4
37
où est la perturbation de la permittivité, avec une approximation (n>>δn), Δε=2nδn.
E x, y, z , t A j z exp i j z A j z expi j z et x, y (2.11)
Nous remplaçons E x, y, z, t de l’équation (2.11) dans les équations des modes couplés (2.8)
et (2.9). Ces équations peuvent être simplifiées en deux modes, qui sont décrits comme suit
[II.8]
dR z
iˆ z R z ik z S z (2.12)
dz
dSz
iˆ z Sz ik z * Sz (2.13)
dz
où et , R(z) est le mode
incident , S(z) est le mode réfléchi, ˆ est le coefficient général de couplage propre (DC)
[II.9] et k z représente le coefficient de couplage (AC) qui permet l’échange d’énergie entre
les modes [II.9].
Le coefficient de couplage k z et le coefficient général de couplage propre ̂z sont
deux paramètres importants dans les équations des modes couplés (2.8) et (2.9). Ce sont des
paramètres fondamentaux pour le calcul de la réponse spectral des réseaux de Bragg.
Le coefficient général de couplage propre (DC) est représenté par la formule suivante [II.1]
2 1 d
neff (2.14)
2 dz
38
1 d
où est la variation de la période optique du réseau, et est la phase du réseau [II.9]. La
2 dz
syntonisation initiale (detuning, en anglais) δ est représentée comme suit [P1]
D (2.15)
1 1
2neff
D
où D 2neff est la longueur d’onde conçue de Bragg et δneff est le changement de l’indice
de réfraction.
Le coefficient de couplage k(z) est représenté par l’équation [P1]
k z neff v (2.16)
où neff est le changement de l’indice de réfraction et v est la visibilité de frange.
Il n’y a aucun signal d’entrée incident à partir de la partie droite du réseau (fig (2.3)),
S(+L/2)=0, alors qu’un signal connu est incident à partir de la partie gauche du réseau, R(-
L/2)=1. A partir de ces deux conditions, la condition initiale de réseau s’écrit avec les
équations (2.17) et (2.18). Les coefficients de réflexion et de transmission sont déduits à partir
des conditions initiales et des équations des modes couplés.
Coté gauche:
L
S 2 ?
(2.17)
R L
1
2
Coté droit :
L
R 2 ?
(2.18)
S L
0
2
39
Fig (2.3). Méthode de calcul de la réponse spectrale de réseau [II.1]
L
S
2
(2.19)
L
R
2
r 2 (2.20)
L
ik sinh B z
2
A z (2.22)
i sinh B L B cosh B L
L L
B cosh B z i sinh B z
2 2
A z (2.23)
i sinh B L B cosh B L
40
B k 2 2 k 2 2 (2.24)
B i 2 k 2 k 2 2 (2.25)
La phase de la lumière réfléchie par rapport à celle incidente est obtenue à partir des équations
(2.22) et (2.23), elle s’écrit comme suit [II.26]
tan 1 B coth B L (2.27)
A la longueur d’onde de Bragg le réseau présente un pique de réflectivité rmax, il s’écrit
comme suit [II.26]
rmax r D tanh 2 k L (2.28)
Il est évident qu’à partir de l’équation (2.28) la réflectivité des réseaux de Bragg peut
atteindre 1 si on augmente l’amplitude de la modulation d’indice et la longueur de réseau.
Le délai de groupe et la dispersion d’un réseau sont obtenus à partir de la phase
du coefficient de réflexion [P1].
Le délai de groupe τp (en ps) pour la lumière réfléchie dans un réseau est défini par [P1]
d 2 d
p (2.29)
dw 2c d
La dispersion d p (en ps/nm) est définie par [P1]
2c d 2
dp (2.30)
2 dw
II.6 Méthode des matrices de transfert
Cette méthode a été en premier lieu utilisée par Yamada [II.27] pour analyser les
guides d’onde optiques. Cette méthode est aussi utilisée pour analyser un réseau de Bragg.
Les équations des modes couplés (2.12) et (2.13) sont résolues par la méthode des matrices de
transfert pour un réseau uniforme et non uniforme. La figure (2.4)(a) est la structure idéale de
base de la méthode des matrices pour les deux types de réseau. Nous rappelons que si la
période Λ est constante la variation de l’indice de réfraction sera constante à chaque section,
et donc la matrice de transfert globale du réseau de Bragg est un multiple de toutes les
matrices élémentaires (FM) du réseau [II.9].
Un réseau non uniforme (voir fig (2.4)(b)) peut être divisé en plusieurs sections de
périodes différentes le long de la fibre. L’onde lumineuse incidente se propage à travers
chaque section non uniforme i qui est décrite par la matrice de transfert Fi.
41
Fig (2.4). Diagramme principal de la méthode des matrices de transfert (a) réseau uniforme
(b) réseau à pas variable [II.1]
Pour une structure d’une fibre de Bragg, la matrice Fi est décrite comme suit [II.9]
ˆ k
cosh B dz i sinh B dz i sinh B dz
B
Fi B
(2.31)
ˆ ˆ
i sinh B dz cosh B dz i sinh B dz
B B
où k est défini par l’équation (2.16), ˆ est défini par l’équation (2.14) et B est représentée
par les équations (2.24) et (2.25).
L’ensemble de réseau est représenté par [II.1]
R L / 2 R L / 2
S FM FM 1 ...... Fi ...... F1 S (2.32)
L / 2 L / 2
II.7 Réseaux de Bragg à pas variable
II.7.1 Principe de base d’un réseau à pas variable
Un réseau de Bragg chirpé (à pas variable) est un réseau qui possède une période
variable. Il y a deux variables qui peuvent être changées pour obtenir un réseau chirpé à partir
de l’équation λD=2neffΛ : Une première qui consiste à changer la période de Bragg ; une autre
réside dans le changement de l’indice de réfraction le long de la direction de propagation de la
fibre [P1].
La figure (2.5) montre un réseau de Bragg chirpé linéaire. Dans ce cas, la période du
réseau varie linéairement avec la position. Ceci fait que le réseau réfléchit différentes
longueurs d’onde (ou fréquences) à chaque point à travers le réseau [II.1].
42
Fig (2.5). Réseau de Bragg chirpé linéaire [II.1]
Le réseau à pas variable est modélisé par la théorie des modes couplés. L’indice de
réfraction d’un réseau chirpé a pour expression [II.9]
2z z
n x, y, z n x, y, z n x, y, z cos 2 d (2.33)
0
où est la période de Bragg, et décrit la phase instantanée du réseau chirpé. Il
n’existe pas de solutions analytiques pour les équations des modes couplés pour un réseau
chirpé. Cependant, la méthode des matrices de transfert est utilisée pour résoudre ces
équations.
Changer l’indice de réfraction neff à travers la direction z a le même effet lorsqu’on
change la période le long de cette direction. Cela veut dire que la période optique est
1 d
changée, alors que la période physique du réseau est fixée. Le terme de phase dans
2 dz
l’équation (2.14) représente un changement de la période optique lorsque la période physique
est fixée, il s’écrit comme suit [P3] [P1]
1 d optical bi
F 2 (2.35)
2 dz L
d D
n eff est l’indice effectif de réfraction, est le taux de changement de la longueur d’onde
db
conçue avec la position dans le réseau (variable chirp), F est le paramètre chirp et L est la
longueur du réseau.
Le paramètre bi représente la variation linéaire de l’indice de réfraction à chaque
section du réseau, il s’écrit comme suit [P4]
43
La valeur de bi , initialisée à zéro dans la première section, augmente par l’intermédiaire du
chirp c g1 jusqu’à ce qu’elle atteint un changement maximale de l’indice de réfraction à
l’ultime période du réseau, où M est le nombre de sections.
Le coefficient général de couplage propre (DC) peut aussi s’écrire sous la forme
suivante [II.1]
2neff 2
neff (2.37)
i
i 0 c g i 1 , i=1….M (2.38)
44
Fig (2.6). Spectres de réflectivité de deux réseaux chirpés avec des valeurs opposées en signe
de la variable chirp [P1]
d D
Les valeurs du chirp de l’indice de réfraction c g1 et la variable chirp reconstruits sont
db
0.0337 et -1.3821 nm/cm, respectivement.
45
Fig (2.7). Spectres réfléchis de réseaux de Bragg avec différentes valeurs de la variable chirp
[P4]
d D
gauche (petites longueurs d’onde). Or, si est positive, la longueur d’onde centrale du
db
réseau se déplace vers la partie droite (grandes longueurs d’onde). Nous remarquons aussi que
les deux réponses spectrales ont la même largeur de bande, mais centrées sur deux longueurs
d’onde différentes.
La figure (2.7) montre que la largeur de bande du spectre réfléchi augmente
d
proportionnellement avec la valeur de la variable chirp D . Nous notons que l’augmentation
db
de la largeur de bande spectrale résulte en une réflectivité réduite en même temps. A partir de
la figure (2.7) nous remarquons que plus on augmente la variable chirp, plus les spectres
réfléchis sont décalés par rapport à la longueur d’onde de Bragg.
Fig (2.8). Spectres réfléchis de trois réseaux chirpés linéaires avec différentes longueurs [P1]
Fig (2.9). Spectres réfléchis de trois réseaux chirpés avec différents changements d’indice de
réfraction [P1]
47
L’augmentation de δneff est limitée par la technologie de fabrication des réseaux de
Bragg. Ce paramètre peut seulement être varié dans une certaine bande limitée [II.1].
Fig (2.10). Relation entre la longueur d’onde centrale et la variable chirp pour trois réseaux de
Bragg avec différentes longueurs: L =10000(μm) (ligne solide), L = 15000(μm) (ligne en
tirets), L = 20000(μm) (ligne en points tillés), paramètres physiques: neff =1.447, δneff =
0.0004 et D = 1550 nm [II.1]
48
Fig (2.11). Relation entre la longueur d’onde centrale et le paramètre chirp pour trois réseaux
de Bragg avec différentes longueurs: L =10000(μm) (ligne solide), L = 15000(μm) (ligne en
tirets), L = 20000(μm) (ligne en points tillés), paramètres physiques: neff =1.447, δneff =
0.0004 et D = 1550 nm [II.1]
49
Fig (2.13). Structure générale d’un réseau apodisé [II.1]
k z neff g z v (2.39)
51
Les fonctions d’apodisation sont citées ci-dessous:
Fonction en cosinus surélevé [P1]
L
z 0, L
z
2 , (2.40)
gz rc 1 cos
L
L
2
z
gz exp g 2 , z 0, L (2.41)
L
où g est le paramètre de Gauss.
52
(a) Réseau chirpé non apodisé
Fig (2.18). Spectres réfléchis d’un réseau chirpé (a) non apodisé (b) apodisé en cosinus surélevé
53
(a) Réflectivité
(b) Dispersion
Fig (2.19). Représentation du (a) spectre réfléchi et (b) de la dispersion de quatre réseaux
apodisés gaussiens
dz
Ra (2.44)
LA
54
Fig (2.20). Structure générale d’un réseau de Bragg (a) mono canal (b) multi canal [P1]
L’intervalle L A - dz n’est pas exposé à la lumière UV, il peut être représenté par une
matrice de décalage de phase [P1]
i p i
exp
0
2 (2.45)
FP i
i p i
0 exp
2
4n eff
pi z (2.46)
et z représente l’intervalle L A - dz .
ˆ i p i k i p
cosh B dz i sinh B dz exp
i sinh B dz exp i
B 2 B 2
Fs i (2.47)
ˆ i p i ˆ i p
i sinh B dz exp
cosh B dz i sinh B dz exp i
B 2 B 2
Par conséquent, l’espacement en longueur d’onde est défini comme suit [P4]
2 (2.49)
2n eff L A
Le nombre de canaux d’un réseau échantillonné est défini par le rapport entre la bande
passante totale d’un réseau mono canal (non échantillonné) et l’espacement en longueur
d’onde [II.31].
N ca (2.50)
Nous définissons la bande passante totale d’un réseau par l’expression suivante [II.31]
2
n eff
1 (2.51)
n eff n dz
eff
Avec neff étant le changement de l’indice de réfraction et dz est la longueur d’une section de
réseau. Supposons que >> 1, l’équation (2.51) devient alors [II.31]
neff dz
2 (2.52)
n eff dz
A partir des équations (2.49) et (2.52), le nombre de canaux d’un réseau échantillonné est
décrit par l’équation suivante [II.31]
2 (2.53)
N ca 1
Ra
Un espacement en longueur d’onde entre deux canaux voisins est déterminé par la
période d’échantillonnage. Un espacement fréquentiel d’un réseau de Bragg échantillonné
fonctionnant à 200 GHz ( =1.6 nm), 100 GHz ( =0.8 nm), 50 GHz ( =0.4 nm) et 25
GHz ( =0.2 nm) correspond à une période d’échantillonnage de 0.517 mm, 1.0267 mm,
2.0534 mm et 4.1068 mm, respectivement.
Les figures (2.21), (2.22) et (2.23) présentent les spectres des canaux réfléchis de trois
réseaux de Bragg chirpés apodisés en cosinus surélevé, avec les périodes d’échantillonnage :
4.1068 mm, 2.0534 mm et 1.0267 mm, respectivement.
Nous remarquons à partir des figures (2.21) (2.22) et (2.23) que plus le nombre de
canaux augmente plus la largeur de bande des canaux est réduite. La largeur de bande d’un
canal du réseau échantillonné simulé fonctionnant à 25 GHz est 0.25 nm, alors que la largeur
de bande d’un canal d’un autre réseau simulé fonctionnant à 50 GHz est 0.4 nm.
56
Fig (2.21). Canaux réfléchis d’un réseau échantillonné fonctionnant à 25 GHz [P1]
Fig (2.22). Canaux réfléchis d’un réseau échantillonné fonctionnant à 50 GHz [P1]
Nous rappelons que la largeur de bande d’un spectre est une bande passante calculée à mi
puissance sur une échelle linéaire, elle est aussi calculée à -3 dB sur une échelle
logarithmique.
La figure (2.24) montre les valeurs des piques de réflectivité et des largeurs de bande
des canaux réfléchis d’un réseau de Bragg échantillonné apodisé en cosinus surélevé,
fonctionnant à 100 GHz. Les paramètres physiques du réseau sont récapitulés dans le tableau
II.3.
57
Fig (2.23). Canaux réfléchis d’un réseau échantillonné fonctionnant à 100 GHz [P1]
Fig (2.24). Valeurs des piques de réflectivité et des largeurs de bande d’un réseau
échantillonné fonctionnant à 100 GHz [P1]
Ce réseau fournit cinq canaux centrés sur 1552.34 nm, 1551.52 nm, 1550.7 nm,
1549.88 nm et 1549.05 nm. La valeur du maximum de réflectivité de ces canaux est entre 0.8
et 0.9, et la largeur de bande est à environ 0.64 nm (voir fig (2.24)). Le nombre de canaux
peut être connu sans passer par le traçage de spectre, il est obtenu en utilisant l’équation
(2.53). Etant donné qu’on a le nombre de sections, M, et la longueur totale du réseau, L, nous
pouvons calculer la longueur des sections du réseau par l’équation suivante
L
dz (2.54)
M
58
Tableau II.3 Paramètres physiques du réseau échantillonné [P2]
Fig (2.25). Illustration du décalage périodique de phase dans un réseau échantillonné [II.32]
10
dz 0.033 cm
300
dz 0.033
Ra 0.32
L A 0.10267
Le nombre de canaux est approximativement égal à cinq, il est calculé par l’équation (2.53)
comme suit
2 2
N ca 1 1 5.25
Ra 0.32
59
2 , (2.55)
m m m 0,1,2,.........M d 1
Md
0 D 2
DPP (2.56)
M d 2 n eff M d L A
Les figures (2.26), (2.27) et (2.28) montrent les canaux réfléchis de trois réseaux
échantillonnés avec les rapports de densification : Md=1, Md=2 et Md=3, respectivement.
Nous notons que la période d’échantillonnage est égale à 1.0267 mm (Δf=100 GHz).
Lorsque Md est mis à 1, le réseau multi canal ne présente aucun décalage de phase et
l’espacement entre canaux adjacents DPP est égal à 0 (voir fig (2.26)). m possède des
décalages de phase de {0, π} dans une période 2×LA pour Md = 2 (voir fig (2.27)) et des
décalages de phase de {0, 2π/3, 4π/3} dans une période 3×LA pour Md = 3 (voir fig (2.28)).
60
Fig (2.27). Canaux réfléchis d’un réseau échantillonné pour Md=2
0
L’espacement entre canaux est réduit à DPP 0.4nm et le nombre de canaux
2
se double pour Md = 2 (voir fig (2.27)), alors que ce nombre se triple lorsque l’espacement
entre canaux est réduit à DPP 0 0.2nm pour Md = 3 (voir fig (2.28)).
3
61
II.9.2.3 Technique de changement linéaire de phase
Une autre méthode pour réduire l’espacement entre canaux 0 est introduite, elle
consiste à appliquer un changement linéaire de phase c z à chaque section du réseau.
Cette phase est décrite comme suit [P3]
c ( z ) c p i 1 , i 1....M (2.57)
Fig (2.29). Illustration de la densification des canaux en appliquant un chirp de phase linéaire
62
II.10 Techniques de décalage des canaux
II.10.1 Ajout d’une extra phase
La première méthode proposée consiste en l’ajout d’une extra phase constante e au
décalage de phase p i , l’équation (2.46) devient
4n eff
pi z e (2.59)
Une extra phase permet de déplacer les canaux avec un décalage de t . La relation
entre le déplacement en longueur d’onde t et l’extra phase e est décrite par l’équation
suivante [P4]
e
t (2.60)
2
Les tableaux II.4 et II.5 illustrent un exemple d’un canal déplacé centré auparavant sur
1553.1215 nm. A partir du tableau II.4, nous pouvons remarquer que le canal se déplace vers
le coté droit (vers les hautes longueurs d’onde) de la longueur d’onde 1553.1215 nm tant que
Tableau II.4. Exemple d’un déplacement positif d’un canal (longueur d’onde centrale initiale
1553.1215 nm)
Tableau II.5. Exemple d’un déplacement négatif d’un canal (longueur d’onde centrale initiale
1553.1215 nm)
63
la valeur de l’extra phase ajoutée est inférieure à π. Lorsque cette valeur augmente le décalage
en longueur d’onde t augmente linéairement jusqu’à ce qu’on atteint un déplacement
maximal ( t =0.4 nm) à l’extra phase π.
Dans le cas général, un déplacement en longueur d’onde maximal est obtenu lorsque l’extra
phase est prise comme suit
A partir du tableau II.5, nous remarquons le déplacement du canal vers le coté gauche
(vers les petites longueurs d’onde) de la longueur d’onde 1553.1215 nm tant que la valeur de
l’extra phase ajoutée est supérieure à π. Lorsque cette phase décroit le décalage en longueur
d’onde croit linéairement jusqu’à atteindre un déplacement maximal ( t =-0.4 nm) à l’extra
phase π.
La figure (2.30) fait une démonstration du déplacement des canaux réfléchis d’un
réseau de Bragg échantillonné pour les extra phases: π/3 et 5π/3. Suivant les résultats de
simulation présentés dans les tableaux II.4 et II.5, nous remarquons que la valeur absolue de
l’intervalle de décalage en longueur d’onde t concernant les deux extra phases appliquées
est similaire ( t =0.13 nm). Seulement, l’extra phase π/3 induit un décalage positif des
canaux (vers les hautes longueurs d’ondes) et l’autre extra phase 5π/3 induit un décalage
négatif des canaux (vers les petites longueurs d’ondes).
Fig (2.30). Démonstration du décalage des canaux réfléchis par l’ajout d’extra phase
64
k1
Définissons l’extra phase qui permet de faire un décalage positif e right comme ,
k2
k3
et l’extra phase qui permet de faire un décalage négatif e left comme .
k4
Si la somme de ces deux extra phases est égale à 2n , nous obtenons les mêmes intervalles de
déplacement mais avec de signes différents. Ceci dit est vrai si la condition (2.63) est vérifiée.
k1 k
3 2n (2.62)
k2 k4
k1 k
, 3 (2.63)
k2 k4
neff neff
D 2 neff T 2 neff L (2.64)
T T L L
(b) [II.35]
Fig (2.32). Gradients de température des refroidisseurs correspondants aux valeurs désirées de
la dispersion [II.34]
obtenu par les deux refroidisseurs, la valeur de la dispersion peut augmenter ou diminuer. En
utilisant la relation présentée graphiquement sur la figure (2.32), la valeur de la dispersion
désirée peut être obtenue.
La figure (2.31) (b) présente une autre méthode pour ajuster la longueur d’onde
centrale d’un réseau de Bragg en utilisant une bobine chauffante [II.35].
66
Fig (2.33). Spectre de réflexion d’un réseau ajusté par un dispositif d’échauffement uniforme
avec les différentes valeurs de voltage appliqué: (a) 0 V; (b) 2 V; (c) 4 V; (d) 6 V [II.35]
La longueur d’onde de Bragg est ajustée en variant le voltage appliqué aux rouleaux
de résistance. Le dispositif de chauffage est dit uniforme à cause de son rayon uniforme qui
permet de fournir une distribution thermique uniforme à travers le réseau. La figure (2.33)
présente un spectre de réflexion d’un réseau ajusté par un dispositif d’échauffement uniforme
avec les différentes valeurs de voltage appliqué : 0 V, 2 V, 4 V et 6 V. Nous remarquons que
le spectre réfléchi se déplace vers les hautes longueurs d’onde lorsqu’on augmente la valeur
du voltage.
Lorsque la température change la longueur d’onde de Bragg varie suivant l’équation
[P5]
D D ( n )T (2.65)
1
(2.66)
T
n est le coefficient thermo-optique, il est approximativement égal à 8.6 ×10-6/°C pour une
fibre dopée germanium, il s’écrit comme suit [P5]
1 n eff
n (2.67)
n eff T
67
Tableau II.6. Evolution de l’indice effectif de réfraction, la période de réseau et la longueur d’onde de
Bragg suivant un changement de température appliqué sur un réseau uniforme mono canal
Fig (2.34). Déplacement des canaux réfléchis d’un réseau échantillonné par application de
différents gradients de température
68
longueurs d’onde de 0.005 nm, 0.0831 nm et 0.1763 nm en appliquant les températures : 30
°C, 50 °C et 80 °C, respectivement.
2 1 d
(n eff n eff ) (2.68)
2 dz
z est la phase de traction, elle décrit le chirp de traction le long du réseau. Le champ de
traction axial z est relié à la phase de traction z par l’équation suivante [P6]
d 2
z (2.69)
dz
1 Pe z (2.70)
D
où est le déplacement en longueur d’onde, la valeur de D est prise à 1550 nm dans le cas
L
général, Pe est la constante photo élastique effective et z est la traction axiale
L
appliquée sur le réseau.
L est le tronçon de réseau soumis à la traction et L est la longueur totale de réseau. La
valeur moyenne de Pe est à environ 0.22 [II.33].
Expérimentalement, un réseau de Bragg à pas variable accordable est réalisé en pliant
un tronçon d’une fibre de Bragg pour la rendre sous forme de Ω [II.36].
Fig (2.35). Transformation d’une fibre de Bragg en une forme de Ω à l’aide d’un étage fixe et
un autre variable [II.36]
69
Fig (2.36). Distribution d’une traction longitudinale appliquée sur une fibre de Bragg avec
différents ΔL [II.36]
Fig (2.37). Spectres de réflexion d’un réseau chirpé mono canal avec différents ΔL
[II.36]
Comme il est montré sur la fig (2.35), une extrémité de la fibre est montée sur l’étage
fixe et l’autre extrémité est montée sur un étage horizontalement déplaçable. La fibre peut être
pliée sous forme de Ω lorsque l’étage variable est déplacé vers l’étage fixe [II.36]. La
distribution de la traction avec différents ΔL (étage variable) est représentée sur la fig (2.36).
70
Fig (2.38). Pli d’une fibre sur une de ses extrémités et création d’un bloc d’adaptation d’indice
sur l’autre extrémité [II.37]
Fig (2.39). Spectres réfléchis d’un réseau à pas variable pour différentes longueurs inactives
[II.37]
A partir de cette figure, nous remarquons qu’autour de la position L/4 (20 mm) la distribution
de la traction n’est pas parfaitement linéaire, ceci est dû au non linéarité du profil de délai de
groupe. La non linéarité croit avec l’augmentation de ΔL.
71
Fig (2.40). Spectres réfléchis restitués pour différentes longueurs inactives [II.37]
Les spectres de réflexion avec différents ΔL sont présentés dans la fig (2.37). Lorsque
ΔL croit de 0.25 mm à 12 mm, la largeur de bande du spectre de réseau varie entre 0.347 nm
et 2.133 nm, où la longueur d’onde centrale est maintenue constante avec une tolérance de 5
pm [II.36]. La valeur de la dispersion est variée entre -925.5 ps/nm et -130.1 ps/nm. Nous
notons que la largeur de bande d’un spectre est inversement proportionnelle à la dispersion.
Nous proposons une autre méthode basée sur le même principe de l’exemple
précédent, mais elle est différente dans son système expérimental. La figure (2.38) montre un
schéma général d’un compensateur accordable de la dispersion chromatique.
Le premier processus utilisé pour faire varier la dispersion est de contrôler la longueur
de réseau en faisant un contact entre le bloc d’adaptation d’indice et une face polie à
l’extrémité de la fibre [II.37]. Par conséquent, la lumière ne peut pas atteindre le tronçon de
fibre derrière le bloc. Lorsque la longueur inactive (d) augmente, la largeur de bande de
réseau devient étroite (voir fig (2.39)). Pour récupérer la largeur de bande initiale, l’extrémité
droite de la fibre est pliée dans une direction compressive pour générer un gradient de traction
négatif, alors le spectre s’élargie vers les courtes longueurs d’onde [II.37] (voir fig (2.40)). Ce
réseau accordable est opérationnel dans un régime de dispersion négative, le signal est alors
inséré à partir de l’extrémité des hautes longueurs d’onde du réseau. Puisque la fibre de Bragg
est pliée à partir du coté des hautes longueurs d’onde (coté gauche de la fibre), cette extrémité
ne va pas subir une traction et le spectre correspondant aux hautes longueurs d’onde reste fixe
durant le processus de pli (voir fig (2.40)).
Une traction positive ou une traction négative est obtenue en utilisant le concept
d’étirement axial ou de compression axiale, respectivement. En pratique, un réseau est fixé
dans un tube cylindrique, il peut être soit compressé soit étiré par une vis micromètre
différentielle [II.38]. La fibre de Bragg est couverte de métal suivant un processus de
métallisation décrit dans [II.38], et elle est soudée de ses deux extrémités, à des tubes
72
Fig (2.41). Emballage d’un réseau de Bragg soumis à une traction axiale en utilisant
une vis micromètre [II.38]
Fig (2.42). Exemple d’une traction appliquée sur le réseau de Bragg emballé suivant (a) Un
étirement allant à une compression, (b) Une compression allant à un étirement [II.38].
d’aluminium dans le but d’éviter un éventuel glissement de la fibre vers l’extérieur. La fibre
de Bragg emballée dans le dispositif présenté dans la fig (2.41) fournit un intervalle de
déplacement en longueur d’onde à environ 15 nm [II.38].
L’effet de la traction appliquée sur le réseau de Bragg emballé est présenté sur la fig
(2.42) pour les différentes positions de la vis. Suivant la position relative entre la vis et la
longueur d’onde centrale, la force de compression déplace le spectre de réseau vers les basses
longueurs d’onde (< 1552 nm), or la force d’étirement implique un déplacement du spectre
vers les hautes longueurs d’onde (> 1552 nm).
Le réseau de Bragg accordable emballé dans le dispositif présenté sur la fig (2.41) est
l’élément le plus important dans un multiplexeur optique à insertion/extraction (MIE)
reconfigurable. Ce composant peut être réglé suivant une longueur d’onde insérée ou extraite.
Le principe de base d’un MIE sera présenté dans la section II.11.3.
Prenons un quatrième exemple, où nous allons appliquer théoriquement une traction
uniforme (constante) le long d’un réseau de Bragg à pas variable multi canal, les tractions
appliquées sont : ε1(z)=100 με, ε2(z)=-100 με et ε3(z)=400 με.
73
Fig (2.43). Démonstration d’un déplacement multi canal par application de tractions
uniformes
L’objet de cet exemple est de déplacer les canaux sans affecter leurs largeurs de bande,
pour cette raison nous avons appliqué un profil de traction constant. Ce réseau échantillonné
peut être utilisé comme un filtre multi canal accordable en longueur d’onde dans un système
de transmission multiplexé en longueurs d’onde. La figure (2.43) démontre la capacité de
réglage des canaux de ce réseau en utilisant des tractions constantes.
Les canaux réfléchis, espacés de 0.4 nm, sont décalés par 0.0328 nm, -0.0328 nm et
0.1366 nm pour les tractions appliquées : 100 με, -100 με et 400 με, respectivement. Nous
remarquons que les décalages en longueur d’onde obtenus sont petits, un intervalle de
déplacement plus large peut être atteint en augmentant la force de la traction appliquée.
74
exemples comme : le résonateur en anneau à quatre ports utilisant un coupleur à fibres
optiques [II.39], un miroir de boucle optique non linéaire [II.40], un démultiplexeur à fibre
soudée suivi d’un filtre passe bande à base de réseau de fibres [II.41], une combinaison de
circulateurs optiques (CO) et des coupleurs de multiplexage en longueurs d’onde [II.42], une
combinaison d’un réseau de Bragg chirpé, un réseau apodisé et des circulateurs optiques (CO)
[II.43] [II.44].
Un DEMUX 1×4 et un MUX à insertion/extraction (avec une longueur d’onde insérée
et autre extraite) seront présentés dans les sections suivantes. Ils sont conçus pour fonctionner
à un espacement fréquentiel de 100 GHz, et ils sont simulés par le logiciel OptiSystem. Ces
composants sont obtenus en faisant une combinaison entre les réseaux de Bragg uniformes,
les circulateurs optiques (CO) et les commutateurs optiques (CoO).
(b)
Les spectres des quatre canaux issus des quatre ports de sortie du DEMUX sont
représentés dans la fig (2.45), pour un espacement fréquentiel de 100 GHz. Nous notons que
les puissances optiques correspondantes aux canaux centrés sur les longueurs d’onde
76
(c)
(d)
Fig (2.45). Spectres des canaux sortants d’un DEMUX 1×4 (a) λ1=1553.3288 nm, (b)
λ2=1552.5244 nm, (c) λ3=1551.7208 nm (d) λ4=1550.918 nm (Résultats obtenus par
OptiSystem)
77
démultiplexées sont : -0.3457 dBm (λ1=1553.3288 nm), -1.1518 dBm (λ2=1552.5244 nm), -
1.5815 dBm (λ3=1551.7208 nm) et -1.5647 dBm (λ4=1550.918 nm). Nous remarquons que les
niveaux d’isolation des quatre canaux démultiplexés, par rapport aux canaux adjacents, sont
différents, ceci est dû aux différents trajets optiques empruntés par chaque longueur d’onde
insérée dans le port d’entrée du DEMUX avant qu’elle soit démultiplexée. Les pertes
d’insertion ont été négligées dans nos simulations.
Fig (2.46). Structure interne d’un multiplexeur à insertion/extraction (conçu par OptiSystem)
Fig (2.47). Circuit de test d’un multiplexeur à insertion/extraction (simulé par OptiSystem)
79
Le principe de fonctionnement du MIE présenté dans la fig (2.46) est démontré en
utilisant un circuit de test simulé par OptiSystem présenté dans la fig (2.47). Quatre lasers à
onde continue avec une puissance d’émission de 0 dBm chacun sont insérés dans les quatre
ports d’entrée d’un multiplexeur avec les longueurs d’ondes suivantes :
(a)
(b)
80
(c)
Fig (2.48). Spectres de sortie issus des différents ports du MIE testé (a) Port d’entrée (b) Port
de sortie (c) Port d’extraction (résultats obtenus par OptiSystem)
λ1=1553.33 nm, λ2=1552.52 nm, λ3=1551.72 nm et λ4= 1550.92 nm. Les quatre canaux sont
multiplexés, et ensuite injectés dans le MIE (ADD_DROP_MUX). La longueur d’onde
λ1=1553.33 nm est choisie pour être extraite et une autre longueur d’onde λ5=1550.12 nm est
insérée. L’analyseur WDM (Wavelength Division Multiplexing) et l’analyseur de spectre
optique sont utilisés pour voir les résultats de simulation de notre circuit de test.
La figure (2.48) montre les spectres de sortie issus des différents ports du MIE testé,
en occurrence, les spectres des canaux multiplexés injectés au port d’entrée, le spectre du
canal extrait et les spectres des canaux de sortie avec le canal ajouté. Le niveau d’isolation du
canal rejeté par rapport aux canaux adjacents, dans le cas d’une opération d’extraction, est
supérieur à 9 dB (voir fig (2.48) (c)). Les trois autres spectres des canaux réfléchis
apparaissent avec le spectre du canal extrait à cause des réflexions résiduelles présentes dans
le circulateur et le réseau uniforme (voir fig (2.48) (c)). Le niveau d’isolation du canal extrait
(λ1=1553.33 nm) par rapport au canal ajouté (λ5=1550.12 nm) est 20 dB, ce qui est acceptable
(voir fig (2.48) (b)). La présence du bruit spectral centré sur 1553.33 nm (longueur d’onde
extraite) sur le port de sortie du MIE testé est dûe aux réflexions imparfaites de cette longueur
d’onde centrale dans le réseau de Bragg uniforme.
III. Conclusion
Les réseaux de Bragg ont un rôle très important dans le secteur des
télécommunications par fibre optique. La théorie des modes couplés a été choisie pour
l’analyse et la modélisation de ces réseaux, cet outil est très important car il nous a permis de
81
comprendre les propriétés optiques globales des réseaux. Les fibres de Bragg ont été simulées
par la méthode des matrices de transfert.
Les caractéristiques spectrales des réseaux de Bragg chirpés linéaires ont été analysées
et discutées dans ce chapitre. Les réseaux à pas variable présentent quelques avantages par
rapport aux réseaux uniformes, les réseaux chirpés sont caractérisés par une grande largeur de
bande, un délai de groupe et une dispersion linéaires, etc. Ces propriétés permettent de
retrouver ce type de réseau dans plusieurs applications dans les systèmes de transmission
optique.
L’apodisation des réseaux de Bragg a été aussi introduite dans ce chapitre. Les rides
présentes sur les réponses spectrales sont réduites et les lobes secondaires sont supprimés
grâce à un couplage variable. La linéarité du délai de groupe et de la dispersion peut être
améliorée en utilisant une apodisation des réseaux de Bragg.
Les spectres des canaux des réseaux échantillonnés simulés avec différentes valeurs de
paramètres physiques ont été présentés et analysés. Nous avons présenté aussi dans ce
chapitre quelques techniques de densification des canaux dans le but de concevoir des réseaux
de Bragg échantillonnés permettant de filtrer un nombre maximum de canaux optiques dans
une liaison multiplexée en longueurs d’onde dense (DWDM, en anglais).
Les réseaux de Bragg peuvent être réglés sur des longueurs d’onde prédéfinies en leur
appliquant des tractions ou des gradients de température. Quelques travaux expérimentaux ont
été présentés pour démontrer le principe de fonctionnement d’un réseau de Bragg accordable.
A la fin de ce chapitre, un démultiplexeur 1×4 et un multiplexeur 4×1 à insertion/extraction
basés sur des réseaux de Bragg uniformes, des circulateurs optiques, des combineurs en
puissance et des commutateurs optiques ont été présentés et discutés. Ces dispositifs optiques
ont été conçus et simulés par le logiciel OptiSystem.
82
Chapitre III
84
II. Algorithmes métaheuristiques
Les métaheuristiques constituent une classe de méthodes qui fournissent des solutions
de bonne qualité en un temps raisonnable à des problèmes combinatoires réputés difficiles
pour lesquels on ne connaît pas de méthode classique plus efficace. On appelle
métaheuristiques (du grec, meta=qui englobe) des méthodes conçues pour échapper aux
minimas locaux. Le terme meta s’explique aussi par le fait que ces méthodes sont des
structures générales dont il faut instancier les composants en fonction du problème, par
exemple ; le voisinage, les solutions de départ ou les critères d’arrêt. Les métaheuristiques
sont généralement des algorithmes stochastiques itératifs, qui progressent vers un optimum
global, c'est-à-dire l’extremum global d’une fonction en évaluant une certaine fonction
objectif. Elles se comportent comme des algorithmes de recherche, tentant d’apprendre les
caractéristiques d’un problème afin d’en trouver une approximation. L’intérêt croissant
apporté aux métaheuristiques est tout à fait justifié par le développement des machines avec
des capacités calculatoires énormes, ce qui a permis de concevoir des métaheurstiques de plus
en plus complexes qui ont fait preuve d’une certaine efficacité lors de la résolution de
plusieurs problèmes à caractère difficile [III .1].
Il existe de nombreuses métaheurstiques allant de la simple recherche locale à des
algorithmes plus complexes de recherche globale. Ces méthodes utilisent un haut niveau
d’abstraction, leur permettant d’être adaptées à un large éventail de problèmes d’optimisation
combinatoire. Nous pouvons partager les méthodes heuristiques en deux catégories. Celles qui
permettent de déterminer un minimum local, et celles qui s’efforcent de déterminer un
optimum global [III .1].
On appelle méthode (ou algorithme de recherche) locale celle qui converge vers un
minimum local. Les méthodes de recherche locale, appelées aussi méthodes de recherche par
voisinage, partent d’une solution initiale et, par raffinements successifs, construisent des
suites de solutions de coûts décroissants pour un problème de minimisation. Le processus
s’arrête lorsqu’on ne peut plus améliorer la solution courante ou parce que le nombre maximal
d’itérations (fixé au départ) est atteint. Quoique, ces méthodes ne soient pas complètes (rien
n’assure qu’elles pourront trouver toutes les solutions existantes), ni n’assurent la preuve
d’optimalité, de telles méthodes parviennent très souvent à trouver des solutions de bonne
qualité dans des temps de calcul raisonnables. En effet, elles sont souvent les premières
méthodes testées sur les nouveaux problèmes combinatoires émergeant des applications
réelles et académiques [III .1]. On trouve dans la littérature de nombreuses méthodes locales.
Les plus anciennes et les plus utilisées sont : la méthode de la descente, le recuit simulé, la
recherche tabou, etc.
Contrairement aux méthodes de recherche locale, les méthodes de recherche globales
ont pour objectif d’atteindre un ou plusieurs optima globaux. Ces méthodes sont appelées
également des méthodes à population. Celles-ci sont d’une grande diversité : parmi elles on
retrouve notamment les algorithmes génétiques, les algorithmes à évolution différentielle, la
recherche dispersée, etc [III .1].
D’autre part, on peut partager les métaheuristiques en deux grandes classes : les
métaheuristiques à solution unique (i.e. évoluant avec une seule solution) et celles à solutions
85
Fig (3.1). Classification des métaheuristiques [III .1]
87
Fig (3.2). Les cinq niveaux d’organisation d’un algorithme génétique [III .2]
pi r pi
i 1 i 1
(3.1)
88
Fig (3.3). Probabilité de chaque individu placé sur une roulette de casino
Ensuite, la bille est lancée et s’arrête sur un individu. Les meilleurs individus peuvent ainsi
être tirés plusieurs fois et les plus mauvais ne jamais être sélectionnés. Cela peut être simulé
par l’algorithme suivant [III.3]
Algorithme : Sélection par roulette
- On calcule la somme S1 de toutes les fonctions d’évaluation d’une population;
- On génère un nombre r entre 0 et S1;
- On calcule ensuite une somme S2 des évaluations en s’arrêtant dès que r est dépassé;
- Le dernier individu dont la fonction d’évaluation vient d’être ajoutée est sélectionné.
89
II.1.3.3 Sélection par tournoi
Soit une population de m individus. On forme m paires d’individus. Ensuite, il faut
déterminer un nouveau paramètre, à savoir la probabilité de victoire du plus fort. Cette
probabilité représente la chance que le meilleur individu de chaque paire soit sélectionné. En
principe, cette probabilité doit être grande. Par expérience, une valeur idéale se situe entre
70% et 100%. A partir des m paires, on détermine ainsi m individus pour la reproduction
[III.6].
II.1.3.4 Elitisme
A la création d’une nouvelle population, il y a de grandes chances que les meilleurs
chromosomes soient modifiés, et donc perdus après les opérations d’hybridation et de
mutation. Pour éviter cela, on utilise la méthode élitiste. Elle consiste à copier un ou plusieurs
des meilleurs chromosomes dans la nouvelle génération. Ensuite, on génère le reste de la
population selon l’algorithme de reproduction usuel. Cette méthode améliore
considérablement les algorithmes génétiques, car elle permet de ne pas perdre les meilleures
solutions [III.5].
II.1.4 Croisement
Le phénomène de croisement est une propriété naturelle de l’ADN (Acide
désoxyribonucléique : ADN est le support de l’hérédité car il constitue le génome des êtres
vivants et se transmet en totalité ou en partie lors des processus de reproduction). Il a pour
objectif d’enrichir la diversité de la population en manipulant les composantes des individus,
c’est-à-dire les chromosomes. C’est par analogie avec la notion d’hydridation de Darwin qu’a
été conçu cet opérateur. Classiquement, les croisements sont envisagés avec un couple
d’individus parents et génèrent deux enfants [III.1].
Cet opérateur favorise l’exploration de l’espace de recherche. Considérons deux gènes
A et B pouvant être améliorés par mutation. Il est peu probable que les deux gènes améliorés
A’ et B’ apparaissent par mutation chez un même individu. Par contre, l’opérateur de
croisement permettra de combiner rapidement A’ et B’ dans la descendance de deux parents
portant chacun un des gènes mutants. Il est alors possible que la présence simultanée des deux
gènes produise un individu encore plus adapté. L’opérateur de croisement assure donc le
brassage du matériel génétique et l’accumulation des mutations favorables. En termes plus
concrets, cet opérateur permet de créer de nouvelles combinaisons des paramètres des
composants [III.2].
90
Fig (3.4). Croisement en un point [III.2]
Une fois la génération intermédiaire à moitié remplie, les individus sont aléatoirement
répartis en couples hermaphrodites (Individu morphologiquement male et femelle, soit
alternativement soit simultanément. Par exemple un escargot). Les chromosomes des parents
sont alors copiés et recombinés de façon à former deux descendants possédant des
caractéristiques issues des deux parents. On forme ainsi la génération t + 1 [III.4].
II.1.5 Mutation
Une mutation est l’inversion d’un bit aléatoire dans un chromosome. Un bit à 1
prendra la valeur 0 et réciproquement (voir fig (3.6)). Les mutations jouent le rôle de bruit en
introduisant des innovations et empêchent l’évolution de se figer. Elles assurent en outre une
recherche aussi bien globale que locale, selon le poids et le nombre des bits mutés. Plus un
individu a été muté, plus la recherche s’éloignera du voisinage local de ce même individu non
muté. De plus, elles garantissent mathématiquement que l’optimum global peut être atteint
[III.3].
D’autre part, une population trop petite peut s’homogénéiser à cause des erreurs
stochastiques (aléatoires) : les gènes favorisés par le hasard peuvent se répandre au détriment
des autres. Une telle situation risque d’aboutir à des résultats qui ne seront pas forcément
optimaux. Ce mécanisme de l’évolution, qui existe même en l’absence de sélection, est connu
sous le nom de dérive génétique. Les mutations permettent de contrebalancer cet effet en
introduisant constamment de nouveaux gènes dans la population [III.4].
Il existe de nombreuses méthodes simulant une mutation. Souvent la probabilité de
mutation pm par bit et par génération est fixée entre 0,001 et 0,01. On peut prendre également
91
Fig (3.6). Mutation [III.2]
une valeur pm = 1/Lm1 où Lm1 est le nombre de bits composant un chromosome. Il est possible
d’associer une probabilité différente à chaque gène. Ces probabilités peuvent être fixes ou
peuvent évoluer dans le temps [III.1].
Si, dans un environnement stable, il est préférable d’avoir un taux de mutation faible,
la survie d’une espèce dans un environnement subissant une évolution rapide nécessite un
taux de mutation élevé, permettant une adaptation rapide. Les taux de mutation d’une espèce
dépendent donc de leur environnement. Après divers essais, il semblerait que la meilleure
stratégie à adopter soit la méthode d’auto-adaptation des probabilités de mutation [III.2].
Lors de la genèse, les probabilités de mutation sont supposées égales à 0,1. Cette
valeur semble bonne par expérience. Au cours du déroulement de l’algorithme, les gènes et
les individus ayant des probabilités de mutation trop élevées ont tendance à disparaitre. De
même, les gènes ayant des probabilités de mutation trop faibles ne peuvent pas évoluer
favorablement et tendent à être supplantés. Les probabilités de mutation dépendent donc du
gène considéré et de la taille de la population. De plus, elles évoluent au cours du temps. Il y a
donc auto-adaptation des probabilités de mutation [III.2].
92
II.2.2 Notions
II.2.2.1 Probabilité de Boltzmann
La probabilité de Boltzmann [III.10], notée PT mesure la probabilité de trouver un
système dans une configuration i avec une énergie Ei, à une température T donnée, dans
l’espace des configurations S. Elle est définie par [III.10]
1 E
PT X i exp i (3.2)
ZT KT
où X est une variable stochastique qui désigne l’état actuel du solide, K est appelé la constante
de Boltzmann et Z(T) est appelée fonction de répartition définie par [III.10]
t
Ej
ZT exp (3.3)
j1 KT
f obj
p X , Y exp (3.4)
T
Lorsque la variation Δfobj est négative ou nulle, l’exponentielle est supérieure ou égale à 1 et
la nouvelle configuration est acceptée [III.12]
- Si Δfobj>0, p(X,Y) est comparée à un nombre aléatoire rand [0,1].
- Si rand<p(X,Y) la configuration Y est essayée.
- Sinon elle est rejetée et une autre configuration est essayée.
93
Les configurations ayant une augmentation en Δfobj, c'est-à-dire une dégradation de la
fonction objectif sont donc moins probables pour une température donnée, d’autant moins que
la température est faible [III.12].
II.2.2.4 Algorithme
Kirkpatrick a fait une analogie entre l’optimisation et le phénomène physique de
refroidissement en faisant une correspondance entre arrangements des atomes et paramètres
de conception, énergie et fonction objectif à minimiser, minimum de l’énergie et minimum
global, chaine de Markov et nombre de configurations explorées à température constante
[III.8].
Cependant, le concept de température d’un système physique n’a pas d’équivalent
direct avec le problème à optimiser. Ainsi, le paramètre température T est simplement un
paramètre de contrôle, indiquant le contexte dans lequel se trouve le système, c'est-à-dire le
stade de recherche. Le critère de Metropolis détermine si une nouvelle configuration générée
présente une variation d’une fonction objectif acceptable. Il permet aussi de sortir des
minimas locaux quand la température est élevée [III.13].
L’algorithme commence avec une température initiale élevée et une configuration
initiale prise au hasard. A l’aide d’un déplacement aléatoire, une nouvelle configuration est
générée selon chaque direction. Selon le critère de Metropolis, elle sera acceptée ou rejetée.
Ce processus est répété, à partir du dernier point accepté, un certain nombre de fois jusqu’à
obtenir l’équilibre thermique. Pendant cette phase, le vecteur pas de déplacement p est
périodiquement ajusté pour s’adapter à la fonction, le meilleur point obtenu est désigné
comme optimum courant [III.13].
La température est diminuée progressivement en générant à chaque palier un ensemble
de solutions à partir de l’optimum courant. Au début de l’algorithme, le paramètre de contrôle
T est élevé, ainsi la probabilité p(X,Y) est proche de 1 et presque toutes les dégradations de la
fonction objectif sont acceptables, ce qui favorise l’exploration. Au contraire, quand T
diminue, les remontées sont de plus en plus difficiles et seules de très faibles dégradations
sont acceptées, favorisant ainsi l’intensification. Le processus complet est répété jusqu’à ce
que le critère d’arrêt soit atteint [III.9].
En résumé, le recuit simulé utilise une double dynamique :
1- Recherche de minimas à température fixée avec la chaîne de Markov ,
2- Diminution par étape de la température.
94
II.2.3 Paramètres
La principale difficulté rencontrée dans la résolution d’un problème d’optimisation par
cette méthode est liée à la détermination du schéma de refroidissement. L’ensemble des
paramètres qui gouvernent la convergence de l’algorithme sont [III.10] :
Valeur initiale du paramètre de contrôle T0 (température initiale),
Facteur de réduction de la température rt,
Nombre d’itérations à température constante (longueur de chaine de Markov) Lm,
Taille de voisinage Ns,
Critère d’arrêt.
f
T0 (3.5)
ln Pr
95
f f opt
(3.7)
f
où f opt est la valeur de la fonction objectif à l’optimum courant, f est la moyenne des
fonctions objectives des configurations générées depuis le début de l’optimisation et est la
précision fixée par l’utilisateur.
D’autres méthodes empiriques, faciles à implanter et conduisant à de bon résultats,
parmi lesquels on cite Painton, en 1994 [III.15], qui propose d’arrêter l’algorithme lorsque
après dix décroissances successives de la température, il n’y a pas d’amélioration de la
configuration optimale. Dolan, en 1990 [III.16], propose d’arrêter l’algorithme de recuit
simulé lorsque la température de recuit est inférieure à une valeur choisie.
x iM x im
Pi, j , i 1,...., n j 1,...., L m (3.8)
cj
Algorithme de RTHu
L’algorithme est comme suit [III.27]
1. Construire la matrice des pas en utilisant (3.8) pour chaque direction.
2. Initialisation : soit un point de départ généré aléatoirement X0, la liste tabou est vide et
le meilleur point est le point initial X*=X0, f*=f(X0).
3. Générer un mouvement aléatoire faisable dans les voisinages V(X*,pk) (k=1,…,Na) de
la meilleure solution courante en utilisant les pas non tabou. Soit Xj le meilleur point
dans tous les voisinages, f(Xj) est sa fonction objectif et pj le pas utilisé pour générer
Xj.
4. Si f(Xj)≤f(X*) alors X*= Xj, f*=f(Xj), ajouter le pas pj à la liste tabou LT.
5. Mettre à jour les pas pi en utilisant (3.8) et retourner à l’étape 3.
6. Si tous les pas sont dans la liste tabou et le critère d’arrêt n’est pas satisfaisant, alors
remise à zéro de la liste tabou LT et retour à l’étape 3.
Critère d’arrêt.
L’algorithme s’arrête si :
1. La fonction objectif ne change plus de manière significative ou si
2. Un grand nombre de points est généré sans aucune amélioration de la fonction
objectif.
98
III. Métaheuristiques hybrides
III.1 Introduction
Le recuit simulé, la recherche tabou et l’algorithme génétique ont été déjà combinés
avec succès avec d’autres méthodes dans plusieurs applications. On peut classifier ces
différentes hybridations selon la taxonomie proposée par Talbi [III.29]. Cette classification
permet de comparer les métaheuristiques hybrides de façon qualitative. La taxonomie
comporte deux aspects. Une classification hiérarchique permet d'abord d'identifier la structure
de l'hybridation. Ensuite, une classification générale spécifie les détails des algorithmes
impliqués dans l'hybridation.
100
Dans les paragraphes suivants, nous allons présenter trois types de métaheuristiques
hybrides. L’algorithme génétique, le recuit simulé et la recherche tabou vont être hybridés
avec des méthodes de recherche locale et de voisinage.
x r x x x n 1 (3.12)
101
Fig (3.7). Les points de réflexion, expansion, contraction et rétrécissement d’un
simplexe en deux dimensions [III.44]
xin
où x est le centroϊde des n meilleurs points, i.e, x . Evaluer f(xr). Si f(x1)≤f(xr)
i 1 n
x e x x r x (3.13)
Evaluer f(xe). Si f(xe)<f(xr), alors remplacer xn+1 par le point xe et aller à l’étape 6 ,
sinon remplacer xn+1 par xr et aller à l’étape 6.
4. Contraction : Si f(xr)≥f(xn), alors faire une contraction entre x et le meilleur des deux
points : xn+1 et xr.
4.1 Contraction extérieure : Si f(xn)≤f(xr)<f(xn+1) (i.e, xr est strictement meilleure que
xn+1), alors faire une contraction extérieure. Calculer [III.44] :
x OC x x r x (3.14)
Evaluer f(xoc). Si f(xoc)≤f(xr), alors remplacer xn+1 par xoc et aller à l’étape 6 , sinon
aller à l’étape 5.
4.2 Contraction intérieure : Si f(xr)≥f(xn+1), alors faire une contraction
intérieure. Calculer [III.44]
x iC x x n 1 x (3.15)
Evaluer f(xic). Si f(xic)≤f(xn+1), alors remplacer xn+1 par xic et aller à l’étape 6 , sinon
aller à l’étape 5.
5. Rétrécissement : Evaluer f sur n nouveaux sommets [III.44]
102
x ' i x 1 x i x 1 , i=2,…,n+1 (3.16)
P0 S j : S j x j,i
n 1
i 1 ; x j,i R n , j 1,....., M (3.17)
j
2. Classer les sommets de chaque simplexe S , j=1,2,…, M, tel que
III.[Link] Sélection
Nous allons décrire un ensemble Q P de membres qui vont être mutés à partir de la
population courante P. Pour chaque génération, la taille de Q est la même que celle de P mais un
nombre important des membres dans P sont choisis avec une probabilité plus grande et ils vont être
103
inclus dans Q. On utilise le schéma de Baker appelé ‘sélection par classement linéaire [III.44] pour
sélectionner les nouveaux membres dans Q. Les chromosomes S j P , j=1,2,…,M, sont classés
dans l’ordre des fonctions fitness (objectif) comme dans l’équation (3.19), alors la probabilité
d’intégrer une copie d’un chromosome S j dans l’ensemble Q est calculée par [III.44]
j 1
ps S j
1
M
max max min
M 1
(3.20)
où min 2 max et 1 max 2 . En utilisant ces probabilités, la population est projetée dans
une roulette, où chaque chromosome Sj est représenté par un espace qui proportionnellement
correspond à p s S j . Les chromosomes de l’ensemble Q sont choisis en tournant la roulette
plusieurs fois jusqu’à ce que toutes les positions dans Q soient occupées.
[III.44]
où ri, i=1,….,nc, sont de vecteurs aléatoires de taille inférieure à 1, d est la distance maximale
entre les paires de parents et x k est la moyenne des kème sommets de tous les parents [III.44]
nc
1
xk
nc
x
i 1
i ,k
, k=1,….,n+1 (3.22)
3. Choisir un nombre aléatoire à partir de l’intervalle unitaire [0,1] pour chaque enfant
Ci,i=1,….,nc dans Q. Si ce nombre est inférieur à une valeur prédéfinie de la
probabilité de mutation pm, alors l’enfant est muté. Soit Im l’ensemble de ses enfants
qui ont été mutés.
4. Appliquer la procédure suivante pour chaque enfant C i x ic,k i 1 , i I m . Sélectionner
n 1
104
x im, n i x u x x ci, n i (3.23)
f i
wi m
, i 1,2,......m
f
j1
j
( y i p)
ei , i 1,2,......m (3.25)
yi p
f i f ( y i ) f (p) , i 1,2,......m
Nous pouvons citer un exemple d’un cas linéaire où la fonction f à minimiser est
linéaire, elle est décrite par l’équation suivante
f x c T x b , c R n , b R (3.26)
1 m
( y i p)
v m f i
yi p
f
j1
j
i 1
1 m
( y i p)
m c y T
i p
yi p
(3.27)
f
j1
j
i 1
1 m y i p ( y i p ) T
c
m yi p
f j i 1
j 1
A c
1 m
y i p ( y i p) T
où et A . Nous notons que la matrice A est définie positive
m
yi p
f
j1
j
i 1
tant que x Ax
T
m
( y i p) T x 2
f p v c T A c 0
T
(3.28)
106
III.5.3 Recherche heuristique de motif (RHM)
Dans cette partie, nous allons présenter en détails la nouvelle méthode de recherche
heuristique de motif (RHM). A chaque itération k correspondant au point xk appartenant à
l’ensemble Rn, la RHM utilise la DDA pour générer une direction v au point xk. Si nous
pouvions obtenir un meilleur mouvement à travers la direction v avec une certaine longueur
de pas, alors nous procédons pour l’itération suivante en faisant une mise à jour de l’itération
courante. Dans le cas contraire, la RHM [III.47] utilise un ensemble fini D des directions
d’espace positif dans Rn pour générer un maillage de points. Pour éviter de chercher
aléatoirement dans toutes ces directions, on rétrécit l’ensemble des directions de l’espace
positif D, en utilisant le paramètre de contrôle 1,1 , pour sélectionner seulement les
directions qui font partie de l’angle cos 1 formé par rapport au vecteur v ou –v. Le choix
de ces deux vecteurs est dépendant si v est une direction de descente ou pas. Nous avons alors
deux différents cas [III.48]
1. Si v est une direction de descente, nous réduisons l’ensemble D des directions
d’espace positif pour obtenir un ensemble de direction réduit D pk comme suit
D pk d D : d T v d v (3.29)
D pk d D : d T v d v (3.30)
En utilisant un pas positif avec une longueur suffisamment petite (α>0), nous pouvons
juger si v est une direction de descente ou pas. La condition suivante doit être vérifiée, sinon
la direction de descente n’est pas considérée [III.48].
f x k v f x k (3.31)
107
3. RM : Si f x k v f x k , alors utiliser l’équation (3.29) pour obtenir D pk . Dans le
cas contraire, utiliser (3.30) pour obtenir D pk . Evaluer la fonction f dans ces points de
test (ou points d’essai).
k 1 k .
5. Si la condition d’arrêt est vérifiée, alors arrêter le calcul. Sinon, mettre k=k+1 et
retourner à l’étape 2.
108
refroidissement h 0.5,0.99 et la température minimale. Mettre la température T
égale à la température maximale Tmax.
2. L’itération principale : Répéter la recherche du RS globale (étape 2.1) m1 fois. Si
le nombre de points acceptés est supérieur à mac sur un total de points de test de
m1, alors éviter l’étape de la RHM (étape 2.2) et procéder à l’étape 3.
2.1 La recherche RS globale : Considérons une solution courante xk, générer un
point d’exploration zk aléatoirement dans le voisinage de xk avec un rayon .
Générer un point de test xRS dans le voisinage de la solution actuelle xk.
2.2 La RHM locale : répéter la procédure suivante m2 fois
2.2.1 DDA : Calculer le vecteur v au point xk comme dans (3.24). Si
f x k k v f x k , alors mettre x k 1 x k k v et continuer pour la
prochaine itération de la boucle de la RHM locale.
2.2.2 Recherche de motif (RM) : Si f x k v f x k , alors utiliser
l’équation (3.29) pour obtenir D pk . Dans le cas contraire, utiliser (3.30)
pour obtenir D pk . Evaluer la fonction f dans ces points de test.
la règle k 1 k .
3. Si la longueur de la période (durée), qui correspond aux M itérations de la recherche
du RS globale, n’est pas atteinte, alors aller à l’étape 2.
4. Si le programme de refroidissement est complet (T≤Tmin) ou les valeurs d’une
fonction de deux points de test consécutifs sont très proches l’une à l’autre ou le
nombre d’itérations dépasse 50n, alors aller à l’étape 5. Sinon, réduire la
température par la relation Ti+1=λh×Ti, incrémenter m2 doucement, décrémenter le
rayon r doucement, et aller à l’étape 2.
5. A partir des meilleurs points trouvés, appliquer la méthode modifiée de Nelder-
Mead.
109
Quand cette méthode tombe dans un minimum local, elle continuera la recherche locale à
partir de ce point en acceptant des points non améliorés pour permettre l’exploration des
autres régions dans l’espace des solutions.
La méthode de la Recherche Tabou Dirigée (RTD) présentée dans cette partie diffère
par rapport aux méthodes antérieures liées à la recherche tabou continue dans plusieurs
aspects. Dans la méthode RTD, trois procédures sont utilisées ; l’exploration, la
diversification et l’intensification. Dans la recherche d’exploration, une procédure d’une
nouvelle recherche locale est introduite pour générer des mouvements d’essai (ou de test),
cette procédure est basée sur la méthode de Nelder-Mead [III.51] et la méthode de recherche
de motif (RM) récemment proposée [III.48]. En plus de ça, de nouveaux concepts d’éléments
de recherche tabou appelés ; les régions tabou (RT), les semi-régions tabou (s-RT) et la liste
tabou avec multi-classement (LT) sont introduits. Un autre élément de mémoire appelé liste
des régions visitées (LRV) est aussi introduit comme un outil pour diversifier la recherche
vers des régions non visitées de l’espace des solutions. Finalement, étant donné que l’un des
meilleurs points trouvés par l’intermédiaire des recherches d’exploration et de diversification
est très proche d’un minimum global, la recherche d’intensification est appliquée une nouvelle
fois à l’étape finale pour raffiner les solutions d’élite visitées plus loin.
La section suivante donnera une description détaillée des éléments de mémoire de la
recherche tabou. Nous introduisons, après, les stratégies de recherche locale et de voisinage
utilisées pour générer des mouvements d’essai. La méthode globale de la RTD sera décrite,
ainsi toutes les étapes de l’algorithme seront présentées dans la section finale de cette partie.
Les éléments de la LT sont classés dans un ordre croissant du plus ancien jusqu’au plus
récent, en utilisant les indices de classement I ir , i=1,….,L. Si l’élément le plus récent dans la
LT est tk, alors I kr 1 . Tant dis que si l’élément le plus ancien est tk’, alors I kr' L [III.53].
110
Les éléments de la LT sont aussi classés dans un ordre croissant suivant les valeurs de
leurs fonctions objectif, en utilisant un autre ensemble d’indices de classement I i fv ,i=1,….,L.
Si le meilleur élément dans la LT est tj, alors I jfv 1 , et si le mauvais élément est tj’, alors
I jfv' L . Nous remarquons que les éléments de la LT sont classés suivant deux différents
types de classement, alors nous considérons celle liste comme un ensemble confus et on
associe à ses éléments ti i 1 les valeurs
L
mi max mir , mifv , i=1,…..,L (3.34)
où mir , mifv 0,1 sont les valeurs classées de la recense et les valeurs classées des fonctions
objectif, respectivement, de l’élément ti [III.53].
Les régions tabou (RT) sont définies comme étant des sphères avec un rayon rTR et ses
centres sont les points de la liste tabou (LT), où rTR>0. Pour chaque RT, nous définissons une
semi-RT. C’est une région qui entoure la RT avec un rayon extérieur rSTR, avec rSTR> rTR. Si
une solution de test se retrouve dans une semi-RT, alors une procédure spécifique est
appliquée pour créer des points de test au voisinage pour éviter de retourner aux solutions
visitées antérieures. Nous présentons la procédure suivante pour cet objet [III.53].
1. Considérons un point d’essai x qui appartient à des semi-RTs avec les centres
t1,….,tυ. On calcule le centroїde t des centres de semi-RTs et la distance maximale
dmax entre x et ces centres, par l’intermédiaire des équations suivantes
1
t ti
i 1
(3.35)
d max max i 1,..., x t i (3.36)
Fig (3.8). Démonstration d’une recherche de voisinage à partir d’un point appartenant à une
semi-RT [III.53]
111
2. Construire des directions de recherche au voisinage qui sont parallèles aux axes de
coordonnées. Ces directions sont déterminées par la fonction
où ei R n est le ième vecteur unité dans l’ensemble Rn. Les points d’essai voisins sont
générés à travers ces directions de recherche avec une taille d’un pas
convenable 0 . Dans le cas où 1 , la taille du pas est choisie de telle sorte
qu’elle soit supérieure à d max rTR dans le but d’éviter une génération de points de test
à l’intérieur de la RT.
Cette information est utilisée pour diriger la recherche vers de nouvelles régions lorsque la
procédure courante de la recherche tabou (RTa) échoue de faire une amélioration ou lorsqu’
une procédure de diversification est requise. Dans un schéma de diversification, on essaye de
générer de nouveaux points d’essai à l’extérieur des régions visitées. Cependant, la génération
de points de test proches des régions les plus fréquemment visitées est une action découragée.
Une fonction est introduite pour distinguer entre une région plus fréquemment visitée et
une autre moins fréquemment visitée, cette fonction est définie comme suit [III.53]
1 e 1 (3.38)
où 0,1est une constante. On note que est une fonction strictement progressive et
bornée en haut par la valeur . Dans le paragraphe suivant, nous allons présenter une
procédure qui utilise les informations stockées dans le LRV pour générer de nouvelles
solutions. La procédure permet d’accepter un point d’essai en dehors des régions visitées, et
spécialement en dehors des régions les plus fréquemment visitées [III.53].
1. Générer un point d’essai x aléatoirement dans le domaine de recherche de f.
112
x i
2. Calculer d i , i=1,….,M, où i est définie par (3.38).
1 i
di
Si min 1i M 1 , alors accepter x, sinon, retourner à l’étape 1.
i
Un point x est accepté par la procédure citée au dessus si la condition suivante est vérifiée
x i
1 i (3.39)
i
où i=1,….,M. Cela veut dire qu’aucun point n’est accepté à l’intérieur d’une région visitée
antérieurement. De plus, un point proche des régions les plus fréquemment visitées sera
difficilement accepté. Par conséquent, plus la valeur de est grande, plus la possibilité
d’accepter un point proche des régions les plus fréquemment visitées est faible. Pour éviter un
cycle infini dans notre procédure un nombre prédéfini d’itérations est choisi [III.53].
d’une solution courante x. Cette opération est appelée la recherche de voisinage et les points
d’essai yi ip1 sont appelés les points de voisinage. Après, nous essayons d’améliorer ces
points de voisinage en exécutant une autre procédure de recherche, appelée la recherche
locale, pour générer un ensemble de points de test y p i i 1 appelés points locaux [III.53].
q
Nous notons que la RTD se base aussi dans sa partie exploration sur la recherche heuristique
de motif adaptative (RHMA), qui est la même utilisée par le RSRHM. Cette dernière a été
conçue en utilisant la méthode de direction de descente approximative (DDA). Les détails de
la méthode de Nelder-Mead (NM) et la méthode de recherche heuristique de motif (RHM) ont
été déjà présentés dans les sections III.4.2 et III.5.3, respectivement.
113
stratégies de recherche locale et de voisinage, pour constituer la méthode de recherche tabou
dirigée (RTD).
Dans la méthode RTD, trois procédures de recherche principales sont utilisées ;
Exploration, Diversification et Intensification. Dans la recherche d’exploration, nous utilisons
les stratégies de recherche locale et de voisinage, qui ont été déjà présentées dans la section
III.6.2, pour explorer l’espace de solutions. De plus, les règles imposées de la liste tabou
multi-classée (LT), la région tabou (RT) et les semi-RTs sont appliquées pour éviter de
revisiter des solutions déjà visitées récemment, ou pour éviter de tomber sur des minimas
locaux. Après, la recherche de diversification est requise pour diversifier la recherche sur
d’autres régions de l’espace de solutions, qui n’ont pas été prises en considération dans la
recherche d’exploration. Nous utilisons la liste des régions visitées (LRV) et la procédure
décrite dans la section III.6.1.3 pour gérer la recherche de diversification. Finalement, dans le
but d’explorer des régions proches autour des meilleurs points visités plus loin, la recherche
d’intensification est appliquée pour raffiner ces meilleurs points. Ces procédures de recherche
sont appliquées de telle façon pour donner à la RTD plus de chance à explorer l’espace de
recherche d’une manière plus efficace.
Recherche d’exploration
Cette recherche commence par une solution initiale x0 j à chaque itération j de la
boucle générale. Une stratégie de recherche locale et de voisinage (RNM, RHM) est utilisée
pour générer n points d’essai y i i 1 dans le voisinage d’une solution courante correspondant
n
à l’itération j. Si une meilleure amélioration est trouvée parmi ces points d’essai, on mit à jour
la solution courante et on procède à l’itération suivante de la boucle intérieure. Dans le cas
contraire où x k j est bien meilleure que tous les autres points d’essai de voisinage, la stratégie
continue à générer q points d’essai locaux y ni i 1 , où q=0 ou 1 dans le cas d’une stratégie
q
RNM et q=2 dans le cas d’une stratégie RHM. Alors, la solution courante x k j est mise à jour
pour être la meilleure des points d’essai locaux et de voisinage. La liste tabou (LT) est aussi
mise à jour ; la solution x k j va remplacer l’élément avec qui correspond la valeur de la
fonction fitness la plus petite. Si une nouvelle région est atteinte, alors la liste des régions
visitées (LRV) doit être mise à jour en ajoutant une information sur cette nouvelle région.
114
Fig (3.9). Structure principale de la méthode RTD
La boucle de recherche d’exploration est répétée au maximum linner fois, où linner est un entier
positif prédéfini.
Recherche de diversification
Ce type de recherche est réalisé lorsque le nombre d’itérations de la recherche
d’exploration dépasse la valeur maximale linner, ou lorsque l’exploration échoue d’obtenir une
amélioration dans quelques itérations consécutives. En utilisant la LRV courante, la procédure
décrite dans la section III.6.1.3 est appliquée pour générer un point d’essai x0 j 1 dans une
nouvelle région. Alors, la recherche d’exploration est répétée à nouveau, en commençant par
le point x0 j 1 .
115
IV. Synthèse et optimisation des réseaux de Bragg par la RTD
IV.1 Introduction
Les réseaux de Bragg sont considérés comme des composants fondamentaux dans
plusieurs applications en télécommunications optiques, notamment dans la compensation de
la dispersion chromatique dans une liaison mono canal ou multiplexée en longueur d’onde,
dans les applications de multiplexage et de démultiplexage en longueur d’onde (MDL) et dans
le routage optique. La synthèse d’une structure d’un réseau de Bragg à partir de son spectre de
réflexion est très utile pour la caractérisation et la conception de ce dispositif. Différentes
méthodes numériques ont été démontrées et appliquées dans ce contexte, nous pouvons citer
quelques unes : la technique de la transformée de fourier [III.54], la méthode Gel’fand–
Levitan–Marchenko [III.55], l’algorithme de Nelder-Mead(NM)[III.56], l’algorithme de
décortication par couches (ADC)[III.57], l’algorithme génétique (AG) [III.58][P5],
l’algorithme du recuit simulé (RS)[III.59][III.60] et la recherche tabou (RTa)[III.61][III.62].
L’algorithme de décortication par couches (ADC) utilise la méthode des matrices de
transfert (MMT) pour faire propager les champs électromagnétiques à l’intérieur du réseau à
base de silice à travers de petites couches. Cette méthode considère toutes les réflexions
multiples de la lumière à l’intérieur de la structure. Cependant, l’ADC utilise la phase et
l’amplitude de la réponse spectrale pour synthétiser les paramètres du réseau. La mesure du
spectre de l’intensité réfléchie est faite en utilisant un analyseur de spectre optique, tandis que
la mesure du spectre de phase impose un système expérimental plus complexe contenant un
système de mesure interférométrique [III.63].
L’algorithme génétique (AG), la recherche tabou (RTa) et le recuit simulé (RS)
reconstruisent les paramètres d’un réseau à partir du spectre de l’intensité réfléchie seulement.
Ces algorithmes utilisent une technique d’optimisation pour rechercher les paramètres
caractérisant un réseau, une technique qui peut générer un spectre réfléchi simulé adapté à un
autre supposé mesuré ou expérimental. Un spectre supposé mesuré est calculé à l’aide de la
méthode des matrices de transfert, le choix des paramètres physiques de ce réseau est lié à son
domaine d’application (compensation de la dispersion, multiplexage/démultiplexage, routage
optique). Un spectre supposé mesuré (théorique) peut devenir plus réaliste (expérimental) en
lui ajoutant un bruit gaussien aléatoire, créant des fluctuations sur le spectre en question et
générant quelques lobes secondaires. La deuxième utilité de ce bruit gaussien aléatoire est de
tester la robustesse et la performance des algorithmes de synthèse dans des environnements
générant des bruits.
Les métaheuristiques convergent lentement avec des temps de calcul considérables,
surtout lorsque le nombre de paramètres à reconstruire est grand, ce qui rend leur rentabilité
mauvaise. Citons un travail où son objectif était de tester un algorithme génétique (AG) pour
optimiser 50 paramètres, correspondants à un vecteur chirp qui varie les périodes d’un réseau
de Bragg, distribué sur 50 sections [P5]. Cet algorithme a convergé en six heures de calcul à
l’aide d’un processeur 2.4 GHz Core 2 duo avec 2 Go de RAM. L’AG a donné une fonction
objectif de 1.9, avec un total désaccord (désadaptation) entre les chirps de période originaux
et les chirps de période optimisés (voir figure (3.10)). Les résultats de cet exemple illustrent
l’échec de ce type de métaheuristiques quand le nombre de paramètres à synthétiser est grand.
116
Fig (3.10). Représentation des paramètres chirps reconstruits (ligne solide) et originaux (ligne
en points tillés) (Résultats obtenus en utilisant l’algorithme génétique (AG) dans [P5])
IV.3 Synthèse des paramètres physiques d’un réseau de Bragg mono canal par la
RTD
Dans cette partie, la RTD est appliquée pour synthétiser le spectre de réflexion d’un
réseau de Bragg à pas variable négatif (période décroissante à travers les sections du réseau)
avec une apodisation en cosinus surélevé. Généralement, un réseau à pas variable négatif est
dédié à la compensation de la dispersion chromatique induite par une fibre de transmission.
118
Tableau III.1 Comparaison entre les performances des algorithmes RSRHM AGCS et RTD [P3]
Quatre paramètres
Fonction Objectif 0.3382
0.9212 1733 0.3326
Nombre d’évaluations de la fonction
2488 5 1462
Temps de calcul (minute)
8 5
Six Paramètres
Fonction Objectif 0.7703
1.1004 3507 0.5179
Nombre d’évaluations de la fonction
3090 25 1349
Temps de calcul (minute)
17 11
Initialisation
Tolérance ‘’tol’’
Choisir les bornes supérieures et inférieures Nombre d’Itérations maximal de boucle
Pour chaque paramètre du réseau de Bragg intérieure ’linner ‘’, Nombre d’Itérations
maximal de boucle globale ‘’lmain’’
Initialisation Liste Tabou (LT) vide, Liste des Régions visitées (LRV) vide
Générer des solutions initiales aléatoires pour L, δneff ,Λ0 ,cg , αrc
Fonction objectif
Exploration
Non
Iter1 <linner
Oui
Non
Diversification x0(j+1)
Intensification (NM)
Fig (3.11) Organigramme
Oui
représentant les étapes de la RTD
pour la synthèse d’un réseau de
Iter2<maxit Bragg mono canal [P3]
dist>tol
120
Non
IV.4 Synthèse des paramètres physiques d’un réseau de Bragg multi canal par la
RTD
Dans cette partie, la RTD est appliquée pour synthétiser un réseau de Bragg à pas
variable échantillonné (multi canal). Un tel réseau est généralement utilisé comme un filtre
multi canal, il est aussi utilisé pour compenser la dispersion chromatique de plusieurs canaux
multiplexés en longueur d’onde (MDL).
Ce type de réseau est conçu en introduisant des décalages de phase entre les sections
du réseau, ces décalages sont décrits comme suit [P3]
4neff
pi z (3.40)
où z est la séparation entre deux sections consécutives de réseau. Cette séparation est
décrite comme suit [P3]
z L A dz (3.41)
121
Fig (3.12) Synthèse d’un réseau de Bragg à pas variable avec apodisation en cosinus surélevé
à partir d’un spectre réfléchi supposé mesuré [P3]
4n eff
pi z c (3.42)
L’espacement en longueur d’onde, défini par l’équation (2.49) dans la section II.9.2.1 du
chapitre II, est divisé par deux (réduit à moitié) lorsque la phase c (z ) s’incrémente
linéairement par π multipliée par un nombre impair. Cette phase est décrite comme suit [P3]
c ( z ) c p i 1 , i 1....N (3.43)
L
2
z
gz exp g
2 , z 0, L (3.45)
L
D D ( n )T (3.46)
1
(3.47)
T
124
0.9
Original
Computed
0.8
0.7
Reflected channels (A.U)
0.6
0.5
0.4
0.3
0.2
0.1
0
1550 1551 1552 1553 1554 1555 1556
Wavelength (nm)
Fig (3.14) Synthèse d’un réseau de Bragg multi canal avec une apodisation gaussienne à partir
de canaux réfléchis originaux [P3]
1 neff
n (3.48)
neff T
Tableau III.4 Valeurs des paramètres physiques reconstruits, la longueur d’onde de Bragg et les
coefficients thermo-optique et expansion thermique d’un réseau de Bragg soumis à quatre différentes
températures [P3]
Le paramètre de chirp de période est optimisé pour compenser les déviations entre la
période de réseau reconstruite et la période originale (figure (3.16) (c)). Par exemple, à la
température 30 °C la période initiale Λ0 synthétisée est 535.6062 nm. Cette valeur est
supérieure à la valeur originale à la même température (Λ0originale=535.6015 nm). Pour
minimiser cette erreur entre les deux périodes (optimisée et originale) à chaque section du
réseau, une optimisation est faite par la RTD dont le but est d’incrémenter la valeur du chirp
de période cg à -0.0404 nm. Par conséquent, la valeur de la période synthétisée à chaque
section sera réduite et plus proche de la valeur de période originale.
La figure (3.17) montre un exemple d’un réseau de Bragg reconstruit à partir d’un
spectre en réflexion supposé mesuré à 70 °C. On note que la RTD a convergé en 3 minutes
avec le même type de processeur cité auparavant. En faisant une comparaison entre le spectre
synthétisé à 70 °C et un autre spectre original à 25 °C, nous remarquons que la longueur
d’onde initiale de Bragg est déplacée de 0.732 nm avec un gradient de température de 45 °C.
Les valeurs des paramètres reconstruits à cette température et leurs déviations
correspondantes sont récapitulées dans les Tableaux III.4 et III.5, respectivement. Nous avons
calculé la valeur moyenne des quatre valeurs des coefficients thermo-optique αn et expansion
thermique αΛ citées dans le Tableau III.4.
Fig (3.15) Evolution de la longueur d’onde de Bragg d’un réseau mono canal synthétisé par la
RTD à partir des spectres en réflexion correspondants à quatre différentes températures [P3]
127
(a)
(b)
(c)
Fig (3.16) Evolution des paramètres reconstruits, (a) indice effectif de réfraction (b) période
initiale du réseau, (c) chirp de réseau, d’un réseau mono canal synthétisé par la RTD à partir
des spectres en réflexion correspondants à quatre différentes températures (30 °C, 50 °C, 70
°C et 80 °C) [P3]
Le Tableau III.6 présente une comparaison entre les valeurs moyennes des coefficients
synthétisés par la RTD, des valeurs théoriques [III.64] et d’autres valeurs obtenues par
l’algorithme de Nelder-Mead (NM) [III.56].
128
1
Reference (25°C)
Target (70°C)
0.9
Synthesized (70°C)
0.8
0.7
Reflectivity (A.U)
0.6
0.5
0.4
0.3
0.2
0.1
0
1548 1548.5 1549 1549.5 1550 1550.5 1551 1551.5 1552
Wavelength (nm)
Fig (3.17) Synthèse d’un réseau de Bragg à pas variable apodisé en cosinus à partir d’un
spectre en réflexion supposé mesuré à 70 °C [P3]
Nous notons que dans la référence [III.56] la synthèse a été faite à partir d’un spectre en
réflexion d’un réseau de Bragg uniforme ; seulement l’indice de réfraction et la période de
réseau ont été optimisés. Nous remarquons que toutes les valeurs comparées sont
approximativement égales, ce qui prouve la bonne performance de la RTD proposée.
ε(z)=a0z2+a1z+a2 (3.49)
129
Tableau III.6 Comparaison entre les valeurs moyennes des coefficients thermo-optique et expansion
thermique obtenues par la RTD [P3], l’algorithme de Nelder-Mead [III.56] et d’autres valeurs
théoriques [III.64]
αn αΛ
Par conséquent, une traction constante est appliquée sur le réseau. Les coefficients a1
et a2 sont utilisés en même temps pour créer une traction avec profil linéaire (figure (3.18)
eps= 3z+100). La force du stress appliqué devient plus grande en la comparant avec la
précédente (avec un seul coefficient). Par conséquent, la largeur de bande du réseau devient
un petit peu grande. Ce petit changement est justifié par le fait que nous avons choisi une
petite valeur du coefficient a1 (a1=3). Un changement de la largeur à mi-puissance indique que
les valeurs du délai de groupe et de la dispersion vont aussi changer. Dans le troisième cas,
nous avons appliqué une traction avec un profil non linéaire en utilisant les trois coefficients
a0, a1 et a2. A partir de la figure (3.18), nous remarquons que lorsque le profil de traction suit
la fonction ε(z)= z2+3z+100 la largeur de bande devient plus large en la comparant avec celle
obtenue en utilisant une traction linéaire.
0.9
0.8
0.7
Reflectivity (A.U)
0.6
0.5
0.4
0.3
eps=0
eps=100
0.2
eps=3z+100
0.1 eps=z 2+3z+100
eps=-z 2-3z-100
0
1548 1548.5 1549 1549.5 1550 1550.5 1551 1551.5 1552 1552.5 1553
Wavelength (nm)
130
Fig (3.19) (a) Etirement et (b) compression d’un réseau de Bragg à pas variable négatif
Les spectres réfléchis correspondants aux trois tractions citées au dessus sont déplacés
vers les hautes longueurs d’onde. Nous remarquons que la déformation du spectre se situe à
droite de la longueur d’onde centrale, ceci est dû à la pression positive appliquée sur la fibre.
Supposant que la fibre utilisée est à pas variable négatif (chirp de période négatif), dans ce
type de fibres les hautes longueurs d’onde se réfléchissent plus rapidement (un délai de
groupe court) (figure (3.19) (a)), alors que les courtes longueurs d’onde mettent plus de temps
pour se réfléchir (un délai de groupe long) (figure (3.19) (b)). Sur la figure (3.18), la
déformation du spectre se présente sur les hautes longueurs d’onde parce que la traction a été
appliquée au début de la fibre de Bragg (à partir des hautes longueurs d’onde). Une traction
négative (voir figure (3.19) (b)) permet de déplacer le spectre vers les courtes longueurs
d’onde (figure (3.18) eps=-z2-3z-100), et de créer des déformations sur le spectre à gauche de
la longueur d’onde centrale.
Les tractions linéaires et non linéaires peuvent être appliquées sur un réseau de Bragg
pour utiliser ce composant comme un compensateur accordable de la dispersion dans une
liaison mono canal ou multiplexée en longueurs d’onde. Une condition doit être respectée ; la
longueur d’onde centrale doit être fixée en changeant la valeur de la dispersion. La force de la
traction est utile lorsque la valeur de la dispersion négative d’un réseau de Bragg doit être
ajustée suivant un changement de la valeur de la dispersion positive d’une fibre monomode de
transmission.
Une traction avec un profil constant est utilisée pour déplacer un spectre sans changer
sa largeur de bande. Une traction positive et négative permet de déplacer la réflectivité d’un
réseau en direction des hautes longueurs d’onde et des courtes longueurs d’onde,
respectivement. Cette technique est appliquée pour utiliser un réseau de Bragg comme un
filtre accordable dans les systèmes de transmission multiplexés en longueurs d’onde (MDL).
131
IV.6.2 Reconstruction d’un profil d’une traction appliquée sur un réseau de
Bragg à partir d’un spectre original déformé
Dans cette partie, nous nous intéresserons à la reconstruction d’un profil d’une
contrainte mécanique appliquée sur un réseau de Bragg à pas variable en utilisant la méthode
de la Recherche Tabou Dirigée (RTD), à partir d’un autre spectre déformé supposé mesuré.
Le nombre de sections M du réseau est 15, alors il y aura 15 paramètres à optimiser : ε(z1),
ε(z2),….. ε(z15). L’intervalle de recherche initial pour chaque paramètre synthétisé est fixé à
100 με (εtarget(z)-100<x< εtarget(z)+100). Les paramètres physiques du réseau de Bragg utilisés
dans cette simulation sont les mêmes présentés dans le Tableau III.2.
150
original strain
reconstructed strain
140
130
µstrain
120
110
100
90
0 1 2 3 4 5 6 7 8 9 10
grating length (mm)
220
original strain
200 reconstructed strain
180
160
µstrain
140
120
100
80
60
0 1 2 3 4 5 6 7 8 9 10
grating length (mm)
132
-80
original strain
-100 reconstructed strain
-120
-140
µstrain
-160
-180
-200
-220
-240
0 1 2 3 4 5 6 7 8 9 10
grating length (mm)
Les différents profils des tractions synthétisés et originaux sont représentés dans les
figures (3.20), (3.21) et (3.22). Nous remarquons que l’accord entre eux est relativement bon.
Les résultats de simulation des trois profils appliqués sont récapitulés dans le Tableau III.7.
Nous rappelons qu’une augmentation des valeurs des trois coefficients a1, a2 et a3 dans le cas
d’un profil constant, linéaire et non linéaire fait augmenter par conséquent la différence de
force de la traction appliquée. Les valeurs moyennes de la force de traction pour les trois
différents profils sont : ε1(z)= 114 µε, ε2(z)= 144 µε et ε3(z)=-144 µε. Le Tableau III.7 montre
que l’erreur maximale et moyenne d’une traction est minimale lorsque la valeur moyenne de
la force de traction appliquée est petite.
La simulation faite par la RTD a donné des fonctions objectif entre 1.2282×10-8 et
2.544×10-8. La méthode a convergé en 30 minutes à l’aide un processeur 2.4 GHz Core 2 Duo
avec 2 Go de RAM.
La RTD est aussi utilisée pour la reconstruction d’un profil d’une contrainte
mécanique appliquée sur un réseau de Bragg à pas variable multi canal. Ce réseau
échantillonné est divisé en 45 sections, alors il y aura 45 paramètres du profil de traction à
synthétiser à partir d’un spectre original déformé. A partir de la figure (3.23), nous
remarquons qu’il y a un bon accord entre les canaux réfléchis reconstruits et originaux. La
traction appliquée suit la fonction ε(z)=-3z2-50z-100.
Dans cet exemple, nous avons soumis le réseau échantillonné à une forte traction
négative, la force moyenne de la traction est -423 με. Une forte traction implique que la
déformation sur le spectre est plus significative, or des rides (fluctuations) se présentent à
gauche de la longueur d’onde Bragg.
La méthode hybride a donné une fonction objectif de 1.5064×10-6 après 80007
évaluations. La RTD a convergé en 10 heures à l’aide d’un processeur 2.4 GHz Core 2 Duo
avec 2 Go de RAM. Il est évident que lorsque le nombre de sections devient plus grand, le
calcul sera plus complexe et le temps d’exécution augmente.
133
Tableau III.7 Résultats de simulation obtenus par la RTD pour les différents profils synthétisés [P6]
0.8
Original
Reconstructed
0.7 without applied strain
0.6
Reflectivity (A.U)
0.5
0.4
0.3
0.2
0.1
0
1548 1549 1550 1551 1552 1553 1554 1555 1556 1557
Wavelength (nm)
Fig (3.23). Synthèse d’un réseau de Bragg multi canal avec une traction négative à partir d’un
spectre en réflexion original avec le profil de traction ε(z)=-3z2-50z-100 [P6]
A partir de la figure (3.24), nous remarquons des erreurs entre le profil original et le
profil synthétisé parce que le nombre de paramètres synthétisés est grand. La valeur de la
distance d’erreur maximale est 99.2 με, alors que la distance d’erreur moyenne est 46.74 με.
V. Conclusion
Ce chapitre a été consacré à la synthèse et l’optimisation des réseaux de Bragg mono
canal et échantillonnés. La RTD, l’AGCS et le RSRHM ont été proposés pour reconstruire les
paramètres physiques de ces réseaux. Nous rappelons que la RTD et le RSRHM appartiennent
à la classe des métaheuristiques hybrides de bas niveau à relais, tandis que l’AGCS appartient
à la classe des méthodes hybrides de bas niveau co-évolutionnaires. Une comparaison entre
ces trois métaheuristiques hybrides, en termes de la fonction objectif, du nombre
d’évaluations et du temps de calcul, a permis de démontrer que la RTD est la méthode la plus
134
0
original strain
-100 reconstructed strain
-200
-300
-400
µstrain
-500
-600
-700
-800
-900
-1000
0 1 2 3 4 5 6 7 8 9 10
grating length (mm)
Fig (3.24). Profil de traction reconstruit et original pour un réseau échantillonné (ε(z)=-3z2-
50z-100) [P6]
avantageuse grâce à sa rapidité de convergence et grâce aussi à la bonne qualité des solutions
obtenues. Dans sa partie exploration, la RTD utilise les stratégies de recherche au voisinage et
locale. La recherche heuristique par motif (RHM), qui se base sur la méthode de la direction
de descente approximative (DDA), permet de générer des points d’essai au voisinage d’une
solution initiale. Si une meilleure solution est trouvée par cette stratégie, la recherche de
Nelder-Mead continue à générer des points d’essai locaux à partir de cette solution. Des
éléments de mémoire sont introduits pour stocker les meilleurs points visités en faisant une
mise à jour de la liste tabou (LT). Les solutions sont classées dans cette liste suivant un ordre
croissant de l’élément le plus ancien jusqu’à l’élément le plus récent, et selon aussi la valeur
de leurs fonctions objectif. La liste des régions visitées (LRV) contient les informations
historiques concernant les régions visitées antérieurement. Une région visitée est caractérisée
par une sphère d’un certain rayon, et une fréquence de visites de cette région. Ces
informations sont utilisées pour diriger la recherche vers de nouvelles régions lorsque la
procédure d’exploration échoue de faire une amélioration, or une procédure de diversification
doit être exécutée. Une procédure d’intensification est appliquée dans la dernière étape de la
RTD pour améliorer les meilleurs points stockés dans la liste tabou (LT).
A la fin de ce chapitre, la recherche tabou dirigée (RTD) a été utilisée pour estimer les
coefficients thermo-optique et expansion thermique d’une fibre de Bragg. Une contrainte
thermique est appliquée sur un réseau de Bragg échantillonné pour déplacer ses canaux vers
d’autres longueurs d’onde désirées. Un gradient de température est aussi appliqué aux deux
extrémités d’une fibre de Bragg échantillonnée pour changer la largeur de bande de chaque
canal sans déplacer sa longueur d’onde centrale, par conséquent la valeur de la dispersion de
chaque canal est modifiée. Cette caractéristique est utilisée pour une compensation accordable
135
de la dispersion chromatique dans une liaison multiplexée en longueur d’onde. La RTD a été
utilisée aussi pour reconstruire des profils de tractions positives et négatives, appliquées sur
une fibre de Bragg, à partir des spectres en réflexion déformés. L’objectif de l’application des
ces contraintes mécaniques est le même que celui d’une contrainte thermique ; un
déplacement des longueurs d’onde centrales des canaux et le changement des valeurs de leurs
dispersions correspondantes.
Les distances d’erreur, entre les paramètres synthétisés et les paramètres originaux,
minimisées par la RTD sont acceptables. Cette méthode a donné de très bons résultats malgré
qu’elle a été testée dans un environnement contenant un bruit gaussien aléatoire. Cette
métaheuristique hybride est avantageuse par rapport à d’autres métaheuristiques combinées
surtout lorsque le nombre de paramètres à reconstruire est très grand.
136
CONCLUSION GENERALE
Conclusion générale
Les réseaux de Bragg ont un rôle très important dans le secteur des
télécommunications optiques. Nous avons choisi ces composants optiques grâce aux
nombreux avantages qu’ils présentent, parmi lesquels leur structure est simple et par
conséquent leur coût de fabrication est faible. Ces réseaux peuvent aussi être accordés en
longueur d’onde et leur gamme de sélectivité est large, ils sont flexibles et leur adaptation est
parfaite avec les fibres mono mode largement utilisées. En plus de ces caractéristiques, les
fibres de Bragg sont des composants de base de quelques types de multiplexeurs et de
démultiplexeurs optiques, elles sont utilisées aussi dans le filtrage que ça soit statique ou
accordable. L’objectif de notre travail était d’analyser, de synthétiser et d’optimiser différents
types de réseaux de Bragg.
Les réseaux de Bragg uniformes, à pas variable, apodisés et échantillonnés ont été
modélisés par la théorie des modes couplés et simulés en utilisant la méthode des matrices de
transfert. Une attention particulière a été faite sur le réseau de Bragg multi canal, car il permet
un filtrage adaptatif. En premier lieu, nous avons présenté quelques techniques pour densifier
le nombre de canaux dans un réseau échantillonné dans le but de démultiplexer ou compenser
la dispersion chromatique de plusieurs canaux multiplexés dans une fibre mono mode. En
second lieu, quelques méthodes de modification des paramètres physiques (indice de
réfraction, période de réseau) d’un réseau de Bragg ont été introduites en appliquant des
contraintes thermiques et mécaniques sur la fibre de Bragg. Nous avons aussi conçu et simulé
un multiplexeur à insertion/extraction basé sur des réseaux de Bragg, des circulateurs optiques
et des commutateurs optiques en utilisant le logiciel OptiSystem.
Dans la dernière partie de cette thèse, une synthèse et une optimisation des réseaux de
Bragg mono canal et multi canal a été faite en utilisant des métaheuristiques hybrides. La
recherche tabou dirigée (RTD), l’algorithme génétique avec codage de simplexe (AGCS) et le
recuit simulé avec recherche heuristique de motif (RSRHM) ont été proposés pour
reconstruire les paramètres physiques des réseaux de Bragg. Une comparaison entre les
performances de ces trois algorithmes de différents concepts a été faite en fonction de leur
fonction objectif, du nombre d’évaluations et du temps de calcul, et elle a prouvé que la RTD
est la méthode la plus appropriée car elle a donné la distance d’erreur moyenne la plus petite
en un temps court par rapport aux autres algorithmes. Cette méthode est aussi préférée dans
les systèmes à grandes dimensions.
La recherche tabou dirigée (RTD) est une méthode basée sur des éléments de mémoire
qui permettent de stocker les meilleurs points visités en faisant une mise à jour de la liste
tabou (LT). Ces points sont classés dans cette liste suivant leur nouveauté et suivant aussi la
valeur de la distance d’erreur obtenue. La liste des régions visitées (LRV) mémorise toutes les
espaces qui ont été déjà visités. Si une mauvaise solution est stockée dans une sphère de rayon
donné la recherche tabou évitera de rechercher des solutions dans cette zone et diversifie la
recherche vers de nouvelles régions lorsque la procédure d’exploration échoue de faire une
amélioration. Les meilleurs solutions stockées dans la LT sont intensifiées à la fin de la RTD
en utilisant une méthode de recherche locale appelée la recherche de Nelder-Mead (RNM).
Nous rappelons qu’une méthode de recherche globale basée sur la direction de descente
approximative (DDA) est utilisée dans la partie exploration de la RTD.
138
La recherche tabou dirigée (RTD) a été utilisée pour reconstruire les coefficients
thermo-optique et expansion thermique d’un réseau de Bragg mono canal en faisant une
synthèse des paramètres physiques d’une fibre de Bragg à partir d’un spectre réfléchi supposé
mesurer. Un gradient de température a été appliqué théoriquement aux deux extrémités d’un
réseau de Bragg pour déplacer la longueur d’onde centrale sans modifier la largeur de bande
du spectre. Le réseau conçu peut être utile pour un démultiplexage accordable en longueurs
d’onde. Cet algorithme hybride a été aussi utilisé pour reconstruire des profils de tractions
positives et négatives de la fibre de Bragg à partir de spectres déformés originaux. La RTD a
pu reconstruire 45 paramètres avec des distances d’erreur acceptables dans un environnement
contenant un bruit gaussien aléatoire. Les résultats de simulation obtenus durant cette thèse
ont montré que l’algorithme de la recherche tabou hybride proposé est une méthode très
prometteuse qui peut être utilisée pour concevoir d’autres types de filtres optiques, notamment
les fibres optiques à cristaux photoniques utilisées pour la compensation de la dispersion
chromatique.
139
ANNEXE 1
140
for i=1:n;
w(i)=df(i)/sum(abs(df));
u(i,:)=-(S(i,:)-x);
u(i,:)=u(i,:)/norm(u(i,:));
end
for j=1:n; v=v+w(j)*u(j,:); end
v=v/norm(v);
S(n+1,:)=x+(0.5*rand)*Edg*v;
S(n+2,:)=x+(0.5+0.5*rand)*Edg*v;
for j=1:2
jn=n+j;
fv(jn)=feval(f,S(jn,:));
end
FCount=FCount+2;
[fvs,is]=sort(fv);
Stmp=S(is,:); S=Stmp; fv=fvs;
xnew=S(1,:);
fnew=fv(1);
else
xnew=x1; fnew=f1;
end
% TL and VRL Setting
inTL=1; TL(inTL,:)=x; fTL(inTL)=fx; rec(inTL)=1;
inVRL=1; VRL(inVRL,:)=x; frVRL(inVRL)=1;
x=xnew; fx=fnew;
Xold=x; Fold=fx;
%%%%%%%%%%%%%%%%%%%%%% Main Loop %%%%%%%%%%%%%%%%%%%%%%%%%
while(inNoImprv < NoImprv & itrc < MaxItrs)
[FM1,XM1,FCount1,TL1,fTL1,rec1,VRL1,frVRL1]=APS(f,n,x,fx,Edg,TL,fTL,rec,VRL
,frVRL);
TL=TL1;
fTL=fTL1;
rec=rec1;
VRL=VRL1;
frVRL=frVRL1;
FCount=FCount+FCount1;
%%%%%%%%%%%%%%% Diversification%%%%%%%%%%%%%%%%%%%%%
accpt=0;
while accpt==0
r=rand(1,n);
xnw=L+(U-L).*r;
for j=1:inVRL
dxV(j)=norm(xnw-VRL(j,:))*(1+0.25*(1-exp(-0.25*(frVRL(j)-1))));
end
if min(dxV)>= H*dVR
accpt=1;
end
end
x=xnw;
fx=feval(f,x);
FCount=FCount+1;
% Update TL
[sortrec,iTLr]=sort(rec);
TLrtmp=TL(iTLr,:); TL=TLrtmp;
fTLtmp=fTL(iTLr); fTL=fTLtmp;
rec=sortrec;
inTL=length(fTL);
for j=1:inTL
rec(j)=j;
141
end
[sortfTL,iTL]=sort(fTL);
TLtmp=TL(iTL,:); TL=TLtmp;
rectmp=rec(iTL); rec=rectmp;
fTL=sortfTL;
fTL(inTL)=fx;
TL(inTL,:)=x;
rec(inTL)=inTL+1;
inVRL=length(frVRL);
VRL(inVRL+1,:)=x; frVRL(inVRL+1)=1;
if FM1 <= Fold
inNoImprv = 0;
else
inNoImprv = inNoImprv+1;
end
Xold=XM1; Fold=FM1;
itrc=itrc+1;
end
%%%%%%%%%%%%% Intensification Search %%%%%%%%%%%%%%%%
[fTLs,ifTL]=sort(fTL);
TLtmp=TL(ifTL,:);
TL=TLtmp; fTL=fTLs;
XMin=TL(1,:); XM=TL;
FMin=fTL(1);
% Applying N-M method starting from Best Point
if n<5
maxitt=100*n;
budget=200*n;
else
maxitt=5000*n;
budget=10000*n;
end
for j=1:1
[x1,FCoun1,f1]=NeMe(n,f,XM(j,:),maxitt,budget);
if f1 < FMin
XMin1=x1(:,1);
FMin=f1;
end
FCount=FCount+FCoun1;
end
%%%%%%%%%%%%%%%%%%%%%%%%% END %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Fibre simulation parameters
neff = 1.447; %core index of photosensitive fibre
l = 10e6;%length of grating in micrometers
delnbar=4e-4;
M=15;
walD = 1550; %design wavelength
wal1 = 1548;
wal2 = 1552;
step =100;
delz1 = l/M;
c = 2.99793e8; %Speed of light
lmin=0;
zel= [lmin:(l-lmin)/M:l];
wal = [wal1:(wal2-wal1)/step:wal2];
wal3 = [wal1:(wal2-wal1)/(step-1):wal2];
wal4 = [wal1:(wal2-wal1)/(step-2):wal2];
alpha=40;
H=rand(1)*1e-5;
G=80;
142
m=2;
B=0.03%;
raised=0.5;
for (r = 1:step+1)
w = wal(r);
f1_2 = [1 0; 0 1];
for (i = 1:M)
z=zel(i);
%%%%%%%%%%%%%%%%%CGP%%%%%%%%%%%%%%%%
A=535.6;
cg=-0.0335;
CGP=x1(3)*(i-1);
CGP1(i)=x1(2)+CGP;
g_z=x1(4)*(1+cos(pi*(z-l/2)/l));
sig1=((2*pi*neff)*(1/w)-pi/CGP1(i))+((2*pi*x1(1))/w);
gammab1=((((g_z*pi*x1(1))/w)^2-((sig1)^2)))^.5;
f11_2=(cosh(gammab1*delz)-sqrt(-1)*(sig1/gammab1)*sinh(gammab1*delz));
f22_2=conj(f11_2);
f21_2=(((sqrt(-1)*g_z*pi*x1(1))/(gammab1*w))*sinh(gammab1*delz));
f12_2=conj(f21_2);
f1_2=f1_2*[f11_2 f12_2;f21_2 f22_2];
end
r3_2(r) = f1_2(2,1)/f1_2(1,1); %amplitude reflection coefficient
R3_2(r) = (abs(r3_2(r)))^2; %power reflection coefficient of single grating
end
%%%%%%%%%%%%%%standard deviation%%%%%%
dcg=cg-x1(3)
dA=A-x1(2)
d1n=delnbar-x1(1)
draised=raised-x1(4)
figure(1), hold on,
plot(wal, R3_2, 'b--')
grid
xlabel('Wavelength (nm)')
ylabel('Computed Reflectivity')
function g=object4(x)
close all
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%COMPUTED
REFLECTIVITY%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Fibre simulation parameters
neff = 1.447; %core index of photosensitive fibre
l = 10e6;%length of grating in micrometers
delnbar=4e-4;
M=15;
walD = 1550; %design wavelength
wal1 = 1548;
wal2 = 1552;
step =100;
c = 2.99793e8; %Speed of light
lmin=0;
zel= [lmin:(l-lmin)/M:l];
wal = [wal1:(wal2-wal1)/step:wal2];
wal3 = [wal1:(wal2-wal1)/(step-1):wal2];
wal4 = [wal1:(wal2-wal1)/(step-2):wal2];
alpha=40;
H=rand(1)*1e-5;
143
G=80;
m=2;
B=0.03%;
A=535.6;
raised=0.5;
for (r = 1:step+1)
w = wal(r);
f1_2 = [1 0; 0 1];
for (i = 1:M)
z=zel(i);
%%%%%%%%%%%%%%%%%CGP%%%%%%%%%%%%%%%%
A=535.6;
cg=-0.0335;
CGP=x(3)*(i-1);
CGP1(i)=x(2)+CGP;
g_z=x(4)*(1+cos(pi*(z-l/2)/l));
sig1=((2*pi*neff)*(1/w)-pi/CGP1(i))+((2*pi*x(1))/w);
gammab1=((((g_z*pi*x(1))/w)^2-((sig1)^2)))^.5;
f11_2=(cosh(gammab1*delz)-sqrt(-1)*(sig1/gammab1)*sinh(gammab1*delz));
f22_2=conj(f11_2);
f21_2=(((sqrt(-1)*g_z*pi*x(1))/(gammab1*w))*sinh(gammab1*delz));
f12_2=conj(f21_2);
f1_2=f1_2*[f11_2 f12_2;f21_2 f22_2];
end
r3_2(r) = f1_2(2,1)/f1_2(1,1); %amplitude reflection coefficient
R3_2(r) = (abs(r3_2(r)))^2; %power reflection coefficient of single grating
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%MEASURED
REFLECTIVITY%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for (r = 1:step+1)
w = wal(r);
f1 = [1 0; 0 1];
for (i = 1:M)
z=zel(i);
%%%%%%%%%%%%%%%%%CGP%%%%%%%%%%%%%%%%
A=535.6;
cg=-0.0335;
CGP=cg*(i-1);
CGP1(i)=A+CGP;
g_z=raised*(1+cos(pi*(z-l/2)/l));
sig1=((2*pi*neff)*(1/w)-pi/CGP1(i))+((2*pi*delnbar)/w);
gammab1=((((g_z*pi*delnbar)/w)^2-((sig1)^2)))^.5;
f11=(cosh(gammab1*delz)-sqrt(-1)*(sig1/gammab1)*sinh(gammab1*delz));
f22=conj(f11);
f21=(((sqrt(-1)*g_z*pi*delnbar)/(gammab1*w))*sinh(gammab1*delz));
f12=conj(f21);
f1=f1*[f11 f12;f21 f22];
end
r3(r) = f1(2,1)/f1(1,1); %amplitude reflection coefficient
R3(r) = (abs(r3(r)))^2; %power reflection coefficient of single grating
end
%%%%%%%%%%%%%%%%%%%%%%PLOTS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(1), hold on,
plot(wal, R3, 'r')
grid
xlabel('Wavelength (nm)')
ylabel('measured Reflectivity')
g=sum((R3-R3_2).^2);
144
Liste des figures
Fig (1.1). Evolution de la capacité de transmission d’un lien optique: (a) lien optique point à
point simple; (b) lien multiplexé spatialement en quatre canaux; (c) Lien multiplexé
en quatre longueurs d’onde (6)
Fig (1.2). Atténuation de la lumière dans une fibre de silice en fonction de la longueur d’onde
(8)
Fig (1.3). Définition de quelques paramètres caractérisant un composant optique (10)
Fig (1.4). Réponses spectrales de deux canaux correspondants à un filtre diélectrique à
couches minces (Thin Film Filter) et un réseau de Bragg (Fiber Grating) (11)
Fig (1.5). Délai de groupe d’un coupleur interférométrique fonctionnant à 100 GHz (11)
Fig (1.6). Dispersion d’un réseau de Bragg (FBG) et d’un filtre diélectrique à couches minces
(TFF) (12)
Fig (1.7). Implémentation d’un lien optique multiplexé en longueurs d’onde (13)
Fig (1.8). Structure d’un coupleur en étoile de base (13)
Fig (1.9). Spectres de réflectivité d’un multi couches conçu avec 5, 10 et 15 périodes, centrés
sur la longueur d’onde 1550 nm, avec nB=1.45 et nE=2.1 (15)
Fig (1.10). Réponse spectrale d’un filtre passe bande avec différents nombres de cavités (17)
Fig (1.11). Schéma d’un filtre à couches minces avec trois ports (17)
Fig (1.12). Schéma d’un composant en zigzag (18)
Fig (1.13). Coupleur mono canal à trois ports (conçu par JDSU) (19)
Fig (1.14). Démultiplexeur à quatre canaux basé sur quatre coupleurs cascadés à base de
filtres de couches minces (19)
Fig (1.15). Performance de deux différents coupleurs à division de bande (a) Coupleur 4×1
(b) Coupleur 4×0 (20)
Fig (1.16). Structure d’un démultiplexeur à 40 canaux basé sur cinq filtres 8×0 à division de
bande combinés avec des filtres mono canal (21)
Fig (1.17). Réponse spectrale d’un démultiplexeur à 40 canaux fonctionnant à 100 GHz
(conçu par CIRCA) (22)
Fig (1.18). Principe d’un démultiplexage par diffraction dans un réseau optique : les longueurs
d’onde λ1, λ2, λ3 venants d’une fibre de transmission sont diffractées dans
différentes directions. Le réseau travaille en (a) réflexion et en (b) transmission
(22)
Fig (1.19). Calcul des angles d’ordre de diffraction (23)
Fig (1.20). Configuration d’un démultiplexeur basé sur un réseau de diffraction avec une
correction d’une aberration sphérique causée par un miroir concave sphérique. La
correction est faite avec une combinaison de deux matériaux d’indices n et N
présentant la même valeur de dispersion (Configuration faite par Stimax®) (24)
Fig (1.21). Réponse spectrale de trois canaux adjacents avec un espacement inter-canal de 100
GHz (documentation Stimax®) (24)
Fig (1.22). Structure d’un démultiplexeur à base de réseau de guides d’onde assemblés (25)
Fig (1.23). Représentation de la région de propagation libre (RPL) entrée/sortie (25)
Fig (1.24). Structure d’un interféromètre de Mach-Zehnder (26)
145
Fig (1.25). Concept de base d’un filtre optique accordable (27)
Fig (1.26). Coupleur directionnel multi-électrodes fabriqué sur un crystal LiNbO3 (29)
Fig (1.27). Principe de fonctionnement des filtres multi réseaux accordables (30)
Fig (2.1). Le changement de l’indice de réfraction dans un réseau de Bragg (34)
Fig (2.2). Diagramme illustrant les propriétés d’un réseau de Bragg (35)
Fig (2.3). Méthode de calcul de la réponse spectrale de réseau (40)
Fig (2.4). Diagramme principal de la méthode des matrices de transfert (a) réseau uniforme
(b) réseau à pas variable (42)
Fig (2.5). Réseau de Bragg chirpé linéaire (43)
Fig (2.6). Spectres de réflectivité de deux réseaux chirpés avec des valeurs opposées en signe
de la variable chirp (45)
Fig (2.7). Spectres réfléchis de réseaux de Bragg avec différentes valeurs de la variable chirp
(46)
Fig (2.8). Spectres réfléchis de trois réseaux chirpés linéaires avec différentes longueurs (47)
Fig (2.9). Spectres réfléchis de trois réseaux chirpés avec différents changements d’indice de
réfraction (47)
Fig (2.10). Relation entre la longueur d’onde centrale et la variable chirp pour trois réseaux de
Bragg avec différentes longueurs: L =10000(μm) (ligne solide), L = 15000(μm)
(ligne en tirets), L = 20000(μm) (ligne en points tillés), paramètres physiques: neff
=1.447, δneff = 0.0004 et D = 1550 nm (48)
Fig (2.11). Relation entre la longueur d’onde centrale et le paramètre chirp pour trois réseaux
de Bragg avec différentes longueurs: L =10000(μm) (ligne solide), L =
15000(μm) (ligne en tirets), L = 20000(μm) (ligne en points tillés), paramètres
physiques: neff =1.447, δneff = 0.0004 et D = 1550 nm (49)
Fig (2.12). Rélfectivité d’un réseau uniforme (49)
Fig (2.13). Structure générale d’un réseau apodisé (50)
Fig (2.14). Fonction en cosinus surélevé ( rc =0.5) (50)
Fig (2.15). Fonction Gaussienne ( g =10) (50)
Fig (2.16). Fonction en sinc (71)
Fig (2.17). Fonction de Blackman (B=0.2) (51)
Fig (2.18). Spectres réfléchis d’un réseau chirpé (a) non apodisé (b) apodisé en cosinus surélevé
(53)
Fig (2.19). Représentation du (a) spectre réfléchi et (b) de la dispersion de quatre réseaux
apodisés gaussiens (54)
Fig (2.20). Structure générale d’un réseau de Bragg (a) mono canal (b) multi canal (55)
Fig (2.21). Canaux réfléchis d’un réseau échantillonné fonctionnant à 25 GHz (57)
Fig (2.22). Canaux réfléchis d’un réseau échantillonné fonctionnant à 50 GHz (57)
Fig (2.23). Canaux réfléchis d’un réseau échantillonné fonctionnant à 100 GHz (58)
Fig (2.24). Valeurs des piques de réflectivité et des largeurs de bande d’un réseau
échantillonné fonctionnant à 100 GHz (58)
Fig (2.25). Illustration du décalage périodique de phase dans un réseau échantillonné (59)
Fig (2.26). Canaux réfléchis d’un réseau échantillonné pour Md=1 (60)
Fig (2.27). Canaux réfléchis d’un réseau échantillonné pour Md=2 (61)
146
Fig (2.28). Canaux réfléchis d’un réseau échantillonné pour Md=3 (61)
Fig (2.29). Illustration de la densification des canaux en appliquant un chirp de phase linéaire
(62)
Fig (2.30). Démonstration du décalage des canaux réfléchis par l’ajout d’extra phase (64)
Fig (2.31). Deux exemples expérimentaux illustrant deux méthodes de réglage de la
dispersion (a) Gradient thermique contrôlé par deux refroidisseurs
électrothermiques [II.34] longueur d’onde de Bragg (b) Anneau chauffant
uniforme (66)
Fig (2.32). Gradients de température des refroidisseurs correspondants aux valeurs désirées de
la dispersion (66)
Fig (2.33). Spectre de réflexion d’un réseau ajusté par un dispositif d’échauffement uniforme
avec les différentes valeurs de voltage appliqué: (a) 0 V; (b) 2 V; (c) 4 V; (d) 6 V
(67)
Fig (2.34). Déplacement des canaux réfléchis d’un réseau échantillonné par application de
différents gradients de température (68)
Fig (2.35). Transformation d’une fibre de Bragg en une forme de Ω à l’aide d’un étage fixe et
un autre variable (69)
Fig (2.36). Distribution d’une traction longitudinale appliquée sur une fibre de Bragg avec
différents ΔL (70)
Fig (2.37). Spectres de réflexion d’un réseau chirpé mono canal avec différents ΔL (70)
Fig (2.38). Pli d’une fibre sur une de ses extrémités et création d’un bloc d’adaptation d’indice
sur l’autre extrémité (71)
Fig (2.39). Spectres réfléchis d’un réseau à pas variable pour différentes longueurs inactives
(71)
Fig (2.40). Spectres réfléchis restitués pour différentes longueurs inactives (72)
Fig (2.41). Emballage d’un réseau de Bragg soumis à une traction axiale en utilisant une vis
micromètre (73)
Fig (2.42). Exemple d’une traction appliquée sur le réseau de Bragg emballé suivant (a) Un
étirement allant à une compression, (b) Une compression allant à un étirement
(73)
Fig (2.43). Démonstration d’un déplacement multi canal par application de tractions
uniformes (74)
Fig (2.44). Schéma d’un DEMUX 1×4 (conçu par OptiSystem) (75)
Fig (2.45). Spectres des canaux sortants d’un DEMUX 1×4 (a) λ1=1553.3288 nm, (b)
λ2=1552.5244 nm, (c) λ3=1551.7208 nm (d) λ4=1550.918 nm (Résultats obtenus
par OptiSystem) (77)
Fig (2.46). Structure interne d’un multiplexeur à insertion/extraction (conçu par OptiSystem)
(79)
Fig (2.47). Circuit de test d’un multiplexeur à insertion/extraction (simulé par OptiSystem)
(79)
Fig (2.48). Spectres de sortie issus des différents ports du MIE testé (a) Port d’entrée (b) Port
de sortie (c) Port d’extraction (résultats obtenus par OptiSystem) (81)
Fig (3.1). Classification des métaheuristiques (86)
Fig (3.2). Les cinq niveaux d’organisation d’un algorithme génétique (88)
147
Fig (3.3). Probabilité de chaque individu placé sur une roulette de casino (89)
Fig (3.4). Croisement en un point (91)
Fig (3.5). Croisement en deux points (91)
Fig (3.6). Mutation (92)
Fig (3.7). Les points de réflexion, expansion, contraction et rétrécissement d’un simplexe en
deux dimensions (102)
Fig (3.8). Démonstration d’une recherche de voisinage à partir d’un point appartenant à une
semi-RT (111)
Fig (3.9). Structure principale de la méthode RTD (115)
Fig (3.10). Représentation des paramètres chirps reconstruits (ligne solide) et originaux (ligne
en points tillés) (Résultats obtenus en utilisant l’algorithme génétique (AG) dans
[P5]) (117)
Fig (3.11) Organigramme représentant les étapes de la RTD pour la synthèse d’un réseau de
Bragg mono canal (120)
Fig (3.12) Synthèse d’un réseau de Bragg à pas variable avec apodisation en cosinus surélevé
à partir d’un spectre réfléchi supposé mesuré (122)
Fig (3.13) Comparaison entre la période de réseau reconstruite et originale (122)
Fig (3.14) Synthèse d’un réseau de Bragg multi canal avec une apodisation gaussienne à partir
de canaux réfléchis originaux (125)
Fig (3.15) Evolution de la longueur d’onde de Bragg d’un réseau mono canal synthétisé par la
RTD à partir des spectres en réflexion correspondants à quatre différentes
températures (127)
Fig (3.16) Evolution des paramètres reconstruits, (a) indice effectif de réfraction (b) période
initiale du réseau, (c) chirp de réseau, d’un réseau mono canal synthétisé par la
RTD à partir des spectres en réflexion correspondants à quatre différentes
températures (30 °C, 50 °C, 70 °C et 80 °C) (128)
Fig (3.17) Synthèse d’un réseau de Bragg à pas variable apodisé en cosinus à partir d’un
spectre en réflexion supposé mesuré à 70 °C (129)
Fig (3.18) Spectres réfléchis correspondants à quatre profils de traction (130)
Fig (3.19) (a) Etirement et (b) compression d’un réseau de Bragg à pas variable négatif (131)
Fig (3.20). Profil de traction reconstruit et original (ε(z)= 3z+100) (132)
Fig (3.21). Profil de traction reconstruit et original (ε(z)=z2+3z+100) (132)
Fig (3.22). Profil de traction reconstruit et original (ε(z)=-z2-3z-100) (133)
Fig (3.23). Synthèse d’un réseau de Bragg multi canal avec une traction négative à partir d’un
spectre en réflexion original avec le profil de traction ε(z)=-3z2-50z-100 (134)
Fig (3.24). Profil de traction reconstruit et original pour un réseau échantillonné (ε(z)=-3z2-
50z-100) (135)
148
Liste des tableaux
149
Liste des symboles et abréviations
Δν : Déviation en fréquence
Δλ : Largeur spectrale
λ : Longueur d’onde
AR : Anti Réflexion
HR : Haute Réflexion
R : Coefficient de réflexion
150
GRIN : Lentilles à GRadient d’Indice
θ : angle d’incidence
I0 : Intensité lumineuse
c : Vitesse de lumière
Δf : Bande de fréquences
ßj : Constante de propagation
t
K mn z : Coefficient de couplage transverse entre les modes n et m
ˆ : Coefficient général de couplage propre (DC)
151
k z : Coefficient de couplage (AC)
δ : Syntonisation initiale
v : Visibilité de frange
r : Réflectivité
: Phase du réseau
d D
: Variable chirp
db
F : Paramètre chirp
g : Paramètre de Gauss
B : Paramètre de Blackman
LA : Période d’échantillonnage
δλ : Espacement inter-canal
152
Md : Rapport de densification des canaux
cp : Chirp de phase
e : Extra phase
T : Gradient de température
n : Coefficient thermo-optique
z : Phase de traction
AG : Algorithmes Génétiques
PT : Probabilité de Boltzmann
K : constante de Boltzmann
T : Température
p(X,Y) : Probabilité
T0 : Température initiale
Ns : Taille de voisinage
153
Lm : Nombre d’itérations à température constante
LT : Liste Tabou
ρ: Coefficient de réflexion
χ : Coefficient d’expansion
γ : Coefficient de contraction
σ : Coefficient de rétrécissement
RM : Recherche de Motif
RT : Régions Tabou
154
BIBLIOGRAPHIE
Références du chapitre I
[I.1] G. R. Hill, ‘‘Wavelength domain optical network techniques’’, Proc. IEEE, vol. 78, p.
121 (1990).
[I.2] T. Miki, ‘‘Optical transport networks’’ Proc. IEEE, vol. 81, p. 1594 (1993).
[I.3] K. Nosu, ‘’Optical FDM Network Technologies’’, Artech House, Boston, 1997
[I.4] M. J. O’Mahoney, ‘‘Optical multiplexing in fiber networks: Progress in WDM and
OTDM’’, IEEE Commun Mag, vol. 33, p. 82 (1995).
[I.5] W. J. Goralski, ‘’SONET: A Guide to Synchronous Optical Network’’, McGraw-Hill,
New York, 1997
[I.6] C. A. Siller, Jr., and M. Shafi, eds., ‘’SONET/SDH”, IEEE Press, New York, 1996
[1.7] Gerd E. KeiserA, ‘’Review of WDM Technology and Applications’’, Optical Fiber
Technology, vol. 5, p. 3-39 (1999).
[I.8] G. E. Keiser, ‘’Optical Fiber Communications’’, 2nd. ed., McGraw-Hill, New York,
1991; 3rd. ed., 1999.
[I.9] Telecommunication Standardization Sector of the International Telecommunication
Union _ITU-T., Place des nations, CH-1211 Geneva 20, Switzerland (http//:[Link])
[I.10] ITU-T Recommendation G.692, ‘‘Optical interfaces for multichannel systems with
optical amplifiers’’,1997.
[I.11] A. Srivastava, ‘‘Wide Bandwidth High Capacity Systems’’, OFC’99, paper FC4.
[I.12] A. Srivastava, L. Zhang, Y. Sun, J. W. Sulhoff, and C. Wolf, ‘‘System Margin
Enhancement with Raman Gain in Multi-Span WDM Transmission’’, OFC’99, paper FC2
[I.13] J. Kani, et al., ‘‘Trinal-wavelength-band WDM transmission over dispersion shifted
fibre’’, Electron Lett, vol. 321 (1999).
[I.14] Bruce Nyman, Mark Farries, and Calvin Si, ‘’ Technology Trends in Dense WDM
Demultiplexers’’, Optical Fiber Technology vol. 7, p. 255-274 (2001).
[I.15] A. Hill and D. Payne, ‘‘Linear Crosstalk in Wavelength Division Multiplexed Optical
Fiber Transmission Systems’’, J. Light. Tech, vol. 643 (1985).
[I.16] P. E. Green, Jr, ‘‘Optical networking update’’, IEEE J. Select. Areas Commun, vol. 14,
p. 764 (1996).
[I.17] Markus K. Tilsch, Robert B. Sargent and Charles A. Hulse, ‘’Dielectric Multilayer
Filters’’, white paper
[I.18] H.A. Macleod, ‘’Thin-Film Optical Filters’’, Institute of Physics Publishing, Dirac
House, Temple Back, Bristol BS1 6BE, UK, 2001.
[I.19] A. Thelen, ‘’Design of Interference Coatings’’, McGraw-Hill Book Company, New
York, 1989.
[I.20] P. Baumeister, ‘’Optical Coating Technology’’, SPIE Optical Engineering Press,
Bellingham, WA, 2004.
[I.21] J.D. Rancourt, ‘’Optical Thin Films: User Handbook’’, SPIE Optical Engineering
Press, Bellingham, WA, 1996.
[I.22] P. Yeh, ‘’Optical Waves in Layered Media’’, John Wiley and Sons, Hoboken, NJ, 2005.
157
[I.23] A. Thelen, “Antireflection Coatings, in Design of Interference Coatings’’,Chap. 4
McGraw-Hill Book Company, New York, 1989.
[I.24] J.A. Dobrowolski, A.V. Tikhonravov, M.K. Trubetskov, B. T. Sullivan, and P.G. Verly,
“Optimal single-band normal-incidence antireflection coatings”, Appl. Opt. vol. 35, p. 644–
658 (1996).
[I.25] P. Baumeister, “Reflection reducing coatings, in Optical Coating Technology’’,Chap. 4
SPIE Optical Engineering Press, Bellingham, WA, 2004.
[I.26] U. Schallenberg, U. Schulz, and N. Kaiser: “Multicycle AR coatings: a theoretical
approach”, Proc. SPIE, vol. 5250, p. 357–366 (2004).
[I.27] Stevenson: “High performance antireflection coatings for telecommunications”, Proc.
SPIE, vol. 5527, p. 79–92 (2004).
[I.28] W. J. Tomlinson, “Applications of GRIN-rod lenses in optical fiber communication
systems”, Appl. Opt. vol. 19, p. 1127–1139 (1980).
[I.29] Y.C. Si, G. S. Duck, J. Ip, and N. Teitelbaum, “Narrow band filter and method of
making same”, U.S. Patent, vol. 5, p. 612-824 (March 18, 1997).
[I.30] M.A. Scobey and D. E. Spock, “Passive DWDM components using microplasma
optical interference filters”, Optical Fiber Communication Conference, 1996 (OFC’96), p.
242–243
[I.31] K. Nosu, H. Ishio, and K. Hashimoto, “Multireflection optical multi/demultiplexer
using interference filters”, Electron. Lett, vol. 15, p. 414–415 (1979).
[I.32] M. Scobey, D. E. Spock, M. [Link], and J.A. Beattie, “EDFA gain equalization using
interference films”, NFOEC 1996, Technical Proceedings, Denver, CO, p. 969–972 (1996).
[I.33] M.A. Scobey, W. J. Lekki , and T.W. Geyer, “Filters create thermally stable, passive
multiplexers”, Laser Focus World, vol. 33, p. 111–116 ( 1997).
[I.34] Jean-Pierre Laude, ‘’ Diffraction Gratings WDM Components ‘’, white paper.
[I.35]M. C. Hutley, ‘’Diffraction gratings’’, Academic Press, London, 1982.
[I.36] J. P. Laude and J. Flamand, ‘’Herstellung von Beugungsgittern für Spektrometrie und
Optoelektronik,“ Feinwerktechnik und Meßtechnik, Carl Hanser Verlag, München, vol. 94, p.
5 (1986).
[I.37] J. M. Lerner, J. Flamand, J. P. Laude, G. Passereau, and A. Thevenon, “Diffraction
gratings ruled and holographic – A review”, Proc. SPIE 240, vol. 14, p. 82-88 (1980).
[I.38]J. Cordelle, J. P. Laude, R. Petit, and G. Pieuchard, “Réseaux classiques, réseaux
holographiques. Les questions d’efficacité et leurs conséquences en spectroscopie
instrumentale”, Nouv. Rev. d’Optique Appliquée, vol. 3, p. 149-154 (1970).
[I.39] E. G. Loewen and E. Popov, ‘’Diffraction gratings and applications’’, Marcel Dekker,
New York, Basel, 1997.
[I.40] [Link]: “A multiplexer/coupler with tapered graded-index glass fibers and a grin
rod lens”, Proc. 10th Europ. Conf. Opt. Commun. (ECOC’84), Berlin, Germany, p. 164-165
(1984).
[I.41] J. P. Laude, J. Flamand, J. C. Gautherin, D. Lepere, P. Gacoin, F. Bos, and J. Lerner,
“STIMAX, a grating multiplexer for monomode or multimode fibers”, Proc. 9th Europ. Conf.
Opt. Commun. (ECOC’83), Geneva, Switzerland, p. 417 (1983).
158
[I.42] J. P. Laude, ‘’DWDM Fundamentals, Components, and Applications’’, ISBN: 1-58053-
177-6, Artech House, Boston, London, 2002.
[I.43] Research by Nippon Telegraph and Telephone Corporation, “400-channel arrayed
waveguide grating with 25GHz spacing”, 2002.
[I.44] Hiroaki Yamada, Kazumasa Takada and Seiko Mitachi, “Crosstalk Reduction in a 10
GHz Spacing Arrayed-Waveguide Grating by Phase-Error Compensation”, Journal of
Lightwave Technol, vol. 16, no. 3 (March 1998).
[I.45] K. A. McGreer, “Arrayed Waveguide Gratings for Wavelength Routing”, University of
Manitoba and TRLabs, IEEE commications Magazine (December 1998).
[I.46] ‘’Arrayed Waveguide Gratings, Application and Design’’, book Chapter 2, white paper.
[I.47] C. Huang, H. Luo, S. Xu, and P. Chen, ‘‘Ultra-low loss, temperature-insensitive 16-
channel 100 GHz dense wavelength division multiplexers based on cascaded all-fiber
unbalanced Mach-Zehnder structure’’, in Conference on Optical Fiber Communications,
Tech. Digest Series, paper TuH2 (1999).
[I.48] Y. P. Li, and C. Henry, ‘‘Silicon Optical Bench Technology’’, in Optical Fiber
Telecommunications IIIB, I. Kaminow and T. Koch, Eds., Academic, NY, 1997.
[I.49] X. Charles et al ,’’Diffraction-Grating-Based (De)Multiplexer Using Image Plane
Transformations’’, IEEE Journal of selected topics in quantum electronics, vol. 8,
(november/december 2002).
[I.50] H. Kobrinski and K.-W. Cheung, ‘‘Wavelength tunable optical filters: Applications and
technology’’, IEEE Commun. Mag, vol. 27, p. 53 (1989).
[I.51] D. Brooks and S. Ruschin, ‘‘Integrated electro-optic multielectrode tunable filter’’, J.
Lightwave Technol, vol. 13, p. 1508 (1995).
[I.52] E. L. Wooten, R. L. Stone, E. W. Miles, and E. M. Bradley, ‘‘Rapidly tunable
narrowband wavelength filter using LiNbO3 unbalanced Mach]Zehnder interferometers’’, J.
Lightwave Technol, vol. 14, p. 2530 (1996).
[I.53] K. Oda, N. Yakato, T. Kominato, and H. Toba, ‘‘A 16-channel frequency selection
switch for optical FDM distribution systems’’, IEEE J. Selected Areas Commun, vol. 8, p.
1132 (1990).
[I.54] J. Stone and L. W. Stulz, ‘‘High-performance fibre Fabry-Pérot filters’’, Electron. Lett,
vol. 27, p. 2239 (1991).
[I.55] M. Born and E. Wolf, ‘’Principles of Optics’’, 6th ed., Pergamon, New York, 1980.
[I.56] M. Zirngibl, C. H. Joyner, and B. Glance, ‘‘Digitally tunable channel-dropping
filterrequalizer based on waveguide grating router and optical amplifier integration’’, IEEE
Photon. Technol. Lett, vol. 6, p. 513 (1994).
[I.57] D. Y. Al-Salameh, M. T. Fatehi, W. J. Gartner, S. Lumish, B. L. Nelson, and K. K.
Raychaudhuri, ‘‘Optical networking,’’ Bell Labs Technol. J., vol. 3, p. 39 (1998).
[I.58] O. Ishida, H. Takahashi, and Y. Inoue, ‘‘Digitally tunable optical filters using AWG
multiplexers and optical switches’’, J. Lightwave Technol, vol. 15, p. 321 (1997).
[I.59] A. Sneh and K. M. Johnson, ‘‘High-speed tunable liquid crystal filter for WDM
networks’’, J. Lightwa¨e Technol, vol. 14, p. 1067 (1996).
159
[I.60] P.-L. Chen, K.-C. Lin, W.-C. Chuang, Y.-C. Tzeng, K.-Y. Lee, and W.-Y. Lee,
‘‘Analysis of a liquid crystal Fabry-Perot etalon filter: A novel model’’, IEEE Photon.
Technol. Lett, vol. 9, p. 467 (1997).
[I.61] J.-P. Weber, B. Stoltz, H. Sano, M. Dasler, O. O¨berg, and J. Walz, ‘‘An integratable
polarizationindependent tunable filter for WDM systems: The multigrating filter’’, J.
Lightwave Technol, vol. 14, p. 2719 (1996).
[I.62] D. A. Smith, J. E. Baran, J. J. Johnson, and K.-W. Cheung, ‘‘Integrated-optic
acoustically tunable filters for WDM networks’’, IEEE J. Selected Areas Commun, vol. 8, p.
1151 (1990).
[I.63]T. Nakazawa, M. Doi, S. Taniguchi, Y. Takasu, and M. Seino, ‘‘Ti:LiNbO3 AOTF for
0.8 nm channel-spaced WDM systems’’, OFC ’98 (1998).
160
Références du chapitre II
[II.1] J. ZHAO, ‘‘An Object-oriented Simulation Program for Fibre Bragg Gratings’’, Master
of engineering, Rand Afrikaans University of Johannesburg, South Africa, 2001.
[II.2] C. R. Giles, “Lightwave applications of fibre Bragg gratings”, Journal of Lightwave
Technology, vol. 15, p. 1391-1404 (1997).
[II.3] C. R. Giles and V. Mizrahi, “Low-loss add/drop multiplexer for WDM lightwave
network”, in Tenth International Conference on Integrated Optics and Optical Fiber
Communication, IOOC95 Technical Digest, paper ThC2-1, p. 66-67 (1995).
[II.4] K.O. Hill, Y. Fujii, D.C. Johnson, and B.S. Kawasaki, “Photosensitivity in optical fibre
waveguides: application to reflection filter fabrication”, Applied Physics Letters, vol. 32, p.
647-649 (1978).
[II.5] G. Meltz, W. W. Morey, and W. H. Glenn, “Formation of Bragg gratings in optical
fibres by a transverse holographic method”, Optics Letters, vol. 14, p.823-825 (1989).
[II.6] A. Othonos and K. Kalli, “Fibre Bragg gratings: fundamentals and applications in
telecommunications and sensing”, Artech House, 1999.
[II.7] C. R. Giles “Lightwave application of fiber Bragg gratings”, Journal of Lightwave
Technology, vol.15, p. 1391-1404 (1997).
[II.8] A. W. Snyder and J. D. Love, “Optical waveguide theory”, Chapman and Hall, London,
p. 542, 1983.
[II.9] T. Erdogan, “Fibre grating spectra”, Journal of Lightwave Technology, vol.15, p. 1277-
1294 (1997).
[II.10] A. D. Kersey, M. A. Davis, H. J. Patrick, M. Leblanc, K. P. Koo, C. G. Askins, M. A.
Putnam, and E. J. Friebele, “Fibre grating sensors”, Journal of Lightwave Technology, vol.15,
, p. 1442-1463 (1997).
[II.11] R. Kashyap, “Photosensitive optical fibres: devices and applications”, Optical Fibre
Technology, p.17-34 (1994).
[II.12] S. J. Spammer, P. L. Swart, and A. A. Chtcherbakov, “Merged Sagnac- Michelson
interferometer for distributed disturbance detection”, Journal of Lightwave Technology,
vol.15, p. 972-976 (1997).
[II.13] G. A. Ball, W. W. Morey, and W. H. Glenn, “Standing-wave monomode erbium fiber
laser”, IEEE Photonics Technology Letters, vol. 3, p. 613- 615 (1991).
[II.14] A. Hamakawa, T. Kato, G. Sasaki, and M. Higehara, “Wavelength stabilization of 1.48
um pump laser by fiber grating”, in Proc. ECOC’96, Oslo, Norway, paper MoC3.6 (1996).
[II.15] P. C. Becker, N. A. Olsson, J. R. Simpson, P. C. Becker, and P. Becker, “Erbium-
doped fiber amplifiers : fundamentals and technology’’, Optics and Photonics series,
Academic Press, 1999.
[II.16] J. A. R. Williams, I. Bennion, K. Sugden, and N. J. Doran, “Fiber dispersion
compensation using a chirped in-fire Bragg grating”, Electronics Letters, vol. 30, p. 985-987
(1994).
[II.17] C. R. Giles and J. M. P Delavaux, “Repeaterless bidirectional transmission of 10 Gb/s
WDM channels”, in ECOC’95, Brussels, paper PD2 (1995).
161
[II.18] R. Kashyap, R. Wyatt, and R. J. Campbell, “Wideband gain flattened erbium fiber
amplifier using a photosensitive fiber blazed grating”, Electronics Letters, vol. 29, p. 154-156
(1993).
[II.19] C. R. Giles and V. Mizrahi, “Low-loss add/drop multiplexers for WDM lightwave
networks”, in Proc. IOOC’95, Hong Kong, paper ThC2-1 (1995).
[II.20] B. H. Lee, Y. Chung, and U. Paek, "Fiber comb filters based on fiber gratings",
COOC99, 5F/13/1999, p. 19-20 (1999).
[II.21] N. G. R. Broderick, D. Taverner, D. J. Richardson, M. Isben, and R. I. Laming,
“Optical pulse compression in fibre Bragg gratings”, Physical Review Letters, vol.79, p.
4566-4569 (1997).
[II.22] Y. Tohmori, F. Kano, H. Ishii, Y. Yoshikuni, and Y Kondo, "Wide tuning with narrow
linewidth in DFB lasers with superstructure grating (SSG)", Electronic Letters, vol.29, p.
1350-1351 (1993).
[II.23] M. Prabhu, N. S. Kim, L. Jianren, J. Xu, K. Ueda, "Highly-efficient ultrabroadband
supercontinuum generation centered at 1484nm using Raman fiber laser", Photonics West,
LASE2001, San Jose, USA, 2001.
[II.24] A. Boskovic, M. J. Guy, S. V. Chernikov, J. R. Taylor And R. Kashyap “All fiber
diode-pumped, femtosecond chirped pulse amplification system”, Electronics Letters, vol. 31,
, p. 877-879 (1995).
[II.25] «Optimisation, écriture et démonstration d’encodeurs à réseaux de Bragg multiples
pour le CDMA » Hafedh Ben Jaâfar, thèse de doctorat, Université de Laval, 2001.
[II.26] Y. Chen and S. Jian, “An introduction to lightwave technology”, China Railway
Publishing, p. 248 (2000).
[II.27] M. Yamada and K. Sakuda, "Analysis of almost-periodic distributed feedback slab
waveguide via a fundamental matrix approach", Applied Optics, vol. 26, p. 3474-3478 (1987).
[II.28] A. Othonos, “Fibre Bragg gratings”, Review of Scientific Instruments, vol. 68, p.
4309-4341 (1997).
[II.29] G. Booch, “Object-oriented development”, IEEE Transactions on Software
Engineering, vol. 12, p. 211 (1986).
[II.30] P. C. Becker, N. A. Olsson, J. R. Simpson, P. C. Becker, and P. Becker, “Erbium-
doped fiber amplifiers : fundamentals and technology’’, Optics and Photonics series,
Academic Press, 1999.
[II.31] SUN, C-c, LI, C-y, X-n, YU, ‘’Study on sampled chirped fiber gratings’’, Optics
communication, vol. 218, p. 297-302 (2003).
[II.32] I. Narvuz, N. F. Guler, ”A novel technique for optical dense comb filters using
sampled fiber Bragg gratings”, Optical Fiber Technology, vol. 14, p. 114-118 (2008).
[II.33] N. Mohammad “ Analysis and Development of a Tunable Fiber Bragg Grating Filter
Based on Axial Tension/Compression”, Master thesis, University of Saskatchewan,
Saskatoon, Canada, 2005.
[II.34] A. F. Abas,’’Chromatic Dispersion Compensation in 40 Gbaud Optical Fiber WDM
Phase-Shift-Keyed Communication Systems’’, Master thesis, Paderborn, Germany , 2006.
162
[II.35] Z. Zhao et al ‘’Adjustment of the central wavelength and the chirp of fiber Bragg
grating separately with external heaters’’, Optics Communications, vol. 242, p. 135–139
(2004).
[II.36] J. Zhou et al ‘’Tunable chirped fiber Bragg grating based on the D-shaped fiber’’,
Optics Communications, vol. 281, p. 2077–2082 (2008).
[II.37] J. Kwon et al ‘’Tunable dispersion compensation with fixed center wavelength and
bandwidth using a side-polished linearly chirped fiber Bragg grating’’, Optical Fiber
Technology, vol. 11, p. 159–166 (2005).
[II.38] Monica L. Rocha et al ‘’Mechanical tuning of a Fiber Bragg grating for optical
networking applications’’, Journal of Microwaves and Optoelectronics, vol. 4 (June 2005).
[II.39] Davies PA, Abd-El-Hamid G, ‘’Four-port fibre-optic ring resonator’’, Electron Lett
vol. 24, p. 662–663 (1988).
[II.40] Blow KJ, Doran NJ, Nelson BP, ‘’Demostration of the nonlinear fibre loop mirror as
an ultrafast alloptical demultiplexer’’, Electron Lett, vol. 26, p. 962–964 (1990).
[II.41] Mizrahi V, Erdogan R, DiGiovanni DJ, Lemaire PJ, MacDonald WM, Kosiski SG,
Cabot S, Sipe JE, ‘’Four channel fibre grating demultiplexer’’, Electron Lett, vol. 30, p. 780–
781 (1994).
[II.42] Pan JJ, Shi Y, ‘’Dense WDM multiplexer and demultiplexer with 0.4nm channel
spacing’’, Electron Lett, vol. 34, p. 74–75 (1998).
[II.43] R. Romero et al , ‘’Multiplexers and Demultiplexers Based on Fibre Bragg Gratings
and Optical Circulators for DWDM Systems’’, M.M. Freire, P. Lorenz, M.M.-O. Lee (Eds.):
HSNMC 2003, LNCS 2720, Springer-Verlag Berlin Heidelberg, p. 442–451 (2003).
[II.44] R. Romero et al, ‘’Chirped fibre Bragg grating based multiplexer and demultiplexer for
DWDM applications’’, Optics and Lasers in Engineering, vol. 43, p. 987–994 (2005)
[II.45] Orlov S. S, Yariv A, and Van Essen S, "Coupled-mode analysis of fiber-optic add-drop
filters for dense wavelength division multiplexing", Opt. Lett, vol. 22, p. 688-690 (1997).
[II.46] Biloudeau F, Johnson D. C, Theriault S, Malo B, Albert J, and Hill K. O, "An all-fiber
densewavelength- division multiplexer/demultiplexer using photoimprinted Bragg gratings",
IEEE Phot. Tech. Lett ,vol. 7, p. 388-390 (1995).
[II.47] Dong L, Hua P, Birks T. A, Reekie L, and Russell P. St. J, "Novel add/drop filters for
wavelengthdivision- multiplexing optical fiber systems using a Bragg grating assisted
mismatched coupler", IEEE Phot. Tech. Lett, vol. 8, p. 1656-1658 (1996).
[II.48] Park K. N, Lee Y. T, Kim M. H, Lee K. S, and Won Y. H, "All-fiber drop-pass with
fiber Bragg gratings", IEEE Phot. Tech. Lett, vol. 10, p. 555-557 (1998).
[II.49] Augustsson T, "Bragg grating-assisted MMI -coupler for add-drop multiplexing",
Journal of Lightwave Technology, vol. 16, p. 1517-1522 (1998).
163
Références du chapitre III
[III.1] Abdesslem layeb ,’’ Introduction aux métaheuristiques ‘’, Chapitre de thèse (2009).
[III.2] Selvaraj Ramkumar, ‘’Algorithmes génétiques’’, article (2007).
[III.3] S. Amédée et R. Francois-Gérard, ‘’Algorithme génétiques’’ TE de fin d’année,
Tutorat de Mr Philippe Audebaud, 2004.
[III.4] Yann Ollivier, ‘’ Introduction aux algorithmes génétiques’’, article (Octobre 2000).
[III.5] Thomas Vallée et Murat Yıldızoglu, ‘’ Présentation des algorithmes génétiques et de
leurs applications en économie’’, article (décembre 2003).
[III.6] Jean-Marc Alliot et Nicolas Durand, ‘’Algorithmes génétiques’’, article (Mars 2005).
[III.7] Omessaad Hajji, ‘’Contribution au développement des méthodes stochastiques.
Application à la conception des dispositifs électrotechniques’’, Thèse de Doctorat, Université
des Sciences et Technologies de Lille, 2003.
[III.8] S. Kirkpatrick et al, ‘’Optimization by Simulated Annealing’’, Science, vol. 220, p.
671-680 (1983).
[III.9] R.E. Randelman et G.S. Grest, ‘’N-City travelling Salesman Problem: Optimization by
Simulated Annealing’’, Journal of Optimization Theory and Application, vol. 4, p. 303-320
(1969).
[III.10] E. Aarts and J. Korst, ‘’Simulated Annealing and Boltzmann Machines’’, John
Wiley&Sons, New York, 1990.
[III.11] N. Metropolis et al, ‘’Equation of state calculations by fast computing machines’’,
journal Chem. Phys, vol. 27, p. 108761092 (1953).
[III.12] Laurent Miclo, ‘’Une étude des algorithmes de recuit simulé sous-admissibles’’,
Annales de la faculté des sciences de Toulouse, 6ème série, tome 4, No4 , p. 819-877 (1995).
[III.13] Sébastien Loustau, ‘’Le recuit siumlé : analyse théorique et applications’’ Projet 2A,
2005.
[III.14] D. Vanderbilt and S.G. louie, ‘’A Monte Carlo Simulated annealing Approch to
Optimization over Continuous Variables’’, Journal of Computational physics, vol. 56, p. 259-
271 (1984).
[III.15] L. Painton and U.M. Diwekar, ‘‘Synthesizing Optimal Design Configurations for a
Brayton Cycle Power Plant”, Comput. Chem. Eng, vol 18, p. 369-381 (1994).
[III.16] W. B. Dolan et al, ‘’Algorithmic efficiency of simulated annealing for heat exchanger
network design”, Comp. Chem. Eng, vol. 14, p. 1039-1050 (1990).
[III.17] F. Glover, ‘’Future paths for integer programming and links to artificial intelligence’’,
Computers & Operations Research, vol. 13, p. 533-549 (1986).
[III.18] D.E. Goldberg, ‘’Genetic algorithms in Search, Optimization and Machine
Learning’’, Addison Wesley, 1989.
[III.19] D.E. Goldberg, ‘’Sizing Populations for Serial and Parallel Genetic Algorithms’’, Int.
Conf. genetic algorithm (1989).
[III.20]F. Glover and M. Laguna, ‘’Tabu Search’’, Colin R. Reeves (ed.), Modern Heuristic
Techniques for combinatorial Problems, Blackwell Scientific Publications, Oxford, p. 70-150
(1993).
164
[III.21] S. Tsubakitani and J. R. Evans, ‘’Optimizing tabu list size for travelling salesman
problem’’, Computers Ops. Res, vol. 25, p. 91-97 (1998).
[III.22] M. Gendreau and G. Laporte, ‘’A Tabu Search Heuristic for The heterogenous Fleet
Vehicle Routing Problem’’, Computers & Operations Research, vol. 26, p.1153-1173 (1999).
[III.23] M. Widmer and A. Hertz, ‘’A new approach for solving the flow shop sequencing
problem’’, Europ. J. Oper. Res, vol 41, p. 186-193 (1989).
[III.24] C. C. Shyur and T.C. Lu, ‘’Applying Tabu Search to Spare Capacity Planning for
Network Restoration’’, Computers & Oper. Res, vol. 26, p.1175-1194 (1999).
[III.25] F. Glover, ‘’Tabu Search-Part II’’, ORSA Journal on Computing, vol 2, p. 4-32
(1990).
[III.26] Michel Gendreau, ‘’ An Introduction to Tabu Search’’, article sur
[Link] juillet 2002.
[III.27] N . Hu, “ Tabu Search Method With Random moves for Globally Optimal’’,
International journal for Numerical Methods in Engineering, vol. 35, p. 1055-1070 (1992).
[III.28] Sébastien Noël, ‘’ Mathématiques hybrides pour la résolution des systèmes
d’ordonnoncement de voitures dans une chaîne d’assemblage automatique’’, mémoire de
maitrise en informatique, Université du Québeq, Montréal, 2007.
[III.29] Talbi, E.-G, ‘’A Taxonomy ofHybrid Metaheuristics’’, Journal of Heuristics, vol. 8, p.
541-564 (2002).
[III.30] Martin, O.c. et S.W. Otto, ‘’Combining Simulated Annealing with Local Search
Heuristics’’, Annals of Operations Research, vol. 63, p. 57-75 (1996).
[III.31] Salcedo-Sanz, S., Y. Xu et X. Yao, Hybrid Meta-Heuristics Algorithmsfor Task
Assignment in Heterogeneous Computing Systems. Computers and Operations Research, vol.
33, p. 820-835 (2006).
[III.32] Fleurent, C et al, ‘’Genetic Hybridsfor the Quadratic Assignment Problem’’,
DlMACS Series in Mathematics and Theoretical Computer Science, vol. 16, p. 190-206
(1994).
[III.33] Azimi, Z.N, ‘’Bybrid Beuristicsfor Examination Timetabling Problem’’, Applied
Mathematics and Computation, vol. 163, p. 705-733 (2005).
[III.34] Tanese, R , ’’Parallel genetic Algrithms for a Hypercube’’, in Proceedings of the
Second International Conference on Genetic Algorithms, , MIT, Cambridge, MA, USA, p.
177-183 (1987).
[III.35] Hammami, M. et K. Ghédira, ‘’COSATS: A new Cooperation Model between
Simulated Annealing and Tabu Searchfor the K-Graph Partitioning Problem’’, in The Fourth
IEEE International Workshop on Soft Computing as Transdisciplinary Science and
Technology, Mororan, Japan, p. 863-873 (2005).
[III.36] Taillard, É.D, ‘’Parallel Iterative Search Methods for Vehicle Routing Problems.
Networks’’, vol. 23, p. 661-673 (1993).
[III.37] Bachelet, V., Z. Hafidi, P. Preux et E.-G. Talbi, ‘’Diversitying Tabu Search by
Genetic Algorithms’’, in INFORMS'98 on Operations Research and Management Sciences
meeting, , Montréal, Canada (1998).
165
[III.38]Krueger. M, ‘’Méthodes d'analyse d'algorithmes d'optimisation stochastiques à l'aide
d'algorithmes génétiques’’, École Nationale Supérieure des Télécommunications, Paris,
France, 1993
[III.39] Shahookar, K. et P. Mazumder, ‘’A Genetic Approach to Standard Cell Placement
Using Meta-Genetic Parameter Optimization’’, IEEE Transactions on Computer Aided
Design, vol. 9, p. 500-511 (1990).
[III.40] Salcedo-Sanz, S., Y. Xu et X. Yao, ‘’Hybrid Meta-Heuristics Algorithms for Task
Assignment in Heterogeneous Computing Systems’’, Computers and Operations Research,
vol. 33, p. 820-835 (2006).
[III.41] J. E. Baker, ‘’Adaptive selection methods for genetic algorithms’’, In: J. J. Grefen-
stette (Ed.), Proceedings of the First International Conference on Genetic Algorithms,
Lawrence Erlbaum Associates, Hillsdale, MA, pp. 101-111 (1985).
[III.42] C. T. Kelley, ‘’Detection and remediation of stagnation in the Nelder-Mead algorithm
using a succinct decrease condition’’, SIAM J. Optim, vol. 10, p. 43-55 (1999).
[III.43] C. T. Kelley, ‘’Iterative Methods for Optimization’’, Frontiers Appl. Math. 18, SIAM,
Philadelphia, PA, 1999.
[III.44] A. Hedar, M. Fukushima,’’ Minimizing multimodal functions by simplex coding
genetic algorithm’’, Optimization Methods and Software, vol. 18 265–282 (2003).
[III.45] K. I. M. McKinnon, ‘’Convergence of the Nelder-Mead simplex method to a non
stationary point’’, SIAM J. Optim, vol. 9, p. 148-158 (1999).
[III.46] T.G. Kolda, R.M. Lewis and V. Torczon, ‘’Optimization by direct search: New
perspectives on some classical and modern methods’’, SIAM Review, vol. 45, p. 385–482
(2003).
[III.47] V. Torczon,’’ On the convergence of pattern search algorithms’’, SIAM J. Optim, vol.
7, p. 1–25 (1997).
[III.48] A. Hedar, M. Fukushima, ‘’Heuristic pattern search and its hybridization with
simulated annealing for nonlinear global optimization’’, Optimization Methods and Software,
vol. 19, p. 291–308 (2004).
[III.49] R. Chelouah, P. Siarry, ‘’Tabu search applied to global optimization’’, European
Journal of Operational Research, vol. 123, p. 256-270 (2000).
[III.50] F. Franze, N. Speciale, ‘’A tabu-search-based algorithm for continuous multiminima
problems’’, International Journal for Numerical Engineering, vol. 50, p. 665-680 (2001).
[III.51] J. A. Nelder, R. Mead, ‘’A simplex method for function minimization’’, The
Computer Journal, vol. 7, p. 308-313 (1965).
[III.52] F. Glover, M. Laguna, ‘’Tabu Search, Kluwer’’, Academic Publishers, MA, USA,
1997.
[III.53] A. Hedar, M. Fukushima, ‘’Tabu Search directed by direct search methods for
nonlinear global optimization’’, Optimization Methods and Software, vol. 170, p. 329–349
(2006).
[III.54] K. A. Winick and J. E. Roman, ‘’Design of corrugated waveguide filters by Fourier-
transform techniques’’, IEEE J. Quantum Electron, vol. 26, p. 1918–1929 (1990).
166
[III.55] E. Peral, J. Capmany, and J. Marti, ‘’Iterative solution to the Gel’Fand–Levitan–
Marchenko coupled equations and application to synthesis of fiber gratings’’, IEEE J.
Quantum Electron, vol. 32, p. 2078–2084 (1996).
[III.56] F. Lhommé, C. Caucheteur, K. Chah, M. Blondel, and P. Mégret, ‘’Synthesis of fiber
Bragg grating parameters from experimental reflectivity: a simplex approach and its
application to the determination of temperature-dependent properties’’, Applied Optics, vol.
44, p. 2078–2084 (2005).
[III.57] J. Skaar, L. Wang, and T. Erdogan, ‘’On the synthesis of fiber Bragg gratings by layer
peeling’’, IEEE J. Quantum Electron, vol. 37, p. 165–173 (2001).
[III.58] F. Casagrande, P. Crespi, A. M. Grassi, A. Lulli, R. P. Kenny, and M. P. Whelan,
‘’From the reflected spectrum to the properties of a fiber Bragg grating: a genetic algorithm
approach with application to distributed strain sensing’’, Applied Optics, vol. 41, p. 5238–
5244 (2002).
[III.59] P. Dong, J. Azana, and A. G. Kirk, ‘’Synthesis of fiber Bragg grating parameters from
reflectivity by means of a simulated annealing algorithm’’, Optics Communications, vol. 228
p. 303–308 (2003).
[III.60] I. Navruz, N. F. Guler, Optimization of reflection spectra for phase-only sampled fiber
Bragg gratings, Optics Communications, vol. 271, p. 119–123 (2007).
[III.61] R. T. Zheng, N. Q. Ngo, P. Shum, S.C. Tjin, and L. N. Binh, ‘’A staged continuous
tabu search algorithm for the global optimization and its applications to the design of fiber
Bragg gratings’’, Computational Optimization and Applications, vol 30, p. 319–335 (2005).
[III.62] N. Q. Ngo, R. T. Zheng, J. H. Ng, S.C. Tjin, and L. N. Binh, ‘’Optimization of fiber
Bragg grating using a hybrid optimization algorithm’’, Journal of Lightwave Technology, vol.
25, p. 799–802 (2007).
[III.63] S. Y. Huang, M. M. Ohn, M. LeBlanc, and R. M. Measures, ‘’Continuous arbitrary
strain profile measurements with fiber Bragg gratings’’, Smart Mater, Struct, vol 7, p. 248–
256 (1998).
[III.64] S. Yin, “Distributed fiber optic sensors,” in Fiber Optic Sensors, F. T. S. Yu and S.
Yin, eds. Marcel Dekker, New York, p. 183–233, 2002.
167
Publications de l’auteur
[P.1] F. Karim et O. Seddiki ,’’ Numerical analysis of raised cosine sampled chirped Bragg
grating for dispersion compensation in dense wavelength division multiplexing systems’’,
International Journal of Commnunications, Issue 1, volume 3, p. 9-16, 2009. ISSN 1998-
4480.
[P.2] F. Karim et O. Seddiki,’’ Study on raised-cosine sampled chirped fiber Bragg grating
for dispersion compensation applications’’, Computers and Simulation in Modern Science,
WSEAS press, p. 97-99, 2009. ISSN 1790-2769.
[P.3] F. Karim et O. Seddiki, ‘’ Synthesis of chirped apodized fiber Bragg grating parameters
using direct tabu search algorithm: Application to the determination of thermo-optic and
thermal expansion coefficients’’, Optics Communications, volume 283, p. 2109–2116, 2010.
ISSN 0030-4018.
[P.4] F. Karim et O. Seddiki,’’ Theoretical analysis of sampled Bragg grating with a chirp in
sampling and grating period for dispersion compensation applications’’, Journal of Optical
Communications, volume 31, p. 18-21, 2010. ISSN 0173-4911.
[P.5] F. Karim et O. Seddiki,’’ Reconstruction of raised cosine sampled grating period chirp
distribution using genetic algorithm’’, Journal of Optical Communications, volume 31, p. 14-
17, 2010. ISSN 0173-4911.
[P.6] F. Karim et O. Seddiki,’’ Direct tabu search algorithm for the fiber Bragg grating
distributed strain sensing’’, Journal of Optics, volume 12, 095401 (8 pages), 2010. ISSN
2040-8978.
156