FR 2954515 A1 20110624 FR 0906269 A 20091223 Procédé et dispositif de goniométrie de signaux interférant un signal utile modulé par plusieurs porteuses La présente invention concerne un procédé et un dispositif de goniométrie de signaux interférant un signal utile modulé par plusieurs porteuses. Elle s'applique en particulier à la goniométrie de signaux parasites présents sur la même bande de fréquence qu'un signal de type OFDM, c'est à dire un signal modulé par répartition de fréquences orthogonales. Lorsque plusieurs émissions sont reçues sur un même canal de fréquence, seule l'émission de plus forte amplitude peut généralement faire l'objet d'une goniométrie. Or, de nombreuses applications requièrent de pouvoir déterminer la direction d'arrivée de toutes les émissions reçues sur un même canal de fréquence, y compris les émissions reçues avec un faible niveau de puissance, de l'ordre par exemple de 20 ou 25 dB au dessous de l'émission de puissance la plus élevée. A titre d'illustration, on souhaite déterminer l'origine de signaux perturbant un programme de télévision émis sur une forme d'onde de type DVB (« Digital Video Broadcasting »). Une première catégorie de méthodes connues permettant d'effectuer une goniométrie en présence de plusieurs sources de signaux regroupe les méthodes à haute-résolution et les méthodes à super résolution. Parmi ces méthodes, on peut citer MUSIC, ESPRIT, la méthode du maximum de vraisemblance, ou encore la méthode CAPON. Ces méthodes ne font aucune hypothèse sur la forme d'onde utilisée par les émissions reçues et permettent de séparer N-1 sources de signaux (dans la pratique N/2 sources) avec N antennes associées à N voies de réception. Or, il est parfois souhaitable d'effectuer la goniométrie de multiples sources en utilisant seulement deux voies de réception. Par ailleurs, les méthodes à haute résolution ne sont pas toujours applicables du fait de leur besoin en ressources de calcul. Une seconde catégorie de méthodes connues concerne la goniométrie sur des signaux interférant une émission à forme d'onde série comportant des séquences connues. Ce type de méthode exploite la connaissance de la forme d'onde des signaux reçus (ou d'une partie des signaux reçus), de manière à améliorer les performances de la goniométrie. Un exemple d'une telle méthode est décrit dans un brevet européen de Thales référencé sous le numéro de publication EP0991952. Cette méthode s'applique cependant plus particulièrement à des formes d'onde à modulation linéaire (par exemple les formes d'onde de type GSM ou UMTS). Or, il est souhaité de pouvoir effectuer une goniométrie sur des signaux interférant des signaux modulés par multiporteuses, comme dans les émissions de type DVB ou WiMax. Dans la suite du texte, on désigne par « symbole OFDM » le signal utile transmis correspondant à la sortie de l'opération de transformée de Fourier après numérisation du signal reçu, et par « symbole », les symboles élémentaires, n-PSK (modulation par déplacement de phase) ou QAM (modulation d'amplitude en quadrature) par exemple, modulant les sous-porteuses d'un symbole OFDM. Les symboles élémentaires 391 et symboles OFDM 381, 382 sont illustrés en figure 3b. Parmi les symboles, on distingue les symboles connus, permettant la synchronisation, l'apprentissage du canal de propagation et la démodulation, et les symboles de données transportant l'information utile. Un but de l'invention est de déterminer la direction d'arrivée de signaux interférant un signal modulé sur plusieurs fréquences porteuses et comportant des symboles connus, quelque soit les amplitudes relatives des signaux interférents et dudit signal modulé. A cet effet, l'invention a pour objet un procédé de goniométrie pour déterminer la direction d'arrivée de signaux interférant un signal utile à répartition de fréquences orthogonales comportant des symboles élémentaires connus, les signaux comprenant le signal utile et les signaux interférents étant reçus par un réseau d'antennes alimentant au moins deux voies de réception, le procédé comprenant une étape de numérisation des signaux reçus sur chaque voie de réception, une étape de synchronisation sur les symboles dudit signal utile, une étape de démodulation par transformée de Fourier des signaux numérisés, le procédé étant caractérisé en ce qu'il comprend pour un ou plusieurs symboles issus de l'opération de transformée de Fourier après numérisation du signal reçu par exemple pour une forme d'onde de type OFDM pour un ou plusieurs symboles OFDM, les étapes suivantes : o estimer sur plusieurs sous-porteuses la contribution des signaux interférents au signal reçu à partir desdits symboles 35 élémentaires connus ; o effectuer une goniométrie sur les signaux interférents à partir des estimations effectuées lors de l'étape précédente. Selon une mise en oeuvre du procédé selon l'invention, l'étape d'estimation de la contribution Yk des signaux interférents comprend, au moins pour chaque sous-porteuse du signal reçu modulée avec des symboles connus, pour chaque voie de réception k, une première sous-étape d'estimation de la contribution hk du canal de propagation utile reçu, effectuée à partir du signal numérisé puis démodulé par transformée de Fourier Zk et des symboles connus utilisés pour moduler le signal utile, puis une deuxième sous-étape, effectuée sous-porteuse par sous-porteuse, de retranchement au signal numérisé puis démodulé par transformée de Fourier Zk d'une estimation de la contribution du signal utile, calculée à partir de ladite estimation de hk et des symboles connus ou démodulés. Selon une mise en oeuvre du procédé selon l'invention, le signal utile à répartition de fréquences orthogonales utilise la technologie MIMO (« Multiple Input Multiple Output ») ; ce signal est émis à partir de plusieurs voies d'émission distinctes, chacune desdites voies d'émission utilisant ses propres symboles connus pour moduler le signal utile émis, l'étape d'estimation de la contribution des signaux interférents comprenant pour chaque voie d'émission 1, une première sous-étape d'estimation de la contribution hk,, du signal utile pour au moins chaque sous-porteuse f du signal reçu modulée avec des symboles connus, pour chaque voie de réception k, à partir du signal numérisé puis démodulé par transformée de Fourier Zk et des symboles connus utilisés pour moduler le signal utile sur ladite voie d'émission 1, une deuxième sous-étape combinant le signal numérisé puis démodulé par transformée de Fourier Zk, les estimations hk,, avec les symboles connus de chacune des voies d'émission correspondantes pour produire une estimation de la contribution des signaux interférents. Le procédé de goniométrie peut comprendre une étape d'estimation des symboles transmis sur les sous-porteuses de données à partir du vecteur signal z reçu sur l'ensemble des K voies de réception et de h(, le vecteur canal de propagation reçu sur l'ensemble des K voies de réception. Avantageusement, le procédé selon l'invention peut être exécuté pour des signaux interférents dont la puissance est inférieure à la puissance reçue du signal utile, l'écart de puissance pouvant être supérieur à 20 ou 25 dB. Le procédé selon l'invention peut notamment être appliqué à des goniométries de signaux interférant des émissions de type DVB ou WiMax. Un avantage du procédé selon l'invention est qu'il permet d'améliorer la précision et la sensibilité d'une goniométrie effectuée sur les signaux de forme d'onde parallèle, par exemple sur les signaux de type 1 o OFDM. D'autres caractéristiques apparaîtront à la lecture de la description détaillée donnée à titre d'exemple et non limitative qui suit faite en regard de dessins annexés qui représentent : 15 • la figure 1, un premier exemple de mise en oeuvre du procédé selon l'invention ; • la figure 2, un deuxième exemple de mise en oeuvre du procédé selon l'invention ; • la figure 3a, un diagramme temps-fréquence illustrant la répartition 20 des sous-porteuses pilotes dans une trame de type DVB-T ; • la figure 3b, une représentation d'un symbole OFDM et de symboles élémentaires, • la figure 4, un exemple de système sur lequel le procédé selon l'invention peut être mis en oeuvre si un seul signal interférent est 25 présent ; • la figure 5, un schéma d'un système sur lequel peut être exécuté le procédé selon l'invention avec une méthode à haute résolution ; • la figure 6, un synoptique illustrant les étapes d'un procédé selon l'invention appliqué à un système de type MIMO. 30 La figure 1 présente un premier exemple de procédé selon l'invention. Des signaux 121 sous forme d'ondes électromagnétiques sont reçus et transmis à un récepteur radioélectrique comportant K voies de réception. Ces signaux reçus comportent à la fois un signal utile modulé 35 selon une forme d'onde de type OFDM et également un signal interférant ce signal utile sur au moins une partie de sa bande de fréquence. Le signal utile comprend des symboles connus émis sur certaines fréquences porteuses, qualifiées par la suite de « porteuses pilotes ». Les autres fréquences porteuses du signal utile OFDM, appelées par la suite « porteuses de données » sont modulées avec des données. La première étape 101 du procédé selon l'invention transpose le signal reçu 121 en bande de base et numérise le signal reçu 121 sur chaque voie de réception selon des techniques connues de l'Homme du métier. La deuxième étape 102 effectue une synchronisation du récepteur sur le signal numérisé Xk obtenu en sortie de la première étape 101. La synchronisation peut notamment s'effectuer en utilisant le procédé décrit dans le brevet européen de la société Thales et référencé sous le numéro de publication EP0665665 qui permet d'effectuer une prise de synchronisation en présence d'interférence. La troisième étape 103 est une étape de démodulation par transformée de Fourier du signal numérisé. Sur chaque symbole OFDM d'indice n et pour chaque voie de réception k, une transformée de Fourier (FFT) est appliquée sur le signal numérisé Xk comme suit : zk(n, t) = FFT(Xk(((n-1)*(NFFT+Ng) +Ng :n*(NFFT+Ng))), où le terme zk(n, f) représente le symbole démodulé, f étant l'indice de la sous-porteuse, NFFT la taille de la FFT qui est un paramètre de la forme d'onde du signal utile, Ng étant le temps de garde. La quatrième étape 104 est une étape d'estimation de la contribution hk(f) du canal de propagation utile, sur chaque porteuse f correspondant à un symbole connu, et pour chaque voie de réception k. Par exemple, pour une voie de réception k et la sous-porteuse d'indice f, cette estimation hk(t) est effectuée via le critère suivant : hk(f) = rdZk(f) / rdd(f) où r (f)_Ed(n,f).d*(n,f) et rdzk(f)_Ezk(n,f).d*(n,f) n n d(n, f) étant un symbole connu émis sur la sous-porteuse pilote d'indice f du symbole OFDM d'indice n. La cinquième étape 105 estime la contribution des signaux interférents en retranchant aux signaux reçus la partie estimée du signal utile sur chaque sous-porteuse f. Ainsi, pour chaque symbole OFDM d'indice n, pour chaque voie de réception k et pour chaque porteuse correspondant à un symbole connu, est effectuée l'estimation suivante : Yk(n,f) = zk(n,f) û hk(f)*d(n,f) où le terme yk(n,f) représente la contribution des signaux interférents sur le symbole d'indice n reçu sur la voie de réception k sur la porteuse f, le terme zk(n,f) représentant le symbole d'indice n du signal reçu démodulé par transformée de fourier, par exemple par FFT (transformée de Fourier rapide). Comme illustré en figure 1, le procédé selon l'invention permet d'isoler la contribution des signaux interférents pour effectuer leur goniométrie 106a comme si le signal utile OFDM était absent. Outre une goniométrie des signaux interférents, une goniométrie 106b du signal utile OFDM peut être effectuée avec une précision et une sensibilité accrue. La figure 2 illustre un deuxième exemple de mise en oeuvre du procédé selon l'invention. Dans le cas où les sous-porteuses pilotes sont réparties parmi les sous-porteuses de données, on peut estimer la contribution du canal utile hk(f) sur les sous-porteuses de données en effectuant une interpolation à partir de l'estimation du canal utile effectuée sur les sous-porteuses pilotes. Par exemple, dans le cas d'une forme d'onde de type DVB, une sous- porteuse sur trois est une sous-porteuse pilote, ce qui permet d'effectuer une interpolation et d'estimer la contribution du canal utile sur les sous-porteuses de données intercalées entre les sous-porteuses pilotes. Dans d'autres formes d'ondes, par exemple les formes d'onde de type WiMAX, un symbole placé en début de trame ne comporte que des sous-porteuses pilotes, les autres symboles comportant quasiment uniquement des sous-porteuses de données (et parfois quelques sous-porteuses pilotes pour assurer le suivi du canal). La figure 3a illustre, à travers un diagramme temps-fréquence, la répartition des sous-porteuses pilotes dans une trame de type DVB-T. Les sous-porteuses pilotes sont représentées par des disques pleins 301, tandis que les sous-porteuses de données sont représentées par des cercles 302. Dans un signal de forme d'onde DVB-T, les sous-porteuses pilotes sont insérées comme suit : • des « continuai pilots » ou sous-porteuses pilotes fixes 311 sont présentes au même emplacement (c'est-à-dire à la même fréquence) sur chaque symbole, • des « scattered pilot » ou sous-porteuses pilotes 312 réparties sont 5 régulièrement espacées toutes les 12 sous-porteuses et leur fréquence varie d'un symbole à l'autre (décalages de 3 sous- porteuses pour une périodicité de 4 symboles). Les autres sous-porteuses sont des porteuses de données. Comme illustré en figure 2, une étape optionnelle 110 du procédé 10 consiste à estimer les symboles transmis sur les sous-porteuses de données. Soit d(n,f) le symbole transmis sur la sous-porteuse f du symbole OFDM n. Le symbole d(n,f) peut être estimé, par exemple, à partir de z(n,f), le vecteur signal reçu sur l'ensemble des K voies de réception et de h(f), le vecteur canal de propagation reçu sur l'ensemble des K voies de réception, 15 de la façon suivante : d[n,f]= h[n,f]"l(h[n,frh[n,f])z[n,f] avec : z[n,f] = [zl[n,f],z2[n,f],...,zK[n,f]]T h[f] = [hi[f],h2[f],...,hK[f]]T Le symbole d(n,f) peut également être estimé en utilisant une méthode 20 décrite dans le brevet français référencé sous le numéro de publication FR2854290 qui permet la démodulation d'une forme d'onde OFDM en présence d'interférences. L'intérêt d'effectuer l'étape optionnelle 110 est qu'elle permet de pouvoir utiliser également les symboles de données pour estimer la 25 contribution des signaux interférents : en sortie de l'étape 110, l'étape 105 d'estimation de la contribution des signaux interférents peut être effectuée sur toutes les sous-porteuses f par la formule : yk(n,f) = zk(n,f) ù hk(f)*d(n,f) où d(n,f) représente, suivant la porteuse f et le symbole OFDM n , soit un 30 symbole de données, démodulé lors de l'étape 110, soit un symbole pilote. Si un seul signal interférent est présent, le procédé selon l'invention permet de déterminer sa direction d'arrivée en utilisant, par exemple, une méthode de goniométrie telle que l'interférométrie, la corrélation vectorielle, ou la formation de voies. La figure 4 présente un premier exemple de système sur lequel le procédé selon l'invention peut être mis en oeuvre si un seul signal interférent est présent. Le système 400 comprend, dans l'exemple, cinq antennes 401, 402, 403, 404, 405 de réception réparties aux coins d'un polygone régulier. Toutes les antennes 401, 402, 403, 404, 405 sont raccordées à un module de commutation 410 permettant de sélectionner les signaux reçus par deux antennes choisies parmi les cinq 401, 402, 403, 404, 405 et d'orienter les signaux sélectionnés vers un module traitement 420 permettant, par exemple, d'effectuer des calculs d'interférométrie. Le module de traitement 420 comprend une chaîne formée respectivement d'un récepteur 421 à deux voies, d'un convertisseur analogique numérique 422, et d'un calculateur 422. Pour le système 400, le procédé selon l'invention peut sélectionner, par exemple, les couples d'antennes suivants : (première antenne 401 ; troisième antenne 403) (deuxième antenne 402 ; quatrième antenne 404) (troisième antenne 403 ; cinquième antenne 405) (quatrième antenne 404 ; première antenne 401) et (cinquième antenne 405 ; deuxième antenne 402). L'ensemble des couples d'antennes est noté C par la suite. Un exemple d'exécution du procédé selon l'invention pour le système 400 est maintenant décrit pour une forme d'onde de type DVB-T, c'est-à-dire DVB Terrestre, avec les paramètres suivants : mode 8 K, une bande de 8 MHz et un temps de garde de 1/8, ce qui correspond à 6817 sous-porteuses, Nfft=8192, Ng=1024, Fe=64/7 MHz. Selon la première étape 101 du procédé, les signaux reçus xk et xm sur les deux antennes k et m sont numérisés et translatés en bande de 25 base. La durée d'acquisition correspond à 8 symboles OFDM. Ensuite, la deuxième étape 102 du procédé est exécutée. Un critère de synchronisation est déterminé sur une durée correspondant à quatre symboles OFDM en utilisant la méthode décrite dans le brevet européen référencé sous le numéro de publication EP0665665. La séquence 30 de synchronisation est formée à partir des sous-porteuses pilotes (signal temporel provenant de la transformée de Fourier inverse des porteuses pilotes dispersées, les autres porteuses étant mises à zéro, en ajoutant l'intervalle de garde Ng) et ce sur quatre symboles OFDM consécutifs. Le critère de synchronisation peut s'écrire de la façon suivante : c(n) = d r (n)RXx (n)rXd (n) Illl avec : N-1 N-1 N-1 d 2 = Id(k)12 rxd (n) = E x(n + k)d * (k) RXX (n) = E x(n + k)x" (n + k) k=0 k=0 k=0 N = 4 * (Nfft+N9') x(n) = [xk(n) ; Xm(n)]T On note maintenant n0 la position de synchronisation choisie de façon à prendre en compte l'ensemble des trajets de propagation dans l'intervalle de garde. Selon la troisième étape 103 du procédé, on effectue la 10 transformée de Fourier sur les signaux reçus sur les voies k et m sur quatre symboles suivant la position de synchronisation n0 : Zk(n, f) = FFT(xk(n0+(n-1)*(NFFT+Ng) +Ng : nO+n*(NFFT+Ng))) zm(n, f) = FFT(xm(n0+(n-1)*(NFFT+Ng) +Ng : n0+n*(NFFT+Ng))) Selon la quatrième étape 104 du procédé, la contribution du canal 15 de propagation hk(f) et hm(f) sur chaque porteuse f est déterminée pour la voie de réception k et pour la voie de réception m. Ainsi, sur chacune des porteuses pilotes dispersées (1 porteuse sur 12), on effectue l'estimation de canal par la formule hk(f) = rdzk(f) / rdd(f) 20 où rdd(f)=Ed(p+4n,f) d*(p+4n,f) n=0 et rd (f) = Ezk(p+4n,f) d*(p+4n,f) n=0 avec p = 1 à 4 pour tenir compte de la périodicité de la position des sous-porteuses pilotes dispersées sur quatre symboles. Et f = 1 +3*(p-1)+12*I 25 A l'issue de cette opération, les coefficients hk(f) sont donc estimés toutes les trois sous-porteuses. Ensuite, une interpolation est mise en oeuvre pour obtenir l'estimation des coefficients hk(f) sur toutes les sous-porteuses. L'étape 110 d'estimation des symboles transmis sur les sous-porteuses de données est exécutée, comme expliquée plus haut. Sur chaque 30 sous-porteuse de données, la décision des symboles est réalisée en tenant compte de la constellation (QPSK, 16-QAM ou 64-QAM). Selon la cinquième étape 105, on estime la contribution des signaux interférents sur les voies de réception k et m, selon la méthode expliquée plus haut, en utilisant l'ensemble des sous-porteuses (sous-porteuses de données et sous-porteuses pilotes). Ensuite, aux fins d'effectuer une goniométrie par interférométrie, on effectue les étapes suivantes : pour chaque couple d'antennes (k ; m), sur les signaux yk(n,f) et ym(n,f) issus du procédé selon l'invention : • dans un premier temps, on estime, pour chaque couple d'antennes (k ;m), des différences de phase entre les signaux reçus par ces antennes k, m, différences de phase notées 4mesurékm, entre les signaux xk(n) et xm(n) reçus sur les deux antennes k et m. Cette différence de phase s'estime en utilisant les signaux estimés yk(n,f) et ym(n,f) issus du procédé selon l'invention, par exemple, via la formule suivante : 8 6817 Dresuré(k,m) = phase ( Il yk (n, f) ym (n, f) ) n=1 f =l • dans un deuxième temps, on recherche la direction d'arrivée (azimut et/ou site) correspondant au jeu des cinq différences de phase mesurées : pour chaque direction d'arrivée (DOA), les différences de phase théoriques entre les antennes k et m sont notées .àehé°k m(DOA) et sont obtenues par des méthodes connues de l'homme du métier faisant intervenir si nécessaire une calibration du réseau d'antennes. Ensuite, un algorithme recherche la direction d'arrivée qui minimise l'écart entre les différences de phase théoriques et les différences de phases mesurées sur chaque couple c d'antenne (k,m) de l'ensemble C. Par exemple, DOA = ArgMin { E 0`ché° (DOA) - Qmesuré 12 } Le procédé selon l'invention rend donc possible, dans un système comportant un réseau de K antennes, K >2, associé à deux voies de réception, la détermination de la direction d'arrivée d'un signal interférent en présence d'un signal utile OFDM. La goniométrie du signal utile peut également être effectuée par l'algorithme d'interférométrie décrit ci-dessus pour effectuer la goniométrie des signaux interférents en exploitant les différences de phase de hk(f) et hm(f), les estimations du canal de propagation utile obtenues sur les voies k et m 6817 4mesu e-k m = phase [ E hk (f) (f) ] f=1 En outre, si le signal utile est reçu suivant plusieurs trajets de propagation, il est aussi possible d'effectuer la goniométrie de chaque trajet utile : on effectue dans un premier temps une FFT inverse pour estimer la réponse impulsionnelle temporelle du canal de propagation comme suit : gk(n) = FFT-1(hk(f)) La goniométrie est effectuée dans un second temps pour chaque trajet 10 détecté en utilisant les coefficients de gk(n) pour calculer les différences de phase : Orsurék m (n) = phase [gk (indice(n)) g;,, (indice(n) ] avec n = 1, .., nombre de trajets détectés indice : tableau contenant les retards des trajets détectés 15 La détection des trajets peut s'effectuer par exemple à partir du critère de synchronisation ou d'après la valeur de lgk(n)12, qui correspond à la puissance du trajet n. Si plusieurs signaux interférents sont présents, l'exécution d'une méthode à haute résolution, placée en sortie de la cinquième étape 105 20 d'estimation des signaux interférents, permet d'effectuer leur goniométrie en présence d'un signal utile OFDM avec de meilleures performances que si le procédé selon l'invention n'était pas mis en oeuvre. Un système sur lequel peut être exécuté le procédé selon l'invention avec une méthode à haute résolution est schématiquement 25 représenté en figure 5. Le système 500 comporte un réseau antennaire pentagonal, chaque antenne 501, 502, 503, 504, 505 du réseau alimentant un module de traitement 510 configuré pour exécuter le procédé de goniométrie selon l'invention avec une méthode à haute résolution. Grâce au procédé selon l'invention, il est possible de séparer plus de sources que de 30 signaux à nombre fixé d'antennes. Le module de traitement 510 comprend une chaîne formée respectivement d'un récepteur 511 à cinq voies, d'un convertisseur analogique numérique 512, et d'un calculateur 512. Par exemple, l'algorithme MUSIC est mis en oeuvre à partir de l'estimation de la matrice de corrélation des signaux interférents (respectivement des signaux utiles) pour estimer les directions d'arrivée des signaux interférents (respectivement des signaux utiles). Le réseau d'antennes comporte cinq capteurs 501, 502, 503, 504, 505 et cinq voies de réception numérisées en parallèle. Pour le signal DVB-T de l'exemple précédent, les étapes 101 à 105 sont effectuées comme dans l'exemple précédent pour les cinq voies de réception en parallèle. La matrice de corrélation des signaux interférents est calculée à partir des signaux y(f). 8 6817 RYY= EEy[n,f]'y[n,f') n=1 f=1 Avec y[n,t] = [y1[n,t], Y2[n,t], ..., YK[n,t]]T La matrice de corrélation du signal utile est calculée à partir du vecteur h des estimées de canal : 8 6817 Ruu = Eh[n,f]•h[n,f]H) n=1 f=1 Avec h[n,t] = [hl[n,t],h2[n,t],...,hK[n,t]]T La figure 6 illustre le procédé selon l'invention appliqué à un système de type MIMO. Le signal utile est émis sur plusieurs antennes émettrices. Dans le cas d'un système de type MIMO, acronyme anglo-saxon pour « Multiple Input, Multiple Output », la quatrième étape 104 et la cinquième étape 105 du procédé sont modifiées de la façon suivante pour tenir compte des signaux émis par les différentes antennes d'émission (les autres étapes restent inchangées). ^ La quatrième étape 104 est répétée 204, 204', 204" pour chacune des M voies d'émission utilisées pour émettre le signal utile, de sorte que l'on estime pour chaque voie de réception k, les M canaux de propagation utile hk,m(f) correspondant à chacune des voies d'émission. ^ La cinquième étape 105 estime les contributions des signaux interférents obtenus pour chaque voie de réception et pour au 30 moins toutes les sous-porteuses ayant été modulées avec des symboles connus par la formule suivante yk(n,f) = zk(n,f) û somme_sur_m hk,m(f)*dm(n,f) où dm(n,f) correspond au symbole élémentaire utile émis depuis l'antenne m sur la porteuse f et le symbole OFDM n. En général, pour une porteuse f correspondant à un symbole élémentaire connu émis par l'une des antennes d'émission, les autres antennes d'émission n'émettent aucune donnée (porteuse mise à zéro), de sorte que la somme écrite ci-dessus ne comporte qu'un seul terme si on prend en compte les sous-porteuses f correspondant à un symbole élémentaire connu émis par l'une des antennes d'émission. L'étape 110 de démodulation des symboles élémentaires peut également être effectuée. Dans ce cas, le terme yk(n,f) peut être estimée sur ~o l'ensemble des sous-porteuses (symboles élémentaires de données et connus). REVENDICATIONS1. Procédé de goniométrie pour déterminer la direction d'arrivée de signaux interférant un signal utile à répartition de fréquences orthogonales comportant des symboles élémentaires connus, les signaux comprenant le signal utile et les signaux interférents étant reçus par un réseau d'antennes alimentant au moins deux voies de réception, le procédé comprenant une étape de numérisation (101) des signaux reçus sur chaque voie de réception, une étape de synchronisation (102) sur les symboles dudit signal utile, une étape de démodulation (103) par transformée de Fourier des signaux numérisés, le procédé étant caractérisé en ce qu'il comprend, pour un ou plusieurs symboles issus de ladite transformée de Fourier, les étapes suivantes : o estimer (104, 204, 105, 205) la contribution des signaux interférents au signal reçu à partir desdits symboles élémentaires connus ; o effectuer une goniométrie (106) sur les signaux interférents à partir des estimations effectuées lors de l'étape précédente. 2. Procédé de goniométrie selon la revendication 1, caractérisé en ce que l'étape (104, 105) d'estimation de la contribution Yk des signaux interférents comprend, au moins pour chaque sous-porteuse du signal reçu modulée avec des symboles connus, pour chaque voie de réception k, une première sous-étape (104) d'estimation de la contribution hk du canal de propagation utile reçu effectuée à partir du signal zk numérisé puis démodulé par transformée de Fourier et des symboles connus utilisés pour moduler le signal utile, puis une deuxième sous-étape (105), effectuée sous-porteuse par sous-porteuse, de retranchement au signal numérisé puis démodulé par transformée de Fourier Zk d'une estimation de la contribution du signal utile, calculée à partir de ladite estimation de hk et des symboles connus ou démodulés . 3. Procédé de goniométrie selon la revendication 1 ou 2, le signal utile OFDM étant émis à partir de plusieurs voies d'émission distinctes, chacune desdites voies d'émission utilisant ses propres symbolesconnus pour moduler le signal utile émis, caractérisé en ce que l'étape (204, 204', 204", 205) d'estimation de la contribution des signaux interférents comprend pour chaque voie d'émission I, une première sous-étape (204, 204', 204") d'estimation de la contribution hk,, du signal utile pour au moins chaque sous-porteuse f du signal reçu modulée avec des symboles connus, pour chaque voie de réception k, à partir du signal Zk numérisé puis démodulé par transformée de Fourier et des symboles connus utilisés pour moduler le signal utile sur ladite voie d'émission I, une deuxième sous-étape (205) combinant le signal Zk numérisé puis ~o démodulé par transformée de Fourier, les estimations hk,, avec les symboles connus de chacune des voies d'émission correspondantes pour produire une estimation de la contribution des signaux interférents. 4. Procédé de goniométrie selon l'une des revendications 2 et 3, 15 caractérisé en ce que l'estimation (104) de la contribution du signal utile est effectuée en déterminant le rapport suivant : hk(f) = rdZk(f) I rdd(f), où rdd (.f) = E d(n, f) . d * (n, .f) et n rdZk (f) = E zk (n, f) . d * (n, f) , d(n, f) étant un symbole connu émis sur la n sous-porteuse d'indice f du symbole OFDM d'indice n. 20 5. Procédé de goniométrie selon l'une des revendications 2 à 4, caractérisé en ce qu'il comprend une étape (110) d'estimation des symboles transmis sur les sous-porteuses de données à partir du vecteur signal z reçu sur l'ensemble des K voies de réception et de h(f), le vecteur canal 25 de propagation reçu sur l'ensemble des K voies de réception. 6. Procédé de goniométrie selon l'une des revendications précédentes, caractérisé en ce que la puissance reçu des signaux interférents est inférieure à la puissance reçue du signal utile, l'écart de puissance étant 30 compris entre 20 et 25 dB.