FR 2765757 A1 19990108 FR 9708547 A 19970701 PROCÉDÉ ET DISPOSITIF DE MODULATION D'UN SIGNAL MULTIPORTEUSE DE TYPE OFDM/OQAM, ET PROCÉDÉ ET DISPOSITIF DE DÉMODULATION CORRESPONDANTS. 1 - Domaine de l'invention 1-1 Généralités Le domaine de l'invention est celui de la transmission ou de la diffusion de données numériques, ou de données analogiques et échantillonnées, destinées à être reçues notamment par des mobiles. Plus précisément, l'invention concerne la mise en oeuvre de signaux multiporteuses de type OFDM/OQAM. En d'autres termes, l'invention s'applique aux signaux de densité 2, voire plus. On sait que la modulation multiporteuse présente de nombreux intérêts, notamment lorsqu'elle est associée à un codage correcteur d'erreur et à un entrelacement. On a ainsi retenu la technique COFDM (Coded Orthogonal Frequency Division Multiplexing) pour la norme de diffusion audionumérique européenne (DAB) et pour la norme de diffusion de télévision numérique de terre (DVB-T). La technique COFDM présente un système d'égalisation particulièrement simple, à savoir l'utilisation d'un intervalle de garde, encore appelé préfixe cyclique. Ce dernier permet d'assurer un bon comportement face aux échos, mais au prix d'une perte relativement importante en efficacité spectrale. Ce problème est notamment discuté dans la demande de brevet française FR95 05455 (dans laquelle la modulation COFDM est appelée OFDM/QAM). Pour pallier ce problème, ce document de brevet présente une nouvelle technique d'implémentation de modulations multiporteuses de type OFDM/OQAM. On notera que les différents types de modulation discutés par la suite sont désignés de façon légèrement différente dans ce document antérieur et dans la présente demande de brevet. La table ci-dessous établit la correspondance: FR95 05455: Document présent: OFDM/QAM OFDM/QAM/OFDM OFDM/OQAM OFDM/OQAM/NYQUIST OFDM/OMSK OFDM/OQAM/MSK OFDM/ICYFA OFDM/OQAM/IOTA Le terme OQAM fait référence à la définition en langue anglaise Offset Quadratic Amplitude Modulation . Ceci exprime le fait que pour de telles modulations, un offset temporel existe entre l'émission de la partie en phase et celle de la partie en quadrature d'un symbole QAM. 1-2 Applications L'invention trouve des applications dans de très nombreux domaines, notamment lorsque une haute efficacité spectrale est souhaitée et que le canal est fortement non-stationnaire. Une première catégorie d'applications concerne la radiodiffusion numérique terrestre, qu'il s'agisse d'image, de son et/ou de données. En particulier, l'invention peut s'appliquer à la diffusion synchrone, qui génère intrisèquement des trajets multiples de longue durée. Elle s'applique également avantageusement à la diffusion vers des mobiles. Une autre catégorie d'applications concerne les radiocommunications numériques. L'invention peut trouver notamment des applications dans des système s de communication numérique vers des mobiles à haut débit. 2 - Rappels 2-1 Canal de transmission En environnement radiomobile, l'onde émise subit des réflexions multiples, et le récepteur reçoit donc une somme de versions retardées du signal émis. Chacune de ces versions est atténuée et déphasée de façon aléatoire. Le récepteur étant supposé en mouvement, l'effet Doppler agit également sur chaque trajet. La conjugaison de ces effets se traduit par un canal non stationnaire présentant des évanouissements profonds à certaines fréquences (canal sélectif en fréquence). Pour les applications qui retiennent notre attention, la bande de transmission est supérieure à la bande de cohérence du canal (bande pour laquelle la réponse fréquentielle du canal peut être considérée comme constante, sur une durée donnée), et il apparaît donc des évanouissements dans la bande, c'est à dire qu'à un instant donné, certaines fréquences sont fortement atténuées. 2-2 Description d'une modulation multiporteuse Une modulation multiporteuse est avant tout une modulation numérique, c'est-à-dire un procédé de génération d'un signal électromagnétique à partir de l'information numérique à transmettre. L'originalité, et l'intérêt, d'une telle modulation est de découper la bande limitée allouée au signal en sous-bandes. Sur ces sous-bandes, choisies de largeur inférieure à la bande de cohérence du canal, le canal peut être considéré comme constant pendant la durée de transmission d'un symbole, choisie inférieure au temps de cohérence du canai. L'information numérique à transmettre pendant cette durée est alors répartie sur chacune de ces sous-bandes. Cela présente notamment deux intérêts: - diminuer la rapidité de modulation (c'est-à-dire augmenter la durée symbole), sans modifier le débit transmis, - modéliser simplement l'action du canal sur chacune des sous-bandes multiplieur complexe. On notera que, en réception, un système peu complexe de correction des données reçues (division complexe par le canal estimé) permet de récupérer l'information émise sur chacune des porteuses de façon satisfaisante, sauf pour les porteuses ayant subi un évanouissement profond. Dans ce cas, si aucune mesure de protection de l'information n'est prise, les données véhiculées par ces porteuses seront perdues. Un système multiporteuse prévoit donc que la génération du signal électrique soit précédée de traitements numériques des données (codage correcteur d'erreurs et entrelacement). Le texte du brevet n" FR95/05455 décrit en détail les deux types de modulation multiporteuse existants. On rappelle ici succinctement leurs caractéristiques. 2-2-2 Notations Espacement entre deux porteuses adjacentes du multiplex de porteuses : v0. Espacement temporel entre deux symboles multiporteuses émis (temps symbole): 0. 2-2-3 La fonction prototype Le filtre de mise en forme de chacune des porteuses du multiplex est le même, il correspond à la fonction prototype caractérisant la modulation multiporteuse. Soit g(t) cette fonction prototype, le signal émis, à chaque instant nT0, sur la mème sous-bande de fréquence centrale vm, est am,nei#m,ne2i#vmtg (t-nt0). En bande basse, l'expression du signal émis, centré autour de la fréquence Mv0, est donc: Les fonctions eitPmfne2i2DnV tg(t~nTo) sont appelées les translatées "tempsfréquence" de g(t). Pour retrouver l'information transmise par chacune des sousporteuses, il faut choisir g(t) de sorte que ses translatées "temps-fréquence" soient séparables. Pour être sûr de cela, on impose que ces translatées soient orthogonales, au sens d'un produit scalaire défini sur l'ensemble des fonctions d'énergie finie (qui est un espace de Hilbert au sens mathématique). Cet espace admet deux produits scalaires possibles, qui sont: - le PS complexe (xly) I y > = x(t)y* (t)dt R - le PS réel (xly) = Ste(Jx(t)y (t)dtssJ R On définit ainsi deux types de modulation multiporteuse: - type complexe, ou encore OFDM/QAM: la fonction g(t) choisie garantit une orthogonalité de ses translatées au sens complexe (Exemple : OFDM, encore nommé OFDM/QAM/OFDM). Dans ce cas, #m,n = O et les données a" n sont complexes. - de type réel, ou encore OFDM/OQAM : la fonction g(t) choisie garantit une orthogonalité de ses translatées au sens réel (Exemples OFDM/OQAM/NYQUIST, OFDM/OQAM/MSK, OFDM/OQAM/IOTA). Dans ce cas, #m,n = (#/2)*(m+n) et les données a",n sont réelles. 2-24 Densité du réseau "temps-fréquence" Ces modulations étant destinées à transmettre de hauts débits, elles seront associées à des efficacités spectrales assez grandes, de l'ordre de 4 bits/s/Hz (Télévision Numérique). Le "mapping" des bits issus du codeur correcteur d'erreur sera ainsi de type QAM. Pour une modulation OFDM/QAM, les parties réelle et imaginaire d'un complexe issu de la constellation QAM sont transmises simultanément, tous les temps symbole Ts. Dans le cas d'une modulation de type OFDM/OQAM, elles sont transmises avec un décalage temporel (Offset QAM) d'un demi temps symbole (T,/2). Pour une même bande de transmission et un même nombre de sous-porteuses, il faudra donc, pour transmettre un même débit, que le rythme d'émission de symboles multiporteuses de type OFDM/OQAM soit deux fois plus rapide que celui de symboles multiporteuses de type OFDM/QAM. Ces deux modes de transmission de l'information sont caractérisés par la densité du réseau temps-fréquence d = l/(v0 lo) Les modulations multiporteuses de type OFDMIOQAM correspondent à une densité d = 2, celles de type OFDMtQAM correspondent à une densité d = 1. On peut donc noter qu'une modulation multiporteuse est caractérisée par: - la densité du réseau "temps-fréquence" sur lequel elle est définie, - la fonction prototype. La mise en oeuvre d'une modulation OFDMtOQAM de densité 2, et de la démodulation correspondante, demande une puissance de calcul importante, et une grande capacité de mémorisation. Cela sous-entend donc que les appareils correspondants seront relativement complexes et coûteux. 3 - Objectifs de l'invention L'invention a notamment pour objectif de pallier, ou réduire, ces inconvénients. Plus précisément, un objectif de l'invention est de fournir des techniques de modulation et de démodulation de signaux multiporteuses qui soient simples et peu coûteuses à mettre en oeuvre, par rapport aux approches directes. En d'autres termes, l'invention a pour objectif de fournir de telles techniques de modulation et de démodulation, qui limitent le nombre d'opérations à effectuer et la capacité de mémorisation nécessaire. 4 - Description de l'invention Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints selon l'invention à l'aide d'un procédé de modulation d'un signal multiporteuse de densité l/(v,.z,)=2, formé de symboles successifs, comprenant chacun M échantillons à émettre, et constitués d'un jeu de 2M fréquences porteuses orthogonales au sens réel, l'écart entre deux fréquences porteuses voisines valant v0 et l'écart entre les instants d'émission de deux symboles consécutifs, ou temps symbole, valant 0, chacune desdites fréquences porteuses étant modulée selon une même fonction prototype de modulation g(t) présentant une longueur de troncature de 2Lto, procédé comprenant, à chaque temps symbole, les étapes suivantes: - obtention d'un jeu de 2M coefficients complexes représentatifs de données àémettre; - calcul de 2LM combinaisons linéaires à partir desdits 2M coefficients complexes obtenus, les coefficients de pondération utilisés dans cesdites combinaisons étant représentatifs de ladite fonction prototype g(t), de façon à obtenir 2LM coefficients; - sommation desdits 2LM coefficients pondérés dans des emplacements mémoire prédéterminés respectifs d'une mémoire comprenant 2LM emplacements mémoire représentant 2L groupes de M sommes partielles distinctes, de façon à former progressivement, dans lesdits emplacements mémoire, sur une durée de 2LIo, M échantillons à émettre; - émission desdits échantillons à émettre. Ainsi, selon l'invention, les données à traiter sont stockées après la pondération, et non avant. I1 est ainsi possible de réduire la capacité de mémoire nécessaire, ainsi que le nombre d'opérations effectuées. Les échantillons à émettre sont construits progressivement, dans chaque champ de mémorisation. Selon un mode de réalisation avantageux de l'invention, un échantillon à émettre à l'instant jx0 + kzdM, noté 5k+jM s'écrit où: C0j à C,,j sont les 2M coefficients complexes générés entre les instants jTo et (j+l)o; ak q et I3k.q sont lesdits coefficients de pondération. Dans le cas d'une modulation OFDM/OQAM, on aura généralement: &alpha;k,q7 q = 0 pour q impair; ssk,q = pour q pair. Le nombre d'opérations effectuées est donc encore réduit de moitié. Dans un mode de réalisation préférentiel de l'invention, le procédé comprend, pour la génération d'un symbole d'indice j formé de M échantillons, les étapes suivantes: - obtention de 2M entrées réelles am,n représentatives d'un signal source; -pré-modulation de chacune desdites entrées réelles, produisant 2M coefficients complexes; -transformation de Fourier inverse desdits 2M coefficients complexes, produisant 2M coefficients transformés complexes C0,j à - pour chacun des M couples (Ckj, C+M)j) dedits coefficients transformés, calcul de 2L couples pondérés, les coefficients de pondération étant représentatives de ladite fonction prototype; - addition du résultat de chacune desdites 2LM valeurs pondérées au contenu de 2LM zones mémoire distinctes, de façon à construire progressivement les échantillons à émettre constituants les symboles j, (j+l), (j+2),... (j+2L-1); - émission de M échantillons correspondant aux M plus anciens contenus desdites zones mémoire puis mise à zéro du contenu de cesdites M zones mémoire. En général, lesdites étapes seront mises en oeuvre au rythme xJM des échantillons. Le contrôle des moyens de mémorisation est très simple. Ainsi, ladite étape d'émission peut être suivie d'une étape de mise à jour desdits emplacements mémoire, comprenant: - un décalage physique du contenu de chacun desdits emplacements mémoire, si ces derniers sont des éléments d'un registre à décalage ; ou - une mise à jour des adresses d'écriture et de lecture desdits emplacements mémoire, si ces derniers sont des éléments d'une mémoire RAM. Selon une caractéristique avantageuse de l'invention, lesdits coefficients représentatifs de données à émettre sont obtenus par la mise en oeuvre d'une transformation mathématique comprenant les étapes suivantes: - application d'une transformation de Fourier inverse réelle; - permutation circulaire du résultat de cette transformée inverse de NV2 coefficients vers la gauche; - multiplication par in de chacun desdits coefficients. Ainsi, il est possible d'obtenir des coefficients transformés complexes, à partir d'une PP li à entrées réelles. A nouveau, cela permet de limiter le nombre d'opérations effectuées. On peut encore simplifier les calculs, en modifiant légèrement l'équation du signal centré autour de la fréquence Mv,, de façon qu'il s'écrive: L'invention concerne également les dispositifs de modulation mettant en oeuvre un tel procédé de modulation. Selon un mode de réalisation particulier, ce dispositif comprend notamment: - des moyens de transformation mathématique délivrant lesdits coefficients représentatifs de données à émettre au rythme To/2M et dans l'ordre suivant (COJS CM+Ij). (CMIj, C2M-1 j); - 2LM-M emplacements mémoire de type RAM à écriture et lecture simultanées; - N multiplieurs complexes fonctionnant au rythme NT2LM, N valant 1, 2, 4,...ou 2L. De cette façon, la place mémoire est encore réduite. L'invention concerne encore un procédé de démodulation d'un signal reçu, correspondant à un signal émis multiporteuse de densité l/(v0.T0) = 2. Selon ce procédé, on reconstruit une estimation des 2M données réelles émises à un temps symbole donné, à l'aide des étapes suivantes: - échantillonnage dudit signal reçu à la fréquence d'échantillonnage to/M, délivrant M échantillons complexes reçus; - mémorisation de chacun desdits M échantillons complexes reçus dans un emplacement prédéterminé d'une mémoire d'entrée comprenant 2ML emplacements complexes, dans laquelle ont été préalablement mémorisés (2L-1)M échantillons reçus pendant les 2L-1 temps symbole précédents; - multiplication des 2ML valeurs contenues dans ladite mémoire d'entrée par des coefficients représentatifs de ladite fonction prototype; - repliement temporel, par sommation de 2M séries de L résultats de multiplication, de façon à obtenir 2M valeurs complexes ; - traitement desdites 2M valeurs complexes pour former lesdites estimations des 2M données réelles émises. Avantageusement, les 2M valeurs complexes issues de l'étape de repliement temporel entre les instants (j+2L-l)To et (j+2L)To s'écrivent: où: rk'+(j+q)M représente l'échantillon reçu à l'instant k'x0+(j+q')X0/M; a 'kq et ss'k,q sont lesdits coefficients de pondération. Le plus souvent, les calculs seront simplifiés, du fait que: - &alpha;'k,q' = 0 pour q' impair; - ss'k,q' = 0 pour q' pair. Selon un mode de réalisation préférentiel, ladite étape de traitement comprend les étapes suivantes - application d'une transformation mathématique inverse de celle effectuée lors de la modulation sur lesdites 2M valeurs complexes, délivrant 2M valeurs transformées; - correction de distorsions de phase et/ou d'amplitude dues au canal de transmission; - extraction de la partie réelle desdites valeurs complexes transformées. En général, lesdites étapes sont mises en oeuvre au rythme TM des échantillons. L'invention concerne également les dispositifs de démodulation mettant en oeuvre ce procédé. Ces dispositifs comprennent: - des moyens d'échantillonnage dudit signal reçu; - des moyens de mémorisation temporaire des échantillons complexes échantillonnés, comprenant 2ML emplacements complexes; - des moyens de multiplication desdits échantillons mémorisés par des coefficients de pondération représentatifs de ladite fonction prototype; - des moyens de repliement temporel, assurant la sommation de L résultats de pondération, de façon à obtenir 2M valeurs complexes; - des moyens de traitement desdites valeurs complexes, délivrant une estimation de 2M données réelles émises à chaque temps symbole. I1 est possible de réduire encore la place mémoire nécessaire dans dispositif, à l'aide de: - moyens de mémorisation comprenant 2ML-M emplacements mémoire complexes de type RAM à écriture et lecture simultanée; - N multiplieurs complexes fonctionnant au rythme NT2LM, où N vaut 1, 2, 4... ou 2L ; - moyens de transformation mathématique fonctionnant au rythme To/2Ms dont les entrées, R0j à R2M1j sont lues dans l'ordre (Rojs RMj). (R, RM+1j),... (RM-1,j, R2M-1 j) Enfin, l'invention concerne plus généralement un procédé de filtrage délivrant des séries de M valeurs complexes de sortie à intervalles réguliers, à partir de 2L séries de 2M valeurs complexes d'entrée, lesdites M valeurs complexes correspondant à une somme pondérée de 2L desdites valeurs complexes d'entrée à traiter, comprenant les étapes suivantes, pour chaque série de valeurs complexes d'entrée: - calcul de 2LM combinaisons linéaires à partir desdits 2M coefficients complexes obtenus, les coefficients de pondération étant issus de 2L filtres complexes ou réels de taille M, de façon à obtenir 2LM coefficients; - sommation de chacune des valeurs pondérées dans un emplacement mémoire prédéterminé, parmi un ensemble de 2ML emplacements mémoire contenant chacun une somme partielle, de façon à former progressivement lesdites valeurs de sortie, dans lesdits emplacements mémoire, sur une durée correspondant à la réception 2L séries de valeurs complexes d'entrée. Le terme filtrage doit bien sûr être pris ici dans son sens général de traitement, ou calcul, effectué sur des données. Ce traitement, qui comprend le calcul d'une somme pondérée, est réalisé progressivement. 5 - Description d'un mode de réalisation préférentiel 5-1 Liste des figures D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à 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, et des dessins annexés, parmi lesquels: - la figure 1 illustre de façon générale et simplifiée le procédé de modulation de l'invention (étape j), permettant de générer M échantillons; - la figure 2 illustre la construction progressive des M échantillons, dans le cas où L=4; pour la forme d'onde Iota; - la figure 3 illustre de façon plus détaillée le fonctionnement du procédé de modulation de l'invention, pour les instants j-l à j+2. - la figure 4 précise l'amorçage de la procédure de modulation de la figure 3, dans le cas où L = 4. - la figure 5 est un schéma de principe d'un circuit de IFT;T complexe connu en soi; - les figures 6A à 6C illustrent des architectures optimisées, mettant en oeuvre des FIFOs, et utilisant respectivement un seul multiplieur (figure 6A), L multiplieurs (figure 6B) ou 2L multiplieurs (figure 6C) ; ce procédé, dans le cas où une structure FIFO est choisie. la figure 7 présente un mode de réalisation optimisé de la FFT inverse, utilisant une FET à entrées réelles; - la figure 8 illustre le fonctionnement du procédé de démodulation de l'invention, dans le cas L=4; - la figure 9 présente le cas général de la démodulation, qui se déduit directement de la figure 8; - la figure 10 illustre une architecture de démodulateur correspondant; - les figures 1 1 et 12 présentent deux mises en oeuvre du filtrage de réception, dans le cas d'une structure FIFO, utilisant respectivement L et 2L multiplieurs. 5-2 Notations Espacement interporteuses : v0. . Durée intersymboles : 0. Densité du réseau: l/(vot0) = 2. .Bande allouée au signal: W= 2Mvo. . Fréquence d'échantillonnage: fe = l/Te = M/o. Longueur de troncature de la fonction prototype : 2lao. En théorie, la fonction prototype est de support temporel et/ou fréquentiel infini. Cependant, pour implémenter le filtre numérique correspondant, il faut tronquer cette fonction. C'est le cas en OFDM/OQAM/NYQUIST (support temporel infini) et en OFDM/OQAM/IOTA (supports temporel et fréquentiel infinis). Typiquement, pour la fonction Iota, L = 4 au minimum. Pour un échantillonnage à Te défini ci-dessus, le filtre numérique sera alors de longueur 2LM coefficients réels. indices des échantillons: Pour être cohérent avec les formules données par la suite, on note: la fonction prototype g(t), le signal émis s(t) ou le signal émis r(t) peuvent être substitués à x(t). 5-3 Algorithme de modulation 5-3-1 Principe Le signal bande basse centré autour de la fréquence Mvo = fe/2, s'écrit Un échantillon du signal s'écrit donc Avec les notations introduites ci-dessus (Formule (2)), on a, après calcul, Etant donnée la décroissance rapide de la fonction prototype, seuls ses échantillons indicés de 0 à 2ML-1 sont considérés comme non nuls. Nous devrons donc avoir 0 L'équation (5) devient Enfin, en posant q = j - n, nous obtenons la formule dont découle l'algorithme de modulation: 2L-1 ZM-I li ni (P+9M) S+jM = a ", ( j ç ) i t i i ) (~l)m(j-y+L)e e 2 M M q=O m=O ~= Rcmodulaon FFT Inverse 7 Ponderation par la fonction protcnpe cjr) 0 avec JEZ En notant Cm.j q les valeurs d'entrées pré-modulées, soit Cm,j~q = = m.(j-q) I m+(j-q) (-1)m(j-q+L) l'équation (7) devient: Cette dernière formule conduit à un algorithme de modulation en trois étapes principales: - Pré-modulation des données grâce à une multiplication complexe simple. - Transformée de Fourier inverse (par l'algorithme IFFT). - Filtrage par la fonction prototype. Ici, comme dans la suite de ce document, on entend par filtrage une opération de pondération des résultats de 2L FFT inverses par certaines valeurs de la fonction prototype, suivie d'une opération de sommation de ces coefficients pondérés. 11 s'agit en d'autres termes d'une combinaison linéaire. On présente ci-dessous deux modes de mise en oeuvre possibles. Le second est optimal et fait plus précisément l'objet de l'invention. Bien qu'en pratique, on travaille au rythme échantillon xo/Ms nous garderons la structure de bloc de M échantillons pour décrire ces modes de mise en oeuvre avec plus de clarté. 5-3-2 Architecture directe Au vu de la formule (8), il faut effectuer 2L'FEr de taille 2M complexes pour générer M échantillons (correspondant à la durée d'un symbole multiporteuse lo) Cependant, le résultat d'une IFFT intervient sur le calcul de 2L blocs consécutifs d'échantillons. Pour calculer les M échantillons courants, il ne faut donc calculer que la Transformée de Fourier inverse des 2M dernières données entrées dans le modulateur, les résultats des 2La1 autres IFl;T ayant été calculés aux étapes précédentes, et stockés en mémoire. L'algorithme de modulation comprend donc les étapes principales suivantes: - Pré-modulation des 2M entrées réelles, délivrant 2M complexes. - Transformée de Fourier inverse de taille 2M complexes (algorithme IFFT). - Stockage du résultat en réactualisant un buffer de taille 2L*2M complexes (contenant les résultats des 2L FFr inverses impliquées dans le calcul). - Filtrage de 2LM éléments du buffer de stockage par la fonction prototype. - Emission des M échantillons complexes ainsi calculés. La taille mémoire requise est donc: - une RAM de taille 2L*2M complexes (buffer d'entrée). - une ROM de taille ML+l réels (coefficients de pondération). (La fonction prototype étant choisie symétrique, sur les 2ML coefficients de pondération, seuls ML+ 1 sont distincts). Cette première façon de procéder laisse entrevoir un gaspillage de mémoire RAM. La deuxième architecture proposée montre qu'il est possible de réduire de plus que de moitié la taille de mémoire RAM nécessaire. Cette réduction s'accompagne d'une réduction du nombre d'opérations, et donc d'une augmentation de la vitesse de traitement. 5-3-3 Architecture optimisée n est possible d'optimiser l'architecture du modulateur, selon l'invention, tant au niveau du filtrage par la fonction prototype que de la FFT inverse. En effet, en analysant la formule (8), on peut vérifier que pour chacune des 2L bP1 impliquées dans le calcul du bloc courant de M échantillons, seuls M points sur 2M sont utilisés. I1 est alors possible de réduire de moitié la mémoire RAM nécessaire stockant les données servant aux calculs des différents blocs d'échantillons après le filtrage, et non avant. De plus, la structure particulière des données complexes en entrée de la FFT inverse (am,n (- 1)m(n+L) im+n) permet d'utiliser un algorithme de FFT inverse à entrées réelles. Afin de préciser ce procédé, nous allons développer la formule (8): Posons: Nous avons: avec 0 L'équation (10) traduit la construction de M valeurs complexes à partir de 2ML valeurs complexes. Elle peut s'écrire de façon plus générale: Dans le mode de réalisation décrit, # 0 si q est impair # gk+qM Si q est impair &alpha; k,q = gk+qM Si q est pair et ssk,q # 0 si q est pair Pour générer M échantillons selon cet algorithme de modulation, on peut donc procéder ainsi que cela est illustré en figure 1 (étape j): - Pré-modulation 11 des 2M entrées réelles. - Transformée de Fourier inverse 12 des 2M données complexes ainsi obtenues, de façon à générer des Ck.j et C+M.j - Pondération 13 (correspondant à l'application de la fonction prototype) du résultat de la Transformée de Fourier inverse par la fonction prototype : L opérations de pondération parallèles. Les L vecteurs de pondération, de taille 2M, ont pour coefficients: [g0,.., g2M-1], [g2M,..., g4M-1],..., [g2LM-2M,.., g2LM-1]. - Addition 14 de ces résultats de pondération au buffer de sortie de taille 2ML complexes. - Décalage 15 du buffer de sortie avec émission de M échantillons, correspondants aux M plus anciennes valeurs contenues dans le buffer. Un échantillon du signal à émettre représente une somme de 2L résultats de t pondérés. Chaque bloc de M cases mémoire consécutives du buffer de sortie contient M sommes partielles de 2L-m,1 termes chacune, où mbXOX varie de 1 à 2L (ioc =2L correspond au bloc "tout à zéro", dû à l'opération de décalage du buffer (étape 14) à l'instant (j-l)T0). Les 2ML éléments issus des L opérations de pondération parallèles sont ici ajoutés aux 2ML éléments du buffer. Après cette opération, les M sommes partielles du buffer correspondant à mblOz = 1 sont complétées, c'est à dire que les M échantillons courants sont calculés et peuvent donc être émis. Ce fonctionnement est décrit par la figure 2, dans le cas où L=4. Chaque ligne illustre la situation du buffer de construction des données à émettre, à un instant donné. n faut 2L temps symbole consécutifs pour construire progressivement un échantillon à émettre. La forme d'onde 2L représentée correspond à la fonction Iota. Elle est représentée par les 2L vecteurs de coefficients 22 [gk] à [gk+7M], où l'indice k varie deOàM-l. A chaque instant, les 2M coefficients en entrée sont multipliés (23) par les coefficients 22, puis additionnnés (24) chacun à une somme partielle. Les M sommes partielles complétées à l'étape 15 sont émises, le contenu du buffer est décalé de M cases mémoire (de sorte à assurer le bon ordre du calcul des M prochains échantillons) et M zéros sont insérés dans les M cases mémoire vacantes. La diagonale 25 illustre ainsi le calcul de Sk+jM. La figure 3 illustre de façon plus détaillée le fonctionnement de cet algorithme, pour les instants j-l à j+2. Si l'on considère l'instant j, les coefficients arnj à transmettre alimentent le module 31 de pré-modulation, qui fournit à la Wr inverse 32 les coefficients cmj. La FFT inverse délivre les C et Ck+Mj (l'indice k variant de 0 à M-l) qui sont soumis à la pondération 33 (L opérations de pondération en parallèle) pour délivrer les résultats 34, qui sont sommés dans un buffer de sortie 35. On a indiqué, sur la figure 3, le contenu exact de ces buffers de sortie. La figure 4 précise l'amorçage de cette procédure de modulation, dans le cas oùL=4. L'architecture du modulateur correspondant à l'algorithme ci-dessus présenté devra donc comporter: une ROM de taille ML+1 réels, contenant les coefficients du filtre. une RAM de taille 2ML complexes correspondant au buffer de sortie. un circuit de Wr complexe (réalisant une Wr inverse) de taille 2M. Pour augmenter la vitesse de traitement, on parallélisera les opérations de pondération en utilisant L RAMs de taille 2M associées à L multiplieurs, ou même 2L RAMs de taille M complexes associées à 2L multiplieurs, au lieu d'une RAM de 2ML complexes. La complexité du circuit de modulation est donc: - pour le filtrage: Afin de réaliser le filtrage, nous multiplions les résultats de la Wr inverse complexe par les 2ML coefficients de la fonction prototype, en réalisant L opérations de pondération en parallèle, dont nous ajoutons les résultats au buffer de sortie. Etant donné que les coefficient de la fonction prototype sont réels, nous avons (2x2ML) multiplications réelles et 2ML additions complexes ou 4ML additions réelles. La taille du buffer de sortie est alors 2ML mémoires complexes ou 4ML mémoires réelles. - pour la transformation lWr: Les résultats indiqués dans le tableau ci-dessous concernent un circuit de IFFT complexe conventionnel, dont le schéma de principe est donné figure 5. Elle comprend un buffer d'entrée 51, recevant 2M entrées, une unité de calcul 52, alimentée par des coefficients mémorisés en mémoire ROM 53 et des valeurs calculées stockées dans une mémoire RAM 54, qui délivre 2M sorties. Un module de contrôle pilote ces divers éléments. Ce tableau donne une estimation des opérations et matériels nécessaires pour la partie de modulation: Modulation Additions Multiplications RAM ROM FIFO (Réelle) (Réelle) (Réelle) (Réelle) (Réelle) Pré-modulation x x x x x L FFT inverse 6M(l+log2M) 4M(l+log2M) 8M 1 2M x Filtrage 4ML ; 4ML ~ 4ML ~ ML+1 x On notera qu'aucune multiplication ni addition n'est nécessaire pour la prémodulation car la multiplication complexe simple à effectuer à cette étape se traduit, au niveau de l'architecture, par des permutations de parties réelles et imaginaires, ainsi que des changements de signe. Selon une variante de l'invention, un gain supplémentaire de place mémoire peut être obtenu. n est en effet possible de n'utiliser que (2L-1) RAMs de M complexes (stockage total : 2ML-M complexes au lieu de 2ML). Pour ce faire, à l'étape j, on lit les M échantillons à émettre dans la RAM correspondante, et on écrit, au fur et à mesure, les M complexes Ck+M,jgk+(2L-1)M aux mêmes adresses. Pour réaliser cette opération de filtrage, on pourra utiliser, au choix, des RAMs ou des FIFOs. Les figures 6A à 6C illustrent ce procédé, dans le cas où une structure FIFO est choisie. Dans le cas de la figure 6A, un seul multiplieur 61 est mis en oeuvre. n assure la multiplication des données délivrées par la FFT inverse par les coefficients de pondération, et alimente un additionneur 62, qui reçoit également la sortie de la mémoire FIFO 63 contenant 2ML-M complexes. Cette FIFO 63 est alimentée par le résultat de l'addition 62. Un module de commande 64 permet de diriger la sortie de la FIFO vers l'extérieur, pour délivrer les M complexes prêts à émettre 65. Il est possible d'utiliser N multiplieurs en parallèle, où N = 1,2,4,...2L. =1,2,4,.. .2L. Ainsi, dans le cas de la figure 6B, on met en oeuvre L(=4) multiplieurs en parallèle 611 à 614. lls sont alimentés alternativement par l'un ou l'autre des coefficients de pondération qui leur sont associés. Ils alimentent chacun un additionneur 621 à 624 qui reçoit également des données des 2L mémoires FIFO 631 à 63 comprenant chacune M complexes. La FIFO 631 délivre les M sorties complexes. Des moyens de sélection 66 permettent de sélectionner la FIFO à prendre en compte à chaque instant. La figure 6C présente le cas de la mise en oeuvre de 2L (=8) multiplieurs. Dans ce cas, la structure ne nécessite plus la présence de moyens de commande. Les 2L mémoires FIFO 63t à 63s sont chacune alimentées par un multiplieur 611 à 618, associé à son propre coefficient de pondération, et associées à un additionneur 621 à 628. n est à noter que la réduction de la place mémoire nécessaire au buffer de sortie est intéressante uniquement si les deux conditions suivantes sont respectées: - l'algorithme utilisé pour réaliser la Wr inverse range ses sorties dans l'ordre "optimal" Coj, CMj,CMj,C1,j,CM+1J,...,CM-1j,C2M-1 et fonctionne au rythme V(2M). En effet, dans le cas de sorties rangées dans l'orde "bit-reverse" ou même dans un ordre "naturel" C0,j, C1j,...C2M-2j, C2M-1j, le gain de mémoire proposé, dû à une écriture et une lecture simultanée dans le buffer de sortie, nécessite alors de réordonner les sorties de la FFT inverse dans l'ordre "optimal", ce qui nécessite une mémorisation en sortie de FFT. L'un dans l'autre, le gain de mémoire sera alors négligeable, voire nul. - les multiplieurs utilisés, pour éviter tout stockage des sorties de la Wr inverse, fonctionnent "rapidement" : le rythme de fonctionnement de N multiplieurs en parallèle doit être égal à Nxo/(2LM), pour N = 1,2,4,...2L. Selon une autre variante de l'invention, il est possible d'optimiser la Wr inverse. Etant donné le caractère complexe "particulier" des entrées (am,n,(-1)m(n+L)im+n) de la FFT inverse à l'émission, il est possible d'utiliser un algorithme de FFT à entrées réelles. Il est connu que déphaser de r les entrées xm d'une FFT de taille 2M, revient à appliquer une permutation circulaire à ses sorties Yk de M/2 vers la gauche. En appliquant ce résultat, il apparaît clairement que l'étape de pré-modulation: (am,nim+n (-1 )m(n+L)) suivie de la m;T inverse complexe peut être réalisée, ainsi que cela est illustré en figure 7, selon les étapes suivantes: - Wr inverse réelle 71 des données (-l)m(n+L) amn. - Permutation circulaire 72 de M/2 des sorties. - Multiplication 73 par in. Un tel algorithme permet de réduire de moitié la place mémoire nécessaire à la FFT, ainsi que le nombre d'opérations. La figure 5 représente ces trois opérations. On peut noter que l'opération de multiplication par (-1)m(n+L) a été omise. En effet, elle peut être évitée. A l'émission, pour générer le signal en bande basse s(t) (équation 1) de façon numérique, il faut multiplier les am,n par (-1)m(n+L) . En réception, I'estimation des données nécessitent à nouveau cette multiplication, comme on le verra par la suite. Etant donné que le retrait de ce terme multiplicatif n'agit en rien sur l'orthogonalité des translatées "temps-fréquence", il est possible de s'affranchir de cette multiplication. Cela revient alors à générer le signal bande basse, centré autour de la fréquence Mvo = fk2, suivant: 6 - Démodulation Le procédé de démodulation doit permettre de récupérer l'information utile émise grâce aux échantillons du signal reçu en réception. Nous supposons ici que le canal "délai-Doppler" (cas le plus général) de fonction de transfert T(f,t) est parfaitement estimé, et qu'il est localement assimilé à un canal multiplicatif complexe Tm,j = Pm je mJ . Etant donné l'orthogonalité des translatées "temps-fréquence" de la fonction prototype, I'information émise à l'instant jT0, sur la porteuse m est ainsi estimée: En pratique, nous travaillons sur les versions échantillonnées à T0 = - du M W signal reçu, la formule de démodulation devient donc: En reprenant les notations données par la formule (2) et en tenant compte du nombre limité de coefficients représentatifs de la fonction prototype (2ML), nous obtenons comme formule de démodulation: PondétT,tion par la I e mrJ (i)ttl+i(l)/tt(j+L) t -L-l fonction prototype -'7i,-771 Pm.j -CRe t A=O r+j +rs7sgA+rq7M Correction de phase et d'amplitude . (14) FKr complexe [p= jM+2qM+k avec 40 IQO La formule (14) suggère cinq étapes pour l'algorithme rapide de démodulation: - Pondération des échantillons reçus par la fonction prototype. - Repliement temporel. - Wr l complexe de taille 2M. - Correction de phase et d'amplitude. - Extraction de la partie réelle. La solution proposée selon l'invention pour minimiser la place mémoire occupée par le filtrage (pondération et repliement temporel) en réception est la suivante: - Insérer les M échantillons reçus dans un buffer d'entrée de taille 2ML complexes. - Multiplier les données de ce buffer par les coefficients représentatifs de la fonction prototype. - Sommer les résultats de ces multiplications (Repliement temporel). - Appliquer une Transformée de Fourier directe aux 2M valeurs complexes ainsi obtenues. - Corriger le résultat de cette FFF en phase et en amplitude. - Extraire la partie réelle. I1 est à noter que cette technique est indépendante de la façon dont le signal a été construit à l'émission. Elle peut s'appliquer à la réception de tout type de signaux multiporteuses OFDM/OQAM. Afin d'illustrer le fonctionnement de cet algorithme, nous décomposons les signaux de taille 2M complexes en 2 sous-signaux de M complexes chacun, de la façon suivante: L'entrée de la FFT sera quant à elle notée k variant de 0 à 2M-1. L'estimation des données émises débutera après un retard de (2L-l)T0. n faut, en effet, que tous les échantillons reçus comprenant la donnée am,j soient stockés dans le buffer d'entrée avant de calculer Les notations (15.1) et (15.2) ci-dessus permettent d'écrire, dans le cas L=4, les entrées de la FFT nécessaires à l'estimation de am,0 ces équations sous la forme (avec k' = RR,0 = r*,g*, + + rk'+2Mgk'+2M+ + rk'+4Mgk'+4M + +rk'+6Mgk'+6M Rk'+M,0= rk'+Mgk'+M + + rk'+3Mgk'+3M + + rk'+5Mgk'+5M + + De même, les entrées de la FFT correspondant à l'estimation des arn. (m=0 à 2M-1) sont ainsi construits: Rk'1 = rk'+Mgk' + + rt'+3Mgt'+2M + + rk'+5Mgk'+4M + + R. k+M.I = = rt'+2Mgt'+M + + rt'+4Mgk'+3M + + rk'+6Mgk'+5M + + La figure 8 illustre le fonctionnement de l'architecture proposée, dans le cas L=4. Le cas général est illutré par la figure 9, et se déduit directement de la figure 8. Sur cette dernière figure, les M échantillons 81 reçus à un instant, sont stockés dans un buffer d'entrée 82. A chaque temps symbole, les données contenues dans ce buffer 82 sont multipliées (83) par les coefficients de pondération 84, représentatifs de la forme d'onde 85 (IOTA dans l'exemple), puis additionnées (86) pour réaliser le repliement. Les données Rkfl correspondantes alimentent la FFT 87, effectuées sur 2M échantillons complexes. On effectue ensuite une correction de phase 88, puis une extraction 89 de la partie réelle. Enfin, on réalise un décalage du contenu du buffer d'entrée 82. La figure 8 présente le contenu de ce buffer pour huit instants successifs, correspondant à la production des sorties âmj à âme+7. Les formules précédentes représentent le cas du démodulateur associé à un modulateur OFDM/OQAM de densité 2. Cependant, cette dernière architecture resterait applicable au cas de la restitution de 2M complexes à partir de M complexes, découlant du cas général de modulation illustré par la formule (11). La formule générale associée serait: avec k' E M-1 .,M - 1) (16) L'algorithme ci-dessus requiert les moyens suivants: - une RAM de taille 2ML complexes (buffer d'entrée). - une ROM de taille (ML+1) réels (coefficients du filtre numérique) . - 1, L ou 2L multiplieurs complexes, selon le degré de parallélisation des opérations de pondération. - un circuit FFT complexe de taille 2M. La figure 10 illustre l'architecture proposée. Le buffer d'entrée, pouvant contenir 2ML complexes, reçoit M échantillons à chaque T0. On effectue la pondération par la fonction prototype par les multiplications 102, puis le repliement, à l'aide de deux additionneurs 1031 et 1032, qui alimentent un buffer 104 de 2M valeurs alimentant la Wr complexe 105. En sortie de la FT;T 105, on effectue une correction de phase et d'amplitude 106, puis l'on sélectionne la partie réelle 107, pour fournir les 2M réels transmis. De même qu'à l'émission, il est possible de paralléliser les opérations de pondérations (multiplications complexes) en utilisant 2L buffers de taille M complexes associés à L, ou 2L multiplieurs, plutôt qu'un seul de taille 2ML. Ces aspects sont illustrés respectivement par les figures 11 et 12. Le fonctionnement se déduit directement de celui décrit en relation avec les figures 6B et 6C pour l'émission. I1 est possible de n'utiliser que (2L-1) RAMs de M complexes (stockage total : 2ML-M complexes au lieu de 2ML) pour stocker les échantillons reçus. n faut pour cela, à l'étape j, lire le kè" (k=0..M-l) échantillon reçu à l'étape (j-(2L1)), et écrire le kème échantillon courant à la même adresse. Pour réaliser cette opération, on pourra, au choix, utiliser des RAMs ou des FIFOs. Les figures 11 et 12 mettent en oeuvre ce procédé, dans le cas d'une structure FIFO. A nouveau, ainsi qu'à l'émission, il est à noter que la réduction de la place mémoire nécessaire au buffer d'entrée est intéressante uniquement si les deux conditions suivantes sont respectées: - l'algorithme utilisé pour réaliser la FFT fonctionne au rythme T avec des entrées se présentant dans l'ordre "optimal" R0j, RMj, R1j, RM+1j, RM+1j,...,R2M-1j. Dans le cas contraire, une mémorisation à l'entrée de la Wr serait nécessaire, et le gain de mémoire sera a lors négligeable, voire nul. - les multiplieurs utilisés, également pour éviter tout stockage des entrées de la Wr, fonctionnent "rapidement" : le rythme de fonctionnement de N multiplieurs en parrallèle doit être égal à N#0/(2LM), pour N = 1,2,4,...2L. REVENDICATIONS 1. Procédé de modulation d'un signal multiporteuse de densité l/(v0.T0)=2, formé de symboles successifs, comprenant chacun M échantillons à émettre, et constitués d'un jeu de 2M fréquences porteuses orthogonales au sens réel, l'écart entre deux fréquences porteuses voisines valant v0 et l'écart entre les instants d'émission de deux symboles consécutifs, ou temps symbole, valant T0, chacune desdites fréquences porteuses étant modulée selon une même fonction prototype de modulation g(t) présentant une longueur de troncature de 2LIo, caractérisé en ce qu'il comprend, à chaque temps symbole, les étapes suivantes: - obtention d'un jeu de 2M coefficients complexes représentatifs de données àémettre; - calcul de 2LM combinaisons linéaires à partir desdits 2M coefficients complexes obtenus, les coefficients de pondération utilisés dans cesdites combinaisons étant représentatifs de ladite fonction prototype g(t), de façon à obtenir 2LM coefficients; - sommation desdits 2LM coefficients pondérés dans des emplacements mémoire prédéterminés respectifs d'une mémoire comprenant 2LM emplacements mémoire représentant 2L groupes de M sommes partielles distinctes, de façon à former progressivement, dans lesdits emplacements mémoire, sur une durée de 2LTO, M échantillons à émettre; - émission desdits échantillons à émettre. 2. Procédé de modulation selon la revendication 1, caractérisé en ce qu'un échantillon à émettre à l'instantjT0 + kTOIM, noté soit s'écrit: où: C0j à C.1j sont les 2M coefficients complexes générés entre les instants jTO et (j+1)TO; et et sskq sont lesdits coefficients de pondération. 3. Procédé de modulation selon la revendication 2, caractérisé en ce que: - = O pour q impair; Bk.q = pour q pair. 4. Procédé de modulation selon la revendication 3, caractérisé en ce qu'il comprend, pour la génération d'un symbole d'indice j formé de M échantillons, les étapes suivantes: - obtention de 2M entrées réelles am représentatives d'un signal source; -pré-modulation de chacune desdites entrées réelles, produisant 2M coefficients complexes; -transformation de Fourier inverse desdits 2M coefficients complexes, produisant 2M coefficients transformés complexes C0j à C,,j - pour chacun des M couples (C, C(k+M)J) dedits coefficients transformés, calcul de 2L couples pondérés, les coefficients de pondération étant représentatives de ladite fonction prototype; - addition du résultat de chacune desdites 2LM valeurs pondérées au contenu de 2LM zones mémoire distinctes, de façon à construire progressivement les échantillons à émettre constituants les symboles j, (j+l), (j+2),... (j+2L-l); - émission de M échantillons correspondant aux M plus anciens contenus desdites zones mémoire puis mise à zéro du contenu de cesdites M zones mémoire. 5. Procédé de modulation selon l'une quelconques des revendications 1 à 4, caractérisé en ce que lesdites étapes sont mises en oeuvre au rythme TM des échantillons. 6. Procédé de modulation selon l'une quelconque des revendications 1 à 5, caractérisé en ce que ladite étape d'émission est suivie d'une étape de mise à jour desdits emplacements mémoire, comprenant: - un décalage physique du contenu de chacun desdits emplacements mémoire, si ces derniers sont des éléments d'un registre à décalage ; ou - une mise à jour des adresses d'écriture et de lecture desdits emplacements mémoire, si ces derniers sont des éléments d'une mémoire RAM. 7. Procédé de modulation selon l'une quelconque des revendications 1 à 6, caractérisé en ce que lesdits coefficients représentatifs de données à émettre sont obtenus par la mise en oeuvre d'une transformation mathématique comprenant les étapes suivantes: - application d'une transformation de Fourier inverse réelle; - permutation circulaire du résultat de cette transformée inverse de M/2 coefficients vers la gauche; - multiplication par i" de chacun desdits coefficients 8. Procédé de modulation selon l'une quelconque des revendications 1 à 7, caractérisé en ce que le signal centré autour de la fréquence Mvo s'écrit: 9. Dispositif de modulation d'un signal multiporteuse de densité l/(v0.T0) = 2, formé de symboles successifs, comprenant chacun M échantillons à émettre, et constitués d'un jeu de 2M fréquences porteuses orthogonales au sens réel, l'écart entre deux fréquences porteuses voisines valant v0 et l'écart entre les instants d'émission de deux symboles consécutifs valant T0, chacune desdites fréquences porteuses étant modulée selon une même fonction prototype de modulation g(t) présentant une longueur de troncature de 2LTO, caractérisé en ce qu'il comprend: - des moyens de mémorisation temporaire de 2L groupes de M sommes partielles; - des moyens de pondération de 2M coefficients complexes représentatifs de données à émettre par des coefficients de pondération représentatifs de ladite fonction prototype g(t); - des moyens de sommation des coefficients pondérés dans des emplacements mémoire prédéterminés respectifs desdits moyens de mémorisation temporaire, de façon à former progressivement, sur une durée de 2LTO, lesdits échantillons à émettre. 10. Dispositif de modulation selon la revendication 9, caractérisé en ce qu'il comprend: - des moyens de transformation mathématique délivrant lesdits coefficients représentatifs de données à émettre au rythme #0/2M et dans l'ordre suivant (C0j, CM+1j),..(CM-1j, C2M.1j); - 2LM-M emplacements mémoire de type RAM à écriture et lecture simultanées; - N multiplieurs complexes fonctionnant au rythme NtJ2LM, N valant 1, 2, 4,...ou 2L. 11. Procédé de démodulation d'un signal reçu, correspondant à un signal émis multiporteuse de densité 1/(VO.TO) = 2, formé de symboles successifs, représentés chacun par M échantillons à émettre, et constitués d'un jeu de 2M fréquences porteuses orthogonales au sens réel, l'écart entre deux fréquences porteuses voisines valant v0 et l'écart entre les instants d'émission de deux symboles consécutifs, ou temps symbole, valant T0, chacune desdites fréquences porteuses étant modulée selon une même fonction prototype de modulation g(t) présentant une longueur de troncature de 2LTO, caractérisé en ce que l'on reconstruit une estimation des 2M données réelles émises à un temps symbole donné, à l'aide des étapes suivantes: - échantillonnage dudit signal reçu à la fréquence d'échantillonnage zdM, délivrant M échantillons complexes reçus; - mémorisation de chacun desdits M échantillons complexes reçus dans un emplacement prédéterminé d'une mémoire d'entrée comprenant 2ML emplacements complexes, dans laquelle ont été préalablement mémorisés (2L-1)M échantillons reçus pendant les 2L-1 temps symbole précédents; - multiplication des 2ML valeurs contenues dans ladite mémoire d'entrée par des coefficients représentatifs de ladite fonction prototype; - repliement temporel, par sommation de 2M séries de L résultats de multiplication, de façon à obtenir 2M valeurs complexes; - traitement desdites 2M valeurs complexes pour former lesdites estimations des 2M données réelles émises. 12. Procédé de démodulation selon la revendication 11, caractérisé en ce que les 2M valeurs complexes issues de l'étape de repliement temporel entre les instants (j+2L-1)To et (j+2L)To s'écrivent: où: rk. &alpha;'k,q' et ss'k,q' sont lesdits coefficients de pondération. 13. Procédé de démodulation selon l'une quelconque des revendications 11 et 12, caractérisé en ce que: - CC'kq = 0 pour q' impair; - ss'k,q' = 0 pour q' pair. 14. Procédé selon l'une quelconque des revendications 11 à 13, caractérisé en ce que ladite étape de traitement comprend les étapes suivantes - application d'une transformation mathématique inverse de celle effectuée lors de la modulation sur lesdites 2M valeurs complexes, délivrant 2M valeurs transformées; - correction de distorsions de phase et/ou d'amplitude dues au canal de transmission; - extraction de la partie réelle desdites valeurs complexes transformées. 15. Procédé de démodulation selon l'une quelconque des revendications 11 à 14, caractérisé en ce que lesdites étapes sont mise en oeuvre au rythme T 16. Dispositif de démodulation d'un signal reçu, correspondant à un signal émis multiporteuse de densité l/(vo.o) = 2, formé de symboles successifs, représentés chacun par M échantillons à émettre, et constitués d'un jeu de 2M fréquences porteuses orthogonales au sens réel, l'écart entre deux fréquences porteuses voisines valant v0 et l'écart entre les instants d'émission de deux symboles consécutifs, ou temps symbole, valant #0, chacune desdites fréquences porteuses étant modulée selon une même fonction prototype de modulation présentant une longueur de troncature de 2LTO, caractérisé en ce qu'il comprend: - des moyens d'échantillonnage dudit signal reçu; - des moyens de mémorisation temporaire des échantillons complexes échantillonnés, comprenant 2ML emplacements complexes; - des moyens de multiplication desdits échantillons mémorisés par des coefficients de pondération représentatifs de ladite fonction prototype; - des moyens de repliement temporel, assurant la sommation de L résultats de pondération, de façon à obtenir 2M valeurs complexes; - des moyens de traitement desdites valeurs complexes, délivrant une estimation de 2M données réelles émises à chaque temps symbole. 17. Dispositif de démodulation selon la revendication 16, caractérisé en ce que lesdits moyens de traitement comprennent: - des moyens de transformation mathématique inverse de celle effectuée lors de la modulation sur lesdites 2M valeurs complexes; - des moyens de correction de distorsions de phase et/ou d'amplitude dues au canal de transmission; - des moyens d'extraction de la partie réelle des valeurs transformées. 18. Dispositif de démodulation selon l'une quelconque des revendications 16 et 17, caractérisé en ce qu'il comprend: - des moyens de mémorisation comprenant 2ML-M emplacements mémoire complexes de type RAM à écriture et lecture simultanée; - N multiplieurs complexes fonctionnant au rythme Nzd2LM, où N vaut 1, 2, 4... ou 2L ; - des moyens de transformation mathématique fonctionnant au rythme T RMj), (Rlj, RM+1j),... (RM-1,1', R2M-1j). 19. Procédé de filtrage délivrant des séries de M valeurs complexes de sortie délivrées à intervalles réguliers, à partir de 2L séries de 2M valeurs complexes d'entrée, lesdites M valeurs complexes correspondant à une somme pondérée de 2L desdites valeurs complexes d'entrée à traiter, caractérisé en ce qu'il comprend les étapes suivantes, pour chaque série de valeurs complexes d'entrée: - calcul de 2LM combinaisons linéaires à partir desdits 2M coefficients complexes obtenus, les coefficients de pondération étant issus de 2L filtres complexes ou réels de taille M, de façon à obtenir 2LM coefficients; - sommation de chacune des valeurs pondérées dans un emplacement mémoire prédéterminé, parmi un ensemble de 2ML emplacement mémoire contenant chacun une somme partielle, de façon à former progressivement lesdites valeurs de sortie, dans lesdits emplacements mémoire, sur une durée correspondant à la réception de 2L séries de valeurs complexes d'entrée.