FR 2850501 A1 20040730 FR 0300983 A 20030129 PROCEDE ET DISPOSITIF DE RECEPTION D'UN SIGNAL RADIO La présente invention concerne les techniques de radiocommunication numérique utilisant un accès multiple à répartition par codes (CDMA, " CodeDivision Multiple Access "). Un signal CDMA traité par un récepteur a pour expression, après filtrage et transposition en bande de base: u y(t)= Z yu(t)+n(t) (1) u=1 o n(t) est un bruit additif et U est le nombre de canaux multiplexés sur la porteuse CDMA, dont les contributions yu(t) ont pour expression générale: yU (t) =bu.s ju(t-iT) (2) o: * bu est le symbole d'information de rang i transmis sur le u-ième canal; * sP(t) est un code généralisé donné par la convolution de la réponse impuisionnelle du u-ième canal avec la portion correspondant au symbole 15 bu du code d'étalement cu affecté au canal. Le nombre U correspond aux nombres d'utilisateurs si chaque utilisateur considéré dispose d'un seul canal. Il peut toutefois y avoir plusieurs canaux par utilisateur. Les codes d'étalement cu sont des séquences d'échantillons discrets 20 appelés " chips ", à valeurs réelles ( 1) ou complexes ( 1 j), ayant une cadence de chips donnée. Les symboles bu sont également à valeurs réelles (_ 1) ou complexes ( 1 j). La durée d'un symbole sur un canal est un multiple de la durée du chip, le rapport entre les deux étant le facteur d'étalement Q du canal. Dans certains systèmes, le facteur d'étalement peut varier d'un canal à un autre. Dans un tel cas, on considère un facteur d'étalement commun Q égal au plus grand commun diviseur (PGCD) des U facteurs d'étalement Qu. Un symbole sur le canal u est alors considéré comme la concaténation de QU/Q symboles consécutifs bu dont les valeurs sont identiques. La durée de la réponse généralisée sU(t) correspond à Q+W-1 chips 5 si W désigne la longueur de la réponse impulsionnelle exprimée en nombre de chips. En échantillonnant à la cadence des chips le signal CDMA y(t) reçu pour un bloc de n symboles sur chacun des canaux, le récepteur obtient des échantillons complexes qu'on peut modéliser par un vecteur Y de nxQ+W-1 10 composantes: Y=A.b+N (3) o: * b désigne un vecteur-colonne de taille nxU, qu'on peut décomposer en bT =(bT,b2T,...,bnT), (.)T représentant l'opération de transposition, les 1 2 n _ 2bu' vecteurs b. étant de taille U pour 1 En d'autres termes, les matrices Al s'écrivent: 2Ai = (çjQ" uç2.. O) (4) avec: u = MU.HU (5) I i I o Mu est une matrice de Toeplitz de taille (nxQ+W-1) x (nxQ+W-Q) obtenue à partir des valeurs cU(q) des chips du code d'étalement cu du u-ième canal pendant la durée du i-ième bit du bloc: c () O... 0 c (2) c (1) cU(2) . 0 MU= c (Q) . c(1) (6) 0 cU(Q) cU(2) 0... O cU(Q) et Hu est un vecteur-colonne de taille (n-1)xQ+W qui, lorsque les U canaux 5 sont reçus de façon synchronisée, contient (i-1)xQ composantes à zéro, suivies par les W échantillons de la réponse impulsionnelle du u-ième canal relative au i-ième symbole bu, suivies par (n-i)xQ autres composantes à zéro. Les décalages temporels de réception selon les différents canaux, en nombres de chips, se traduisent par des décalages correspondants des W échantillons 10 de la réponse impulsionnelle des canaux le long du vecteur HU. Le récepteur le plus couramment utilisé, dit récepteur en râteau ou "rake", utilise un ou plusieurs filtres adaptés pour estimer la valeur des symboles transmis sur chaque canal à partir d'une estimation de la réponse impulsionnelle du canal selon un ou plusieurs trajets de propagation. L'opération effectuée par de tels récepteurs revient à effectuer le produit matriciel: Z = *.Y (7) o * est la transposée conjuguée d'une estimation = ( 1, 2, ..., An) de la matrice A, les matrices Ai découlant des réponses impulsionnelles estimées en 20 appliquant les relations (4) et (5). Les nxU composantes Zu du vecteur Z sont des estimations souples respectives des nxU symboles bu du vecteur b. Si les décodages effectués en aval admettent des estimations souples en entrée, on peut utiliser directement les composantes du vecteur Z. Sinon, le signe de ces composantes est pris pour former les estimations dures des symboles. Le récepteur à filtre adapté est optimal lorsque les codes généralisés 5 (vecteurs nu) sont orthogonaux deux à deux, c'est-à-dire lorsque la matrice de corrélation R = A*.A est diagonale. En général, les systèmes adoptent des codes d'étalement orthogonaux deux à deux et ayant de bonnes propriétés d'autocorrélation, ce qui permet de vérifier cette condition en première approximation. Toutefois, lorsqu'on prend en compte la réponse impulsionnelle du canal, la condition d'orthogonalité n'est plus remplie. L'approximation cidessus devient mauvaise particulièrement en présence de trajets multiples de propagation. Certains récepteurs corrigent a posteriori les estimations souples des 15 symboles issues du récepteur à filtres adaptés en tenant compte de l'interférence entre les utilisateurs et/ou de l'interférence entre symboles sur un même canal, ce qui améliore sensiblement les performances. Une telle méthode, basée sur un algorithme dit MFPIC ("Matched Filter Parallel Interférence Cancellation"), est décrite dans WO 01/99301. Elle a pour 20 avantage de ne pas trop accroître la complexité globale des calculs par rapport au récepteur "rake" traditionnel. Cependant, ces corrections a posteriori font suite à une optimisation du système (3) reposant sur l'approximation ci-dessus. Cet algorithme MFPIC appartient à la classe des algorithmes de détection multi-utilisateurs (MUD, "Multi-User Detection") qui offrent de 25 meilleures performances que le simple récepteur à filtres adaptés. Il existe dans cette classe des algorithmes plus précis que le MFPIC, en ce qu'ils peuvent prendre en compte des termes de la matrice R plus éloignés de sa diagonale. Un bon exemple d'algorithme MUD utilisable dans un tel contexte est 30 l'algorithme SDP décrit dans l'article de M. Abdi, et ai., "Semidefinite Positive Relaxation of the Maximum-Likelihood Criterion Applied to Multiuser Detection in a CDMA Context", IEEE Signal Processing Letters, Vol. 9, No. 6, juin 2002, pp. 165-167. Les meilleures performances de ces algorithmes MUD sont obtenues au prix d'une complexité fortement augmentée. En général, leur complexité est 5 plus que linéaire en la taille nxU du problème à résoudre, ce qui les rend très coteux en comparaison avec des récepteurs plus classiques tels que le "rake" ou le MFPIC. Un but de la présente invention est de trouver un bon compromis entre les performances et la complexité d'un récepteur CDMA. L'invention propose ainsi un procédé de traitement d'un signal reçu par l'intermédiaire d'une interface radio, comportant des contributions de plusieurs canaux multiplexés par des codes d'étalement respectifs. Ce procédé comprend les étapes suivantes: la/ estimer des paramètres de réponse des canaux multiplexés; lb! calculer des estimations souples de symboles transmis sur les canaux multiplexés, en fonction du signal reçu et des paramètres de réponse estimés; /ci répartir les symboles dont des estimations souples viennent d'être calculées entre un premier ensemble de symboles vérifiant un critère 20 de confiance appliqué auxdites estimations souples et un second ensemble de symboles ne vérifiant pas le critère de confiance; Idl déterminer un signal modifié en retranchant du signal soumis au précédent calcul d'estimations souples des contributions estimées correspondant respectivement aux symboles du premier ensemble; et 25 le! calculer de nouvelles estimations souples des symboles du second ensemble seulement, en fonction du signal modifié et des paramètres de réponse estimés. Les symboles pour lesquels le premier calcul procure des estimations suffisamment fiables ne sont plus à estimer dans le second calcul. Leur 30 contribution au signal est retranchée pour prendre en compte l'interférence qu'ils engendrent pour les autres symboles. On peut employer divers algorithmes de détection pour les calculs - 6 successifs d'estimations souples. Dans une réalisation préférée du procédé, l'étape /e! est exécutée selon un algorithme de détection de nature plus complexe que l'étape /b/, notamment selon un algorithme de détection multi-utilisateurs. Le procédé 5 profite alors des performances d'un algorithme complexe de type MUD, mais en appliquant celui-ci à un problème de taille plus petite que celui faisant l'objet du premier calcul. Le procédé permet ainsi de réaliser un compromis entre la performance de la détection des symboles et la complexité des calculs requis. 10 Dans la conception ou la configuration du récepteur, on peut privilégier la performance ou la complexité dans ce compromis, par un réglage approprié du critère de confiance employé. Il est d'autre part possible de réitérer une fois ou plus la séquence des étapes Ici à /e!. Le nombre plus ou moins grand d'itérations et/ou la sévérité du 15 critère de confiance appliqué à chaque itération (ce critère de confiance peut varier d'une itération à la suivante) sont des paramètres qui peuvent aussi être réglés en fonction du compromis performances/complexité recherché. Un autre aspect de la présente invention se rapporte à un dispositif de traitement d'un signal reçu par l'intermédiaire d'une interface radio, comportant 20 des contributions de plusieurs canaux multiplexés par des codes d'étalement respectifs. Ce dispositif comprend des moyens d'estimation de paramètres de réponse des canaux multiplexés, des premiers moyens de calcul d'estimations souples de symboles transmis sur les canaux multiplexés, en fonction du signal reçu et des paramètres de réponse estimés, des moyens de répartition des 25 symboles entre un premier ensemble de symboles vérifiant un critère de confiance appliqué aux estimations souples produites par les premiers moyens de calcul et un second ensemble de symboles ne vérifiant pas le critère de confiance, des moyens de détermination d'un signal modifié en retranchant du signal reçu des contributions estimées correspondant respectivement aux 30 symboles du premier ensemble, et des seconds moyens de calcul de nouvelles estimations souples des symboles du second ensemble seulement, en fonction du signal modifié et des paramètres de réponse estimés. Un tel dispositif peut notamment être incorporé à une station de base d'un système de radiocommunication CDMA. Lorsque plusieurs itérations sont prévues, le dispositif comprend en outre des seconds moyens de répartition des symboles du second ensemble 5 entre un premier sous-ensemble de symboles vérifiant un second critère de confiance appliqué aux estimations souples produites par les seconds moyens de calcul et un second sous-ensemble de symboles ne vérifiant pas le second critère de confiance, des moyens de détermination d'un second signal modifié en retranchant du signal modifié qui a été soumis aux seconds moyens de 0 - calcul des contributions estimées correspondant respectivement aux symboles du premier sous-ensemble, et des troisièmes moyens de calcul de nouvelles estimations souples des symboles du second sous-ensemble seulement, en fonction du second signal modifié et des paramètres de réponse estimés. Un autre aspect de la présente invention se rapporte à un programme 15 d'ordinateur à installer dans un récepteur de radiocommunication, le programme comprenant des instructions pour mettre en oeuvre un procédé tel que défini ci-dessus lors d'une exécution du programme par une unité de traitement de signal du récepteur. D'autres particularités et avantages de la présente invention 20 apparaîtront dans la description ci-après d'exemples de réalisation non limitatifs, en référence aux dessins annexés, dans lesquels: - la figure 1 est un schéma synoptique d'un exemple de dispositif de réception selon l'invention; et - la figure 2 est un organigramme d'un exemple de procédé selon 25 l'invention. Le dispositif représenté sur la figure 1 fait partie de l'étage de réception d'une station de radiocommunication apte à communiquer avec plusieurs stations distantes 1. Les canaux montants utilisés par ces stations distantes 1 sont 30 multiplexés par la technique CDMA, de sorte que le signal radio capté par l'antenne 2, ramené en bande de base, peut être représenté sous la forme 8 (1)-(2) pour U canaux multiplexés en provenance de V stations (1 La station incorporant le dispositif est par exemple une station de base d'un système de radiocommunication cellulaire de troisième génération de type UMTS (" Universal Mobile Telecommunication System "). Sur la figure 1, l'unité 3 représente schématiquement les modules effectuant de façon classique les pré-traitements de réception du signal (amplification, filtrages, conversion en bande de base, échantillonnage à la fréquence des chips). Cette unité 3 délivre des blocs Y de nxQ+W-1 échantillons, correspondant à des blocs de n symboles émis simultanément sur 10 les U canaux. Si les blocs de n symboles se succèdent sans interruption sur les canaux, il y a un recouvrement de W échantillons (chips) entre les blocs Y successifs, correspondant à la durée de la réponse impulsionnelle. Les blocs de signal reçu Y sont fournis à un module 4 qui estime les réponses impulsionnelles H. des U canaux CDMA multiplexés, à l'aide de 15 corrélations avec les codes d'étalement produits par un générateur de codes pseudo-aléatoires 5 (1 Le premier terme soustrait dans (8) correspond à de l'interférence intersymboles sur le u-ième canal tandis que le second terme soustrait correspond à de l'interférence entre canaux. Les estimations sfiu produites par le module 8 sont des approximations des "softbits", à un coefficient multiplicatif 4/a2 près, a désignant la puissance du bruit additif capté, dont une estimation est 10 classiquement rendue disponible par le module de sondage 4. Le signe de sfiu forme une estimation dure du bit bu, tandis que sa valeur absolue mesure la vraisemblance de cette estimation. Selon l'invention, ces estimations souples sfiu sont examinées par un module 10 de manière à identifier un ensemble F de symboles dont les 15 estimations sont considérées comme les plus fiables. A titre d'exemple, le module 10 trie les estimations sfiu dans l'ordre des valeurs absolues décroissantes et place dans l'ensemble F les K % des symboles dont les estimations sont classées les premières. Une autre possibilité est de placer dans l'ensemble F les symboles dont les estimations 20 sfiu ont une valeur absolue supérieure à un seuil de confiance p qui peut être pris proportionnel à 4/a2. Le seuil p ou le pourcentage K peut être rendu adaptatif, par exemple en fonction du rapport signal-sur-bruit observé par le récepteur. Pour les symboles de l'ensemble F ainsi déterminé, le dsipositif 25 délivrera les estimations souples sfiu produites par le module 8 ou des 10 estimations dures bu constituées par leurs signes, selon les besoins des circuits de traitement situés en aval. Les symboles qui ne font pas partie de cet ensemble vont faire l'objet d'un second calcul d'estimations souples qui utilise de préférence un algorithme MUD. Les colonnes de la matrice des codes généralisés sont triées par un module 11. Chaque colonne iu associée à un symbole de l'ensemble F est extraite de pour être multipliée par l'estimation dure bu de ce symbole (multiplieur 12), ce qui donne une estimation de la contribution de l'interférence 10 causée par ce symbole sur les autres. Cette contribution bi.QPu est soustraite I I du vecteur Y représentant le signal d'entrée par le soustracteur 13. Une telle soustraction est effectuée pour chaque symbole de l'ensemble F: y, y - bu 2Y (9) (iu)EF n U = ZE b'u.W +N- E buhu (10) i=lu=1 (iu)eF ubhY.Q+N (11) (i,u)eF La relation (11) montre que la prise en compte de des contributions des symboles bien estimés réduit la taille du système à traiter de nxU à nxU - card(F). La réduction de taille est très appréciable si suffisamment de symboles ont pu être bien estimés lors de la première passe. Ceci facilite le 20 recours à un algorithme de nature plus complexe pour la seconde estimation des symboles restants, notamment à un algorithme MUD. Cette seconde estimation est effectuée par le module de détection 14 de la figure 1, auquel est soumis le signal modifié Y'. Chaque vecteur colonne avec (i,u) E F est enlevé de la matrice par le module 11 pour former une 25 matrice de codes généralisés réduite d'après laquelle le module 15 obtient la matrice de corrélation réduite R = . . Le vecteur Y' et les matrices réduites - 11 et R constituent les données d'entrée de l'algorithme MUD mis en oeuvre par le module 14, qui est par exemple l'algorithme SDP précédemment cité. Les estimations souples sf1u (ou dures) produites par le module MUD 14 sont finalement combinées aux estimations déterminées de manière 5 suffisamment fiable par l'algorithme MFPIC au cours de la première passe, pour reconstituer l'ensemble des données démodulées. Les estimations souples sfju produites par ce module MUD 14 peuvent aussi faire l'objet d'un tri entre estimations fiables et estimations peu fiables, ces dernières faisant alors l'objet d'un nouveau calcul d'estimations. Ce 10 processus de tri et de réestimation peut être réitéré un certain nombre de fois. Un tel procédure itérative est illustrée par la figure 2. Les deux premières étapes 20 et 21, préalables aux itérations, consistent en l'estimation classique des réponses des canaux et des codes généralisés Qu (opérations des modules 4 et 6 de la figure 1). A l'étape 22, 15 l'ensemble E des symboles à estimer au cours de la prochaine itération est initialisé au produit cartésien {1,2,...,n} x {1,2,...,U} correspondant à la totalité des symboles transmis au titre du bloc courant sur les canaux à traiter. Dans chaque itération, la première étape 23 consiste à assembler les vecteurs colonne Qu associés aux symboles de l'ensemble E pour former la 20 matrice de codes qui sera utilisée pour la détection suivante, et à calculer la matrice de corrélation R = A.A. La détection est opérée à l'étape suivante 24 sur la base du bloc Y à l'aide des matrices et R. A titre d'exemple, l'algorithme employé à l'étape 24 est le MFPIC dans la première itération et le SDP dans chaque itération suivante. Les estimations souples obtenues à l'étape 24 sont soumises au critère de confiance pour construire, à l'étape 25, l'ensemble F des index (i,u) des symboles les plus fiables, qui sont par exemple les index (i,u) de l'ensemble E tels que | sf1u | >p. - 12 Si toutes les estimations sont jugées suffisamment fiables (F = E lors du test 26), la procédure de détection est terminée et les estimations souples sf*u qui ont été calculées peuvent être délivrées à l'étape 27. Le seuil p peut varier au fur et à mesure des itérations. Il peut notamment décroître pour 5 rendre le critère de confiance de moins en moins sévère. Pour limiter la procédure à un nombre maximum X d'itérations, on peut fixer p = 0 pour la X-ième itération. Si certaines estimations sont insuffisamment fiables (F # E au test 26), les estimations dures bu des symboles de F sont déterminées à l'étape 28. 10 Les contributions estimées bhP de ces symboles sont soustraites du bloc Y à l'étape 29 (opération des modules 12 et 13 de la figure 1). Pour initialiser l'itération suivante, l'ensemble E des symboles à estimer est mis à jour à l'étape 30 en lui retirant les symboles correctement estimés de l'ensemble F. Une procédure comme celle illustrée par la figure 2 peut être mise en 15 oeuvre en programmant un processeur de traitement de signaux numériques prévu dans le récepteur radio. - 13 REVENDICATIONS 1. Procédé de traitement d'un signal (Y) reçu par l'intermédiaire d'une interface radio, comportant des contributions de plusieurs canaux multiplexés par des codes d'étalement respectifs, ce procédé comprenant les étapes suivantes: la/ estimer des paramètres de réponse des canaux multiplexés; /b! calculer des estimations souples de symboles transmis sur les canaux multiplexés, en fonction du signal reçu et des paramètres de réponse estimés; Ici répartir les symboles dont des estimations souples viennent d'être calculées entre un premier ensemble de symboles vérifiant un critère de confiance appliqué auxdites estimations souples et un second ensemble de symboles ne vérifiant pas le critère de confiance; /d/ déterminer un signal modifié en retranchant du signal soumis au 15 précédent calcul d'estimations souples des contributions estimées correspondant respectivement aux symboles du premier ensemble; et le! calculer de nouvelles estimations souples des symboles du second ensemble seulement, en fonction du signal modifié et des paramètres de réponse estimés. 2. Procédé selon la revendication 1, dans lequel l'étape le! est exécutée selon un algorithme de détection de nature plus complexe que l'étape lb!. 3. Procédé selon la revendication 1 ou 2, dans lequel l'étape /e! est exécutée selon un algorithme de détection multi-utilisateurs. 4. Procédé selon l'une quelconque des revendications précédentes, dans lequel on réitère au moins une fois la séquence des étapes Ici à lel. 5. Procédé selon la revendication 4, dans lequel le critère de confiance varie d'une itération de la séquence des étapes Ici à le! à la suivante. 14 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel le critère de confiance est exprimé comme une proportion des estimations souples les plus faibles. 7. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel le critère de confiance est exprimé comme un seuil de confiance à atteindre par les estimations souples des symboles. 8. Dispositif de traitement d'un signal (Y) reçu par l'intermédiaire d'une interface radio, comportant des contributions de plusieurs canaux multiplexés par des codes d'étalement respectifs, ce dispositif comprenant des moyens 10 (4-6) d'estimation de paramètres de réponse des canaux multiplexés, des premiers moyens (7-9) de calcul d'estimations souples de symboles transmis sur les canaux multiplexés, en fonction du signal reçu et des paramètres de réponse estimés, des moyens (10) de répartition des symboles entre un premier ensemble de symboles vérifiant un critère de confiance appliqué aux 15 estimations souples produites par les premiers moyens de calcul et un second ensemble de symboles ne vérifiant pas le critère de confiance, des moyens (12-13) de détermination d'un signal modifié (Y') en retranchant du signal reçu des contributions estimées correspondant respectivement aux symboles du premier ensemble, et des seconds moyens (14) de calcul de nouvelles 20 estimations souples des symboles du second ensemble seulement, en fonction du signal modifié et des paramètres de réponse estimés. 9. Dispositif selon la revendication 8, dans lequel les seconds moyens de calcul (14) sont agencés pour appliquer un algorithme de détection de nature plus complexe que les premiers moyens de calcul (7-9). 10. Dispositif selon la revendication 8 ou 9, dans lequel les seconds moyens de calcul (14) sont agencés pour appliquer un algorithme de détection multi-utilisateurs. il. Dispositif selon l'une quelconque des revendications 8 à 10, comprenant en outre des seconds moyens de répartition des symboles du 30 second ensemble entre un premier sous-ensemble de symboles vérifiant un 15 second critère de confiance appliqué aux estimations souples produites par les seconds moyens de calcul et un second sous-ensemble de symboles ne vérifiant pas le second critère de confiance, des moyens de détermination d'un second signal modifié en retranchant du signal modifié qui a été soumis aux 5 seconds moyens de calcul des contributions estimées correspondant respectivement aux symboles du premier sous- ensemble, et des troisièmes moyens de calcul de nouvelles estimations souples des symboles du second sous-ensemble seulement, en fonction du second signal modifié et des paramètres de réponse estimés. 12. Programme d'ordinateur à installer dans un récepteur de radiocommunication, le programme comprenant des instructions pour mettre en ceuvre un procédé selon l'une quelconque des revendications 1 à 7 lors d'une exécution du programme par une unité de traitement de signal du récepteur.