Brevet FR2783119A1 — publié le 20000310

PROCEDE D'EGALISATION NUMERIQUE, ET RECEPTEUR DE RADIOCOMMUNICATION METTANT EN EUVRE UN TEL PROCEDE La présente invention concerne l'égalisation numérique des signaux. Elle trouve une application importante dans le domaine des radiocommunications. Le procédé s'applique lorsqu'on reçoit un signal issu d'un émetteur par l'intermédiaire d'un canal de transmission entre émetteur et récepteur, dont la réponse est connue ou a été préalablement estimée. Un problème principal qui se pose alors est celui du compromis entre les performances de l'égaliseur et sa complexité. Une estimation complète, selon le maximum de vraisemblance, de tous les symboles discrets composant le signal émis est possible, par exemple en employant l'algorithme de Viterbi (voir G.D. Forney Jr.: " The Viterbi Algorithm ", Proc. of the IEEE, Vol. 61, No. 3, mars 1973, pages 268-278). Néanmoins, dès que la réponse impulsionnelle des canaux devient longue ou que le nombre de valeurs discrètes possibles des symboles devient important, la complexité exponentielle de ces méthodes les rend impraticables. On considère le cas d'un canal de radiocommunication servant à la transmission d'un signal composé de séquences ou trames successives de n symboles dk (l de type QPSK (quaternary phase shift keying)... Après conversion en bande de base, numérisation et filtrage adapté, un vecteur Y du signal reçu reflétant les symboles émis sur la durée d'une trame a pour expression: 2 - rO 00 0 d1 YN,1 Y l À d Y2 rl r0 0 d2 YN,2 r1 r0 0 rO = = =A.D + YN (1) Y Yk rw r1 dk + YN,k 0 rw 0 0 rw 0YL 0 0 rW dn YN,L o W+1 est la longueur, en nombre de bits, de la réponse impulsionnelle estimée du canal, r = (r0,rl,...,rw) est la réponse impulsionnelle estimée du canal, les rq étant des nombres complexes tels que rq=O si q W, Yk est le k-ième échantillon complexe reçu avec l signal par l'émetteur et du filtrage de réception. La matrice A de taille Lxn a une structure de type Toeplitz le long de sa diagonale principale, c'est-à-dire que si ai.j désigne le terme situé à la i-ième ligne et à la j-ième colonne de la matrice A, alors ai+l,j+ l=ai,j pour 1 (structure de matrice-bande); et ail=ril pour 1 La relation matricielle (1) exprime que le signal reçu Y est l'observation, affectée d'un bruit additif, du produit de convolution entre la réponse impulsionnelle du canal et les symboles émis. Ce produit de convolution peut encore s'exprimer par sa transformée en Z: Y(Z) = R(Z).D(Z) + YN(Z) (2) o D(Z), Y(Z), R(Z) et YN(Z) sont les transformées en Z - 3 - respectives des symboles émis, du signal reçu, de la réponse impulsionnelle du canal et du bruit: E INCORPORER Equation.2 ECO (3) O INCORPORER Equation.2 WLE (4) E INCORPORER Equation.2 OW] (5) Une solution classique pour résoudre un système tel que (1) est la méthode dite de forçage à zéro (" zero forcing "), suivant laquelle on détermine le vecteur E INCORPORER Equation.2 [LIOZF à n composantes continues qui minimise l'erreur quadratique E INCORPORER Equation.2 ED.. Une discrétisation des composantes du vecteur E INCORPORER Equation.2 ELEZF relative à chaque canal intervient ensuite, souvent par le biais d'un décodeur de canal. La solution E INCORPORER Equation.2 lDzF au sens des moindres carrés est donnée par: DzF=(AHA)-IAHY, o AH désigne la matrice transposée conjuguée de A. On est alors ramené au problème de l'inversion de la matrice hermitienne définie positive AHA. Cette inversion peut être réalisée par divers algorithmes classiques, d'une manière directe (méthodes de Gauss, de Cholesky...) ou par des techniques itératives (algorithmes de Gauss-Seidel, du gradient...). L'erreur d'estimation D-DzF est égale à (AHA)-lAHYN, ce qui montre que la solution obtenue est affectée d'un bruit de variance: Y2 = E(|D-DZF) = N0x Trace[(AHA)1] (6) o No est la densité spectrale de puissance du bruit. On voit qu'il se produit une amplification du bruit (" noise enhancement ") quand la matrice AHA est mal conditionnée, c'est-à-dire quand elle a une ou plusieurs valeurs propres proches de 0. Cette amplification du bruit est le principal inconvénient des méthodes de résolution classiques. Dans la pratique, les cas de mauvais conditionnement de la -4- matrice AHA sont fréquents, particulièrement en présence de trajets multiples de propagation. On connaît un moyen relativement simple de remédier en partie à cet inconvénient, en acceptant dans la solution un résidu d'interférence, c'est-à-dire en adoptant non pas la solution optimale au sens des moindres carrés, mais la solution: DMMSE=(AHA+ No0)-lAHY, o N0 désigne une estimation de la densité spectrale du bruit, que le récepteur doit alors calculer. Cette méthode est connue sous le nom de MMSE (minimum mean square error) Elle permet de diminuer la variante d'estimation par rapport à la méthode de " zero forcing ", mais en introduisant un biais. Les méthodes de " zero forcing " et analogues reviennent à opérer un filtrage inverse du signal reçu par un filtre, modélisant la fonction de transfert 1/R(Z), calculé par une certaine approximation (quadratique dans le cas du " zero forcing "). Lorsqu'une ou plusieurs racines du polynôme R(Z) (équation (5)) sont situées sur le cercle unité, le filtre inverse théorique présente des singularités telles qu'il ne peut pas être estimé par une approximation satisfaisante. Dans le cas de l'approximation quadratique, ceci correspond à la divergence de la variance de l'erreur a2 lorsque la matrice AHA a une valeur propre nulle (relation (6)). Ce problème n'est pas rencontré dans les méthodes telles que l'algorithme de Viterbi qui prennent intrinsèquement en compte la nature discrète des symboles, mais qui requièrent une puissance de calcul très supérieure pour les systèmes de taille importante. La présente invention a pour but de proposer un procédé d'égalisation procurant un bon compromis entre la fiabilité des estimations et la complexité de l'égaliseur. Un autre but est d'obtenir un égaliseur nécessitant une puissance de calcul raisonnable et capable de traiter, avec des performances comparables à celles d'un égaliseur de Viterbi, des signaux dont les symboles ont un nombre d'états relativement élevés et/ou des signaux reçus suivant un canal de réponse impulsionnelle relativement étalée. L'invention propose ainsi un procédé d'égalisation numérique, pour estimer des symboles discrets d'un signal transmis à partir d'échantillons numériques d'un signal reçu par l'intermédiaire d'un canal de transmission représenté par une réponse impulsionnelle finie de W+l coefficients, W étant un entier plus grand que 1. Ce procédé comprend les étapes suivantes: - déterminer les W racines dans le plan complexe de la transformée en Z de la réponse impulsionnelle du canal; - répartir les W racines en un premier ensemble de W-p racines et un second ensemble de p racines, p étant un entier plus grand que 0 et plus petit que W, les racines du second ensemble étant plus proches du cercle unité que celles du premier ensemble selon un critère de distance déterminé dans le plan complexe; - obtenir un signal intermédiaire en appliquant au signal reçu une première méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z, formée par un polynôme en Z-1 de degré W-p, a pour racines les W-p racines du premier ensemble; et obtenir des estimations des symboles discrets du signal transmis en appliquant au signal intermédiaire une seconde méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z, formée par un polynôme en Z-1 de degré p, a pour racines les p racines du second ensemble. La " première méthode d'égalisation " sera généralement choisie de façon à traiter les symboles inconnus comme des variables continues. Elle conduit alors -- 6 -- à une opération semblable à un filtrage inverse dont la fonction de transfert serait d'une forme approchant l'expression 1/RS(Z), o RS(Z) désigne le polynôme en Z-1 de degré W-p ayant pour racines les W-p racines les plus éloignées du cercle unité. Elle peut notamment être du type " zero forcing ". Cette opération ne génère qu'une amplification du bruit réduite, puisque les racines de la fonction de transfert en Z associée sont relativement éloignés du cercle unité. 1Q Pour les p racines les plus proches du cercle unité, on adopte des mesures permettant de s'affranchir ou de limiter l'incidence du problème de l'amplification du bruit. On peut choisir une méthode MMSE ou analogue comme " seconde méthode d'égalisation ". Toutefois, cette seconde méthode tiendra avantageusement compte de la nature discrète des symboles inconnus. Elle pourra notamment reposer sur un algorithme à treillis, tel que l'algorithme de Viterbi, dont la mise en oeuvre est courante dans les égaliseurs de canal quand la taille du système n'est pas trop grande. La seconde méthode d'égalisation est généralement d'une mise en oeuvre plus complexe que la première. Dans chaque cas particulier, le choix du nombre p permet de rechercher le meilleur compromis entre la fiabilité des estimations, qui fait préférer les valeurs élevées de p, et la complexité de l'égaliseur, qui fait préférer les valeurs faibles de p. Un autre aspect de la présente invention se rapporte à un récepteur de radiocommunication, comprenant: - des moyens de conversion pour produire des échantillons numériques à partir d'un signal radio reçu par l'intermédiaire d'un canal de transmission représenté par une réponse impulsionnelle finie de W+l coefficients, W étant un entier plus grand que 1; - des moyens de mesure de la réponse impulsionnelle - 7 - du canal; des moyens de calcul des W racines dans le plan complexe de la transformée en Z de la réponse impulsionnelle mesurée; - des moyens de répartition des W racines en un premier ensemble de W-p racines et un second ensemble de p racines, p étant un entier plus grand que O et plus petit que W, les racines du second ensemble étant plus proches du cercle unité que celles du premier ensemble selon un critère de distance déterminé dans le plan complexe; - un premier étage d'égalisation pour obtenir un signal intermédiaire en appliquant au signal reçu une première méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z, formée par un polynôme en Z-1 de degré W-p, a pour racines les W-p racines du premier ensemble; et - un second étage d'égalisation pour obtenir des estimations de symboles discrets d'un signal transmis sur le canal en appliquant au signal intermédiaire une seconde méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z, formée par un polynôme en Z-1 de degré p, a pour racines les p racines du second ensemble. D'autres particularités et avantages de la présente invention 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 récepteur de radiocommunication selon l'invention; - la figure 2 est un organigramme montrant un mode de réalisation du procédé selon l'invention; et - la figure 3 est un diagramme illustrant les performances du procédé. Le récepteur représenté sur la figure 1 comporte un étage radio 1 qui reçoit le signal radio capté par -8- l'antenne 2 et le convertit en bande de base. Le signal en bande de base est numérisé par un convertisseur analogique-numérique 3, puis fourni à un filtre de réception 4. Le filtre 4 assure un filtrage adapté à la mise en forme des signaux par l'émetteur. Il délivre un signal numérique à raison d'un échantillon complexe par symbole émis. Ce signal numérique est fourni à un démodulateur comprenant d'une part un module 6 de synchronisation et d'estimation de canal, et d'autre part un égaliseur 7. La synchronisation et l'estimation de canal sont par exemple effectuées de manière classique à l'aide d'une séquence de synchronisation incluse par l'émetteur dans chaque trame de signal. La détection de cette séquence, connue du récepteur, permet d'une part de synchroniser le récepteur par rapport à la structure temporelle des trames émises, et d'autre part d'estimer la réponse impulsionnelle r = (r0,r1, .,rw) du canal sur lequel les trames sont transmises. La réponse impulsionnelle calculée..DTD: par le module 6 est fournie à l'égaliseur 7. L'égaliseur 7 fonctionne par exemple conformément à l'organigramme représenté sur la figure 2 pour traiter chaque trame synchronisée du signal reçu, se présentant (Y1 sous la forme d'un vecteur Y= y, avec L=n+ W en reprenant les notations précédentes. Le module d'estimation de canal 6 ayant fourni les W+1 coefficients complexes rq de la réponse impulsionnelle estimée du canal, la première étape 10 consiste à rechercher les W racines de la transformée en Z de cette réponse impulsionnelle, donnée par l'équation (5). Diverses méthodes classiques de recherche de racines complexes d'un polynôme peuvent être utilisées à l'étape 10. On pourra à cet égard se reporter à l'ouvrage de 9 - E. DURAND: " Solutions Numériques des Equations Algébriques; Tome I: Equations du Type F(x)=0 ", Editions Masson, 1960. Les W racines complexes ainsi trouvées al,a2,....,W sont ensuite ordonnées de façon à pouvoir les répartir en deux ensembles, l'un contenant les W-p racines les plus éloignées du cercle unité, et l'autre les p racines les plus proches du cercle unité. Pour cela, une distance 8q est calculée à l'étape 11 pour chacune-des racines aq (l 1 A l'étape 12, les racines aq de la fonction de transfert R(Z) sont triées dans l'ordre des distances décroissantes: 81 2 82 2...> &. On sépare alors les W-p premières racines al,...,Caw_p, qui sont les plus éloignées du cercle unité, des p racines restantes awp+l,...,aw. A l'étape 13, l'égaliseur 7 développe un polynôme en z-1 défini par: w-p W-p RS(z) = H (1- q.Z1) = E sq. q (8) q=l q=0 Ceci permet de déterminer les coefficients sq de la fonction de transfert RS(Z) associée à la réponse impulsionnelle s = (s0,sl,...,swp) d'un canal virtuel, qui correspondrait au canal de transmission estimé avec élimination des contributions les plus proches des zones de singularité. On peut alors procéder à une première égalisation 14 revenant à effectuer un filtrage inverse approchant la fonction de transfert 1/RS(Z). Plusieurs implémentations - 10 - peuvent être retenues pour effectuer ce filtrage inverse. On peut notamment effectuer une égalisation par " zero forcing " comme indiqué précédemment. Au sujet de ces méthodes, on pourra se reporter à l'ouvrage de J.G. Proakis: " Digital Communications " McGraw-Hill, 2e édition, 1989. Le filtrage inverse 14 produit un signal intermédiaire sous la forme d'un vecteur Y' de L'=n+p échantillons Y''... fyfL. Dans le cas d'une méthode de " zero forcing ", le vecteur Y' est obtenu par la relation matricielle: Y' = (A' H A')-1 A H y (9) Dans l'expression (9), A' désigne une matrice de n+W lignes et n+p colonnes ayant une structure de Toeplitz, formée à partir des coefficients sq du polynôme RS(Z): s O O O S1 So 0 s1 so '. 0 s1 À 0 A= sw-p s1 (10) 0 sw_p 0 0 sw_p 0. 0 oswp Grâce au tri des racines aq, les valeurs propres de la matrice A' H A' sont relativement éloignées de 0. En variante, on pourrait réaliser le filtrage inverse en mettant en cascade W-p cellules de filtrage correspondant chacune à l'inverse d'une fonction de transfert Rq (Z) = 1 - caqZ1, pour l - 11 - I q* 2 -2 m zm - 1 + q.- Z- + Sq. +.. + q. Z +.. (11) Rq (Z) Le développement (11) est causal, et stable puisque le domaine de convergence contient le cercle unité. La cellule de filtrage inverse peut donc être réalisée sous forme transverse ou sous forme récursive. Si czq > 1, on peut développer 1/Rsq (Z) sous la forme: 1 - -2 22 m zm S - -q. Z. 1 + aq. Z + aq.z +...+ q.Z +... (12) Rq (Z)) Ce développement (12) est anti-causal et stable. Pour la réalisation de la cellule de filtrage inverse, on tronque le développement (12) et on adopte une implémentation sous forme transverse. L'anti-causalité provoque un retard correspondant à la longueur de la réponse retenue. On note que les développements (11) et (12) justifient le critère de distance au cercle unité ôq utilisé conformément à la relation (7). A l'étape 15, l'égaliseur 7 développe un polynôme de degré p en Z-1, dont les racines correspondent aux p racines de R(Z) les plus proches du cercle unité, tel que R(Z) = RS(Z).RI(Z): w p RI(Z) = ro. H (1 - q z1) = 1 tq.Zq (13) q=W-p+1 q=0 Les coefficients complexes t définissent la réponse q impulsionnelle d'un autre canal de transmission virtuel, dont l'égalisation par une méthode de type " zero forcing " ou analogue poserait des problèmes d'amplification du bruit. Le signal intermédiaire Y' est alors soumis à une égalisation selon une autre méthode, sur la base de la réponse impulsionnelle t = (t0,tl,..., tp). Cette seconde égalisation 16 est avantageusement effectuée à l'aide d'un - 12 - treillis de Viterbi (voir l'article précité de G.D. Forney Jr., ou l'ouvrage précité de J.G. Proakis). Le second étage d'égalisation 16 produit les estimations dk des symboles de la trame (1 d'éventuelles erreurs de transmission. La figure 3 illustre les performances du procédé ia dans le cas de la transmission d'une trame de signal selon le format du système radiotéléphonique cellulaire européen GSM, en remplaçant la modulation binaire de type GMSK par une modulation de phase à huit états (modulation 8-PSK). La réponse impulsionnelle du canal était tronquée à cinq temps bits (W=4). La figure 3 montre la dépendance entre le taux d'erreur binaire BER, exprimé en %, et le rapport Eb/N0 entre l'énergie par bit et la densité spectrale du bruit, exprimé en décibels. Le BER est celui observé dans les estimations des symboles après le désentrelacement et le décodage de canal effectués conformément aux méthodes employées dans le GSM. La courbe I montre les résultats procurés par la méthode de " zero forcing " pur, c'est-à- dire dans le cas limite o p=0. La courbe II montre le résultat théorique qui serait obtenu en égalisant le canal purement avec l'algorithme de Viterbi (cas limite o p=W). Dans la pratique, le treillis correspondant devrait comporter 84=4096 états, de sorte que l'égaliseur de Viterbi serait irréalisable avec les techniques actuelles. L'écart entre les courbes I et II illustre la supériorité de l'algorithme de Viterbi qui délivre les estimations selon le maximum de vraisemblance. Les courbes III et IV montrent les résultats obtenus par le procédé selon l'invention, respectivement dans les cas o p=l et p=2. On voit l'amélioration très sensible des résultats déjà obtenue pour la valeur p=l par rapport - 13 - au " zero forcing " pur. A titre indicatif, l'égalisation d'une trame de signal GSM par l'algorithme de Viterbi pur dans les conditions de la figure 3 requerrait de l'ordre de 8,45 millions d'opérations en virgule flottante, soit environ 1,83 Gflops, alors que la mise en oeuvre de la présente invention dans les mêmes conditions requiert de l'ordre de 19000 opérations en virgule flottante (,4,2 Mflops) dans le cas o p=l, y compris la recherche des racines de R(Z) et le filtrage inverse 1/RS(Z) par la méthode " zero forcing ". Ce nombre est de l'ordre de 129000 opérations (=28 Mflops) dans le cas o p=2, ce qui reste compatible avec la puissance des processeurs de traitement de signal (DSP) actuellement disponibles. - 14 - R E V E N D I C A T I ON S 1. Procédé d'égalisation numérique, pour estimer des symboles discrets (dk) d'un signal transmis à partir d'échantillons numériques (yk) d'un signal reçu par l'intermédiaire d'un canal de transmission représenté par une réponse impulsionnelle finie de W+1 coefficients (ro,rl,...,rw), W étant un entier plus grand que 1, comprenant les étapes suivantes: déterminer les W racines (al,a2,.. ,Cw) dans le plan complexe de la transformée en Z (R(Z)) de la réponse impulsionnelle du canal; - répartir les W racines en un premier ensemble de W-p racines (al,..., W_p) et un second ensemble de p racines (Cw_p+l,...,aw), p étant un entier plus grand que 0 et plus petit que W, les racines du second ensemble étant plus proches du cercle unité que celles du premier ensemble selon un critère de distance déterminé dans le plan complexe; - obtenir un signal intermédiaire (Y') en appliquant au signal reçu (Y) une première méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z (RS(Z)), formée par un polynôme en Z-1 de degré W-p, a pour racines les W-p racines du premier ensemble; et - obtenir des estimations (dk) des symboles discrets du signal transmis en appliquant au signal intermédiaire une seconde méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z (RI(Z)), formée par un polynôme en Z-1 de degré p, a pour racines les p racines du second ensemble. 2. Procédé selon la revendication 1, dans lequel la première méthode d'égalisation produit le signal - 15 - intermédiaire sous la forme d'un vecteur Y' de n+p échantillons (Y',,Y'n+ p) obtenu selon la relation: Y' = (A'H A)-1 Ai H y o Y est un vecteur formé de n+W échantillons (yk) du signal reçu, et A' est une matrice de n+W lignes et n+p colonnes ayant une structure de Toeplitz formée à partir des coefficients (Sq) dudit polynôme en Z-1 de degré W-p (Rs(Z)) 3. Procédé selon la revendication 1 ou 2, dans lequel la seconde méthode d'égalisation comporte la mise en oeuvre d'un algorithme de Viterbi. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le critère de distance au cercle unité, utilisé pour répartir les W racines al,...,aw de la transformée en Z (R(Z)) de la réponse impulsionnelle du canal entre les premier et second ensembles, s'exprime par une distance Èq de la forme Èq = 1 - q si |q de la forme 8q = 1 - 1/|qq si q > 1, pour 1 5. Récepteur de radiocommunication, comprenant: - des moyens de conversion (1,3,4) pour produire des échantillons numériques (Yk) à partir d'un signal radio reçu par l'intermédiaire d'un canal de transmission représenté par une réponse impulsionnelle finie de W+1 coefficients (r0,r1,...,rw), W étant un entier plus grand que 1; - des moyens (6) de mesure de la réponse impulsionnelle du canal; - des moyens de calcul des W racines (al,U2,...,aw) dans le plan complexe de la transformée en Z (R(Z)) de la réponse impulsionnelle mesurée; - 16 - - des moyens de répartition des W racines en un premier ensemble de W-p racines (al,...,Wp) et un second ensemble de p racines (,W_p+i,...,-W), p étant un entier plus grand que O et plus petit que W, les racines du second ensemble étant plus proches du cercle unité que celles du premier ensemble selon un critère de distance déterminé dans le plan complexe; un premier étage d'égalisation pour obtenir un signal intermédiaire en appliquant au signal reçu (yk) une première méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z (RS(Z)), formée par un polynôme en Z-1 de degré W-p, a pour racines les W-p racines du premier ensemble; et - un second étage d'égalisation pour obtenir des estimations (dk) de symboles discrets d'un signal transmis sur le canal en appliquant au signal intermédiaire une seconde méthode d'égalisation sur la base d'une réponse impulsionnelle finie dont la transformée en Z (RI(Z)), formée par un polynôme en Z-1 de degré p. a pour racines les p racines du second ensemble. 6. Récepteur selon la revendication 5, dans lequel le premier étage d'égalisation est agencé pour produire le signal intermédiaire sous la forme d'un vecteur Y' de n+p échantillons (Y'1,'.. Y'n+p) obtenu selon la relation: Y' = (A'H A')-' A'H y o Y est un vecteur formé de n+W échantillons (yk) du signal reçu, et A' est une matrice de n+W lignes et n+p colonnes ayant une structure de Toeplitz formée à partir des coefficients (Sq) dudit polynôme en Z-1 de degré W-p (RS(Z)). - 17 - 7. Récepteur selon la revendication 5 ou 6, dans lequel le second étage d'égalisation est agencé pour mettre en oeuvre un algorithme de Viterbi. 8. Récepteur selon l'une quelconque des revendications 5 à 7, dans lequel les moyens de répartition des racines utilisent, pour répartir les W racines entre les premier et second ensembles, un critère de distance au cercle unité s'exprimant par une distance q de la forme 3q = 1 - |q si cql 6q = 1 - 1/|aq si |q > 1, pour l1