FR 3120437 A1 20220909 FR 2101987 A 20210302 Cette centrale de navigation inertielle comprend une unité d’acquisition de données de navigation (1) et un premier étage de navigation (2) apte à calculer des valeurs de navigation hybrides par fusion de données de navigation . Elle comporte en outre un deuxième étage de navigation (3) apte à calculer des valeurs de navigation certifiées indépendamment des valeurs de navigation hybrides, lesdites valeurs de navigation certifiées étant certifiées avec une limite d’erreur de géolocalisation Figure pour l’abrégé : Fig 1 Procédé et centrale de calcul de données de navigation inertielle La présente invention concerne les centrales de navigation inertielle et concerne plus particulièrement un procédé et une centrale de calcul de données de navigation inertielle. Etat de la technique antérieure Les centrales de navigation inertielle sont par exemple, mais non exclusivement, destinées à être utilisées pour la géolocalisation d’aéronefs, de bateaux, de véhicules terrestres, … pour la conduite de tirs, ou la désignation de cibles. De nombreux travaux ont été effectués pour l’amélioration de la précision des centrales inertielles de navigation. Les travaux d’amélioration sont fondés essentiellement sur l’amélioration des précisions de mesure des capteurs inertiels, tels que des accéléromètres et des gyroscopes, ou sur la fusion ou l’hybridation des données inertielles avec des données provenant d’autres capteurs ou avec des modèles comportementaux complémentaires. La fusion ou l’hybridation de données inertielles et de données provenant d’autres capteurs permet souvent d’améliorer considérablement les précisions typiques de géolocalisation d’une centrale inertielle de navigation. En revanche, l’hybridation des données ne permet souvent pas de calculer un indicateur statistique de précision des erreurs de géolocalisation fiable. Pour qu’un tel indicateur de précision d’erreurs de géolocalisation soit fiable, il faut pour cela connaître la statistique des erreurs des capteurs considérés. Ceci est généralement le cas pour les capteurs inertiels dont les fabricants peuvent garantir des statistiques d’erreurs en production. Mais ce n’est en revanche pas nécessairement le cas en ce qui concerne les autres capteurs utilisés dans les hybridations de données. En effet, ces techniques de fusion de données s’appuient souvent sur des modèles comportementaux du véhicule, en l’espèce de l’aéronef, ou de la scène observée, ou encore de bases de données limitées, ce qui ne permet pas d’associer une statistique d’erreur fiable à ce type de modèle. Par ailleurs, les techniques de fusion de données, bien que performantes, sont généralement non sécurisées et ne permettent pas de garantir une statistique d’erreurs dans des cas de dysfonctionnements, qui ne sont pas dus à des pannes des capteurs, mais à des emplois non robustes, c’est-à-dire en dehors des plages de fonctionnement sécurisées. Dans ces cas, il ne suffit pas d’associer plusieurs centrales inertielles redondantes car il s’agit d’un cas dysfonctionnel commun au système de navigation intégré dans le véhicule. Par ailleurs, les données de géolocalisation des centrales de navigation peuvent être utilisées dans des chaînes fonctionnelles critiques, que ce soit par exemple pour la conduite autonome d’un véhicule, tel qu’un drone, un aéronef, un train, un véhicule automobile, un bateau, un sous-marin, …, ou encore la conduite d’un tir de canon ou la désignation de cible dans un cas d’emploi militaire. Dans ces cas, il est nécessaire de pouvoir certifier la centrale de navigation avec un niveau de criticité suffisant en face du risque de sûreté associé, et de pouvoir fournir une limite statistique des données de géolocalisation avec une probabilité d’avoir une erreur de navigation en dehors de cette limite, qui soit compatible avec les exigences de sûreté de la mission. On connaît, dans l’état de la technique, et en particulier dans le domaine de l’aéronautique, des centrales inertielles de navigation utilisant une hybridation avec radionavigation et assurant la fusion de données inertielles avec des données de géolocalisation provenant de systèmes GPS, GLONASS, … Les statistiques d’erreurs sont définies dans le document GLOBAL POSITIONING SYSTEM STANDARD POSITIONING SERVICE PERFORMANCE STANDARD [DoD & GPS NavSTAR – 2020], en ce qui concerne le GPS, par exemple, et sont utilisées pour définir des rayons de protection d’erreurs de position de la navigation inertielle hybridée GPS, même en cas de panne d’un ou de plusieurs satellites. On pourra également se référer au document US 5, 760, 737 ou au document WO 2012 031940 qui décrivent des exemples d’applications de cette technique d’hybridation avec radionavigation. Mais ces techniques de calcul de rayon de protection des erreurs de positions ne sont valables qu’en cas d’emploi avec un système de positionnement global GPS et ne traitent pas des problématiques de type erreur de multitrajet au sol, leurrage, voire cas d’emploi sans signal GNSS (brouillage permanent, …). Le but de l’invention est donc de pallier les inconvénients précités et de fournir un procédé et une centrale de calcul de données de navigation inertielle qui soient capables de délivrer des données de géolocalisation certifiées en fournissant des données de géolocalisation avec un indicateur de précision associé à une probabilité d’erreur. Un autre but de l’invention est d’utiliser des données provenant d’une navigation inertielle hybride non certifiée dans une chaîne fonctionnelle de navigation inertielle critique. L’invention a donc pour objet, selon un premier aspect, un procédé de calcul de données de navigation inertielle, dans lequel : -on acquiert des données de navigation, et - on calcule des valeurs de navigation hybrides à partir de ces données de navigation. Selon ce procédé, on calcule des valeurs de navigation certifiées, indépendamment des valeurs de navigation hybrides, lesdites valeurs de navigation certifiées étant certifiées avec une limite d’erreur de géolocalisation. Dans un mode de mise en œuvre, on calcule, pour les valeurs de navigation certifiées, une première limite d’erreur de géolocalisation admissible correspondant à une géolocalisation avec une probabilité d’erreurs prédéterminée. L’invention permet ainsi de calculer, en parallèle d’une navigation inertielle hybride performante mais non certifiée et sans limite de protection fiable des erreurs, fondée sur l’hybridation des données inertielles, une navigation inertielle plus simple et certifiée, qui permet de calculer une limite de protection fiable des données de navigation inertielle hybride, sans avoir à certifier cette navigation hybride. Par limite de protection, on entend une limite statistique des données de géolocalisation avec une probabilité d’avoir une erreur de navigation en dehors de cette limite. Les valeurs de navigation hybrides sont certifiées en calculant une limite de protection basée sur une comparaison avec une navigation inertielle sécurisée et certifiée, non impactée par des capteurs externes ou modèles mathématiques comportementaux incohérents. Ainsi, les valeurs de navigation hybride peuvent être intégrées et exploitées dans une chaîne fonctionnelle critique certifiée. Dans un mode de mise en œuvre, on calcule, pour les valeurs de navigation hybride, une deuxième limite d’erreur de géolocalisation admissible à partir de la première limite d’erreur de géolocalisation et d’un écart entre les valeurs de géolocalisation hybrides et les valeurs de géolocalisation certifiées. Selon une autre caractéristique du procédé selon l’invention, on calcule un statut de validité des valeurs de navigation hybrides par comparaison avec des valeurs de seuil respectives. Dans un mode de mise en œuvre, lors du calcul des valeurs de navigation certifiées, on utilise un filtre de Kalman pour le calcul de covariances d’erreur à partir des données de navigation. Par exemple, le filtre de Kalman est un filtre de Kalman invariant. Selon encore une autre caractéristique du procédé selon l’invention, on procède à un recalage à vitesse nulle d’un étage de navigation inertielle apte à calculer les valeurs de navigation certifiées. On peut prévoir que les valeurs de navigation hybrides soient calculées à partir de données inertielles issues de capteurs, notamment de trois gyroscopes et de trois accéléromètres, dans trois directions de l’espace, et de données additionnelles, notamment de modèles de déplacement. L’invention a également pour objet une centrale de navigation inertielle, comprenant une unité d’acquisition de données de navigation et un premier étage de navigation apte à calculer des valeurs de navigation hybrides par fusion de données de navigation. Cette centrale comporte en outre un deuxième étage de navigation apte à calculer des valeurs de navigation certifiées, indépendamment des valeurs de navigation hybrides, lesdites valeurs de navigation certifiées étant certifiées avec une limite d’erreur de géolocalisation. D’autres buts, caractéristiques et avantages de l’invention apparaîtront à la lecture de la description suivante, donnée à titre d’exemple nullement limitatif et faite en référence aux dessins annexés sur lesquels : est un schéma de l’architecture matérielle d’une centrale de navigation inertielle conforme à l’invention ; illustre les principales étapes d’un procédé de calcul de données de navigation inertielle conforme à l’invention montre un schéma illustrant les première et deuxième données de navigation avec leur limite d’erreur de géolocalisation admissible sous la forme de cercles de protection limite de l’erreur de position ; et Exposé détaillé d’au moins un mode de réalisation de l’invention On se réfère tout d’abord à la qui illustre schématiquement l’architecture d’une centrale de navigation inertielle conforme à l’invention. Cette centrale comporte une unité UMI 1 d’acquisition de données navigation inertielles comprenant des capteurs de mesure, par exemple trois gyroscopes et trois accéléromètres délivrant des données inertielles dans trois directions de l’espace, un premier étage de navigation hybridée 2 destiné à calculer des valeurs de navigation inertielle hybrides à partir des données de navigation inertielles délivrées par l’unité d’acquisition 1 et par des capteurs complémentaires ou modèles comportementaux 4 fournissant des données de navigation complémentaires, telles que des données de géolocalisation GPS, des données de distances parcourues, un modèle de déplacement du véhicule,…et un deuxième étage de navigation 3 fournissant des valeurs de navigation sécurisées. Les premier et deuxième étages de navigation inertielle sont chacun constitués par un calculateur dûment programmé pour mettre en œuvre les fonctions qui vont être décrites ci-après. Ces calculateurs peuvent toutefois être formés par des partitions de calcul d’un même calculateur. Le premier et le deuxième étages de calcul assurent précisément le calcul de valeurs de données de navigation Le premier étage de navigation inertielle 2 assure un calcul des données de navigation hybrides, par exemple la position, la vitesse et l’altitude, en trois dimensions par hybridation, c’est-à-dire par fusion des données délivrées par l’unité 1 d’acquisition de données inertielles et des données d’aide externes fournies par les capteurs ou les modèles comportementaux 4. Ce premier étage de calcul comporte un calculateur 5, qui assure le calcul de ces données de navigation hybride. Comme cela sera décrit en détail par la suite, il reçoit également des données de recalage de position et d’initialisation P. Le deuxième étage de calcul comporte un calculateur 6, qui assure le calcul des données de navigation certifiées indépendamment des données de navigation hybrides fournies par l’unité d’acquisition 1 et un module de calcul additionnel 7 destiné au calcul de limites de protection des erreurs de géolocalisation. Le calculateur 6 reçoit, en entrée, des données de réglage de position, de recalage à vitesse nulle et d’initialisation P’ . Le calcul de la navigation hybride est une fonction permettant en premier lieu de calculer une solution de géolocalisation inertielle (position 3D, vitesse 3D, attitudes 3D) à l’aide de la fusion des mesures de l’UMI et de données d’aides externes ou modèles comportementaux du porteur (récepteur radionavigation, odomètre et modèle de déplacement véhicule, capteur loch (pour un porteur marine), recalage à vitesse nulle (ZUPT), capteur optique, lidar etc…) par exemple à l’aide de techniques telles que l’intégration d’une navigation inertielle couplé à un filtre de Kalman qui permet de calculer des covariances d’erreurs à partir des caractéristiques des capteurs inertiels et des capteurs d’aide ou encore à l’aide d’autres techniques telles que du filtrage particulaire ou réseaux de neurones et intelligence artificielle.. Le calcul de la navigation hybride est une fonction permettant en second lieu d’estimer la matrice 3x3 de rotation de passage du repère géographique local « North-East-Down » [N] au repère « Body » [B] lié au véhicule : la notation désignant une estimation de la grandeur , étant donc entaché d’erreur Il s’agit encore d’une fonction permettant l’estimation des composantes de vitesses du véhicule dans le repère géographique local « North-East-Down » [N], l’estimation des coordonnées de position horizontale du véhicule exprimées selon le modèle ellipsoïde WGS84, par exemple la longitude estimée et la latitude estimée , l’estimation de l’altitude du véhicule exprimée par rapport au géoïde. Le calcul de la navigation hybride assure également le calcul d’une matrice de covariance des erreurs de géolocalisation si un filtre de Kalman est utilisé, s’appuyant sur la modélisation des erreurs de géolocalisation. Des erreurs de géolocalisation de la navigation hybride sont en outre définies par modélisation mathématique. Il s’agit, par exemple d’erreurs d’angle, dans 3 dimensions, de la matrice suivante, projetée dans le repère Body [B] : , définie par les relations suivantes : Avec la matrice antisymétrique 3x3 du vecteur exprimé dans le repère [B] ; Et où : « logSO3 » est la carte logarithmique du groupe des matrices de rotation SO(3) vers les vecteurs de rotations R3 (R au sens de l’ensemble des réels) et « expSO3 » est la carte exponentielle des vecteurs de rotation R3 du groupe des matrices de rotation SO(3), SO(3) étant le groupe spécial orthogonal des matrices carrées de dimension 3 Les estimations des erreurs de vitesses tridimensionnelles, entachées d’erreur, délivrées par le premier étage de navigation hybridée 2, et exprimées dans le repère géographique local [N] par , sont définies par la relation suivante : Les erreurs de position horizontale bidimensionnelles exprimées dans le plan « North-East » du repère géographique local [N] : et , sont définies par les relations suivantes : On peut également définir l’erreur de position horizontale radiale comme étant la norme de ces deux composantes Avec le rayon terrestre selon le modèle WGS 84 dans la direction Nord Avec le rayon terrestre selon le modèle WGS 84 dans la direction Est Les erreurs d’altitude délivrées par le premier étage de navigation inertielle hybridée 2 qui correspondent à l’erreur de position selon l’axe vertical du repère géographique local [N] sont définies par Cette fusion de données permet d’améliorer la précision de la géolocalisation en réduisant les erreurs présentées ci-dessus, d’estimer et de compenser des erreurs des capteurs inertiels, d’estimer et de compenser des erreurs des aides externes, et d’estimer un indicateur statistique de performance, par exemple une matrice de covariance des erreurs d’un filtre de Kalman. En ce qui concerne le deuxième étage de navigation 3, le calcul de la navigation certifiée, appelée sécurisée, est une fonction permettant de calculer une solution de géolocalisation inertielle (position 3D, vitesse 3D, attitudes 3D) à l’aide de la fusion des mesures de l’unité UMI 1 et de données d’aides externes sécurisées, par exemple des données de recalage à vitesse nulle à l’arrêt du véhicule, ou recalage ZUPT, sécurisé par un indicateur d’arrêt transmis au système supposé être au niveau adéquat de sureté de fonctionnement ; des données de recalage de position ponctuel ou tout recalage utilisant un capteur intègre pour lequel dont on sait associer une probabilité de défaillance statistique et dont on sait borner les erreurs, hors défaillance, statistiquement et indépendamment des cas d’emplois considérés. Ces calculs de fusion de données sont réalisés grâce à l’intégration d’une navigation inertielle couplé à un filtre de Kalman qui permet de calculer des covariances d’erreurs à partir des caractéristiques des capteurs inertiels. Le deuxième étage de calcule de navigation 3 assure le calcul de l’estimation de la matrice 3x3 de rotation de passage du repère géographique local « North-East-Down » [N] au repère « Body » [B] lié au véhicule : La notation désignant une estimation de la grandeur , étant donc entaché d’erreur Il assure en outre l’estimation des composantes tridimensionnelles de la vitesse du véhicule dans le repère géographique local « North-East-Down » [N], ainsi que l’estimation des coordonnées de position horizontale du véhicule exprimées selon le modèle ellipsoïde WGS84 par exemple la longitude estimée et la latitude estimée Le deuxième étage de calcule de navigation 3 calcul en outre l’estimation de l’altitude du véhicule exprimée par rapport au géoïde, ainsi que le calcul d’une matrice de covariance des erreurs de géolocalisation fournie par le filtre de Kalman s’appuyant sur la modélisation des erreurs de géolocalisation. Cette matrice « » de dimension 9x9 représente l’estimation des covariances des 9 erreurs de géolocalisation choisies, à savoir: Les erreurs d’angle tridimensionnelles de la matrice projetées dans le repère Body [B] : , définies par Les erreurs de vitesse tridimensionnelles exprimées dans le repère géographique local [N] : , définies par Les erreurs de position horizontale bidimensionnelles exprimées dans le plan « North-East » du repère géographique local [N] : et , définies par : et Les erreurs d’altitude qui correspond à l’erreur de position selon l’axe vertical du repère géographique local [N] : , définies par Il est également possible de définir l’erreur de position horizontale radiale comme étant la norme de ces deux composantes : Avec : le rayon terrestre selon le modèle WGS 84 dans la direction Nord, et le rayon terrestre selon le modèle WGS 84 dans la direction Est Cette fusion de données permet d’améliorer la précision de la géolocalisation en réduisant les erreurs présentées ci-dessus, d’estimer et de compenser des erreurs des capteurs inertiels, d’estimer un indicateur statistique de performance, par exemple une matrice de covariance des erreurs d’un filtre de Kalman, qui n’est pas entaché d’erreurs liées à des mauvais cas d’emplois des aides externes puisqu’aucune aide externe complexe et non sécurisée n’est utilisée. Cet indicateur est uniquement entaché d’erreurs liées à une défaillance du système mais cette défaillance est connue et faible et correspond à la probabilité de défaillance non détectée du système. Cette géolocalisation est simple, sécurisée mais tout de même performante, à l’aide des recalages ZUPT ou de position ponctuels, pour permettre un rayon de protection sécurisé exploitable en opération. Ces recalages ponctuels permettent de maintenir une précision de géolocalisation sécurisée acceptable pour la mission. Le deuxième étage de navigation 3 assure par ailleurs le calcul de la limite de protection des données de géolocalisation et de validité des données hybrides. Ce calcul permet en particulier de définir une limite statistique des données de géolocalisation avec une probabilité d’avoir une erreur de navigation en dehors de cette limite. Cette limite de protection, notée « PL », des erreurs des données de géolocalisation de la navigation sécurisée est calculée en utilisant la matrice de covariance « Psafe » des erreurs de géolocalisation du filtre de navigation sécurisé (matrice de covariance fournie par le filtre de Kalman). Elle contient les variances des erreurs de géolocalisation et leurs corrélations entre elles. Cette matrice de covariance est le résultat de la propagation des caractéristiques d’erreurs des senseurs inertiels (statistique gaussienne des bruits, erreurs de biais, erreurs de facteurs d’échelle des senseurs, erreurs de mésalignements entre les senseurs, …) vers les erreurs de géolocalisation et des prises en compte des recalages vitesse nulle (ZUPT) ou de position absolue par le filtre de navigation La probabilité statistique d’être en en dehors de cette limite de protection en l’absence de panne d’un des senseurs inertiels, notée « Proba » est par exemple : Proba = 1e-5 Dans ce qui suit, des exemples de calcul de limite de protection vont être détaillés. En premier lieu, la limite de protection peut être définie par un cercle de protection limite de l’erreur de position horizontale : Où est la valeur propre maximale de la matrice de covariance 2x2 « » des erreurs de position horizontale Nord et Est du filtre de Kalman de la navigation sécurisée est l’inverse de la fonction cumulative de distribution de la variable aléatoire gaussienne normée de dimension N et décorrélée pour une probabilité Par exemple, pour N = 2, ln(x) = logarithme népérien évalué en x Pour N = 2, Proba = 1e-5et Qinv(1e-5 , 2) soit environ 4.7985 En l’absence de panne de la centrale de navigation, cette limite HPL(Proba) ainsi calculée est supérieure à avec une probabilité au moins égale à « 1-Proba » En l’absence de panne de la centrale de navigation, ce rayon de protection limite HPL(Proba) ainsi calculé, centré sur la position horizontale vraie du véhicule doit contenir avec une probabilité au moins égale à « 1-Proba » la position calculée par la navigation sécurisée. La limite de protection peut également être appliquée à l’erreur de position verticale (altitude) : Où est la variance de l’erreur d’altitude du filtre de Kalman de la navigation sécurisée et est l’inverse de la fonction cumulative de distribution de la variable aléatoire gaussienne normée de dimension N et décorrélée pour une probabilité Par exemple, pour N = 1, Où erfcinv est l’inverse de la fonction d’erreur complémentaire de Gauss Pour N = 1, Proba = 1e-5 et Qinv(1e-5 , 1) ≈ 4.4172 En l’absence de panne de la centrale de navigation, cette limite ZPL(Proba) ainsi calculée est supérieure à avec une probabilité au moins égale à « 1-Proba » La limite de protection des erreurs peut également être appliquée aux composantes de vitesse nord, Est et verticale dans le repère géographique local : Où est la variance de l’erreur de vitesse selon l’axe nord du repère géographique local du filtre de Kalman de la navigation sécurisée et est l’inverse de la fonction cumulative de distribution de la variable aléatoire gaussienne normée de dimension N et décorrélée pour une probabilité En l’absence de panne de la centrale de navigation, cette limite VnPL(Proba) ainsi calculée est supérieure à avec une probabilité au moins égale à « 1-Proba ». De même pour respectivement les erreurs de vitesse Est et verticale, on a : La limite de protection des erreurs d’angle tridimensionnelles projetées dans le repère Body [B], qui est un repère d’intérêt pour ces erreurs dans un cas de désignation de cible notamment, de la matrice , s’écrit : est la variance du filtre de Kalman de la navigation sécurisée de la composante selon l’axe x du repère Body [B] des erreurs d’angle de la matrice et est l’inverse de la fonction cumulative de distribution de la variable aléatoire gaussienne normée de dimension N et décorrélée pour une probabilité En l’absence de panne de la centrale de navigation, cette limite CbxPL(Proba) ainsi calculée est supérieure à avec une probabilité au moins égale à « 1-Proba » De même, pour respectivement les erreurs d’angle projetée dans le repère Body [B] de la matrice selon les axes YB et ZB, on a : Le deuxième étage de navigation certifiée 3 assure par ailleurs le calcul d’une limite de protection, noté « PL_hybrid », des erreurs des données de géolocalisation de la navigation hybride en utilisant les limites de protection PL des erreurs des données de géolocalisation de la navigation sécurisée associés à la probabilité « Proba » tel que décrit précédemment, et l’écart entre les données de géolocalisation de la navigation hybride et de la navigation sécurisée. Diverses limites de protection, qui s’appuient sur les calculs des limites de protection définies pour les données de géolocalisation sécurisée, peuvent être calculées. On peut en premier lieu calculer un cercle de protection limite de l’erreur de position hybride dans le plan horizontal du repère géographique local : HPL_hybrid(Proba) = HPL(Proba) + ∆Pos Où : HPL(Proba) est le cercle de protection limite des erreurs de position horizontales de la navigation sécurisée présenté auparavant, et est la distance entre les deux positions sécurisées et hybrides dans le plan horizontal du repère géographique local [N], Avec et En l’absence de panne de la centrale de navigation, HPL_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba » En effet, par construction (inégalité de Cauchy Schwartz), Or, il y a une probabilité supérieure ou égale à « 1-Proba », par choix de HPL(Proba), que . On sait donc garantir par le choix de HPL_hybrid(Proba) qu’il y a une probabilité supérieure ou égale à « 1-Proba » que En l’absence de panne de la centrale de navigation, ce rayon de protection limite HPL_hybrid(Proba) ainsi calculé, centré sur la position horizontale vraie du véhicule doit contenir avec une probabilité au moins égale à « 1-Proba » la position calculée par la navigation hybride. On peut en second lieu calculer une limite de protection de l’erreur de position verticale (altitude) hybride : ZPL_hybrid(Proba) = ZPL(Proba) + ∆h Où HPL(Proba) est le cercle de protection limite des erreurs de position horizontales de la navigation sécurisée présenté auparavant est la distance entre les deux altitudes sécurisée et hybride. En l’absence de panne de la centrale de navigation, ZPL_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba » En effet, par construction (inégalité de Cauchy Schwartz), Or, il y a une probabilité supérieure ou égale à « 1-Proba », par choix de ZPL(Proba), que . On sait donc garantir par le choix de ZPL_hybrid(Proba) qu’il y a une probabilité supérieure ou égale à « 1-Proba » que On peut également calculer une limite de protection des erreurs de vitesse hybride Nord, Est et Verticale dans le repère géographique local. Par exemple, par le même principe que décrit ci-dessus, on peut définir, pour la limite de protection des erreurs de vitesse hybride Nord une limite de protection de l’erreur , à partir de la relation : VnPL_hybrid(Proba) = VnPL(Proba) + ∆Vn Où: VnHPL(Proba) est la limite de protection de l’erreur de vitesse selon l’axe nord du repère géographique local [N], et est l’écart entre les deux vitesses nord sécurisée et hybride En l’absence de panne de la centrale de navigation, Vn_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba » On peut également définir une limite de protection de l’erreur : VePL_hybrid(Proba) = VePL(Proba) + ∆Ve Où : VeHPL(Proba) est la limite de protection de l’erreur de vitesse selon l’axe Est du repère géographique local [N] ,et = écart entre les deux vitesses Est sécurisée et hybride En l’absence de panne de la centrale de navigation, Ve_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba » On peut encore définir une limite de protection de l’erreur : VzPL_hybrid(Proba) = VzPL(Proba) + ∆Vz Où : VzHPL(Proba) est la limite de protection de l’erreur de vitesse selon l’axe vertical du repère géographique local [N] , et est l’écart entre les deux vitesses verticales sécurisée et hybride En l’absence de panne de la centrale de navigation, Vz_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba » On peut encore calculer une limite de protection des erreurs d’angle tridimensionnelles de la matrice projetée dans le repère Body [B]. A titre d’exemple, la limite de protection de l’erreur est : Où : CbxPL_hybrid(Proba) est la limite de protection de la composante selon XB des erreurs d’angle de la matrice et En l’absence de panne de la centrale de navigation, CbxPL_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba » En effet, on a : Soit encore Et donc Or pour des angles d’amplitude faible, par exemple de norme inférieure à 10e-3 radian, on a, au premier ordre : On en déduit pour des angles d’amplitude faible, par exemple de norme inférieure à 10e-3 radian, au premier ordre : Et donc Par unicité entre matrice antisymétrique d’un vecteur et ce même vecteur, en en déduit : On a alors l’inégalité suivante pour la composante selon XB : Or, il y a une probabilité supérieure ou égale à « 1-Proba », par choix de CbxPL(Proba), que On sait donc garantir par le choix de CbxPL_hybrid(Proba) qu’il y a une probabilité supérieure ou égale à « 1-Proba » que De même, les limites de protection des erreurs et sont calculées. Par exmple, en ce qui concerne , cette limite est calculée à partir de la relation suivante: Où : CbyPL_hybrid(Proba) est la limite de protection de la composante selon YB des erreurs d’angle de la matrice , et En l’absence de panne de la centrale de navigation, CbyPL_hybrid(Proba) ainsi calculé est supérieur à avec une probabilité supérieure ou égale à « 1-Proba ». La limite CbzPL_hybrid(Proba) est calculée de manière similaire. En référence à la , qui illustre les principales phases d’un procédé de calcul de données de navigation inertielles conforme à l’invention, selon une première étape 10, l’unité d’acquisition de données inertielles 1 acquiert les données en provenance d’au moins trois gyroscopes et d’au moins trois accéléromètres inertiels qui mesurent les rotations inertielles et les forces spécifiques inertielles dans les trois directions de l’espace. De même, au cours de cette étape 10, les données provenant des autres capteurs 4 sont également acquises. On notera que les capteurs de l’unité d’acquisition de données inertielles 1 sont développés et certifiés en fonction des contraintes de sécurité requises, en fonction de la criticité des évènements redoutés. Lors de l’étape 11 suivante, le premier étage de navigation 2 délivre les premières données de navigation par fusion des données inertielles issues de l’unité d’acquisition de données inertielles 4 et de toutes autres données de navigation complémentaires provenant des capteurs additionnels 4. On notera qu’à l’issue de cette étape 11 de calcul, les données de navigation fournies à l’issue de cette première étape de calcul sont relativement précises mais sont difficiles à certifier. Lors de l’étape 12 suivante, le deuxième étage de navigation 3 assure le calcul des données de navigation sécurisées à partir, notamment, des covariances d’erreur et à partir des caractéristiques des capteurs inertiels, comme décrit précédemment. On notera que ces étapes d’acquisition et de calcul sont réalisées en permanence. Toutefois, il est réalisé des phases de recalage de position, par exemple lors du passage à des endroits connus. On met également en œuvre des phases de recalage à vitesse nulle afin en particulier de recaler les données de vitesse. Lors de l’étape 13 suivante, on procède à un calcul de limites de protection qui correspond à un niveau d’erreur de protection, pour les premier et deuxième étages de navigation inertielle. En se référant également à la , sur laquelle P désigne une position réelle, et P1 et P2 désignent la position calculée par le premier et le deuxième étage de navigation 2 et 3, respectivement, au cours de cette étape, on calcule, comme indiqué précédemment, la limite de protection des erreurs des données de géolocalisation de la navigation sécurisée issue du deuxième étage de navigation inertielle. Cette limite de protection, notée HPL, constitue ici un cercle de protection limite de l’erreur de géolocalisation centrée sur la position P2 calculée par le deuxième étage de navigation 3. Lors de l’étape 14 suivante, on calcule la limite de protection des erreurs des données de géolocalisation calculées par le premier étage de navigation 2, tel qu’exposé précédemment. Cette limite de protection HPL_hybride forme un cercle de protection limite de l’erreur de position fournie par le premier étage de navigation 2, centré sur la position P1 géolocalisée par ce premier étage de navigation et dont le rayon est constitué par la somme de l’écart entre les données de géolocalisation fournies par les premiers et deuxièmes étages de navigation et de la limite de protection HPL calculée lors de l’étape 13 précédente. Lors de l’étape 15 suivante, on calcule un statut de validité V des premières données de navigation hybrides par comparaison avec des valeurs de seuil respectives. Tout d’abord, un test de cohérence des données de géolocalisation hybride est préalablement effectué afin de s’assurer de la cohérence de la donnée, tel que la validité déclarée des données, le format de la donnée, l’intervalle de la donnée, etc…. Dans le cas contraire un statut d’invalidité de la donnée est émis. Par ailleurs, les écarts mesurés entre la navigation hybride et la navigation optimale calculés précédemment sont également comparés à des seuils afin de lever des statuts d’anomalie vers l’utilisateur et l’informer de ne plus considérer les données de navigation hybride. On peut par exemple utiliser une fraction des limites de protections « PL » des erreurs de la navigation sécurisée comme seuil de levée d’anomalie. Le calcul de validité peut être du type : Si test de cohérence des données de géolocalisation hybride est invalide Ou si ∆Pos > Seuil_pos Ou si ∆h > Seuil_Z Ou si ∆Vn > Seuil_Vn Etc… Alors les données de la navigation hybride sont invalides Sinon les données de la navigation hybride sont valides Fin Si On notera que les calculs mis en œuvre au sein de la centrale de navigation sont réalisés au sein de deux calculateurs séparés, ou de partitions d’un calculateur, assurant respectivement le calcul des données de navigation hybride et le calcul des données de navigation sécurisée, ce qui permet de ne pas contraindre les couches basses logicielles et matérielles de l’organe de calcul de la navigation hybride qui n’est pas dans la chaine fonctionnelle critique de l’invention. La solution proposée permet d’utiliser les données d’une navigation hybride non certifiée quelconque dans une chaine fonctionnelle critique grâce à la sécurisation de ces données par une partition logicielle ou matérielle certifiée qui calcule pour ces données de navigation non certifiée une limite de protection basée sur la comparaison avec une navigation sécurisée, non impactée par des capteurs externes ou modèles mathématiques comportementaux incohérents et certifiée. Ainsi, toute technique d’hybridation de la navigation (filtre de Kalman complexe, filtre particulaire, réseau de neurone, intelligence artificielle, …) peut être utilisée sans avoir à certifier ces traitements, tout en garantissant un niveau d’intégrité des données grâce au calcul des limites de protection fiables et certifiées. Les données de la navigation sécurisée ainsi que leurs limites de protection sont également disponibles en solution de secours, utilisées par exemple lorsque l’écart entre les données de géolocalisation de la navigation hybride et sécurisée est supérieur à un seuil par exemple. Procédé de calcul de données de navigation inertielle, dans lequel : -on acquiert des données de navigation ; et -on calcule des valeurs de navigation hybrides à partir de ces données de navigation, caractérisé en ce que : on calcule des valeurs de navigation certifiées, indépendamment des valeurs de navigation hybrides, lesdites valeurs de navigation certifiées étant certifiées avec une limite d’erreur de géolocalisation. Procédé selon la revendication 1, dans lequel on calcule, pour les valeurs de navigation certifiées, une première limite d’erreur de géolocalisation admissible correspondant à une géolocalisation avec une probabilité d’erreur prédéterminée. Procédé selon la revendication 2, dans lequel on calcule, pour les valeurs de navigation hybrides, une deuxième limite d’erreur de géolocalisation admissible à partir de la première limite d’erreur de géolocalisation et d’un écart entre les valeurs de navigation hybrides et les valeurs de navigation certifiées. Procédé selon l’une quelconque des revendications 1 à 3, dans lequel on calcule un statut de validité (V) des valeurs de navigation hybrides par comparaison avec des valeurs de seuil respectives. Procédé selon l’une quelconque des revendications 1 à 4 dans lequel, lors du calcul des valeurs de navigation certifiées, on utilise un filtre de Kalman pour le calcul de covariances d’erreur à partir des données de navigation. Procédé selon la revendication 5, dans lequel le filtre de Kalman est un filtre de Kalman invariant. Procédé selon l’une quelconque des revendications 1 à 5, dans lequel on procède à un recalage à vitesse nulle de l’étage de navigation inertielle apte à calculer les valeurs de navigation certifiées. Procédé selon l’une quelconque des revendications 1 à 7, dans lequel les valeurs de navigation hybrides sont calculées à partir des données inertielles issues de capteurs de mesure, notamment d’au moins trois gyroscopes et d’au moins trois accéléromètres, dans trois directions de l’espace, et de données additionnelles, notamment de modèles de déplacement. Centrale de navigation inertielle, comprenant une unité d’acquisition de données de navigation (1) et un premier étage de navigation (2) apte à calculer des valeurs de navigation hybrides par fusion de données de navigation, caractérisée en ce qu’elle comporte en outre un deuxième étage de navigation (3) apte à calculer des valeurs de navigation certifiées, indépendamment des valeurs de navigation hybrides, lesdites valeurs de navigation certifiées étant certifiées avec une limite d’erreur de géolocalisation.