Brevet FR3054034A1 — publié le 20180119

- 1 - «Procédé et dispositif de caractérisation d'une source optique» Domaine technique La présente invention concerne un procédé de caractérisation d'une source optique. Elle concerne aussi un dispositif de caractérisation d'une source optique, et un système comprenant un tel dispositif. Un tel dispositif permet à un utilisateur de caractériser une source optique, par exemple une longueur d'onde d'un rayonnement émis par cette source optique. Etat de la technique antérieure On connaît différents dispositifs selon l'état de l'art pour mesurer une longueur d'onde d'un rayonnement émis par une source optique. Ces dispositifs nécessitent d'explorer une plage de longueurs d'onde possibles afin de déterminer quelle est la longueur d'onde, parmi celles de la plage, qui est émise par la source. Cette exploration peut se faire typiquement : - avec une solution coûteuse, et exploration par barrette de photodiodes. Des solutions existantes proposent des configurations de Fabry-Perot sans pièces mécaniques mobiles, mais avec des barrettes comme détecteur. L'aspect sans pièces mobiles est un avantage certain pour la robustesse et la fiabilité. Par contre, l'usage de barrettes de détection est un inconvénient au niveau du coût, surtout lorsqu'il s'agit d'aller mesurer des longueurs d'onde au-delà de 1100 nm et que le silicium n'est plus utilisable. Il faut alors s'orienter vers des matériaux plus onéreux comme l'InGaAs pour les longueurs d'onde télécom ou encore du PbSe ou MCT pour le moyen infrarouge. Les barrettes coûtent typiquement plusieurs milliers d'euros. - avec une solution lente, et exploration par Fabry Perot accordable. Des solutions à Fabry-Perot accordables n'utilisent pas de barrettes mais seulement une ou plusieurs photodiodes. Cependant, pour réaliser l'exploration souhaitée et obtenir 3054034 -2- l'information souhaitée (position du pic de transmission du FabryPerot), il faut déplacer un des miroirs du Fabry Perrot. Ce balayage ralentit fortement la mesure. Il faut en effet faire une acquisition pour chaque position du miroir, au lieu d'acquérir tout 5 l'interférogramme simultanément dans le cas de l'utilisation de détecteurs barrettes. Pour un temps de réponse typique de 0,1 ms et une mesure de 100 positions de miroirs le temps minimal d'une mesure est alors de 10 ms limitant le taux de rafraîchissement à 100 Hz (sans compter les temps de traitement) 10 Le but de la présente invention est de proposer un procédé et/ou un dispositif de caractérisation d'une source optique qui permette une caractérisation plus économique et/ou plus rapide que l'état de l'art. 15 Exposé de l'invention Cet objectif est atteint avec un procédé de caractérisation d'une source optique, ladite source optique comprenant une cavité fixe ayant un intervalle spectral libre, ledit procédé comprenant : - une génération d'un premier rayonnement par la source selon un 20 état de la source ou différents états explorés de la source, - une réception d'au moins une partie de ce premier rayonnement par au moins un capteur, chaque capteur comprenant un étalon Fabry Perot suivi d'une photodiode, - une mesure d'un signal par chaque capteur et pour chaque état 25 exploré de la source, en fonction de l'au moins une partie du premier rayonnement reçu par ce capteur à cet état exploré de la source, - à partir des signaux mesurés et pour chaque état exploré de la source, un calcul d'une première donnée représentative de la 30 longueur d'onde du premier rayonnement, le calcul comprenant, pour chaque état exploré de la source, une sélection d'une valeur sélectionnée de la première donnée parmi plusieurs valeurs possibles, ladite sélection comprenant une élimination des valeurs de la première donnée ne correspondant pas à une constante 3054034 - 3 - modulo l'intervalle spectral libre de la cavité fixe exprimé selon les unités de la première donnée. Chaque capteur ne comprend de préférence pas de barrette de photodiodes. 5 L'étalon Fabry Perot de chaque capteur ne comprend de préférence aucune partie mobile. La première donnée illustrant la longueur d'onde du premier rayonnement peut être : - la longueur d'onde du premier rayonnement, ou 10 - la fréquence temporelle du premier rayonnement, ou - la période temporelle du premier rayonnement, ou - le nombre d'onde du premier rayonnement. La sélection peut éliminer en outre les valeurs de la première donnée situées en dehors d'un intervalle de zone spectrale du premier 15 rayonnement. Dans ce cas : - chaque étalon Fabry Perot peut avoir une courbe de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, ayant de préférence une période supérieure ou égale à l'intervalle de zone spectrale du premier 20 rayonnement, et/ou - la réception d'au moins une partie du premier rayonnement comprend de préférence une réception par plusieurs capteurs comprenant chacun un étalon Fabry Perot suivi d'une photodiode, chaque étalon Fabry Perot ayant une courbe de transmission de 25 l'intensité du premier rayonnement en fonction de la longueur d'onde du premier rayonnement, l'assemblage des courbes de transmission des différents étalons Fabry Perot formant une courbe globale de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier 30 rayonnement, ayant de préférence une période supérieure ou égale à l'intervalle de zone spectrale du premier rayonnement. L'au moins un capteur comprend de préférence plusieurs capteurs, les différents étalons Fabry Perot des différents capteurs comprenant de préférence aux moins un couple de deux étalons Fabry Perot dont : 3054034 -4- - un premier étalon Fabry Perot ayant une première courbe de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, et - un deuxième étalon Fabry Perot ayant une deuxième courbe de 5 transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, de sorte que : - une même longueur d'onde du premier rayonnement correspond à la fois à un maximum de pente de la première courbe de 10 transmission et à une pente de la deuxième courbe de transmission inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente, et/ou - une même longueur d'onde du premier rayonnement correspond à la fois à un maximum de pente de la deuxième courbe de 15 transmission et à une pente de la première courbe de transmission inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente. Le procédé selon l'invention peut comprendre, avant la sélection éliminant les valeurs de la première donnée ne correspondant pas à la 20 constante modulo l'intervalle spectral libre de la cavité fixe, une détermination ou une connaissance de la constante. La détermination de la constante peut comprendre, dans le cas d'une génération du premier rayonnement par la source selon différents états explorés de la source : - un calcul préliminaire de la première donnée pour chaque état 25 exploré de la source, et - à partir du premier calcul préliminaire de la première donnée pour l'ensemble de différents états explorés de la source, une détermination de la constante. La sélection d'une valeur sélectionnée de la première donnée peut 30 comprendre, après l'élimination des valeurs de la première donnée ne correspondant pas à une constante modulo l'intervalle spectral libre de la cavité, une sélection finale de la valeur sélectionnée, de préférence par une méthode des moindres carrés. La source peut être un oscillateur paramétrique optique ou un laser. 3054034 - 5 - La source peut être une source au moins doublement résonnante, comprenant deux cavités partageant un même cristal de génération de rayonnement, dont: 5 - une cavité mobile générant le premier rayonnement, et - la cavité fixe générant un deuxième rayonnement. Le procédé selon l'invention peut comprendre en outre un calcul, pour chaque état exploré de la source, d'une deuxième donnée représentative de la longueur d'onde du deuxième rayonnement en fonction du calcul de la 10 première donnée représentative de la longueur d'onde du premier rayonnement. Selon une autre particularité de l'invention, il est proposé un dispositif de caractérisation d'une source optique, ledit dispositif comprenant : 15 - au moins un capteur, agencé pour recevoir au moins une partie d'un premier rayonnement généré par une source selon un état de la source ou différents états explorés de la source, chaque capteur comprenant un étalon Fabry Perot suivi d'une photodiode, chaque capteur étant agencé pour mesurer un signal pour chaque 20 état exploré de la source en fonction de l'au moins une partie du premier rayonnement reçu par ce capteur à cet état exploré de la source, - une unité de calcul, agencée et/ou programmée pour calculer, à partir des signaux mesurés par chaque capteur et pour chaque 25 état exploré de la source, une première donnée représentative de la longueur d'onde du premier rayonnement, l'unité de calcul étant agencée et/ou programmée pour réaliser une sélection, pour chaque état exploré de la source, d'une valeur sélectionnée de la première donnée parmi plusieurs valeurs possibles, ladite 30 sélection comprenant une élimination des valeurs de la première donnée ne correspondant pas à une constante modulo un intervalle spectral libre exprimé selon les unités de la première donnée. 3054034 - 6 - L'unité de calcul peut être en outre agencée et/ou programmée pour : - commander dans quel état se trouve la source et/ou recevoir des informations sur l'état de la source ou les différents états explorés de 5 la source, et - associer ou synchroniser chaque état exploré de la source avec le signal mesuré par chaque capteur pour cet état exploré de la source. Chaque capteur ne comprend de préférence pas de barrette de photodiodes. 10 L'étalon Fabry Perot de chaque capteur ne comprend de préférence aucune partie mobile. La première donnée illustrant la longueur d'onde du premier rayonnement peut être : - la longueur d'onde du premier rayonnement, ou 15 - la fréquence temporelle du premier rayonnement, ou - la période temporelle du premier rayonnement, ou - le nombre d'onde du premier rayonnement. L'unité de calcul peut être agencée et/ou programmée pour réaliser la sélection en éliminant en outre les valeurs de la première donnée situées en 20 dehors d'un intervalle de zone spectrale du premier rayonnement. Dans ce cas : - chaque étalon Fabry Perot peut avoir une courbe de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, ayant de préférence une période 25 supérieure ou égale à l'intervalle de zone spectrale du premier rayonnement. - l'au moins un capteur peut comprendre plusieurs capteurs comprenant chacun un étalon Fabry Perot suivi d'une photodiode, chaque étalon Fabry Perot ayant une courbe de transmission de 30 l'intensité du premier rayonnement en fonction de la longueur d'onde du premier rayonnement, l'assemblage des courbes de transmission des différents étalons Fabry Perot formant une courbe globale de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier 3054034 - 7 - rayonnement, ayant de préférence une période supérieure ou égale à l'intervalle de zone spectrale du premier rayonnement. L'au moins un capteur comprend de préférence plusieurs capteurs, les différents étalons Fabry Perot des différents capteurs comprenant de 5 préférence aux moins un couple de deux étalons Fabry Perot dont : - un premier étalon Fabry Perot ayant une première courbe de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, et - un deuxième étalon Fabry Perot ayant une deuxième courbe de 10 transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, de sorte que : - une même longueur d'onde du premier rayonnement correspond à la fois à un maximum de pente de la première courbe de 15 transmission et à une pente de la deuxième courbe de transmission inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente, et/ou - une même longueur d'onde du premier rayonnement correspond à la fois à un maximum de pente de la deuxième courbe de 20 transmission et à une pente de la première courbe de transmission inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente. L'unité de calcul peut être agencée et/ou programmée pour, avant 25 d'effectuer la sélection éliminant les valeurs de la première donnée ne correspondant pas à la constante modulo l'intervalle spectral libre de la cavité fixe, déterminer ou avoir une connaissance de la constante. L'unité de calcul peut être, pour déterminer la constante, agencée et/ou programmée pour : 30 - calculer de manière préliminaire la première donnée pour chaque état exploré de la source, et - à partir du calcul préliminaire de la première donnée pour l'ensemble de différents états explorés de la source, déterminer la constante. 3054034 - 8 - L'unité de calcul peut être agencée et/ou programmée pour réaliser la sélection comprenant, après l'élimination des valeurs de la première donnée ne correspondant pas à une constante modulo l'intervalle spectral libre, une 5 sélection finale de la valeur sélectionnée, de préférence par une méthode des moindres carrés. L'unité de calcul peut être agencée et/ou programmée pour calculer, pour chaque état exploré de la source, une deuxième donnée représentative de la longueur d'onde du deuxième rayonnement en fonction du calcul de la première donnée représentative de la longueur d'onde du 10 premier rayonnement. Suivant encore un autre aspect de l'invention, il est proposé un système, comprenant : 15 - une source optique comprenant une cavité fixe ayant un intervalle spectral libre, ladite source étant agencée pour générer un premier rayonnement selon un état de la source ou selon différents états explorés de la source, et - un dispositif de caractérisation de cette source selon l'invention. 20 La source peut être un oscillateur paramétrique optique ou un laser. La source peut être une source au moins doublement résonnante, comprenant deux cavités partageant un même cristal de génération de rayonnement, dont: - une cavité mobile générant le premier rayonnement, et 25 - la cavité fixe générant un deuxième rayonnement. Description des figures et modes de réalisation D'autres avantages et particularités de l'invention apparaîtront à la lecture de la description détaillée de mises en oeuvre et de modes de 30 réalisation nullement limitatifs, et des dessins annexés suivants : - la figure 1 est une vue schématique d'un mode de réalisation préféré de système selon l'invention 100 comprenant un mode de réalisation préféré de dispositif selon l'invention 10, 3054034 - 9 - - la figure 2 illustre les courbes de transmission des trois Fabry-Perots 5 du dispositif 10 dans la zone spectrale 15 d'émission de la source 1 du dispositif 10 ; les bruits de mesure sont schématisés par des traits verticaux sur ces courbes, 5 - la figure 3 illustre une séquence de sauts de modes émis par la source 1 du système 100 lors du déplacement d'un miroir 7a de la source 1, - la figure 4 illustre, pendant cette séquence, les valeurs du nombre d'onde de la figure 3 modulo l'intervalle spectral libre (ISL) de la cavité fixe 2 de la source 1, 10 - la figure 5 illustre différentes valeurs d'un nombre d'onde - calculées selon l'invention avant la prise en compte du modulo de l'ISL (ronds) et - des modes réellement émis (points) - la figure 6 illustre, dans le cas des valeurs de la figure 5, la 15 différence entre les mesures calculées et les nombres d'onde effectivement émis, - la figure 7 illustre la différence entre les mesures calculées ou simulées après restriction des nombres d'onde possiblement émis et les nombres d'onde effectivement émis : 20 - avant le calcul du modulo (ronds) et - après le calcul du modulo (étoiles), - la figure 8 illustre l'erreur moyenne observée sur mille séquences de sauts de modes avec un bruit de mesure aléatoire en utilisant la connaissance de l'ISL (courbe 82) et sans utiliser la 25 connaissance de l'ISL (courbe 81), et - la figure 9 illustre l'écart-type à 3x6 de l'erreur observée sur mille séquences de sauts de modes avec un bruit de mesure aléatoire en utilisant la connaissance de l'ISL (courbe 92) et sans utiliser la connaissance de l'ISL (courbe 91). 30 Ces modes de réalisation étant nullement limitatifs, on pourra notamment considérer des variantes de l'invention ne comprenant qu'une sélection de caractéristiques décrites ou illustrées par la suite isolées des autres caractéristiques décrites ou illustrées (même si cette sélection est 3054034 -10- isolée au sein d'une phrase comprenant ces autres caractéristiques), si cette sélection de caractéristiques est suffisante pour conférer un avantage technique ou pour différencier l'invention par rapport à l'état de la technique antérieure. Cette sélection comprend au moins une caractéristique de 5 préférence fonctionnelle sans détails structurels, et/ou avec seulement une partie des détails structurels si cette partie uniquement est suffisante pour conférer un avantage technique ou à différencier l'invention par rapport à l'état de la technique antérieure. 10 On va tout d'abord décrire, en référence aux figures 1 et 2, un système 100 selon l'invention comprenant un dispositif 10 selon l'invention. Le système 100 comprend une source optique 1 comprenant une cavité fixe 2 délimitée par deux miroirs 2a, 2b et ayant un intervalle spectral libre, ladite source étant agencée pour générer, selon un état de la 15 source ou selon différents états explorés de la source, un premier rayonnement 3 qui dépend de l'état de la source lors de la génération de ce rayonnement 3. L'intervalle spectral libre correspond à la distance (en longueur d'onde, fréquence, période, nombre d'onde, etc...) entre deux pics de 20 résonnance successifs de la cavité 2. La source 1 est typiquement un oscillateur paramétrique optique (OPO) ou un laser. La source 1 est une source au moins doublement résonnante, comprenant deux cavités partageant un même cristal 9 de génération de 25 rayonnement, dont: - une cavité mobile 7, délimitée par deux miroirs 7a, 7b et générant le premier rayonnement 3, et - la cavité fixe 2 générant un deuxième rayonnement 8. La source 1 doublement résonante comprend le cristal non linéaire 9 30 et trois miroirs 7a et 2a et 2b/7b dont deux (2a et 2b/7b) déposés directement sur le cristal 9 et un miroir mobile 7a. Un de ces trois miroirs 2b/7b est commun aux deux cavités 2 et 7. Dans la présente description, différents « états » de la source 1 se distinguent par une différence d'un paramètre physique de la source 1 qui 3054034 -11- influe sur la longueur d'onde du premier rayonnement 3 émis par la source 1. Dans la présente description de la figure 1, nous allons considérer le cas non limitatif où les différents états de la source 1 se différencient par 5 différentes positions du miroir 7a. On explore plusieurs états de la source 1 de manière à faire varier la longueur d'onde du premier rayonnement 3 (même si différents états de la source 1 peuvent correspondre à la même longueur d'onde du premier rayonnement 3) mais en conservant l'intervalle spectral libre de la cavité 2. 10 On va par la suite considérer le cas particulier mais non limitatif d'une source 1 qui est un oscillateur paramétrique optique doublement résonnant. La source 1 est par exemple une source de référence X-FLR8 fabriquée par Blue Industry and Science. Le système 100 comprend en outre le dispositif 10 selon l'invention 15 de caractérisation de cette source 1. Le dispositif 10 comprend au moins un capteur 4 (de préférence au moins deux ou idéalement au moins trois capteurs 4), chaque capteur 4 étant agencé pour recevoir au moins une partie du premier rayonnement 3 généré par la source 1 selon un état de la source ou différents états 20 explorés de la source 1, chaque capteur 4 comprenant un étalon Fabry Perot 5 suivi d'une photodiode 6. Chaque capteur 4 est agencé pour mesurer un signal pour chaque état exploré de la source 1 en fonction de l'au moins une partie du premier rayonnement 3 reçue par ce capteur 4 à cet état exploré de la source. 25 Le dispositif 10 comprend au moins un miroir 71, 72 ; 73 semi réfléchissant et/ou dichroïque agencé pour répartir le rayonnement 3 vers les différents capteurs 4. L'intensité lumineuse de chacune des parties du rayonnement 3 arrivant juste en amont de chacun des capteurs 4, c'est à dire : 30 - entre le miroir (ou dichroïque) 71 et le Fabry Perot 51 pour le premier capteur, - entre le miroir (ou dichroïque) 72 et le Fabry Perot 52 pour le deuxième capteur, 3054034 - 12 - - entre le miroir (ou dichroïque) 73 et le Fabry Perot 53 pour le troisième capteur, est connue soit par calibration initiale (typiquement lors de la fabrication du dispositif 10 ou en usine) ou en utilisant : 5 - un capteur optionnel de référence 66 (sans Fabry Perot) agencé pour mesurer une intensité d'une partie du rayonnement 3 sans Fabry Perot, et - une connaissance des coefficients de réflexion et de transmission de chacun des miroirs 71, 72, et 73. Chaque photodiode 6, 66 est par exemple une photodiode de 10 référence G6854-01 fabriquée par Hamamatsu. Chaque étalon 5 est par exemple un étalon de la gamme OP-7423 fabriqué par Light machinery. Chaque capteur 4 ne comprend pas de barrette de photodiodes, mais comprend une unique photodiode 6. 15 L'étalon Fabry Perot 5 de chaque capteur 4 ne comprend aucune partie mobile. Le dispositif 10 comprend une unité de calcul 11. L'unité 11 ne comprend que des moyens techniques, de préférence des moyens électroniques et/ou informatiques et/ou logiciels. 20 L'unité de calcul 11 est agencée et/ou programmée pour : - commander dans quel état se trouve la source 1 et/ou recevoir des informations sur l'état de la source 1 ou les différents états explorés de la source 1, et - associer ou synchroniser chaque état exploré de la source 1 avec le 25 signal mesuré par chaque capteur 4 pour cet état exploré de la source 1. Ainsi le signal de chaque photodiode 6 donne une information sur la transmission du rayonnement 3 à travers chaque Fabry Perot. L'unité de calcul 11 est agencée et/ou programmée pour calculer, à partir des signaux mesurés par chaque capteur 4 et pour chaque état 30 exploré de la source 1, une première donnée (par exemple ks) représentative de la longueur d'onde du premier rayonnement 3 pour chaque état exploré de la source 1, l'unité de calcul 11 étant agencée et/ou programmée pour réaliser une sélection, pour chaque état exploré de la 3054034 -13- source 1 d'une valeur sélectionnée de la première donnée (par exemple ks) parmi plusieurs valeurs possibles. La première donnée illustrant la longueur d'onde du premier rayonnement est : 5 - la longueur d'onde du premier rayonnement, ou - la fréquence temporelle du premier rayonnement, ou - la période temporelle du premier rayonnement, ou - comme décrit dans cet exemple illustratif mais non limitatif, le nombre d'onde ks du premier rayonnement 3. 10 Plus exactement, l'unité de calcul 11 est agencée et/ou programmée: - pour connaitre une constante «b» ; L'unité de calcul 11 est, dans le cas d'une génération du premier rayonnement 3 par la source 1 selon différents états explorés de la source 1, agencée et/ou programmée pour déterminer la constante «b», et pour cela est 15 agencée et/ou programmée pour : o calculer de manière préliminaire la première donnée ks pour chaque état exploré de la source 1, en éliminant des valeurs de la première donnée ks situées en dehors d'un intervalle 15 de zone spectrale du premier rayonnement 3 20 et en réalisant une sélection d'une valeur de la première donnée ks pour chaque état exploré de la source 1, de préférence par une méthode des moindres carrés o à partir du calcul préliminaire de la première donnée ks pour l'ensemble de différents états explorés de la source 1, 25 déterminer la constante «b». - pour calculer, à partir des signaux mesurés par chaque capteur 4 et pour chaque état exploré de la source 1, la première donnée ks représentative de la longueur d'onde du premier rayonnement 3 pour chaque état exploré de la source 1, en effectuant une 30 sélection d'une valeur sélectionnée de la première donnée ks parmi plusieurs valeurs possibles, l'unité de calcul 11 étant agencée et/ou programmée : o pour réaliser une élimination (qui fait partie de la sélection) des valeurs de la première donnée ks situées en 3054034 - 14 - dehors de l'intervalle 15 de zone spectrale du premier rayonnement 3 o pour réaliser une élimination (qui fait partie de la sélection) des valeurs de la première donnée ks ne 5 correspondant pas à la constante «b» modulo l'intervalle spectral libre de la cavité fixe 2 exprimé selon les unités de la première donnée k' o puis pour réaliser, après les deux précédentes éliminations, une sélection finale (qui fait partie de la sélection) de la 10 valeur sélectionnée, de préférence par une méthode des moindres carrés. L'unité de calcul 11 est agencée et/ou programmée pour calculer, pour chaque état exploré de la source 1, une deuxième donnée kc représentative de la longueur d'onde du deuxième rayonnement 8 en 15 fonction du calcul de la première donnée ks représentative de la longueur d'onde du premier rayonnement 3. La deuxième donnée kc illustrant la longueur d'onde du deuxième rayonnement 8 est : - la longueur d'onde du deuxième rayonnement, ou 20 - la fréquence du deuxième rayonnement, ou - comme décrit dans cet exemple illustratif mais non limitatif, le nombre d'onde kc du deuxième rayonnement 8. La manière dont est agencée et/ou programmée l'unité de calcul 11 sera décrite ci-dessous plus en détails en référence à un mode de 25 réalisation de procédé selon l'invention mis en oeuvre par le dispositif 10. En référence à la figure 2, chaque étalon Fabry Perot 5 a une courbe de transmission 12, 13, 14 de l'intensité du premier rayonnement 3, en fonction de la longueur d'onde du premier rayonnement 3, ayant une période supérieure ou égale à l'intervalle 15 de zone spectrale du premier 30 rayonnement 3. Chaque étalon Fabry Perot 5 a une courbe de transmission 12, 13, 14 de l'intensité du premier rayonnement 3 en fonction de la longueur d'onde du premier rayonnement 3, l'assemblage des courbes de transmission des différents étalons Fabry Perot 5 formant une courbe globale 16 de 3054034 -15- transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement 3, ayant une période supérieure ou égale à l'intervalle 15 de zone spectrale du premier rayonnement. L'au moins un capteur 4 comprend plusieurs capteurs 4, les différents 5 étalons Fabry Perot 5 des différents capteurs 4 comprenant aux moins un couple de deux étalons Fabry Perot 5 (de préférence uniquement des couples d'étalons Fabry Perot), chaque couple comprenant : - un premier étalon Fabry Perot ayant une première courbe de transmission 13 de l'intensité du premier rayonnement 3, en 10 fonction de la longueur d'onde du premier rayonnement 3, et - un deuxième étalon Fabry Perot ayant une deuxième courbe de transmission 12 de l'intensité du premier rayonnement 3, en fonction de la longueur d'onde du premier rayonnement 3, de sorte que : 15 - une même longueur d'onde du premier rayonnement 3 correspond à la fois à un maximum de pente (c'est-à-dire de la dérivée première) 130 de la première courbe de transmission 13 et à une pente 120 de la deuxième courbe de transmission 12 inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente 130, 20 et/ou - une même longueur d'onde du premier rayonnement 3 correspond à la fois à un maximum de pente 121 de la deuxième courbe de transmission 12 et à une pente 131 de la première courbe 13 de transmission inférieure ou égale, en valeur absolue, à 10% de ce 25 maximum de pente 121. On remarque que cela est valable pour tous les couples 12, 13 ou 13, 14 ou 14, 12 des courbes de transmission de tous les étalons Fabry Perot 5 du dispositif 10. On remarque en outre que les maximums de pente 121, 130 de 30 même signe (soit pente croissante soit pente décroissante) des différentes courbes de transmission 12, 13, 14 sont régulièrement espacés par un même intervalle de séparation 17 à l'intérieur de l'intervalle 15. Ainsi, dans ce mode de réalisation, on utilise des étalons Fabry-Perot 5 (sans pièces mécaniques mobiles) sans imager la figure d'interférences 3054034 - 16 - sur un détecteur barrette. On ne fait qu'un seul point de mesure par étage de Fabry-Perot 5. Pour améliorer la précision on réalise un accord en longueur d'onde dans toute la bande de gain d'émission de la source 1 en ne modifiant 5 qu'une seule cavité 7 de la source cohérente 1. Ainsi, on sait que toutes les longueurs d'onde sont identiques modulo l'ISL de la cavité 2 de la source 1. Cela donne un supplément d'information qui permet de remonter à chacune des longueurs d'onde émises des faisceaux 3, 8. Ce mode de réalisation a les avantages suivants : 10 - très compact : la solution peut faire moins de 10 cm de long. - très rapide : en utilisant une électronique adaptée 11 on peut faire une mesure à chaque impulsion laser (> 100 kHz). L'invention peut fonctionner aussi avec une source 1 continue avec une mesure au maximum de l'électronique 11 de mesure et de 15 traitement des photodiodes 6. L'invention n'est pas limitée par le temps de réponse d'un piézoélectrique qui déplacerait une partie d'un Fabry Perot accordable. - peu coûteux : il ne comprend pas de composants chers : les étalons Fabry Perot non accordables 5, les photodiodes 6 sans 20 barrette sont très économiques par rapport à un système à barrette On va maintenant décrire, en référence aux figures 1 à 9, un mode de réalisation de procédé selon l'invention mis en oeuvre par le dispositif 10 25 et par le système 100. Ce procédé comprend les étapes 1) à 8) décrites ci-après: 1) une génération du premier rayonnement 3 (et du deuxième rayonnement 8) par la source 1 selon un état de la source ou différents états explorés de la source 1. La source cohérente 1 est un OPO 30 doublement résonnant. Un faisceau laser de pompe de nombre d'onde kp traverse le cristal 9 non linéaire à polarisation périodiquement inversée. Le cristal 9 rayonne deux nouveaux faisceaux : le premier rayonnement 3 (aussi appelé signal) de nombre d'onde ks et le deuxième rayonnement 8 3054034 - 17 - (aussi appelé complémentaire) de nombre d'onde kc. Les nombres d'onde sont reliés par la loi de conservation de l'énergie qui implique que : kp = ks + kc Comme schématisé à la figure 1, le rayonnement 3 résonne dans la 5 cavité mobile 7 formée par les miroirs 7a et 7b. Le rayonnement 8 résonne dans la cavité fixe 2 entre les miroirs 2a et 2b. La cavité mobile 7 a un intervalle spectral libre (ISL) Aks d'environ 0,26 cm-1 (qui dépend de la position du miroir 7a); la cavité fixe 2 a un intervalle spectral libre (ISL) Akc de 0,25 cm-1. 10 Le miroir 7a est mobile, par exemple monté sur un élément piézoélectrique. L'acceptance spectrale (dite bande de gain paramétrique) typique de la source 1 est de 5 cm-1. C'est-à-dire que le cristal non linaire 9 ne peut rayonner que sur une zone du spectre de 5 cm-1 si les paramètres suivants 15 sont fixés : nombre d'onde du laser de pompe, température du cristal 9, pas du réseau d'inversion de polarisation , etc... Le nombre de modes possiblement émis pour une position fixe de la cavité complémentaire est donc de 5/0,25 = 20 modes. Nous supposerons les nombres d'onde suivants pour chacun des 20 faisceaux : - nombre d'onde du rayonnement pompe : 9 394,5 cm-1 ; - nombre d'onde « signal » ou du rayonnement 3 : de 6 662,5 cm-1 à 6 667,5 cm-1 (cet intervalle étant nommé, dans la présente description, l'intervalle 15 de zone spectrale du premier rayonnement); 25 - nombre d'onde « complémentaire » ou du rayonnement 8: de 2 727 cm-1 à 2 732 cm-1. 2) une commande, par l'unité 11, de l'état dans lequel se trouve la source 1 et/ou une réception, par l'unité 11, des informations sur l'état de 30 la source 1 ou sur les différents états explorés de la source 1 ; comme expliqué précédemment, on explore différents états de la source 1 en changeant un paramètre physique de la source 1 qui influe sur la longueur d'onde du premier rayonnement 3 émis par la source, par exemple en changeant la position du miroir 7a. On explore plusieurs états de la source 1 3054034 -18- de manière à faire varier la longueur d'onde du premier rayonnement 3 (même si différents états de la source 1 peuvent correspondre à la même longueur d'onde du premier rayonnement 3) mais en conservant l'intervalle spectral libre de la cavité fixe 2. 5 Chacun des faisceaux 3, 8 a un peigne de mode possible dans sa cavité respectivement 7, 2. L'oscillation a lieu si, au sein de la bande de gain paramétrique, il se trouve un recouvrement suffisant entre deux modes signal et complémentaire, c'est-à-dire si la condition de conservation de l'énergie est respectée. 10 En appliquant une tension sur l'élément piézoélectrique du miroir 7a, la cavité mobile 7 va légèrement bouger de quelques centaines de nanomètres. Le peigne de mode associé va alors se décaler et les coïncidences vont aussi se déplacer au gré des nouveaux recouvrements entre les modes signal du faisceau 1 et complémentaire du faisceau 8. 15 La figure 3 présente un exemple d'une série de sauts de modes ks effectués par la source 1 lors du déplacement du miroir 7a et la figure 4 présente ces mêmes nombres d'onde ks modulo l'ISL de la cavité complémentaire 2. Seul le miroir 7a a bougé. La cavité complémentaire 2 est donc restée 20 fixe. Le peigne des nombres d'onde possiblement émis par la cavité 2 est donc toujours le même. Tous les modes du faisceau 8 sont donc par définition égaux modulo l'ISL de la cavité 2 (nous négligeons ici la dispersion de l'ISL). Par conservation de l'énergie, ceux du faisceau 3 le sont aussi. C'est pourquoi, tous les nombres d'onde du faisceau 8 sont 25 reliés par l'équation 1 et par conservation de l'énergie, ceux du faisceau 3 sont reliés par l'équation 2 : kc = a [ISLc] équation 1 ks = b [ISLc] équation 2 où chacune parmi « a » et « b » est une constante, plus exactement un 30 nombre réel supérieur ou égal à zéro, la constante « b » étant celle déjà mentionnée précédemment dans la présente description, [ISLc] signifiant « modulo l'ISL de la cavité fixe ISLc», ISLc étant l'ISL de la cavité fixe exprimée avec les unités d'un nombre d'onde (inverse d'une longueur). De manière générale, on notera que ce principe de l'invention est valable et 3054034 - 19 - généralisable pour toute source 1 dont un milieu à gain (cristal laser, cristal non linéaire...) est contenu dans une cavité résonnante. Ces deux équations pourraient naturellement être exprimées en fréquence temporelle ou période temporelle ou longueur d'onde. Par exemple, dans le cas de 5 l'équation 2 : f, = b' [ISLc'] équation 2' t, = b" [ISLc"] équation 2" À, = b"' [ISLcm] équation 2"' 10 où : - chacune parmi « b' » « b" » et « b"' » est une constante, plus exactement un nombre réel supérieur ou égal à zéro - f, est la fréquence temporelle du rayonnement 3 - [ISLc'] signifiant « modulo l'ISL de la cavité fixe ISLc'», ISLc' étant 15 l'ISL de la cavité fixe exprimée avec les unités d'une fréquence temporelle (inverse d'un temps) - t, est la période temporelle du rayonnement 3 - [ISLc"] signifiant « modulo l'ISL de la cavité fixe ISLc"», ISLc" étant l'ISL de la cavité fixe exprimée avec les unités d'une période temporelle (un 20 temps) - À, est la longueur d'onde du rayonnement 3 - [ISLcm] signifiant « modulo l'ISL de la cavité fixe ISLc"'», ISLc"' étant l'ISL de la cavité fixe exprimée avec les unités d'une longueur d'onde (une longueur) 25 3) une réception d'au moins une partie du premier rayonnement 3 par l'au moins un capteur 4, et même par les plusieurs capteurs 4 dans le présent mode de réalisation. 30 4) une mesure d'un signal par chaque capteur 4 et pour chaque état exploré de la source 1, en fonction de l'au moins une partie du premier rayonnement 3 reçu par ce capteur 4 à cet état exploré de la source 1. Les Fabry-Perot 51, 52, 53 sont dimensionnés de sorte que leur intervalle spectral libre (ISL) soit respectivement de 961,29 pm, 961,54 pm 3054034 -20- et 961,79 dam. Leur transmission en fonction du nombre d'onde dans la bande de gain paramétrique est décrite par la figure 2. Sur la figure 2, a aussi été simulée l'intensité lue par chaque photodiode 6 pour chacun des modes de la figure 3 à laquelle a été ajouté un bruit 5 blanc de 3 % de l'intensité maximum à 3x6, où 6 représente l'écart type du bruit. 5) une association ou synchronisation, par l'unité 11, de chaque état exploré de la source 1 avec le signal mesuré par chaque capteur 4 pour cet 10 état exploré de la source 1, 6) une connaissance a priori de la constante «b» ; dans le cas d'une génération du premier rayonnement 3 par la source 1 selon différents états explorés de la source 1, la connaissance de la constante «b» peut provenir 15 d'une détermination de cette constante «b», la détermination de cette constante «b» comprenant: a. un calcul préliminaire de la première donnée ks représentative de la longueur d'onde du premier rayonnement 3 pour chaque état exploré de la source 1, le calcul comprenant, pour chaque 20 état exploré de la source 1, une sélection d'une valeur de la première ks donnée parmi plusieurs valeurs possibles, ladite sélection comprenant : i. une élimination des valeurs de la première donnée ks situées en dehors de l'intervalle 15 de zone spectrale du 25 premier rayonnement 3. L'objectif est donc désormais de retrouver à partir des informations des photodiodes 6 les nombres d'onde ks du rayonnement 3 émis par la source 1 dans sa zone d'émission possible : de 6 662,5 cm-1 à 6 667,5 cm-1 soit 5 cm-1. Pour atteindre une 30 résolution de 2.10-3 cm-1, on échantillonne la zone d'émission 15 en 2500 points espacés de 2.10-3 cm-1. À chacun de ces points correspond un triplet de transmissions des trois photodiodes 6. 3054034 -21- ii. une sélection de la valeur sélectionnée, de préférence par une méthode des moindres carrés. Par la méthode des moindres carrés, on compare chacun des triplets mesurés aux 2500 triplets théoriques selon l'équation 5 3 mi = V[(Tithi - Tm1)2 + (rth2- Tm2)2 + (rth3- Tm3)2] équation 3 où Tith, correspond à la transmission théorique du i-ème Fabry-Perot (sur les 3 possibles) du j-ème triplet (sur les 2500 possibles), et Tm, correspond à la transmission mesurée du i-ème Fabry-Perot. Pour chacun des triplets mesurés on retient la valeur de nombre d'onde ks du faisceau 3 10 parmi les 2500 possibilités, qui minimise la valeur de « rn, ». Ces nombres d'onde calculés ks sont présentés par à la figure 5. Malheureusement, le bruit de mesure perturbe le calcul et dégrade le résultat. La figure 6 présente la différence point à point entre le calcul et les nombres d'onde effectivement émis. 15 b. à partir du calcul préliminaire de la première donnée ks pour l'ensemble de différents états explorés de la source 1, une détermination de la constante «b». En effet, à partir de ces mesures bruitées de la figure 5, nous pouvons tout de même calculer précisément le « b » présenté à l'équation 2. En effet, 20 le bruit de mesure est alors moyenné sur l'ensemble des 250 acquisitions. S'il s'agit d'un bruit blanc gaussien, le bruit est alors diminué d'un facteur supérieur à 15 (-V250). Afin de limiter l'effet d'un nombre aberrant nous utilisons la médiane (med) plutôt que la moyenne. Nous avons ainsi : b = med(ksmesurés [ISLc]) équation 4 25 Dans notre exemple, b = 0,1404 cm-1. Comme nous l'avons vu plus haut, nous savons maintenant que la source 1 ne peut émettre qu'à des nombres d'onde qui respectent l'équation 2 avec la valeur « b » que nous venons de calculer (équation 4). Une autre alternative à la détermination de b serait de l'enregistrer (par 30 exemple en usine) dans le dispositif 10, mais cela ne permet pas de prendre en compte les variations de b dues par exemple aux variations de température. 3054034 -22- 7) à partir des signaux mesurés et pour chaque état exploré de la source 1, un calcul de la première donnée ks représentative de la longueur d'onde du premier rayonnement 3 pour chaque état exploré de la source 1, 5 le calcul comprenant, pour chaque état exploré de la source 1, une sélection d'une valeur sélectionnée de la première donnée ks parmi plusieurs valeurs possibles, ladite sélection comprenant : a. une élimination des valeurs de la première donnée ks situées en dehors de l'intervalle 15 de zone spectrale du premier 10 rayonnement 3, comme à l'étape 6.a.i précédente b. une élimination des valeurs de la première donnée ks ne correspondant pas à une constante «b» modulo l'intervalle spectral libre de la cavité fixe 2 exprimé selon les unités de la première donnée 15 c. puis, après les deux précédentes éliminations, une sélection finale de la valeur sélectionnée, de préférence par une méthode des moindres carrés. On réutilise la méthode des moindres carrés en limitant les solutions possibles aux seuls nombres d'onde de la bande de gain paramétrique (zone 20 spectrale 15) respectant l'équation 2. Ainsi contrairement à l'étape 6), nous n'avons plus le choix entre 2500 triplets possibles mais entre 19 modes possibles respectant la condition de modulo. Si b = 0, nous avons alors 20 modes possibles. On refait le calcul des « m, » de l'équation 3 pour les 25 seuls 19 triplets possibles. La valeur de nombre d'onde retenue pour une mesure donnée est celle qui minimise « m, » parmi les 19 valeurs calculées. On obtient alors la figure 7 qui présente, en plus des résultats 60 (ronds éparpillés) présentés à la figure 6, les résultats 70 issus de ces nouveaux moindres 30 carrés (étoiles quasiment toutes alignées). La figure 7 est éloquente, le bruit est quasiment nul et seul un point de mesure est faux à un ISL près (mesure n° 196). Un filtre éliminant les points de mesures à proximité d'un saut de modes peut l'éliminer. Un tel filtre est connu pour la spectroscopie ; en effet, les points proches des sauts 3054034 -23- de modes reflètent une source qui n'est pas parfaitement pure spectralement. Ils sont donc filtrés pour garantir la qualité de la mesure indépendamment du mesureur de nombre d'onde. Dans l'exemple décrit, on obtient alors les résultats présentés au 5 tableau ci-dessous : Sans modulo (points 60) Avec modulo (points 70) Erreur moyenne (x10-3 cm-1) 22,9 0,04 Écart-type (3x6, x10-3 cm-1) 19,2 0,0 Pourcentage de mesures 70,4 69,2 restantes après filtrage (%) Les résultats sont nettement meilleurs en utilisant la connaissance de l'ISL de la cavité fixe 2. Nous pouvons aussi voir que le filtrage ne supprime 10 pas beaucoup plus de points de mesures lors qu'il est appliqué au premier calcul (sans utilisation du modulo) plutôt qu'au second. Dans les deux cas il reste environ 70 % des acquisitions. La figure 8 et la figure 9 présentent les résultats suivants: - erreur moyenne au premier calcul (référence 81, sans utilisation du 15 modulo) et au second calcul (référence 82, avec utilisation du modulo de l'équation 2), et - écart-type au premier calcul (référence 91, sans utilisation du modulo) et au second calcul (référence 92, avec utilisation du modulo de l'équation 2), 20 pour mille simulations de séquence de sauts de modes, le bruit étant aléatoire d'une séquence à l'autre. 8) Dans ce mode de réalisation de procédé selon l'invention, la cavité mobile 7 génère le premier rayonnement 3 et la cavité fixe 2 génère le 25 deuxième rayonnement 8. Ce procédé comprend en outre un calcul, pour chaque état exploré de la source, de la deuxième donnée kc représentative de la longueur d'onde du deuxième rayonnement 8 en fonction du calcul de la première donnée ks représentative de la longueur d'onde du premier rayonnement 3. Connaissant : 30 - la valeur de kp de la pompe, connue et maîtrisée, et qui a priori ne change pas lors des changements d'états de la source 1, et 3054034 -24- - les différentes valeurs de ks calculées pour les différents états explorés de la source 1, alors on calcul les différentes valeurs de kc calculées pour les différents états explorés de la source 1, avec l'équation : 5 kc = kp - ks Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. 10 Ainsi, dans des variantes combinables entre elles des modes de réalisations précédemment décrits : - le procédé selon l'invention ne calcule pas nécessairement la constante « b », qui peut être déjà connue (par exemple par une calibration en usine ou autre). Dans ce cas, le procédé ou dispositif selon l'invention peut être 15 mis en oeuvre en explorant qu'un seul état de la source 1, c'est-à-dire avec un seul point 70 de la figure 7 ; - la source 1 peut ne comprendre qu'une seule cavité fixe 2, sans cavité mobile 7. Dans ce cas, la cavité 2 peut comprendre le cristal 9 entre les miroirs 2a et 2b fixes. L'intervalle spectral libre ISLc de la cavité 2 est 20 toujours fixe ; on change l'état de la source 1 par exemple en changeant la température du cristal 9, et le dispositif et le procédé selon l'invention sont mis en oeuvre cette fois ci sur la base de l'équation suivante : kc = b [ISLc], la constante b étant connue ou déterminée sur le même principe que précédemment décrit, et la première donnée étant kc; 25 - le dispositif et le procédé selon l'invention peut être basé sur des mesures, par les capteurs 4, sur le signal de la cavité 2 plutôt que sur le signal de la cavité 7, dans un cas avec un miroir 2a ou 2b mobile et de préférence le miroir 7a fixe, par exemple en utilisant un gaz à la place d'un cristal 9 ; 30 - la première donnée n'est pas nécessairement le nombre d'onde ; on peut aussi mettre en oeuvre le dispositif et le procédé selon l'invention avec comme première donnée une fréquence temporelle, une période temporelle, ou une longueur d'onde du rayonnement reçu par l'au moins un capteur 4. 3054034 -25- Bien entendu, les différentes caractéristiques, formes, variantes et modes de réalisation de l'invention peuvent être associées les unes avec les autres selon diverses combinaisons dans la mesure où elles ne sont pas incompatibles ou exclusives les unes des autres. En particulier toutes les 5 variantes et modes de réalisation décrits précédemment sont combinables entre eux. 10 REVENDICATIONS1. Procédé de caractérisation d'une source optique (1), ladite source optique comprenant une cavité fixe (2) ayant un intervalle spectral libre, ledit procédé comprenant : une génération d'un premier rayonnement (3) par la source selon un état de la source ou différents états explorés de la source, une réception d'au moins une partie de ce premier rayonnement par au moins un capteur (4), chaque capteur comprenant un étalon Fabry Perot (5) suivi d'une photodiode (6), une mesure d'un signal par chaque capteur et pour chaque état exploré de la source, en fonction de l'au moins une partie du premier rayonnement reçu par ce capteur à cet état exploré de la source, à partir des signaux mesurés et pour chaque état exploré de la source, un calcul d'une première donnée représentative de la longueur d'onde du premier rayonnement, le calcul comprenant, pour chaque état exploré de la source, une sélection d'une valeur sélectionnée de la première donnée parmi plusieurs valeurs possibles, ladite sélection comprenant une élimination des valeurs de la première donnée ne correspondant pas à une constante (b) modulo l'intervalle spectral libre de la cavité fixe exprimé selon les unités de la première donnée. 2. Procédé selon la revendication 1, caractérisé en ce que chaque capteur est du type sans barrette de photodiodes. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que toutes les parties de l'étalon Fabry Perot de chaque capteur sont immobiles. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la première donnée illustrant la longueur d'onde du premier rayonnement est : - la longueur d'onde du premier rayonnement, ou 3054034 - 27- - la fréquence du premier rayonnement, ou - le nombre d'onde du premier rayonnement. 5. Procédé selon l'une quelconque des revendications précédentes, 5 caractérisé en ce que ladite sélection élimine en outre les valeurs de la première donnée situées en dehors d'un intervalle (15) de zone spectrale du premier rayonnement. 6. Procédé selon la revendication 5, caractérisé en ce que chaqùe étalon 10 Fabry Perot a une courbe de transmission (12, 13, 14) de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, ayant une période supérieure ou égale à l'intervalle de zone spectrale du premier rayonnement. 15 7. Procédé selon la revendication 5 ou 6, caractérisé en ce que la réception d'au moins une partie du premier rayonnement comprend une réception par plusieurs capteurs (4) comprenant chacun un étalon Fabry Perot (5) suivi d'une photodiode (6), chaque étalon Fabry Perot ayant une courbe de transmission de l'intensité du 20 premier rayonnement en fonction de la longueur d'onde du premier rayonnement, l'assemblage des courbes de transmission des différents étalons Fabry Perot (5) formant une courbe globale (16) de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, ayant une période 25 supérieure ou égale à l'intervalle de zone spectrale du premier rayonnement. 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'au moins un capteur comprend plusieurs 30 capteurs, les différents étalons Fabry Perot des différents capteurs comprenant aux moins un couple de deux étalons Fabry Perot dont : - un premier étalon Fabry Perot ayant une première courbe de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, et 3054034 -28- - un deuxième étalon Fabry Perot ayant une deuxième courbe de transmission de l'intensité du premier rayonnement, en fonction de la longueur d'onde du premier rayonnement, de sorte que : 5 - une même longueur d'onde du premier rayonnement (3) correspond à la fois à un maximum de pente (130) de la première courbe de transmission (13) et à une pente (120) de la deuxième courbe de transmission (12) inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente (130), et/ou 10 - une même longueur d'onde du premier rayonnement (3) correspond à la fois à un maximum de pente (121) de la deuxième courbe de transmission (12) et à une pente (131) de la première courbe de transmission (13) inférieure ou égale, en valeur absolue, à 10% de ce maximum de pente (121). 15 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend, avant la sélection éliminant les . - orrespon an pas a a cons an e (b) modulo l'intervalle spectral libre de la cavité fixe, une 20 détermination ou une connaissance de la constante (b). 10. Procédé selon la revendications 9, caractérisé en ce que la détermination de la constante (b) comprend, dans le cas d'une génération du premier rayonnement (3) par la source selon différents 25 états explorés de la source : - un calcul préliminaire de la première donnée pour chaque état exploré de la source, et - à partir du premier calcul préliminaire de la première donnée pour 30 l'ensemble de différents états explorés de la source, une détermination de la constante (b). 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la sélection d'une valeur 35 sélectionnée de la première donnée comprend, après l'élimination des 3054034 - 29- valeurs de la première donnée ne correspondant pas à une constante (b) modulo l'intervalle spectral libre de la cavité, une sélection finale de la valeur sélectionnée, de préférence par une méthode des moindres carrés. 5 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la source est un oscillateur paramétrique optique ou un laser. 10 13. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la source est une source au moins doublement résonnante, comprenant deux cavités partageant un même cristal (9) de génération de rayonnement, dont: 15 - une cavité mobile (7) générant le premier rayonnement (3), et - la cavité fixe (2) générant un deuxième rayonnement (8). 14. Procédé selon la revendication 13 comprend en outre un calcul, pour chaque état exploré de la source, 20 d'une deuxième donnée représentative de la longueur d'onde du deuxième rayonnement en fonction du calcul de la première donnée représentative de la longueur d'onde du premier rayonnement. 15. Dispositif (10) de caractérisation d'une source optique (1), ledit 25 dispositif comprenant : au moins un capteur (4), agencé pour recevoir au moins une partie d'un premier rayonnement (3) généré par une source selon un état de la source ou différents états explorés de la source, chaque capteur comprenant un étalon Fabry Perot (5) suivi d'une 30 photodiode (6), chaque capteur étant agencé pour mesurer un signal pour chaque état exploré de la source en fonction de l'au moins une partie du premier rayonnement reçu par ce capteur à cet état exploré de la source, 3054034 - 30- une unité de calcul (11), agencée et/ou programmée pour calculer, à partir des signaux mesurés par chaque capteur et pour chaque état exploré de la source, une première donnée représentative de la longueur d'onde du premier rayonnement, 5 l'unité de calcul étant agencée et/ou programmée pour réaliser une sélection, pour chaque état exploré de la source, d'une valeur sélectionnée de la première donnée parmi plusieurs valeurs possibles, ladite sélection comprenant une élimination des valeurs de la première donnée ne correspondant pas à une constante (b) 10 modulo un intervalle spectral libre exprimé selon les unités de la première donnée. 15 20 25 30