FR 2721778 A1 19951229 FR 9407984 A 19940623 Procédé d'estimation d'une erreur de phase résiduelle sur les échantillons d'un signal numérique démodulé, et procédé de correction correspondant, Le domaine de l'invention est celui de la réception de signaux numériques, et notamment, mais non exclusivement, de signaux transmis à l'aide d'une pluralité de fréquences porteuses émises simultanément. En d'autres termes, I'invention s'applique en particulier aux récepteurs de signaux transmis selon la technique du multiplexage par répartition en fréquence (en anglais : Frequency Division Multiplex (FDM)), et par exemple aux récepteurs du système COFDM (Coded Orthogonal Frequency Division Multiplex (multiplexage de fréquences orthogonales codées)), mis en oeuvre notamment dans le cadre du projet européen Eurêka 147 1DABE (Digital Audio Broadcasting (diffusion audionumérique)). De façon classique, la réception d'un signal numérique consiste essentiellement à démoduler un signal transmis, à l'aide d'un oscillateur local de transposition de fréquence, puis à récupérer après filtrage et échantillonnage des valeurs numériques en phase et en quadrature. On effectue ensuite un décodage canal de ces valeurs numériques, en fonction du procédé de codage utilisé à l'émission. Si la transmission met en oeuvre plusieurs porteuses, une opération intermédiaire de récupération du symbole affecté à chaque porteuse est réalisée, par filtrage sélectif ou par l'application d'une transformation mathématique adéquate. On connaît de nombreux moyens de contrôle de la fréquence d'oscillation de l'oscillateur local, et notamment les dispositifs de contrôle automatique de fréquence (CAF). Un exemple de CAF, applicable aux signaux multiporteuses est par exemple décrit dans la demande de brevet française FR 93 02198, aux noms des mêmes déposants. Toutefois, de légères variations de fréquences peuvent apparaître dans l'oscillateur local, que ne peuvent combattre les techniques de CAF. Ces variations, appelées bruit de phase, se traduisent entre autres par un léger décalage de phase sur les échantillons à décoder. Ainsi que cela apparaîtra plus précisément par la suite, ce décalage de phase n'est pas pris en compte par les démodulateurs. Les échantillons fournis aux décodeurs de canal peuvent donc ne pas correspondre aux valeurs émises, et, par conséquent, perturber le décodage de canal. Pour limiter ce problème, le décodage de canal met en général en oeuvre des techniques de décodage à décision douce, telles que la méthode de Viterbi. Ces techniques permettent en général de reconstituer la séquence de données source, même si quelques valeurs démodulées sont fausses. Toutefois, dans certaines situations, le décodage de canal peut se trouver perturbé. C'est par exemple le cas si une série importante de valeurs fausses est reçue, du fait d'un bruit de phase trop élevé. Par ailleurs, plus le risque d'erreur dû au bruit de phase est élevé, plus le décodage de canal doit mettre en oeuvre des traitements complexes, donc coûteux, tant sur le plan du matériel que sur celui du temps de traitement. L'invention a notamment pour objectif de pallier ces différents inconvénients de l'état de la technique. Plus précisément, un objectif essentiel de l'invention est de fournir un procédé d'estimation de l'erreur de phase sur les échantillons d'un signal numérique, après leur démodulation. Un autre objectif de l'invention est également de fournir un tel procédé, permettant de limiter, voire de supprimer, le bruit de phase induit par l'oscillateur local. Encore un autre objectif de l'invention est de fournir un tel procédé, permettant de limiter les traitements de décodage source, notamment en ce qui concerne la correction des erreurs. Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints selon l'invention à l'aide d'un procédé d'estimation d'une erreur de phase résiduelle sur les échantillons de données d'un signal numérique démodulé, comprenant les étapes suivantes - décision dure de la valeur de chacun desdits échantillons, par association à chacun desdits échantillons de l'élément de données le plus probable, parmi une constellation d'éléments de données formant un alphabet de modulation - détermination d'un décalage de phase entre chacun desdits échantillons et son élément de données associé - estimation d'une information représentative de ladite erreur de phase résiduelle, à partir du décalage de phase correspondant à au moins deux échantillons (typiquement : 50 échantillons). En d'autres termes, I'invention permet de déterminer le décalage de phase subi par chaque échantillon démodulé, du fait des légers décalages de l'oscillateur local (bruit de phase), par une analyse "brute" de ces échantillons, avant le décodage proprement dit. Cela permet d'une part de compenser ce bruit de phase (ce que ne permet pas l'estimation de la réponse du canal de transmission), et d'autre part d'agir sur l'oscillateur local. De façon avantageuse, ladite étape d'estimation comprend une première étape de pondération desdits décalages de phase, en fonction d'un coefficient de confiance associé à chacun desdits échantillons, délivré par un module estimateur de canal (réponse du canal et/ou puissance du bruit). Ainsi, moins l'échantillon considéré a été perturbé lors de la transmission, plus la confiance accordée à l'estimation est grande. Le coefficient de confiance est par exemple le module au carré de l'estimation du canal utilisée pour la démodulation de l'échantillon. Préférentiellement, ladite étape d'estimation comprend également une seconde étape de pondération desdits décalages de phase, en fonction d'une information représentative de la puissance de bruit globale, délivré par un module d'évaluation du spectre de bruit. Cela permet d'augmenter sensiblement la valeur de l'estimation du décalage de phase, lorsque le niveau de bruit est important, de façon à compenser le fait que des erreurs sur les décisions dures induisent des estimations de décalage de phase ayant un signe erroné, et, en conséquence, une valeur absolue de l'erreur de phase résiduelle estimée trop faible. Selon un mode de réalisation avantageux de l'invention, ladite étape d'estimation comprend une étape de calcul de l'équation où: restim est ladite information représentative de l'erreur de phase résiduelle; a2 est la puissance de bruit globale f est une fonction croissante et telle que f(0) = ; est le décalage déterminé pour la porteuse k Pk est le coefficient de confiance associé à la porteuse k (correspondant sensiblement au rapport signal/bruit). Préférentiellement, de façon notamment à limiter les calculs dans ladite étape de détermination d'un décalage de phase, ledit décalage de phase est assimilé à sa partie imaginaire, sa partie réelle étant systématiquement fixée égale à 1. On vérifie en effet aisément que seule cette partie imaginaire est réellement significative, la partie réelle exacte étant toujours proche de 1. De façon avantageuse, ladite étape d'estimation comprend une étape de correction, consistant à sensiblement accentuer ladite erreur de phase résiduelle, I'accentuation étant croissante avec ladite erreur de phase. Cette étape de correction a pour but de compenser les sous-estimations du décalage de phase qui apparaissent lorsque ce décalage est important. Dans le cas où l'invention s'applique à un signal constitué d'une pluralité de porteuses émises simultanément, ladite estimation prend préférentiellement en compte au moins deux échantillons correspondant à des porteuses différentes. Comme on le verra par la suite, l'invention est particulièrement adaptée aux signaux de ce type, et par exemple aux signaux COFDM. De façon préférentielle, ladite étape d'estimation comprend une troisième étape de pondération, consistant à accorder une confiance sensiblement inversement proportionnelle à l'éloignement de l'élément de données considéré par rapport au centre de la constellation. En effet, plus l'élément de données est éloigné du centre, plus la différence de phase qui le sépare de son voisin est faible, ce qui peut conduire, dans les situations les plus défavorables, à de fausses décisions dures, et donc à de fausses estimations du déphasage. En revanche, à proximité du centre de la constellation, ce risque est beaucoup plus limité. Un mode de pondération simplifié peut consister à ne prendre en compte que des symboles proches du centre de la constellation. L'estimation de l'erreur résiduelle selon l'invention peut être utilisé de nombreuses manières, à l'intérieur de la chaîne de réception, et en particulier pour au moins une des opérations appartenant au groupe comprenant: - correction de l'effet de ladite erreur de phase résiduelle sur lesdits échantillons - contrôle d'un oscillateur local de démodulation du signal reçu; - correction des signaux démodulés en phase et en quadrature correspondant au signal reçu; - contrôle de traitements de décision douce sur la valeur desdits échantillons. Ainsi, en particulier, I'invention s'applique à un procédé de correction de l'effet d'une erreur de phase résiduelle sur les échantillons de données d'un signal numérique démodulé, comprenant les étapes suivantes - obtention d'une information (estim représentative de ladite erreur de phase résiduelle, par le procédé d'estimation d'une erreur de phase résiduelle selon le procédé décrit précédemment ; et - rotation desdits échantillons de eicestiml délivrant des échantillons corrigés. Avantageusement, lesdites étapes d'obtention et de rotation sont répétées au moins une fois sur lesdits échantillons corrigés, de façon à optimiser la correction. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante de plusieurs modes de réalisation de l'invention, donnés à titre d'exemples illustratifs et non limitatifs, et des dessins annexés dans lesquels - la figure 1 est un schéma synoptique simplifié d'un récepteur mettant en oeuvre le procédé de l'invention; - la figure 2 illustre les états que peuvent prendre les symboles d'une modulation à 16 états, et l'effet du décalage de phase; - la figure 3 présente, sous la forme d'un organigramme simplifié, le procédé d'estimation de l'invention, tel que mis en oeuvre dans la figure 1;; - la figure 4 présente un exemple d'une fonction f de pondération en fonction de la puissance du bruit, mise en oeuvre dans le procédé de la figure 3 - la figure 5 illustre un exemple d'une fonction g de correction de l'effet des forts déphasages - la figure 6 illustre un procédé de correction d'une erreur de phase résiduelle, mettant en oeuvre le procédé de la figure 3, et mis en oeuvre dans le récepteur de la figure 1. Comme indiqué précédemment, I'invention s'applique à la détection, notamment en vu de sa correction, d'une erreur de phase résiduelle apparaissant lors de la démodulation d'un signal numérique, que ce signal soit transmis selon une technique monoporteuse ou multiporteuse, quels que soient les types de modulation (par exemple MDP4, MDP8, 16QAM, 64QAM,...) et de démodulation (différentielle ou cohérente) mis en oeuvre. Le mode de réalisation préférentiel décrit par la suite s'applique au système COFDM. L'homme du métier saura aisément en déduire l'application à d'autres systèmes. On rappelle rapidement que le système COFDM, décrit en particulier dans le brevet français FR-86 09622 déposé le 2 juillet 1986 et dans le document "Principes de modulation et de codage canal en radiodiffusion numérique vers les mobiles" (M. Alard et R. Lassalle ; Revue de l'U.E.R, n" 224, août 1987, pp. 16 & 90), repose notamment sur la combinaison de deux principes essentiels. Le premier principe est la répartition de l'information à transmettre sur un grand nombre de porteuses (par exemple 512 ou 1024, sur une bande de 8 MHz) modulées chacune à un faible débit, afin de réduire l'effet de sélectivité du canal de transmission. Le second principe du système COFDM consiste à corréler par un procédé de codage des éléments d'information consécutifs et à les transmettre en des points distants du domaine temps-fréquence (technique de l'entrelacement en temps et en fréquence). L'éloignement de ces points est choisi de façon à assurer l'indépendance statistique des perturbations qu'ils risquent de subir lors de la transmission. Le codage employé est par exemple de type convolutif. Il est avantageusement associé à un décodage mettant en oeuvre un algorithme de Viterbi en décision douce. Plus précisément, le codeur convolutif délivre des valeurs complexes Ck appartenant à un alphabet de modulation. Le choix de l'alphabet spécifie le type de modulation utilisé. Par exemple, pour une modulation à 4 états de phase (MDP4), l'alphabet utilisé est (1 + i, 1 - i, -1 + i, -1 - i}. De nombreux autres types de modulation peuvent être utilisés, tels que les modulations MDP8, 16QAM, 64QAM ou les modulations par codage en treillis selon la méthode d'Ungerboeck. Chaque porteuse de chaque symbole est modulée par l'une de ces valeurs complexes Ck, k représentant l'indice de la porteuse. On sait que le canal de transmission induit des distorsions, qui correspondent à une multiplication de chaque valeur Ck par un complexe Hk, représentant la réponse du canal pour la keme porteuse. Par ailleurs, il apparaît toujours du bruit lors de la transmission, qui s'ajoute au signal transmis. Si l'on désigne par Nk le bruit affectant la porteuse k, la valeur complexe modulant cette porteuse à la réception peut s'écrire Yk = Hk.Ck + Nk Plusieurs méthodes d'estimation de Hk sont connues, qui délivrent une valeur estimée Hkest.Par exemple, dans le cas de la démodulation différentielle, on détermine au rang j (où j représente la dimension temporelle) un estimateur simplifié du canal déduit du rang j-1 Ho zest = Yj,k/Cj,k Nk est une variable gaussienne aléatoire, dont la variance ak (la puissance de bruit) peut varier d'une porteuse à l'autre, du fait des interférences. Cette variance est déterminée de façon connue, par un module d'estimation du bruit. Ainsi, le décodeur de canal est alimenté par des métriques déterminées à partir des échantillons corrigées Zk suivants Zk = Yk/Hkest avec une confiance associée correspondant sensiblement à La figure 1 est un schéma synoptique simplifié d'un récepteur mettant en oeuvre le procédé de l'invention. Le signal reçu 1 1 est tout d'abord transposé en bande de base, par le multiplieur 12 contrôlé par l'oscillateur local 13, qui délivre la fréquence de transposition fo. Puis le signal transposé est filtré par un filtre passe-bas 14, et converti en un signal numérique par un échantillonneur 15. Ensuite, un module 16 de génération des composantes en phase et en quadrature délivre les voies 17, et 17Q. Dans le cas du COFDM, les voies I et Q sont soumises à une transformation F.F.T. 18, qui délivre les échantillons 19 correspondant à chaque porteuse Yk. Un module 110 d'estimation de la réponse du canal fournit l'estimation 111 du canal, HkeSt, qui divise (112) les échantillons 19, pour délivrer les décisions brutes 113 sur la valeur complexe de l'échantillon reçu, soit sensiblement Zk = Ck + Nk/Hk Les éléments décrits jusqu'ici sont classiques d'une chaîne de démodulation COFDM. Ensuite, selon les techniques connues, les échantillons bruités 113 sont directement soumises à un module 125 de décodage de canal. Pourtant, ces échantillons peuvent être perturbées par un décalage de phase résiduel, ainsi que cela est illustré par la figure 2. Cette figure 2 illustre les 16 états 21l à 2116 que peuvent prendre les symboles d'une modulation à 16 états (16QAM par exemple) dans l'espace des complexes défini par les deux axes 22I et 22Q. Dans la pratique, les symboles reçus ne correspondent pas toujours exactement à ces symboles originaux 211 à 2116. En effet, de légères variations de fréquence de l'oscillateur local 13, appelées bruit de phase, peuvent provoquer une rotation d'un angle de déphasage 4 pour toutes les porteuses, qui n'est pas corrigée par l'estimation du canal Hkest Dans ce cas les échantillons reçus bruités Zk 23l à 2316 (représentés par des nuages de façon à symboliser le bruit) subissent une rotation de (, ou, en d'autres termes, sont multipliés par ejc. On constate aisément sur cette figure 2 qu'une telle situation peut conduire à des erreurs de démodulation. Ainsi, par exemple, l'élément 24, appartenant au nuage 232, correspond au symbole 211. Pourtant, il est plus proche du symbole original 21l, auquel il sera assimilé dans le cas des méthodes classiques. La sensibilité au phénomène augmente avec la taille de l'alphabet. L'invention permet d'éviter de telles erreurs de décodage. Pour cela, le récepteur (figure 1) comprend un module supplémentaire 114 d'estimation de l'erreur résiduelle 115 (, alimenté par les décisions dures Zk. Le procédé mis en oeuvre par ce module est décrit par la suite, en relation avec la figure 3. Avantageusement, ce module 114 tient compte d'un coefficient de confiance 119 pk2, délivré par le module 110 d'estimation du canal, et d'une information ak2 de puissance globale de bruit délivré par un module classique 121 d'estimation du bruit. L'erreur résiduelle ( 115 peut être utilisée de plusieurs façons, simultanément ou non. Notamment, il est possible de corriger l'erreur résiduelle sur les échantillons bruts 113, en les multipliant (122) par e-jC 123, après les avoir soumis à un retard 124 compensant le temps de calcul du module 114. Il s'agit d'une caractéristique essentielle de l'invention. Ainsi, le décodage de canal 125 est effectué sur des échantillons 126 non soumis, ou faiblement, à l'erreur de phase résiduelle. Comme on le verra par la suite, le traitement peut être réitéré, pour limiter encore l'effet de cette erreur de phase résiduelle. Avantageusement, le décodage canal 125 tient compte du rapport 126 a2k/ P2k Par ailleurs, l'erreur résiduelle peut servir à agir sur la démodulation des symboles suivants, après un filtrage passe-bas 116. Cela peut être fait selon deux manières alternatives - en agissant directement (117) sur l'oscillateur local 13 de façon à corriger le décalage de phase ; ou - en multipliant (118) les voies en quadrature 17I et 17Q par de. La figure 3 présente, sous la forme d'un organigramme simplifié, le procédé d'estimation de l'invention, tel que mis en oeuvre dans le module 114. Ce procédé comprend tout d'abord une étape 31 de décision dure, consistant à associer à l'échantillon Zk la valeur Pk la plus probable de l'alphabet de modulation utilisé, c'est-à-dire, en référence à la figure 2, la valeur 21i la plus proche. Ensuite, on détermine (32) le décalage de phase Ck entre l'échantillon Zk et la valeur associée Pk, à savoir: Pour simplifier les calculs, le décalage de phase peut être limité à sa partie imaginaire, ce qui revient à déterminer Après une étape 33 facultative d'optimisation et de correction, détaillée par la suite, on réalise (34) une moyenne sur plusieurs décalages de phase correspondant à des porteuses distinctes d'un même symbole (cas d'un signal multiporteuse) et/ou à plusieurs symboles successifs dans le temps, pour obtenir l'information 35 représentative de l'erreur résiduelle. L'étape 33 de correction peut comprendre quatre sous-étapes, qui peuvent être mis en oeuvre ou non indépendamment en fonction des besoins, du type du signal, du niveau de perturbation du canal, de l'application (réception dans les mobiles par exemple), etc... Le traitement peut donc comprendre une première sous-étape 36 de pondération en fonction d'un coefficient de confiance Pk2 délivré par l'estimateur de canal. En effet, il est souhaitable d'attacher une plus grande importance aux valeurs pour lesquels le canal n'est pas trop perturbé. Une seconde sous-étape 37 de pondération consiste à tenir compte de la puissance du bruit ak2, fournie par un module d'estimation du spectre de bruit. Cette pondération consiste à pondérer par P2k = I Hk l 2/a2k Par ailleurs, ainsi qu'on le constate aisément sur la figure 2, les symboles proches du centre de la constellation sont moins sujet à une erreur de décision que les symboles éloignés, à déphasage égal. En conséquence, une étape d'optimisation 38 complémentaire consiste à accorder une confiance supérieure à ces symboles proches du centre. Un mode de réalisation simple de cette pondération peut être de ne prendre en compte que les symboles proche du centre. Enfin, de même que le bruit, les forts déphasages introduisent des erreurs dans les décisions dures et donc des sous-estimations du déphasage. Cela peut être compensé, après moyennage 34, par une étape 39 de correction de biais. Cette étape 39 peut par exemple consister à associer au décalage ( la valeur g(C), la fonction g augmentant les décalages élevés (en valeur absolue). La figure 5 illustre une telle fonction g, dans le cas d'une modulation 64QAM. On constate qu'un déphasage calculé de 5 se voit associé à la valeur 7,5". On peut également tenir compte du bruit global a2, en multipliant le décalage de phase calculé par une fonction f croissante et telle que f(0) = 1. La figure 4 présente un exemple d'une telle fonction, variant linéairement de 1à 1,5 quand ak2 varie de 0 à ok2max. Cela a pour effet d'augmenter la valeur de (k quand le bruit est important, et donc de compenser le fait que des erreurs sur les décisions dures (dues au bruit) entraîne l'apparition de décalages estimés ayant un signe opposé au décalage réel, et donc une réduction de la valeur moyenne. Dans la pratique, on peut donc calculer directement la valeur suivante Comme indiqué précédemment, l'invention permet de corriger le décalage de phase calculé selon cette méthode. Un tel procédé de correction est illustré en figure 6. On calcule (61) tout d'abord l'erreur de phase résiduelle (, selon le procédé de la figure 3, puis on effectue une rotation de correction 62 correspondante aux décisions dures. Cela revient à multiplier les valeurs Zk par eIC. Avantageusement, ce traitement peut être réitéré (63) au moins une fois. ( étant plus petit la seconde fois, il apparaît moins d'erreurs sur les décisions dures. En conséquence, l'estimation sera meilleure. Dans ce cas, les étapes de correction peuvent être adaptées à chaque itération. En particulier, lors d'une seconde passe, il sera préférable de prendre en compte toutes les porteuses (suppression de l'étape 38) de façon à limiter l'effet du bruit. REVENDICATIONS 1. Procédé d'estimation d'une erreur de phase résiduelle sur les échantillons (Zk) de données d'un signal numérique démodulé, caractérisé en ce qu'il comprend les étapes suivantes - décision dure (31) de la valeur de chacun desdits échantillons, par association à chacun desdits échantillons de l'élément de données le plus probable, parmi une constellation d'éléments de données formant un alphabet de modulation - détermination (32) d'un décalage de phase entre chacun desdits échantillons et son élément de données associé - estimation (33, 34) d'une information représentative de ladite erreur de phase résiduelle, à partir du décalage de phase correspondant à au moins deux échantillons. 2. Procédé selon la revendication 1, caractérisé en ce que ladite étape d'estimation (33) comprend une première étape de pondération (36) desdits décalages de phase, en fonction d'un coefficient de confiance associé à chacun desdits échantillons, délivré par un module estimateur de canal. 3. Procédé selon la revendication 2, caractérisé en ce que ladite étape d'estimation (33) comprend une étape de calcul de l'équation où: (estim est ladite information représentative de l'erreur de phase résiduelle a2 est la puissance de bruit globale f est une fonction croissante et telle que f(0) = ; est le décalage déterminé pour la porteuse k Pk est le coefficient de confiance associé à la porteuse k. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que, dans ladite étape de détermination (32) d'un décalage de phase, ledit décalage de phase est assimilé à sa partie imaginaire, sa partie réelle étant systématiquement fixée égale à 1. S. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ladite étape d'estimation (33) comprend une étape de correction (37), consistant à sensiblement accentuer ladite erreur de phase résiduelle, l'accentuation étant croissante avec ladite erreur de phase. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il est mis en oeuvre sur un signal constitué d'une pluralité de porteuses émises simultanément, et en ce que ladite estimation prend en compte au moins deux échantillons correspondant à des porteuses différentes. 7. Procédé selon la revendication 6, caractérisé en ce que ladite étape d'estimation (33) comprend une troisième étape de pondération (38), consistant à accorder une confiance sensiblement inversement proportionnelle à l'éloignement de l'élément de données considéré par rapport au centre de la constellation. 8. Application du procédé selon l'une quelconque des revendications 1 à 7 à au moins une des opérations appartenant au groupe comprenant: - correction de l'effet de ladite erreur de phase résiduelle sur lesdits échantillons - contrôle d'un oscillateur local de démodulation du signal reçu - correction des signaux démodulés en phase et en quadrature correspondant au signal reçu - contrôle de traitements de décision douce sur la valeur desdits échantillons. 9. Procédé de correction de l'effet d'une erreur de phase résiduelle sur les échantillons de données d'un signal numérique démodulé, caractérisé en ce qu'il comprend les étapes suivantes - obtention (61) d'une information esnm représentative de ladite erreur de phase résiduelle, par le procédé d'estimation d'une erreur de phase résiduelle selon l'une quelconque des revendications 1 à 7; - rotation (62) desdits échantillons de e-j#estim, délivrant des échantillons corrigés. 10. Procédé selon la revendication 9, caractérisé en ce que lesdites étapes d'obtention (61) et de rotation (62) sont répétées (63) au moins une fois sur lesdits échantillons corrigés.