Brevet FR2648233A1 — publié le 19901214

La présente invention a tout d'abord pour objet un procédé de détermination de la vitesse par rapport à l'air d'un hélicoptère en vue de son pilotage, hélicop- tère d'axe longitudinal OX, d'axe transversal OY et d'axe vertical OZ, pourvu d'un rotor de sustentation dont le pas collectif P0, le pas cyclique longitudinal PX et le pas cyclique transversal Py sont commandables, ainsi que d'un rotor anticouple dont le pas Pz est com- mandable, procédé dans lequel on détermine les deux composantes VX et Vy de cette vitesse sur les axes OX et OY, respectivement, en fonction des valeurs commandées desdits pas P0, PX, Py et PZ, des valeurs des angles d'inclinaison OX et y de l'hélicoptère autour des axes Èx y OY et OX, respectivement, et de la valeur d'accélération Z de l'hélicoptère selon l'axe OZ. La vitesse ainsi déterminée est affichée sur le tableau de bord de l'hélicoptère, à l'usage du pilote. En effet, la vitesse de l'hélicoptère par rapport à l'air, souvent appelée vitesse-air, est la vitesse propre de l'hélicop- tère par rapport à la masse d'air dans laquelle il se déplace. En l'absence de vent, la vitesse air est la même que la vitesse par rapport au sol, souvent appelée vitesse- sol. En présence de vent, la vitesse-sol est égale à la somme de la vitesse-air et de la vitesse du vent. La vitesse-air est importante pour des raisons de sécurité, car c'est elle qui conditionne, en vol, l'équilibre de l'hélicoptère. I1 est donc primordial que le pilote con- naisse à tout instant la vitesse-air avec la meilleure précision possible. Il existe des anémomètres qui permettent de mesurer cette vitesse. Toutefois les anémomètres classiques, qui repo- sent sur une mesure de pression à l'aide d'un tube de Pitot, par exemple, ne donnent des résultats valables que si la vitesse est supérieure à un certain seuil, couram- ment de l'ordre de 45 noeuds, par exemple. Les anémome- tres à laser ne souffrent pas de cet inconvénient, mais ils sont d'un prix de revient et d'un poids qui rend leur emploi difficile à bord d'un hélicoptère. C'est pourquoi un procédé du type défini ci-dessus, qui ne nécessite pas, pour sa mise en oeuvre, l'embarquement d'un appareillage lourd et spécifique, puisqu'il utilise des grandeurs nor- malement disponibles à bord d'un hélicoptère, est appré- ciable. On connait déjà un procédé du type défini ci-dessus, par la demande française N 2 613 078, au nom de la demande- resse. Ce procédé est basé sur une modélisation théorique du fonctionnement de l'hélicoptère, modélisation qui permet d'établir des formules qui relient les composantes VXet Vy de la vitesse aux pas P0, PX, Py et PZ' aux angles d'inclinaison ex et ey et à l'accélération yZ. Cependant, les valeurs de certains des paramètres interve- nant dans ces formules étant difficilement accessibles, on prévoit d'effectuer, pour chaque cas pratique, au moins un vol d'étalonnage permettant une détermination expérimen- tale des valeurs des coefficients qui interviennent dans les formules théoriques. A cet effet, on fait passer l'hélicoptère par un certain nombre de configurations de vol différentes, pour lesquelles on se donne les moyens d'estimer la vitesse-air, à l'aide de mesures directes ou indirectes, et pour lesquelles on mesure Pot PX' PY' PZ' ex. et Z. En identifiant la valeur de la vitesse-air estimée et celle fournie par la formule théorique de cal- cul pour chaque configuration de vol, on obtient un système d'équations dont les inconnues sont les coefficients que l'on veut déterminer. On résout ce système par la méthode des moindres carrés, par exemple. Ainsi, dans le procédé connu, la forme, ou structure, des formules qui permettent le calcul de la vitesse-air sont établies de façon théo- rique, tandis que les valeurs des coefficients qui interviennent dans ces formules sont déterminées de façon expérimentale. Toutefois, si ce procédé connu donne, dans certains cas, des résultats satisfaisants, il n'en est pas toujours ainsi, vraisemblablement parce que le modèle théorique du fonctionnement qui lui sert de base ne correspond pas, malgré sa complexité, à la réalité. Il serait sans doute possible d'affiner ce modèle théorique pour un type d'hélicoptère bien déterminé, mais dans ce cas, les cal- culs effectués pour ce type d'hélicoptère ne seraient pas valables pour un autre type d'hélicoptère, ce qui rend le procédé connu long et coûteux. La présente invention vise à pallier ces inconvénients en procurant un procédé qui donne, quel que soit le type d'hélicoptère, des résultats toujours précis, sans qu'il soit nécessaire d'effectuer, pour chaque type différent d'hélicoptère, une modélisation théorique différente. A cet effet, la présente invention a pour objet un procédé du type défini ci-dessus, caractérisé par le fait que: - au cours d'une première étape, on effectue un vol d'éta- lonnage faisant passer l'hélicoptère par une pluralité de N configurations de vol différentes, pour chacune des- quelles on mesure, d'une part, une pluralité de grandeurs permettant d'estimer les vitesses V et Vy, et d'autre part les valeurs des pas Po Px, PY et PZ' des angles &X et y et de l'accélération yZ, - au cours d'une deuxième étape, on estime lesdites vites- ses VX et Vy lors de chaque configuration de vol, on calcule une première matrice des coefficients de corréla- tion de tous les couples possibles de variables d'un premier ensemble de variables comprenant la vitesse estimée /VX et des variables construites à partir des pas Po PX' Py et PZ, des angles eX et ey et de l'accéléra- tion yZ, et on calcule une deuxième matrice de façon comparable à la première mais relative à un deuxième ensemble de variables comprenant la vitesse estimée Vy et lesdites variables construites, - au cours d'une troisième étape, on détermine, à l'aide de ladite première matrice, la variable construite du premier ensemble la plus correlée à la vitesse estimée Vx, A on soustrait de cette vitesse estimée VX ladite variable construite la plus correlée multipliée par un coefficient multiplicatif choisi pour que la vitesse estimée ainsi modifiée soit totalement decorrélée de ladite variable construite la plus corrélée, on calcule l'erreur commise en admettant que la valeur modifiée est nulle, on calcule une première matrice modifiée comme on a calculé ladite première matrice mais en remplaçant la vitesse estimée par la vitesse estimée modifiée, et on recommence ainsi, Jusqu'à avoir soustrait de la vitesse estimée 'x un nombre nx de variables construites du premier ensemble, suffisant pour que l'erreur commise, en admettant que la valeur modifiée est nulle, soit inférieure à un seuil, de façon à disposer, pour la valeur de la vitesse estimée VX, d'une formule de calcul en fonction de n variables construites x du premier ensemble, et on procède de même pour la deuxiè- me matrice, de façon à disposer d'une formule semblable pour y, en fonction de n variables construites du deu- y y xième ensemble, et, - au cours d'une quatrième étape, en vol normal, on mesure les valeurs des pas Po, PX, Py et PZ, des angles X et y et de l'accélération yz, on calcule les vitesses Vx et Vyen appliquant les formules déterminées au cours de la troisième étape, et en réponse aux valeurs mesurées au cours de la quatrième étape, et on affiche, à l'usage du pilote de l'hélicoptère, les valeurs VX et Vy calculées. Dans le procédé de la présente invention, les formules qui permettent de calculer V et V en vol normal sont x y établies de façon expérimentale, en tenant compte des va- riables qui sont les plus corrélées aux vitesses estimées VX et Vy sur l'ensemble des configurations de vol du vol d'étalonnage. Il en résulte que, quel que soit le type et les particularités de l'hélicoptère pour lequel on applique le procédé, les formules de calcul de Vx et Vy collent à la réalité sans qu'il soit nécessaire de changer quoi que ce soit au procédé. Ce résultat remarquable est atteint par la demanderesse notamment car elle applique au cours de la troisième étape et indépendamment à chacune des matrices de coefficients de corrélation, un algorithme qui n'est autre que l'algorithme connu sous le nom de régression pas à pas (en anglais "stepwise regression"). De façon surprenante, les résultats obtenus sont très sa- tisfaisants, alors que les variables considérées ne sont pas indépendantes. Par ailleurs, avec le procédé de l'in- vention, les formules de calcul de VX et de Vy ne risquent pas d'tétre inutilement complexes, puisque l'on ne tient compte que des variables les plus corrélées avec les composantes VX et Vy. Cette caractéristique permet, en vol normal, d'avoir un temps de calcul ramené à ce qui est juste nécessaire. Avantageusement, les variables construites au cours de la deuxième étape à partir des pas Po, PX, Py et PZ, des angles &X et ey et de l'accélération yZ comprennent: - des variables simples du premier ordre, à savoir POP PX' Py PZ et yZ, - des variables combinées du premier ordre, à savoir P x- et P, et, - des variables simples du second ordre, à savoir les carrés et les produits deux à deux des variables Po PX' Py, PZ &X' Oy et YZ. Dans ce cas, l'expérience montre que l'erreur ne varie plus de façon significative au delà d'un nombre nx ou n de variables égal à 15, et parfois même bien avant que Y cette valeur ne soit atteinte. Le choix des variables combinées P -'X et Py- est justifié du fait que les x x modélisations théoriques, ainsi que l'expérience, mon- trent que ce terme est prépondérant. Il est donc avanta- geux, du point de vue de la durée des calculs effectués au cours de la troisième étape, de le prévoir parmi les variables construites à partir de Po0, Px, Py, PZ, eX' ey et Yz. Avantageusement encore, à la fin de la troisième étape, on considère comme des inconnues les n et n coefficients multiplicatifs des formules de calcul de Vx et de V, respectivement, on reporte dans ces formules les résultats des mesures effectuées au cours de la première étape, afin de former au moins un système d'équations linéaires dont les n et n coefficients multiplicatifs sont les incon- x y nues, on résout ce système et on corrige, dans les formu- les de calcul de Vx et de Vy, les valeurs des nx et ny coefficients multiplicatifs, respectivement, en les rem- plaçant par celles ainsi déterminées. Dans ce cas, l'application de l'algorithme de "stepwise regression" aux matrices de coefficients de corrélation, on matrice de corrélation, sert uniquement à déterminer la structure des formules, les valeurs numériques des coefficients employés dans ces formules étant détermi- nées par résolution d'un système d'équations obtenus en reportant, dans ces formules, les résultats des mesures effectuées au cours du vol d'étalonnage. Ceci permet en général d'améliorer la précision relative aux valeurs des coefficients employés dans les formules. De plus, ceci permet, lorsque par exemple on doit changer un capteur sur un hélicoptère dont les formules de calcul étaient connues, de faire un nouvel étalonnage de ces formules sans avoir à repasser par l'étape des matrices de corré- lation et de la régression pas à pas. En effet, lorsque l'on se contente de changer un capteur sur un hélicoptère, il n'y a aucune raison pour que la structure des formules de calcul de VX et de Vy en soit affectée, et seule la valeur numérique des coefficients qui interviennent dans ces calculs risque de se trouver changée. Dans un tel cas, il est suffisant de procéder à un étalonnage comme indi- qué ci-dessus, sans reprendre le processus permettant la détermination de la structure des formules. Dans une première mise en oeuvre du procédé de l'invention, on forme un premier et un deuxième systèmes d'équations linéaires, en identifiant, pour une pluralité de configu- rations de vol prises parmi les N configurations de vol, la vitesse estimée VX ou Vy obtenue au cours de la deuxième étape, à la quantité obtenue par la formule de calcul correspondante. Cette première mise en oeuvre est assez simple, mais néces- site, comme cela sera mieux compris dans la suite, de disposer de valeurs estimées VX et Vy assez proches de la réalité. Il peut en 8tre ainsi si on a pris soin d'embar- quer, pendant le vol d'étalonnage, un anémomètre laser qui permet de mesurer directement la vitesse-air, par exemple. Dans une deuxième mise en oeuvre du procédé de l'inven- tion, on mesure, au cours de la première étape, la vi- tesse par rapport au sol et les angles d'attitude de l'hélicoptère pour chacune des configurations de vol, ainsi que l'instant de datation de celle-ci, et, pour former le système d'équations linéaires, on choisit parmi les configurations de vol, des paires de configurations de vol d'instants de datation suffisamment proches pour admettre que le vent est le même dans ces deux configura- tions de vol, on transpose dans le repère terrestre, et pour chaque configuration de vol, le résultat de la mesure de la vitesse par rapport au sol et les formules de cal- cul des composantes de la vitesse par rapport à l'air, et on prend pour équations à résoudre celles obtenues en identifiant la variation mesurée de la vitesse par rapport au sol de l'hélicoptère entre les deux configurations d'une paire, et la variation calculée à l'aide des for- mules de calcul, de la vitesse par rapport à l'air de l'hélicoptère. En utilisant ce procédé, on n'utilise, pour le calcul des coefficients des formules de calcul, que les mesures de vitesse sol et d'attitude effectuées au cours du vol d'étalonnage, mesures dont le résultat est connu avec une bonne précision. Dans une troisième mise en oeuvre du procédé de l'invention, on mesure, au cours de la première étape, la vitesse par rapport au sol et les angles d'attitude de l'hélicoptère pour chacune des configurations de vol, ainsi que l'ins- tant de datation de celle-ci, et, pour former le système d'équations linéaires, on admet que le vent suit, dans le repère terrestre, une loi polynomiale en fonction du temps de coefficients inconnus,on transpose la vitesse du vent dans le repère hélicoptère pour chaque configuration de vol, et on identifie, pour une pluralité de configurations de vols prises parmi les N configurations de vol, le résul- tat de la mesure de la vitesse par rapport au sol à la somme de la vitesse du vent et de la vitesse par rapport à l'air obtenue par les formules de calcul de Vx et Vly, de façon à déterminer en même temps les n et n coeffi- x y cients multiplicatifs et les coefficients inconnus de la loi polynomiale du vent. Dans ce cas, le temps de calcul est un peu plus long, car le système comprend davantage d'inconnues, mais la préci- sion du résultat est accrue. Avantageusement encore, au cours de la première étape, on détermine les N configurations de vol différentes afi-n qu'un nombre minimal d'entre elles appartienne à chacun des sous-domaines du domaine des configurations possibles de l'hélicoptère en vol normal, chaque sous- domaine étant défini par un ensemble de plages de valeurs des grandeurs mesurées pour chaque configuration de vol. On assure ainsi une bonne répartition des configurations de vol pour lesquelles on effectue les mesures, afin que les matrices de corrélation soient bien représentatives de l'ensemble des phénomènes susceptibles de se produire en vol normal. Avantageusement encore, on mesure, au cours de la première étape, la vitesse par rapport au sol et les angles d'atti- tude de l'hélicoptère, ainsi que le module de la vitesse par rapport à l'air de l'hélicoptère, on admet que le vent moyen dans le repère terrestre est resté constant pendant le vol d'étalonnage, on forme une pluralité d'équations, chaque équation identifiant, pour une configuration de vol et dans le repère terrestre, le module mesuré de la vitesse par rapport à l'air et le module calculé de la différence entre la vitesse par rapport au sol et la vitesse du vent, on soustrait membre à membre, et deux à deux, ces équa- tions, on résout le système d'équations ainsi formé, pour calculer la vitesse du vent moyen pendant le vol d'étalonnage, et, au cours de la deuxième étape, on transpose la vitesse du vent moyen ainsi calculée dans le repère hélicoptère pour estimer la vitesse-air comme somme de la vitesse par rapport au sol et de la vitesse du vent moyen. Dans ce cas, il n'est pas nécessaire d'embarquer un ané- momètre laser, par exemple, pendant le vol d'étalonnage. En effet, avec le procédé ci-dessus il est possible d'es- timer les composantes VX et Vy de la vitesse-air au cours du vol d'étalonnage à partir de la seule connaissance du module de cette vitesse, des trois composantes de la vitesse sol et des angles d'attitude au cours de ce vol. Or le module de la vitesse-air est mesurable, au moins lorsqu'il est supérieur à un certain seuil, par le capteur anémobarométrique normalement prévu sur l'hélicoptère, et les composantes de la vitesse sol, l'attitude et le cap sont des grandeurs normalement disponibles à bord d'un hélicoptère. Le procédé ci-dessus repose sur l'hypothèse que le vent, mesuré dans le repère terrestre, reste le même pendant toute la durée du vol d'étalonnage. Les com- posantes du vent moyen sont alors calculées à partir de l'ensemble des configurations de vol pour lesquelles on dispose d'une mesure valable du module de la vitesse-air. Ensuite, pour estimer les composantes de la vitesse air dans chaque configuration de vol, on utilise les composan- tes de la vitesse sol mesurée dans cette configuration et celle du vent moyen déterminée comme ci-dessus. Avantageusement encore, on mesure, au cours de la quatrième Il étape en vol normal, à l'aide de moyens anémobarométri- ques, la vitesse de l'hélicoptère par rapport à l'air, et en réponse à la valeur de la vitesse anémobarométrique ainsi mesurée, on sélectionne, comme valeur de la vitesse par rapport à l'air, le résultat de la mesure, celui du calcul par les formules de calcul déterminées au cours de la troisième étape, ou une combinaison des deux. Dans ce cas, on utilise les résultats de la mesure anémo- barométrique lorsque celle-ci est valable, c'est-à-dire lorsque la vitesse air est suffisamment élevée, par exemple. La présente invention a également pour objet un dispositif pour la mise en oeuvre du procédé de l'invention, caracté- risé par le fait qu'il comprend: - des moyens pour mesurer, au cours d'une première étape au cours de laquelle est effectué un vol d'étalonnage faisant passer l'hélicoptère par une pluralité de N confi- gurations de vol différentes, d'une part une pluralité de grandeurs permettant d'estimer les vitesses Vx et Vy, et d'autre part les valeurs des pas Po PX' Py et PZ,' des angles eX' ty et de l'accélération Z - des moyens pour estimer, au cours d'une deuxième étape, lesdites vitesses VX et Vy lors de chaque configuration de vol, pour calculer une première matrice des coefficients de corrélation de tous les couples possibles de variables 3C d'un premier ensemble de variables comprenant la vitesse estimée V et des variables construites à partir des pas x Po PX, Py et PZ, des angles &X et &y et de l'accéléra- tion yZ, et pour calculer une deuxième matrice de façon comparable à la première mais relative à un deuxième ensemble de variables comprenant la vitesse estimée Vy et lesdites variables construites, - des moyens pour déterminer, au cours d'une troisième étape, et à l'aide de ladite première matrice, la variable construite du premier ensemble la plus corrélée à la vitesse estimée VX, pour soustraire de cette vitesse es- timée VX ladite variable construite la plus corrélée multipliée par un coefficient multiplicatif choisi pour que la vitesse estimée ainsi modifiée soit totalement decorrelée de ladite variable construite la plus correlée, pour calculer l'erreur commise en admettant que la valeur modifiée est nulle, pour calculer une première matrice modifiée comme a été calculée ladite première matrice mais en remplaçant la vitesse estimée VX par la vitesse esti- mée modifiée, et pour recommencer ainsi jusqu'à ce que soient soustraites de la vitesse estimée VX un nombre nx de variables construites du premier ensemble suffisant pour que l'erreur commise, en admettant que la valeur modifiée est nulle, soit inférieure à un seuil, de façon à disposer pour la valeur de la vitesse estimée VX d'une formule de calcul, en fonction de n variables construites du premier ensemble, et pour procéder de même pour la deuxième matrice de façon à disposer d'une formule sembla- ble pour V., en fonction de n variables construites du y deuxième ensemble, et, - des moyens pour mesurer, au cours d'une quatrième étape, en vol normal, les valeurs des pas Pot PX' Py et PZ, des angles eX et 6y et de l'accélération yZ, des moyens pour calculer les vitesses VX et Vy à l'aide des formules dé- terminées au cours de la troisième étape et en réponse aux valeurs mesurées au cours de la quatrième étape, et des moyens pour afficher, à l'usage du pilote de l'hélicoptère, les valeurs VX et Vy calculées. La présente invention sera mieux comprise à la lecture de la description suivante de la mise en oeuvre préférée du procédé de l'invention, et de la forme de réalisation préférée du dispositif de l'invention, et de quelques variantes, faite en se référant aux dessins annexés, sur lesquels, - la figure 1 représente schématiquement un hélicoptère en cours de vol, et, - la figure 2 représente un schéma par blocs du disposi- tif de l'invention, embarqué à bord de l'hélicoptère de la figure 1. En se référant à la figure 1, un hélicoptère 1 est repré- santé en cours de vol. L'hélicoptère i est animé, par rapport à la masse d'air dans laquelle il se déplace, d'une vitesse propre, ou vitesse air, représentée par un vecteur V. Il est animé, par rapport au sol, d'une vitesse représentée par un vecteur W qui est ici différente de sa vitesse propre T'du fait de l'existence d'un vent qui entralne la masse d'air dans laquelle se déplace l'hélicoptère 1 à une vitesse, par rapport au sol, représentée par un vecteur vent W. Le procédé et le dispositif de l'invention permettent la détermination, à bord de l'hélicoptère 1, et à partir des informations délivrées par ses instruments de bord, de sa vitesse propre, ou vitesse air V. De façon connue, un repère OXYZ est lié à l'hélicoptère 1, l'axe OX étant un axe longitudinal dirigé vers l'avant de l'hélicoptère 1, l'axe OY étant un axe transversal dirigé vers la droite de son pilote, et l'axe OZ un axe dirigé vers le bas, ici confondu avec l'axe du moyeu de son rotor. De même, le repère terrestre comprend ici un axe N diri- gé vers le Nord, un axe E dirigé vers l'Est, et un axe V dirigé verticalement vers le bas. Ici, on s'intéresse essentiellement aux deux composantes de la vitesse air selon les axes OX et OY, notées VX et Vy respectivement. De façon connue, l'hélicoptère 1 est pourvu d'un rotor de sustentation pourvu de pales dont le pas collectif P0, le pas cyclique longitudinal PX et le pas cyclique trans- versal Py sont commandables. Le pas ou angle de calage d'une pale est l'angle que fait la corde de portance nulle de son profil avec le plan de rotation des pales. Le pas collectif P0 représente l'angle de calage moyen de l'en- semble des pales, étant entendu qu'il est prévu de pouvoir faire varier de façon cyclique, au cours de chaque tour du rotor, l'angle de calage instantané de chaque pale autour de sa valeur moyenne. Les pas cycliques PX et Py corres- pondent aux variations cycliques de l'angle de calage. Les pas Po, PX et Py sont commandés par le pilote de l'hé- licoptère qui agit à cet effet sur la position et l'orien- tation d'un plateau cyclique qui commande lui-même le pas de chacune des pales. Une action sur le pas collectif P0 fait simplement monter ou descendre l'hélicoptère, tandis qu'une action sur les pas cycliques PX et Py fait basculer le disque rotor, c'est-à-dire le disque décrit par les pales et entratne un déplacement de l'hélicoptère dans le plan horizontal. Par ailleurs, l'hélicoptère est pourvu d'un rotor anticouple, dont le pas PZ est lui aussi commandable par le pilote de l'hélicoptère. Il est connu, notamment par la demande française N 2 613 078, de déterminer les composantes VX et Vy de la vitesse air grâce à des formules de calcul qui les expriment en fonction des valeurs des pas P0o PX Py et PZcommandées par le pilote, des valeurs des angles d'in- clinaisons eX et y de l'hélicoptère 1 autour des axes OY et OX respectivement, ainsi que de la valeur de l'accé- lération yZ de l'hélicoptère 1 selon l'axe OZ. Dans la mesure o les valeurs de P0o PX, PYI PZ ' e y et Z sont facilement disponibles à bord de l'hélicoptère, en sortie de capteurs de déplacement, inclinomètres et accélé- romètres appropriés, un tel procédé est de mise en oeuvre facile. Le procédé de l'invention prévoit, comme le procédé connu, d'établir des formules permettant le calcul de VX et de Vy à partir de la connaissance de P0o PX' PY PZ' XI y et IZ. Néanmoins, il diffère de ce procédé connu par le fait qu'il prévoit d'établir la structure de ces for- mules non pas à partir d'un modèle théorique du fonctionne- ment de l'hélicoptère, mais à partir de données expérimen- tales recueillies au cours d'un vol d'étalonnage, comme cela va maintenant être décrit. Le procédé de l'invention commence donc par une première étape au cours de laquelle un vol d'étalonnage est effec- tué. Au cours de ce vol, l'hélicoptère passe par une pluralité de N configurations de vol différentes. Pour chacune d'entre elles on mesure d'une part les valeurs de POP PX, PYs PZ, &X ey et Z et d'autre part une pluralité de grandeurs permettant d'estimer les composantes VX et Vy de la vitesse air. Dans le cas o l'on prévoit d'embarquer un anémomètre laser au cours du vol d'étalonnage, cette pluralité de grandeurs comprend simplement les valeurs des composantes VX et Vy directement fournies par l'anémomètre laser. Cependant, compte tenu du coût et du poids d'un tel anémomètre laser, il en est rarement ainsi. Dans ce cas, le procédé de l'invention prévoit de mesurer, en vue de l'estimation ultérieure des vitesses VX et Vy.., les gran- deurs suivantes: - les trois composantes S., Sy et SZ de la vitesse sol 7 dans le repère OXYZ, - les trois angles d'attitude, de l'hélicoptère 1, c'est- à-dire les angles de cap t, de tangage e, et de roulis et, - le module V de la vitesse air V. On notera ici que toutes ces grandeurs sont normalement disponibles à bord de l'hélicoptère 1, les composantes SX, Syet Sz en sortie d'un dispositif de mesure de la vitesse sol, les angles ', 9 et en sortie de la centrale d'attitude et le module V en sortie d'un capteur anémo- barométrique de type connu, comprenant par exemple un tube de Pitot. Toutefois, et comme cela a déjà été signalé, le résultat de la mesure faite par le capteur anémobarométri- que n'est valable que si le module V est supérieur à un seuil déterminé VS, valant par exemple 45 noeuds. Avant que ne soit expliqué comment, au cours de la deuxième étape du procédé de l'invention, on estime les composantes VX et Vy à partir de SX, S SZ, t ' 8, 'f et V, les critères de choix des configurations de vol pour lesquelles sont effectuées, au cours de la première étape, les mesures, sont maintenant exposés. Ces critères de choix obéissent à deux principes: - les configurations de vol retenues doivent représenter des configurations d'équilibre stable de l'hélicoptère 1, - la répartition des configurations de vol retenues doit être représentative de l'ensemble des conditions d'utili- sation de l'hélicoptère en vol normal. En conséquence, sur le nombre total de configurations de vol qui pourraient être retenues au cours du vol d'éta- lonnage, un certain nombre seulement le sont. Le nombre total de configurations qui pourraient être retenues est lié à la cadence de mesure des différents capteurs et appareils qui délivrent les grandeurs mesurées. Dans le cas décrit, ces capteurs et appareils sont comman- dés pour que, toutes les 250 ms, un nouveau jeu de para- mètres PO' Px' Py, PZ XI eyt Z' Sx' Sy, Sz,, e , et V soit disponible. Chaque jeu de paramètres disponible est donc mémorisé à titre provisoire, et ne le sera à titre définitif que s'il satisfait à un certain nombre de critères. C'est ainsi que les dérivées par rapport au temps des expressions PX_-X et Py- ey sont calculées en permanence, et que tout jeu de paramètre disponible à un instant ou l'une de ces dérivées est supérieure à un seuil qui lui est propre n'est pas retenu. En effet, il est connu que, en première approximation, les vitesses VX et Vy sont proportionnelles aux expressions PX- eX et Py- ey respec- tivement. On évite ainsi de retenir des configurations de vol relatives aux phases de fortes accélérations et de régime transitoire, qui conduiraient à de mauvais résultats. De plus, afin d'éliminer les configurations de vol instables, on ne retient que celles pour lesquelles un nombre déterminé de jeux de paramètres successifs stables a été retenu, nombre par exemple égal à quatre. A cet effet, la valeur moyenne de chacun des paramètres sur les quatre jeux successifs est calculée, ainsi que l'écart type de ce paramètre, qui est comparé à un seuil propre à ce paramètre. Si les écart types de tous les paramètres sont inférieurs à leurs seuils respectifs, alors une confign- ration de vol est retenue. Dans cette configuration, la valeur retenue pour chaque paramètre est la valeur moyenne calculée sur les quatre jeux successifs. Enfin, afin d'assurer une bonne répartition des configu- rations de vol retenues dans tout le domaine des configu- rations possibles de l'hélicoptère 1 en vol normal, il est prévu de diviser ce domaine en un certain nombre de sous-domaines, et d'assurer, dans chacun de ces sous- domaines, un nombre de configurations de vol retenues compris entre une valeur maximale et une valeur minimale. Ici, et à titre d'exemple, les sous-domaines sont au nombre de 12, et ils sont définis de la façon suivante. Tout d'abord, 8 sous-domaines, numérotés de 1 à 8, cor- * respondent à des vitesses air faibles, c'est-à-dire dont le module V est inférieur au seuil VS en deçà duquel la mesure de ce module par le capteur anémobarométrique est entachée d'erreur. Les sous-domaines 1 à 8 correspon- dent aux critères suivants, les vitesses étant exprimées en noeuds: sous-domaine 1, ou de montée: Sz -7,5 -2,5 - sous-domaine 4, ou d'hybridation: 3 +7,5 -2,5 SX -2,5 Ensuite, 4 sous-domaines, numérotés de 9 à 12, correspon- dent à des vitesses air élevées, c'est-à-dire dont le module V est supérieur au seuil VS au delà duquel la mesure de ce module par le capteur anémobarométrique est valable. Les sous-domaines 9 a 12 correspondent aux critères suivants: - sousdomaine 9, ou de palier: -2,5 sous-domaine elle est susceptible d'être classée, en fai- sant des essais sur les valeurs des paramètres V, SX, Sy et SZ, ainsi que sur la valeur de la dérivée par rapport au temps de l'angle de cap t, pour cette configuration de vol. Lorsque le sous-domaine dans lequel la configuration de vol à classer est susceptible d'être classée est détermi- né, on regarde si le nombre maximal de configuration classées dans ce sous-domaine est atteint. Si tel est le cas, le sous-domaine est plein, et la configuration à classer est rejetée. Si le sous-domaine n'est pas plein, l'écart entre la confi- guration à classer et chacune des configurations déjà classées est calculé. Si cet écart est supérieur à un seuil pour chaque configuration déjà classée, alors la configu- ration à classer est effectivement classée dans son sous- domaine. L'écart calculé ici est le carré des différences entre les deux valeurs de la composante SX et les deux valeurs de la composante Sy Le classement comprend essentiellement la mémorisation des paramètres P0O PX' Py' Pz' eX sx sY, s, V, Y, et q ainsi que de données de datation de la configuration de vol. Comme données de datation, on choisit ici le rang de classement de la configuration et l'intervalle de temps qui la sépare de la configuration de rang précédent, que cette configuration soit ou non classée dans le même sons-domaine. Ainsi on peut rétablir si nécessaire la chronologie du classement et connaître l'instant auquel les données relatives à toute configura- tion de vol classée ont été mesurées. Au cours du vol d'étalonnage, le pilote commande l'héli- coptère 1 pour permettre le classement d'un nombre suffi- sant de points de vol dans chaque sous-domaine. Il est prévenu dès qu'un point est classé dans un sous-domaine, et il est informé sur l'état de remplissage des divers sous-domaines. Dès que le nombre minimal de points est atteint dans chacun des sous-domaines 1 à 12, le vol d'étalonnage est terminé. La deuxième étape du procédé peut alors commencer. Celle- ci est effectuée en temps différé, après la première étape d'acquisition des valeurs de PoP PX? PY Pz ' x Y' ?Z' SX, Sy', S V, v, a et q pour les N configurations de vols régulièrement réparties dans le domaine des configu- rations possibles pour l'hélicoptère 1. La première tâche effectuée au cours de la deuxième étape concerne l'estimation des composantes VX et Vy de la vi- tesse air pour chacune des configurations de vol classées, à partir de la seule connaissance des valeurs de SX, Sy' sz, y,, e et pour cette configuration. Cette estima- tion repose ici sur l'hypothèse que le vent, exprimé dans le repère terrestre, est resté constant au cours de la durée du vol d'étalonnage. Les composantes WN, WE et WV de ce vent, exprimé dans le repère terrestre, sont alors --4 déterminées. La connaissance de la vitesse du vent W ainsi déterminée, et de la vitesse sol S mesurée pour chaque configuration de vol permet d'estimer la vitesse sol V pour chaque configuration de vol, compte tenu de la relation connue: s = V' V (1) Pour déterminer les composantes du vent moyen, au cours du vol d'étalonnage, on utilise les configurations de vol classées dans le sous- domaine 9, car, dans ce sous-domaine, la mesure du module V de la vitesse air est valable et la composante Vz est proche de zéro. On suppose que le nombre de configurations de vol classées dans ce sous-domaine est égal à n. On calcule d'abord, pour chaque configuration de vol du sous-domaine 9 les trois composantes SN' SE et SV du vec- teur Sdans le repère terrestre, à partir de ses composantes SX, Sy et SZ mesurées dans le repère hélicoptère. Une telle transposition se fait à partir de formules de change- ment de repère bien connues de l'homme du métier, qui font intervenir les angles d'attitude À, $ et f. On dispose donc, pour chaque configuration de vol, de rang i, de quatre grandeurs scalaires qui sont: - le module V. du vecteur V, et, -les trois composantes SNi, SEi et SVi du vecteur S dans le repère terrestre. Or, la relation (1) peut s'écrire: S -W = v (2) La relation (2), écrite en module et à l'instant ti, four- nit l'équation: (Si - WN)2 + (SEi -)2 (S - w) 2= V 2 (3) L'équation (3) traduit le fait que, dans le repère terres- tre, l'extrémité du vecteur tse trouve sur une sphère SP cen:_ee au por.: de cmor=or=ees _,, 5_, et S. e: de r=von V. On dispose donc de n équa:ionr te--es que la re'as:on 3 (3), chacune d'entre e' les étant e:a:Àve à un istan- ce rang i différent. Si l'on soustrait membre à memboe l'équa: on rea-_ve à instant t. de cei; e re!a-ve à 1'nsza== tj, on obtient I l'équation (S,- W2 (S)2 (s E - (s.-.- (4 + (Svj _ Wv)2 _ (S_ _ w.)2 Vj+ (s- v) - (s - =v2 _V2 J i qui devient 2(SNi - Sj) W + 2 (S. - S Wj)+ 2 (Sv. - SVj) V. = $N2 2 2 2 2 2 2 S NJ S+ S E _ 5 S + VJ i L'équation (5) traduit le fait que, dans le repère terr es- tre, l'extrémité du vecteur -se trouve sur le plan radical de la sphère SPi et de la sphère SPJ, c'est-à- dire le plan qu contient le cercle selon lequel ces den= sphères se coupent, lorsque celles-ci ont au moine un point commun, ce qui est, en principe, le cas ici. En écrAvant toutes les équations telles que (3) pour tous les couples i, J, i étant dAffirent de J, on obtient un système de n (n-l)/2 équations à trois inconnues YNI VE *t W7. On résout alors ce système par la méthode des moindres carrés pour déterminer les trois inconnues WN, YE et WV, qui sont les trois composantes du vent moyen à déterminer. Dans le procédé ci-dessus, on comprend que c'est le fait de résoudre le système de n (n-l)/2 équations par régres- sion linéaire qui introduit un moyennage sur les compo- santes WN, WE et W. En effet, si l'on se place dans une situation théorique idéale, o le vecteur W est rigoureusement constant et o les mesures de Vi, SNi, SEi et Svi sont exactes, il suffit de trois équations telles que l'équation (5) pour détermi- ner WN, WE et WN. Dans ce cas, les trois sphères définies par ces trois équations ont un point commun et un seul dont les trois coordonnées sont les trois composantes recherchées. Dans ce cas idéal, la multiplication des ins- tants de mesure permet de former d'autres groupes de trois sphères dont le point commun est toujours le mame. Toutefois, dans les cas pratiques o le vecteur 'varie légèrement entre deux instants de mesure et o les mesures de Vi, SNi' SEi et Svi sont entachées d'erreurs, on comprend que les n sphères dtéquation (5) n'ont plus, à proprement parler, un point commun et que ce point commun éclate, en fait, en un nuage de points, nuage de volume plus ou moins large à l'intérieur duquel passent toutes les sphères. Dans ce cas, le procédé ci-dessus apporte une solution au problème posé par la détermination des coordonnées moyen- nes de ce volume. De plus, le passage par les équations (5) des plans radicaux, et l'application de la méthode des moindres carrés au système d'équations (5) conduit à un résultat qui se révèle tout à fait satisfaisant, à l'égard de la précision du résultat comme à l'égard du temps de calcul. Dans le cas particulier o la composante du vent vertical est faible, on peut simplifier le procédé de l'inven- tion, afinde gagner du temps de calcul pour la résolu- tion du système d'équations. En effet, si le vent vertical est faible, on peut admettre que la composante VVi selon l'axe vertical du repère terrestre, de la vitesse propre V de l'hélicoptère vaut: Vvi = V. sin e. (6) Vi expression dans laquelle A est l'angle de tangage mesuré au même instant t. que le module V. de la vitesse i 1 propre V. La projection de la relation (2) sur l'axe vertical du repère terrestre fournit l'équation: SVi - WVi = Vi sin e () ce qui permet immédiatement la détermination de WVi et de la composante W comme moyenne de WVi, selon la V relation i wv = (l/i) Svk - Vk sine k (8) k=l La projection de la relation (2) sur le plan horizontal du repère terrestre fournit l'équation: (s w - WN (S, s - = r cc s ( 9r) L' équa;ion (9) taldu= e fa-= que, dans le plan hor:- zon:a., l'exreme le a projecó=on du vecteur W se trouve sur un cercle C. cenré au poi=: de coordcnées (SN., SEi) et de ra30on V cos >.. C= dlsose donc de n éaua::ors::es qoe la rea:tion (9). Si on soustrait membre à membre 1'équation relative à 1':nstant t. de celle re:ative à l'inszant tj, on obtient l'équation: 2 _>(S _ 2 (sNj - WN) - (s - %2 + (SEi - w2 (sEi -) = V cos2 e - V2 cos2 i (10) qui devient 2(si - sj) w + 2 (s - sE) () 2s -ss2 - s 2 2 = SNi - SNj + Si - 5j +V cos2 L'téquation (11) traduit le fait que, dans le plan hori- zontal, la projection de l'est,émité du vecteur W se trouve sur l'axe radical du cercle Ci et du cercle Cj c'est-à-dire de la droite qui contient les deux points d'intersection de ces deux cercles. On obtient alors, comme précédemment, un système de n(n-1)/2 équations, mais à deu= inconnues WN et yE seule- ment, système que l'on résout par la méthode des moindres carrs. Lorsque les composantes W, WE et Wz du vent moyen au cours du vol d'étalonnage ont été déterminées à partir des données relatives aux configurations de vol classées dans le sous-domaine 9, on calcule les composantes WXi et WYi du vent moyen dans le repère OXYZ de l'hélicoptère 1 pour chacune des N configurations classées dans l'ensem- ble des sous-domaines 1 à 12. A cet effet, on utilise les formules connues de changement de repère, avec les valeurs mesurées ti' a, iet ii des angles d'attitude pour chaque configuration. L'application de la relation (2) dans le repère OXYZ permet de déterminer les valeurs estimées Vxi et Vyi de la vitesse air par les formules: V.= S.- w. Vxi Sxi -xi _1 (12) V =S -WI VYi Yi WYi On aborde alors la deuxième tache de la deuxième étape du procédé de l'invention, qui concerne la formation de deux matrices de corrélation, une relative à la vitesse estimée VX, et l'autre relative à la vitesse estimée V. A cet effet, on calcule d'abord, pour chacune des N confi- gurations de vol, les valeurs d'un certain nombre de va- riables construites à partir des pas Po, PX, Py et PZ, des angles eX et ey, et de l'accélération IZ. Ici, ces variables construites comprennent: 5 variables simples du premier ordre, à savoir Po,' PX' Py' Pz et Yz, - 2 variables combinées du premier ordre, à savoir p X-X et P -r -, et, - 28 variables simples du second ordre, à savoir les carrés et les produits deux à deux des variables P0ot PX' PY Pz, X, 8ey et z. Les 28 variables simples du second ordre comprenant tous 2 2 2 les carrés PO2' PX' ' Z et tous les produits deux à deux P0PX PoY, o eyyz On forme alors un premier ensemble d'ici 36 variables comprenant les 35 variables construites précédentes et la vitesse estimée VX. Ces 36 variables prennent chacune M valeurs, chacune de ces M valeurs étant relative à une des M configurations de vol classées dans les sousdomaines 1 à 8 au cours du vol d'étalonnage. Il est donc possible de calculer, pour chacun des couples possibles de ces 36 variables, un coefficient de corrélation qui traduit le fait que les variations des valeurs de chaque variable du couple, sur les M configurations de vol, sont plus ou moins corrélées. La formule de calcul d'un tel coefficient à partir des M valeurs différentes de chaque variable est bien connue de l'homme de métier. Comme une façon habi- tuelle de présenter les résultats de ces calculs est de former un tableau d'ici 36 lignes sur 36 colonnes, d'af- fecter chacune des 36 variables à chacune des 36 lignes et à chacune des 36 colonnes, et de placer le coefficient de corrélation de chaque couple au croisement de la ligne et de la colonne correspondant à ce couple, on dit que l'on calcule la matrice des coefficients de corrélation, ou matrice de corrélation. On forme également un deuxième ensemble comprenant les variables construites précédemment et la vitesse esti- mée Vy,, pour lequel on calcule la matrice de corrélation. Toutefois, pour la vitesse estimée Vy, on tient compte, pour le calcul de la matrice de corrélation, des N confi- gurations de vol classées dans les sous-domaines 1 à 12, et non pas seulement de M configurations de vol classées dans les sous-domaines 1 à 9. En effet, comme cela a été indiqué, les sous-domaines 9 A 12 correspondent à un module V de la vitesse air supérieur au seuil V$ par exemple égal à 45 noeuds. Or, dans ce cas, non seulement la mesure de ce module est valable, mais de plus, la disposi- tion du tube de Pitot est telle que l'on peut admettre que Vx v (13) Il en résulte que, en pratique, lorsque le module V est supérieur au seuil VS, on ne détermine pas VX par le calcul, mais on utilise le résultat de la mesure anémo- barométrique. Il est donc inutile de prendre en compte les données relatives aux sous-domaines 9 à 12 pour éta- blir les formules de calcul de la vitesse VX, puisque pour des configurations comparables survenant en cours de vol normal, on n'utilisera pas la formule de calcul de Vx. Par contre, on tient compte de la totalité des sous-domaines pour la vitesse Vy, car on détermine toujours cette vites- se par le calcul. Les deux matrices de corrélation étant ainsi construites, on détermine, au cours de la troisième étape, les formules qui permettront de calculer, en vol normal les composantes VX et Vy en fonction des grandeurs Pot PX, PY' PZ AX' 8yet z. A cet effet, et en considérant d'abord le cas de la compo- sante Vx, on détermine d'abord la variable construite la plus corrélée à la vitesse estimée Vx. Pour cela on déter- mine, à l'aide de la matrice de corrélation relative à VX, la variable construite dont le coefficient de corréla- tion avec Vx est le plus élevé. En principe, et comme la théorie et l'expérience le confirment, cette première va- riable la plus corrélée s'avère 8tre la variable (P.- eX) o On soustrait alors de la vitesse estimée Vx la variable la plus corrélée, ici (PX- eX) multiplié par un coeffi- cient multiplicatif a1 choisi pour que la vitesse estimée modifiée soit totalement décorrélée de (Px- GX). Si on appele Vxm la vitesse estimée modifiée, on peut écrire: Vxm = Vx ' al (Px x) (14) La valeur du coefficient ai est liée à celle du coefficient de corrélation et aux valeurs moyennes de VX et de (Px- X) par une loi bien connue de l'homme de métier. On calcule également l'erreur commise en admettant que la valeur modifiée Vxm est nulle, c'est-à-dire que la vitesse estimée Vx est égale au produit de la variable la plus correlée (Px- eX) par le coefficient multiplicatif ai. On calcule alors une matrice de correlation modifiée, en remplaçant la vitesse estimée Vx par la vitesse estimée modifiée V... On détermine alors la variable la plus cor- rélée à la vitesse estimée modifiée Vxm, ou deuxième variable la plus corrélée à la vitesse Vx. Supposons ici que cette variable soit la variable PxPz. On soustrait alors PxPz, multipliée par un coefficient multiplicatif a2, pour que la variable estimée modifiée une deuxième fois, se trouve totalement décorrélée de la variable PxPz La variable modifiée une deuxième fois vaut alors: 3o Vxm - al(PX- eX) - a2PxPz (15) On calcule l'erreur commise en admettant que cette valeur modifiée est nulle, et ainsi de suite. On recommence ainsi jusqu'à avoir soustrait de la vitesse estimée Vx un nombre nx de variables, ici au plus égal à , suffisant pour que l'erreur commise soit inférieure à un seuil donné. En admettant que l'erreur est nulle, on dispose alors d'une formule qui permet de calculer la vitesse estimée VXen fonction des nx variables les plus corrélées, c'est- à-dire en fonction de POt PX, PY' PZ eX, ey et yz. Pour un type d'hélicoptère, cette formule de calcul pourra être de la forme: II 2 VX al(PX eX) + a2PxPz +...+anx 2 (16) tandis que pour un autre type, elle pourra ftre de la forme: vX at î(Px-ex) +..+a' Py z (16') et ainsi de suite. On notera que l'algorithme défini ci-dessus est celui connu sous le nom de régression pas à pas, ou, selon la terminologie anglo-saxonne plus répandue, "stepwise regression". En appliquant le même algorithme à la matrice de correla- tion relative à la vitesse estimée V^y, on aboutit à une formule semblable, par exemple de la forme: Vy= b1( y- Y + 2PY X + +bnyPZ () ou encore de la forme: V b(- b'29X +b'ypz (1') y=bl(P -e)+ b'1 +..+b, (l?')Z On utilise alors les formules (16) et (17), ou, le cas échéant, les formules (16') et (171) pour calculer, au cours de la quatrième et dernière étape du procédé de l'invention, les valeurs des vitesses VX et Vy en fonc- tion des valeurs de POP PX, PY' PZ, eX, ey et YZ mesu- rées en permanence au cours du vol normal et pour afficher les valeurs calculées VX et Vy devant le pilote de l'hé- licoptère 1. Cependant si le procédé tel qu'il vient d'être décrit est très efficace pour déterminer la forme, ou la structure des formules de calcul, il se peut dans certains cas qu'il conduise à des valeurs numériques des coefficients multi- plicatifs tels que ai et bi qui manquent de précision. i i Par structure des formules de calcul, il faut entendre ici le nombre et la nature des variables construites qu'elles font intervenir, indépendamment des valeurs numériques des coefficients multiplicatifs affectant ces variables dans les formules. Par exemple, les formules (16) et (16') sont de structure différentes dans la mesure o la variable PY eX intervient dans la formule (16) et non dans la formule (16'). Pour obtenir des valeurs plus précises pour les coefficients multiplicatifs ai et bi, le procédé de l'invention prévoit donc, à titre de perfectionnement, d'utiliser les formules de calcul dont la structure est déterminée comme cela a été décrit, mais dQnt les valeurs des coefficients multi- plicatifs ai et b. peuvent être calculés indépendamment 1 1 de façon plus précise, afin de corriger leurs valeurs obtenues précédemment. A cet effet, il est possible, à la fin de la troisième étape, lorsque la structure de chacune des deux formules de calcul est connue, de considérer les coefficients multiplicatifs ai et b. comme des inconnues. En repor- tant, dans la formule donnant VX, et pour chacune des M configurations de vol appartenant aux sous-domaines 1 à 8, les résultats des mesures relatives à P0o PX, PY' PZ' PX' e et yZ effectuées au cours de la première étape, et en identifiant la quantité ainsi obtenue à la vitesse VX estimée au cours de la deuxième étape, on obtient un sys- tème de M équations linéaires dont les n inconnues sont x les coefficients multiplicatifs a.. Une résolution de ce système par la méthode des moindres carrés permet de déterminer avec une bonne précision les valeurs des coef- ficients a.. On corrige alors les valeurs de ces coeffi- cients dans la formule destinée au calcul de V en vol x normal, en remplaçant leurs anciennes valeurs par celles ainsi déterminées. On procède de même pour la détermination des coefficients bi intervenant dans la formule de calcul de Vy, en consi- dérant cette fois l'ensemble des N configurations de vol des sous-domaines 1 à 12. I1 est possible d'améliorer, selon une première variante du procédé de l'invention, le calcul des coefficients multiplicatifs ai et bi, sans avoir à faire l'hypothèse que le vent reste constant pendant toute la durée du vol d'étalonnage, hypothèse qui, on le rappelle, est nécessai- re pour l'estimation des vitesses VX et Vy dans chaque configuration de vol du vol d'étalonnage. Selon cette première variante, on choisit, parmi-les M configurations de vol des domaines 1 à 8, des paires de configurations de vol d'instants de datation suffisamment proches pour que l'on puisse admettre que le vent est resté le même entre ces deux configurations. A cet effet, après avoir fixé un intervalle de temps maximal, entre deux configurations, intervalle de temps en deçà duquel on peut admettre que le vent est resté le même, on part, par exemple, de la dernière confign- ration de vol classée, et on remonte systématiquement la suite des configurations de vols classées, en formant des paires appartenant toutes deux aux sous-domaines 1 à 8, séparées par un intervalle de temps inférieur à l'in- tervalle de temps maximal admis. Afin d'éviter de prendre plusieurs fois la même configuration de vol dans des paires différentes, on marque toute configuration prise dans une paire. Pour chacune des deux configurations de chaque paire, par exemple la configuration de rang k et celle de rang $, on transpose dans le repère terrestre, les résultats k et des mesures de la vitesse de l'hélicoptère par rapport au sol, grace aux formules connues de changement de repère. De même, on transpose dans le repère terrestre les ré- _ - sultats Vk et VL donnés par les formules de calcul de la vitesse air. Comme le vent est supposé constant entre les deux confi- gurations, on peut écrire, dans le repère terrestre: Sk: k (18) La formule (18) traduit le fait que la variation mesurée de la vitesse par rapport au sol de l'hélicoptère entre les deux configurations de rang i et j est égale à la va- riation calculée, à l'aide de formules de calcul, de la vitesse air de l'hélicoptère 1, puisque l'on suppose que le vent n'a pas varié entre ces deux configurations. Le report, dans l'équation (18), des resultats mesurés dans chaque configuration de vol, les-coefficients ai et b. étant considérés comme des inconnues, permet d'établir, pour chaque paire de configurations, une équation linéaire. La résolution du système d'équations linéaires ainsi obtenu par la méthode des moindres carrés permet la détermination des coefficients a. et b. et permet donc I 1. l'identification simultanedes formules de calcul de VX et de Vy sur les sous-domaines 1 à 8. Pour prendre en compte les points des sous-domaines 9 à 12 pour l'identification de la formule de calcul de Vy un processus analogue peut être repris sur l'ensemble des domaines 1 à 12, en prenant pour vitesse VX celle donnée par la formule de calcul pour les points des sous- domaines 1 à 8, et celle mesurée par le capteur anémo- barométrique pour les points des sous-domaines 9 à 12. Selon une deuxième variante, on peut, au lieu de former des paires de points, admettre que le vent suit, dans le repère terrestre, une loi polynomiale en fonction du temps, de coefficients inconnus, soit, par exemple, en supposant que la composante verticale du vent est nulle: m WN L ai ti i=o (19) m WE =-1 Ri ti i=o En transposant la vitesse W -dans le repère hélicoptère OXYZ, on obtient alors, pour chaque configuration de vol classée, une équation linéaire en identifiant le résultat de la mesure de la vitesse S par rapport au sol à la somme de la vitesse du vent, déduite des relations (19), et de la vitesse par rapport à l'air, obtenue par les formules de calcul. On applique ainsi simplement la relation (1): s = X + v (1) On peut ainsi obtenir un système d'équations linéaires, dont les inconnues sont les coefficients a., bi, ainsi que les coefficients inconnus ci et 2i de la loi poly- nomiale du vent. Naturellement, des simplifications sont possibles si le vol d'étalonnage se produit dans un plan parallèle au sol, car les formules de passage du repère terrestre au repère hélicoptère et inversement ne dépendent plus que de l'angle de cap Y. De même, lorsque les structures des formules de calcul ont été déterminées et que l'on désire simplement effec- tuer une nouvelle détermination des coefficients multi- plicatifs a. et bi., il est possible d'effectuer un vol 1 1 d'étalonnage simplifié. Comme cela a été indiqué, on prévoit, en réponse au ré- sultat de la mesure du module V de la vitesse air effec- tué par le capteur anémobarométrique, de sélectionner, comme valeur de la vitesse VX par rapport à l'air, le résultat de la mesure si celui-ci est supérieur au seuil VS, ou celui du calcul par la formule de calcul, si le résultat de la mesure est inférieur au seuil VS, ou encore une combinaison des deux, dans les cas intermédiaires. Ceux-ci correspondent ici au sous-domaine 4, appelé pour cette raison sous-domaine d'hybridation. Le procédé de l'invention est mis en oeuvre gràce au 3? dispositif qui est maintenant décrit en référence à la figure 2. L'élément principal du dispositif est un microprocesseur 15, ici du type commercialisé par la société INTEL sous la référence 8086, pourvu d'un circuit de mémoire 2. Le circuit de mémoire 2 comprend une mémoire de type ROM de 64 k octets, dans laquelle est stocké le programme de commande du microprocesseur, une mémoire de type RAMl de 16 k octets et une mémoire de type EEPROM de 16 k octets, ces deux dernières mémoires étant destinées à mémoriser les données manipulées par le microprocesseur 15. Le microprocesseur 15 est raccordé par un bus 11 aux instruments de bords de l'hélicoptère 1. * Ces instruments comprennent quatre capteurs de déplace- ment 3, 4, 5 et 6, montés sur les commandes de pas collec- tif, de pas longitudinal, de pas transversal et de pas du rotor anticouple, respectivement, qui délivrent les grandeurs Pop PX, Py, et PZ sous forme numérique, respec- tivement. Les instruments de bord comprennent également un dispo- sitif accélérométrique 7, de type connu, qui mesure les trois accélérations yX, y et Z de.l'hélicoptère 1 le long des axes OX, OY et OZ. Le dispositif accélérométrique 7 est pourvu ici d'une unité de calcul qui calcule les angles d'inclinaison X et ey d'après les formules connues: X = x YZ (20) et O =. /. (21) Si ces angles n'étaient pas calculés dans le dispositif accélérométrique 7, ils pourraient être mesurés par des inclinomètres de type connu. Le bottier accélérométrique ? délivre les grandeurs DX, y et yZ sous forme numérique. Les instruments de bord comprennent également un disposi- tif 8 de mesure des trois composantes SX, Sy et Sz de la vitesse par rapport au sol, dans le repère OXYZ. Un tel dispositif, qui comprend par exemple un radar Doppler, délivre les grandeurs SX, Sy et Sz sous forme numérique. Il est également prévu une centrale d'attitude 9 qui mesure les angles de cap Y, de tangage & et de roulis f de l'hélicoptère 1, et les délivre sous forme numérique. Enfin, un capteur anémobarométrique 10 mesure le module V de la vitesse air et délivre cette grandeur sous forme numérique. Par ailleurs, il est prévu un indicateur de vol station- naire 12, du type connu à aiguilles croisées, qui permet au pilote, de prendre connaissance de la vitesse air de composantes VX et Vy calculées par le microprocesseur 15. Un tel indicateur comprend deux aiguilles perpendiculaires mobiles en translation perpendiculairement à leurs direc- tions longitudinales respectives. La première est entrai- née en un déplacement proportionnel à la composante V et la deuxième en un déplacement proportionnel à la composante Vy. Lorsque les composantes VX et Vy sont nulles, les aiguilles se croisent au centre du cadran. Lorsque les composantes VX et Vy sont non nulles elles se croisent en un point qui s'écarte du centre du cadran de façon repré- sentative du vecteur de composantes VX et Vy. Le programme qui commande le microprocesseur 15 est mémorisé dans la mémoire 2. Il est prévu pour que le microprocesseur 15 exécute toutes les tâches qui viennent d'être exposées, et délivre à l'indicateur 11 les données VX et Vy en vol normal, calculées à l'aide des formules de calcul déterminées à la suite du vol d'étalonnage. La programmation du microprocesseur 15 pour qu'il effec- tue les t:ches décrites est évidemment à la portée de l'homme du métier, et ne sera donc pas davantage décrite. A titre d'exemple, pour un nombre N de points de l'ordre de 300, le temps de calcul nécessaire pour la détermina- tion de la structure des formules de calcul, et pour la détermination des valeurs de leurs coefficients multipli- catifs est de l'ordre de 20 à 40 minutes. Lorsque, comme cela a été décrit, le microprocesseur 15 commande les instruments de bord pour disposer d'un jeu de paramètres toutes les 250 millisecondes, la durée du vol d'étalonnage est de l'ordre de une heure, toujours dans l'hypothèse ou l'on prévoit un nombre N de configu- rations classées de l'ordre de 300. Revendications 1. Procédé de détermination de la vitesse par rapport à l'air d'un hélicoptère en vue de son pilotage, hélicop- tère (1) d'axe longitudinal OX, d'axe transversal OY et d'axe vertical OZ, pourvu d'un rotor de sustentation dont le pas collectif P0, le pas cyclique longitudinal PXet le pas cyclique transversal Py sont commandables, ainsi que d'un rotor anticouple dont le pas PZ est com- mandable, procédé dans lequel on détermine les deux composantes V.. et V de cette vitesse sur les axes OX et 1% Y OY, respectivement, en fonction des valeurs commandées desdits pas Po0 PXt Py et PZ, des valeurs des angles d'inclinaison eX et 8y de l'hélicoptère (1) autour des axes OY et OX, respectivement, et de la valeur d'accélé- ration yZ de l'hélicoptère selon l'axe OZ, procédé caractérisé par le fait que: - au cours d'une première étape, on effectue un vol d'étalonnage faisant passer l'hélicoptère par une plura- lité de N configurations de vol différentes, pour chacune desquelles on mesure, d'une part, une pluralité de gran- deurs (Sx, Sy, SZ, e,.,V) permettant d'estimer les vitesses VX et Vy, et d'autre part les valeurs des pas P0, PX, Py et PZ, des angles &X et ey et de l'accéléra- tion 7z, - au cours d'une deuxième étape, on estime lesdites vi- tesses VXet Vy lors de chaque configuration de vol, on calcule une première matrice des coefficients de corréla- tion de tous les couples possibles de variables d'un premier ensemble de variables comprenant la vitesse estimée VX et des variables construites à partir des pas P0ot PX, Pyet PZ' des angles eX et ey et de l'accélération yZ, et on calcule une deuxième matrice de façon comparable à la premiere mais relative-à un deuxième ensemble de variables comprenant la vitesse estimée Vy et lesdites variables construites, - au cours d'une troisième étape, on détermine, à l'aide de ladite première matrice, la variable construite du premier ensemble la plus corrélée à la vitesse estimée Vx, on soustrait de cette vitesse estimée V ladite va- X riable construite le plus corrélée mut:ipliée par un coefficient multiplicatif choisi pour que la vitesse estimée ainsi modifiée soit totalement décorrélée de ladite variable construite la plus corrélée, on calcule l'erreur commise en admettant que la valeur modifiée est nulle, on calcule une première matrice modifiée comme on a calculé ladite première matrice mais en remplaçant la vitesse estimée VX par la vitesse estimée modifiée, et on recommence ainsi jusqu'à avoir soustrait de la vitesse estimée VX un nombre nx de variables construites du pre- mier ensemble, suffisant pour que l'erreur commise, en admettant que la valeur modifiée est nulle, soit inférieure à un seuil, de façon à disposer, pour la valeur de la vitesse estimée VX, d'une formule de calcul en fonction de n variables construites du premier ensemble, et on procède de même pour la deuxième matrice, de façon à disposer d'une formule semblable pour V , en fonction de n variables construites du deuxième ensemble, et, y - au cours d'une quatrième étape, en vol normal, on mesure les valeurs des pas Po, PX, Py et PZ, des angles eX et &y et de l'accélération 7Z, on calcule les vitesses VX et Vy en appliquant les formules déterminées au cours de la troisième étape, et en réponse aux valeurs mesurées au cours de la quatrième étape, et on affiche, à l'usage du pilote de l'hélicoptère (1), les valeurs VX et Vy cal- culées. 2. Procédé selon la revendication 1, dans lequel les variables construites au cours de la deuxième étape à partir des pas Pop PX Py et Pz' des angles X et e et de l'accélération]Z comprennent: - des variables simples du premier ordre, à savoir Pot PXP PY' Pz et RZi - des variables combinées du premier ordre, à savoir Px'- et pY-Oy, et, - des variables simples du second ordre, à savoir les carrés et les produits deux à deux des variables PO, PX' PYI PZ'.XI t yet yZ. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel, à la fin de la troisième étape, on considère comme des inconnues les n et n coefficients multiplica- x y tifs des formules de calcul de V et de Vy, respectivement, x on reporte dans ces formules les résultats des mesures effectuées au cours de la première étape, afin de former au moins un système d'équations linéaires dont les nx et n coefficients multiplicatifs, sont les inconnues, on Y résout ce système et on corrige, dans les formules de cal- cul de VX et de V, les valeurs des n et ny coefficients multiplicatifs, respectivement, en les remplaçant par celles ainsi déterminées. 4. Procédé selon la revendication 3, dans lequel on résout le système d'équations linéaires par la méthode des moindres carrés. 5. Procédé selon l'une des revendications 3 ou 4, dans lequel, on forme un premier et un deuxième systèmes d'équations linéaires, en identifiant, pour une pluralité de configurations de vol prises parmi les N configurations de vol, la vitesse estimée VX ou Vy obtenue au cours de la deuxième étape, à la quantité obtenue par la formule de calcul correspondante. 6. Procédé selon l'une des revendications 3 ou 4, dans lequel on mesure, au cours de la première étape, la vi- tesse par rapport au sol et les angles d'attitude de l'hélicoptère pour chacune des configurations de vol, ainsi que l'instant de datation de celle-ci, et, pour former le système d'équations linéaires, on choisit parmi les configurations de vol, des paires de configurations de vol d'instants de datation suffisamment proches pour admettre que le vent est le même dans ces deux configu- rations de vol, on transpose dans le repère terrestre, et pour chaque configuration de vol, le résultat de la mesure de la vitesse par rapport au sol et les formules de calcul des composantes de la vitesse par rapport à l'air, et on prend pour équations à résoudre celles obte- nues en identifiant la variation mesurée de la vitesse par rapport au sol de l'hélicoptère entre les deux confi- gurations d'une paire, et la variation calculée à l'aide des formules de calcul, de la vitesse par rapport à l'air de l'hélicoptère. 7. Procédé selon l'une des revendications 3 ou 4, dans lequel on mesure, au cours de la première étape, la vites- se par rapport au sol et les angles d'attitude de l'héli- coptère pour chacune des configurations de vol, ainsi que l'instant de datation de celle-ci, et, pour former le système d'équations linéaires, on admet que le vent suit, dans le repère terrestre, une loi polynomiale en fonction du temps, de coefficients inconnus, on transpose la vites- se du vent dans le repère hélicoptère pour chaque confi- guration de vol, et on identifie, pour une pluralité de configurations de vols prises parmi les N configurations de vol, le résultat de la mesure de la vitesse par rapport au sol à la somme de la vitesse du vent et de la vitesse par rapport à l'air obtenue par les formules de calcul de V_ et Vy, de façon à déterminer en même temps les n et n coefficients multiplicatifs et les coefficients inconnus de la loi polynomiale du vent. 8. Procédé selon l'une des revendications 1 à 7, dans lequel, au cours de la première étape, on détermine les N configurations de vol différentes afin qu'un nombre minimal d'entre elles appartienne à chacun des sous- domaines du domaine des configurations possibles de l'héli- coptère en vol normal, chaque sous-domaine étant défini par un ensemble de plages de valeurs des grandeurs mesu- rées pour chaque configuration de vol. 9. Procédé selon la revendication 8, dans lequel, au cours de la première étape, on mesure, à l'aide de moyens anémo- barométriques, le module de la vitesse de l'hélicoptère par rapport à l'air, et chaque sous-domaine est également défini par des plages pour la vitesse anémobarométrique ainsi mesurée. 10. Procédé selon l'une des revendications 8 on 9, dans lequel on détermine les N configurations de vol pour que l'écart entre deux quelconques d'entre elles appartenant à un même sous-domaine soit supérieur à un seuil. 11. Procédé selon l'une des revendications 1 à 10, dans lequel, on mesure, au cours de la première étape, la vi- tesse par rapport au sol et les angles d'attitude de l'hé- licoptère, ainsi que le module de la vitesse par rapport à l'air de l'hélicoptère, on admet que le vent moyen dans le repère terrestre est resté constant pendant le vol d'étalonnage, on forme une pluralité d'équations, chaque équation identifiant, pour une configuration de vol et dans le repère terrestre, le module mesuré de la vitesse par rapport à l'air et le module calculé de la différence entre la vitesse par rapport au sol et la vitesse du vent, on soustrait membre à membre, et deux à deux, ces équa- tions, on résout le système d'équations ainsi formé, pour calculer la vitesse du vent moyen pendant le voi d'étalon- nage, et, au cours de la deuxième étape, on transpose la vitesse du vent moyen ainsi calculée dans le repère héli- copt-re pour estimer la vitesse air comme somme de la vitesse par rapport au sol et de la vitesse du vent moyen. 12. Procédé selon l'une des revendications 1 à 11, dans lequel on mesure, au cours de la quatrième étape en vol normal, à l'aide de moyens anémobarométriques, la vitesse de l'hélicoptère par rapport à l'air, et en réponse à la valeur de la vitesse anémobarométrique ainsi mesurée, on sélectionne, comme valeur de la vitesse par rapport à l'air, le résultat de la mesure, celui du calcul par les formules de calcul déterminées au cours de la troisième étape, ou une combinaison des deux. 13. Dispositif pour la mise en oeuvre du procédé selon l'une des revendications 1 à 12, pour la détermination de la vitesse par rapport à l'air d'un hélicoptère en vue de son pilotage, hélicoptère (1) d'axe longitudinal OX, d'axe transversal OY et d'axe vertical OZ, pourvu d'un rotor de sustentation dont le pas collectif P0, le pas cyclique longitudinal PX et le pas cyclique transversal Py sont commandables, ainsi que d'un rotor anticouple dont le pas Pz est commandable, dispositif de détermination des deux composantes V et Vy de cette vitesse sur les axes OX et OY, respectivement, en fonction des valeurs commandées desdits pas Pol PX' Py et PZ, des valeurs des angles d'in- clinaison e X et y de l'hélicoptère autour des axes OY et OX respectivement, et de la valeur de l'accélération Z de l'hélicoptère selon l'axe OZ, dispositif caractérisé par le fait qu'il comprend: - des moyens (3-10) pour mesurer, au cours d'une première étape au cours de laquelle est effectué un vol d'étalon- nage faisant passer l'hélicoptère par une pluralité de N configurations de vol différentes, d'une part une plura- lité de grandeurs (SX s, Sy$ 7 ' V) permettant d'estimer les vitesses V, et Vy, et d'autre part les valeurs des pas P., P", P et Pz, des angles ex. , et de l'accélération 7Z, - des moyens (15, 2) pour estimer, au cours d'une deuxième étape, lesdites vitesses VX et Vy lors de chaque configu- ration de vol, pour calculer une première matrice des coefficients de corrélation de tous les couples possibles de variables d'un premier ensemble de variables comprenant la vitesse estimée Vx et des variables construites à par- tir des pas P0, PX, Py et PZ, des angles PX et gy et de l'accélération 7Z, et pour calculer une deuxième matrice de façon comparable à la première mais relative à un deuxième ensemble de variables comprenant la vitesse esti- mee Vy et lesdites variables construites, - des moyens (15, 2) pour déterminer, au cours d'une troi- sième étape, et à l'aide de ladite première matrice, la variable construite du premier ensemble la plus corrélée à la vitesse estimée VX, pour soustraire de cette vitesse estimée V ladite variable construite la plus corrélée multipliée par un coefficient multiplicatif choisi pour que la vitesse estimée ainsi modifiée soit totalement décorrélée de ladite variable construite la plus corrélée, pour calculer l'erreur commise en admettant que la valeur modifiée est nulle, pour calculer une première matrice modifiée comme a été calculée ladite première matrice mais en remplaçant la vitesse estimée VX par la vitesse estimée modifiée, et pour recommencer ainsi jusqu'à ce que soient soustraites de la vitesse estimée V., un nombre n de variables construites du premier ensemble suffisant X pour que l'erreur commise, en admettant que la valeur modifiée est nulle, soit inférieure à un seuil, de façon à disposer pour la valeur de la vitesse estimée VX d'une formule de calcul, en fonction de n variables construites x du premier ensemble, et pour procéder de même pour la deu- xième matrice de facon à disposer d'une formule semblable pour V , en fonction de n variables construites du deu- y xième ensemble, et, - des moyens (3, 7) pour mesurer, au cours d'une qua:rième 13 étape, en vol normal, les valeurs des pas Pot PX' Py et Pz, des angles eX et ey et de l'accélération IZ, des moyens (15,2) pour calculer les vitesses VX et Vy à l'aide des formules déterminées au cours de la troisième étape et en réponse aux valeurs mesurées au cours de la quatrième étape, et des moyens (12) pour afficher, à l'usage du pilote de l'hélicoptère (1), les valeurs VX et Vy calculées.