FR 2813465 A1 20020301 FR 0011160 A 20000829 La présente invention concerne de manière générale une méthode de détermination de la fonction de transfert spatiale d'un canal de transmission au moyen d'un réseau d'antennes. Plus précisément, lorsque la propagation du signal entre un émetteur et un réseau d'antennes de réception se fait suivant une pluralité de trajets, la présente invention permet d'estimer, pour chaque trajet, la direction d'arrivée conjointement à l'atténuation et la rotation de phase subies par le signal le long de ce trajet. La présente invention s'applique plus particulièrement à l'estimation de canal et de direction d'arrivée en téléphonie mobile. La propagation d'un signal entre un émetteur et un récepteur, notamment entre un téléphone mobile et une station de base peut se faire selon des trajets multiples. Chaque trajet est caractérisé par un retard, une direction d'arrivée, une atténuation et une rotation de phase. il est connu d'utiliser un réseau d'antennes pour faire de la formation de faisceau et/ou de l'annulation d'interférence selon une ou plusieurs directions. Le traitement d'antenne consiste en une pondération des sorties des différentes antennes par des coefficients complexes avant sommation, les coefficients étant choisis de manière à obtenir le diagramme d'antenne équivalente souhaité. On peut ainsi former un faisceau dans la direction d'arrivée du signal utile tout en plaçant des zéros dans les directions d'arrivée des signaux interférents. La plupart des techniques de formation de faisceau requièrent toutefois la connaissance préalable de la direction d'arrivée du signal. La formation de faisceau a été appliquée à la téléphonie mobile, notamment pour diriger un faisceau de réception d'une station de base vers un terminal mobile (liaison montante). La station de base est alors équipée d'une antenne adaptative (dite " antenne intelligente ") capable de pointer dans la direction d'un trajet de propagation issu d'un terminal mobile. Outre la direction d'arrivée du signal utile, il est souvent nécessaire de connaître les caractéristiques du canal de transmission. Ainsi lorsqu'un système de téléphonie mobile fonctionne selon un mode d'accès à répartition par codes (CDMA pour Code Division Multiple Access) et exploite la diversité en réception grâce à un récepteur RAKE, les coefficients complexes utilisés dans les différentes branches du récepteur sont obtenus à partir de l'atténuation et de la rotation de phase subies le long des différents trajets. I est connu d'utiliser une antenne adaptative exploitant la diversité de trajets. Pour ce faire, les systèmes de l'état de la technique utilisent un récepteur RAKE par antenne et effectuent la formation de faisceau à partir des signaux issus des différents récepteurs. On peut également envisager de pointer plusieurs faisceaux selon les directions d'arrivée des différents trajets et d'effectuer ensuite la combinaison des différents trajets dans un récepteur RAKE. Dans un cas comme dans l'autre on procède indépendamment à l'estimation de la direction d'arrivée du signal utile et des caractéristiques du canal de transmission. Cette approche présente le double inconvénient de nécessiter deux étapes d'estimation et de ne pas exploiter la totalité de l'information disponible à la seconde étape d'estimation. Le but de la présente invention est de remédier aux inconvénients précités et de proposer une méthode d'estimation conjointe de la direction d'arrivée et des caractéristiques du canal de transmission. Plus particulièrement, lorsque le canal de transmission comporte plusieurs trajets de propagation, l'invention permet d'estimer la direction d'arrivée de chaque trajet et le coefficient complexe (atténuation et rotation de phase) caractérisant chaque trajet. L'invention est définie par une méthode d'estimation de canal et de direction d'arrivée d'un signal émis par un émetteur et reçu par un réseau d'antennes après s'être propagé selon au moins un trajet, la méthode comprenant, pour chaque trajet, une première étape d'estimation des déphasages (dû) des signaux reçus par les différentes antennes du réseau, une seconde étape d'estimation de l'angle d'arrivée (0) du signal ainsi que de la rotation de phase (v) subie par le signal le long dudit trajet à partir desdits déphasages et une troisième étape d'estimation de l'atténuation (a) subie par le signal le long dudit trajet à partir des valeurs estimées (v, a) de la rotation de phase et de l'angle d'arrivée. Avantageusement, la première étape comprend, pour chaque signal reçu (xe), la minimisation d'une première distance entre une pluralité de valeurs dudit signal (xt(t)) prises à une pluralité d'instants (t) et les valeurs prises aux mêmes instants d'un signal pilote (b(t)) ayant subi le déphasage et l'atténuation à estimer. Selon une caractéristique de l'invention, la seconde étape comprend une levée d'ambiguïté des déphasages, la levée d'ambiguïté opérant de proche en proche d'une antenne à la suivante à partir d'une antenne de référence du réseau. Avantageusement, la levée d'ambiguïté utilise une relation affine entre les déphasages et le rang de l'antenne dans le réseau. Selon une autre caractéristique de l'invention, la rotation de phase (v) et l'angle d'arrivée (0) sont estimés à partir d'une régression linéaire sur les valeurs estimées des déphasages. Avantageusement, la rotation de phase (v) et l'angle d'arrivée (O) sont estimés en minimisant une seconde distance J(v, q() = (v + ( -1).(- -) o =2'r.cos(o).d/; o d est le pas du réseau, A la longueur d'onde du signal et L le nombre d'antennes du réseau. Selon une autre caractéristique de l'invention, l'on procède à une nouvelle estimation (e) des déphasages à partir de la valeur estimée de la rotation de phase (v) et de celle de l'angle d'arrivée (). Avantageusement, l'atténuation (a) est estimée en minimisant une troisième distance entre une pluralité de valeurs dudit signal (xe(t)) prises à une pluralité d'instants (t) et les valeurs prises aux mêmes instants d'un signal pilote (b(t)) déphasé des valeurs de déphasage (e') nouvellement estimées et ayant subi l'atténuation (a) à estimer, la distance étant calculée sur l'ensemble des antemnnes du réseau. L'invention est également définie par un dispositif de réception de signal comprenant un réseau d'antennes et des moyens d'estimation adaptés à mettre en ceuvre les étapes de la méthode d'estimation exposée ci-dessus. Selon une caractéristique de l'invention, le dispositif comprend en sortie de chaque antenne une pluralité de filtres adaptés aux différents trajets de propagation du signal et une pluralité de moyens d'estimation, chaque moyen d'estimation étant associé à un trajet (i) et recevant les sorties des filtres adaptés au trajet correspondant. Selon une autre caractéristique de l'invention, le dispositif comprend des moyens de formation de voie, chaque moyen de formation de voie étant associé à un trajet (i) et recevant du moyen d'estimation associé audit trajet l'estimation Si de l'angle d'arrivée dudit trajet pour former une voie dans la direction correspondante. En outre, chaque moyen de formation de voie peut également recevoir, des moyens d'estimation autres que le moyen d'estimation associé, les estimations (È i., i'ri) des angles d'arrivée des autres trajets pour placer des zéros dans les directions correspondantes. Avantageusement, le dispositif comprend une pluralité de moyens de multiplication complexe, chaque moyen de multiplication complexe étant associé à un trajet (i) et multipliant la sortie d'un moyen de formation de voie par le coefficient complexe &e coefficient d'atténuation fournies par le moyen d'estimation associé audit trajet. Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, apparaîtront plus clairement à la lecture de la description faite en relation avec les figures jointes, parmi lesquelles: La Fig. 1 représente de manière schématique une propagation multi-trajet entre une source et un réseau d'antennes; La Fig. 2 représente un diagramme illustrant une procédure de levée d'ambiguïté utile à un mode de réalisation de l'invention; La Fig. 3 représente un dispositif de réception selon un mode de réalisation de l'invention; Nous nous plaçons dans l'hypothèse illustrée en Fig. 1 o un signal est émis par une source et reçu par un réseau d'antennes. Le réseau illustré est uniforme et linéaire mais, bien entendu, d'autres modes de réalisation peuvent être envisagés, notamment des réseaux sectoriels présentant des antennes élémentaires arrangées selon un arc de cercle. Le signal se propage selon une pluralité p de trajets, chaque trajet i=l..p arrivant sur le réseau selon un angle d'incidence Oi. L'ensemble des trajets constitue le canal de transmission entre la source émettrice et le réseau récepteur. Si plusieurs sources j=l..J émettent chacune un signal Sj, il y aura J canaux de transmission distincts, chaque canal étant constitué d'une pluralité de trajets. Dans le contexte d'un système de télécommunication mobile, les canaux de transmission considérés correspondront à des liaisons montantes entre des terminaux mobiles, dénommés encore ci-après utilisateurs, et la station de base. Lorsque plusieurs liaisons montantes sont simultanément établies avec une même station de base, un canal de transmission donnée peut être perturbé par les signaux des autres utilisateurs (interférence multiutilisateur). En outre, le signal d'un utilisateur donné s'étant propagé selon différents trajets arrive sur le réseau après des temps de propagation différents. Les canaux de transmission étant généralement dispersifs en temps du fait de leur sélectivité en fréquence, leurs réponses impulsionnelles sont de durée non négligeable, ce qui crée des interférences entre les signaux issu du même utilisateur mais ayant suivi des trajets distincts (interférence multi-trajet). Dans certains systèmes, notamment ceux utilisant un accès CDMA par séquences directes (DS-CDMA), on arrive à isoler dans une certaine mesure chaque trajet de chaque utilisateur. Pour ce faire, on exploite les bonnes propriétés d'autocorrélation et d'intercorrélation des séquences d'étalement (spreading codes) en filtrant le signal reçu par un banc de filtres adaptés. Les interférences multi-utilisateur et multi-trajet ne sont toutefois pas complètement éliminées et une interférence résiduelle subsiste. Si les utilisateurs sont suffisamment nombreux, le théorème central limite peut être appliqué et l'interférence résiduelle vue comme un bruit blanc additif gaussien. L'invention présuppose que le récepteur puisse, au moins grossièrement, séparer au niveau de chaque antenne les différents trajets d'un signal issu d'un utilisateur. Ce sera le cas pour un système DS-CDMA, si l'on utilise en sortie de chaque antenne une batterie de filtres adaptés à ces différents trajets. Ce sera encore le cas si le canal de transmission entre l'utilisateur et le récepteur est essentiellement constitué d'un seul trajet. Chaque trajet i, i=l..p, est en réalité la combinaison d'un grand nombre de trajets élémentaires ayant subi des réflexions et des diffusions produites par l'environnement du récepteur et de l'émetteur. A chacun de ces trajets élémentaires sont associées une rotation de phase et une atténuation qui suivent des lois uniformes. D'après le théorème central limite, le signal si(t) provenant du trajet i peut alors s'écrire: s, (t) = a, (t).ejrv ( . bi (t) + ni (t) (1) o a, (t) est l'atténuation associée au trajet i obéissant à une loi de Rayleigh, v (t) est la rotation de phase, uniformément répartie sur] - r; r], b (t) est l'information émise et n, (t) est un bruit blanc additif gaussien incluant l'interférence résiduelle. Le canal est variant dans le temps. Un des paramètres permettant de quantifier cette variation est le temps de cohérence (At)c, qui garantit que deux réponses du canal à deux instants séparés de plus de (At), sont non corrélées. Si les antennes employées dans le réseau ne sont pas omni-directionnelles, le coefficient ca (t) intègre également le gain d'antenne. Nous considérons maintenant le canal de transmission entre un utilisateur et le récepteur. Nous supposerons que le récepteur puisse isoler un trajet i du canal, soit que le canal ne comporte qu'un seul trajet, soit que l'on utilise un filtrage adapté au trajet i. Un seul trajet étant considéré, nous pouvons omettre l'indice i pour simplifier les notations. Nous supposerons également que le signal arrive sur un réseau linéaire uniforme de L antennes sous un angle d'incidence (t). En notant xt le signal reçu au niveau de l'antenne e, ú = (1,...,L) on peut écrire: xt (t) = a(t).e' (t). b(t) + n, (t) (2) avec f( O = t)+( e-1) . ( t)="t)+2 rre-D)cos(6(Ot) (3) o d est le pas du réseau, est la longueur d'onde, b(t) est le symbole émis par l'utilisateur, n. (t) est le bruit blanc additif gaussien au niveau de l'antenne e incluant le bruit thermique, l'interférence multi-trajet et l'interférence multiutilisateur. La réponse du canal est ainsi modélisée par un coefficient d'atténuation a(t) et une rotation de phase v(t). Le déphasage total l'antenne ú et l'antenne de référence 1 sous l'angle d'arrivée 6(t). L'invention propose une estimation conjointe de v(t), (t), a(t) à partir des signaux x(t) et d'une information a priori sur symbole émis b(t). Le symbole b(t) peut être un symbole pilote parfaitement connu du récepteur ou encore un symbole issu d'une opération de décodage antérieure à décisions souples ou dures et présentant un degré de fiabilité suffisant. Nous supposerons v(t), (t), a(t)constants sur une fenêtre d'observation de longueur T choisie faible devant le temps de cohérence (At)c. La longueur T sera choisie d'autant plus courte que la vitesse du mobile est élevée et choisie d'autant plus longue que le rapport signal sur bruit est faible. Elle peut évoluer dans le temps en fonction des conditions de transmission. On cherche 0, a et v tels que la distance quadratique D soit minimale. D= 1(t)-.e.b( (4) La resolution d'un tel problème nécessite la résolution d'un système non linéaire complexe. Le problème est résolu en le scindant en trois étapes. Etape 1: on estime T 2 D = x, (t)- a.e'.b(t)| (5) t=! Etape 2: on estime 9 et v à partir d'une régression linéaire sur les valeurs estimées ô. Etape 3: à partir des valeurs estimées v et È on estime finalement a en L minimisant la distance quadratique D=IDt. t=1 Dans la première étape on cherche à estimer le déphasage du signal au niveau de chaque antenne, c'est-à-dire la valeur ót de ót qui minimise D, : Dt = xt(t) - a.ej.b(t)2 =, (() - a.eje. b(t)). (xt (l) - a.l eJ.b'(t)) t=1 = ( |X!(t)12 + a2,I|b()12 - 2.Re(x1 (i).a.e'jG.b*(I))) t:l (6) Si l'on prend la partie réelle dans l'équation (6), on obtient: x,(t).a.e-'ó. b (t)= a.(xR (t) + j.xt (t)).(bR (t) - j.b' (t)).(cos(ót) - j.sin(,)) = a( [t ().b (t) + x (t).b' (t)]+ j.[xf (t).b (t)-x,(t).b,(t)) (cos(t) /sin(ó)) (7) donc Re(x, (t).a.e-,' ú (t)) = a( [,R (t).bR (t)+ x (t).b' (t).cos(,) + [x! (t).bR(t) - x (t).b' (t)].sin(et)) (8\ En utilisant les équations (6) et (8) on peut calculer la dérivée de la distance D, selon 4!. On obtient: aDj = T a[Re(x, (t).a.e-j4etb* (t)) = -2a.- ( - [x (t).b (t) + xc (t).b' (t)].sin(ó,) + [xc (t).b (t)- X (t).b (t)].cos(,)) *=1 (9) On cherche la valeur 4, de 4t qui annule la dérivée de l'équation (9). 4e doit donc vérifier l'équation suivante: x kl (t). bR (t)-xc ( XR (t)] cos(,)= [x (t).b (t) + x (t).b'(t)].sin() t=l t=l (10) Pour simplifier les écritures définissons les sommes suivantes: Se=ôx! (t),bR (t) *=1 A I T S, fi xc (t).b' (t) t=l t=l En remplaçant les sommes ci-dessus dans l'équation (10), on obtient la valeur estimée du déphasage 4,: Dans la seconde étape, on estime a et v à partir des valeurs estimées sû des déphasages dû. Le calcul des valeurs estimées È et v sera fait dans le cas simple d'un réseau linéaire uniforme mais peut s'étendre à un réseau quelconque. Sur la première antenne (e = 1) prise comme antenne de référence, la rotation de phase ne dépend pas de la direction d'arrivée: 4j = v. On peut donc provisoirement choisir v égal à 41. L'équation (10) a néanmoins deux solutions distinctes sur ] - ,r;iz], qui correspondent à un maximum et un minimum de la fonction D, (4,) et diffèrent d'un angle z. Ces deux valeurs estimées sont: 2 e [- 72>; 2 J et 0 = 41 + or (12) 2 2 Le minimum correspond à la valeur pour laquelle la dérivée seconde 2D avec I = 1 est positive. ao 2_ [(s;3 + S4.( 1 ___ = 2a.[(St + 5).cos(;) + (s; - St).sin(C,)] (13) ó est celle des deux valeurs fournies par l'équation (10) qui correspond à une valeur positive de la formule (13) pour e = 1. Considérons maintenant les autres antennes ( > 1). La valeur de ó est liée à 0 et v par les égalités suivantes: = v + ót = v + (e -1).p = v + 2r d (e - 1) cos0 (14) L'équation (10) a plusieurs solutions. Soit so l'unique solution telle que K [2'1. Les autres solutions sont = ao + k.zr avec k entier relatif 2'2] D'après la relation (14) on doit toujours avoir k _ 2d(ú l+ d o E]-2 X(-1)+v;27cX (e-1)+v] (15) On doit donc rechercher toutes les solutions possibles qui vérifient les conditions (14) et (15), en remplaçant v par sa valeur estimée déjà disponible ót. Si d/ = 1/2, la longueur du segment de l'équation (15) est 2(e - 1). Il existe donc au plus 2( - 1) valeurs possibles de t pour d / tous les t avec k pair sont des minima et tous les t avec k impair des maxima. De même, si possibles de [, seules (ú - 1) valeurs sont des minima. Il reste donc à choisir un minimum parmi les différents minima obtenus. La procédure de levée d'ambigulté est illustrée en Fig. 2. Pour = 2 et d/ 2, on choisit parmi tous les minima possibles celui qui conduit à la plus faible variation de l'estimation de suivante en partant de l'antenne de référence (=l. On dispose ainsi d'un ensemble de L valeurs estimées t pour e =1...,L} à partir desquelles les valeurs estimées Q et ó sont ensuite obtenues par régression linéaire. Plus précisément on cherche les valeurs 0 et 0 qui minimisent l'écart quadratique J: )= (v + (17) Ces deux valeurs sont obtenues par annulation des dérivées partielles de la fonction j: av _ i _L avec,:, (e, e f- i= -1) Lú=1L j (,)2.",(-.(+ e1.-t=.≤> (î-1)2.r ( e1). ó-Oe1 (19) aveco 1).ó =_-1.ó t(e- _) 1)2 En remplaçant l'équation (18) dans l'équation (19), on obtient: (ei)2 V =(e-_).4_-( -_.eO-.e_ -1--1 récepteur RAKE. Nous avons considéré jusqu'à présent un seul trajet du canal de transmission. Si le canal présente plusieurs trajets i=l..p, la méthode d'estimation conjointe doit être appliquée à chaque trajet i. Le canal peut alors être caractérisé par un ensemble de triplets de valeurs estimées (t, j i), i=l. .p, chaque triplet caractérisant lui-même un trajet. La Fig. 3 illustre un dispositif de réception selon un mode de réalisation de l'invention. Le dispositif comprend un réseau d'antennes (3001)... (300L). Chaque antenne (300t) est reliée à une batterie de filtres adaptés (310 i), i=l..p, chaque filtre (310 i) étant adapté à un trajet i. Les sorties des différents filtres (310 i) relatifs à un même trajet i sont reliées à un module d'estimation conjointe (330i) mettant en oeuvre la méthode d'estimation décrite ci-dessus et fournissant un triplet de valeurs estimées ( t, j i, O i) L'estimation est effectuée en parallèle pour tous les trajets. Les signaux de sortie relatifs à un même trajet i=l..P sont également dirigés vers des moyens de formation de voies (3201,..,320p). La valeur estimée ai fournie par le module d'estimation conjointe au formateur de voie (320i) permet à ce dernier de pointer son faisceau de réception dans la direction d'arrivée du trajet i. Dans une variante de réalisation non représentée, chaque formateur de voie (320i) reçoit également les valeurs estimées Èi, fournies par les autres estimateurs (330i,), i'Éi, des trajets i' du même utilisateur. Ceci permet au formateur de voie (320i) de placer des zéros dans le diagramme de réception pour les directions concernées (ôi.. i',i) de manière à mieux séparer les signaux reçus des différents trajets. Les valeurs estimées i et â i sont transmises à un multiplicateur complexe (340i) qui multiplie le signal de sortie du formateur de voie (3200) par le coefficient complexe e-i, conjugué du coefficient complexe du canal. Les sorties des P multiplicateurs sont alors sommées dans un sommateur (350) pour fournir une sortie de type MRC (Maximum Ratio Combining). Ce dispositif de réception permet de réaliser à la fois une formation de voies dans les directions des différents trajets et une combinaison de type MRC comme dans un filtre de Rake. Il est clair pour l'homme du métier que les opérations représentées par les différents modules de la Fig. 3 peuvent être effectuées par un processeur programmé à cet effet ou par une pluralité de processeur dédiés. Bien que l'invention ait été essentiellement décrite dans le cadre de la réception d'un signal émis par un terminal mobile et reçu par l'antenne adaptative d'une station de base, il est clair pour l'homme du métier que l'invention ne se limite pas à cet exemple d'application mais s'étend à tout type d'application requérant l'estimation de direction d'arrivée d'un signal utile et des caractéristiques d'un canal de transmission. REVENDICATIONS 1) Méthode d'estimation de canal et de direction d'arrivée d'un signal émis par un émetteur et reçu par un réseau d'antennes après s'être propagé selon au moins un trajet, caractérisée en ce qu'elle comprend, pour chaque trajet, une première étape d'estimation des déphasages (e,) des signaux reçus par les différentes antennes du réseau, une seconde étape d'estimation de l'angle d'arrivée (6) du signal ainsi que de la rotation de phase (v) subie par le signal le long dudit trajet à partir desdits déphasages et une troisième étape d'estimation de l'atténuation (a) subie par le signal le long dudit trajet à partir des valeurs estimées (v, È) de la rotation de phase et de l'angle d'arrivée. 2) Méthode d'estimation selon la revendication 1, caractérisée en ce que la première étape comprend, pour chaque signal reçu (xe), la minimisation d'une première distance entre une pluralité de valeurs dudit signal (xe(t) ) prises à une pluralité d'instants (t) et les valeurs prises aux mêmes instants d'un signal pilote (b(t)) ayant subi le déphasage et l'atténuation à estimer. 3) Méthode d'estimation selon la revendication 2, caractérisée en ce que le déphasage $, du signal xt reçu par une antenne ú est estimé par: t = Arc tgl ' 4 S- o A T S =Yxi(t).bR(t) t=l A T t::l St=ZXIR(t) b(t)l A T t=l T Se3=-xTR(t).bR(t) t=1 o xR(t) et x'(t) sont respectivement la partie réelle et la partie imaginaire de la valeur du signal reçu x" à l'instant t, bR(t)et b, (t) la partie réelle et la partie imaginaire de la valeur du signal pilote à l'instant t et T une longueur de fenêtre temporelle. 4) Méthode d'estimation selon l'une des revendications précédentes, caractérisée en ce que la seconde étape comprend une levée d'ambiguïté des déphasages, la levée d'ambiguïté opérant de proche en proche d'une antenne à la suivante à partir d'une antenne de référence du réseau. ) Méthode d'estimation selon la revendication 4, caractérisée en ce que la levée d'ambiguïté utilise une relation affine entre les déphasages et le rang de l'antenne dans le réseau. 6) Méthode d'estimation selon l'une des revendications précédentes, caractérisée en ce que la rotation de phase (v) et l'angle d'arrivée (0) sont estimés à partir d'une régression linéaire sur les valeurs estimées des déphasages. 7) Méthode d'estimation selon la revendication 6, caractérisée en ce que la rotation de phase (v) et l'angle d'arrivée () sont estimés en minimisant une seconde distance J(v,o) (v+(e-1).q-t 4 o óest la valeur estimée du déphasage du signal t=l reçu par l'antenne de rang e, ó =21r. cos(>.d/I o d est le pas du réseau, i la longueur d'onde du signal et L le nombre d'antennes du réseau. 8) Méthode d'estimation selon l'une des revendications précédentes, caractérisée en ce que l'on procède à une nouvelle estimation (ôi) des déphasages à partir de la valeur estimée de la rotation de phase (v) et de celle de l'angle d'arrivée (a), 9) Méthode d'estimation selon la revendication 8, caractérisée en ce que l'atténuation (a) est estimée en minimisant une troisième distance entre une pluralité de valeurs dudit signal (xt(t)) prises à une pluralité d'instants (t) et les valeurs prises aux mêmes instants d'un signal pilote (h(t)) déphasé des valeurs de déphasage (s) nouvellement estimées et ayant subi l'atténuation (a) à estimer, la distance étant calculée sur l'ensemble des antennes du réseau. ) Dispositif de réception de signal comprenant un réseau d'antennes (300i) et des moyens d'estimation (330i) adaptés à mettre en oeuvre les étapes de la méthode selon l'une des revendications précédentes. 11) Dispositif de réception selon la revendication 10, caractérisé en ce qu'il comprend en sortie de chaque antenne une pluralité de filtres adaptés (310i) aux différents trajets de propagation du signal et une pluralité de moyens d'estimation (330i), chaque moyen d'estimation étant associé à un trajet (i) et recevant les sorties des filtres adaptés (310i) au trajet correspondant. 12) Dispositif de réception selon la revendication 11, caractérisé en ce qu'il comprend des moyens de formation de voie (320i), chaque moyen de formation de voie (320i) étant associé à un trajet (i) et recevant du moyen d'estimation (330i) associé audit trajet l'estimation Ai de l'angle d'arrivée dudit trajet pour former une voie dans la direction correspondante. 13) Dispositif de réception selon la revendication 12, caractérisé en ce que chaque moyen de formation de voie (320i) reçoit également des moyens d'estimation (330i') autres que le moyen d'estimation associé (330i) les estimations (0s, i',i) des angles d'arrivée des autres trajets pour placer des zéros dans les directions correspondantes. 14) Dispositif de réception selon la revendication 12 ou 13, caractérisé en ce qu'il comprend en outre une pluralité de moyens de multiplication complexe (340i), chaque moyen de multiplication complexe étant associé a un trajet (i) et multipliant la sortie d'un moyen de formation de voie (320i) par le coefficient complexe Èi.e-e' ou V iet x i sont les valeurs estimées de la rotation de phase et du coefficient d'atténuation fournies par le moyen d'estimation associé audit trajet.