Brevet FR2671923A1 — publié le 19920724

i Dispositif de démodulation cohérente de données numériques entrelacées en temps et en fréquence, à estimation de la réponse fréquentielle du canal de transmission et seuillage, et émetteur 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 brevets français 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-temps 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. L'intérêt de la démodulation différentielle réside essentiellement dans sa simplicité de mise en oeuvre et son absence d'inertie après un évanouissement profond C'est cette approche qui a été utilisée pour valider les principes généraux du système COFDM. 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 des conditions de réception spécifiques aux systèmes de diffusion vers les mobiles en milieu perturbé, 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 performances Ceci est d'autant plus vrai dans le cas d'évanouissements profonds et rapides rencontrés lorsque la fréquence porteuse ou la vitesse du véhicule augmente. En d'autres termes, la démodulation cohérente est en principe plus performante que la démodulation différentielle, mais nécessite que le dispositif de récupération de porteuse puisse fournir une bonne estimation de la réponse fréquentielle du canal à tout instant. On connaît, par la demande de brevet FR 90 01491, du 06 02 90, au nom des mêmes déposants, un procédé de diffusion à entrelacement temps- fréquence permettant une démodulation cohérente, en prévoyant l'insertion parmi les éléments d'information utiles à transmettre d'éléments de référence de valeur et de position dans l'espace fréquence f temps t connues des récepteurs. L'idée de base de ce procédé consiste donc à 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 à la réception. 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. Plus précisément, l'estimation de la réponse du canal est obtenue par un filtrage d'interpolation, soit par convolution circulaire, ainsi que cela est décrit dans la demande de brevet FR 90 01491 déjà mentionnée, soit par transformation de Fourier Cette seconde solution présente l'avantage de nécessiter, à qualité égale, un nombre d'opérations inférieur à la première. Toutefois, on constate que ces méthodes n'apportent pas en pratique des résultats satisfaisants En effet, si le gain en faveur de la démodulation cohérente par rapport à la démodulation différentielle est théoriquement de 3 d B lorsque la réponse du canal est estimée de manière parfaite, il n'est en pratique que de l'ordre de 0,5 d B. Ce mauvais résultat est dû essentiellement au fait que l'estimation de la réponse du canal est fortement entachée de bruit, ce qui nuit à la qualité de l'interpolation. L'invention a notamment pour objectif de pallier cet inconvénient de l'état de la technique. Plus précisément, l'invention a pour objectif de fournir un dispositif de démodulation cohérente d'un signal numérique entrelacé en temps et en fréquence, présentant un gain important par rapport à la démodulation différentielle, de l'ordre de 2 d B. L'invention a donc pour objectif de fournir un tel dispositif, dans lequel l'influence du bruit est réduite, et donc les résultats de l'interpolation améliorés. Un objectif particulier de l'invention est de fournir un tel dispositif, ne nécessitant qu'une adaptation simple et peu onéreuse des décodeurs connus, et utilisant des moyens de calcul et des informations déjà présents dans ces décodeurs. Un autre objectif de l'invention est de fournir un tel dispositif, ne nécessitant pas d'adaptation du signal à émettre. De façon complémentaire, l'invention propose également un procédé de diffusion, permettant, sans l'ajout d'éléments de référence supplémentaires, une interpolation de meilleure qualité à la réception. Le dispositif et le procédé de l'invention peuvent avantageusement être mis en oeuvre conjointement Toutefois, ils sont indépendants, et peuvent être utilisés l'un sans l'autre. Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints selon l'invention à l'aide d'un dispositif de démodulation cohérente d'un signal numérique du type constitué d'éléments numériques répartis dans l'espace temps- fréquence et émis sous la forme de symboles constitués d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu desdits éléments numériques et diffusées simultanément, ledit signal numérique comprenant des éléments de référence, de valeur et de position dans ledit espace temps-fréquence connues dudit dispositif de démodulation, le dispositif comprenant des moyens d'estimation par transformation de Fourier de la réponse fréquentielle du canal de transmission à tout instant, réalisant la transformation des échantillons reçus correspondant à des éléments de référence du domaine fréquentiel dans le domaine temporel, la multiplication dans le domaine temporel desdits échantillons transformés par une fenêtre temporelle rectangulaire (fa) et la transformation inverse après ladite multiplication des échantillons obtenus du domaine temporel dans le domaine fréquentiel, et lesdits moyens d'estimation comprenant des moyens de seuillage desdits échantillons dans ledit domaine temporel, assurant l'élimination systématique des échantillons inférieurs à un certain seuil. De cette façon, les éléments de référence reçus avec une faible puissance, donc les éléments les plus perturbés par le bruit, ne sont pas pris en compte. Avantageusement, lesdits moyens de seuillage tiennent compte pour le calcul dudit seuil de la valeur a 2 de la puissance du bruit affectant ledit canal de transmission. Dans un mode de réalisation préférentiel, ledit seuil a une valeur comprise entre 5 a et 6 a. De façon avantageuse lesdits moyens de seuillage tiennent également compte pour le calcul dudit seuil d'une estimation de la réponse impulsionnelle dudit canal de transmission. Dans un mode de réalisation simplifié de l'invention, le seuil peut être fixe. Préférentiellement, lesdits moyens de seuillage sont placés en amont desdits moyens de fenêtrage Toutefois, ces moyens peuvent également être placés en aval. De façon avantageuse, ladite transformation des échantillons reçus est une transformation de format égal au nombre M d'éléments de référence par symboles émis, et en ce que ladite multiplication par une fenêtre temporelle réalise simplement l'ajout d'une série de (N M) zéros aux M éléments de référence transformés. L'invention concerne également un procédé de démodulation cohérente d'un signal numérique du type constitué d'éléments numériques répartis dans l'espace temps- fréquence, comprenant une étape d'estimation par transformation de Fourier de la réponse fréquentielle du canal de transmission à tout instant, comprenant les étapes de: la transformation des échantillons reçus correspondant à des éléments de référence du domaine fréquentiel dans le domaine temporel, la multiplication dans le domaine temporel desdits échantillons transformés par une fenêtre temporelle rectangulaire fn, la transformation inverse après ladite multiplication des échantillons obtenus du domaine temporel dans le domaine fréquentiel, et comprenant de plus une étape de seuillage desdits échantillons dans ledit domaine temporel, assurant l'élimination systématique des échantillons inférieurs à un certain seuil. De façon complémentaire, l'invention présente un procédé de diffusion d'un signal numérique du type constitué d'éléments numériques répartis dans l'espace temps-fréquence et émis sous la forme de symboles constitués d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu desdits éléments numériques et diffusées simultanément, ledit signal numérique com- prenant des éléments de référence, de valeur et de position dans ledit espace temps-fréquence connues dudit dispositif de démodulation, comprenant des moyens d'émission desdites porteuses attribuant sélectivement aux porteuses portant des éléments de référence une puissance supérieure à celle utilisée pour les porteuses portant des éléments d'information utile. En effet, l'objectif de l'invention est de faciliter la distinction des éléments de référence par rapport au bruit induit par le canal de transmission Cela peut donc être réalisé à la réception, par seuillage, mais également à l'émission, en augmentant la puissance des éléments de référence Il sera alors bien sûr plus facile de les distinguer à la réception. Ces deux moyens sont clairement indépendants Toutefois, ils peuvent être avantageusement mis en oeuvre simultanément. 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 présente un schéma synoptique des moyens de seuillage d'un dispositif de démodulation cohérente selon l'invention, mettant en oeuvre une transformation de Fourier et un seuillage des échantillons; les figures 2 A et 2 B présentent respectivement deux exemples d'estima- tions théoriques de la réponse impulsionnelle du canal dans l'hypothèse o toutes les porteuses sont des porteuses de référence, et dans le cas de l'insertion d'une porteuse de référence toutes les R porteuses; la figure 3 illustre un exemple d'estimation de la réponse impulsionnelle du canal obtenue en pratique, dans le cas d'un dispositif ne mettant pas en oeuvre les moyens de seuillage de l'invention; la figure 4 est le schéma synoptique des moyens de seuillage d'un dispositif tel que représenté en figure 1, dans le cas ou le seuil dépend du niveau de bruit et de la réponse impulsionnelle du canal de transmission; la figure 5 présente une comparaison des résultats obtenus dans le cas d'une démodulation différentielle et dans le cas de démodulations cohérentes avec un élément de référence sur quatre porteuses et un élément de référence sur huit porteuses. Les différents aspects du mode de réalisation qui va être décrit ci- après concernent plus particulièrement la réception de radiodiffusion sonore numérique vers des mobiles. Toutefois, il est clair que le principe du dispositif de démodulation cohérente de signaux numériques à haut débit de l'invention peut s'appliquer à tout type de récepteurs, dès lors que la diffusion des données se fait sous la forme de données numériques entrelacées en temps et en fréquence et que les données comprennent des éléments de référence. Ce dispositif s'applique notamment, mais non exclusivement, à la réception de signaux émis selon le procédé décrit dans la demande de brevet FR 90 01491. 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. C'est notamment le cas du procédé de diffusion COFDM Selon ce procédé, le signal émis est constitué par une suite de symboles de modulation formant un multiplex de N porteuses orthogonales. Soit {ff} l'ensemble des fréquences porteuses considérées, avec: fk = k/ts, k = O à N-1 o t 6 représente la durée allouée à un symbole de modulation. On peut alors définir la base de signaux élémentaires suivante: Y (t) avec k = O à N-1, j = à + o: f'i k (t) = gk (t jts) avec: valeurs dans un alphabet fini, représentant le signal de données émis. Le signal OFDM associé s'écrit alors: + N-1 x(t) = Re le 2 i( E CJ,k Y j,k(t))l j= k=O Pour s'affranchir de tout problème de sélectivité en fréquence du canal, un intervalle de garde de durée A est insérée devant chaque signal Ytk(t), afin d'absorber le brouillage intersymboles. ts représente désormais la durée du signal utile, Ts = t 6 + A, la durée du symbole, et A la durée de l'intervalle de garde. On définit alors les signaux émis par la relation Y-j,k(t) = gk(t-j Ts) avec: - reçu sur chaque porteuse k et à chaque instant j. De façon à permettre une démodulation cohérente, 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: Hik =/o,k ei:jk pour tous les instants j et toutes les fréquences k. Une façon avantageuse de le mettre en oeuvre consiste à utiliser certaines porteuses judicieusement équi-réparties dans le domaine temps- fréquence, comme fréquences pilotes de référence de phase et d'amplitude Ceci est tout à fait possible avec le système COFDM, du fait de la nature bidimensionnelle du signal émis. Cela permet d'obtenir une estimation des Hjk pour certains instants et certaines fréquences correspondant à l'insertion de ces références. L'estimation de la réponse du canal Hk à tout instant j T 8 et pour toutes les fréquences fk est ensuite obtenue par filtrage d'interpolation. Ce filtrage numérique peut être réalisé classiquement par un produit de convolution du signal d'entrée et de la réponse impulsionnelle du filtre Chaque valeur de sortie est alors remplacée par la somme pondérée de ces voisins. Il existe une autre façon avantageuse d'effectuer ce filtrage, qui utilise la propriété selon laquelle la Transformée de Fourier d'un produit de convolution est égale au produit des transformées L'opération nécessite alors une Transformée de Fourier Rapide (TFD) directe (respectivement inverse), un fenêtrage (pondération) du signal que l'on veut interpoler, et une TFD inverse (respecti- vement directe). L'invention concerne plus particulièrement cette seconde méthode En effet, une des propriétés intéressantes de la TFD est que la transformée d'un produit de convolution est égale au produit des transformées Il apparaît donc que le nombre d'opérations à effectuer par cette méthode est inférieur à celui nécessaire au filtrage RIF (filtre à réponse impulsionnelle finie), pour un rendement équivalent Par ailleurs, elle n'impose que peu de modification des décodeurs, des moyens de calcul d'une TFD étant déjà présents. L'insertion des porteuses de référence, à raison d'une référence toutes les R porteuses, permet au récepteur de disposer d'une estimation, entachée de bruit, de la réponse fréquentielle sous- échantillonnée du canal représentée par: H(\) pour \ 2 = n R n: 0 (N- 1)/R avec: R 1: nombre de porteuses entre deux références Il est donc nécessaire de déterminer le produit de convolution suivant, qui correspond au signal de sortie du filtrage de H(,)) par un filtre à réponse impulsionnelle de support fini A^ N-1 H(k) = z H'()) F(k,) ( 1) O=o avec H'(,) = H si) si = n R H'(,,)) = O sinon et F(k) : coefficients du filtre passe-bas d'interpolation. Les N échantillons estimés H(k) sont obtenus par une somme pondérée des échantillons H(/) voisins. Etant donnés les suites de N points H'()) (réponse du canal dans le domaine fréquentiel) et F(\) (réponse du filtre dans le domaine fréquentiel), dont les Transformées de Fourier inverses sont respectivement h'(n) et f(n), la Transformée de la convolution circulaire ( 1) peut s'écrire h(n) = h'(n) f(n) avec 'i(n) TFD q sur N points de H(k) h'(n) TFD -1 sur N points de H'(k) f(n) TFD -1 sur N points de F(k) Cette méthode d'interpolation nécessitera donc les trois opérations successives suivantes: TFD inverse pour obtenir les h'(n) et f(n) à partir des H'(k) et F(k) (passage du domaine fréquentiel au domaine temporel). Produit de h'(n) par f(n). TFD directe pour obtenir les H(k) à partir des h(n) (passage du domaine temporel au domaine fréquentiel). Il est souhaitable de choisir pour N une puissance de 2 afin de réaliser facilement les opérations de TFD. La figure 1 présente le schéma d'un dispositif de démodulation mettant en oeuvre de tels moyens d'estimation de la réponse du canal. Le signal reçu échantillonné yn est démodulé de façon classique, par une transformation de Fourier discrète (TFD) 11, produisant, dans le domaine fréquentiel les échantillons: Yk = Hk À Ck + Nk, k variant de O à N 1. Les échantillons définitifs Ck sont obtenus par projection, dans un module 12, des valeurs Yk sur les estimations Hk de la réponse fréquentielle du canal. L'interpolation de cette réponse fréquentielle est obtenue de la façon r",J r UO r. suivante Un module 13 réalise l'extraction des N/R échantillons Hk (Ho, Hr' r- rv. H 2 r Hn r) correspondant aux éléments de référence de l'ensemble des échantillons Yk, et l'insertion de (N N/R) éléments fictifs nuls entre ces références, de façon à obtenir N éléments, selon la séquence suivante: Ho, 0, 0, Hr, 0,, 0, Hnr. Par une transformation inverse 14 à N points, on obtient dans le domaine il temporel les valeurs h'n correspondant à l'interpolation dans le domaine temporel du sous-échantillonnage Hk de la réponse fréquentielle du canal Plus précisément, la suite des N valeurs h'n obtenues après la TFD -1 14 constitue une estimation de la réponse impulsionnelle du canal. Si on appelle h(n) la transformée inverse calculée sur N points de H(k), on montre que, en ne tenant pas compte du bruit, h'(n) obtenue à partir de la TFD inverse de H'(k) peut s'exprimer de la manière suivante: h'(n) = h(n) + h(n + N/R) + h(n + ( 2 N/R)) + h(n + (R 1)N/R) o les indices sont pris modulo N. Les figures 2 A et 2 B présentent respectivement deux exemples d'estima- tions de la réponse impulsionnelle du canal h(n) (correspondant à l'hypothèse o toutes les porteuses sont des porteuses de référence), et h'(n) (correspondant à un sous-échantillonnage de la réponse fréquentielle du canal, c'est-à-dire à l'utilisation d'une porteuse de référence toutes les R porteuses). Il apparaît clairement, en figure 2 B, que l'insertion d'éléments fictifs nuls entre les éléments de référence entraîne une répétition de la réponse impulsion- nelle De façon à obtenir une estimation 51 correspondant à celle de la figure 2 A, et donc éliminer la partie 52, il est nécessaire de fenêtrer temporellement. On montre que l'estimation h'(n) de la réponse impulsionnelle du canal ne présente pas de recouvrement, si et seulement si l'étalement max des différents échos vérifie: A max d'un signal, appliqué ici au cas particulier d'un signal complexe dont la Trans- formée de Fourier est dissymétrique. Il est donc nécessaire d'appliquer une fenêtre f(n) 15 de pondération, définie par: f(n) = 1 pourn:0,,(N/R -1) f(n) = O pour n: N/R, N 1 Les N échantillons fournis par la fenêtre 15 sont ensuite transformés par une transformation de Fourier discrète 16 à N points, fournissant dans le domaine fréquentiel l'estimation Hk de la réponse du canal. Dans un autre mode de réalisation, il est possible de ne pas insérer d'échantillons nuls parmi les échantillons Hk extraits par le module 13 d'extraction. Dans ce cas, le module 14 effectue une transformation de Fourier inverse seulement sur N/R points Cette seconde solution présente l'avantage de nécessiter un nombre d'opérations moins important. De plus, dans ce cas, le fenêtrage 15 correspond à un simple ajout aux N/R termes h'n d'une série de N N/R échantillons nuls L'application sur le signal h'n d'une fenêtre temporelle rectangulaire suivie de la TFD, peut être interprétée comme un échantillonnage parfait du signal (si les conditions de Shannon sont réunies) dans le domaine fréquentiel En augmentant ainsi la durée de l'enregistrement, en le portant de N/R à N avec des échantillons com- plémentaires nuls, on obtient une analyse spectrale plus fine. Selon les modes de réalisation, on peut choisir de limiter le nombre d'opérations, en adoptant cette seconde méthode, ou de privilégier l'utilisation d'un format de FFT standard, en mettant en oeuvre la première méthode présentée. Enfin, il est clair que l'ordre des transformations TFD -1 et TFD peut être inversé Il est tout à fait possible que le module 14 effectue une transformation directe, et le module 16 une transformation inverse. La figure 3 est un exemple d'estimation de la réponse impulsionnelle du canal, tel qu'on peut l'obtenir en pratique On obtient ainsi un tableau de N/R valeurs complexes que l'on peut graduer en retard Ce tableau comporte M raies distinctes auxquelles s'ajoute un bruit complexe Gaussien centré. Cette estimation de la réponse du canal est donc fortement entachée de bruit ce qui nuit à la qualité de l'interpolation. Aussi, si le gain en faveur de la démodulation cohérente par rapport à la démodulation différentielle est théoriquement de 3 d B lorsque la réponse du canal est estimée de manière parfaite, il n'est en pratique que de l'ordre de 0,5 d B. Selon une caractéristique essentielle de l'invention, la réponse impulsion- nelle du canal est traitée de façon à limiter l'influence du bruit Il est ainsi prévu un module 17 du seuillage de cette réponse, assurant l'élimination systématique de tout échantillon au-dessous d'un certain seuil Ainsi que cela est précisé par la suite, ce seuil peut être fixe ou adaptatif. Il est clair en effet que toute raie 21 inférieure au niveau de bruit 22 est totalement inexploitable L'invention propose donc de placer sur ce signal un seuil 23, éliminant tout signal dont le module est inférieur à ce seuil, et en particulier le bruit 22, et ne conservant donc que les raies utiles 24 A, 24 B, 24 c. Le seuillage est fait avantageusement avant l'opération de fenêtrage, notamment dans le cas de la seconde solution décrite En effet, il y a alors moins d'échantillons à traiter Toutefois, il peut également être réalisé entre les opérations de fenêtrage et de transformation directe. Ce procédé s'avère particulièrement performant avec des signaux du type des signaux COFDM En effet, l'information utile se trouve répartie sur un nombre relativement réduit de raies Elle est donc en très grande majorité au-dessus du seuil, et est conservée En revanche, l'essentiel du bruit est éliminé. Avantageusement, on choisit un seuil variable, notamment en fonction du niveau de bruit. La figure 4 présente un schéma synoptique de tels moyens de seuillage, selon un mode de réalisation préférentiel de l'invention. Le seuillage 17 est donc effectué sur les échantillons h'n, en fonction d'un seuil variable 32 Un module 33 de calcul détermine la valeur du seuil, en tenant compte de l'estimation U 2 de la puissance du bruit fournit par un module 34 d'estimation. Il est à noter que, dans les décodeurs COFDM, cette information a 2 est déjà disponible Le dispositif de l'invention ne nécessite donc pas d'importants moyens de traitement, mais exploite au contraire les moyens et les informations présents dans chaque décodeur, qu'il s'agisse du calcul de la FFT ou de l'estima- tion de î 2. Cette estimation peut par exemple être obtenue selon le procédé décrit dans le brevet FR 88 15216, déposé le 18 11 88 au nom des mêmes déposants, qui consiste à profiter de l'absence de signal pendant un temps de symbole, pour effectuer une analyse spectrale du bruit. On a pu constater que, dans les exemples de mise en oeuvre présentés plus loin, le seuil optimal se situe entre 5 a et 6 a, a étant l'écarttype du bruit. Le module 33 de calcul du seuil peut également tenir compte d'une estimation de la réponse impulsionnelle du canal, et en particulier dunombre de raies significatives En effet, plus il y a de raies, plus la puissance est répartie. Cette information est fournie par un module 35 d'estimation de la réponse impulsionnelle A nouveau, le module 35 est déjà présent dans les décodeurs COFDM, et utilisé pour la synchronisation. Le module 17 de seuillage peut par exemple être un comparateur, ou un circuit de polarisation. De nombreuses autres configurations peuvent bien sûr être mises en oeuvre Ainsi, on peut adjoindre aux moyens décrits un module de décision, ne mettant en oeuvre le seuillage que si le seuil dépasse une certaine valeur de référence. Il est également possible de choisir un seuil fixe, non adaptatif, qui peut alors être égal à environ 5 fois un "a critique"' as Dans ce cas, a 3 peut par exemple caractériser une puissance de bruit correspondant à un TEB d'environ 104. On présente ci-dessous quelques résultats numériques obtenus avec un dispositif selon l'invention, en relation avec la figure 5. Dans cet exemple, on utilise la technique de modulation COFDM Le nombre de porteuses N du multiplex est 512 Chaque symbole de longueur T 3 = ps à une durée utile t 3 = 64 Us Chaque porteuse est modulée en modulation de phase à 4 états. On sait que si on utilise une référence toutes les R porteuses, l'informa- tion sur la réponse du canal sera préservée tant que L Zmax est l'étalement maximum de la réponse impulsionnelle du canal. Cette réponse impulsionnelle peut être modélisée par une distribution exponentielle de la forme: i P( 6) = efi pour Z> O avec Co: moyenne et écart-type des retards. La figure 5 compare trois types de démodulation démodulation différentielle 41; démodulation cohérente avec une référence sur 4 porteuses 42; démodulation cohérente avec une référence sur 8 porteuses 43; pour un écart-type des retards Z% égal à 1 As. Ces courbes tiennent compte de la perte en efficacité en puissance ( 10 log 1/R) due à l'insertion des porteuses de référence (courbe paramétrée en Energie par bit utile sur No). On constate donc qu'on obtient, à l'aide du dispositif de l'invention un gain compris entre 1,6 d B (pour R = 4) et 2 d B (pour R = 8) par rapport à la démodulation différentielle. En d'autres termes, le traitement vis à vis du bruit de l'estimation de la réponse impulsionnelle du canal avant interpolation, permet de se situer à environ 1 d B de la courbe obtenue en démodulation cohérente avec une estimation du canal parfaite. Il est clair que le dispositif de l'invention présente encore un intérêt plus important lorsque le nombre d'états de constellation croît. Par ailleurs, l'invention propose un procédé de diffusion améliorée, permettant également de faciliter une démodulation cohérente On a vu en effet qu'il était nécessaire de bien distinguer les éléments de référence du bruit. Une autre solution pour obtenir ce résultat est d'augmenter à l'émission la puissance des porteuses de référence par rapport aux porteuses transportant de l'information Le nombre de porteuses de référence N/R étant relativement faible par rapport au nombre total N, la diminution incidente de l'efficacité en puissance du système est donc minime. Bien sûr, le procédé de diffusion et le dispositif de démodulation de l'invention peuvent être mis en oeuvre simultanément, pour améliorer encore la qualité de démodultation Ils peuvent également être utilisés indépendamment. REVENDICATIONS 1 Dispositif de démodulation cohérente d'un signal numérique du type constitué d'éléments numériques répartis dans l'espace temps-fréquence et émis sous la forme de symboles constitués d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu desdits éléments numériques et diffusées simultanément, ledit signal numérique comprenant des éléments de référence, de valeur et de position dans ledit espace temps-fréquence connues dudit dispositif de démodula- tion, dispositif caractérisé en ce qu'il comprend des moyens d'estimation par transforma- tion de Fourier de la réponse fréquentielle du canal de transmission à tout instant, réalisant la transformation ( 14) des échantillons reçus correspondant à des éléments de référence du domaine fréquentiel dans le domaine temporel, la multiplication ( 15) dans le domaine temporel desdits échantillons transformés par une fenêtre temporelle rectangulaire (f%) et la transformation ( 16) inverse après ladite multiplication des échantillons obtenus du domaine temporel dans le domaine fréquentiel, et en ce que lesdits moyens d'estimation comprennent des moyens ( 17) de seuillage desdits échantillons dans ledit domaine temporel, assurant l'élimination systémati- que des échantillons inférieurs à un certain seuil. 2 Dispositif selon la revendication 1, caractérisé en ce que lesdits moyens de seuillage tiennent compte pour le calcul dudit seuil de la valeur a 2 de la puissance du bruit affectant ledit canal de transmission. 3 Dispositif selon la revendication 2, caractérisé en ce que ledit seuil a une valeur comprise entre 5 a et 6 u. 4 Dispositif selon l'une quelconque des revendications 1 à 3, caractérisé en ce que lesdits moyens ( 17) de seuillage tiennent compte pour le calcul dudit seuil d'une estimation de la réponse impulsionnelle dudit canal de transmission. 5 Dispositif selon la revendication 1, caractérisé en ce que ledit seuil est fixe. 6 Dispositif selon l'une quelconque des revendications 1 à 5, caractérisé en ce que lesdits moyens ( 17) de seuillage sont placés en amont desdits moyens ( 15) de fenêtrage. 7 Dispositif selon l'une quelconque des revendications 1 à 6, caractérisé en ce que ladite transformation ( 14) des échantillons reçus est une transformation de format égal au nombre M d'éléments de référence par symboles émis, et en ce que ladite multiplication ( 15) par une fenêtre temporelle réalise simplement l'ajout d'une série de (N M) zéros aux M éléments de référence transformés. 8 Procédé de démodulation cohérente d'un signal numérique du type constitué d'éléments numériques répartis dans l'espace temps-fréquence et émis sous la forme de symboles constitués d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu desdits éléments numériques et diffusées simultanément, ledit signal numérique comprenant des éléments de référence, de valeur et de position dans ledit espace temps-fréquence connues dudit dispositif de démodula- tion, caractérisé en ce qu'il comprend une étape d'estimation par transformation de Fourier de la réponse fréquentielle du canal de transmission à tout instant, comprenant les étapes de: la transformation ( 14) des échantillons reçus correspondant à des éléments de référence du domaine fréquentiel dans le domaine temporel, la multiplication ( 15) dans le domaine temporel desdits échantillons transformés par une fenêtre temporelle rectangulaire f 1,, la transformation inverse ( 16) après ladite multiplication des échantillons obtenus du domaine temporel dans le domaine fréquentiel, et comprenant de plus une étape de seuillage ( 17) desdits échantillons dans ledit domaine temporel, assurant l'élimination systématique des échantillons inférieurs à un certain seuil. 9 Procédé de diffusion d'un signal numérique du type constitué d'éléments numériques répartis dans l'espace temps-fréquence et émis sous la forme de symboles constitués d'un multiplex de N fréquences porteuses orthogonales modulées par un jeu desdits éléments numériques et diffusées simultanément, ledit signal numérique comprenant des éléments de référence, de valeur et de position dans ledit espace temps-fréquence connues dudit dispositif de démodula- tion, caractérisé en ce qu'il comprend des moyens d'émission desdites porteuses attribuant sélectivement aux porteuses portant des éléments de référence une puissance supérieure à celle utilisée pour les porteuses portant des éléments d'information utile. Emetteur d'un signal numérique, mettant en oeuvre le procédé de 0 diffusion de la revendication 9. A: