Brevet FR2988835A1 — publié le 20131004

Procédé de détermination d'un état de crédibilité de mesures de capteurs d'un aéronef et système correspondant La présente invention concerne un procédé de détermination d'un état de crédibilité de mesures de capteurs d'un aéronef lors d'un vol dudit aéronef, lesdits capteurs comprenant des capteurs de pression statique, de pression totale et d'incidence. Le pilotage d'un aéronef repose sur la connaissance de paramètres de vol de celui-ci, telles que sa vitesse par rapport à l'air ambiant, son altitude et son incidence. Ces paramètres sont déterminés au moyen de sondes situées sur la carlingue de l'aéronef. De manière connue, ces sondes comprennent des capteurs de pression statique, des sondes de Pitot pour la mesure de la pression totale, des sondes d'incidence montées sur un dispositif à palette ou pneumatique, et des sondes de température totale. Ces sondes sont alors raccordées à des moyens de détermination de la grandeur correspondante. Notamment, un anémomètre détermine la vitesse de l'aéronef par rapport à l'air à partir des mesures de pressions totale et statique, et un altimètre détermine l'altitude de l'aéronef à partir de mesures de la pression statique. Ces mesures sont ensuite regroupées et affichées sur un dispositif d'affichage qui constitue une source d'information centrale à partir de laquelle le pilotage de l'aéronef est effectué. De manière connue, les sondes d'incidence et de pression se présentent sous la forme de palettes et de tubes protubérants de la peau de l'aéronef. Elles sont ainsi exposées à des facteurs météorologiques ou mécaniques qui peuvent en altérer le fonctionnement, en particulier par obstruction d'orifices de ces sondes par du givre ou des poussières ou insectes, ou par blocage des dispositifs à palette. De telles défaillances conduisent à la production de mesures erronées, et notamment à l'affichage d'incidences, de vitesses et/ou d'altitudes fausses, pouvant conduire le pilote à effectuer des manoeuvres inappropriées. Par exemple, de fausses informations de pilotage peuvent conduire au décrochage de l'aéronef ou à la perte de son contrôle du fait d'une vitesse trop élevée. Afin de minimiser les conséquences de telles défaillances, les régulations en aéronautique imposent aux constructeurs de disposer de moyens de mesures de ces caractéristiques redondants. Ainsi, les aéronefs comprennent généralement au moins une sonde de secours identique à chaque sonde susceptible de faillir. Toutefois, cette solution ne donne pas entière satisfaction. En effet, les sondes de secours existantes sont pour la plupart de type protubérant, et par conséquent présentent les mêmes risques de défaillance que les sondes qu'elles sont destinées à suppléer. Ainsi, en cas de défaillance des sondes de mesure, aucune mesure fiable n'est fournie au pilote. En outre, le pilote ne dispose d'aucun moyen pour vérifier la fiabilité des informations fournies par les sondes. Un but de l'invention est donc de détecter d'éventuels dysfonctionnements de capteurs de l'aéronef pour être à même d'avertir le pilote de ces dysfonctionnements et de lui fournir des informations alternatives fiables. A cette fin, l'invention a pour objet un procédé du type précité, caractérisé en ce qu'il comprend les étapes suivantes : - déterminer une information de vitesse de l'aéronef, à partir de mesures de pressions statique Ps et totale PT desdits capteurs de pression statique et de pression totale, - déterminer un premier coefficient de portance Ez de l'aéronef à partir d'une mesure d'incidence a dudit capteur d'incidence et de ladite information de vitesse, - déterminer une masse estimée (h de l'aéronef, - déterminer si une équation de sustentation de l'aéronef est vérifiée, à partir de ladite masse estimée C-1-7 , d'une information représentative du facteur de charge de l'aéronef, de ladite information de vitesse et dudit premier coefficient de portance de l'aéronef, avec une première tolérance prédéterminée, - activer un état de crédibilité optimale, dans lequel les mesures desdits capteurs de pression statique, de pression totale et d'incidence sont jugées fiables, si ladite équation de sustentation est vérifiée avec ladite première tolérance, - activer au moins un état de crédibilité non optimale, dans lequel les mesures d'au moins un capteur parmi lesdits capteurs de pression statique, de pression totale et d'incidence sont jugées non fiables, si ladite équation de sustentation n'est pas vérifiée avec ladite première tolérance. Suivant des modes particuliers de réalisation, le procédé de détermination comporte l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - ladite information de vitesse est un nombre de Mach Ma, et l'étape consistant à déterminer si l'équation de sustentation de l'aéronef est vérifiée comprend la vérification de l'égalité 0,7SPsMa2Ez = nziiig avec la première tolérance déterminée, où S représente une surface de référence de l'aéronef, ledit premier coefficient de portance az est une projection selon un axe Z choisi, avantageusement un axe Z perpendiculaire à l'axe longitudinal de l'aéronef 1 défini notamment comme l'axe moyen de poussée des moteurs, d'un coefficient de portance de l'aéronef selon un axe orthogonal au vecteur vitesse de l'aéronef et d'un coefficient de traînée parallèle à ce vecteur vitesse, nz désigne ladite information représentative du facteur de charge de l'aéronef selon l'axe Z, et g l'accélération de la pesanteur ; - ledit premier coefficient de portance est déterminé à partir d'une table de correspondance fournissant une estimation de la valeur dudit premier coefficient de portance à partir de ladite information de vitesse, de la mesure d'incidence a et d'une configuration de vol de l'aéronef ; - la détermination de la masse estimée de l'aéronef comprend les étapes suivantes : - détermination d'une première masse de l'aéronef, à partir de caractéristiques de l'aéronef déterminées avant ou pendant le décollage de l'aéronef, - détermination d'une deuxième masse de l'aéronef lors du vol de l'aéronef, à partir d'une équation de sustentation de l'aéronef exprimant la masse de l'aéronef en fonction d'une information représentative du facteur de charge de l'aéronef, d'un deuxième coefficient de portance de l'aéronef, d'une information de vitesse de l'aéronef et d'une pression statique d'une masse d'air traversée par l'aéronef, déterminées à partir de mesures de capteurs de l'aéronef lors du vol de l'aéronef, - évaluation de la masse estimée, à partir d'une masse calculée tenant compte desdites première et deuxième masses, par exemple à partir de la moyenne pondérée desdites première et deuxième masses ; - lesdites première et deuxième masses sont des masses de l'aéronef à un instant initial au parking, et l'étape d'évaluation de ladite masse estimée comprend la détermination de ladite masse calculée à partir desdites première et deuxième masses, la détermination d'une masse de carburant consommé depuis ledit instant initial, ladite masse estimée étant évaluée en soustrayant ladite masse de carburant consommé de ladite masse calculée ; - le procédé comprend une étape de détermination d'une altitude auxiliaire de l'aéronef à partir d'au moins une mesure d'altitude géographique obtenue indépendamment d'un capteur de pression, notamment issue d'un capteur de positionnement par satellite, et une étape de mise à disposition d'un équipage de l'aéronef de ladite altitude auxiliaire, au moins lorsque l'état de crédibilité non optimale est activé ; - ladite altitude auxiliaire est déterminée à partir d'une altitude ZGPS déterminée indépendamment d'un capteur de pression de l'aéronef, sous la forme : Zp** (t) = ZGPS (t)- AZ** (j) où AZ**(j) est un terme correctif ajusté périodiquement, destiné à corriger un décalage entre l'atmosphère standard et l'atmosphère réelle de vol de l'aéronef, ledit terme correctif prenant une valeur figée lorsque l'état de crédibilité non optimale est activé ; - le procédé comprend, au moins lorsque l'état de crédibilité non optimale est activé, une étape de détermination d'au moins une information de vitesse auxiliaire à partir d'au moins une mesure d'un capteur de position géographique de l'aéronef, et une étape de mise à disposition d'un équipage de l'aéronef de ladite information de vitesse auxiliaire ; - l'état de crédibilité non optimale est choisi parmi : - un état de crédibilité forte, dans lequel les mesures d'au moins un capteur parmi les capteurs de pressions statique et de pression totale sont jugées non fiables et les mesures dudit capteur d'incidence sont jugées fiables, et - un état de crédibilité faible dans lequel les mesures du capteur d'incidence sont jugées non fiables ; - au moins dans l'état de crédibilité forte, le procédé comporte une étape de détermination d'une information de vitesse auxiliaire de crédibilité forte, sur la base de la mesure d'incidence du capteur d'incidence, en fonction d'un estimateur de pression statique déterminé à partir d'au moins une mesure obtenue indépendamment d'un capteur de pression de l'aéronef ; - ladite information de vitesse auxiliaire de crédibilité forte est déterminée périodiquement à chaque instant à partir d'une information de vitesse auxiliaire de crédibilité forte obtenue à un instant précédent, avantageusement à l'aide d'une série par récurrence ; - le procédé comprend une estimation de la vitesse du vent, ladite vitesse du vent étant déterminée périodiquement à partir d'informations de vitesses issues d'un capteur de position géographique ; - au moins dans l'état de crédibilité faible, le procédé comporte une étape de détermination d'une information de vitesse auxiliaire de crédibilité fiable, comprenant les étapes suivantes : - détermination d'informations de vitesse de l'aéronef par rapport au sol à partir de mesures d'un capteur de position géographique de l'aéronef, - estimation d'une vitesse du vent, - détermination de ladite vitesse auxiliaire de crédibilité faible par rapport à l'air, à partir desdites informations de vitesse (de l'aéronef par rapport au sol et de ladite estimation de la vitesse du vent ; - le procédé comprend une étape de détermination d'un état de fiabilité de mesures dudit capteur d'incidence, comprenant les étapes suivantes : - déterminer un estimateur de pression statique à partir d'au moins une mesure obtenue indépendamment d'un capteur de pression de l'aéronef, - déterminer une information de vitesse auxiliaire, sur la base dudit estimateur de pression statique, - déterminer un deuxième coefficient de portance de l'aéronef, à partir d'une mesure d'incidence a dudit capteur d'incidence et de ladite information de vitesse auxiliaire, - déterminer si l'équation de sustentation de l'aéronef est vérifiée, à partir de ladite masse estimée rn , de ladite information de vitesse auxiliaire et dudit deuxième coefficient de portance de l'aéronef, avec une deuxième tolérance prédéterminée, - confirmer l'état de état de crédibilité faible, dans lequel les mesures du capteur d'incidence sont jugées non fiables, si ladite équation de sustentation n'est pas vérifiée avec ladite deuxième tolérance prédéterminée, - activer l'état de crédibilité forte, dans lequel les mesures d'au moins un capteur parmi les capteurs de pression statique et de pression totale sont jugées non fiables et les mesures dudit capteur d'incidence sont jugées fiables, si ladite équation de sustentation est vérifiée avec ladite tolérance prédéterminée ; - le procédé comprend, dans l'état de crédibilité non optimale, une étape d'affichage d'un message d'alerte à destination d'un équipage de l'aéronef, indiquant que les mesures d'au moins un capteur choisi parmi lesdits capteurs de pression statique, de pression totale et d'incidence ne sont pas fiables. L'invention a également pour objet un système de détermination d'un état de crédibilité de mesures de capteurs d'un aéronef lors d'un vol dudit aéronef, lesdits capteurs comprenant des capteurs de pression statique, de pression totale et d'incidence, caractérisé en ce qu'il comprend : - des moyens configurés pour déterminer une information de vitesse de l'aéronef, à partir de mesures de pressions statique Ps et totale PT desdits capteurs de pression statique et de pression totale, - des moyens configurés pour déterminer un premier coefficient de portance de l'aéronef à partir d'une mesure d'incidence a dudit capteur d'incidence et de ladite information de vitesse, - des moyens configurés pour déterminer une masse estimée de l'aéronef, - des moyens configurés pour déterminer si une équation de sustentation de l'aéronef est vérifiée, à partir de ladite masse estimée fil , d'une information représentative du facteur de charge, de ladite information de vitesse et dudit coefficient de portance de l'aéronef, avec une première tolérance prédéterminée, - des moyens configurés pour activer un état de crédibilité optimale, dans lequel les mesures desdits capteurs de pression statique, de pression totale et d'incidence sont jugées fiables, si ladite équation de sustentation est vérifiée avec ladite première tolérance, et - des moyens configurés pour activer au moins un état de crédibilité non optimale, dans lequel les mesures d'au moins un capteur parmi lesdits capteurs de pression statique, de pression totale et d'incidence sont jugées non fiables, si ladite équation de sustentation n'est pas vérifiée avec ladite première tolérance. L'invention a également pour objet un procédé de détermination d'un état de crédibilité de mesures d'au moins un capteur d'incidence d'un aéronef lors d'un vol dudit aéronef, caractérisé en ce qu'il comprend les étapes suivantes : - déterminer un estimateur de pression statique Ps** d'une masse d'air traversée par l'aéronef à partir d'au moins une mesure d'altitude obtenue indépendamment d'un capteur de pression, notamment issue d'un capteur de positionnement par satellites, - déterminer une information de vitesse auxiliaire M** de l'aéronef, sur la base dudit estimateur de pression statique Ps**, - déterminer un premier coefficient de portance az** de l'aéronef à partir d'une mesure d'incidence a dudit capteur d'incidence et de ladite information de vitesse M**, - déterminer une masse estimée m de l'aéronef, - déterminer si une équation de sustentation de l'aéronef est vérifiée, à partir d'une information représentative du facteur de charge nz, de ladite masse estimée , de ladite information de vitesse M** et dudit coefficient de portance Ez** de l'aéronef, avec une première tolérance prédéterminée, - activer au moins un état de crédibilité fiable, dans lequel les mesures dudit capteur d'incidence sont jugées fiables, si ladite première équation de sustentation est vérifiée avec ladite première tolérance, - activer un état de crédibilité faible, dans lequel les mesures dudit capteur d'incidence sont jugées non fiables, si ladite première équation de sustentation n'est pas vérifiée avec ladite première tolérance. Suivant des modes particuliers de réalisation, le procédé de détermination d'un état de crédibilité de mesures d'au moins un capteur d'incidence comporte l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - ladite information de vitesse auxiliaire est un nombre de Mach M**, et l'étape consistant à déterminer si l'équation de sustentation de l'aéronef est vérifiée comprend la vérification de l'égalité 0,7SPs M**262 = nZing , où S représente une surface de référence de l'aéronef, ledit coefficient de portance est une projection selon un axe Z choisi, avantageusement un axe Z perpendiculaire à l'axe longitudinal de l'aéronef défini comme un axe moyen de poussée des moteurs, d'un coefficient de portance de l'aéronef selon un axe orthogonal au vecteur vitesse de l'aéronef et d'un coefficient de traînée parallèle à ce vecteur vitesse, nz désigne ladite information représentative du facteur de charge de l'aéronef selon l'axe Z, et g l'accélération de la pesanteur, - ledit coefficient de portance az** est déterminé à partir d'une table de correspondance fournissant une estimation de la valeur dudit coefficient de portance Cz** à partir de ladite information de vitesse auxiliaire, de la mesure d'incidence a et d'une configuration de vol de l'aéronef, - le procédé comprend une estimation de la vitesse du vent, ladite vitesse du vent étant déterminée périodiquement à partir d'informations de vitesses indépendamment d'un capteur de pression lorsque l'état de crédibilité fiable est activé, et ladite vitesse du vent ayant une valeur figée lorsque l'état de crédibilité faible est activé, - l'étape de détermination d'au moins une information de vitesse auxiliaire de l'aéronef comprend les étapes suivantes : - détermination d'informations de vitesse de l'aéronef par rapport au sol à partir de mesures d'altitude obtenues indépendamment d'un capteur de pression, - estimation d'une vitesse du vent, - détermination de ladite information de vitesse de l'aéronef, par rapport à l'air, à partir desdites informations de vitesse de l'aéronef par rapport au sol et de ladite estimation de la vitesse du vent ; - le procédé comprend, au moins lorsque l'état de crédibilité faible est activé, une étape de mise à disposition d'un équipage de l'aéronef de ladite information de vitesse auxiliaire, - le procédé comprend, au moins lorsque l'état de crédibilité faible est activé, une étape de détermination d'une altitude auxiliaire de l'aéronef à partir d'au moins une mesure d'altitude obtenue indépendamment d'un capteur de pression, et une étape de mise à disposition d'un équipage de l'aéronef de ladite altitude auxiliaire, au moins lorsqu'un état de crédibilité non optimale est activé, - ladite altitude auxiliaire Zp** est déterminée à partir d'une altitude ZGF,s( ) déterminée par un capteur d'altitude de l'aéronef, sous la forme : Zp**(t). Gps (t)- AZ- (j) où AZ (j) est un terme correctif ajusté périodiquement, destiné à corriger un décalage entre l'atmosphère standard et l'atmosphère réelle de vol de l'aéronef, ledit terme correctif prenant une valeur figée lorsqu'un état de crédibilité non optimale est activé, - le procédé comprend, dans l'état de crédibilité faible, une étape d'affichage d'un message d'alerte à destination d'un équipage de l'aéronef, indiquant que les mesures dudit capteur d'incidence ne sont pas fiables, - l'état de crédibilité fiable est choisi parmi : - un état de crédibilité optimale, dans lequel les mesures desdits capteurs de pression statique, de pression totale et d'incidence sont jugées fiables, - un état de crédibilité forte, dans lequel les mesures d'au moins un capteur parmi les capteurs de pression statique et de pression totale sont jugées non fiables et les mesures dudit capteur d'incidence sont jugées fiables ; - le procédé comprend une étape de détermination d'un état de crédibilité de mesures desdits capteurs de pression statique et de pression totale, comprenant les étapes suivantes : - déterminer une information de vitesse de l'aéronef, à partir de mesures de pressions statique Ps et totale PT desdits capteurs de pression statique et de pression totale, - déterminer un deuxième coefficient de portance de l'aéronef à partir d'une mesure d'incidence a dudit capteur d'incidence et de ladite information de vitesse, - déterminer si une deuxième équation de sustentation de l'aéronef est vérifiée, à partir de ladite masse estimée m, de ladite deuxième information de vitesse et dudit deuxième coefficient de portance de l'aéronef, avec une deuxième tolérance prédéterminée, - activer l'état de crédibilité optimale si lesdites première et deuxième équations de sustentation sont vérifiées, avec lesdites première et deuxième tolérances prédéterminées respectivement, - activer l'état de crédibilité forte, si ladite première équation de sustentation est vérifiée avec ladite première tolérance et ladite deuxième équation de sustentation n'est pas vérifiée avec ladite deuxième tolérance ; - dans l'état de crédibilité forte, le procédé comporte une étape de détermination d'une information de vitesse auxiliaire de crédibilité forte, sur la base de la mesure d'incidence du capteur d'incidence, en fonction d'une valeur de pression statique déterminée à partir d'au moins une mesure obtenue indépendamment d'un capteur de pression, notamment à partir d'un capteur de positionnement par satellites de l'aéronef, - l'information de vitesse auxiliaire de crédibilité forte est déterminée périodiquement à chaque instant à partir d'une information de vitesse auxiliaire de crédibilité forte obtenue à un instant précédent, avantageusement à l'aide d'une série par récurrence. L'invention a également pour objet un système de détermination d'un état de crédibilité de mesures d'au moins un capteur d'incidence d'un aéronef lors d'un vol dudit aéronef, caractérisé en ce qu'il comprend : - des moyens configurés pour déterminer un estimateur de pression statique d'une masse d'air traversée par l'aéronef à partir d'au moins une mesure d'altitude obtenue indépendamment d'un capteur de pression, - des moyens configurés pour déterminer une information de vitesse auxiliaire de l'aéronef, sur la base dudit estimateur de pression statique, - des moyens configurés pour déterminer un premier coefficient de portance de l'aéronef à partir d'une mesure d'incidence dudit capteur d'incidence et de ladite information de vitesse, - des moyens configurés pour déterminer une masse estimée d7 de l'aéronef, - des moyens configurés pour déterminer si une équation de sustentation de l'aéronef est vérifiée, à partir d'une information représentative du facteur de charge, de ladite masse estimée, de ladite information de vitesse et dudit coefficient de portance de l'aéronef, avec une première tolérance prédéterminée, - des moyens configurés pour activer au moins un état de crédibilité fiable, dans lequel les mesures dudit capteur d'incidence sont jugées fiables, si ladite première équation de sustentation est vérifiée avec ladite première tolérance, - des moyens configurés pour activer un état de crédibilité faible, dans lequel les mesures dudit capteur d'incidence sont jugées non fiables, si ladite première équation de sustentation n'est pas vérifiée avec ladite première tolérance.35 L'invention a également pour objet un procédé de détermination d'un estimateur d'altitude d'un aéronef en vol à un instant de détermination, à partir d'au moins une information d'altitude instantanée obtenue indépendamment d'un capteur de pression audit instant de détermination, caractérisé en ce qu'il comprend les étapes suivantes : - détermination, lors du vol de l'aéronef, d'au moins un terme correctif d'altitude, à partir d'au moins un écart d'altitude entre une information d'altitude antérieure obtenue indépendamment d'un capteur de pression de l'aéronef, déterminée en au moins un instant antérieur audit instant de détermination et une altitude pression issue d'au moins une mesure de pression statique d'un capteur de pression statique de l'aéronef, déterminée audit instant antérieur, - détermination de ladite information d'altitude instantanée indépendamment d'un capteur de pression audit instant de détermination, - calcul dudit estimateur d'altitude par correction de ladite information d'altitude instantanée avec ledit terme correctif d'altitude. Suivant des modes particuliers de réalisation, le procédé de détermination d'un estimateur d'altitude comporte l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - le procédé comprend la détermination périodique, en une pluralité d'instants antérieurs audit instant de détermination, d'un écart d'altitude entre une information d'altitude antérieure obtenue indépendamment d'un capteur de pression et une altitude pression issue d'au moins une mesure de pression statique d'un capteur de pression statique de l'aéronef ; - ledit terme correctif d'altitude est égal à l'écart d'altitude déterminé à l'instant antérieur le plus proche dudit instant de détermination ; - le procédé comprend la détermination en chaque instant antérieur d'un état de fiabilité de la valeur de l'écart d'altitude déterminée audit antérieur, comprenant l'estimation de la fiabilité de mesures dudit capteur de pression statique ; - la détermination d'un état de fiabilité de la valeur de l'écart d'altitude déterminée à un instant antérieur comprend la détermination d'une différence entre ledit écart d'altitude et l'écart d'altitude déterminé à l'instant antérieur immédiatement précédent, la valeur de l'écart d'altitude étant jugée non fiable si ladite différence est supérieure à un seuil prédéterminé ; - la détermination d'un état de fiabilité de la valeur de l'écart d'altitude comprend l'actionnement par un pilote d'une commande pour indiquer que ledit écart d'altitude est fiable ; - ledit terme correctif d'altitude est égal à l'écart d'altitude jugé fiable déterminé à l'instant antérieur le plus proche dudit instant de détermination ; - le procédé comprend une étape de détermination d'un estimateur d'altitude auxiliaire obtenu au moyen d'une table d'atmosphère standard à partir d'une pression QNH local et d'une information d'altitude instantanée obtenue indépendamment d'un capteur de pression de l'aéronef, indiquant une altitude de l'aéronef au-dessus de l'isobare du QNH ; - ledit procédé comprend la détermination, à partir dudit estimateur d'altitude, d'un estimateur de pression statique d'une masse d'air traversée par l'aéronef, - ledit estimateur de pression statique est déterminé au moyen d'une table d'atmosphère standard. L'invention a également pour objet un système de détermination d'un estimateur d'altitude d'un aéronef en vol à un instant de détermination, à partir d'au moins une information d'altitude instantanée obtenue indépendamment d'un capteur de pression audit instant de détermination, caractérisé en ce qu'il comprend les étapes suivantes : - des moyens configurés pour déterminer, lors du vol de l'aéronef, au moins un terme correctif d'altitude, à partir d'au moins un écart d'altitude entre une information d'altitude antérieure obtenue indépendamment d'un capteur de pression et une altitude pression issue d'au moins une mesure de pression statique d'un capteur de pression statique de l'aéronef, déterminé en au moins un instant antérieur audit instant de détermination, - des moyens configurés pour déterminer ladite information d'altitude instantanée indépendamment d'un capteur de pression audit instant de détermination, - des moyens configurés pour calculer ledit estimateur d'altitude par correction de ladite information d'altitude instantanée avec ledit terme correctif d'altitude. L'invention a également pour objet un procédé de détermination de caractéristiques d'une masse d'air traversée par un aéronef lors d'un vol dudit aéronef, caractérisé en ce qu'il comprend une estimation d'une vitesse horizontale du vent autour dudit aéronef comprenant : - la détermination d'un estimateur de pression statique de ladite masse d'air à partir d'au moins une mesure d'altitude obtenue indépendamment d'un capteur de pression, - l'estimation d'une vitesse vraie de l'aéronef par rapport à ladite masse d'air sur la base dudit estimateur de pression statique, - la détermination d'une vitesse horizontale et d'une vitesse verticale de l'aéronef par rapport au sol à partir d'informations de position géographique issues d'un capteur de position géographique de l'aéronef, - la détermination de ladite vitesse horizontale du vent, à partir de ladite vitesse horizontale et de ladite vitesse verticale de l'aéronef par rapport au sol, et de ladite vitesse vraie. Suivant des modes particuliers de réalisation, le procédé de détermination de caractéristiques de ladite masse d'air comporte l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement 10 possibles : - la détermination de l'estimateur de pression statique comprend la détermination d'une altitude estimée de l'aéronef à partir de ladite mesure d'altitude, ledit estimateur de pression statique étant déterminé à partir de ladite altitude estimée au moyen d'une table d'atmosphère standard ; 15 - la détermination de ladite vitesse vraie de l'aéronef par rapport à ladite masse d'air comprend : - la détermination d'un nombre de Mach de l'aéronef, sur la base dudit estimateur de pression statique et d'une mesure d'incidence de l'aéronef par au moins un capteur d'incidence de l'aéronef, 20 - l'estimation d'une température statique de ladite masse d'air à partir d'une mesure d'une température totale de ladite masse d'air par au moins un capteur de température totale dudit aéronef, et dudit nombre de Mach de l'aéronef, - la détermination de ladite vitesse vraie de l'aéronef par rapport ladite masse d'air à partir dudit nombre de Mach et de ladite estimation de température 25 statique ; - ledit nombre de Mach est déterminé périodiquement à chaque instant à partir d'un nombre de Mach obtenu à un instant précédent, avantageusement à l'aide d'une série par récurrence ; - le procédé comprend la détermination répétée d'un estimateur d'écart de 30 température standard, calculé comme un écart entre ladite estimation de température statique et une température statique attendue à ladite altitude estimée de l'aéronef ; - le procédé comprend la détermination d'un état de fiabilité de mesures dudit capteur d'incidence, l'activation d'au moins un état de crédibilité fiable si les mesures dudit capteur d'incidence sont jugées fiables, et l'activation d'un état de crédibilité faible si 35 les mesures sont jugées non fiables ; - dans l'état de crédibilité faible, la détermination de ladite vitesse vraie de l'aéronef par rapport à ladite masse d'air comprend : - l'estimation d'une incidence auxiliaire de l'aéronef, indépendamment d'un capteur d'incidence, à partir d'une équation de sustentation de l'aéronef, - l'estimation de ladite vitesse vraie comme un quotient entre une accélération de l'aéronef et une vitesse angulaire de cet aéronef, sur la base de ladite incidence auxiliaire et de sa dérivée ; - lorsque l'état de crédibilité faible est activé, la valeur dudit estimateur d'écart de température standard est figée à la dernière valeur de cet estimateur déterminée dans l'état de crédibilité fiable ; - le procédé comprend, dans l'état de crédibilité faible, une étape d'estimation d'une température statique auxiliaire comme une somme de ladite température statique attendue à ladite altitude estimée de l'aéronef et de ladite valeur figée dudit estimateur d'écart de température standard. L'invention a également pour objet un système de détermination de caractéristiques d'une masse d'air traversée par un aéronef lors d'un vol dudit aéronef, caractérisé en ce qu'il comprend des moyens configurés pour estimer une vitesse horizontale du vent autour dudit aéronef comprenant : - des moyens pour déterminer un estimateur de pression statique de ladite masse d'air à partir d'au moins une mesure d'altitude obtenue indépendamment d'un capteur de pression, - des moyens pour estimer une vitesse vraie de l'aéronef par rapport à ladite masse d'air sur la base dudit estimateur de pression statique, - des moyens pour déterminer une vitesse horizontale et une vitesse verticale de l'aéronef par rapport au sol à partir d'informations de position géographique issues d'un capteur de position géographique de l'aéronef, - des moyens pour déterminer ladite vitesse horizontale du vent, à partir de ladite vitesse horizontale et de ladite vitesse verticale de l'aéronef par rapport au sol, et de ladite vitesse vraie. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins suivants, sur lesquels : - la figure 1 illustre de manière schématique un aéronef dans des conditions de vol classiques, auquel est appliquée l'invention ; - la figure 2 est un schéma synoptique général représentant le système selon l'invention ; - la figure 3 est un schéma synoptique illustrant le procédé selon l'invention ; - la figure 4 est un schéma synoptique illustrant les étapes mises en oeuvre pour évaluer la fiabilité de mesures d'incidence selon un mode de réalisation ; - la figure 5 est un schéma synoptique représentant des étapes mises en oeuvre pour évaluer la fiabilité de mesures de pression selon un mode de réalisation ; - la figure 6 illustre un mode de représentation par un dispositif d'affichage auxiliaire d'un premier ensemble d'informations relatives au vol ; - la figure 7 illustre un mode de représentation par un dispositif d'affichage auxiliaire d'un deuxième ensemble d'informations relatives au vol ; - la figure 8 illustre un mode de représentation par un dispositif d'affichage auxiliaire d'un troisième ensemble d'informations relatives au vol. La Figure 1 représente, de manière schématique, un aéronef 1 en vol auquel est appliqué le procédé selon l'invention. L'aéronef 1 est représenté sur la Figure 1 par son seul centre de gravité. Son axe longitudinal est orienté selon un axe X, qui forme avec l'horizontale A un angle 0 appelé assiette de l'aéronef. Il se déplace par rapport à l'air selon un vecteur vitesse Vair , qui forme avec l'horizontale A un angle 7,,, appelé pente air de l'aéronef. L'angle a entre l'axe longitudinal X de l'aéronef 1 et son vecteur vitesse est appelé angle d'incidence. Ces angles vérifient ainsi la relation : 0 = + yak. La vitesse propre Vp de l'aéronef, qui est sa vitesse par rapport à l'air dans un plan horizontal, est liée à sa vitesse Vs par rapport au sol dans ce plan horizontal par le triangle des vitesses selon la relation : Vp = Vs - W , où W désigne le vecteur vitesse du vent dans le plan horizontal. Le système 2 de détermination de paramètres de l'aéronef 1 selon l'invention est représenté de manière schématique sur la Figure 2. Le système 2 est propre à déterminer des caractéristiques de l'aéronef 1 ou de l'air ambiant lors du vol, et à évaluer la fiabilité des mesures de capteurs de l'aéronef 1. Ce système 2 comprend un calculateur central 3 muni d'une mémoire 4. La mémoire 4 comprend notamment des tables de correspondance et des abaques fournissant la valeur d'une caractéristique de vol de l'aéronef en fonction de la valeur d'une ou plusieurs autres caractéristiques. En particulier, la mémoire 4 comprend une table d'atmosphère type donnant l'altitude pression de l'aéronef, notée Zp, en fonction de la pression statique autour de l'aéronef. Cette altitude pression est définie comme l'altitude qu'aurait un aéronef dans une atmosphère standard s'il était à cette pression statique. La mémoire 4 comprend également une table de correspondance fournissant la valeur d'un coefficient de portance de l'aéronef à partir de valeurs du nombre de Mach M et d'incidence a de l'aéronef. Cette table de correspondance est issue d'essais en vol préalablement réalisés. La mémoire 4 comprend par ailleurs des abaques, telles qu'un abaque permettant de déterminer la masse de l'aéronef 1 en fonction de la portance et de la traînée pendant la phase initiale de décollage, et de l'accélération en résultant. Ce système 2 comprend par ailleurs plusieurs capteurs 5, et notamment un capteur de pression statique 5a, un capteur de pression totale 5b, un capteur de température 5c et un capteur d'incidence 5d. Le capteur de pression statique 5a est propre à mesurer la pression statique Ps, c'est-à-dire la pression atmosphérique au niveau de l'aéronef. Le capteur de pression totale 5b est par exemple une sonde Pitot. Il est propre à mesurer la pression totale PT, somme de la pression dynamique Pdy,, et de la pression statique Ps. Le capteur de température 5c est propre à mesurer la température totale, notée TAT pour « Total Air Temperature », correspondant à la température d'impact de l'air dans la sonde destinée à la mesure. Cette température est constante sur tout le filet d'air qui arrive sur la sonde. La température totale est liée à la température statique de l'air, notée SAT pour « Static Air Temperature », et correspondant à la température que mesurerait un thermomètre à l'arrêt dans la masse d'air, par la relation : SAT = TAT TAT (1 4. y-1 m2\ (1+ 0.2M2) (1) 2 où M désigne le nombre de Mach de l'aéronef et y 1.4 est le coefficient adiabatique de l'air. Le capteur d'incidence 5d est propre à déterminer l'incidence a de l'aéronef 1. Tous les capteurs sont reliés au calculateur 3. Le système 2 de détermination comprend également un altimètre 6, un anémomètre 7, un machmètre 8 et un accéléromètre 9. L'altimètre 6 est relié au capteur de pression statique 5a et au calculateur 3. Il est propre à déterminer l'altitude pression Zp de l'aéronef 1 par rapport à un niveau de référence à partir de la mesure de la pression statique Ps. A cette fin, l'altimètre 6 utilise par exemple une table d'atmosphère type comprenant des valeurs tabulées de pressions statiques en fonction de l'altitude. L'anémomètre 7 est relié aux capteurs de pression statique 5a et totale 5b et au calculateur 3. Il est propre à déterminer, à partir des pressions statique Ps et totale PT, la pression dynamique Pdyn = PT Ps et à déduire de cette pression dynamique la vitesse conventionnelle Vc de l'aéronef par rapport à l'air, à partir de la relation : 101325 2 -3,5 PdYn = 1+ V 661,471 -1 (2) où 101325 correspond à la pression atmosphérique au sol, 661,471 est la vitesse du son au sol en noeuds. La vitesse conventionnelle est ainsi donnée par 1 ( P jird +1 -1 101d3e 25 0,2 Vc = 661,471 (3) Cette vitesse conventionnelles Vc est celle qui aurait produit la même pression dynamique Paya en volant dans l'atmosphère standard au sol. L'équivalent de vitesse EV peut être déduit de l'altitude pression Zp et du nombre de Mach M par : -1 1 poEV2 =- 2 pa2 M (4) 2 2 où po désigne la masse volumique de l'air au sol, a la vitesse du son et p la masse volumique de l'air. Lorsque le nombre de Mach de l'aéronef est faible (M Le machmètre 8 est relié aux capteurs de pression statique 5a et totale 5b et au calculateur 3. Il est propre à déduire de la pression totale PT et de la pression statique Ps le nombre de Mach de l'aéronef. On notera ainsi Ma le nombre de Mach tel que déterminé par le machmètre 8. L'accéléromètre 9 comprend une centrale à inertie. Il est propre à déterminer un vecteur d'accélération:ide l'aéronef, et notamment ses composantes Jx et Jz suivant l'axe longitudinal X de l'aéronef et l'axe de lacet Z de l'aéronef 1 respectivement. Le système 2 de détermination comprend par ailleurs un capteur de position géographique, avantageusement un capteur d'altitude, tel qu'un capteur de position par satellites, par exemple un capteur GPS 10. Ce capteur GPS 10 est propre à déterminer la position de l'aéronef 1, notamment son altitude exprimée de manière classique au-dessus du géoïde de référence WGS 84 (pour World Geodetic System), appelée altitude GPS ZGps. A partir de cette position, le système 2 est propre à estimer la vitesse horizontale GS de l'aéronef 1 par rapport au sol, la vitesse Vz(GPS) de l'aéronef 1 selon l'axe vertical B et une altitude pression reconstituée, notée Zp**. Le capteur GPS 10 est relié au calculateur 3. Le calculateur 3 est propre à déterminer, à partir d'une valeur de Mach M de l'aéronef, sa vitesse par rapport à l'air, appelée vitesse vraie et notée TAS pour « True Air Speed », suivant la relation : TAS = A/1\lyR -SA T (5) où y est le coefficient adiabatique de l'air et R la constante universelle des gaz parfaits. Le calculateur 3 est également propre à déterminer la norme de la vitesse propre V de l'aéronef, norme de la projection horizontale de la vitesse vraie TAS, à partir de la relation : TAS 2 = vr2 +VZ (GPS) (6) Le système 2 de détermination comprend par ailleurs des moyens 11 de détermination de la quantité de carburant contenue dans les réservoirs de l'aéronef 1. Les moyens 11 sont reliés au calculateur 3. Ces moyens comprennent par exemple des jauges de réservoirs et des débitmètres. Les jauges de réservoirs sont propres à mesurer la quantité, notamment la masse de carburant contenue dans chaque réservoir, permettant au calculateur 3 de déterminer la masse de carburant dans l'aéronef 1. Les débitmètres sont propres à mesurer le débit massique de carburant fourni à chaque moteur, permettant au calculateur 3 de déduire la masse FU de carburant consommée, et à partir d'une mesure initiale de la masse de carburant, la masse de carburant restant. L'aéronef 1 comprend des moyens 12 d'interface homme-machine. Ces moyens 12 sont reliés au système 2 de détermination. Ils sont propres à présenter des informations à destination de l'équipage et notamment du pilote, et à recevoir des instructions du pilote destinées en particulier au système 2 de détermination. L'aéronef 1 comprend ainsi des instruments de bord, propres à présenter au pilote des informations relatives au vol de l'aéronef 1. Notamment, ces instruments comprennent des moyens d'affichage classiques, propres à afficher les caractéristiques du vol issues des mesures de pression et d'incidence par les capteurs de pression et d'incidence. L'aéronef 1 comprend également un dispositif d'affichage auxiliaire 14. Le dispositif d'affichage auxiliaire 14 comprend des moyens 15 d'affichage sélectif de différentes évaluations de caractéristiques du vol, en fonction d'un état de crédibilité affecté aux mesures des capteurs, et notamment des capteurs d'incidence et de pression. Le dispositif d'affichage auxiliaire 14 est complété en outre par des moyens 16 d'affichage de messages d'alerte destinés à informer le pilote de l'état de crédibilité affecté aux mesures des capteurs, et en particulier à l'avertir lorsque les mesures d'un ou plusieurs capteurs ne sont pas fiables. Ces messages peuvent être des signaux visuels et/ou sonores, par exemple des signaux lumineux. Avantageusement, ces signaux sont des messages textuels ou symboliques et sont inclus dans le dispositif d'affichage des pannes de l'avionique résidente. L'aéronef 1 comprend par ailleurs une interface 17 de saisie, par exemple des boutons de commande et un clavier, permettant au pilote de donner des instructions au système de détermination ou d'entrer des valeurs numériques de paramètres de vol. On a représenté sur la Figure 3 un exemple de mise en oeuvre du procédé selon un mode de réalisation de l'invention, pour la surveillance de la crédibilité d'informations de pilotage fournies par les sondes et appareils de mesure de cet aéronef, et en particulier des valeurs de pression et d'incidence mesurées. Le procédé comprend une étape 21 de détermination d'une masse au parking initiale mpj0 de l'aéronef 1, en vue de déterminer une masse estimée de l'aéronef 1 à tout instant lors de son vol. Cette masse au parking initiale mp,0 est élaborée en estimant, selon différentes méthodes, plusieurs masses mp, au parking de l'aéronef 1 et en déterminant la masse au parking initiale mp,o de l'aéronef 1 à partir des estimations mp, ainsi obtenues. L'étape 21 comprend ainsi avantageusement la détermination de trois masses au parking mpl, mp2, mp3, et le choix comme valeur de masse au parking initiale mp,0 de la valeur médiane de ces trois masses (Io= 1, 2 ou 3). L'étape 21 comprend une phase 23 de détermination de la première masse au parking mpi de l'aéronef 1, à un instant initial t1 avant le décollage de l'aéronef 1. Cette masse au parking est déterminée par estimation et sommation de la masse de base de l'aéronef 1, la masse de carburant et la masse du chargement. Cette masse n'est généralement pas exacte. Notamment, la mesure de la masse de carburant contenue dans les réservoirs n'est pas une mesure exacte. On estime l'erreur faite sur cette masse mpi à quelques pourcents. En outre, la masse de base et la masse du chargement peuvent subir des erreurs dues à la mesure, à la transmission d'informations, au pilote et/ou à la saisie dans le système. L'étape 21 comprend une phase 25 de détermination de la deuxième masse au parking mp2 de l'aéronef 1, réalisée lors du décollage, à un instant t2 choisi par exemple 3 secondes après la fin du signal d'enfoncement des trains "Weight off wheels". Cette deuxième masse au parking mp2 est déterminée à partir de caractéristiques de l'aéronef à l'instant t2, notamment la poussée nominale des moteurs au décollage dans des conditions de tarage, de température et de pression saisies par le pilote, et la traînée théorique de l'aéronef dans la configuration "train sorti + volets décollage" dans des conditions mesurées de vitesse et d'incidence de l'aéronef. La phase 25 comprend la détermination d'une masse m2 de l'aéronef 1 à l'instant t2 lors du décollage, à laquelle est ajoutée la masse de carburant FU(t2) consommée depuis le décollage. La phase 25 comprend à cette fin la détermination des forces de poussée F et de traînée T exercées sur l'aéronef 1 à cet instant t2. La phase 25 comprend en particulier la mesure de la température totale TAT par le capteur de température 5c et la détermination de la température statique SAT par le calculateur 3 à partir de cette mesure. Elle comprend en outre la mesure de l'incidence a par le capteur d'incidence 5d, la mesure de la pression statique Ps par le capteur de pression statique 5a, de l'accélération longitudinale Jx et selon l'axe Z Jz de l'aéronef 1 par l'accéléromètre 9, et d'une information de vitesse de l'aéronef. Cette information de vitesse est avantageusement un nombre de Mach indépendant des pressions mesurées et de l'incidence a, noté M**. Ce nombre de Mach M** est déterminé à partir du triangle des vitesses, selon la relation : GS+1/z(GPS)k -1/17(t21 .\1,,R - SAT(t2) dans laquelle : - GS et Vz(GPS) sont respectivement la vitesse horizontale de l'aéronef 1 par rapport au sol et la vitesse de l'aéronef 1 selon l'axe vertical B, dérivée de l'altitude pression GPS par rapport au temps, déterminées à partir de mesure du capteurs GPS 10 à l'instant t2; - k est un vecteur unitaire ascendant parallèle à l'axe vertical B ; - W(t2) est un vecteur vitesse du vent au décollage saisi par le pilote ; et SAT(t2)est la température statique au décollage. A cet instant t2, les roues de l'aéronef 1 ne sont plus en contact avec le sol, de telle sorte que la force de traînée T est déterminée à partir des valeurs du nombre de Mach M** et d'incidence a déterminées. (t2 (7) La force de traînée T est par exemple déterminée à partir de l'équation de propulsion de l'aéronef donnée par : T = 0,7SPsM**2 ax(ce,M* *,conf) (8) où désigne une estimation du coefficient de traînée de l'aéronef, fonction de l'incidence a, du nombre de Mach M** et de la configuration de vol de l'aéronef 1, notée « conf », relative notamment aux volets sortis (ici configuration premier segment). Le coefficient de traînée ax est par exemple estimé à partir d'une table de correspondance stockée dans la mémoire 4 du calculateur 3, fournissant une estimation de la valeur de ax à partir du nombre de Mach M** et de l'incidence a de l'aéronef 1. Cependant, ce coefficient de traînée dépend faiblement du nombre de Mach, et il peut être avantageux, pour simplifier le procédé, de déterminer ax en fonction de l'incidence a et de la configuration seulement, voire de fixer la valeur de ax. La force de poussée F, due aux moteurs, est une fonction connue de la température et de la pression. La masse m2 de l'aéronef 1 au décollage est liée à l'accélération Jx, à la traînée T et à la force de poussée F par la relation : F +T = x x m2 où Fx et TX sont les projections de la force de poussée F et de la traînée T respectivement sur l'axe X. La masse m2 est par exemple déterminée à partir de l'accélération, de la traînée et de la force de poussée, au moyen d'un abaque préalablement établi et stocké dans la mémoire 4. L'abaque est établi en déterminant par pesée de manière précise la masse exacte de l'aéronef pour une série de vols d'essais, puis en mesurant pour chaque vol d'essai l'accélération Jx correspondante, avec différentes configurations de l'aéronef. Cette accélération est mesurée avec un accéléromètre. La deuxième masse au parking mp2 est déterminée à partir de la masse au décollage m2 et d'une mesure par les débitmètres de la masse de carburant consommée entre les instants t1 et t2. La deuxième masse au parking mp2 est ainsi égale à : mp2 = m2 (t2 )+ FU(t2) (10) où FU(t2)désigne la masse de carburant consommé entre les instants t1 et t2 (9) L'étape 21 comprend en outre une phase de détermination 26 de la troisième masse au parking mp3 de l'aéronef 1, également réalisée lors du décollage. La détermination 26 de la troisième masse au parking mp3 comprend la détermination d'une masse m3 de l'aéronef 1 à l'instant t2 lors du décollage, à laquelle est ajoutée la masse de carburant F02) consommée depuis le décollage. La masse m3 est calculée selon une équation de sustentation de l'aéronef mettant en correspondance sa masse effective à l'instant t2, le facteur de charge nz, l'incidence mesurée a, le nombre de Mach estimé M**, et la pression statique mesurée Ps. L'équation de sustentation de l'aéronef est donnée de manière générale par : nzm g =0,7SPsM26z(a,M,conf) (11) où m est la masse de l'aéronef, S représente une surface de référence de l'aéronef, nz le facteur de charge de l'aéronef 1 selon l'axe Z, M le nombre de Mach de l'aéronef et az un coefficient de portance de l'aéronef 1 selon l'axe Z. La phase 26 comprend ainsi la détermination du facteur de charge nz, de la pression statique Ps, du nombre de Mach et du coefficient de portance de l'aéronef 1 à l'instant t2. Le facteur de charge nz à l'instant t2 est alors déterminé par le calculateur 3 à partir de l'expression : nz(t2)= jz(t2) (12) où Jz est l'accélération selon l'axe Z déterminée par l'accéléromètre 9 à l'instant t2. Le coefficient de portance az (t2) estimé à l'instant t2 est une projection selon l'axe Z perpendiculaire à l'axe longitudinal de l'aéronef 1 du coefficient de portance Cz selon un axe orthogonal au vecteur vitesse de l'aéronef 1, et du coefficient de traînée Cx parallèlement à ce vecteur vitesse. Le coefficient de portance 6z vérifie donc : 6z = Cx sin a+Cz cos a. Il comprend en général majoritairement un terme de sustentation. Le coefficient de portance az est fonction de l'incidence a, du nombre de Mach M et de la configuration de vol de l'aéronef 1, relative notamment aux volets sortis. A l'instant t2 considéré, il s'agit de la configuration premier segment (train sorti et braquage de volets choisi pour le décollage). Le coefficient de portance az est par exemple estimé à partir d'une table de correspondance stockée dans la mémoire 4 du calculateur 3, fournissant une estimation de la valeur de az à partir du nombre de Mach et de l'incidence a de l'aéronef 1. Cette table de correspondance est préalablement déterminée en effectuant une série de vols d'essais pour un modèle d'aéronef donné. Différentes configurations d'incidence, de facteur de masse, de charge et de Mach, et différentes configurations d'aéronefs sont balayées pour déterminer dans chaque cas le coefficient C. . Le coefficient de portance est ainsi estimé au moyen de la table de correspondance à partir du nombre de Mach M**(t2) déterminé à l'instant t2 et d'une valeur d'incidence a(t2) mesurée par le capteur d'incidence à cet instant t2, dans la configuration premier segment. Il est ainsi noté az- (t2) : az (t2) = az(a(t2),M (t2),1er segment). La masse m3 de l'aéronef 1 à l'instant t2 est ainsi déduite de l'équation de sustentation selon l'expression : m3 (t2 077SPsA4**2 (t2 )6; (t2 ) (13) nz La masse au parking mp3 est déduite par ajout de la masse de carburant consommé FU(t2) entre l'instant t1 et l'instant t2, déterminée à partir de mesures des débitmètres. mp3 = m3 (t2)+ FU(t2) (14) Les trois masses au parking mn mp2 et mp3 ainsi déterminées sont généralement différentes. L'étape 21 comprend ainsi une phase 27 de détermination d'une valeur de masse au parking initiale de l'aéronef 1, notée mp,o, à partir des trois masses au parking mo, mp2 et mp3 préalablement déterminées. Cette masse au parking mpo est avantageusement égale à la médiane des masses mpi, mp2 et mp3. Ainsi, si une de ces masses a une valeur aberrante, cette valeur aberrante est exclue et n'intervient pas dans la valeur de la masse au parking mp,o. Lors du vol de l'aéronef 1, des paramètres de vol sont déterminés. Ces paramètres comprennent des caractéristiques de l'air ambiant et du vol de l'aéronef 1, et sont déterminés de manière continue ou périodique à des instant t, du vol, lors d'une étape 28. Cette étape 28 comprend une phase de détermination 29 de caractéristiques de l'air ambiant, notamment de la pression statique Ps et des températures statique SAT et totale TAT. L'étape 28 comprend également une phase de détermination 31 de caractéristiques de vol de l'aéronef 1, et en particulier de son incidence a par le capteur d'incidence 5d, de l'altitude pression Zp par l'altimètre 6, du nombre de Mach Ma par le machmètre 8 et de l'accélération J par l'accéléromètre 9, à l'instant t, considéré. Le nombre de Mach Ma constitue une deuxième information de vitesse, en plus de l'équivalent de vitesse EV. L'étape 28 comprend par ailleurs une phase de détermination 33 du coefficient de portance az et du facteur de charge nz de l'aéronef 1 selon l'axe Z à partir des valeurs du nombre de Mach Ma, d'incidence a et d'accélération J déterminées lors de la phase 31. Comme précédemment décrit, le facteur de charge nz est déterminé à partir de l'expression : nz (15) où Jz est l'accélération selon l'axe Z déterminée lors de la phase 31. Par ailleurs, le coefficient de portance az est estimé à partir de la table de correspondance stockée dans la mémoire 4 du calculateur 3, fournissant une estimation de la valeur de az à partir des valeurs du nombre de Mach Ma et d'incidence a déterminées lors de la phase 31. L'étape 28 comprend en outre une phase de détermination 39 d'une masse estimée m de l'aéronef 1 lors de son vol. La détermination 39 de la masse estimée in comprend une phase de détermination 41 d'une quatrième masse au parking de l'aéronef 1, lors du vol de l'aéronef 1, à partir de l'équation de sustentation de l'aéronef. Cette quatrième masse est ainsi appelée masse de sustentation msp. Cette masse de sustentation est déterminée en estimant la masse instantanée de l'aéronef 1 à plusieurs instants t distincts lors du vol, avantageusement en continu, en déduisant à partir de chacune des masses instantanées ainsi évaluées une masse au parking de l'aéronef 1, et en effectuant une moyenne pondérée de ces masses au parking. La masse instantanée de l'aéronef 1 à chaque instant t, notée m4(t), est déduite de l'équation de sustentation à partir de la pression statique Ps mesurée lors de la phase 29, du nombre de Mach Ma déterminé lors de la phase 31, et du coefficient de portance az et du facteur de charge nz déterminés lors de la phase 33, selon l'expression : 0,7SPsMa 2 az(a,Ma,conf) ma (t)= nzg (16) Une masse au parking m4' est déduite par ajout de la masse de carburant consommée FU(t) entre l'instant t1 et l'instant t, déterminée à partir de mesures des débitmètres. Cependant, cette valeur de masse au parking m4' dépend de mesures qui pourraient être erronées, notamment en cas de turbulence, ou de défaillances d'instruments de mesure. On détermine ainsi lors de la phase 41 la masse de sustentation msp comme une moyenne pondérée des masses au parking m4' déterminées entre l'instant t2 de décollage et l'instant t,, cette moyenne étant réalisée seulement sur les instants t lors desquels les mesures sont considérées comme fiables. La masse de sustentation msp est ainsi déterminée à un instant tv à partir de l'expression : fzm (t)(m4(t)+ FU(t)) msp (t' )= t2 (17) prm (t) 12 Xm (t) désigne ici une fonction caractéristique de masse prenant la valeur 1 ou 0 selon des critères géométriques et dynamiques indépendants visant à caractériser le vol rectiligne horizontal de l'aéronef 1. La valeur de la fonction xm(t) à chaque instant t est déterminée par le calculateur 3. Notamment, lorsque l'aéronef 1 est dans une phase de vol stabilisée et calme, c'est-à-dire lorsque la vitesse verticale Vz(GPS) issue du capteur GPS 10 est inférieure à un seuil prédéterminé, la valeur absolue de l'inclinaison (issue des centrales à inertie) est inférieure à 5°, le facteur de charge n z est voisin de 1 et une énergie du facteur de charge, définie comme sa variance, est proche de zéro, xn,(t) prend la valeur 1. Au contraire, lorsque le vol est agité ou non stabilisé en ligne droite, ou lorsque le calculateur 3 détecte une défaillance d'un capteur ou instrument de mesure, comme décrit ci-après, il fixe la valeurzm (t) = 0 . La masse de sustentation msp n'a donc pas une valeur figée, mais est ajustée lors du vol. La masse estimée m de l'aéronef 1 est alors déterminée lors d'une étape 43 à partir des valeurs de masse au parking mpo et msp. La masse estimée est ainsi du type /Ti = f //valeurs 00, m sp) La masse estimée Fri est par exemple une somme pondérée de la masse médiane moo déterminé lors de l'étape 21 et de la masse de sustentation msp. Les coefficients de pondération sont par exemple calculés à partir de la fonctionzr,(t), de telle sorte que le coefficient attribué à la masse de sustentation soit proportionnel à la fraction de temps lors de laquelle cette masse a pu être déterminée. La masse estimée di est dans cet exemple donnée par : le)(t,)=[(1- p(t,))m pi° + p(t,)msp(tv)1- FU(t,) (18) où la fonction p(t,) est définie par : fi, fx.(t) p(t,). t, -t Ainsi, lorsque le vol est stabilisé depuis longtemps, p(t,) tend vers 1 et la masse estimée m tend vers la masse de sustentation msp dont on a soustrait la masse de carburant consommé. On remarque en outre que tout biais dans la valeur FU(t,), résultant par exemple d'une erreur affectant la débitmétrie de l'aéronef, n'a d'influence que sur le terme (1- p(0)moo, qui est d'autant plus faible que p(tv) est proche de 1, c'est-à-dire que la majorité du vol s'effectue en croisière. Ainsi, lorsque le vol est stabilisé depuis longtemps, un biais dans la valeur FU(tv) n'aura quasiment aucune incidence sur la valeur de la masse estimée m . L'étape 28 comprend par ailleurs une phase de détermination 45 de paramètres de vols auxiliaires, destinés à permettre le calcul de tests des paramètres de vols issus des capteurs ou instruments de mesure, et à remplacer ces paramètres en cas d'identification d'une défaillance d'un ou plusieurs de ces capteurs ou instruments. Ces paramètres auxiliaires sont donc déterminés de manière continue ou périodique à partir de l'instant t2 de décollage de l'aéronef. L'indice n de ces paramètres correspond à la nième valeur déterminée de ces paramètres lorsque la détermination est périodique. La phase 45 comprend notamment la détermination 47 de paramètres de vol à partir de données issues du capteur GPS 10, la détermination 49 d'une première information de vitesse auxiliaire de l'aéronef 1, indépendante des pressions mesurées, la détermination 50 d'une température statique auxiliaire et l'estimation 51 de la vitesse du vent. La connaissance de l'altitude pression de l'aéronef (ou de la pression statique environnant l'aéronef) est d'importance primordiale pour tout procédé de surveillance analytique des paramètres de vol, visant en particulier à détecter une corruption éventuelle des valeurs de pressions statique et totale mesurées. Ainsi, selon l'invention, (19) une information d'altitude issue d'un capteur de position géographique est en permanence disponible lors du vol de l'aéronef 1. Cette information d'altitude est avantageusement issue du capteur GPS 10, l'altitude étant exprimée de manière classique au-dessus du géoïde de référence WGS 84. La phase 47 comprend ainsi la détermination d'une altitude mesurée par le GPS et notée ZGPS, d'un estimateur d'altitude pression, également appelé altitude pression reconstituée et noté Zp**, et d'une pression statique reconstituée notée Ps**, déterminée à partir de l'altitude pression reconstituée Zp** au moyen d'une table d'atmosphère standard. Deux modes sont à distinguer. Selon un premier mode, appelé mode primaire, l'altitude pression reconstituée Zp** est déterminée en retranchant de l'altitude GPS ZGPS un terme correctif d'altitude afin de tenir compte du décalage entre l'atmosphère standard et l'atmosphère réelle de vol. Ce décalage peut être dû à une pression au niveau du sol inférieure, respectivement supérieure, à la pression standard de 1013 bar, générant un décalage des courbes isobares vers le bas, respectivement vers le haut, et à un gradient de température non standard entre le sol et l'aéronef, modifiant l'espacement des isobares. Ce décalage est déterminé à des intervalles DT constants, par exemple égaux à 30 secondes, par comparaison entre l'altitude pression reconstituée Zp** et l'altitude pression Zp issue de l'altimètre 6. Ces écarts constituent une suite de valeurs AZ(j), où j est un entier naturel. Les instants successifs de mesure de ces écarts seront notés t(j). Ainsi : 32(i)= ZGps(t(i))- p(t(i)) (20) Le terme correctif d'altitude est calculé à partir des écarts AZ(i) comme une suite AZ**(j) donnée par la formule AZ**(j) = AZ(i), où i est un indice inférieur ou égal à j, choisi selon l'algorithme suivant : - Si aucune incohérence des paramètres anémo-clino-barométriques n'a été déclarée, c'est-à-dire si les mesures des capteurs de pression statique 5a, de pression totale 5b et d'incidence 5d n'ont pas été jugés non fiables, et si AZ(j)E [AZ(j -1)- ST; AZ(j -1)+ ST] l'indice i est égal à l'indice j, c'est-à-dire A**Z(j). AZ(j). 3.**Z est donc ainsi automatiquement recalé. ST est un seuil de tolérance. La vérification de l'appartenance de AZ(j) à l'intervalle [AZ(j -1)- ST;AZ(j -1)+ ST] revient à analyser analytiquement les sauts de OZ(j) de manière à refuser un recalage aberrant en cas de corruption de l'altitude pression Ze(t(j)). Le seuil de tolérance ST peut par exemple prendre la valeur de 100 ft + 5%IVzI. Le premier terme de 100 pieds est destiné à couvrir un gradient exceptionnel d'isobare en altitude (générateur du vent géostrophique) ainsi que la fluctuation statistique de la valeur d'altitude issue des GPS ; le second terme égal à 5 pourcents de la vitesse verticale VZ(GPS) est destiné à couvrir, en cas de changement d'altitude, un écart de température d'au plus par rapport au standard 30`C en valeur absolue. ST peut avantageusement être adapté suite à l'analyse de résultats d'essais en vol. A titre d'exemple, un vent géostrophique de 200 noeuds au pôle correspond à un gradient d'isobare de 40 pieds en 30 secondes de vol. - Si aucune incohérence des paramètres anémo-clino-barométriques n'a été déclarée, et si AZ(j)« [AZ(j -1)- ST;AZ(j -1)+ ST], la valeur de l'indice i reste égale à sa valeur précédente, et A-Z n'est pas recalé. De plus, le processus de recalage de la suite d'écarts d'altitude AZ**(i) sur la suite AZ est arrêté jusqu'à une relance éventuelle par le pilote comme cela sera décrit plus bas. - Lorsque au contraire les paramètres anémo-clino-barométriques déterminés à partir des capteurs de pression statique 5a, de pression totale 5b et d'incidence 5d sont jugés incohérents donc non fiables, comme décrit ci-après, le processus de recalage de la suite d'écarts d'altitude AZ**(j) sur la suite AZ est arrêté. Ainsi, la valeur i est figée au dernier indice "pertinent", c'est-à-dire correspondant à un instant passé du vol pour lequel l'incohérence ne s'était pas encore manifestée. Ainsi i Cependant, si aucune information disponible à bord ne vient accréditer la validité de ces mesures de pression statique, l'altitude GPS de l'aéronef 1 au-dessus du géoïde et la pression statique environnant l'aéronef 1 sont évaluées selon un mode secondaire, en fonction du QNH local. Le QNH est un code international destiné au calage de l'altimètre de telle sorte que celui-ci indique l'altitude topographique du terrain où est élaboré le QNH lorsque l'aéronef se trouve au sol sur ce terrain. Le QNH local, obtenu ponctuellement par prévision météorologique, ou bien par des moyens radiocom ou datalink, est alors inséré dans ce mode secondaire par le pilote, et le calculateur 3 extrait la pression Pst statique théorique dans une table d'atmosphère standard, correspondant à une altitude pression située à une altitude ZGPS au-dessus de l'isobare du QNH. Dans ce mode secondaire, la pression statique reconsituée Ps** est égale à la pression statique théorique Pst. Il est à noter que cette estimation de la pression statique est d'autant meilleure que l'avion vole près de la surface du sol ou de la mer. Avantageusement, l'entrée dans ce mode secondaire de reconstitution de la pression statique peut se faire par appui sur la même commande que la commande de relance de AZ** en mode primaire. Le duplexage de cette commande peut être réalisé implicitement par l'affichage d'un calage altimétrique par le pilote sur le dispositif d'affichage auxiliaire 14 : si le calage demandé est « standard », l'actionnement de la commande de recalage relancera le processus de recalage automatique selon le mode primaire ; si le calage est de type manuel (QNH affiché), typiquement sous la surface de transition, la commande de recalage activera le mode secondaire de reconstitution de la pression statique. La première vitesse auxiliaire de l'aéronef 1 déterminée lors de la phase 49 est indépendante de la pression statique Ps issue du capteur de pression statique 5a. Elle est avantageusement extrapolée en temps réel à partir de la valeur courante de l'incidence mesurée. Cette vitesse est par exemple un nombre de Mach, appelé Mach fort et noté Ce nombre de Mach est déterminé périodiquement à partir de l'équation de sustentation, sous la forme d'une série par récurrence, selon la relation : [0.7S - P: - 5;(tv) 2 nziflg où : - n est un indice temporel correspondant à un instant tv, (22) - Ps** est la pression statique reconstituée déterminée lors de la phase 47 à l'instant tv, nz est le facteur de charge déterminé lors de la phase 31 à l'instant tv, tri est la masse estimée lors de la phase 39 à l'instant tv, et a*z )= az (a(tv),M,4,conf) est un coefficient de portance estimé à l'instant tv. Le coefficient de portance estimé C*z (tv) est déterminé à partir d'une table de correspondance, fournissant la valeur de a*z ) à partir de la valeur de l'incidence a à l'instant tv déterminée lors de la phase 31, et du dernier Mach fort déterminé, c'est-à-dire le Mach fort Mn*_1 d'indice n-1. Ce coefficient de portance estimé Ô*z(tv) diffère ainsi du coefficient de portance Cz déterminé lors de la phase 33 en ce qu'il ne fait pas intervenir le nombre de Mach issu du Machmètre 8 à partir de la pression statique mesurée par le capteur de pression statique 5a et de la pression totale mesurée par le capteur de pression totale 5b. Le Mach fort Mn ne dépend donc pas des pressions mesurées à l'instant tv. La série Mn* est initialisée à partir d'un premier terme M; . Ce premier terme est une valeur initiale fiable, provenant de mesures indépendantes des capteurs de pressions. La fréquence de détermination du Mach fort est par exemple comprise entre 1 et 10Hz, de préférence 4Hz. Le Mach fort dépend cependant de l'incidence a. En cas de défaillance du capteur d'incidence 5d détectée, cette vitesse Mn* n'est donc plus calculée. Le Mach fort Mn ainsi calculé permet de déduire une estimation de température statique SAT, indépendante des mesures de pression provenant des capteurs 5a ou 5b, un estimateur d'écart de température standard AISA,* , estimant l'écart de la température SAT,* par rapport à la température statique SATed attendue à l'altitude pression reconstituée Zp** de l'aéronef, et une estimation de la vitesse horizontale du vent. La température statique SAT,: est obtenue lors de la phase 50 à partir de la température totale TAT issue du capteur 5c, supposée fiable, et du Mach fort Mn* , en remplaçant dans l'équation (1) le nombre de Mach M par le Mach fort Mn* suivant l'expression : = SAT,* TAT(tv) TAT(tv) (1+ 7-1A4.2) (1-1-o.2n4,72) (23) 2 " L'estimateur d'écart de température standard AISA,* estimant l'écart de la température SAT, par rapport à la température statique SATstd attendue à l'altitude pression GPS Zp** de l'aéronef, est déterminé selon l'expression : AISA,* = SATstd - SAT; (24) la température SATstd étant par exemple déterminée à partir d'une table d'atmosphère standard. Cet estimateur d'écart de température standard n'est pas utilisé tant que l'incidence mesurée par le capteur d'incidence 5d est jugée fiable et que le Mach fort peut être calculé. Cependant, il permet, lorsque cette incidence n'est plus jugée fiable, d'estimer une température statique en corrigeant une température statique SATstd obtenue à partir de l'altitude pression GPS Zp**, comme décrit ci-après. La vitesse horizontale du vent par rapport au sol est déterminée lors de la phase 51 au moyen d'un estimateur de vent horizontal W; , évaluant périodiquement la vitesse du vent en force et en direction. Cet estimateur de vent est avantageusement déterminé à partir du triangle des vitesses selon lequel la vitesse Vs de l'aéronef par rapport au sol, dans un plan horizontal, est la somme de sa vitesse propre Vp et de la vitesse du vent, soit : 17v: = vs> - vp* (25) La vitesse horizontale Vs de l'aéronef par rapport au sol, c'est-à-dire sa vitesse dans un plan horizontal au sol, est estimée à partir de mesures du capteur GPS 10 par exemple. Ainsi Vs = GS . La vitesse propre Vp de l'aéronef est la composante horizontale de sa vitesse par rapport à l'air, c'est-à-dire la composante horizontale de sa vitesse vraie. Une vitesse vraie de l'aéronef, appelée vitesse vraie forte et notée TAS,* , est déterminée à partir du Mach fort Mn* , suivant une relation analogue à la relation (5) : TAS., = Mn yR - SAT; (26) En négligeant la composante verticale de la vitesse du vent, la vitesse de l'aéronef 1 selon l'axe vertical B par rapport au sol est égale à sa vitesse verticale par rapport à l'air. Ainsi, une valeur de la norme de la vitesse propre, appelée vitesse propre forte et notée V; , est déterminée à partir de la vitesse vraie forte à partir de la relation : TASa* 2 = Ve:2 (GPS)2 (27) où vz (GPS) est la vitesse de l'aéronef 1 selon l'axe vertical B par rapport au sol, déterminée à partir de mesures du capteur GPS 10 par exemple. La vitesse propre V; de l'aéronef est donc estimée par : Vp = uhlTASa* 2 - VZ (GPS)2 (28) où le vecteur uh est un vecteur unitaire horizontal portant la projection horizontale du vecteur vitesse par rapport à l'air, estimé par l'intermédiaire de l'incidence a. A partir de l'hypothèse selon laquelle le dérapage (i.e. l'angle entre le vent relatif et le cap avion) est nul et l'incidence a de l'aéronef donnée par le capteur d'incidence 5d est fiable, l'estimateur de vent 14/; horizontal est ainsi déterminé à partir de l'expression : Wr = GS-LZNIM.a2 -(yRSAT;)-Vz(GPS)2 (29) Lors du vol de l'aéronef 1, le calculateur 3 met en oeuvre plusieurs tests permettant de déterminer un état de crédibilité des valeurs de paramètres de vol déterminées au moyen de mesures de capteurs. En particulier, ces tests permettent de déterminer la cohérence des informations de vitesse, de pression et d'incidence, telles que le nombre de Mach Ma, la pression statique Ps et l'incidence a. Le procédé comprend ainsi une étape 70 de détermination d'un état de crédibilité des informations de pression statique, de pression totale et d'incidence, au moyen d'un premier test. Ce test 70 est un test de surveillance analytique, mis en oeuvre en continu lors du vol, tant qu'un dysfonctionnement n'a pas été détecté Le test 70 est basé sur l'équation de sustentation. Il est réalisé en déterminant si les paramètres de vol mesurés vérifient effectivement cette équation de sustentation. Cette équation fait en effet intervenir la pression statique Ps, le nombre de Mach Ma et l'incidence a, via le coefficient de portance en axes avion az . Le test 70 comprend ainsi une phase 72 de détermination de la quantité : )= nz(t+ 0)iii(tvki dans laquelle les paramètres Ps, Ma, az sont ceux déterminés lors de l'étape 28 à un instant tv immédiatement antérieur à l'instant tvk de la détermination 72. L'indice k 0,7SPs(tv)Ma2(tvez Mtv (30) ),Ma(t,),conf) indique qu'il s'agit du Orne test 70 mis en oeuvre par le calculateur 3. nz est déterminé à un instant (t,' + 0), où 0 est un déphasage temporel déterminé une fois pour toutes dans tout ou partie du domaine de vol, afin de maximiser statistiquement une corrélation entre l'incidence oe(t') et le facteur de charge nz (t' +0) sur l'aéronef considéré, comme décrit par la suite. L'équation de sustentation est vérifiée si T1=1 en l'absence de toute erreur dans les mesures et dans la loi az . Ainsi, si les valeurs de la pression statique Ps, du nombre de Mach Ma et de l'incidence a issues de mesures des capteurs sont correctes et si la loi az est correcte, la quantité T1 doit être égale à 1. Un léger écart par rapport à cette valeur peut néanmoins être toléré, en raison notamment du bruit des mesures, de l'écart de la masse estimée m par rapport à la masse réelle de l'aéronef, et de l'imprécision potentielle de la valeur du coefficient de portance az . Le test 70 comprend ainsi une étape 73 lors de laquelle le calculateur 3 vérifie l'appartenance de la valeur T1 à un intervalle restreint autour de la valeur 1, noté 11- + ei est un nombre prédéterminé définissant l'écart toléré, par exemple égal à 0.1, correspondant à une erreur de 10%. Si T1 E 11-e1;1+ ei[, le premier test 70 est positif, et les informations de vitesse, de pression et d'incidence sont jugées de crédibilité optimale 74. Dans l'état de crédibilité optimale 74, les paramètres de vol utilisés directement ou indirectement pour le pilotage de l'aéronef sont donc : l'altitude pression Zp ; la température statique SAT ; le nombre de Mach Ma ; la vitesse vraie TAS ; la vitesse indiquée IAS. La vitesse vraie de l'aéronef est déterminée à partir du Mach Ma, suivant la relation : TAS = Mo\hR - SAT (31) La vitesse indiquée de l'aéronef, notée IAS pour « Indicated Airspeed », est la vitesse directement issue des mesures de pressions, traduites par l'anémomètre 7. Au cours d'une longue phase de vol stabilisée, la fonction p(t) converge vers la valeur 1, de sorte que la masse estimée m(t) converge vers la masse map(t)-FU(t), qui converge elle-même vers la masse m4 (t) issue de l'équation de sustentation. Ainsi, la quantité Ti(tti,r) converge naturellement vers la valeur 1 lors d'un vol stabilisé en croisière. Ceci est obtenu même si la loi az comporte un biais ou si la mesure du capteur d'incidence est biaisée. Dans l'état de crédibilité optimale 74, le dispositif d'affichage auxiliaire 12 affiche un ensemble de caractéristiques de vol issues des mesures de capteurs de l'aéronef 1, déterminées notamment à partir des pressions statique et totale et de l'incidence mesurées. Par exemple, le dispositif d'affichage auxiliaire 12 affiche l'altitude pression Z,,, le nombre de Mach Ma et la vitesse indiquée IAS. Le dispositif d'affichage auxiliaire 12 affiche également la pente air yair et l'incidence a mesurées. Le dispositif d'affichage auxiliaire 12 affiche en outre une vitesse indiquée de décrochage de l'aéronef 1, notée IAS,, sur la base du facteur de charge n, courant. Il s'agit donc de la vitesse IAS minimale qui peut être atteinte par l'aéronef 1 sans décrocher à la valeur du facteur de charge courante. Cette vitesse de décrochage IAS, est déterminée par le calculateur 3 en fonction de l'incidence as de décrochage, de la vitesse indiquée IAS et de l'incidence a courantes, selon la relation : = IAS (32) où ao est l'incidence de portance nulle dans la configuration actuelle de l'aéronef. Elle est généralement négative, sa valeur absolue pouvant atteindre quelques degrés au braquage de volets le plus élevé. Cette incidence ao est tabulée suivant la configuration de l'aéronef. Si au contraire Ti ;1 + T1 correspondant à un instant Cl , le test 70 est négatif, ce qui signifie que l'une au moins des informations de pression statique, de pression totale ou d'incidence est erronée. Le calculateur 3 estime alors que les mesures de vitesse, d'altitude barométrique et d'incidence ne sont plus fiables, et active dans une étape de transition 75 un état de crédibilité faible 76. Tant que cet état de crédibilité faible 76 est activé, les paramètres de vol déterminés par le calculateur 3 à un instant tv et destinés à l'aide au pilotage sont indépendants des valeurs de pression statique, de pression totale et d'incidence issues de capteurs à cet instant tv, même s'ils peuvent dépendre de valeurs antérieures mesurées dans l'état de crédibilité optimale. Lors de l'étape de transition 75, le calculateur 3 fige la valeur de l'écart de température standard AISA,* et du décalage d'altitude pression AZ-(j) aux dernières valeurs fiables connues. Ces dernières valeurs fiables connues sont celles déterminées un certain laps de temps avant le test 70 précédent, à un instant noté tf : tf est le dernier instant auquel les paramètres de vol ont été considérés comme fiables. On choisira par exemple t, tyk1 _ 30s . Notamment, la mesure de température totale TAT par le capteur de température totale 5c est jugée non fiable, de telle sorte que le recalage de l'écart de température standard AISA,* à partir des relations (23) et (24) ne peut être effectuée. L'écart de température standard AISA,* prend donc la valeur constante A/SAff- = AISA,* (tf). La température statique est alors estimée à partir l'altitude pression reconstituée Zp** selon la relation : SAT- = SATstd(Zp)+ A/SAtf** (33) SAT- est appelée température statique faible. La température statique SAT,* n'est quant à elle plus déterminée. Le décalage d'altitude pression AZ- (j), utilisé pour la détermination 47 de l'altitude pression reconstituée Zp- , prend par ailleurs la valeur constante AZtf- = AZ (tf) (la valeur de i est figée, comme décrit ci-dessus). Cependant, comme décrit précédemment, cette valeur AZtf peut être modifiée à tout moment par action manuelle du pilote sur la commande de relance, par exemple en basse altitude avant l'atterrissage si le pilote vérifie que son altitude pneumatique indiquée au calage QNH est cohérente avec son altitude géoïde donnée par un instrument de secours ou sa hauteur radiosonde augmentée de l'altitude du terrain. En outre, la valeur de l'incidence n'est pas fiable, de telle sorte que l'estimateur de vent ne peut plus être déterminé selon la relation (29) ci-dessus. Cet estimateur, alors noté W- , prend alors une valeur constante W;* = Mt, (tf). Cependant, la valeur de l'estimateur de vent peut être mise à jour ponctuellement, de manière automatique, lorsque la trajectoire de l'aéronef présente une évolution suffisante, selon un critère décrit ci-après. La mise à jour ou recalage de l'estimateur de vent horizontal W en crédibilité faible est réalisée suivant une méthode similaire à la méthode de détermination de l'estimateur de vent horizontal Wn* en crédibilité forte, en substituant à l'incidence a mesurée par le capteur d'incidence une incidence reconstituée notée aR , et à la vitesse vraie forte TAS,* une vitesse vraie estimée notée TÂS . L'évaluation de l'incidence reconstituée aR est réalisée à partir de l'équation de sustentation de l'aéronef, en restant dans un domaine de vol dans lequel le coefficient de portance de l'aéronef peut être exprimé comme une fonction affine de l'incidence aR , dont les coefficients sont prédéterminés, selon l'expression : 5; (aR A4) = 2(M)+P(M)'aR (34) valable dans un domaine restreint d'incidence [ao(M),ai (M)]. Ainsi, l'incidence reconstituée est donnée par : az - 2(M) nzifig 2(M) «R - (35) ,u (Ad - 0,7SP;* M **2,u(M) ,u (M) De même, la dérivée de l'incidence reconstituée aR est exprimée en fonction de la dérivée du facteur de charge nz, par : g = tri nz (36) p(M)0,7SPs** M *2 ) Les valeurs de 2(M) et de ,u(M) sont déterminées à partir des tables donnant la valeur du coefficient de portance en fonction de l'incidence et du nombre de Mach. La vitesse vraie estimée TÂS est alors estimée comme un quotient entre une accélération de l'aéronef et une vitesse angulaire de cet aéronef, sur la base de l'incidence reconstituée aR et de sa dérivée. Cette estimation est valable sous réserve que l'évolution de la trajectoire de l'aéronef soit suffisamment nette, selon le critère décrit ci-après, et en supposant que le dérapage est nul. Le critère utilisé pour estimer si l'évolution de la trajectoire est suffisante repose sur la formule suivante qui exprime, par le biais de la relation fondamentale de la dynamique, que la projection de l'accélération de l'avion dans le référentiel terrestre, projetée dans le plan de l'aéronef, est supérieure à une valeur déterminée, avantageusement 0.4g : gn + g > 0,40 g (37) où h' est le facteur de charge et gn +g l'accélération de l'aéronef. Ainsi, l'évolution de la trajectoire de l'aéronef est jugée suffisante si : (cos 0 sin yoi )2 + (nz - cos çai cos 0)2 > 0,16 (38) où °désigne l'assiette de l'aéronef et çoi son angle de gîte. On remarque notamment qu'il suffit d'avoir (cos °sin qyi )2 > 0,16 , donc d'incliner suffisamment l'avion (quelles que soient la poussée, la traînée et la portance) pour obtenir une évolution de trajectoire "suffisante". Typiquement, lorsque l'aéronef 1 est en trajectoire proche du vol rectiligne, le calcul de la vitesse vraie estimée sera déclaré invalide et l'estimateur de vent horizontal W'' ne sera pas mis à jour. A l'inverse, lorsque l'évolution de la trajectoire est estimée suffisante, la vitesse vraie estimée TÂS est estimée, ce qui permet de recaler l'estimateur de vent horizontal W selon l'expression : W** =GS-t7hIITÂS2 - Vz (GPS)2 (39) où vh est un vecteur unitaire horizontal portant la projection horizontale du vecteur vitesse air estimé par l'intermédiaire de l'incidence reconstituée aR. Lors de l'étape de transition 75, le calculateur 3 fixe par ailleursxm (t) = 0 , de sorte que les valeurs de pression statique, de Mach et d'incidence ne sont plus utilisées pour le calcul de la masse estimée. Ainsi, la masse estimée bi à un instant tt, > tf est donnée par : FU[tf;tv] (40) où FU[tf;ti désigne la masse de carburant consommée entre tfet tv. L'estimation de la vitesse horizontale du vent et de l'écart de température standard permet, en complément de l'altitude pression reconstituée Zp**, de calculer des estimations de vitesse de l'aéronef 1 à partir de sa vitesse GS par rapport au sol. Ces estimations sont par exemple une estimation du nombre de Mach de l'aéronef, appelé Mach faible et noté M**, et une estimation d'un équivalent de vitesse, appelé équivalent de vitesse faible et noté EV-. La vitesse GS de l'aéronef 1 par rapport au sol est par exemple déterminée par le calculateur 3 à partir du capteur GPS 10. Elle peut également être déterminée à partir de l'accélération de l'aéronef telle que mesurée par des centrales à inertie IRS (Inertial Reference System) de l'aéronef 1, éventuellement en hybridant les informations reçues du capteur GPS et des centrales à inertie. La vitesse GS de l'aéronef 1 et l'estimation de la vitesse horizontale du vent W- , égale à la valeur figée Wt*f* ou à une valeur recalée, permettent de déterminer un vecteur vitesse propre de l'aéronef, appelé vitesse propre faible et noté Vp donné par : Vp = GS-W- (41) La vitesse verticale Vz (GPS) de l'aéronef 1 selon l'axe vertical B par rapport à l'air est déterminée à partir de mesures du capteur GPS 10, en négligeant le vent vertical. Une vitesse vraie de l'aéronef, appelée vitesse vraie faible et notée TAS**, est alors déterminée comme une norme du vecteur vitesse de l'aéronef par rapport à l'air, somme de la vitesse propre et de la vitesse verticale, selon la relation : TAS (t). GS+1/z(GPS)k-W (42) dans laquelle : GS est la vitesse horizontale de l'aéronef 1 par rapport au sol à l'instant t ; Vz(GPS) désigne la vitesse de l'aéronef 1 selon l'axe vertical B à l'instant t ; est un vecteur unitaire ascendant parallèle à l'axe vertical B ; W est l'estimateur de vent figé ou recalé. En fonction de l'altitude pression estimée Zp et de l'écart de température standard A/SAtf** figé, on estime par ailleurs la température statique environnante SAT- =SATstd(Zp**)+AISAf , et la vitesse estimée du son .\/),R. SAT** . Le Mach faible M** est ensuite déterminé à partir de la vitesse vraie faible et de la vitesse du son estimée, par la relation : M(t> tf). TAS *(t) .\/,,R - SAT- c'est-à-dire GS+Vz(GPS)k -wf*f* .NlyR - SAT- Le Mach faible M** est ainsi indépendant de toute valeur de pression mesurée postérieurement à tf. Le calcul périodique du Mach fort Mn est interrompu, car ce calcul dépend de l'incidence a. m (t > tf ). (43) (44) En outre, un équivalent de vitesse faible, noté EV**, est déterminé à partir du Mach faible M**et de l'altitude pression reconstituée Z*; selon la relation : 1 poEV**2 = 2 yf's M (45) 2 L'équivalent de vitesse faible EV** est ainsi obtenu à partir du Mach faible M et de la pression statique reconstituée P;* par la relation : EV** = Y's - M , la pression Po statique reconstituée P;* étant elle-même déterminée à partir de l'altitude pression reconstituée Zp** selon une table de l'atmosphère standard. Dans l'état de crédibilité faible 76, les paramètres de vol utilisés directement ou indirectement pour le pilotage de l'aéronef, à un instant tv, sont donc : l'altitude-pression GPS Zp** = ZGPS OZ ,f ; **10 la température statique SAT = SAT std(Z p )+ AIS& ; - l'estimateur de vent horizontal W figé à la valeur W;* ou recalé ; - le Mach faible M ; - la vitesse vraie faible TAS** ; - l'équivalent de vitesse faible EV**. 15 Ainsi, en crédibilité faible, les seules mesures en temps réel utilisées sont celle issues du capteur GPS. Pour affiner le test 70 et discriminer les erreurs sur la pression statique Ps ou sur la pression totale PT, dues à un dysfonctionnement des sondes, d'une erreur sur la mesure de l'incidence a, due à un dysfonctionnement du capteur d'incidence 5d, le 20 calculateur 3 effectue un deuxième test 78. Ce test 78 est effectué dans la suite du test 70, pendant un délai d limité après t ; par exemple de 5 secondes. Ce deuxième test 78 est destiné à déterminer la crédibilité de l'incidence a telle que mesurée lors de la phase 31. Ce test 78 d'incidence est réalisé en déterminant si l'équation de sustentation est 25 vérifiée, en remplaçant dans cette équation la pression statique Ps mesurée par le capteur de pression statique 5a et le nombre de Mach Ma déterminé par le machmètre 8 à partir de mesures des capteurs de pression statique et totale, par des paramètres de vols auxiliaires correspondants. Ces paramètres de vols auxiliaires correspondants sont la pression statique 30 reconstituée Ps**, déterminée conformément à la phase 45, et le Mach faible M-déterminé lors de l'étape 76. En outre, le coefficient de portance, noté az**, est estimé à partir de la valeur de l'incidence a à l'instant tv déterminée lors de la phase 31, et du Mach faible M**. Le test 78 comprend la détermination 80 du rapport T2 vérifiant : T2 (tvkl + d > tv > tvkl )= nz (tv + 0)iii(tv où a()= Cz (cr(tv),M (tv)conf). Cette étape 80 est suivie d'une étape 82 lors de laquelle le calculateur 3 vérifie l'appartenance de T2 à l'intervalle 11- e2;1+ ez [, où E2 est un paramètres de tolérance prédéterminé, par exemple égal à el. Le rapport T2 ne dépend pas de la pression statique et du nombre de Mach mesurés. Ainsi, un écart de T2 par rapport à la valeur 1 théorique attendue incrimine spécifiquement la mesure de l'incidence a, car il s'agit de la seule mesure postérieure à tvid autre que des mesures issues du capteurs GPS, qui sont quant à elles supposées fiables. Si Tz « ]1- e2;1+ ez [, le calculateur 3 identifie une défaillance du capteur d'incidence 5d, n'excluant pas une défaillance des capteurs de pression. Le calculateur 3 maintient ainsi l'état de crédibilité faible 76. Le dispositif d'affichage auxiliaire 14 affiche alors des caractéristiques de vol auxiliaires ne faisant pas intervenir les mesures issues des capteurs de pression 5a, 5b et d'incidence 5d, ces mesures étant jugées peu fiables. Ainsi, chacune des caractéristiques de vol faisant intervenir une de ces mesures dans l'état de crédibilité optimale est remplacée par une caractéristique homologue, i.e. représentative de la même caractéristique de vol, sur le dispositif d'affichage auxiliaire 14. Par exemple, l'altitude pression Zp, le nombre de Mach Ma et la vitesse indiquée IAS sont remplacés sur le dispositif d'affichage auxiliaire 12 par des caractéristiques homologues, qui sont respectivement l'altitude pression reconstituée Zp** ou l'altitude ZGPS, selon le choix du pilote, le nombre de Mach faible M**, et l'équivalent de vitesse faible EV**. En outre, la pente air yak et l'incidence a ne sont plus affichées, et disparaissent donc du dispositif d'affichage auxiliaire 14. Par ailleurs, la vitesse indiquée de décrochage IASS est remplacée par un équivalent de vitesse de décrochage faible EV: sur la base du facteur de charge n, courant. L'équivalent de vitesse de décrochage faible EVS est déterminé par le 0,7SP;* M**2 (tv)az* (tv) (46) calculateur 3 en fonction du facteur de charge n, courant et d'un équivalent de vitesse de décrochage EVso sous une accélération de 1g, selon la relation : EV;* = EVso V nz (47) L'équivalent de vitesse de décrochage EVso est obtenu à partir d'une table de décrochage en fonction de la masse estimée in de l'aéronef 1. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité faible. Ce message indique ainsi au pilote que les mesures des capteurs de pression et d'incidence ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Si au contraire T2 E e2;1+ e2[, et sous réserve qu'un test complémentaire évaluant une cohérence dynamique de l'incidence mesurée avec le facteur de charge, détaillé ci-après, n'ait pas fourni un résultat négatif pendant une durée prédéterminée, par exemple 5 minutes, le calculateur 3 exclut une défaillance du capteur d'incidence 5d, ce qui implique une défaillance d'au moins un des capteurs de pression statique et totale. Le calculateur 3 lève ainsi l'état de crédibilité faible, et active dans une étape 83 un état de crédibilité forte 84. Lors de la transition 83, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité forte. Ce message indique que les mesures des capteurs de pression ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Dans cet état de crédibilité forte 84, les paramètres de vol déterminés à un instant tv sont indépendants des valeurs de pression statique et totale issues de capteurs à cet instant tv, mais peuvent dépendre de l'incidence issue du capteur d'incidence 5d. Lors de la transition 83 d'un état de crédibilité faible 76 à un état de crédibilité forte 84, le calculateur 3 réinitialise le calcul de la suite des Mach forts, interrompu dans l'état de crédibilité faible, à partir d'un terme initial M; correspondant à la dernière valeur du nombre de Mach faible avant la transition 83. Le terme initial gest ainsi remis à jour à chaque fois que l'état de crédibilité forte 84 est activé. La température statique SAT; est à nouveau évaluée, suivant l'expression : TAT(tv) SAT; = (48) + 0.2Mn*2) Par ailleurs, le calculateur 3 reprend la détermination de l'estimateur d'écart de température standard AISAn* , interrompue dans l'état de crédibilité faible en raison de l'indisponibilité du nombre de Mach fort Mn* . Cet estimateur est déterminé conformément à l'étape 50 ci-dessus. Le calculateur 3 reprend également le calcul de l'estimateur de vent W,*, , évaluant périodiquement la vitesse du vent par rapport au sol, conformément à la phase 51 ci- dessus. Ce calcul fait en effet appel au nombre de Mach fort et à la température statique forte SAT; , à nouveau disponibles. La vitesse vraie forte TAS,* est déterminée à partir du nombre de Mach fort Mn* et de la température statique SAT; selon l'expression (23) ci-dessus. En outre, un équivalent de vitesse fort, noté EV,* , est déterminé à partir du Mach fort Mn et de l'altitude pression reconstituée Zp** selon la relation : 2 2 * poEV'2 = -1 2Ps 114,7*2 (49) Ainsi, les paramètres de vol utilisés directement ou indirectement pour le pilotage de l'aéronef 1 en crédibilité forte, sont : l'altitude pression GPS 4,* = Z -GPS AZ la température statique SAT; ; l'incidence mesurée a ; le Mach fort Mn* ; la vitesse vraie forte TAS,* ; l'équivalent de vitesse fort EV,,* . Dans l'état de crédibilité forte 84, le dispositif d'affichage auxiliaire 14 affiche des caractéristiques de vol auxiliaires, ne faisant pas intervenir les mesures issues des capteurs de pression. Ainsi, chacune des caractéristiques de vol faisant intervenir une mesure d'un capteur de pression dans l'état de crédibilité optimale est remplacée par une caractéristique homologue sur le dispositif d'affichage auxiliaire 14, indépendante de toute mesure de capteurs de pression réalisée dans l'état de crédibilité forte. Par exemple, l'altitude pression Zp est remplacée sur le dispositif d'affichage auxiliaire 14 par une caractéristique homologue, qui est l'altitude pression reconstituée Zp** ou l'altitude ZGPS, selon le choix du pilote. Le nombre de Mach Ma est remplacé par le nombre de Mach fort Mn* , et la vitesse indiquée IAS est remplacée par l'équivalent de vitesse fort EV; . La pente air yak et l'incidence a mesurées restent affichées. Le dispositif d'affichage auxiliaire 12 affiche en outre un équivalent de vitesse de décrochage fort noté EV; , déterminé par le calculateur 3 de manière analogue à la vitesse indiqué de décrochage IASs, en fonction de l'équivalent de vitesse fort EV; , de l'incidence as de décrochage et de l'incidence a courantes, selon la relation : EV; = EV,*, a-ao es - eo (50) La détermination de ces caractéristiques de vol et de l'air ambiant permet ainsi de fournir des informations au pilote de l'aéronef tout en s'affranchissant des pressions statique et totale potentiellement erronées mesurées par les capteurs. Comme précédemment décrit, dans l'état de crédibilité forte, la température totale mesurée par le capteur 5c et l'incidence mesurée par le capteur 5d sont jugées fiables, ce qui permet, grâce au calcul de l'altitude pression reconstituée Zp** , de déterminer un nombre de Mach (Mach fort) et une vitesse vraie (vitesse vraie forte) de l'aéronef. Ces informations de vitesse permettent de déterminer une image de l'atmosphère, et en particulier un estimateur de vent horizontal et un estimateur d'écart de température standard. A l'inverse, dans l'état de crédibilité faible, la température totale mesurée par le capteur 5c et l'incidence mesurée par le capteur 5d sont jugées non fiables. L'altitude pression reconstituée 4, l'estimateur de vent horizontal et l'estimateur d'écart de température standard sont alors figés aux dernières valeurs fiables connues, l'altitude pression reconstituée Zp** et l'estimateur de vent horizontal pouvant cependant être recalés. En outre, la valeur figée de l'écart de température standard AISA'* (telle que déterminés dans l'état de crédibilité forte) permet d'estimer une température statique à partir de l'altitude pression reconstituée Z. . La connaissance de cette image simplifiée de l'atmosphère permet alors de déterminer un nombre de Mach (Mach faible) et un équivalent de vitesse (équivalent de vitesse faible) de l'aéronef. Quel que soit l'état de crédibilité forte ou faible en vigueur dans le système, le calculateur 3 met en oeuvre des tests complémentaires pour évaluer, de manière périodique, la crédibilité des informations fournies par les capteurs de pression et/ou d'incidence, détecter une éventuelle évolution de cette crédibilité et confirmer ou modifier un état de crédibilité faible ou fort. Dans l'état de crédibilité faible 76, les mesures issues des capteurs de pression 5a et 5b et d'incidence 5d sont jugées peu fiables. Pour réévaluer la crédibilité de ces mesures, le calculateur 3 évalue de manière répétée la fiabilité de la mesure d'incidence a. Si cette mesure est jugée non fiable, le calculateur 3 maintient l'état de crédibilité faible. Si au contraire cette mesure est jugée fiable, le calculateur 3 lève l'état de crédibilité faible, et active un état de crédibilité forte 84, dans lequel seules les mesures de pression sont jugées peu fiables. Dans l'état de crédibilité forte 84, les mesures issues des capteurs de pression 5a, 5b sont jugées peu fiables, tandis que les mesures du capteur d'incidence 5d sont estimées fiables. Si les mesures d'incidence sont jugées fiables lors de l'évaluation par les tests complémentaires, le calculateur 3 maintient l'état de crédibilité forte 84. Si au contraire les mesures d'incidence sont jugées peu fiables lors de l'évaluation par les tests complémentaires, le calculateur 3 active l'état de crédibilité faible 76. Les étapes mises en oeuvre par le calculateur 3 pour évaluer la fiabilité des mesures du capteur d'incidence 5d dans l'état de crédibilité forte 84 ou faible 76 sont illustrées sur la Figure 4. Cette réévaluation repose sur deux tests de crédibilité de l'incidence. Un premier test 90 mesure la cohérence dynamique de l'incidence mesurée a avec le facteur de charge nz. Une modification de l'incidence a entraîne une variation du facteur de charge n2, avec un léger retard noté 0. Si l'on considère un échantillon de temps, lors duquel l'incidence varie de da, suffisamment court pour que les valeurs m et PsMa2puissent être considérés comme constantes, la dérivation de l'équation de sustentation sur cet échantillon de temps donne : dnz = 0,7spsma2 `ta6z da (51) mg)^ aa aaz est constante sur l'échantillon de aa Si l'on considère de plus que la grandeur temps, on obtient dnz = K.doe , où K est une constante. Par conséquent, si la valeur de l'incidence mesurée par le capteur d'incidence 5d est juste, cette incidence doit être une fonction affine du facteur de charge nz sur un échantillon de temps court. Le test 90 repose ainsi sur l'évaluation de la covariance entre l'incidence a et le facteur de charge n2, mesurant la corrélation entre ces deux paramètres. Le test 90 comprend une étape 92 de détermination de l'applicabilité du test 90, une étape 94 de calcul d'un coefficient de corrélation entre le facteur de charge et l'incidence, et une étape 96 de comparaison de ce coefficient de corrélation à une valeur seuil prédéterminée. Le test 90 ne peut être effectué que si l'incidence et le facteur de charge varient suffisamment sur l'échantillon de temps considéré. La détermination 92 de l'applicabilité du test 90 comprend l'évaluation d'une variance du facteur de charge nz, mesurant la dispersion du facteur de charge nz autour de sa moyenne sur un intervalle de temps choisi, et appelé énergie empirique brute. L'intervalle de temps est par exemple égal à 2.5 secondes. L'énergie empirique brute en, du facteur de charge est alors donnée par : 1 100 (52) e2nz 100 k = flz (to +_k 40 + çb)-nz =1 où nz , moyenne du facteur de charge nz sur l'intervalle [to + 0;t0 +2.5+0], est donné par : ___ 1 100 nz = Enz(to +-k 40 +0) (53) 100 k=1 0 est le déphasage entre le facteur de charge et l'incidence et to un instant de mesure choisi. Oest par exemple égal à 0.01 seconde. Ce déphasage 0 est dépendant d'un type d'avion. Il est déterminé une fois pour toutes dans tout ou partie du domaine de vol afin de maximiser en moyenne, dans la partie du domaine de vol considérée, la corrélation entre a(t) et nz(t+0). Lors de l'étape 92, le calculateur 3 compare ensuite la valeur de l'énergie empirique brute en, du facteur de charge à une valeur seuil minimale enzmin prédéterminée. La valeur seuil en,min est par exemple comprise entre 0,01 et 0,02, avantageusement égale à 0,015. Si enz L'étape 94 de calcul d'un coefficient de corrélation entre le facteur de charge et l'incidence comprend le calcul d'une covariance empirique brute entre l'incidence a et le facteur de charge n, sur un échantillon glissant des valeurs mesurées. La covariance empirique brute Covo(nz,a) entre nz (t + 0)et a(t) à l'instant t=to est déterminée suivant l'expression : 1 100 f ( ( Covo(nz,a)=-E nz to +- 0)-nz to 100 40 k=1 \ 40 1 100 k nz = En4t0+- 40 +0 où 100 k=1 100 . 1 v cer t_ 100 ° 40 A partir de ces grandeurs, le calculateur 3 détermine lors de l'étape 94 un coefficient de corrélation brute Corro(nz,a)entre le facteur de charge et l'incidence, selon l'expression : Corro(nz,a)=Covo(nz,a) en,ea où eoe désigne l'énergie empirique brute de l'incidence a dans l'intervalle de temps 100( k - 2 [to ; to + 2,5], donnée par : e2a = 10010° + 40)- a - k=1 Le coefficient de corrélation est par définition compris dans l'intervalle [- 1;+ 1]. Plus sa valeur est proche de 1, plus l'incidence et le facteur de charge sont corrélés. Le déphasage 0 est tabulé pour comprendre une pluralité de valeurs dépendant notamment du domaine de vol considéré. Puis, dans l'étape 96, le calculateur 3 compare le coefficient de corrélation déterminé lors de l'étape 94 à un seuil de corrélation prédéterminé Corrmin Si Corro(nz,a) ço,' et/ou Pa > Prra x , le test 100 n'est pas poursuivi. Le calculateur 3 réitère alors l'étape 102 jusqu'à ce que le test 100 puisse être effectué. Si ço çomax et P É Pniax , le test 100 peut être poursuivi et le calculateur 3 passe à l'étape 106. Lorsque çon,ax et Pa É Pmax, c'est-à-dire lorsque la pente air et l'inclinaison de l'aéronef sont quasiment nulles, si on néglige la composante verticale du vent, l'assiette Ode l'aéronef est sensiblement égale à son incidence. Lors de l'étape 106, le calculateur 3 évalue la fiabilité du capteur d'incidence 5d en comparant l'assiette O, déterminée au moyen de la centrale à inertie, à l'incidence a mesurée par le capteur d'incidence 5d. Par exemple, le calculateur 3 détermine la différence (0 - a)et vérifie que cette différence est proche de la valeur 0, c'est-à-dire est qu'elle est inférieure à un seuil e3 en valeur absolue, où e3 est un nombre prédéterminé définissant l'écart toléré. On prendra par exemple e3 = 1° . Si 10-a > e3 , le test 100 est négatif. Le calculateur 3 réitère alors l'étape 80 tant qu'un résultat positif n'a pas été obtenu. Si IO - al 5_ e3, le test 100 est positif. Ce test 100 permet de valider la valeur absolue de l'incidence en vol rectiligne horizontal. Il permet ainsi d'éliminer des erreurs constantes sur la valeur de l'incidence. Les tests 90 et 100 sont répétés en permanence par le calculateur 3. Les deux tests 90 et 100 ne sont pas toujours en mesure de donner un résultat de manière concomitante. En effet, le test 90 de corrélation entre l'incidence et le facteur de charge est applicable lorsque le facteur de charge varie, tandis que le test 100 de comparaison entre l'incidence et l'assiette est applicable lorsque l'aéronef est en palier rectiligne. Les résultats positifs à ces deux tests ne sont donc pas nécessairement obtenus ensemble. Ainsi, lorsque le calculateur 3 obtient un résultat positif à l'un des tests 90, 100, il active 107 un état intermédiaire et mémorise ce résultat pendant une durée prédéterminé At, en attente d'un résultat positif ou négatif à l'autre de ces tests 100, 9. La durée At est par exemple égale à 5 min (300 s). Si un résultat positif est obtenu pour le deuxième de ces tests 100, 90 pendant cette durée At, le calculateur 3 considère que les mesures du capteur d'incidence 5d sont fiables. A noter que si ces résultats sont obtenus durant l'intervalle [tvki,tvki + di, c'est-à- dire pendant le délai d après un test 70 négatif, ils ne sont pas pris en compte. Si l'état de crédibilité en cours est l'état de crédibilité forte 84, le calculateur maintient cet état. Si l'état de crédibilité en cours est l'état de crédibilité faible 76, le calculateur 3 lève cet état de crédibilité faible 76 et active l'état de crédibilité forte 84, précédemment décrit. Notamment, la transition entre ces deux états est similaire à la transition 83 décrite précédemment. Le test permanent 70 est alors réactivé dès l'entrée en crédibilité forte. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche des caractéristiques de vol auxiliaires adaptées à cet état de crédibilité forte 84. Notamment, l'altitude pression GPS Z; reste affichée, mais le nombre de Mach faible M** est remplacé par le nombre de Mach fort Mn* , et l'équivalent de vitesse faible EV** est remplacé par l'équivalent de vitesse fort EV,* En outre, la pente air yak et l'incidence a sont à nouveau affichées par le dispositif d'affichage auxiliaire. Enfin, l'équivalent de vitesse de décrochage fort EV; est à nouveau déterminé à partir de mesures d'incidence issues du capteur d'incidence. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité forte. Ce message indique au pilote que les mesures des capteurs de pression 5a, 5b ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Si aucun résultat positif n'est obtenu pour le deuxième de ces tests 100, 90 pendant la durée At, le calculateur 3 considère que les mesures du capteur d'incidence 5d ne sont pas fiables. Si l'état de crédibilité en cours est l'état de crédibilité faible 76, le calculateur maintient cet état de crédibilité faible 76 et réitère les tests 90 et 100. Si l'état de crédibilité en cours est l'état de crédibilité forte 84, le calculateur 3 active l'état de crédibilité faible 76. Notamment, lors de la transition 75 vers l'état de crédibilité faible 76, les valeurs de l'estimateur de vent W,*, , de l'écart de température standard AISA,* et du décalage d'altitude pression AZ* (j) sont figées aux dernières valeurs fiables connues, qui sont par exemple celles déterminées à un instant tf antérieur au test 90 ou 100 négatif. L'altitude pression reconstituée Zp** et l'estimateur de vent horizontal peuvent cependant être recalés comme décrit ci-dessus. Le dispositif d'affichage auxiliaire 14 affiche par ailleurs des caractéristiques de vol auxiliaires adaptées à cet état de crédibilité faible 76. Notamment, l'altitude pression GPS Zp** reste affichée, mais le nombre de Mach fort Mn* est remplacé par le nombre de Mach faible M**, et l'équivalent de vitesse fort EV; est remplacée par l'équivalent de vitesse faible EV**. En outre, la pente air y, , r et l'incidence a ne sont plus affichées, et disparaissent du dispositif d'affichage auxiliaire 14. Enfin, l'équivalent de vitesse de décrochage faible Ey;* est obtenu à partir d'une table de décrochage en fonction du facteur de charge n, courant et d'un équivalent de vitesse de décrochage EVso sous une accélération de 1g. L'équivalent de vitesse de décrochage EVso est obtenu à partir d'une table de décrochage en fonction de la masse estimée iii de l'aéronef 1. Par ailleurs, le dispositif d'affichage auxiliaire 14 affiche un message d'alerte destiné à informer le pilote de l'état de crédibilité faible. Ce message indique au pilote que les mesures des capteurs de pression et d'incidence ne sont pas fiables, et que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. La combinaison des tests 90 et 100, mesurant respectivement la corrélation entre l'incidence et le facteur de charge et la cohérence entre l'incidence et l'assiette en vol horizontal, permet ainsi d'évaluer la fiabilité des mesures d'incidence par le capteur d'incidence 5d. Comme précédemment décrit, l'état de crédibilité faible 76 est activé lorsque le premier test 70 et le deuxième test 78 s'avèrent successivement négatifs. L'état de crédibilité faible 76 est également activé de manière systématique au décollage, et pendant chaque changement de configuration de l'aéronef 1. Après le décollage, l'état de crédibilité faible 76 est levé et l'état de crédibilité forte 84 activé dès que la configuration « deuxième segment » est atteinte, sous réserve que le test 90 de corrélation entre l'incidence et le facteur de charge ne soit pas négatif. De même, dès que le changement de configuration de l'aéronef 1 est achevé, l'état de crédibilité faible 76 est levé et l'état de crédibilité forte 84 établi, sous réserve que le test 90 de corrélation entre l'incidence et le facteur de charge ne soit pas négatif. Par ailleurs, dans l'état de crédibilité forte 84, le calculateur 3 réévalue la fiabilité des capteurs de pression 5a, 5b. Ces tests sont mis en oeuvre de manière permanente dans l'état de crédibilité forte 84. Pour réévaluer la fiabilité des capteurs de pression dans l'état de crédibilité forte, le calculateur 3 met en oeuvre deux tests 110, 112, comme représenté sur la Figure 5. Dans un premier test 110 des capteurs de pression, le calculateur 3 compare le nombre de Mach Ma obtenu par le machmètre 8 à un instant donné à partir des valeurs de pression statique et totale au nombre de Mach fort Mn obtenu au même instant, qui ne dépend pas de ces pressions. En effet, si les valeurs de pression mesurées sont correctes, ces deux vitesses devraient être sensiblement égales. Ainsi, lors du test 110, le calculateur 3 détermine le nombre de Mach Ma conformément à la phase 31 ci-dessus, et le nombre de Mach fort Mn* , conformément à l'étape 49 ci-dessus. Le calculateur 3 détermine alors le rapport Ma, et vérifie que ce Mn rapport est proche de la valeur 1, c'est-à-dire est compris dans un intervalle /1- e4;1+e4[, où e4 « 1 , est un nombre prédéterminé définissant l'écart toléré. Si Ma e -(4;1+ e4[, le test 110 est négatif. Les valeurs des pressions mesurées Mn sont jugées non fiables, et l'état de crédibilité forte est maintenu. Si Ma E - e4;1+ e4[, le test est positif. Mn Le deuxième test 112 des capteurs de pression est complémentaire au premier test 110. Lors de ce test 112, le calculateur 3 compare l'altitude pression Zp obtenue à un instant donné à partir de la mesure de la pression statique à l'altitude pression GPS Zp** obtenue au même instant, qui ne dépend pas de la mesure de la pression statique. Si la valeur de pression statique mesurée est correcte, ces deux altitudes pression devraient être sensiblement égales. Ainsi, lors du test 112, l'altimètre 6 fournit au calculateur 3 l'altitude pression Zp conformément à la phase 31 ci-dessus, et le calculateur 3 détermine l'altitude pression GPS Zp**, conformément à l'étape 47 ci-dessus. Le calculateur 3 détermine alors la différence Z*p* -Zp , et vérifie que ce rapport est proche de zéro, c'est-à-dire est compris dans un intervalle e5;+ e5[, où E5 est un nombre prédéterminé définissant l'écart toléré. Si Zp - Zp e5;+ e5 [ , le test 112 est négatif. La valeur de la pression statique mesurée est jugée non fiable, et l'état de crédibilité forte est maintenu. Si Zp* - Zp E e5 ;+ e5 [ , le test 112 est positif. Dans une étape 113, le calculateur détermine si deux tests 110 et 112 sensiblement concomitants, par exemple effectués dans un intervalle de temps de l'ordre de la seconde, sont positifs. Si c'est le cas, les mesures des capteurs de pression sont jugées fiables. Le calculateur 3 lève donc l'état de crédibilité forte et active l'état de crédibilité optimale 74, décrit ci-dessus. Dans le cas contraire, l'état de crédibilité forte 84 est maintenu. Ainsi, lors d'un changement de configuration ou au décollage, une fois l'état de crédibilité faible levé et l'état de crédibilité forte établi, l'état de crédibilité optimale peut être atteint si les tests 110 et 112 s'avèrent positifs. Ces deux tests 110 et 112 permettent ainsi de réévaluer la crédibilité des mesures de pression lors du vol, lorsque ces mesures ont été jugées peu crédibles, et de fournir au pilote des informations issues de ces mesures lorsque les capteurs de pression fonctionnent à nouveau. L'utilisation de deux tests complémentaires garantit par ailleurs une meilleure sécurité. On a représenté sur la Figure 6 un dispositif d'affichage auxiliaire 14 selon un mode de réalisation de l'invention. Comme précédemment décrit, le dispositif d'affichage auxiliaire 14 est relié au système 2 et en particulier au calculateur 3 de l'aéronef 1, et reçoit des informations et des ordres de commande de celui-ci. Le dispositif 14 d'affichage auxiliaire comprend des moyens 15, propres à afficher des évaluations des valeurs de vitesse, d'altitude, d'incidence et de pente air, choisies en fonction d'estimations de la fiabilité des capteurs d'incidence et de pression et transmises par le calculateur 3. Ces moyens 15 forment un dispositif de visualisation. Le dispositif d'affichage auxiliaire 14 affiche, dans l'état de crédibilité optimale, des caractéristiques de vol dites « principales » de l'aéronef, telles qu'une vitesse, une altitude et un nombre de Mach, obtenues à partir de mesures des capteurs de pression statique 5a, de pression totale 5b et d'incidence 5d. Lorsque les mesures d'au moins un de ces capteurs sont jugées non fiables, chaque caractéristique de vol « principale » obtenue à partir de mesures de ce capteur est remplacée sur le dispositif d'affichage par une caractéristique de vol « auxiliaire », homologue à la caractéristique principale remplacée, i.e. représentative de la même caractéristique de vol, et indépendante de toute mesure du capteur jugé non fiable effectuée lorsque les mesures de ce capteur sont jugées non fiables. Ainsi, la vitesse indiquée IAS est remplacée par une vitesse (l'équivalent de vitesse forte EV,* ou faible EV- ), l'altitude pression Zp est remplacée par une altitude (l'altitude pression GPS Zp** ou l'altitude au-dessus du géoïde de référence ZGPS), et le nombre de Mach Ma est remplacé par un nombre de Mach (nombre de Mach fort Mn* ou faible M**). Les Figure 6, 7 et 8 représentent ainsi des informations projetées sur le dispositif 15, affichées sous la forme de symboles, lors de la mise en oeuvre du procédé précédemment décrit, lorsque l'aéronef 1 est dans l'état de crédibilité optimale 74, dans l'état de crédibilité forte 84 et dans l'état de crédibilité faible 76 respectivement. Ces symboles comportent un indicateur 122 de commande de vol, affichant un symbole 124 de maquette de l'aéronef, occupant une position constante, qui matérialise une projection à l'infini de l'axe longitudinal X de l'aéronef 1, et une ligne d'horizon artificiel 126, au centre d'une échelle graduée de pente 128. La position de la ligne 126 par rapport au symbole 124 représente l'assiette Ode l'aéronef 1, cette assiette étant indiquée sur l'échelle graduée de pente 128 en regard de la ligne d'horizon artificiel 126. L'indicateur 122 comprend également un symbole 130 de vecteur vitesse, indiquant la direction du vecteur vitesse de l'aéronef par rapport à l'air. Ce symbole 130 est disposé le long de l'échelle graduée de pente 128. L'écart sur l'échelle graduée de pente 128 entre le symbole 124 de maquette et le symbole 130 de vecteur vitesse est égal à l'incidence a de l'aéronef 1. Avantageusement, les symboles de maquette 124 et de vitesse 130 ne sont affichés que dans les états de crédibilité optimale 74 ou forte 84. La pente air et l'incidence de l'aéronef 1 ne sont ainsi pas affichées dans l'état de crédibilité faible 76. Le dispositif 15 affiche par ailleurs un indicateur de vitesse 132. Cet indicateur 132 comprend une échelle de vitesse 134 graduée, représentée sous la forme d'un segment s'étendant entre deux points fixes, et un symbole de vitesse 136, disposé en regard de l'échelle de vitesse 134. Le symbole de vitesse 136 indique une vitesse de l'aéronef 1. Le symbole 136 a par exemple une forme de chevron. Un symbole 138, indiquant sous forme numérique la valeur de la vitesse, est ajouté en regard du symbole de vitesse 136. La vitesse indiquée dépend de l'état de crédibilité établi par le calculateur 3. La vitesse affichée au moyen du symbole de vitesse 136 est par exemple la vitesse indiquée IAS dans l'état de crédibilité optimale 74, l'équivalent de vitesse forte EV; dans l'état de crédibilité forte 84 et l'équivalent de vitesse faible EV- dans l'état de crédibilité faible 76. Un symbole 141, disposé en bas de l'échelle de vitesse 134 graduée, indique un nombre de Mach de l'aéronef. Le type de nombre de Mach affiché dépend de l'état de crédibilité courant. Dans l'état de crédibilité optimale 74, il s'agit du nombre de Mach Ma issu du machmètre 8, dans l'état de crédibilité forte 84, le nombre de Mach fort Mn* est affiché, et dans l'état de crédibilité faible 76, le nombre de Mach faible M** est affiché. Les symboles 138 et 141 sont munis d'une identification représentative du type de vitesse et de nombre de Mach affiché, avantageusement un code couleur, de manière à informer le pilote du type de vitesse dont il dispose. Le code couleur choisi indique de manière intuitive quels types de vitesse et de nombre de Mach sont affichés. Par exemple, les symboles 138 et 141 sont verts encadrés de vert dans l'état de crédibilité optimale 74, i.e. quand la vitesse indiquée IAS et le nombre de Mach Ma sont affichés, verts encadrés de jaune dans l'état de crédibilité forte 84, i.e. quand l'équivalent de vitesse forte EV,* et le nombre de Mach fort Mn* sont affichés, et jaunes encadrés de jaune dans l'état de crédibilité faible 76, i.e. quand l'équivalent de vitesse faible EV- et le nombre de Mach faible M** sont affichés. Comme représenté sur la Figure 6, dans l'état de crédibilité optimale, le dispositif 15 affiche en outre, superposé à l'indicateur de vitesse 132, un ruban 139a indiquant la vitesse indiquée IASa de décrochage de l'aéronef 1. Le ruban 139a est par exemple de couleur rouge. Dans l'état de crédibilité forte 84, un ruban 139a identique indique l'équivalent de vitesse de décrochage fort EV; . Dans l'état de crédibilité faible 76, le ruban 139a est remplacé par un ruban rouge discontinu 139b, comme représenté sur la Figure 8. Le ruban 139b indique l'équivalent de vitesse de décrochage faible EV; décrit ci-dessus. Le dispositif 15 affiche également un indicateur d'altitude 140. Cet indicateur 140 comprend une échelle d'altitude 142 graduée, représentée sous la forme d'un segment s'étendant entre deux points fixes, et un symbole d'altitude 144, disposé en regard de l'échelle d'altitude 142. Le symbole d'altitude 144 indique une altitude de l'aéronef 1. Le symbole 144 a par exemple une forme de chevron. Un symbole 146, indiquant sous forme numérique la valeur de l'altitude, est ajouté en regard du symbole d'altitude 144. L'altitude indiquée dépend de l'état de crédibilité établi par le calculateur 3. Dans l'état de crédibilité optimale 74, l'altitude affichée au moyen du symbole d'altitude 144 est par exemple l'altitude pression Zp ou, au choix du pilote, une altitude corrigée du calage, notée 4. Le calage est par exemple un calage QNH ou QFE. Le QFE est un code international destiné au calage de l'altimètre par rapport à un terrain donné de telle sorte que celui-ci indique une altitude nulle lorsque l'aéronef se trouve au sol sur ce terrain. Dans l'état de crédibilité forte 84 ou faible 76, l'altitude affichée au moyen du symbole d'altitude 144 est l'altitude pression GPS Zp** ou, au choix du pilote, l'altitude au-dessus du géoïde de référence ZGps, donnée par le capteur GPS. Ces choix entre Zp ou Zc et Zp** ou ZGps peuvent être effectués au moyen de l'interface de saisie 17, par exemple par actionnement d'un bouton dédié. Le dispositif 15 indique alors, dans un emplacement réservé 150, sur quelle référence l'altitude est calée. Par exemple, dans l'état de crédibilité optimale 74, si l'altitude affichée est l'altitude pression Zp, le calage « STD » est affiché dans l'emplacement 150, comme représenté sur la Figure 6. Si l'altitude affichée est une altitude corrigée du calage 4, « QNH » ou « QFE » est affiché dans l'emplacement 150, selon le calage choisi. Dans l'état de crédibilité forte 84 ou faible 76, si l'altitude pression GPS Zp** est affichée, le calage « STD » est affiché dans l'emplacement 150, comme représenté sur les Figure 7. Si l'altitude ZGps au-dessus du géoïde de référence est affichée, le calage «GEO » est affiché dans l'emplacement 150, comme représenté sur la Figure 8. Le symbole 146 est muni d'une identification représentative du type d'altitude affichée, avantageusement un code couleur. Le code couleur choisi indique de manière intuitive quel type de vitesse est affiché. Par exemple, le symbole 146 est vert encadré de vert dans l'état de crédibilité optimale 74, i.e. quand l'altitude Zp ou 4 est affichée, vert encadré de jaune lorsque l'hybridation lâche de Zp** sur Zp est en cours, et jaune encadré de jaune lorsque l'hybridation lâche de Zp** par Zp est arrêtée. Le symbole 146 encadré de jaune signifie donc que l'altitude pression GPS Zp** ou l'altitude au-dessus du géoïde de référence ZGps est affichée, et la couleur du symbole lui-même indique si le processus d'hybridation lâche de Zp** sur Zp est en cours ou non. Le dispositif d'affichage auxiliaire 14 est complété en outre par des moyens 16 d'affichage de messages d'alerte, comme décrit précédemment. Ces moyens 16 comprennent par exemple une fenêtre d'affichage 152, propre à afficher des messages textuels destinés à informer le pilote de l'état de crédibilité affecté aux mesures des capteurs, et de la fiabilité des mesures affichées par les instruments de bord classiques. Cette fenêtre 152 est par exemple intégrée à un système d'alerte de l'équipage (ou CAS pour Crew Alerting System). Par exemple, dans l'état de crédibilité forte, la fenêtre 152 affiche un message indiquant au pilote que les mesures des capteurs de pression 5a, 5b ne sont pas fiables, et un message indiquant que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Ces messages sont par exemple du type « STATIC &/or total pressure input fail » et « M/IAS unreliable (caution) » respectivement, comme illustré sur la Figure 7. Dans l'état de crédibilité faible 76, la fenêtre 152 affiche un message indiquant au pilote que les mesures des capteurs de pression 5a, 5b et d'incidence 5d ne sont pas fiables, et un message indiquant que les vitesses affichées par les instruments de bord classiques sont probablement incorrectes. Ces messages sont par exemple du type « AOA input fail » et « M/IAS unreliable (caution) » respectivement, comme illustré sur la Figure 8.Le dispositif d'affichage auxiliaire selon l'invention permet donc de fournir au pilote, à tout moment, des informations de pilotage fiables. En effet, lorsque les mesures d'au moins un capteur parmi les capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d) sont jugées non fiables, chaque caractéristique de vol déterminée à partir de mesures de ces capteurs est remplacée sur le dispositif d'affichage auxiliaire par une caractéristique de vol auxiliaire. Chaque caractéristique de vol auxiliaire est indépendante de toute mesure de capteur effectuée à un instant auquel ce capteur est jugé non fiable. Les caractéristiques de vol auxiliaires peuvent cependant dépendre de mesures du ou de chaque capteur jugé non fiable effectuées préalablement, lorsque ce ou ces capteurs étaient jugé(s) fiables. Le dispositif d'affichage auxiliaire selon l'invention permet également de prévenir l'équipage en cas de suspicion de défaillance d'un ou plusieurs capteurs. Le procédé et le système selon l'invention permettent ainsi de fournir à l'équipage d'un aéronef une évaluation de la fiabilité des caractéristiques de vol issues de mesures de capteurs de l'aéronef, et de l'alerter en cas de possible défaillance. En outre, ce procédé et ce système permettent de réévaluer en permanence la fiabilité des capteurs pour détecter une éventuelle remise en état de fonctionnement de capteurs jugés non fiables ou au contraire une défaillance d'un capteur auparavant jugé fiable. Le procédé et le système selon l'invention permettent également de fournir à l'équipage des caractéristiques de vol auxiliaires indépendantes des mesures des capteurs jugés non fiables. Ainsi, même en cas de défaillance de capteurs de pression et/ou d'incidence, l'équipage dispose de valeurs fiables des caractéristiques de vol. Il devra être compris que les exemples de réalisation présentés ci-dessus ne sont pas limitatifs. Notamment, la vérification de l'équation de sustentation lors des étapes 70 et 78 peut être réalisée selon une autre forme que celle précédemment décrite. Notamment, elle peut faire intervenir un autre type d'information de vitesse que le nombre de Mach Ma. En outre, les caractéristiques de vol auxiliaires peuvent être déterminées selon d'autres expressions, et éventuellement à partir d'autres types de capteurs qu'un capteur GPS, par exemple un autre type de capteur de position par satellites, un radar Doppler, ou une centrale à inertie. Le coefficient de portance peut également être projeté sur un autre axe que l'axe Z représenté, par exemple selon un axe perpendiculaire au vecteur vitesse. Par ailleurs, le dispositif d'affichage auxiliaire peut afficher d'autres informations auxiliaires, telles que la température statique ou l'estimateur de vent horizontal. Selon une variante, l'aéronef 1 comprend plusieurs capteurs du même type redondants, par exemple plusieurs capteurs de pressions statique et totale et/ou plusieurs capteurs d'incidence. Selon cette variante, le calculateur teste par exemple la fiabilité de chacun des capteurs d'un même type séparément. Ainsi, si un seul capteur d'un type donné est jugé fiable, le calculateur 3 utilise les mesures issues de ce capteur, et avertit le pilote d'une défaillance du ou des autre(s) capteur(s) du même type. Si tous les capteurs d'un même type sont non fiables, le calculateur 3 active un état de crédibilité forte ou faible suivant les cas, conformément au procédé décrit ci-dessus. Bien entendu, d'autres modes de réalisation peuvent être envisagés. En outre, les caractéristiques techniques des modes de réalisation et variantes mentionnées ci-dessus peuvent être combinées entre elles. REVENDICATIONS1.- Procédé de détermination d'un état de crédibilité de mesures de capteurs (5) d'un aéronef (1) lors d'un vol dudit aéronef (1), lesdits capteurs (5) comprenant des capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d), caractérisé en ce qu'il comprend les étapes suivantes : - déterminer (31) une information de vitesse (Ma) de l'aéronef (1), à partir de mesures de pressions statique Ps et totale PT desdits capteurs de pression statique (5a) et de pression totale (5b), - déterminer (33) un premier coefficient de portance az de l'aéronef (1) à partir d'une mesure (31) d'incidence a dudit capteur d'incidence (5d) et de ladite information de vitesse (Ma), - déterminer (39) une masse estimée/ri de l'aéronef (1), - déterminer (70) si une équation de sustentation (T1) de l'aéronef (1) est vérifiée, à partir de ladite masse estimée m , d'une information représentative du facteur de charge (nz) de l'aéronef (1), de ladite information de vitesse (Ma) et dudit premier coefficient de portance (Ez ) de l'aéronef (1), avec une première tolérance prédéterminée (e1), - activer un état de crédibilité optimale (74), dans lequel les mesures desdits capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d) sont jugées fiables, si ladite équation de sustentation est vérifiée avec ladite première tolérance (e1), - activer (75, 83) au moins un état de crédibilité non optimale (76; 84), dans lequel les mesures d'au moins un capteur parmi lesdits capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d) sont jugées non fiables, si ladite équation de sustentation n'est pas vérifiée avec ladite première tolérance (e1). 2.- Procédé de détermination selon la revendication 1, caractérisé en ce que ladite information de vitesse est un nombre de Mach Ma, et en ce que l'étape (70) consistant à déterminer si l'équation de sustentation (T1) de l'aéronef (1) est vérifiée comprend la vérification de l'égalité 0,7SPsMa2az = n z ing avec la première tolérance déterminée (e1), où S représente une surface de référence de l'aéronef (1), ledit premier coefficient de portance az est une projection selon un axe Z choisi, avantageusement un axe Z perpendiculaire à l'axe longitudinal de l'aéronef (1) défini notamment comme un axe moyen de poussée des moteurs, d'un coefficient de portance (Ci) de l'aéronef (1) selon un axe orthogonal au vecteur vitesse de l'aéronef (1) et d'un coefficient de traînée (Cx)parallèle à ce vecteur vitesse, nz désigne ladite information représentative du facteur de charge de l'aéronef (1) selon l'axe Z, et g l'accélération de la pesanteur. 3.- Procédé de détermination selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que ledit premier coefficient de portance (6f) est déterminé à partir d'une table de correspondance fournissant une estimation de la valeur dudit coefficient de portance (5g ) à partir de ladite information de vitesse (Ma), de la mesure d'incidence a et d'une configuration de vol de l'aéronef (1). 4.- Procédé de détermination selon l'une quelconque des revendications 1 à 3, caractérisé en ce que la détermination (39) de la masse estimée (e7) de l'aéronef (1) comprend les étapes suivantes : - détermination (21) d'une première masse (mp10) de l'aéronef (1), à partir de caractéristiques de l'aéronef (1) déterminées avant ou pendant le décollage de l'aéronef (1), - détermination (41) d'une deuxième masse (map) de l'aéronef (1) lors du vol de l'aéronef (1), à partir d'une équation de sustentation de l'aéronef (1) exprimant la masse de l'aéronef (1) en fonction d'une information représentative du facteur de charge (ni) de l'aéronef (1), d'un deuxième coefficient de portance (5f) de l'aéronef (1), d'une information de vitesse (Ma) de l'aéronef (1) et d'une pression statique (Ps) d'une masse d'air traversée par l'aéronef (1), déterminées à partir de mesures de capteurs (5a, 5b, 5c, 5d) de l'aéronef (1) lors du vol de l'aéronef (1), - évaluation (39) de la masse estimée (dl), à partir d'une masse calculée tenant compte desdites première et deuxième masses, par exemple à partir de la moyenne pondérée desdites première et deuxième masses.. 5.- Procédé de détermination selon la revendication 4, caractérisé en ce que lesdites première (mp10) et deuxième (msp) masses sont des masses de l'aéronef (1) à un instant initial (t1) au parking, et en ce que l'étape (39) d'évaluation de ladite masse estimée (tri) comprend la détermination de ladite masse calculée à partir desdites première et deuxième masses, la détermination d'une masse (FU) de carburant consommé depuis ledit instant initial (t1), ladite masse estimée (m) étant évaluée en soustrayant ladite masse (FU) de carburant consommé de ladite masse calculée. 6.- Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend une étape (47) de détermination d'une altitude auxiliaire (Zp**) de l'aéronef (1) à partir d'au moins une mesure d'altitude géographique obtenue indépendamment d'un capteur de pression, notamment issue d'un capteur de positionnement par satellite, et une étape de mise à disposition d'un équipagede l'aéronef (1) de ladite altitude (Zp**) auxiliaire, au moins lorsqu'un état de crédibilité non optimale est activé. 7.- Procédé de détermination selon la revendication 6, caractérisé en ce que ladite altitude auxiliaire (Zp**) est déterminée à partir d'une altitude ZGps déterminée indépendamment d'un capteur de pression de l'aéronef (1), sous la forme : Zp** (0= ZGPS (t) AZ (j) où AZ**(j) est un terme correctif ajusté périodiquement, destiné à corriger un décalage entre l'atmosphère standard et l'atmosphère réelle de vol de l'aéronef (1), ledit terme correctif prenant une valeur figée lorsqu'un état de crédibilité non optimale est activé. 8.- Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend, au moins lorsque l'état de crédibilité non optimale (76; 84) est activé, une étape (49, 76) de détermination d'au moins une information de vitesse auxiliaire (TAS,* , M;, ,TAS** , M* * ) à partir d'au moins une mesure d'un capteur (10) de position géographique de l'aéronef (1), et une étape de mise à disposition d'un équipage de l'aéronef de ladite information de vitesse auxiliaire (TAS;,, M;; , TAS**, M**). 9.- Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce que l'état de crédibilité non optimale (76; 84) est choisi parmi : - un état de crédibilité forte (84), dans lequel les mesures d'au moins un capteur parmi les capteurs de pressions statique (5a) et de pression totale (5b) sont jugées non fiables et les mesures dudit capteur d'incidence (5d) sont jugées fiables, et - un état de crédibilité faible (76) dans lequel les mesures du capteur d'incidence (5d) sont jugées non fiables. 10.- Procédé de détermination selon la revendication 9, caractérisé en ce que, au moins dans l'état de crédibilité forte, le procédé comporte une étape (49) de détermination d'une information de vitesse auxiliaire (Mn*) de crédibilité forte, sur la base de la mesure d'incidence (a) du capteur d'incidence (5d), en fonction d'un estimateur de pression statique (Ps**) déterminé à partir d'au moins une mesure (ZGps) obtenue indépendamment d'un capteur (10) de pression de l'aéronef (1). 11.- Procédé de détermination selon la revendication 10, caractérisé en ce que ladite information de vitesse auxiliaire (M'*) de crédibilité forte est déterminée périodiquement à chaque instant à partir d'une information de vitesse auxiliaire decrédibilité forte (Mn*) obtenue à un instant précédent, avantageusement à l'aide d'une série par récurrence. 12.- Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend une estimation (51) de la vitesse (W; , W ) du vent, ladite vitesse du vent (1/n) étant déterminée périodiquement à partir d'informations de vitesses issues d'un capteur (10) de position géographique. 13.- Procédé de détermination selon la revendication 12, caractérisé en ce que, au moins dans l'état de crédibilité faible, le procédé comporte une étape (76) de détermination d'une information de vitesse auxiliaire (TAS** , M** ) de crédibilité fiable, comprenant les étapes suivantes : - détermination d'informations de vitesse (GS,Vz(GPS)) de l'aéronef (1) par rapport au sol à partir de mesures d'un capteur (10) de position géographique de l'aéronef (1), - estimation (51) d'une vitesse (W ) du vent, - détermination de ladite vitesse auxiliaire (TAS* * ,M ) de crédibilité faible par rapport à l'air, à partir desdites informations de vitesse (GS,Vz(GPS)) de l'aéronef (1) par rapport au sol et de ladite estimation de la vitesse (W ) du vent. 14.- Procédé de détermination selon l'une quelconque des revendications 9 à 13, caractérisé en ce qu'il comprend une étape de détermination d'un état de fiabilité de mesures dudit capteur d'incidence (5d), comprenant les étapes suivantes : - déterminer un estimateur de pression statique (Ps**) à partir d'au moins une mesure (ZGps) obtenue indépendamment d'un capteur (10) de pression de l'aéronef, - déterminer (76) une information de vitesse auxiliaire (M**), sur la base dudit estimateur de pression statique (Ps**), - déterminer un deuxième coefficient de portance (Cz ) de l'aéronef (1), à partir d'une mesure (31) d'incidence a dudit capteur d'incidence (5d) et de ladite information de vitesse auxiliaire (M**), - déterminer (78) si l'équation de sustentation de l'aéronef (1) est vérifiée, à partir de ladite masse estimée m , de ladite information de vitesse auxiliaire ( M**) et dudit deuxième coefficient de portance (a z) de l'aéronef (1), avec une deuxième tolérance prédéterminée (e2),- confirmer (75) l'état de état de crédibilité faible (76), dans lequel les mesures du capteur d'incidence (5d) sont jugées non fiables, si ladite équation de sustentation n'est pas vérifiée avec ladite deuxième tolérance prédéterminée, - activer (83) l'état de crédibilité forte (84), dans lequel les mesures d'au moins un capteur parmi les capteurs de pression statique (5a) et de pression totale (5b) sont jugées non fiables et les mesures dudit capteur d'incidence (5d) sont jugées fiables, si ladite équation de sustentation est vérifiée avec ladite tolérance prédéterminée (e2). 15.- Procédé de détermination selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend, dans l'état de crédibilité non optimale (76 ; 84), une étape d'affichage d'un message d'alerte à destination d'un équipage de l'aéronef (1), indiquant que les mesures d'au moins un capteur choisi parmi lesdits capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d) ne sont pas fiables. 16.- Système (2) de détermination d'un état de crédibilité de mesures de capteurs d'un aéronef (1) lors d'un vol dudit aéronef (1), lesdits capteurs comprenant des capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d), caractérisé en ce qu'il comprend : - des moyens configurés pour déterminer une information de vitesse (Ma) de l'aéronef (1), à partir de mesures de pressions statique Ps et totale PT desdits capteurs de pression statique (5a) et de pression totale (5b), - des moyens configurés pour déterminer un premier coefficient de portance (az ) de l'aéronef (1) à partir d'une mesure (31) d'incidence a dudit capteur d'incidence (5d) et de ladite information de vitesse (Ma), - des moyens configurés pour déterminer une masse estimée t'il de l'aéronef (1), - des moyens configurés pour déterminer si une équation de sustentation de l'aéronef (1) est vérifiée, à partir de ladite masse estimée m , d'une information représentative du facteur de charge (nz), de ladite information de vitesse (Ma) et dudit coefficient de portance (6z ) de l'aéronef (1), avec une première tolérance prédéterminée (Ci), - des moyens configurés pour activer un état de crédibilité optimale (74), dans lequel les mesures desdits capteurs de pression statique (5a), de pression totale (5b) et d'incidence (5d) sont jugées fiables, si ladite équation de sustentation est vérifiée avec ladite première tolérance (e1), - des moyens configurés pour activer au moins un état de crédibilité non optimale (76; 84), dans lequel les mesures d'au moins un capteur parmi lesdits capteurs depression statique (5a), de pression totale (5b) et d'incidence (5d) sont jugées non fiables, si ladite équation de sustentation n'est pas vérifiée avec ladite première tolérance (c1).5