FR 2983666 A1 20130607 FR 1161061 A 20111201 PROCEDE D'ESTIMATION D'UN CANAL RADIOELECTRIQUE La présente invention concerne un procédé d'estimation d'un canal de propagation radioélectrique entre un émetteur et un récepteur. Elle concerne 5 également un équipement récepteur et un programme d'ordinateur correspondants. L'invention se rapporte plus particulièrement aux systèmes de transmission multiporteuses, par exemple de type à multiplexage par division en fréquences orthogonales, plus connu sous l'acronyme anglais OFDM (« Orthogonal 10 Frequency Division Multiplexing »). Ce type de transmission est de plus en plus fréquemment utilisé. Il a notamment été adopté dans les systèmes LTE (« Long Term Evolution »), TEDS (« TETRA Enhanced Data Service »), DAB (« Digital Audio Broadcasting ») et DVB-T (« Digital Video Broadcasing-Terrestrial »). 15 Dans les systèmes de transmission OFDM, les données sont en général organisées en trames. Dans chaque trame, certains symboles, répartis dans le plan temps-fréquence sont insérés parmi les informations utiles destinées au récepteur. Ces symboles, appelés "symboles pilotes" sont connus de l'émetteur et du récepteur. Ils sont utilisés à des fins de synchronisation et d'estimation du 20 canal de propagation. Un canal de propagation est défini comme étant le chemin qu'emprunte un signal émis depuis un émetteur vers un récepteur. Ce canal comprend une pluralité de trajets entre l'émetteur et le récepteur, chaque trajet étant caractérisé notamment par un retard et une atténuation. L'effet de ce canal multi-trajets, dans 25 le cas d'un canal non sélectif en fréquence, sur un symbole de données émis est généralement modélisé simplement par un coefficient multiplicatif complexe que le récepteur cherche à estimer afin de récupérer le symbole émis avec le minimum possible d'erreurs. L'estimation de canal est définie comme étant la détermination de ces coefficients pour tous les symboles émis. 30 Conventionnellement, le récepteur commence par estimer le canal à la position des symboles pilotes puis il effectue une interpolation en temps et en fréquence afin d'estimer le canal sur l'ensemble de la trame. Le canal estimé est généralement modélisé par un vecteur colonne, dit vecteur canal, ayant un nombre de lignes égal au nombre de symboles dans une trame. Cependant, cette estimation de canal n'est pas optimale car elle est très sensible aux variations du canal et ne permet pas de réduire efficacement les effets liés au bruit. Le document FR2814011 décrit un procédé d'estimation de canal ne présentant pas les inconvénients des procédés conventionnels. Cependant, le procédé décrit est très contraignant en termes de complexité. La présente invention vise à améliorer la situation. A cet effet, l'invention concerne tout d'abord un procédé d'estimation d'un canal de propagation radioélectrique entre un émetteur et un récepteur, ledit émetteur émettant un signal comprenant des trames de symboles répartis en temps et en fréquence parmi lesquels certains symboles, dits symboles pilotes, sont connus du récepteur, caractérisé par le fait qu'il comporte les étapes mises en oeuvre dans le récepteur de : - calcul d'une matrice MF de covariance du canal dans le domaine fréquentiel; - calcul d'une matrice MT de covariance du canal dans le domaine temporel; - décomposition des matrices MF et MT de covariance en vecteurs propres selon les relations MF = WFHNFWF et MT = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matrice MF; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MF; - WT est une matrice de vecteurs propres de la matrice MT ; - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MT ; - calcul du produit de Kronecker des matrices de vecteurs propres WF et WT pour obtenir une matrice de vecteurs propres globale W ; - calcul d'une matrice diagonale de valeurs propres globale N égale au 30 produit de Kronecker des matrices de valeurs propres NF et NT ; et - estimation du canal de propagation à l'aide des symboles pilotes et des matrices de vecteurs propres W et de valeurs propres N en utilisant le critère du maximum a posteriori ou de la minimisation de l'erreur quadratique. On entend par symbole pilote tout symbole connu du récepteur, c'est-à-dire à la fois les symboles pilotes, tels que désignés dans les normes des systèmes OFDM connus comme LTE et TEDS, et/ou les symboles de synchronisation et/ou les symboles déjà démodulés et décidés par le récepteur dans une étape antérieure de traitement. La présente invention exploite ainsi la propriété du canal de propagation radio selon laquelle les contraintes que doit respecter le canal dans le domaine fréquentiel sont indépendantes des contraintes qu'il doit respecter dans le domaine temporel. Cela permet d'exprimer la matrice de covariance globale du canal comme le produit de Kronecker des deux matrices de covariance MF et MT correspondant respectivement au domaine fréquentiel et au domaine temporel. Le produit de Kronecker de deux matrices A et B, la matrice A ayant pour composantes (ag), où i est un entier compris entre 1 et m et j est un entier compris entre 1 et n, est le produit matriciel noté AOB et défini par l'expression suivante : C111" - - - - - - Clin" AOB= a iniB - - - - - - a innB L'exploitation de la séparabilité de la matrice de covariance selon l'invention est très favorable en terme de complexité. Ainsi, pour une trame de n symboles, avec nf sous-porteuses et nt temps-symboles, la matrice de covariance globale est de dimension (n,n) avec n=nf.n, . Cette matrice est égale au produit de Kronecker d'une matrice de covariance MF de dimension (lof) avec une matrice de covariance MT de dimension (nt,nt). Il suffit donc de mémoriser (nf2+nt2) valeurs au lieu de nf2.nt2 valeurs. La matrice W des vecteurs propres étant le produit de Kronecker des deux matrices WF et WT de vecteurs propres, il suffit également de (nf.nf+nt.nt) = nf2+nt2 valeurs pour la représenter. De plus, la matrice diagonale N des valeurs propres étant le produit de Kronecker des deux matrices diagonales NF et NT, il suffit de nf+nt valeurs pour la représenter. Le procédé de la présente invention permet ainsi de réduire de manière significative la complexité de l'estimation de canal, particulièrement en termes de mémoire de stockage nécessaire. Selon une réalisation, le critère utilisé est le critère du maximum a posteriori et l'étape d'estimation du canal comprend les sous-étapes de : - calcul d'une matrice P selon la relation P = T.W.EH.E.WH. dans laquelle E est une matrice diagonale comprenant les symboles pilotes aux positions des symboles pilotes, les autres éléments de la diagonale étant nuls ; - décomposition de la matrice P en vecteurs propres selon la relation P = XHQX dans laquelle X est une matrice de vecteurs propres de la matrice P et Q est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice P ; et - estimation du canal selon la relation C = WH..XH.(Q+a2.1)1.X.T.W.EH.R dans laquelle C est un vecteur représentant le canal, I est la matrice identité, R est un vecteur comprenant les symboles reçus et a 2représente la puissance du bruit dans le canal. Plus précisément, la matrice est une matrice diagonale dont les termes de la diagonale sont les racines carrées des termes de la diagonale de la matrice diagonale N. Lorsque les symboles pilotes sont régulièrement répartis en temps et en 20 fréquence, la matrice P peut être écrite sous la forme d'un produit de Kronecker de deux matrices relatives respectivement au domaine fréquentiel et au domaine temporel. Dans ce cas, la complexité de l'estimation de canal est encore réduite. Selon une réalisation préférée, l'étape d'estimation du canal comprend les sous étapes de : 25 - extraction d'une matrice W' à partir de la matrice W et d'une matrice N' à partir de la matrice N, la matrice N' comprenant un nombre déterminé n' de valeurs propres de la matrice N et la matrice W' comprenant les vecteurs propres associés à ces n' valeurs propres; et - estimation du canal de propagation à l'aide des symboles pilotes et des 30 matrices de vecteurs propres W' et de valeurs propres N' en utilisant le critère du maximum a posteriori ou de minimisation de l'erreur quadratique. De préférence, les n' valeurs propres sont les plus grandes valeurs propres de la matrice N . Cela permet de réduire la complexité de l'estimation du canal tout en minimisant la perte de performances grâce au fait que les vecteurs propres conservés dans la matrice W' représentent la plus large part de l'énergie, représentée par la somme des valeurs propres conservées dans la matrice N'. Avantageusement, les n' vecteurs propres de la matrice W' sont les produits de Kronecker d'un nombre limité n'f de vecteurs propres de la matrice WF et d'un nombre limité n't de vecteurs propres de la matrice WT et les n' valeurs propres de la matrice N' sont les produits des valeurs propres de la matrice NF correspondant aux n'f vecteurs propres de la matrice WF et des valeurs propres de la matrice NT correspondant aux n't vecteurs propres de la matrice WT Dans ce cas, la matrice W' est séparable en temps et en fréquence. Il suffit alors de (n'f.nf+n't.nt) valeurs pour la représenter. Ainsi, la complexité du procédé est encore plus réduite. De préférence, n' est inférieur ou égal au nombre de symboles pilotes dans chaque trame. Avantageusement, le critère utilisé est le critère du maximum a posteriori et l'étape d'estimation du canal comprend les sous-étapes de : - calcul d'une matrice selon la relation r =N' dans laquelle E est une matrice diagonale comprenant les symboles pilotes aux positions des symboles pilotes, les autres éléments de la diagonale étant nuls ; - décomposition de la matrice P' en vecteurs propres selon la relation = X/HQ/X/ dans laquelle X' est une matrice de vecteurs propres de la matrice P' et Q' est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice P' ; et - estimation du canal selon la relation C=WH.-N/W.X/H.(Q/-Ecy2.1)1.X'.11/.W.EH.R dans laquelle C est un vecteur représentant le canal, I est la matrice identité, R est un vecteur comprenant les symboles reçus et 6 2représente la puissance du bruit dans le canal. Plus précisément, la matrice1T' est une matrice diagonale dont les termes de la diagonale sont les racines carrées des termes de la diagonale de la matrice diagonale N'. Lorsque les symboles pilotes sont régulièrement répartis en temps et en fréquence, la matrice P' peut être écrite sous la forme d'un produit de Kronecker de deux matrices relatives respectivement au domaine fréquentiel et au domaine temporel. Dans ce cas, la complexité de l'estimation de canal est encore réduite. En variante, le critère utilisé est le critère de la minimisation de l'erreur quadratique et le canal est estimé selon la relation c = w'H.N1N' .R dans laquelle C est un vecteur représentant le canal, I est la matrice identité, R est un vecteur comprenant les symboles reçus et a 2 représente la puissance du bruit dans le canal. Cette variante permet d'obtenir des résultats proches de ceux obtenus par l'application du critère du maximum a posteriori. Elle est de complexité équivalente dans le cas où la matrice W' est séparable en temps et en fréquence. Avantageusement, le signal est un signal multiporteuses, notamment un signal OFDM. L'invention peut également s'appliquer à des systèmes monoporteuses. L'invention concerne également un équipement récepteur apte à recevoir un signal transmis à travers un canal de propagation radioélectrique, ledit signal comprenant des trames de symboles répartis en temps et en fréquence parmi lesquels certains symboles, dits symboles pilotes, sont connus dudit récepteur, caractérisé par le fait qu'il comporte des moyens de : - calcul d'une matrice MF de covariance du canal dans le domaine fréquentiel; - calcul d'une matrice MT de covariance du canal dans le domaine temporel; - décomposition des matrices MF et MT de covariance en vecteurs propres selon les relations MF = WFHNFWF et MT = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matriceMF ; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MF; - WT est une matrice de vecteurs propres de la matrice MT ; - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MT ; - calcul du produit de Kronecker des matrices de vecteurs propres WF et WT pour obtenir une matrice de vecteurs propres globale W ; - calcul d'une matrice diagonale de valeurs propres globale N égale au produit de Kronecker des matrices de valeurs propres NF et NT ; et - estimation du canal de propagation à l'aide des symboles pilotes et des matrices de vecteurs propres W et de valeurs propres N en utilisant le critère du maximum a posteriori ou de la minimisation de l'erreur quadratique. L'invention concerne aussi un programme d'ordinateur comportant des instructions pour la mise en oeuvre du procédé de l'invention lorsque le programme est exécuté par au moins un processeur. Les organigrammes des figures 2 à 4 illustrent schématiquement le déroulement de ce programme d'ordinateur, pour des modes de réalisation préférés de l'invention. Des modes de réalisation de l'invention vont maintenant être décrits de façon plus précise mais non limitative en regard des dessins annexés sur lesquels : - la figure 1 est un schéma d'un équipement récepteur selon l'invention ; - la figure 2 est un organigramme illustrant le fonctionnement du procédé d'estimation de canal selon un premier mode de réalisation de l'invention; - la figure 3 est un organigramme illustrant le fonctionnement du procédé d'estimation de canal selon un deuxième mode de réalisation de l'invention ; et - la figure 4 est un organigramme illustrant le fonctionnement du procédé d'estimation de canal selon un troisième mode de réalisation de l'invention. La figure 1 représente un équipement récepteur 2 d'un système de transmission OFDM, comme par exemple un terminal mobile d'un système radio mobile professionnel PMR (« Professional Mobile Radio ») de type TEDS. Le récepteur 2 comprend une antenne 4 pour la réception d'un signal OFDM transmis depuis un émetteur 6 à travers un canal de propagation radioélectrique 8. Le signal OFDM émis par l'émetteur 6 est organisé en trames de symboles répartis en temps et en fréquence parmi lesquels certains symboles, dits symboles pilotes, sont connus du récepteur 2 et sont stockés dans une mémoire 10 dudit récepteur 2 . Chaque trame comprend ainsi n symboles avec nf sous-porteuses et nt temps-symboles, n étant égal au produit de met de nt. Le récepteur 2 comprend des moyens de conversion temps-fréquence 12 du signal reçu pour le passer du domaine temporel au domaine fréquentiel. Ces moyens de conversion temps-fréquence 12 mettent en oeuvre une transformation de Fourier rapide FFT (« Fast Fourier Transform »). Le récepteur 2 comprend également un module d'estimation de canal 14 utilisant les symboles pilotes stockés dans la mémoire 10 et sa connaissance des contraintes physiques du canal 8 dans les domaines temporel et fréquentiel pour estimer le canal de propagation 8. D'une part, dans le domaine fréquentiel, l'étalement en fréquence du canal 8, nommé étalement Doppler, dû aux réflexions sur des obstacles proches, est limité. Cet étalement est compris entre -FD et +FD, où FD est la fréquence Doppler maximale donnée par la relation FD - .F dans laquelle est la vitesse du D c pl récepteur 2, c est la vitesse de la lumière, et Fp est la fréquence porteuse. Les composantes du spectre de fréquence du canal de propagation 8 selon l'axe des 20 fréquences sont donc comprises entre ces limites -FD et +FD. D'autre part, l'étalement temporel du canal 8, dû aux réflexions sur les obstacles lointains, est limité. Cet étalement temporel dépend de la bande de fréquences utilisée et de l'environnement. A titre d'exemples, à une fréquence porteuse de 400 MHz, en environnement urbain, l'étalement temporel est de 25 l'ordre de 5 ps alors qu'en environnement montagneux, cet étalement est de l'ordre de 15 ps. Les composantes de la réponse temporelle du canal 8 sont donc comprises entre des limites fixes pour des conditions d'environnement données. Les limites du spectre de fréquences et de la réponse temporelle du canal 8 sont connues du récepteur 2 et sont stockées dans la mémoire 10. 30 Le module d'estimation de canal 14 comprend des moyens de calcul 16 de matrices MF et MT de covariance du canal respectivement dans le domaine fréquentiel et dans le domaine temporel à partir respectivement des limites de l'étalement temporel et de l'étalement fréquentiel stockées dans la mémoire 10. Le module d'estimation de canal 14 comprend également des moyens de décomposition 18 des matrices MF et MT de covariance en vecteurs propres et en valeurs propres. Ces moyens sont aptes à calculer une matrice de covariance M globale du canal à partir des matrices MF et MT . Le module d'estimation de canal 14 comprend en outre des moyens d'estimation 20 du canal de propagation à l'aide des symboles pilotes stockés dans la mémoire 10 et des vecteurs propres et des valeurs propres obtenus par les moyens de décomposition 18 en utilisant le critère du maximum a posteriori ou de la minimisation de l'erreur quadratique. Le récepteur 2 comprend également un décodeur 22 fournissant une estimation des symboles émis à partir de l'estimation de canal effectuée par le module d'estimation de canal 14. La structure du récepteur de l'invention ayant été décrite, le fonctionnement du procédé d'estimation de canal mis en oeuvre dans le module d'estimation de canal 14 va être détaillé. De manière générale, le signal reçu par le récepteur 2 s'écrit sous la forme R=E.0+B où R est un vecteur de dimension n constitué des symboles reçus, C est un vecteur de dimension n représentant le canal de propagation 8 de dimension n, E est une matrice diagonale de dimension (n,n) constituée des symboles émis dans une trame et B est un vecteur de dimension n représentant le bruit du canal 8. On considère conventionnellement que le bruit du canal est un bruit Gaussien de variance, ou puissance, 62. La probabilité de recevoir le vecteur R R-EC 2 connaissant le vecteur C est alors égale à e 21:7 2 En outre, pour un canal radiomobile terrestre variant selon une loi de Rayleigh, _CH .M-1.0 la probabilité du canal est égale à e 2 où M est la matrice de covariance du canal et où la notation X" indique qu'il s'agit d'une matrice X conjuguée et 30 transposée. En ce qui concerne le choix de 62 et de M, il faut noter que l'on peut prendre pour M la matrice de covariance du canal normalisée, c'est-à-dire le canal de puissance moyenne unitaire. Dans ce cas, 62 représente alors l'inverse du rapport signal sur bruit. Il est possible de se fixer une valeur de rapport signal sur bruit cible, cette valeur ne changeant pas quelle que soit la valeur effective de la puissance de bruit et de la puissance de signal utile. D'autres solutions sont possibles comme, par exemple, d'estimer le rapport signal sur bruit au fur et à mesure de la réception et de la démodulation des informations reçues, par exemple sur les symboles pilotes qui sont connus de l'émetteur et du récepteur; Dans ce cas, il est possible d'adapter dynamiquement le rapport signal sur bruit dans le processus d'estimation du canal. L'organigramme de la figure 2 illustre un premier mode de réalisation du procédé d'estimation de canal de l'invention, mettant en oeuvre le critère du maximum a posteriori ou MAP. La résolution du problème d'estimation du canal au sens du MAP revient à maximiser la probabilité que le canal soit égal à C, connaissant le vecteur R, tout en répondant à ses contraintes physiques exprimées dans la matrice de covariance M. R-EC 2 _CH.N1-1.0 Cela revient à maximiser la probabilité e 3:72 .e 2 , ce qui revient à minimiser l'opposé du logarithme de cette expression égal à R-ECO2 +CH.M-1.0 a 2 La solution C est donc celle pour laquelle l'expression suivante, notée [1], obtenue en annulant le gradient par rapport à C de l'expression précédente, est vérifiée : (EH.E+152.M-1).0 =EH.R [1]. Seuls les symboles pilotes étant connus du récepteur 2, des zéros sont placés dans la matrice E aux emplacements des autres symboles inconnus du récepteur. Or, les contraintes du canal 8 dans le domaine temporel sont indépendantes de ses contraintes dans le domaine fréquentiel. La matrice de covariance M du canal est donc séparable en temps et en fréquence. La présente invention tire profit de cette séparabilité de la matrice de covariance M en l'exprimant comme le 2 9 83666 11 produit de Kronecker de deux matrices de covariance MF et MT , la matrice MF exprimant les contraintes du canal dans le domaine fréquentiel et la matrice MT exprimant les contraintes du canal dans le domaine temporel. Ainsi, à l'étape 30, les matrices MF et MT de covariance du canal 5 respectivement dans le domaine fréquentiel et dans le domaine temporel sont calculées de manière conventionnelle par les moyens de calcul 16 en utilisant respectivement les limites d'étalement temporel et Doppler du canal 8 stockées dans la mémoire 10. Ces matrices MF et MT sont stockées dans la mémoire 10. A l'étape 32, les matrices MF et MT de covariance sont décomposées en 10 vecteurs propres et en valeurs propres par les moyens de décomposition 18 selon les relations MF = WFHNFWF et MT = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matriceMF ; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MF; 15 - WT est une matrice de vecteurs propres de la matrice MT ; et - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MT . A l'étape 34, les moyens de décomposition 18 calculent le produit de Kronecker des matrices de vecteurs propres WF et WT pour obtenir une matrice 20 de vecteurs propres globale W de dimension (n, n). Ils calculent également une matrice diagonale de valeurs propres globale N de dimension (n, n) contenant les produits des valeurs propres des matrices NF et NT . Les matrices WF , NF, WT et NT sont également stockées dans la mémoire 10. La matrice de covariance globale M du canal 8 est alors égale à M=WHNW 25 L'expression [1] peut donc s'écrire : WHE+a . 2 W).C=EH.R En posant C=WH.,N.b où b est un vecteur de dimension n, et où est une matrice diagonale dont chaque terme de la diagonale est la racine carrée du terme correspondant de la diagonale de la matrice diagonale N , l'expression [1] s'écrit : + a 2413 = [2] 30 où I est la matrice identité. A l'étape 36, les moyens d'estimation de canal 20 calculent la matrice P =,N.W.EH.E.WH. et décomposent cette matrice P en vecteurs propres et valeurs propres selon la relation P = XHQX dans laquelle X est une matrice de vecteurs propres de la matrice P et Q est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice P. A l'étape 38, les moyens d'estimation de canal 20 estiment le canal 8 à partir de la relation 2 selon l'expression [3] suivante : C = WH ...N.XH .(Q +a 2.1)-1.X.,N.W.EH.R L'organigramme de la figure 3 illustre un deuxième mode de réalisation du procédé d'estimation de canal de l'invention, mettant également en oeuvre le critère MAP. Sur cette figure, les étapes 30 à 34 sont identiques à celles de la figure 2. La description de ces étapes n'est donc pas répétée. A l'étape 40, les moyens d'estimation de canal 20 extraient une matrice diagonale N' à partir de la matrice N , la matrice N' comprenant un nombre déterminé n' des valeurs propres les plus grandes de N et déterminent une matrice W' comprenant les vecteurs propres associés à ces n' valeurs propres les plus grandes. Les matrices N' et W' sont stockées dans la mémoire 10. Généralement, les valeurs propres de la matrice N décroissent rapidement. Ainsi, la somme des valeurs propres conservées dans la matrice N' représente une large part de l'énergie du canal 8. La matrice W' comprend ainsi les vecteurs propres les plus représentatifs pour exprimer les contraintes du canal. L'expression [1] peut donc s'écrire : (EHE+62.W'H.N'-1.W')C =EH.R En posant C =W'H.,/lVT'.b où b est un vecteur de dimension n', l'expression [1] s'écrit 252.1)13-1/1VP.W'.EH.R [4] où I est la matrice identité. A l'étape 42, les moyens d'estimation de canal 20 calculent une matrice P' selon la relation r= et décomposent cette matrice P' en vecteurs propres et valeurs propres selon la relation P'=X'HQ'X' dans laquelle X' 30 est une matrice de vecteurs propres de la matrice P' et Q' est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice P'. A l'étape 44, les moyens d'estimation de canal 20 estiment le canal 8 à partir de la relation [4] selon l'expression [5] suivante C .X'H .(Q'±ci 2 .I)-1.X'..N/1VT' [5]. La matriceW'étant de dimension (n',n) et la matrice N' étant de dimension (n',n'), ce mode de réalisation permet des calculs moins complexes que le premier mode de réalisation. Dans une réalisation préférée, les n' vecteurs propres de la matrice W' sont les produits de Kronecker d'un nombre limité n'f de vecteurs propres de la matrice WF et d'un nombre limité n't de vecteurs propres de la matrice WT et les n' valeurs propres de la matrice N' sont les produits des valeurs propres de la matrice NF correspondant aux n'f vecteurs propres de la matrice WF et des valeurs propres de la matrice NT correspondant aux n't vecteurs propres de la matrice WT Dans ce cas, la matrice W' est séparable en temps et en fréquence. On peut alors la mémoriser sous la forme du produit de Kronecker de 2 matrices et elle correspond alors à une dimension totale de [(n'f,nf)+(n't,nt)], ce qui permet de réduire la complexité du procédé. L'organigramme de la figure 4 illustre un troisième mode de réalisation du procédé d'estimation de canal de l'invention, mettant en oeuvre le critère de minimisation de l'erreur quadratique. Il s'agit, selon ce mode de réalisation, de minimiser l'erreur quadratique e_Ecr, en imposant que le canal soit une combinaison linéaire d'un nombre limité de vecteurs propres de la matrice de covariance M. On utilise ainsi dans ce cas les matrices W'et N' telles que déterminées à l'étape 40 du deuxième mode de réalisation. Ainsi, sur la figure 4, les étapes 30 à 40 sont identiques à celles de la figure 3. La description de ces étapes n'est donc pas répétée. A l'étape 50, les moyens d'estimation de canal 20 estiment le canal 8 selon la 30 relation C = WH.N/.(N/ +a 2.I) .E.W/H .R [6] où la notation X" indique qu'il s'agit d'une matrice X conjuguée et transposée. Cette expression permet d'obtenir des résultats proches de ceux de l'expression [5] du deuxième mode de réalisation. Lorsque la matrice W' est séparable en temps et en fréquence, la complexité de ce mode de réalisation est équivalente à celle du deuxième mode de réalisation. Bien entendu, d'autres modes de réalisation peuvent encore être envisagés. Plus particulièrement, le calcul des matrices de covariance peut être réalisé de manière dynamique afin de prendre en compte les variations des contraintes du canal. Il est également possible de prendre en compte d'autres paramètres dans le calcul des matrices de covariance notamment le profil des retards, etc. REVENDICATIONS1. Procédé d'estimation d'un canal de propagation radioélectrique (8) entre un émetteur (6) et un récepteur (2), ledit émetteur (6) émettant un signal comprenant des trames de symboles répartis en temps et en fréquence parmi lesquels certains symboles, dits symboles pilotes, sont connus du récepteur (2), caractérisé par le fait qu'il comporte les étapes mises en oeuvre dans le récepteur (2) de : - calcul (30) d'une matrice MF de covariance du canal (8) dans le domaine fréquentiel; - calcul (30) d'une matrice MT de covariance du canal (8) dans le domaine temporel; - décomposition (32) des matrices MF et MT de covariance en vecteurs propres selon les relations MF = WFHNFWF et MT = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matrice MF ; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MF; - WT est une matrice de vecteurs propres de la matrice MT ; - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MT ; - calcul (34) du produit de Kronecker des matrices de vecteurs propres WF et WT pour obtenir une matrice de vecteurs propres globale W ; - calcul (34) d'une matrice diagonale de valeurs propres globale N égale au produit de Kronecker des matrices de valeurs propres NF et NT ; et - estimation du canal de propagation (8) à l'aide des symboles pilotes et des matrices de vecteurs propres W et de valeurs propres N en utilisant le critère du maximum a posteriori ou de la minimisation de l'erreur quadratique. 2. Procédé selon la revendication 1, dans lequel le critère utilisé est le critère du maximum a posteriori et l'étape d'estimation du canal comprend les sous-étapes de :- calcul (36) d'une matrice P selon la relation P = dans laquelle E est une matrice diagonale comprenant les symboles pilotes aux positions des symboles pilotes, les autres éléments de la diagonale étant nuls ; - décomposition (36) de la matrice P en vecteurs propres selon la relation P = XHQX dans laquelle X est une matrice de vecteurs propres de la matrice P et Q est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice P; et - estimation (38) du canal (8) selon la relation C = WH .,N.X1/.(Q + a2.I) .R dans laquelle C est un vecteur représentant le canal, I est la matrice identité, R est un vecteur comprenant les symboles reçus et a 2 représente la puissance du bruit dans le canal. 3. Procédé selon la revendication 1, dans lequel l'étape d'estimation du canal comprend les sous étapes de : - extraction (40) d'une matrice W' à partir de la matrice W et d'une matrice N' à partir de la matrice N , la matrice N' comprenant un nombre déterminé n' de valeurs propres de la matrice N et la matrice W' comprenant les vecteurs propres associés à ces n' valeurs propres; et - estimation (44 ; 50) du canal de propagation (8) à l'aide des symboles 20 pilotes et des matrices de vecteurs propres W' et de valeurs propres N' en utilisant le critère du maximum a posteriori ou de la minimisation de l'erreur quadratique. 4. Procédé selon la revendication 3, dans lequel les n' valeurs propres sont les 25 plus grandes valeurs propres de la matrice N. 5. Procédé selon la revendication 3 ou 4, dans lequel les n' vecteurs propres de la matrice W' sont les produits de Kronecker d'un nombre limité n'f de vecteurs propres de la matrice WF et d'un nombre limité n't de vecteurs propres de la 30 matrice WT et les n' valeurs propres de la matrice N' sont les produits des valeurs propres de la matrice NF correspondant aux n'f vecteurs propres de la matriceWF et des valeurs propres de la matrice NTcorrespondant au n't vecteurs propres de la matrice WT 6. Procédé selon l'une des revendications 3 à 5, dans lequel n' est inférieur ou 5 égal au nombre de symboles pilotes dans chaque trame. 7. Procédé selon l'une des revendications 3 à 6, dans lequel le critère utilisé est le critère du maximum a posteriori et l'étape d'estimation du canal comprend les sous-étapes de : 10 - calcul (42) d'une matrice P' selon la relation r =N' dans laquelle E est une matrice diagonale comprenant les symboles pilotes aux positions des symboles pilotes, les autres éléments de la diagonale étant nuls ; - décomposition (42) de la matrice P' en vecteurs propres selon la 15 relation P/ = X/HQ/X/ dans laquelle X' est une matrice de vecteurs propres de la matrice P' et Q' est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice P'; et - estimation (44) du canal selon la relation C = W/H.-N/W.X'11.(Q/-Fa2.1)1.X'.1T'.W'.E11.R dans laquelle C est un vecteur 20 représentant le canal, I est la matrice identité, R est un vecteur comprenant les symboles reçus et a 2 représente la puissance du bruit dans le canal. 8. Procédé selon la revendication 6, dans lequel le critère utilisé est le critère de la minimisation de l'erreur quadratique et le canal est estimé (50) selon la relation 25 C = WH.N'.(N/-E62.1)1.(W'.EH.E.WH)1W'.E11.R dans laquelle C est un vecteur représentant le canal, I est la matrice identité, R est un vecteur comprenant les symboles reçus et a 2 représente la puissance du bruit dans le canal. 9. Procédé selon l'une des revendications précédentes, dans lequel le signal est 30 un signal multiporteuses, notamment un signal OFDM. 10. Equipement récepteur (2) apte à recevoir un signal transmis à travers un canal de propagation radioélectrique (8), ledit signal comprenant des trames de symboles répartis en temps et en fréquence parmi lesquels certains symboles, dits symboles pilotes, sont connus dudit récepteur (2), caractérisé par le fait qu'il comporte des moyens de : - calcul (16) d'une matrice MF de covariance du canal (8) dans le domaine fréquentiel; - calcul (16) d'une matrice MT de covariance du canal (8) dans le domaine temporel; - décomposition (18) des matrices MF et MT de covariance en vecteurs propres selon les relations MF = WFHNFWF et MT = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matriceMF ; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MF; - WT est une matrice de vecteurs propres de la matrice MT ; - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice MT ; - calcul (18) du produit de Kronecker des matrices de vecteurs propres WF et WT pour obtenir une matrice de vecteurs propres globale W ; - calcul (18) d'une matrice diagonale de valeurs propres globale N égale au produit de Kronecker des matrices de valeurs propres NF et NT ; et - estimation (20) du canal de propagation à l'aide des symboles pilotes et des matrices de vecteurs propres W et de valeurs propres N en utilisant le critère du maximum a posteriori ou de la minimisation de l'erreur quadratique. 11. Programme d'ordinateur comportant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 9 lorsque le programme est exécuté par au moins un processeur.