Brevet FR2822568A1 — publié le 20020927

La présente invention concerne une méthode et un dispositif d'égalisation pour voie descendante d'un système de télécommunication de type MC-CDMA. Plus particulièrement, la présente invention concerne une méthode et un dispositif d'égalisation de type GMMSE. L'accès multiple à répartition par codes sur porteuses multiples (MC-CDMA pour Multi-Carrier Code Division Multiple Access) combine la modulation OFDM (Orthogonal Frequency Division Multiplex) et la technique d'accès multiple CDMA. Cette technique d'accès multiple a été proposée pour la première fois par N. Yee et al. dans l'article intitulé Multicarrier CDMA in indoor wireless radio networks paru dans Proceedings of PIMRC '93, Vol. 1, pages 109-113,1993. Les développements de cette technique ont été passés en revue par S. Hara et al. dans l'article intitulé Overview of multicarrier CDMA publié dans IEEE Communication Magazine, pages 126-133, Décembre 1997. Contrairement à la méthode DS-CDMA (Direct Sequence Code Division Multiple Access) dans laquelle le signal de chaque utilisateur est multiplié dans le domaine temporel pour étaler son spectre fréquentiel, la signature multiplie ici le signal dans le domaine fréquentiel, chaque élément de la signature multipliant le signal d'une sous-porteuse différente. Plus précisément la Fig. 1 représente la structure d'un émetteur MC-CDMA pour un utilisateur k donné. Soit dk(i) le ième symbole à émettre de l'utilisateur k où dk(i) appartient à l'alphabet de modulation. Le symbole dk(i) est tout d'abord multiplié en 110 par une séquence d'étalement ou signature de l'utilisateur, notée ck(t) , constituée de N chips , chaque chip étant de durée Tc, la durée totale de la séquence d'étalement correspondant à une période symbole T. Les résultats de la multiplication du symbole dk(i) par les différents chips sont convertis par le convertisseur série/ parallèle 120 en un bloc de L symboles, où L est en général un multiple de N. Nous considérerons, pour des raisons de simplification de présentation, que L=N. Le bloc de L symboles est ensuite soumis à une transformée de Fourier rapide inverse (IFFT) dans le module 130 avant d'être transmis au convertisseur parallèle / série 140. Afin d'éviter l'interférence inter-symbole, un intervalle de garde, de longueur supérieure à la durée de la réponse impulsionnelle du canal de transmission, est adjoint au symbole MC-CDMA. Cet intervalle est obtenu par l'adjonction (non représentée) d'un suffixe choisi identique au début dudit symbole. Le symbole ainsi obtenu est amplifié en 150 avant d'être transmis sur le canal de l'utilisateur. On voit onc que la méthode MC-CDMA s'analyse en un étalement dans le domaine spectral (avant IFFT) suivi d'une modulation OFDM. En pratique, une station de base transmet les données pour un utilisateur k sous forme de trames de I symboles, chaque symbole dk(i) étant étalé par une signature réelle Ck(t) de durée égale à la période symbole T , telle que c)=Osi[0,T[. Le signal modulé à l'instant t=i.T+.T, peut alors s'écrire, si l'on omet les intervalles de garde entre symboles MC-CDMA : où vk est l'amplitude du signal transmis pour l'utilisateur k, supposée constante pour un bloc de transmission. Un récepteur MC-CDMA pour un utilisateur k donné a été illustré schématiquement en Fig. 2. Le signal reçu démodulé est échantillonné à la fréquence chip et les échantillons appartenant à l'intervalle de garde sont supprimés (suppression non représentée). Le signal obtenu peut s'écrire : où K est le nombre d'utilisateurs, cke=ck(.T) , et où hkp(i) traduit la réponse du canal de l'utilisateur k à la fréquence de la sous-porteuse du symbole MC-CDMA émis au temps i. T et où b(t) est le bruit reçu. Si l'on considère la voie descendante, les canaux de transmission ont des caractéristiques identiques et l'on peut écrire hk#=h#. On se limitera par la suite à l'étude de la voie descendante. Les échantillons obtenus par échantillonnage à la fréquence chip sont parallélisés dans un convertisseur série/ parallèle 210 avant de subir une FFT dans le module 220. Les échantillons dans le domaine fréquentiel, en sortie de 220, sont égalisés et désétalés par la signature de l'utilisateur k . Pour ce faire, les échantillons du domaine fréquentiel sont multipliés par les coefficients qe,c:e dans les multiplieurs 2300,.., 230L-),puis sommés en 240 pour fournir un symbole de sortie dk(i) . Différentes possibilités d'égalisation ont été envisagées dans l'état de la technique : - la combinaison MRC (Maximum Ratio Combining) définie par l'emploi des coefficients q#=h#* où .* est l'opération de conjugaison; - la combinaison EGC (Equal Gain Combining) définie par l'emploi des coefficients qf=e-JfPt où h(=pt.eJ9' ; - le forçage à zéro ZF (Zéro Forcing) où qf=hi1(i); - le forçage à zéro ZF avec seuil (Th) où ='(ï)si (/)) (Minimum Mean Square Error) sur chacune des porteuses : qe= Iifl he + 2 où #h##2+# #2est la variance du bruit sur la porteuse. On rappelle qu'en MC-CDMA , la présence d'une période de garde permet de négliger l'interférence inter-symbole. L'égalisation peut être donc simplement réalisée au moyen d'une simple multiplication par un coefficient complexe, porteuse par porteuse. Le récepteur illustré en Fig. 2 décode les données d'un utilisateur k sans prendre en compte l'interférence due aux autres utilisateurs. Il est dénommé pour cette raison mono-utilisateur ou SU (Single User). Les techniques de détection multi-utilisateur sont connues, notamment dans les systèmes de télécommunication CDMA. Elles ont pour caractéristique commune de prendre en compte les interférences générées par les autres utilisateurs. Une technique de détection multi-utilisateur ou MUD (Multi User Detection) pour MC-CDMA a été présentée dans l'article de J-Y. Beaudais, J. F. Hélard et J. Citerne intitulé A novel linear MMSE détection technique for MC-CDMA publié dans Electronics Letters, Vol. 36, N 7, pages 665-666, le 30 Mars 2000. La méthode d'égalisation proposée n'opère plus porteuse à porteuse mais symbole MC-CDMA par symbole MC-CDMA, en prenant en compte l'ensemble des porteuses. Elle est appelée pour cette raison égalisation GMMSE (Global Minimum Mean Square Error) ou, de manière équivalente, égalisation M-MMSE (Matrix Minimum Mean Square Error). Elle a pour but de minimiser l'erreur quadratique moyenne entre les symboles estimés dk(i) et les symboles transmis dk(i). Un récepteur MC-CDMA à égalisation GMMSE pour un utilisateur k (encore appelée dans ce cas per user MMSE) a été représenté en Fig. 3. Il diffère de celui de la Fig. 2 en ce que l'égalisation est réalisée au moyen d'une multiplication 331 par une matrice Q des signaux des différentes porteuses. Après désétalement par multiplication en 3320,...,332L-1 avec la signature de l'utilisateur k et sommation en 340, le signal obtenu est multiplié en 360 par le niveau vk d'émission relatif à l'utilisateur k pour fournir un symbole estimé #k (nous omettrons, par souci de simplification, l'indice temporel i). Un récepteur MC-CDMA à égalisation GMMSE, estimant les symboles dk transmis à destination des différents utilisateurs, a été représenté en Fig. 4. Les modules 410,420 et 431 sont identiques aux modules 310,320 et 331. Le désétalement est effectué en 432, au moyen de la multiplication par la matrice CH où C est la matrice des N codes possibles (.H désigne le transposé conjugué). Plus précisément, la matrice C est la matrice dont les colonnes sont constituées par les N codes possibles ck, k=0,..,N 1 , chaque code ck étant un vecteur dont les composantes forment la séquence c#k.Le vecteur en sortie de 432 est soumis à une multiplication en 433 par la matrice V" où V=Diag(vo, .... vN-I) est la matrice diagonale représentant les amplitudes vk des signaux émis par la station de base à destination des N utilisateurs possibles. Lorsque la station de base sert moins de N utilisateurs, autrement dit lorsque le système ne travaille pas à pleine charge, les niveaux d'émission correspondants aux codes non utilisés sont nuls. Le signal MC-CDMA en sortie de l'étape de FFT, 420, peut s'écrire sous la forme d'un vecteur Y de dimension N dont les composantes correspondent aux différentes porteuses. En reprenant les notations ci-dessus, on peut exprimer Y en fonction du vecteur D des symboles dk transmis à destination des différents utilisateurs : Y=HCVD+b=AD+b (3) où A=HCV et où b est le vecteur des composantes de bruit. La matrice diagonale H représente la réponse fréquentielle du canal aux différentes sous-porteuses, soit H=Diag (ho, hl, 1 hN-,). Le vecteur Y est un observable à partir duquel on peut estimer le vecteur D par filtrage de Wiener. Si D est le vecteur des symboles estimés pour les différents utilisateurs d'après l'équation de Wiener-Hopf : f)=RDYR-' Y=(RDDAHARDDAH+ab 2.1 N'y (4) où RDY, Ryy, RDD sont respectivement la matrice de covariance de Y et D et les matrices d'autocovariance de Y et de D. Si l'on suppose que les symboles transmis dk sont décorrélés, on a RDD=IN et donc : Î)=AH(AA il +O"b 2.1 N 1 y (5) Puisque AH=VHCHHH on déduit de (5) que la matrice d'égalisation Q peut se mettre sous la forme : Q=H"(AA"+6h.IN)-' (6) Comme l'indique (6), l'opération d'égalisation passe par le calcul de la matrice AAH et par l'inversion de la matrice AAH+#b2.IN (ou de manière équivalente la résolution d'un système linéaire de rang N). Ces calculs consomment des ressources du terminal mobile et, ce, d'autant qu'ils doivent être effectués fréquemment, notamment en cas de variations rapides du canal de transmission. Or, l'on sait que les ressources des terminaux constituent précisément la partie critique d'un système de télécommunication mobile. Le problème à la base de la présente invention est de proposer une méthode d'égalisation GMMSE pour voie descendante d'un système de télécommunication MC-CDMA qui ne mobilise pas d'importantes ressources calculatoires. Ce problème est résolu par l'invention définie comme un méthode d'égalisation pour voie descendante d'un système de télécommunication transmettant des symboles MC-CDMA sur une pluralité de fréquences porteuses, dans laquelle une pluralité de liaisons est établie entre un émetteur et une pluralité de récepteurs pour transmettre une pluralité de signaux codés, chaque liaison utilisant un code d'accès distinct parmi N codes possibles, et dans laquelle, pour au moins récepteur, un vecteur (Y) représentant les composantes du signal reçu par le récepteur sur les différentes fréquences porteuses est soumis à une étape de filtrage adapté pour fournir un vecteur observable (Z, Z') et que ledit vecteur observable est utilisé pour effectuer une estimation desdits symboles transmis selon un critère de minimisation de l'erreur quadratique moyenne. Avantageusement, ladite estimation comprend le calcul de la matrice AHA+#b2IN avec A=HCV où H est une matrice donnant les caractéristiques du canal de transmission aux différentes fréquences porteuses, C est la matrice des codes possibles, V est une matrice diagonale donnant les niveaux d'émission des signaux codés au moyen desdits N codes possibles, 66 est la variance de bruit et IN est la matrice unité de taille NxN. Ladite estimation pourra comprend une étape de résolution d'un système d'équations linéaires fllâ--Z où 2=AhA+ Selon un mode de réalisation, on calcule la matrice B=CH#H#2C où C est la matrice des codes possibles et où H est une matrice donnant les caractéristiques du canal de transmission aux différentes fréquences porteuses. On pourra sélectionner à partir des niveaux d'émission, réels ou estimés, des signaux codés au moyen desdits N codes possibles, un sous-ensemble de codes parmi les N codes possibles. Dans ce cas, la matrice B pourra être avantageusement réduite à une matrice B' de taille PxP en ne retenant que les P lignes et P colonnes correspondant aux codes sélectionnés et la matrice C des codes possibles à une matrice C' de taille NxP ne en retenant que les P colonnes correspondant aux codes sélectionnés. La matrice sera alors calculée sous la forme B'=C"" H C' et où H est une matrice donnant les caractéristiques du canal de transmission aux différentes fréquences porteuses. Si P codes parmi sont sélectionnés, la matrice V=Diag(vo,..,vN~) ,où les v, sont les niveaux réels ou estimés, pourra être réduite à une matrice V' de taille PxP ne comportant plus que les niveaux associés aux dits codes sélectionnés, ladite estimation comprenant une étape de calcul de la matrice Q'=V'HB ,V,+ ab 2.1@ où ai est la variance du bruit sur chaque porteuse et Ip la matrice unité de taille PxP. Dans ce cas, ladite estimation pourra comprendre une étape de résolution d'un système d'équations linéaires SZ'D'=Z' où Z' est la réduction du vecteur observable Z aux P composantes correspondant aux codes sélectionnés et où D'est le vecteur des symboles estimés sur les différentes liaisons utilisant lesdits codes sélectionnés. L'invention est également définie par un dispositif d'égalisation pour un récepteur de système de télécommunication de type MC-CDMA, comprenant des moyens pour mettre en oeuvre la méthode d'égalisation telle qu'exposée ci-dessus. L'invention concerne enfin un récepteur pour système de télécommunication de type MC-CDMA comprenant un tel dispositif d'égalisation. 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 de manière schématique la structure d'un émetteur MCCDMA connu de l'état de la technique ; La Fig. 2 représente de manière schématique la structure d'un récepteur MCCDMA mono-utilisateur connu de l'état de la technique; La Fig. 3 représente de manière schématique la structure d'un récepteur MCCDMA utilisant une égalisation per user MMSE, connu de l'état de la technique; La Fig. 4 représente de manière schématique la structure d'un récepteur MCCDMA utilisant une égalisation GMMSE, connu de l'état de la technique; La Fig. 5 représente de manière schématique la structure d'un récepteur MCCDMA utilisant une méthode d'égalisation selon un premier mode de réalisation de l'invention ; La Fig. 6 représente de manière schématique la structure d'un récepteur MCCDMA utilisant une méthode d'égalisation selon un second mode de réalisation de l'invention. Un première idée pour simplifier l'opération d'égalisation est de tirer profit du fait que le système ne travaille généralement pas à pleine charge. Lorsque le nombre P d'utilisateurs effectivement présents est inférieur au nombre N de codes possibles, on construit une matrice V' de taille PxP obtenue à partir de la matrice V en éliminant les éléments relatifs aux codes non utilisés. On peut supposer, sans perte de généralité que ces codes portent les indices 0,..P-1, d'où V'=Diag(vo,,,,vp-I)' La matrice C peut être également réduite à une matrice C' de taille NxP ne comportant que les colonnes qui correspondent aux codes utilisés. La matrice A, produit des matrices H, C et V se réduit alors également à une matrice A' de taille NxP : A'=HC'V' et le vecteur Y s'écrit : Y=A'D'+b où D' est le vecteur D réduit aux composantes correspondants aux P utilisateurs présents. L'équation (5) peut alors se mettre sous la forme : D'=AtH(A'A'H+U2 b l'y (7) Cependant, l'équation (7) n'apporte malheureusement pas de simplification substantielle des calculs car elle fait intervenir la matrice A'A'H qui est encore de taille NxN bien que A' soit de taille réduite NxP. La taille de la matrice à inverser (ou de manière équivalente le rang du système d'équations linéaires à résoudre) reste donc fonction du nombre N de codes possibles, c'est-à-dire encore du nombre maximal d'utilisateurs que peut supporter le système. Nous avons supposé jusqu'ici qu'un utilisateur n'employait qu'un seul code. De manière plus générale, nous ferons l'hypothèse qu'une pluralité de liaisons peuvent être établies entre l'émetteur et un utilisateur, chaque liaison employant un code distinct. Selon cette définition, un canal de transmission entre l'émetteur et un utilisateur peut supporter plusieurs liaisons. Ainsi, N devient le nombre maximal de liaisons que peut autoriser le système et P est le nombre de liaisons effectivement établies. L'idée générale à la base de l'invention est de trouver une formulation de l'opération d'égalisation susceptible de simplification lorsque le nombre P de liaisons (de codes utilisés) est inférieur au nombre maximal N. Considérons le vecteur Y en sortie de la FFT, on a vu que Y=AD+b où D est le vecteur des symboles à estimer. Soit Z=AHY, le résultat du filtrage adapté de Y : sait que Z est un résumé exhaustif de Y. On peut prendre comme observable le vecteur Z au lieu de Y et estimer le vecteur D à partir de cet observable au moyen d'un filtrage de Wiener. Le vecteur estimé D s'écrit donc : #=RDZRZZ-1Z (8) Or Roz=E(DyH A}=(RooA H +ROb}A=A H A et de même R =AHARDDAHA+AHRbbA=(AAH+c762IN)AHA d'où D=(AHA+66IN)-'Z=(AHA+ablN)-lA"Y (9) Lorsque le système ne travaille pas à pleine charge, il est possible de réduire la matrice A à une matrice A' et de réécrire l'équation (9) sous la forme réduite : D=(A,a A,+6b Ip)-A,H Y (10) où Ip est la matrice unité de taille PxP. On remarque que cette nouvelle formule d'estimation ne nécessite plus que l'inversion d'une matrice de taille PxP (ou de manière équivalente la résolution d'un système linéaire de P équations à P inconnues) au lieu d'une matrice de taille NxN comme précédemment. L'estimation des symboles s'en trouve considérablement simplifiée et, ce, d'autant que la charge du système est plus faible. Nous allons maintenant montrer que, selon un mode de réalisation de l'invention, le calcul de la matrice A'HA' figurant dans l'équation (10) peut être davantage simplifié lorsque la matrice des codes C possède certaines propriétés. En effet, considérons tout d'abord la matrice AHA . Celle-ci peut encore s'exprimer: AH A=VH(CHHHHC)V=VH(CHIHI2C)V=VHBV (11) Ou B=C''III ZC . Supposons tout d'abord que le code employé soit un code de Walsh-Hadamard (WH). On rappelle que le code WH peut être généré comme les colonnes d'une matrice C2p obtenue par récurrence: La matrice est réelle, symétrique, et égale à sa propre inverse. De plus, les séquences sont orthogonales. Les éléments de C peuvent s'exprimer sous la forme : c!I=(-l)M où est le produit scalaire des représentations binaires de / et/ : (i,j)= X jk k On en déduit la propriété suivante : C;Ckj=Ck,lEJ)j (13) où i#j représente l'addition bit à bit des entiers i et/. Eneffet : 'J}nkn+ )nkn= I(in+ jn)kn= 'JJiEB j)nkn et donc k,iy+k,=k,i j) . n n n n Il résulte de (13) que le produit conjugué de deux codes WH est encore un code WH. Supposons maintenant que les codes utilisés soient des codes de Fourier, définis par : cr,k = exp(-j 2nik/N). La matrice du code est symétrique. C ki Ckl =Cki-i (14) où "j-i" désigne ici la soustraction modulo N. Il en résulte là aussi que le produit conjugué de deux codes de Fourier est encore un code de Fourier. Cette propriété peut être généralisée à une famille de codes obtenus par rotation des codes précédents, encore appelés codes tournés, tels que présentés dans l'article de A. Bury et al. intitulé Diversity comparison of spreading transforms for multi-carrier spread spectrum transmission methods publié dans Electronic Wireless 2000 Proceedings, Dresden Germany. La matrice d'un code tourné est définie par : Cr = C. D (r) où D est une matrice diagonale définie par le vecteur r(ro,..,rN 1)T soit D=Diag(ro,..,m-I) où r,=exp(i(p,). Il est important de noter que coder N symboles à destination des N utilisateurs au moyen de la matrice d'un code tourné est équivalent à effectuer une rotation de chacun de ces symboles par un angle puis à appliquer la matrice du code initial. Autrement dit, un code tourné d'un code initial peut se ramener au code initial au prix d'une rotation de la constellation pour chacun des utilisateurs. Si la matrice C du code initial est relative à un code de Fourier, alors pour le code de Fourier tourné : cf/kJ=r'4rrJ.ckjxp( avec Ày=ffJj -ffJ,-ffJj-1 De même, pour un code de Hadamard tourné : Ck,h.j=r; ,ck".rj'ck,j =exp( ffJj -ffJ,),ck,j(f),=Àyek,j(f)1 (17) avec Ày = Dans ce cas les coefficients bij de la matrice B=CH#H#2C peuvent s'exprimer par : On calcule tout d'abord, par exemple, la première ligne de la matrice B, c'est-àdire les coefficients b0,J, puis on en déduit les autres lignes par permutation et multiplication : Si l'on dispose d'une table de f0-1f, et, le cas échéant de #iJ/#-,J, le calcul de la #0, J matrice B ne nécessite que le calcul d'une seule de ses lignes. A titre d'exemple, l'expression (20) s'écrit dans le cas d'un code WH : biJ=b0,l#J. Si l'on utilise par exemple un code WH, un code WH tourné ou encore un code de Fourier, un code de Fourier tourné, la première ligne de la matrice B peut être avantageusement obtenue par simple transformée du vecteur p=( ho Z, .., hN 2)T. En effet, dans le cas d'un code WH, la première ligne est obtenue par transformée de Hadamard de p: et dans le cas d'un code de Fourier, la première ligne est obtenue par transformée de Fourier de p: Avantageusement, la transformée de Fourier sera effectuée au moyen d'une FFT. Si des versions tournées de ces codes sont utilisées, les coefficients de la première ligne de B se déduisent de (21) et (22) au moyen d'une simple multiplication par un coefficient complexe. De manière générale, si le code employé vérifie la propriété (18), notamment pour un code de Hadamard, Hadamard tourné, Fourier ou Fourier tourné, on calculera avantageusement une ligne de la matrice B (ou une colonne puisqu'elle est hermitienne) et on en déduira les autres lignes (respectivement, les autre colonnes) au moyen d'opérations de permutation et le cas échéant de multiplication. Considérons maintenant la matrice A'HA' figurant dans le calcul d'estimation (10). La matrice A' peut encore s'écrire A'=(HC)'V' où (HC) ' est la réduction de la matrice HC à ses P premières colonnes. On en déduit une première expression de A'HA' : A'" A'=V'"(C"IHI2C)'V'=V'"B'V' (23) Selon cette variante, on effectue d'abord de calcul de la matrice B=CH#H#2C en utilisant la méthode simplifiée, exposée plus haut, puis l'on procède à la réduction de dimension. Alternativement, la matrice A'HA' peut être obtenue directement à partir de A'=HC'V' . Une seconde expression de A'HA' s'en déduit: A,H A'=V'"(CI"IHI2C')V' (24) Selon cette variante, on réduit d'abord la matrice des codes C, puis on calcule la matrice C'H#H#2C'. Lorsque le nombre de codes P est faible, on optera pour le calcul de A'HA' selon la formule (24). Ce calcul conduit à environ (NP2+2P2+NP) opérations. Si, en revanche, la réduction du nombre de codes n'est pas très importante, on optera pour le calcul selon la formule (23), dès lors qu'il existe un algorithme rapide pour l'application de la transformée C. Ce calcul est d'autant plus simple que, dans la plupart des cas, les codes de la matrice C vérifient la propriété (18). Si l'on utilise par exemple un code WH, la première ligne de la matrice B est obtenue par une transformée de Hadamard du vecteur p=(ho2,..,hN j2)T nécessitant Nlog2(N) opérations. Le calcul complet de A'HA' requiert (Nlog2(N)+2p2) opérations. Cette opération sera donc plus efficace dès que log2(N) sera plus petit que P2. Cette conclusion s'applique également lorsque le code utilisé est un code de Fourier. Avantageusement, on pourra passer de manière dynamique du calcul de A'HA' selon (23) au calcul selon (24) et réciproquement, suivant la valeur de P. Nous avons supposé jusqu'à présent que la réduction des matrices consistait à se limiter aux composantes relatives aux codes utilisés. En fait, une réduction plus poussée peut être obtenue sur la base de l'un des critère suivants : (I) on sélectionne les P-1(P entier donné tel que P Th où Th est un seuil prédéterminé. Avantageusement, le seuil Th sera adaptatif et dépendra du niveau d'émission pour l'utilisateur considéré. On sélectionnera alors les codes cm tels que : v. >,u. vk avec, par exemple, 0 liaisons), on prendra en compte le code d'indice k,5e{kl,..,kà } associé au niveau vks le plus faible et l'on sélectionnera les codes cm tels que : vm> . vk# Bien entendu une combinaison des critères précédents est également envisageable. En outre, en complément à la sélection à partir des niveaux d'émission, on peut sélectionner parmi ceux retenus selon ce critère, seuls les codes les plus intrinsèquement brouilleurs (c'est-à-dire hormis toute considération de niveau d'émission) pour le code considéré. En effet, on sait, par exemple de l'article de D. Mottier et D. Castelain intitulé 'A spreading séquence allocation procédure for MCCDMA tansmission systems' publié dans VTC' Fall 2000, Sept. 2000 inclus ici par référence, que, si les coefficients hk du canal sont corrélés, pour un code d'utilisateur donné, certains codes, à niveau d'émission égal, sont plus brouilleurs que d'autres. Plus précisément, l'interférence entre deux codes c, et cJ est liée aux propriétés de la séquence produit Wy, définie par wy{k]=cJk] c/fk]. Nous appellerons ce critère de sélection le critère de brouillage intrinsèque. La Fig. 5 illustre un égaliseur GMMSE selon un premier mode de réalisation de l' invention. Les modules 510,520 sont identiques aux modules 410 et 420 de la Fig. 4. Le vecteur v=(vo, ...,vN-1)T des niveaux d'émission, reçus ou estimés par le récepteur de l'utilisateur k considéré, est soumis en 570 à un module de sélection mettant en oeuvre un critère de sélection selon (I) ou (II), éventuellement en combinaison avec le critère de brouillage intrinsèque. Le module de sélection 570 détermine un ensemble Fk #{0,..,N-1} d'indices de codes à sélectionner et le transmet aux multiplexeurs 551 , 561 et 571. Le multiplexeur 571 sélectionne les composantes de v relatives aux indices de Fk pour fournir le vecteur de taille réduite v'. De même, les multiplexeurs 551 et 561 fournissent respectivement les matrices de taille réduite B' et C. Le module 550 est un filtre adapté fournissant Z'=A'"Y . Notons qu'en toute rigueur l'expression filtrage adapté ne convient qu'au cas où la sélection retient tous les indices des codes utilisés. Nous l'étendrons au cas où les codes faiblement interférents ne sont pas sélectionnés. Le module de calcul matriciel 590 calcule ensuite la matrice fl~~V,IIB,V,+ U. 2.1 et la transmet au module d'estimation 580. Ce dernier résout le système linéaire de P équations à P inconnues dk : SL'D'=Z' . La Fig. 6 illustre un égaliseur GMMSE selon un second mode de réalisation de l'invention. Les modules 610,620, 650,651, 661,670, 671,680, 690 sont respectivement identiques aux modules 510,520, 550, 551, 561,570, 571,580, 590 de la Fig. 5. Comme précédemment, le module 650 effectue un filtrage adapté sur vecteur Y et fourni le vecteur Z'=A'"Y. L'égaliseur comprend ici deux branches possibles pour le calcul de B' selon la valeur de P : Si P2>log2N, le calcul est effectué dans la branche supérieure comprenant les modules 660 et 661. Le vecteur p=(lhoI2,..,lhN-II2)T issu d'un estimateur de canal est soumis en 660 à une transformée (Hadamard, Fourier) correspondant au code utilisé pour fournir la première ligne de la matrice B. Puis, à l'aide d'une table donnant la permutation des indicesf0-1f, le module 660 génère au moyen de (20) les autres lignes de la matrice B. Le multiplexeur 661 réduit ensuite la matrice B à une matrice B' de taille PxP en fonction de l'ensemble Fk des indices à sélectionner, fourni par le module de sélection 670. En revanche, si P2#log2N, le calcul est effectué dans la branche inférieure. Le module 662 reçoit du module 651 la matrice réduite C' et calcule directement la matrice B' selon B'=C'H#H#2C'. Le module 690 de calcul matriciel reçoit le vecteur de taille réduite v' du multiplexeur 671 et la matrice B' de la branche supérieure ou de la branche inférieure, selon le cas. Il calcule la matrice fl=v",B'V'+ Crb '.1 et la transmet au module 680. Ce dernier résout alors le système linéaire de P équations à P inconnues dk : #'#'=Z'. Bien que le dispositif de réception MC-CDMA utilisant la méthode d'égalisation selon l'invention ait été représenté sous la forme de modules fonctionnels, il va de soi que tout ou partie de ce dispositif peut être réalisé au moyen d'un processeur unique dédié ou programmé pour remplir l'ensemble des fonctions décrites ou encore sous la forme d'une pluralité de processeurs dédiés ou programmés pour remplir chacun une partie desdites fonctions. REVENDICATIONS 1) Méthode d'égalisation pour voie descendante d'un système de télécommunication transmettant des symboles MC-CDMA sur une pluralité de fréquences porteuses, dans laquelle une pluralité de liaisons est établie entre un émetteur et une pluralité de récepteurs pour transmettre une pluralité de signaux codés, chaque liaison utilisant un code d'accès distinct parmi N codes possibles la méthode étant caractérisée en ce que, pour au moins récepteur, un vecteur (Y) représentant les composantes du signal reçu par le récepteur sur les différentes fréquences porteuses est soumis à une étape de filtrage adapté pour fournir un vecteur observable (Z, Z') et que ledit vecteur observable est utilisé pour effectuer une estimation desdits symboles transmis selon un critère de minimisation de l'erreur quadratique moyenne. 2) Méthode d'égalisation selon la revendication 1, caractérisée en ce que ladite estimation comprend le calcul de la matrice AHA+#b2IN avec A=HCV où H est une matrice donnant les caractéristiques du canal de transmission aux différentes fréquences porteuses, C est la matrice des codes possibles, V est une matrice diagonale donnant les niveaux d'émission des signaux codés au moyen desdits N codes possibles, #b2 est la variance de bruit et IN est la matrice unité de taille NxN. 3) Méthode d'égalisation selon la revendication 2, caractérisée en ce que ladite estimation comprend une étape de résolution d'un système d'équations linéaires 52D=Z où SZ=A"A+6bIN , Z est ledit vecteur observable et D est le vecteur des symboles estimés pour les différentes liaisons. 4) Méthode d'égalisation selon la revendication 1, caractérisée en ce que l'on calcule la matrice B=CH#H#2C où C est la matrice des codes possibles et où H est une matrice donnant les caractéristiques du canal de transmission aux différentes fréquences porteuses. 5) Méthode d'égalisation selon la revendication 1, caractérisée en ce que l'on sélectionne à partir des niveaux d'émission, réels ou estimés, des signaux codés au moyen desdits N codes possibles, un sous-ensemble de codes parmi les N codes possibles. 6) Méthode d'égalisation selon la revendication 4 et 5, caractérisée en ce que l'on réduit la matrice B à une matrice B' de taille PxP en ne retenant que les P lignes et P colonnes correspondant aux codes sélectionnés. 7) Méthode d'égalisation selon la revendication 5, caractérisée en ce que l'on réduit la matrice C des codes possibles à une matrice C' de taille NxP ne en retenant que les P colonnes correspondant aux codes sélectionnés. 8) Méthode d'égalisation selon la revendication 7, caractérisée en ce que l'on calcule la matrice B'=C'H#H#2C' et où H est une matrice donnant les caractéristiques du canal de transmission aux différentes fréquences porteuses. 9) Méthode d'égalisation selon la revendication 6 ou 8, caractérisée en ce que P codes parmi N ayant été sélectionnés, la matrice V=Diag(v,,..,v,-,) ,où les v, sont les niveaux réels ou estimés, est réduite à une matrice V' de taille PxP ne comportant plus que les niveaux associés aux dits codes sélectionnés, ladite estimation comprenant une étape de calcul de la matrice Q!=VnB'Y+al.IP où ai est la variance du bruit sur chaque porteuse et Ip la matrice unité de taille PxP. 10) Méthode d'égalisation selon la revendication 9, caractérisée en ce que ladite estimation comprend une étape de résolution d'un système d'équations linéaires #'#'=Z' où Z' est la réduction du vecteur observable Z aux P composantes correspondant aux codes sélectionnés et où D' est le vecteur des symboles estimés sur les différentes liaisons utilisant lesdits codes sélectionnés. 11) Dispositif d'égalisation pour un récepteur de système de télécommunication de type MC-CDMA, caractérisé en ce qu'il comprend des moyens pour mettre en oeuvre la méthode d'égalisation selon l'une des revendications précédentes. 12) Récepteur pour système de télécommunication de type MC-CDMA, caractérisé en ce qu'il comprend un dispositif d'égalisation selon la revendication 11.