FR 2947052 A1 20101224 FR 0903042 A 20090623 PROCEDE DE MESURE DE VITESSE D'UN AERONEF, PAR ANEMOMETRIE LASER DOPPLER Le domaine de l'invention est celui des procédés de mesure de vitesse d'un aéronef par anémométrie laser Doppler au moyen d'au moins un Lidar (de l'anglo-saxon "Light Detection and Ranging") dans lequel l'émission est réalisée par exemple par un laser (de l'anglo-saxon, "Light Amplification by Stimulated Emission of Radiation") à émission continue. Ce type de procédé consiste à mesurer le décalage de fréquence entre une onde de référence émise et une onde rétrodiffusée par les aérosols naturels de l'air, utilisés comme traceurs de vitesse, situés sur l'axe de visée. Ce décalage est représentatif de la projection, sur l'axe de visée, de la vitesse de l'aéronef relativement à l'air. Pour effectuer des mesures de référence, on focalise généralement sur une zone de focalisation choisie suffisamment loin de l'aéronef pour que la masse d'air qui s'y trouve soit aussi peu perturbée que possible par l'aéronef. Dans ce cas, l'onde rétrodiffusée est issue principalement des traceurs se trouvant dans la zone de focalisation. L'information utile, portée par le décalage de fréquence Doppler, est obtenue en effectuant une détection de type cohérent. Une onde issue d'une source de lumière cohérente est séparée en deux. Une première onde appelée � faisceau émis � est envoyée dans la zone de focalisation et une seconde onde appelée � oscillateur local � constitue une référence pour la détection cohérente du décalage de fréquence Doppler. Les aérosols présents naturellement dans l'atmosphère rétrodiffusent la lumière du faisceau émis. Le faisceau émis est rétrodiffusé à une fréquence décalée d'un décalage Doppler Afi par rapport à la fréquence du faisceau émis. Le signal rétrodiffusé est mélangé avec l'onde de référence pour former un battement hétérodyne sur la surface photosensible d'un détecteur. L'écart Doppler Afi, sur un axe de visée Xi, est l'écart entre la fréquence fi du signal rétrodiffusé et la fréquence F de l'oscillateur local. Il permet de déduire une mesure de la projection, sur l'axe de visée Xi, de la vitesse relative de l'aéronef par rapport au milieu ambiant, c'est-à-dire par rapport à la masse d'air se trouvant dans la zone de focalisation. L'expression qui lie ces deux grandeurs est la suivante: 2 vi=Afi.X/2 (A) vi étant la projection sur la ligne de visée Xi, du vecteur vitesse de l'aéronef par rapport à la masse d'air se trouvant sur l'axe de visée ou dans la zone de focalisation. étant la longueur d'onde du faisceau émis dans le milieu. L'écart Doppler Afi sur une ligne de visée Xi est obtenu classiquement par l'analyse du signal rétrodiffusé. Parmi les méthodes utilisables, l'analyse spectrale par transformée de Fourier est généralement utilisée. Le vecteur vitesse d'un aéronef est calculé à partir des composantes vi du vecteur vitesse suivant trois axes de visée non coplanaires, les composantes du vecteur vitesse vi étant obtenues par l'équation (A) à partir des écarts Doppler Mi des signaux rétrodiffusés sur les axes de visée concernés. Ce vecteur vitesse est exprimé sous la forme d'une Vitesse Vraie généralement nommée TAS (de l'anglo- saxon "True Air Speed") et des angle d'incidence et de dérapage généralement nommées AOA ( de l'anglo-saxon, "Angle Of Attack") et SSA (de l'ango-saxon, "Side Slip Angle") La fréquence du signal rétrodiffusé est déterminée à partir de la fréquence de la raie présente sur le spectre du signal rétrodiffusé. Cependant, lorsqu'un aéronef se déplace dans une masse d'air inhomogène en termes de vitesse, c'est-à-dire lorsque qu'il existe plusieurs masses d'air se déplaçant à des vitesses différentes, il n'est pas possible de déterminer de manière fiable les composantes du vecteur vitesse sur les axes de visée car il est parfois impossible de déterminer une fréquence unique sur les spectres des signaux rétrodiffusés sur les axes de visée. C'est le cas lorsqu'on est en présence de nuages lointains, de turbulences, d'échos de sol ou de perturbations de type vortex de sillage. Les procédés de détermination de la vitesse de l'aéronef connus présentent l'inconvénient, soit de déterminer une vitesse dont la valeur peut être fausse (le pilote ignorant que la valeur qui lui est fournie n'est pas fiable), soit de ne pas déterminer de vecteur vitesse (la vitesse de l'aéronef est alors indisponible). Les inconvénients précités sont spécifiques aux lasers qui émettent continûment un spectre étroit. Un procédé basé sur un laser pulsé ne présente pas ces inconvénients mais ce procédé est plus coûteux et moins précis, pour les mesures de vitesse, qu'un procédé basé sur un laser à émission continue car les spectres issus des signaux rétrodiffusés sont spectralement plus étalés. Le but de l'invention est de proposer un procédé de détermination de la vitesse d'un aéronef par rapport à une masse d'air qui soit précis, peu coûteux, fiable, ou avec une fiabilité connue et qui présente une bonne disponibilité, même lorsque la masse d'air dans laquelle se trouve l'avion est inhomogène. A cet effet l'invention a pour objet un procédé de mesure de la vitesse d'un aéronef à un instant donné, par anémométrie laser Doppler, ledit procédé comprenant les étapes suivantes : • acquisition d'un signal rétrodiffusé sur trois premiers 20 axes de visée orientés, • calcul du spectre du signal rétrodiffusé sur chacun des trois premiers axes de visée, caractérisé en ce que le procédé comprend en outre les étapes suivantes : 25 • acquisition d'un signal rétrodiffusé supplémentaire sur au moins un axe de visé supplémentaire non confondu avec un des trois premiers axes de visée, • calcul du spectre du signal rétrodiffusé sur au moins un axe de visée supplémentaire, 30 • calcul de plusieurs vecteurs vitesse candidats, • calcul des projections des vecteurs candidats sur au moins un axe de visée, appelé axe de calcul, • calcul de la cohérence des vecteurs candidats en fonction de leur(s) projection(s) sur le(s) axe(s) de calcul et en fonction du (ou des) spectre(s) du signal (ou des signaux) rétrodiffusé(s) sur le(s) axe(s) de calcul , • choix du vecteur vitesse mesuré parmi les vecteurs vitesse candidats en faisant une analyse de vraisemblance à partir des cohérences desdits vecteurs vitesse candidats. Le procédé a également pour objet une des caractéristiques suivantes prises seules ou en combinaison : - l'étape de calcul de la cohérence d'un vecteur candidat 10 comprend les étapes suivantes: • projection du vecteur candidat sur le(s)dit(s) axe(s) de calcul, • calcul de la cohérence individuelle entre le vecteur candidat et le signal rétrodiffusé sur chaque axe de calcul à partir de 15 la projection du vecteur vitesse candidat sur le(s)dit(s) axe(s) de calcul et le spectre du signal rétrodiffusé sur le(s)dit(s) axe(s) de calcul, • calcul du produit des cohérences individuelles, s'il y a plusieurs axes de calcul; 20 -l'analyse de vraisemblance consiste à choisir, parmi plusieurs vecteurs candidats : • soit le vecteur vitesse candidat qui a la cohérence la plus élevée, lorsque l'écart entre la cohérence la plus élevée et les cohérences des autres vecteurs vitesse candidats est supérieure à un 25 seuil d'écart de cohérence prédéterminé, • soit une combinaison de plusieurs vecteurs vitesse candidats comprenant le vecteur vitesse candidat dont la cohérence est la plus élevée ainsi que au moins un vecteur dont la cohérence présente, avec la cohérence la plus élevée, un écart inférieur à un 30 seuil d'écart de cohérence prédéterminé - on attribue un indicateur de confiance au vecteur vitesse mesuré; - préalablement à l'analyse de vraisemblance on rejette le(s) vecteur(s) candidat(s) qui a (ont) une cohérence inférieure à un seuil de 35 cohérence prédéterminé; - préalablement à l'analyse de vraisemblance on pondère la cohérence des vecteurs candidats en fonction d'au moins un vecteur vitesse mesuré à un instant précédent; -préalablement à l'analyse de vraisemblance on pondère la cohérence des vecteurs candidats en fonction d'une estimation du vecteur vitesse réalisée à partir d'au moins un vecteur vitesse mesuré à un instant précédant l'instant donné et éventuellement en fonction d'un intervalle de confiance attribué à l'estimation du vecteur vitesse à l'instant donné; - préalablement à l'analyse de vraisemblance on pondère la cohérence des vecteurs candidats en fonction d'informations inertielles; -comprend une étape d'identification de raies, à partir d'un critère de sélection de raies prédéterminé, sur les spectres des signaux rétrodiffusés sur les axes de visée; - on détermine les vecteurs candidats en: • sélectionnant 3 axes de base parmi les axes de visée, l' (ou les) axe(s) non sélectionné(s) étant appelé(s) axe(s) de redondance, • construisant, à partir des raies identifiées sur les axes de base, de vecteurs vitesse candidats dont les composantes sur les trois 20 axes de base sont égales aux composantes des raies identifiées sur les trois axes de base, et en ce que I' (ou les) axe(s) de calcul est (ou sont) le(s) axe(s) de redondance. I' (ou les) axe(s) de calcul est (ou sont) le(s) axe(s) de visée et en ce 25 en ce qu'on construit les vecteurs candidats en: • construisant des ensembles de raies présentant une dimension (M) au moins égale à 4 et égale au nombre d'axes de visée sur lesquels on a identifié au moins une raie, chaque raie d'un ensemble est prise sur un axe de visée différent des autres raies dudit 30 ensemble, • construisant, à partir de chaque ensemble de raies, un vecteur candidat minimisant les écarts avec ledit ensemble; - on construit un vecteur candidat, à partir de chaque ensemble, en mettant en oeuvre une méthode de résolution aux moindres 35 carrés; - il comprend une étape d'attribution à chaque projection du vecteur candidat sur un axe de calcul, d'un signal attendu centré sur la projection du vecteur candidat sur ledit axe et en ce que la cohérence individuelle d'un vecteur candidat sur un axe de calcul est calculée à partir du spectre du signal rétrodiffusé sur l'axe de calcul et du signal attendu; - la cohérence individuelle entre le vecteur candidat et le signal rétrodiffusé sur ledit axe de calcul est égale à la corrélation entre le signal attendu sur l'axe de redondance et le spectre du signal rétrodiffusé sur l'axe de redondance Xi; -la cohérence individuelle entre le vecteur candidat et le signal rétrodiffusé sur ledit axe de calcul est calculée par une méthode des moindre carrés; - les axes de calcul sont les axes de visée et en ce qu'on obtient les vecteurs candidats en : • déterminant un volume d'incertitude à partir d'au moins un vecteur vitesse de l'aéronef calculé antérieurement, • découpant le volume d'incertitude en vecteurs candidats; e volume d'incertitude est calculé au moyen d'un filtre de 20 Kalman; - le volume d'incertitude dépend d'informations sur la dynamique de l'aéronef; -la cohérence individuelle entre un vecteur candidat et le signal rétrodiffusé sur un axe de visée est calculée à partir d'un seuil de 25 puissance attendu Ps et de la densité spectrale de puissance Ri du spectre du signal rétrodiffusé sur l'axe d'indice i à la fréquence de la projection du vecteur candidat d'indice j sur l'axe d'indice i à partir de la formule suivante : N Ri COHij = fJ exp(ù PS) , - au moins un vecteur vitesse calculé antérieurement est 30 calculé en mettant en oeuvre le procédé selon l'une quelconque des revendications 1 à 14; Le procédé selon l'invention permet d'améliorer la fiabilité de la mesure de la vitesse d'un aéronef lorsque celui-ci se trouve dans un environnement inhomogène en termes de vitesse. Plus particulièrement, le procédé selon l'invention permet d'améliorer la fiabilité de la mesure de la vitesse d'un aéronef dans un environnement nuageux ou lorsqu'au moins un axe de visé est orienté vers le sol. Ce procédé permet en outre d'augmenter la disponibilité de la mesure de vitesse dans un environnement inhomogène en termes de vitesse. II ne dégrade pas la précision sur la mesure de vitesse lorsque l'environnement est homogène. II permet en outre au système utilisateur (pilote ou système de pilotage automatique) recevant l'information sur la vitesse de l'aéronef de choisir de se servir de l'information de vitesse ou non, selon le niveau de fiabilité de l'information qui lui est transmise. Le procédé selon le second mode de réalisation de l'invention permet d'améliorer la fiabilité de la mesure de la vitesse d'un aéronef dans un environnement perturbé par des perturbations du type vortex. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels : • la figure 1 est une représentation schématique d'un aéronef se 20 déplaçant dans une masse d'air, • les figures 2a, 2b, 2c et 2d représentent respectivement des exemples de spectres de signaux rétrodiffusés en fonction de leur fréquence selon un premier, un deuxième, un troisième et un quatrième axes de visée, • la figure 3a est une représentation schématique du gain de la 25 fonction de focalisation du faisceau laser en fonction de la distance à l'aéronef; la figure 3b est une représentation schématique d'un exemple de coefficient de rétrodiffusion d'une masse d'air en fonction de sa distance à un aéronef, • la 4a est une représentation schématique d'un spectre d'un signal 30 rétrodiffusé en fonction de sa fréquence, sur un axe de visée sur lequel il existe une perturbation de type turbulence, • la figure 4b est une représentation schématique d'un spectre d'un signal rétrodiffusé en présence d'une couche nuageuse portée par un vent plus fort, • la figure 5 est une représentation schématique de deux vecteurs vitesse candidats obtenus à partir de 3 axes de mesure ainsi que de leurs projections sur un axe de visée supplémentaire; sur la figure 5, on a également représenté les spectres des signaux rétrodiffusés représentés sur les figures 2a à 2d, en fonction de la vitesse, sur les axes de visée. D'une figure à l'autre, les mêmes éléments sont repérés par les mêmes références. La figure 1 illustre un aéronef 1 équipé d'un laser 2 émettant continûment des faisceaux, appelés faisceaux émis 4 à une longueur d'onde X. L'aéronef est animé d'une vitesse V relativement à une masse d'air. Le procédé de mesure de la vitesse est basé sur une technique d'anémométrie laser Doppler. A cet effet, on focalise les faisceaux émis 4 sur chacun des axes Xi dans une zone de focalisation ZOF;. Une zone de focalisation est choisie de manière à ce que le mouvement de la masse d'air dans cette zone soit indépendant de celui de l'avion ou suive un modèle prédéfini. La zone de focalisation ZOF; sur un axe de visée est une zone située à une distance moyenne de l'aéronef appelée distance de focalisation d; et dont l'épaisseur, sur l'axe de visée, est appelée épaisseur de la zone de focalisation Ti. Sur la figure 1, on a représenté les zones de focalisation ZOF1, ZOF2, ZOF3, ZOF4. Pour plus de clarté, on n'a représenté sur la figure 1 que les épaisseurs TI et T2 des axes XI et X2 ainsi que les distances de focalisations associées di et d2. La distance de focalisation est choisie en fonction de la taille et de la dynamique de l'aéronef. En variante, on ne focalise pas les faisceaux laser. Le procédé utilise généralement la même optique pour l'émission et la réception mais on peut tout à fait envisager d'utiliser des optiques d'émission et de réception différentes. Les faisceaux sont émis par un unique laser mais on peut envisager d'envoyer des faisceaux provenant chacun d'un laser différent. On envoie des faisceaux 4 selon trois premiers axes de visée X1, X2, X3 et selon au moins un axe supplémentaire X4. Dans le mode de réalisation représenté sur les figures 1 et 5, on a utilisé uniquement 4 axes de visée. Dans un souci de simplicité, on a représenté les 4 axes de visée sur un espace à deux dimensions mais ces axes ne sont pas coplanaires. Chaque axe de visé Xi est orienté. Les axes de visée ne sont pas confondus entre eux. i est l'indice d'un axe de visée Xi. i= 1,...,N. N est le nombre d'axes de visée Xi utilisés dans le procédé. Les aérosols présents dans l'air sur un axe de visée rétrodiffusent le faisceau émis et produisent un signal rétrodiffusé sur cet axe de visée. On capte les signaux rétrodiffusés au moyen d'un détecteur 5 lors d'une étape d'acquisition des signaux rétrodiffusés. Sur les figures 2a, 2b, 2c, 2d, on a représenté un exemple de spectres de signaux rétrodiffusés sur le premier axe de visée X1 (figure 2a), sur le second axe de visée X2 (figure 2b), sur le troisième axe de visée X3 (Figure 2c) ainsi que sur le quatrième axe de visée X4 (figure 2d). Le spectre FTi(f) du signal rétrodiffusé sur l'axe de visé Xi s'exprime en fonction de la fréquence f de ce signal. Il s'agit par exemple de la transformée de Fourier du signal. Lorsque le signal rétrodiffusé est issu d'aérosols se déplaçant tous à la même vitesse, le signal rétrodiffusé présente une fréquence unique. Le spectre de ce signal présente une unique raie de puissance importante, comme visible sur les figures 2b, 2c et 2d. Toutefois, comme représenté sur la figure 2a, lorsque plusieurs types d'aérosols présentant des vitesses différentes sont présents sur l'axe de visée, l'énergie du signal rétrodiffusé est répartie sur plusieurs fréquences. Avec les procédés de l'art antérieur utilisant l'information provenant uniquement de trois axes de visée, la composante du vecteur vitesse sur l'axe de visée concerné a plusieurs valeurs possibles et on ne sait pas identifier quelle est la composante utile. C'est par exemple le cas lorsque la masse d'air présente dans la zone de focalisation forme un vortex dans lequel la vitesse des aérosols est répartie sur un intervalle de vitesse. On a représenté schématiquement sur la figure 4a un exemple de spectre du signal issu de la rétrodiffusion par une perturbation de type vortex. La majorité de la puissance du signal est étalée sur une plage de fréquence finie. Il n'est pas possible d'identifier avec précision la fréquence du signal. On ne peut pas en déduire de façon fiable et précise une composante du vecteur vitesse sur l'axe de visée. Sur la figure 4b, on a représenté schématiquement un spectre d'un signal rétrodiffusé en présence d'une couche nuageuse située à distance de l'aéronef et portée par un vent plus fort que la masse d'air située dans la zone de focalisation. L'intensité du signal est répartie sur deux gaussiennes de même hauteur. II n'est pas possible de reconnaitre le signal rétrodiffusé par la couche nuageuse du signal utile sur des critères liés à l'allure spectrale. Le même problème se pose lorsque l'axe de visée est dirigé vers le sol et qu'un écho sol vient s'ajouter au signal utile. Pour un spectre présentant deux raies d'intensité importante, comme représenté sur la figure 2a, une solution serait de discriminer la raie utile (issue de la masse d'air dans la zone de focalisation) de la raie parasite (issue du nuage ou du sol) en se basant sur la valeur du spectre ou sur la fréquence de celle-ci mais ces méthodes ne sont pas fiables. En effet, on ne sait pas déterminer lequel du nuage ou de la masse d'air circule le plus vite. De plus, le gain de la fonction de focalisation GF(x) du faisceau incident ne présente pas une sélectivité spatiale suffisante pour que le signal issu de la rétrodiffusion par des porteurs situés en dehors de la zone de focalisation soit suffisamment atténué pour ne pas être capté par le détecteur. Sur la figure 3a, on a représenté le gain de la fonction de focalisation GF (x) sur un axe de focalisation en fonction de la distance x à un aéronef. Le gain de la fonction de focalisation GF (x), comme visible sur la figure 3a, est une gaussienne maximum à la distance de focalisation. La puissance IT du signal rétrodiffusé par une masse d'air située à une distance x de l'aéronef est donnée par : IT(x) = GF(x) *p(x) I3(x) est le coefficient de rétrodiffusion de la masse d'air à la distance x. Le coefficient de rétrodiffusion p(x) d'une masse d'air dépend des caractéristiques des traceurs et de leur nombre dans la masse d'air. Sur la figure 3b, on a représenté schématiquement un exemple de variation du coefficient de rétrodiffusion de la masse d'air entourant un avion. La masse d'air a un faible coefficient de rétrodiffusion autour de la distance d de focalisation, et un plateau ayant un fort coefficient de rétrodiffusion à distance importante de l'aéronef. Ce profil est caractéristique de la présence d'un nuage sur l'axe de visé concerné à une distance importante d'un aéronef. Même si un nuage ou le sol qui présentent des coefficients de rétrodiffusion importants se trouvent en dehors de la zone de focalisation, à une distance où la fonction de focalisation est faible, la contribution (en intensité) au signal rétro-diffusé est non négligeable. Il n'existe pas de relation biunivoque entre l'intensité du signal rétro-diffusé et la distance aux traceurs contributeurs. On ne peut pas se baser sur la puissance du spectre pour différencier deux raies d'un spectre du type de celui de la figure 2a. Ainsi, dans le procédé selon l'invention, on se propose d'utiliser la redondance apportée par l'acquisition d'un signal rétrodiffusé sur au moins ~o un axe de visée supplémentaire, c'est-à-dire sur au moins un quatrième axe de visée, pour calculer le vecteur vitesse de l'aéronef. Dans la suite de la description, on s'intéressera tout d'abord au cas où on utilise un seul axe de visée supplémentaire puis on généralisera au cas où on a plusieurs axes de visée supplémentaires. 15 Après acquisition des signaux sur les 4 axes de visée, on calcule le spectre de chaque signal rétrodiffusé sur chacun de ces axes. On calcule les vecteurs vitesse candidats ÿc et leurs cohérences respectives qui sont fonction des projections des vecteurs candidats sur au moins un axe de visée, appelé axe de calcul X1, X2, X3, X4 et en fonction des spectres 20 FT1(V), FT2(V), FT3(V), FT4(V) des signaux rétrodiffusés sur le (ou les) axe(s) de calcul respectifs X1, X2, X3, X4. Le vecteur vitesse V de l'aéronef est choisi parmi les vecteurs vitesse candidats ÿc en faisant une analyse de vraisemblance de leurs cohérences COH. On notera que les spectres de puissance en vitesse FT1(V), 25 FT2(V), FT3(V), FT4(V) des signaux rétrodiffusés sur les axes de visée respectifs X1, X2, X3, X4 sont directement proportionnels aux spectres FT1(f), FT2(f), FT3(f), FT4(f).représentés en fréquence sur les figures 2a à 2d. 30 Dans un premier mode de réalisation selon l'invention, le vecteur vitesse V est choisi en suivant le procédé décrit ci-après. On sélectionne sur chacun des axes de visée X1, X2, X3, X4, une ou plusieurs raies. Une raie est identifiée lorsqu'elle vérifie un critère de sélection de raie prédéterminé. Un critère de sélection de raie peut consister à identifier une raie lorsque la 35 valeur du spectre du signal dépasse un seuil S prédéterminé. Le spectre est filtré par un opérateur adapté au signal attendu, par exemple, par une gaussienne dont les caractéristiques correspondent à la turbulence attendue de la masse d'air. Une première raie est identifiée à la fréquence à laquelle la valeur du spectre filtré est maximale et supérieure au seuil S. La raie suivante est déterminée par le même procédé après avoir éliminé du spectre une plage de fréquence se trouvant autour de la fréquence de la première raie identifiée. Le processus de détection de raies s'arrête lorsqu'il n'y a plus de maximum dépassant le seuil S. D'autres principes de détection de raies sont possibles afin de s'adapter au mieux des caractéristiques du matériel. Une raie est par exemple identifiée lorsqu'un pourcentage minimum de l'énergie du spectre est confiné dans un intervalle de fréquence restreint. Sur les exemples de spectres représentés sur les figures 2a à 2d, on obtient deux raies identifiées R11 et R12 sur le premier axe XI dont les valeurs respectives du spectre FTRI1, FTR12 (ou puissance transportée à ces fréquences) sont supérieures au seuil prédéterminé S. Ces raies sont centrées sur des fréquences de raies f1 1, f12 qui correspondent respectivement à des composantes de vitesse V1 1 et V12 sur le premier axe (voir figure 5.) On associe un écart type E11 et E12 aux raies identifiées. On a identifié une unique raie R21, R31, R41 sur chacun des deuxième X2, troisième X3 et quatrième X4 axes de visée, centrées sur des fréquences de raie respectives f2l, f3l, f41 qui correspondent respectivement à des composantes de vitesse V21, V31 et V41. Les valeurs du spectre pour ces raies sont respectivement FTR21, FTR31, FTR41. On associe aux raies identifiées des écarts types respectifs E21, E31 et E41 qui sont fonction de la forme des raies correspondantes. On sélectionne trois axes, appelés axes de base. Les axes de base sont sélectionnés à partir d'un critère de sélection prédéterminé. Les axes de visée non sélectionnés sont appelés axes de redondance. A partir des axes de base, on construit des vecteurs vitesse candidats dont les composantes sur les trois axes de base sont les composantes en vitesse associées aux raies identifiées sur ces axes de base. Autrement dit, on construit à partir des raies identifiées sur les axes de base les vecteurs à 3 dimensions possibles qui présentent sur chaque axe de base, une composante en vitesse égale à la composante en vitesse associée à une raie identifiée sur l'axe de base considéré. Une raie sur chacun des 3 axes de base conduit à un seul vecteur candidat. Si un axe possède 2 raies, et deux axes possèdent une raie, on obtient deux vecteurs candidats. Si deux axes possèdent 2 raies et un axe possède une raie, on obtient 4 vecteurs candidats et ainsi de suite. De préférence, le critère de sélection des trois axes de base consiste à choisir les axes de visée dont le spectre présente le moins de raies de sorte à produire un nombre minimum de vecteurs vitesse candidats ce qui permet de limiter le coût du procédé en termes de calculs. Comme représenté sur la figure 5, en prenant les exemples de spectres représentés sur les figures 2a, 2b, 2c, à partir des raies identifiées R11, et R12 sur le premier axe de visée XI, de la raie R21 sur X2 et de la raie R31 sur X3, on construit deux vecteurs vitesse candidats Vcl et Tlc2. Les axes XI et X2 sont situés dans le plan de la feuille contrairement aux axes X3 et X4. Les projections des vecteurs candidats sur les axes X3 et X4 sont orthogonales mais les traits pointillés représentant ces projections sur la figure 5 ne forment pas des angles droits avec ces axes car ils n'appartiennent pas au plan de la feuille. On projette ensuite les vecteurs candidats sur des axes de calcul. Dans le mode de réalisation représenté sur la figure 5, l'axe de calcul est l'axe de redondance X4, on projette les deux vecteurs candidats T-/.cl et T7c.2 sur l'axe X4, on a également représenté des signaux attendus S41(V) et S42(V) sur l'axe de redondance X4, des deux vecteurs vitesse candidats Vcl et Tc2 . On associe en effet à chaque projection du vecteur candidat d'indice j sur un axe de redondance d'indice i un signal attendu Sii(V) sur l'axe de redondance. Le signal attendu sur un axe de redondance Xi est centré sur la projection du vecteur candidat sur ledit axe. Autrement dit la puissance du signal attendu est maximale au niveau de la projection du vecteur candidat sur l'axe de redondance considéré. On peut associer à ce signal attendu une forme et une puissance dépendant de la puissance et de la forme des raies à partir desquelles le vecteur vitesse candidat est construit. Par exemple, comme représenté sur la figure 5, le signal attendu représenté en traits pointillés est une gaussienne centrée autour de la projection du vecteur candidat sur ledit axe supplémentaire. L'écart type Ec du signal attendu est estimé en fonction des écarts types des raies candidates identifiées sur les axes de mesure. On calcule ensuite la cohérence COH des vecteurs candidats. On appelle cohérence individuelle COHi entre un vecteur candidat et le signal rétrodiffusé sur un axe de calcul d'indice i, la grandeur caractérisant la concordance entre le signal attendu, issu de la projection du vecteur candidat sur l'axe de calcul considéré, et le spectre du signal rétrodiffusé sur l'axe de calcul considéré. On donnera par la suite plusieurs exemples non limitatifs de calcul de la cohérence. La cohérence individuelle entre le vecteur vitesse candidat d'indice j et un axe de calcul Xi est par exemple égale à la corrélation CORij entre le signal attendu sur l'axe de calcul et le spectre du signal rétrodiffusé sur l'axe de redondance Xi. La corrélation CORij est alors calculée de la manière suivante : CORij = f Sig (V) * FTi(V)dV Ci-dessus, comme dans la suite du texte, les calculs de cohérence sont réalisés en fonction de la vitesse mais on pourrait de façon équivalente les écrire en fonction de la fréquence f d'un signal rétrodiffusé sur cet axe, cette fréquence étant reliée à la vitesse et à la fréquence de l'oscillateur local par l'équation A. Dans le premier mode de réalisation les axes de calcul étant les axes de redondance, on calcule les cohérences individuelles sur les axes de redondance. Dans l'exemple représenté sur la figure 5, dans lequel on a un axe de redondance X4 la corrélation COR41 est importante car le signal attendu associé à ijcl se superpose avec une raie d'intensité importante du spectre du signal rétrodiffusé sur l'axe X4. A l'inverse, la valeur de la corrélation COR42 est très faible car le signal attendu associé à Vc2 se superpose avec une zone du spectre du signal rétrodiffusé sur l'axe X4 dont l'intensité est très faible. En variante, la cohérence individuelle entre un vecteur candidat et un axe de calcul est calculée par la méthode des moindres carrés à partir de la projection du vecteur candidat d'indice j sur l'axe de calcul d'indice i et le spectre du signal rétrodiffusé sur cet axe. Dans un premier exemple, on calcule la somme SCji des carrés des écarts, appelée somme des résidus, entre le signal attendu du vecteur candidat d'indice j et le spectre du signal rétrodiffusé sur l'axe de calcul d'indice i. SC4 j = f (S4j (V) ù FT4(V))2 dV 0 La cohérence individuelle est alors : sc,,u COHj = e sco où SCO représente le résidu attendu ou son estimée, il s'agit d'une valeur prédéterminée dépendant de l'application. II est par exemple déterminé par un compromis entre les performances représentées par la probabilité de détection de la vitesse et la probabilité de fausses alarmes. Plus on choisit une probabilité de fausse alarmes faible, plus la probabilité de détection de la vitesse est faible. Dans un second exemple moins coûteux que le premier en termes de calcul on calcule la somme des résidus à partir de la fréquence (ou vitesse) de la projection vecteur candidat d'indice j sur l'axe d'indice i et de la fréquence moyenne (ou vitesse moyenne) du spectre du signal rétrodiffusé sur l'axe de calcul d'indice i. La somme des résidus est égale au carré de la différence entre la fréquence (ou vitesse) de la projection vecteur candidat d'indice j sur l'axe d'indice i et la fréquence (ou vitesse) moyenne du signal rétrodiffusé sur l'axe d'indice i. La cohérence individuelle est alors : (fùfpi)2 COHj = e Fso fi est fréquence moyenne du spectre du signal rétrodiffusé sur axe d'indice i fpj est la fréquence de la projection du vecteur candidat d'indice j sur l'axe d'indice i. Fs0 est le résidu attendu exprimé en hertz (ou en m.s-1 si on fait le calcul en termes de vitesse.) Dans le cas où on a un unique axe de redondance, la cohérence individuelle est égale à la cohérence ou cohérence globale COHj du vecteur candidat d'indice j. On fait une analyse de vraisemblance à partir des cohérences COH de chaque vecteur vitesse candidat afin de déterminer le vecteur vitesse de l'aéronef, appelé vecteur vitesse mesuré. Le vecteur vitesse mesuré est le vecteur qui a le maximum de vraisemblance. Avant l'analyse de vraisemblance, on rejette les vecteurs candidats qui ont une cohérence inférieure à un seuil de cohérence prédéterminé. Si tous les vecteurs candidats sont rejetés, la mesure est déclarée invalide ou une mesure dégradée est fournie avec un indicateur de défaut adapté. Sinon, pour les vecteurs candidats non rejetés, on fait une analyse de vraisemblance classique dont un exemple non limitatif est donné ci-après. Lorsqu'un vecteur candidat T7.c a une cohérence COH plus élevée que les autres vecteurs candidats et que l'écart entre la cohérence COH de ce vecteur vitesse candidat et la cohérence COH la plus proche est supérieure à un seuil prédéterminé, le vecteur candidat dont la cohérence COH est la plus élevée est le vecteur vitesse dont la vraisemblance est maximale. Lorsque l'écart entre la cohérence COH la plus élevée et la cohérence COH la plus proche, est inférieure à un seuil d'écart de cohérence prédéterminé, le vecteur vitesse mesuré résulte de la moyenne des vecteurs candidats suivants : le vecteur vitesse candidat qui a la cohérence la plus élevée et les vecteurs vitesse candidats qui ont une cohérence s'écartant de la cohérence la plus élevée d'une valeur inférieure au seuil d'écart de cohérence. La moyenne est pondérée en fonction de la cohérence COH de chacun desdits vecteurs vitesse candidats sélectionnés. En variante, on ne rejette pas de vecteurs candidats, on fait 25 l'analyse de vraisemblance pour tous les vecteurs candidats. Dans le mode de réalisation représenté sur la figure 5, le vecteur vitesse mesuré est le premier vecteur vitesse candidat Tic1. Avantageusement, on pondère la cohérence des vecteurs candidats en fonction de critères supplémentaires pris isolément ou de façon 30 combinée avant d'effectuer l'analyse de vraisemblance. Par exemple, on tient compte des vecteurs vitesse mesurés à un ou plusieurs instants précédents en pondérant, préalablement à l'analyse de vraisemblance, la cohérence COH des vecteurs vitesse candidats en fonction d'un ou de plusieurs vecteurs mesurés antérieurement. A cet effet, 35 on pondère par exemple la cohérence d'un vecteur candidat en fonction de l'écart entre celui-ci et le (ou les) vecteur(s) mesuré(s) antérieurement. On favorise un vecteur candidat qui est proche des vecteurs vitesse mesurés à l'instant précédent par rapport à un vecteur candidat éloigné en multipliant sa cohérence par un facteur de pondération plus important. Dans un autre exemple, on pondère la cohérence d'un vecteur candidat au moyen d'un facteur F de pondération dépendant d'une estimation du vecteur vitesse Î/e à l'instant donné réalisée à partir d'un ou de plusieurs vecteurs vitesse mesurés antérieurement et éventuellement d'un intervalle de confiance I attribué à l'estimé. On estime par exemple, le vecteur vitesse à l'instant donnée au moyen d'un filtre de Kalman. (VcùVe)2 Le facteur de pondération est par exemple égal à F = e Dans le cas d'un traitement différé, on peut en outre pondérer la cohérence des vecteurs candidats en fonction des vecteurs vitesse mesurés postérieurement. On peut en outre se baser sur des vecteurs vitesse calculés à partir d'autres types d'informations (par exemple, des informations inertielles ou des informations issues de satellites) pour pondérer les cohérences. On associe un indicateur de qualité au vecteur vitesse mesuré V . Cet indicateur de qualité correspond à la probabilité que la précision et/ou l'intégrité de la mesure de vitesse soient dans un intervalle de performance prédéterminé. L'indicateur de qualité est, par exemple, basé sur la valeur de la cohérence du vecteur candidat sélectionné. Plus la cohérence du vecteur candidat est importante, meilleur est l'indicateur de qualité. Dans le premier mode de réalisation, lorsqu'il y a plus d'un axe de redondance, les axes de calcul sont les axes de redondance. On projette les vecteurs candidats sur tous les axes de redondance. On calcule la cohérence COHji individuelle des vecteurs vitesse candidats sur chacun des axes de redondance j. La cohérence COHj d'un vecteur candidat, appelée cohérence globale est le produit des cohérences individuelles. Dans un deuxième mode de réalisation de l'invention, on sélectionne, comme dans le premier mode de réalisation, des raies sur chacun des axes de visée X1, X2, X3, X4. On construit ensuite tous les ensembles possibles de raies à m dimensions. La dimension M est au moins égal à 4 et égal au nombre d'axes présentant au moins une raie. Chaque raie d'un ensemble est prise sur un axe de visée différent des autres raies de l'ensemble. On cherche à former des ensembles de N raies (correspondant à une raie sur chaque axe) mais si un axe des dépourvu de raies, on cherche à former toutes les combinaisons de N-1 raies. Ce mode de réalisation n'est pas réalisable si moins de quatre axes présentent au moins une raie identifiée. Dans ce cas, la mesure de vitesse est indisponible. Dans le cas où on a 4 axes de visée dont 3 comprennent une raie et 1 comprend un raie, on forme deux ensembles de raies à 4 dimensions. A partir de chaque ensemble de raies, on construit un vecteur candidat à 3 dimensions par une méthode classique de résolution aux moindres carrés. On passe ainsi d'un espace à m dimensions à un vecteur à 3 dimensions qui minimise les écarts entre l'ensemble M et le vecteur à 3 dimension. On utilise par exemple une méthode classique de résolution aux moindres carrés. Dans ce mode de réalisation, les axes de calcul sont les axes de visée. On projette les vecteurs candidats sur tous les axes de visée et on attribue éventuellement, à chaque projection, un signal attendu. On calcule ensuite la cohérence COH des vecteurs candidats qui sont les produits, effectués sur tous les axes de visée, des cohérences individuelles des vecteurs sur chaque axe de visée. Les cohérences individuelles sont calculées de la même manière que dans le premier mode de réalisation. Le vecteur vitesse mesuré est ensuite obtenu en faisant une analyse de vraisemblance éventuellement précédée d'une pondération des cohérences similaire à l'analyse effectuée dans le premier mode de réalisation. Dans un troisième mode de réalisation selon l'invention, on se base sur le fait que le vecteur vitesse de l'avion ne varie pas brusquement. Les vecteurs vitesse candidats ÿc à un instant donné t se trouvent dans un volume d'incertitude dont le contour est défini en fonction du vecteur vitesse calculé à l'instant précédent t - At. On estime par exemple le volume d'incertitude à partir d'au moins une mesure de vitesse précédente par propagation, par exemple à partir d'un filtre de Kalman. En variante, l'estimation du volume d'incertitude est réalisée en outre à partir d'informations sur la dynamique du porteur (par exemple des informations inertielles qui fournissent des mesures peu bruitées.) On obtient ainsi une meilleure estimée du volume d'incertitude. On découpe le volume d'incertitude en vecteurs vitesse candidats au moyen par exemple de l'une des deux méthodes classiques suivante : - tirage aléatoire dans le volume d'incertitude (méthode dite � particulaire �.) Cette méthode conduit à un très grand nombre de vecteurs candidats. - espacement régulier (méthode à grille). L'espacement des vecteurs candidats dépend de la turbulence de l'air. Par exemple, on choisit la turbulence minimale attendue afin de bien couvrir le domaine d'incertitude. Cette méthode conduit à un moins grand nombre de vecteurs candidats que la précédente. Dans ce mode de réalisation les axes de calcul sont les axes de visée. On projette les vecteurs candidats sur tous les axes de visée. On calcule une cohérence COH de chaque vecteur candidat. La cohérence d'un vecteur candidat est égal au produit effectué sur les N axes de visée, des cohérences individuelles. La cohérence individuelle COHij entre un vecteur candidat d'indice j et spectre du signal rétrodiffusé un axe de visée d'indice i est calculée à partir du spectre du signal rétrodiffusé sur l'axe i et de la projection du vecteur candidat j sur cet axe i. Elle est donnée par la formule suivante. N Ri COHij = exp(--) Ps Ps est seuil de puissance attendu. Ri est la densité spectrale de puissance du signal rétrodiffusé sur l'axe d'indice i à la fréquence (ou vitesse) de la projection du vecteur d'indice j sur l'axe d'indice i. On fait ensuite une analyse de vraisemblance, avec éventuellement une pondération préalable des cohérences, similaire à celle qui a été effectuée dans le premier mode de réalisation pour déterminer le vecteur vitesse mesuré. Avantageusement, le troisième mode de réalisation peut être initialisé par un des deux premiers modes de réalisation. Dans ce cas, les mesures de vitesse effectuées aux instants précédents sont calculées à partir d'un des procédés selon le premier ou le deuxième mode de réalisation. Ces mesures initialisent le filtre de propagation pour l'estimation du volume d'incertitude à un instant donné. Le troisième mode de réalisation n'utilisant pas de seuil sur les spectres des signaux sur les différents axes, fonctionne avec des signaux plus faibles mais au prix d'une puissance de calcul plus grande, que les premiers modes de réalisation. La puissance de calcul est d'autant plus grande que l'estimée précédente est mauvaise. REVENDICATIONS1. Procédé de mesure de la vitesse d'un aéronef (1) à un instant donné, par anémométrie laser Doppler, ledit procédé comprenant les étapes suivantes : • acquisition d'un signal rétrodiffusé sur trois premiers axes de visée (XI ,X2, X3), • calcul du spectre du signal rétrodiffusé sur chacun des trois premiers axes de visée (X1, X2, X3), caractérisé en ce que le procédé comprend en outre les étapes suivantes : • acquisition d'un signal rétrodiffusé supplémentaire sur au moins un axe de visé supplémentaire (X4) non confondu avec un des trois premiers axes de visée (XI, X2, X3), • calcul du spectre du signal rétrodiffusé sur au moins un axe de visée supplémentaire (X4), • calcul de plusieurs vecteurs vitesse candidats (T-/>c ), • calcul des projections des vecteurs candidats sur au moins un axe de visée (X1, X2, X3, X4), appelé(s) axe(s) de calcul, • calcul de la cohérence des vecteurs candidats (T7c) en fonction de leur(s) projection(s) sur l'(ou les) axe(s) de calcul et en fonction du (ou des) spectre(s) du signal (ou des signaux) rétrodiffusé(s) sur le(s) axe(s) de calcul (X1, X2, X3, X4), • choix du vecteur vitesse mesuré (V) parmi les vecteurs vitesse candidats (17c) en faisant une analyse de vraisemblance à partir des cohérences desdits vecteurs vitesse candidats (Vc ). 2. Procédé de mesure de la vitesse d'un aéronef (1) selon la revendication précédente, caractérisé en ce que l'étape de calcul de la cohérence d'un vecteur candidat comprend les étapes suivantes: • projection du vecteur candidat sur le(s)dit(s) axe(s) de calcul, • calcul de la cohérence individuelle entre le vecteur candidat et le signal rétrodiffusé sur chaque axe de calcul à partir dela projection du vecteur vitesse candidat sur le(s)dit(s) axe(s) de calcul et le spectre du signal rétrodiffusé sur le(s)dit(s) axe(s) de calcul • calcul du produit des cohérences individuelles, s'il y a plusieurs axes de calcul. 3. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications précédentes, caractérisé en ce que l'analyse de vraisemblance consiste à choisir, parmi plusieurs vecteurs candidats : • soit le vecteur vitesse candidat qui a la cohérence la plus élevée, lorsque l'écart entre la cohérence la plus élevée et les cohérences des autres vecteurs vitesse candidats est supérieure à un seuil d'écart de cohérence prédéterminé, • soit une combinaison de plusieurs vecteurs vitesse candidats comprenant le vecteur vitesse candidat dont la cohérence est la plus élevée ainsi que au moins un vecteur dont la cohérence présente, avec la cohérence la plus élevée, un écart inférieur à un seuil d'écart de cohérence prédéterminé. 4. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on attribue un indicateur de confiance au vecteur vitesse mesuré (V ). 5. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications précédentes, caractérisé en ce que préalablement à l'analyse de vraisemblance on rejette le(s) vecteur(s) candidat(s) qui a (ou ont) une cohérence inférieure à un seuil de cohérence prédéterminé. 6. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications précédentes, caractérisé en ce que préalablement à l'analyse de vraisemblance on pondère la cohérence des vecteurs candidats en fonction d'au moins un vecteur vitesse mesuré à un instant précédant l'instant donné.35 7. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications précédentes, caractérisé en ce que préalablement à l'analyse de vraisemblance on pondère la cohérence des vecteurs candidats en fonction d'une estimation du vecteur vitesse réalisée à partir d'au moins un vecteur vitesse mesuré à un instant précédant l'instant donné et éventuellement en fonction d'un intervalle de confiance attribué à l'estimation du vecteur vitesse à l'instant donné. 8. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications précédentes, caractérisé en ce que préalablement à l'analyse de vraisemblance on pondère la cohérence des vecteurs candidats en fonction d'informations inertielles. 9. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 2 à 8, caractérisé en ce qu'il comprend une étape d'identification de raies, à partir d'un critère de sélection de raies prédéterminé, sur les spectres des signaux rétrodiffusés sur des axes de visée. 10. Procédé de mesure de la vitesse d'un aéronef (1) selon la revendication précédente, caractérisé en ce qu'on détermine les vecteurs candidats en: • sélectionnant 3 axes, appelés axes de base, parmi les axes de visée (X1, X2, X3, X4), I' (ou les) axe(s) non sélectionné(s) étant appelé(s) axe(s) de redondance, • construisant, à partir des raies identifiées sur les axes de base, de vecteurs vitesse candidats (Pc) dont les composantes sur les trois axes de base sont égales aux composantes des raies identifiées sur les trois axes de base, et en ce que I' (ou les) axe(s) de calcul est (ou sont) le(s) axe(s) de redondance. 11. Procédé de mesure de la vitesse d'un aéronef (1) la revendication 9, caractérisé en ce que I' (ou les) axe(s) de calcul est (ou sont) le(s)axe(s) de visée et en ce en ce qu'on construit les vecteurs candidats en: • construisant des ensembles de raies présentant une dimension (M) au moins égale à 4 et égale au nombre d'axes de visée sur lesquels on a identifié au moins une raie, chaque raie d'un ensemble est prise sur un axe de visée différent des autres raies dudit ensemble, • construisant, à partir de chaque ensemble de raies, un vecteur candidat minimisant les écarts avec ledit ensemble. 12. Procédé de mesure de la vitesse d'un aéronef (1) selon la revendication précédente caractérisé en ce qu'on construit un vecteur candidat, à partir d'un ensemble de raies, en mettant en oeuvre une méthode de résolution aux moindres carrés. 13. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 9 à 12, caractérisé en ce qu'il comprend une étape d'attribution à chaque projection du vecteur candidat sur un axe de calcul, d'un signal attendu centré sur la projection du vecteur candidat sur ledit axe de calcul et en ce que la cohérence individuelle d'un vecteur candidat sur un axe de calcul est calculée à partir du spectre du signal rétrodiffusé sur l'axe de calcul et à partir du signal attendu. 14. Procédé de mesure de la vitesse d'un aéronef (1) selon la revendication la revendication précédente, caractérisé en ce que la cohérence individuelle entre le vecteur candidat et le signal rétrodiffusé sur ledit axe de calcul est égale à la corrélation entre le signal attendu sur l'axe de calcul et le spectre du signal rétrodiffusé sur l'axe de redondance Xi. 15. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 9 à 13, caractérisé en ce que la cohérence individuelle entre le vecteur candidat et le signalrétrodiffusé sur ledit axe de calcul est calculée par une méthode des moindre carrés. 16. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 1 à 8, caractérisé en ce que les axes de calcul sont les axes de visée et en ce qu'on obtient les vecteurs candidats en :déterminant un volume d'incertitude à partir d'au moins un vecteur vitesse de l'aéronef calculé antérieurement à l'instant donné, et en découpant le volume d'incertitude en vecteurs candidats. 17. Procédé de mesure de la vitesse d'un aéronef (1) selon la revendication précédente, caractérisé en ce que le volume d'incertitude est calculé au moyen d'un filtre de Kalman. 18. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 16 à 17, caractérisé en ce que le volume d'incertitude dépend d'informations sur la dynamique de l'aéronef. 19. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 15 à 18, caractérisé en ce que la cohérence individuelle entre un vecteur candidat et le signal rétrodiffusé sur un axe de visée est calculée à partir d'un seuil de puissance attendu Ps et de la densité spectrale de puissance Ri du spectre du signal rétrodiffusé sur l'axe d'indice i à la fréquence de la projection du vecteur candidat d'indice j sur l'axe d'indice i à partir de la formule suivante : ~-N7 Ri COHij =11 exp(ù ) Ps 20. Procédé de mesure de la vitesse d'un aéronef (1) selon l'une quelconque des revendications 16 à 19 caractérisé en ce que le(s) vecteur(s) vitesse calculé antérieurement est (ou sont) calculé(s) en mettant en oeuvre le procédé selon l'une quelconque des revendications 1 à 14.