FR 3079308 A1 20190927 FR 1800238 A 20180322 I Dispositif multi-antennes pour la réjection de multi-trajets dans un système de navigation par satellite et procédé associé Le domaine de l'invention concerne les systèmes de navigation 5 par satellite et plus particulièrement les dispositifs de radiolocalisation par calcul de temps de propagation des signaux émis par les satellites. On entend ici par « système de navigation par satellite » tout système dédié à la navigation large zone, comme par exemple les systèmes 10 existants GNSS (« Global Navigation Satellite System ») appelés GPS, GLONASS ou GALILEO, ainsi que tous leurs équivalents et dérivés. L'homme du métier connaît bien le principe de localisation des systèmes de navigation par satellite. Le signal radiofréquence émis par un satellite est codé et on utilise le temps mis par ce signal pour atteindre le récepteur à 15 localiser pour déterminer la distance entre ce satellite et ce récepteur, appelée pseudo-distance. La précision des systèmes de navigation par satellite est affectée par un certain nombre d'erreurs. Ces erreurs peuvent se décliner en deux catégories : les contributions globales et les contributions locales. On peut citer pour les contributions globales, les erreurs liées au 20 passage des ondes électromagnétiques dans l'ionosphère et les erreurs liées aux satellites (erreurs d'orbite et erreurs d'horloge). En ce qui concerne les contributions locales, on peut citer les erreurs liées au passage des ondes électromagnétiques dans la troposphère, les erreurs de réflexion des signaux, les erreurs liées aux interférences, les erreurs dues aux zones 25 blanches et le bruit des récepteurs. En environnement urbain comme illustré par la figure 1, le dispositif de radiolocalisation est particulièrement affecté par le phénomène de réflexion des signaux sur les obstacles urbains comme par exemple les 30 façades d'immeuble. La figure 1 représente un véhicule automobile circulant en environnement urbain et un satellite émettant des signaux de 2 radiolocalisation. Les multi-trajets MP sur les signaux utiles LOS constituent les défauts prépondérants pour la localisation. En effet, les signaux mutitrajets MP introduisent un biais dans l'estimation du temps de propagation des signaux ce qui provoque des erreurs de localisation du récepteur. Il est important de pouvoir supprimer ou estimer ces multi-trajets pour améliorer la précision de localisation des dispositifs de localisation. La demande de brevet internationale WO 2012/025306 décrit un dispositif de réception de signaux de radionavigation par satellite apte à 10 estimer et supprimer les multi-trajets affectant le signal reçu. Un tel dispositif comprend plusieurs antennes ou capteurs et plusieurs corrélateurs. Un inconvénient de cette solution est qu'elle met en oeuvre un procédé d'estimation des paramètres des signaux multi-trajets qui suppose que tous les capteurs du dispositif sont identiques et en particulier, qu'ils 15 présentent tous des diagrammes de rayonnement identiques. Cependant, dans la réalité, ce constat n'est pas vérifié et, au contraire, les diagrammes de rayonnement sont différents, notamment en raison de phénomènes de couplage. L'absence de prise en compte de ces différences peut conduire à rendre inopérantes la solution décrite dans la demande 20 précitée. L'invention propose un perfectionnement du procédé et du dispositif décrits dans la demande WO 2012/025306 pour prendre en compte les différences de diagrammes de rayonnement des différentes antennes du 25 dispositif de réception. Ainsi, l'invention a pour objet un procédé d'estimation des paramètres de signal utile et de signaux multi-trajet originaires d'un signal de radiolocalisation émis par un satellite, au moyen d'un dispositif de localisation 30 comprenant au moins deux capteurs aptes à recevoir ledit signal, le procédé comprenant les étapes de : 3 - Corréler le signal reçu par lesdits capteurs avec un code local au moyen de corrélateurs, Construire, pour chaque capteur, une fonction d'intercorrélation échantillonnée du signal reçu avec le code local, 5 Déterminer une fonction d'intercorrélation spatio-temporelle à partir de la concaténation des fonctions d'intercorrélation obtenues à l'étape précédente pour chaque capteur, Estimer des paramètres représentatifs du signal utile et des signaux multi-trajet par application d'un algorithme de maximum de vraisemblance, lesdits paramètres représentatifs incluant au moins une amplitude complexe estimée indépendamment pour chaque capteur. Selon un aspect particulier de l'invention, lesdits paramètres représentatifs incluent le retard de propagation, la fréquence Doppler et la 15 direction d'arrivée d'un signal. Selon un aspect particulier de l'invention, l'étape d'estimer les paramètres représentatifs du signal utile et des signaux multi-trajet est exécutée itérativement pour chaque signal multi-trajet et le signal utile et comprend les sous-étapes de : 20 Soustraire, à la fonction d'intercorrélation spatio-temporelle, le ou les modèles paramétriques de signaux estimés aux itérations précédentes, - Estimer les paramètres représentatifs d'un signal multi-trajet ou du signal utile par application d'un algorithme du maximum de 25 vraisemblance au résultat de la soustraction précédente. Selon un aspect particulier de l'invention, l'étape d'estimer les paramètres représentatifs d'un signal multi-trajet ou du signal utile comprend les sous-étapes de : - Estimer la direction d'arrivée du signal en exploitant conjointement 30 tous les capteurs, 4 - Estimer l'amplitude complexe du signal indépendamment pour chaque capteur, Projeter l'amplitude complexe estimée sur le sous-espace défini par la direction d'arrivée estimée. 5 L'invention a aussi pour objet un dispositif de localisation apte à discriminer un signal utile de signaux multi-trajets, le dispositif comprenant au moins deux capteurs de signaux de radiolocalisation émis par un satellite, chacun desdits capteurs étant connecté à au moins deux lignes de traitement 10 des signaux, chaque ligne de traitement comprenant un corrélateur et au moins une ligne à retard apte à introduire un retard correspondant à un multiple de la durée d'intégration des corrélateurs, le dispositif comprenant des moyens de traitement des signaux configurés pour mettre en oeuvre un algorithme d'estimation du maximum de vraisemblance pour estimer des paramètres représentatifs du signal utile et des signaux multi-trajet, à partir des signaux pris en sortie des lignes à retard, lesdits paramètres représentatifs incluant au moins une amplitude complexe estimée indépendamment pour chaque capteur. Selon un aspect particulier de l'invention, un corrélateur d'une 20 première ligne de traitement est espacé d'une fraction de la période symbole du code d'étalement par rapport au corrélateur d'une seconde ligne de traitement. Selon un aspect particulier de l'invention, l'algorithme d'estimation du maximum de vraisemblance est de type « SAGE ». 25 D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : - La figure 1, un mobile recevant les signaux de radiolocalisation émis 30 par un satellite dans un environnement urbain, 5 - La figure 2, un schéma d'une architecture de dispositif de localisation selon l'invention, La figure 3, plusieurs diagrammes illustrant une fonction de corrélation mesurée sur chaque voie de réception pour différents codes 5 d'étalement PRN, - La figure 4, un organigramme décrivant les étapes principales du procédé selon l'invention. On reprend ici la description du dispositif de localisation décrit dans la 10 demande WO 2012/025306. Le dispositif de localisation selon la présente invention est basée sur le dispositif décrit antérieurement dans la demande précitée et apporte des améliorations par rapport à celui-ci. La figure 2 décrit une architecture d'un dispositif de localisation selon 15 l'invention. Le dispositif de localisation selon l'invention comporte une pluralité de capteurs A1, A2, An. Ces capteurs sont des antennes élémentaires conçues pour capter les signaux de radiolocalisation émis par les satellites d'un système de navigation tel que le système GPS ou Galileo. Le réseau de capteurs A1, A2, An permet d'obtenir l'information directionnelle par échantillonnage spatial du front d'onde reçu. Le dispositif de localisation comporte également des moyens de traitement des signaux reçus par l'ensemble des capteurs A1, A2 à An. Des 25 lignes de traitement de signal sont connectées en sortie de chaque capteur. Nous décrivons dans ce paragraphe la chaîne de traitement du signal comprise entre la sortie du capteur A1 et les moyens de calcul mettant en oeuvre l'algorithme du maximum de vraisemblance. La présence de l'étage RF, de descente en fréquence et d'un éventuel convertisseur analogique-30 numérique en sortie de chaque antenne n'étant pas utile à la compréhension de l'invention, ces éléments ne sont pas représentés sur la figure 2, le lecteur 6 pourra les considérer comme faisant partie des capteurs. Le capteur Al est relié à une première ligne comportant successivement un corrélateur C11 (constitué d'un multiplicateur et d'un intégrateur), puis une première ligne de retard R111 et une mième ligne de retard Rllm, les dites lignes étant 5 connectées en parallèle. La sortie des lignes de retard est connectée à des moyens de calcul pouvant mettre en oeuvre un algorithme de maximum de vraisemblance. L'ensemble des corrélateurs permet d'obtenir une estimation de la 10 fonction d'intercorrélation entre le code local (généré dans le récepteur) et le code reçu afin de pouvoir estimer les retards des différents signaux. Les lignes à retard vont échantillonner le signal post-corrélation avec une période d'échantillonnage de Tint afin de pouvoir estimer les fréquences Doppler des différents signaux. 15 A la suite du traitement de corrélation, des lignes de retard R111 et R1 1 m sont connectées en sortie de l'intégrateur. Ces lignes de retard sont disposées en parallèles. L'écart introduit par chaque ligne de retard correspond toujours à un multiple de la durée d'intégration du corrélateur. A 20 titre d'exemple non limitatif, la première ligne de retard R111 n'introduit pas de retard, un retard de Tint est introduit par la seconde ligne de retard R112 et la mième ligne de retard R11 m introduit un retard de (m-1)Tint. La ligne de traitement peut comporter un nombre m de lignes de retard en parallèle, chaque ligne de retard étant espacée d'une durée Tint correspondant à la durée d'intégration des corrélateurs. Les lignes de retard peuvent être élaborées par tout moyen de traitement du signal en analogique comme numérique. De même, le capteur Al est relié à une seconde ligne de traitement comportant successivement un corrélateur C12 espacé d'une fraction de la période symbole du code d'étalement (durée notée Te sur la figure 2) par rapport au corrélateur C11 et m lignes de retard R121 à R12m agencées de la même façon que la première ligne de traitement. Le capteur 7 Al peut être relié à des troisième et quatrième lignes de traitement et voire plus qui ne sont pas présentées sur la figure 2. Les lignes de traitement comportant les corrélateurs Cl 1 et C12 sont identiques. Plus généralement, les structures des lignes de traitement en sortie de chaque antenne sont 5 identiques entre elles. Selon un mode de réalisation plus simple, les lignes de traitement peuvent ne pas comporter de lignes de retard. Dans ce cas, les sorties des corrélateurs sont directement connectées aux moyens de calcul pP mettant en oeuvre un algorithme de maximum de vraisemblance. 10 Le dispositif de localisation comporte une pluralité de capteurs Al à An et chaque capteur est connecté à une ligne de traitement telle que décrite précédemment à la suite du capteur Al. Les capteurs et les lignes de traitement, comportant les lignes de retard, forment un réseau spatio- 15 temporel multi-corrélateurs. Les lignes de retard sont reliées directement sur de multiples entrées à des moyens de calcul et de traitement pP mettant en oeuvre un algorithme de maximum de vraisemblance. Plusieurs types d'algorithme de maximum de vraisemblance peuvent être mis en oeuvre. De préférence, un algorithme 20 de type SAGE traite les signaux issus des lignes de retard dans le but d'estimer les paramètres caractéristiques des signaux reçus (direction incidente en azimut et élévation, retard, Doppler par exemple). L'estimation de ces paramètres permet de discriminer le signal utile des signaux multitrajets. L'algorithme de maximum de vraisemblance est apte à traiter les 25 signaux provenant du système implémenté selon une architecture spatio-temporelle multi-corrélateurs. Dans la présente invention, contrairement aux hypothèses prises dans la demande WO 2012/025306, on suppose que les capteurs Al, A2, An sont 30 tous différents et présentent ainsi des diagrammes de rayonnement différents. 8 Un exemple de modules de fonctions de corrélation mesurées au moyen d'un dispositif tel que celui de la figure 2 est illustré sur la figure 3 pour un dispositif comprenant quatre antennes, pour différentes séquences 5 de codes d'étalement PRN au même instant. Sur chacun des six diagrammes de la figure 3, on a représenté les fonctions de corrélation obtenues pour chacune des quatre voies de réception. On remarque sur ces figures que les fonctions de corrélation mesurées sont différentes pour chaque voie de réception et sont aussi différentes selon le code d'étalement 10 PRN utilisé. Il résulte de ces deux diagrammes que les gains diffèrent entre les différentes antennes d'un dispositif en fonction de la direction d'arrivée des signaux. Ainsi, l'application de la méthode d'estimation des paramètres des signaux multi-trajets, telle que décrite dans la demande WO 2012/025306, 15 est inopérante car elle conduit à faire virtuellement apparaitre des signaux non reçus. En effet, cette méthode considère que les gains d'antennes sont identiques pour chacune des voies de réception. La figure 3 identifie en outre sous forme de symboles « + » les niveaux théoriques reçus, prenant en compte des mesures en chambre 20 anéchoïque, des diagrammes de rayonnement des antennes. On peut y vérifier que ces niveaux ne se superposent pas aux pics de corrélation des signaux mesurés et que le classement des voies n'est pas non plus reproduit. Enfin, cette figure identifie sous forme de symboles « * » les niveaux 25 estimés des signaux reçus au moyen de la présente invention. On peut y vérifier cette fois que les niveaux estimés correspondent bien aux niveaux reçus sur les quatre voies de réception, ce qui valide l'apport de l'invention. On détaille à présent les étapes du calcul des paramètres des signaux 30 (trajet direct et multi-trajets) reçus par le dispositif de la figure 2. 9 L'expression y(t) du signal ramené en bande de base en sortie d'antennes identiques peut se mettre sous la forme : L-1 y (t) =E a(0 , 1)x y1 x exp(2»-c.v/.t)x c(t - r) + b(t) 1=0 où l'indice « / » correspond au nombre de trajets reçus pour un même signal 5 (trajet direct et échos), le numéro '0' étant affecté au trajet direct, a(0,) représente le vecteur directionnel du réseau d'antennes en fonction des angles d'arrivée 0 et ço de chaque signal. n représente l'amplitude complexe du signal, vi sa fréquence Doppler, ri son retard de propagation, c(.) le code d'étalement et enfin b(.) représente le vecteur bruit thermique en sortie de 10 chaque antenne. Ce signal traverse ensuite les corrélateurs représentés sur la figure 2. L'expression X du signal reçu post-corrélation pour l'architecture STAP multi-corrélateurs proposée peut se mettre sous la forme vectorisée suivante : L-1 X=EX/(Nri)+b',pjit i=o avec '11, =Vi,t9i,go,,rrovrif On rappelle que m représente le nombre de lignes à retard par ligne de traitement, p désigne le nombre total de lignes de traitements par antenne et n représente le nombre d'antennes. bmp,Tint représente le vecteur bruit thermique en sortie de chaque antenne 20 intégré sur une durée Tint, IP représente le vecteur des paramètres à estimer pour le trajet 1 et X, Xi et bmp, Tint sont de dimension mxnxp. On note rd et vri les retards et décalages Doppler relatifs du trajet 'I ' par rapport au trajet direct. Le corrélateur est asservi sur le trajet direct et commet une erreur Er sur le code (retard compté positivement) par rapport au 25 trajet de référence d'indice '0' et ev sur la fréquence du trajet direct. 7, 15 10 désigne l'amplitude complexe du signal post corrélation. La notation OT représente la transposée d'un vecteur, (.)H représente la transposée-conjuguée d'un vecteur et l'opérateur ce représente le produit de Kronecker. 5 La matrice de corrélation temporelle fL, du signal est construite de la façon suivante : Dans un premier temps, les sorties multi-corrélateurs sont concaténées dans un vecteur colonne pour reconstruire une fonction d'intercorrélation échantillonnée du signal reçu avec le code local. 10 Ces fonctions d'intercorrélation obtenues pour différents instants post-corrélation sont concaténées toujours dans un vecteur colonne de manière à tracer l'évolution temporelle des fonctions d'intercorrélation. Cette évolution temporelle permettra de caractériser les Doppler relatifs des échos. 15 1-1 c(e r r rl ev -1/ r1)= r(s, -r + +(I-=-1) r(er -z- e - T2 -v r(er -r 2 r(er - r 2 r(e -(P-1) e )x exp[- jrc.(e, ,.1).Tint] 2 x exp[-d)lint} -1/ Te) 2 e) x exp[- v 2je-.(e, - ).2.Tint 2 )x exp[- r1).2.Tint 2jg.(ev -v 2 e (1=1-1-2T)xexpE- 2j1r.(e, 2 e avec r() la fonction d'intercorrélation d'un code reçu avec sa réplique locale. Les décalages Doppler relatifs sont supposés constants sur la durée 20 de traitement m. Tint, les retards relatifs évoluent comme : yr/ rrlm = r /int où fc désigne la fréquence porteuse du signal. , Jc 11 Selon l'architecture spatio-temporelle multi-corrélateurs proposée sur la figure 2, la matrice d'autocorrélation temporelle î, du bruit post-corrélation est égale à celle du code local pré-corrélation et est indépendante du Doppler, elle s'écrit : r(0) r (Te) r(0) A A r (p .Te) r (T e) r (T e) = M r (T e) r(0) r (T e) r(0) M r (T e) r (T e) r(p.Te) A A r(Te) r(0) En supposant le bruit blanc spatialement, la matrice d'autocorrélation spatio-temporelle du bruit post corrélation s'écrit : E(b Np jint X b 11\11j = C = X I mn 141, Enfin, on note itcH 15 kif (er - 2-,,,ev-vd)=14cH (Sr - T ri,Cv -Vri)Xkin Oftp-11 soit : 5 10 12 (p -1) \H x x exp[- 2pr.(s, - V,. r(er + -Te) 2 (19- l) r(er r1,1 Te) 2 r(e, r(e, 11,2 C rr/ ev Vr/ ) = x x exp[- 2pr.(e, ] -T r(er -rri,N +-Te) 2 x x exp[- 2j21-.(E, (p -1) Te) r(er -rd,N 2 Le principe de l'algorithme SAGE consiste à décomposer le signal reçu sur l'ensemble des trajets et à estimer, par maximum de vraisemblance, 5 les paramètres de chaque trajet. Pour l'architecture et le modèle de signal proposé, l'opposé de la fonction Log-vraisemblance A1(,) pour le trajet 1 s'écrit à une constante et un coefficient multiplicatif prés : 10 v1)=(X1(y1)-5(Ir xC-1x(Xl(m11)-50 où la notation A désigne l'estimée d'une grandeur. On montre alors que maximiser la vraisemblance pour le trajet revient à maximiser le terme 15 À/(y/) définit par : 2 (01,i) 0 itici (er -rd,ev -y À1(411)= mx cH (er -Tr1,ev-VrI)X11C(er -rd, v ri) 13 La phase « E-Step » de l'algorithme SAGE consiste à isoler un trajet particulier, et la phase « M-Step » consiste à estimer les paramètres du trajet par maximum de vraisemblance. On boucle alors sur l'ensemble des trajets 5 puis on reboucle itérativement jusqu'à convergence de l'algorithme. Le critère de convergence porte généralement sur la norme de la mise à jour du vecteur des paramètres à estimer. On peut également ajouter un critère sur l'amplitude du trajet à estimer pour le distinguer du bruit. 10 La mise en équation des phases « E-Step » et « M-Step » en mode poursuite est la suivante. Lors de la phase « E-Step », le modèle paramétrique des signaux estimés aux itérations précédentes est soustrait aux signaux reçus de la façon suivante : L-1 15 *4 = X- E x, /.0 Le modèle paramétrique de signal est mis à jour de la façon suivante pour prendre en compte des amplitudes du signal reçu potentiellement différentes sur les différentes antennes de réception : = , x,(wi)=7, RcVr -rd,ev-vrI) 20 Sous cette dernière forme, l'amplitude complexe n'est plus un scalaire mais un vecteur de dimension égale au nombre de voies de réception du dispositif selon l'invention. Lors de la phase « M-Step », les paramètres du signal multi-trajet 25 suivant (parmi les L-1 signaux multi-trajets et le signal utile) sont estimés de la façon suivante. En premier lieu, on recherche le retard et le Doppler relatifs du signal à estimer. Afin de réduire la charge calculatoire, on construit dans un premier 14 temps le vecteur Y1 qui exploite les estimées des DOA (direction d'arrivée des signaux), puis on maximise les vraisemblances pour le retard et le Doppler relatif au voisinage des solutions de l'itération précédente : 5 Yro = ia(br, 0/ )11 ttrio = arg max ( litCH (er -12r1)XY1012 itCH (er rrl,êv «i>s r1)X1-1C(8 rl v -12 d) r10 = arg max 12 Vrl Yloi kif (êr î rl e v - V r1)XîiC (esr ri e v - V d) Ensuite, on recherche les directions d'arrivée (DOA) du signal à 10 estimer sur les différents capteurs. Afin de réduire la charge calculatoire, on construit dans un premier temps le vecteur Zio qui exploite les estimées des retards et Doppler relatifs, puis on maximise les vraisemblances pour les DOA au voisinage des solutions de l'itération précédente : =fr (s) 'cf rl é'1, eri x 51/. 15 Ô, = arg max ia(0, , x Zio "04 = argmaxla(Ô, , ça 1)H x Zio Enfin, on détermine l'amplitude complexe post-corrélation du signal à estimer pour chaque capteur en supposant les amplitudes différentes, a priori, pour chaque capteur. 20 Pour cela, on effectue une estimation indépendante des amplitudes complexes sur les différentes antennes : frm kc,(êT -erjr X ÎLtemp = (e^ I.% -"; "s* c r rl V - rl)" rl f". V - ri) 15 Puis, on projette le vecteur résultant sur le sous-espace défini par la direction d'arrivée estimée : 5 al =a(él,01)-Fi _templ x jemp)X al aH xa Où l'opérateur « - » désigne un produit de vecteurs terme à terme et ou l'amplitude complexe estimée à la dernière étape n'est plus un scalaire mais un vecteur de dimension le nombre de voies de réception du dispositif 10 selon l'invention. Autrement dit, on conserve, dans l'estimation indépendante de l'amplitude complexe, la composante correspondant à la direction d'arrivée estimée à l'étape précédente. Ainsi, on privilégie la direction d'arrivée estimée. 15 En pratique, il est entendu que I ne sera pas à recalculer numériquement à la volée lors du traitement par le récepteur afin de ne pas alourdir la charge calculatoire. Au pire, des valeurs numériques pourront être pré-calculées et stockées, voire, sous certaines hypothèses, une solution analytique pourra être proposée. 20 Par exemple, si on néglige l'évolution du retard en fonction du Doppler relatif (ce qui revient à négliger une variation de l'ordre de 6.1040 s pour un Doppler de 100 Hz et une durée d'observation de 10 ms) et que le terme (Er -rd) est multiple du pas de temps entre les corrélateurs, alors chaque 25 bloc (terme correspondant à 1 ligne à retard post-corrélation) de la matrice se ramène à un Dirac multiplié par un terme de phase : 16 8p (k) x exp[- 2 j71-.(e, - yr, ) int] T p (k) x exp[- 2 j .(e, - v ).2.Tint ] k(k.Te,s, 8p (k) x exp{- 2 ft r .(e v - v ri) .m .7' int] avec : 5 p (k) = [0 A 0 1 0 A Or A k-1 k k+1 A /21 2 2 La figure 4 schématise, sur un organigramme, les étapes principales du procédé d'estimation des paramètres des signaux multi-trajets et du signal utile, selon l'invention. 10 Une première étape 401 du procédé consiste à corréler le signal reçu par les capteurs avec un code local au moyen des corrélateurs du dispositif. Une deuxième étape 402 du procédé consiste à construire, pour chaque capteur, une fonction d'intercorrélation échantillonnée du signal reçu avec le code local. 15 Une troisième étape 403 du procédé consiste à déterminer une fonction d'intercorrélation spatio-temporelle à partir de la concaténation des fonctions d'intercorrélation obtenues à l'étape précédente pour chaque capteur. Puis, on estime, de façon itérative, les paramètres représentatifs du 20 signal utile et des signaux multi-trajet par application d'un algorithme de maximum de vraisemblance. Pour cela, à chaque itération, on soustrait lors de l'étape 404, à la fonction d'intercorrélation spatio-temporelle, le ou les modèles paramétriques de signaux dont les paramètres représentatifs ont été estimés aux itérations 25 précédentes. 17 Puis, on estime à l'étape 405 les paramètres représentatifs d'un signal multi-trajet ou du signal utile par application d'un algorithme du maximum de vraisemblance au résultat de la soustraction précédente. Cette étape 405 comprend au moins les sous-étapes suivantes : 5 - Estimer la direction d'arrivée du signal en exploitant conjointement tous les capteurs, Estimer l'amplitude complexe du signal indépendamment pour chaque capteur, Projeter l'amplitude complexe estimée sur le sous-espace défini 1 o par la direction d'arrivée estimée. L'invention permet de discriminer le signal utile des signaux multitrajets en prenant en compte les caractéristiques différentes des voies de réception du dispositif. Cet avantage est obtenu de par l'estimation 15 indépendante de l'amplitude complexe du signal reçu par chaque capteur. L'invention s'applique aux stations sol de référence ou d'observation des constellations de satellites de navigation. Cette invention peut être utilisée dans des terminaux mobiles pour améliorer la solution de navigation dans des environnements urbains, REVENDICATIONS1. Procédé d'estimation des paramètres de signal utile et de signaux multitrajet originaires d'un signal de radiolocalisation émis par un satellite, au moyen d'un dispositif de localisation comprenant au moins deux capteurs aptes à recevoir ledit signal, le procédé comprenant les étapes de : - Corréler (401) le signal reçu par lesdits capteurs avec un code local au moyen de corrélateurs, Construire (402), pour chaque capteur, une fonction d'intercorrélation échantillonnée du signal reçu avec le code local, - Déterminer (403) une fonction d'intercorrélation spatio-temporelle à partir de la concaténation des fonctions d'intercorrélation obtenues à l'étape précédente pour chaque capteur, Estimer (404,405) des paramètres représentatifs du signal utile et des signaux multi-trajet par application d'un algorithme de maximum de vraisemblance, lesdits paramètres représentatifs incluant au moins une amplitude complexe estimée indépendamment pour chaque capteur. 2. Procédé d'estimation selon la revendication 1 dans lequel lesdits paramètres représentatifs incluent le retard de propagation, la fréquence Doppler et la direction d'arrivée d'un signal. 3. Procédé d'estimation selon l'une des revendications précédentes dans lequel l'étape d'estimer les paramètres représentatifs du signal utile et des signaux multi-trajet est exécutée itérativement pour chaque signal multi-trajet et le signal utile et comprend les sous-étapes de : Soustraire (404), à la fonction d'intercorrélation spatio-temporelle, le ou les modèles paramétriques de signaux estimés aux itérations précédentes, 19 Estimer (405) les paramètres représentatifs d'un signal multi-trajet ou du signal utile par application d'un algorithme du maximum de vraisemblance au résultat de la soustraction précédente. 5 4. Procédé d'estimation selon la revendication 3 dans lequel l'étape d'estimer (405) les paramètres représentatifs d'un signal multi-trajet ou du signal utile comprend les sous-étapes de : Estimer la direction d'arrivée du signal en exploitant conjointement tous les capteurs, io Estimer l'amplitude complexe du signal indépendamment pour chaque capteur, Projeter l'amplitude complexe estimée sur le sous-espace défini par la direction d'arrivée estimée. 15 5. Dispositif de localisation apte à discriminer un signal utile de signaux multi-trajets, le dispositif comprenant au moins deux capteurs (Al ,A2) de signaux de radiolocalisation émis par un satellite, chacun desdits capteurs étant connecté à au moins deux lignes de traitement des signaux, chaque ligne de traitement comprenant un corrélateur (C11 ,C12) 20 et au moins une ligne à retard (R11m) apte à introduire un retard correspondant à un multiple de la durée d'intégration des corrélateurs, le dispositif comprenant des moyens de traitement des signaux (pP) configurés pour mettre en oeuvre un algorithme d'estimation du maximum de vraisemblance pour estimer des paramètres représentatifs du signal 25 utile et des signaux multi-trajet, à partir des signaux pris en sortie des lignes à retard, lesdits paramètres représentatifs incluant au moins une amplitude complexe estimée indépendamment pour chaque capteur. 6. Dispositif de localisation selon la revendication 5 dans lequel un 30 corrélateur (C11) d'une première ligne de traitement est espacé d'une 20 fraction de la période symbole du code d'étalement par rapport au corrélateur (C12) d'une seconde ligne de traitement. 7. Dispositif de localisation selon l'une des revendications 5 ou 6 dans 5 lequel l'algorithme d'estimation du maximum de vraisemblance est de type « SAGE ». 10