FR 3149083 A1 20241129 FR 2305087 A 20230523 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 ), présentant N sorties pour la fourniture d’un système multiphasé (S) de signaux déphasés entre eux de T/N, et un processeur (30) de traitement du système multiphasé (S) pour une estimation (SG) du déphasage. Le processeur (30) est pourvu de moyens (34) de filtrage du système multiphasé (S) conçus pour estimer un état de ce système multiphasé à partir d’une mesure, sachant que : l’état à estimer comporte des valeurs constantes d’amplitudes et de décalages en amplitude des 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 ; la mesure comporte les signaux tels que fournis par le coupleur multimodal (24). 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 un processeur de traitement du système multiphasé pour en extraire 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. Il est connu d’appliquer une transformation de Clarke (ou équivalent) 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. 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, 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 précède nécessairement l’estimation du déphasage. Le traitement du système multiphasé pour en extraire une estimation du déphasage entre les deux signaux d’entrée procède donc en deux étapes successives. 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, 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. Mais cette deuxième solution présente le même inconvénient que la précédente, à savoir un calibrage qui précède nécessairement l’estimation du déphasage et un traitement du système multiphasé en deux étapes successives. Elle a cependant le mérite de montrer qu’il est possible de procéder au calibrage 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 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 ; un processeur de traitement du système multiphasé pour en extraire une estimation du déphasage entre les deux signaux d’entrée ; dans lequel le processeur de traitement est pourvu de moyens de filtrage du système multiphasé pour estimer un état de ce systè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. Ainsi, l’estimation du déphasage se fait par un filtrage adaptatif, comportant ce déphasage en tant que paramètre d’état variable à estimer et faire évoluer au cours du temps, et en même temps que le calibrage puisque les valeurs d’amplitudes et de décalages en amplitude des N signaux du système multiphasé sont également intégrées dans la configuration des moyens de filtrage en tant que paramètres d’état constants à estimer. Il est en outre astucieux de prévoir un paramètre d’état constant d’aide à l’estimation de la valeur variable pour faciliter le suivi des évolutions au cours du temps de cette dernière. Il en résulte un moyen simple et efficace pour fusionner le calibrage et l’estimation de déphasage. On observe par ailleurs qu’il n'est généralement pas nécessaire non plus d’opérer un déroulement de phase à l’issue d’un tel filtrage puisque la continuité du déphasage est implicitement intégrée dans le modèle d’état des moyens de filtrage. De façon optionnelle, les moyens de filtrage comportent un filtre récursif non linéaire de type Kalman étendu ou sans parfum. De façon optionnelle également, lequel les moyens de filtrage comportent un filtre du système multiphasé à deux fonctions de transfert : une fonction de transition 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 calculant une estimation de la mesure à partir de la valeur prédite ou déjà estimée de l’état du système multiphasé ; le filtre étant en outre défini de sorte que : l’état du système multiphasé se présente sous la forme d’un vecteur dont les coefficients comportent les valeurs constantes d’amplitudes et de décalages en amplitude des N signaux du système multiphasé, la valeur variable relative au déphasage à estimer et la valeur constante d’aide à l’estimation de la valeur variable ; et la mesure se présente sous la forme d’un vecteur à N coefficients dont les N coefficients comportent respectivement les N signaux du système multiphasé tels que fournis par le coupleur multimodal. De façon optionnelle également, la fonction de transition est définie de telle sorte que la valeur prédite du vecteur d’état est égale à la somme : de la valeur déjà estimée 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 d’un vecteur de bruit de traitement. De façon optionnelle également, les moyens de filtrage comportent un filtre de Kalman sans parfum et la fonction d’observation est définie de telle sorte que l’estimation de la mesure est égale à la somme : d’un vecteur à N coefficients dont les N coefficients sont respectivement les N signaux du système multiphasé exprimés en fonction des coefficients de la valeur déjà estimée du vecteur d’état ; et d’un vecteur de bruit de mesure. De façon optionnelle également et en variante, les moyens de filtrage comportent un filtre de Kalman étendu et la fonction d’observation est définie de telle sorte que le vecteur de mesure est égal à la somme : du produit matriciel d’une matrice jacobienne des N signaux du système multiphasé exprimés en fonction des coefficients de la valeur déjà estimée du vecteur d’état avec la valeur déjà estimée du vecteur d’état ; et d’un vecteur de bruit de mesure. De façon optionnelle également, le signal périodique traité par interférométrie est un signal sinusoïdal de période 2π et le coupleur multimodal est à N=3 sorties pour la fourniture d’un système triphasé de trois signaux sinusoïdaux p 1 , p 2 et p 3 , c’est-à-dire formés pour être déphasés entre eux de 2π /3, exprimables sous la forme suivante : ; ; ; où a 1 , a 2 et a 3 sont les coefficients du vecteur d’état relatifs aux valeurs constantes d’amplitudes des trois signaux sinusoïdaux du système triphasé, o 1 , o 2 et o 3 sont les coefficients du vecteur d’état relatifs aux valeurs constantes de décalages en amplitude des trois signaux sinusoïdaux du système triphasé, et θ est le coefficient du vecteur d’état comportant la valeur variable relative au déphasage à estimer. 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. En fin, 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 une estimation du déphasage entre les deux signaux d’entrée ; dans lequel le traitement comporte un filtrage du système multiphasé pour estimer un état de ce système 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 N signaux du 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 N signaux du système multiphasé tels que fournis par le coupleur multimodal. 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 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 , et 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 . 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. 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 peut être définie angulairement à 2π, 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 pourrait par exemple procéder à un traitement en deux temps, i.e. calibrage puis estimation du déphasage par transformation de Clarke recalée, tel qu’enseigné dans le document Halir et al précité ou dans le document WO 2022/238170 A1 également précité. Mais conformément aux principes généraux de la présente invention et de façon avantageuse, il comporte des moyens de filtrage du système triphasé S pour estimer de façon récursive un état de ce système triphasé S à partir d’une mesure. Il s’agit donc d’un filtrage de type Kalman. Ces moyens de filtrage sont définis et configurés 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. 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 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ù la calibration est réalisée 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’une précédente calibration/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 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 d’un système triphasé et échantillonné temporellement, le vecteur d’état à un instant d’échantillonnage k peut donc 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 multiphasé, 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 multiphasé : , 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 : . 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é 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 : , 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 multiphasé 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é, 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. 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 filtrage précités sous la forme du module logiciel ou programme d’ordinateur 34 apte à s’exécuter cette fois-ci sur les n systèmes triphasés fournis. Le fonctionnement du système interférométrique 10, ou celui équivalent du système interférométrique 40, va maintenant être détaillé 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 électriques 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). Enfin, au cours d’une deuxième et dernière étape 104 de procédé d’estimation de déphasage, le système triphasé ou ensemble de 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. 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 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 signaux de mesure fournie par le transducteur 60. La figure 5 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 4, 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é. 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 de combiner les fonctions de calibrage et d’estimation de déphasage en les réalisant conjointement et simultanément par un filtrage spécifique et astucieux. 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 . 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 ; un processeur (30) de traitement du système multiphasé (S) pour en extraire 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 ) ; caractérisé en ce que le processeur de traitement (30) est pourvu de moyens (34) de filtrage du système multiphasé (S) conçus pour estimer un état de ce systè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 ). Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 1, dans lequel les moyens de filtrage (34) comportent un filtre récursif non linéaire de type Kalman étendu ou sans parfum. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 1 ou 2, dans lequel les moyens de filtrage (34) comportent un filtre du système multiphasé (S) à deux fonctions de transfert : une fonction de transition calculant une valeur prédite à partir d’une valeur déjà estimée de l’état du système multiphasé (S) ; et une fonction d’observation calculant une estimation de la mesure à partir de la valeur prédite ou déjà estimée de l’état du système multiphasé (S) ; le filtre étant en outre défini de sorte que : l’état du système multiphasé (S) se présente sous la forme d’un vecteur dont les coefficients comportent les valeurs constantes d’amplitudes et de décalages en amplitude des N signaux du système multiphasé (S), la valeur variable relative au déphasage à estimer et la valeur constante d’aide à l’estimation de la valeur variable ; et la mesure se présente sous la forme d’un vecteur à N coefficients dont les N coefficients comportent respectivement les N signaux du système multiphasé (S) tels que fournis par le coupleur multimodal (24 ; 24 1 … 24 n ). Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 3, dans lequel la fonction de transition est définie de telle sorte que la valeur prédite du vecteur d’état est égale à la somme : de la valeur déjà estimée 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 d’un vecteur de bruit de traitement. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 3 ou 4, dans lequel les moyens de filtrage (34) comportent un filtre de Kalman sans parfum et la fonction d’observation est définie de telle sorte que l’estimation de la mesure est égale à la somme : d’un vecteur à N coefficients dont les N coefficients sont respectivement les N signaux du système multiphasé (S) exprimés en fonction des coefficients de la valeur déjà estimée du vecteur d’état ; et d’un vecteur de bruit de mesure. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon la revendication 3 ou 4, dans lequel les moyens de filtrage comportent un filtre de Kalman étendu et la fonction d’observation est définie de telle sorte que le vecteur de mesure est égal à la somme : du produit matriciel d’une matrice jacobienne des N signaux du système multiphasé (S) exprimés en fonction des coefficients de la valeur déjà estimée du vecteur d’état avec la valeur déjà estimée du vecteur d’état ; et d’un vecteur de bruit de mesure. Dispositif (22 ; 22 1 … 22 n ) d’estimation d’un déphasage selon l’une quelconque des revendications 3 à 6, dans lequel le signal périodique (s 0 ) traité par interférométrie est un signal sinusoïdal de période 2π et le coupleur multimodal (24 ; 24 1 … 24 n ) est à N=3 sorties pour la fourniture d’un système triphasé de trois signaux sinusoïdaux p 1 , p 2 et p 3 , c’est-à-dire formés pour être déphasés entre eux de 2π /3, exprimables sous la forme suivante : ; ; ; où a 1 , a 2 et a 3 sont les coefficients du vecteur d’état relatifs aux valeurs constantes d’amplitudes des trois signaux sinusoïdaux du système triphasé, o 1 , o 2 et o 3 sont les coefficients du vecteur d’état relatifs aux valeurs constantes de décalages en amplitude des trois signaux sinusoïdaux du système triphasé, et θ est le coefficient du vecteur d’état comportant la valeur variable relative au déphasage à estimer. 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 7, 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 du système multiphasé (S) pour en extraire 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 ) ; caractérisé en ce que le traitement comporte un filtrage (104) du système multiphasé (S) pour estimer un état de ce système multiphasé (S) à partir d’une mesure, ce filtrage (104) é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 N signaux du système multiphasé (S), une 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.