FR 3131008 A1 20230623 FR 2113624 A 20211217 Procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes (2), comprenant des étapes de : estimation (202) d’une matrice de covariance de signaux acquis par le réseau d’antennes (2) ; calcul (204) d’un vecteur propre normalisé de la matrice de covariance ; corrélation (205) du vecteur propre normalisé avec une première table de référence de sorte à produire un premier spectre de corrélation comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée ; corrélation (206) du vecteur propre normalisé avec une deuxième table de référence de sorte à produire un deuxième spectre de corrélation comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée ; construction (312) d’une troisième table de référence par combinaison linéaire de la première table de référence et de la deuxième table de référence avec respectivement une composante de polarisation de l’onde électromagnétique dans la première direction et une composante de polarisation de l’onde électromagnétique dans la deuxième direction ; corrélation (314) du vecteur propre normalisé avec la troisième table de référence de sorte à produire un troisième spectre de corrélation comprenant des indices de corrélation associés à différentes directions d’arrivée ; identification (316) d’une direction d’arrivée associée à un indice de corrélation maximal du troisième spectre de corrélation. Figure pour l’abrégé : Fig. 3 procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes DOMAINE DE L'INVENTION La présente invention concerne un procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes. ETAT DE LA TECHNIQUE On connaît de l’état de la technique des procédés d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes. Certains de ces procédés sont associés à un réseau d’antennes toutes polarisées verticalement ou horizontalement. Or, lorsque la polarisation de l’onde incidente diffère des polarisations verticales ou horizontales, l’estimation de la direction d’arrivée est biaisée, ce qui est d’autant plus vrai lorsque la composante de polarisation croisée des antennes de réception ne peut être maîtrisée, lors notamment d’une intégration sur porteur. Le choix d’un procédé dépend donc de la capacité du réseau antennaire à fonctionner sur l’ensemble des polarisations linéaires voire elliptiques. En particulier, le document EP2435847 décrit un procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes polarisées différemment, qui comprend des étapes conduisant à l’élaboration d’une fonction de coût C(Xj). Cette fonction de coût est le résultat de la somme de deux termes au carré. Cette fonction peut être représentée par une courbe, avec Xj en abscisse et C(Xj) en ordonnée. Il est identifié une valeur X 0 associée à un maximum de cette courbe, la valeur X 0 constituant une estimée de la direction d’arrivée de l’onde électromagnétique. Toutefois, ce procédé n’estime pas conjointement la direction d’arrivée et la polarisation de l’onde incidente. Un but de la présente divulgation est d’obtenir une estimée plus précise de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes. Il est à cet effet proposé, selon un premier aspect, un procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes mis en œuvre par ordinateur, le procédé comprenant des étapes de : - estimation d’une matrice de covariance de signaux acquis par le réseau d’antennes, les signaux étant issus de l’onde électromagnétique, - calcul d’un vecteur propre normalisé de la matrice de covariance, le vecteur propre étant associé à une valeur propre de la matrice de covariance, - corrélation du vecteur propre normalisé avec une première table de référence de sorte à produire un premier spectre de corrélation comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée, la première table de référence comprenant des réponses normalisées du réseau d’antennes à une première onde électromagnétique de référence polarisée dans une première direction respectivement associées à différentes directions d’arrivée de la première onde électromagnétique de référence, - corrélation du vecteur propre normalisé avec une deuxième table de référence de sorte à produire un deuxième spectre de corrélation comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée, la deuxième table de référence comprenant des réponses normalisées du réseau d’antennes à une deuxième onde électromagnétique de référence polarisée dans une deuxième direction respectivement associées à différentes directions d’arrivée de la deuxième onde électromagnétique de référence, la deuxième direction étant différente de la première direction, - construction d’une troisième table de référence par combinaison linéaire de la première table de référence et de la deuxième table de référence avec respectivement une composante de polarisation de l’onde électromagnétique dans la première direction et une composante de polarisation de l’onde électromagnétique dans la deuxième direction, - corrélation du vecteur propre normalisé avec la troisième table de référence de sorte à produire un troisième spectre de corrélation comprenant des indices de corrélation associés à différentes directions d’arrivée, - identification d’une direction d’arrivée associée à un indice de corrélation maximal du troisième spectre de corrélation. Un avantage du procédé proposé et qu’il propose conjointement d’estimer conjointement la direction d’arrivée ainsi que la polarisation de l’onde incidente, ce qui permet, avec le choix du réseau d’antennes adéquat, une intégration de ce dernier très près du véhicule porteur. Le procédé selon le premier aspect peut comprendre également les caractéristiques suivantes, prises seules ou en combinaison lorsque cela est techniquement possible. De préférence, le procédé selon le premier aspect comprend en outre des étapes de : - sélection d’une direction d’arrivée associée à un indice de corrélation maximal du premier spectre de corrélation ou du deuxième spectre de corrélation, - construction d’une table de recherche en polarisation à partir d'une première réponse normalisée du réseau d’antennes à la première onde électromagnétique de référence, la première réponse normalisée étant associée dans la première table de référence à la direction d’arrivée sélectionnée, et à partir d’une deuxième réponse normalisée du réseau d’antennes à la deuxième onde électromagnétique de référence, la deuxième réponse normalisée étant associée dans la deuxième table de référence à la direction d’arrivée sélectionnée, - corrélation du vecteur propre normalisé avec la table de recherche en polarisation, de sorte à obtenir un quatrième spectre de corrélation comprenant des indices de corrélation associés à différentes valeurs complexes, - ajustement de la composante de polarisation de l’onde électromagnétique dans la deuxième direction à une valeur complexe associée à un indice de corrélation maximal du quatrième spectre de corrélation, puis normalisation d’un vecteur de polarisation formé par les composantes de polarisation de l’onde électromagnétique dans la première direction et la composante de polarisation de l’onde électromagnétique dans la deuxième direction, la construction de la troisième table de référence étant mise en œuvre après la normalisation du vecteur de polarisation. De préférence, le premier spectre de corrélation comprend un premier indice de corrélation maximal et deuxième spectre de corrélation comprend un deuxième indice de corrélation maximal, et la direction d’arrivée sélectionnée à l’étape de sélection est : - une direction d’arrivée associée au premier indice de corrélation maximal dans le premier spectre de corrélation, lorsque le premier indice de corrélation maximal est supérieur au deuxième indice de corrélation maximal, - une direction d’arrivée associée au deuxième indice de corrélation maximal dans le deuxième spectre de corrélation, lorsque le premier indice de corrélation maximal n’est pas supérieur au deuxième indice de corrélation maximal. De préférence, le procédé selon le premier aspect comprend la création d’un espace de recherche comprenant différentes valeurs pour une variable complexe, et dans lequel la construction de la table de recherche en polarisation comprend une combinaison linéaire de la première réponse et de la deuxième réponse avec respectivement la composante de polarisation de l’onde électromagnétique dans la première direction et un poids dépendant du module de la deuxième composante de polarisation de l’onde électromagnétique et d’une des différentes valeurs pour la variable complexe, la combinaison linéaire étant répétée pour chacune des différentes valeurs. De préférence, le procédé selon le premier aspect comprend des étapes de : - mise à jour de la direction d’arrivée sélectionnée à la valeur de la direction d’arrivée identifiée, - mise à jour de l’espace de recherche en un nouvel espace de recherche incluant la composante de polarisation de l’onde électromagnétique dans la deuxième direction du vecteur de polarisation obtenue à l’issue de l’étape de normalisation, puis - répétition des étapes de sélection, construction d’une table de recherche en polarisation, corrélation du vecteur propre normalisé avec la table de recherche en polarisation, ajustement, construction d’une troisième table de référence, corrélation du vecteur propre normalisé avec la troisième table de référence, et identification d’une direction d’arrivée. De préférence, l’espace de recherche forme une spirale dans le plan complexe, et la mise à jour de l’espace de recherche comprend une dilatation de la spirale, de sorte que le nouvel espace de recherche forme une nouvelle spirale dans le plan complexe qui est dilatée par rapport à la spirale. De préférence, le procédé selon le premier aspect comprend en outre une comparaison entre l’indice de corrélation maximal du troisième spectre de corrélation et un seuil prédéfini, ladite répétition n’étant mise en œuvre que si l’indice de corrélation maximal du troisième spectre de corrélation est inférieur au seuil prédéfini. De préférence, la première direction et la deuxième direction sont orthogonales. Il est également proposé, selon un deuxième aspect, une mémoire lisible par ordinateur stockant des instructions exécutables par l’ordinateur pour l’exécution des étapes du procédé selon le premier aspect. Il est également proposé, selon un troisième aspect, un dispositif d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes, le dispositif d’estimation comprenant au moins un processeur configuré pour mettre en œuvre le procédé selon le premier aspect. Il est également proposé, selon un quatrième aspect, un système comprenant : - un réseau d’antennes configuré pour acquérir des signaux issus d’une onde électromagnétique, - un dispositif d’estimation de la direction d’arrivée de l’onde électromagnétique au réseau d’antennes sur la base des signaux, le dispositif d’estimation étant conforme au troisième aspect. Il est également proposé un porteur mobile, tel qu’un aéronef, comprenant un système selon le quatrième aspect. DESCRIPTION DES FIGURES D’autres caractéristiques, buts et avantages de l’invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : La illustre de façon schématique un système selon un mode de réalisation. La est un organigramme d’étapes d’un procédé susceptible d’être mis en mis en œuvre par le système de la . La est un organigramme d’étapes d’un procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes, selon un mode de réalisation. La est un diagramme représentant des valeurs propres et des indices associés, ces valeurs propres étant impliquées dans le procédé de la . La figure 5 représente deux courbes formées par des spectres de corrélation , calculés au cours d’une mise en œuvre du procédé de la . La est un organigramme détaillant un mode de réalisation d’une étape du procédé de la . La représente un espace de recherche pour une variable complexe p, dans le plan complexe. La figure 8a représente une courbe formée par un spectre calculé au cours d’une itération du procédé de la , ce spectre étant fonction de la variable complexe p. La figure 8b représente les deux courbes de la figure 5 ainsi qu’une courbe formée par une autre spectre de corrélation calculés au cours d’une itération du procédé de la . La figure 9a représente une courbe formée par un spectre calculé au cours d’une autre itération du procédé de la , ce spectre étant fonction de la variable complexe p. La figure 9b représente les deux courbes de la figure 5 ainsi qu’une courbe formée par une autre spectre de corrélation calculé au cours d’une autre itération du procédé de la . Sur l’ensemble des figures, les éléments similaires portent des références identiques. Procédé d’estimation de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes (2) mis en œuvre par ordinateur, le procédé comprenant des étapes de : - estimation (202) d’une matrice de covariance de signaux acquis par le réseau d’antennes (2), les signaux étant issus de l’onde électromagnétique, - calcul (204) d’un vecteur propre normalisé ( ) de la matrice de covariance, le vecteur propre étant associé à une valeur propre de la matrice de covariance, - corrélation (205) du vecteur propre normalisé ( ) avec une première table de référence ( ) de sorte à produire un premier spectre de corrélation ( ) comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée, la première table de référence ( ) comprenant des réponses normalisées du réseau d’antennes (2) à une première onde électromagnétique de référence polarisée dans une première direction respectivement associées à différentes directions d’arrivée de la première onde électromagnétique de référence, - corrélation (206) du vecteur propre normalisé ( ) avec une deuxième table de référence ( ) de sorte à produire un deuxième spectre de corrélation ( ) comprenant des indices de corrélation respectivement associés à différentes directions d’arrivée, la deuxième table de référence ( ) comprenant des réponses normalisées du réseau d’antennes (2) à une deuxième onde électromagnétique de référence polarisée dans une deuxième direction respectivement associées à différentes directions d’arrivée de la deuxième onde électromagnétique de référence, la deuxième direction étant différente de la première direction, - construction (312) d’une troisième table de référence ( ) par combinaison linéaire de la première table de référence ( ) et de la deuxième table de référence ( ) avec respectivement une composante de polarisation (P1) de l’onde électromagnétique dans la première direction et une composante de polarisation (P2) de l’onde électromagnétique dans la deuxième direction, - corrélation (314) du vecteur propre normalisé ( ) avec la troisième table de référence ( ) de sorte à produire un troisième spectre de corrélation ( ) comprenant des indices de corrélation associés à différentes directions d’arrivée, - identification (316) d’une direction d’arrivée associée à un indice de corrélation maximal du troisième spectre de corrélation ( ). Procédé selon la revendication 1, comprenant en outre des étapes de : - sélection (210) d’une direction d’arrivée associée à un indice de corrélation maximal du premier spectre de corrélation ( ) ou du deuxième spectre de corrélation ( ), - construction (304) d’une table de recherche en polarisation à partir d'une première réponse normalisée du réseau d’antennes (2) à la première onde électromagnétique de référence, la première réponse normalisée étant associée dans la première table de référence ( ) à la direction d’arrivée sélectionnée, et à partir d’une deuxième réponse normalisée du réseau d’antennes (2) à la deuxième onde électromagnétique de référence, la deuxième réponse normalisée étant associée dans la deuxième table de référence ( ) à la direction d’arrivée sélectionnée, - corrélation (306) du vecteur propre normalisé ( ) avec la table de recherche en polarisation, de sorte à obtenir un quatrième spectre de corrélation ( ) comprenant des indices de corrélation associés à différentes valeurs d’une variable complexe, - ajustement (310) de la composante de polarisation de l’onde électromagnétique dans la deuxième direction (P2) à une valeur complexe (p2) associée à un indice de corrélation maximal du quatrième spectre de corrélation ( ), puis normalisation d’un vecteur de polarisation formé par les composantes de polarisation (P1) de l’onde électromagnétique dans la première direction et la composante de polarisation (P2) de l’onde électromagnétique dans la deuxième direction, la construction (312) de la troisième table de référence ( ) étant mise en œuvre après la normalisation du vecteur de polarisation. Procédé selon la revendication 2, dans lequel le premier spectre de corrélation ( ) comprend un premier indice de corrélation maximal et deuxième spectre de corrélation ( ) comprend un deuxième indice de corrélation maximal, et dans lequel la direction d’arrivée sélectionnée à l’étape de sélection (210) est : - une direction d’arrivée associée au premier indice de corrélation maximal dans le premier spectre de corrélation ( ), lorsque le premier indice de corrélation maximal est supérieur au deuxième indice de corrélation maximal, - une direction d’arrivée associée au deuxième indice de corrélation maximal dans le deuxième spectre de corrélation ( ), lorsque le premier indice de corrélation maximal n’est pas supérieur au deuxième indice de corrélation maximal. Procédé selon l’une des revendications 2 et 3, comprenant la création (302) d’un espace de recherche comprenant différentes valeurs pour la variable complexe, et dans lequel la construction (304) de la table de recherche en polarisation comprend une combinaison linéaire de la première réponse et de la deuxième réponse avec respectivement la composante de polarisation de l’onde électromagnétique dans la première direction et un poids dépendant du module de la deuxième composante de polarisation de l’onde électromagnétique et d’une des différentes valeurs pour la variable complexe, la combinaison linéaire étant répétée pour chacune des différentes valeurs. Procédé selon la revendication 4, comprenant les étapes de : - mise à jour de la direction d’arrivée sélectionnée à la valeur de la direction d’arrivée identifiée, - mise à jour de l’espace de recherche en un nouvel espace de recherche incluant la composante de polarisation (P2) de l’onde électromagnétique dans la deuxième direction du vecteur de polarisation obtenue à l’issue de l’étape de normalisation, puis - répétition des étapes de sélection (303), construction (304) d’une table de recherche en polarisation, corrélation (306) du vecteur propre normalisé ( ) la table de recherche en polarisation, ajustement (310), construction (312) d’une troisième table de référence , corrélation (314) du vecteur propre normalisé ( ) avec la troisième table de référence ( ), et identification (316) d’une direction d’arrivée. Procédé selon la revendication 5, dans lequel l’espace de recherche forme une spirale dans le plan complexe, et la mise à jour de l’espace de recherche comprend une dilatation de la spirale, de sorte que le nouvel espace de recherche forme une nouvelle spirale dans le plan complexe qui est dilatée par rapport à la spirale. Procédé selon l’une des revendications 5 et 6, comprenant en outre une comparaison entre l’indice de corrélation maximal du troisième spectre de corrélation et un seuil prédéfini, ladite répétition n’étant mise en œuvre que si l’indice de corrélation maximal du troisième spectre de corrélation est inférieur au seuil prédéfini. Procédé selon l’une des revendications précédentes, dans lequel la première direction et la deuxième direction sont orthogonales. Mémoire (8) lisible par ordinateur stockant des instructions exécutables par l’ordinateur pour l’exécution des étapes du procédé selon l’une des revendications précédentes. Dispositif d’estimation (4) de la direction d’arrivée d’une onde électromagnétique à un réseau d’antennes (2), le dispositif d’estimation (4) comprenant au moins un processeur (6) configuré pour mettre en œuvre le procédé selon l’une des revendications 1 à 8. Système (1) comprenant : - un réseau d’antennes (2) configuré pour acquérir des signaux issus d’une onde électromagnétique, - un dispositif d’estimation (4) de la direction d’arrivée de l’onde électromagnétique au réseau d’antennes (2) sur la base des signaux, le dispositif d’estimation (4) étant conforme à la revendication 10. Porteur mobile, tel qu’un aéronef, comprenant un système (1) selon la revendication 11.