FR 2977429 A1 20130104 FR 1155867 A 20110630 Procédé de réception d'un signal multiporteuse de type OQAM, dispositif de réception et programme d'ordinateur correspondants. 1. Domaine de l'invention Le domaine de l'invention est celui des communications numériques. Plus précisément, l'invention concerne la transmission de signaux multiporteuses de type OFDM/OQAM (en anglais « Orthogonal Frequency Division Multiplexing / Offset Quadrature Amplitude Modulation ») ou BFDM/OQAM (en anglais « Biorthogonal Frequency Division Multiplexing / OQAM »), pour lesquels les porteuses sont mises en forme par un filtre prototype. Encore plus précisément, l'invention concerne la synchronisation des récepteurs, et notamment l'estimation des décalages en phase et/ou en temps introduits par le canal de transmission et affectant de tels signaux multiporteuses,. L'invention trouve notamment des applications dans le domaine des communications sans fil (DAB, DVB-T, WLAN, optique non guidée, etc) ou filaires (xDSL, PLC, optique, etc). 2. Art antérieur Les techniques de transmission à porteuses multiples présentent de nombreux avantages, notamment dans le contexte de canaux multi-trajets. Ainsi, les modulations de type OFDM sont particulièrement bien adaptées pour contrer les effets des évanouissements dans les canaux sélectifs en fréquence. Toutefois, ces modulations OFDM présentent l'inconvénient de générer un signal multiporteuse présentant une mauvaise localisation fréquentielle, et nécessitent donc l'introduction d'un intervalle de garde dans le domaine temporel pour limiter les interférences. Or l'insertion d'un tel intervalle de garde engendre une diminution de l'efficacité spectrale du signal multiporteuse. Des solutions alternatives ont alors été proposées pour limiter les interférences tout en s'affranchissant de l'insertion d'un intervalle de garde. Ces techniques reposent sur la mise en forme du signal par des filtres (pour un signal discrétisé) ou des fonctions (pour un signal continu), dits prototypes, permettant une meilleure localisation fréquentielle grâce à des propriétés d'orthogonalité restreintes au corps des réels. Il s'agit par exemple des modulations de type OFDM/OQAM ou BFDM/OQAM, classiquement utilisées pour les communications radiofréquences, telles que décrites notamment dans les documents "Analysis of OFDM/OQAM systems based on the filterbank theory", P. Siohan et N. Lacaille, Proc. GLOBECOM'99, Rio de Janeiro, Brazil, Dec. 1999, pp. 2279-2284, et "Design of BFDM/OQAM systems based on Biorthogonal modulated filter Banks", C. Siclet et P. Siohan, Proc. GLOBECOM'00, San Francisco, USA, Nov. 2000, pp. 701-705. La figure 1 illustre plus précisément le schéma d'un modulateur OFDM/OQAM 10 délivrant un signal multiporteuse de type OFDM/OQAM. Un tel signal peut être représenté, en bande de base et en temps discret, sous la forme suivante : M-1n=+00 j27rm(k-D/2) s[k] = 1 amnh[k-nN]e M m=0 n=-.0 avec : - am n un symbole de données à valeur réelle à transmettre sur une sous-porteuse m à l'instant n ; - M le nombre de fréquences porteuses ; - N = M / 2 un décalage temporel discret ; - h[n] le filtre prototype utilisé par le modulateur, de longueur L ; - D le délai introduit par le filtre prototype ; - 0m n un terme de phase choisi de manière à réaliser une alternance partie réelle- partie imaginaire permettant l'orthogonalité, par exemple égal à (n+m) +e7Lmn,avec s={0,+1} ; - k correspond à l'instant kTs, avec Ts la fréquence d'échantillonnage. Dans le cas d'une modulation OFDM/OQAM orthogonale, le délai D est tel que D = L -1, avec L la longueur du filtre prototype. Dans le cas d'une modulation BFDM/OQAM biorthogonale, le délai D peut être choisi d'une manière plus flexible, et peut être tel que DL-1. Comme illustré en figure 1, les symboles de données réels am n subissent un pré-traitement ou une pré-modulation 11, permettant notamment d'assurer une quadrature en temps et en fréquence des porteuses du signal multiporteuse. Plus précisément, au cours de ce pré-traitement, les symboles de données réels am n sont multipliés par un premier terme de phase en 7r / 2 , permettant d'assurer un déphasage en temps et en fréquence des porteuses du signal multiporteuse, et par un deuxième terme permettant de tenir compte de la longueur du filtre prototype. Les symboles de données obtenus en sortie de ce module de pré-traitement, notés am n, peuvent s'exprimer sous la forme suivante : .7t D-N n -j -m am,n = am,ne 2 2Ne 2 N Ces symboles de données sont ensuite convertis du domaine fréquentiel vers le domaine temporel, en utilisant classiquement une transformée de Fourier inverse discrète (en anglais IDFT pour «Inverse Discrete Fourier Transform») de taille M, dans un module IDFT 12. Les symboles transformés obtenus en sortie du module IDFT 12, notés u0,n à u2N-1,n sont alors filtrés par le filtre prototype h[n] 13, puis sur-échantillonnés et décalés pour obtenir le signal multiporteuse OFDM/OQAM s[k]. Plus précisément, le filtre prototype h[n~ peut s'exprimer sous sa forme polyphase, comprenant M composantes polyphases Gl (z) , définies par : Gl(z)=1h[l+nM]z-n . n Cependant, malgré le gain qu'offrent les modulations de type OQAM (OFDM/OQAM ou BFDM/OQAM) par rapport aux modulations de type OFDM au niveau de l'efficacité spectrale, du fait qu'elles permettent de s'affranchir de l'insertion d'un intervalle de garde ou d'un préfixe cyclique, les modulations OFDM sont souvent préférées car leur mise en oeuvre est plus simple. En particulier, au niveau de la réception, la présence d'un préfixe cyclique pour les modulations OFDM rend l'opération de synchronisation des récepteurs peu complexe, du fait de la corrélation qu'un tel préfixe crée entre le début et la fin de chaque symbole multiporteuse OFDM. L'absence d'un tel intervalle de garde pour les modulations multiporteuses de type OQAM rend l'opération de synchronisation plus complexe, et nécessite soit l'introduction de données supplémentaires, dites pilotes, qui se présentent sous forme de préambule en début de trame ou sous forme de données isolées réparties dans la trame, soit des calculs coûteux dans le cas d'une synchronisation aveugle. En particulier, dans le cas d'un système de transmission utilisant des pilotes, l'opération de synchronisation est mise en oeuvre en effectuant une corrélation du signal multiporteuse reçu dans le domaine temporel. Cette corrélation peut, par exemple, être créée à partir de pilotes déterministes groupés sous forme de préambule, ou encore à partir de pilotes répartis basés sur des séquences pseudo-aléatoire. Le document «A data-aided symbol timing estimation algorithm for OFDM/OQAM systems » de T. Fusco, A. Petrella, et M. Tanda propose notamment une technique de synchronisation temporelle d'un récepteur OFDM/OQAM, basée sur l'utilisation de préambule. Bien que cette technique offre de bons résultats en termes de performance et de complexité, elle impose l'utilisation d'un préambule spécifique, selon lequel les porteuses impaires sont nulles, pour établir une relation entre les échantillons du signal multiporteuse. En particulier, cette relation est obtenue pour un système de transmission non-causal, et donc difficilement implémentable. Les sorties de l'étape de transformée de Fourier inverse, pour ce préambule spécifique, sont donc liées comme celles d'une transformée de Fourier inverse avec des entrées réelles et, par suite, il faut un filtre prototype de longueur impaire pour conserver ces relations après filtrage. Un autre inconvénient est que la durée du préambule augmente avec la longueur du filtre prototype, ce qui engendre une diminution de l'efficacité spectrale pour des filtres prototypes longs. Il existe donc un besoin pour une nouvelle technique de transmission de signaux 35 multiporteuses de type OFDM/OQAM ou BFDM/OQAM, permettant de simplifier l'opération de synchronisation des récepteurs par rapport aux techniques de l'art antérieur. 3. Exposé de l'invention L'invention propose une solution nouvelle qui ne présente pas l'ensemble de ces inconvénients de l'art antérieur, sous la forme d'un procédé de réception d'un signal multiporteuse de type OQAM mettant en oeuvre une étape d'estimation, dans le domaine temporel, d'un décalage en temps et/ou en phase du signal multiporteuse, mettant en oeuvre au moins un estimateur pour l'estimation du décalage en temps, dit estimateur temporel, et/ou au moins un estimateur pour l'estimation du décalage en phase, dit estimateur de phase. Selon l'invention, le signal multiporteuse comprenant au moins un préambule, l'estimateur temporel et/ou l'estimateur de phase tient compte des coefficients d'un filtre prototype utilisé en émission pour mettre en forme le préambule. Ainsi, l'invention repose sur une approche nouvelle et inventive de la synchronisation des récepteurs utilisés pour démoduler des signaux multiporteuses de type OFDM/OQAM ou BFDM/OQAM, appelés plus généralement signaux multiporteuses OQAM, présentant une complexité réduite par rapport aux techniques de l'art antérieur. En particulier, l'utilisation de tels estimateurs permet de réduire de moitié la complexité de l'étape d'estimation, un estimateur selon l'invention pouvant être réduit à effectuer une somme sur M/4 termes alors qu'un estimateur selon l'art antérieur tel que décrit dans le document Fusco précité notamment nécessite d'effectuer une somme sur M/2 termes, avec M le nombre de sous- porteuses d'un symbole multiporteuse OQAM du signal multiporteuse. De plus, une telle technique de réception permet de tenir compte de la causalité du filtre prototype. Un récepteur mettant en oeuvre une telle technique de réception est ainsi physiquement réalisable, puisqu'il est basé sur l'utilisation d'un filtre prototype causal. Or dans le document Fusco précité, le schéma de modulation proposé est non causal. Autrement dit, le paramètre de délai D n'est pas pris en considération, et par suite, les relations de conjugaison à la sortie de l'étape de transformation fréquence/temps sont identiques à celles d'une transformation fréquence/temps avec des entrées réelles si les porteuses impaires (ou paires) sont modulées par une valeur nulle. En particulier, la durée du préambule du signal multiporteuse dépend de l'estimateur selon l'invention ainsi que de la longueur de filtre prototype utilisé en émission. Afin d'éviter un allongement de la durée du préambule, il est possible de commuter entre deux filtres prototypes. En d'autres termes, différents filtres prototypes peuvent être utilisés en émission, l'un pour mettre en forme le préambule et l'autre pour mettre en forme la partie utile du signal multiporteuse. Ainsi, il est possible d'utiliser un filtre prototype très court pour mettre en forme le préambule utilisé pour la synchronisation (par exemple de longueur L = M , permettant de générer un préambule formé d'un unique symbole multiporteuse OQAM par exemple), et un filtre prototype plus long pour la partie utile. Il est également possible d'utiliser des filtres prototypes plus longs pour mettre en forme le préambule, par exemple de longueur L = qM , avec q un entier supérieur à 1. Dans ce cas, il peut être nécessaire d'ajouter un ou plusieurs colonnes de symboles de préambules nuls dans le préambule, avant les données utiles. Selon une caractéristique particulière de l'invention, l'étape d'estimation tient également compte de la position d'au moins une paire de valeurs maximales obtenues en émission pour le signal multiporteuse, ladite position étant définie par rapport aux sorties d'une étape de transformation du domaine fréquentiel vers le domaine temporel d'un ensemble de symboles de données destinés à former le préambule, dits symboles de préambule. En effet, les inventeurs de la présente demande de brevet ont mis en évidence des relations particulières de conjugaison entre les différentes sorties de l'étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission (par une transformation de Fourier rapide inverse par exemple). Du fait de la nature du signal multiporteuse OQAM et de ces relations de conjugaison, on met en évidence une symétrie dans le signal multiporteuse en sortie du modulateur, que l'on exploite en réception. L'estimation des décalages en phase et/ou en temps selon l'invention peut donc être effectuée plus rapidement et plus simplement. En particulier, le procédé de réception selon l'invention comprend une étape de suppression des échantillons temporels du signal multiporteuse présentant une amplitude inférieure à un seuil prédéterminé, mise en oeuvre préalablement à l'étape d'estimation. De cette façon, il est possible de supprimer les échantillons du signal multiporteuse reçu présentant une faible valeur, ces échantillons étant fortement perturbés par le canal de transmission. On conserve alors uniquement les échantillons du signal multiporteuse reçu correspondant aux valeurs crêtes obtenues en sortie du module de transformée fréquence/temps, ces échantillons étant peu affectés par le bruit et le canal. Connaissant les positions de ces valeurs crêtes ou maximales, il est possible de définir des estimateurs simplifiés pour l'étape d'estimation. Selon un autre aspect de l'invention, l'estimateur temporel met en oeuvre l'une des équations suivantes : M/4-1 h(k)h(M/2-1-k)r(k+2)r(M/2-1-k+z) k=0 zLSl = arg max 2 et M/4-1 h(M-1-k)h(M/2+k)r(M/2+k+2)r(M-1-k+z) k=0 et/ou l'estimateur de phase met en oeuvre l'une des équations suivantes : M/4-1 ~LSt = 2 L h(k)h(M /2 -1- k)r(k + 2)r(M /2 -1- k + II) et k=0 t-Ls2 = arg max 2 {M/4-1 h(M -1- k)h(M /2+ k)r(M /2+ k + 2)r(M -1- k + 2) k=0 avec : zl gl ou zLS2 l'estimation du décalage en temps en sortie de l'estimateur temporel ; ~LSt ou 0LS2 l'estimation du décalage en phase en sortie de l'estimateur de phase ; LH l'argument d'un nombre complexe ; M le nombre de sous-porteuses d'un symbole OQAM du signal multiporteuse ; k un entier tel que 0 k M / 4 -1 ; h(k) les coefficients dudit filtre prototype utilisé en émission ; r(k) le signal multiporteuse reçu, tel que r(k) = s(k - 2)ei(b + b(k) ; s(k) le signal multiporteuse émis ; b(k) un bruit blanc gaussien ; 2 le décalage en temps du signal multiporteuse reçu par rapport au signal multiporteuse émis ; 4 le décalage en phase du signal multiporteuse reçu par rapport au signal multiporteuse émis. Les estimateurs proposés reposent donc sur la connaissance du filtre prototype utilisé à l'émission pour mettre en forme le préambule, et permettent de diminuer la complexité de calcul des estimations du décalage en temps et/ou en fréquence du fait des relations de symétrie obtenues en sortie du module de transformation fréquence/temps du signal multiporteuse OQAM. Selon un premier exemple, pour des symboles de préambule pmn n définis en entrée de l'étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : pmn=~l2 pour m=4p et m=1+4p et pmn=-~l2 pour m=2+4p et m=3+4p avec pmn n un symbole de préambule associé à une sous-porteuse d'indice m à l'instant n, p un entier, et 0 m M -1, l'estimateur temporel met en oeuvre l'une des équations suivantes : zLSl = arg max l r(2)r(M / 2 -1+1-1 , z 2LS2 =argmaxlr(M/2+2)r(M-1+2)l. z Selon un deuxième exemple, pour des symboles de préambule pmn n définis en entrée de l'étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en 1 e:PLS2 _ - 7 émission par : P2m,n ="e l2 P2m+1,n = J~ / 2pour 0 2M, il est possible d'ajouter des colonnes de données nulles (formées de symboles de prémambule nuls aux instants n +1, n+ 2 , etc) pour améliorer les performances. L'ajout d'au moins colonne de symboles de préambule nuls est en revanche nécessaire pour les estimateurs temporels selon la deuxième équation (îLs2 )- Selon un troisième exemple, pour des symboles de préambule pm n définis en entrée de l'étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : pour 0 2LS1 = arg max lr(2)r(M / 2 -1+ 2)l , 2 2LS2 =argmaxlr(M/2+2)r(M-1+2)l, 2 zLS3 =argmax1r(2)r(M-l+2)+r(M/2-1+2)r(M/2+2)l. 2 Selon un cinquième exemple, pour des symboles de préambule p,n n définis en entrée de ladite étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : 1pm,n =1 si mod(m,4) = 0 pour 0mM-1 pm,n = 0 sinon l'estimateur temporel met en oeuvre l'une des équations suivantes : 2LSt =argmaxlr(2)r(M/2-1+2)+r(M/4-1+2)r(M/4+21 2 2LS2 = argmaxIr(M / 2+2)r(M-1+2)+r(3M / 4-1+2)r(3M / 4+2)~ 2 2LS3 = arg max lr(2)r(M -1+ 2) + r(M / 4 -1+ 2)r(3M / 4 + 2) 2 +r(M/4+2)r(3M/4-1+2)+r(M/2-1+2)r(M/2+2)~ Selon ces trois derniers exemples, on considère un système de transmission OQAM mettant en oeuvre un nombre de porteuses M multiple de 4 et un délai D tel que D = qM -1, ce qui permet d'obtenir d'autres relations particulières entre les sorties de l'étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission, et donc de simplifier encore les estimateurs. A nouveau, comme noté en relation avec les deux premiers exemples, il peut être nécessaire d'ajouter une ou plusieurs colonne de symboles de préambule nuls dans le préambule, avant la partie utile, tels que pm,n+i = 0 pour 0 m M -1 et i =1,2,.... L'ajout de ces symboles de préambule nuls est toujours facultatif pour les estimateurs temporels selon la première équation (îLsi ), mais nécessaire pour les estimateurs temporels selon la deuxième équation (zLs2 ) et la troisième équation (2Ls3 ). Dans un autre mode de réalisation, l'invention concerne un dispositif de réception d'un signal multiporteuse de type OQAM, comprenant des moyens d'estimation, dans le domaine temporel, d'un décalage en temps et/ou en phase du signal multiporteuse. Les moyens d'estimation comprennent au moins un estimateur pour l'estimation du décalage en temps, dit estimateur temporel, et/ou au moins un estimateur pour l'estimation du décalage en phase, dit estimateur de phase. Selon l'invention, le signal multiporteuse comprenant au moins un préambule, l'estimateur temporel et/ou l'estimateur de phase tient compte des coefficients d'un filtre prototype utilisé en émission pour mettre en forme le préambule. Un tel dispositif de réception est notamment adapté à mettre en oeuvre le procédé de réception décrit précédemment. Il s'agit par exemple d'un récepteur mettant en oeuvre un démodulateur OFDM/OQAM ou BFDM/OQAM. Un tel récepteur pourra bien sûr comporter les différentes caractéristiques relatives au procédé de réception selon l'invention, qui peuvent être combinées ou prises isolément. Ainsi, les caractéristiques et avantages de ce récepteur sont les mêmes que ceux du procédé de réception, et ne sont donc pas détaillés plus amplement. L'invention concerne également un programme d'ordinateur comportant des instructions pour la mise en oeuvre d'un procédé de réception tel que décrit précédemment lorsque le programme est exécuté par un processeur. En effet, le procédé selon l'invention peut être mis en oeuvre de diverses manières, notamment sous forme câblée ou sous forme logicielle. 4. 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 particulier, donné à titre de simple exemple illustratif et non limitatif, et des dessins annexés, parmi lesquels : - la figure 1, décrite en relation avec l'art antérieur, présente un schéma de modulation OFDM/OQAM classique ; - la figure 2 illustre un système de transmission mettant en oeuvre une technique de réception selon un mode de réalisation de l'invention ; - la figure 3 représente l'amplitude des sorties du module de transformation fréquence/temps pour un préambule donné ; - les figures 4 et 5 illustrent l'effet des filtres prototypes sur les valeurs pics obtenues en sorties du module de transformation fréquence/temps pour des préambules spécifiques ; - les figures 6 et 7 illustrent les fonctions de coût de deux exemples d'estimateurs selon l'invention ; - les figure 8 et 9 montrent les performances obtenues en utilisant des exemples d'estimateurs selon l'invention ; - la figure 10 présente la structure simplifiée d'un récepteur mettant en oeuvre une technique de réception selon un mode de réalisation de l'invention. 5. Description d'un mode de réalisation de l'invention 5.1 Principe général L'invention se place dans le contexte des systèmes de transmission mettant en oeuvre une modulation de type OFDM/OQAM ou BFDM/OQAM et mettant en oeuvre des préambules, et propose une technique permettant de simplifier l'opération de synchronisation des récepteurs. On note que ces systèmes de transmission pouvant être implémentés sous forme de bancs de filtres, ils sont également appelés FBMC/OQAM (en anglais «Filter Bank Multicarrier/OQAM »), ou plus simplement OQAM. Le principe général de l'invention repose sur l'utilisation, côté réception, d'un estimateur pour le décalage en phase et/ou en temps du signal multiporteuse, tenant compte du filtre prototype utilisé en émission pour mettre en forme le préambule. Un tel filtre prototype peut être très court et générer, à l'émission, un préambule localisé sur un unique symbole multiporteuse OQAM. En particulier, si l'on considère un canal mutli-trajets, les estimateurs proposés selon l'invention permettent d'estimer le retard maximal introduit par le canal, et donc de compenser correctement les effets du canal. 5.2 Exemple de mise en oeuvre 5.2.1 Système de transmission On présente, en relation avec la figure 2, un système de transmission mettant en oeuvre une technique de réception selon l'invention. Plus précisément, les opérations mises en oeuvre côté émission, à savoir les opérations de pré-traitement, conversion du domaine fréquentiel vers le domaine temporel et filtrage, du modulateur 10 sont classiques, et ont d'ores et déjà été décrites en relation avec la figure 1. Nous notons toutefois que dans le cas d'un système de transmission mettant en oeuvre un préambule, les symboles de données entrants à transmettre sur une sous-porteuse m à l'instant n peuvent être notés cmn au lieu de am n , et être à valeur complexe pour les symboles de préambule au lieu d'être à valeur réelle pour les symboles de données utiles. Par exemple, les symboles de données aux instants n = 0 et n =1 sont des symboles de préambule, notés pm n , et les symboles de données aux instants suivants sont des symboles de données utiles, notés am,n Le signal multiporteuse FBMC/OQAM en sortie du modulateur 10 peut donc s'exprimer sous la forme suivante : j m(k- s[k] cm nh[k - nN]e MD/2) _ m=0 avec cm,n = Pm,n des symboles de données à valeur complexe pour la mise en forme du préambule (préférentiellement à valeur réelle ou imaginaire pure), et cm,n = am,n des symboles de données à valeur réelle pour la mise en forme de la partie utile. A nouveau, on rappelle que des filtres prototypes distincts peuvent être utilisés pour mettre en forme le préambule et la partie utile du signal multiporteuse s [k] . Le passage dans le canal de transmission peut déformer le signal multiporteuse ski, en introduisant un bruit additif blanc gaussien, un décalage en temps 2Ts (également appelé retard), et/ou un décalage en phase 4 du signal multiporteuse. Le signal multiporteuse reçu peut donc s'exprimer sous la forme suivante : r[k] = s[k +b[k] avec b[k] un bruit blanc additif gaussien complexe circulaire centré de densité spectrale 62 . b Côté réception, on cherche à estimer ces paramètres de décalage en temps 2 et décalage en phase 4 au cours d'une étape d'estimation 21, puis à compenser le signal multiporteuse reçu au cours d'une étape de compensation 22, afin de synchroniser le récepteur. En sortie de l'étape de compensation 22, on obtient un signal synchronisé, qui peut s'exprimer sous la forme suivante : rsynchro [ k] = r [k + Île Ce signal synchronisé peut alors être démodulé de manière classique, lors d'une étape de démodulation 23, afin d'obtenir une estimation des symboles de données émis cm,n 5.2.2 Estimation des décalages en temps et en phase On décrit ci-après plus en détail l'étape d'estimation 21 des paramètres de décalage en temps 2 et décalage en phase 4 . Les inventeurs de la présente demande de brevet, également inventeurs de la demande de brevet français FR 1151590 déposée le 28 février 2011 au nom du même Demandeur, ont mis en évidence dans la demande de brevet précitée des relations particulières entre les différentes sorties des modules de transformation fréquence/temps 12 et de filtrage polyphase 13 mis en oeuvre en émission. Plus précisément, ils ont démontré que les sorties du module de transformation 11 M-1 n=+00 fréquence/temps sont conjuguées deux à deux, et les composantes polyphases du filtre prototype sont para-conjuguées deux à deux. Il est donc possible d'utiliser cette symétrie pour réutiliser une partie des résultats des multiplications survenant à des instants de filtrage différents, et réduire ainsi la complexité du filtrage. Il est également possible, selon l'invention en objet, d'utiliser ces relations pour réduire la complexité des estimateurs mis en oeuvre lors de l'étape d'estimation 21. Ainsi, en choisissant un délai tel que D = qM -1, avec q un entier, les inventeurs ont démontré qu'on obtient les relations suivantes en sortie module de transformation fréquence/temps 12, quels que soient les symboles de donnés entrants (données utile ou de préambule) : "k ,n = (-l)n "Al* /2-k-1,n "Al* uM/2+k,n = ~-l)n -k-1,n pour 0 /c M / 4 -1 pour 0 /c M / 4 -1 avec : um n un symbole transformé associé à la sortie d'indice m de l'étape de transformation fréquence/temps 12 à un instant n ; - * l'opérateur conjugué. Par conséquent, si on se place à l'instant n = 0 et quel que soit le type de préambule (aléatoire ou déterministe), on aura pour 0 k M / 4 -1 et un filtre prototype de longueur L = qM , avec q un entier tel que q 1 : h[M/2-1-k]s[k]=h[k]s*[M/2-1-k] et h[M-1-k]s[M/2+k]=h[M/2+k]s*[M-1-k] En l'absence de bruit, c'est-à-dire si l'on considère r[k]=s[k-z]ei(b, ces relations peuvent également s'écrire : h[M/2-1-k]r[k+2]=h[k]r*[M/2-1-k+T]ei20 (1) h[M-1-k]r[M/2+k+2]=h[M/2+k]r*[M-1-k+2]ei20 (2) pour 0/c M/4-1. carrés » fréquence/temps suivante Il est possible de définir deux estimateurs se basant sur une mesure de type « moindres (en anglais LS ou « Least Square ») entre les sorties du module de transformation 12 vérifiant une propriété de conjugaison. Le premier estimateur est basé sur la relation (1) et consiste à minimiser la fonction de coût : M/4-1 2 1/11(îM= ~h[M/2-1-k]r[k+î]-h[k]r [M/2-1-k+î]e~2~~ k=0 En développant cette équation, on obtient : M/4-1 M/4-1 t~l(î,¢)= h2[M/2-1-k]~r[k+î]~2+ h2[k]lr[M/2-1-k+-èr k=0 k=0 M/4-1 -29i e-2io h[k]h[M/2-1-k]r[k+î]r[M/2-1-k+î] k=0 ~1(î~~)= a(î)-a(î,~), avec : M/2-1 a(-)= h2[M /2-1-k]Ir[k+ Ir et k=0 M/4-1 e-2i(b 1 h[k]h[M/2-1-k]r[k+î]r[M/2-1-k+î] k=0 On peut constater que la quantité a(î) est indépendante de l'estimation du décalage en phase et faiblement liée à l'estimation du décalage temporel î . Par suite, le premier estimateur LS revient à maximiser conjointement la fonction : (zLS~LS) = argmaxP(î,¢) Pour obtenir l'estimation du décalage en phase Q~LS , on fixe l'estimation du décalage temporel î dans l'équation précédente, et on fait varier l'estimation du décalage en phase 0 , ce qui aboutit à l'estimateur de phase suivant : M/4-1 h(k)h(M / 2 -1- k)r(k + 2)r(M / 2-1- k + z) k=0 où L{.} désigne l'argument d'un nombre complexe. Par conséquent, en se basant sur les équations précédentes, on obtient l'estimateur suivant pour le décalage temporel : M/4-1 h(k)h(M/2-1-k)r(k+2)r(M/2-1-k+2) k=0 On note que pour cet estimateur temporel, le préambule utilisé peut être égal à un seul symbole multiporteuse FBMC/OQAM (de durée M/2) parce que les relations sur lesquelles se base 2LSt = arg max 2 cet estimateur (relations entre les premiers M/2 échantillons) restent valables et ne sont pas affectées par le symbole multiporteuse FBMC/OQAM suivant, qui est décalé de M/2 échantillons par rapport au préambule. Cette remarque s'applique même pour une longueur du filtre prototype supérieure à M. Le deuxième estimateur est basé sur la relation (2) et consiste à minimiser la fonction de coût suivante : M/4-1 2 IP-2(îM= ~h[M-1-k]r[M/2+k+î]-h[M/2+k]r [M-1-k+î]e~2~~ k=0 En suivant le même développement que ci-dessus, on obtient les estimateurs de phase et temporel suivants : {M/4-1 OLS2 = 2 L h(M -1- k)h(M / 2 + k)r(M / 2 + k + 2)r(M -1- k + 2) k=0 M/4-1 h(M-1-k)h(M/2+k)r(M/2+k+2)r(M-1-k+z) (4) k=0 On note que pour cet estimateur temporel, le préambule comprend de manière préférentielle un symbole multiporteuse FBMC/OQAM nul avant les données utiles. Le préambule utilisé peut donc être égal à deux symboles multiporteuses FBMC/OQAM. On constate que les deux estimateurs de phase et les deux estimateurs temporels proposés nécessitent une connaissance a priori sur le filtre prototype utilisé en émission, mais sur une plage de corrélation plus courte (M / 2 échantillons au lieu de NI). La complexité de calcul des estimateurs est donc diminuée de moitié. 5.2.3 Exemples de préambules et d'estimateurs temporels simplifiés Afin d'améliorer l'estimation de la synchronisation temporelle, on propose ci-après deux types de préambules et les estimateurs associés, permettant de diminuer encore la complexité de mise en oeuvre de la synchronisation. Ainsi, une simplification au niveau de calcul peut être faite pour les estimateurs proposés ci-dessus et pour des préambules spécifiques, aboutissant à des estimateurs plus performants. Cette simplification est due à la présence de valeurs pics à la sortie du module de transformation fréquence/temps, encore appelées valeurs crêtes ou valeurs maximales, pour le préambule. La position de ces valeurs pics dépend du préambule utilisé. La performance des estimateurs est liée à la position de ces valeurs pics par rapport aux coefficients du filtre prototype utilisé en émission. Un choix optimal de la position de ces valeurs pics est fait de sorte qu'elles soient filtrées par des coefficients de filtre assez importants. Pour ce faire, on propose différents préambules, permettant de remplir au moins l'un des t-L g2 = arg max 2 critères suivants : - la quantité de données contenues dans le préambule, et donc la longueur du préambule, doit être réduite au maximum, afin d'augmenter l'efficacité spectrale ; - la possibilité de réaliser de bonnes performances, i.e. la capacité à estimer précisément le retard 2 , qui dépend également de la conception de l'estimateur ; - la possibilité d'être utilisé pour d'autres fonctionnalités, comme l'estimation du canal de transmission. Dans les exemples ci-après, en se basant sur ces critères, les inventeurs ont proposé deux principaux types de préambule, avec pour chacun plusieurs variantes, et les estimateurs correspondants. Les choix de préambules et d'estimateurs sont évalués en considérant un filtre prototype optimisé pour la localisation temps-fréquence et de longueur L = M, noté TFL. Toutefois, l'invention ne se limite pas à l'utilisation d'un tel filtre prototype, et le principe proposé s'applique de manière équivalente à des filtres prototypes présentant une longueur arbitraire, notamment un multiple de M, avec M pair. A) Premier type de préambule On présente tout d'abord deux nouveaux préambules particulièrement intéressants, car ils peuvent être utilisés pour d'autres fonctions, comme l'estimation de canal. Selon un premier exemple, on considère un préambule émis à l'instant n présentant la structure suivante : pmn=~l2 pour m=4p et m=1+4p pmn=-~l2 pour m=2+4p et m=3+4p avec pmn un symbole de préambule associé à une sous-porteuse d'indice m à l'instant n, p un entier, et 0 m M - 1 . Selon un deuxième exemple, on considère un préambule émis à l'instant n présentant la structure suivante : P2m,n = ~ l 2 25 pour 0 30 multiporteuse, afin de conserver les propriétés de conjugaison. En utilisant l'estimateur temporel 'Tisi défini à l'équation (3), le préambule peut être formé d'une seule colonne de symboles de préambule (pour n = 0 par exemple), puisque les20 relations de conjugaison obtenues entre les premiers M / 2 échantillons du signal multiporteuse sont valables même si la colonne suivante (pour n =1 par exemple) appartient à la partie utile du signal multiporteuse. En revanche, si l'on utilise l'estimateur temporel zLS2 défini à l'équation (4), le préambule doit être formé d'au moins deux colonnes de symboles de préambule, dont une première définie par les équations précédentes avec p2m,n = / 2 et p2m+l,n = J~ / 2, et au moins une deuxième définie par pm,n+1 = 0 Un préambule selon le deuxième exemple permet d'obtenir des valeurs pics aux sorties d'indices M / 4 -1 et M / 4 de l'étape de transformation fréquence/temps plutôt qu'aux sorties d'indices 0 et M / 2 -1 selon le premier exemple. Ceci présente un intérêt, car si l'on considère un préambule selon le premier exemple, après filtrage par le coefficient h[0] du filtre prototype, la valeur pic de la sortie d'indice 0 va être atténuée, et par suite, la multiplication des deux valeurs pics aux sorties d'indices 0 et M / 2 -1 dans l'équation de l'estimateur temporel 'Tisi définie à l'équation (3) va être perturbée par le bruit et le canal multi-trajets. Il est donc souhaitable de proposer ce deuxième exemple de préambule, qui permet de décaler la position des valeurs pics. Les préambules selon ces deux premiers exemples présentent donc des paires de valeurs pics en sortie de l'étape de transformation fréquence/temps en émission à des positions différentes. Ainsi, le préambule selon le premier exemple présente des valeurs maximales aux sorties du module de transformation fréquence/temps 12 d'indices 0, M / 2 -1, M / 2 et M -1. Le préambule selon le deuxième exemple présente des valeurs maximales aux sorties du module de transformation fréquence/temps 12 d'indices M / 4 -1, M / 4 , 3M / 4 -1 et 3M / 4 . En utilisant le fait que des paires de valeurs pics peuvent être générées par le modulateur FBMC/OQAM à des positions particulières en choisissant un préambule approprié et en tenant compte des propriétés de symétrie/conjugaison intrinsèques au modulateur, il est possible simplifier encore l'opération d'estimation des décalages en phase et en temps du signal multiporteuse, en réduisant la complexité de calculs associés aux estimateurs. Plus précisément, on peut déduire que les valeurs pics du signal reçu contribuent majoritairement au calcul du retard en utilisant les estimateurs précédents définis aux équations (3) et (4). Ces valeurs pics sont les moins affectées en présence de bruit ou d'un canal multi-trajet. Par suite, on peut rendre ces estimateurs plus robustes et moins complexes en négligeant le calcul de produits des échantillons qui ne correspondent pas à des valeurs pics. On obtient alors des estimateurs LS simplifiés ou réduits. Par suite, ayant une connaissance a priori sur les valeurs pics du signal émis, on peut obtenir, pour chaque préambule, les estimateurs LS réduits correspondants aux estimateurs précédents. Ainsi, si l'on considère le préambule selon le premier exemple, le premier estimateur réduit correspondant à l'estimateur TLsi défini à l'équation (3) peut s'exprimer sous la forme suivante : zLsl = arg max l r(2)r(M / 2 -1 + 2)~ , puisque la première paire de valeurs pics est obtenue aux 2 sorties d'indices 0 et M / 2 -1, et le deuxième estimateur réduit correspondant à l'estimateur îLs2 défini à l'équation (4) peut s'exprimer sous la forme suivante : 1-Ls2 = arg max l r(M / 2+ 2)r(M -1+ 2)~ , puisque la deuxième paire de valeurs pics est 2 obtenue aux sorties d'indices M / 2 et M - 1 . Si l'on considère le préambule selon le deuxième exemple, le premier estimateur réduit correspondant à l'estimateur îLsi défini à l'équation (3) peut s'exprimer sous la forme suivante : 2LSi = arg maxlr(M / 4 -1 + 2)r(M / 4 + T)1, puisque la première paire de valeurs pics est 2 obtenue aux sorties d'indices M / 4 -1 et M / 4 ; et le deuxième estimateur réduit correspondant à l'estimateur îLs2 défini à l'équation (4) peut s'exprimer sous la forme suivante : 2Ls2 = arg max l r(3M / 4 -1 + 2)r( 3M / 4 + 2)~ , puisque la deuxième paire de valeurs pics est 2 obtenue aux sorties d'indices 3M / 4 -1 et 3M / 4 . B) Deuxième type de préambule On présente ci-après trois autres exemples de préambules formés d'une colonne de symboles de préambule dont la moitié des porteuses présente une valeur nulle, et éventuellement d'une ou plusieurs colonnes de symboles de préambule dont toutes les porteuses présentent une valeur nulle, dits symboles nuls. De tels préambules sont particulièrement intéressants du fait qu'ils permettent d'obtenir une propriété de conjugaison supplémentaire. Ainsi, si l'on se place dans le contexte d'un système de transmission FBMC/OQAM 25 mettant en oeuvre un nombre M de porteuse multiple de 4 et un délai D = qM -1 (avec D = L -1 dans le cas orthogonal et D L -1 dans le cas bi-orthogonal), et que l'on suit le raisonnement décrit dans la demande de brevet FR 1151590 susmentionnée, il est possible de montrer que les sorties de l'étape de transformation fréquence/temps appliquée symboles de préambule sont liées comme suit : uk n = (-l)n uM/2-k-1,n uM/2+k,n «-l)n uM-k-1,n uk n = (-l)n uM-k-1,n um n un symbole transformé associé à la sortie d'indice m de l'étape de transformation fréquence/temps 12 à un instant n ; - * l'opérateur conjugué. On considère par exemple que la première colonne de symboles de préambule est émise à un instant n pair. Tenant compte de cette nouvelle propriété de conjugaison, un nouvel estimateur peut être défini. En effet, on a pour 0 k M / 4 -1 la relation suivante, si le filtre prototype est 10 symétrique et de longueur L = M : s[k]=s*[M-1-k] pour 0kM/2-1. On note que si l'on considère un filtre prototype de longueur L = qM , il est souhaitable d'ajouter 2q -1 colonnes de symboles de préambule nuls en fin de préambule pour séparer le préambule (après modulation) des données utiles. La relation ci-dessus reste donc valable pour 15 0 On propose ci-après un cinquième exemple de préambule, permettant de limiter la valeur maximale des paires de pics obtenus en sortie de l'opération de transformation fréquence/temps, afin de limiter le PAPR (en anglais « Peak to Average Power Ratio », en français « facteur de crête ») du signal multiporteuse émis. Le préambule selon ce cinquième exemple, à l'instant n, présente par exemple la structure suivante : {19m,n = 1 19m,n = 0 si mod(m,4) = 0 pour 0mM-1. sinon Dans ce cas, les valeurs maximales sont obtenues sur les sorties d'indice 0, M / 4 -1, M / 4,M / 2 -1, M / 2 et M -1, comme illustré en figure 3 pour un nombre de porteuses M =128 , où l'axe des abscisses représente les indices des sorties du module de transformation fréquence/temps, et l'axe des ordonnées représente l'amplitude du préambule en sortie du module de transformation fréquence/temps. Un tel préambule permet donc de distribuer les valeurs maximales sur plus de sorties du module de transformation fréquence/temps, avec des amplitudes moins élevées, et donc de réduire le PAPR du signal émis. Comme déjà indiqué en relation avec les deux premiers exemples de préambule, il peut être nécessaire d'ajouter dans ces trois derniers exemples une ou plusieurs colonnes de symboles de préambule nuls dans le préambule après modulation, avant la partie utile, afin de conserver les propriétés de conjugaison après filtrage. Ainsi, si l'on utilise l'estimateur temporel zLs2 défini à l'équation (4) ou l'estimateur temporel 2Ls3 défini à l'équation (5), le préambule doit être formé d'au moins deux colonnes de symboles de préambule, dont une première définie par les équations précédentes, et au moins une deuxième définie par pm,n+1 = 0 pour 0 m M -1. Plus précisément, pour utiliser un filtre prototype de longueur L = qM , il est souhaitable d'insérer 2q -1 colonnes de symboles de préambules nuls dans le préambule, pour séparer la première colonne de symboles de préambule définie par les équations précédentes et les données utiles. En revanche, pour utiliser l'estimateur temporel 2LSt défini à l'équation (3), le préambule peut n'être formé que d'une seule colonne de symboles de préambule (pour n = 0 par exemple). Comme décrit précédemment en relation avec les deux premiers exemples de préambule, il est également possible de réduire/simplifier les estimateurs proposés et d'augmenter leurs performances en tenant compte de la position des paires de valeurs pics générées par le modulateur FBMC/OQAM. Ainsi, si l'on considère le préambule selon le troisième exemple, le premier estimateur réduit correspondant à l'estimateur 2Ls3 défini à l'équation (3) peut s'exprimer sous la forme suivante : zLSl = arg maxlr(M / 4 -1+ 2)r(M / 4 + 2)I , puisque la première paire de valeurs pics est 2 obtenue aux sorties d'indices M / 4 -1 et M / 4 , le deuxième estimateur réduit correspondant à l'estimateur 2Ls2 défini à l'équation (4) peut s'exprimer sous la forme suivante : zLs2 = arg max l r(3M / 4 -1 + 'Or( 3M / 4 + z)l , puisque la deuxième paire de valeurs pics est z obtenue aux sorties d'indices 3M / 4 -1 et 3M / 4 , et le troisième estimateur réduit correspondant à l'estimateur 2Ls3 défini à l'équation (5) peut s'exprimer sous la forme suivante : zLs3 = arg max I r(M / 4 -1 + 2)r( 3M / 4 + 2) + r(M / 4 + 2)r(3M / 4 -1 + 2)l , puisque les z deux paires de pics sont obtenues aux sorties d'indices M / 4 -1, M / 4 , 3M / 4 -1 et 3M / 4 . Si l'on considère le préambule selon le quatrième exemple, le premier estimateur réduit correspondant à l'estimateur rLsi défini à l'équation (3) peut s'exprimer sous la forme suivante : 20 2LS1 = arg max lr(2)r(M / 2 -1 + 2)l , puisque la première paire de valeurs pics est obtenue aux 2 sorties d'indices 0 et M / 2 -1 ; le deuxième estimateur réduit correspondant à l'estimateur C^LS2 défini à l'équation (4) peut s'exprimer sous la forme suivante : 2LS2 = arg max lr(M / 2+ 2)r(M -1+ 2)~ , puisque la deuxième paire de valeurs pics est 2 obtenue aux sorties d'indices M / 2 et M -1, et le troisième estimateur réduit correspondant à l'estimateur 2LS3 défini à l'équation (5) peut s'exprimer sous la forme suivante : 2LS3 = arg max l-1 + 2) + r(M / 2 -1 + 't-)r(M / 2+2)1 , puisque les deux paires de 2 pics sont obtenues aux sorties d'indices 0, M / 2 -1, M / 2 et M - 1 . Si l'on considère le préambule selon le cinquième exemple, le premier estimateur réduit correspondant à l'estimateur iLSi défini à l'équation (3) peut s'exprimer sous la forme suivante : 2LSi = arg max ~r(2)r(M / 2 -1+r) + r(M / 4 -1 + 2)r(M /4+2)1 puisque les valeurs pics 2 sont obtenues aux sorties d'indices 0, M / 2 -1, M / 4 -1 et M / 4 ; le deuxième estimateur réduit correspondant à l'estimateur 2LS2 défini à l'équation (4) peut s'exprimer sous la forme suivante : 2LS2 = arg rnax Ir(M / 2 + 2)r(M -1 + 2) + r(3M / 4 -1 + 2)r(3M / 4 + 2)~ puisque les 2 valeurs pics sont obtenues aux sorties d'indices M / 2 , M -1, 3M / 4 -1, 3M / 4 , et le troisième estimateur réduit correspondant à l'estimateur 2LS3 défini à l'équation (5) peut s'exprimer sous la forme suivante : 2LS3 = arg nax l r(2)r(M -1 + 2) + r(M / 4 -1 + 2)r(3M / 4 + 2) 2 +r(M /4+ 2)r(3M /4-1+ 2) + r(M /2-1+ 2)r(M / 2+ 2)l On constate donc que les estimateurs simplifiés/réduits proposés selon l'invention présentent de bonnes performances par rapport aux estimateurs non simplifiés. En effet, en présence de bruit et d'un canal multi-trajet, les échantillons du signal multiporteuse à faibles amplitudes vont être perturbés et par suite diminuer la performance de l'estimateur. On propose donc selon au moins un mode de réalisation de l'invention d'éliminer ces échantillons et de ne conserver que les valeurs pics du signal multiporteuse, qui sont elles moins affectées par le bruit et le canal. Ceci conduit à l'obtention d'estimateurs réduits qui présentent de meilleures performances que les estimateurs non simplifiés, tels que définis aux équations (3), (4) et (5), car ces estimateurs simplifiés permettent de réduire le nombre des valeurs crêtes autour de la valeur à estimer. 5.3 Génération des préambules Bien entendu, d'autres préambules pourraient être proposés, permettant d'obtenir une ou plusieurs paires de valeurs pics en sortie du module de transformation fréquence/temps à des positions spécifiques. De tels préambules présentent de bonnes performances lorsqu'ils sont utilisés par un récepteur selon l'invention, mais pourraient également être utilisés par un récepteur classique quelconque. On propose ci-après deux exemples d'algorithmes permettant de générer de tels préambules. On suppose pour ce faire que le nombre de porteuses M est une puissance de 2 telle que M 8, que la longueur du filtre prototype est L = M et que, de même, le nombre de symboles de préambule p,n n non nul est limité à M, soit une seule colonne de symboles de préambule à un instant donné. Comme déjà indiqué, si un filtre prototype plus long est utilisé, il suffit d'ajouter des colonnes de symboles de préambule nuls dans le préambule pour le séparer des données utiles. Sous ces hypothèses, seuls les symboles de préambules p,n o correspondant à la colonne de symboles de préambule à l'instant n = 0 , également appelés pilotes, sont à prendre en compte, et le signal en sortie du module de transformation fréquence/temps peut s'écrire, en ne considérant en entrée et sortie que des séquences discrètes et en s'affranchissant de l'indice temporel n, sous la forme suivante : M-1 j2tcm(k-(M-1)/2) j-m uk = pme m e 2 m=o qui après simplification donne : M-1 tc tcm 2tcmk M-1 rn uk pme 2 e~ M ej m = p m=o m=o 7r 7r m mk avec On k = - 2 m + M + 27r . La détermination des préambules peut être mise en oeuvre en : - recherchant les valeurs des pilotes produisant une ou deux paires de valeurs pics uk à valeurs maximales (en sortie du module de transformation fréquence/temps) ; - supposant, dans un premier temps que les pilotes sont soit réels purs, soit imaginaires purs. Du fait des relations de symétrie en sortie du module de transformation fréquence/temps précédemment exposées, i.e. uM/2-1-k = uk et uM_1-k = uM/2+k pour 0 k M / 4 -1, les pics se produisent toujours par paire. Par ailleurs, il est souhaitable que ces paires de pics apparaissent en nombre limité pour éviter que leur énergie soit dispersée sur un nombre trop important de sorties uk du module de transformation fréquence/temps. Par exemple, si l'on cherche à obtenir une paire de valeurs pics d'amplitude maximale en ko et (M / 2 -1- ko), il faut déterminer les valeurs de pilotes pm donnant le maximum de l'expression maxpm luko l2 . Sachant que L l uk l2 - l uko l2 + 1 k,kk0l uk l2 , le choix de la fonction objective devrait en même temps garantir qu'en dehors des positions des valeurs pics, les valeurs des sorties uk du module de transformation fréquence/temps présentent une faible énergie. Les deux algorithmes proposés ci-après permettent de déterminer des préambules purement réels produisant respectivement une ou deux paires de valeurs pics. Ils peuvent bien entendu s'étendre à la génération de préambules purement imaginaires ou de préambules mixtes, combinant des symboles de préambule purement réels et purement imaginaires, ou même complexes. 5.3.1 Détermination d'une paire unique de valeurs pics Comme indiqué ci-dessus, obtenir une valeur crête (pic) à la sortie d'indice ko du module de transformation fréquence/temps revient à maximiser l'expression l uko l Or d'après l'équation simplifiée de la valeur uk proposée ci-dessus, on constate que toutes les phases d'indices m et M - m , pour 1 m M / 2 -1, sont symétriques par rapport à l'axe imaginaire. En effet, les quantités ejeo'k0 et e1em'2'k° sont réelle et imaginaire pures respectivement, et sont égales à 1 et j(-1)k0 en supposant que M = 2q avec q 3 entier. En effet, le terme de la phase d'indice M - m est tel que : - z(M-m)+27r (M - m)ko BM-rn,ko --2(M - m) + M M BM-rn,ko = 2m+~- M -M mk° = Par suite on peut écrire l'équation simplifiée de la valeur uk de la manière suivante : M/2-1 -j-rn j~cyn j2tcink0 M/2-1 j-m -j2irmko uko = po +j(-1) pM/2 + prne k 2 e M e M - pM-,tee 2 e m m=1 m=1 k uko = po + j(-1)k pM/2 + { Jem, ko -Jem, ko l - pM_rne 1 M/2-1 rn=1 Si l'on suppose que les symboles de préambules pm ont des valeurs provenant d'une constellation issue d'une modulation de phase à deux états (en anglais BPSK pour «Binary phase-shift keying »), le problème revient à déterminer les signes des symboles de préambules pm . Or, les couples de symboles de préambule pm et pm_m qui contribuent à la somme dans l'équation ci-dessus ont des arguments Om k et (-0m ,k +7L) qui les positionnent de manière symétrique par rapport à l'axe imaginaire. Donc en fonction de leurs signes respectifs, le vecteur résultant associé à chacun de ces symboles de préambule est soit réel pur, soit imaginaire pur. La procédure de maximisation peut donc être mise en oeuvre par l'algorithme suivant : 1) préciser la valeur désirée de la position ko pour la valeur pic ; 2) choisir les signes des symboles de préambule po et pM/2 ; 3) calculer, pour chaque indice m tel que 1 m M / 2 -1, la valeur de Om ko et comparer sa projection sur l'axe des réels (cosinus) et sa projection sur l'axe des imaginaires (sinus), ce qui revient à tester la condition suivante : si 1cos(Om,ko )1 >- 1sin(0m ko alors bool(m) = 0 sinon bool(m) = 1, où l'opérateur bool représente une variable à deux états ; 4) si bool(m) = 0 , affecter des signes opposés aux symboles de préambule pm et Pm_m, c'est-à-dire pm = -Pm-m Ce signe est choisi tel que Sign(pm) * Sign(cos (0m ko» = Sign(po) . Si bool(m) =1, affecter le même signe aux symboles de préambule pm et pm_m, c'est-à-dire pm = pm_m . Ce signe est choisi tel que Sign(pm) * Sign(sin (0m ko» = Sign«-1)k0 pm/2). A titre d'exemple, on suppose que l'on cherche à obtenir une valeur pic à la position d'indice ko = M / 4 . On considère des symboles de préambule tels que po = pm 2 =1. La phase associée est telle que 0m ko = zm / M , ce qui implique que Icos(0m ko )l Isin(0m ko )l pour 1 conduisant à l'obtention de deux paires de valeurs pics en sortie du module de transformation fréquence/temps. On suppose ici que les symboles de préambule présentent une valeur telle que pm = ±1, et on cherche à déterminer deux paires de valeurs pics séparées de M/2 échantillons. Les deux premières valeurs pics à maximiser sont alors localisées aux positions k0 et M / 2 + k0 , et les deux autres valeurs se déduisent à partir des relations de symétrie. La sortie du module de transformation fréquence/temps correspondant à l'indice M / 2 + k0 peut s'exprimer sous la forme suivante : M-1 zrn 27rm(M/2+4)) M_t zmn 2zmko -Jrn J J m -J- -J J "Al /2+ko - terne 2 e M e m -(-1) terne 2 e m e m m=0 m=0 En comparant cette équation et l'équation simplifiée de la valeur uk proposée en relation avec le premier algorithme, on peut remarquer que l'optimisation (c'est-à-dire la maximisation) des deux sorties d'indices kp et M / 2 + k0 aboutit à une contradiction au niveau des signes. En effet, si l'on fait des choix d'optimisation identiques pour uko et uM/2+ko (ce qui revient à fixer les signes de pp et pM /2 ), on peut remarquer que l'optimisation des indices pairs de m aboutit à la même distribution des signes pour uko et uM/2+ko . Par contre, l'optimisation des indices impairs aboutit à deux choix qui ont des signes opposés. Par suite, si l'on gagne sur l'un des vecteurs résultant de uko , on perd sur le même vecteur résultant de uM /2+ko et vice versa. Pour obtenir des valeurs pics de même amplitude pour les valeurs uko et uM/2+ko , il faut donc conserver l'optimisation des indices pairs et choisir correctement les vecteurs résultants des indices impairs de uko et uM /2+ko , de façon à maximiser les résultantes finales pour lesquelles on impose d'avoir la même amplitude. Pour ce faire, on peut tout d'abord remarquer que la phase de l'indice m (pair ou impair) et celle de l'indice M / 2+ m sont décalées de 7L / 2 . En effet, on a : B - 7s M/2+m +7z(M/2+m)+2~(M/2+m)k0 M/2+m,ko - - 2 ( ) M M 9M / 2+m, kp = - Ir Ir 7CM m + + + 7Lk0 + 27Lk = + zk0 + 9m, k Donc la différence de phase entre BM/2+m ko et Bm ko est ±7L / 2, avec un signe qui dépend de la parité de kp . Ceci implique que si l'on maximise pm pour optimiser les sorties uko d'indices pairs, le vecteur résultant imaginaire est égal, en norme, à la résultante réelle. Il faut désormais choisir la manière optimale pour ajouter les vecteurs résultants des indices impairs, sachant que les vecteurs résultants de uko sont opposés à ceux résultants de uM /2+ko , comme indiqué précédemment. Dans les deux cas, on peut dire que si l'optimisation du couple (m, M - m) a une projection maximale sur l'axe des réels, l'optimisation du couple (M / 2+ m, M / 2 - m) a, nécessairement, une projection maximale sur l'axe des imaginaires et vice versa. Or on rappelle que les vecteurs résultants de ces deux couples ont la même amplitude. La solution optimale est donc d'obtenir un vecteur résultant des indices impairs qui est perpendiculaire à celui des indices pairs. Ceci peut être réalisé en choisissant, par exemple, le couple (pm, pM_m) qui maximise le vecteur résultant réel de uko et le couple (pM/2+m, pM/2-m) qui maximise le vecteur résultant imaginaire de uM/2+ko La procédure d'optimisation peut donc être mise en oeuvre par l'algorithme suivant : 1) préciser la valeur de la position ko pour la valeur pic et déduire la valeur de la position M / 2 + ko ; 2) choisir les signes des symboles de préambule po et pm /2 ; 3) calculer, pour chaque indice m tel que 1 m M / 2 -1, la valeur de 0m ko et comparer sa projection sur l'axe des réels (cosinus) et sa projection sur l'axe des imaginaires (sinus), ce qui revient à tester la condition suivante : si Icos(0m ko >- Isin(0m ko )1 alors bool(m) = 0 sinon bool(m) = 1, où l'opérateur bool représente une variable à deux états ; 4) si bool(m) = 0 : si m est impair : affecter des signes opposés aux symboles de préambule pm et pm_m - Ce signe est choisi tel que Sign(pm) * Sign (cos (0m ko» _ -Sign(po) dans un premier cas, ou Sign ( pm ) * Sign (cos (0m ko ) ) = Sign( po ) dans un deuxième cas ; si m est pair : affecter des signes opposés aux symboles de préambule pm et pM _ m , tels que Sign ( pm ) * Sign (cos (0m ko» = Sign( po ) ; si bool(m) =1 : affecter des signes opposés aux symboles de préambule pm et pM_m . Ce signe est choisi tel que Sign(pm) * Sign(sin(0m ko» = Sign ((-1)k0 pm/2) dans le premier cas, ou Sign (pm)*Sign (sin (0mko»=-Sign ((-l)k°pM/2) dans le deuxième cas. On note que plusieurs choix sont possibles sur les directions finales des résultantes impaires. A titre d'exemple, on suppose que l'on cherche à obtenir une valeur pic à la position d'indice ko = M / 4 et en respectant les différentes étapes de l'algorithme proposé ci-dessus, on obtient le même préambule que pour le premier algorithme, en inversant les signes des indices impairs pour 1 m M / 4 -1 et pour 3M / 4 + 1 m M -1. A nouveau, comme pour le premier algorithme, si l'on considère que les symboles de préambule peuvent également prendre des valeurs imaginaires, il est possible d'obtenir d'autres préambules encore plus performants, du fait de la présence de degrés de liberté supplémentaires. 5.4 Effet du filtre prototype Si l'on considère à nouveau le schéma du modulateur FBMC/OQAM des figures 1 et 2, on constate que les symboles de données obtenus en sortie du module de pré-traitement 11, et entrant dans le module de transformation fréquence/temps 12, dépendent du délai D introduit par le filtre prototype. Or pour avoir les trois relations de conjugaison présentées en relation avec le deuxième type de préambule, à savoir : uk,n «-l)n uM/2-k-1,n pour 0 k M / 4 -1 uM/2+k,n «-l)n uM-k-1,n pour 0 k M / 4 -1 uk,n = (-l)n uM-k-1,n pour 0 k M / 2 -1 il faut que le délai D soit tel que D = qM -1, avec q un entier supérieur ou égal à 1 (ce qui correspond à L = qM pour le cas OFDM/OQAM). De plus, pour pouvoir utiliser le troisième estimateur 2Ls3 défini à l'équation (5), il faut que la longueur de préambule soit égale à la longueur de filtre prototype. Par suite et à titre d'exemple, pour une longueur L = qM , le préambule doit être composé de 2q colonnes de symboles de préambule (correspondant aux instants n, n +1, ..., n+ 2q -1). Comme décrit précédemment, la première colonne de symboles de préambule comprend des pilotes spécifiques, et les autres 2q -1 colonnes de symboles de préambule des valeurs nulles. Les inventeurs de la présente demande de brevet ont constaté que le filtre prototype utilisé en émission influe directement sur la performance des estimateurs. En effet, pour conserver une amplitude importante des valeurs pics obtenues en sortie du module de transformation fréquence/temps, ces valeurs doivent être filtrées (i.e. multipliées) par des coefficients de filtrage qui ne sont pas négligeables par rapport aux autres coefficients. Ainsi, si le pic obtenu doit être filtré par un coefficient de faible valeur, il est souhaitable d'augmenter la taille du préambule pour décaler la position de ce pic. Il est bien entendu possible d'utiliser le premier estimateur îLSt défini à l'équation (3) et le deuxième estimateur îLS2 défini à l'équation (4) pour des filtres prototypes de longueur L = qM en utilisant les préambules définis précédemment, éventuellement en ajoutant une seule colonne de symboles de préambule nuls. Toutefois, il est bien connu que les coefficients les plus significatifs de la plupart des filtres prototypes localisés en temps et en fréquence sont distribués sur les M/2 coefficients qui sont autour du centre de ces filtres prototypes. Par suite, en utilisant à titre d'exemple les préambules selon les deuxième ou troisième exemples ci-dessus et les estimateurs réduits correspondants, on peut profiter de l'une des deux paires de valeurs pics obtenues aux positions {M / 4 -1,M / 41 ou 13M / 4 -1,3M / 41 en sortie du module de transformation fréquence/temps, en ajoutant q -1 colonnes de symboles de préambule nuls au préambule et en filtrant le préambule par les coefficients h[qM/2-M/4-1] et h[qM/2-M/4]. A titre d'exemple, la figure 4 illustre les valeurs pics obtenues aux positions {M / 4 -1,M / 41 et 13M / 4 -1,3M / 41 en sortie du module de transformation fréquence/temps, avant filtrage, en utilisant un préambule Pn émis à l'instant n selon le troisième exemple décrit ci-dessus, à savoir : 1P2m,n = l 2 P2m+1,n = 0 On considère un filtre prototype 40, illustré en pointillé, présentant une longueur L = M . Les valeurs pics vont donc être filtrées par les coefficients {h[M / 4 -1],h[M / 4],h[3M / 4 -1],h[3M / 4]} du filtre prototype, qui ont des amplitudes assez importantes pour conserver des valeurs pics une fois le signal émis. On note que les M / 2 premiers échantillons du signal multiporteuse émis (correspondant au préambule dont la première paire de valeurs crêtes fait partie) préservent une relation de conjugaison (à un coefficient près) même si le symbole multiporteuse OQAM suivant émis à l'instant n + 1, noté Sn + 1, n'est pas nul. Par suite, le premier estimateur îLSt défini à l'équation (3) peut être utilisé. En revanche, les M/2 échantillons suivant perdent les relations de conjugaison si le symbole Sn+ 1 n'est pas nul, et le deuxième estimateur 2LS2 défini à l'équation (4) ne peut pas être appliqué. La figure 5 illustre un autre exemple des valeurs pics obtenues aux positions en sortie du module de transformation fréquence/temps, avant filtrage, en utilisant un préambule Pn émis à pour 0M/2-1 l'instant n selon le troisième exemple décrit ci-dessus, et considérant un filtre prototype 50, illustré en pointillé, présentant une longueur L = qM = 2M . On peut remarquer que pour les M/2 premiers échantillons du signal multiporteuse émis (correspondant au préambule dont la première paire de valeurs crêtes correspondant aux positions 1M / 4 -1,M / 41 fait partie), la première paire de valeurs pics va être affaiblie à cause du filtrage par des coefficients à faibles amplitudes. Il est donc préférable d'utiliser la deuxième paire de valeurs pics (correspondant aux positions 13M / 4 -1,3M / 41), qui est filtrée par un couple de coefficients d'amplitudes importantes {h[3M / 4 -11,h[3M / 4]}, pour estimer le décalage temporel. En particulier, pour que la deuxième paire de valeurs pics ne soit pas perturbée par les symboles multiporteuses OQAM émis aux instants suivants, le symbole multiporteuse OQAM émis à l'instant n + 1, noté Sn + 1, est de préférence nul. On note toutefois que dans le cas illustré en figure 5, même si le symbole Sn+1 n'est pas nul, la perturbation engendrée par les échantillons de ce symbole Sn+1 n'est pas trop gênante car ces échantillons, qui vont s'ajouter aux valeurs pics, sont filtrés par des coefficients faibles. Plus généralement, on constate que la première paire de valeurs pics est préférentiellement utilisée dans le calcul de l'estimateur réduit pour q impair, et que la deuxième paire de valeurs pics est préférentiellement utilisée dans le calcul de l'estimateur réduit pour q pair. Par suite, il est possible de généraliser la réduction du premier estimateur défini à l'équation (3) et du deuxième estimateur défini à l'équation (4), quelque soit q, entier supérieur ou égal à 1, par l'équation suivante : zLS4 =arginaxlr(qM/2-M/4-l+z)r(gM/2-M/4+2)l (6). z La figure 6 illustre la fonction de coût V4 de cet estimateur en utilisant un filtre prototype de type TFL et de longueur L = 4M , pour un rapport signal à bruit de 15dB, un nombre de porteuses M =128 , un retard 2 / Ts =10 et un préambule selon le troisième exemple décrit ci-dessus. On constate sur cette figure 6 la présence de deux pics distants de M / 2 échantillons. Le premier pic est dû à la multiplication des deux valeurs maximales obtenues en sortie du module de transformée fréquence/temps localisées aux positions 13M / 4 -1,3M / 41, et filtrées par les coefficients du filtre prototype {h[qM / 2 - M / 4 -11,h[gM / 2 - M / 411. Le deuxième pic est dû à la multiplication des deux valeurs maximales obtenues en sortie du module de transformée fréquence/temps localisées aux positions 1M / 4 -1,M / 41, et filtrées par les coefficients du filtre prototype {h[qM / 2+ M / 4 -11,h[gM / 2+ M / 411, ce qui est équivalent à 2+ M / 2 dans le quatrième estimateur 2LS4 défini à l'équation (6). On note que la relation h[qM/2+M/4]s[gM/2+M/4-1] = h[qM/2+M/4-1]s*[qM/2+M/4] reste valable même si l'on ajoute q -1 colonnes de symboles de préambule nuls dans le préambule, car la qième colonne n'aura pas d'effet sur les valeurs pics {M / 4 -1,M / 41 de préambule après filtrage, ces valeurs pics étant situées dans la première moitié des sorties du module de transformation fréquence/temps. Il est également possible de profiter des deux paires de valeurs pics obtenues en sortie du module de transformation fréquence/temps aux positions {M -1,01 ou {M / 2 -1,M / 21, en utilisant à titre d'exemple le préambule selon le quatrième exemple ci-dessus et les estimateurs réduits correspondants Dans ce cas, les deux paires de valeurs pics obtenues aux positions précitées doivent être filtrées par des coefficients de filtrage {h[qM / 2 - 11,h[gM / 211, et il faut ajouter q colonnes de symboles de préambuleM nuls dans le préambule. Plus généralement, on constate que la paire de valeurs pics correspondant aux positions {M -1,0} est préférentiellement utilisée dans le calcul de l'estimateur réduit pour q pair, et que la paire de valeurs pics correspondant aux positions {M / 2 -1,M / 21 est préférentiellement utilisée dans le calcul de l'estimateur réduit pour q impair. Par suite, il est possible de généraliser la réduction du premier estimateur défini à l'équation (3) et du deuxième estimateur défini à l'équation (4), quelque soit q, entier supérieur ou égal à 1, par l'équation suivante : 2LS5 = arg inax Ir(gM / 2 -1 + r)r(gM / 2+ r)l (7) r La figure 7 illustre la fonction de coût vis de cet estimateur en utilisant un filtre prototype de type TFL et de longueur L = 4M , pour un rapport signal à bruit de 15dB, un nombre de porteuses M =128 , un retard r / Ts =10 et un préambule selon le quatrième exemple décrit ci- dessus. Si l'on ajoute q -1 colonnes de symboles de préambule nuls dans le préambule, on obtient pratiquement les mêmes performances en termes de coût, car le symbole multiporteuse OQAM utile suivant le préambule, correspondant à la colonne de symboles de données utiles à l'instant q, qui va s'ajouter au préambule, sera filtré par les Mpremiers coefficients de filtre prototype, qui ont généralement une faible amplitude par rapport aux Mcoefficients suivants. Les performances des estimateurs selon l'invention sont présentées en relation avec les figures 8 et 9. Plus précisément, la figure 8 propose une comparaison entre le premier estimateur îLSt défini à l'équation (3) et le troisième estimateur 2Ls3 défini à l'équation (5), avant et après simplification en utilisant les préambules selon le deuxième type de préambule présenté. Le premier estimateur îLSi réduit en utilisant le préambule selon le quatrième exemple est noté zLSl-R - Le troisième estimateur 2Ls3 réduit en utilisant le préambule selon le troisième exemple est noté 2Ls3-R Les courbes obtenues illustrent la performance des estimateurs après simplification en terme d'erreur quadratique moyenne RMSE (en anglais «Root Mean Square Error») par rapport au rapport signal à bruit SNR (en anglais « Signal to Noise »). La performance des estimateurs réduits est notamment due à l'amélioration de leur robustesse contre l'effet d'un canal multi-trajets et du bruit gaussien, en éliminant les échantillons les plus perturbés. L'erreur quadratique moyenne s'exprime ici sous la forme suivante : 1 1 re RMSE (z/T) =- -~(1- -îi)2 M re i=1 où re est le nombre de réalisations, îi la valeur estimée du retard 2 à la ième itération et T le temps symbole (avec T = M en discret). Selon la figure 8, le premier estimateur 2LS1_R réduit en utilisant le préambule selon le quatrième exemple est le plus performant. La figure 9 illustre les performances des mêmes estimateurs, en terme de RMSE en fonction de nombre de porteuses, pour un rapport signal à bruit de 10 dB et en présence d'un canal de transmission de type Ma. Les estimateurs montrent toujours un gain après réduction. Le premier estimateur 2LS1_R réduit en utilisant le préambule selon le quatrième exemple est toujours le plus performant, du fait que les échantillons qui entrent dans le calcul de l'estimateur ont les amplitudes les plus importantes après filtrage, et, par suite, sont les moins perturbés lors de la transmission. Au niveau de la complexité, les premier et deuxième estimateurs définis par les équations (3) et (4), qui utilisent la corrélation entre la première ou la deuxième moitié des échantillons du préambule après filtrage sont moins complexes, en termes de calcul opératoire, que le troisième estimateur défini par l'équation (5) utilisant la corrélation entre tous les échantillons de préambule. Après simplification de ces estimateurs, en exploitant la présence des valeurs pics après l'étape de transformation fréquence/temps, la complexité de ces estimateurs devient négligeable. 5.5 Structure du récepteur On présente finalement, en relation avec la figure 10, la structure simplifiée d'un récepteur comprenant un démodulateur OFDM/OQAM ou BFDM/OQAM mettant en oeuvre une technique de réception selon un exemple décrit ci-dessus. Un tel récepteur comprend une mémoire 101 comprenant une mémoire tampon, une unité de traitement 102, équipée par exemple d'un microprocesseur µP , et pilotée par le programme d'ordinateur 103, mettant en oeuvre le procédé de réception selon l'invention. A l'initialisation, les instructions de code du programme d'ordinateur 103 sont par exemple chargées dans une mémoire RAM avant d'être exécutées par le processeur de l'unité de traitement 102. L'unité de traitement 102 reçoit en entrée un signal multiporteuse r [k] . Le microprocesseur de l'unité de traitement 102 met en oeuvre les étapes du procédé de réception décrit précédemment, selon les instructions du programme d'ordinateur 103, pour effectuer une estimation des retards en temps et en phase affectant le signal multiporteuse reçu r [k] . Pour cela, le récepteur comprend, outre la mémoire tampon 101, des moyens d'estimation comprenant au moins un estimateur estimateur temporel, et/ou au moins un estimateur de phase tenant compte des coefficients d'un filtre prototype utilisé en émission pour mettre en forme le préambule. Ces moyens sont pilotés par le microprocesseur de l'unité de traitement 102. REVENDICATIONS1. Procédé de réception d'un signal multiporteuse de type OQAM mettant en oeuvre une étape d'estimation (21), dans le domaine temporel, d'un décalage en temps et/ou en phase dudit signal multiporteuse, ladite étape d'estimation (21) met en oeuvre au moins un estimateur pour l'estimation dudit décalage en temps, dit estimateur temporel, et/ou au moins un estimateur pour l'estimation dudit décalage en phase, dit estimateur de phase, caractérisé en ce que, ledit signal multiporteuse comprenant au moins un préambule, ledit estimateur temporel et/ou ledit estimateur de phase tient compte des coefficients d'un filtre prototype utilisé en émission pour mettre en forme ledit préambule. 2. Procédé de réception selon la revendication 1, caractérisé en ce que ladite étape d'estimation (21) tient également compte de la position d'au moins une paire de valeurs maximales obtenues en émission pour ledit signal multiporteuse, ladite position étant définie par rapport aux sorties d'une étape de transformation du domaine fréquentiel vers le domaine temporel (12) d'un ensemble de symboles de données destinés à former ledit préambule, dits symboles de préambule. 3. Procédé de réception selon la revendication 1, caractérisé en ce que ledit estimateur temporel met en oeuvre l'une des équations suivantes : M/4-1 h(k)h(M/2-1-k)r(k+2)r(M/2-1-k+T) k=0 zLSf = arg max T et M/4-1 h(M-1-k)h(M/2+k)r(M/2+k+2)r(M-1-k+T) k=0 et/ou ledit estimateur de phase met en oeuvre l'une des équations suivantes : {M/4-1 h(k)h(M /2-1-k)r(k+2)r(M /2-1-k+T) et k=0 {M/4-1 h(M-1-k)h(M/2+k)r(M/2+k+T)r(M-1-k+T) k=0 avec : 25 TLSt ou zLs2 l'estimation dudit décalage en temps en sortie dudit estimateur temporel ; ~LSt ou OLs2 l'estimation dudit décalage en phase en sortie dudit estimateur de phase ; LH l'argument d'un nombre complexe ; M le nombre de sous-porteuses d'un symbole OQAM dudit signal multiporteuse ; k un entier tel que 0 k M / 4 -1 ; TLs2 = arg max T OLst = 1 ~LS2 _ - 1h(k) les coefficients dudit filtre prototype utilisé en émission ; r(k) le signal multiporteuse reçu, tel que r(k) = s(k - 2)ei(b + b(k) ; s(k) le signal multiporteuse émis ; b(k) un bruit blanc gaussien ; 2 le décalage en temps du signal multiporteuse reçu par rapport au signal multiporteuse émis ; 4 le décalage en phase du signal multiporteuse reçu par rapport au signal multiporteuse émis. 4. Procédé de réception selon la revendication 2, caractérisé en ce que, pour des symboles de préambule pmn n , définis en entrée de ladite étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : pm,n =~12 pour m=4p et m=1+4p et pmn= -NE l 2 pour m=2+4p et m=3+4p avec pmn n un symbole de préambule associé à une sous-porteuse d'indice m à l'instant n, p un entier, et 0 m M -1, ledit estimateur temporel met en oeuvre l'une des équations suivantes : zLSl = arg max l r(2)r(M / 2 -1 + z)l , z 15 2Ls2 =argmaxlr(M/2+2)r(M-1+2)l. z 5. Procédé de réception selon la revendication 2, caractérisé en ce que, pour des symboles de préambule pmn n , définis en entrée de ladite étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : 192m,n = ~ l 2 192m+1,n = J~ / 2 20 avec pmn un symbole de préambule associé à une sous-porteuse d'indice m à l'instant n, et J2 = -1 ledit estimateur temporel met en oeuvre l'une des équations suivantes : zLS1 =argmaxlr(M/4-1+2)r(M/4+2)~, z zLs2 =argmaxlr(3M/4-l+z)r(3M/4+z)l. z 25 6. Procédé de réception selon la revendication 2, caractérisé en ce que, pour des symboles de préambule pmn n , définis en entrée de ladite étape de transformation du domaine fréquentiel vers le pour 0 7. Procédé de réception selon la revendication 2, caractérisé en ce que, pour des symboles de préambule p,n n , définis en entrée de ladite étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : {p2m,n = (-1)m / 2 P2m+1,n = 0pour 0 8. Procédé de réception selon la revendication 2, caractérisé en ce que, pour des symboles de préambule p,n n , définis en entrée de ladite étape de transformation du domaine fréquentiel vers le domaine temporel mise en oeuvre en émission par : pm,n =1 si mod(m,4) = 0 - 0 sinon Pour 0 m M -1 hm,n - avec p,n n un symbole de préambule associé à une sous-porteuse d'indice m à l'instant n, ledit estimateur temporel met en oeuvre l'une des équations suivantes : 2LSt = argmaxIr(2)r(M / 2-1+2)+r(M / 4-1+2)r(M / 4+2)1 2 2LS2 = arg naxlr(M / 2+2)r(M-1+2)+r(3M / 4-1+2)r(3M / 4+2)l 2 2LS3 = arg nax l r(2)r(M -1+ 2) + r(M / 4 -1+ 2)r(3M / 4 + 2) 2 +r(M / 4 + 2)r(3M / 4 -1+ 2) + r(M /2 -1+ 2)r(M / 2+ 2)l 9. Dispositif de réception d'un signal multiporteuse de type OQAM, comprenant des moyens d'estimation (21), dans le domaine temporel, d'un décalage en temps et/ou en phase dudit signal multiporteuse, lesdits moyens d'estimation (21) comprenant au moins un estimateur pour l'estimation dudit décalage en temps, dit estimateur temporel, et/ou au moins un estimateur pour l'estimation dudit décalage en phase, dit estimateur de phase, caractérisé en ce que, ledit signal multiporteuse comprenant au moins un préambule, ledit estimateur temporel et/ou ledit estimateur de phase tient compte des coefficients d'un filtre prototype utilisé en émission pour mettre en forme ledit préambule. 10. Programme d'ordinateur comportant des instructions pour la mise en oeuvre d'un procédé selon la revendication 1 lorsque ledit programme est exécuté par un processeur.