FR 2846165 A1 20040423 FR 0212900 A 20021017 Procédé d'estimation de la phase dans un système de communication numérique et boucle à verrouillage de phase Domaine technique de l'invention La présente invention concerne le domaine des communications numériques io et notamment un procédé d'estimation de la phase dans un système de communication numérique et une boucle à verrouillage de phase. Etat de la technique Dans un système de communication numérique comportant un émetteur et un récepteur, le signal numérique à transmettre - en général une suite de symboles- est converti avant transmission sous la forme d'un signal analogique à 20 temps continu qui est alors transmis par un milieu de propagation physique, l'air ou tout autre milieu physique de propagation. Lorsque le signal est reçu par le récepteur, il est alors traité et converti sous une forme numérique au moyen d'un échantillonnage approprié, lequel intervient en général à une fréquence fe que l'on cherche en général à rendre synchrone par rapport à la fréquence f, à l'émission 25 des symboles. Malheureusement les horloges se trouvant dans les circuits oscillateurs équipant l'émetteur et le récepteur ne sont jamais synchrones et il faut alors compenser toute dérive en fréquence entre ces horloges de manière à pouvoir traiter convenablement le signal reçu et en extraire les symboles émis. Cette dérive en fréquence affectant les oscillateurs à l'émission et à la réception engendre un déphasage parasite sur les observations du signal en sortie du démodulateur complexe situé dans le récepteur. D'autres facteurs contribuent à accentuer ce déphasage parasite. En premier lieu, on peut citer le temps que prend ST - 02-GR2-004 -2 le signal numérique pour traverser le milieu de propagation. En second lieu tout mouvement de l'émetteur par rapport au récepteur est source d'effet Doppler et tend, lui également, à introduire un déphasage perturbateur. En se référant au modèle en bande de base, les observations Yk en sortie du démodulateur complexe situé dans le récepteur peuvent être exprimées par la formule suivante: sur: Yk = ak e + Il k Avec ak correspondant aux symboles émis, k étant le déphasage parasite et llk étant un bruit supplémentaire. Des techniques sont déjà connues - basées sur des circuits estimateurs de 1 5 phase - pour estimer ce déphasage parasite ik et permettre sa correction. Les estimateurs de phase les plus sophistiqués, qui traitent simultanément des pans entiers d'observations reçues s'appuient sur des traitements numériques extrêmement lourds: le filtrage particulaire, les méthodes de Monte-Carlo par 20 chaîne de Markov etc... En pratique, de telles techniques s'avèrent irréalisables à mettre en oeuvre sur un plan pratique en raison de la puissance de calcul considérable devant être mises en jeu. Pour cette raison, on oppose à la complexité de ces estimateurs de phase la 25 simplicité de mise en oeuvre des boucles à verrouillage de phase qui traitent séquentiellement les observations reçues les unes après les autres et non par pans entiers d'observations reçues. Classiquement une boucle à verrouillage de phase est basée sur un algorithme numérique itératif permettant d'approcher une valeur estimée de la phase. Si traditionnellement, les boucles à verrouillage de phase sont 30 réalisées au moyen de circuits analogiques, le même traitement est aujourd'hui entièrement numérique Il est à noter que ce traitement numérique dépend étroitement du type de modulation envisagée. ST - 02-GR2-004 -3 A titre d'exemple, considérons le cas d'une modulation à déplacement de phase à deux états (MDP2), connue dans la littérature anglo-saxonne sous l'appellation Binaîy Phase Shift Keying BPSK. Dans cette modulation MDP2 ou BPSK, les symboles ak qui sont transmis sont égaux à -1 ou +1. Du fait du 5 déphasage parasite évoqué précédemment, on n'obtient plus -1 et +1 en sortie du démodulateur complexe, mais ces mêmes valeurs affectées d'un déphasage. Une boucle à verrouillage de phase bien connue pour prendre en charge la correction de ce déphasage est la boucle dite de COSTAS qui prend appui sur l'utilisation d'un algorithme du gradient associé à une fonction de cot J donnée par la formule: J(p) = E (|y2ke i2( - 2 Dans laquelle E étant l'opérateur Espérance mathématique. L'application de l'algorithme du gradient à la variable (p permet de faire converger l'algorithme vers la phase estimée: (Pk =( k-1 - 7 aJ ( P) IP 1 =D k-1 La boucle de Costas est finalement obtenue en supprimant l'espérance mathématique: c'est l'algorithme du gradient stochastique pour minimiser la fonction de cot J. On connaît d'autres formules pour d'autres types de modulation et 25 notamment la modulation en quadrature, également dénommée Quadrature Amplitude Modulation à quatre états ou 4QAM. D'une manière générale, quelle que soit le type de modulation employé, les boucles de verrouillages de phase sont construites sur une formule générale du type: ST - 02-GR2-004 -4 Typiquement, comme on le voit sur la formule précédente, toutes les boucles consistent à calculer la phase (P k en fonction de l'élément précédent (P k-i et d'une fonction F des deux éléments Yk et p k-i. Il est à noter que, dans cette formule, en fonction du type de déphasage à corriger, on peut avoir recours à une correction 5 sophistiquée au niveau du paramètre y et notamment prévoir la mise en oeuvre d'un filtre correcteur du second ordre (proportionnel intégral), voire d'un ordre supérieur. Toutes les boucles à verrouillage de phase connues - reprenant classiquement le modèle analogique traditionnellement connu - présentent une io même limitation. L'évaluation de la phase (9 k est essentiellement basée à partir de la valeur précédente (à k-1 et à partir d'une fonction d'une ou plusieurs observations passées. De ce fait l'estimation de la phase reste imparfaite et, par suite, sa correction. Il convient d'améliorer le modèle de la boucle à verrouillage de phase afin d'accroître la précision de l'estimation et l'efficacité de la correction. Exposé de l'invention Un objet de l'invention est une nouvelle structure de boucle à verrouillage de phase permettant d'accroître la précision de l'estimation de la phase par rapport à une boucle classique.. Un autre objet de l'invention consiste à proposer un procédé d'estimation de phase pour un récepteur numérique parfaitement adapté au traitement d'un signal numérique dans un récepteur doté d'un système de correction d'erreurs. Selon l'invention, ces objets sont réalisés au moyen d'un procédé d'estimation de phase et d'un dispositif de boucle à verrouillage de phase. L'invention propose un Procédé d'estimation d'une phase dans un système de communication numérique comportant les étapes suivantes: ST - 02-GR2-004 - réception et mémorisation d'un bloc d'observations Yk; - exécution d'au moins une boucle à verrouillage de phase sur une séquence prédéterminée d'observations extraites dudit bloc. L'invention propose également un Dispositif de boucle à verrouillage de phase pour un récepteur numérique comportant: - des moyens pour recevoir et stocker des observations par bloc; - une première boucle de verrouillage de phase permettant de générer une première valeur intermédiaire; - une seconde boucle de verrouillage de phase permettant de générer une seconde valeur intermédiaire; - des moyens pour dériver desdites première et seconde valeur intermédiaire une valeur estimée de la phase. En réalisant un traitement par bloc au niveau de la boucle à verrouillage de phase - lequel traitement par bloc n'intervient classiquement qu'en aval, au niveau des codes correcteurs d'erreur, on améliore sensiblement l'estimation de la phase tout en conservant la grande simplicité de mise en oeuvre des boucles 20 traditionnelles à verrouillage de phase. On peut en particulier briser le lien chronologique entre les observations et l'algorithme itératif qui est mis en oeuvre par la boucle à verrouillage de phase. L'invention consiste à rapprocher deux concepts a priori antinomiques, à savoir le 25 traitement par bloc et la boucle de phase. Dans l'approche classique, on conçoit la boucle de phase comme un traitement itératif, prenant appui sur les observations passées pour en dériver une estimation de la phase à un instant donné. La présente invention transgresse radicalement cette approche et réalise un procédé dans lequel, une fois un bloc d'observations mémorisées, l'estimation de la phase peut 30 être effectuée à partir d'une ou plusieurs boucles de phase opérant sur toute séquence quelconque d'observations au sein du même bloc. De préférence, on réalise une première boucle à verrouillage de phase sur un premier ensemble d'observations dudit bloc et une seconde boucle à verrouillage de ST - 02-GR2-004 -6 phase sur une autre séquence d'observations extraites du même bloc. Le résultat de chacune des boucles est combiné pour fournir une estimation affinée de la phase à chaque instant d'échantillonnage. Comme on le voit, en combinant le traitement par bloc et la simplicité des boucles à verrouillage de phase, on peut sensiblement améliorer la précision de l'estimation de phase et, par suite, la correction. En effet, le traitement par bloc permet en effet de tenir compte, pour le calcul de 9 k non seulement des observations antérieures mais également des observations postérieures au sein i0 d'un même bloc. La précision de la boucle de verrouillage de phase en est alors considérablement améliorée. Bien que le procédé introduise un retard de traitement lié à la mémorisation d'un bloc complet d'observations avant l'estimation de la phase pour chacune des 1 5 observations Yk considérées, ce retard n'est pas préjudiciable et est parfaitement adapté à l'usage de certains codes correcteurs d'erreurs, comme les turbos codes ou les codes en blocs. Préférentiellement le traitement par bloc sera effectué au moyen d'une 20 première boucle de phase opérant dans le sens chronologique des observations, et d'une seconde boucle de phase opérant en sens inverse. Dans un mode de réalisation préféré, la seconde boucle à verrouillage de phase est initialisée avec une valeur qui est celle fournie par la première boucle en 25 fin d'itération. Description des dessins D'autres caractéristiques, buts et avantages de l'invention apparaîtront à la lecture de la description et des dessins ci-après, donnés uniquement à titre d'exemples non limitatifs. Sur les dessins annexés: ST - 02-GR2-004 La figure 1 illustre le schéma de principe du procédé d'estimation d'une phase selon l'invention. La figure 2 illustre l'effet d'un déphasage pour une modulation BPSK. La figure 3 illustre le mode de réalisation préféré basé sur l'emploi de deux boucles à verrouillage de phase opérant en sens contraire. La figure 4 illustre l'application du procédé de la figure 1 à une modulation 4i0 PSK ou 4-QAM. La figure 5 illustre l'application de l'invention à une modulation de type BPSK à deux états. Description d'un mode de réalisation préféré En relation avec la figure 1, on décrit à présent le procédé général 20 d'estimation de la phase, lequel convient pour tout type de modulation, qu'il s'agisse d'une modulation BPSK, 4-QAM etc.... Le procédé comporte une première étape 11 au cours de laquelle un ensemble de n+1 observations ( par hypothèse Y0 à Yn) en sortie du démodulateur complexe est réceptionné et mémorisé au sein d'un bloc de n+1 observations. Typiquement, on pourra envisager un bloc de 1000 à 10 000 25 observations et l'homme du métier choisira un nombre en fonction de l'application et du type de modulation envisagée. Une fois la mémorisation du bloc effectuée, on réalise une première boucle à verrouillage de phase qui va opérer sur une séquence prédéterminée 30 d'observations du bloc. Typiquement, on pourra considérer toute séquence quelconque et plus particulièrement une séquence chronologique d'observations Yk, suivant l'ordre de réception des données. La réalisation concrète d'une boucle à verrouillage de phase ne pose aucune difficulté en soi et, pour des raisons de simplicité d'exposé, celle-ci ne sera pas développée plus avant. On se bornera à ST - 02-GR2-004 dire que, en fonction du type de modulation utilisé, on mettra en oeuvre un algorithme itératif suivant la formule générale rappelée précédemment. P k = k-1 - ' F ( Yk, 9 k-1) avec k = 1 à n Dans une étape 12, on procède à l'initialisation de la première boucle à verrouillage de phase. Typiquement, cette initialisation permet de fixer les premières valeurs de l'algorithme itératif, et notamment la première valeur p o. De préférence, io lorsque la phase présente une continuité d'un bloc à l'autre, on initialisera d'une manière avantageuse la première boucle en prenant en considération la dernière valeur estimée au bloc précédent. Pour une boucle d'ordre supérieur à 1, plusieurs paramètres devront être fixés et, à nouveau, on pourra avantageusement tirer profit de toute continuité exploitable. Dans une étape 13, on exécute la première boucle à verrouillage de phase de manière à construire la séquence 9p o. p 1, 9 2, p 3..., (p,, comme cela est illustré dans la figure 3 par la flèche allant de la gauche vers la droite. Dans une étape 14, on initialise à présent une seconde boucle à verrouillage de phase allant dans le sens opposé à la première (représentée sur la figure 2 par la flèche dirigée de la droite vers la gauche). De préférence, on initialise la première valeur de la seconde boucle, à savoir (p', avec la dernière valeur numérique 9 n calculée par la première boucle à verrouillage de phase. Dans une étape 15, on exécute la seconde boucle à verrouillage de phase de manière à construire la séquence 9' n-1, 9' n-2, 9' 2, p' 1...., 9' o calculée dans le sens inverse à précédemment. (p' k = (P k+l - ry F ( Yk,9'k+1) avec k = n-1 à O Dans une étape 16, on combine les résultats intermédiaires produits par la première et par la seconde boucle à verrouillage de phase de manière à produire la phase estimée 9p" k suivant la formule ST - 02-GR2-004 ( k = G (k, (pk) dans laquelle la fonction G est adaptée au type de modulation considéré. 5 Dans un mode de réalisation particulier, on choisira G de manière à générer une somme pondérée sous la forme: (P k = Ax wk + Bx p' k De préférence, on pourra d'ailleurs choisir des coefficients variables Ak et Bk de manière à privilégier l'une des deux boucles à verrouillage de phase en fonction de k. En effet, on pourra avantageusement choisir les " poids " de la somme pondérée de manière à privilégier la première boucle dans la partie droite du bloc de la figure 2 et, inversement, ajouter plus de poids à la seconde boucle dans la 1 5 partie la plus à gauche du bloc. Ainsi, on privilégie toujours la boucle ayant effectué le plus d'itérations et pourra prétendre à une plus grande précision dans le calcul de la phase. En relation avec la figure 4, on décrit à présent l'application du procédé 20 général de la figure 1 à une modulation de type 4-BPSK ou 4- QAM. Dans une modulation de ce type, on transmet les symboles complexes ak qui sont respectivement -1, +1, -i et +i. Le procédé débute par une étape 41 consistant, comme précédemment, en la réception et la mémorisation d'un bloc d'observations reçues du démodulateur complexe. Dans une étape 42, on procède à l'initialisation de la première boucle à verrouillage de phase adaptée pour une modulation 4-QAM ( k = (P k-1 + y Img (Y4k e -i4ç(k-1) 30 dans laquelle IMG correspond à la partie imaginaire du nombre complexe. ST - 02-GR2-004 - 10 Il est à noter que, en pratique, on pourra prévoir l'emploi d'un filtre numérique sophistiqué pour le calcul du facteur y en fonction de modèle du déphasage que l'on souhaitera corriger. Pour des déphasages simples, on pourra s'accommoder d'un correcteur proportionnel simple alors que dans des cas plus complexes, on 5 pourra avantageusement recourir à un correcteur intégral, voire à un filtre d'ordre supérieur. Pour l'initialisation de la boucle à verrouillage de phase, on donnera la première valeur de (po. De préférence, le filtre y pourra être réalisé au moyen d'un io filtre numérique d'un ordre 2, par exemple du second ordre suivant la formule ciaprès: 7 = Yi + Y2/(1 +Z-1) ]5 et que l'on initialisera en tenant compte de tous les facteurs de continuité exploitables. Dans une étape 43, on exécute la première boucle à verrouillage de phase de manière à construire toute la séquence (p o. (p 1, (P 2, (à 3 (p n Dans une étape 44, on initialise à présent une seconde boucle à verrouillage de phase de COSTAS avec la dernière valeur calculée par la première boucle et, dans une étape 45 on exécute cette boucle pour construire la séquence (P' n-, q' n-2 p2, (P' 1 ---, (p' o calculée dans le sens inverse à précédemment. 25 (Pk (P k+1 + y img (Y4k e -i49(k+l) k=n-1...0 Les deux résultats intermédiaires sont alors combinés dans une étape 46 de manière à générer l'estimation du déphasage (P" k de préférence suivant la formule 30 (P k = G (k, 9 k) Dans une étape 47, le procédé achève le traitement du bloc courant. ST - 02-GR2-004 - 1i1 Pour illustrer la généralité de l'invention, on va montrer à présent en référence à la figure 5 que l'invention peut contribuer a améliorer considérablement l'estimation de phase dans une modulation de type BPSK. Le procédé démarre, une fois encore, par la réception et la mémorisation d'un bloc de d'observations Yk dans une étape 51. Dans une étape 52, on procède à l'initialisation de la première boucle à verrouillage de phase adaptée pour une modulation BPSK. D'une manière particulièrement avantageuse et nouvelle, on combine cette première boucle à verrouillage de phase avec l'exploitation d'informations statistiques portant sur les symboles ak transmis dans cette modulation BPSK. On réalise à cet effet une 15 boucle à verrouillage de phase suivant la formule: (P k = (1 k-1 + 'Y Img (Yk e -ip(k-1) th[ Lk /2 + 2/a2 Re( Yk e ip(k)]) dans laquelle: th: étant l'opérateur tangeante hyperbolique, Re: étant l'opérateur renvoyant la partie réelle d'un nombre complexe, Ca2: étant la variance du bruit; et Lk = Ln [ p(ak =1) / p(ak =-)] et Ln étant le logarithme népérien. p(ak =1) est la probabilité pour que le symbole ak soit égal à +1 et p(ak =-1) est la probabilité pour que le symbole ak soit égal à -1, ST - 02-GR2-004 - 12 Comme précédemment, on pourra prévoir l'emploi d'un filtre numérique du second ordre - voire d'un ordre supérieur - pour le facteur y s Dans une étape 53, on exécute la première boucle à verrouillage de phase de manière à construire la séquence (p 0, 9p 1, (p 2, ( 3 -..., (P Dans une étape 54, on initialise à présent une seconde boucle à verrouillage de phase construite comme précédemment, en tenant compte d'informations i o statistiques portant sur les symboles. avec la dernière valeur calculée par la première boucle et, dans une étape 55 on exécute cette boucle pour construire la séquence (p' n-, (1' n-2, (p' 2 9' 1...., (P' o calculée dans le sens inverse à précédemment. (k = () k+1 + y Jmg (yk e -iq(k Il th[ Lk /2 + 2/2 Re( Yk e -i](k+) Deux variantes de réalisation sont encore possibles. Si la puissance de bruit reste faible, on peut approximer la fonction tangente hyperbolique par la fonction signe. On obtient alors la formule suivante typique d'une boucle à remodulation 20 (désignée sous l'appellation anglo-saxonne décision feeback loop) au terme additif Lk /2 près. On peut noter que, dans le cas contraire lorsque la puissance de bruit est 25 importante, l'opérateur tangente hyperbolique peut s'approximer par la fonction identité et l'on retombe alors sur le formule d'une boucle de COSTAS (au terme additif Lk /2 près). Dans une étape 56, on combine les deux résultats intermédiaires obtenues 30 des deux boucles à verrouillage de phase de manière à générer l'estimation du déphasage (9" k de préférence suivant la formule (13 k = G ( 9k ( 'k) ST - 02-GR2-004 - 13 Comme précédemment, le traitement du bloc courant s'achève dans une étape 57. Les deux exemples précédents montrent que le procédé suivant l'invention s'applique à tout type de modulation et à tout type de boucle à verrouillage de phase. Il est clair que l'homme du métier pourra réaliser toutes les adaptations appropriées, voire même combiner avantageusement le traitement par bloc de la 1o boucle de phase réalisé très en amont avec tous les traitements ultérieurs réalisés en aval, qu'il s'agisse de traitement de codes correcteurs d'erreur, turbocodes etc... ST - 02-GR2-004 - 14 Revendications 1. Procédé d'estimation d'une phase dans un système de communication numérique comportant les étapes suivantes: - réception et mémorisation d'un bloc d'observations Yk; - exécution d'au moins une boucle à verrouillage de phase sur une séquence 1o prédéterminée d'observations extraites dudit bloc. 2. Procédé d'estimation d'une phase dans un système de transmission numérique selon la revendication 1 caractérisé en ce qu'il comporte: exécution d'une première boucle de verrouillage de phase sur lesdites observations dans le but de générer une première valeur intermédiaire; exécution d'une seconde boucle de verrouillage de phase sur lesdites observations dans le but de générer une seconde valeur intermédiaire; - combinaison desdites première et seconde valeur intermédiaire pour générer une 20 valeur estimée de la phase. 3. Procédé selon la revendication 2 caractérisé en ce que ladite première boucle est exécutée sur la séquence d'observations dans l'ordre chronologique d'apparitions, et en ce que la seconde séquence est exécutée sur la séquence inverse. 4. Procédé selon la revendication 3 caractérisé en ce que la seconde boucle à verrouillage de phase est initialisée à partir de la dernière valeur calculée par ladite première boucle à verrouillage de phase. 5. Procédé selon la revendication 4 caractérisé en ce qu'il comporte les étapes suivantes: - réception et mémorisation d'un bloc d'observations Yk avec k variant de O à n; ST - 02-GR2-004 - 15 - initialisation d'une première boucle de phase à partir des observations Yk reçues; - exécution de ladite première boucle de phase suivant une formule: k = ( k-1 - 7 F ( Yk, q k-1) avec k = 1 à n avec F étant une fonction adaptée au type de modulation considérée initialisation d'une seconde boucle de phase à partir des observations Yk avec k variant de n à O; -exécution de ladite seconde boucle à verrouillage de phase suivant une formule: (P'k (P'k+1 - 7 F( Yk, 9'k+1) avec k = n-1 à O - combinaison des résultats produits par lesdites première et seconde boucle de manière à générer une valeur estimée de la phase. 6. Procédé selon l'une quelconque des revendications précédentes caractérisé en ce que la modulation est de type BPSK avec une boucle à verrouillage de phase définie par ( k = k- + Img (Yk e-i(k-) th[ Lk /2 + 2/G2 Re( Yk e-i(k-)]) dans laquelle: th: étant l'opérateur tangente hyperbolique, Re: étant l'opérateur renvoyant la partie réelle d'un nombre complexe, 2: étant la variance du bruit; et Lk= Ln [ p(ak=1) / p(ak =-1)] et Ln étant le logarithme népérien, p(ak =1) la probabilité pour que le symbole ak soit égale à +1 et p(ak =-1) la probabilité pour que le symbole ak soit égale à -1. ST - 02-GR2-004 - 16 7. Procédé selon la revendication 6 caractérisé en ce que ledit facteur y est réalisé au moyen d'un filtre numérique d'ordre 2 ou supérieur. 8. Dispositif de boucle à verrouillage de phase pour un récepteur numérique comportant: - des moyens pour recevoir et stocker des observations par bloc; - une première boucle de verrouillage de phase permettant de générer une première valeur intermédiaire; io - une seconde boucle de verrouillage de phase permettant de générer une seconde valeur intermédiaire; - des moyens pour dériver desdites première et seconde valeur intermédiaire une valeur estimée de la phase. 9. Dispositif de boucle à verrouillage de phase selon la revendication 8 caractérisé en ce que ladite première et seconde boucle à verrouillage de phase sont réalisées au moyen de la formule: 9 k = 9 k-1 - y F( Yk, p k-1) avec k = 1 à n ou (P k = 9ki - y F( Yk, 9k+1) avec k = n-1 à 0 avec F étant une fonction adaptée au type de modulation considérée 10. Dispositif selon la revendication 9 caractérisé en ce que première valeur 25 calculée par ladite seconde boucle est déterminée par le dernier calcul réalisé par la dite première boucle à verrouillage de phase. ST - 02-GR2-004