FR 3150298 A1 20241227 FR 2306538 A 20230622 Ce dispositif (22) d’estimation d’un déphasage, entre deux signaux d’entrée (s 1 , s 2 ) issus d’un même signal périodique (s 0 ) de période T traité par interférométrie, comporte un coupleur multimodal (24) recevant les deux signaux d’entrée (s 1 , s 2 ) et fournissant un système multiphasé (S), une unité (30) de traitement du système multiphasé (S) pour un calibrage et une estimation (SG) du déphasage, et des moyens (32) de stockage de paramètres de calibrage. L’unité de traitement (30) est configurée (34) pour tester si le système multiphasé (S) remplit des conditions prédéterminées requises pour réaliser une estimation des paramètres de calibrage et pour : réaliser l’estimation des paramètres de calibrage et leur enregistrement dans les moyens de stockage (32) si les conditions prédéterminées requises sont remplies ; extraire, des moyens de stockage (32), les paramètres de calibrage d’une estimation précédente si les conditions prédéterminées requises ne sont pas remplies. Figure pour l’abrégé : Fig. 1 Dispositif et procédé d’estimation d’un déphasage interférométrique, dispositif électronique de caractérisation d’un fluide La présente invention concerne un dispositif d’estimation d’un déphasage entre deux signaux d’entrée issus d’un même signal périodique traité par interférométrie, un dispositif électronique de caractérisation d’un fluide à transducteur comportant un tel dispositif d’estimation de déphasage, et un procédé d’estimation d’un déphasage entre deux signaux d’entrée issus d’un même signal périodique traité par interférométrie. Un dispositif électronique de caractérisation d’un fluide, liquide ou gazeux, par interférométrie comporte au moins un capteur conçu pour interagir avec le fluide et un transducteur traitant, par interférométrie, les signaux issus de cette interaction pour fournir une séquence d’au moins un signal de mesure représentative d’une interaction du fluide avec chaque capteur. Il est bien entendu que l’interaction avec le fluide, la fourniture de la séquence de signaux de mesure et son éventuel traitement sont des opérations qui peuvent s’effectuer en parallèle au cours du temps. En variante, le traitement peut être exécuté a posteriori . Ce type de dispositif est parfois appelé capteur multivarié (de l’anglais « multivariate sensor »), notamment lorsqu’il comporte une pluralité de capteurs. Il peut être envisagé pour des applications variées d’estimations olfactives et est généralement, dans ce cas, qualifié de « nez électronique » (en milieu gazeux) ou « langue électronique » (en milieu liquide). Il est alors utilisé pour détecter, discriminer, identifier et quantifier des composés organiques volatils dans un fluide gazeux à analyser, ou des composés présents dans un liquide. Il peut être exploité dans différents domaines industriels tels que : l’industrie du parfum, par exemple pour comparer, étudier ou concevoir des compositions olfactives pures ou mélangées, la protection de l’environnement, en particulier pour détecter les pollutions odorantes ou la surveillance de la qualité d’environnements plus ou moins confinés, la surveillance de sites industriels présentant un risque de contamination par des matériaux volatils ou solubles dans un solvant potentiellement dangereux ou odorants, la santé, par exemple pour proposer un substitut d’odorat à des personnes souffrant d’anosmie ou pour détecter des marqueurs biologiques volatils tels que les émanations d’une activité microbiologique infectieuse, l’industrie alimentaire, par exemple pour détecter des contaminations dans une chaîne de fabrication et/ou distribution de nourriture, n’importe quel autre domaine industriel dans lequel le contrôle d’un produit odorant quelconque peut s’avérer utile. L'invention s’applique plus particulièrement à un dispositif d’estimation d’un déphasage entre deux signaux d’entrée issus d’un même signal périodique de période T traité par interférométrie, comportant : un coupleur multimodal à deux entrées, pour la réception des deux signaux d’entrée, et à N sorties, N ≥ 3, pour la fourniture d’un système multiphasé de N signaux formés pour être déphasés entre eux de T/N ; et une unité de traitement du système multiphasé pour en extraire des paramètres de calibrage et une estimation du déphasage entre les deux signaux d’entrée. Un tel dispositif d’estimation de déphasage est notamment particulièrement avantageux pour un dispositif électronique de caractérisation d’un fluide par interférométrie parce qu’il permet de déterminer sans aucune ambiguïté modulo T l’évolution au cours du temps du décalage de phase, caractéristique de l’interaction du fluide avec tout capteur adapté, entre les deux signaux qui lui sont fournis. Mais il est applicable plus généralement à tout système produisant deux signaux, issus d’un même signal périodique de période T traité par interférométrie, dont on veut déterminer le décalage de phase à tout instant sans ambiguïté modulo T. En général le signal périodique est de forme sinusoïdale et de période exprimée angulairement, d’où une période T généralement égale à 2π. Il est connu d’appliquer une transformation de Clarke (ou toute autre transformation équivalente telle qu’une transformation de Concordia ou Fortescue) au système multiphasé pour en obtenir théoriquement une représentation circulaire centrée à partir de laquelle il est ensuite très aisé d’extraire une valeur du déphasage à estimer à chaque instant par simple application de la fonction Arc tangente notée Arctan dans le cas de signaux sinusoïdaux dont la périodicité est égale à 2π lorsqu’elle est exprimée angulairement. Il reste ensuite souvent nécessaire de procéder à une opération de déroulement de phase (de l’anglais « phase unwrapping ») pour traiter les discontinuités de modulo. Mais dans la réalité, le système multiphasé présente des défauts qui sont dus aux imperfections du traitement par interférométrie, du coupleur multimodal, d’un éventuel capteur en sortie du coupleur multimodal (par exemple un capteur photographique CCD lorsqu’il s’agit d’interférométrie optique) et éventuellement du générateur de signal périodique dont sont issus les deux signaux d’entrée. Ces défauts se présentent notamment sous la forme de différences entre les valeurs d’amplitudes des N signaux qui le constituent : par couplage multimodal et éventuelle capture optique, ces N signaux devraient présenter la même amplitude, mais ce n’est pas le cas. Les défauts observés se présentent en outre généralement sous la forme de décalages en amplitude des N signaux par rapport à l’axe nul autour duquel ils devraient évoluer : par génération du signal périodique initial, couplage multimodal et éventuelle capture optique, ces N signaux devraient tous évoluer autour de l’axe nul, ou au moins avec un même décalage par rapport à cet axe nul, mais ce n’est pas le cas non plus. Il en résulte que la transformation de Clarke d’un tel système multiphasé imparfait fournit sur une séquence, non pas un cercle centré, mais une ellipse décentrée. Ainsi, une application de la fonction Arctan sur le résultat de la transformation de Clarke à chaque instant produit une estimation biaisée du déphasage. Une première solution pour résoudre ces défauts du système multiphasé est apportée dans l’article de Halir et al, intitulé « Direct and sensitive phase readout for integrated waveguide sensors », publié dans IEEE Photonics Journal, volume 5, numéro 4, en août 2013. Elle consiste à estimer des paramètres de calibrage par changement de coordonnées à partir d’un recalage de l’ellipse obtenue par transformation de Clarke sur un cercle centré, sur une séquence suffisamment longue, à enregistrer ce recalage en mémoire en tant que paramètres de calibrage, puis à appliquer ce recalage par changement de coordonnées à chaque instant sur la transformation de Clarke du système multiphasé, avant de procéder au déroulement de phase. Cette approche est avantageuse parce que le calibrage ne nécessite que les données du système multiphasé lui-même, ce qui permet de le qualifier de calibrage aveugle. En revanche, l’estimation des paramètres du calibrage aveugle peut s’avérer très imprécise, notamment lorsque les mesures exploitées pour le calibrage ne permettent pas de reconstituer avec une certaine précision l’ellipse à partir de laquelle les paramètres de calibrage sont calculés. Or rien n’est prévu pour s’en assurer et le fait de prendre des mesures sur une séquence suffisamment longue ne garantit rien. Une deuxième solution pour résoudre ces défauts du système multiphasé, qui prend la forme d’un perfectionnement relativement proche de la première, est apportée par le document de brevet WO 2022/238170 A1. Elle consiste à estimer les paramètres de calibrage à partir des valeurs d’amplitudes et de décalages en amplitude du système multiphasé sur une séquence suffisamment longue, à enregistrer ces paramètres en mémoire, puis à appliquer ce calibrage par traitement à chaque instant du système multiphasé avant de lui appliquer la transformation de Clarke et l’opération de déroulement de phase. Il s’agit là encore avantageusement d’un calibrage aveugle. Mais cette deuxième solution présente le même inconvénient que la précédente, à savoir une estimation des paramètres de calibrage qui peut s’avérer très imprécise, notamment lorsque les mesures exploitées pour le calibrage ne permettent pas d’évaluer les valeurs d’amplitudes et de décalages en amplitude du système multiphasé. Or rien n’est prévu pour s’en assurer et le fait de prendre des mesures sur une séquence suffisamment longue ne garantit rien non plus dans ce cas. Elle a cependant le mérite de montrer qu’il est possible de procéder au calibrage aveugle directement à partir des valeurs d’amplitudes et de décalages en amplitude du système multiphasé. Il peut ainsi être souhaité de prévoir un dispositif d’estimation de déphasage à calibrage aveugle qui permette de s’affranchir d’au moins une partie des problèmes et contraintes précités. Il est donc proposé un dispositif d’estimation d’un déphasage entre deux signaux d’entrée issus d’un même signal périodique de période T traité par interférométrie, comportant : un coupleur multimodal à deux entrées, pour la réception des deux signaux d’entrée, et à N sorties, N ≥ 3, pour la fourniture d’un système multiphasé de N signaux formés pour être déphasés entre eux de T/N ; une unité de traitement du système multiphasé pour en extraire des paramètres de calibrage et une estimation du déphasage entre les deux signaux d’entrée ; et des moyens de stockage des paramètres de calibrage ; dans lequel l’unité de traitement est configurée pour tester si les N signaux du système multiphasé remplissent des conditions prédéterminées requises pour réaliser une estimation des paramètres de calibrage et pour : réaliser l’estimation des paramètres de calibrage et leur enregistrement dans les moyens de stockage si les conditions prédéterminées requises sont remplies ; et extraire des moyens de stockage, les paramètres de calibrage d’une estimation précédente si les conditions prédéterminées requises ne sont pas remplies. Ainsi, en prévoyant un test, quel qu’il soit, sur des conditions prédéterminées identifiées comme étant requises pour réaliser l’estimation des paramètres du calibrage aveugle, et en prévoyant également de réutiliser les paramètres de calibrage d’une estimation précédente qui a rempli les conditions du test lorsque celui-ci n’est pas favorable pour un nouveau calibrage, on note une amélioration de ce dernier. Cette proposition est à la fois simple et astucieuse. Il est ensuite à la portée de l’homme du métier d’identifier quelles sont les conditions prédéterminées à appliquer pour réaliser le test en fonction du calibrage proposé. Par exemple, si le calibrage proposé est celui du document Halir et al précité, celles-ci peuvent porter sur la possibilité de réaliser le recalage elliptique, tandis que si c’est celui du document de brevet WO 2022/238170 A1, elles doivent plutôt porter sur la capacité à retrouver les valeurs d’amplitudes et de décalages en amplitude du système multiphasé. D’autres conditions prédéterminées peuvent être adaptées à d’autres méthodes de calibrage envisageables. De façon optionnelle, les conditions prédéterminées requises que l’unité de traitement est configurée pour tester comportent le fait que chaque signal du système multiphasé fourni résulte de suffisamment de valeurs de déphasage différentes pour qu’elles couvrent au moins la période T. De façon optionnelle également, le test réalisé par l’unité de traitement comporte une vérification que chaque signal du système multiphasé fourni varie suffisamment en amplitude pour en déduire ses valeurs extrémales pouvant être atteintes. De façon optionnelle également, les conditions prédéterminées requises que l’unité de traitement est configurée pour tester comportent le fait qu’une transformation de Clarke, Concordia ou Fortescue du système multiphasé fourni produise suffisamment de valeurs de déphasage différentes pour définir une unique solution d’ellipse par ajustement elliptique. De façon optionnelle également, l’unité de traitement est pourvue de moyens de filtrage du système multiphasé conçus pour estimer un état de ce sytème multiphasé à partir d’une mesure, ces moyens de filtrage étant plus précisément définis par le fait que : l’état à estimer comporte des valeurs constantes d’amplitudes et de décalages en amplitude des N signaux du système multiphasé, une valeur variable relative au déphasage à estimer et une valeur constante d’aide à l’estimation de la valeur variable ; et la mesure comporte les N signaux du système multiphasé tels que fournis par le coupleur multimodal ; et les paramètres de calibrage comportent les valeurs constantes d’amplitudes et de décalages en amplitude des N signaux du système multiphasé. De façon optionnelle également, l’unité de traitement est configurée pour : réaliser une transformation de Clarke, Concordia ou Fortescue du système multiphasé et en déduire des paramètres d’une ellipse ajustée sur cette transformation de Clarke, Concordia ou Fortescue ; procéder à un calibrage de cette ellipse ajustée par transformation paramétrée en cercle recentré, dans lequel les paramètres de calibrage comportent les paramètres de la transformation ; puis estimer le déphasage entre les deux signaux d’entrée après application des paramètres de calibrage à la transformation de Clarke du système multiphasé. De façon optionnelle également, l’unité de traitement est configurée pour : procéder à un calibrage des N signaux périodiques du système multiphasé à partir de valeurs d’amplitudes et de décalages en amplitude de ces N signaux périodiques, dans lequel les paramètres de calibrage comportent une combinaison de ces valeurs d’amplitudes et de décalages en amplitude ; réaliser une transformation de Clarke, Concordia ou Fortescue du système multiphasé après calibrage et en déduire des paramètres d’un cercle centré ajusté sur cette transformation de Clarke, Concordia ou Fortescue ; puis estimer le déphasage entre les deux signaux d’entrée à partir de cette transformation de Clarke, Concordia ou Fortescue. Il est également proposé un dispositif électronique de caractérisation d’un fluide, comportant : au moins un capteur conçu pour interagir avec le fluide ; et un transducteur conçu pour fournir, en interaction avec ledit au moins un capteur, une séquence d’au moins un signal de mesure représentative d’une interaction du fluide avec chaque capteur ; dans lequel le transducteur est un système de traitement d’un signal périodique par interférométrie et le dispositif électronique de caractérisation comporte au moins un dispositif d’estimation d’un déphasage entre deux signaux d’entrée issus dudit signal périodique tel que défini précédemment. De façon optionnelle, le transducteur est un système d’amplification de variation d’indice optique par interférométrie de Mach-Zehnder à interférences multimodales, le coupleur multimodal étant alors à sorties optiques et associé à un capteur optique. Il est également proposé un procédé d’estimation d’un déphasage entre deux signaux d’entrée issus d’un même signal périodique de période T traité par interférométrie, comportant : un couplage multimodal des deux signaux d’entrée pour la fourniture d’un système multiphasé de N signaux formés pour être déphasés entre eux de T/N où N ≥ 3 ; un traitement du système multiphasé pour en extraire des paramètres de calibrage et une estimation du déphasage entre les deux signaux d’entrée ; et un enregistrement des paramètres de calibrage dans des moyens de stockage ; dans lequel le traitement comporte un test sur des conditions prédéterminées requises pour réaliser une estimation des paramètres de calibrage et est conçu pour exécuter sélectivement l’une ou l’autre des opérations suivantes : l’estimation des paramètres de calibrage et leur enregistrement dans les moyens de stockage si les conditions prédéterminées requises sont remplies par les N signaux du système multiphasé ; et l’extraction, des moyens de stockage, de paramètres de calibrage d’une estimation précédente si les conditions prédéterminées requises ne sont pas remplies par les N signaux du système multiphasé. L’invention sera mieux comprise à l’aide de la description qui va suivre, donnée uniquement à titre d’exemple et faite en se référant aux dessins annexés dans lesquels : la représente schématiquement la structure générale d’un système interférométrique comportant un dispositif d’estimation de déphasage selon un mode de réalisation de l’invention, la représente schématiquement la structure générale d’un système interférométrique comportant une pluralité mutualisée de dispositifs d’estimation de déphasage tels que celui de la , la illustre les étapes successives d’un procédé interférométrique incluant une estimation de déphasage selon un mode de réalisation de l’invention, la reprend les étapes successives du procédé interférométrique de la selon une première variante de réalisation, la reprend les étapes successives du procédé interférométrique de la selon une deuxième variante de réalisation, la représente schématiquement la structure générale d’un dispositif électronique de caractérisation d’un fluide selon un mode de réalisation de l’invention, la illustre un exemple de diagramme temporel dans lequel se superposent plusieurs signaux de réponse, ou sensorgrammes, pouvant être obtenus par le dispositif électronique de la , dans des conditions maîtrisées de mesure multivariée, la illustre, sous la forme d’un diagramme circulaire, un exemple de signature normée pouvant être calculée par le dispositif électronique de la à partir de sensorgrammes tels que ceux de la , la illustre les étapes successives d’un procédé de caractérisation d’un fluide pouvant être mis en œuvre par le dispositif électronique de la , et la illustre schématiquement l’amélioration apportée par un procédé interférométrique incluant une estimation de déphasage à calibrage aveugle selon l’invention, notamment selon la première variante de réalisation de la . Le système interférométrique 10 représenté schématiquement sur la comporte un interféromètre 12, par exemple un interféromètre optique de Mach-Zehnder. Un tel interféromètre 12 reçoit un unique signal optique cohérent périodique s 0 de période T, notamment un signal optique sinusoïdal s 0 de pulsation 2πf où f est la fréquence porteuse, engendré et émis, dans un guide d’onde, par exemple par une source laser 14 de lumière cohérente. La période T peut être définie temporellement, auquel cas elle vaut l’inverse de la fréquence porteuse, ou angulairement, auquel cas elle vaut 2π pour un signal sinusoïdal. Dans l’interféromètre 12, le signal périodique s 0 est conduit vers deux branches optiques distinctes généralement de même longueur, l’une constituant une branche optique de référence 16, l’autre une branche optique d’interaction 18. La première, la branche optique de référence 16, produit un premier signal optique périodique s 1 dont le déphasage par rapport à s 0 dépend de la longueur parcourue dans cette branche optique de référence 16. La deuxième, la branche optique d’interaction 18, produit un deuxième signal optique périodique s 2 dont le déphasage par rapport à s 0 dépend non seulement de la même longueur parcourue dans cette branche optique d’interaction 18, mais en outre d’une interaction optique avec un dispositif d’interaction 20 disposé sur ce parcours. Le déphasage éventuel entre s 1 et s 2 dépend donc directement de l’interaction optique avec le dispositif 20. Pour une application de caractérisation de fluide, liquide ou gazeux, le dispositif d’interaction 20 est un capteur conçu pour interagir avec le fluide, comportant notamment un site réactif sur lequel sont greffés des composants réactifs aptes à interagir par adsorption/désorption avec des composés susceptibles d’être présents dans le fluide. Par ces interactions, l’indice de réfraction dans la branche optique d’interaction 18 est localement modifiable et variable au cours du temps, ce qui engendre le déphasage précité entre s 1 et s 2 , ainsi que sa variation au cours du temps qui produit un signal temporel de déphasage généralement appelé sensorgramme. D’autres applications sont envisageables et tout dispositif d’interaction 20 susceptible d’avoir un impact sur le déphasage du signal s 2 convient. Plus généralement, selon les applications, d’autres systèmes interférométriques peuvent convenir, autres qu’un interféromètre optique de Mach-Zehnder. Il est par exemple possible de prévoir un interféromètre de Fizeau, de Fabry-Pérot, de Jamin, de Ramsey-Bordé, de Michelson, etc., ou même une combinaison d’interféromètres. En sortie de l’interféromètre 12, les deux signaux périodiques s 1 et s 2 pourraient être recombinés pour interférer entre eux et produire un unique signal en bande de base porteur de l’information de déphasage. Mais cette dernière est alors ambiguë, de sorte qu’ils sont avantageusement fournis en tant que signaux d’entrée d’un dispositif 22 d’estimation de leur déphasage qui va maintenant être détaillé. Ce dispositif 22 comporte un coupleur multimodal 24 à deux entrées pour la réception des deux signaux d’entrée s 1 et s 2 , et à N sorties, N ≥ 3, pour la fourniture d’un système multiphasé de N signaux en bande de base formés pour être déphasés entre eux de T/N. Dans le cas de signaux optiques sinusoïdaux, dont la période T peut être définie angulairement à 2π comme indiqué précédemment, et pour un coupleur multimodal à trois sorties tel que celui illustré à titre de simple exemple non limitatif dans la , un système triphasé de trois signaux optiques sinusoïdaux s’ 1 , s’ 2 et s’ 3 déphasés entre eux de 2π/3 est fourni. Le dispositif 22 d’estimation de déphasage comporte en outre, dans le cas de la fourniture d’un système triphasé optique, un capteur optique associé au coupleur multimodal tel qu’une caméra munie de photorécepteurs, par exemple une caméra prenant la forme d’un capteur photographique CCD (de l’anglais « Charge Coupled Device ») 26, pour la fourniture d’un système triphasé électrique S de trois signaux sinusoïdaux p 1 , p 2 et p 3 déphasés entre eux de 2π/3. Ces trois signaux prennent la forme suivante : ; ; . Dans ce système triphasé S, a 1 , a 2 et a 3 sont des valeurs d’amplitudes constantes mais a priori différentes, ce qui illustre le fait que les trois signaux p 1 , p 2 et p 3 sont d’amplitudes différentes à cause des imperfections de la source laser 14, de l’interféromètre 12, du coupleur multimodal 24 et/ou du capteur CCD 26. De même, o 1 , o 2 et o 3 sont des valeurs de décalages en amplitude constantes mais a priori différentes, ce qui illustre le fait que les trois signaux p 1 , p 2 et p 3 sont de décalages en amplitude différents à cause des imperfections de la source laser 14, de l’interféromètre 12, du coupleur multimodal 24 et/ou du capteur CCD 26. Quant à θ(t), il s’agit du signal temporel de déphasage à estimer qui est codé comme s’il s’agissait de la fréquence commune des trois signaux p 1 , p 2 et p 3 du système triphasé S. Le système triphasé S est fourni en entrée d’un système informatique 28 du dispositif 22 d’estimation de déphasage, ce système informatique 28 étant spécifiquement configuré pour un traitement du système triphasé S afin d’en extraire une estimation du déphasage entre les deux signaux d’entrée s 1 et s 2 , ainsi que son évolution au cours du temps, sous la forme d’un signal temporel SG de type sensorgramme par exemple en contexte de caractérisation de fluide. Le système informatique 28 tel que représenté schématiquement sur la , comporte une unité de traitement 30 (i.e. un processeur) associée de façon classique à une mémoire 32 (par exemple une mémoire RAM ou autre). Ce système peut par exemple être mis en œuvre dans un dispositif informatique tel qu’un ordinateur classique comportant un processeur associé à une ou plusieurs mémoires pour le stockage de fichiers de données et de programmes d’ordinateurs dont les instructions sont destinées à être exécutées par le processeur. Leurs fonctions pourraient aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés. Ainsi, en variante, le dispositif informatique mettant en œuvre le système informatique 28 pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques (sans programme d’ordinateur) pour la réalisation des mêmes actions. Ce système informatique 28 peut par exemple procéder à un traitement en deux temps, i.e. calibrage puis estimation du déphasage par transformation de Clarke ou équivalente recalée, tel qu’enseigné dans le document Halir et al précité ou dans le document WO 2022/238170 A1 également précité. Ce traitement est défini et configuré sous la forme par exemple d’instructions d’un module logiciel ou programme d’ordinateur 34 enregistré en mémoire 32 et exécutable par l’unité de traitement 30. Plus précisément, le calibrage réalisé par le traitement précité est défini pour extraire directement ou indirectement, du système triphasé S lui-même, des paramètres de calibrage tels que mentionnés précédemment et pour les enregistrer en mémoire 32. Il s’agit en d’autres termes d’un calibrage aveugle. Selon l’enseignement du document Halir et al, l’unité de traitement 30 peut être configurée à l’aide du programme d’ordinateur 34 pour : réaliser une transformation de Clarke, Concordia ou Fortescue du système triphasé S et en déduire des paramètres d’une ellipse ajustée sur cette transformation de Clarke, Concordia ou Fortescue, procéder au calibrage de cette ellipse ajustée par transformation paramétrée en cercle recentré, dans lequel les paramètres de calibrage comportent les paramètres de la transformation, puis estimer le déphasage entre les deux signaux d’entrée après application des paramètres de calibrage à la transformation de Clarke du système triphasé S. Selon l’enseignement du document WO 2022/238170 A1, l’unité de traitement 30 peut être configurée à l’aide du programme d’ordinateur 34 pour : procéder directement au calibrage des N signaux périodiques du système triphasé S à partir de valeurs d’amplitudes et de décalages en amplitude de ces N signaux périodiques, dans lequel les paramètres de calibrage comportent une combinaison de ces valeurs d’amplitudes et de décalages en amplitude ; réaliser une transformation de Clarke, Concordia ou Fortescue du système triphasé S après calibrage et en déduire des paramètres d’un cercle centré ajusté sur cette transformation de Clarke, Concordia ou Fortescue, puis estimer le déphasage entre les deux signaux d’entrée à partir de cette transformation de Clarke, Concordia ou Fortescue. Mais en variante et de façon avantageuse, le système informatique 28 comporte des moyens de filtrage du système triphasé S pour estimer de façon récursive un état de ce système à partir d’une mesure. Il s’agit donc d’un filtrage de type Kalman. Ces moyens de filtrage peuvent également être définis et configurés sous la forme d’instructions du module logiciel ou programme d’ordinateur 34 enregistré en mémoire 32 et exécutable par l’unité de traitement 30. Ils sont plus précisément astucieusement configurés de sorte que l’état à estimer comporte les valeurs constantes d’amplitudes a 1 , a 2 et a 3 et de décalages en amplitude o 1 , o 2 et o 3 des trois signaux sinusoïdaux p 1 , p 2 et p 3 du système triphasé S, une variable θ relative au déphasage θ(t) à estimer et une valeur constante δθ d’aide à l’estimation et au suivi de la variable θ, tandis que la mesure comporte les trois signaux sinusoïdaux du système triphasé S tels que fournis par le coupleur multimodal 24 et le capteur CCD 26. En optimisant le modèle d’état par comparaison avec la mesure selon un critère de minimisation d’erreur, par exemple la minimisation d’une erreur quadratique moyenne de façon bien connue en filtrage récursif, le calibrage aveugle se réalise, par estimation des paramètres constants a 1 , a 2 et a 3 , o 1 , o 2 et o 3 et δθ, en même temps que l’estimation de déphasage, par suivi temporel de la variable θ. En tant que paramètre constant d’aide à l’estimation et au suivi de θ, δθ est destiné à converger vers la valeur constante égale à 0. Comme indiqué précédemment, ce suivi temporel tel que modélisé par filtrage assure une continuité qui permet de se passer d’un déroulement de phase. Expérimentalement, on observe que l’estimation du déphasage θ présente de bien meilleurs résultats par le filtrage récursif proposé que dans l’état de la technique où le calibrage est réalisé préalablement à l’estimation et séparément. Les résultats sont par ailleurs d’autant meilleurs que l’initialisation des paramètres a 1 , a 2 , a 3 , o 1 , o 2 , o 3 , θ et δθ est pertinente. C’est le cas lorsque l’on dispose déjà de valeurs a 1 , a 2 et a 3 , o 1 , o 2 , o 3 résultant d’un précédent calibrage/estimation et lorsque l’on initialise θ et δθ à 0. Par défaut, les paramètres a 1 , a 2 , a 3 , o 1 , o 2 et o 3 peuvent être initialisés à des valeurs prédéfinies quelconques. Concrètement et à titre d’exemple non limitatif, les moyens de filtrage peuvent se présenter sous la forme d’un filtre à deux fonctions de transfert : une fonction de transition F calculant une valeur prédite à partir d’une valeur déjà estimée de l’état du système multiphasé, et une fonction d’observation H calculant une estimation de la mesure à partir de la valeur prédite ou déjà estimée de l’état du système multiphasé. Ce filtre est alors avantageusement défini d’une façon générale de sorte que : l’état du système multiphasé se présente sous la forme d’un vecteur dont les coefficients comportent les 2N valeurs constantes d’amplitudes et de décalages en amplitude des N signaux sinusoïdaux du système multiphasé, la valeur variable relative au déphasage à estimer et la valeur constante d’aide à l’estimation de la valeur variable, c’est-à-dire par exemple un vecteur à 2N+2 coefficients, et la mesure se présente sous la forme d’un vecteur à N coefficients dont les N coefficients comportent respectivement les N signaux sinusoïdaux du système multiphasé tels que fournis par le coupleur multimodal 24 et le capteur CCD 26. Pour simplifier les notations, dans le cas N = 3 du système triphasé S échantillonné temporellement, le vecteur d’état à un instant d’échantillonnage k peut s’écrire sous la forme : La fonction de transition F permet donc de modéliser le lien entre X k et X k+1 , par exemple selon l’équation générale : , où V k est un vecteur modélisant un bruit de traitement. Plus précisément, un modèle simple de lien entre X k et X k+1 peut s’exprimer de la façon suivante : , où V k est de mêmes dimensions que X k et est un vecteur de bruit blanc gaussien. Avantageusement, au vu de la définition qui est donnée du vecteur d’état dans la présente invention, la fonction de transition est définie de la façon suivante : . Il en résulte que la valeur prédite X k+1 du vecteur d’état est égale à la somme de la valeur déjà estimée X k du vecteur d’état, dans laquelle le coefficient comportant la valeur variable relative au déphasage à estimer est mis à jour par ajout de la valeur du coefficient comportant la valeur constante d’aide à l’estimation de la valeur variable, et du vecteur de bruit de traitement V k . La fonction d’observation H, calculant une estimation de la mesure à partir de la valeur prédite ou déjà estimée de l’état du système triphasé S, permet de modéliser le lien entre X k ou X k+1 et cette estimation Y k+1 de la mesure, par exemple selon l’équation générale suivante dans le cas d’une estimation de la mesure à partir de la valeur déjà estimée de l’état du système triphasé S : , où N k est un vecteur modélisant un bruit de mesure. Plus précisément, un modèle simple de lien entre X k et Y k+1 peut s’exprimer de la façon suivante : , où N k est de mêmes dimensions que Y k +1 et est un vecteur de bruit blanc gaussien. Avantageusement et selon un premier mode de réalisation dans lequel les moyens de filtrage 34 comportent un filtre de Kalman sans parfum UKF (de l’anglais « Unscented Kalman Filter »), la fonction d’observation est définie de la façon suivante pour le système triphasé S : . Il en résulte que l’estimation de la mesure Y k+1 est égale à la somme d’un vecteur à N=3 coefficients, dont les trois coefficients sont respectivement les trois signaux du système triphasé S exprimés en fonction des coefficients de la valeur déjà estimée X k du vecteur d’état, et du vecteur de bruit de mesure N k . Avantageusement également et selon un deuxième mode de réalisation dans lequel les moyens de filtrage 34 comportent un filtre de Kalman étendu EKF (de l’anglais « Extended Kalman Filter »), la fonction d’observation peut être définie de la façon suivante pour le système triphasé S : , où . Il en résulte que l’estimation de la mesure Y k+1 est égale à la somme du produit matriciel d’une matrice jacobienne des trois signaux du système triphasé S exprimés en fonction des coefficients de la valeur déjà estimée X k du vecteur d’état avec la valeur déjà estimée X k du vecteur d’état, et du vecteur de bruit de mesure N k . Conformément aux principes généraux du filtrage récursif, l’estimation des paramètres de l’état du système triphasé S, c’est-à-dire l’estimation des coefficients du vecteur d’état dans le cas d’une représentation vectorielle et matricielle du filtrage, et le suivi de leur éventuelle évolution au cours du temps se fait par comparaison de la mesure et de son estimation, plus précisément par minimisation d’une erreur d’estimation en adaptant les paramètres en conséquence. Il en résulte que les paramètres permettant le calibrage automatique des mesures fournies sont estimés conjointement au paramètre de déphasage, d’où la fourniture d’un moyen simple et efficace pour fusionner le calibrage et l’estimation de déphasage, cette dernière produisant le signal temporel SG de type sensorgramme par exemple en contexte de caractérisation de fluide. On notera que le filtrage précité prévoit également que les paramètres de calibrage a 1 , a 2 et a 3 , o 1 , o 2 , o 3 du système triphasé S, estimés en même temps que le déphasage θ, soient enregistrés en mémoire 32. Conformément aux principes généraux de la présente invention, et quel que soit le mode de réalisation choisi pour le calibrage et l’estimation du déphasage entre les deux signaux d’entrée s 1 et s 2 , notamment parmi les trois enseignements précités, l’unité de traitement 30 exécutant le programme d’ordinateur 34 est configurée pour tester si les N signaux d’un système multiphasé remplissent des conditions prédéterminées requises pour réaliser l’estimation des paramètres de calibrage et pour : réaliser l’estimation des paramètres de calibrage et leur enregistrement en mémoire 32 si les conditions prédéterminées requises sont remplies, et extraire, de la mémoire 32, les paramètres de calibrage d’une estimation précédente si les conditions prédéterminées requises ne sont pas remplies. Par exemple, si le calibrage est réalisé selon l’enseignement du document Halir et al, c’est-à-dire par recalage sur un cercle centré de l’ellipse obtenue par transformation de Clarke ou équivalente d’un système multiphasé, les conditions prédéterminées peuvent porter sur la possibilité supposée de réaliser un ajustement elliptique fiable permettant justement d’obtenir une telle ellipse. Il est à la portée de l’homme du métier d’en concevoir une implémentation. Mais selon un mode de réalisation avantageux et astucieux, ces conditions prédéterminées requises peuvent comporter un nombre minimal prédéfini de valeurs de déphasage suffisamment différentes produites par transformation de Clarke ou équivalente du système multiphasé, par exemple selon un modèle prédéfini de distance entre elles et/ou de partitionnement, pour définir et obtenir une unique solution d’ellipse par ajustement elliptique. En d’autres termes, si ce nombre minimal est atteint, alors l’ajustement elliptique est considéré comme possible, fiable et de solution unique, de sorte que le calibrage par ajustement elliptique et recalage sur un cercle centré, l’enregistrement en mémoire 32 des paramètres de calibrage ainsi obtenus, puis l’estimation de déphasage peuvent être réalisés sur la base du système multiphasé fourni. En revanche, si ce nombre minimal n’est pas atteint, les paramètres de calibrage d’un calibrage précédent sont extraits de la mémoire 32 et appliqués au système multiphasé fourni pour réaliser l’estimation de déphasage par calcul Arctan et déroulement de phase de la transformation de Clarke du système multiphasé. Par exemple également, si le calibrage est réalisé selon l’enseignement du document WO 2022/238170 A1 ou du filtrage récursif détaillé ci-dessus, c’est-à-dire directement à partir des valeurs d’amplitudes et de décalages en amplitude d’un système multiphasé, les conditions prédéterminées peuvent porter sur la capacité supposée à retrouver ces valeurs. Il est à la portée de l’homme du métier d’en concevoir une implémentation. Mais selon un mode de réalisation avantageux et astucieux, ces conditions prédéterminées requises peuvent comporter le fait que chaque signal du système multiphasé fourni résulte de suffisamment de valeurs de déphasage différentes, par exemple selon un modèle prédéfini de distance entre elles et/ou de partitionnement, pour qu’elles couvrent au moins la période T, c’est-à-dire la période 2π dans le cas de signaux sinusoïdaux. Il est alors possible de reconstituer, au moins par interpolation et pour chaque signal du système multiphasé, un unique signal sinusoïdal passant au plus près de ce signal sur au moins une période angulaire parcourue, de sorte qu’il devient possible d’en déduire ses valeurs extrémales, donc les valeurs de son amplitude et de son décalage en amplitude. En variante, les conditions prédéterminées requises peuvent comporter le fait que chaque signal du système multiphasé fourni varie suffisamment rapidement en amplitude pour en déduire ses valeurs extrémales pouvant être atteintes. En d’autres termes, cela signifie qu’il comporte une portion temporelle de forme sinusoïdale sur au moins une période angulaire permettant, au moins par interpolation sinusoïdale également, d’estimer ses valeurs extrémales. Ainsi, si au moins l’interpolation sinusoïdale est possible, alors l’estimation des valeurs extrémales à partir desquelles on peut estimer l’amplitude et le décalage en amplitude de chaque signal du système multiphasé fourni est considérée comme possible, fiable et de solution unique, de telle sorte que le calibrage, l’enregistrement en mémoire 32 des paramètres de calibrage ainsi obtenus, puis l’estimation de déphasage peuvent être réalisés sur la base du système multiphasé fourni (soit successivement dans le cas du document WO 2022/238170 A1, soit simultanément dans le cas du filtrage récursif détaillé ci-dessus). En revanche, si l’estimation des valeurs extrémales des signaux du système multiphasé fourni n’est pas possible, les paramètres de calibrage d’un calibrage précédent sont extraits de la mémoire 32 et appliqués au système multiphasé fourni pour réaliser l’estimation de déphasage : par transformation de Clarke du système multiphasé, calcul Arctan et déroulement de phase dans le cas du document WO 2022/238170 A1 ; par filtrage récursif en imposant les paramètres d’amplitudes et de décalages en amplitude du calibrage précédent pour le filtrage récursif détaillé ci-dessus. Le système interférométrique 40 représenté schématiquement sur la comporte une pluralité mutualisée de dispositifs d’estimation de déphasage tels que celui de la . Les éléments communs repris du système interférométrique 10 portent la même référence. Ainsi, ce système 40 comporte une pluralité d’interféromètres 12 1 … 12 n qui reçoivent tous le signal optique sinusoïdal s 0 engendré et émis par la source laser 14 et fournissent chacun deux signaux d’entrée, notés respectivement s 1,1 , s 1,2 … s n,1 , s n,2 , à l’un d’une pluralité de dispositifs 22 1 … 22 n d’estimation de leurs déphasages respectifs, avec n ≥ 2. Chaque interféromètre 12 i , 1≤i≤n, comporte avantageusement un dispositif d’interaction qui lui est propre dans sa branche optique d’interaction de sorte que les n déphasages obtenus sont respectivement propres aux n interféromètres 12 1 … 12 n . Les n dispositifs 22 1 … 22 n d’estimation de déphasages comportent respectivement n coupleurs multimodaux 24 1 … 24 n à deux entrées, pour la réception des n couples de signaux d’entrée s 1,1 , s 1,2 … s n,1 , s n,2 , et à trois sorties pour la fourniture respective de n systèmes triphasés de trois signaux sinusoïdaux s’ 1,1 , s’ 1,2 , s’ 1, 3 … s’ n,1 , s’ n,2 , s’ n, 3 . Dans le cas de la fourniture de systèmes triphasés optiques, les n dispositifs 22 1 … 22 n d’estimation de déphasages partagent le même capteur photographique CCD 26 associé aux n coupleurs multimodaux 24 1 … 24 n qui lui fournissent les n systèmes triphasés de trois signaux sinusoïdaux s’ 1,1 , s’ 1,2 , s’ 1, 3 … s’ n,1 , s’ n,2 , s’ n, 3 , pour la fourniture d’un ensemble S de n systèmes triphasés électriques dans chacun desquels les trois signaux sinusoïdaux sont déphasés entre eux de 2π/3. Ces 3n signaux prennent la forme suivante : ; ; ; … ; ; . Les n dispositifs 22 1 … 22 n d’estimation de déphasages partagent enfin le même système informatique 28 spécifiquement configuré pour un traitement de l’ensemble S des n systèmes triphasés fournis afin d’en extraire une estimation des n déphasages respectifs entre les couples de signaux d’entrée s 1,1 , s 1,2 … s n,1 , s n,2 , ainsi que leurs évolutions respectives au cours du temps, sous la forme d’un ensemble de n signaux temporels SG i , 1≤i≤n, par exemple constitué de n sensorgrammes en contexte de caractérisation de fluide. Le fonctionnement du système informatique 28 de la est inchangé par rapport à ce qui a été détaillé en référence à la , moyennant une simple adaptation pour traiter simultanément les n systèmes triphasés fournis au lieu d’un seul. C’est la raison pour laquelle il porte la même référence 28 et comporte la même unité de traitement 30 associée à la mémoire 32 dans laquelle sont enregistrés les moyens de traitement précités (conformément à l’enseignement du document Halir et al, du document WO 2022/238170 A1 ou du filtrage récursif détaillé ci-dessus) sous la forme du module logiciel ou programme d’ordinateur 34 apte à s’exécuter cette fois-ci sur l’ensemble S des n systèmes triphasés fournis, ainsi que les paramètres de calibrage estimés. Le fonctionnement du système interférométrique 10, ou celui équivalent du système interférométrique 40, va maintenant être détaillé d’une façon générale en référence à la . Au cours d’une étape préalable d’interférométrie 100, le signal optique sinusoïdal s 0 engendré et émis par la source laser 14 traverse l’interféromètre 12 ou la pluralité d’interféromètres 12 1 … 12 n ., pour la fourniture des deux signaux d’entrée s 1 , s 2 ou des 2n signaux d’entrée s 1,1 , s 1,2 … s n,1 , s n,2 . Au cours d’une première étape 102 de procédé d’estimation d’un déphasage entre les deux signaux d’entrée s 1 , s 2 ou entre les 2n signaux d’entrée s 1,1 , s 1,2 … s n,1 , s n,2 deux à deux, ces signaux d’entrée subissent un couplage multimodal, éventuellement combiné avec une capture CCD, pour l’obtention du système triphasé S de trois signaux sinusoïdaux déphasés entre eux de 2π/3 (ou plus généralement d’un système multiphasé de N signaux déphasés entre eux de T/N où N ≥ 3) ou de l’ensemble S de n systèmes triphasés dans chacun desquels les trois signaux sinusoïdaux sont déphasés entre eux de 2π/3 (ou plus généralement n systèmes multiphasés dans chacun desquels les N signaux sont déphasés entre eux de T/N où N ≥ 3). Au cours d’une étape de test suivante 104, l’unité de traitement 30 exécute le programme d’ordinateur 34 pour tester si les conditions prédéterminées requises pour réaliser une estimation des paramètres de calibrage à partir du système triphasé S ou de l’ensemble S de n systèmes triphasés fourni sont remplies. Ces conditions prédéterminées ont été détaillées précédemment et ne seront pas rappelées ici. On notera simplement que l’étape 104 peut nécessiter un traitement du système triphasé S ou de l’ensemble S de n systèmes triphasés tel qu’une transformation de Clarke ou équivalente dans le cas où l’enseignement du document Halir et al est suivi pour réaliser le calibrage et l’estimation de phase. Si les conditions prédéterminées sont remplies à l’issue de ce test 104, le procédé d’estimation de déphasage passe : à une étape 106 de calibrage, incluant l’estimation des paramètres de calibrage, et d’estimation de déphasage(s) entre les deux signaux d’entrée s 1 , s 2 ou entre les 2n signaux d’entrée s 1,1 , s 1,2 … s n,1 , s n,2 deux à deux, et à une étape 108 d’enregistrement en mémoire 32 des paramètres de calibrage estimés à l’étape 106. Etant donnée que l’étape 108 dépend du succès de l’étape de test 104 et du fait que l’étape 106 procède effectivement à l’estimation des paramètres de déphasage(s) par traitement direct ou indirect du système triphasé S ou de l’ensemble S de n systèmes triphasés fourni, une fonction logique « ET » symbolisée par la référence « & » s’interpose entre les étapes 106 et 108 de la pour illustrer la vérification de ces deux conditions. Si les conditions prédéterminées ne sont pas remplies à l’issue du test 104, le procédé d’estimation de déphasage passe à une deuxième étape de test 110 au cours de laquelle il est vérifié si des paramètres de calibrage d’une estimation précédente sont disponibles en mémoire 32. Si c’est le cas, le procédé passe à une étape 112 d’extraction de ces paramètres de calibrage de la mémoire 32 et le procédé passe à l’étape 106 au cours de laquelle l’estimation des paramètres de calibrage par traitement direct ou indirect du système triphasé S ou de l’ensemble S de n systèmes triphasés fourni est remplacée par une application des paramètres de calibrage extraits à l’étape 112. Sinon, le procédé passe directement à l’étape 106 au cours de laquelle l’estimation des paramètres de calibrage par traitement direct ou indirect du système triphasé S ou de l’ensemble S de n systèmes triphasés fourni est remplacée par une application de paramètres de calibrage définis par défaut pour un fonctionnement en mode dégradé, par exemple un fonctionnement sans calibrage lorsque les paramètres de calibrage définis par défaut sont des valeurs arithmétiques neutres. En variante, si aucun paramètre de calibrage n’est disponible en mémoire 32 lors de l’étape de test 106, le procédé d’estimation de déphasage pourrait être stoppé. De façon plus spécifique, lorsque le calibrage du système triphasé S ou de l’ensemble S de n systèmes triphasés fourni et l’estimation de déphasage(s) sont séparés, c’est-à-dire notamment lorsque l’enseignement du document Halir et al ou du document WO 2022/238170 A1 est appliqué, le fonctionnement de la peut être précisé conformément au schéma de la . Sur cette figure, les étapes 100, 102, 104, 106, 108, 110 et 112 de la restent inchangées, mais l’étape 106 de calibrage et d’estimation de déphasage(s) peut être subdivisée en trois sous-étapes 106A, 106B et 106C. La sous-étape 106A n’est exécutée que si les conditions prédéterminées sont remplies à l’issue de l’étape de test 104. Selon l’enseignement du document Halir et al, elle peut consister à déduire, de la ou les transformation(s) de Clarke ou équivalente(s) réalisée(s) sur le système triphasé S ou sur l’ensemble S de n systèmes triphasés à l’étape de test 104, les paramètres d’une ou plusieurs ellipse(s) ajustée(s) sur cette ou ces transformation(s), puis à procéder à un calibrage de cette ou ces ellipse(s) ajustée(s) par transformation(s) paramétrée(s) en cercle(s) recentré(s), dans lequel les paramètres de calibrage comportent les paramètres de la ou des transformation(s). En variante, selon l’enseignement du document WO 2022/238170 A1, elle peut consister à procéder à un calibrage direct des trois signaux périodiques du système triphasé, ou des 3n signaux périodiques de l’ensemble S de n systèmes triphasés, à partir des valeurs d’amplitudes et de décalages en amplitude de ces trois ou 3n signaux périodiques, dans lequel les paramètres de calibrage comportent une combinaison de ces valeurs d’amplitudes et de décalages en amplitude. La sous-étape 106B est exécutée suite à l’étape 106A ou suite à l’étape 112. Elle consiste à appliquer les paramètres de calibrage, obtenus respectivement à la sous-étape 106A ou à l’étape 112, à l’ellipse ou aux ellipses ajustée(s) (selon Halir et al), ou au système triphasé S ou à l’ensemble S de n systèmes triphasés (selon WO 2022/238170 A1). La sous-étape 106C est exécutée suite à l’étape 106B ou suite à l’étape 110. Selon l’enseignement du document Halir et al, elle consiste à estimer le ou les déphasage(s) sur la ou les ellipse(s) ajustée(s) et calibrée(s), puis à procéder au(x) déroulement(s) de phase(s). Selon l’enseignement du document WO 2022/238170 A1, elle consiste à procéder à une transformation de Clarke ou équivalente du système triphasé S ou de l’ensemble S de n systèmes triphasés calibré, à estimer le(s) déphasage(s) sur le(s) cercle(s) ajusté(s), puis à procéder au(x) déroulement(s) de phase(s). De façon plus spécifique également, lorsque le calibrage du système triphasé S ou de l’ensemble S de n systèmes triphasés fourni et l’estimation de déphasage(s) sont simultanés, c’est-à-dire notamment lorsque l’enseignement du filtrage récursif détaillé ci-dessus est appliqué, le fonctionnement de la peut être précisé conformément au schéma de la . Sur cette figure, les étapes 100, 102, 104, 106, 108, 110 et 112 restent inchangées, mais l’étape 106 de calibrage et d’estimation de déphasage(s) est exécutée conformément à une variante 106’ ou 106’’ selon que les conditions prédéterminées requises sont remplies ou pas à l’issue de l’étape de test 104. Conformément aux deux variantes, le système triphasé ou ensemble de n systèmes triphasés (ou plus généralement multiphasé(s)) S est filtré pour en extraire une estimation du déphasage entre les deux signaux d’entrée s 1 , s 2 ou entre les 2n signaux d’entrée s 1,1 , s 1,2 … s n,1 , ce filtrage étant plus précisément du type à estimer un état de ce (ou chaque) système triphasé (ou plus généralement multiphasé) à partir d’une mesure, ce filtrage étant plus précisément défini par le fait que : l’état à estimer comporte des valeurs constantes d’amplitudes et de décalages en amplitude des trois signaux du (ou de chaque) système triphasé (ou plus généralement des N signaux du (ou de chaque) système multiphasé), une variable relative au déphasage à estimer et une valeur constante d’aide à l’estimation de la valeur variable ; et la mesure comporte les trois signaux du (ou de chaque) système triphasé (ou plus généralement les N signaux du (ou de chaque) système multiphasé) tels que fournis par le couplage multimodal 102. Selon la variante 106’, les paramètres de calibrage sont estimés en même temps que l’estimation de déphasage(s) par filtrage récursif. Ils sont par exemple initialisés à des valeurs par défaut. Mais ils peuvent aussi être initialisés à des valeurs de paramètres de calibrage d’une estimation précédente s’il en existe en mémoire 32. Selon la variante 106’’, les paramètres de calibrage sont imposés à des valeurs par défaut, suite à l’étape 110, ou à des valeurs de paramètres de calibrage d’une estimation précédente, suite à l’étape 112, ce qui apporte une contrainte sur le filtrage récursif réalisé, celui-ci n’étant alors plus autorisé à faire converger les valeurs constantes d’amplitudes et de décalages en amplitude des trois signaux du (ou de chaque) système triphasé (ou plus généralement des N signaux du (ou de chaque) système multiphasé) vers d’autres valeurs que celles des paramètres de calibrage imposés. Bien entendu, l’étape 108 ne peut être ensuite réalisée que si c’est la variante 106’ de l’étape 106 qui est exécutée. Le système interférométrique 40 de la peut être implémenté dans un dispositif électronique 50 de caractérisation olfactive d’un fluide tel que celui représenté schématiquement sur la . Il s’agit d’un exemple non limitatif de dispositif de caractérisation de fluide selon un mode de réalisation possible de la présente invention pour une application olfactive non limitative d’identification d’odeurs. Il comporte une chambre de mesure 52 destinée à recevoir un fluide, par exemple un gaz tel que de l’air ambiant. Pour ce faire, il comporte un dispositif d’aspiration 54 conçu pour aspirer l’air situé à l’intérieur de la chambre de mesure 52 et le faire sortir à l’extérieur. Il comporte en outre une entrée d’air 56 pouvant être sélectivement fermée pour garder l’air ambiant dans la chambre de mesure 52 ou bien ouverte pour permettre l’évacuation de l’air ambiant de la chambre de mesure 52 et son renouvellement par l’activation du dispositif d’aspiration 54. Il est ainsi doté de moyens de maîtrise des flux entrant et sortant. Dans sa chambre de mesure 52, Le dispositif électronique de caractérisation olfactive 50 comporte plusieurs capteurs, notamment des capteurs olfactifs 58, répartis respectivement sur autant de sites réactifs, par exemple une soixantaine, conçus pour interagir avec des composés susceptibles d’être présents dans la chambre de mesure 52 lorsque le dispositif 50 est placé à proximité d’un fluide à analyser émetteur de ces composés, notamment lorsque l’entrée d’air 56 est à proximité du fluide considéré. Les composés émis sont généralement des composés organiques volatils mais la présente invention ne se limite pas à de tels composés. Chaque capteur olfactif 58 est lui-même par exemple un biocapteur conçu pour interagir avec les composés d’une famille particulière de composés organiques volatils. En pratique, chaque capteur olfactif 58 peut comporter une molécule, telle qu’un peptide immobilisé sur un substrat ou un polymère recouvrant une surface, complémentaire des composés de la famille associée à ce capteur olfactif 58. L’idée générale est de fonctionnaliser des sites réactifs à l’aide de capteurs olfactifs 58 (i.e. biocapteurs, polymères, nanotubes de carbones, etc.) de telle sorte qu’ils adsorbent et désorbent de manière différenciée les composés organiques volatils, de former une réponse d’interaction moléculaire différenciée des capteurs olfactifs, et d’amplifier la réponse sous la forme d’une séquence S de signaux électriques de mesure à l’aide d’un dispositif de transduction physique. En variante, le dispositif électronique de caractérisation olfactive 50 pourrait être adapté pour être mis en contact avec tout autre fluide, liquide ou gazeux, que de l’air ambiant. Il pourrait également selon une version particulièrement simple ne pas comporter le dispositif d’aspiration 54 et l’entrée d’air 56, voire même la chambre de mesure 52. Dans cette version simple, les capteurs olfactifs 58 sont alors susceptibles d’être directement mis au contact du fluide à analyser sans maîtrise de flux. Les capteurs olfactifs 58 sont associés à au moins un transducteur 60 avec lequel ils interagissent. Ce transducteur 60 est arrangé et configuré pour mesurer tout changement de propriété physique provoqué par une interaction des capteurs olfactifs 18 avec le fluide à analyser. Il fournit le signal électrique S de mesure qui se présente par exemple sous la forme d’une séquence de signaux électriques de mesure et caractérise ce fluide puisque cette séquence est représentative des composés organiques volatils avec lesquels les capteurs olfactifs 58 peuvent interagir dans la chambre de mesure 52. Plus précisément, le transducteur 60 est un système de traitement d’un signal périodique, i.e. le signal périodique s 0 précédemment défini dans le système interférométrique 40, par interférométrie. Il comporte plus précisément la pluralité d’interféromètres 12 1 … 12 n qui reçoivent tous le signal périodique s 0 et présentent chacun l’un des capteurs 58 en tant que dispositif d’interaction de leur branche d’interaction. Il y a donc autant d’interféromètres 12 1 … 12 n dans le transducteur 60 que de capteurs 58. A titre d’exemple non limitatif, le transducteur 60 est comme évoqué précédemment un système d’amplification de variation d’indice optique par interférométrie de Mach-Zehnder, par exemple selon une technologie à matrice d’interféromètres de Mach-Zehnder, dite technologie MZI (de l’anglais « Mach-Zehnder Interferometer »), et même plus précisément une technologie MZI à interférences multimodales, dite technologie MZI/MMI (de l’anglais « Multi Mode Interference »), puisque le transducteur 60 comporte en outre les n coupleurs multimodaux 24 1 … 24 n à deux entrées, pour la réception des n couples de signaux d’entrée s 1,1 , s 1,2 … s n,1 , s n,2 , et à trois sorties pour la fourniture respective de n systèmes triphasés de trois signaux sinusoïdaux s’ 1,1 , s’ 1,2 , s’ 1,3 … s’ n,1 , s’ n,2 , s’ n,3 , ainsi que le capteur photographique CCD 26 mutualisé pour fournir l’ensemble S des n systèmes triphasés en tant que signal électrique de mesure. Ce transducteur 60 est donc configuré pour mesurer tout changement d’un indice de réfraction dû à une interaction du fluide étudié avec au moins l’un quelconque des capteurs olfactifs 58 grâce au décalage de phases détectable entre le bras de référence et le bras de détection de l’interféromètre sur lequel ce capteur olfactif quelconque est disposé. Le dispositif électronique de caractérisation olfactive 50 comporte en outre plusieurs modules fonctionnels qui vont être décrits ci-dessous. Dans l’exemple décrit, ces modules sont de nature logicielle. Ainsi, le dispositif 50 comporte le système informatique 28 précédemment décrit, celui-ci comportant l’unité de traitement 30 et la mémoire associée 32 dans laquelle plusieurs programmes d’ordinateurs ou plusieurs fonctions d’un même programme d’ordinateur sont enregistrés, dont notamment le module logiciel ou programme d’ordinateur 34. Ces programmes d’ordinateurs comportent des instructions conçues pour être exécutées par l’unité de traitement 30 afin de réaliser les fonctions des modules logiciels. Ils sont présentés comme distincts, mais cette distinction est purement fonctionnelle. Ils pourraient tout aussi bien être regroupés selon toutes les combinaisons possibles en un ou plusieurs logiciels. Leurs fonctions pourraient aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés, tels que des circuits numériques. Ainsi, en variante, le système informatique 28 pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques (sans programme d’ordinateur) pour la réalisation des mêmes fonctions. Le dispositif électronique de caractérisation olfactive 50 comporte ainsi tout d’abord un module logiciel 62, destiné à être exécuté par l’unité de traitement 30, de commande du dispositif d’aspiration 54 (s’il est prévu dans de dispositif 50), de l’entrée d’air 56 (si elle est également prévue dans de dispositif 50) et du transducteur 60. Il comporte en outre de façon optionnelle mais avantageuse un module logiciel 64, destiné à être exécuté par l’unité de traitement 30, de sélection, parmi les capteurs olfactifs 58 du dispositif électronique de caractérisation olfactive 50, d’un sous-ensemble de capteurs sensibles à des composants volatils caractéristiques d’une empreinte olfactive recherchée. Ces composants volatils caractéristiques peuvent varier d’une application ou d’un fluide étudié à l’autre de sorte que la sélection de capteurs olfactifs réalisée par le module logiciel 64 peut varier elle aussi et être paramétrée. Le sous-ensemble sélectionné comporte par exemple M ≥ 1 capteur(s) olfactif(s), notamment avantageusement plusieurs capteurs olfactifs (M ≥ 2). Le dispositif électronique de caractérisation olfactive 50 comporte en outre le module logiciel ou programme d’ordinateur 34, déjà défini et décrit en détail, destiné à être exécuté par l’unité de traitement 30, pour extraire M sensorgrammes respectivement représentatifs des interactions des M capteurs olfactifs sélectionnés avec les composés organiques volatils concernés à partir des valeurs propres à ces M capteurs olfactifs dans la séquence S de systèmes triphasés fournie par le transducteur 60. La figure 7A illustre ainsi les diagrammes temporels superposés d’une soixantaine de sensorgrammes obtenus sur une durée d’environ 190 secondes selon un protocole de mesure bien maîtrisé, impliquant une commande du dispositif d’aspiration 54 et de l’entrée d’air 56, dans lequel : les capteurs olfactifs 58 sont tout d’abord exposés à un environnement fluidique de référence à fluide porteur sans présence des composés cibles d’un fluide à analyser au cours d’un premier état de référence identifiable par une première portion PH1 des sensorgrammes, ils sont ensuite exposés au fluide à analyser au cours d’un deuxième état analytique d’adsorption déclenché par une injection maîtrisée de ce fluide dans la chambre de mesure 52, ce deuxième état étant identifiable par une deuxième portion PH2 des sensorgrammes, et ils sont finalement exposés de nouveau à l’environnement fluidique de référence au cours d’un troisième état final de désorption par une évacuation maîtrisée du fluide à analyser hors de la chambre de mesure 52, ce troisième état étant identifiable par une troisième portion PH3 des sensorgrammes. De retour à la figure 6, le dispositif électronique de caractérisation olfactive 50 comporte en outre un module logiciel 66 optionnel, destiné à être exécuté par l’unité de traitement 30, pour opérer un traitement préalable éventuel sur les M sensorgrammes fournis par le module logiciel 34. Ce traitement préalable comporte par exemple un filtrage passe-bas implémenté sous la forme d’un filtre numérique à réponse impulsionnelle finie ou infinie. Il s’agit de filtrer le bruit de mesure à hautes fréquences dans les signaux bruts tels que fournis par le module logiciel 34. Un filtre de Butterworth à réponse impulsionnelle finie du premier ordre et fréquence de coupure normalisée à 0,45 (i.e. valeur du ratio entre la fréquence de coupure et la fréquence d’échantillonnage égale à 0,45) convient. Ce traitement préalable comporte par exemple en outre un calcul de norme au sens du document de brevet WO 2020/141281 A1 sur les M sensorgrammes filtrés ou non pour l’obtention de M sensorgrammes filtrés et/ou normés . Le dispositif électronique de caractérisation olfactive 50 comporte en outre un module logiciel 68, destiné à être exécuté par l’unité de traitement 30, pour obtenir de façon bien connue et non détaillée une caractérisation ou signature SIG de la composition du fluide à analyser à partir des M sensorgrammes ou . Cette caractérisation ou signature SIG peut prendre la forme d’une signature olfactive normée telle qu’illustrée dans la sous la forme d’un diagramme circulaire. On notera que ce module peut procéder en deux temps : tout d’abord l’obtention d’une première signature intermédiaire, puis la transformation de cette première signature intermédiaire par normalisation. La signature SIG illustrée dans la comporte soixante-quatre composantes. Le module logiciel 68 peut donc optionnellement mais avantageusement être en outre programmé pour transformer la signature SIG, normée ou non, en une autre signature SIG’ simplifiée par réduction de composantes. Une méthode par analyse discriminante linéaire LDA (de l’anglais « Linear Discriminant Analysis »), analyse en composantes principales PCA (de l’anglais « Principal Component Analysis »), analyse en composantes indépendantes ICA (de l’anglais « Independent Component Analysis »), auto-encodeur, etc., convient. Pour une signature SIG à soixante-quatre composantes telle que celle de la , une signature SIG’ simplifiée à deux ou trois composantes peut être obtenue. Un procédé de caractérisation olfactive d’un fluide correspondant à l’exécution des modules logiciels 34 et 62 à 68 va maintenant être détaillé conformément à la succession d’étapes de la . Au cours d’une première étape optionnelle 200 réalisée par exécution du module logiciel 64, une partie des capteurs olfactifs 58 est sélectionnée en fonction d’une empreinte olfactive recherchée. Au cours d’une étape suivante 202, le dispositif électronique de caractérisation olfactive 50 est approché d’un fluide de manière à obtenir une interaction du fluide avec chacun des capteurs olfactifs 58 sélectionnés. A cette occasion, le module logiciel 62 est exécuté. Au cours d’une étape suivante 204, le transducteur 60 fournit une séquence S de signaux de mesure en interaction avec les capteurs olfactifs 58. Au cours d’une étape suivante 206 réalisée par exécution du module logiciel 34, l’unité de traitement 30 fournit les M sensorgrammes représentatifs d’une interaction du fluide avec chaque capteur olfactif 58 sélectionné après avoir réalisé un calibrage aveugle et une estimation de déphasages. Au cours d’une étape optionnelle suivante 208 réalisée par exécution du module logiciel 66, l’unité de traitement 30 fournit les M sensorgrammes . Enfin, au cours d’une dernière étape 210 réalisée par exécution du module logiciel 68, l’unité de traitement 30 fournit la signature SIG ou SIG’ à partir des M sensorgrammes ou . Il apparaît clairement qu’un dispositif d’estimation de déphasage tel que l’un de ceux décrits précédemment en référence aux figures 1 et 2 permet d’améliorer le calibrage aveugle en prévoyant une étape de test sur des conditions requises pour réaliser l’estimation des paramètres de calibrage ainsi qu’une solution de repli vers les paramètres de calibrage d’une estimation précédente si le test échoue. Il en résulte une meilleure estimation du déphasage. Un tel dispositif trouve en outre une application particulièrement avantageuse dans la caractérisation d’un fluide par l’utilisation d’un dispositif électronique tel que celui illustré dans la . A titre d’illustration et afin de démontrer l’intérêt d’un tel calibrage amélioré dans le dispositif d’estimation de déphasage, une expérimentation a été menée sur un système triphasé résultant d’une interférométrie présentant une amplitude de déphasage d’environ 2 radians, ce qui est loin de couvrir la période 2π requise. Le système triphasé S tel que fourni est présenté en haut à gauche de la . Il présente 3 signaux p 1 , p 2 , p 3 dépendant du temps et présentant des variations d’amplitude A (exprimée sans unité, en niveaux de gris par exemple). En calibrant ce système triphasé S à l’aide de paramètres de calibrage d’une estimation précédente, par exemple en suivant les étapes 110, 112, 106B et 106C du procédé de la , il est possible d’obtenir le système triphasé calibré S’ dans lequel les signaux p’ 1 , p’ 2 , p’ 3 sont corrigés en amplitudes et décalages en amplitude. Bien que ce calibrage ne soit pas idéal, parce qu’il n’est pas réalisé sur le système triphasé S fourni lui-même, il reste avantageux de le prévoir en cas d’échec du test 104 parce que l’on remarque bien que le sensorgramme SG’ obtenu à partir de S’ est plus proche de la réalité de l’expérimentation que le sensorgramme SG obtenu à partir de S. En effet, S présente des variations de déphasage entre 0 et 0,7 radians, ce qui ne correspond pas à l’expérimentation conduite, alors que S’ évolue bien entre 0 et 2 radians. On notera par ailleurs que l’invention n’est pas limitée aux modes de réalisation décrits précédemment. Il apparaîtra en effet à l'homme de l'art que diverses modifications peuvent être apportées aux modes de réalisation décrits ci-dessus, à la lumière de l'enseignement qui vient de lui être divulgué. Dans la présentation détaillée de l’invention qui est faite précédemment, les termes utilisés ne doivent pas être interprétés comme limitant l’invention aux modes de réalisation exposés dans la présente description, mais doivent être interprétés pour y inclure tous les équivalents dont la prévision est à la portée de l'homme de l'art en appliquant ses connaissances générales à la mise en œuvre de l'enseignement qui vient de lui être divulgué. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage entre deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) issus d’un même signal périodique (s 0 ) de période T traité par interférométrie, comportant : un coupleur multimodal (24 ; 24 1 … 24 n ) à deux entrées, pour la réception des deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ), et à N sorties, N ≥ 3, pour la fourniture d’un système multiphasé (S) de N signaux formés pour être déphasés entre eux de T/N ; une unité (30) de traitement du système multiphasé (S) pour en extraire des paramètres de calibrage et une estimation (SG ; SG i ) du déphasage entre les deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) ; et des moyens (32) de stockage des paramètres de calibrage ; caractérisé en ce que l’unité de traitement (30) est configurée (34) pour tester si les N signaux du système multiphasé (S) remplissent des conditions prédéterminées requises pour réaliser une estimation des paramètres de calibrage et pour : réaliser l’estimation des paramètres de calibrage et leur enregistrement dans les moyens de stockage (32) si les conditions prédéterminées requises sont remplies ; et extraire, des moyens de stockage (32), les paramètres de calibrage d’une estimation précédente si les conditions prédéterminées requises ne sont pas remplies. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 1, dans lequel les conditions prédéterminées requises que l’unité de traitement (30) est configurée (34) pour tester comportent le fait que chaque signal du système multiphasé (S) fourni résulte de suffisamment de valeurs de déphasage différentes pour qu’elles couvrent au moins la période T. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 2, dans lequel le test réalisé par l’unité de traitement (30) comporte une vérification que chaque signal du système multiphasé (S) fourni varie suffisamment en amplitude pour en déduire ses valeurs extrémales pouvant être atteintes. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon l’une quelconque des revendications 1 à 3, dans lequel les conditions prédéterminées requises que l’unité de traitement (30) est configurée (34) pour tester comportent le fait qu’une transformation de Clarke, Concordia ou Fortescue du système multiphasé (S) fourni produise suffisamment de valeurs de déphasage différentes pour définir une unique solution d’ellipse par ajustement elliptique. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon l’une quelconque des revendications 1 à 4, dans lequel l’unité de traitement (30) est pourvue de moyens (34) de filtrage du système multiphasé (S) conçus pour estimer un état de ce sytème multiphasé (S) à partir d’une mesure, ces moyens de filtrage (34) étant plus précisément définis par le fait que : l’état à estimer comporte des valeurs constantes d’amplitudes et de décalages en amplitude des N signaux du système multiphasé (S), une valeur variable relative au déphasage à estimer et une valeur constante d’aide à l’estimation de la valeur variable ; et la mesure comporte les N signaux du système multiphasé (S) tels que fournis par le coupleur multimodal (24 ; 24 1 … 24 n ) ; et dans lequel les paramètres de calibrage comportent les valeurs constantes d’amplitudes et de décalages en amplitude des N signaux du système multiphasé (S). Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon l’une quelconque des revendications 1 à 4, dans lequel l’unité de traitement (30) est configurée pour : réaliser une transformation de Clarke, Concordia ou Fortescue du système multiphasé (S) et en déduire des paramètres d’une ellipse ajustée sur cette transformation de Clarke, Concordia ou Fortescue ; procéder à un calibrage de cette ellipse ajustée par transformation paramétrée en cercle recentré, dans lequel les paramètres de calibrage comportent les paramètres de la transformation ; puis estimer le déphasage entre les deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) après application des paramètres de calibrage à la transformation de Clarke du système multiphasé (S). Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon l’une quelconque des revendications 1 à 4, dans lequel l’unité de traitement (30) est configurée pour : procéder à un calibrage des N signaux périodiques du système multiphasé (S) à partir de valeurs d’amplitudes et de décalages en amplitude de ces N signaux périodiques, dans lequel les paramètres de calibrage comportent une combinaison de ces valeurs d’amplitudes et de décalages en amplitude ; réaliser une transformation de Clarke, Concordia ou Fortescue du système multiphasé (S) après calibrage et en déduire des paramètres d’un cercle centré ajusté sur cette transformation de Clarke, Concordia ou Fortescue ; puis estimer le déphasage entre les deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) à partir de cette transformation de Clarke, Concordia ou Fortescue. Dispositif électronique (50) de caractérisation d’un fluide, comportant : au moins un capteur (58) conçu pour interagir avec le fluide ; et un transducteur (60) conçu pour fournir, en interaction avec ledit au moins un capteur (58), une séquence d’au moins un signal de mesure (S) représentative d’une interaction du fluide avec chaque capteur (58) ; dans lequel le transducteur (60) est un système de traitement d’un signal périodique (s 0 ) par interférométrie et le dispositif électronique de caractérisation (50) comporte au moins un dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage entre deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) issus dudit signal périodique (s 0 ) selon l’une quelconque des revendications 1 à 7. Dispositif électronique (50) de caractérisation d’un fluide selon la revendication 8, dans lequel le transducteur (60) est un système d’amplification de variation d’indice optique par interférométrie de Mach-Zehnder à interférences multimodales, le coupleur multimodal (24 ; 24 1 … 24 n ) étant alors à sorties optiques et associé à un capteur optique (26). Procédé d’estimation d’un déphasage entre deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) issus d’un même signal périodique (s 0 ) de période T traité par interférométrie, comportant : un couplage multimodal (102) des deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) pour la fourniture d’un système multiphasé (S) de N signaux formés pour être déphasés entre eux de T/N où N ≥ 3 ; un traitement (104, 106, 110, 112) du système multiphasé (S) pour en extraire des paramètres de calibrage et une estimation (SG ; SG i ) du déphasage entre les deux signaux d’entrée (s 1 , s 2 ; s 1 ,1 , s 1,2 … s n,1 , s n,2 ) ; et un enregistrement (108) des paramètres de calibrage dans des moyens de stockage (32) ; caractérisé en ce que le traitement (104, 106, 110, 112) comporte un test (104) sur des conditions prédéterminées requises pour réaliser une estimation des paramètres de calibrage et est conçu pour exécuter sélectivement l’une ou l’autre des opérations suivantes : l’estimation (106 ; 106A ; 106’) des paramètres de calibrage et leur enregistrement (108) dans les moyens de stockage (32) si les conditions prédéterminées requises sont remplies par les N signaux du système multiphasé (S) ; et l’extraction (112), des moyens de stockage (32), de paramètres de calibrage d’une estimation précédente si les conditions prédéterminées requises ne sont pas remplies par les N signaux du système multiphasé (S).