FR 2961897 A1 20111230 FR 1002675 A 20100625 Filtre de navigation pour système de navigation par corrélation de terrain La présente invention concerne un filtre de navigation. Elle s'applique notamment au domaine de la navigation de porteurs, et plus particulièrement à des systèmes de navigation opérant notamment suivant des procédés de corrélation de terrain, la navigation par corrélation de terrain étant désignée par le sigle NCT. La navigation par corrélation de terrain ou NCT, constitue un moyen particulier de navigation, pouvant s'appliquer à une grande diversité de véhicules porteurs, par exemple des aéronefs, des sous-marins, des missiles 1 o autoportés, etc. Il existe trois moyens principaux connus visant à répondre aux besoins de navigation de porteurs. Le premier principal moyen connu comprend les techniques de navigation inertielle. Le deuxième principal moyen connu comprend les techniques de radionavigation. Le troisième 15 principal moyen connu comprend les techniques de navigation par corrélation de terrain. La navigation inertielle consiste à exploiter des informations fournies par des centrales inertielles. Le fonctionnement d'une centrale inertielle se fonde sur le principe de relativité d'Einstein-Galilée, postulant qu'il est 20 possible, sans l'aide de signaux extérieurs à un porteur, de mesurer d'une part la vitesse de rotation du porteur par rapport à un référentiel inertiel, par exemple défini par un repère géocentrique associé à des étoiles fixes, et d'autre part la force spécifique appliquée au porteur : typiquement son accélération dans le référentiel inertiel, diminuée de l'accélération de la 25 pesanteur. D'une manière typique, un système de navigation inertielle, communément désigné INS selon le sigle correspondant à la terminologie anglaise "Inertial Navigation System", est un dispositif permettant de mesurer ces deux grandeurs au moyen de capteurs tels que des gyromètres et des accéléromètres, communément au nombre de trois de chaque type, disposés 30 selon trois axes orthogonaux, l'ensemble de ces capteurs formant une unité de mesure inertielle, communément désignée "IMU" selon le sigle correspondant à la terminologie anglaise "Inertial Measurement Unit". L'intégration temporelle des données d'accélération, et la projection dans le repère de navigation à partir des données de vitesse de rotation, permettent la détermination de la position et de la vitesse du porteur par rapport à la terre, en connaissance d'un état initial de ces données. Cependant, un inconvénient lié à l'intégration temporelle est que l'erreur associée aux données ainsi déterminées, est une fonction croissante du temps. Cette erreur croît plus que linéairement, typiquement de façon exponentielle, la variation de l'erreur étant communément désignée dérive de la centrale inertielle. II est ainsi nécessaire, pour des applications nécessitant une navigation de précision, d'hybrider les mesures inertielles avec d'autres mesures de position et/ou de vitesse et/ou d'attitude du porteur fournies par des capteurs complémentaires, tels que des baroaltimètres, des odomètres, des sondes de Pitot, etc., dans le but de diminuer la dérive de la centrale inertielle. De tels capteurs fournissent des informations sur l'état cinématique du porteur sans nécessiter aucun accès à des signaux extérieurs ou des cartes embarquées, et sont communément désignés capteurs de bas niveau. La radionavigation consiste à exploiter les signaux issus de balises émettrices de signaux radioélectriques, pour en tirer une information de positionnement du porteur vis-à-vis de ces balises. Une technique de radionavigation largement utilisée est la technique de géo-positionnement par satellite, communément désignée selon le sigle GNSS correspondant à la terminologie anglaise "Global Navigation Satellite System", et dont un des représentants est la technique de GPS, correspondant à la terminologie anglaise "Global Positioning System". Un des inconvénients propres aux techniques de radionavigation, est lié au fait que la réception des signaux en provenance des balises n'est pas garantie en tout lieu à tout moment, et peut notamment être affectée par l'environnement géophysique du porteur, ainsi que par le bruit électromagnétique environnant, des techniques de brouillage pouvant notamment mettre en cause le fonctionnement d'un dispositif de radionavigation. En outre, les balises émettrices étant maintenues par des opérateurs, l'intégrité des données de radionavigation qui en sont issues est fortement dépendante de la bonne volonté de ces derniers. La radionavigation, et notamment le géo-positionnement par satellite et la navigation inertielle sont par exemple des techniques de navigation complémentaires, et une hybridation des deux techniques peut s'avérer très performante en pratique. La navigation inertielle constitue en effet un très bon estimateur local de positionnement dérivant à long terme, et le géopositionnement par satellite étant peu fiable sur une courte durée en raison des inconvénients précités, mais ne présentant pas de dérive. Cependant, dans les applications les plus critiques, et notamment pour les applications militaires, il est essentiel de recourir à d'autres sources d'information de position et/ou de vitesse et/ou d'attitude du porteur pour assurer une hybridation avec une technique de navigation inertielle. II est notamment souhaitable que ces sources alternatives permettent des mesures de position et/ou de vitesse et/ou d'attitude du porteur qui soient autonomes, non sujettes au brouillage, et discrètes. La navigation par corrélation de terrain ou NCT consiste à exploiter des mesures de données géophysiques délivrées par un capteur adéquat, à des données de référence propres à un terrain de couverture de la navigation. Les capteurs sont ainsi utilisés conjointement à une carte de référence du terrain, également désignée carte embarquée. Ces capteurs permettent la lecture d'une donnée caractéristique sur le terrain, et la navigation par corrélation de terrain consiste à comparer ces valeurs avec les données de la carte embarquée, la carte embarquée étant un relevé a priori des valeurs de ces données sur la zone de navigation considérée, obtenu par des moyens adéquats, et désigné pour la suite filière de production des données. La navigation par corrélation de terrain est particulièrement adaptée pour une hybridation avec une technique de navigation inertielle, et permet de pallier les insuffisances de la radionavigation. Bien sûr, il est possible, pour des performances optimales, de recourir à un système de navigation permettant une hybridation des trois techniques de navigation précitées. D'une manière générale, tout système de navigation impliquant une corrélation de terrain comprend ainsi une pluralité de capteurs embarqués compris dans la centrale inertielle, ainsi que le capteur de terrain, une carte embarquée représentant la meilleure connaissance possible sur la réalité des données géophysiques que le capteur embarqué doit mesurer, et un filtre de navigation. Le filtre de navigation permet d'arbitrer en temps réel, entre les informations fournies par la centrale inertielle et celles fournies par la comparaison entre les mesures fournies par le capteur de terrain et la carte embarquée. L'arbitrage réalisé par le filtre en fonction de sa connaissance a priori des erreurs sur les mesures fournies. Cette connaissance est regroupée dans des modèles d'erreur. Les modèles d'erreur sont relatifs à la centrale inertielle, les erreurs de la centrale inertielle étant plus ou moins grandes selon la qualité de l'équipement ; les modèles d'erreur concernent également le capteur de terrain, ainsi que la carte embarquée, les erreurs de cette dernière étant plus ou moins grande selon la qualité de la filière de production des données. Les modèles d'erreur des équipements proviennent des informations fournies par les fabricants, et / ou bien de mesures effectuées via des études spécifiques. Les modèles d'erreur des cartes embarquées sont fournis par les producteurs de données. Un aspect essentiel de la navigation est le caractère stochastique des phénomènes considérés. En effet, les capteurs produisent des erreurs suivant des modèles stochastiques, et la connaissance des données géophysiques étant peu maîtrisée, la résolution du problème de la navigation par un filtrage rend les performances de navigation intrinsèquement stochastiques. Ainsi, le filtre utilisé dans un système de navigation peut être considéré comme un estimateur d'un processus stochastique, c'est-à-dire comme le dispositif qui donne à tout instant l'état dynamique du porteur modélisé comme une variable aléatoire. Un premier exemple de système de navigation impliquant une corrélation de terrain est fondé sur la technique de navigation altimétrique. Cette technique consiste à faire naviguer un porteur aérien à l'aide d'une centrale inertielle, d'un capteur de terrain de type radioaltimètre ou scanneur laser multi-faisceau, mesurant la distance du porteur au terrain dans une ou plusieurs direction(s) donnée(s), et d'une carte embarquée de type Modèle Numérique de Terrain ou MNT, recensant les altitudes de points du sol sur une grille régulière géolocalisée. Un deuxième exemple de système de navigation impliquant une corrélation de terrain est fondé sur la technique de navigation bathymétrique. Cette technique consiste à faire naviguer un porteur marin ou sous-marin à l'aide d'une centrale inertielle, d'un capteur de terrain de type sondeur bathymétrique mono-faisceau ou multi-faisceau, mesurant la distance du porteur au fond de la mer dans une ou plusieurs direction(s) donnée(s), et d'une carte embarquée de type carte bathymétrique, recensant les altitudes de points du fond marin sur une grille régulière géolocalisée. Un troisième exemple de système de navigation impliquant une corrélation de terrain est fondé sur la technique de navigation gravimétrique. Cette technique consiste à faire naviguer un porteur aérien, marin ou sous-marin à l'aide d'une centrale inertielle, d'un capteur de terrain de type gravimètre ou accéléromètre, mesurant le champ de pesanteur local ou son anomalie, et d'une carte embarquée de type carte d'anomalie gravimétrique, recensant les valeurs de l'anomalie du champ de pesanteur en des points du globe sur une grille régulière normalisée. Un quatrième exemple de système de navigation impliquant une corrélation de terrain est fondé sur la technique de navigation par vision. Cette technique consiste à faire naviguer un porteur aérien à l'aide d'une centrale inertielle, d'un capteur de terrain de type caméra embarquée qui délivre des images du paysage survolé à une fréquence donnée dans le domaine visible ou infrarouge, et de deux cartes embarquées, une carte embarquée de type Orthoimage géolocalisée, c'est-à-dire une image ré-échantillonnée de telle sorte qu'ont été supprimés les effets du relief, c'est-à-dire pour laquelle l'échelle est la même en tous les points, ainsi qu'une carte embarquée de type MNT. Dans le cadre de systèmes de navigation impliquant une corrélation de terrain, les concepteurs sont notamment confrontés à un certain nombre de problèmes techniques énoncés ci-dessous : il est nécessaire de définir un système de navigation permettant d'atteindre une qualité souhaitée de navigation selon un ensemble de critères de performance déterminés, par exemple garantissant une erreur moyenne de positionnement inférieure à un seuil donné, ceci à moindre coût ; il est nécessaire de déterminer des modèles d'erreur les plus fidèles possible pour la centrale inertielle, le capteur de terrain et la carte embarquée ; il est nécessaire de définir les missions d'un porteur, notamment en terme de trajectoire de consigne, lors d'une phase de préparation de mission, afin de déterminer une 30 35 trajectoire optimale le long de laquelle la qualité du signal délivré par le capteur de terrain soit maximale, la trajectoire optimale devant également être définie en fonction d'autres critères de performance de la mission du porteur et de contraintes opérationnelles liées au théâtre de la mission. La phase de préparation de mission doit par exemple se fonder sur un critère de navigabilité qui soit pertinent, c'est-à-dire représentatif de la richesse du signal délivré par le capteur de terrain ; il est nécessaire de définir un filtre de navigation performant, robuste et capable de prendre en considération au mieux tous les modèles d'erreur relatifs aux différents composants du système, c'est-à-dire l'erreur de la centrale inertielle, du capteur de terrain et de la carte embarquée. 15 La présente invention a principalement pour objet de résoudre le problème technique précité, relatif à la définition d'un filtre de navigation performant. Selon des techniques connues de l'état de la technique, les filtres de navigation employés dans des systèmes de NCT sont des filtres de 20 navigation de type filtre de Kalman étendu, communément désigné par le sigle EKF correspondant à la terminologie anglaise "Extended Kalman Filter". Ces filtres sont notoirement peu robustes en cas d'une faible information remontant du terrain, présentant des cas pouvant entraîner une divergence du filtre. Des solutions connues permettant d'éviter ces divergences, 25 consistent à coupler des filtres EKF avec des algorithmes de recalage par bloc. Cependant, d'une manière typique, des phases de recalage par blocs peuvent durer de l'ordre d'une vingtaine de secondes à une minute, laps de temps pendant lequel aucune information sur le comportement du système n'est remontée. De tels effets de "silence" peuvent avoir un impact fâcheux 30 sur la navigation d'un porteur, qui peut parcourir ainsi jusqu'à 15 km, pour une vitesse de déplacement égale à 250 m/s, sans aucun moyen d'apprécier la qualité de sa navigation. Des filtres non linéaires plus génériques que les filtres de Kalman sont connus de l'état de la technique. Il s'agit des filtres particulaires, qui pallient 35 les défauts des filtres EKF. Néanmoins, les filtres particulaires présentent 10 des effets de dégénérescence au bout d'un temps relativement long de navigation. De ce fait, les filtres particulaires n'ont jusqu'ici jamais été exploités en pratique dans des applications de navigation par corrélation de terrain. Un but de la présente invention est de pallier au moins les inconvénients précités, en proposant un filtre de navigation offrant de bonnes caractéristiques de robustesse, de bonnes performances de convergence, et ne présentant pas l'inconvénient liés aux effets de silence précités. 10 A cet effet, la présente invention a pour objet un filtre de navigation pour système de navigation par corrélation de terrain délivrant une estimation de l'état cinématique d'un porteur associé à une matrice de covariance à un instant discrétisé k à partir d'une pluralité de données comprenant les 15 mesures restituées par au moins un capteur de terrain, le modèle associé au capteur de terrain, les données d'une carte embarquée, un modèle d'erreur de la carte embarquée, les mesures restituées par une centrale inertielle et un modèle de la centrale inertielle, le filtre de navigation étant caractérisé en ce qu'il comprend un premier filtre dit filtre de convergence et un second filtre 20 dit filtre de poursuite, le filtre de navigation comprenant des moyens de basculement sélectionnant pour le calcul de l'état cinématique du porteur et de la matrice de covariance, l'un ou l'autre desdits filtres de convergence et filtre de poursuite, en fonction de la comparaison d'un indice de qualité calculé à partir de la matrice de covariance restituée par le filtre de 25 navigation à l'instant précédent k-1, à une valeur seuil prédéterminée, le filtre de navigation restituant les valeurs calculées par le filtre sélectionné. Dans un mode de réalisation de l'invention, l'indice de qualité peut être égal à une norme des termes d'une sous-matrice de la matrice de covariance. 30 Dans un mode de réalisation de l'invention, l'indice de qualité peut être égal à une norme des termes de la sous-matrice de covariance comprenant les termes de position du porteur dans un plan horizontal. Dans un mode de réalisation de l'invention, l'indice de qualité peut être égal à la somme des carrés des termes diagonaux de la sous-matrice de5 covariance comprenant les termes de position du porteur dans un plan horizontal. Dans un mode de réalisation de l'invention, le filtre de convergence peut être un filtre de type particulaire. Dans un mode de réalisation de l'invention, le filtre de convergence peut être un filtre de type particulaire marginalisé. Dans un mode de réalisation de l'invention, le filtre de poursuite peut être un filtre de type filtre de Kalman étendu désigné suivant l'acronyme EKF. Dans un mode de réalisation de l'invention, le filtre de poursuite peut être un filtre du type désigné suivant l'acronyme UKF. Dans un mode de réalisation de l'invention, le filtre de poursuite peut être associé à un dispositif de rejet de mesures aberrantes, rejetant les mesures si le rapport entre l'innovation et l'écart type du bruit de mesure dépasse une deuxième valeur de seuil prédéterminée. Dans un mode de réalisation de l'invention, les moyens de basculement restituent en sortie du filtre de navigation les données restituées par ledit filtre de convergence lorsque le nombre de mesures consécutives rejetées par le dispositif de rejet est supérieur à une troisième valeur de seuil prédéterminée. Dans un mode de réalisation de l'invention, le filtre de navigation peut en outre comprendre des moyens de comparaison de l'écart type local des mesures de navigation par corrélation de terrain à une quatrième valeur de seuil prédéterminée, le filtre de navigation se basant uniquement sur les mesures fournies par le baro-altimètre compris dans la centrale inertielle lorsque l'écart type est inférieur à ladite quatrième valeur de seuil prédéterminée. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description, donnée à titre d'exemple, faite en regard des 30 dessins annexés qui représentent : la figure 1, un schéma illustrant de manière synoptique la structure d'un système de navigation impliquant une corrélation de terrain ; la figure 2, un schéma illustrant de manière synoptique la structure d'un filtre de navigation selon un exemple de réalisation de l'invention ; la figure 3, un logigramme illustrant le fonctionnement d'un filtre de navigation selon un exemple de réalisation de l'invention. La figure 1 présente un schéma illustrant de manière synoptique la structure d'un système de navigation impliquant une corrélation de terrain. Un système de navigation 1 comprend un bloc de navigation 10, un filtre de navigation 11 et une carte embarquée 12. Le bloc de navigation 10 peut comprendre une centrale inertielle 101 et un ou une pluralité de capteurs de terrain 102. La centrale inertielle 101 peut notamment comprendre un calculateur inertiel 1010, recevant des données en provenance d'un modèle de pesanteur 1011, d'une IMU 1012 et d'un ou une pluralité de capteurs de bas niveau 1013. Le filtre de navigation 11 reçoit des données de position, de vitesse et d'attitude du porteur, en provenance du calculateur inertiel 1010. Le filtre de navigation 11 reçoit des mesures de données géophysiques en provenance des capteurs de terrain 102. Egalement, le filtre de navigation 11 accède aux données contenues dans la carte embarquée 12. Le filtre de navigation 11 peut être mis en oeuvre dans des dispositifs de calcul adéquats, et restitue des estimations de l'état cinématique du porteur. Le filtre de navigation 11 est également apte à appliquer des corrections aux paramètres de configuration du calculateur inertiel 1010 et des capteurs de terrain 102 et de la carte embarquée 12. D'une manière typique, le filtre de navigation 11 peut par exemple corriger des biais des capteurs de terrain 102, ou encore des dérives de la centrale inertielle 101 ou bien des paramètres du modèle d'erreur de la carte embarquée 12. La carte embarquée 12 peut par exemple être formée par un ensemble de différentes cartes de différentes natures correspondant à chacun des capteurs de terrain impliqués, données stockées dans une mémoire. La solution proposée par la présente invention est basée sur la collaboration d'un premier filtre appelé par la suite « filtre de convergence », par exemple de type particulaire, offrant de meilleures performances dans les phases de convergence ou de divergence, associé à un second filtre appelé par la suite « filtre de poursuite », par exemple issu de la famille des filtres de Kalman offrant de bonnes performances dans les phases de poursuite d'un porteur, c'est-à-dire dans les phases où l'erreur de positionnement est assez faible,. Ainsi que cela est décrit précédemment, il est possible, selon des méthodes en elles-mêmes connues de l'état de la technique, de recourir à un algorithme de corrélation par blocs lorsque le modèle de mesure « verticale » est licite. Selon un tel algorithme, un profil de mesure est stocké puis comparé avec le terrain dans la zone d'incertitude jusqu'à obtenir une corrélation satisfaisante. Lorsque l'incertitude sur les paramètres est suffisamment faible, il est alors possible de recourir à un filtre de type EKF, qui prend le relais de l'algorithme de corrélation par blocs. Cependant de telles méthodes présentent l'inconvénient de générer des effets de silence dans certaines situations. Un filtre de navigation selon la présente invention peut être illustré par le schéma en synoptique présenté par la figure 2. En référence à la figure 2, un filtre de navigation 11 fournissant en sortie une estimation de l'état cinématique d'un porteur, par exemple une estimation de la position, de la vitesse et de l'attitude de ce dernier, reçoit en entrée une pluralité de données, comprenant les mesures réalisées par le capteur de terrain, les mesures restituées par la centrale inertielle, le modèle de capteur de terrain, les données de la carte embarquée, le modèle d'erreur de carte, par exemple de type Gauss-Markov, le modèle de la centrale inertielle. Le filtre de navigation 11 peut comprendre un premier filtre désigné filtre de convergence 21, par exemple de type filtre particulaire, et un second filtre désigné filtre de poursuite 22, par exemple de type filtre de Kalman. Plus précisément, le filtre de poursuite 22 peut être du type désigné EKF ou UKF respectivement selon les terminologies anglaises « Extended Kalman Filter » et "Unscented Kalman Filter".35 La navigation optimale peut être considérée comme un problème de filtrage non linéaire. Le système dynamique de filtrage, discrétisé dans le domaine temporel, peut être formalisé par la relation suivante : Xk =fk(Xk-1)+Wk (1),où: Y k =h k(X k)+V k Xk désigne l'état, d'une dimension d, du système dynamique à une étape temporelle k ; fk(Xk_1) désigne le modèle de transition d'état appliqué à l'état dynamique du système à l'étape temporelle précédente k-1 ; Wk désigne le bruit de process, considéré comme issu d'une distribution 1 o normale de moyenne nulle et de covariance Qk : Wk -:N(0, Qk) ; Mathématiquement, la mesure de la position de la plate-forme par un capteur géophysique se représente par la seconde équation du système d'équations formulé dans la relation (1) ci-dessus, dans laquelle h est une fonction a priori non linéaire de la position et Wk est une erreur de mesure 15 généralisée dépendant à la fois : - de l'erreur de mesure du capteur proprement dit, - de l'incertitude associée à la cartographie de la donnée géophysique ; Yk désigne l'observation de l'état dynamique du système à l'étape 20 temporelle k ; hk désigne le modèle d'observation ; Vk désigne le bruit d'observation, issu d'une distribution de type bruit blanc Gaussien, de moyenne nulle et de covariance Rk : Vk N(0, Rk) . Le bruit de process Wk et le bruit d'observation Vk sont supposés indépendants. L'hybridation de cette mesure et de l'équation d'évolution, la première 25 équation du système d'équations posé par la relation (1) ci-dessus, revient mathématiquement à un problème de filtrage non linéaire. Le problème posé est alors d'estimer à tout instant t la distribution de probabilité conditionnelle de l'état connaissant l'ensemble des mesures. Xk désignant l'état de la plate-forme, le filtrage optimal du processus 30 aléatoire rassemblant l'ensemble des variables aléatoires, {Xk}, à partir de la suite des variables aléatoires, {yk}, revient à calculer la probabilité conditionnelle exprimée suivant la relation suivante : p(X(k)ly(t),t On peut néanmoins classer les filtres dans les deux catégories de filtres de poursuite, adaptés aux phases où les incertitudes sur la position sont faibles, ou filtres de convergence, adaptés aux phases de grande incertitude. Pour chacune de ces deux catégories de filtres, il est possible de sélectionner : - les filtres de la famille Kalman pour la phase de poursuite, - les filtres de type particulaire pour la phase de convergence. Quelques approximations du filtrage non linéaire sur la base de filtres de la famille Kalman et de la famille des filtres particulaires sont décrites ci- après. Pour la suite, il est considéré le filtrage d'un processus X satisfaisant au modèle d'évolution d'état décrit par la relation suivante : Xk+1 - J (k,Xk )+Vk (7), où vk est un bruit blanc. Un processus d'observation Yk est relié à Xk par la relation suivante : Yk =h(k,Xk)+Wk (8), où Wk est un bruit blanc, pouvant être modélisé comme indépendant de vk comme dans les hypothèses du système formulé par la relation (1) précédemment. Lorsque les fonctions x->f(k,x) et x--+h(k,x) sont linéaires en x, et lorsque les bruits d'état et d'observation sont des bruits gaussiens, la distribution de probabilité conditionnelle lak (dx) est gaussienne, et ses deux premiers moments sont donnés par des équations récursives qui sont les équations du filtre de Kalman. L'avantage majeur offert par le filtre de Kalman est que la performance du filtre peut être calculée hors ligne en résolvant une équation de type Ricatti. Dans le cas non linéaire général, il est nécessaire d'approximer les équations permettant de gérer les étapes de prédiction et de correction du filtre. Il est à observer que, quelle que soit l'approximation choisie, la performance de la méthode de filtrage ne pourra d'une manière générale être appréciée que par simulation Monte Carlo. L'approche la plus élémentaire, basée sur des filtres de type EKF, repose sur la linéarisation récursive, par approximation de Taylor, des fonctions x-->f(k,x) et x-*h(k,x) autour de l'estimée de l'état, qui peut alors être propagé par un filtre de Kalman, comme dans le cas linéaire. II est possible de remplacer l'approximation de Taylor par une approximation stochastique, ou approximation d'Hermite, qui conduit au filtre de type UKF. Dans un tel filtre, des points, dénommés sigma-points sont utilisés, et les sigma-points sont associés à des poids, issus de formules de quadrature numérique, pour représenter les distributions de probabilité conditionnelle en évitant toute linéarisation. Il est également possible d'approximer la loi conditionnelle pk(dx) par un mélange de gaussiennes, suivant la relation suivante : n=N ,uk (dx) = E pk Gauss; (dx) (9), où n=1 Gauss désigne une loi Normale de moyenne Mk et de covariance Sk . Cette approximation conduit à une combinaison de N filtres de Kalman. Les avantages procurés par ces techniques de la famille Kalman sont liés à leur simplicité de mise en oeuvre, au faible coût en terme de temps de calcul, et aux bonnes performances dans les phases de poursuite. Les inconvénients sont liés à la mauvaise prise en compte de la multi-modalité ou d'autres formes d'ambiguïté, même temporaires, liées souvent à de grandes incertitudes, et à la difficulté intrinsèque à linéariser une fonction d'observation définie seulement ponctuellement, par lecture sur une carte numérique. Le filtrage particulaire, décrit ci-après, vise à apporter une solution à certains de ces inconvénients en approximant les calculs bayésiens d'intégrales des étapes de prédiction et de correction par des méthodes de Monte Carlo. Il est à observer que l'approximation des intégrales bayésiennes peut conduire à différentes mises en oeuvre. L'approche de type filtrage particulaire consiste à approcher numériquement le filtre bayésien, c'est-à-dire la distribution de probabilité conditionnelle pk(dx) à l'aide de particules, qui sont autant d'hypothèses pour représenter l'état caché X, éventuellement affectées de poids exprimant la pertinence de chaque particule pour représenter l'état caché. Une approximation peut alors se formuler suivant la relation suivante : ll N P[Xk E Cl~.xlYo:kJ wksei (10), où i=1 les particules k forment une grille aléatoire qui suit les évolutions de la suite des états cachés, et se concentre sur les zones d'intérêt, au niveau 15 desquelles la vraisemblance est forte. Dans l'étape de prédiction, on utilise donc le modèle d'évolution des erreurs d'estimation inertielle, des biais inertiels et des éventuels hyperparamètres, pour permettre aux particules d'explorer l'espace d'état et de former des hypothèses pertinentes. 20 Dans l'étape de correction, on évalue la vraisemblance de chaque particule, qui permet de quantifier la pertinence de chacune de ces hypothèses réalisées pour l'état caché avec les observations mesurées par le capteur de terrain et avec les estimations de position délivrées par la centrale, et on met ainsi à jour le poids de chaque particule. A chaque 25 itération, ou bien seulement lorsqu'un critère permettant de quantifier le déséquilibre dans la répartition des poids des particules dépasse un seuil déterminé, les particules sont redistribuées en fonction de leurs poids respectifs, avec pour effet de multiplier les particules les plus prometteuses, ou au contraire d'éliminer les particules les moins prometteuses. 30 Ce mécanisme de mutation / sélection a pour effet de concentrer automatiquement les particules, c'est-à-dire la puissance de calcul disponible, dans les régions d'intérêt de l'espace d'état. Une variante communément utilisée pour les filtres particulaires consiste à exploiter la nature conditionnellement linéaire gaussienne du problème pour propager un système de particules dans un espace d'état de dimension réduit : cette méthode est dite de Rao-Blackwell ou de marginalisation. II est possible d'employer cette méthode pour des applications du filtrage particulaire à la navigation inertielle. Dans le cadre de la navigation inertielle, l'équation d'état est usuellement linéaire avec des bruits gaussiens, et l'équation d'observation est non linéaire avec des bruits qui ne sont pas nécessairement gaussiens, mais ne mettent en jeu que certaines composantes de l'état : essentiellement les erreurs d'estimation inertielle en position, ces composantes sont par exemple désignées par le terme d'états « non linéaires » XNL et les autres composantes de l'état par le terme « d'états linéaires » XL. II est à remarquer en premier lieu que la loi conditionnelle de XL sachant les observations passées et sachant les états non linéaires passés XNL ne dépend en fait que des états non linéaires passés, et que cette loi est gaussienne, donnée par un filtre de Kalman dont les observations sont les états non linéaires. Il est à remarquer ensuite que la loi jointe des états non linéaires peut s'exprimer à l'aide d'innovations. La méthode de Rao-Blackwell tire avantage de ces remarques et approche la loi conditionnelle jointe des états linéaires et des états non linéaires passés sachant les observations passées à l'aide d'une approximation particulaire, sur un espace d'état de dimension réduite, de la loi conditionnelle des états XNL passés sachant les observations passées, et associe à chaque particule, qui représente donc une trajectoire possible des états « non linéaires », un filtre de Kalman pour prendre en compte la loi conditionnelle de XL sachant les états XNL passés. Par rapport à une approximation particulaire classique, la méthode de Rao-Blackwell réduit la variance de l'erreur d'approximation et n'utilise l'approximation particulaire que sur un espace d'état de dimension réduite, ce qui augmente d'autant l'efficacité de l'algorithme en minimisant son coût calculatoire. Les principaux avantages de l'approche particulaire résident dans la prise en compte immédiate des non linéarités, des contraintes portant sur l'état caché ou sur son évolution, des statistiques de bruit non gaussiennes, etc., et dans leur grande facilité de mise en oeuvre : il suffit en effet de savoir simuler des réalisations indépendantes du modèle d'état caché, et de savoir calculer la fonction de vraisemblance en un point quelconque de l'espace d'état, cette dernière restriction pouvant d'ailleurs être levée. De plus les grandes incertitudes sur l'état produisent des multi-modalités qui sont naturellement prises en compte par ce type d'approximation, ce qui est caractéristique de la phase de convergence. Les principaux inconvénients de cette approche résident essentiellement dans un temps de calcul plus important, mais qui se réduit essentiellement à la répétition, sur un échantillon de grande dimension de l'ordre de quelques centaines à quelques milliers, de tâches par ailleurs très simples de complexité comparable voire inférieure à la complexité d'un filtre de Kalman étendu, ainsi que dans l'introduction d'aléa supplémentaires pour la simulation (lors de l'étape de prédiction) et pour la redistribution. Il est donc indiqué de réduire à tout instant au minimum le nombre de particules et l'aléa supplémentaire introduit. La présente invention propose de combiner astucieusement les deux familles de filtres, Kalman et particulaire, en sélectionnant le meilleur 20 algorithme dans chacune des phases de convergence et de poursuite. Pour une meilleure compréhension de l'invention, l'algorithme de poursuite de type UKF est développé ci-après, et choisi à titre d'exemple non limitatif de l'invention. 25 Le filtrage de type UKF permet de s'affranchir de la linéarisation des fonctions f et h, en obtenant directement à partir de f et h des approximations numériques de la moyenne et de la covariance de l'état dynamique du système. A cette fin, un ensemble de points désignés "sigma-points" ou a-points est élaboré. Les fonctions f et h sont évaluées pour les a-points. Le 30 nombre de a-points créés dépend de la dimension d du vecteur d'état dynamique. Les a-points sont alors propagés dans le modèle de prédiction non linéarisé. La moyenne et la covariance sont alors calculées par moyenne pondérée à partir des a-points, chaque a-point X; étant associé à un poids co. 35 x =x±S ,/d+2 (11), où i~ Sx désigne la racine carrée de la matrice de covariance, Px , obtenue par décomposition de Cholesky, et ei est le i-ème vecteur de base. Les poids associés w(md) ,..., w(dm) pour l'estimation de la moyenne et w(e d ,..., wd~) pour l'estimation de la variance sont donnés par la relation suivante : w(m) _ 2 0 d+2 w(c) d+2+(1-a2+,Q) (12),OÙ: (m) = w(c) = 1 , i 0 2 est un coefficient d'échelle défini par : 2= a2 (d + K) - d , a, /3, K étant des paramètres à ajuster. En utilisant le paramétrage admissible a =1,,(3 = 0 qui rend identiques les poids pour le calcul de la moyenne et de la covariance, il est alors possible de fixer 2= K. On considère alors les poids : 1 w0 d+~ et w i =w =2(d+2) (13). Ainsi, pour l'étape de prédiction, propre à tout système de filtrage, les a-points x_ d,..., xd peuvent être définis par la relation : x0 =Xk_1 et x+i -k-1 ±Sk-1 .ei d+2 (14), où : Xk_1 désigne l'estimation a posteriori de l'état dynamique à l'instant k-1 ; Sk _ 1 est la racine carrée de la matrice de covariance de l'état à La construction des a-points proposée par Rudolph van der Merwe peut par exemple être utilisée. Selon cette construction, les a-points x_ d,..., xd peuvent être formulés suivant la relation suivante : x0 = x ±t l'instant k-1, Pk -1, obtenue par décomposition de Cholesky et ei est le ième vecteur de base. Le vecteur moyen d'état dynamique à l'instant k peut être formulé par la relation suivante : d Xk = E w fk(xi) (15). i = -d La matrice de covariance moyenne à l'instant k peut être formulée par la relation suivante : d -7 Pk =Qk+ E w fk(xi)-Xk)1(fk(xi)-X kJ (16). -d J a - Pour l'étape de correction, il est possible de définir les a-points 1 o x- d ,..., xd par la relation : T x+i =Xk ±Sk ei(17),où: Sk est la racine carrée de la matrice de covariance moyenne, Pk , 15 obtenue par décomposition de Cholesky et ei est le i-ème vecteur de base. La prédiction du vecteur moyen de mesure peut être formulée par la relation : k = E d cohk(xi) (18). i = -d La prédiction de la matrice de covariance moyenne "k à l'instant k 20 peut être formulée par la relation suivante : d ( (hk(x.) T "k=Rk + E wl hk(xi)-k ) -k (19). \ La matrice de corrélation peut être formulée par la relation suivante : d \ T Ck = E co xi - Xk I - (hk (x) - k I (20). i=-d / \ / L'estimation de l'état dynamique du système à l'instant k ainsi que la covariance associée peuvent alors être formulées par la relation suivante : L'algorithme peut être initialisé avec : ~ = EMXO j et PO = E[(XD - X0 )- (X0 - X0 ~J (22). Le paramètre d'échelle 2, représentatif de la dispersion des a-points autour de la valeur moyenne, peut par exemple être fixé à 0,9. Lorsque les incertitudes sont plus importantes, par exemple dans le cas où les performances de la centrale inertielle sont insuffisantes ou lorsque les informations remontant du terrain ne sont pas assez discriminantes, le filtre de poursuite 22 peut ne plus être en mesure de traiter correctement le problème de filtrage, car la distribution conditionnelle de l'état sachant la mesure peut alors devenir multi modale. Il s'agit alors d'utiliser un filtre de convergence plus apte à ramener les fortes incertitudes à des valeurs plus faibles. Une telle fonction de filtre de convergence est particulièrement bien traitée par les filtres de type particulaire. La suite de la description décrit en détails la technique de Rao-Blackwellisation, ou de marginalisation, à titre d'exemple non limitatif de l'invention, cette technique permettant d'optimiser le filtre particulaire utilisé dans l'invention pendant la phase de convergence. Dans le cas où la forme particulière des équations inertielles permet de Rao-Blackwelliser le filtre particulaire, l'état du système peut alors être décomposé en : - des composantes dites non linéaires désignées XNL, à savoir l'erreur inertielle en position et en attitude, qui apparaissent de manière non linéaire dans les équations d'observation, {Îk +Ck -k-1 Ykk PP -Ck - - 1 "k CT k k k 30 - des composantes dites linéaires désignées XL, à savoir l'erreur inertielle en vitesse et les biais, qui apparaissent de manière linéaire dans les équations d'état. Du point de vue mathématique, le cadre général du problème de filtrage est explicité ci-après. Le système d'évolution d'état peut être formulé suivant la relation suivante : f XkL =FkL'LXkL-1 +FkL,NXkN-1 +WkL - XkN FkN'LXkL_1 + FkN,NXkN_1 + WkN = Wk désigne le bruit de process, décomposé en composantes linéaires et non linéaires, considéré comme issu d'une distribution normale de (TITI, \ moyenne nulle et de matrice de covariance Qk : Wk = k -,N(O,Qk), la ~Wk i (23), où : matrice de covariance s'écrivant Qk = QkL QkL,N N,L N Qk Qk / Le vecteur d'observation peut se formuler suivant la relation suivante, ou "équation de mesure" : Y k hkN (xkN )+hkLxk +Vk (24), où : Vk désigne le bruit d'observation, issu d'une distribution de type bruit blanc Gaussien, de moyenne nulle et de covariance Rk : Vk -N(0, Rk) . L'estimation de la densité de probabilité a posteriori de l'état p(Xk ,Xk L /Yk) peut alors, en utilisant le théorème de Bayes, être formulée suivant la relation suivante : 25 p(Xk,XNL/Yk)=p(Xk /XNL,Yk)-p(XNL/Yk) (25). Kalman Particulaire20 La partie non linéaire de l'état dynamique du système est représentée à l'instant k par les particules , i=1,...,N, associées aux poids col et aux prédicteurs de Kalman mk pour la partie linéaire. L'initialisation peut être réalisée à partir de la matrice de covariance 5 initiale de l'état, en disposant les particules sur une grille régulière recouvrant l'incertitude à plus ou moins trois sigmas et à zéro pour la partie linéaire. La marginalisation peut alors se dérouler de la façon suivante, décomposée en différents items, pour chaque valeur de i de 1 à N : - Simulation de deux vecteurs aléatoires gaussiens indépendants Ul et Wl , centrés et de matrice de covariance PL/Net Qk ; et on pose : N i {~N N - Fk mk-1 +u»+ Jk k-1 + de sorte que le vecteur aléatoire moyenne Fk -mk-1+f(1) FN PL~N FN + N k k-1 (k) Qk 20 - Définition du vecteur moyenne suivant la relation : mil = FLmi + L (ei v LIN (- FNm - f ~) (27), k k-1 k k-1 / k k-1 k k-1 k k k-1 k k-1 et de matrice de covariance N Pklk 1 - \Fk L Pk-1 (Fk L )* + Qk )- Kkllk 1 (FNPLIN(FLY kk1 + Qk L ) (28), et de la matrice de covariance : où la matrice de gain K1 est définie par la relation suivante : KLIIN (FLPkII (FN)* +Qk N)FkNPkII (FkN)* +Q 1-1 klk-1 k (29). 10 15 k est gaussien, de vecteur 25 - Définition du poids, formulé par la relation suivante : (( r i LIN * kt/ ~k °~ q(Yk - hk \k )- Hkmkl k 1 HkPkl k-1Hk + Rk k-1 - Normalisation du poids. - Si la taille effective de l'échantillon devient inférieure à un seuil prédéterminé, les particules sont redistribuées selon un échantillonnage multinomial. - Définition du vecteur moyenne formulé suivant la relation suivante : l l mi =mk1k 1 +KkIN(ÿk -hk( k/-Hkmklk 11 (31), et de la matrice de covariance suivant la relation : pLIN =Pklk- pLiN 1 où la matrice de gain ((est définie par la relation : N Kk I N = Pk~k- Hk (HkFkl k 1Hk + Rk Y (33). Dans le but de rendre le filtre plus robuste, une étape de post- 20 régularisation peut avantageusement être réalisée après l'étape de correction particulaire par un noyau gaussien de paramètre p de régularisation, par exemple fixé à 0,2. Le filtre de navigation 11 proposé permet de tirer profit du meilleur des 25 deux filtres de convergence et de poursuite, par exemple particulaire et de Kalman 21, 22, selon le schéma présenté dans la figure 3 décrit ci-après. A un instant k du temps discrétisé, le filtre de navigation 11 reçoit une mesure, c'est-à-dire une observation de l'état dynamique du système Yk, ainsi que cela est illustré par une première étape 301 sur la figure. 30 La première étape 301 est suivie d'une étape de comparaison 302, lors de laquelle la valeur d'un indice de qualité calculé à partir de la matrice de covariance restituée par le filtre de navigation à l'instant discrétisé précédent k-1, est comparée à une valeur seuil prédéterminée. L'indice de qualité et la valeur seuil sont précisés en détails ci-après. 10 15 A l'issue de l'étape de comparaison 302, des moyens de basculement compris dans le filtre de navigation 11, permettent de forcer l'un ou l'autre des deux filtres de convergence et de poursuite 21, 22 à réaliser les calculs de l'état dynamique du système et de la matrice de covariance associée, pour l'instant k, ainsi que cela est illustré sur la figure par les étapes intermédiaires 303 à 306. Le filtre de navigation 11 restitue alors pour l'instant k, une estimation de l'état du système dynamique et la matrice de covariance associée ainsi que cela est représenté par la dernière étape 307 sur la figure, tels que restitués par le filtre 21, 22 sélectionné. Par exemple, si le filtre de navigation 11 restituait à l'instant k-1 les données calculées par le filtre de poursuite 22, et qu'à l'étape de comparaison 302, l'indice de qualité est inférieur à la valeur seuil prédéterminée, alors pour l'instant k, les moyens de basculement forcent le filtre de poursuite 22 à réaliser les calculs, et le filtre de navigation 11 restitue pour l'instant k l'état du système et la covariance tels que calculés par le filtre de poursuite 22. Si le filtre de navigation 11 restituait à l'instant k-1 les données calculées par le filtre de poursuite 22, et qu'à l'étape de comparaison 302, l'indice de qualité est supérieur à la valeur seuil prédéterminée, alors pour l'instant k, les moyens de basculement forcent le filtre de convergence 21 à réaliser les calculs, et le filtre de navigation 11 restitue pour l'instant k l'état du système et la covariance tels que calculés par le filtre de convergence 21. Ainsi, pour ce dernier exemple, si à l'instant k+1, l'indice de qualité reste supérieur à la valeur seuil prédéterminée, alors pour l'instant k+1, les moyens de basculement forcent le filtre de convergence 21 à réaliser les calculs, et le filtre de navigation 11 restitue également pour l'instant k+1 l'état du système et la covariance tels que calculés par le filtre de convergence 21. Si en revanche l'indice de qualité devient inférieur à la valeur seuil prédéterminée, alors pour l'instant k+1, les moyens de basculement forcent le filtre de poursuite 22 à réaliser les calculs, et le filtre de navigation 11 restitue pour l'instant k+1 l'état du système et la covariance tels que calculés par le filtre de poursuite 22. L'indice de qualité peut par exemple être la norme des termes d'une sous-matrice de covariance. Par exemple, l'indice de qualité peut être défini comme une norme des termes de la matrice de covariance relatifs à la position du porteur projetée dans un plan horizontal (c'est-à-dire sur les termes en x et y). Plus précisément, la norme peut être définie comme la somme des carrés des termes diagonaux de la sous-matrice de covariance comprenant les termes de position dans un plan horizontal. Ainsi, par exemple, si la trajectoire du porteur passe sur une longue zone de terrain plat, le filtre de convergence 21, de type particulaire, peut prendre le relais du filtre de poursuite 22, de type filtre de Kalman, afin de permettre une reconvergence dès la réapparition d'un signal exploitable pour la NCT. Par ailleurs, il est avantageusement possible d'adjoindre au filtre de poursuite 22, de type filtre de Kalman, un dispositif de rejet de mesures aberrantes à partir de l'analyse de l'innovation du filtre. L'innovation du filtre est définie comme l'écart entre la mesure obtenue et la mesure prédite (mesure effectuée en l'état prédit). Par exemple, si l'innovation dépasse une deuxième valeur seuil prédéterminée par rapport à l'écart type de la mesure prévue, c'est-à-dire l'écart type du bruit de mesure Vk la mesure peut être rejetée. Ce mode de réalisation permet de traiter les cas où l'erreur de carte ponctuelle est très grande et ne satisfait pas au modèle prévu pour la zone. Avantageusement encore, dans le but de permettre la détection de la divergence du filtre de poursuite 22, de type filtre de Kalman, une variable de type compteur peut être introduite, cette variable représentant le nombre de mesures consécutives qui ont été rejetées par le filtre de navigation 11. Lorsque ce nombre dépasse une troisième valeur de seuil prédéterminée, le filtre de navigation 11 suppose que le filtre de poursuite 22, de type filtre de Kalman, est « perdu », et peut automatiquement basculer en mode de filtrage particulaire via le filtre de convergence 21, en augmentant d'un certain facteur (supérieur à 1), par exemple un facteur trois, les incertitudes sur la position, afin de tenter de retrouver la position du porteur. Avantageusement encore, afin d'éviter au filtre de navigation 11 de confondre un véritable signal avec de l'erreur, par exemple dans le cas d'un terrain très plat, et ainsi de rendre le filtrage plus robuste, une quatrième valeur de seuil minimal peut être introduite sur l'écart type local des mesures impliquées dans la NCT. Si cet écart type est inférieur à la quatrième valeur de seuil, le filtre de navigation 11 peut alors utiliser uniquement les mesures de capteurs de bas niveau, dans le cas contraire il peut utiliser l'ensemble de la mesure pour réaliser la NCT. REVENDICATIONS1- Filtre de navigation (11) pour système de navigation par corrélation de terrain délivrant une estimation de l'état cinématique d'un porteur associé à une matrice de covariance à un instant discrétisé k à partir d'une pluralité de données comprenant les mesures restituées par au moins un capteur de terrain (102), le modèle associé au capteur de terrain (102), les données d'une carte embarquée (12), un modèle d'erreur de la carte embarquée (12), les mesures restituées par une centrale inertielle (101) et un modèle de la centrale inertielle (101), le filtre de navigation (11) étant caractérisé en ce qu'il comprend un premier filtre dit filtre de convergence (21), et un second filtre dit filtre de poursuite (22), le filtre de navigation (11) comprenant des moyens de basculement sélectionnant pour le calcul de l'état cinématique du porteur et de la matrice de covariance, l'un ou l'autre desdits filtres de convergence (21) et filtre de poursuite (22), en fonction de la comparaison (302) d'un indice de qualité calculé à partir de la matrice de covariance restituée par le filtre de navigation (11) à l'instant précédent k-1, à une valeur seuil prédéterminée, le filtre de navigation (11) restituant les valeurs calculées par le filtre (21, 22) sélectionné. 2- Filtre de navigation (11) selon la revendication 1, caractérisé en ce que l'indice de qualité est égal à une norme des termes d'une sous-matrice de la matrice de covariance. 3- Filtre de navigation (11) selon la revendication 1, caractérisé en ce que l'indice de qualité est égal à une norme des termes de la sous-matrice de covariance comprenant les termes de position (x,y) du porteur dans un plan horizontal. 4- Filtre de navigation (11) selon la revendication 1, caractérisé en ce que l'indice de qualité est égal à la somme des carrés des termes diagonaux de la sous-matrice de covariance comprenant les termes de position (x,y) du porteur dans un plan horizontal. 5- Filtre de navigation (11) selon l'une quelconque des revendications précédentes, caractérisé en ce que le filtre de convergence (21) est un filtre de type particulaire. 6- Filtre de navigation (11) selon l'une quelconque des revendications précédentes, caractérisé en ce que le filtre de convergence (21) est un filtre de type particulaire marginalisé. 7- Filtre de navigation (11) selon l'une quelconque des revendications précédentes, caractérisé en ce que le filtre de poursuite (22) est un filtre de type filtre de Kalman étendu EKF. 8- Filtre de navigation (11) selon l'une quelconque des revendications précédentes, caractérisé en ce que le filtre de poursuite (22) est un filtre de type UKF. 9- Filtre de navigation (11) selon l'une quelconque des revendications précédentes, caractérisé en ce que le filtre de poursuite (22) est associé à un dispositif de rejet de mesures aberrantes du filtre de navigation (11), rejetant les mesures si le rapport entre l'innovation et l'écart type du bruit de mesure dépasse une deuxième valeur de seuil prédéterminée. 10- Filtre de navigation (11) selon la revendication 9, caractérisé en ce que les moyens de basculement restituent en sortie du filtre de navigation (11) les données restituées par ledit filtre de convergence (21) lorsque le nombre de mesures consécutives rejetées par le dispositif de rejet est supérieur à une troisième valeur de seuil prédéterminée. 11- Filtre de navigation (11) selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre des moyens de comparaison de l'écart type local des mesures de navigation par corrélation de terrain à une quatrième valeur de seuil prédéterminée, le filtre de navigation (11) se basant uniquement sur lesmesures fournies par le baro-altimètre (1013) compris dans la centrale inertielle (101) lorsque l'écart type est inférieur à ladite troisième valeur de seuil prédéterminée.