FR 2781626 A1 20000128 FR 9809578 A 19980723 La présente invention concerne un démodulateur dit QPSK destiné à reconstruire deux signaux binaires véhiculés par deux porteuses de même fréquence mais en quadrature de phase. La présente invention concerne plus spécifiquement un estimateur de l'erreur de la fréquence sur laquelle une boucle à verrouillage de phase du démodulateur est réglée. La figure 1 représente sous forme de constellation des couples de valeurs, ou symboles, correspondant aux valeurs suc- cessives des deux signaux binaires reconstruits par le démodulateur. Cette constellation représente des points dont chacun a pour coordonnées les valeurs I et Q des deux signaux binaires à un instant donné. Théoriquement, les points sont confondus avec des points P1 à P4 de coordonnées respectives (1, 1), (-1, 1), (-1, -1), et (1, -1). Toutefois, comme cela est représenté dans le premier quadrant de la constellation, des points qui devaient être confondus avec le point Pl forment un nuage de points autour de ce point Pl à cause de diverses erreurs et de bruit dans la transmission. Un démodulateur QPSK est généralement formé d'une bou- cle à verrouillage de phase (PLL) destinée à régler la fréquence d'un oscillateur local sur la fréquence de porteuse par une analyse des points de la constellation. Afin de filtrer le bruit, la PLL a une bande passante très faible, de l'ordre du millième de la fréquence de transmission, ou fréquence symbole, des signaux binaires. De ce fait, la PLL a aussi une plage de capture étroite, qui de plus diminue lorsque le bruit augmente. Cette plage est de l'ordre de 0,5 % de la fréquence symbole avec un rapport signal sur bruit de 3dB. La fréquence de porteuse est généralement imparfaitement connue et le décalage peut être de plusieurs fois la fréquence symbole. Ainsi, si on ne prend pas de mesures particulières, la PLL ne parvient à se régler sur la fréquence de porteuse. Lorsque la fréquence de la PLL est mal réglée, la cons- tellation de points reçus tourne, commroe cela est illustré par des flèches en figure 1, à une vitesse égale à l'erreur de fréquence. Une première solution classique pour trouver une fré- quence située dans la plage de capture de la PLL consiste à régler la fréquence de la PLL sur des valeurs successives choi- sies jusqu'à ce que l'on détecte un verrouillage de la PLL. Une détection de verrouillage de PLL est une opération relativement lente qui nécessite l'analyse d'un nombre élevé de symboles pour chaque fréquence essayée. Le nombre de symboles à analyser croît avec le bruit. Dans les transmissions par satellite, on tolère généralement un rapport signal sur bruit aussi faible que 3 dB. Dans ce cas, il faut analyser plusieurs dizaines de milliers de symboles pour détecter le verrouillage, ce qui ralentit considé- rablement la cadence des essais de fréquence successifs et donc la vitesse à laquelle le démodulateur est réglé sur la fréquence de porteuse. D'autres solutions utilisent un estimateur d'erreur de fréquence destiné à ramener la fréquence du démodulateur en une seule tentative à l'intérieur de sa plage de capture. Un exemple d'estimateur de fréquence est décrit dans l'article intitulé "Frequency Detectors for PLL Acquisition in Timing and Carrier Recovery" par David G. Messerschmitt, IEEE Transactions on Comtmunications, volume 27, N 9, septembre 1979. Cet estimateur exploite certaines propriétés du signal et fournit une valeur en principe proportionnelle à l'erreur de fréquence. Toutefois, ce type d'estimateur fournit la position du centre de gravité du spectre du signal et ne fonctionne donc convenablement que si ce spectre est symétrique. Dans la plupart des situations réelles, le spectre n'est pas toujours symétrique du fait des caractéristiques des éléments de transmission du signal, qui fluctuent par rapport aux caractéristiques théoriques souhaitées. Ainsi, la valeur de fréquence fournie par l'estima- teur fluctue en pratique de plusieurs centièmes de la fréquence symbole et tombe souvent hors de la page de capture de la PLL. Un détecteur d'erreur de fréquence est décrit dans l'article intitulé "New Phase and Frequency Detectors for Carrier Recovery in PSK and QAM Systems", par Hikmet Sari et Said Moridi, IEEE.R Transactions on Communications, volume 36, N 9, septembre 1988. Ce détecteur analyse l'évolution de l'erreur de phase de la PLL. En effet, dans le cas d'une erreur de fréquence, la phase évolue linéairement en fonction du temps, ceci modulo x/2 dans un démodulateur QPSK. En absence de bruit, la phase varie en dent de scie entre -x/4 et 7/4, et la fréquence correspond à la pente de chacune des rampes de la dent de scie. Le système pro- posé fournit l'erreur de phase mesurée lorsqu'elle est comprise entre un seuil positif et un seuil négatif, et fournit la der- nière valeur mesurée entre ces seuils sinon. Toutefois, en présence de bruit important, les points de chaque rampe sont superposés à des erreurs aléatoires pouvant souvent excéder impossible d'avoir un résultat fiable. Un objet de la présente invention est de prévoir un es- timateur d'erreur de fréquence d'un démodulateur QPSK, qui soit particulièrement fiable, même en présence de bruit important. Cet objet est atteint grâce à un procédé d'estimation de l'erreur de fréquence d'un démodulateur destiné à reconstruire deux signaux binaires véhiculés sur deux porteuses de même fré- quence mais en quadrature de phase, comprenant les étapes consistant à former des vecteurs ayant pour composantes les cou- ples successifs de valeurs des deux signaux binaires; à appliquer à chaque vecteur une transformation qui multiplie par quatre son angle au moins lorsqu'il est égal à un multiple des vecteurs transformés. Selon un mode de réalisation de la présente invention, la transformation est une rotation de quatre fois l'angle du vecteur. Selon un mode de réalisation de la présente invention, l'erreur de fréquence est fournie par le calcul de la dérivée de l'angle du vecteur moyen. Selon un mode de réalisation de la présente invention, la dérivée est fournie sous forme de différence, modulo x/2, de deux valeurs d'angle successives. Selon un mode de réalisation de la présente invention, la transformation est une transformation linéaire par morceaux qui fait correspondre à chaque vecteur un vecteur image situé sensiblement dans le même quadrant du plan image qu'un vecteur image obtenu par multiplication par quatre de l'angle du vecteur. Selon un mode de réalisation de la présente invention, on attribue une phase de O aux vecteurs image situés dans les deuxième et troisième quadrants du plan image, une phase +1 aux vecteurs image situés dans le premier quadrant, et une phase -1 aux vecteurs image situés dans le quatrième quadrant, une infor- mation d'erreur de fréquence étant fournie par la différence des valeurs de phase attribuées au vecteur image courant et au vec- teur image précédent, laquelle différence est fournie modula 2 si elle est positive et modulo -2 si elle est négative. Ces objets, caractéristiques et avantages, ainsi que d'autres de la présente invention seront exposés en détail dans la description suivante de modes de réalisation particuliers faite à titre non-limitatif en relation avec les figures jointes parmi lesquelles: la figure 1, précédemment décrite, représente une cons- tellation correspondant à des symboles reconstitués par un démodulateur QPSK; la figure 2 représente un schéma de principe d'un esti- mateur d'erreur de fréquence selon l'invention; la figure 3 illustre l'image d'un cercle de la constel- lation de la figure 1 par une transformation simplifiée permettant d'obtenir un résultat équivalent à celui que l'on obtiendrait en utilisant le circuit de la figure 2; et la figure 4 représente un exemple de filtre de moyenne utilisé dans le circuit de la figure 2. Afin d'estimer l'erreur de fréquence d'un démodulateur QPSK, on propose selon la présente invention, de calculer la phase moyenne des symboles et de dériver cette phase moyenne afin de trouver l'erreur de fréquence. On pourrait pour cela utiliser un détecteur de phase classique, en effectuant en outre un filtrage passe-bas sur les valeurs de phase. Ceci aurait pour conséquence un lissage du bruit permettant de retrouver une évolution proche d'une dent de scie d'amplitude n/2. Toutefois, en présence de bruit important, un grand nombre de valeurs de phase instantanées dépassent les limites d'amplitude de la dent de scie et engendrent des contributions erronées au calcul de moyenne. Par exemple, si la dent de scie varie entre -x/4 et x/4, une phase instantanée de valeur -x/4 - s est transformée en une phase de valeur + iz/4 - s, d'o il résulte que cette phase instantanée est comptabilisée avec une erreur de i/2 dans le calcul de moyenne. Dans un milieu fortement bruité, la probabilité pour avoir de telles valeurs hors limites est si élevée que la valeur moyenne est trop entachée d'erreurs pour être exploitable. Pour éviter ces inconvénients, on prévoit selon l'invention, de calculer la moyenne sur des symboles qui ont subi une rotation de quatre fois leur angle de départ tout en conser- vant sensiblement leur module. Avec une telle transformation, deux symboles quelcon- ques qui ont la même position modulo ir/2 dans la constellation d'origine ont la même position absolue dans la constellation image. En d'autres termes, toute erreur de t/2 que l'on pourrait faire en considérant les phases dans la constellation d'origine du fait qu'on ne sait pas lequel des quatre symboles (I, Q) a été émis, est supprimée dans la constellation image grâce à une telle transformation. En effet, tout symbole affecté d'une erreur de n/2 est confondu avec le symbole correct dans la constellation image. La figure 2 représente un schéma de principe d'un exem- ple de circuit mettant en oeuvre le procédé susmentionné. Les couples de valeurs I et Q des deux signaux binaires reconstruits par le démodulateur sont fournis à un convertisseur cartésien- polaire 10 qui calcule le module R du vecteur de coordonnées (I, Q) et l'angle 0 de ce vecteur par rapport à l'axe horizontal (I à la figure 1). L'angle O est multiplié par 4 en 12 avant d'être fourni, avec le module R, à un convertisseur polaire-cartésien 14 qui calcule un couple de coordonnées (x, y) correspondant. Les moyennes xl et yl des évolutions des coordonnées x et y sont fournies par des filtres passe-bas respectifs 16 et 17. Un cir- cuit 19 fournit le quart de l'argument de chaque vecteur de composantes (xl, yl) fournies par les filtres 16 et 17. Ainsi, le circuit 19 fournit la phase moyenne p correspondant à 1 'angle de la moyenne des vecteurs de composantes (I, Q), dépourvue des erreurs de r/2 qui pouvaient être provoquées par des symboles ayant une erreur de phase hors des limites de -K/4 et 7/4. Un dérivateur 21 fournit l'erreur de fréquence F recherchée à partir de la phase moyenne p. Cette erreur est proportionnelle à la différence d entre deux angles p successifs, modulo n/2, c'est-à-dire que si la différence d est, par exemple, supérieure à 7/4, on prend la valeur d - c/2 qui est comprise entre -X/4 et x/4. L'erreur de fréquence F peut être utilisée directement pour corriger la valeur de consigne de fréquence d'une boucle à verrouillage de phase 23. On remarque que le circuit de la figure 2 ne calcule pas directement la moyenne des angles des vecteurs de composantes (x, y), mais calcule d'abord la moyenne de chacune des compo- santes x et y, c'est-à-dire un vecteur moyen, puis calcule l'angle de ce vecteur moyen. Cette façon de procéder constitue un aspect important de l'invention. En effet, elle permet de lever des ambiguités lorsque des symboles sont reçus près des axes I et Q de la constellation de la figure 1. Par exemple, on suppose que la constellation reçue est décalée de -x/4 par rapport à la constellation théorique, c'est- à-dire que les points correspondant aux points théoriques P1 à P4 se trouvent sur les axes I et Q. On dira alors que la phase de la constellation est nulle. Il est très probable de recevoir un symbole de coordonnées (1, E) correspondant au point Pl, et un symbole de coordonnées (e, 1), correspondant au point P2, o E est une faible composante parasite liée au bruit. Ces points sont à des angles respectifs E et 7/2 - E. La multiplication par quatre de ces angles fournit respectivement 4E et 2x 4E. La moyenne de ces angles, divisée par 4, fournit la valeur x/4 counre estimation de la phase, ce qui est incorrect puisque la phase est nulle dans cet exemple. Par contre, les composantes des vecteurs après transformation étant respectivement (1, 4E) et (1, -4E), le vecteur moyen correspondant, de composantes (1, 0), a un angle nul, ce qui est la valeur exacte de la phase de la constellation revue. On remarquera en outre que, si on n'applique pas la transformation selon l'invention, ni le calcul de l'angle moyen, ni le calcul du vecteur moyen, ne fournissent une phase correcte (X/4 dans les deux cas). Le fait de calculer d'abord le vecteur moyen présente également l'avantage de donner un effet prépondérant aux vecteurs de grand module qui ont une probabilité plus importante d'avoir un angle correct que les vecteurs de faible module. En effet, une même composante de bruit affecte beaucoup moins l'angle d'un vecteur de module élevé que l'angle d'un vecteur de faible module. Le circuit de la figure 2 est particulièrement complexe à réaliser en pratique du fait des convertisseurs cartésien- polaire et polaire-cartésien qui doivent réaliser des calculs trigonométriques. Il en va de même pour le circuit 19 qui fournit l'angle du vecteur de composantes (xl, yl). Selon un mode de réalisation de la présente invention, on propose un circuit fournissant un résultat final équivalent à celui de la figure 2, mais qui soit plus simple à réaliser. Pour cela, on utilise une transformation linéaire qui ne transforme plus un cercle en un cercle comme le ferait le circuit de la figure 2, mais qui multiplie par quatre les angles multiples de ic/4. Un exemple de telle transformation est la suivante: en notant Sup = maxdII, IQI), Inf = min(lIIl, IQI), et Sgn = signe(IQ(IIl-IQI)), si 4Inf les composantes des vecteurs image. La figure 3 représente l'image obtenue selon cette transformation d'un cercle de rayon 1 de la constellation de la figure 1. Les images des points théoriques Pl à P4 se trouvent sur l'axe des x, du côté négatif. Les vecteurs d'angle compris entre, sensiblement, 30 et 60 (modulo 90 ) ont des images dans les deuxième et troisième quadrants de la figure 3. Les vecteurs d'angle compris entre 0 et sensiblement 30 (modulo 90 ) ont des images dans le premier quadrant, tandis que les vecteurs d'angle compris entre sensiblement 60 et 90 (modulo 90 ) ont des images dans le quatrième quadrant. Comme dans le circuit de la figure 2, les filtres passe-bas 16 et 17 fournissent les moyennes xl et yl des compo- santes x et y obtenues après cette transformation. D'autres transformations linéaires par morceaux pour- ront être envisagées. Elles devront toutefois, comme la transformée ci-dessus, multiplier par quatre au moins les angles multiples de 7/4. Les circuits 19 et 21 de calcul d'angle et de dériva- tion de la figure 2 pourraient être conservés. Toutefois, le circuit de calcul d'angle 19 doit effectuer un calcul d'arc- tangente, ce qui est relativement complexe et incompatible avec la simplification de calcul obtenue grâce à la transformée linéaire par morceaux qui précède. Selon un mode de réalisation, on propose d'utiliser seulement la polarité de l'angle image moyen pour fournir une indication d'erreur de fréquence. Plus exactement, si le point image moyen (xl, yl) est situé dans le deuxième ou troisième quadrant de la figure 3, on affecte une valeur nulle à la phase p. Si le point image moyen est situé dans le premier quadrant, on affecte la valeur 1 à la phase p. Finalement, si le point image moyen est situé dans le quatrième quadrant, on affecte la valeur -1 à la phase p. La dérivation effectuée par le circuit 21 consiste à soustraire la valeur précédente de la phase p à la valeur courante, le résultat étant fourni modulo 2 si la dif- férence est positive et modulo -2 si la différence est négative (en d'autres termes, la différence est annulée si elle est égale à 2 ou à -2). Ainsi, lorsque la constellation reçue tourne à vitesse constante, cette dérivation fournit une succession de valeurs 1 ou -1 (selon le sens de rotation), séparées par un même nombre de valeurs nulles, le nombre de valeurs nulles augmentant lorsque la vitesse de rotation diminue. Cette dérivation ne fournit donc pas directement une valeur de fréquence que l'on peut appliquer en consigne à la PLL 23, mais des incréments positifs ou négatifs que l'on applique à cette consigne, les incrémennts étant d'autant plus rapprochés que la vitesse de rotation de la constellation reçue est élevée. Les opérations décrites ci-dessus seront généralement effectuées sur des signaux numériques, par exemple de manière câblée. Le fait que le dernier mode de réalisation décrit ne nécessite que des calculs linéaires permet d'utiliser un circuit câblé particulièrement simple ne devant pas réaliser d'opérations trigonométriques. La figure 4 représente un exemple numérique des filtres 16 et 17. La composante x ou y est fournie à un additionneur 30 qui reçoit également la sortie xl ou yl du filtre par l'intermé- diaire d'un élément 32 de retard d'un échantillon. La sortie xl ou yl du filtre est fournie par un soustracteur 34 qui calcule la différence entre la sortie de l'additionneur 30 et la sortie de l'élément de retard 32, multipliée par un coefficient a. De bons résultats sont obtenus avec a = 1/8. REVEtNDICATICNS 1. Procédé d'estimation de l'erreur de fréquence d'un démodulateur destiné à reconstruire deux signaux binaires (I, Q) véhiculés sur deux porteuses de même fréquence mais en quadrature de phase, caractérisé en ce qu'il comprend les étapes suivantes: - former des vecteurs ayant pour composantes les cou- ples successifs de valeurs (I, Q) des deux signaux binaires; - appliquer à chaque vecteur une transformation qui multiplie par quatre son angle au moins lorsqu'il est égal à un multiple de x/4 et conserve sensiblement son module; et - calculer la moyenne (xl, yl) des vecteurs transfor- mes. 2. Procédé d'estimation d'erreur de fréquence selon la revendication 1, caractérisé en ce que la transformation est une rotation de quatre fois l'angle (0) du vecteur. 3. Procédé d'estimation d'erreur de fréquence selon la revendication 1, caractérisé en ce que l'erreur de fréquence est fournie par le calcul de la dérivée de l'angle (p) du vecteur moyen. 4. Procédé d'estimation d'erreur de fréquence selon la revendication 3, caractérisé en ce que la dérivée est fournie sous forme de différence, modulo x/2, de deux valeurs d'angle successives. 5. Procédé d'estimation d'erreur de fréquence selon la revendication 1, caractérisé en ce que la transformation est une transformation linéaire par morceaux qui fait correspondre à chaque vecteur un vecteur image situé sensiblement dans le même quadrant du plan image qu'un vecteur image obtenu par multiplica- tion par quatre de l'angle du vecteur. 6. Procédé d'estimation d'erreur de fréquence selon la revendication 2 ou 5, caractérisé en ce que l'on attribue une phase de O aux vecteurs image situés dans les deuxième et troi- sième quadrants du plan image, une phase +1 aux vecteurs image situés dans le premier quadrant, et une phase -1 aux vecteurs image situés dans le quatrième quadrant, une information d'erreur de fréquence étant fournie par la différence des valeurs de phase attribuées au vecteur image courant et au vecteur image précé- dent, laquelle différence est fournie modula 2 si elle est positive et modulo -2 si elle est négative.