Brevet FR2818468A1 — publié le 20020621

La présente invention concerne une méthode et un dispositif de détection multiutilisateur. Plus particulièrement, la présente invention concerne une méthode et un dispositif de détection multi-utilisateur à maximum de vraisemblance pour un système de télécommunication DS-CDMA (Direct Séquence Code Division Multiple Access). Dans un système de télécommunication mobile DS-CDMA, la séparation des communications en provenance ou à destination des différents utilisateurs est réalisée en multipliant chaque symbole complexe d'un utilisateur par une séquence d'étalement qui propre à ce dernier, encore appelée pour cette raison signature de l'utilisateur. La fréquence d'étalement (chip rate) étant supérieure à la fréquence des symboles, le signal transmis par chaque utilisateur est distribué (ou étalé) dans l'espace des fréquences. Le rapport entre la bande occupée par le signal étalé et la bande occupée par le signal d'information est appelé facteur d'étalement. A la réception, la séparation d'un utilisateur donné est obtenue grâce à un filtrage adapté à la signature correspondante. Lorsque le canal de transmission présente une pluralité de trajets de propagation, la sortie du filtrage adapté comporte autant de pics de corrélation. Chaque trajet du canal peut être modélisé par un coefficient multiplicatif complexe et un retard. Les signaux s'étant propagés selon les différents trajets peuvent être alignés et combinés au moyen de coefficients complexes conjugués des coefficients de trajet, réalisant ainsi un filtrage adapté au canal de transmission. Pour simplifier la terminologie nous engloberons dans l'expression générale filtrage adapté à l'utilisateur k à la fois l'opération de filtrage adapté à la signature de l'utilisateur k et celle de filtrage adapté au canal de transmission. Pour combattre l'interférence entre signaux à destination (liaison descendante) ou en provenance (liaison montante) des différents utilisateurs, il a été proposé des méthodes de détection multi-utilisateur et notamment des méthodes de détection itératives telles que celles connues sous le nom de PIC (Parallel Interference Cancellation) et SIC (Serial Interference Cancellation). Elles sont fondées sur l'itération d'un cycle d'élimination des interférences comportant l'estimation des symboles émis, l'évaluation des interférences et leur soustraction des signaux reçus. Bien que performantes, ces méthodes ne sont pas optimales dans la mesure où elles ne fournissent pas une estimation au sens du maximum de vraisemblance des symboles transmis par les différents utilisateurs. Une méthode de détection multi-utilisateur à maximum de vraisemblance inspirée de l'algorithme de Viterbi a été proposée par S. Verdu dans un article intitulé Minimum probability of error for asynchronous Gaussian multiple access channels publié dans IEEE Transactions on Information Theory, pages 85-96, Janvier 1986 mais sa complexité est prohibitive car elle varie de manière exponentielle avec le nombre d'utilisateurs . Plus récemment il a été proposé par L. Brunel et al., dans un article intitulé Euclidian space lattice decoding for joint detection in CDMA system publié dans Proceedings of ITW, page 129, Juin 1999, une méthode de détection multi-utilisateur à maximum de vraisemblance utilisant une représentation par réseau de points. Selon cette méthode on détermine un vecteur caractéristique du signal reçu représentant une statistique suffisante pour la détection au maximum de vraisemblance des symboles transmis par les différents utilisateurs. On montre sous certaines conditions que le vecteur caractéristique peut être représenté comme le point d'un réseau perturbé par un bruit. La détection consiste alors à rechercher le point du réseau le plus proche du point correspondant au vecteur reçu. Cependant, la dimension du réseau à utiliser étant de 2.K où K est le nombre d'utilisateurs le nombre de points reste encore très élevé. Pour simplifier la détection, il a été proposé de limiter la recherche du plus proche voisin aux points du réseau appartenant à une sphère centrée autour du point reçu. Nous exposerons ci-après cette méthode de détection simplifiée dite méthode de détection par sphères : Nous plaçons dans le contexte d'un système de télécommunication mobile à accès multiple avec étalement de spectre par séquences directes (DS-CDMA) comprenant K utilisateurs communiquant de manière synchrone avec une station de base. Soit dk(i) le symbole complexe émis par l'utilisateur k à l'instant i. Ce symbole appartient à la constellation de modulation Ak utilisée par l'utilisateur k, que l'on appellera encore alphabet de symboles de l'utilisateur k Chaque utilisateur k transmet un bloc de N symboles avec une amplitude du signal ak . Les symboles sont étalés par une signature complexe sk (t) = SR (t) + j.sk (t) de durée égale à la période symbole T : Les K symboles complexes dk(i)=dk(i)+j.dk(i) transmis à l'instant i sont placés dans un vecteur ligne de valeurs réelles d2(i) défini comme : Du fait de la synchronisation des utilisateurs, le signal modulé correspondant peut s'écrire, en fonction du temps t : Nous supposons que le canal est un canal idéal à bruit blanc additif gaussien. Soit 7i=St+rjt le signal reçu au temps t et 17, un bruit gaussien complexe de moyenne nulle et dont les composantes ont une variance No. Soit les vecteurs lignes y(0=(),..j)) et y2(i(y; (iy; (i...,y;(i),y(i)) où yk(i)=yk (i)+j.yk(t est la sortie complexe à l'instant i du filtre adapté à l'utilisateur k : T T avec Rlk= fSI(t)Sk*(t)dt=Rt1+ j.RA pour Ie,.e = 1,...,K et m(i)= .sk*(t-i.1)dt 0 0 On notera R(i) la matrice d'autocorrélation des séquences d'étalement. Si l'on décompose les éléments complexes de (3) en leurs parties réelles et imaginaires, on obtient : Soient A2=Diaai,ai,...,o,a/:) et R2 la matrice de taille 2K x 2K telle que : L'équation (4) peut alors se mettre sous forme matricielle : où M2 est une matrice réelle de taille 2K x 2K définie par M2=A2R2 et où le vecteur de bruit n2(f(n(i(/){(i)r . {i\n'K(i)) a pour matrice de covariance jV0R2 . Nous démontrerons ci-après que y2(i), tel que donné par l'équation (6), peut être représenté comme un point d'un réseau A2 de dimension 2.K, de matrice génératrice M2 corrompu par un bruit n2. Nous appellerons réseau réel de points A de dimension # tout ensemble de vecteurs de R#vérifiant : x=bvl+mv2+...+b,.-vK où bEZ,tü=1,...,x et où {VI, V2," " ", V K} est une base sur R" . Un exemple de réseau de points de dimension 2 a été représenté en Fig. 1. Les points du réseau forment un sous-groupe abélien additif de RK ,c'est d'ailleurs le plus petit sous-groupe de R" contenant les vecteurs {v),V2,...,v} et un Zmodule de R" . Ces vecteurs de base forment les lignes de la matrice génératrice G du réseau. On peut donc écrire x = bG où b=(b,...,laZ" . (7) La région délimitée par les vecteurs de base est appelée parallélotope fondamental et son volume, noté vol(A) ou det(A), est nommé volume fondamental. Ce volume fondamental n'est autre que le module du produit vectoriel des # vecteurs de base et est donc égal à #det(G)# où det désigne le déterminant. S'il existe plusieurs choix possibles pour la matrice génératrice d'un même réseau, il n'existe par contre qu'une unique valeur pour le volume fondamental. La région de Voronoï V ou cellule de Dirichlet d'un point x appartenant au réseau est l'ensemble des points de R# plus proches de x que de tout autre point du réseau. Le volume de cette région est égal au volume fondamental. Le rayon d'empilement p du réseau est le rayon de la plus grande sphère inscrite dans la région de Voronoï et le rayon de recouvrement celui de la plus petite sphère circonscrite à cette même région. Le rayon d'empilement est donc le rayon des sphères dont l'empilement constitue le réseau de points et le rayon de recouvrement est celui des sphères les plus petites qui, centrées sur les points du réseau, permettent de recouvrir tout l'espace R#. La densité du réseau est le rapport entre le volume de la sphère de rayon p et le volume fondamental. Enfin, le coefficient d'erreur (kissing number) #(#) du réseau est le nombre de sphères tangentes à une même sphère dans l'empilement ou, en d'autres termes, le nombre de voisins d'un point du réseau, situés à la distance minimale dEmun = 2p . Considérons à nouveau l'équation (6). Les nombres complexes Supposons par exemple que les composantes dk (i) et dlk\i) soient des symboles de modulation PAM d'ordre M : Si l'on effectue la transformation : df(i)=1(dfV)tM-l) et di(i)=1(dl(i)tM-l) soit encore de manière vectorielle : les composantes df(i) et dV) sont des éléments de Z et par suite d'2V) est un vecteur de Z2K. De manière générale, s'il existe une transformation affine transformant les composantes duif) et dkl(i) en des éléments de Z, on peut représenter le vecteur d'2 (i) par un vecteur de Z2K De manière similaire, on effectue la transformation correspondante sur Y2(i) c'est-à-dire : Moyennant cette transformation, que nous supposerons implicite dans la suite, le vecteur d2(1)M2(i) appartient alors à un réseau de points A2 de dimension 2.K tel que défini par la relation (7) avec G=M2(i) . Le vecteur Y2(i) peut donc être considéré comme un point du réseau A2 corrompu par un bruit n2(i). Si l'on suppose que les composantes du vecteur de bruit n2(i) sont des variables aléatoires indépendantes gaussiennes centrées, le problème de la détection au sens du maximum de vraisemblance des symboles émis par les différents utilisateurs se ramène à la recherche du point z2 du réseau A2 tel que sa distance à Y2(i) soit minimale. En réalité, les composantes du vecteur de bruit n2(i) sont corrélées et la matrice de covariance de n2(i) est NoR2 . Afin de se ramener au cas décorrélé il faut effectuer préalablement au décodage une opération de blanchiment du bruit. La matrice R étant hermitienne, la matrice d'autocorrélation R2 est symétrique définie positive et peut donc faire l'objet d'une factorisation de Cholesky : R2=W2W2I (13) où W2 est une matrice triangulaire inférieure de taille 2Kx2K. On définit un vecteur d'observation blanchie : 2(i)=y2(i)W2T-' (14) ainsi qu'un nouveau réseau de points S22 constitué des vecteurs de composantes (xt(i),XI(i),...,xf(i1x (i)) avec xi(i)=X2(i)W2'~] où X2(z) est un vecteur de composantes (xi i.xi (i),. rxK(ICK(1)) appartenant à A2. On peut facilement montrer, qu'après blanchiment, la matrice de covariance du bruit filtré m(i)wr1 est égale à Nol2K où 12K est la matrice identité de dimension 2K. La détection comprend donc une première étape de blanchiment du vecteur d'observation suivie d'une étape de recherche du plus proche voisin au sein du réseau de points #2. Pour réduire le nombre de points à tester, comme illustré en Fig. 1, on peut limiter la recherche à une sphère centrée autour du point #2. En pratique, le choix du rayon de la sphère résulte d'un compromis : il ne doit pas être trop grand pour ne pas conduire à un nombre de points trop élevé et suffisamment grand pour inclure au moins le plus proche voisin. La Fig. 2 représente schématiquement un dispositif de détection multi-utilisateur utilisant une méthode de détection par sphères. Le signal reçu n est filtrée par une batterie de filtres adaptés à chacun des utilisateurs, 2101,...,210K. Les composantes réelles et imaginaires du vecteur d'observation Y2(i) en sortie des filtres adaptés sont transmis à une unité de calcul matriciel effectuant l'opération de blanchiment spectral selon la relation (14). Les composantes réelles et imaginaires du vecteur blanchi #2(i) sont ensuite transmises à une unité de détection par sphères recherchant le plus proche voisin du point reçu au sein du réseau #2 de dimension 2.K. Les coordonnées du plus proche voisin donnent directement les composantes réelles et imaginaires des symboles estimés #k(i) pour les différents utilisateurs. On peut montrer que la méthode de détection par sphères telle qu'exposée cidessus possède une complexité en O(K6) ce qui peut s'avérer très pénalisant lorsque le nombre d'utilisateurs est élevé. Le but de la présente invention est de proposer sous certaines conditions une simplification de la méthode de détection par sphères. A cette fin, l'invention est définie par une méthode de détection d'une pluralité de symboles (dk(i)) transmis par ou pour une pluralité K d'utilisateurs, chaque symbole appartenant à une constellation de modulation et faisant l'objet d'un étalement spectral par une séquence d'étalement, ladite méthode comprenant une étape de filtrage adapté pour fournir un vecteur complexe (y(0,y(ï)) caractéristique dudit signal reçu, ledit vecteur complexe étant décomposé en un premier vecteur (y (y (0 ) et un second vecteur (y (/),y (i) ) et au moins les plus proches voisins des premier et second vecteurs étant recherchés au sein d'un réseau de points (A,Q) généré par lesdites constellations de modulation, les symboles transmis étant estimés à partir des composantes desdits plus proches voisins. Avantageusement, les séquences d'étalement (sk(t)) sont constitués de multiples réels (sk(t) ) d'un même coefficient complexe (a) . Selon un mode de réalisation de l'invention, la recherche est limitée à un premier ensemble de points du réseau appartenant à une première zone prédéterminée (#R) autour du premier vecteur et à un second ensemble de points du réseau appartenant à une seconde zone prédéterminée (#I) autour du second vecteur. De même la recherche peut être limitée à un premier ensemble de points du réseau appartenant à une première zone prédéterminée (En) autour de l'origine et à un second ensemble de points du réseau appartenant à une seconde zone prédéterminée (SI) autour de l'origine. Lesdites première et seconde zones prédéterminées sont par exemple des sphères. Avantageusement, la recherche du plus proche voisin du premier vecteur est effectuée sur une pluralité de composantes de ce dernier, la recherche étant limitée pour chacune desdites composantes à un intervalle défini par une borne inférieure et une borne supérieure, lesdites bornes étant choisies de manière à ce que ledit intervalle ne comprenne pas des points relatifs à des symboles ne pouvant appartenir à la constellation de modulation. De même la recherche du plus proche voisin du second vecteur peut être effectuée sur une pluralité de composantes de ce dernier, la recherche étant limitée pour chacune desdites composantes à un intervalle défini par une borne inférieure et une borne supérieure, lesdites bornes étant choisies de manière à ce que ledit intervalle ne comprenne pas des points relatifs à des symboles ne pouvant appartenir à la constellation de modulation. Préalablement à la recherche du plus proche voisin, le premier vecteur (yR(i)) est avantageusement soumis à un traitement matriciel visant à substantiellement décorréler les différentes composantes de bruit de ce dernier. De même, préalablement à la recherche du plus proche voisin, le second vecteur (yI(i)) peut être soumis à un traitement matriciel visant à substantiellement décorréler les différentes composantes de bruit de ce dernier. Selon une variante de réalisation de l'invention, ladite étape de recherche est étendue à celle d'un premier ensemble de points plus proches voisins dudit premier vecteur, dits premiers voisins, et à un second ensemble de points plus proches voisins dudit second vecteur, dits seconds voisins, et les symboles transmis sont estimés de manière souple à partir des symboles générant lesdits premiers et second voisins et des distances séparant lesdits premiers voisins dudit premier vecteur d'une part et lesdits seconds voisins dudit second vecteur d'autre part. Selon un mode de réalisation particulier de l'invention, on détermine à partir des symboles estimés les contributions de chaque utilisateur aux signaux obtenus par l'étape de filtrage adapté et, pour un utilisateur k donné, on élimine en sortie de l'étape de filtrage adapté les contributions des autres utilisateurs correspondant aux symboles déjà estimés. Alternativement, on détermine les contributions de chaque utilisateur au signal reçu à partir des symboles estimés et, pour un utilisateur k donné, on élimine en entrée de l'étape de filtrage adapté les contributions des autres utilisateurs correspondant aux symboles déjà estimés. Si les symboles desdits K utilisateurs sont transmis de manière synchrone, ledit réseau de points sera préférablement de dimension K. Si les symboles desdits K utilisateurs sont transmis de manière asynchrone et se propagent selon une pluralité de trajets, la dimension du réseau sera préférablement égale au nombre de symboles des différents utilisateurs pouvant interférer et non encore estimés. La présente invention est également définie par un dispositif de détection d'une pluralité de symboles (dk(i transmis par ou pour une pluralité K d'utilisateurs, chaque symbole appartenant à une constellation de modulation et faisant l'objet d'un étalement spectral par une séquence d'étalement, le dispositif comprenant des moyens pour mettre en oeuvre la méthode exposée ci-dessus. Ce dispositif pourra en particulier utilisé dans un récepteur de système de télécommunication mobile DS-CDMA . Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, apparaîtront plus clairement à la lecture de la description suivante faite en relation avec les figures jointes, parmi lesquelles : La Fig. 1 représente un réseau de points utile à la méthode de détection employée dans le récepteur illustré en Fig. 2; La Fig. 2 représente de manière schématique la structure d'un récepteur DSCDMA multi-utilisateur utilisant une méthode de détection par sphères; La Fig. 3 représente de manière schématique la structure d'un dispositif de détection multi-utilisateur selon un premier mode de réalisation de l'invention. La Fig. 4 représente de manière schématique la structure d'un dispositif de détection multi-utilisateur selon un second mode de réalisation de l'invention. L'idée à la base de l'invention est d'abaisser la dimension du réseau de points en utilisant des séquences d'étalement de type particulier pour les différents utilisateurs. Considérons à nouveau un système de télécommunication DS-CDMA avec K utilisateurs synchrones. Si nous choisissons des séquences d'étalement sk(t) à valeurs réelles, les termes imaginaires de la matrice R2 et par suite de la matrice M2(i) sont nuls. Par suite, on peut modéliser le système par un réseau de points réel A de dimension K et de matrice génératrice M(i) : où y'),d"(/),n(0 (resp. yl(i),dl(i),n1(i)) sont les vecteurs constitués des parties réelles (resp. des parties imaginaires) des composantes de y(i),d(i),n(i) ; T M(i)=AR(i) où R(i) est la matrice constituée par les coefficients Rlk = fSf(I)Sk(t)dt et A o est le vecteur des amplitudes des K utilisateurs. Les vecteurs d'observation yR(i) et yl(i) appartiennent à RK. Après transformation éventuelle selon une relation du type de celle de (12), les vecteurs yR(i) et yI(i) peuvent être considérés comme des points d'un réseau A de matrice génératrice M(i) corrompus par du bruit. On montre aisément que les vecteurs de bruit n R (j) et n' (i) ont tous deux pour matrice de covariance No.R(i). R(i) étant une matrice symétrique définie positive on peut la factoriser selon une décomposition de Cholesky : R = WWT où W est une matrice réelle triangulaire inférieure de taille KxK. Afin de décorréler les composantes de bruit, les vecteurs d'observation réels yR(i)et yI(i) sont tout d'abord soumis à une opération de blanchiment : Dans un deuxième temps, on recherche les plus proches voisins des vecteurs yR(i) et #I(i) appartenant au réseau de points Q constitué des vecteurs xl=x(i)w,- où X(i) appartient à A. On peut montrer facilement, qu'après blanchiment, les matrices de covariance des bruits filtrés nR (i)W T-1 et n7(/)\Vr sont toutes deux égales à N0Ik où IK est la matrice identité de dimension K. On voit donc que le fait d'utiliser des signatures réelles conduit à une recherche de deux plus proches voisins dans un même réseau de dimension K alors que dans le cas général, complexe, le décodage nécessite une recherche dans un réseau de dimension 2K. En fait on peut facilement montrer que le résultat précédent s'étend à tout jeu de signatures portées un même nombre complexe, c'est-à-dire telles que : sk(t)=ask(t) où a est un nombre complexe et sf(t) est réel. Ce sera bien entendu en particulier le cas si toutes les signatures sont imaginaires. La recherche du plus proche voisin pour le vecteur #I(i) se faisant selon le même principe que pour #R(i), nous n'exposerons que la première. Cette recherche consiste à déterminer le point x minimisant la métrique : où Y'=x+n' et x=(xi,...,xx) est un point appartenant au réseau Q. Alternativement, on notera que le vecteur yR(i) n'a pas besoin d'être blanchi si l'on utilise une métrique basée sur la matrice de covariance : On désignera par la suite, pour des raisons de simplification, z le vecteur d'observation blanchi (y"(0) ou non (yl?(i et Il.11 la métrique intervenant dans la relation (19) ou (20). Les points du réseau Q sont constitués des vecteurs x tels que x=bG où G est la matrice génératrice du réseau et b=la,...,1), les composantes b, appartiennent à l'anneau des entiers Z . Le détecteur restreint avantageusement son calcul de métrique aux points qui sont situés à l'intérieur d'une zone de la constellation située autour du point reçu, préférentiellement à l'intérieur d'une sphère de rayon donné #C centrée sur le point reçu z. Seuls les points du réseau situés à une distance quadratique inférieure à C du point reçu sont donc considérés pour la minimisation de la métrique. En pratique, le décodeur effectue la minimisation suivante : Pour ce faire, le décodeur recherche le plus petit vecteur w dans l'ensemble translaté z-#. On peut exprimer les vecteurs z et w comme suit : Il est important de noter que p et sont des vecteurs réels. Comme w=z-x où x appartient au réseau S2, on a la relation , = p, - b, pour i=1,...,K avec Le vecteur w est un point du réseau dont les coordonnées #I sont exprimées dans le repère translaté centré sur le point reçu z. Le vecteurw appartient à une sphère de rayon quadratique C centrée en 0 si : Dans le nouveau système de coordonnées défini par , la sphère de rayon quadratique C centrée en y est donc transformée en une ellipse centrée à l'origine. La factorisation de Cholesky de la matrice de Gram r=GGT donne r=AeT , où A est une matrice triangulaire inférieure d'éléments ##. Il faut noter que si le vecteur y a été blanchi, il n'y a pas lieu d'effectuer cette factorisation car la matrice génératrice de Q vaut AW et est donc déjà triangulaire inférieure. Toutefois, dans le cas où l'on n'a pas procédé au blanchiment préalable, la décomposition de Cholesky est nécessaire. Dans tous les cas, on peut écrire : En posant on obtient : En s'intéressant tout d'abord à la plage de variations possibles de #K, puis en rajoutant une à une les composantes, on obtient les K inégalités suivantes, qui définissent tous les points à l'intérieur de l'ellipse : On peut montrer que les inégalités (26) imposent aux composantes entières de b de satisfaire à : où #x# est le plus petit entier supérieur au réel x et #x#est le plus grand entier inférieur au réel x . Le décodeur possède K compteurs internes, à savoir un compteur par dimension, chaque compteur comptant entre une borne inférieure et supérieure comme indiqué en (27), étant entendu qu'à chaque compteur est associé un couple de bornes particulier. En pratique, ces bornes peuvent être mises à jour de manière récursive. Avantageusement, on liste toutes les valeurs du vecteur b pour lesquelles le point correspondant du réseau x = bG est situé en deçà de la distance quadratique C du point reçu. Les points du réseau situés en dehors de la sphère considérée ne sont pas testés. Avantageusement, pour chaque composante 1,..,K les bornes supérieure et inférieure de la recherche seront ajustées ne manière à ne pas contenir des points qui sont certainement en dehors de la constellation. Ainsi les compteurs ne perdent pas de temps à parcourir des points qui, de toutes façons, ne sont pas des solutions. Par exemple, si tous les utilisateurs emploient une même constellation de modulation PAM de taille M, les bornes de recherche ne pourront sortir de l'intervalle [0,M-1]. De plus, la recherche à l'intérieur de la sphère peut être considérablement accélérée en mettant à jour le rayon #C avec la dernière norme euclidienne calculée #w#. Finalement, on sélectionne comme meilleur point x celui associé à la plus petite norme #w#. De manière à être sûr que le décodeur trouve au moins un point du réseau, on choisit avantageusement un rayon de recherche supérieur au rayon de recouvrement du réseau. On peut par exemple le prendre égal à la borne supérieure de Rogers : où Vk est le volume d'une sphère de rayon unité dans l'espace réel. La Fig. 3 illustre schématiquement un dispositif de détection multi-utilisateur selon un premier mode de réalisation de l'invention. Le signal reçu est tout d'abord filtré par une batterie de filtres adaptés aux différents utilisateurs, 310i,.. ,310k. Le vecteur d'observation en sortie des filtres adaptés est décomposé en un vecteur d'observation réel yR(i)=(yR(i),...,yK(t) et un vecteur d'observation imaginaire yI(i)=(ylI(i),..,ykI(i)). Après transformation éventuelle du type de (12) (non représentée), les vecteurs y/l(i) et yl(i) subissent un blanchiment spectral en 320 et 321 afin de décorréler les échantillons de bruit. Les vecteur blanchis yR(i) et #I(i) font ensuite l'objet d'une recherche de plus proche voisin comme décrit plus haut dans les détecteurs par sphères 330 et 331. Le point trouvé par le détecteur 330 donne (moyennant éventuellement une transformation inverse à celle de (12)) les composantes réelles des symboles estimés pour les K utilisateurs. De même, le détecteur 331 donne les composantes imaginaires de ces symboles estimés. Au lieu de fournir directement des symboles de la constellation, le récepteur peut être adapté pour fournir les symboles sous forme de décisions souples. Dans ce cas, la recherche à l'intérieur de la sphère de détection ne se limite plus au plus proche voisin mais est étendue à une pluralité des plus proches voisins du point relatif au signal reçu. Plus précisément, soit #R et #l les sphères centrées respectivement autour de R(i) et y7(/). On associe à tout couple d'"'"'(i)=(vm,vm) de points voisins appartenant à (ER,El) tel que les K composantes d#'m(ï) appartiennent aux constellations de modulation des utilisateurs, une probabilité a posteriori pm,m , probabilité que le vecteur dk'"'(1) défini par ce point ait été émis, étant donné l'observation ni). Soit 0 l'ensemble de ces couples. On définit un symbole souple d'un utilisateur k comme le Mk-uplet (r,..., iMk) où Mk est le cardinal de la constellation de modulation de l'utilisateur k et où #j est la probabilité que le symbole sj ait été émis. On a : Les probabilités a posteriori pm,m pourront par exemple être exprimées en fonction des distances Am et ,m séparant les vecteurs YR (i) et Y'(i) de vm et vm L'équation (2) présupposait que les signaux des différents utilisateurs étaient synchrones. Lorsque cette hypothèse n'est pas valide, un symbole étalé d'un utilisateur k à un instant donné peut interférer avec deux symboles étalés successifs d'un autre utilisateur k'. Si nous supposons que la dispersion des retards de transmission #k des différents utilisateurs est inférieure à une période symbole T, le symbole d'un utilisateur k émis à l'instant i, dk(i) peut interférer avec les symboles dk'(i-l) et dk'(i+l) d'un utilisateur k'. On peut supposer sans perte de généralité que 0 vecteurs: le vecteur dp=(dl(z},...,dk-I(i),ìk(i-l ),dk+l(i-1 ),...,ìK(i-I)) relatif aux symboles déjà détectés, le vecteur (iF=(d1(1-Hl),...,dk-I(1+1),dk(1),dk+(1',..,dK(1))relatif aux symboles futurs et le vecteur d'observation y(i). On peut montrer que le vecteur d'observation complexe peut s'écrire sous la forme d'une contribution passée et d'une contribution future : où n est un vecteur de bruit de matrice de covariance N0RF, Rp et RF sont respectivement la matrice de corrélation des signatures avec les signatures passées et la matrice de corrélation des signatures avec les signatures futures. Lorsque les signatures sont à valeurs réelles (ou plus généralement sont des multiples réels d'un même nombre complexe) l'équation (30) peut se décomposer sous la forme de deux équations ne faisant intervenir que des vecteurs réels , de manière analogue aux équations (15) et (16) : Au lieu d'effectuer la détection par sphères sur les vecteurs d'observation yR et yI celle-ci opérera sur les vecteurs débarrassés des interférences dues aux symboles passés déjà estimés à savoir : :y R ~âpR APp et y'-dr'ARr . Alternativement, au lieu de procéder à l'élimination soustractive des interférences au niveau du vecteur d'observation, cette élimination peut être envisagée de manière équivalente, en amont, sur les signaux étalés . Les symboles estimés sont réétalés spectralement et les contributions des différents utilisateurs au signal reçu sont soustraites une à une, à la volée (au fur et à mesure de leur réétalement), en entrée des filtres adaptés. Les entrées des filtres adaptés recevront donc les signaux épurés des contributions passées à savoir n-Ik(t) où : Le premier terme de droite dans l'équation (33) représente la contribution des utilisateurs qui sont en avance sur l'utilisateur k et pour lesquels par conséquent une estimation des symboles ?*#(/) est déjà disponible. Le second terme représente la contribution des utilisateurs qui sont en retard sur l'utilisateur k: la contribution due aux symboles dk(i) ne peut bien entendu être évaluée à l'instant où le symbole dk(i) est désétalé. En revanche, la contribution due aux symboles précédents dk,(i-1) peut être évaluée puisque une estimation des symboles d/(i-1) est déjà disponible. Il importe de noter que, contrairement au cas synchrone où l'étape de détection par sphères estimait simultanément tous les symboles #k(i), dans le cas asynchrone, les symboles des différents utilisateurs sont estimés les uns après les autres suivant leur ordre d'arrivée. La Fig. 4 illustre un second mode de réalisation de l'invention mettant en oeuvre une élimination soustractive des interférences pour des utilisateurs asynchrones . Pour simplifier la représentation, seule la branche de détection concernant un utilisateur k a été illustrée. La somme Ik(t) symbolise la somme des contributions des autres utilisateurs et est soustraite en 405k au signal n en entrée du filtre adapté 410k. La détection par sphères pour l'utilisateur k démarre à l'instant où le symbole dk(i) est désétalé. Il faut noter qu'à cet instant les symboles interférents non encore estimés des autres utilisateurs n'auront fait l'objet que d'un désétalement partiel (sauf bien entendu si les symboles interférents sont synchrones). Le module 430k (resp. 431k) opère sur les composantes en sortie de 420k (resp. 421k). Alternativement les sorties de tous les modules 430k (resp. 431k) pourront être démultiplexées pour être traitées par un module commun 440 travaillant K fois plus rapidement. A la différence de la Fig. 3, ici seule la kième sortie des modules 430, et 431k est utilisée pour fournir le symbole estimé z(/). Le symbole dk(i) est ensuite réétalé et la contribution mdk(i)sk(t-iT -Tk) de l'utilisateur k est évaluée puis soustraite en entrée de chacun des filtres adaptés 41Ok' via / (. Si les canaux de transmission des différents utilisateurs sont du type multi-trajet, le problème de l'élimination des interférences est plus complexe puisqu'il faut considérer les interférences entre tous les trajets de tous les utilisateurs. Le signal reçu peut s'écrire : où Pk est le nombre de trajets du canal de transmission de l'utilisateur k, #pk=#k+#pk est le retard cumulé du retard (#k) à l'émission de l'utilisateur k et du retard (#pk) de propagation le long du trajet p du canal k et Cpk est le facteur multiplicatif complexe associé à ce trajet. On supposera à nouveau que 05z Les filtres adaptés 410k effectuent alors un filtrage adapté à la signature de l'utilisateur k et au canal de transmission k par combinaison MRC (Maximum Ratio Combining) des signaux des différents trajets. Plus précisément, le filtre 310k effectue l'opération suivante : Du fait de la dispersion des retards, la durée à prendre en compte pour l'interférence entre utilisateurs est de 2T et le symbole d'un utilisateur pourra interférer avec deux symboles consécutifs non encore estimés d'un autre utilisateur. La détection conjointe porte alors sur tous les symboles interférents non encore estimés, les symboles interférents estimés étant, quant à eux, exploités pour l'élimination soustractive. La complexité du détecteur par sphères est plus élevée que dans le cas asynchrone mono-trajet puisque la dimension du réseau à prendre en compte est de 2K-1au lieu de K. Enfin, il faut encore noter que l'estimation des contributions interférentes dans le cas asynchrone mono-trajet ou multi-trajet peut être réalisée soit à partir de symboles de constellation (décisions dures) soit à partir des symboles souples des différents utilisateurs. Bien que certains modes de réalisation de l'invention aient été représentés sous forme de modules fonctionnels, il est clair que le dispositif selon l'invention peut être réalisé sous forme d'un processeur programmé pour exécuter les différentes fonctions représentées ou sous la forme d'une pluralité de processeurs dédiés aptes à mettre en oeuvre une ou plusieurs de ces fonctions. REVENDICATIONS 1) Méthode de détection d'une pluralité de symboles (dk(i)) transmis par ou pour une pluralité K d'utilisateurs, chaque symbole appartenant à une constellation de modulation et faisant l'objet d'un étalement spectral par une séquence d'étalement, ladite méthode comprenant une étape de filtrage adapté (310I,..,310k) pour fournir un vecteur complexe (y(i),#(i)) caractéristique dudit signal reçu, caractérisée en ce que ledit vecteur complexe est décomposé en un premier vecteur (yR(I,yR(1)) et un second vecteur (y (/),y (/)) et qu'au moins les plus proches voisins des premier et second vecteurs sont recherchés (330,331) au sein d'un réseau de points (A,Q) généré par lesdites constellations de modulation, les symboles transmis étant estimés à partir des composantes desdits plus proches voisins. 2) Méthode de détection selon la revendication 1, caractérisée en ce que les séquences d'étalement (sk(t)) sont constitués de multiples réels (sk0(t)) d'un même coefficient complexe (#). 3) Méthode de détection selon la revendication 1 ou 2, caractérisée en ce que la recherche est limitée à un premier ensemble de points du réseau appartenant à une première zone prédéterminée (#R) autour du premier vecteur et à un second ensemble de points du réseau appartenant à une seconde zone prédéterminée (Si) autour du second vecteur. 4) Méthode de détection selon la revendication 1 ou 2, caractérisée en ce que la recherche est limitée à un premier ensemble de points du réseau appartenant à une première zone prédéterminée (#R) autour de l'origine et à un second ensemble de points du réseau appartenant à une seconde zone prédéterminée (SI) autour de l'origine. 5) Méthode de détection selon la revendication 3 ou 4, caractérisée en ce que lesdites première et seconde zones prédéterminées sont des sphères. 6) Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que la recherche du plus proche voisin du premier vecteur est effectuée sur une pluralité de composantes de ce dernier, la recherche étant limitée pour chacune desdites composantes à un intervalle défini par une borne inférieure et une borne supérieure, lesdites bornes étant choisies de manière à ce que ledit intervalle ne comprenne pas des points relatifs à des symboles ne pouvant appartenir à la constellation de modulation. 7) Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que la recherche du plus proche voisin du second vecteur est effectuée sur une pluralité de composantes de ce dernier, la recherche étant limitée pour chacune desdites composantes à un intervalle défini par une borne inférieure et une borne supérieure, lesdites bornes étant choisies de manière à ce que ledit intervalle ne comprenne pas des points relatifs à des symboles ne pouvant appartenir à la constellation de modulation. 8) Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que, préalablement à la recherche du plus proche voisin, le premier vecteur (yR(i)) est soumis à un traitement matriciel (320) visant à substantiellement décorréler les différentes composantes de bruit de ce dernier. 9) Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que, préalablement à la recherche du plus proche voisin, le second vecteur (yI(i)) est soumis à un traitement matriciel (321) visant à substantiellement décorréler les différentes composantes de bruit de ce dernier. 10) Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que ladite étape de recherche est étendue à celle d'un premier ensemble de points plus proches voisins dudit premier vecteur, dits premiers voisins, et à un second ensemble de points plus proches voisins dudit second vecteur, dits seconds voisins, et que les symboles transmis sont estimés de manière souple à partir des symboles générant lesdits premiers et second voisins et des distances séparant lesdits premiers voisins dudit premier vecteur d'une part et lesdits seconds voisins dudit second vecteur d'autre part. 11) Méthode de détection selon l'une des revendications précédentes, caractérisée en ce que l'on détermine à partir des symboles estimés les contributions de chaque utilisateur aux signaux obtenus par l'étape de filtrage adapté et que, pour un utilisateur k donné, l'on élimine en sortie de l'étape de filtrage adapté les contributions des autres utilisateurs correspondant aux symboles déjà estimés. 12) Méthode de détection selon l'une des revendications 1 à 10, caractérisée en ce que l'on détermine (340) les contributions de chaque utilisateur au signal reçu à partir des symboles estimés et que, pour un utilisateur k donné, l'on élimine en entrée de l'étape de filtrage adapté les contributions des autres utilisateurs correspondant aux symboles déjà estimés. 13) Méthode de détection selon l'une des revendications 1 à 10, caractérisée en ce que, les symboles desdits K utilisateurs étant transmis de manière synchrone, ledit réseau de points est de dimension K. 14) Méthode de détection de selon la revendication 11 ou 12, caractérisée en ce que, les symboles desdits K utilisateurs étant transmis de manière asynchrone et se propageant selon une pluralité de trajets, la dimension du réseau est égale au nombre de symboles des différents utilisateurs pouvant interférer et non encore estimés. 15) Dispositif de détection d'une pluralité de symboles (dk(/)) transmis par ou pour une pluralité K d'utilisateurs, chaque symbole appartenant à une constellation de modulation et faisant l'objet d'un étalement spectral par une séquence d'étalement, le dispositif comprenant des moyens pour mettre en oeuvre la méthode revendiquée selon l'une des revendications précédentes. 16) Récepteur pour système de télécommunication mobile DS-CDMA comprenant un dispositif de détection selon la revendication 15.