FR 2604316 A2 19880325 FR 8613271 A 19860923 Procédé et installation de communication de donnéep numériciues vers des mobiles L'invention concerne la diffusion de données sous forme numérique en milieu perturbé par des parasites et/ou par l'existence de trajets multiples. Elle concerne plus particulièrement les procédés et installations permettant la radiodiffusion sonore numérique, ce qui n'exclut pas la transmission d'autres types d'informations. Mais l'adaptation à la radio diffusion sonore interdit les procédures de répétition et, en conséquence, implique une grande robustesse du système de transmission. La demande de brevet FR 86 09622 décrit et revendique un procédé de diffusion de données provenant d'une source qui comprend la mise des données sous forme d'échantillons numériques, le codage des bits des échan tillons à l'aide d'un code redondant conduisant à des symboles de modulation A, l'émission desdits symboles en les répartissant en parallèle dans l'espace fréquence 1-temps _t avec entrelacement temporel et émission simultanée sur M voies de diffusion à des fréquences différentes (M étant supérieur à 1). En mettant en oeuvre ce procédé, deux éléments A successifs ne sont jamais émis à des instants successifs à la même fréquence. on voit que l'on répartit le débit numérique total d'un programme donné sur M voies de diffusion, c'est-à-dire sur M porteuses. Si l'on doit diffuser simultanément L programmes, ayant chacun M porteuses, le récepteur à donc à traiter M porteuses parmi N = L.M porteuses pour assurer la réception d'un programme. A titre d'exemple représentatif de la radiodiffusion sonore numérique vers des mobiles, on peut chercher à transmettre simultanément L programmes en stéréo phonie à un débit de l'ordre de 2 x 200 kbits/s avec M porteuses par programme dans un canal télévision de 8 MHz. Dans le mode particulier de mise en oeuvre décrit à titre d'exemple dans la demande de brevet principal, un démodulateur élémentaire muni d'un filtre d'entrée ajustable pour permettre la sélection doit être prévu sur chacune des M voies affectées à un même pro gramme, ce qui se traduit par une complexité excessive lorsque le nombre de voies est élevé. on connait par ailleurs depuis longtemps les problèmes posés par la transmission numérique dans canal sélectif, c'est-à-dire dans lequel le débit numérique à transmettre est supérieur à l'inverse de l'étalement de la réponse impulsionnelle du canal (habituellement d'environ 5 ps). Lorsqu'on souhaite obtenir une efficacité spectrale élevée, on utilise alphabet de symboles ayant un spectre étroit et une fonction d'autocorrélation large, c' est-à-dire un codage connu sous le terme OFDM, abréviation de la désignation anglo-saxonne *Orthogonal Frequency Division Multi- plexing". Cette utilisation implique, à la réception, l'utilisation de transformées de Fourier discrètes. En conséquence, elle n'a été mise en oeuvre jusqu'ici que pour des débits et un nombre de canaux relativement faibles, par exemple pour la transmission de données dans la bande 60-104 kHz, comme solution de remplacement aux techniques de codage à réponse partielle décrites dans la Recommandation V37 du C.C.I.T.T. On pourra trouver une description de ce procédé dans divers documents, et notamment dans "Data transmission by frequency division multiplexing using the discrete Fourier transformu par S. S. Weinstein et autres, IEEE Com., Vol. Com. 19, No. 5, octobre 1971 et dans mAn orthogonally multiplexed QAM system using the discrete Fourier transformu par B. Hirosaki dans IEEE COM., Vol. COM. 29, No. 7, juillet 1981. L'invention vise à permettre la radiodiffusion sonore numérique vers des récepteurs, éventuellement mobiles, en mettant en oeuvre le procédé de la demande de brevet principal, mais sans exiger, à la réception, la présence d'un nombre de démodulateurs égal au nombre de voies élémentaires, et ce tout en permettant d'at teindre une efficacité spectrale élevée, de l'ordre de 1 bit/Hz. Dans ce but, l'invention propose notamment un procédé caractérisé en ce que les symboles sont émis par modulation d'amplitude de porteuses en quadrature- et en ce qu'on démodule simultanément, à la réception, les M voies, et elles seules, correspondant à un même pro gramme par application de la transformée de Fourier discrète simultanément à ces voies. L'invention propose également une installation de diffusion de données permettant de mettre en oeuvre le procédé ci-dessus défini. Cette installation comprend : - un émetteur destiné à diffuser plusieurs programmes, ayant : un codeur convolutionnel fournis sant, à partir des échantillons de données à diffuser, des symboles numériques pour chacun des programmes ; des moyens d'entrelacement des symboles dans le temps et de distribution des symboles entre plusieurs voies à des fréquences différentes, les fréquences affectées aux différents programmes étant intercalées, et - un récepteur comprenant des moyens de désen- trelacement et des moyens de décodage et de sélection de programme par transformée de Fourier discrète appliquée uniquement aux voies correspondant à un programme déter miné. L'invention sera mieux comprise à la lecture de la description qui suit de modes particuliers de réalisation donnés à titre d'exemples non limitatifs. La description se réfère aux dessins qui l'accompagnent, dans lesquels - la Figure 1 est un schéma de principe montrant une répartition possible des fréquences utilisées pour la transmission de L programmes, chacun sur M fré- quences, - les Figures 2A et 2B sont des schémas de principe montrant la réalisation de transformées de Fourier discrètes, respectivement avec L et M 4 et vec L et M - la Figure 3 est un synoptique montrant la constitution de principe d'un récepteur de mise en oeuvre de l'invention. On remarquera tout d'abord que, dans un canal de diffusion de télévision de 8 MHz, il est possible de loger N = 256 porteuses espacées de 31,25 kHz. Les fréquences allouées au premier programme sont choisies parmi les M fréquences fl,l ; f2,1 ; ... ; fM,l régu lièrement réparties-dans le spectre et entrelacées avec les fréquences affectées aux autres programmes pour accroître l'immunité aux bruits présentant un étalement en fréquence restreint. Dans la pratique, et pour des problèmes de raideur de flanc du filtre de réception, on sera généralement amené à adopter un nombre de porteuses effectivement utilisées inférieur à 256 et, par exemple, à adopter 14 porteuses pour chacun de 16 programmes, ce qui correspond à 224 fréquences seulement parmi les 256 disponibles. La mise en oeuvre de l'invention implique, d'une part, le choix d'un procédé de modulation particulier, d'autre part, la démodulation simultanée de toutes les porteuses correspondant à un même programme (mais de ces porteuses seulement) par transformée de Fourier dis crète. Ces deux éléments de l'invention seront succes sivement décrits. Les données sont transmises en modulation d'amplitude de porteuses en quadrature, c'est-à-dire par le procédé couramment dénommé QAM, dont un mode de réalisation est décrit dans l'article de Weinstein et autres déjà mentionné. Suivant l'invention, chaque sym bole élémentaire s(t) véhicule 2N bits sous la forme : ailleurs : s(t) := 0 Dans la formule (1) ci-dessus, fk f 0 Le temps à constitue un intervalle de garde entre symboles. Si on désigne par T. la durée d'un symbole, le signal émis peut s'écrire S(t) s(t-p T S). Les données a(k) et b(k) de la formule (1) dépendent alors de l'indice p. Si le canal de transmission a une fonction de transfert H(f), en supposant que sa réponse impulsion- nelle a un support inférieur à A, on peut écrire que la contribution du canal au signal est h(t) 0, pour t et t > à. En d'autres termes, l'interférence entre symboles d'ordre p-1 et d'ordre p n'affecte ce dernier que sur un temps au plus égal à A. On pourra négliger l'intervalle de temps [-à,o] qui a pour seul but d'absorber l'inter férence entre symboles en protégeant l'intervalle utile [0,NT]. Le signal x(t) reçu par le récepteur pour p = 0, pendant l'intervalle de temps [0,NTI s'écrit N-1 x(t) E p k [a(k) cos (2w fk+wk) b(k) sin (2w fk t+Wk)] k=O (2) Dans la formule (2), Pk exp (jwk ) est la réponse du canal à la fréquence f k : Wk est le déphasage intro duit par la ligne et Pk est un facteur multiplicatif. Si maintenant on pose X(k) [a(k) + i b (k)] p k exp (jWk) on a pour 0 t La démodulation doit permettre de retrouver les valeurs successives X(k) à partir du signal reçu x(t). Pour simplifier, on considèrera simplement la démodula tion du symbole élémentaire s d'ordre p = 0, la même opération devant être répétée à intervalles de temps T s L'intervalle de temps [-A,O] pourra être négligé. Il faut seulement remarquer au passage que à doit rester faible par rapport à T . puisque sa présence correspond à une consommation inutile de puissance émise. La démodulation est effectuée en projetant le signal x(t) reçu sur deux porteuses en quadrature à la fréquence f 0+1/2T, c'est-à-dire à la fréquence centrale du canal (T étant la fréquence bit). Le signal complexe -Z x(t) sur les deux porteuses en quadrature peut alors s'écrire : En échantillonnant le signal 'x"(t) à la fréquence bit 1/T, on obtient alors des échantillons successifs _W1, x(nT) L'équation (5) peut être mise sous une forme différente en faisant le changement de variable ,(_,)n %' x(n) IN] x (nT) On trouve alors c'est-à-dire : [x(n)] TFD_ 1 [X(k)] que l'on peut aussi bien écrire Conformément à l'invention, il n'y a pas, à la réception, démodulation par application de la transfor mée de Fourier discrète aux N points complexes, repré sentant autant d'échantillons, pendant la durée T s d'un symbole. Si le signal émis est composé de L programmes de M porteuses chacun, le résultat recherché et obtenu se limite à la transformée de Fourier discrète sur M porteuses parmi les N. La sélection du programme est donc utilisée en tant que moyen de décimation de la transformée de Fourier discrète. Pour réaliser de façon efficace cette trans formée de Fourier discrète ou TFD, on utilise l'un des algorithmes de transformée de Fourier rapide (ou TFR) bien connus. Si on utilise les notations courantes en' TFR, on peut poser : W exp (-2jw/n) et son calcul peut être décomposé en écrivant k = lr+s, où s désigne le numéro du programme traité parmi 1 et peut varier de s 0 à s 1-1 et où r est le numéro d'ordre de la porteuse et varie de 0 à m-1. n ml+m, c'est-à-dire que n peut être identifié par son quotient entier 1 par M et par le reste M. On voit qu'un programme donné correspond à une valeur particulière de s et correspond à M porteuses espacées de 1/MT, chacune affectée d'une valeur de r. On peut alors écrire Ken fonction de s et de r sous la forme : Le terme en , pouvant être désigné par q(s,m), devrait être calculé à l'aide de M transformées de Fourier discrètes sur L points si la sélection de programme ne s'effectuait qu'après transformation. Mais, du fait qu'il suffit de retenir un seul programme à la fois, la fonction q(s,m) n'est plus fonction que de m car elle correspond à une valeur unique de s. Le calcul des différentes valeurs q(s,m) est obtenu au moyen de N multiplications complexes seulement. Certaines de ces multiplications peuvent de plus être triviales. Connaissant q, on obtient X par M multiplica tions par les facteurs Wms et par une transformée de Fourier discrète effectuée sur M points seulement. Cette dernière opération peut être effectuée par l'un des algorithmes courants de transformée de Fourier rapide si M a une valeur telle qu'il est factorisable. La Figure 2A donne un exemple de mise en oeuvre de processus pour M 4 porteuses et L programmes. Dans ce cas, où L est une puissance de deux, l'échelonnement du calcul apparaît directement. La première étape fait intervenir des Opapillonsu sur toutes les m valeurs x échantillonnées. Chaque papillon" est un opérateur effectuant une addition, une soustraction et une multi plication complexe sur les deux valeurs d'entrée et fournit deux valeurs de sortie. Dans l'exemple donné en Figure 2A, N/2 = 16 papillons sont nécessaires dans cette première phase et fournissent des résultats qui sont rangés à des adresses déterminées dans une mémoire vive à écriture et lecture. Au cours de la seconde étape, l'invention prévoit de faire intervenir les papillons non pas sur tous les résultats de la première étape, mais uniquement sur ceux qui occupent les posi tions mémoires correspondant au programme à recevoir. Huit opérateurs seulement sont alors nécessaires (et seront d'ailleurs matériellement les mêmes que ceux utilisés dans la première étape). On progresseainsi dans le calcul en n'effectuant que les opérations con tenues dans le cadre en traits épais de la Figure 2A jusqu'à obtenir les quatre valeurs nécessaires, ici X 21 X 10, X18 et X26. On constate que le nombre total d'opé rateurs requis n'est plus que de 16 + 8 + 4 + 2 + 2 N. Parmi les algorithmes de transformée de Fourier rapide, il sera avantageux d'en utiliser qui sont basés sur des valeurs de N égales à des puissances de 2. Cela permet d'utiliser des structures très répétitives, telles que radix 2 ou 4, bien adaptées à la mise en oeuvre au moyen d'opérateurs câblés. Dans le cas illustré en Figure 2B, on réalise une TFD partielle sélectionnant un programme parmi L programmes correspondant chacun à M porteuses. De nouveau, il n est pas nécessaire de filtrer à l'entrée les fréquences utiles. Cette sélection s'effectue du fait même de la décimation au cours de la TFD. Si l'on revient à la formule (7), la méthode classique de calcul impliquerait le calcul de tous les q(s,m) par M six TFO sur cinq points. Mais une seule valeur de s étant à prendre en considération, par exemple s les valeurs requises de q(3,m) sont obtenues par N multiplications complexes, au cours de six TFD sur chacune cinq points, indiquée par la pre mière étape de calcul sur la Figure 2B. Les M résultats obtenus sont multipliés par les facteurs Wm3, dits "twiddlen. Enfin, il reste à effectuer une TFD sur M points seulement. Dans le cas de la Figure 2B où M est factorisable et s'écrit M 2 x 3, l'opération est effectuée par deux TFR successives, d'abord sur deux points, puis sur trois points. Les calculs nécessaires à la sélection du programme sont de nouveau uniquement ceux contenus dans le cadre en trait épais. Le résultat de la transformée de Fourier dis crète partielle sera constitué par M nombres complexes X(k) correspondant à la formule (6) ci-dessus, avec k = Lr+s où r est compris entre 0 et N-1 et s est le numéro du programme que l'on souhaite recevoir. En d'autres termes : X(k) [a(k)+ib(k)] Pk exp ('(Pk) où a et b sont les termes déjà mentionnés dans l'équation (1) et Pk exp (iwk) est la réponse du canal à la fréquence f k* Les variations de p k dans le temps sont toujours lentes par rapport à la durée T s d'un symbole, ce qui permet de récupérer la phase wk par des processus clas siques d'apprentissage et d'effectuer une démodulation cohérente. On pourrait également utiliser une démodulation différentielle, bien qu'elle soit moins performante lorsqu'on travaille sur quatre phases. Par contre, la démodulation différentielle donne des résultats satis faisants lorsqu'on travaille sur deux phases, mais l'emploi de deux phases seulement a un effet défavorable sur l'efficacité spectrale. on voit qu'après multiplication par exp _(jWk), on retrouve les données émises a et b affectées d'un facteur multiplicatif Pk- Pour retrouver les données réelles à partir du résultat de la multiplication, on met en oeuvre un processus d'évaluation et de décision. Pour améliorer les performances, il est souhaitable d'utiliser des codes correcteurs et de mettre en oeuvre une décision "douce" du type déjà décrit dans la demande de brevet FR 86 09622. L'entrelacement temps-fréquence, déjà défini dans la demande de brevet FR 86 09622, permet de placer le décodeur dans des conditions optimales d'indépendance entre les échantillons reçus. Cet entrelacement en blocs, convolutif ou pseudo-aléatoire, permet, en asso ciation avec le décodeur, d'éliminer les erreurs groupées dues à de l'évanouissement sur une fréquence donnée. La mise en oeuvre du procédé au niveau de l'émetteur ne soulève aucun problème particulier de réa lisation. L'émetteur utilisé peut avoir une constitution similaire à celle déjà décrite dans la demande FR 86 09622, si ce n'est qu'il met en -oeuvre une modulation d'amplitude sur porteuses en quadrature. Le récepteur peut comporter des étages d'entrée classique, fournis sant un signal à fréquence intermédiaire. Le décodeur 10 peut avoir la constitution de principe montrée en Figure 3. Ce décodeur 10 comprend un filtre de canal 12, constitué par un filtre passe-bande dont la largeur correspond à l'ensemble du spectre occupé par les N porteuses de diffusion. Le signal fourni par le filtre 12 est projeté sur deux porteuses en quadrature à la fréquence centrale du canal. Cette fréquènce est fournie par un oscillateur 14 qui attaque directement l'une des entrées d'un multi- plieur 161 et, par l'intermédiaire d'un déphaseur 18, -l'entrée d'un multiplieur 1612. Les deux multiplieurs reçoivent le signal de sortie du filtre 12. La sortie de chaque multiplieur attaque un convertisseur analogique- numérique 20A ou 2012 fournissant les échantillons d'entrée au circuit 22 de calcul de TFR. Au circuit 22 est associé un circuit 24 de sélection de programme qui détermine les adresses, dans la mémoire stockant les résultats intermédiaires contenue*dans le circuit 22, où se trouvent les'échantillons sur lesquels le calcul doit être poursuivi. Les échantillons fournis par le circuit 22 sont appliqués à un circuit à apprentissage 26 d'estimation de la phase et de démodulation : ce circuit peut être constitué à partir de composants intégrés disponibles dans le commerce, tels que les circuits TMS 320. Enfin, les données 28, quantifiées par exemple sur quatre bits, sont appliquées à un circuit de désentre- lacement, puis à un décodeur de Viterbi 32. Les circuits 30 et 32 peuvent avoir une constitution similaire à celle déjà décrite dans la demande FR 86 09622. Le démodulateur, constitué des éléments 14 à 20, peut être d'un type connu. Il faut d'ailleurs noter que l'on peut soit asservir l'oscillateur de transposition 14, soit effectuer l'opération de transposition après conversion analogique-numérique par multiplication complexe en bande de base effectuée immédiatement avant la TFR, cette dernière solution étant avantageuse dans le cas d'une intégration à grande échelle. A titre d'exemple de documents décrivant ces appareils de démo- dulation, on peut citer l'article de Baccetti et autres, "Full Digital Adaptive Equalization in 64QAM Radio Systems" dans les comptes rendus de la réunion ICC 86 à Toronto. On pourra également trouver une description d'un appareil comportant un égaliseur adaptatif dans "A Stepped Square 256QAM for Digital Radio System* par Toshihiko Ryu et Coll., publié dans les comptes rendus de la même réunion sous la référence 46-6, pp. 1477- 1481. REVENDICATIONS 1. Procédé de diffusion de données provenant d'une source, comprenant la mise des données sous forme d'échantillons numériques, le codage des bits des échan tillons à l'aide d'un code redondant conduisant à des symboles de modulation A, l'émission desdits symboles en les répartissant en parallèle dans l'espace fréquence f- temps _t avec entrelacement temporel et émission simulta née sur M voies de diffusion à des fréquences diffé- rentes, caractérisé en ce que les symboles sont émis par modulation d'amplitude de porteuses en quadrature et en ce qu'on démodule simultanément, à la réception, les M voies, et elles seules, par transformation de Fourier discrète simultanément à ces voies. 2. Procédé selon la revendication 1, caractérisé en ce qu'on diffuse simultanément un nombre L de pro-_ grammes de données supérieur à 1, les M voies correspon dant à un même programme étant intercalées en fréquence entre les voies correspondant aux autres programmes. 3. Procédé selon la revendication 2, caractérisé en ce que l'on effectue la démodulation par transformée de Fourier rapide, à l'aide d'opérateurs appliqués, dans une première étape, à la totalité des échantillons de données reçus, puis dans les étapes suivantes aux seuls échantillons nécessaires à la réception des M voies correspondant au programme choisi. 4. Installation permettant la mise en oeuvre du procédé selon la revendication 1, caractérisée en ce qu'elle comprend : - un émetteur destiné à diffuser plusieurs programmes, ayant : un codeur convolutionnel fournis sant, à partir des échantillons de données à diffuser, des symboles numériques pour chacun des programmes ; des moyens d'entrelacement des symboles dans le temps et de distribution des symboles entre plusieurs voies à des fréquences différentes, les fréquences affectées aux différents programmes étant intercalées, et - Un récepteur comprenant des moyens de désen- trelacement et des moyens de décodage et de sélection de programme par transformée de Fourier discrète appliquée uniquement aux voies correspondant à un programme déter miné.