FR 2658016 A1 19910809 FR 9001491 A 19900206 i Procédé de diffusion de données numériques, notamment pour la radiodiffusion à haut débit vers des mobiles, à entrelacement tempsfréquence et démodulation cohérente, et récepteur correspondant. Le domaine de l'invention est celui de la diffusion de données numériques destinées à être reçues notamment par des mobiles circulant en milieu urbain, c'est-à-dire dans des conditions de multipropagation (processus de Rayleigh) engendrant des phénomènes d'évanouissement (fading), et en présence de parasites et de brouillage Plus généralement, l'invention concerne la diffusion de signaux numériques dans des canaux affectés de trajets multiples dont les caractéristiques varient dans le temps. L'invention s'applique plus particulièrement, mais non exclusivement, au système de radiodiffusion sonore numérique dénommé système COFDM (Coding Orthogonal Frequency Division Multiplex (Multiplexage de Fréquences Or- thogonales Codées)), tel que décrit dans les demandes de brevets françaises 86 09622 du 2 juillet 86, et 86 13271 du 23 septembre 86 au nom des mêmes déposants. Ce système de diffusion numérique est basé sur l'utilisation conjointe d'un dispositif de codage de canal et d'un procédé de modulation par multiplexage de fréquences orthogonales. Le procédé de modulation proprement dit de ce système connu permet de s'affranchir des problèmes liés à la sélectivité en fréquence du canal Il consiste à assurer la répartition d'éléments numériques constitutifs du signal de données dans l'espace fréquence f temps t et à émettre simultanément des jeux d'éléments numériques sur N voies de diffusion parallèles au moyen d'un multiplex de fréquences utilisant des porteuses orthogonales En particulier, ce type de modulation permet d'éviter que deux éléments successifs du train de données soient émis à la même fréquence. Le processus de codage connu vise dans son ensemble à permettre le traitement des échantillons issus du démodulateur pour absorber l'effet des variations d'amplitude du signal reçu, dues au processus de Rayleigh Le codage est avantageusement un codage convolutif, éventuellement concaténé à un codage du type Reed-Solomon Le décodage peut être avantageusement un décodage de Viterbi à décision douce. De façon connue, les éléments numériques codés sont de plus entrelacés, en temps et en fréquence, de façon à maximiser l'indépendance statistique des échantillons vis-à-vis du processus de Rayleigh et du caractère sélectif du canal. La démodulation du signal reçu peut être différentielle ou cohérente. On sait que la démodulation cohérente offre, en théorie, une meilleure résistance au bruit que la démodulation différentielle, permettant d'obtenir un gain en performance d'environ 3 d B Toutefois, il apparaît clairement que, dans les conditions de réception spécifiques du procédé, il est particulièrement difficile d'extraire du signal modulé une référence de phase et d'amplitude pour chaque porteuse du multiplex L'erreur faite sur l'estimation des porteuses entraîne donc, dans le cas de la démodulation cohérente, une dégradation sensible des performan- ces Ceci est d'autant plus vrai dans le cas d'évanouissement profonds et rapides rencontrés lorsque la fréquence porteuse ou la vitesse du véhicule augmente. La démodulation différentielle constitue une autre solution dont l'intérêt essentiel réside dans sa simplicité de mise en oeuvre et son absence d'inertie après un évanouissement profond C'est donc cette seconde approche qui a été utilisée, pour valider les principes généraux du système COFDM, au détriment d'une dégraduation des performances d'environ 3 d B par rapport à la démodulation différentielle. L'invention a pour objectif de pallier ces inconvénients, et d emontrer qu'il existe cependant une possibilité élégante de mettre en oeuvre une démodulation cohérente. Plus précisément, l'invention a pour but de fournir un procédé de diffusion numérique à entrelacement temps-fréquence permettant une démodulation cohérente, en facilitant l'extraction de référence de phase et d'amplitude. Un autre objectif de l'invention est de limiter les pertes de capacité de transmission d'un tel procédé. Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints par un procédé de diffusion de données numériques, notamment pour la radiodiffusion à haut débit vers des mobiles, du type assurant la répartition desdites données sous forme d'éléments numériques dans l'espace fréquence f temps t et l'émission de symboles constitués chacun d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu des éléments numériques, et diffusées simultanément, certains desdits éléments numériques étant des éléments de référence de valeur et de position dans ledit espace fréquence f temps t connues des récepteurs. Ainsi, il est possible de déterminer, par interpolation, une référence de phase et d'amplitude pour chacun des éléments numériques, et donc de réaliser une démodulation cohérente. Il est nécessaire, pour effectuer une interpolation efficace, de disposer d'un nombre suffisant d'éléments de référence. De façon avantageuse, ces éléments de référence sont répartis dans l'espace fréquentiel de telle façon que l'écart fréquentiel Z 3 F entre deux desdits éléments de référence est inférieur ou égal à 1/( 2 Tm), 'Tm étant la valeur maximale de l'étalement de la réponse impulsionnelle du canal de transmission, et, dans l'espace temporel, de façon que l'écart temporel LT entre deux desdits éléments de référence est inférieur ou égal à l/( 2 ff 2;),, étant la fréquence Doppler maximale déterminée selon l'équation fÉma Lu = fo v/c, o fo est la fréquence porteuse maximale, v est la vitesse de déplacement maximale desdits mobiles, c est la célérité de la lumière. Ces éléments peuvent être placés de façon quelconque dans l'espace tempsfréquence Cependant, on les disposera de préférence de façon qu'ils soient équirépartis dans l'espace fréquence f temps t, et par exemple, en quinconce dans ledit espace fréquence f temps t, par rapport à la dimension fréquentielle et/ou à la dimension temporelle. L'invention concerne également des récepteurs susceptibles d'exploiter ce procédé. Avantageusement, de tels récepteurs comprennent des moyens d'analyse desdits éléments de référence, de façon à déterminer par interpolation la réponse du canal de transmission en phase et/ou en amplitude pour chacun desdits éléments numériques. Selon un mode de réalisation avantageux, lesdits moyens d'analyse comprennent: des premiers moyens de mémorisation des données avant filtrage d'interpolation, des moyens de filtrage d'interpolation selon l'axe fréquentiel, des seconds moyens de mémorisation des résultats dudit filtrage selon l'axe fréquentiel, des moyens de filtrage d'interpolation selon l'axe temporel, pour chacune desdites fréquences porteuses. De façon préférentielle, lesdits moyens de filtrage d'interpolation selon l'axe fréquentiel comprennent un filtre à réponse impulsionnelle finie et lesdits moyens de filtrage d'interpolation selon l'axe temporel comprennent une batterie de N filtres à réponse impulsionnelle finie, fonctionnant indépendamment sur chacune desdites N fréquences porteuses. Ces récepteurs peuvent alors avantageusement comprendre des moyens de démodulation cohérente, par projection des porteuses modulées par lesdits éléments numériques, cette projection s'effectuant sur les résultats fournis par lesdits moyens d'analyse. 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 d'exemple illustratif et non limitatif, et des dessins annexés dans lesquels: la figure 1 est un schéma synoptique d'une chaîne d'émission-réception mettant en oeuvre le procédé de l'invention; la figure 2 est une représentation du multiplex de fréquences porteuses dans le domaine temps-fréquence, comprenant des éléments de référence et des éléments d'information; la figure 3 illustre le fonctionnement du filtrage d'interpolation selon l'invention sur l'axe fréquentiel la figure 4 illustre le fonctionnement du filtrage d'interpolation selon l'invention sur l'axe temporel; la figure 5 illustre un exemple de réalisation du module d'estimation des porteuses selon l'invention. les figures 6 A, 6 B et 6 C illustrent différents modes de mise en place des éléments de référence. Les différents aspects du mode de réalisation qui va être décrit ci-après concernent plus particulièrement la radiodiffusion sonore numérique vers des mobiles. Toutefois, il est clair que le principe de diffusion numérique à haut débit de l'invention peut s'appliquer à tout type de communication, notamment dans des canaux soumis au processus de Rayleigh. Dans l'application de radiodiffusion sonore numérique, un objectif peut être par exemple de transmettre seize programmes stéréophoniques dans une bande de fréquence de 8 M Hz de largeur, avec un débit numérique de l'ordre de 250 kbits après compression par programme stéréophonique. Une chaîne de transmission telle que décrite dans les demandes de brevet mentionnées en préambule est présentée en figure 1. Chacun des N ( 16) canaux C O à C n 1 subit en parallèle un codage 10, puis un entrelacement temps-fréquence 11 sur une voie séparée, avant d'être soumis collectivement à un processus 12 de multiplexage temporel et de modulation OFDM. Le codage 10 est avantageusement du type convolutif L'entrelacement temps-fréquence 11 a pour objectif de brasser les éléments numériques de chaque canal afin de leur conférer une indépendance maximale vis-à-vis des parasites et du brouillage du canal de diffusion 13. La modulation OFDM, consiste à moduler des symboles constitués chacun d'un multiplex de fréquences orthogonales diffusées simultanément sur N voies. Cette opération peut être réalisée par une transformation de Fourier sur la séquence numérique codée et entrelacée de chaque canal Ci. A titre d'exemple, dans une bande de fréquence de 8 M Hz, on peut définir 512 fréquences porteuses séparées de 15 625 Hz Parmi celles-ci 448 sont utilisables, après élimination de la fréquence centrale du spectre et des porteuses latérales ( 1/8 e du spectre) pour tenir compte des contraintes de filtrage. La chaîne de réception comporte les étapes classiques de filtrage et de transposition en bande de base 14, de démodulation 15, de désentrelacement en fréquence 16 et de décodage 17 du canal désentrelacé, avantageusement du type décodage de Viterbi à décision douce. Le signal émis est constitué par une suite de symboles de modulation formant un multiplex de N porteuses orthogonales. Soit {fk} l'ensemble des fréquences porteuses considérées, avec: fk = f O + k/ts,k = O à N-1 o ts représente la durée allouée à un symbole de modulation. On peut alors définir la base de signaux élémentaires suivante jk (t) avec k = O à N-1, j = à + 00: i Yj,k (t) = gk (t -j Ts) avec: O valeurs dans un alphabet fini, représentant le signal de données émis. Le signal OFDM associé s'écrit alors: +co N-1) x(t) = Re ( N Cjk t"j,k(t) j=.00 k=O Pour s'affranchir de tout problème de sélectivité en fréquence du canal, un intervalle de garde de durée LS est inséré devant chaque signal 'jk(t) afin d'absorber le brouillage intersymboles. ts représente donc la durée du signal utile, ila durée de l'intervalle de garde et Ts = ts + 4, la durée du symbole. Le canal de transmission peut être modélisé par la relation L J Y 3 k = H k C k + N k o Hk est la réponse du canal à la fréquence fk et à l'instant j Ts. Njk est un bruit gaussien complexe. Le dispositif de récupération de porteuse utilisé en démodulation cohérente doit pouvoir fournir une estimation de la réponse du canal à tout instant et pour toutes les fréquences porteuses, sous la forme jk P Jk ek o: 8 jk est l'amplitude de la réponse du canal, -jk est la phase de la réponse du canal. L'idée de base de l'invention consiste à utiliser certaines porteuses judicieusement réparties dans le domaine temps-fréquence comme fréquences pilotes de référence de phase et/ou d'amplitude. En d'autres termes, on insère parmi les données à transmettre et à des emplacements prédéterminés des mots fixes, servant de référence en amplitude et/ou en phase de réception. La figure 2 illustre un exemple de répartition de ces fréquences pilotes dans le domaine temps-fréquence Parmi les éléments d'information 21 sont introduits des éléments de référence 22, de position dans le domaine temps- fréquence et de valeur connues des récepteurs, permettant l'estimation de la réponse du canal. Celle-ci peut être ensuite obtenue à tout instant j Ts et pour toutes les fréquences fk par filtrage d'interpolation, si l'emplacement entre les fréquences pilotes prend en compte les caractéristiques du milieu de transmission à trajets multiples. la première caractéristique de ce milieu est la variation en fréquence de la fonction de transfert du canal On sait que la bande de cohérence fréquentielle (C f)c est directement déterminée par l'inverse de l'étalement de la réponse impulsionnelle du canal Tm: (Z-f)C z /Tm. On pourra donc obtenir une bonne estimation de la réponse du canal à toutes les fréquences fk si l'écart a F entre deux fréquences pilotes vérifie la relation suivante: LF c: célérité de la lumière. Ainsi, on aura une bonne estimation de la réponse du canal à tous les instants j Ts si les symboles de référence apparaissent avec une périodicité T telle que: AT 1 I( 2 fnzar) = (L\ t)c/2. La fonction d'interpolation permettant d'obtenir ces estimations peut notamment être réalisée par des techniques de filtrage numérique à réponse impulsionnelle finie Toutefois, toutes les techniques de filtrage connues peuvent être envisagées, et par exemple le filtrage récursif. Dans un premier temps, ce filtrage d'interpolation est réalisée selon l'axe fréquentiel, sur chacun des symboles portant des échantillons de référence, c'est- à-dire, dans l'exemple de la figure 2, sur les lignes 0, 3, 6 La figure 3 représente le principe de reconstruction par interpolation de la réponse du canal A partir des échantillons de référence reçus 31 A 4 31 B' 31 C, correspondant à l'émission à un mot de référence identique, et espacés de b F, un filtre de fonction de transfert H(f) à réponse impulsionnelle finie permet de reconstituer les échantillons manquants 32 A, 32 B, 32 c, distants de fsj, par interpolation. Pour cela, le filtre de fonction de transfert H(f) doit avoir une fréquence d'échantillonnage fs 1 égale à ( 1/ts) Le nombre de coefficients Q de ce filtre doit d'autre part vérifier la relation: Q > > (àF/fsl). Dans un second temps, un deuxième filtrage d'interpolation utilisant les résultats précédents est effectué, sur l'axe temporel, pour chaque fréquence. Ce filtrage est réalisé par une batterie de N filtres en parallèle, chacun travaillant indépendamment sur une fréquence fk. La figure 4 illustre cette interpolation dans le domaine temporel Les échantillons de référence reçus 41 A' 41 B, 41 c permettent, par l'intermédiaire du filtre de fonction de transfert G(F) à réponse impulsionnelle finie, la détermination des échantillons manquants 42 A, 42 B, 42 C espacés de Ts. La fréquence d'échantillonnage de chacun des filtres de fonction de transfert G(F) doit être: fs 2 = ( 1/Ts) et leur nombre P de coefficients doit vérifier P >> (AT/Ts). La figure 5 présente un mode de réalisation du module d'estimation des porteuses selon l'invention. Les données 51 issues du module de désentrelacement en fréquence 52, qui réalise avantageusement une transformation de Fourier, sont dirigées d'une part vers un module 53 d'estimation des porteuses, et d'autre part vers un module 54 de projection sur les références. Ce module 54 réalise, à l'aide de l'estimation des porteuses 55, une démodulation cohérente des données 51 Après démodulation, le récepteur réalise les opérations de désentrelacement 56 et de décodage 57. Dans un mode de réalisation particulier, l'estimation des porteuses peut être réalisée de la façon suivante: les données 51 sont stockées dans une mémoire 58 Lorsque le nombre voulu de données est disponible en mémoire, un premier filtrage 59 est réalisé, selon l'axe fréquentiel Les données issues de ce filtrage sont stockées dans une seconde mémoire 60, puis filtrées salon l'axe temporel Ce filtrage temporel est effectué par N filtres indépendants 61 A, 61 B, 61 C agissant chacun sur une fréquence porteuse. Ces différents modules 58, 59, 60, 61 A, 61 B, 61 C sont pilotés par un séquenceur 62, qui assure la bonne synchronisation des opérations. Cet exemple de réalisation n'est pas limitatif de l'invention De nombreu- ses autres possibilités apparaîtront à l'homme du métier Notamment, il est possible de réaliser d'abord le filtrage temporel, puis ensuite le filtrage fréquentiel. Avantageusement, les éléments de référence sont placés régulièrement Les figures 6 A, 6 B et 6 C illustrent différentes possibilités d'implantation La figure 6 A présente le cas d'une implantation en quinconce, par rapport aux fréquences, la figure 6 B une implantation en quinconce par rapport au temps, et la figure 6 C en quinconce par rapport aux dimensions. A titre d'exemple d'application, on présente ci-dessous les caractéristiques du procédé de l'invention dans le cas de la modulation COFDM telle qu'elle est actuellement définie. Le nombre de porteuses du multiplex est N = 448 Chaque symbole est de longueur Ts = 80 gs, et a donc une durée utile de ts = 64 gs L'écart entre les porteuses fs est donc de 15625 Hz. Si on considère un étalement de la réponse impulsionnelle du canal Tm = 8 gs, ce qui correspond au cas extrême d'un écart de 2,5 km entre un signal et son écho simple, la bande de cohérence fréquentielle (ti F)c est de 125 k Hz. L'écart entre deux éléments de référence doit donc être tel quetn Fp Une référence toutes les quatre porteuses permet donc d'obtenir une bonne interpolation sur l'axe fréquentiel. Avec une fréquence centrale f O de 1 G Hz et une vitesse du véhicule de km/h, la fréquence Doppler maximale fmax est d'environ 100 Hz La bande de cohérence temporelle du canal (f)c étant de 10 ms, on en déduit TP En choisissant par exemple,8 TP = 0 8 ms, on a besoin d'une référence tous les 10 symboles: b Tp 8 ( 00 = 10 Ts 80 Dans ce cas, les éléments de référence n'occupent donc qu'un quarantième de la ressource, ce qui représente une perte relativement faible. il REVENDICATIONS 1 Procédé de diffusion de données numériques, notamment pour la radiodiffusion à haut débit vers des mobiles, du type assurant la répartition desdites données sous forme d'éléments numériques dans l'espace fréquence f temps t et l'émission de symboles constitués chacun d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu des éléments numériques, et diffusées simultanément, certains desdits éléments numériques étant des éléments de référence ( 22), de valeur et de position dans ledit espace fréquence f temps t connues des récepteurs. 2 Procédé selon la revendication 1, caractérisé en ce que l'écart fréquentiel 1 F entre deux desdits éléments de référence ( 22) est inférieur ou égal à 1/( 2 Tm), Tm étant la valeur maximale de l'étalement de la réponse impulsion- nelle du canal de transmission. 3 Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que l'écart temporel, T entre deux desdits éléments de référence ( 22) est inférieur ou égal à 1/( 2 f 112 a), f,7 a étant la fréquence Doppler maximale déterminée selon l'équation furax = f O v/c, o: f O est la fréquence porteuse maximale, v est la vitesse de déplacement maximale desdits mobiles, c est la célérité de la lumière. 4 Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que lesdits éléments de référence ( 22) sont équirépartis dans l'espace fréquence f temps t, et par exemple, en quinconce dans ledit espace fréquence f- temps t, par rapport à la dimension fréquentielle et/ou à la dimension temporelle. 5 Récepteur d'émissions diffusées selon le procédé de l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il comprend des moyens ( 53) d'analyse desdits éléments de référence ( 22), de façon à déterminer par interpola- tion la réponse du canal de transmission en phase et/ou en amplitude pour chacun desdits éléments numériques. 6 Récepteur selon la revendication 5, caractérisé en ce que lesdits moyens ( 53) d'analyse comprennent: 4 ' 2658016 des premiers moyens ( 58) de mémorisation des données avant filtrage d'une interpolation, des moyens ( 59) de filtrage d'interpolation selon l'axe fréquentiel, des seconds moyens ( 60) de mémorisation des résultats dudit filtrage selon l'axe fréquentiel, des moyens ( 61 A, 61 B, 61 c) de filtrage d'interpolation selon l'axe temporel, pour chacune desdites fréquences porteuses. 7 Récepteur selon la revendication 6, caractérisé en ce que lesdits moyens ( 59) de filtrage d'interpolation selon l'axe fréquentiel comprennent un filtre à réponse impulsionnelle finie. 8 Récepteur selon l'une quelconque des revendications 6 ou 7, caractérisé en ce que lesdits moyens ( 61 A, 61 B, 61 C) de filtrage d'interpolation selon l'axe temporel comprennent une batterie de N filtres à réponse impulsionnelle finie, fonctionnant indépendamment sur chacune desdites N fréquences porteuses. 9 Récepteur selon l'une quelconque des revendications 5 à 8, caractérisé en ce qu'il comprend des moyens de démodulation cohérente ( 54), par projection des porteuses modulées par lesdits éléments numériques, cette projection s'effectuant sur les résultats fournis par lesdits moyens ( 53) d'analyse.