La présente invention concerne un procédé et un système de réception radio haut débit pour recevoir des signaux dans un système radio utilisant la technique CDMA (Code Division Multiple Access).

Elle s'applique notamment, mais non exclusivement aux systèmes de téléphonie mobile haut débit (au moins 2 Mbits/s), tels que les systèmes UMTS (Universal Mobile Telecommunications System) en Europe. Ces systèmes sont conçus pour offrir une large gamme de services ayant des débits et des facteurs d'étalement différents. Ces services incluent un mode de transmission par paquet haut débit en liaison descendante (vers les terminaux mobiles) appelé HSDPA (High Speed Downlink Packet Access), qui est caractérisé par l'abandon du contrôle de puissance en boucle fermée, et une adaptation de lien utilisant des modulations à constellations variables (QPSK, MAQ16, MAQ64) et un faible facteur d'étalement.

Dans les systèmes de téléphonie mobile, l'information est transmise en utilisant une technique d'accès multiples. Certains systèmes utilisent une division en fréquence (FDMA) et temporelle (TDMA), les utilisateurs du réseau étant distingués les uns des autres à la fois sur la base de la fréquence respective utilisée, les informations à transmettre étant réparties dans des tranches de temps attribuées respectivement à chaque utilisateur. Dans les systèmes basés sur une division de codes (CDMA), les utilisateurs communiquent entre eux en utilisant la même bande de fréquences radio. Pour pouvoir distinguer les utilisateurs les uns des autres, un code d'étalement est attribué respectivement à chacun d'entre eux pour toute la durée d'une connexion, ce code servant à étaler en bande de base le spectre du signal à transmettre. Pour reconstituer l'information transmise les récepteurs doivent utiliser le même code pour effectuer l'opération inverse de l'étalement. Comparée aux autres méthodes d'accès multiples, cette technique présente l'avantage d'être plus souple en terme d'accès et de débit que l'on peut faire varier en jouant sur le facteur d'étalement.

Dans les transmissions radio, la forme du milieu entre l'émetteur et le récepteur d'un signal radio interfère avec les transmissions effectuées, et entraîne une propagation suivant des trajets multiples dus aux réflexions qui se produisent en

différents points du canal radio (en particulier en milieu urbain). Il en résulte que les composantes d'un même signal transmis arrivent au récepteur avec des puissances et des retards différents.

Pour les récepteurs CDMA, on utilise un récepteur en râteau, appelé récepteur Rake, qui permet de reconstituer le signal transmis à partir des différentes composantes reçues ayant suivi des chemins de propagation différents. Ces récepteurs sont basés sur la reconstitution d'un profil de délais ou modèle équivalent de canal radio. A cet effet, on transmet avec l'information une séquence de symboles pilotes connue des récepteurs, et sur la base de cette connaissance, les récepteurs effectuent une estimation (réponse impulsionnelle représentant l'ensemble des trajets du canal radio) du canal radio par lequel le signal reçu a été transmis. Un filtre adapté est décalé sur le signal reçu, par exemple d'une demi unité de code d'étalement, pendant que la puissance reçue est mesurée. Cette technique permet de construire un graphe de réponse impulsionnelle donnant des informations sur la puissance et les retards dus à la propagation par trajets multiples des composantes du signal reçu par un canal radio donné.

Si la technique CDMA apparaît bien adaptée aux services temps réel à faibles débits, elle semble mal adaptée aux services hauts débits par paquets. En effet, les performances du récepteur Rake reposent sur les propriétés d'inter-et d'auto- corrélation des séquences d'étalement, et ces propriétés sont d'autant meilleures que la séquence d'étalement est longue et donc que le facteur d'étalement est grand. Or plus le débit augmente, plus le facteur d'étalement diminue. La séquence d'étalement devient plus courte et donc les propriétés d'inter-et d'auto-corrélation des séquences d'étalement se dégradent, conduisant à l'apparition d'interférences entre symboles d'un même signal transmis. Il en résulte que les performances du récepteur Rake s'écroulent pour des facteurs d'étalement inférieurs à 8, en particulier lorsque la modulation utilisée présente un grand nombre d'états.

Une étude de la dégradation des performances du récepteur Rake à cause des interférences entre symboles a montré que l'utilisation d'une technique d'égalisation est nécessaire lorsque le facteur d'étalement est inférieur à 16 (voir [1]"On the Rake receiver performance"H. Boujemaa, M. Siala, VTC 2000 Fall, Boston USA).

Il s'avère donc que les récepteurs de type Rake apparaissent seuls très mal adaptés pour répondre aux besoins de la téléphonie mobile haut débit.

Or la mise en oeuvre de techniques de détection et de décodage optimales s'avère à l'heure actuelle quasi impossible, car elle entraîne une complexité de calcul très importante, en particulier dans le cas où les canaux de transmission présentent une réponse impulsionnelle trop longue comme en milieu urbain.

Par ailleurs, dans la technique TDMA (Time Division Multiple Access), il existe plusieurs techniques de détection et de décodage sous-optimales.

On connaît par exemple la technique utilisant un égaliseur de type LMMSE (Linear Minimum Mean Square Error) permettant de réduire non seulement les interférences entre symboles, mais également les Interférences entre utilisateurs.

Pour plus de détails, on peut se référer aux documents : - [2] "Linear receivers for the DS-CDMA downlink exploiting orthogonality of spreading sequences"de I. Ghauri et D. T. M Slock, dans Proc. 32th Asilomar Conf. on Signals, Systems and Comp., Asilomar, CA, Nov. 1-4 1998, et - [3]"Interference Suppression in CDMA Downlink through Adaptive Channel Equalization", de M. Heikkilâ, P. Komulainen, et J. Lilleberg, dans VTC 99 Fall, Tokyo, Japan.

Il est à noter que dans les liaisons descendantes, les interférences entre utilisateurs sont dues au fait que la propagation dans les canaux s'effectue suivant des trajets multiples, étant donné que les séquences d'étalement entre les différents utilisateurs sont orthogonales. Cette solution s'avère très efficace pour réduire les interférences entre utilisateurs, par comparaison avec les récepteurs de type Rake. Toutefois, cette solution ne permet pas de réduire sensiblement les interférences entre symboles du fait des caractéristiques linéaires de l'égaliseur LMMSE.

On a également proposé de disposer en sortie d'un récepteur Rake, un égaliseur à estimation de la séquence la plus probable (MLSE : Maximum Likelihood Sequence Estimation). Pour plus de détails concernant cette technique, on peut se référer par exemple au document [4] "Joint multipath combining and MLSE equalization (Rake-MLSE Receiver) for WCDMA systems"de S. Tantikovit, et A. U. H. Sheikh, dans VTC 2000 Spring, Tokyo, Japan.

Cette solution est optimale du point de vue de la détection de séquence, et

proche de la solution optimale sur le plan de la détection d'erreurs sur les symboles transmis. Toutefois, la complexité de cette solution augmente d'une manière exponentielle avec l'étalement des retards dans un même canal et avec la taille de la constellation de la modulation employée. Elle n'est donc pas applicable à tous les services UMTS. En outre, cette solution ne prend pas en compte les dégradations de performances dues à une mauvaise estimation de canal et de codage de canal. De plus, elle ne propose pas d'algorithme de sortie souple ou pondérée.

On connaît par ailleurs la technique d'égalisation DDFSE (Delayed Decision Feedback Sequence Estimation) qui consiste à réduire la complexité des états du treillis en utilisant la technique dite de"traitement par survivant" (Per Survivor Processing). Pour plus de détails concernant cette technique, on peut se référer par exemple au document [5] "Delayed Decision-Feedback Sequence Estimation"de A. Duel-Hallen et C. Heegard, dans IEEE Transactions on Communications, vol. 37, pp. 428-436, May 1989.

Dans les systèmes TDMA, cette technique présente l'inconvénient d'être sensible à la propagation d'erreurs, ce qui nécessite un préfiltrage. Dans les systèmes CDMA, cette technique apparaît inapplicable étant donné que le modèle équivalent de canal en sortie du récepteur Rake varie à chaque symbole transmis, du fait qu'il dépend du code d'étalement qui change à chaque symbole.

Une méthode de détection et de décodage par itération conçue pour la technique TDMA, appelée"turbo-détection"å également été décrite dans l'article [6] "Iterative Correction of Intersequential Interference : Turbo-equalization"de C.

Douillard, M. Jezequel, C. Berrou, A. Picart, P. Didier et A. Glavieux, publié dans European Transactions on Telecommunications, vol. 6, p. 507 à 511, Sept.

1995. Dans cette technique de détection et de décodage, on utilise un égaliseur à maximum de vraisemblance à entrées et sorties pondérées, appelé SISO MLSE, et le processus de décodage est de type Viterbi également à entrées et sorties pondérées, désigné par SOVA. Ce dernier a fait l'objet d'une publication intitulée [7] "A low Complexity Soft Output Viterbi Decoder Architecture", ICC'93 p. 733 à 740, Genève, Suisse, Mai 1993.

Cette technique de détection et de décodage a fait l'objet de développements importants qui ont conduit à la mise en oeuvre de détecteurs optimaux à maximum de probabilité a posteriori (MAP). Pour plus de détails sur ces détecteurs, on peut se reporter aux articles intitulés :

- [8]"Optimal Decoding of Linear Codes for Minimizing Symbol Error Rate" publié par L. R. Bahl, J. Cocke, F. Jelinek et J. Raviv dans IEEE Transactions on Information Theory, vol. IT-20, p. 284-287, Mars 1994 ; - [9]"Iterative Equalization and Decoding in Mobile Communications Systems", publié par G. Baush, H. Khorram et J. Hagenauer dans Proc.

EPMCC'97, p. 307-312, Bonn, Allemagne, Sept. 1997.

Toutefois, cette solution n'a pas été transposée aux systèmes CDMA, et en outre elle introduit une complexité démesurée au niveau du récepteur, de l'ordre de ML, M étant le nombre de points de la constellation de modulation et L le nombre d'échos pris en compte dans le canal de propagation. De plus, elle n'aborde pas le problème de l'estimation de canal.

Enfin, dans la publication [10] "Turbo-Equalization over Frequency Selective Channel"-International Symposium on Turbo-Codes, Brest, France, Sept.

1997, on a également proposé une technique itérative de détection de symboles et de décodage de canal, appelée"turbo-égalisation"substantiellement différente de la technique de turbo-détection précitée, et qui suppose une estimation bruitée du canal de transmission. Toutefois, par rapport à la technique de turbo-détection, cette technique introduit des dégradations de performances, dépendant fortement de la technique d'égalisation employée à la première itération. A ce sujet, on peut se référer à la publication [11]"Joint Equalization and Decoding : Why Choose the Iterative Solution ? "de A.

Roumy, I. Figalkow et D. Pirez, dans IEEE VTC'1999 Fall, Amsterdam, Netherlands, Sept. 1999.

Cette technique n'est pas transposable aux systèmes CDMA, car elle repose sur des techniques de filtrage qui ne peuvent pas être appliquées lorsque le canal varie indépendamment d'un symbole transmis à l'autre.

La présente invention a pour but de résoudre ces problèmes en vue de proposer un procédé de réception et une structure de récepteur, relativement simples, présentant des performances proches des performances optimales, tout en utilisant une modulation d'ordre élevé et un facteur d'étalement faible. Cet objectif est atteint par la prévision d'un procédé de réception d'un signal, transmis sur un canal de transmission à trajets multiples, suivant une technique d'étalement de spectre avec un faible facteur d'étalement, ce signal étant

transmis sous la forme de séquences de symboles binaires codés comportant des symboles pilotes prédéfinis et des symboles de données, multipliée par une séquence d'étalement, ce procédé comprenant une étape de détermination d'une estimation du canal de transmission à l'aide des symboles pilotes reçus.

Selon l'invention, ce procédé comprend les étapes consistant à : - reconstituer le signal transmis à l'aide de l'estimation de canal à partir des signaux reçus et transmis par les multiples trajets du canal de transmission, - déterminer un modèle équivalent de canal tel qu'il est vu à partir du signal reconstitué, à partir de l'estimation de canal, - réduire les interférences entre les symboles du signal transmis reconstitué, ces interférences résultant d'un facteur d'étalement faible, à l'aide d'un détecteur DDFSE à base de treillis à nombre d'états réduit et utilisant le modèle équivalent de canal, et délivrant en sortie des valeurs estimées des symboles codés reçus, - désentrelacer les symboles codés, et - décoder les valeurs estimées des symboles codés et désentrelacés, pour restituer les symboles de données transmis.

Avantageusement, les valeurs estimées des symboles codés reçus, obtenues après réduction d'interférence et après décodage sont des valeurs pondérées ou souples.

Selon une particularité de l'invention, les étapes de réduction des interférences entre symboles et de décodage sont incluses dans un premier processus itératif dans lequel les symboles codés et désentrelacés obtenus à une itération n sont réestimés lors du décodage en fonction des symboles de données obtenus à l'issue du décodage et après correction d'erreurs, la différence entre les symboles codés réestimés obtenus à l'itération n et les symboles codés, désentrelacés, obtenus à la même itération est réentrelacée, puis appliquée en entrée du détecteur DDFSE et retranchée aux symboles codés obtenus à l'itération n+1 en sortie du détecteur DDFSE.

Avantageusement, l'estimation de canal est améliorée par la méthode des moindres carrés (LS).

De préférence, l'estimation de canal est améliorée par l'algorithme de l'erreur quadratique moyenne minimum (MMSE).

Selon une autre particularité de l'invention, l'estimation de canal est effectuée dans un processus itératif dans lequel les symboles codés et désentrelacés obtenus à une itération n sont réestimés lors du décodage, en fonction des symboles de données obtenus à l'issue du décodage et après correction d'erreurs, les symboles codés réestimés étant entrelacés, une estimation de canal étant effectuée sur la base des symboles codés réestimés et entrelacés, un modèle équivalent de canal étant déterminé à partir de l'estimation de canal, l'estimation de canal et le modèle équivalent de canal déterminés à une itération n étant utilisés respectivement pour reconstituer le signal transmis et réduire les interférences entre symboles à l'itération suivante n+1.

L'invention concerne également un système de réception d'un signal, transmis sur un canal de transmission à trajets multiples, suivant une technique d'étalement de spectre avec un faible facteur d'étalement, ce signal étant transmis sous la forme d'une séquence de symboles binaires codés comportant des symboles pilotes prédéfinis et des symboles de données, multipliée par une séquence d'étalement, ce système comprenant un récepteur en râteau pour reconstituer le signal transmis à l'aide d'une estimation de canal, à partir des signaux reçus et transmis par les multiples trajets du canal de transmission, et des moyens d'estimation de canal à partir de symboles pilotes reçus, pour délivrer une estimation du canal de transmission au récepteur en râteau.

Selon l'invention ce système comprend en outre : - des moyens de modélisation de canal pour déterminer un modèle équivalent de canal tel qu'il est vu en sortie du récepteur en râteau, en fonction de l'estimation de canal, - des moyens de réduction d'interférences entre symboles reçus, comportant un détecteur DDFSE à base de treillis à nombre d'états réduit, pour réduire les interférences entre symboles reçus à l'aide du modèle équivalent de canal, et restituer des valeurs estimées des symboles codés reçus, - des moyens de désentrelacement pour désentrelacer les valeurs estimées des symboles codés reçus, et - des moyens de décodage pour décoder les valeurs estimées et désentrelacées et fournir les symboles de données transmis.

Avantageusement, les valeurs estimées et décodées délivrées par les moyens de réduction d'interférences entre symboles et les moyens de décodage sont des valeurs pondérées ou souples.

Selon une particularité de l'invention, ce système comprend en outre : - des moyens pour réestimer les symboles codés en fonction des symboles de données décodés et après correction d'erreurs, - des premiers moyens de soustraction pour soustraire aux symboles codés réestimés, les symboles codés estimés et désentrelacés, et obtenir une séquence de symboles codés réestimés extrinsèque, - des premiers moyens d'entrelacement pour entrelacer la séquence de symboles codés réestimés extrinsèque, et - des seconds moyens de soustraction pour soustraire de la séquence de symboles reçus et estimés par les moyens de réduction à une itération suivante, la séquence de symboles codés réestimés extrinsèque.

Selon une autre particularité de l'invention, ce système comprend en outre : - des seconds moyens d'entrelacement pour entrelacer la séquence de symboles codés réestimés en sortie des moyens de décodage, - des second moyens d'estimation de canal pour fournir une estimation de canal de transmission à partir de la séquence entrelacée de symboles codés réestimés, aux moyens de détermination d'un modèle équivalent de canal et au récepteur en râteau.

Un mode de réalisation préféré de l'invention sera décrit ci-après, à titre d'exemple non limitatif, avec référence aux dessins annexés dans lesquels La figure 1 représente schématiquement un émetteur classique conçu pour émettre des signaux selon la technique CDMA ; La figure 2 représente schématiquement un système de réception selon l'invention ; Les figures 3 et 4 représentent schématiquement deux variantes préférées du système de réception montré sur la figure 2.

Les figures 5 et 6 illustrent sous la forme de courbes de taux d'erreur de bit en fonction du rapport signal sur bruit, les performances des récepteurs représentés sur les figures 2 à 4 ;

La figure 7 représente schématiquement un mode de réalisation d'un estimateur de canal selon l'invention.

La figure 1 représente un émetteur CDMA selon l'art antérieur, un tel émetteur comprenant une source 1 de signaux fournissant des séquences de symboles binaires notées ui° _ {ul,..., uzoT de TO symboles, un encodeur de canal 2 qui délivre une séquence codée notée clo Cl,..., C T Chaque symbole de donnée un = {un,1, ..., un,k0)T contient ko bits, tandis que chaque symbole cn = {cn,1, ..., cn,n0}T contient no bits. Les bits codés sont entrelacés par un entrelaceur 3 et complétés pour correspondre à un format de transmission prédéterminé, c'est-à-dire des trames de longueur-c incluant des symboles pilotes pour permettre au récepteur d'effectuer une estimation de canal. Les bits résultants sont groupés en symboles de type ak = (ak akq) contenant q bits avant d'être appliqués à un modulateur 4 de type M-PSK (Phase-Shift Keying) d'ordre M, MAQ16 ou MAQ64, qui délivre un symbole modulé correspondant s (k).

Ensuite, le signal s (k) est multiplié en 5 à une séquence d'étalement prédéfinie c (q) pour la transmission en question, le signal résultant traversant un filtre 6 de type en racine carrée de cosinus surélevé également appelé en racine de Nyquist, avec un facteur d'occupation spectrale (ou roll-off) de 0,22.

Sur la figure 2, le récepteur selon l'invention comprend un récepteur 11 en râteau de type Rake qui utilise une estimation de canal fournie par un estimateur de canal 10, le récepteur Rake étant suivi d'un dispositif de réduction d'interférences entre symboles 12, auquel est appliqué un modèle équivalent 16 du canal de transmission. La sortie du dispositif de réduction d'interférences 12 est connectée à un décodeur 15, par l'intermédiaire d'un désentrelaceur 14 effectuant l'opération inverse de l'entrelacement 3. Par ailleurs, l'estimation de canal fournie par l'estimateur 10 est appliquée au modèle de canal 16 pour déterminer un modèle équivalent de canal.

Selon l'invention, le dispositif 12 de réduction d'interférences entre symboles est du type DDFSE (Delayed Decision Feed-back Sequence Estimation) à entrée et sortie pondérées, c'est-à-dire de la forme d'un logarithme ln [p (l)/p (0)] du rapport de probabilités p (l) et p (0) que le signal soit respectivement au niveau logique 1 et 0. De même, le décodeur 15 présente également des entrées et sorties pondérées.

Dans le cas d'une propagation multi-trajets, le signal reçu à l'entrée du récepteur à étalement de spectre au temps t présente la forme suivante : L 7) = () (Â ;) e, (-- ()) + ) (1) 1=1 k dans laquelle ek (t)-EN-'e (kN + q) g (t-qT,) est la forme d'onde pour le symbole modulé s (k), e (q) est la séquence d'étalement, N est le facteur d'étalement, g (t) est la fonction de transfert du filtre en racine de Nyquist 6, Te et T, sont respectivement les périodes dite"chip"et de symbole, L est le nombre de trajets dans le canal, hl (t) et zl (t) sont respectivement l'amplitude complexe et le retard du 1-ième trajet, et (t) est le bruit blanc gaussien avec une densité spectrale de puissance No.

Si d (i) est le résultat du produit des symboles transmis par la séquence d'étalement, c (i), s L , L J représentant la fonction"partie entière". U-K Le signal reçu peut également être écrit de la manière suivante : r (t) = E d (i) h' (t-iT,) (2) dans laquelle : L h' (t) = E hl (t + iTC) g (ttz (t + iTC)) (3) 1=1 Conformément au théorème de Nyquist, le signal reçu r (t) est échantillonné à deux fois la fréquence de chip pour obtenir un vecteur r (i) des échantillons empilés, qui est utilisé pour estimer le symbole de chip i.

Le vecteur r (i) présente la forme suivante : dans laquelle Ml et M2 représentent la longueur de h'(t) en multiple de Ter - H(i) = [hM1+M2(i),###,h0(i),h1(i),###,hM1+M2(i)] (5)

[hi((j-M1)Tc),hi((j-M1)Tc+Tc/2),###<BR> <BR> <BR> <BR> <BR> ###,hi(M2Tc),01,2j]T, 0#j#M1+M2,<BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> hj 2j,hi(-M1Tc),hi(-M1Tc+Tc/2),###<BR> <BR> <BR> <BR> <BR> <BR> ###,hi((M2-j)Tc)]T, 0#j#M1+M2, et d (i) =[d(i-M1-M2),..., d (i),###,d(i + Mi + Mf (8) Le signal réduit (après opération inverse de l'étalement) sur la branche j (trajet j du canal de transmission) et pour le symbole k peut être écrit de la manière suivante : dans laquelle dk est le produit de la séquence d'étalement par les symboles transmis et zji = #j-#i)/Tc.

A partir des résultats du document [1], on peut facilement montrer que la sortie ôk du récepteur Rake 11 peut être représentée de la manière suivante : L L' 0 k i (10) =1 dans laquelle : L j=l (11) j=l gl (k) est la 1-ième amplitude du modèle équivalent en sortie du récepteur Rake 11, (2L'+1) est le nombre d'échos du modèle équivalent 16.

Si l'on suppose que les retards de trajet sont espacés d'un multiple de la période de chip Tc, les paramètres du modèle équivalent 16 sont donnés par les formules suivantes :

Les courbes représentées sur la figure 5 donnent le taux d'erreur binaire BER en fonction du rapport signal sur bruit Eb/No, dans le cas idéal (courbe Cl) et en sortie du récepteur Rake (courbe C2).

Ces courbes de même que celles de la figure 6, ont été obtenues par simulation avec un facteur d'étalement de 4, un canal de transmission EQ-4 avec quatre trajets dont les retards sont séparés respectivement d'une période chip, chaque trajet présentant une forme gaussienne complexe circulaire ou présentant une atténuation de Rayleigh. Le code de sortie est un code systématique récursif à 16 états, de rapport 1/2 et de polynômes de génération générant une séquence préencodée c qui est envoyée à un entrelaceur pseudo-aléatoire et divisé en trames.

La comparaison des courbes Cl et C2 montre que les performances d'un tel récepteur sont très mauvaises.

Si comme cela est préconisé dans les références [2] et [3], on dispose un égaliseur LMMSE en amont du récepteur Rake, pour réduire la dégradation due

aux interférences entre symboles résultant d'un facteur d'étalement faible, l'estimation du i-ième symbole de chip est de la forme suivante : dans laquelle d est la variance de la séquence de chip.

La courbe C3 sur la figure 5 qui illustre les performances obtenues avec un tel égaliseur, montre que cette solution ne permet pas d'améliorer sensiblement les performances du récepteur Rake.

Pour améliorer sensiblement les performances du récepteur Rake 11 selon l'invention, on dispose en sortie de celui-ci, un dispositif 12 de réduction d'interférences entre symboles conçu autour d'un détecteur DDFSE sous- optimal basé sur un treillis à nombre d'états réduit, un tel détecteur utilisant un modèle équivalent de canal 16.

Avant d'être appliqués au détecteur DDFSE, les échantillons sont retardés de manière à rendre causal le modèle équivalent de canal 16. Ils peuvent alors être représentés par la formule suivante : dans laquelle le vecteur h (k) = gl-L (k), ce vecteur représentant le vecteur des coefficients de canal [ho (k),..., @ h2Lk)] T.

Ensuite, le détecteur DDFSE travaille sur un treillis à nombre d'états réduit (voir document [8] ), par rapport à la technique BCJR qui applique un critère MAP (à maximum de probabilité a posteriori) à un treillis complet à Q2L états, Q étant le nombre de points de la constellation de modulation PSK, et (2L41) étant le nombre de trajets du modèle équivalent 16.

Un tel treillis est en fait une machine à états finis, étalés dans le temps et dont les transitions entre états ne dépendent que de l'état précédent, le nombre d'états à chaque instant étant constant. Dans un tel treillis, une section représente l'ensemble des transitions entre les états correspondants à deux instants successifs.

Dans le détecteur DDFSE, le treillis est ainsi réduit à un nombre d'états Que, vr étant un nombre entier positif appelé mémoire réduite et choisi tel que vr < 2L' dans le cas du détecteur DDFSE.

On dit généralement qu'une séquence d'entrée cq de treillis se termine par un sous-état s si cq se termine par la sous-chaîne s = an-#n +1. A une profondeur n, r l'espace des sous-états Sn coincide avec l'espace complet Sn des états du treillis BCJR si vr = 2L'. Si or < 2L', Sn est réduit à un sous-ensemble comprenant tous les sous-états possibles s dérivés de tous les états : #Sn#=Q#r,#n#[0,#] et #Bn#=Q#r+1,#n#[1,#] (20) La notation employée dans ces relations s'applique à la définition d'un sous- treillis T (S, B) sur lequel s'applique l'algorithme DDFSE.

Dans chaque section et pour toutes les transitions, le calcul de métrique de branche implique une convolution de la réponse impulsionnelle en temps discret du canal à la séquence de 2L'+1 symboles déjà estimés. Seuls les premiers #r+1 symboles estimés pour cette séquence étant disponibles sur la transition en cours de traitement et sur le sous-état de sous-treillis de départ avec lequel elle est connectée.

A chaque indice temporel ne [l, T] et pour tous les indices de bit je [l, q], l'algorithme optimal symbole par symbole BCJR fourni les logarithmes des rapports de probabilité a posteriori, selon la formule suivante : dans laquelle h est une estimation (ou ré-estimation) du vecteur transverse des coefficients de canal (si possible converti à phase minimale), et Y1 est une séquence observée de longueur #. Dans la dérivée qui suit, le conditionnement par # est implicite et omis pour simplifier les expressions.

Si l'on effectue une marginalisation sur les séquences de symboles d'entrée à bits marqués, l'équation (21) peut être réécrite de la manière suivante :

avec Ak représentant des quantités non négatives, le logarithme du rapport de probabilité a posteriori , (ans) peut être évalué à l'aide de la formule suivante : <BR> <BR> <BR> <BR> <BR> # min {-ln p (a1#,y1#)}- min {- ln p (a1#,y1#)} (2<BR> <BR> a1#,an,j=0 a1#,an,j=1 dans laquelle {-lnp(a1#,y1#)} représente le coût métrique, correspondant au bruit, du chemin dans le treillis, associé à la séquence d'entrée a1 et à la séquence reçue Y1. Du fait de la réduction du treillis, le dispositif DDFSE 12 évalue la quantité (-lnp(a1#,y1#)} d'une manière sous-optimale sur la base de l'algorithme PSP dit"par survivant"consistant à sélectionner un seul survivant par noeud.

Pour un sous-treillis T (SB) donné et une branche métrique particulière, on note p (b) le coût métrique du meilleur chemin commençant au sous-état 0 à la profondeur 0, et se terminant au sous-état 0 à la profondeur a (si l'on prend en compte les symboles de queue de la séquence), et passant par la branche b E Bn de la section n. On suppose également que chaque branche b (E Bn contient trois champs : un champ sous-état de départ b-E : Sn-,, un champ sous-état d'arrivée b+ E Sn et un champ libellé modélisant un symbole d'entrée à label de bit pour le code d'interférence entre symboles convolutionnel de taux 1 variant en fonction du temps à l'instant n. La sortie du dispositif DDFSE 12 peut être représentée de la manière suivante : <BR> <BR> <BR> <BR> <BR> min µn#(b) - min µn# (b)<BR> <BR> <BR> b#Bn,bj#=0 b#Bn,bj#=1 Le coût métrique µn#(b) considéré dans la formule précédente peut toujours être décomposé en une somme de trois termes : llnZ (b) = µn-1# (b-)+#n(b)+µn#(b+) (26) dans laquelle ju (y) représentant la métrique accumulée vers l'avant du meilleur sous-chemin partant du sous-état 0 e So, et se terminant au sous-état s e S,,, est calculé de manière récursive à l'aide de la formule suivante : <BR> <BR> <BR> <BR> Abn (s) = min {µn-1# (b-)+#n(b)}<BR> <BR> b#Bn-1,B+=8 avec les conditions limites suivantes :

µ0#(0)=0 et µ0#(s)=#,#s#0 (28) et dans laquelle, un (s) représentant la métrique accumulée vers l'arrière du meilleur sous-chemin partant du sous-état s E Sn, et se terminant au sous-état 0 e Sh est calculé de manière récursive à l'aide de la formule suivante : <BR> <BR> <BR> µn#(s) = min {µn+1#(b+)+#n+1(b)}<BR> <BR> b#Bn+1,b-=s (29) avec les conditions limites suivantes : f (0) = 0 et µ##(s)=#,#s#0 (30) La métrique de branche an (b) basée sur l'algorithme PSP et utilisée par le dispositif DDFSE 12 est exprimée de la manière suivante : Dans cette équation, le symbole complexe sn pénétrant dans le code d'interférence entre symboles au temps n résulte de la simple redéfinition du libellé de branche bV. La séquence complexe de symboles {sn-#r,...,sn-1} est déduite simplement du sous-état b-, tandis que la séquence estimée de symboles {#n-2L',...,#n-#v-1} est obtenue en parcourant en arrière le chemin"survivant"qui se termine en b-et en redéfinissant les libellés des branches le constituant. Les chemins"survivants"sont supposés être mémorisés dans une matrice glissante de remontée (traceback) de profondeur 2L'.

La courbe C4 sur la figure 5 qui illustre les performances de cette solution a été obtenue avec une complexité de treillis réduit DDFSE de quatre états (v, = 1).

Cette courbe montre que même avec un facteur d'entrelacement faible (= 4), et un réduction importante de du nombre d'états du treillis, cette solution permet de se rapprocher sensiblement de la solution idéale, et en particulier de la solution illustrée par la courbe C5, qui met en oeuvre une estimation de canal de type MLSE (Maximum Likelihood Sequence Estimation) qui est inapplicable à la technique CDMA du fait de sa complexité.

Pour améliorer davantage les performances d'un tel récepteur, l'invention propose comme représenté sur la figure 3, de relier la sortie du dispositif 12 à l'entrée positive d'un comparateur 13 dont la sortie est connectée au

désentrelaceur 14. La sortie du désentrelaceur 14 est connectée à une entrée recevant le logarithme de rapport de probabilité intrinsèque des bits codés du décodeur 15 qui fournit sur une première sortie le logarithme de rapports de probabilité a posteriori sur les bits de données transmis, et sur une seconde sortie le logarithme du rapport de probabilité a posteriori sur les bits codés qui est déterminé par une opération inverse au décodage portant sur les bits de données décodés après correction d'erreurs. La seconde sortie du décodeur 15 est rebouclée via un comparateur 18 et un entrelaceur 17 sur l'entrée négative du comparateur 13 et sur une entrée du dispositif DDFSE 12. Par ailleurs, la sortie du désentrelaceur 14 est appliquée également sur l'entrée négative du comparateur 18.

Le code de canal peut être décodé grâce au décodeur 15 de manière optimale en utilisant l'algorithme BCJR. Le dispositif 12 de réduction d'interférences délivre des rapports de probabilité a posteriori sur la valeur des bits an,j des symboles an constituant la séquence a1, à l'aide des logarithmes des rapports de probabilité a priori de ceux-ci provenant du décodeur (valant 0 à la première itération) et en tenant compte de la séquence reçue Y1 et d'une valeur estimée (ou réévaluée) h (n) du vecteur (équivalent) des coefficients de canal à l'instant n.

Les rapports de probabilité a posteriori approximés sur les bits X (ajj) peuvent être divisés en deux parties selon la formule suivante : #'(an,j)=#a(an,j)+#e(an,j) Après désentrelacement par le désentrelaceur 14, la séquence complète de logarithmes de rapports de probabilité extrinsèque devient une séquence de logarithmes de rapports de probabilité intrinsèque portant sur les bits des symboles codés, qui est appliquée au décodeur 15. D'une manière analogue, à la sortie du décodeur 15, chaque logarithme de rapport de probabilité a posteriori portant sur le bit codé S (cn, j) peut être décomposé en une partie a priori et une partie extrinsèque. Cette dernière peut être calculée en soustrayant bit à bit dans le comparateur 18, le logarithme du rapport a priori Ba (cnJ) à la sortie du décodeur, du logarithme correspondant du rapport a posteriori B (Cn ;) #e(cn,j)=#(cn,j)-#a(cn,j) (33) La séquence de logarithmes de rapports de probabilité extrinsèque portant sur les bits codés, en sortie du décodeur 15 est ré-entrelacée par l'entrelaceur 17 et renvoyée vers le décodeur 15 à la suite de la détection suivante de N séquences

de logarithmes de rapport de probabilité a priori portant sur les bits de symboles. En répétant cette procédure un certain nombre de fois, on obtient une augmentation importante du rapport signal sur bruit au niveau des bits de donnée des séquences reçues.

Du fait de sa structure Viterbi très régulière, et de ses bonnes performances au regard de sa complexité modérée, le détecteur DDFSE dont est muni le dispositif 12 apparaît parfaitement adapté à la turbo-détection telle qu'elle est effectuée dans le récepteur représenté sur la figure 3.

Pour la première itération, le fonctionnement du système représenté sur la figure 3 est identique à celui qui a été décrit en référence à la figure 2. A la seconde itération et aux suivantes, l'équation (31) donnant la métrique de branche sen (b) utilisée par le dispositif DDFSE 12 présente un terme de plus : Elle est calculée une seule fois durant le traitement récursif vers l'avant, puis mémorisée.

Le logarithme de la probabilité a priori In Pr (b=b) sur la branche b E Bn dans l'équation (31) correspond exactement au logarithme de la probabilité a priori du libellé bv qu'elle porte, si bien que : ln Pr (b = b) = In Pr (b'= b') = In Pr (a,,, = by) (35) Si l'on suppose une décorrélation parfaite entre les logarithmes des probabilités à priori sur les bits de symbole a,,, j après réentrelacement de la séquence de logarithmes de rapport de probabilités extrinsèques provenant du code Co, on obtient : Finalement, en utilisant les équations (25) et (36), la sortie ' (j) du dispositif DDFSE 12 appliquée au symbole anj peut être décomposée en une somme de deux termes logarithmiques.

'(anj) j)=#a(an,j)+#'e(an,j) où : Pr(an,j=1) #a(an,j)=ln<BR> Pr(an,j=0) représente le logarithme du rapport a priori portant sur le bit anj fourni par le décodeur 15, et où :

avec : le second terme de cette équation représentant le logarithme du rapport de probabilité extrinsèque sur le bit an, j apportée par tous les autres bits des symboles à libellés de bit de la séquence ai tout au long du processus de décodage.

Il convient d'insister sur le fait que dans le cas où vr= 2L', l'algorithme effectué par le dispositif DDFSE 12, devient formellement équivalent à un algorithme de type Min-Log-BCJR appliqué sur la totalité du treillis de canal. Si le traitement ne portait que sur un treillis à états réduits, les séquences estimées obtenues à partir de l'historique des chemins et impliquées dans les dérivées des métriques de branche, introduisent une dégradation de performance, due à un possible effet de propagation d'erreur. Toutefois, il apparaît que les canaux équivalents à la sortie du récepteur Rake 11 n'introduisent pas de propagation d'erreur significative dans la structure du dispositif de réduction DDFSE 12. Il en résulte que le choix de v, = 1 est suffisant dans la plupart des cas.

Il s'avère que les performances des estimations de canal conventionnelles qui sont obtenues par des traitements de corrélation et de calcul de moyenne, appliqués à une séquence de symboles pilotes, se dégradent avec des facteurs d'étalement faibles du fait des interférences entre symboles.

Pour illustrer ce fait, la figure 6 montre des courbes de taux d'erreur binaire en fonction du rapport signal sur bruit en sortie du décodeur 15, obtenues avec différentes solutions. Sur cette figure, la courbe C6 correspond au cas idéal. La courbe C7 qui correspond au cas où l'on utilise une estimation de canal conventionnelle, montre que cette solution présente des performances relativement mauvaises, assez éloignées du cas idéal.

La présente invention propose d'améliorer la qualité de l'estimation de canal

réalisée dans les systèmes représentés sur les figures 2 à 4, en utilisant la connaissance de la structure des interférences entre symboles. A cet effet, on utilise l'estimateur de canal 10 représenté sur la figure 7. Cet estimateur de canal comprend un estimateur de canal conventionnel 30 dont la sortie est reliée à un correcteur 31 d'estimation de canal utilisant la méthode MMSE (méthode de l'erreur quadratique moyenne minimum) ou LS (méthode des moindres carrés), délivrant une estimation de canal qui est utilisée par le dispositif 12 et par le dispositif 16 de modélisation de canal.

Un exemple d'estimateur de canal conventionnel est par exemple décrit dans la référence [1].

Si l'on suppose que les retards de chemins sont espacés d'un multiple de la période de chip Tc et que l'étalement des délais est inférieur à la période de symbole Ts, les estimations de canal conventionnelles sont obtenues par la formule suivante : <BR> <BR> <BR> <BR> <BR> h= =Mh+n (41)<BR> <BR> / dans laquelle : M Mji]0#j,i#L-1, (42) Mii=l, O<i<L-1, (43) p étant le nombre de symboles de la séquence pilote, h = (hl,..., étant les coefficients de canal parfaits (sans bruit), n étant l'estimation de bruit de canal qui est supposé avoir une variance égale à NolEpilot et Epilot étant l'énergie de

symboles pilotes. L'estimation de canal MMSE se déduit de l'estimation de canal conventionnelle par la formule suivante : hMMSE = LHh (46) <BR> <BR> dans laquelle<BR> L=argmin##MMSE-h#. (47 En utilisant l'équation (18), on peut en déduire : -MMSE = MH No-ih (48) Epilot dans laquelle M est la matrice transposée conjuguée de la matrice M.

Une estimation par la méthode LS peut également être effectuée : hLS = (MHM)-1 MHh. (49) Toutefois, cette méthode d'estimation ne tient pas compte de la puissance de bruit et ainsi dégrade les performances en termes de rapport de signal sur bruit, par rapport à l'estimation MMSE. Comme cela apparaît sur la figure 6, les simulations montrent qu'une estimation de canal classique, basée seulement sur la corrélation avec des symboles pilotes donne de mauvaises performances avec un facteur d'étalement faible (courbe C7). La structure des interférences entre symboles doit être prise en compte, ce qui est possible grâce aux méthodes MMSE et LS. Les performances obtenues par la méthode MMSE sont illustrées par la courbe C9 de la figure 6 qui montre une sensible amélioration par rapport aux méthodes conventionnelles.

Selon une variante avantageuse de l'invention représentée sur la figure 4, le récepteur selon l'invention comporte une boucle itérative de détection permettant d'effectuer une réestimation de canal. Cette boucle comprend un entrelaceur 19 connecté en sortie du décodeur 15 et dont la sortie est reliée par l'intermédiaire d'un comparateur à seuil 20 à un dispositif 21 d'estimation itérative de canal. Le comparateur 20 permet de transformer la sortie pondérée ou souple du décodeur 15 en une sortie dite"dure", c'est-à-dire égale à 0 ou 1 selon que la valeur pondérée est supérieure ou non à un seuil prédéterminé.

Par ailleurs, le signal appliqué en entrée du récepteur Rake 11 est également appliqué au dispositif 21 d'estimation itérative de canal et à un dispositif d'estimation de canal 10 utilisé pour effectuer une première estimation de canal

à l'aide de la séquence pilote, ces dispositifs fournissant des estimations de canal, par l'intermédiaire de commutateurs respectifs 23,24 au dispositif de modélisation de canal 16 conçu pour déterminer un modèle équivalent de canal qui est appliqué en entrée du dispositif DDFSE 12.

Ce récepteur effectue un premier décodage sur la base des estimations de canal déterminées par le dispositif 10 à l'aide des symboles pilotes, et appliquées au dispositif 16 de modélisation de canal (commutateurs 23,24 respectivement fermé et ouvert). Ensuite les codes estimés provenant du décodeur de canal 15 sont utilisés par le dispositif 21 pour ré-estimer le canal, ceci afin de réduire le bruit d'estimation de canal pour l'itération suivante, les estimations de canal ainsi déterminées étant appliquées au dispositif 16 de modélisation (commutateurs 23,24 respectivement ouvert et fermé).

Bien entendu, le processus itératif d'estimation de canal représenté sur la figure 4 peut également être appliqué au système représenté sur la figure 2, c'est-à- dire, à un récepteur ne comportant pas la boucle itérative de détection représentée sur la figure 3.

Les courbes C8 et C10 sur la figure 6 illustrent les performances obtenues grâce à la solution représentée sur la figure 4 utilisant une méthode d'estimation de canal itérative, la courbe C8 correspondant au cas où l'on utilise un estimateur de canal conventionnel et la courbe C10 un estimateur de canal avec correction MMSE. Ces deux courbes montrent que cette méthode permet d'améliorer les performances, et dans le cas d'une estimation de canal avec correction MMSE, de se rapprocher sensiblement du cas idéal illustré par la courbe C6.