FR 2954024 A1 20110617 FR 0958910 A 20091214 MÉTHODE D'ESTIMATION DE PARAMÈTRES OFDM PAR ADAPTATION DE COVARIANCE DESCRIPTION 5 DOMAINE TECHNIQUE La présente invention concerne une méthode d'estimation de paramètres de modulation OFDM (Orthogonal Frequency Division Multiplex). Elle trouve à s'appliquer particulièrement dans les systèmes de 10 télécommunication opportunistes (Cognitive Radio). ÉTAT DE LA TECHNIQUE ANTÉRIEURE La modulation OFDM est bien connue de l'état de la technique et est employée dans de nombreux systèmes de 15 télécommunication tels que DVB-T, ADSL, Wi-Fi (IEEE 802 a/g), WiMax (IEEE 802.16). Elle présente l'avantage d'une bonne efficacité spectrale et d'une bonne protection contre les évanouissements sélectifs en fréquence. 20 On rappelle que dans un système OFDM, les symboles d'information à transmettre sont groupés par blocs de N symboles, où N est généralement une puissance de 2, les symboles OFDM étant obtenus en effectuant une IFFT (Inverse Fast Fourier Transform) sur lesdits blocs de 25 symboles d'information. Généralement, un préfixe cyclique est rajouté au début de chaque symbole OFDM pour absorber l'interférence inter-symbole ou ISI et faciliter l'égalisation à la réception. Le signal OFDM constitué par ces symboles OFDM peut être 30 éventuellement ensuite translaté en fréquence. De manière générale, le signal émis par un système OFDM peut être représenté en bande de base par : N-1 -2in NT (t-DT,-kT, ) sa(t) N LLak,ne ga(t-kT) k n=0 où E est la puissance du signal, N est le nombre de porteuses du multiplex OFDM, akn sont les symboles d'information relatifs au bloc k, appartenant à un alphabet de modulation M -aire, typiquement BPSK, QPSK 10 ou QAM, 1/1'e est le débit des symboles d'information où Te est le temps « bribe » ou « chip », T est la durée totale du symbole OFDM avec T =(N + D)T, où D est la taille du préfixe cyclique exprimée en nombre de bribes, ga(t) est la réponse impulsionnelle du filtre de 15 mise en forme des symboles OFDM, de support temporel [O,T], destinée à apodiser le spectre du signal. On a représenté de manière schématique un signal OFDM en Fig. 1. Il est constitué d'une séquence de symboles OFDM, chaque symbole présentant une durée 20 totale (N + D)T dont une durée utile NT et un intervalle de garde de durée Tprefix = DTe , contenant le préfixe cyclique. On rappelle que, de manière classique, le préfixe cyclique est une copie de la fin du symbole OFDM à l'intérieur de l'intervalle de garde. 25 Dans certains systèmes OFDM, les préfixes cycliques sont simplement omis, autrement dit les parties utiles des symboles sont séparées par des intervalles de garde « vides ». Cette technique de transmission permet aussi (1)5 3 d'éliminer l'interférence intersymbole mais rend plus complexe l'égalisation du signal. Après propagation dans le canal de transmission, le signal OFDM reçu par le récepteur peut s'exprimer par : y, (t) = h O sa(t)+b(t) (2) où h(Dsa est la convolution entre le signal OFDM émis, sa(t) est la réponse impulsionnelle du canal de transmission h(t), et b(t) est une fonction aléatoire décrivant le bruit. On supposera que la longueur de la réponse impulsionnelle du canal est inférieure à la durée de l'intervalle de garde, de sorte que l'on pourra faire abstraction de l'interférence intersymbole (ISI). La Fig. 2 représente de manière schématique la structure d'un récepteur OFDM. Après éventuelle démodulation en bande de base, le signal reçu est échantillonné en 210 à la fréquence chip, puis les échantillons sont soumis à une conversion série/parallèle en 220 pour former des blocs de N + D échantillons. Les D premiers échantillons correspondant à l'intervalle de garde sont rejetés et le bloc des N échantillons restants correspondant à la partie utile du symbole OFDM est soumis à une FFT en 230. Les symboles démodulés obtenus sont ensuite soumis à une conversion série en 240. 4 En définitive, en supposant que le récepteur est bien synchronisé en temps et en fréquence, les symboles démodulés peuvent s'exprimer par : Zkn =hnakn+bkn (3) où hn est un coefficient complexe qui dépend de la réponse impulsionnelle du canal de transmission, et bk,n est une variable aléatoire représentant un échantillon de bruit. Le bon fonctionnement de ce récepteur nécessite une synchronisation précise en temps et en fréquence. En effet, on comprend qu'une mauvaise synchronisation en temps entraînera un glissement temporel progressif de la fenêtre de troncation et une mauvaise synchronisation en fréquence, une rotation en phase des échantillons, pouvant être représentée par un facteur multiplicatif e2`n°fa'e où 4f est l'offset de fréquence entre la fréquence de démodulation du récepteur et la fréquence porteuse du multiplex OFDM. La synchronisation temporelle et fréquentielle du récepteur est généralement réalisée grâce à l'acquisition d'une séquence d'apprentissage. Le bon fonctionnement du récepteur suppose également que l'on connaisse les paramètres de modulation du signal OFDM. Par paramètres de modulation, ou plus simplement paramètres OFDM, on entend ici la durée utile NTT, la durée de l'intervalle de garde DT', et/ou la période de répétition (N+D)TT de ces symboles, ou bien encore l'inverse de ces valeurs. On notera à cet égard que 1/NT représente l'espacement entre sous-porteuses et 1/(N + D)T, la fréquence symbole. Bien souvent le récepteur ne connaît pas a priori les paramètres de modulation OFDM et doit donc procéder 5 à leur estimation à l'aveugle avant de pouvoir démoduler le signal. On connaît par ailleurs les systèmes radio opportuniste dont on trouvera une description dans la thèse fondatrice de J. Mitola intitulée « Cognitive radio : an integrated agent architecture for software defined radio », Royal Institute of Technology, Stockholm, PhD Dissertation, 8 Mai 2000. Lorsqu'un utilisateur primaire d'un tel système utilise une modulation OFDM (WiFi, Wi-Max, LTE, WRAN, 802.22, DVB-T etc.), un utilisateur secondaire, encore dénommé opportuniste, doit être en mesure de détecter si un signal OFDM est présent ou non dans une bande donnée, et par conséquent d'estimer un paramètre OFDM d'un tel signal. Différentes méthodes ont été proposées pour estimer à l'aveugle les paramètres OFDM d'un signal. Ces méthodes exploitent, pour la plupart, la présence du préfixe cyclique dans le signal OFDM et les propriétés de cyclostationnarité qui en découlent. Les estimateurs proposés sont basés sur la fonction d'autocorrélation du signal OFDM. On trouvera notamment un exemple d'une telle méthode d'estimation dans l'article de P. Liu et al. intitulé « A blind time- parameters estimation scheme for OFDM in multi-path channel », publié dans Proc. 2005 Int'l Conference on 6 Information, Communications and Signal Processing, vol. 1, pp. 242-247, 23-26 Sept. 2005. Ces méthodes d'estimation nécessitent cependant d'acquérir un nombre élevé de symboles OFDM pour effectuer le calcul de la fonction d'autocorrélation. En outre, ces méthodes ne fonctionnent pas ou mal lorsque le signal OFDM est dépourvu de préfixe cyclique. Plusieurs méthodes ont déjà été envisagées avec succès pour estimer de manière aveugle des paramètres OFDM en absence de préfixe cyclique. Elles sont respectivement basées sur un critère à maximum de vraisemblance, un filtrage adapté et une minimisation de l'entropie en sortie d'un récepteur OFDM adaptatif. On trouvera une description de ces méthodes dans les articles de A. Bouzegzi et al., respectivement intitulés « Maximum Likelihood based methods for OFDM intercarrier characterization », IEEE International Symposium on Personal Indoor and Mobile Radio Communications (PIRMC), 2008 ; « Matched filter based algorithm for blind recognition of OFDM systems, IEEE Vehicular Technology Conference, Sept. 2008 ; et « A fourth-order based algorithm for characterization of OFDM signals », IEEE Workshop on Signal Processing and Advances in Wireless Communication (SPAWC), Juillet 2008. Toutefois, à faible rapport signal sur bruit les méthodes précitées sont inopérantes et conduisent à des erreurs d'estimation des paramètres OFDM. En outre, elles nécessitent de connaître le rapport signal sur 7 bruit, ce qui peut s'avérer délicat, en particulier dans un contexte opportuniste. Aucune des méthodes précitées ne permet de déterminer conjointement les paramètres OFDM et le rapport signal sur bruit. Or, une telle méthode d'estimation permettrait à un récepteur de déterminer rapidement l'état du canal OFDM et de retourner à l'émetteur une information dite CSI (Channel State Information) le décrivant. L'information CSI est avantageusement utilisée par l'émetteur pour effectuer un ordonnancement de fréquences, une adaptation de liaison, etc. Le but de la présente invention est par conséquent de proposer une méthode d'estimation de paramètres de modulation d'un signal OFDM, avec un taux de succès élevé, même à faible rapport signal sur bruit et ce, que le signal soit ou non dépourvu de préfixe. Un but subsidiaire de la présente invention est de permettre à un récepteur d'effectuer une estimation conjointe des paramètres OFDM et du rapport signal sur bruit. EXPOSÉ DE L'INVENTION La présente invention est définie par une méthode d'estimation d'au moins un paramètre de modulation d'un signal OFDM, ledit signal étant échantillonné pendant une fenêtre temporelle pour fournir une séquence d'échantillons, dans laquelle : (a) on regroupe des échantillons consécutifs de ladite séquence, en fonction d'une 8 valeur testée du dit paramètre, pour former une pluralité de vecteurs ; (b) on estime la matrice de covariance des vecteurs ainsi obtenus ; (c) on calcule la matrice de covariance théorique que l'on obtiendrait à partir d'un signal OFDM dont ledit paramètre de modulation aurait ladite valeur testée ; (d) on mesure une distance entre la matrice de covariance et la matrice de covariance théorique ; et en ce que l'on itère les étapes (a), (b), (c) et (d), la valeur estimée du paramètre étant celle réalisant la distance la plus faible à l'étape (d). Selon une variante de réalisation, la matrice de covariance théorique est obtenue en fonction de ladite valeur testée du paramètre de modulation, d'une valeur testée de la puissance du signal et d'une valeur testée de la puissance de bruit, les étapes (c) et (d) étant alors itérées pour une pluralité de valeurs testées de puissance de signal et de bruit. Avantageusement, la méthode d'estimation peut également fournir une estimation du rapport signal sur bruit à partir des valeurs de puissance de signal et de bruit minimisant ladite distance conjointement avec la valeur estimée de paramètre de modulation. La matrice de covariance théorique peut être calculée au moyen de : a'6b~=à FeFY+(VI30 où e, à',', à sont respectivement les valeurs testées du paramètre de modulation, de la puissance de signal et de la puissance de bruit ; le est la matrice unité de taille N x N où est le nombre de sous-porteuses du 5 multiplex OFDM, connu ou supposé, et la matrice t- est une matrice de taille PxN avec P=[(NTe+DTe)/Te~ dont les éléments sont donnés par : Fa ùù2 nnmNTù2 nnNT, m,n N N où NTe est la durée utile des symboles OFDM, connue ou testée, DTe est la durée du préfixe des symboles OFDM, connue ou testée, et Te est la période d'échantillonnage du signal OFDM. 15 La matrice de covariance quant à elle peut être estimée par : Kù1 H Ryy (Ô') ùLykYk K k=o T 20 où Yk =(yk,o,...ykP_1) est le kième vecteur de ladite pluralité de vecteurs, yk p = ya (pT +DTe+k( NTe+DTC)) où NTe est la durée utile des symboles OFDM, connue ou testée, DTe est la durée du préfixe des symboles OFDM, connue ou testée, Te est la période d'échantillonnage 25 du signal OFDM et K= avec P=L(NT,+DT,)/Te J et M 10 est le nombre total d'échantillons dans la fenêtre temporelle. La distance entre la matrice de covariance et la matrice de covariance théorique peut être calculée au moyen de : _ 2 ` l C O M E T ( , a'gb)= R> ( 6 ) - R> ( ,da'd b) F où II'IIF est la norme de Frobenius. Alternativement, la distance entre la matrice de covariance et la matrice de covariance théorique peut être calculée au moyen de : _ K 2 JMCOMET( '~a'ab)ùP2 RYYùRYY( ~da,db)M F où II'IIF est la norme de Frobenius. Le paramètre de modulation peut être la durée totale d'un symbole OFDM, la durée utile d'un symbole OFDM, la durée du préfixe d'un symbole OFDM, le nombre de sous-porteuses d'un symbole OFDM. Ladite méthode d'estimation peut être également appliquée à un ensemble des paramètres de modulation constitué de la durée utile d'un symbole OFDM, de la durée de préfixe d'un symbole OFDM et du nombre de sous-porteuses du multiplex OFDM, les valeurs estimées de ces paramètres étant alors les valeurs testées minimisant conjointement ladite distance. 11 BRÈVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention fait en référence aux figures jointes parmi lesquelles : La Fig. 1 représente de manière schématique un signal OFDM ; La Fig. 2 représente de manière schématique un récepteur OFDM connu de l'état de la technique ; La Fig. 3 représente un organigramme de la méthode d'estimation d'un paramètre de modulation OFDM selon l'invention ; La Fig. 4 représente un organigramme de la méthode d'estimation de paramètres de modulation OFDM selon un premier mode de réalisation de l'invention ; La Fig. 5 représente un premier exemple de fonction de coût utilisée dans la méthode de la Fig. 4; La Fig. 6 représente un second exemple de fonction de coût utilisée dans la méthode de la Fig. 4; La Fig. 7 représente un premier exemple de fonction de coût normalisée utilisée dans un second mode de réalisation de l'invention ; La Fig. 8 représente un second exemple de fonction de coût normalisée utilisée dans le second mode de réalisation de l'invention ; La Fig. 9 donne le taux de bonne estimation d'un paramètre de modulation OFDM en fonction du rapport signal sur bruit, dans le cas d'une estimation conjointe ; 12 La Fig. 10 donne le taux de bonne estimation de la variance du signal en fonction de ce même rapport, dans le cas d'une estimation conjointe La Fig. 11 compare plusieurs méthodes d'estimation d'un paramètre OFDM en fonction du rapport signal sur bruit, dans le cas d'une synchronisation parfaite du récepteur et pour une faible durée de préfixe ; La Fig. 12 compare plusieurs méthodes d'estimation d'un paramètre OFDM en fonction du rapport signal sur bruit, dans le cas d'une synchronisation parfaite du récepteur et en absence de préfixe ; La Fig. 13 compare plusieurs méthodes d'estimation d'un paramètre OFDM en fonction du rapport signal sur bruit, en absence de synchronisation du récepteur et en absence de préfixe. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Nous considérerons dans la suite le cas d'un système de télécommunication utilisant une modulation OFDM. Nous supposerons que le signal OFDM émis a la forme donnée par l'expression (1) et que les symboles d'information peuvent être représentés par des variables aléatoires indépendantes, identiquement distribuées et prenant leurs valeurs dans un alphabet de modulation M -aire. Après démodulation en bande de base, le signal reçu est échantillonné à une fréquence 1 {~ /eT , Te supérieure à la largeur de la bande OFDM considérée 30 pour satisfaire au critère de Nyquist. Ainsi, l'on est 13 certain que la période d'échantillonnage Te est inférieure à la période chip T. Soit T° la durée de la fenêtre d'observation, le signal à temps discret obtenu est composé de M° = T° échantillons (où Lx] désigne le plus grand entier inférieur à x). Si le récepteur connaît les paramètres NT, DT et N du signal OFDM, les échantillons peuvent être regroupés par symboles OFDM de taille P=L(NT +DT )/T J et l'expression (2) peut être alors réécrite sous forme vectorielle: yk = Gak + b (4) où ak est le k-ième vecteur de taille N des symboles transmis, b est un vecteur de taille P, représentant le bruit, supposé additif et gaussien et yk est un vecteur de taille P constitué de P échantillons consécutifs ykp = ya(pT +DT +k(NT +DT )) , p=O,...,Pû1. Les P échantillons de bruit sont des variables aléatoires indépendantes entre elles et du signal utile ; G = HF où H est la matrice de taille N> Of représentative de la réponse fréquentielle du canal et F est une matrice de taille PxN dont les éléments sont donnés par : -2i7Gnm é -2ir~nDT F = e NT e NT m,n NIN (5) 14 Par souci de simplification, mais sans perte de généralité, on supposera dans la suite que le canal de transmission est sans évanouissement, autrement dit que la matrice H est égale à la matrice identité à un coefficient multiplicatif près. La relation (4) peut alors s'écrire sous la forme: yk = Fak + b (6) Considérons maintenant le cas où le récepteur opère en aveugle. Les paramètres NT, DT et N, caractérisant le signal OFDM sont alors des inconnues que l'on peut regrouper sous la forme d'un paramètre à estimer 9 =(NT,,DT,,N). Les échantillons du signal OFDM reçus par le récepteur sont regroupés sous forme de vecteurs de taille P=[(NTc+DTc)/Te], notés yak, avec: 3ëk =Faaak+bak (7) où, de manière similaire à la relation (6), a4k est le vecteur de taille N des symboles transmis, bak est un vecteur de taille P représentant les échantillons de bruit additif gaussien, et Fa est la matrice de taille 25 PxN dont les éléments sont donnés par : Fe = 1 e -27GnmNT e -2nn F NT m,n N (8) 15 On comprendra que yak est un vecteur de P variables aléatoires dont la matrice de covariance vaut . ,Q,6b)=E(yaya )=aa2t +ab21 (9) où 6b et 6b représentent respectivement la variance des symboles transmis et du bruit, et la est la matrice 10 unité de taille Nx N. Avec ces conventions, le rapport a,/6b n'est autre que le rapport signal sur bruit. La matrice de covariance dépend de manière générale du paramètre 9 et des puissances respectives du signal utile et du bruit. 15 L'idée à la base de la présente invention est d'estimer au moins un paramètre de modulation du signal OFDM en utilisant une méthode d'adaptation de la covariance. 20 La Fig. 3 représente de manière générale un organigramme de la méthode d'estimation d'un paramètre de modulation OFDM selon l'invention. A l'étape 310, on échantillonne le signal OFDM en bande de base pendant une fenêtre temporelle de largeur 25 donnée. On obtient ainsi une séquence d'échantillons. A l'étape 320, on regroupe par paquets des échantillons de la séquence ainsi obtenue, en fonction d'une valeur testée, 9, du paramètre à estimer, chaque paquet étant formé par des échantillons consécutifs de5 16 ladite séquence et étant représenté par un vecteur noté Yék A l'étape 330, on estime la matrice de covariance R,(6) des vecteurs ainsi obtenus. A l'étape 340, on calcule ensuite la matrice de covariance théorique R06) que l'on obtiendrait à partir d'un signal OFDM dont ledit paramètre de modulation aurait ladite valeur testée O. A l'étape 350, on calcule une distance entre la 10 matrice de covariance en 330 et la matrice de covariance théorique obtenue en 340. Les étapes 320 à 350 sont répétées pour une pluralité de valeurs possibles du paramètre O. la valeur estimée, 0, de ce paramètre est obtenue en 360 15 comme celle minimisant la distance calculée à l'étape 350. L'homme du métier comprendra que l'on obtient la valeur 0 qui conduit à la matrice de covariance estimée la plus proche de la matrice de covariance théorique. 20 La Fig. 4 représente de manière plus détaillée la méthode d'estimation d'au moins un paramètre de modulation OFDM selon un premier mode de réalisation de l'invention. 25 A l'étape 410, on échantillonne le signal reçu en bande de base, y, pendant une fenêtre temporelle de largeur To. La période d'échantillonnage T est choisie suffisamment faible pour respecter le critère de Nyquist, compte tenu de la largeur de bande OFDM 17 d'intérêt. On obtient ainsi une séquence de Mo échantillons. A l'étape 420, on choisit une valeur à tester du paramètre 6, c'est-à-dire des valeurs de durée utile NTe, de durée de préfixe DTe, et de nombre de sous-porteuses N. On choisit également une valeur de puissance du signal transmis P = d- et une valeur de puissance de bruit P1,= 6b. Ce choix peut être aveugle ou semi-aveugle si le récepteur connaît l'une des valeurs précitées. On choisit aussi une valeur JcoMET arbitrairement grande. A l'étape 430, on regroupe des échantillons consécutifs de la séquence pour former k vecteurs yk, k=O,...,Kù1, avec K=[ p, de taille P=[(NT,+DTc)/Te] . Les P composantes yk p , p = 0,...,Pù1 du vecteur Yk sont donc obtenues par : ykp =ya(PT +DTe+k(NTe+DTe)) (10) La valeur testée du paramètre de modulation, 6, détermine la taille de ces vecteurs ainsi que leurs points de départ dans la séquence. A l'étape 440, on estime la matrice de covariance RYY(6) en calculant : /_ K-1 Ryy (e ) ùLYkYk \ K k=o T avec Yk =(Yko,....ykp-1) A l'étape 460, on détermine la distance, au sens de la norme de Frobenius, entre la matrice de covariance Ryy e estimée à partir de (11) et la matrice Ryy(9) théorique obtenue préalablement en 450 par : Ryy(8,6a,6b)=a FeF;+6bl (12) où 6â est une valeur testée de la puissance du signal et 6b est une valeur testée de la puissance de bruit. On rappelle que la norme de Frobenius d'une matrice A est définie par A F =Tr(AHA) où Tr(.) est la fonction trace, autrement dit on calcule la valeur de la fonction de coût définie par : z JCOMET(e,6a'6b)= MRùyy (8)ùR,,(6,da,6b) F On compare cette valeur à la valeur JCOMET et si 18 (13) min `~COMET 8(,da l 6 b A l'étape 470, on vérifie si l'on a parcouru un ensemble ordonné prédéfini de paramètres (0,6a,6b). A défaut, on passe à l'élément suivant de cet ensemble en 475 et l'on retourne à l'étape 430. Cet ensemble ordonné est par exemple le produit cartésien d'ensembles élémentaires IIe XII6a XII6b où H est un ensemble de valeurs possibles du paramètres prévues par une norme OFDM, et F4,, F4, sont des intervalles de valeurs de puissance. La relation d'ordre sur cet ensemble peut être l'ordre lexicographique. A l'étape 480, on estime les valeurs des paramètres 0, 6a, 6b à l'aide de . (0,6a,6b ) = argminJcoMET (0Ida1db ) 8,6a,6b autrement dit, on récupère les valeurs de ces paramètres dans la zone mémoire précitée. On comprend que la méthode d'estimation 20 représentée en Fig. 4 permet d'estimer conjointement les paramètres de modulation du signal OFDM (représentés par 0) ainsi que les niveaux de puissance de signal et de bruit. Cette méthode permet par conséquent non seulement de déterminer NT, Ln', et N 25 mais aussi par exemple le rapport signal sur bruit G,2/6b2. Il est clair pour l'homme du métier que si certains paramètres de modulation sont déjà connus du récepteur, les valeurs connues de ces paramètres (14)15 20 remplacent leurs valeurs testées (9,6a,6b) dans le calcul des expressions (12) et (13). Dans ce cas, la dimension correspondante du pavé n'est pas parcourue en 475. On a illustré en Fig. 5 un exemple de la fonction NT, de coût JcoMET en fonction de la variable réduite dans le cas où seul le paramètre NT, est supposé inconnu, les autres paramètres prenant les valeurs 10 connues, DT, =DT , N=N, da =6a 6b= 6b. Le rapport signal sur bruit 6â/6b est ici égal à 10 dB. On remarque que la fonction de coût JcoMET est convexe et présente un minimum prononcé en NT, = NT , ce qui permet une estimation précise de la longueur utile du symbole 15 OFDM. De manière similaire, si la longueur utile était connue et si l'on recherchait par exemple la longueur de préfixe ou le nombre de sous-porteuses, la fonction de coût aurait encore un minimum prononcé pour la valeur réelle correspondante. 20 La Fig. 6 représente un autre exemple de fonction de coût dans les mêmes conditions que la Fig. 5, à la différence près que le rapport signal sur bruit est ici égal à 0 dB, autrement dit 6a= 6b. On remarque que la 25 fonction de coût présente toujours un minimum pour NT,= NT mais qu'en revanche elle n'est plus convexe. De manière plus générale, on constate que la fonction de coût perd sa propriété de convexité aux faibles NT niveaux de rapport signal sur bruit. Cette absence de convexité peut être préjudiciable à la recherche du minimum de la fonction, par exemple si la cette recherche n'est pas effectuée au moyen d'un balayage systématique mais d'une descente de gradient. Dans le cas d'un faible rapport signal sur bruit, on préfère utiliser un second mode de réalisation décrit ci-après. Le second mode de réalisation met en oeuvre une fonction de coût normalisée définie par : JCOMET (ô,6a9 b JMCOMET ( a, ~ b ) ù E(ICOMET (Ô96a9db» où E(JcoMET(e,6a,6b)) désigne la moyenne de la fonction de coût. On peut montrer que pour les faibles valeurs de rapport signal sur bruit, c'est-à-dire pour 6a « 6b, cette moyenne vaut : _ PZ 4 E (J COMET,àa,àb) K 6b Etant donné que 6b est une constante, la fonction de coût normalisée peut se réduire à : JMCOMET ,6a,6b COMET laalab (17) La méthode d'estimation selon le second mode de réalisation de l'invention diffère de celle illustrée (15) 25 (16) 22 en Fig. 4 en ce que la fonction de coût JMCOMET est utilisée en lieu et place de JCOMET . La description des étapes de traitement ne sera donc pas reprise ici. On a représenté en Fig. 7 un premier exemple de fonction de coût normalisée, dans les mêmes conditions que celles de la Fig. 5. De même, on a représenté en Fig. 8 un premier exemple de fonction de coût normalisée, dans les mêmes conditions que celles de la Fig. 6. On remarque que la fonction de coût normalisée est convexe et présente un minimum prononcé dans les deux cas, en particulier à faible niveau de rapport signal sur bruit, cf. Fig. 8. La méthode d'estimation selon le second mode de réalisation est par conséquent robuste vis-à-vis du bruit. La Fig. 9 représente le taux de bonne estimation (ou encore d'identification) du paramètre NT en fonction du rapport signal sur bruit, dans le cas d'une estimation conjointe de NT, 6Q et 6b. On remarque que le taux de bonne estimation est égal à 1 avec la fonction de coût normalisée JMCOMET dès que le rapport signal sur bruit est supérieur à -8 dB. En revanche, ce taux de bonne estimation n'est atteint avec la fonction de coût JcoMET que pour un rapport signal sur bruit supérieur à 2 dB. La Fig. 10 représente le taux de bonne estimation de la variance 6Q en fonction du rapport signal sur bruit, dans le cas d'une estimation conjointe de NT, 6Q et 6b. On remarque là encore qu'un taux de bonne 23 estimation égal 1 est atteint dès -8 dB pour JMCOMET et seulement à partir de 2 dB pour la fonction de coût JCOMET ' Les Figs. 11, 12 et 13 comparent les performances de différentes méthodes d'estimation d'un paramètre OFDM, ici le paramètre NT, en termes de taux de bonne estimation, en fonction du rapport signal sur bruit. On a désigné par I et II les courbes relatives aux méthodes d'estimation connues, basées respectivement sur la corrélation et la corrélation cyclique du signal reçu, telles que décrites dans l'article de A. Bouzezgi et al. intitulé « A second order statistics based algorithm for blind recognition of OFDM based systems » publié dans IEEE Global Telecommunications Conference, Nov. 2008, et par III et IV les courbes relatives respectivement au second et au premier modes de réalisation de l'invention. On a pris comme base de comparaison un signal OFDM de N=64 sous-porteuses, de durée utile NT =3.2 µs et de rapport Tc/Te=2. On a supposé en Fig. 11 que la durée du préfixe était égale au 1/32 ième de la durée utile du symbole OFDM (c'est-à-dire D/N =1/32). Dans les Figs. 12 et 13, on supposé que les symboles OFDM étaient dépourvus de préfixe (1)=0). Enfin, dans les Figs. 11 et 12, on a supposé que le récepteur était parfaitement synchronisé sur le signal OFDM, sans offset temporel (début de la fenêtre d'échantillonnage calé sur le début d'un symbole OFDM) ni fréquentiel (fréquence de démodulation en bande de base égale à la fréquence de la porteuse). En revanche, en Fig. 13, le récepteur 24 présente un offset temporel et un offset fréquentiel aléatoires par rapport au signal OFDM. On notera que lorsque la durée du préfixe est courte ou lorsque le préfixe est absent les méthodes d'estimation du paramètre temporel NT selon l'invention conduisent à des taux de bonne estimation sensiblement plus élevés qu'avec les méthodes connues de l'art antérieur. Ceci est d'ailleurs d'autant plus vrai que le préfixe est court. En cas d'absence totale de préfixe, les méthodes traditionnelles sont totalement inopérantes, cf. Fig. 12, alors que les méthodes selon l'invention continuent à donner de très bons résultats, et ce, même lorsque le récepteur n'est pas synchronisé, cf. Fig. 13. La méthode d'estimation selon l'invention peut être mise en oeuvre dans un récepteur OFDM ou un récepteur d'utilisateur secondaire dans un système radio opportuniste, à l'aide de circuits dédiés ou de modules logiciel exécutés par un microprocesseur, de manière connue en soi. REVENDICATIONS1. Méthode d'estimation d'au moins un paramètre de modulation d'un signal OFDM, ledit signal étant échantillonné pendant une fenêtre temporelle pour fournir une séquence d'échantillons, caractérisée en ce que : (a) on regroupe (310) des échantillons consécutifs de ladite séquence, en fonction d'une valeur testée (9) du dit paramètre, pour former une pluralité de vecteurs ; (b) on estime (320) la matrice de covariance ( Ryyle) ) des vecteurs ainsi obtenus ; (c) on calcule (330) la matrice de covariance théorique (Ryy(6)) que l'on obtiendrait à partir d'un signal OFDM dont ledit paramètre de modulation aurait ladite valeur testée ; (d) on mesure une distance (340) entre la 20 matrice de covariance et la matrice de covariance théorique ; et en ce que l'on itère les étapes (a), (b), (c) et (d), la valeur estimée du paramètre étant celle réalisant la distance la plus faible à l'étape (d). 25 2. Méthode d'estimation selon la revendication 1, caractérisée en ce que la matrice de covariance théorique est obtenue en fonction de ladite valeur testée du paramètre de modulation, d'une valeur 30 testée de la puissance du signal (à) et d'une valeur 26 testée de la puissance de bruit (6b), les étapes (c) et (d) étant alors itérées pour une pluralité de valeurs testées de puissance de signal et de bruit. 3. Méthode d'estimation selon la revendication 2, caractérisée en ce qu'elle fournit en outre une estimation du rapport signal sur bruit à partir des valeurs de puissance de signal et de bruit minimisant ladite distance conjointement avec la valeur estimée de paramètre de modulation. 4. Méthode d'estimation selon la revendication 2 ou 3, caractérisée en ce que la matrice de covariance théorique est calculée au moyen de : a,6b)=a FeF- +6bl où 9, 6a, 6b sont respectivement les valeurs testées du paramètre de modulation, de la puissance de signal et 20 de la puissance de bruit ; Ië est la matrice unité de taille N x N où est le nombre de sous-porteuses du multiplex OFDM, connu ou supposé, et la matrice Fe est une matrice de taille PxN avec P=[(NT,+DT,)/Te~ dont les éléments sont donnés par : -2i7Lnm é -2i7Ln~T` e NT`e NT` V1~ 25où NTe est la durée utile des symboles OFDM, connue ou testée, DTe est la durée du préfixe des symboles OFDM, connue ou testée, et Te est la période d'échantillonnage du signal OFDM. 5. Méthode d'estimation selon l'une quelconque des revendications précédentes, caractérisée en ce que la matrice de covariance est estimée par : K-1 Ryy \e f ùLykYk / g k=o H T où Yk =(yk,o'•••yk,P_1) est le kième vecteur de ladite pluralité de vecteurs, ykp =(pT +DTe+k(NTC+DTe)) où NTe est la durée utile des symboles OFDM, connue ou testée, 15 DTe est la durée du préfixe des symboles OFDM, connue ou testée, Te est la période d'échantillonnage du signal OFDM et K=L P avec P=[(NTe+DTe)/Te] et M est le nombre total d'échantillons dans la fenêtre temporelle. 20 6. Méthode d'estimation selon la revendication 5, caractérisée en ce que la distance entre la matrice de covariance et la matrice de covariance théorique est calculée au moyen de : 10 25 ~COMET(~, a'gb)= MR-> (6)ùR> (,da 'db)Foù II.IIE est la norme de Frobenius. 7. Méthode d'estimation selon la revendication 5, caractérisée en ce que la distance entre la matrice de covariance et la matrice de covariance théorique est calculée au moyen de : _ K 2 `lMCOMET( ' a'ab)ùP2 RYYùRYY( d a , b)M F où II•IIE est la norme de Frobenius. 8. Méthode d'estimation selon l'une des revendications précédentes, caractérisée en ce que ledit paramètre de modulation est la durée totale d'un symbole OFDM. 9. Méthode d'estimation selon l'une des revendications 1 à 7, caractérisée en ce que ledit paramètre de modulation est la durée utile d'un symbole OFDM. 10. Méthode d'estimation selon l'une des revendications 1 à 7, caractérisée en ce que ledit paramètre de modulation est la durée du préfixe d'un symbole OFDM. 11. Méthode d'estimation selon l'une des revendications 1 à 7, caractérisée en ce que ledit paramètre de modulation est le nombre de sous-porteuses d'un symbole OFDM. 29 12. Méthode d'estimation selon l'une des revendications 1 à 7, caractérisée en ce qu'elle est appliquée à un ensemble des paramètres de modulation constitué de la durée utile d'un symbole OFDM, de la durée de préfixe d'un symbole OFDM et du nombre de sous-porteuses du multiplex OFDM, les valeurs estimées de ces paramètres étant alors les valeurs testées minimisant conjointement ladite distance.10