FR 2733869 A1 19961108 FR 9505455 A 19950502 Signal multiporteuse numérique, procédé de construction d'un tel signal et procédés d'émission et de réception correspondants. 1. Domaine de l'invention Domaine général 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 des signaux produits à l'aide de nouvelles modulations, ainsi que les techniques de modulation et de démodulation correspondantes. Depuis de nombreuses années, on cherche à construire des modulations adaptées à des canaux fortement non-stationnaires, tels que les canaux de transmission vers des mobiles. Dans de tels canaux, le signal émis est affecté d'évanouissements et de trajets multiples. Les travaux menés par le CCE1T dans le cadre du projet européen EUREKA 147 (DAB: Digital Audio Broadcasting, ou Diffusion Audionumérique) ont montré l'intérêt, pour ce type de canaux, des modulations multiporteuses, et en particulier de l'OFDM (Orthogonal Frequency Division Multiplexing). L'OFDM a été retenu dans le cadre de ce projet européen comme base de la norme DAB. Cette technique est également envisagée comme modulation pour la diffusion de programmes de télévision. Cependant, on constate un certain nombre de limitations (précisées par la suite) lorsqu'on aborde le problème de modulations codées à haute efficacité spectrale, telles que celles requises pour les applications de TV numérique. JIL Applications possibles 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 nonstationnaire. 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 intrinsè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èmes de communication numérique vers des mobiles à haut débit, dans le cadre par exemple de l'UMTS (projet RACE). Elle peut également être envisagée pour des réseaux locaux radio à haut débit (type HIPERLAN). Une troisième catégorie d'applications est celle des transmissions sous-marines. Le canal de transmission en acoustique sous-marine est fortement perturbé du fait de la faible vitesse de transmission des ondes acoustiques dans l'eau. Ceci conduit à un étalement important des trajets multiples et du spectre Doppler. Les techniques de modulation multiporteuses sont donc bien adaptées à ce domaine, et tout particulièrement les techniques objet de la présente invention. 2. Etatdelatechnique 2.1. Remarques théoriques sur la représentation des signaux Avant de présenter les signaux selon l'invention, on décrit ci-dessous les signaux connus. Cette description repose sur une approche générale des signaux multiporteuses définie par les inventeurs, et nouvelle en elle-même. Cette généralisation n'a en effet aucun équivalent dans l'état de la technique, et n'est nullement évidente pour l'homme du métier. Elle doit donc être considérée comme une partie de l'invention, et non comme appartenant à l'état de la technique. On s'intéresse à des signaux réels (une grandeur électrique par exemple), à énergie finie, et fonction du temps. Les signaux peuvent donc être représentés par des fonctions réelles de L2 (R). En outre, ces signaux sont à bande limitée w et leur spectre est contenu dans [fc-2stc+2]s fc étant la "fréquence porteuse" du signal. On peut donc de façon équivalente représenter un signal réel a(t) par le son enveloppe complexe s(t) avec: s(t) = e''9=' F,[a](t) (1) OÙ FA désigne le filtre analytique. Le signal s(t) appartient à un sous espace vectoriel (caractérisé par la limitation de bande à + ) de l'espace des fonctions complexes d'une variable réelle et de carré bande sommable L2 (R). On peut définir cet espace vectoriel de deux façons différentes, selon qu'on le construit sur le corps des complexes ou le corps des réels. A chacun de ces espaces, on peut associer un produit scalaire à valeur dans C ou dans R et construire un espace de Hilbert. On appellera H l'espace de Hilbert construit sur le corps des complexes et HR l'espace de Hilbert construit sur le corps des réels. Les produits scalaires correspondants s'écrivent: dans le cas de H (2) dans le cas de HR (3) Les normes associées sont évidemment identiques dans les deux cas: 2.2. Principes généraux de l'OFDM Les principes généraux de l'OFDM sont par exemple présentées dans le brevet français FR-86 09622 déposé le 2 juillet 1986. L'idée de base de cette technique est de transmettre des symboles codés comme des coefficients de formes d'ondes élémentaires confinées autant que possible dans le plan temps-fréquence, et pour lesquels le canal de transmission peut être considéré comme localement stationnaire. Le canal apparaît alors comme un simple canal multiplicatif caractérisé par la distribution du module des coefficients, qui suit une loi de Rice ou de Rayleigh. On assure ensuite la protection contre les évanouissements à l'aide d'un code utilisable en décision pondérée, en association à un entrelacement en temps et en fréquence qui garantit que les symboles intervenant dans la maille minimale du code soient dans toute la mesure du possible affectés par des évanouissements indépendants. Cette technique de codage avec entrelacement dans le plan temps-fréquence est connue sous le nom de COFDM. Elle est par exemple décrite dans le document [23] (voir annexe 1 (pour simplifier la lecture, la plupart des références de l'état de la technique sont listées dans cette annexe 1. Celle-ci, ainsi que les annexes 2 et 3 doivent bien sûr être considérées comme des éléments à part entière de la présente description)). n existe essentiellement deux types de modulations OFDM connues. Les appellations utilisées dans la littérature étant souvent ambiguës, nous introduirons ici des appellations nouvelles plus précises tout en rappelant la correspondance avec la littérature existante. Nous utiliserons l'appellation générique OFDM, suivi d'un suffixe précisant le type de modulation à l'intérieur de cette famille. 2.3. OFDM/OAM mol2. Principes théoriques Une première catégorie de modulations est constituée d'un multiplex de porteuses QAM (Quadrature Amplitude Modulation), ou éventuellement en QPSK (Quadrature Phase Shift Keying) dans le cas particulier de données binaires. Nous désignerons par la suite ce système sous le nom OFDM/QAM. Les porteuses sont toutes synchronisées, et les fréquences porteuses sont espacées de l'inverse du temps symbole. Bien que les spectres de ces porteuses se recouvrent, la synchronisation du système permet de garantir l'orthogonalité entre les symboles émis par différentes porteuses. Les références [1] à [7] donnent un bon aperçu de la littérature disponible sur ce sujet. Pour plus de simplicité dans l'écriture, et selon l'approche nouvelle de l'invention, on représentera les signaux par leur enveloppe complexe décrite ci-dessus. Dans ces conditions, l'équation générale d'un signal OFDM/QAM s'écrit: Les coefficients am,n prennent des valeurs complexes représentant les données émises. Les fonctions xm.n( t) sont des translatées dans l'espace temps-fréquence d'une même fonction prototype x(t): Si Iti (6) ailleurs xm,n(t)=ei(2#mv0@+#)x(t-n#0) avec v0#0 = 1 (7) # étant une phase quelconque, que l'on peut arbitrairement fixer à 0.La fonction x(t) est centrée, c'est-à-dire que ses moments d'ordre 1 sont nuls, soit: #t|x(t)|2dt=#f|X(f)|2df = 0, (8) X(f) désignant la transformée de Fourier de x(t). Dans ces conditions, on observe que: ItIXm,n(t) 2dt = (9) #f|Xm,n(f)| df = mv0 Les barycentres des fonctions de base forment donc un réseau du plan temps fréquence engendré par les vecteurs (T0, 0) et (0, v0), ainsi que cela est illustré en figure 1. Ce réseau est de densité unité, c'est à dire que v0T0 = 1. On pourra se reporter à l'article [9] pour une discussion plus détaillée sur ce sujet. La fonction prototype x(t) a ceci de particulier que les fonctions {xn,n} sont orthogonales entre elles, et plus précisément constituent une base hilbertienne de L2 (R), soit: 1 si(m,n) = (m',n') = # (10) 0 sinon Projeter un signal sous cette base équivaut simplement à découper le signal en séquences de durée 0 et à représenter chacune de ces séquences par le développement en série de Fourier correspondant. Ce type de décomposition constitue un premier pas vers une localisation à la fois en temps et en fréquence, par opposition à l'analyse de Fourier classique, qui assure une localisation fréquentielle parfaite avec une perte totale de l'information temporelle. Malheureusement, si la localisation temporelle est excellente, la localisation fréquentielle est beaucoup moins bonne, du fait de la décroissante lente de X(f). Le théorème de Balian-Low-Coifman-Semmes (voir [9], p 976) montre d'ailleurs que si l'on appelle X la transformée de Fourier de x, tx(t) et fX(f) ne peuvent être simultanément de carré sommable. 2.3.2. L'OFDM/OAM avec intervalle de garde D'une manière générale, on peut caractériser la tolérance d'une modulation OFDM vis-à-vis des trajets multiples et de l'étalement Doppler par un paramètre mesurant de manière globale la variation du niveau d'interférence entre symboles (lES) en fonction d'un décalage en temps ou en fréquence. La justification de ce concept est donnée en annexe 2. Ce paramètre de tolérance est appelé # et est défini par la relation: 5 1/4XtAf (11) avec: En vertu de l'inégalité de Heisenberg, 5 ne peut pas dépasser l'unité. Compte tenu du théorème de Balian-Low-Coifman-Semmes cité précédemment, le paramètre 5 vaut 0 pour l'OFDM/QAM. Il s'agit là d'un défaut important de la modulation OFDM/QAM telle que décrite ci-dessus. Ceci se caractérise en pratique par une forte sensibilité aux erreurs temporelles, et par conséquent aux trajets multiples. Ce défaut peut être contourné par l'utilisation d'un intervalle de garde décrit par exemple dans [5]. Il s'agit là d'un artifice consistant à prolonger la fenêtre rectangulaire de la fonction prototype. La densité du réseau des symboles de base est alors strictement inférieure à l'unité. Cette technique est possible du fait que l'on retrouve à l'intérieur d'un symbole prolongé par un intervalle de garde une infinité de versions translatées du symbole initial. Bien sûr, ceci ne fonctionne que parce que la fonction prototype est une fenêtre rectangulaire. En ce sens, l'OFDM(QAM avec intervalle de garde constitue un point singulier unique. La modulation OFDM/QAM avec intervalle de garde est à la base du système DAB. Cet intervalle de garde permet de limiter l'interférence entre symboles, au prix d'une perte de performance, puisqu'une partie de l'information émise n'est pas réellement utilisée par le récepteur, mais sert seulement à absorber les trajets multiples. Ainsi, dans le cas du système DAB, où l'intervalle de garde représente 25% du symbole utile, la perte est de 1 dB. En outre, il existe une perte supplémentaire, due au fait que pour obtenir une efficacité spectrale globale donnée, il faut compenser la perte due à l'intervalle de garde par un meilleure efficacité du code utilisé. Cette perte est marginale dans le cas du système DAB, parce que l'efficacité spectrale est faible. Par contre, si l'on vise une efficacité spectrale globale de 4 bits/Hz, il faut utiliser un code à 5 bits/Hz, soit selon le théorème de Shannon une perte de l'ordre de 3 dB. La perte globale est donc dans ce cas d'environ 4 dB. 2.3.3. Autres systèmes OFDM/OAM On peut imaginer d'autres systèmes de type OFDM/QAM. Malheureusement, aucune modulation QAM filtrée (c'est à dire utilisant une mise en forme conventionnelle de type demi-Nyquist (ou, plus exactement, racine carrée de Nyquist )), ne vérifie les contraintes d'orthogonalité requises. Les fonctions prototypes connues vérifiant les critères d'orthogonalité requis sont: - la fenêtre rectangulaire; - le sinus cardinal. Ces deux exemples sont triviaux, et apparaissent duaux l'un de l'autre par transformée de Fourier. Le cas de la fenêtre rectangulaire correspond à l'OFDM/QAM sans intervalle de garde. Le cas du sinus cardinal correspond à un multiplex fréquentiel classique (c'est-à-dire dont les porteuses ont des spectres disjoints) avec un roll-off de 0%, ce qui constitue un cas asymptotique difficilement réalisable en pratique. Dans chacun de ces cas, on observe que la fonction prototype est parfaitement bornée, soit en temps, soit en fréquence, mais possède une décroissance médiocre (en 1/t ou 1/f) dans le domaine dual. Le théorème de Balian-Low-Coifman-Semmes laisse d'ailleurs peu d'espoir que des solutions satisfaisantes puissent exister. Comme indiqué précédemment, ce théorème démontre que tx(t) et fX(f) ne peuvent être simultanément de carré sommable. On ne peut donc pas espérer trouver une fonction x(t) telle que x(t) et X(f) décroissent simultanément avec un exposant inférieur à -3/2. Ceci n'exclut d'ailleurs pas que puissent exister des fonctions satisfaisantes aux yeux d'un ingénieur. Néanmoins, un article récent [10] traitant de ce sujet exhibe un autre exemple de fonction prototype ayant les propriétés requises. L'allure de la fonction prototype proposée dans cet article est très éloignée de ce que l'on peut souhaiter en terme de concentration temporelle. II est donc probable qu'il n'existe pas de solution satisfaisante de type OFDM/QAM. En conclusion, 1'OFDM/QAM, correspondant à l'utilisation d'un réseau de densité 1 et à des coefficients am.n complexes ne peut être mis en pratique que dans le cas d'une fenêtre temporelle rectangulaire et de l'utilisation d'un intervalle de garde. L'homme du métier cherchant d'autres modulations est donc conduit à se tourner vers les techniques décrites ci-dessous sous le nom d'OFDM/OQAM. 2.4. OFDM/OOAM Une deuxième catégorie de modulations utilise en effet un multiplex de porteuses OQAM (Offset Quadrature Amplitude Modulation). Nous désignerons par la suite ce système sous le nom OFDM/OQAM. Les porteuses sont toutes synchronisées, et les fréquences porteuses sont espacées de la moitié de l'inverse du temps symbole. Bien que les spectres de ces porteuses se recouvrent, la synchronisation du système et le choix des phases des porteuses permet de garantir l'orthogonalité entre les symboles émis par différentes porteuses. Les références [11-18] donnent un bon aperçu de la littérature disponible sur ce sujet. Pour plus de simplicité dans l'écriture, on représentera les signaux sous leur forme analytique. Dans ces conditions, l'équation générale d'un signal OFDM/OQAM s'écrit: Les coefficients am,n prennent des valeurs réelles représentant les données émises. Les fonctions xm,n( t) sont des translatées dans l'espace temps-fréquence d'une même fonction prototype x(t): si m + n est pair (15) si m + n est impair avec v0#0=1/2, # étant une phase quelconque que l'on peut arbitrairement fixée égale à 0. Les barycentres des fonctions de base forment donc un réseau du plan tempsfréquence engendré par les vecteurs (X0, 0) et (0, v0), tel qu'illustré en figure 2. Ce réseau est de densité 2. Les fonctions xn,.n(t) sont orthogonales au sens du produit scalaire dans R. Dans les approches connues, la fonction prototype est bornée en fréquence, de telle sorte que le spectre de chaque porteuse ne recouvre que celui des porteuses adjacentes.En pratique, les fonctions prototypes considérées sont des fonctions paires (réelles ou éventuellement complexes) vérifiant la relation suivante: X(f) = 0 si|f|#v0 (16) |X(f)|2+|X(f - v0)|2 = 1/v0 si 0#f#v0 Un choix possible pour x(t) est la réponse impulsionnelle d'un filtre demi-Nyquist 100% roll-off, soit Si If ailleurs Lorsqu'on observe x(t) et sa transformée de Fourier, on note que X(f) est à support borné et que x(t) décroît en t-2, c'est à dire un résultat notablement meilleur que la limite théorique découlant du théorème de Balian-Low-Coifman-Semmes.Les formes d'onde élémentaires sont mieux localisées dans le plan temps fréquence que dans le cas de I'OFDM/QAM, ce qui confere à cette modulation un meilleur comportement en présence de trajets multiples et de Doppler. Comme précédemment, on peut définir le paramètre 5 mesurant la tolérance de la modulation au délai et au Doppler. Ce paramètre 5 vaut 0.865. 3. Inconvénients des systèmes connus Ces systèmes connus présentent de nombreux inconvénients et limites, notamment dans les canaux très perturbés, et lorsque une haute efficacité est requise. 3. 1. OFDM/QAM Le problème principal du système OFDM/QAM est qu'il nécessite impérativement l'utilisation d'un intervalle de garde. Comme indiqué précédemment, ceci engendre une perte d'efficacité notable lorsque l'on vise de hautes efficacités spectrales. De plus, les symboles émis sont mal concentrés dans le domaine fréquentiel, ce qui limite également les performances dans des canaux fortement non-stationnaires. En particulier, cet étalement rend difficile l'utilisation d'égaliseurs. 3.2. OFDM/OOAM A l'inverse, les performances fréquentielles de l'OFDMtOQAM sont plutôt satisfaisantes et le problème de la perte liée à l'intervalle de garde ne se pose pas. En revanche, la réponse impulsionnelle de la fonction prototype a une décroissance temporelle relativement lente, soit en l/x2. Ceci implique deux types de difficultés. Tout d'abord, la forme d'onde peut difficilement être tronquée sur un intervalle de temps court, ce qui implique un traitement complexe au niveau du récepteur. En outre, ceci complique également d'éventuels systèmes d'égalisation. En d'autres termes, l'efficacité des techniques OFDM/OQAM est supérieure à celle de l'OFDM/QAM, mais ces techniques s'avèrent plus complexes à mettre en oeuvre, et donc coûteuses, en particulier dans les récepteurs. 4. Présentation de l'invention 4.1 Obiectifs de l'invention L'invention a notamment pour objectif de pallier ces différents inconvénients et limitations de l'état de la technique. Ainsi, un objectif de l'invention est de fournir un signal numérique destiné à être transmis ou diffusé vers des récepteurs, qui permette d'obtenir de meilleures performances dans des canaux non stationnaires, et tout particulièrement dans des canaux fortement non stationnaires. L'invention a également pour objectif de fournir un tel signal, permettant d'obtenir une haute efficacité spectrale. Un autre objectif de l'invention est de fournir un tel signal, qui évite les inconvénients de l'OFDM/QAM liés à l'intervalle de garde, tout en conservant une réponse temporelle de la fonction prototype aussi concentrée que possible, de façon notamment à simplifier le traitement au niveau du récepteur. L'invention a également pour objectif de fournir un tel signal, permettant la réalisation de récepteurs à complexité et coût limités, en particulier en ce qui concerne la démodulation et l'égalisation. Un objectif complémentaire de l'invention est de fournir des émetteurs, des récepteurs, des procédés de transmission ou de diffusion, des procédés de réception et des procédés de construction, c'est-à-dire de définition, d'une modulation correspondants à un tel signal. 4.2 Caractéristiques principales de l'invention Ces objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints selon l'invention par un signal multiporteuse destiné à être transmis vers des récepteurs numériques, notamment dans un canal de transmission non-stationnaire, correspondant au multiplexage en fréquence de plusieurs porteuses élémentaires correspondants chacunes à une série de symboles, deux symboles consécutifs étant séparés d'un temps symbole T0, signal dans lequel, d'une part, l'espacement v0 entre deux porteuses voisines est égal à la moitié de l'inverse du temps symbole 0, et dans lequel, d'autre part, chaque porteuse subit un filtrage de mise en forme de son spectre présentant une largeur de bande strictement supérieure à deux fois ledit espacement entre porteuses v0.Ce spectre est choisi de façon que chaque élément de symbole soit concentré autant que possible à la fois dans le domaine temporel et dans le domaine fréquentiel. Notamment, un tel signal peut répondre à l'équation suivante: où : am,n est un coefficent réel représentatif du signal source, choisi dans un alphabet de modulation prédéterminé; m est un entier représentant la dimension fréquentielle; n est un entier représentant la dimension temporelle; t représente le temps; X,,(t) est une fonction de base, translatée dans l'espace temps-fréquence d'une même fonction prototype x(t) paire prenant des valeurs réelles ou complexes, soit où (p est un paramètre de phase arbitraire, la transformée de Fourier X(f) de la fonction x(t) ayant un support s'étendant au delà de l'intervalle [-v0,v0], et où lesdites fonctions de base {xn} sont orthogonales entre elles, la partie réelle du produit scalaire de deux fonctions de base différentes étant nulle. Le symbole "+" indique que xn,n(t) peut prendre indifféremment un signe négatif ou positif. I1 ne signifie pas, bien sûr, que x,l, n(t) prend les deux valeurs. Ainsi, l'invention repose sur un système de modulation qui utilise des fonctions prototypes aussi concentrées que possible dans le plan temps-fréquence. L'intérêt de cette approche est de disposer d'une modulation produisant un signal évitant les inconvénients de l'OFDM/QAM liés à l'intervalle de garde, tout en conservant une réponse temporelle de la fonction prototype aussi concentrée que possible, de façon à simplifier le traitement au niveau du récepteur. En d'autres termes, l'invention a pour objet de nouveaux systèmes de modulation construits comme l'OFDM/OQAM sur un réseau orthogonal de densité 2, sans pour autant que la fonction prototype soit à support borné en fréquence. Parmi les modulations proposées, on trouve soit des modulations utilisant des fonctions prototypes à support borné en temps, soit des fonctions prototypes qui ne sont bornés ni en temps ni en fréquence, mais qui présentent par compte des propriétés de décroissance rapide à la fois en temps et en fréquence, et une concentration quasi-optimale dans le plan tempsfréquence. De tels signaux ne sont nullement évidents pour l'homme du métier, au vu de l'état de la technique. Comme indiqué précédemment, il existe fondamentalement deux modes de construction de modulations de type OFDM. Le premier mode de construction connu utilise un réseau de densité 1. Cette première solution utilise une base de décomposition des signaux ou tout signal est découpé en intervalles, chaque intervalle étant ensuite décomposé sous forme de série de Fourier. C'est la solution OFDM/QAM. La littérature donne peu d'exemples de solutions alternatives construites sur le même réseau, et les résultats obtenus sont d'un faible intérêt pratique [10]. En outre, la technique OFDM/QAM est la seule qui puisse bénéficier de la méthode de l'intervalle de garde. La solution OFDM/QAM est donc un point singulier qui ne permet pas d'extensions. Le deuxième mode de construction connue (OFDM/OQAM) utilise un réseau de densité 2. L'orthogonalité entre symboles centrés sur une même fréquence ou sur des fréquences adjacentes est garantie par une mise en forme des signaux de type demi Nyquist et par le choix adéquat de la phase des signaux. Enfin, l'orthogonalité au delà des fréquences adjacentes est garantie par le fait que les supports fréquentiels sont disjoints. Par conséquent, la construction de nouvelles modulations ne vérifiant pas cette propriété n'est pas évidente. Toutes les variantes de l'invention décrites ci-dessous présentent l'avantage d'utiliser une fonction prototype, soit limitée dans le domaine temporel, soit à décroissance rapide, de telle sorte que la fonction puisse être aisément tronquée. Selon une première variante, ladite fonction prototype x(t) est une fonction paire, nulle en dehors de l'intervalle [-#0,#0], et vérifiant la relation: x(t)=0 si|t|##0 Ix(t)12 +Ix(t-r0)12 = 1/TO si O t De façon avantageuse, ladite fonction prototype x(t) est définie par: si |t|##0 ailleurs Dans ce premier cas (appelé par la suite OFDM/MSK), les performances en terme de résistance au Doppler et aux trajets multiples sont équivalentes à 1'OFDM/OQAM, et la réalisation du récepteur est simplifiée. Selon une deuxième variante de l'invention, ladite fonction prototype x(t) est caractérisée par l'équation: la fonction y(t) étant définie par sa transformée de Fourier Y(f): où G(f) est une fonction gaussienne normalisée du type: G(f) = (2a)"4e2 a étant un paramètre réel strictement positif. k varie de -oo à +oo. De façon avantageuse, le paramètre a est égal à l'unité. La modulation correspondante est appelée par la suite OFDM/IOTA. Dans ce cas, la fonction prototype correspondante, notée 3, est identique à sa transformée de Fourier. La réalisation du récepteur est plus simple que dans le cas de l'OFDM/OQAM, bien que légèrement plus complexe que dans le cas précédent, mais les performances sont sensiblement supérieures. L'invention concerne également un procédé de transmission d'un signal numérique, notamment dans un canal de transmission non-stationnaire, comprenant les étapes suivantes: - codage canal d'un signal numérique à transmettre, délivrant des coefficients numériques réels a",n choisis dans un alphabet prédéterminé; - construction d'un signal s(t) répondant à l'équation définie ci-dessus; - émission d'un signal ayant pour enveloppe complexe ledit signal s(t) vers au moins un récepteur. De façon avantageuse, un tel procédé comprend de plus une étape d'entrelacement en fréquence et/ou en temps, appliquée aux éléments binaires formant ledit signal numérique à transmettre ou aux coefficients numériques a,,. Cela permet d'assurer des performances optimales dans des canaux nonstationnaires. L'invention concerne également les émetteurs d'un tel signal. L'invention concerne encore un procédé de réception d'un signal tel que décrit cidessus, qui comprend les étapes suivantes: - réception d'un signal ayant pour enveloppe complexe un signal r(t) correspondant au signal s(t) à l'émission; - estimation de la réponse du canal de transmission, comprenant une estimation de la réponse en phase g,n et de la réponse en amplitude pmn; - démodulation dudit signal r(t), comprenant les étapes suivantes:: - multiplication dudit signal r(t) par la fonction prototype x(t) - repliement de la forme d'onde filtrée modulo 2#0; - application d'une transformée de Fourier (FFr) , - correction de la phase #m,n induite par le canal de transmission ; - correction de la phase correspondant au terme im+n; - sélection de la partie réelle du coeeficient obtenu am,n correspondant au coefficient am,n émis pondéré par la réponse en amplitude #m,n du canal de transmission. De façon préférentielle, ce procédé de réception comprend une étape de désentrelacement en fréquence et/ou en temps desdits coefficients numériques réels #m,n et, éventuellement, des valeurs correspondantes #m,n de la réponse de l'amplitude du canal, ledit désentrelacement étant inverse de un entrelacement mis en oeuvre à l'émission, et/ou une étape de décodage en décision pondérée adapté au codage canal mis en oeuvre à l'émission. L'invention concerne également les récepteurs correspondants. Enfin, l'invention concerne également un procédé préférentiel de construction d'une fonction prototype x(t) pour un signal tel que décrit ci-dessus, comprenant les étapes suivantes: - sélection d'une fonction gaussienne G(f) normalisée du type: G(f) = (2a)H4 e~t - détermination de ladite fonction prototype x(t) telle que: la fonction y(t) étant définie par sa transformée de Fourier Y(f): Ce procédé permet notamment de définir la fonction prototype S, décrite dessus. 5. Description de modes de réalisation particuliers de l'invention 5.1. Liste des figures - la figure 1 illustre un réseau de densité 1, correspondant à celui mis en oeuvre dans le cas de la modulation connue OFDM/QAM; - la figure 2 illustre un réseau de densité 2, correspondant à celui mis en oeuvre dans le cas de la modulation connue OFDM/OQAM et dans le cas de l'invention; - les figures 3A à 3D, 4C à 4D, 5A à 5D, 6A à 6D et 7A à 7D illustrent respectivement les modulations connues OFDM/QAM (3), OFDM/QAM avec intervalle de garde (4), OFDM/OQAM (5) et les modulations de l'invention OFDM/MSK (6) et OFDM/IOTA (7), selon les aspects suivants . A: la fonction prototype x(t); . B : la transformée de Fourier en linéaire de la fonction prototype; .C : le module de la fonction d'ambiguïté en linéaire (telle que définie dans l'annexe 2); . D : la fonction d'intersymbole (telle que définie dans l'annexe 2); - la figure 7E montre en échelle logarithmique la décroissance du signal OFDM/IOTA; - la figure 8 la fonction d'ambiguïté d'une fonction gaussienne; - la figure 9 est un schéma synoptique d'un émetteur (et du procédé d'émission correspondant) utilisable selon l'invention; - la figure 10 est un schéma synoptique d'un récepteur (et du procédé de réception correspondant) utilisable selon l'invention - la figure 11 illustre plus précisément le procédé de démodulation mis en oeuvre dans le récepteur de la figure 10. L Principes théoriques des signaux selon l'invention Toutes les fonctions de base de l'OFDM/OQAM définies en (15) peuvent se réécrire sous la forme: Les barycentres des fonctions de base forment donc un réseau du plan temps fréquence engendré par les vecteurs (r0, 0) et (0, v0) (voir figure 2). Ce réseau est de densité 2, c'est à dire que v0#0 = 1/2. Comme indiqué en [16], ces fonctions constituent une base hilbertienne de HR. Afin de simplifier l'écriture, nous omettrons par la suite les inversions de signe. D'une manière générale, on cherche les conditions sur x(t) pour que la famille {xm.n} constitue une base hilbertienne de HR. On imposera que x(t) soit une fonction paire. Le produit scalaire de xm,n et de xi ni peut s'écrire: soit, en posant t'= t-(n+(n+n')#0/2 et #'0=(n-n')#0: = #e[i(m-m')+(n-n')i(m-m')(n+n') J(fonction paire + i fonction impaire)] L'orthogonalité est donc obtenue si le coefficient de l'intégrale est un nombre imaginaire pur. L'analyse de ce coefficient montre qu'il suffit pour cela que m-m' ou que n-n' soit impair. Le réseau peut donc se décomposer en quatre sous-réseaux, tel que cela apparaît sur la figure 2 ({m pair, n pair},{m pair, n impair impair n pair},{m impair, n impair}) qui sont orthogonaux entre eux (toute fonction de l'un des sous-réseaux est orthogonale à toute fonction d'un autre sous-réseau). Une condition suffisante pour que {xm,n} constitue une base hilbertienne est donc que: R = 0 # m-m' pair, # n-n'pair, (m,n)#(m',n') (21) I1 suffit donc de trouver une fonction x(t) paire telle que les fonctions du type: soient orthogonales entre elles au sens du produit scalaire dans R. D'ailleurs, si tel est le cas, ces fonctions sont aussi orthogonales au sens du produit scalaire dans C, pour des raisons de symétrie analogues à celles évoquées ci-dessus. Une autre façon d'exprimer cette condition est d'utiliser la fonction d'ambiguïté de x [19]: n suffit alors de trouver une fonction x(t) paire telle que: A (2nlO, 2mv0) = O ,#(m,n)# (0,0) (24) Si l'on compare le problème ainsi posé à celui de trouver une base hilbertienne au sens du produit scalaire dans C, les contraintes d'orthogonalité sont nettement moins fortes, puisque le réseau concerné est deux fois moins dense. En effet, les fonctions de base sont centrées sur les points du réseau {2mv0,2n#0}, c'est à dire un réseau de densité 1/2. On voit donc apparaître ici de façon intuitivement évidente les raisons de l'inapplicabilité du théorème de Balian-Low-Coifman-Semmes. Dans le cas de l'OFDM/OQAM, l'orthogonalité des fonctions X2m.2n(t) entre elles est obtenue par deux contraintes de natures différentes. En effet, si m # m', est nul parce que ces fonctions ont des spectres disjoints. Par ailleurs, (x2m 2n | x2m 2" ) est nul parce que X(f) a une mise en forme de type demi Nyquist. Comme le montre l'abondante littérature déjà citée, I'homme du métier considère qu'il est impératif de vérifier ces deux contraintes. En particulier, il estime que la fonction prototype doit être à support borné en fréquence. 5.3. Principes généraux de l'invention L'invention repose sur une approche tout à fait nouvelle des signaux multiporteuses du type OFDM/OQAM, selon laquelle l'orthogonalité est obtenue non plus par le respect des deux contraintes mentionnées ci-dessus, mais par une définition spécifiques des fonctions prototypes. En d'autres termes, l'invention a pour objet de nouveaux signaux, basés sur des systèmes de modulation construits comme l'OFDM/OQAM sur un réseau orthogonal de densité 2, sans pour autant que la fonction prototype soit à support borné en fréquence. Le principe utilisé est de construire des réseaux orthogonaux de densité 1/2, puis d'en déduire des réseaux de densité 2 par un choix judicieux des phases des signaux. De très nombreux signaux peuvent être construits selon la technique de l'invention. On donne ci-dessous deux exemples non limitatifs de tels signaux, appelés respectivement OFDM/MSK et OFDM/IOTA. Une méthode particulière pour construire de tels signaux est également donnée, à titre d'exemple non limitatif, en annexe 3. Cette méthode fait bien sûr partie de l'invention, et n'a été renvoyée en annexe que pour simplifier la lecture de la présente description. 5.4. La modulation OFDM/MSK Nous considérons ici une nouvelle modulation construite selon la même équation générique que la modulation OFDM/OQAM (équations 14 et 15), mais à partir d'une fonction prototype différente. Nous la baptiserons OFDM/MSK parce que chaque porteuse est modulée en MSK [20]. La fonction prototype s'écrit: Si t En fait, on constate a posteriori que cette modulation peut être considérée comme duale de l'OFDM/OQAM, puisqu'elle correspond à un échange des axes temps et fréquence. L'intérêt essentiel de cette modulation par rapport à l'OFDM'OQAM est que la fonction prototype est strictement limitée dans le temps, ce qui simplifie notablement l'implémentation du récepteur, puisque le nombre de coefficients du filtre d'entrée est considérablement réduit.Par ailleurs, les performances en présence de trajets multiples sont inchangées, le paramètre 5 étant identique. 5.5. La modulation IOTA La modulation OFDM/IOTA résulte en revanche d'une approche totalement nouvelle et originale dans le domaine du traitement de signal que nous baptisons transformée IOTA (pour "Isotropic Orthogonal Transform Algorithm"), et décrite en annexe 3. 5.5.1. Equation du signal Nous considérons ici une nouvelle modulation construite selon la même équation générique que la modulation OFDM/OQAM (équations 14 et 15), mais à partir d'une fonction prototype différente. Nous la baptiserons OFDM/IOTA en raison du choix de la fonction prototype. La fonction prototype s'écrit: désignant la fonction IOTA définie en annexe 3. On notera que la méthode de construction donnée en annexe 3 permet d'obtenir une infinité de solutions, la fonction IOTA constituant une solution remarquable. Les fonctions de base de la modulation OFDM/IOTA s'écrivent donc: Le signal émis s'écrit donc: avec: 5.5.2. Commentaires des figures et avantages liés à la décroissance rapide Afin de mieux mettre en évidence, de façon visuelle, les avantages de l'invention, on présente pour chaque modulation discutée précédemment: .A la fonction prototype x(t); . B la transformée de Fourier en linéaire de la fonction prototype; . C le module de la fonction d'ambiguïté en linéaire (telle que définie dans l'annexe 2); .D la fonction d'intersymbole (telle que définie dans l'annexe 2). Les vues de la fonction d'ambiguïté (figures indicées C) permettent de juger du confinement dans le plan temps-fréquence de la fonction prototype. Les vues de la fonction d'intersymbole (figures indicées D) permettent d'apprécier la sensibilité d'une modulation au délai et au Doppler. Les erreurs de phase ne sont pas considérées, puisque toutes les modulations sont équivalentes sur ce plan. Les figures 3A à 3D concernent le cas connu de l'OFDM/QAM classique. Le défaut principal de cette modulation n'est pas, comme pourrait le faire penser la réponse en fréquence de la fonction prototype, la faible décroissance du niveau des lobes secondaires. En fait, la sensibilité de l'OFDM aux erreurs fréquentielles n'est que légèrement supérieure à celle des autres modulations considérées. L'IES a par contre une statistique différente, qui se traduit par une fermeture horizontale de l'oeil équivalente à celle d'une modulation à roll-off nul. n existe donc des traces, certes improbables, mais qui peuvent créer des erreurs systématiques en l'absence de codage. Ce détail est "inesthétique", mais sans conséquence réelle en présence de codage. Par contre, cette faible décroissance fait que l'énergie d'IES se répartit sur un grand nombre de symboles voisins, ce qui rend très délicate toute tentative d'égalisation. Paradoxalement, le véritable problème vient de la limitation brutale de la réponse temporelle, qui correspond à une fonction d'ambiguïté triangulaire selon cet axe. Ceci donne une fonction d'intersymbole avec une très forte sensibilité aux erreurs temporelles la pente est verticale et le paramètre 5 est donc nul. C'est ce qui justifie l'utilisation d'un intervalle de garde. Les figures 4C et 4D concernent l'OFDM/QAM avec un intervalle de garde (les fonction prototype et transformée de Fourier sont identiques à celles de 1'OFDM/QAM illiustrées en figures 3A et 3C. L'utilisation d'un intervalle de garde crée une zone plate au niveau de la fonction d'ambiguïté. En fait, on devrait plutôt parler dans ce cas de crossambiguïté. On retrouve évidemment cette partie plate au niveau de la fonction d'intersymbole, ce qui donne une immunité aux erreurs temporelles. Les figures représentent le cas d'un intervalle de garde 0,25 T0. Au niveau des erreurs fréquentielles, les propriétés sont les mêmes que celles de l'OFDM standard. Le coût de l'intervalle de garde est admissible lorsqu'on s'intéresse à des modulations à faible efficacité spectrale. n devient par contre rédhibitoire si l'on cherche une efficacité spectrale élevée : prenons par exemple un intervalle de garde égal au quart du symbole utile. Dans ces conditions, il faut pour atteindre une efficacité nette de 4 bits/s/Hertz un sytème de modulation et de codage ayant une efficacité brute de S bits/s/Hertz, soit une perte de 3 dB par rapport à la capacité limite de Shannon. Encore faut-il ajouter à cette perte la perte supplémentaire de 1 dB due à la puissance "inutilement" émise dans l'intervalle de garde. Au total, ce sont donc 4 dB qui sont perdus par rapport à l'optimum. Les figures SA à SD présentent le cas de l'OFDM/OQAM. La réponse temporelle de 1'OFDM/OQAM présente une meilleure allure que celle de l'OFDM/QAM. Néanmoins la décroissance temporelle n'est qu'en l/t2. La fonction d'ambiguïté s'annule sur un réseau de densité 1/2. La sensibilité aux erreurs en fréquence est supérieure à celle aux erreurs temporelles. Le paramètre 5 vaut 0.8765. Les figures 6A à 6D concernent le premier mode de réalisation de l'invention, I'OFDM/MSK. On vérifie qu'elle présente des propriétés strictement identiques à celles de l'OQAM en inversant les échelles temporelles et fréquentielles. Le paramètre 5 est inchangé. Enfin, les figures 7A à 7D présentent la modulation OFDM/IOTA. Celleci présente une décroissance rapide (au sens mathématique du terme) en temps et en fréquence, ce qui permet d'envisager l'égalisation dans les meilleurs conditions possibles. Elle présente par ailleurs une symétrie parfaite par rapport à ces deux axes. Sa fonction d'intersymbole est quasiment idéale. D'une manière générale son comportement se rapproche de celui de la gaussienne. Le paramètre 5 vaut 0.9769. On peut comparer la fonction d'ambiguïté de la fonction S (figure 7C) à celle d'une gaussienne, telle qu'illustrée en figure 8. L'allure générale de ces deux fonctions est très similaire au niveau du sommet. Elles different par contre à la base. La figure 7E montre en échelle logarithmique la décroissance en temps du signal IOTA. On observe que l'amplitude du signal décroît linéairement en échelle logarithmique (en temps et en fréquence, bien sûr, puisque les deux aspects sont identiques), soit de façon exponentielle en échelle linéaire. Cette propriété permet donc dans une réalisation pratique de tronquer la forme d'onde et de limiter ainsi la complexité du récepteur. 5.6. Principe d'un émetteur La figure 9 présente un synoptique simplifié d'un émetteur d'un signal selon l'invention. Le procédé d'émission s'en déduit directement. On considère une source binaire à haut débit (typiquement quelques dizaines de Mégabits/s). Par source binaire, on entend une série d'éléments de données correspondant à un ou plusieurs signaux 91 source de tous types (sons, images, données) numériques ou analogiques échantillonnés. Ces données binaires sont soumises à un codage canal 92 binaire à binaire adapté à des canaux évanouissants. On pourra par exemple utiliser un code en treillis (Trellis Coded Modulation), concaténé éventuellement avec un code de Reed-Solomon. Plus précisément, si l'on souhaite une efficacité spectrale de 4 bits/Hertz, on peut utiliser un code de rendement 2/3 associé à une modulation 8AM, prenant 8 niveaux d'amplitude. Ensuite, conformément au principe exposés dans le brevet FR-88 15216, on répartit (93) ces données codées dans l'espace temps-fréquence de façon à apporter la diversité nécessaire, et à décorréler l'évanouissement (fading) de Rayleigh affectant les symboles émis. Plus généralement, on effectue un premier codage binaire à binaire, un entrelacement en temps et en fréquence et un codage binaire à coefficients, communément appelé mapping . n est clair que l'entrelacement peut être indifféremment effectué avant ou après le mapping, selon les besoins et les codes utilisés. A l'issue de cette opération de codage, on dispose des symboles réels à émettre amn. Le principe de réalisation du modulateur 94 OFDM/MSK ou OFDM/IOTA est analogue à celui d'un émetteur OFDM/OQAM. Seule la forme d'onde prototype differe. On pourra se référer à [15] pour une description détaillée du système de modulation. Pour construire le signal à émettre, on regroupe les symboles de même rang n, et l'on calcule: Cette opération peut être avantageusement réalisée sous forme numérique par une transformée de Fourier rapide (for) 1 ) portant sur tous les symboles de même rang n, suivi d'une multiplication de la forme d'onde résultante par la fonction prototype IOTA, et enfin d'une addition des symboles de rangs différents (sommation selon l'indice n). Le signal complexe ainsi généré est alors converti sous forme analogique 98, puis transposé à la fréquence finale par un modulateur 99 à deux voies en quadrature (modulateur I & ), et enfin amplifié 910 avant d'être émis 911. 5.7. Principe d'un récepteur La figure 10 illustre de façon schématique un récepteur d'un signal selon l'invention (ainsi que le procédé de réception correspondant). Le récepteur OFDM/MSK ou OFDM/IOTA est sensiblement analogue à celui adapté à la modulation OFDM/OQAM. Les étages d'entrée sont conventionnels. Le signal est préamplifié 101, puis converti en fréquence intermédiaire 102 afin de réaliser le filtrage canal 103. Le signal en fréquence intermédiaire est ensuite converti en bande de base sur deux voies en quadrature 105. En outre, on réalise les fonctions de correction automatique de gain (CAG) 104, qui contrôle la préamplification 101. Une autre solution consiste à transposer le signal en fréquence intermédiaire sur une fréquence porteuse basse, de façon à échantillonner le signal sur une seule voie, la représentation complexe étant alors obtenue par filtrage numérique. Alternativement, le signal RF peut être transposé directement en bande de base (conversion directe), le filtrage canal étant alors réalisé sur chacune des deux voies I & . Dans tous les cas, on peut se ramener à une représentation discrète du signal de l'enveloppe complexe correspondant au signal reçu. Afin de détailler le traitement numérique en bande de base, nous considérerons une modulation de type multiporteuse caractérisée par l'équation de l'enveloppe complexe du signal émis: Soit un canal de transmission caractérisé par sa fonction de transfert variable T(f,t) (voir annexe 2). L'enveloppe complexe du signal reçu r(t) s'écrit: Le démodulateur estime (106) la fonction de transfert T(f,t) par des moyens classiques, qui peuvent par exemple utiliser un réseau de référence de porteuses explicites selon le brevet FR-90 01491. Pour démoduler le signal proprement dit (107), on assimile localement le canal à un canal multiplicatif caractérisé par une amplitude et une phase correspondant à la valeur de T(f,t) pour l'instant et la fréquence considérée.Pour estimer am,n(t), le signal reçu est donc assimilé au signal: r(t) = JS(f)T(mv0,nr0)e2df = T(mv0,nr0) s(t) (33) On posera: Le démodulateur effectue donc le traitement suivant: Dans le cas d'un canal stationnaire de fonction de transfert #ei#, on retrouve évidemment: #m,n = pam,n (36) En pratique, le traitement 107 est effectué sous forme numérique, selon le procédé illustré en figure 11. Le récepteur fonctionne de façon analogue à un récepteur OFDM/OQAM [13-16]. Il effectue les traitements suivants: - multiplication 111 dudit signal reçu r(t) par la fonction prototype x(t) 112; - #repliement# 113 de la forme d'onde filtrée modulo 2to; - application 114 d'une transformée de Fourier (Fr); - correction 115 de la phase e en fonction de l'estimation du canal 116, comprenant par exemple une estimation Pm.n de la réponse en amplitude et une estimation #m,n de la réponse en phase du canal de transmission; - correction 117 de la phase correspondant au terme im+n, les éléments de données étant alternativement en phase et en quadrature;; - sélection 118 de la partie réelle du coceficient obtenu a, @ correspondant au coefficient am,n émis pondéré par la réponse en amplitude Pm,n du canal de transmission. Cet algorithme permet donc de calculer globalement tous les coefficients d'un indice n donné. L'ordre de grandeur de la complexité correspondante est approximativement le double de celle de l'algorithme utilisé pour l'OFDM/QAM. Les coefficients ainsi obtenus sont ensuite désentrelacés 108, symétriquement à l'entrelacement mis en oeuvre à l'émission, puis décodés 109, avantageusement selon une technique de décodage à décision douce, mettant par exemple en oeuvre un algorithme du type de l'algorithme de Viterbi. Si le décodage de canal tient compte de l'estimation de la réponse de l'amplitude du canal Pm n les valeurs correspondantes sont également désentrelacées 110. Par ailleurs, le désentrelacement est bien sûr effectué avant ou après le mapping, selon le moment ou l'entrelacement a été mis en oeuvre à l'émission. ANNEXE 1: REFERENCES [1] M.L. Doeltz, E.T. Heald and D.L. Martin, Binary data transmission techniques for linear systems Proceedings of the IVRE, pp. 656-661, May 1957. [2] R.R. Mosier, A data transmission system using pulse phase modulation IRE Conv. Rec. Ist Nat'l Conv Military Electronics (Washington, D.C., June 17-19, 1957) pp. 233-238. [3] G.A. Franco and G. Lachs, An orthogonal coding technique for communnications 1961 IRE Internat'l Conv. Rec., vol. 9, pp. 126-133. [4] H.F. Harmuth, On the transmission of information by orthogonal time functions AIEE Trans. (Communications and Electronics) vol. 79, pp. 248-255, July 1960. [5] S.B. Weinstein and Paul M. Ebert, Data transmission by frequency-division multiplexing using the discrete Fourier transform IEEE Trans. Commun., vol. COM-19, pp. 628-634, Oct. 1971. [6] L.J. Cimini, Analysis and simulation of a digital mobile channel using orthogonal frequency division multiplexing, IEEE Trans. Commun., vol. COM-33, pp. 665-675, July 1985. [7] E.F. Casas and C. Leung, OFDM for data communication over mobile radio FM channels - Part I: Analysis and experimental results, IEEE Trans. Commun., vol. 39, pp. 783-793, May 1991. [8] E.F. Casas and C. Leung, OFDM for data communication over mobile radio FM channels - Part li : Performance improvement, IEEE Trans. Commun., vol. 40, pp. 680-683, April 1992. [9] I. Daubechies, The wavelet transform, time-frequency localization and signal analysis, IEEE Trans. Inform. Theory, vol. IT-36, pp. 961-1005, Sept. 1990. [10] H.E. Jensen, T. Hoholdt, and J. Justesen, Double series representation of bounded signals, IEEE Trans. Inform. Theory, vol. IT-34, pp. 613-624, July 1988. [1 1] R.W. Chang, Synthesis of band-limited orthogonal signals for multi-channel data transmission, Bell Syst. tech. J., vol. 45, pp. 1775-1796, Dec. 1966. [12] B.R. Saltzberg, Performance of an efficient parallel data transmission system, IEEE Trans. Commun. Technol., vol. COM-15, pp. 805-811, Dec. 1967. [13] R.W. Chang, A theorical study of performance of an orthogonal multiplexing data transmission schème, IEEE Trans. Commun. Technol., vol. COM-16, pp. 529-540, Aug. 1968. [14] B. Hirosaki, An analysis of automatic equalizers for orthogonally multiplexed QAM systems, IEEE Trans. Commun., vol. COM-28, pp. 73-83, Jan. 1980. [15] B. Hirosaki, An orthogonally multiplexed QAM system using the discrete Fourier transform, IEEE Trans. Commun., vol. COM-29, pp. 982-989, July. 1981. [16] B. Hirosaki, A maximum likehood receiver for an orthogonally mulitplexed QAM system, IEEE Journal on Selected Areas in Commun., vol. SAC-22, pp. 757-764, Sept 1984. [17] B. Hirosaki, S. Hasegawa, and A. Sabato, Advanced group-band modem using orthogonally multiplexed QAM technique, IEEE Trans. Commun., vol. COM-34, pp. 587-592, June. 1986. [18] John A. C. Bingham, Multicarrier modulation for data transmission : An idea whose time has come, IEEE Communications Magazine, pp. 5-14, May 1990. [19] P.M. WOODWARD, Probability and information theory with application to Radar, Pergamon Press, London 1953. [20] F. Amoroso and J. A. Kivett, Simplified MSK signalling technique, IEEE Trans. Commun., vol. COM-25, pp. 433441, April 1977. [21] P. A. Bello, Characterization of randomly time-variant linear channels, IEEE Trans. Commun. Systems, pp. 360-393, Dec. 1964. [22] P.M. WOODWARD, Probability and information theory with application to Radar , Pergamon Press, London 1953. [23] M. ALARD et R. LASSALLE Principes de modulation et de codage canal en radiodiffusion numérique vers les mobiles Revue de l'U. E. R, nO 224, août 1987, pp. 168-190. ANNEXE 2 1. Modélisation du canal 1.1. Modèle général On peut considérer un canal dispersif comme un système linéaire ayant une réponse impulsionnelle variable dans le temps. n existe deux façons de définir cette réponse impulsionnelle. On s'inspirera largement des conventions proposées dans [21]: la réponse impulsionnelle à l'entrée (Input Delay Spread Function) g(t,)définie par: où s(t) et r(t) représentent respectivement les signaux émis et reçus # la réponse impulsionnelle en sortie (Output Delay Spread Function) h(t,T) définie par: On a évidemment h(t,) = g(t + T,) h(t,T) représente la réponse impulsionnelle du canal à l'instant t. Munis de ces conventions, nous pouvons définir les fonctions caractéristiques suivantes: # la fonction d'étalement délai-Doppler U(,v) (Delay-Doppler Spread Function) caractérisée par: # la fonction d'étalement Doppler-délai V(v,l) (Doppler-Delay Spread Function) caractérisée par: On a tout simplement: V(v,#)=ei2#v#U(#,v) # la fonction de transfert variable (Time-Variant Transfert Function) T(f,t) caractérisée par: On retrouve donc la même équation que dans le cas d'un canal stationnaire, la différence étant simplement que la fonction de transfert devient variable dans le temps. Cette fonction de transfert T(f,t) est la transformée de Fourier bidimensionnelle de U(T,V), soit: Dans tous les cas, nous considérerons que U(T,V) est à support borné, ce qui permet de représenter la fonction de transfert T(f,t) par un réseau de valeurs discrètes en vertu du théorème d'échantillonnage. 1.2. Le modèle délai-Doppler statique Le modèle délai-Doppler est défini par l'équation: Cette équation fait apparaître le canal comme une somme de canaux élémentaires caractérisés par une amplitude, une phase, un offset temporel et un offset fréquentiel. Aussi est-il légitime de s'intéresser au comportement des différentes modulations en présence de ce type de canal, que nous baptiseront délai-Doppler statique. L'équation du canal s'écrit alors sous la forme simplifiée suivante: r(t) = Aeios(t- z)e'2" 2. Performances de l'OFDM dans les canaux non-stationnaires 2.1. Cas général Considérons une modulation multiporteuse OFDM de type quelconque (OFDM/QAM, OFDM/OQAM ou OFDM/IOTA) caractérisée par l'équation générique: ak étant des variables réelles, E étant un réseau bidimensionnel de densité 2 dans l'espace temps-fréquence, les fonctions xk(t) étant des translatées en temps et en fréquence d'une même fonction prototype x(t), et qui constituent une base hilbertienne de L2(R). On notera qu'aucune hypothèse n'est faite sur la structure du réseau E. Dans le cas particulier de l'OFDM/QAM, ce réseau se décompose en deux sous-réseaux co-localisés avec des phases en quadrature. L'opération de démodulation s'écrit: étant une phase estimée par le démodulateur et r(t) l'enveloppe complexe du signal reçu. On peut donc écrire: Or: On en déduit que: La valeur optimale de est celle qui maximise le coefficient an, soit: Bien que générales, les équations ci-dessus ne sont guère manipulables. Elles montrent cependant que le signal utile et l'intersymbole apparaissent comme des intégrations de la fonction d'ambiguïté pondérée par la fonction d'étalement délai Doppler. 2.2. Cas du canal statique Si l'on s'intéresse à un canal de type délai-Doppler statique, caractérisé par une phase 0, un retard X et un offset v (on normalisera à 1 l'amplitude A), on effectue la démodulation de manière similaire en introduisant dans l'estimateur un paramètre de phase . Le résultat de cette opération s'écrit: Le signal démodulé s'écrit donc finalement: Le deuxième terme représente l'interférence entre symboles (IES). Si l'on considère les données ak comme des variables aléatoires indépendantes de variance a2, la variance I de l'IES s'écrit:: Or les coefficients ck sont les coefficients de la décomposition de la fonction ei(#-#)e-2i@v(#+#)xn(@ + T), de norme unité, sur la base hilbertienne des fonctions xk(t). On a donc: En d'autres termes, la variance du signal reçu est constante et se répartit entre le signal "utile" cnan et l'IES, de variance I = (1 - cn2)a2.Le calcul du coefficient cn donne: Or la fonction d'ambiguïté de xn s'écrit: AXn(#,v) = e2i#(vn#-#nv)Ax(#,v) Finalement, on peut écrire: On considérera que la phase de démodulation # s'écrit sous la forme #opt+ ##, où @pt est la phase de démodulation qui minimise l'IES, c'est-à-dire qui maximise cn, soit: Sept = 8 + rvT + 2#(#nv-vn#) Alors, la variance de l'IES s'écrit simplement:: i = (1 - (e[Ax(T,v)e])2)a2 Lorsque la fonction prototype est paire, (ce qui correspond au cas de la méthode de construction de bases hilbertiennes décrite dans le texte principal), la fonction d'ambiguïté est réelle, et on a donc: I = (1 - AX2(T,V)COS2 i ))a2 Ce résultat est tout à fait remarquable, puisqu'il démontre que la sensibilité au délai et au Doppler de toute modulation multiporteuse ne dépend que de la fonction d'ambiguïté de sa fonction prototype. On appellera par la suite fonction d'intersymbole (par abus de langage, pour fonction d'interférence entre symboles) la fonction dans le cas général, qui représente la valeur quadratique moyenne de l'intersymbole normalisée par la valeur quadratique moyenne des données dans le cas d'une estimation de phase optimale. 3. Analvse comparative des différents types d'OFDM 3.1. Limites théoriques Nous nous intéressons ci-après aux propriétés de la fonction d'intersymbole. On constate que la sensibilité d'une modulation multiporteuse est directement liée au comportement de la fonction d'ambiguïté de la fonction prototype correspondante au voisinage de (0,0). Le problème posé est tout à fait analogue aux problèmes d'incertitude rencontrés dans le domaine radar, et on pourra se référer à la littérature abondante sur le sujet (voir par exemple [22]).Sans perte de généralité, on peut choisir la fonction x(t), par une translation temporelle et fréquentielle adéquate, de telle sorte que ses moments d'ordre un soient nuls, soit: Dans ces conditions, on vérifie aisément que les dérivées partielles du premier ordre s'annulent: On peut caractériser le comportement de la fonction d'ambiguïté autour de (0,0) à partir des dérivées partielles de second ordre: On posera Considérons le développement de Taylor-Young de la fonction d'ambiguïté en (0,0): Ax(d#,dv)=1-2#2(#t2dv2+#f2d#2)+ dvd#+o(dv2+d#2) On en déduit le développement de Taylor-Young de la variance de l'intersymbole I=(l-(#e[Ax(#,v)])2 cos2 ##)#2, soit: : I(d#,dv,d#)=#2[4#2(#t2dv@+#f2d#2)-2 dvd#+d#2+o(dv2+d#2+d#2)] On en déduit que la fonction d'intersymbole Is admet à l'origine un cône tangent d'équation: L'intersection de ce cône avec le plan z = 1 (intersymbole maximum) délimite une surface de contour elliptique dont l'aire # peut être considérée comme une mesure de la sensibilité au délai et au Doppler. Lorsque x est nul, cette ellipse a pour axes de symétrie les axes temporel et fréquentiel, et s'étend de +1/ 2Mf selon l'axe temporel et #1/ 2Xt selon l'axe fréquentiel. On a donc: #=1/4##t#f En vertu de l'inégalité de Heisenberg, # ne peut pas dépasser l'unité. Ce résultat se généralise au cas où x est différent de 0.Considérons la fonction y(t) obtenue en multipliant la fonction x(t) par une wobulation: On peut donc écrire: Il est donc toujours possible d'annuler lly en choisissant ss de façon appropriée. Or l'opération de multiplication par une wobulation réalise un simple changement d'axes de la fonction d'ambiguïté associée, avec conservation des aires. On en déduit le paramètre # est donc toujours compris entre 0 et 1. Ce résultat est extrêmement important, puisqu'il permet de comparer les performances de toutes les MCM dans des canaux dispersifs à partir d'un paramètre unique. On constate donc que ces performances ne dépendent que de la concentration de la fonction prototype associée. L'optimum est atteint virtuellement par la gaussienne, mais cet optimum est inaccessible, puisque les gaussiennes ne permettent pas de construire une base hilbertienne. ANNEXE 3 1. Introduction Cette annexe donne une méthode de construction de fonction prototypes vérifiant les critères d'orthogonalité requis. La méthode permet d'obtenir une infinité de fonctions, parmi lesquelles une solution particulière (appelée fonction IOTA) possédant la particularité d'être identique à sa transformée de Fourier. 2. Fonction d'ambiguïté Ce chapitre rappelle les principales propriétés de la fonction d'ambiguïté, et décrit différents opérateurs agissant sur cette fonction. 2.1 Rappels sur la fonction d'ambiguïté 2.1.1 Définitions Soit une fonction x(t) et sa transformée de Fourier X(f). On peut lui associer ses produits temporel et fréquentiel définis respectivement par: Yx(t,T) = x(t+ r/2) x'(t- r/2) @x(f,v)=X(f+v/2)X*(f-v/2) La transformée de Wigner-Ville et la fonction d'ambiguïté de x sont alors données par: 2.1.2. Propriétés de symétrie de la fonction d'ambiguïté Soit une fonction x(t).On notera respectivement par x et x les fonctions définies de la manière suivante: x@(t)=x(-t) x(t) = x*(-t) On a alors les relations: soit, en posant u = -t: On en conclut en particulier que si une fonction x est paire, c'est à dire que x = x-, sa fonction d'ambiguïté est réelle. Par ailleurs, on notera la relation suivante: En combinant ces deux relations, on obtient: A(r,v) = A(r,-v) 2.1.3. Fonction d'ambiguïté et transformée de Fourier On peut réécrire la définition de la fonction d'ambiguïté de la façon suivante: ou encore :Ax (T, v) = Ax(v,#) 2.1.4.Fonction d'ambiguïté et translation temps fréquence Considérons une fonction translatée d'une fonction prototype x(t) quelconque, soit: La fonction d'ambiguïté associée s'écrit: soit, en posant u = t - r, 2.2 Orthogonalité et fonction d'ambiguïté 2.2.1 Cas général On considère deux fonctions translatées d'une même fonction x(t), soit: Le produit scalaire de ces deux fonctions s'écrit: soit, en posant u=t-(#k+#k')/2: 3. Bases hilbertiennes sur des réseaux orthogonaux 3.1. Principes généraux de construction On considère un ensemble de fonctions {Xm,n} défini par: On cherche les conditions sur x(t) pour que cet ensemble {Xm,n} constitue une base hilbertienne de HR.On imposera que x(t) soit une fonction paire, dont la fonction d'ambiguïté Ax est donc réelle. Le produit scalaire dans R de Xm,n et de Xm',n' n' peut s'écrire On notera la relation de congruence modulo 2: (m-m')+(n-n')+(m-m')(n+n')#1-(m-m'+1)(n-n'+1) Par conséquent, si (m,n) #(m',n') modulo 2, le produit scalaire est nul. Le réseau {Xm,n} peut donc se décomposer en quatre sous-réseaux caractérisés par: {m pair, n pair},{m pair, n impair},{m impair n pair},{m impair, n impair}. L'orthogonalité entre fonctions appartenant à des sous-réseaux différents est donc automatique, et ne dépend pas des propriétés de la fonction prototype, dès l'instant où celle-ci est paire. I1 reste ensuite à garantir que les fonctions d'un même sous-réseau soient orthogonales entre elles. I1 suffit pour cela que la fonction d'ambiguïté Ax vérifie: Ax(2n#0,2mv0)= O V(m,n) $ (0,0) On constate donc que le problème de la construction de bases hilbertiennes de HR sur un réseau orthogonal de densité 2 se ramène à celui de la construction d'une fonction prototype paire dont la fonction d'ambiguïté s'annule sur un réseau de densité 1/2. 3.2. Méthodes d'orthogonalisation 3.2.1. Onhogonalisation temporelle Définition: Soit une fonction x(t) de transformée de Fourier X(f). On appelle O@ l'opérateur d'orthogonalisation temporelle qui associe à x(t) une fonction y(t) définie par sa transformée de Fourier Y(f): Par construction, on a: soit, par transformée de Fourier inverse: soit encore Ay(2n#0,0)=0 #n#0 et Ay(0,0)=1 On réalise donc bien l'orthogonalisation sur l'axe temporel. On note en outre que cet opérateur normalise y. Soit x une fonction gaussienne et y = Otx. Considérons l'expression: Puisque X est une gaussienne, on peut écrire: X(f + mvO)X (f - mvO) = cm|X(f)|2 où cm est une constante. On en déduit que: #y(f,2mv0)=cm#y(f,0) Par transformée de Fourier inverse, on obtient: Ay(#,2mv0)=cmAy(#,0) Par conséquent: #m,#n##0 Ay(2n#0,2mv0)=0 L'opérateur d'orthogonalisation temporelle Ot orthogonalise donc l'ensemble du réseau, à l'exception de l'axe des fréquences. Théorème 1 Soit x une fonction gaussienne et y = Otx, alors: #m,#n##0 Ay(2n#0,2mv0)=0 3.2.2. Orthogonalisation fréquentielle Définition Soit une fonction x(t). On appelle Of l'opérateur d'orthogonalisation fréquentielle qui associe à x(t) une fonction y(t) définie par: Par construction, on a: soit, par transformée de Fourier: soit encore Ay(0, 2mv0)=0 #m#0 et Ay(0,0)=l On réalise donc bien l'orthogonalisation sur l'axe fréquentiel. On note en outre que cet opérateur normalise y. Soit x une fonction gaussienne et z = Ofy, avec y = Otx. Considérons l'expression: On peut donc écrire: &gamma;z(t,2n#0)-&gamma;y(t,2n#0)P(t) où P(t) est une fonction périodique de période 0, qui admet un développement en série de Fourier du type Par transformée de Fourier, on obtient: Or De plus, par construction, #m#0, Az(0, 2mv0)=0 On a donc finalement: #(m,n)#(0,0), Az(2n#0,2mv0)=0 Ainsi, la fonction d'ambiguïté de z s'annule en dehors de (0,0) pour tous les multiples de 210 et de 2v0, soit un réseau de densité 1/2. Théorème 2 Soit x une fonction gaussienne et z = OfOtX, , alors: #(m,n)#(0,0), Az(2n#0,2mv0)=0 3.3. L'opérateur d'orthogonalisation O Au vu de ce qui précède; il apparaît clairement qu'il existe une échelle tempsfréquence qui symétrise l'écriture des équations: Il suffit pour cela de choisir #0 = v0 = 1/ W. On renormalisera donc les échelles en conséquence, sans que cela nuise à la généralité des démonstrations. 3.3.1. Définition On appelle O l'opérateur d'orthogonalisation qui associe à une fonction x la fonction y définie par: En outre, on appellera par la suite F l'opérateur de transformée de Fourier. 3.3.2. Idempotence de l'opérateur O Soit z = Oy et y = Ox. On peut écrire: On a donc OOx = Ox, ce qui démontre l'idempotence de l'opérateur O. De la même façon, l'opérateur dual F-1OF est également idempotent, puisque F-1OFF-1OF = F-1OOF = F-1OF. 3.3.3. Lemme 1 Soit P une fonction périodique de période 1/ 4 et D une distribution de la forme: Soit x une fonction quelconque: Lemme 1 Soit P une fonction périodique de période 1/ et D une distribution de la forme Soit x une fonction quelconque.On a: D * (P x)=P(D * x) Lemme 2 Soit la fonction y définie par Ya = D * xa , avec x&alpha;=(2&alpha;)1/4 e~ , et D étant une distribution de la forme On peut donc écrire: Considérons la somme: Soit encore, par application du résultat donné en appendice (≈4): puis, en réorganisant les indices et en redéfinissant k comme k + k' + k": On peut donc écrire: avec On peut estimer aisément le coefficient c en réécrivant la relation précédente sous la forme: Soit, par transformée de Fourier: En particulier, on en déduit: ||y&alpha;||2=Ay&alpha;(0,0)=cAx&alpha;(0,0)=c||x&alpha;||2 On a donc finalement:: Lemme 2 Soit la fonction Ya définie par Ya = D * x&alpha;, avec x&alpha;=(2&alpha;)1/4 #-#&alpha;u2, et D étant une distribution de la forme On peut écrire: 3.3.5. Commutativité des opérateurs O et F-1OF Nous allons maintenant démontrer que les opérateurs O et F-1OF commutent lorsqu'ils sont appliqués à une gaussienne. Soit x&alpha;=(2&alpha;)1/4e-#&alpha;u2. Alors Fxa = x1/&alpha; et O xa = PaXa P&alpha; étant défini par la relation: et sa transformée de Fourier D, par: Soit Ya = F-1OFx&alpha;et za =Oy&alpha;. On peut écrire: y&alpha;=F-1OFx&alpha;=F-1Ox1/&alpha;=F-1(P1/&alpha;x1/&alpha;)=D1/&alpha;*x&alpha; Xa et Ya étant de norme unité, on peut écrire, en application du lemme 2: De la même façon, on définit: w1/&alpha;=FOX&alpha;=F(P&alpha;x&alpha;)=D&alpha;*x1/&alpha; On peut écrire: x1/&alpha; @ &alpha; et w" a étant de norme unité, on a, en application du lemme 2: Soit, par transformée de Fourier inverse: F-1OFOx&alpha;=F-1Ow1/&alpha;=D1/&alpha;*(P&alpha;x&alpha;) Or, par application du lemme 1: Dua* (Paxa ) = Pa (D1.&alpha;* xa ) On en déduit que: OF-1OFx&alpha;-F-1OFOx&alpha; Théorème 3 Pour toute fonction gaussienne x, les opérateurs O et F-tOF commutent, soit: OF-1OFx=F-1OFOx Corollaire 1 Soit z&alpha;=OF-1OFx&alpha;, avec x&alpha; = (2&alpha;)1/4 e-#&alpha;u2, alors Fz&alpha;=z1/&alpha; Démonstration:: Fz&alpha;=FF-1OFOx&alpha;=OF-1Ox&alpha;=OF-1OFx1/&alpha; = z1/&alpha; Cas particulier remarquable Fzl-zl Cette fonction particulière confère une parfaite symétrie aux axes temps et fréquence, et constitue donc la fonction prototype de la transformée IOTA (Isotropic Orthogonal Transform Algorithm). On notera cette fonction particulière 9. Corollaire 2 Soit x une fonction gaussienne et z = OF-1OFx, alors Oz = z Démonstration: Oz = OOF-1OFx = OF-1OFx = z Corollaire 3 Soit x une fonction gaussienne et z = OF-1OFx, alors F-1OFz = z Démonstration: F-1OFz=F-1OFF-1OFOx=F-1OOFOx=F-1OFOX=z 3.3.6. Fonction d'ambiguïté des fonctions z Considérons le théorème 2, avec la normalisation #0 = V0 = l/W2. Alors: Of-OetOt=F-1OF Par conséquent, le théorème 2 peut se réécrire: Théorème 4 Soit x une fonction gaussienne et z = F-1OFOx, alors: #(m,n)#(0,0), Az(n#2,m#2)=0 4. Appendice Soit une fonction gaussienne normalisée Xa définie par: x&alpha;(u)=(2&alpha;)1/4e-#&alpha;u2 Le produit x&alpha;(u-a)x&alpha;(u-b) peut done s'écrire: Or, on a l'identité: a+b a-b (u-a)2+(u-b)2=2##u-#2+##2# 2 2 Finalement, on peut écrire: REVENDICATIONS 1. Signal multiporteuse destiné à être transmis vers des récepteurs numériques, notamment dans un canal de transmission non-stationnaire, correspondant au multiplexage en fréquence de plusieurs porteuses élémentaires correspondant chacunes à une série de symboles, deux symboles consécutifs étant séparés d'un temps symbole 0, caractérisé d'une part en ce que l'espacement v0 entre deux porteuses voisines est égal à la moitié de l'inverse du temps symbole 0, et d'autre part en ce que chaque porteuse subit un filtrage de mise en forme de son spectre présentant une largeur de bande strictement supérieure à deux fois ledit espacement entre porteuses v0, et choisi de façon que chaque symbole soit fortement concentré dans le domaine temporel et dans le domaine fréquentiel. 2. Signal selon la revendication 1, caractérisé en ce que son enveloppe complexe répond à l'équation suivante: où: a.n est un coefficient réel, choisi dans un alphabet de modulation prédéterminé; m est un entier représentant la dimension fréquentielle; n est un entier représentant la dimension temporelle; t représente le temps; x",n(t) est une fonction de base, translatée dans l'espace temps-fréquence d'une même fonction prototype x(t) paire prenant des valeurs réelles ou complexes, soit:: où (p est un paramètre de phase arbitraire, la transformée de Fourier X(f) de la fonction x(t) ayant un support s'étendant au delà de l'intervalle [-v0,v0], et où lesdites fonctions de base {xmn) sont orthogonales entre elles, la partie réelle du produit scalaire de deux fonctions de base différentes étant nulle. 3. Signal selon la revendication 2, caractérisé en ce que ladite fonction prototype x(t) est une fonction paire, nulle en dehors de l'intervalle [-T0,T0], et vérifiant la relation: x(t)=0 si|t|##0 # |x(t)|2+|x(t-#0)|2 =l/#0 si 0 #t 4. Signal selon la revendication 3, caractérisé en ce que ladite fonction prototype x(t) est définie par: Si t S ailleurs 5. Signal selon la revendication 2, caractérisé en ce que ladite fonction prototype x(t) est caractérisée par l'équation: la fonction y(t) étant définie par sa transformée de Fourier Y(f): où G(f) est une fonction gaussienne normalisée du type: G(f) = (2a)"4e-t ,a étant un paramètre réel strictement positif. 6. Signal selon la revendication 5, caractérisé en ce que le paramètre a est égal à l'unité. 7. Procédé de transmission d'un signal numérique, notamment dans un canal de transmission non-stationnaire, caractérisé en ce qu'il comprend les étapes suivantes: - codage canal d'un signal numérique à transmettre, délivrant des coefficients numériques réels am,n choisis dans un alphabet prédéterminé; - construction d'un signal s(t) répondant à l'équation suivante: où: . m est un entier représentant la dimension fréquentielle; . n est un entier représentant la dimension temporelle; . t représente le temps; . xm,n(t) est une fonction de base, translatée dans l'espace temps-fréquence d'une même fonction prototype x(t) paire prenant des valeurs réelles ou complexes, soit:: où (p et un paramètre de phase arbitraire, la transformée de Fourier X(f) de la fonction x(t) ayant un support s'étendant au-delà de l'intervalle [-v0,v0], lesdites fonctions de base [xm,n] étant orthogonales entre elles, la partie réelle du produit scalaire de deux fonctions de base différentes étant nulle; - émission d'un signal ayant pour enveloppe complexe ledit signal s(t) vers au moins un récepteur. 8. Procédé selon la revendication 7, caractérisé en ce qu'il comprend une étape d'entrelacement en fréquence et/ou en temps, appliquée aux éléments binaires formant ledit signal numérique à transmettre ou aux coefficients numériques amn issu du codage canal. 9. Procédé de réception d'un signal selon l'une quelconque des revendications 1 à 7, caractérisé en ce qu'il comprend les étapes suivantes: - réception d'un signal ayant pour enveloppe complexe un signal r(t); - estimation de la réponse du canal de transmission, comprenant une estimation de la réponse en phase 6m,n et de la réponse en amplitude Pm n; - démodulation dudit signal r(t), comprenant les étapes suivantes:: - multiplication dudit signal r(t) par la fonction prototype x(t); - repliement de la forme d'onde obtenue modulo 210; - application d'une transformée de Fourier (Wr); - correction de la phase Omn induite par le canal de transmission; - correction de la phase correspondant au terme im+n; - sélection de la partie réelle du coeeficient obtenu am,n correspondant au coefficient am n émis pondéré par la réponse en amplitude #m,n du canal de transmission. 10. Procédé selon la revendication 9, caractérisé en ce qu'il comprend une étape de désentrelacement en fréquence et/ou en temps desdits coefficients numériques réels am.n et, éventuellement, des valeurs correspondantes #m,n de la réponse en amplitude du canal, ledit désentrelacement étant symétrique à un entrelacement mis en oeuvre à l'émission, et/ou une étape de décodage en décision pondérée adapté au codage canal mis en oeuvre à l'émission. 11. Procédé de construction d'une fonction prototype x(t) pour un signal selon l'une quelconque des revendications 2 à 6, caractérisé en ce qu'il comprend les étapes suivantes - sélection d'une fonction gaussienne G(f) normalisée du type: G(f)=(2&alpha;)1/4e-#&alpha;f - détermination de ladite fonction prototype x(t) telle que: la fonction y(t) étant définie par sa transformée de Fourier Y(f):