FR 3022419 A1 20151218 FR 1401367 A 20140617 PROCEDE D'ESTIMATION D'UN CANAL DE PROPAGATION RADIOELECTRIQUE DOMAINE TECHNIQUE 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 é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 OFDM (« Orthogonal Frequency Division Multiplexing >)) ou de type SC-FDMA (Single-Carrier / Frequency Division Multiple Access) ou encore de type FB-MC (Filter Bank Multi-Carrier), ces techniques étant toutes des techniques de transmission de type multiporteuses. ETAT DE L'ART 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 20 (« TETRA Enhanced Data Service »), DAB (« Digital Audio Broadcasting ») et DVB-T (« Digital Video Broadcasting-Terrestrial »). Dans les systèmes de transmission multiporteuses, les données sont en général codées sous forme de symboles organisés en trames envoyées dans un 25 signal par un émetteur à destination d'un récepteur à travers un canal de propagation, chaque trame utilisant Nt sous-porteuses fréquentielles sur chacune desquelles sont émis Nt symboles, soit n nombre total de symboles dans le plan temps-fréquence égal à n = Nt X Nf. 30 Le canal de propagation est défini comme étant le chemin radioélectrique qu'emprunte le signal émis depuis l'émetteur vers le récepteur. De manière générale, le signal reçu par le récepteur s'écrit sous la forme R= E.0 +B où R est un vecteur colonne de dimension n constitué des symboles reçus, E est une matrice diagonale de dimension (n,n) dont la diagonale est constituée des n symboles émis dans une trame, C est un vecteur colonne, dit vecteur canal, de dimension n représentant le canal de propagation et B est un vecteur colonne de dimension n représentant le bruit du canal. On représentera conventionnellement dans la suite les vecteurs R, C et B sous la forme de la concaténation de Nt vecteurs, chacun étant de dimension Nf, le premier vecteur correspondant aux Nf symboles, correspondant au premier symbole de la trame, successivement situées sur Nf sous-porteuses, et les vecteurs suivants correspondants de façon similaire aux Nt symboles successifs. De même, dans la description qui suit, la diagonale de la matrice E est constituée par la concaténation, pour les Nt symboles successifs dans la trame, des Nf valeurs successives du canal de propagation pour les sous-porteuses successives. Le canal de propagation comprend, pour un signal émis par l'émetteur, une pluralité de trajets entre l'émetteur et le récepteur dus notamment à la réflexion et/ou à la diffraction du signal sur des obstacles tels que, par exemple, des collines, des bâtiments, des véhicules etc. Chaque trajet est ainsi caractérisé par un retard et une atténuation et le délai entre le signal reçu par le récepteur sur le trajet radioélectrique le plus court et le signal reçu par le récepteur sur le trajet radioélectrique le plus long constitue l'étalement temporel du canal. L'ensemble des trajets avec leurs puissances et leurs retards constituent le profil temporel du canal. Selon les caractéristiques de l'environnement de propagation, un profil temporel typique de propagation est attendu. Les différents trajets empruntés par le signal peuvent en outre se combiner entre eux de façon destructive, notamment lorsque les réflexions du signal ont lieu sur des obstacles situés à proximité du récepteur et/ou de l'émetteur. Plus précisément, chacun des trajets peut être affecté par un évanouissement ou fading dit « de Rayleigh », connu de l'homme du métier, dont la pseudo-fréquence est égale au double de la fréquence Doppler, cette fréquence Doppler étant déterminée à partir de la vitesse de l'émetteur et/ou du récepteur lorsqu'au moins l'un des deux est mobile et de la fréquence porteuse du signal transmis. La différence entre la fréquence la plus faible et la plus élevée parmi les signaux reçus par le récepteur sur les différents trajets radioélectrique d'un signal monochrome émis par l'émetteur constitue l'étalement fréquentiel du canal. Le spectre de puissance de ce signal reçu est appelé profil fréquentiel du canal. Selon les caractéristiques de l'environnement de propagation, un profil fréquentiel typique de propagation est attendu. Un canal à multi-trajets avec fading de Rayleigh peut donc être sélectif en temps et/ou en fréquence. Afin d'estimer le canal de propagation, certains symboles, répartis dans le plan temps-fréquence, sont insérés dans chaque trame 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 canal de propagation. Dans un système multiporteuses, l'effet du canal de propagation sur un symbole de données, dans le plan temps-fréquence, est généralement modélisé simplement par un coefficient multiplicatif complexe que le récepteur cherche à estimer afin de déterminer le symbole émis avec le minimum possible d'erreur. L'estimation de canal est définie comme étant la détermination de ces coefficients pour tous les symboles émis. De manière connue, le récepteur commence, en général, par estimer le canal à la position des symboles pilotes puis il effectue une estimation du canal sur l'ensemble de la trame. L'estimation du canal sur les positions ne correspondant pas à des symboles pilotes est en général effectuée au moyen d'une interpolation à partir de l'estimation du canal sur les symboles pilotes. Le document FR2814011 décrit un procédé d'estimation de canal permettant 30 de prendre en compte les caractéristiques physiques connues du canal de propagation. Cependant, le procédé décrit est très contraignant en termes de complexité. Afin de résoudre en partie ces inconvénients, la demande de brevet FR2983666 Al déposée par la Demanderesse décrit un procédé d'estimation de canal de complexité réduite permettant de prendre en compte les caractéristiques physiques connues du canal de propagation et mettant notamment 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 à iR-EC112 CH .M1-1.0 maximiser le terme e 20-2 .e 2 , c'est-à-dire à minimiser l'opposé du IIR logarithme de cette expression égal à +CH.M -1C tout en répondant 1 ' aux variations physiques du canal en temps et en fréquence, exprimées au moyen de la matrice de covariance M1 globale du canal. 0-2 représente la variance du bruit sur le canal. La matrice de covariance M1 globale du canal est obtenue à partir de la connaissance des variations physiques du canal, à savoir : - le profil temporel du canal, défini dans une fenêtre temporelle, qui entraine des variations du canal de propagation le long de l'axe des fréquences et à partir duquel peut être calculée une matrice de covariance Mn exprimant donc les variations fréquentielles du canal dues audit profil temporel, et - le profil fréquentiel, défini dans une fenêtre fréquentielle, qui entraine des variations du canal de propagation le long de l'axe du temps et à partir duquel peut être calculée une matrice de covariance M1t exprimant donc les variations temporelles du canal dues audit profil fréquentiel. La matrice de covariance temporelle M1t exprimant la contrainte temporelle du canal peut être obtenue par Transformée de Fourier inverse du profil fréquentiel du canal tandis que la matrice de covariance fréquentielle Mif exprimant la contrainte fréquentielle du canal peut être obtenue par Transformée de Fourier du profil temporel du canal. La matrice de covariance M1 globale du canal est le produit tensoriel de Kronecker des deux matrices de covariance temporelle Mit et fréquentielle Mif : M1 = Mit Mlf Souvent, l'hypothèse théorique considérée pour estimer le canal dans le domaine fréquentiel correspond à un profil temporel du canal qui est symétrique et centré sur une position dite « de synchronisation temporelle du signal ». Une telle position de synchronisation temporelle du signal est généralement prédéterminée de manière connue, par exemple, au moyen d'une séquence de synchronisation temporelle spécifique (ou par d'autres moyens), et est utilisée par le récepteur pour synchroniser temporellement le signal reçu de l'émetteur, c'est-à-dire pour déterminer un instant précis associé à la fenêtre du profil temporel du signal (dans ce cas le milieu de la fenêtre). Dans ce cas, les valeurs propres et les vecteurs propres de la matrice de covariance fréquentielle Mif sont réels. Or, en pratique, le profil temporel du canal n'est pas nécessairement symétrique et centré sur la position de synchronisation temporelle. Il peut être, par exemple, symétrique mais centré sur une position centrale effective temporelle différente de la position dite "de synchronisation temporelle" du signal car la synchronisation temporelle est souvent réalisée sur les premiers trajets de propagation des signaux reçus, qui sont en général les moins atténués, alors que des trajets ultérieurs peuvent être plus atténués et très retardés mais encore significatifs. Dans ce cas, les valeurs propres de la matrice de covariance fréquentielle Mif sont réelles mais les vecteurs propres de la matrice de covariance fréquentielle Mif sont complexes. De même, la contrainte théorique considérée dans le domaine temporel correspond à un profil fréquentiel du canal qui est symétrique et centré sur une position dite « de synchronisation fréquentielle du signal ». Une telle position de synchronisation fréquentielle du signal peut être prédéterminée de manière connue, par exemple, au moyen d'une séquence de synchronisation fréquentielle spécifique (ou par d'autres moyens), et est utilisée par le récepteur pour synchroniser fréquentiellement le signal reçu de l'émetteur, c'est-à-dire pour déterminer une fréquence précise associée à la fenêtre d'étalement du profil fréquentiel du signal (dans ce cas le milieu de la fenêtre). Dans ce cas, les valeurs propres et les vecteurs propres de la matrice de covariance Mit sont réels. Or, en pratique, le profil fréquentiel du canal n'est pas nécessairement symétrique et centré sur la position de synchronisation fréquentielle. Il peut être, par exemple, symétrique mais centré sur une position centrale effective fréquentielle différente de la position dite "de synchronisation fréquentielle" du signal. Dans ce cas, les valeurs propres de la matrice de covariance M1t sont réelles mais les vecteurs propres de la matrice de covariance M1t sont complexes. La résolution du problème d'estimation du canal au sens du MAP avec des vecteurs propres complexes dans les matrices de covariance temporelle Mif et/ou fréquentielle 11.11t rend le récepteur complexe et coûteux car sa mémoire et la capacité de son processeur doivent être importantes afin de traiter les données reçues pour estimer le canal de propagation, ce qui présente un inconvénient important. PRESENTATION GENERALE DE L'INVENTION La présente invention vise à résoudre cet inconvénient en proposant un procédé d'estimation d'un canal de propagation radioélectrique entre un émetteur et un récepteur qui permette de réduire la complexité et le coût d'implémentation du récepteur tout en permettant l'estimation d'un canal dont les contraintes correspondent à celles de la réalité. 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 utilisant chacune Nf sous-porteuses fréquentielles sur chacune desquelles sont émis Nt symboles, où parmi l'ensemble des symboles, certains symboles, dits symboles pilotes, sont connus dudit récepteur, ledit signal étant synchronisé par le récepteur à partir d'une position de synchronisation temporelle et d'une position de synchronisation fréquentielle, le profil temporel du canal étant symétrique et centré sur une position centrale effective temporelle T -moyen distincte de la position de synchronisation temporelle, le profil fréquentiel du canal étant symétrique et centré sur une position centrale effective fréquentielle f -moyen distincte de la position de synchronisation fréquentielle. Le procédé, mis en oeuvre par le récepteur, est remarquable par le fait qu'il comprend les étapes de : - détermination d'une matrice de covariance fréquentielle intermédiaire M2f représentant les variations fréquentielles du canal dues à un profil temporel théorique du canal, symétrique et centré sur la position de synchronisation temporelle, ce profil temporel théorique du canal étant obtenu à partir du profil temporel du canal translaté en temps de (tmoyen) de façon à être symétrique autour de la position de synchronisation temporelle (position zéro temporelle); - détermination d'une matrice de covariance temporelle intermédiaire M2t représentant les variations temporelles du canal dues à un profil fréquentiel théorique du canal, symétrique et centré sur la position de synchronisation fréquentielle, ce profil fréquentiel théorique du canal étant obtenu à partir du profil fréquentiel du canal translaté en fréquence de (-fmoyen) de façon à être symétrique autour de la position de synchronisation fréquentielle (position zéro fréquentielle) ; - calcul d'une matrice de covariance globale intermédiaire M2 selon le produit de Kronecker: M., = M 2tC)M2f ; - calcul d'une matrice E' selon la relation E'. E.TH , où T = 1 ®Tf et dans laquelle E est une matrice diagonale comprenant les symboles pilotes émis dans une 20 trame aux positions des symboles pilotes, les autres symboles étant nuls, Tf est une matrice diagonale (Ni, Nf) de centrage définie par 0 rm où ço= 271-. °Ven et Tt est une matrice diagonale (Ni, Nf T e12" o 0 O eorf-'),° Nt) de centrage définie par Tt 1 0 0 où f nt°Yen ; 0 et 1^1, ei2w O 0 O e 1( Nt-1)w - calcul d'un vecteur colonne représentant un canal intermédiaire C' qui minimise 1,2 la relation R-E.0+ C'H. Mil. C' où R est un vecteur colonne de dimension (Nt x 02 Nf) constitué des symboles reçus, et ; - estimation du canal de propagation selon la relation C = T .C', à partir du vecteur C' calculé. Par les termes « détermination d'une matrice de covariance fréquentielle intermédiaire M2f représentant les variations fréquentielles du canal dues au profil temporel du canal centré sur la position de synchronisation temporelle », on entend qu'on détermine les effets du profil temporel du signal sur les variations fréquentielles du signal lorsque ce profil temporel est centré sur la position de synchronisation temporelle du signal. De même, par les termes « détermination d'une matrice de covariance temporelle intermédiaire M2t représentant les variations temporelles du canal, dues au profil fréquentiel du canal centré sur la position de synchronisation fréquentielle », on entend qu'on détermine les effets du profil fréquentiel du signal sur les variations temporels du signal lorsque ce profil fréquentiel est centré sur la position de synchronisation fréquentielle du signal. Le produit de Kronecker de deux matrices A et B, la matrice A ayant pour composantes (ai1), 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 : ab,13 A OB = - - an,'B Lorsque le profil temporel du canal est symétrique et centré autour de la position centrale effective temporelle tmoyen mais que ladite position centrale effective temporelle Tmoyen ne correspond pas à la position de synchronisation temporelle, la matrice de covariance fréquentielle Mit exprimant la variation du 30 canal le long de l'axe fréquentiel est une matrice dont les vecteurs propres sont complexes. Le procédé selon l'invention permet de définir, à partir du même profil temporel mais centré autour de la position de synchronisation temporelle, par exemple par transformation de Fourier du profil temporel du canal, une matrice de covariance fréquentielle intermédiaire M2f, qui correspond à la matrice de covariance fréquentielle Mif du canal réel mais qui a des vecteurs propres réels. En d'autres termes, la matrice de covariance fréquentielle intermédiaire M2f correspond à la contrainte due au profil temporel du canal mais avec un décalage temporel tel que la contrainte en terme de profil temporel du canal soit symétrique et centrée non pas sur la position centrale effective temporelle Crnoyen mais sur la position de synchronisation temporelle. La contrainte de profil temporel du canal est alors centrée et symétrique autour de la position de synchronisation temporelle. En outre, dans le cas particulier où la position centrale effective temporelle tmoyen correspond à la position de synchronisation temporelle, alors M2f = Mlf. De même, lorsque le profil fréquentiel du canal est symétrique autour de la position centrale effective fréquentielle f -moyen mais que ladite position centrale effective fréquentielle f -moyen ne correspond pas à la position de synchronisation fréquentielle, la matrice de covariance temporelle M1t exprimant la variation du canal le long de l'axe temporel est une matrice dont les vecteurs propres sont complexes. Le procédé selon l'invention permet de définir, à partir du même profil fréquentiel mais centré maintenant autour de la position de synchronisation fréquentielle, par exemple par transformation de Fourier inverse du profil fréquentiel du canal, une matrice de covariance temporelle intermédiaire M2t, qui correspond à la matrice de covariance temporelle M1t du canal réel mais qui a des vecteurs propres réels. En d'autres termes, la matrice de covariance temporelle intermédiaire M2t correspond à la contrainte due au profil fréquentiel du canal mais avec un décalage fréquentiel tel que la contrainte en terme de profil fréquentiel du canal soit symétrique et centrée non pas sur la position centrale effective fréquentielle f -moyen mais sur la position de synchronisation fréquentielle. La contrainte en terme de profil fréquentiel du canal est alors centrée et symétrique autour de la position de synchronisation fréquentielle. En outre, dans le cas 302 2 4 19 10 particulier où la position centrale effective fréquentielle fmoyen correspond à la position de synchronisation fréquentielle, alors M2t = M1t. La matrice de covariance intermédiaire globale M2 qui est le produit de 5 Kronecker de la matrice de covariance temporelle intermédiaire M2t et de la matrice de covariance fréquentielle intermédiaire M2f correspond ainsi à un canal virtuel de propagation C' et possède avantageusement des vecteurs propres réels. On peut utiliser, pour estimer le canal C' correspondant à cette matrice intermédiaire M2, une méthode basée sur le critère du maximum a posteriori 10 (MAP) telle que, par exemple, la méthode décrite dans le document FR2983666 Al ou bien dans le document FR2814011. Le vecteur colonne C correspondant au canal réel est ensuite avantageusement obtenu à partir de la matrice de canal intermédiaire C' et de la matrice de centrage T définie selon l'invention. 15 L'exploitation de la séparabilité de la matrice de covariance intermédiaire M2, dont les vecteurs propres sont réels, est très favorable en termes de complexité. Ainsi, pour une trame de n symboles, avec Nf sous-porteuses et Nt temps-symboles, la matrice de covariance intermédiaire M2 est de dimension (n,n) avec n = Nt x N. Cette matrice est égale au produit de Kronecker d'une matrice de 20 covariance M2t de dimension (Nf, Nf) avec une matrice de covariance M2f de dimension (Nt, Nt). Il suffit donc de mémoriser (Nf2 + Nt2) valeurs au lieu de (Nf2 x Nt2) valeurs. 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, en tenant compte des contraintes 25 physiques réelles du canal de propagation. De préférence, sur la diagonale de la matrice E, les symboles émis sont concaténés de façon à ce que, pour le premier symbole de la trame, on place d'abord les symboles correspondant aux sous-porteuses successives, puis on fait 30 de même pour le deuxième symbole de la trame, et ainsi de suite jusqu'au dernier symbole de la trame. Selon un aspect de l'invention, le procédé comprend une étape préliminaire de détermination de la position de synchronisation temporelle et de la position de synchronisation fréquentielle du signal. Selon un autre aspect de l'invention, le procédé comprend une étape de détermination de la position centrale effective temporelle tmoyen et de la position centrale effective fréquentielle f -moyen du signal. De préférence, l'étape de calcul du vecteur canal intermédiaire C' comprend 10 les sous-étapes de : - décomposition des matrices M2f et M2t de covariance en vecteurs propres selon les relations M21 = WFHNFWF et M2, = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matrice M2f - NF est une matrice diagonale comprenant des valeurs propres 15 associées aux vecteurs propres de la matrice M2f WT est une matrice de vecteurs propres de la matrice M2t ; - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice M2t ; - calcul du produit de Kronecker des matrices de vecteurs propres WF et WT pour 20 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 - calcul du vecteur canal intermédiaire C' à l'aide des symboles pilotes donnés par la matrice E' et des matrices de vecteurs propres W et de valeurs propres N en 25 utilisant le critère du maximum a posteriori. 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 multiporteuses connus comme LTE et TEDS, et/ou les symboles de 30 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 matrice W des vecteurs propres étant le produit de Kronecker des deux matrices WF et WT de vecteurs propres, il suffit de seulement (Nt X Nf Nt x Nt) = Nt2+Nt2 valeurs pour la représenter. De même, la matrice diagonale N des valeurs propres étant le produit de Kronecker des deux matrices diagonales NF et NT, il suffit seulement de Nt-ENt valeurs pour la représenter. Cela permet 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, l'étape de calcul du vecteur canal intermédiaire C' 10 comprend les sous-étapes de : - calcul d'une matrice P selon la relation P = 1/1 ; - 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 15 propres de la matrice P; et - calcul du vecteur C' selon la relation C'= WH.N/11.XH.(Q + dans laquelle I est la matrice identité, R est un vecteur comprenant les symboles reçus et a' représente la puissance du bruit dans le canal. 20 Plus précisément, la matrice N/1V1 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 25 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 de calcul de la matrice C' comprend les 30 sous étapes de : - 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 - calcul de la matrice C' à 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. 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't de vecteurs propres de la matrice WF et d'un nombre limité n't de vecteurs propres de la matrice WF et les n' valeurs propres de la matrice N' sont les produits des valeurs propres de la matriceNF correspondant aux n't vecteurs propres de la matrice WF et des valeurs propres 20 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't x Nf n't x Nt) valeurs pour la représenter. Ainsi, la complexité du procédé est encore plus réduite. 25 De préférence, n' est inférieur ou égal au nombre de symboles pilotes dans chaque trame. Avantageusement, l'étape de calcul de la matrice C' comprend les sous-étapes 30 de : - calcul d'une matrice P' selon la relation P'=-Nr1^1-7.W.E'H.E.W'H.; - décomposition de la matrice P' en vecteurs propres selon la 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 - calcul de la matrice C' selon la relation C'= c72.1)1.X.W.E'H.R dans laquelle I est la matrice identité, R est un vecteur comprenant les symboles reçus et o-2 représentela puissance du bruit dans le canal. Plus précisément, la matrice ,FN-1' 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. Avantageusement, le signal est un signal multiporteuses, notamment un signal 20 OFDM. L'invention peut également s'appliquer à des systèmes monoporteuse. Dans ce cas, on a alors M2f = Mlf = [1], Où [1] représente la matrice identité de dimension (1,1), cette matrice est dégénérée en un scalaire égale à 1. 25 L'invention concerne aussi un équipement récepteur apte à recevoir d'un émetteur au moins un signal transmis à travers un canal de propagation radioélectrique, ledit signal comprenant des trames utilisant chacune Nt sous-porteuses fréquentielles sur chacune desquelles sont émis Nt symboles, où parmi 30 l'ensemble des symboles, certains symboles, dits symboles pilotes, sont connus dudit équipement récepteur, le signal étant synchronisé par le récepteur à partir d'une position de synchronisation temporelle et d'une position de synchronisation fréquentielle, le profil temporel du canal étant symétrique et centré autour d'une yen position centrale effective temporelle tmodistincte de la position de synchronisation temporelle, le profil fréquentiel du canal étant symétrique et centré autour d'une position centrale effective fréquentielle fmoyen distincte de la position de synchronisation fréquentielle. Le procédé, mis en oeuvre dans le récepteur, est remarquable par le fait qu'il comprend : - des moyens de détermination d'une matrice de covariance fréquentielle intermédiaire M2f représentant les variations fréquentielles du canal, dues à un profil temporel théorique du canal, symétrique et centré sur la position de synchronisation temporelle, ce profil temporel théorique du canal étant obtenu à partir du profil temporel du canal translaté en temps de (tmoyen) de façon à être symétrique autour de la position de synchronisation temporelle ; - des moyens de détermination d'une matrice de covariance temporelle intermédiaire M2t représentant les variations temporelles du canal dues à un profil fréquentiel théorique du canal, symétrique et centré sur la position de synchronisation fréquentielle, ce profil fréquentiel théorique du canal étant obtenu à partir du profil fréquentiel du canal translaté en fréquence de (fmoyen) de façon à être symétrique autour de la position de synchronisation fréquentielle ; - des moyens de calcul d'une matrice de covariance globale intermédiaire M2 20 selon le produit de Kronecker: M2 =M2 °M2f ; - des moyens de calcul d'une matrice E' selon la relation E'. E.TH , où T =T,OTf et dans laquelle E est une matrice diagonale comprenant les symboles pilotes émis dans une trame aux positions des symboles pilotes, les autres symboles étant nuls, Tf est une matrice diagonale (Ne, Nf) de centrage définie par 1 0 0 el9 ei2v Tf = 0 0 où ço= 2z r. °Y en et T est une matrice diagonale (Nt, Nf o Nt) de centrage définie par T = 0 ei2v O O O eioh-ow - des moyens de calcul d'un vecteur colonne représentant un canal intermédiaire II42 C' qui minimise la relation R-E.01 2 C". M2 . Ci où R est un vecteur colonne de dimension (Nt x Nt) constitué des symboles reçus, et ; - des moyens d'estimation du canal de propagation selon la relation c = TH.C', à partir du vecteur C' calculé. L'invention concerne aussi un programme d'ordinateur comportant des 10 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 à 3 illustrent schématiquement le déroulement de ce programme d'ordinateur, pour des modes de réalisation 15 préférés de l'invention. DESCRIPTION DES FIGURES 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 : 20 - 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 mise en oeuvre de l'invention ; et où yt =-2,71-.fino-ven ; la figure 3 est un organigramme illustrant le fonctionnement du procédé d'estimation de canal selon un deuxième mode mise en oeuvre de l'invention. DESCRIPTION DETAILLEE DE L'INVENTION Equipement récepteur 2 La figure 1 représente un équipement récepteur 2 d'un système de transmission multiporteuses, comme par exemple un terminal mobile d'un système radio mobile professionnel PMR ( Les caractéristiques du profil fréquentiel et du profil temporel du canal 8 sont connues du récepteur 2 et sont stockées dans la mémoire 10. Le récepteur 2 est aussi configuré pour déterminer et/ou recevoir une position de synchronisation temporelle du signal et une position de synchronisation fréquentielle du signal. La position de synchronisation temporelle peut être prédéterminée de manière connue, par exemple, au moyen d'une séquence de synchronisation temporelle spécifique (ou par tout autre moyen). La position de synchronisation fréquentielle peut être prédéterminée de manière connue, par exemple, au moyen d'une séquence de synchronisation fréquentielle spécifique (ou par tout autre moyen). Le récepteur 2 est aussi configuré pour déterminer une position centrale effective temporelle du signal correspondant au milieu de la fenêtre d'étalement du profil temporel du signal reçu et une position centrale effective fréquentielle du signal correspondant au milieu de la fenêtre d'étalement du profil fréquentiel du signal reçu. lI. Module d'estimation de canal 14 a) Moyens de détermination 16 Le module d'estimation de canal 14 comprend tout d'abord des moyens de détermination 16 configurés pour déterminer : - la matrice de covariance fréquentielle intermédiaire M2f représentant les variations fréquentielles du canal dues au profil temporel du canal centré sur la position de synchronisation temporelle, et - la matrice de covariance temporelle intermédiaire M2e représentant les variations temporelles du canal, dues au profil fréquentiel du canal centré sur la position de synchronisation fréquentielle. La détermination de la matrice de covariance fréquentielle intermédiaire M2f peut être réalisée, par exemple, en appliquant une transformation de Fourier à la fenêtre d'étalement du profil temporel du signal centrée sur la position de synchronisation temporelle du signal. De même, la détermination de la matrice de covariance temporelle 20 intermédiaire M2e peut être réalisée, par exemple, en appliquant une transformation de Fourier inverse à la fenêtre d'étalement du profil fréquentiel du signal centrée sur la position de synchronisation fréquentielle du signal. De plus, les moyens de détermination 16 sont configurés pour calculer une matrice de covariance globale intermédiaire M2 selon le produit de Kronecker : 25 M2 = M2t °M2f - b) Moyens de calcul 18 Le module d'estimation de canal 14 comprend également des moyens de calcul 18 d'une matrice E' selon la relation E', E.TH , où T = T, 0Tf, , dans laquelle 30 E est une matrice diagonale comprenant les symboles pilotes émis dans une trame aux positions des symboles pilotes, les autres symboles étant nuls, Tf est une matrice diagonale (Ne, Ne) de centrage définie par 1 0 0 eie -12ço Tf = 0 O 0 1- où q)= 27r. MOT en et Tt est une matrice diagonale (Ni, Nf o 1 0 0 0 e` 0 0 ei (Nt-Die Nt) de centrage définie par 7; = où yf = -27r. fmo,e. . N, c) Moyens d'estimation 20 Le module d'estimation de canal 14 comprend en outre des moyens d'estimation 20 du canal de propagation 8 configurés pour calculer un vecteur Il R-E1 .CI 11 2 colonne C' qui minimise la relation 2 + CH. M. C' où R est un vecteur colonne de dimension (Nt x Nf) constitué des symboles reçus et estimer le canal selon la relation c = TH.c , à partir du vecteur C' calculé. III. Décodeur 22 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é. IV. Mise en oeuvre De manière générale, le signal reçu par le récepteur 2 s'écrit sous la forme R-E.C+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, 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 5 de variance, ou puissance, 0-2. La probabilité de recevoir le vecteur R R-E 2 connaissant le vecteur canal C est alors égale à e 2e2 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ù M1 est la matrice de covariance 10 globale représentant les contraintes physiques réelles en temps en en fréquence du canal 8 et où la notation XH indique qu'il s'agit d'une matrice X conjuguée et transposée. En ce qui concerne le choix de 02 et de M1, il faut noter que l'on peut prendre 15 pour M1 la matrice de covariance du canal normalisée, c'est-à-dire le canal de puissance moyenne unitaire. Dans ce cas, 02 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 20 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. 25 Le procédé selon l'invention met 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é qu'un canal intermédiaire virtuel soit égal à un vecteur canal intermédiaire C', connaissant le vecteur R, tout en répondant à ses 30 contraintes physiques exprimées dans une matrice de covariance intermédiaire globale notée M2 symétrique et centrée sur une position de synchronisation temporelle et sur une position de synchronisation fréquentielle. Cela revient à maximiser la probabilité e R-E'C' 2 CH .1V12-1.0 20.2 2 .e , ce qui revient à minimiser l'opposé du logarithme de cette expression égal à R - E'C' 2 + cull.m2-1.ce 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 : (E'H.E'+o-2.M = E" .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 M2 du canal est donc séparable en temps et en fréquence en l'exprimant comme le produit de Kronecker de deux matrices de covariance M2f et M2t, la matrice M2f exprimant les contraintes du canal dans le domaine fréquentiel et la matrice M2t exprimant les contraintes du canal dans le domaine temporel. Ainsi, dans une étape 24a/24b, les matrices M2f et M2t de covariance du canal 8, respectivement dans le domaine fréquentiel et dans le domaine temporel, sont déterminées par les moyens de détermination 16 en utilisant respectivement la fenêtre d'étalement du profil temporel et la fenêtre d'étalement du profil fréquentiel du canal 8 stockées dans la mémoire 10. Ces matrices M2f et M2t sont stockées dans la mémoire 10. Selon l'invention, la matrice de covariance intermédiaire M2 est déterminée par les moyens de détermination 16 dans une étape 26 selon le produit de Kronecker : M2 M2r °M2f. " Une matrice E' est ensuite calculée par les moyens de calcul 18 dans une étape 28 selon la relation E' = E.TH , où T=T,OTf et dans laquelle E est une matrice diagonale comprenant les symboles pilotes émis dans une trame aux positions des symboles pilotes, les autres symboles étant nuls. Le vecteur colonne C' qui minimise la relation + CH. M11. C', est 2 ensuite calculé par les moyens d'estimation 20 dans une étape 30a/30b en mettant en oeuvre le critère du maximum a posteriori ou MAP afin d'estimer le canal de propagation 8 selon la relation c = TH.c'. L'étape 30a/30b de calcul du vecteur C' est décrite ci-après dans un premier mode de mise en oeuvre 30a (étapes 32 à 38) et dans un deuxième mode de mise en oeuvre 30b (étapes 32 à 34 et 40 à 44) du procédé selon l'invention. a) Premier mode de mise en oeuvre L'organigramme de la figure 2 illustre un premier mode de mise en oeuvre du procédé d'estimation de canal de l'invention mettant en oeuvre le critère du MAP. Dans une étape 32, les matrices m2f et M2, de covariance sont décomposées en vecteurs propres et en valeurs propres par les moyens de décomposition 18 25 selon les relations MF = WFHNFWF et MT = WTHNTWT , dans lesquelles : WF est une matrice de vecteurs propres de la matrice m2f ; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice m2f ; - W. est une matrice de vecteurs propres de la matrice M2, ; et NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice M2t . 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 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 intermédiaire globale M2 du canal 8 est alors égale à M2 = WHNW. L'expression [1] peut donc s'écrire :(E'H Et+o-2.WH.N-1.W).C=EIH.R En posant C'= w".-/.b où b est un vecteur de dimension n, et oùIN 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 : (72.4 = [2] où I est la matrice identité. A l'étape 36, les moyens d'estimation de canal 20 calculent la matrice P .E' .W H .-J7^1 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/T.XH .(Q + cr2J) [3] . a) Deuxième mode de mise en oeuvre L'organigramme de la figure 3 illustre un deuxième mode de mise en oeuvre du procédé d'estimation de canal de l'invention mettant également en oeuvre le critère MAP. Sur cette figure, les étapes 24, 26, 28 et 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 :(E'H !-1 MP).0 = E'll .R En posant C'= où b est un vecteur de dimension n', l'expression [1] s'écrit (-N = .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 P'=-^11'.W.E'H.E'.W'H. 1.N1' et décomposent cette matrice P' en vecteurs propres et valeurs propres selon la relation P=XIHQIX 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 44, les moyens d'estimation de canal 20 estiment le canal 8 à partir de la relation [4] selon l'expression [5] suivante C'= W'H .(Q'+cs2.1)-1.X1.1'.WEIH .R [5]. La matrice W' é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 10 matrice NF correspondant aux n't vecteurs propres de la matrice w, 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 deux matrices et elle correspond alors à une dimension totale de Rn'f, nf) + (n't, nt)], ce qui permet de 15 réduire la complexité du procédé. Dans les deux modes de mises en oeuvre décrits ci-avant, la matrice du canal réel C est ensuite estimée selon l'invention dans une étape 50 selon la relation C = TH. C' à partir du vecteur C' calculé. 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. 20 25 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 utilisant chacune Nf sous-porteuses fréquentielles sur chacune desquelles sont émis Nt symboles, où parmi l'ensemble des symboles, certains symboles, dits symboles pilotes, sont connus dudit récepteur (2), ledit signal étant synchronisé par le récepteur (2) à partir d'une position de synchronisation temporelle et d'une position de synchronisation fréquentielle, le profil temporel du canal (8) étant symétrique et centré autour d'une position centrale effective temporelle tmoyen distincte de la position de synchronisation temporelle, le profil fréquentiel du canal (8) étant symétrique et centré autour d'une position centrale effective fréquentielle f -moyen distincte de la position de synchronisation fréquentielle, le procédé, mis en oeuvre par le récepteur (2), étant caractérisé par le fait qu'il comprend les étapes de : - détermination (24a) d'une matrice de covariance fréquentielle intermédiaire M2f représentant les variations fréquentielles du canal (8) dues à un profil temporel théorique du canal (8), symétrique et centré sur la position de synchronisation temporelle, ce profil temporel théorique du canal étant obtenu à partir du profil temporel du canal translaté en temps de (-Tmoyen) de façon à être symétrique autour de la position de synchronisation temporelle ; - détermination (24b) d'une matrice de covariance temporelle intermédiaire M2t représentant les variations temporelles du canal (8) dues à un profil fréquentiel théorique du canal (8), symétrique et centré sur la position de synchronisation fréquentielle, ce profil fréquentiel théorique du canal étant obtenu à partir du profil '--moyen) fréquentiel du canal translaté en fréquence de ( f de façon à être symétrique autour de la position de synchronisation fréquentielle; - calcul (26) d'une matrice de covariance globale intermédiaire M2 selon le produit de Kronecker: M2 = M2t 0 M2 f ; - calcul (28) d'une matrice E' selon la relation E', E.T H , où T = Tt 0 T f et dans laquelle E est une matrice diagonale comprenant les symboles pilotes émis dans une trame aux positions des symboles pilotes, les autres symboles étant nuls, Tfest une matrice diagonale (Ni, Nt) de centrage définie par 0 - Tj 0 où v=2A-.2-mown et T, est une matrice diagonale (Nt, Nf 1 0 0 0 e"e ei2Y Nt) de centrage définie par T, = où = -27r. fmoven ; Nt - calcul (30a ; 30b) d'un vecteur colonne représentant un canal intermédiaire C' qui IIR-E1 .C112 fli -1 minimise la relation + C .M2 . C ou R est un vecteur colonne de 0.2 dimension (Nt x Nt) constitué des symboles reçus, et ; - estimation (50) du canal de propagation selon la relation C =T .C' , à partir du vecteur C' calculé. 2. Procédé selon la revendication 1, caractérisé en ce que l'étape de calcul du vecteur canal intermédiaire C' comprend les sous-étapes de : - décomposition (32) des matrices M2f et 102t de covariance en vecteurs propres selon les relations M2f = WFHNFWF et M2, = WTHNTWT , dans lesquelles : - WF est une matrice de vecteurs propres de la matrice M2f ; - NF est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice M2f ; WT est une matrice de vecteurs propres de la matrice M2t ; - NT est une matrice diagonale comprenant des valeurs propres associées aux vecteurs propres de la matrice M2t ; - 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- calcul du vecteur canal intermédiaire C' à 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. 3. Procédé selon la revendication 2, caractérisé en ce que l'étape de calcul du vecteur canal intermédiaire C' comprend les sous-étapes de : - calcul (36) d'une matrice P selon la relation P = 11.W.EH.E.W11.1 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 - calcul du vecteur canal intermédiaire C' selon la relation C= WH.-\ 11.X11.(Q+0-2.1)1.X.I.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 cy2 représente la puissance du bruit dans le canal. 4. Procédé selon la revendication 2, caractérisé en ce que l'étape de calcul de du vecteur canal intermédiaire C' comprend les sous étapes : - d'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 - calcul du vecteur canal intermédiaire C' à 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. 5. Procédé selon la revendication 4, caractérisé en ce que les n' valeurs propres sont les plus grandes valeurs propres de la matrice N. 6. Procédé selon la revendication 4, dans lequel les n' vecteurs propres de la matrice W' sont les produits de Kronecker d'un nombre limité n't 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 5 valeurs propres de la matriceNF correspondant aux n't vecteurs propres de la matrice WF et des valeurs propres de la matrice NT correspondant au n't vecteurs propres de la matrice WT 7. Procédé selon l'une des revendications 4 à 6, dans lequel n' est inférieur ou 10 égal au nombre de symboles pilotes dans chaque trame. 8. Procédé selon l'une des revendications 4 à 7, dans lequel l'étape de calcul du vecteur canal intermédiaire C' comprend les sous-étapes de : - calcul (42) d'une matrice P' selon la relation P'=-N.W'.EH.E.W'H..N dans 15 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 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 20 aux vecteurs propres de la matrice P'; et - calcul du vecteur canal intermédiaire C' selon la relation C =+ cr2.1)1.X'.-^ 1/V1'.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 0-2représente la puissance du bruit dans le canal. 25 9. Procédé selon l'une des revendications précédentes, dans lequel le signal est un signal multiporteuses, notamment un signal OFDM. 10. Equipement récepteur (2) apte à recevoir d'un émetteur (6) au moins un signal 30 transmis à travers un canal de propagation radioélectrique (8), ledit signal comprenant des trames utilisant chacune Nt sous-porteuses fréquentielles sur chacune desquelles sont émis Nt symboles, où parmi l'ensemble des symboles,certains symboles, dits symboles pilotes, sont connus dudit équipement récepteur (2), le signal étant synchronisé par le récepteur (2) à partir d'une position de synchronisation temporelle et d'une position de synchronisation fréquentielle, le profil temporel du canal (8) étant symétrique et centré autour d'une position centrale effective temporelle tmoyen distincte de la position de synchronisation temporelle, le profil fréquentiel du canal (8) étant symétrique et centré autour d'une position centrale effective fréquentielle f -moyen distincte de la position de synchronisation fréquentielle, le procédé, mis en oeuvre dans le récepteur, est caractérisé par le fait qu'il comprend: - des moyens (16) de détermination d'une matrice de covariance fréquentielle intermédiaire M2f représentant les variations fréquentielles du canal (8) dues à un profil temporel théorique du canal (8), symétrique et centré sur la position de synchronisation temporelle, ce profil temporel théorique du canal étant obtenu à partir du profil temporel du canal translaté en temps de ( Ctmoyen) de façon à être symétrique autour de la position de synchronisation temporelle ; - des moyens (16) de détermination d'une matrice de covariance temporelle intermédiaire M2t représentant les variations temporelles du canal (8) dues à un profil fréquentiel théorique du canal (8), symétrique et centré sur la position de synchronisation fréquentielle, ce profil fréquentiel théorique du canal étant obtenu à partir du profil fréquentiel du canal translaté en fréquence de C-( f moyen) de façon à être symétrique autour de la position de synchronisation fréquentielle ; - des moyens (16) de calcul d'une matrice de covariance globale intermédiaire M2 selon le produit de Kronecker: M2 = M2t ° M2 f ; - des moyens (18) de calcul d'une matrice E' selon la relation E=E.TH , où T = T, 0 T f et dans laquelle E est une matrice diagonale comprenant les symboles pilotes émis dans une trame aux positions des symboles pilotes, les autres symboles étant nuls, Tf est une matrice diagonale (Nf, Nf) de centrage définie par0 0 32 0 où q)=27r. rm `n" et T, est une matrice diagonale (Nt, Nf 0 ei2e' Tf = 1 0 0 0 e'w e2" 0 O ei(M-1)y, Nt) de centrage définie par T, = l'Inn 'en OÙ = -271«. ; N, - des moyens (20) de calcul d'un vecteur colonne représentant un canal 1 2 intermédiaire C' qui minimise la relation II R-E 11 ± CH . Mil . C où R est un vecteur c.2 colonne de dimension (Nt x Nf) constitué des symboles reçus, et ; - des moyens (20) d'estimation du canal de propagation (8) selon la relation C = TH .C', à partir du vecteur C' calculé. 11. Programme d'ordinateur comportant des instructions pour la mise en oeuvre du 10 procédé selon l'une des revendications 1 à 9 lorsque le programme est exécuté par au moins un processeur.