FR 3064350 A1 20180928 FR 1700303 A 20170323 1 Procédé de calcul d'une vitesse d'un aéronef, procédé de calcul d'un rayon de protection, système de positionnement et aéronef associés La présente invention concerne un procédé de calcul d'une vitesse d'un aéronef. La présente invention concerne en particulier un procédé de calcul d'une vitesse intègre protégée contre les possibles défaillances d'un système de positionnement par satellites ou GNSS (de l'anglais « Global Navigation Satellite System »), un système de positionnement et un aéronef associés. Les aéronefs sont fréquemment équipés de systèmes de pilotage sophistiqués 10 dans lesquels de nombreuses fonctions sont contrôlées par des systèmes automatiques ou de systèmes d'affichage synthétiques dans lesquels une représentation tridimensionnelle de l'environnement et des informations sur la vitesse de l'aéronef sont affichées. De tels systèmes supposent une très bonne connaissance de la vitesse et de la 15 position de l'aéronef pour assurer la sécurité de l'aéronef et de ses passagers. En particulier, ces mesures doivent non seulement être précises, mais également intègres. Une mesure est dite « précise » si la différence entre la grandeur mesurée et la grandeur réelle est inférieure à un premier seuil prédéfini. Une mesure est dite « intègre » s'il est possible de calculer une probabilité que 20 l'erreur de mesure soit supérieure à un deuxième seuil déterminé, en prenant en compte les causes de défaillances possibles. Par exemple, une intégrité à le ou moins par heure est fréquemment exigée pour des mesures de vitesse d'un aéronef. Le deuxième seuil est fréquemment appelé « rayon de protection ». Lorsque l'erreur de mesure est un vecteur ayant des composantes dans les axes horizontaux et 25 verticaux, comme par exemple un vecteur d'erreur de position ou un vecteur d'erreur de vitesse, il est possible d'associer à cette mesure un rayon de protection horizontal et un rayon de protection vertical. Lorsque le rayon de protection porte sur la position un volume de protection associé est fréquemment défini par le rayon de protection horizontal et le rayon de protection vertical. 30 Des méthodes de calcul de vitesse basées soit sur l'utilisation des mesures de Doppler de la porteuse soit sur la dérivation de la position obtenue par les mesures de pseudo-distance de code données par un système de positionnement par satellite tel que le système GPS (de l'anglais « Global Positioning System ») sont connues. Cependant, le système GPS (stations au sol, satellites en orbite) fournit un 35 engagement d'intégrité sur une position mesurée, mais pas sur la variation temporelle de cette position, qu'elle soit obtenue de l'une ou l'autre manière. En effet, la notion de 3064350 2 défaillance affectant la position GPS est bien définie et quantifiée, alors que la notion de défaillance affectant la vitesse GPS ne l'est pas. Il est donc possible que l'erreur de position reste comprise dans son volume de protection, mais varie rapidement à l'intérieur de ce volume. Il en découle une grande erreur de vitesse même dans des cas où la position reste comprise dans son volume de protection. On connaît du document US2014074397 Al un procédé de calcul d'une vitesse intègre utilisant un système de positionnement par satellite couplé à des mesures inertielles. Cependant, ce procédé est complexe et demande des ressources importantes en termes de calculs puisque les rayons de protection horizontaux et verticaux sont obtenus par l'utilisation d'un filtre de Kalman. Un tel procédé n'est pas intégrable dans des systèmes avioniques possédant des capacités de calcul limitées. Un des buts de l'invention est de proposer un procédé de calcul d'une vitesse d'un aéronef qui soit intègre et précise, tout en nécessitant moins de puissance de calcul que les procédés de calcul de l'état de la technique. A cet effet, l'invention propose un procédé de calcul d'une vitesse intègre d'un aéronef vis-à-vis des défaillances d'un système de positionnement par satellites, comprenant: - l'obtention d'une vitesse mesurée de l'aéronef à partir d'une unité de mesure inertielle de l'aéronef, - l'obtention d'une position mesurée de l'aéronef, associé à un rayon de protection d'intégrité de position, à partir d'un module de positionnement par satellite de l'aéronef, - le calcul, par une boucle de correction, d'une vitesse corrigée, le calcul de la vitesse corrigée comprenant le calcul d'une position calculée par intégration de la vitesse corrigée, et la correction de la vitesse mesurée en fonction d'une différence entre la position calculée et la position mesurée, - le calcul d'un rayon de protection d'intégrité sur la vitesse corrigée, égal à la somme d'une contribution de la pire variation de l'erreur de position mesurée à l'intérieur d'un volume de protection défini par le rayon de protection d'intégrité de position, et d'une contribution des imperfections de l'unité de mesure inertielle. Le procédé de calcul comprend en option une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou selon toutes les combinaisons techniquement possibles : - le procédé comprend, en outre : - la mesure, par l'unité de mesure inertielle, d'au moins une accélération de l'aéronef dans un référentiel lié à l'aéronef, 3064350 3 - la projection de chaque accélération mesurée dans un référentiel géographique pour obtenir au moins une accélération projetée, la projection étant, en outre, fonction d'une accélération de la pesanteur, et - l'intégration de l'accélération projetée pour obtenir la vitesse mesurée. 5 - le calcul d'une vitesse corrigée comprend : - la correction de l'accélération projetée d'une quantité égale au produit entre un gain, appelé deuxième gain, et l'écart de position, et - la correction de la vitesse mesurée d'une quantité égale au produit entre un premier gain et l'écart de position. 10 - le premier gain et le deuxième gain sont constants. - le rayon de protection d'intégrité de la vitesse corrigée est égal à la somme d'un rapport du double du rayon de protection d'intégrité de la position mesurée sur une constante de temps de la boucle de correction, et d'une première fonction de l'erreur maximale de l'accélération projetée. 15 - la première fonction de l'erreur maximale de l'accélération projetée est égale à la somme d'une deuxième fonction et d'un produit, la deuxième fonction étant une fonction d'une dynamique de la trajectoire de l'aéronef et le produit étant un produit d'une erreur maximale de l'accélération projetée et de la constante de temps. - la constante de temps est égale à la racine carrée du rapport entre, au 20 numérateur, le rayon de protection d'intégrité de la position mesurée et, au dénominateur, l'erreur maximale de l'accélération projetée. L'invention concerne également un système de positionnement comprenant une unité de mesure inertielle, un module de géolocalisation par satellite et un calculateur, le calculateur étant configuré pour mettre en oeuvre un procédé de calcul d'une vitesse 25 intègre tel que précédemment décrit. L'invention concerne également un aéronef équipé d'un système de positionnement tel que précédemment décrit. L'invention et ses avantages seront mieux compris à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en référence aux 30 dessins annexés, sur lesquels : - la figure 1 est une représentation schématique d'un exemple de système de positionnement, - la figure 2 est un ordinogramme des étapes d'un procédé de calcul d'un rayon de protection mis en oeuvre par le système de positionnement de la figure 1, et 35 - la figure 3 est une représentation schématique d'un autre exemple de système de positionnement. 3064350 4 Un premier exemple de système de positionnement 10 d'un aéronef est représenté sur la figure 1. Le système électronique de positionnement 10 est configuré pour estimer au moins une position mesurée Pm et une vitesse de l'aéronef, appelée vitesse corrigée. 5 Dans la suite, la vitesse corrigée sera également appelée « vitesse intègre ». Par exemple, le système de positionnement 10 est propre à estimer la position et la vitesse de l'aéronef dans un référentiel géographique. En particulier, le système de positionnement 10 est propre à estimer un vecteur vitesse V de l'aéronef. Il est entendu par « référentiel géographique local » un référentiel comprenant une 10 origine et trois axes. Les axes du référentiel sont fixés par rapport à des directions géographiques. Par exemple, les directions « Nord-Sud », « Est-Ouest » et la direction verticale, associées à une origine qui est la position courante de l'aéronef forment un référentiel géographique local. De préférence, la position mesurée Pm comprend un vecteur comportant une 15 latitude Lat, une longitude Lon et une altitude Alt. En variante, la position mesurée Pm comprend des coordonnées cartésiennes dans un référentiel ECEF. Un référentiel ECEF (de l'anglais « Earth-centered Earth-fixed ») est un référentiel dont l'origine est le centre de la Terre et les directions sont fixes par rapport à la surface de la Terre. . Il est entendu par « vecteur vitesse » un vecteur comprenant trois coordonnées de 20 vitesses définies chacune selon une direction du référentiel géographique. Par exemple, le vecteur vitesse comprend une vitesse verticale Vv et deux vitesses horizontales Ve et Vn. La vitesse verticale Vv est également appelée « composante verticale ». Un vecteur ayant pour coordonnées les deux vitesses horizontales Ve et Vn est 25 appelé « composante horizontale » de la vitesse. Il est à noter que le vecteur vitesse est susceptible d'être exprimé en coordonnées sphériques. Dans ce cas, le vecteur vitesse comprend une amplitude et deux angles. L'amplitude est une norme de la vitesse dans le référentiel. Chaque angle est un angle entre le vecteur vitesse et une direction du référentiel géographique. Par exemple, un 30 angle est l'angle entre le vecteur vitesse et le nord, également appelé « angle de dérive ». Le passage d'une représentation à l'autre pour le vecteur vitesse est obtenue par une conversion géométrique connue. Le système de positionnement 10 est, en outre, configuré pour calculer au moins un rayon de protection d'intégrité. Un rayon de protection d'intégrité est également connu 35 sous l'acronyme IPL (de l'anglais « Integrity Protection Limit »). 3064350 5 Il est entendu, par « rayon de protection d'intégrité » associé à une grandeur calculée ou mesurée, la grandeur calculée ou mesurée étant une estimation d'une grandeur réelle, une valeur telle que la probabilité pour que la différence, en valeur absolue, entre la grandeur calculée ou mesurée et la grandeur réelle soit supérieure ou 5 égale au rayon de protection d'intégrité, sans qu'une alarme ait été émise, est inférieure ou égale à un seuil prédéterminé. Le seuil prédéterminé est alors appelé « probabilité de non-intégrité » (en anglais « probability of hazardously misleading information »). Le seuil prédéterminé est exprimé, par exemple, en heures-1. 10 Ainsi, alors que la précision borne l'erreur vis-à-vis des événements normaux, ou normaux et rares, le rayon de protection d'intégrité définit une borne supérieure de l'erreur vis-à-vis des événements anormaux tels que des défaillances. Les concepts "normaux", "normaux et rares" et "anormaux", tels qu'utilisés dans cette description, sont définis au paragraphe 2.1.1 de l'annexe R du standard 15 aéronautique RTCA DO-229 défini par la « Radio Technical Commission for Aeronautics ». Le rayon de protection d'intégrité associé à une grandeur est exprimé dans la même unité de mesure que cette grandeur. Par exemple, un rayon de protection d'intégrité associé à une vitesse est exprimé en mètres par seconde. 20 Lorsque la vitesse mesurée est intègre, la probabilité que la différence entre la vitesse réelle et la vitesse mesurée soit supérieure ou égale au rayon de protection d'intégrité associé à la vitesse mesurée, sans qu'une alarme ne soit émise, en prenant en compte toutes les défaillances possibles des constituants du système de positionnement, est inférieure ou égale à un seuil prédéterminé, qui est la probabilité de non-intégrité. La 25 probabilité d'une telle défaillance est en général notée sous la forme « 10-x par heure ». Le système de positionnement 10 est configuré pour calculer au moins un rayon de protection d'intégrité de la vitesse de l'aéronef. Le système de positionnement 10 est, en outre, configuré pour émettre une alarme d'intégrité de vitesse AlarmeV. Plus précisément, le système de positionnement 10 est configuré pour calculer au 30 moins un rayon de protection hVIL de la composante horizontale de la vitesse et un rayon de protection vVIL de la composante verticale de la vitesse. La probabilité que l'erreur de vitesse horizontale soit supérieure au rayon de protection horizontal hVIL de la vitesse, sans que l'alarme d'intégrité de vitesse AlarmeV soit émise, est plus petite que la probabilité de non-intégrité. 3064350 6 La probabilité que l'erreur de vitesse verticale soit supérieure au rayon de protection vertical vVIL de la vitesse, sans que l'alarme d'intégrité de vitesse AlarmeV soit émise, est plus petite que la probabilité de non-intégrité. Le système de positionnement 10 comporte une unité de mesure inertielle 15, un 5 module de géolocalisation par satellite 20 et un calculateur 25. L'unité de mesure inertielle 15 est configurée pour mesurer des paramètres inertiels tels que des accélérations, des angles ou encore des vitesses angulaires. L'unité de mesure inertielle 15 comprend des accéléromètres et des gyromètres. En particulier, l'unité de mesure inertielle 15 est configurée pour mesurer au moins 10 une accélération dans un référentiel associé à l'aéronef. Il est entendu par « référentiel associé à l'aéronef » un référentiel comprenant une direction longitudinale (également appelée « axe de roulis »), une direction transverse (également appelée « axe de tangage ») et une direction normale (également appelée « axe de lacet »). Les axes de roulis, de tangage et de lacet sont des axes communément 15 utilisés dans le domaine de l'aéronautique. En particulier, l'unité de mesure inertielle 15 est configurée pour mesurer trois accélérations selon trois directions du référentiel associé à l'aéronef. Selon l'exemple de la figure 1, l'unité de mesure inertielle 15 est configurée pour mesurer une accélération Al de l'aéronef selon la direction longitudinale, une accélération 20 At selon la première direction transversale et une accélération Ano selon la direction normale. L'unité de mesure inertielle 15 est propre à fonctionner en mode « à composants liés » (également connu sous le nom de mode « strapdown »). Le mode strapdown est un mode dans lequel les composants de l'unité de mesure inertielle 15 sont fixés par rapport 25 à l'aéronef et non montés sur des cardans compensant les rotations de l'aéronef. L'unité de mesure inertielle 15 est également configurée pour mesurer au moins un angle d'orientation par rapport au référentiel géographique. En particulier, l'unité de mesure inertielle 15 est configurée pour mesurer trois angles d'orientation par rapport au référentiel géographique local. 30 Il est entendu par « trois angles d'orientation » les trois angles d'Euler. Les angles d'Euler sont définis comme étant les angles de 3 rotations successives permettant d'aligner le référentiel lié à aéronef avec le référentiel géographique local. Les angles d'Euler comprennent le cap C, le roulis Ro et le tangage Tan. Le cap C est, par exemple, un cap magnétique défini par rapport à la direction du pôle Nord magnétique. 3064350 7 Il est à noter que d'autres représentations de l'orientation, telles que les quaternions cosinus directeurs, sont possibles si ces représentations permettent de calculer une matrice de passage du référentiel lié à l'aéronef au référentiel géographique. Par exemple, l'unité de mesure inertielle 15 comprend un module AHRS (de 5 l'anglais « Altitude and Heading Reference System ») comprenant des accéléromètres, des gyromètres, des magnétomètres, et propre à mesurer des angles de roulis Ro et de tangage Tan et un angle de cap magnétique C. Les modules AHRS sont en général moins performants que les modules IRS (de l'anglais « lnertial Reference System ») pendant les phases de vol dynamiques de l'aéronef. 10 Le module de géolocalisation par satellite 20 est configuré pour mesurer une position Pm de l'aéronef par l'exploitation de signaux radio émis par des satellites d'une constellation de satellites d'un système de positionnement par satellite. Le module de géolocalisation par satellite 20 est également appelé « récepteur de radionavigation », ou « récepteur de positionnement par satellite », « Module GNSS » ou encore « module de 15 positionnement GNSS » (d'après l'acronyme GNSS, pour « Global Navigation Satellite System » en anglais). La position Pm mesurée par le module de géolocalisation par satellite 20 est également appelée « position GNSS ». Par exemple, le module de géolocalisation par satellite 20 est un module de positionnement par satellite du système GPS (de l'anglais « Global Positioning System »). 20 En variante, il s'agit un module de positionnement par satellite d'un système GNSS différent du système GPS, par exemple GALILEO, BEIDOU ou encore GLONASS. Selon un mode de réalisation, le module de géolocalisation par satellite 20 est compatible avec un système d'augmentation par satellite ou SBAS (de l'anglais « Satellite-Based Augmentation System »). Un système d'augmentation par satellite est 25 complémentaire des systèmes GNSS et permet d'augmenter la précision de ceux-ci. Des systèmes SBAS ont été développés par différentes agences gouvernementales sous les appellations « EGNOS » en Europe, « MTSAT » ou « MSAS » au Japon, ou encore « WAAS » aux Etats-Unis d'Amérique. Le module de géolocalisation par satellite 20 est, en outre, propre à calculer au 30 moins un rayon de protection d'intégrité de la position mesurée Pm de l'aéronef, associé à un signal d'alarme sur la position AlarmeP. Par exemple, le module de géolocalisation par satellite 20 est propre à calculer un rayon de protection horizontal d'intégrité hPIL et un rayon de protection vertical d'intégrité vPIL de la position mesurée Pm. Le module de géolocalisation par satellite 20 est configuré pour émettre le signal 35 d'alarme sur la position AlarmeP lorsque les rayons de protection hPIL et vPIL de la position ne bornent plus l'erreur de mesure de la position mesurée Pm. 3064350 8 Des techniques d'intégrité connues prennent en compte les défaillances possibles des différents types de systèmes GNSS considéré, par exemple les techniques ABAS (de l'anglais « Aircraft Based Augmentation Systems ») de type RAIM (de l'anglais « Receiver autonomous integrity monitoring »), ou les techniques SBAS. 5 Les rayons de protection horizontal hPIL et vertical vPIL de la position mesurée Pm sont calculés par le module de géolocalisation par satellite 20 en prenant en compte les défaillances du système GNSS. La probabilité que l'erreur de position GNSS horizontale soit supérieure au rayon de protection horizontal hPIL de la position, sans que l'alarme d'intégrité de position 10 AlarmeP ne soit émise, est inférieure à la probabilité de non intégrité de la position. La probabilité que l'erreur de position GNSS verticale soit supérieure au rayon de protection vertical vPIL de la position, sans que l'alarme d'intégrité de position AlarmeP ne soit émise, est inférieure à la probabilité de non intégrité de la position. En particulier, les rayons de protection horizontaux hPIL et verticaux vPIL de la 15 position sont calculés en prenant en compte les événements normaux, normaux et rares, et anormaux pouvant affecter le système GNSS considéré et le module de géolocalisation par satellite 20. Le calculateur 25 est configuré pour estimer au moins une vitesse Vv, Vn, Ve de l'aéronef. 20 Le calculateur 25 est, en outre, configuré pour calculer au moins un rayon de protection d'intégrité hVIL, vVIL associé à une vitesse de l'aéronef, et pour commander l'émission du signal d'alarme d'intégrité de vitesse AlarmeV. De préférence, le calculateur 25 est configuré pour estimer un rayon de protection d'intégrité hVIL pour la composante horizontale de la vitesse de l'aéronef et un rayon de 25 protection d'intégrité vVIL pour la composante verticale de la vitesse de l'aéronef. Le calculateur 25 comporte un module de correction 30, un module 35 de projection 35, un premier module de calcul 40, un deuxième module de calcul 45, une boucle de correction 50, également appelée correcteur ou lisseur, et un troisième module de calcul 55. 30 Le calculateur 25 est, par exemple, un ordinateur comportant une mémoire et un processeur propre à exécuter des instructions de programme d'un produit programme d'ordinateur stocké dans la mémoire. Les modules de correction 30, de projection 35, les premier, deuxième et troisième modules de calcul 40, 45 et 55 et la boucle de correction 50 sont prévus sous la forme d"instruction de programme exécutables par le processeur 35 du calculateur. 3064350 9 En variante, au moins un des modules de correction 30, de projection 35, les premier, deuxième et troisième modules de calcul 40, 45 et 55 et/ou la boucle de correction 50 sont réalisés sous la forme de circuits intégrés dédiés et/ou de circuits logiques programmables. 5 Le module de correction 30 est configuré pour calculer un cap géographique, également appelé « cap corrigé » Cc, à partir du cap magnétique mesuré par l'unité de mesure inertiellel5. Par exemple, le cap corrigé Cc est l'angle entre l'axe de roulis projeté sur l'horizontale et la direction du Nord géographique projetée sur l'horizontale. Par exemple, le module de correction 30 comprend des tables de déclinaison 10 magnétique et est propre à calculer le cap corrigé Cc, selon un algorithme d'interpolation, à partir des tables de déclinaison magnétique. Le module de projection 35 est configuré pour calculer au moins une accélération de l'aéronef dans le référentiel géographique. Par exemple, le module de projection 35 est configuré pour calculer une accélération An de l'aéronef selon la direction nord-sud, une 15 accélération Ae de l'aéronef selon la direction est-ouest et une accélération Av de l'aéronef selon la direction verticale à partir des accélérations Al, At, Ano, des angles de roulis Ro et de tangage Tan, et de l'angle de cap géographique Cc. Le premier module de calcul 40 est configuré pour calculer une accélération g de la pesanteur. Par exemple, le premier module de calcul 40 est propre à calculer 20 l'accélération g de la pesanteur à partir de la latitude Lat et de l'altitude Alt fournies par le module de géolocalisation par satellite 20. Le deuxième module de calcul 45 est configuré pour calculer une accélération Ac de Coriolis à partir de la latitude Lat et de l'altitude Alt fournies par le module de géolocalisation par satellite 20 et de la vitesse Vn, Vv, Ve de l'aéronef. 25 La boucle de correction 50 est configurée pour calculer, à partir d'au moins une accélération An, Ae, Av dans le référentiel géographique, une vitesse Vn, Vv, Ve de l'aéronef. Selon l'exemple de la figure 1, la boucle de correction 50 est configuré pour calculer trois composantes du vecteur vitesse V de l'aéronef. 30 La boucle de correction 50 comporte un premier comparateur 60, un premier intégrateur 65, un deuxième comparateur 70, un deuxième intégrateur 75, un troisième comparateur 80 et un premier multiplicateur 85. Le troisième module de calcul 55 est configuré pour calculer au moins un rayon de protection d'intégrité hVIL, vVIL d'une vitesse Vn, Vv, Ve de l'aéronef. Par exemple, le 35 troisième module de calcul 55 est configuré pour calculer un rayon de protection 3064350 10 d'intégrité hVIL de la composante horizontale de la vitesse de l'aéronef et un rayon de protection d'intégrité vVIL de la vitesse verticale de l'aéronef. Le fonctionnement du système de positionnement va maintenant être décrit en référence à la figure 2, qui représente un ordinogramme des étapes d'un procédé de 5 calcul d'un rayon de protection d'intégrité hVIL, vVIL d'une vitesse de l'aéronef. Le procédé de calcul d'une vitesse intègre comprend une phase 100 de mise en oeuvre d'un procédé de calcul d'une vitesse intègre de l'aéronef et une étape 110 de calcul du rayon de protection hVIL, vVIL. La phase de mise en oeuvre 100 d'un procédé de calcul d'une vitesse de l'aéronef 10 comprend une étape 120 de mesure d'une accélération, une première étape 130 d'intégration, une première étape 140 de correction, une deuxième étape 150 d'intégration et une première étape 160 d'estimation. Au cours de l'étape de mesure 120, l'accélération Al selon la direction longitudinale, l'accélération At selon la direction transversale et l'accélération Ano selon la 15 direction normale sont mesurées par l'unité de mesure inertielle 15. Le tangage Tan, le roulis Ro et le cap C sont également mesurés par l'unité de mesure inertielle 15. Le cap corrigé Cc est calculé à partir du cap magnétique C par le module de correction 30. Les accélérations mesurées Al, At, Ano sont projetées, par le module de projection 20 35, dans le référentiel géographique pour obtenir au moins une accélération Ae, An, Av projetée. Selon l'exemple de la figure 1, l'accélération projetée Ae selon la direction nord-sud, l'accélération projetée Ae selon la direction est-ouest et l'accélération projetée Av selon la direction verticale sont calculées. Il est à noter qu'une erreur maximale iacc est définie qui borne l'accélération 25 projetée. En particulier, une erreur maximale acch 5 est définie pour les accélérations Ae, An - selon des directions horizontales et une erreur maximale6accv est définie pour 1"accélération Av selon la direction verticale. Ainsi, l'erreur maximale 5 -acch borne l'erreur d'accélération projetée sur l'horizontale 30 et l'erreur maximale 6accv borne l'erreur d'accélération projetée sur la verticale. L'erreur maximale 6acc sur l'accélération projetée traduit en particulier la contribution des erreurs de mesure de l'accélération mesurée par l'unité de mesure inertielle 15 et des erreurs de projection. En outre, la position mesurée Pm de l'aéronef est mesurée par le module de 35 géolocalisation par satellite 20. 3064350 11 Lors de la première étape d'intégration 130, au moins une accélération Ae, An, Av dans le référentiel géographique est intégrée par le premier intégrateur 65 pour obtenir une vitesse mesurée V1, V2, V3 de l'aéronef selon une direction du référentiel géographique. 5 Selon l'exemple de la figure 1, au cours de la première étape d'intégration 130, la vitesse V1 mesurée selon la direction nord-sud, la vitesse V2 mesurée selon la direction est-ouest et la vitesse V3 mesurée selon la direction verticale sont calculées. Par exemple, les accélérations Al, At1, At2 mesurées par l'unité de mesure inertielle 15 étant définies dans le référentiel lié à l'aéronef, les vitesses mesurées V1, 10 V2, V3 sont obtenues chacune à partir de l'accélération Ae, An, Av projetée correspondante, de l'accélération de la pesanteur g calculée par le premier module de calcul 40 et de l'accélération de Coriolis Ac calculée par le deuxième module de calcul 45. Par exemple, chaque accélération projetée Ae, An, Av est corrigée, par le premier 15 comparateur 60, en fonction de l'accélération de Coriolis Ac et de l'accélération de la pesanteur g calculées. Chaque vitesse mesurée V1, V2, V3 est alors obtenue par intégration, après correction, de l'accélération projetée Ae, An, Av correspondante. Lors de la première étape 140 de correction, au moins une vitesse corrigée Vv, 20 Vn, Ve est calculée. Plus précisément, une vitesse corrigée Vn selon la direction nord- sud, une vitesse corrigée Ve selon la direction est-ouest et une vitesse corrigée Vv selon la direction verticale sont calculées. Chaque vitesse corrigée Vv, Vn, Ve est obtenue par correction de la vitesse mesurée V1, V2, V3 correspondante par un correcteur du premier ordre. Le correcteur du 25 premier ordre est ici formé par le deuxième comparateur 70, le deuxième intégrateur 75, le troisième comparateur 80 et le premier multiplicateur 85. En particulier, chaque vitesse corrigée Vv, Vn, Ve obtenue par l'ajout à la vitesse mesurée correspondante, par le deuxième comparateur 70, d'une première quantité Q1. Lors de la deuxième étape d'intégration 150, chaque vitesse corrigée Vv, Vn, Ve 30 est intégrée par le deuxième intégrateur 75 pour obtenir une position calculée Pc. A l'étape 160 d'estimation, un écart E entre la position Pm mesurée par le module de positionnement par satellite 20 et la position Pc calculée est estimée par le troisième comparateur 80. La première quantité Q1 est également estimée par le premier multiplicateur 85. 35 La première quantité Q1 est égale au produit entre l'écart E et un premier gain K1. Ainsi, la vitesse mesurée V1, V2, V3 est corrigée en fonction de l'écart E. 3064350 12 Le premier gain K1 est constant. Une constante de temps Tau est définie pour le correcteur du premier ordre de la boucle de correction 50. La constante de temps Tau est égale à l'inverse du premier gain K1. Cela s'écrit mathématiquement : 5 Tau =11 K1 (Equation 1) L'étape 160 d'estimation est suivie de l'étape de calcul 110 du rayon de protection d'intégrité vVIL, hVIL. L'étape 120 de mesure d'une accélération, la première étape 130 d'intégration, la première étape 140 de correction, la deuxième étape 150 d'intégration et la première 10 étape 160 d'estimation sont répétées dans cet ordre avec une période temporelle P. La période temporelle est, par exemple, égale à 20 Hertz (Hz). Au cours de l'étape de calcul 110, le module de géolocalisation par satellite 20 calcule au moins un rayon de protection d'intégrité hPIL, vPIL de la position. Plus précisément, le module de géolocalisation par satellite 20 calcule un rayon de protection 15 d'intégrité horizontal hPIL et un rayon de protection d'intégrité vertical vPIL. Chaque rayon de protection d'intégrité hVIL, vVIL associé à la vitesse corrigée Vn, Ve, Vv est alors calculé en fonction d'un rayon de protection d'intégrité hPIL, vPIL de la position Pm fourni par le module de géolocalisation par satellite 20 et de l'erreur maximale d'accélération M -acc de l'unité de mesure inertielle 15. 20 Chaque rayon de protection d'intégrité hVIL, vVIL sur la vitesse est calculé par un calcul analytique. Il est entendu par « calcul analytique » que chaque rayon de protection d'intégrité hVIL, vVIL sur la vitesse est calculé selon une formule mathématique dépendant de valeurs disponibles à un instant donné. En particulier, le résultat d'un calcul analytique 25 n'est pas obtenu en intégrant en fonction du temps des grandeurs ou à partir d'une estimation en fonction de valeurs passées. Chaque rayon de protection d'intégrité hVIL, vVIL sur la vitesse corrigée est égal à la somme de la contribution de la pire variation de l'erreur de position Pm à l'intérieur du volume de protection défini par les rayons de protection sur la position hPIL et vPIL, et de 30 la contribution des imperfections de l'unité de mesure inertielle 15. La contribution de la pire variation de l'erreur de position, lorsque cette variation est propagée au travers de la boucle de correction 50, de constante de temps Tau, est égale à un rapport, appelé premier rapport, d'un rayon de protection d'intégrité hPIL, vPIL de la position mesurée Pm, multiplié par 2, sur la constante de temps Tau. 35 La contribution des imperfections de l'unité de mesure inertielle est représentée par une fonction, appelée première fonction Vil-inertiel. 3064350 13 Ainsi, chaque rayon de protection d'intégrité de la vitesse hVIL, vVIL est conforme à l'équation : VIL = 2 x PIL + VIL ine'id (équation 2) Tau dans laquelle VIL est le rayon de protection hVIL, vVIL sur la vitesse calculé, PIL 5 est le rayon de protection hPIL, vPIL sur la position correspondant et Tau est la constante de temps de la boucle de correction 50. Afin de prendre en compte la contribution inertielle due à l'erreur maximale d'accélération projetée et intégrée sur la constante de temps Tau lors d'une phase de vol stabilisé, la première fonction VILinertiel est conforme à l'équation : 10 VIL inertiel = acc X Tau (équation 3) dans laquelle 5 - acc est l'erreur maximale sur la mesure d'accélération projetée qui traduit la contribution des erreurs de mesure de l'accélération par l'unité de mesure inertielle 15 et des erreurs de projection lors d'une phase de vol stabilisé. Il est entendu par « phase de vol stabilisé » une période pendant laquelle l'aéronef 15 est quasiment en translation à vitesse constante. La fonction 5 -acc est définie différemment pour le calcul d'un rayon de protection hVIL portant sur la vitesse horizontale et pour le calcul d'un rayon de protection vVIL portant sur la vitesse verticale. L'erreur horizontale d'accélération projetée 5 est calculée en prenant en bacon 20 compte la contribution des accéléromètres bacch et l'erreur d'attitude 3.(ph. L'erreur d'attitude Acph est une borne supérieure de l'erreur de roulis et de tangage lors d'une phase de vol stabilisé. La contribution des accéléromètres bacch est une borne supérieure de l'erreur d'accélération horizontale mesurée par les accéléromètres. L'accélération mesurée par les accéléromètres en phase de vol stabilisé étant 25 proche de la réaction -g de l'aéronef porté par la masse d'air à la pesanteur, l'erreur horizontale d'accélération projetée 5 est, par exemple, conforme à l'équation : batch batch = I batch+ g XIAqj (équation 4) dans laquelle « » désigne une valeur absolue. L'erreur verticale d'accélération projetée 5 - accv est calculée en prenant en compte la 30 contribution des accéléromètres battu et l'influence de l'erreur d'altitude sur l'estimation de la pesanteur. Par exemple, l'erreur verticale d'accélération est conforme à l'équation : 5accv =lbaccv1+ kv x vPIL (équation 5) dans laquelle k, est un coefficient de sensibilité de la pesanteur terrestre à l'altitude. Le coefficient de sensibilité de la pesanteur terrestre à l'altitude 1 REVENDICATIONS1.- Procédé de calcul d'une vitesse intègre d'un aéronef vis-à-vis des défaillances d'un système de positionnement par satellites, comprenant: - l'obtention (120) d'une vitesse mesurée (V1, V2, V3) de l'aéronef à partir d'une unité de mesure inertielle (15) de l'aéronef, - l'obtention d'une position mesurée (Pm) de l'aéronef, associée à un rayon de protection d'intégrité de position (hPIL, vPIL), à partir d'un module de positionnement par satellite (20) de l'aéronef, - le calcul, par une boucle de correction (50), d'une vitesse corrigée (Vn, Ve, Vv), le calcul de la vitesse (Vn, Ve, Vv) corrigée comprenant le calcul (150) d'une position calculée (Pc) par intégration de la vitesse corrigée (Vn, Vv, Ve), et la correction (140) de la vitesse mesurée (V1, V2, V3) en fonction d'une différence (c) entre la position calculée (Pc) et la position mesurée (Pm), - le calcul d'un rayon de protection d'intégrité (vVIL, hVIL) sur la vitesse corrigée (Vn, Ve, Vv), égal à la somme d'une contribution de la pire variation de l'erreur de position mesurée (Pm) à l'intérieur d'un volume de protection défini par le rayon de protection d'intégrité de position (hPIL, vPIL), et d'une contribution des imperfections de l'unité de mesure inertielle (15). 2.- Procédé de calcul d'une vitesse intègre selon la revendication 1, comprenant, en outre : - la mesure, par l'unité de mesure inertielle (15), d'au moins une accélération (AI, Ano, At) de l'aéronef dans un référentiel lié à l'aéronef, - la projection de chaque accélération (Al, Ano, At) mesurée dans un référentiel géographique pour obtenir au moins une accélération projetée (An, Ae, Av), la projection étant, en outre, fonction d'une accélération de la pesanteur (g), et - l'intégration de l'accélération projetée (An, Ae, Av) pour obtenir la vitesse mesurée (V1, V2, V3). 3.- Procédé de calcul d'une vitesse intègre d'un aéronef selon la revendication 2, dans lequel le calcul d'une vitesse corrigée comprend : - la correction de l'accélération projetée d'une quantité égale au produit entre un gain, appelé deuxième gain (K2), et l'écart de position (e), et - la correction de la vitesse mesurée (V1, V2, V3) d'une quantité égale au produit entre un premier gain (K1) et l'écart de position (s). 3064350 20 4.- Procédé de calcul d'une vitesse intègre d'un aéronef selon la revendication 3, dans lequel le premier gain (K1) et le deuxième gain (K2) sont constants. 5 5.- Procédé de calcul d'une vitesse intègre d'un aéronef selon l'une quelconque des revendications 2 à 4, dans lequel le rayon de protection d'intégrité (hVIL, vVIL) de la vitesse corrigée (Vn, Vc, Ve) est égal à la somme d'un rapport du double du rayon de protection d'intégrité (hPIL, vPIL) de la position mesurée (Pm) sur une constante de temps (Tau) de la boucle de correction (50), et d'une première fonction (VIL 1 de merttel, 10 l'erreur maximale de l'accélération projetée (An, Ae, Av). 6.- Procédé de calcul d'une vitesse intègre d'un aéronef selon la revendication 5, dans lequel la première fonction (VlI-inertiel) de l'erreur maximale de l'accélération projetée (An, Ae, Av) est égale à la somme d'une deuxième fonction (Kdyn) et d'un produit, la 15 deuxième fonction (Kdyn) étant une fonction d'une dynamique de la trajectoire de l'aéronef et le produit étant un produit d'une erreur maximale (6accv, 6acch) de l'accélération projetée et de la constante de temps (Tau). 7.- Procédé de calcul d'une vitesse intègre d'un aéronef selon la revendication 6, 20 dans lequel la constante de temps (Tau) est égale à la racine carrée du rapport entre, au numérateur, le rayon de protection d'intégrité (hPIL, vPIL) de la position mesurée (Pm) et, au dénominateur, l'erreur maximale (6accv, bacch) de l'accélération projetée. 8.- Système de positionnement (10) comprenant une unité de mesure inertielle 25 (15), un module de géolocalisation par satellite (20) et un calculateur, le calculateur étant configuré pour mettre en oeuvre un procédé de calcul d'une vitesse intègre selon l'une quelconque des revendications 1 à 7. 9.- Aéronef équipé d'un système de positionnement (10) selon la revendication 8. 30