FR 2768278 A1 19990312 FR 9711526 A 19970911 Procédé d'estimation d'un décalage de phase parasite lors de la réception d'un signal multiporteuse, et récepteur correspondant. Le domaine de l'invention est celui de la réception de signaux multiporteuses. Plus précisément, l'invention concerne l'estimation du décalage de phase introduit par les opérations de démodulation. L'invention s'applique à tous les types de signaux mettant en oeuvre une pluralité de fréquences porteuses, c'est-à-dire aux systèmes mettant en oeuvre des signaux transmis selon la technique de multiplexage par répartition en fréquence (en anglais: Frequency Division Multiplex (FDM)), et par exemple le système COFDM (Coded Orthogonal Frequency Division Multiplex (multiplexage de fréquences orthogonales codées)), mis en oeuvre notamment dans le cadre des projets européens Eureka 147 "DAB" (Digital Audio Broadcasting (diffusion audionumérique)), et RACE dTTb (digital Terrestrial TV broadcasting (diffusion numérique de signaux de télévision terrestre)). Dans de tels systèmes de transmission, les données source à transmettre sont organisées en symboles (constitués d'une ou plusieurs données source) modulant chacun, pendant un intervalle de temps prédéterminé, une fréquence porteuse choisie parmi une pluralité de porteuses. Le signal formé par l'ensemble des porteuses modulées est transmis vers un ou plusieurs récepteurs, qui reçoivent un signal émis pertubé par le canal de transmission. La démodulation consiste généralement, dans son principe, à estimer la réponse du canal de transmission pour chaque porteuse et chaque symbole, puis à diviser le signal reçu par cette estimation pour obtenir une estimation du symbole émis. On connaît de nombreuses techniques de démodulation, qui peut être différentielle ou cohérente. Une technique facilitant la démodulation cohérente à l'aide de symboles de référence (ou "pilotes") connus des récepteurs et régulièrement insérés parmi les symboles utiles est par exemple décrite dans le brevet FR-94 07984 au nom des mêmes déposants. Le signal numérique multiporteuse considéré inclut donc un certain nombre de "pilotes", c'est-à-dire des porteuses modulées par de valeurs connues du récepteur. Soit C l'ensemble des indices de ces portueuses. Pour k appartenant à C, on note Pnk la valeur complexe modulant la porteuse numéro k pendant la durée symbole numéro n. On rappelle que cette valeur est connue du récepteur. En réception, on note Rn,k la valeur observée sur le nième symbole de la porteuse numéro k. Typiquement: Rn,k = Pn,k Hn,k eJ% (+ bruit) (1) o Hnk est la réponse fréquentielle complexe du canal à la fréquence de la porteuse numéro k; On est une phase introduite par la démodulation et que l'on cherche à estimer (On est commune à toutes les porteuses). Hn,k varie lentement dans le temps (Hnk peu différent de Hn+l,k), alors que On peut varier considérablement d'un symbole au suivant si l'oscillateur local du récepteur est mal réglé. Estimer On sert donc par exemple à connaître l'erreur de fréquence de l'oscillateur local afin de la corriger ensuite (Contrôle Automatique de Fréquence, ou CAF). En outre, comme décrit dans les documents de brevet français FR-95 10067 et FR-95 10068, au nom des mêmes déposants que la présente demande de brevet, la connaissance de Qn permet de diminuer le biais introduit par une distorsion blanche en fréquence, et de tenir compte d'une estimation d'une telle distorsion blanche pour la démodulation. En pratique, seules les variations de An sont à estimer: la valeur de 4o, fixe par définition, peut être "incorporée" dans la réponse du canal Hnlk. Comme Hnk est peu différent de Hn+lk' (tn-0n.1) noté n, est classiquement estimé à partir des Rnk suivant la formule: 6. Arg (R Pk)(R P_ k} (2) Cette approche connue s'explique de la façon suivante. Notons Ank le module carré de Hn,k et Qn lesquels An3k varie beaucoup avec k. Finalement, la phase 8n de la somme donne une estimation de (On-4nl) dont la fiabilité augmente avec le nombre de porteuses de C. En pratique, les Pnk ont généralement le même module Q. Pnk peut même être complètement indépendant de n (porteuses continues): dans ce cas, la somme (2) peut être simplifiée car P*nkPn1,k est toujours égale à Q. de la façon suivante: on Arg { E [Rn, R. (2bis) Dans la suite, on gardera la forme de l'équation (2) pour ne pas perdre en généralité, mais la forme (2bis) est évidemment à préférer dans le cas de porteuses continues. Dans certains cas, on cherche à estimer non pas (nn-l) mais (On-(n.-t), t>l1; la méthode reste bien sûr la même. Il suffit de remplacer l'indice n- 1 par n-t. Les inventeurs ont observés que, si cette technique classique est efficace en condition de transmission peu brouillée, elle s'avère déficiente dans certaines conditions, en présence de brouilleurs. En effet, dans ce cas, certaines porteuses sont polluées en permanence par un bruit très fort, parfois nettement plus fort que le signal utile. Dans la somme (2), cela se traduit par quelques termes complètement faux, généralement dotés d'une forte amplitude. Comme n est assez petit (surtout lorsque l'oscillateur local tend vers la bonne fréquence), ces quelques termes sont susceptibles de perturber complètement le résultat de la somme. L'invention a notamment pour objectif de pallier cet inconvénient de l'étude de l'art. Plus précisément, un objectif de l'invention est de fournir une technique permettant d'atténuer fortement l'effet des brouilleurs sur l'estimation du décalage de phase induit par le récepteur. En d'autres termes, l'invention a pour objectif de fournir un procédé, et un récepteur correspondant, qui optimise l'estimation du décalage de phase, notamment lorsque ce dernier est faible. Un autre objectif de l'invention est de fournir une telle technique, qui soit simple à mettre en oeuvre, et qui ne nécessite pas des calculs complexes, ni de moyens spécifiques dans les récepteurs. L'invention a encore pour objectif de founir une telle technique, qui ne nécessite pas de modification du signal à emettre, et qui soit donc compatible avec les techniques d'émission déjà mises en oeuvre. 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'un décalage de phase parasite lors de la réception d'un signal multiporteuse formé d'une succession temporelle de symboles modulant une pluralité de fréquences porteuses, au moins certaines desdites fréquences porteuses d'au moins certains desdits symboles portant des éléments de référence, dont la valeur à l'émission est connue du récepteur effectuant ladite réception, procédé du type comprenant une étape d'estimation de la variation de phase an entre au moins deux symboles portant des éléments de référence, selon lequel la contribution de chacun desdits éléments de référence dans ladite estimation est pondérée par une information représentative du bruit affectant la fréquence porteuse portant ledit élément de référence. Ainsi l'importance attribuée à chaque pilote est fonction des perturbations qui l'affectent, et donc de sa crédibilité. Il est à noter que cette approche nouvelle n'est pas évidente, mais repose notamment sur la formulation du problème exposé plus haut, qui n'avait jamais été faite ou envisagée. On rappelle que d'une façon générale, le décalage de phase est induit notamment, mais non exclusivement, par le récepteur. Ladite information représentative du bruit peut notamment appartenir au groupe comprenant: - la variance du bruit affectant chaque fréquence porteuse (a2k); - l'amplitude du produit Rn,k.R*n-i,k; - le taux d'erreur affectant chaque fréquence porteuse; - une information indiquant que la porteuse n'est pas fiable. Ces informations, et notamment la variance du bruit, sont souvent déjà calculées pour d'autres applications, et notamment le calcul du rapport signal à bruit (S/B) sur chaque porteuse pour optimiser le décodage de Viterbi à entrées pondérées ("soft decoding"). Les calculs à effectuer sont donc très réduits. Ul s'agit simplement de pondérer les calculs habituels en fonction d'une information disponible. Ainsi, dans le cas d'une estimation de de la variance, l'estimation de la variation de phase on est avantageusement calculée de la façon suivante: Arg J [ 11 (3) = Arg [72 (R.,k.P,k)(R _lk.ePnIk) (3) kC ( k o: Pn,k est la valeur de l'élément de référence tel qu'émis; Rn,k est la valeur de l'élément de référence tel que reçu; n représente la dimension temporelle; k représente la dimension fréquentielle; C étant l'ensemble des fréquences porteuses portant des éléments de référence. Lorsque le récepteur calcule une estimation de Pnk = Ank/ca2k, F2k étant une estimation de la variance du bruit et Ank étant une estimation du module carré de la réponse fréquentielle du canal sur la fréquence porteuse k à l'instant n, l'estimation de la variation de phase on est par exemple calculée de la façon suivante: 1 0 n= Arg [pnk(Rn. k.P*.k)(R,k Pk)*] (4) IA kEc k -ck -1,k Selon un mode de réalisation avantageux de l'invention, ladite pondération consiste simplement en un seuillage, seuls les éléments de référence portés par des fréquences porteuses pour lesquelles ladite information représentative du bruit dépasse un seuil prédéterminé étant prises en compte. On vérifie en effet que cette solution, plus expéditive mais plus simple, est suffisante dans la plupart des cas. Dans ce cas, si ladite information représentative du bruit est une estimation de la variance du bruit, l'estimation de la variation de phase An peut être calculée de la façon suivante: =n Arg { [I(ak)(Rnk.P; ,k)(Rn _,k Plk)} (5) o: I(ak) vaut 1 tant que ak est en dessous dudit seuil, et 0 sinon. Si ladite information représentative du bruit est une estimation de Pn,k = Ank/a2k, l'estimation de la variation de phase on est avantageusement calculée de la façon suivante: n=Arg {kEc [J(pn,k) (Rk. P:n,k) (Rnl,k.PnIk)*]} (6) o: J(Pn.k) vaut 1 tant que Pk est au dessus dudit seuil, et 0 sinon. L'invention concerne également les récepteurs mettant en oeuvre un tel procédé. Ces récepteurs comprennent des moyens d'estimation pondérant la contribution de chacun desdits éléments de référence dans ladite estimation par une information représentative du bruit affectant la fréquence porteuse portant ledit élément de référence. Préférentiellement, ladite information représentative du bruit est une estimation de la variance du bruit affectant chaque fréquence porteuse (G2k). D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante d'un mode de réalisation préférentiel de l'invention, donné à titre de simple exemple illustratif et non limitatif de la figure unique, illustant de façon schématique le procédé de l'invention. Comme indiqué en préambule, l'invention s'applique à la réception des signaux multiporteuses, et concerne plus précisément l'estimation d'un décalage de phase induit principalement (mais non exclusivement) par le récepteur. Le procédé de l'invention peut être par exemple être mis en oeuvre pour optimiser la précision du décalage de phase pris en compte dans les documents de brevet français FR-95 10067 et FR-95 10068, au nom des mêmes déposants que la présente demande de brevet. L'invention permet également, très simplement, d'améliorer l'estimation du décalage de phase entre deux symboles consécutifs ou non, par exemple pour agir sur la commande automatique de fréquence (CAF). La figure unique illustre le principe de l'invention. Le signal transmis au travers du canal de transmission est classiquement reçu par des moyens de réception 1 qui effectuent les opérations classiques permettant de retrouver les différentes porteuses émises (amplification, filtrage, FFT...). Toujours classiquement, ont extrait (2) du signal reçu les pilotes reçus Rn,k. Ces valeurs Rn,k alimentent le module 3 de calcul de l'estimation du décalage de phase, qui reçoit également, d'une mémoire interne 4, les valeurs connues Pr,k des pilotes. Par ailleurs, selon l'invention, le module d'estimation 3 reçoit le calcul de la variance (2k, de moyens 5 de calcul de la variance du bruit, destinée également à d'autres applications classiques. Plusieurs techniques classiques sont connues pour calculer Dans une variante de l'invention, la variance (2k est remplacée par la valeur pn k=An.k/{52k. Il est également possible de tenir compte d'autres informations représentatives des brouilleurs, soit en remplacement, soit en complément. Le module d'estimation 3 délivre l'estimation de la variation de phase ân, qui peut être utilisée pour diverses opérations, et notamment pour contrôler le CAF 6. Cette information on tient compte, selon l'invention, d'un pondération fonction de la confiance que l'on peut attribuer à chaque porteuse, du fait des brouilleurs détectés. Dans le cas o le module 5 délivre l'information a2k, le module d'estimation peut donc calcusler: n =Arg { 2 (Rk P.Pk)(Rnlk. n1,k (3) k Selon une version simplifiée, on peut se contenter d'une pondération plus expéditive, mais nettement plus simple: n= Arg { * [I(k)(Rflk. P:k)(Rfl1k.P..Ik) * (5) O I(ck) vaut 1 tant que ak est en dessous d'un certain seuil S, et 0 lorsque sk dépasse S. Dans une mise en oeuvre matérielle, le bloc chargé d'estimer ak délivre un signal de contrôle lorsque ak dépasse S. Le module 3 chargé de calculer la somme (4) inhibe la prise en compte du terme courant sur réception de ce signal. Si le récepteur calcule Pn.k=An,k/(2k au lieu de (a2k, le calcul devient: n= Arg {k [p n,k(RnP.k)(Rn.kP lk)]} (4) Dans la version simplifiée, on a alors: n= Arg {kc [J(P,'k)(Rfk'P'k)(R!-_.k PflIk) (6) O J(Pnk) vaut 1 tant que Pnk est au dessus d'un certain seuil S, et 0 lorsque pn.k est inférieur à S. Dans une mise en oeuvre matérielle, le bloc chargé d'estimer Pn,k délivre un signal de contrôle lorsque Pn,k devient inférieur à S. Le module 3 chargé de calculer la somme (6) inhibe la prise en compte du terme courant sur réception de ce signal. REVENDICATIONS 1. Procédé d'estimation d'un décalage de phase parasite lors de la réception d'un signal multiporteuse formé d'une succession temporelle de symboles modulant une pluralité de fréquences porteuses, au moins certaines desdites fréquences porteuses d'au moins certains desdits symboles portant des éléments de référence, dont la valeur à l'émission est connue du récepteur effectuant ladite réception, procédé du type comprenant une étape d'estimation de la variation de phase 8n entre au moins deux symboles portant des éléments de référence, par analyse desdits éléments de référence, caractérisé en ce que la contribution de chacun desdits éléments de référence dans ladite estimation est pondérée par une information représentative du bruit affectant la fréquence porteuse portant ledit élément de référence. 2. Procédé selon la revendication 1, caractérisé en ce que ladite information représentative du bruit appartient au groupe comprenant: - la variance du bruit affectant chaque fréquence porteuse ((2k); - l'amplitude du produit Rnk. Rn-i,k; - le taux d'erreur affectant chaque fréquence porteuse; - une information externe indiquant que la porteuse n'est pas fiable. 3. Procédé selon la revendication 2, caractérisé en ce que ladite information représentative du bruit est une estimation de la variance du bruit, et en ce que l'estimation de la variation de phase on est calculée de la façon suivante: ân Arg L 2 (Rnk.pnk) (Rflk' nlk)I| o: Pnk est la valeur de l'élément de référence tel qu'émis; Rn,k est la valeur de l'élément de référence tel que reçu; n représente la dimension temporelle; k représente la dimension fréquentielle; C étant l'ensemble des fréquences porteuses portant des éléments de référence. 4. Procédé selon la revendication 2, caractérisé en ce que ladite information représentative du bruit est une estimation de Pn,k = Ank/(r2k, a2k étant une estimation de la variance du bruit et Ank étant une estimation du module carré de la réponse fréquentielle du canal sur la fréquence porteuse k à l'instant n, et en ce que l'estimation de la variation de phase in est calculée de la façon suivante: =Arg EPn,k( nkPn.k)( n-,k n-1,k) o: Pn,k est la valeur de l'élément de référence tel qu'émis; Rn,k est la valeur de l'élément de référence tel que reçu; n représente la dimension temporelle; k représente la dimension fréquentielle; C étant l'ensemble des fréquences porteuses portant des éléments de référence. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ladite pondération consiste en un seuillage, seuls les éléments de référence portés par des fréquences porteuses pour lesquelles ladite information représentative du bruit dépasse un seuil prédéterminé étant prises en compte. 6. Procédé selon la revendication 5, caractérisé en ce queladite information représentative du bruit est une estimation de la variance du bruit, et en ce que l'estimation de la variation de phase 5n est calculée de la façon suivante: o= Arg {kc [I(ak) (R k P:,k) (RnlkPI,k) *]} o: I(sk) vaut 1 tant que ak est en dessous dudit seuil, et 0 sinon; Pnk est la valeur de l'élément de référence tel qu'émis; Rnk est la valeur de l'élément de référence tel que reçu; n représente la dimension temporelle; k représente la dimension fréquentielle; C étant l'ensemble des fréquences porteuses portant des éléments de référence. 7. Procédé selon la revendication 5, caractérisé en ce que ladite information représentative du bruit est une estimation de Pn,k = Ank/&2k, G2k étant une estimation de la variance du bruit et Ank étant une estimation de l'amplitude carrée de la réponse fréquentielle du canal sur la fréquence porteuse k à l'instant n, et en ce que l'estimation de la variation de phase 5n est calculée de la façon suivante: an= Arg {k[J(Pnlk)(Rfl.k.P:,k)(Rn.,I kPnf-lIk)*]} o: J(Pnk) vaut 1 tant que Pk est au dessus dudit seuil, et O sinon; Pnk est la valeur de l'élément de référence tel qu'émis; Rnk est la valeur de l'élément de référence tel que reçu; n représente la dimension temporelle; k représente la dimension fréquentielle; C étant l'ensemble des fréquences porteuses portant des éléments de référence. 8. Récepteur d'un signal multiporteuse formé d'une succession temporelle de symboles modulant une pluralité de fréquences porteuses, au moins certaines desdites fréquences porteuses d'au moins certains desdits symboles portant des éléments de référence, dont la valeur à l'émission est connue du récepteur, du type comprenant des moyens d'estimation de la variation de phase ân entre au moins deux symboles portant des éléments de référence, par analyse desdits éléments de référence, caractérisé en ce que lesdits moyens d'estimation pondèrent la contribution de chacun desdits éléments de référence dans ladite estimation par une information représentative du bruit affectant la fréquence porteuse portant ledit élément de référence. 9. Récepteur selon la revendication 8, caractérisé en ce que ladite information représentative du bruit est une estimation de la variance du bruit affectant chaque fréquence porteuse (G2k).