FR 3091357 A1 20200703 FR 1874398 A 20181231 Description Titre de l'invention : PROCEDE DE DETERNIINATION CONJOINTE DE LA DIRECTION D'ARRIVEE ET DE LA PO- LARISATION D'UN SIGNAL Domaine technique 100011 Le domaine technique de l'invention est celui de la détermination de la direction d'arrivée et de la polarisation d'un signal et plus particulièrement de la détermination conjointe de la direction d'arrivée et de la polarisation d'un signal. 100021 La présente invention concerne un procédé de détermination de la polarisation et de la direction d'arrivée et en particulier un procédé de détermination conjointe de la polarisation et de la direction d'arrivée. La présente invention concerne également un dispositif permettant de mettre en oeuvre le procédé selon l'invention. ARRIERE-PLAN TECHNOLOGIQUE DE L'INVENTION 100031 Les algorithmes de détection classiques détectent des émetteurs radioélectriques en étudiant la puissance du signal reçu par un dispositif d'acquisition comportant une ou plusieurs antennes. En effet, en l'absence d'émetteurs, la valeur de la puissance du signal reçu est plus faible qu'en présence d'un émetteur et donc, en choisissant un seuil adapté, on peut détecter la présence d'un émetteur par seuillage de la puissance du signal reçu. 100041 La détection d'un émetteur est souvent suivie d'une étude visant à déterminer les ca- ractéristiques du signal reçu, comme sa direction d'arrivée ou encore sa polarisation, pour caractériser l'émetteur. 100051 11 existe de nombreux algorithmes permettant de déterminer de manière précise la direction d'arrivée d'un signal reçu. Certains, comme l'algorithme dc formation dc voies, ont besoin de la polarisation du signal reçu, d'autres, comme par exemple l'algorithme ESPRIT ou l'algorithme MUSIC, peuvent déterminer la direction d'arrivée quelle que soit la polarisation du signal. Dans tous les cas, l'analyse de la polarisation et l'analyse de la direction d'arrivée sont réalisées successivement et non conjointement, ce qui mobilise plus de ressources de calcul ct prend plus de temps. [0006] Il existe donc un besoin d'un algorithme permettant de déterminer conjointement et de manière précise la direction d'arrivée et la polarisation d'un signal reçu. Résumé de l'invention [0007] L'invention offre une solution aux problèmes évoqués précédemment, en proposant un algorithme capable de déterminer précisément et conjointement la direction d'arrivée et la polarisation d'un signal reçu. 100081 Un premier aspect de l'invention concerne un procédé de détermination d'une pola- 2 risation et d'une direction d'arrivée parmi un ensemble de directions d'arrivée, d'un signal émis par un émetteur à une fréquence comprise dans une gamme de fréquences, le procédé étant mis en oeuvre par un dispositif comportant une unité de calcul et un réseau antennaire comprenant une pluralité d'antennes configurécs pour réaliser l'acquisition en simultané du signal et à recevoir et transmettre en simultané le signal acquis sous forme numérique à l'unité de calcul, le réseau antennaire possédant un repère pour chaque direction d'arrivée de l'ensemble de directions d'arrivée et pour chaque fréquence de la gamme de fréquences, chaque repère comportant un premier axe et un second axe, le procédé comportant les étapes suivantes : pour chaque antenne du réseau antennaire : - recevoir, numériser le signal et transmettre le signal numérisé à l'unité de calcul ; estimer la fréquence du signal à partir d'au moins un signal numérisé ; regrouper les signaux numérisés au sein d'un vecteur de signaux numérisés H pour chaque direction d'arrivée de l'ensemble de directions d'arrivée : - décomposer chaque signal numérisé sur le repère correspondant à la direction d'arrivée et à la fréquence estimée ; - regrouper la décomposition de chaque signal numérisé suivant le premier axe du repère dans un premier vecteur, regrouper la décomposition de chaque signal numérisé suivant le second axe du repère dans un second vecteur et créer une matrice X comportant une première ligne correspondant au premier vecteur et une seconde ligne correspondant au second vecteur ; déterminer la direction d'arrivée et la polarisation du signal à l'aide du vecteur de signaux numérisés H et de la matrice X de chaque direction d'arrivée de l'ensemble de directions d'arrivée. [0009] Grâce à l'invention, la direction d'arrivée et la polarisation du signal reçu sont dé- terminées conjointement par le biais du vecteur de signaux numérisés H et de la matrice X de chaque direction d'arrivée de l'ensemble de directions d'arrivée. [0010] Outre les caractéristiques qui viennent d'être évoquées dans le paragraphe précédent, le procédé selon un premier aspect de l'invention peut présenter une ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. [0011] Selon un mode de réalisation, l'étape de détermination de la direction d'arrivée et de la polarisation du signal comporte les sous-étapes suivantes : - pour chaque direction d'arrivée de l'ensemble de directions (l'arrivée : - calculer la matrice C se définissant comme : C = (XX*)-' XH*HX'; 3 - calculer la trace de la matrice C; déterminer la direction d'arrivée du signal se définissant comme la direction d'arrivée de la matrice C de trace maximale ; déterminer la polarisation du signal à partir de la matrice C correspondant à la direction d'arrivée du signal. [0012] Ainsi, la direction d'arrivée et la polarisation du signal reçu sont déterminées conjointement par le biais de la matrice C de trace maximale. [0013] Selon un premier exemple de réalisation compatible avec le mode de réalisation précédent, l'étape d'estimation de la fréquence du signal est réalisée en calculant la transformée de Fourier d'au moins un signal numérisé. [0014] Ainsi, l'estimation de la fréquence du signal est rapide. [0015] Selon un deuxième exemple de réalisation non compatible avec le premier exemple de réalisation et compatible avec le mode de réalisation précédent, l'étape d'estimation de la fréquence du signal est réalisée en calculant l'intercorrélation des spectrogrammes de plusieurs signaux numérisés. [0016] Ainsi, l'estimation de la fréquence du signal est très précise. [0017] Selon un mode de réalisation compatible avec les modes et exemples de réalisation précédents, l'ensemble des directions d'arrivée comporte une pluralité de directions d'arrivée séparées deux à deux par un pas angulaire donné. [0018] Ainsi, les directions d'arrivée de l'ensemble de directions d'arrivée sont réparties uniformément sur une zone angulaire donnée. [0019] Selon un mode de réalisation compatible avec les modes et exemples de réalisation précédents, la polarisation peut être linéaire ou circulaire. [0020] Un deuxième aspect de l'invention concerne un dispositif permettant la mise en oeuvre du procédé selon un premier aspect de l'invention, caractérisé en ce qu'il comporte : un réseau antennaire comportant une pluralité d'antennes configurées pour recevoir simultanément et à numériser des signaux pour les transmettre à une unité de calcul sous forme de signaux numérisés ; une unité de calcul configurée pour réaliser des fonctions de calcul et de détection sur des signaux numérisés. [0021] Un troisième aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en oeuvre toutes les étapes du procédé selon un premier aspect de l'invention, exceptée l'étape d'acquisition et de numérisation. [0022] Un quatrième aspect de l'invention concerne un support d'enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en oeuvre toutes les étapes du procédé selon un premier aspect de l'invention, exceptée l'étape d'acquisition et de numérisation. [0023] L'invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Brève description des dessins [0024] Les figures sont présentées à titre indicatif et nullement limitatif de l'invention. La figure 1 montre un schéma synoptique d'un procédé selon un premier aspect de l'invention. La figure 2 montre une représentation schématique d'un dispositif selon un deuxième aspect de l'invention. La figure 3 montre une représentation schématique du champ électrique d'une onde polarisée elliptiquement. La figure 4 montre le spectre d'un signal S. [0025] DESCRIPTION DETAILLEE D'AU MOINS UN MODE DE REALISATION DE L'INVENTION [0026] Sauf précision contraire, un même élément apparaissant sur des figures différentes présente une référence unique. [0027] Un premier aspect de l'invention concerne un procédé de détermination de la pola- risation et de la direction d'arrivée d'un signal émis par un émetteur dans un envi- ronnement proche d'un dispositif selon un deuxième aspect dc permettant la mise en oeuvre du procédé selon un premier aspect de l'invention. [0028] Par « environnement proche », on entend une zone circulaire entourant le dispositif selon un deuxième aspect de l'invention. [0029] Par « émetteur », on entend un dispositif qui émet des signaux sous forme d'ondes électromagnétiques. Les ondes électromagnétiques sont par exemple de type radio ou de type radar. [0030] Le signal peut être caractérisé par une fréquence, une direction d'arrivée et une pola- risation. [0031] La polarisation d'un signal peut être de plusieurs types, comme par exemple dc type linéaire, circulaire droite ou circulaire gauche. Tous ces types de polarisation sont des cas particuliers de polarisation elliptique. [0032] [fig.3] montre une représentation schématique du champ électrique d'une onde polarisée elliptiquement. [0033] La pointe du vecteur du champ électrique E de l'onde trace une ellipse caractérisée par un demi-grand axe a, un demi-petit axe b et un angle d'orientation v. [0034] Les polarisations linéaire, circulaire gauche et circulaire droite sont des polarisations elliptiques avec des valeurs particulières de demi-grand axe a et de demi-petit axe b résumées dans le tableau suivant. [0035] [Tableaux 1] Type de polarisation Valeur de a Valeur de b Polarisation linéaire 1 0 Polarisation circulaire gauche 1 1 v1.2 -,72 Polarisation circulaire droite 1 1 ,,,I7 vt2 100361 Dans le contexte de l'invention, déterminer la polarisation P du signal revient à dé- terminer le type de polarisation et ses paramètres, par exemple une polarisation elliptique caractérisée par les paramètres a, b et W. 100371 Le procédé 100 selon un premier aspect de l'invention comporte six étapes 101, 102, 103, 104, 105, 110. 100381 La dernière étape 110 du procédé 100 comporte par exemple quatre sous-étapes 106, 107, 108, 109. 100391 Ifie.11 montre un schéma synoptique illustrant les étapes 101 à 109 du procédé 100 selon un premier aspect de l'invention, selon un exemple de réalisation. 100401 Ifie.21 montre une représentation schématique du dispositif 200 selon un deuxième aspect de l'invention. 100411 Le dispositif 200 comporte une unité de calcul 201 réalisant les étapes 102 à 109 du procédé 100 selon l'exemple de réalisation illustré à la figure 1 et un réseau antennaire 202 réalisant l'étape 101 d'acquisition et de numérisation du procédé 100 selon un premier aspect de l'invention. [0042] Le réseau antennaire 202 comporte N antennes 203, N étant strictement supérieur à 1. [0043] Sur la figure 2, le réseau antennaire 202 présente six antennes 203 qui captent en même temps le signal S en provenance de l'environnement proche, c'est-à-dire que chaque antenne 203 reçoit un signal S différent niais que l'acquisition est déclenchée simultanément pour toutes les antennes 203. [0044] On entend par « antenne » 203, un dispositif comportant un élément antennaire capable de capter des signaux de son environnement et un convertisseur analogique numérique capable de numériser un signal. [0045] Le signal S reçu par chaque antenne 203 est ensuite numérisé de manière à ce que tous les signaux S reçus par les antennes 203 soient numérisés simultanément. Après numérisation, chaque antenne 203 transmet le signal numérisé obtenu après numérisation du signal S à l'unité de calcul 201. [0046] Le signal numérisé obtenu par la n-eme antenne 203 du réseau antennaire 203 est noté S. Ainsi, dans le cas de la figure 2, le signal numérisé obtenu par la troisième antenne 203 du réseau antennaire est noté S3. 6 [0047] Dans la présente demande, « unité de calcul » renvoie à un dispositif possédant des moyens de stockage 204 et des moyens de calculs 205 sur les signaux numérisés. Des moyens de stockage peuvent par exemple être un disque dur, une mémoire SSD ou une mémoire Flash simple ou à grille, locale ou distante. Des moyens de calcul peuvent être un microprocesseur ou un circuit logique programmable comme un FPGA (pour « Field-Programmable Gate Array »). [0048] Après réception des signaux numérisés, l'unité de calcul 201 met en oeuvre les étapes 102 à 109 du procédé 100 selon l'exemple de réalisation illustré à la figure 1. [0049] La deuxième étape 102 du procédé 100 consiste à estimer la fréquence du signal S. [0050] La deuxième étape 102 du procédé 100 peut par exemple être réalisée en appliquant une transformée de Fourier à l'un des signaux numérisés pour obtenir son spectre. [0051] On entend par « spectre d'un signal». l'amplitude du signal en fonction de la fréquence. [0052] [fig.4] montre une représentation schématique du spectre du signal S de fréquence fs. [0053] La deuxième étape 102 du procédé 100 peut également être réalisée en calculant le spectrogramme de chaque signal numérisé, puis en évaluant r intercorrelation entre les spectrogrammes. [0054] On entend par « spectrogramme », un diagramme temps-fréquence indiquant pour chaque intervalle de temps et chaque spectre, l'amplitude et la phase du signal numérisé. [0055] A l'issue de la deuxième étape 102 du procédé 100, la fréquence fs du signal S est déterminée. [0056] La troisième étape 103 du procédé 100 consiste à rassembler les signaux numérisés au sein d'un vecteur H se définissant comme : [0057] [Math.1] TI (S ) (Sn)ni _ N [0058] Par exemple, si le réseau antennaire 202 comporte six antennes 203, comme dans le cas de la figure 2, le vecteur H s'écrit : [0059] [Math.2] 11(S) [S1 52 S3 54 S:r, S6] [0060] Le vecteur H est un vecteur de dimension N. [0061] Les étapes 104 à 107 du procédé 100 sont ensuite réalisées pour chaque direction d'an-h/6C i d'un ensemble de directions d'arrivées prédéfini. [0062] L'ensemble de directions d'arrivées correspond à un nombre fini de directions d'an-h/6C i comprises dans une zone angulaire d'acquisition du réseau antennaire 202. [0063] Dans l'exemple de la figure 2, les antennes 203 sont uniquement capables d'acquérir les signaux dont la direction d'arrivée i est comprise dans un demi-cercle représenté en 7 pointillés, correspondant à une zone angulaire d'acquisition de 180'. [0064] Les directions d'arrivée i de l'ensemble de directions d'arrivée sont par exemple choisies de manière à être espacées d'un même pas angulaire 0 deux à deux. [0065] Sur la figure 2, l'ensemble de directions d'arrivée comporte quatre directions d'arrivée i espacées deux à deux d'un pas angulaire 0 de 45'. [0066] L'ensemble de directions d'arrivée peut par exemple correspondre à 72 directions d'arrivée i correspondant à des directions d'arrivée réparties sur une zone angulaire de 360 et espacées deux à deux d'un pas angulaire de 5'. [0067] L'ensemble de directions d'arrivée peut par exemple correspondre à 360 directions d'arrivée i correspondant à des directions d'arrivée réparties sur une zone angulaire de 360 et espacées deux à deux d'un pas angulaire de 1'. [0068] A chacune des directions d'arrivée i de l'ensemble des directions d'arrivée est associée une pluralité de repères R, chaque repère R correspondant à une fréquence appartenant à une gamme de fréquence prédéfinie dans laquelle est comprise la fréquence fs du signal S. [0069] La gamme de fréquence est par exemple la gamme de fréquences radio ou encore la gamme de fréquences radar. [0070] Ainsi, chaque direction d'arrivée i de l'ensemble des directions d'arrivée est associée à un unique repère R correspondant à la fréquence estimée f, du signal S. Cet unique repère R comporte un premier axe nommé par la suite axe horizontal et un second axe nommé par la suite axe vertical. [0071] La quatrième étape 104 du procédé 100 consiste à décomposer chaque signal numérisé sur le repère R correspondant à la direction d'arrivée i de l'ensemble de directions d'anivée traitée et à la fréquence estimée fs. [0072] A l'issue de la quatrième étape 104 du procédé 100, une décomposition selon l'axe horizontal et une décomposition selon l'axe vertical sont obtenues pour chaque signal numérisé. [0073] Par exemple, si le réseau antennaire 202 comporte six antennes 203, comme dans le cas de la figure 6, six décompositions selon l'axe horizontal SI" S2" SI" 54" Si" Sr," et six décompositions selon l'axe vertical SI v 5,V Slv 54" Sçv S6v sont obtenues. [0074] La cinquième étape 105 du procédé 100 consiste à rassembler les décompositions selon l'axe horizontal au sein d'un premier vecteur d'une part, et les décompositions selon l'axe vertical au sein d'un second vecteur d'autre part, et de créer une matrice X comportant une première ligne constituée du premier vecteur et une seconde ligne constituée du second vecteur. [0075] Par exemple, si le réseau antennaire 202 comporte six antennes 203, comme dans le cas de la figure 2, la matrice X s'exprime ainsi : [0076] 8 [Math.3] [0077] La matrice X est une matrice comportant deux lignes et N colonnes. [0078] La sixième étape 106 du procédé 100 scion l'exemple de réalisation illustré à la figure 1 consiste à calculer la matrice C se définissant comme : [0079] [Math.4] CCfs, LX(fs, OX(fs, i)) IX(fs, (H(S)XH(S) )X(fs, [0080] Avec X't correspondant à la transformée conjuguée de la matrice X. [0081] La matrice C est une matrice comportant deux lignes et deux colonnes. [0082] A l'issue de la sixième étape 106 du procédé 100 selon l'exemple de réalisation illustré à la figure 1, une matrice C est obtenue pour chaque direction d'arrivée i de l'ensemble des directions d'arrivée. Ainsi, dans le cas de la figure 2, quatre matrices C sont obtenues. [0083] La septième étape 107 du procédé selon l'exemple de réalisation illustré à la figure 1 consiste à calculer la trace de chaque matrice C précédemment obtenues. 100841 On rappelle que la trace d'une matrice correspond à la somme de ces éléments diagonaux. 100851 A l'issue de la septième étape 107 du procédé 100, le calcul de la trace de chaque matrice C calculées à la sixième étape 106 du procédé 100 selon l'exemple de réalisation illustré à la figure 1 sont obtenues. Ainsi, dans le cas de la figure 2, on dispose de quatre calculs de trace, correspondant chacun à une matrice C, correspondant elles-mêmes chacune à une direction d'arrivée i de l'ensemble de directions d'arrivée. 100861 La huitième étape 108 du procédé 100 selon l'exemple de réalisation illustré à la figure 1 consiste à déterminer la direction d'arrivée is du signal S parmi les directions d'arrivée i de l'ensemble des directions d'arrivée. 100871 Pour cela, les traces calculées à la septième étape 107 du procédé 100 sont comparées pour déterminer la trace maximale. 100881 La direction d'arrivée i de l'ensemble des directions d'arrivée correspondant à la trace maximale est la direction d'arrivée i5 du signal S. 100891 La neuvième étape 109 du procédé selon l'exemple de réalisation illustré à la figure 1 consiste alors à déterminer la polarisation P du signal S, c'est-à-dire le type de polarisation et ses paramètres, à partir de la matrice C correspondant à la direction d'arrivée is du signal S. 100901 Pour cela, une matrice de polarisation A est définie. En partant du cas général de la polarisation elliptique comme représentée à la figure 3, la matrice A s'écrit : X(fs, Esf, st st C 040 st 1, 1. S \- i S2 S3 S4 S5 S6 9 [0091] [Math.5] A('), a, b) = E(a, b)Q(q)) [0092] Avec : [0093] [Math.6] QN)= [cos (tp) -si sin (q)) cosn(q)(11))1 [0094] [Math 7] E(a, b)= [a ib 100951 Le vecteur FI dépend de la matrice A selon la relation suivante : 100961 IMath.81 H ( S ) =A(4), a, b)X(fs, i) 100971 Or: 100981 IMath.81 C(fs, i) = (X(fs,i)X(fs, ) IX(fs,i) (H(S)'H(S))X(fs,i)' [0099] Donc, en remplaçant H par son expression en fonction de A. la matrice C s'exprime de la manière suivante : [0100] [Math.9] C(fs,i) = (X(fs, i)X(fs,i) X(fs,i)X(fs,i) A(tp, a, b) A(g), a, b)X(fs, i)X(fs, ) 101011 En simplifiant l'expression précédente, on obtient : 101021 IMath.101 ( tp, a, b) = A ( tp, a, b) tp, a, b)X(fs, i)X(fs, i) [0103] La matrice C dépend donc uniquement des paramètres de polarisation a, b et v et en remplaçant la matrice A par son expression, on obtient : [0104] [Math.11] C ( tP, a, b) = a-cos ( q) )2 + b`sin (w)2 ( a2 - b2) cos ( ) sin (w) + iab (a2 b2 COS ( q.)) sin ( tp ) - iab a2 sin (IP ) 2 ± b2 cos( ) 2 - [0105] On calcule lc quotient c suivant : [0106]C(qia, b ) (0,1 ) ) - C ( tp, a, b ) ( 1,0) ) 2ab_ e - ( C q), a, ) a2 + b2 [0107] Avec im(Y) la partie imaginaire de Y, tr(M) la trace de M, C(0,1) le terme de la matrice C situé à la première ligne et à la deuxième colonne et C(1,0) le terme de la matrice C situé à la deuxième ligne et à la première colonne, car : [0108] [Math.12] tr (C ( tp, a, b ) ) a' ( cos ( )2 -I- sin (tp ) 2 ) b2 ( COS ( ql) 2 + sin ( ) 2 = a2 + b- [0109] Si la valeur absolue de e est inférieure à 0.5, la polarisation est de type linéaire donc a est égal à 1 et b est égal à 0. [0110] L'angle d'ouverture v peut alors être obtenu à l'aide de la formule suivante : [0111] [Math.13] 1 C.( qi a b ) (1 0) C(q) a b)(11) tan- ( ° ' ' - tan-1 (te( ' ' - (I) C( , a, b ) (0,0) C( , a, b ) (0,1) - [0112] Avec re(Y) la partie réelle de Y, tan la fonction arc tangente. C(0,0) le ternie de la matrice C situé à la première ligne et à la première colonne et C(1,1) le ternie de la matrice C situé à la deuxième ligne et à la deuxième colonne, car, comme b est égal à zéro : [0113] [Math.14] f c(tp, a, b)(1,0) (C (q, a, b)(1,1) = a-cos Wsin = tan tp ret C(t), a, b)(0,0) I re C ( tp, a, b)(0,1) akos [0114] Si la valeur absolue de c est supérieure à 0.5, la polarisation est de type circulaire. [0115] Si c est négatif, la polarisation est circulaire gauche et a et b ont les valeurs corres- pondantes dans le tableau précédent. [0116] Si c est positif, la polarisation est circulaire droite et a et b ont les valeurs corres- pondantes dans le tableau précédent. 11 [Revendication 1]