Brevet FR2798542A1 — publié le 20010316

RECEPTEUR A MULTIPLEXAGE PAR REPARTITION EN FREQUENCES ORTHOGONALES AVEC ESTIMATION ITERATIVE DE CANAL ET PROCEDE CORRESPONDANT DESCRIPTION Domaine technique La présente invention a pour objet un récepteur à multiplexage par répartition en fréquences orthogonales avec estimation itérative de canal et un procédé correspondant. Elle trouve une application en radiocommunications et plus particulièrement dans la technique d'accès multiple dite OFDM pour "Orthogonal Frequency Division Multiplexing". Elle peut être appliquée-entre autres- au système radiomobile européen HIPERLAN II. Etat de la technique antérieure La technique OFDM [1] est une technique multiporteuse qui permet de répartir les utilisateurs dans le plan temps-fréquence d'une façon simple. Elle permet en outre de transmettre des signaux à haut débit sans avoir à utiliser d'égaliseur. Cette méthode a été largement utilisée dans des contextes de diffusion tels que le DVB-T ("Digital Video Broadcasting") [2] et le DAB ("Digital Audio Broadcasting") [3]. En contexte radiomobile, l'OFDM est présent dans la norme HIPERLAN II. La technique OFDM est à la fois une technique d'accès multiple et une technique de modulation. Le principe de base de la technique OFDM est de produire un certain nombre de signaux à bande étroite tous orthogonaux entre eux. En prenant les précautions d'usage, on utilise ces propriétés d'orthogonalité pour retrouver les données transmises. La réalisation d'un tel système fait appel à la mise en oeuvre d'une transformée de Fourier inverse à l'émission et d'une transformée de Fourier à la réception. La figure 1 annexée illustre une chaîne de transmission OFDM classique avec un seul capteur. Cette chaîne comprend un circuit 10 de conversion sérieparallèle recevant des symboles A, un circuit 12 de transformation de Fourier inverse 12, des moyens d'émission 14, des moyens de réception 20, un circuit 22 de transformation de Fourier, un convertisseur parallèle-série 24 et enfin un moyen de décision 26 qui restitue les symboles estimés Â. L'émetteur OFDM classique traite ,le flux de données par bloc. Il gère ce flux par séquences de Nt symboles et en réalise la transformée de Fourier inverse. Cela revient à dire que la transformée de Fourier inverse produit Nf sous-porteuses, chacune véhiculant un des symboles de la séquence de départ. Ce bloc, appelé symbole OFDM, contient les symboles de données, mais peut aussi contenir des symboles pilotes qui peuvent être utilisés pour des besoins de synchronisation ou d'estimation de canal. Contrairement au cas des signaux CDMA ("Code Division Multiple Access") ou TDMA ("Time Division Multiple Access"), où un symbole pilote occupe d'emblée l'ensemble de la bande de transmission, la technique OFDM nécessite de réellement répartir des symboles pilotes sur l'ensemble du plan temps-fréquence. Le canal radiomobile emprunté lors d'une communication entre un émetteur et un récepteur est généralement de type multi-trajets avec évanouissements rapides de Rayleigh. Ce phénomène est dû à la conjonction du mouvement du mobile et de la propagation de l'onde radioélectrique selon plusieurs chemins. Le récepteur traite le signal reçu par bloc de symboles OFDM (bloc temps-fréquence). Le signal est reçu sur un réseau de L capteurs, créant L branches de diversité. L'estimation de canal est réalisée sur chacune de ces branches et les résultats sont combinés au sens du MRC ("Maximum Ratio Combining") pour finalement estimer les données émises. Un récepteur à L branches de diversité est représenté sur la figure 2. Il comprend L,capteurs 301, 302, ...,30L, L circuits de transformation de Fourier 321, 322, 32L, L convertisseurs parallèle-série 341, 342, 34L, L circuits d'estimation de canal 361, 362,36L et un additionneur 38 délivrant les symboles estimés Â. Du point de vue du récepteur, après démodulation, le canal affectant un bloc temps-fréquence peut être représenté sous forme d'une matrice temps-fréquence, ou encore d'une surface dans l'espace temps-fréquenceamplitude. Le problème est donc traité dans un espace bidimensionnel, contrairement au TDMA [4] où le problème est unidimensionnel. L'estimation de canal est basée sur l'utilisation des symboles pilotes. Ils permettent d'obtenir directement une estimation du canal aux emplacements des pilotes en vue d'une interpolation pour estimer le canal affectant le reste des symboles. Ces techniques présentent des inconvénients. En effet, le canal vu par le récepteur peut varier de manière significative d'un bloc temps-fréquence à autre. Cette variation est principalement due aux changements de conditions de propagation entre l'émetteur et le récepteur. D'un point de vue physique, le caractère variable du canal peut être caractérisé par le produit BdxTm où Bd représente la largeur de la bande Doppler et Tm l'étalement des retards. Plus le produit BdxTm est grand, plus le canal varie rapidement dans les domaines temporel et fréquentiel. Les méthodes de réception de l'art antérieur ne cherchent pas à optimiser l'estimation du, canal. Elles se contentent de réaliser une estimation du canal aux positions des symboles pilotes puis d'étendre cette estimation aux données par interpolation. Les interpolations sont généralement réalisées de manière linéaire. Trois des méthodes les plus courantes peuvent être évoquées : *La première prend en compte les trois plus proches symboles pilotes du symbole auquel on veut estimer le canal. On calcule le plan passant par les trois symboles pilotes et on en déduit le canal au point considéré. Même en respectant le critère de Nyquist du point de vue des symboles pilotes, c'est-à-dire en utilisant suffisamment de symboles pilotes et en les répartissant de manière à échantillonner correctement le plan temps- fréquence, cette méthode est sensible aux fortes variations de canal et ne permet pas de réaliser une estimation fiable du canal, surtout dans les cas où le produit BdxTm est élevé. *La deuxième méthode est une forme simple de la technique du MMSE ("Minimum Mean Square Error") : elle consiste à rechercher le plan constant moyennant les valeurs du canal au niveau des symboles pilotes et à en déduire les valeurs du canal affectant les données émises. Cette modélisation de canal est bien adaptée à des canaux variant très peu sur le bloc reçu, c'est-à- dire pour des produits BdxTm faibles. Cependant, dès que le canal devient plus sélectif, la modélisation plane montre ses limites et les performances sont dégradées. *La troisième méthode est une autre forme de MMSE, avec un plan recherché non constant. Cette méthode s'adapte donc mieux aux cas où le canal varie lentement, mais est moins adaptée que la deuxième dans le cas de canaux presque constants. Ces trois méthodes sont donc adaptées à des cas très spécifiques de propagation, mais en aucune façon à des canaux de type multitrajet sélectifs en temps et en fréquence. La présente invention a justement pour but de remédier à cet inconvénient. Exposé de l'invention La présente invention a pour but principal d'améliorer les performances des systèmes OFDM existants ou à venir. Cette amélioration, obtenue par une optimisation de l'estimation de canal, permet d'augmenter sensiblement la capacité du système. Cette amélioration est engendrée par une optimisation du fonctionnement du récepteur OFDM dans le cas des évanouissements lents mais également dans le cas plus complexe des évanouissements très rapides. Il est alors possible de contrecarrer les dégradations en performance engendrées par une variation rapide de canal sur le bloc temps-fréquence considéré en réception. L'invention permet de réduire, à qualité de réception constante, le nombre relatif et/ou la puissance des symboles pilotes. Ce but est atteint par la prise en compte de manière optimale d'un nombre arbitraire de symboles pilotes de blocs temps-fréquence consécutifs dans l'estimation de canal. Il est également atteint par le caractère optimal de la prise en compte dans l'estimation de canal (d'une partie ou de la totalité) des symboles de données de ces blocs, qui sont bien sûr plus nombreux que les symboles pilotes. L'invention peut être utilisée quelle que soit la façon dont les symboles pilotes sont introduits dans le flux d'information transmise. Le récepteur de l'invention effectue un traitement bloc par bloc à chaque fois qu'un nombre donné de symboles OFDM est disponible. Sur chaque branche de diversité, le canal multi-trajets est tout d'abord estimé grossièrement au moyen des symboles pilotes associés au bloc reçu et éventuellement à d'autres blocs. Cette estimation a pour but l'initialisation de l'algorithme d'estimation de canal itératif. On traite ensuite l'ensemble des symboles (données et pilotes) pour obtenir l'estimation de canal permettant de générer des sorties souples des symboles de données émis. Les sorties souples obtenues à la fin d'une itération peuvent à nouveau être utilisées conjointement aux symboles pilotes pour apporter une amélioration supplémentaire à l'estimation de canal, et donc améliorer davantage les estimations souples des symboles de données. La technique proposée permet en outre de tenir compte de la structure codée des symboles de donnée et peut être optimisée dans ce sens, ce qui conduit à produire des sorties souples de meilleure qualité. L'estimation du canal multi-trajets repose, d'une part, sur l'utilisation de l'algorithme itératif dit E. M. ("Expectation Maximisation") [5], [6], [7], pour trouver la réalisation de canal la plus vraisemblable conditionnellement au bloc reçu à traiter et au codage du canal éventuellement utilisé. Elle repose également sur la décomposition du canal multi-trajets bidimensionnel sur chaque branche de diversité selon le théorème d'expansion de Karhunen-Loève [8]. Cette décomposition permet une caractérisation souple des variations temporelles des trajets dues à l'effet Doppler ainsi que des variations fréquentielles dues à l'étalement temporel et s'intègre facilement dans l'algorithme E. M. lui-même. De façon plus précise, la présente invention a pour objet un récepteur pour radiocommunications à multiplexage par répartition en fréquences orthogonales (OFDM) comprenant : i) une pluralité de L branches de diversité traitant des blocs de symboles numériques, chaque bloc comprenant des symboles de données et des symboles pilotes répartis dans un bloc bidimensionnel temps-fréquence à Nt intervalles de temps et à Nf intervalles de fréquence, chaque branche de diversité comprenant un capteur radioélectrique, des moyens délivrant un signal de sortie à N composantes constituant les composantes d'un vecteur R#, où désigne le rang de la branche de diversité ( allant de 0 à L-l), ii) un estimateur de canal traitant les L signaux délivrés par les L branches de diversité et délivrant des estimations souples des symboles de données, iii) un organe de décision recevant les estimations souples des symboles de données et délivrant une estimation des symboles de données, ce récepteur étant caractérisé en ce que : (a) l'estimateur de canal traite un vecteur Ce à N composantes caractérisant le canal dans un bloc bidimensionnel temps-fréquence, ces moyens d'estimation étant aptes à définir une base de N vecteurs Bk qui sont les N vecteurs propres normalisés de la matrice de covariance temps-fréquence du canal, ces moyens décomposant chaque vecteur Ce dans cette base, les N coefficients de cette décomposition étant notés Gk# avec k allant de 0 à N-1, les coefficients Gk# définissant, pour chaque branche de diversité, un vecteur G#, qui est une représentation du canal vu en sortie de ladite branche de diversité, (b) l'estimateur de canal traite un nombre fini (D) d'itérations selon un algorithme d'estimation-maximisation (EM) fondé sur le critère de probabilité maximum a posteriori (MAP), les moyens d'estimation étant initialement mis en oeuvre en prenant en compte les symboles pilotes contenus dans le bloc bidimensionnel temps-fréquence considéré et éventuellement des symboles pilotes contenus dans les blocs temps-fréquence voisins, ce qui conduit à une estimation d'ordre 0, l'estimateur de canal prenant en compte, ensuite, les symboles de données pour les autres itérations et ainsi de suite, l'estimateur de canal délivrant finalement, après une dernière itération D, les coefficients optimaux Gk#(D) (k de 0 à (N-1) et de 0 à L-1), définissant chaque branche de diversité, le vecteur Ge représentant le canal. Dans un mode particulier de réalisation Nt=Nf et les blocs temps-fréquences sont carrés. La présente invention a également pour objet un procédé de réception dont les opérations correspondent aux fonctions remplies par les différents moyens du récepteur qui vient d'être défini. Brève description des figures - la figure 1, déjà décrite, montre une chaîne de transmission OFDM classique à un seul capteur ; - la figure 2, déjà décrite, montre un récepteur OFDM classique à plusieurs capteurs et L branches de diversité ; - la figure 3 illustre un récepteur conforme à l'invention ; - la figure 4 illustre le processus d'estimation itérative selon l'invention ; - la figure 5 illustre un mode particulier de réalisation du récepteur de l'invention ; - la figure 6 montre un exemple de répartition des symboles pilotes et des symboles de données dans un bloc temps-fréquence ; - la figure 7 est une représentation des vecteurs propres de la matrice de corrélation du canal à spectre Doppler classique et profil d'intensité multi-trajets exponentiel pour BdTm=10-5 ; - la figure 8 est une représentation des mêmes vecteurs propres mais pour BdTm=10-3 ; - la figure 9 donne les variations du taux d'erreur binaire (TEB) en fonction du rapport Es/Io pour un récepteur conforme à l'invention et pour divers récepteurs classiques, dans le cas d'un produit BdTm égal à 10-5 avec 16 symboles pilotes répartis comme indiqué sur la figure 5 ; - la figure 10 donne ces mêmes variations mais pour un produit BdTm égal à 10-3. Description détaillée de modes particuliers de réalisation Un récepteur conforme à l'invention est représenté schématiquement sur la figure 3. Ce récepteur comprend des moyens déjà représentés sur la figure 2 et qui portent les mêmes références. Mais le récepteur comprend des moyens pour mettre en oeuvre un algorithme d'estimation itératif qui est schématisé par le bloc (E) qui alimente l'organe de décision 42. L'entité élémentaire du signal OFDM correspond à la transformée de Fourier inverse d'une séquence de symboles. Le récepteur de l'invention traite bloc par bloc le signal reçu. La taille du bloc traité ne dépend pas nécessairement du nombre de porteuses du système OFDM et peut prendre en compte tout ou partie d'un ou de plusieurs symboles OFDM. La forme et la taille du bloc traité en réception est libre, de manière à s'adapter au mieux au système. L'estimation de canal est réalisée bloc par bloc. Un bloc est composé de N symboles amn d'énergie Emn et de position bidimensionelle (mF, nT) où F et T sont l'espacement respectivement en fréquence et en temps entre deux symboles adjacents. Ces symboles prennent leurs valeurs dans un alphabet Q de type modulation à déplacement de phase (MDP) arbitraire. Chaque bloc est composé de ND symboles de donnée indexés dans l'ensemble SD et Np symboles pilotes indexés dans l'ensemble Sp. En règle générale, les récepteurs classiques utilisent des symboles pilotes de plus forte puissance que les symboles de données. Cette différence de puissance permet au mieux d'estimer le canal, mais risque d'introduire des interférences entre porteuses et donc de réduire la capacité du système OFDM. La technique d'estimation de canal selon l'invention permet de réaliser une estimation de canal optimale quelle que soit la valeur de la puissance des symboles pilotes ; dans la suite, on notera Ep l'énergie des symboles pilotes et ED l'énergie des symboles de donnée. Le canal multitrajets emprunté par le signal OFDM est composé de plusieurs trajets présentant ou pouvant présenter des variations temporelles et fréquentielles dues à l'effet Doppler. Chaque trajet est caractérisé par une puissance moyenne et un spectre de puissance Doppler (SPD) donnés qui dépendent de l'environnement et de la vitesse du mobile. De plus, les évanouissements subis par chaque trajet peuvent être aussi bien de type Rayleigh que de type Rice. On note Jo(.) la fonction de Bessel de première espèce d'ordre 0. A titre d'exemples, la fonction d'autocorrélation temps-fréquence du canal à spectre de puissance Doppler classique et au profil d'intensité multitrajets exponentiel de puissance moyenne (0,0) vu sur une branche de diversité est donnée par Le récepteur multicapteurs OFDM de l'invention est composé de L capteurs 301, 302, ..., 30L décorrélés spatialement, donnant naissance à L branches de diversité. Sur chacune de ces branches, le signal reçu est en premier lieu démodulé par la transformée de Fourier dans les circuits 321, 322, ..., 32L. On suppose que le signal en sortie de la 1ième branche de diversité associé au symboles amn s'écrit Rmn# = Cmn#amn# + Nmn# où c est le facteur de gain du canal discret de la #ème branche de diversité vu par le symbole amn et N est un bruit blanc ganssien complexe additif de variance No. Les facteurs de gain sont indépendants d'une branche de diversité à l' autre, mais corrélés en temps et en fréquence entre eux sur la même branche. L'invention a pour but d'estimer les facteurs de gain c du canal. Soit ( . )l'opérateur de transposition. Pour des raisons de notation, on introduit la fonction d'indexation 8(k)=(m(k),n(k))entre l'ensemble monodimensionnel {k}k=0N-1 et l' ensemble d'indexation bidimensionnel SD#SP. De plus, pour chaque bloc transmis, on introduit le vecteur signal en sortie du filtre adapté de la #ème branche de diversité Pour s'affranchir de la dépendance de l'amplitude de chaque symbole MDP amn sur son index (m, n), on introduit le vecteur normalisé du bloc transmis : avec A#(k) = a#(k)/#a#(k)#. . Sur ces bases, il est possible de réécrire les composantes du vecteur reçu sur la #ème branche de diversité : où C##(k) est la # (k) ème composante du vecteur canal discret multiplicatif équivalent sur la # ème branche : On cherche à estimer pour chaque bloc et sur chaque branche le vecteur C . Pour l'estimation du canal, selon le critère du Maximum a Posteriori, on utilise une représentation adéquente du canal multitrajets discret pour chaque branche de diversité. Cette représentation est basée sur le théorème de décomposition orthogonale de Karhunen-Loève. Le vecteur représentant le canal multitrajets discret équivalent sur la #ème branche de diversité CI! s'exprime de la façon suivante: où {Bk}k=0N-1 sont les N vecteurs propres normalisés de la matrice d'autocorrélation temps-fréquence du canal discret F=E[C C#*T] et {Gk#}k=0N-1 sont des . variables aléatoires gaussiennes complexes, indépendantes, centrées et de variances égales aux valeurs propres {#k}k=0# de la matrice hermitienne F. La (p,q)ème entrée de la matrice F est donnée par Les vecteurs {G#}l=0L-1, avec G# = (G0#, ... , G0N-1#)Tsont la représentation adéquate du canal discret multiplicatif vu en sortie des L branches de diversité. L'estimation du canal revient donc à estimer {Gk#}k=0,#=0N-1,L-1. Cette estimation est réalisée itérativement par la formule suivante en notant Gp#(d) l'estimation de Gp1 à la dème itération de l'algorithme : où Bp# (k) est la kème composante de Bp et Le codage n'est pas pris en compte dans cette formule. A cet égard, si une partie du vecteur A transmis est codée par un code quelconque (convolutif, en bloc,...) alors, à chaque itération de l'algorithme EM, les probabilités conditionelles discrètes P(A5(k) = A) RÊ1L-1 @{G (d) L-1 peuvent être calculées exactement en utilisant le treillis de ce code et l' algorithme de Bahl [9]. L'initialisation de l'algorithme est réalisée par la projection des symboles pilote reçus sur les N vecteurs propres {Bk}k=0N-1 de la matrice de corrélation F. On utilise donc: comme pème composante de l' estimation initiale G(0). Les vecteurs {Bk}k=0N-1 sont des éléments connus au niveau du récepteur. Ils correspondent à un modèle de canal choisi. On obtient cette famille de vecteurs en calculant la matrice F de corrélation théorique du modèle correspondant et en tirant de celle-ci ses vecteurs propres {Bk}k=0N-1 et ses valeurs propres {#k}k=0N-1 associées. La recombinaison des L branches permettant de tirer profit de la diversité spatiale du récepteur est réalisée lors du traitement itératif des probabilités P(A5(k) = A 1 {R f } :6, {G f Cd)} :6). Après un nombre D choisi d'itérations de l'algorithme, les probabilités P(A8(k) = A 1 {R f } :6, @{G (d L-1 permettent d'obtenir une estimation souple des symboles de donnée. La décision est réalisée en traitant ces sorties souples soit directement par un organe de décision si les données ne sont pas codées, soit par un décodeur si les données sont codées. L'algorithme d'estimation améliorée de canal selon l' invention est illustré sur la figure 4. L'algorithme E. M. est symbolisé par le bloc 50, la sélection de la base adéquate par le bloc 52 qui reçoit le produit BdTm et le calcul des pondérations par le bloc 54 qui reçoit BdTm et No. L'algorithme E. M. du bloc 50 est itératif et comprend une initialisation représentée schématiquement par le bloc 3o et D itérations 3d. La démodulation est représentée par le bloc 58 qui délivre les probabilités P(A # # # {R#}#=0L-1, {G#(d)}#=0L-1). Un récepteur conforme à la description générale de l'invention est représenté schématiquement sur la figure 5. Ce récepteur comprend des moyens déjà représentés sur la figure 2 et qui portent les mêmes références. Mais le récepteur comprend des moyens pour mettre en oeuvre un algorithme d'estimation itératif qui est schématisé par le rebouclage de la sortie de l'additionneur 38 sur les moyens d'estimation 361, 36L. Par ce rebouclage, on applique à ces moyens un signal noté ##(k)(d) où d représente le rang de l'itération. Le récepteur représenté comprend en outre un interrupteur symbolique 40 qui se ferme à la dernière itération D pour alimenter l'organe de décision 42. Il est supposé, dans ce cas particulier de réalisation, que le signal OFDM reçu est composé par des symboles prenant leurs valeurs dans un alphabet Q de type modulation à déplacement de phase à 2 ou 4 états. Les autres caractéristiques du signal reçu telles la forme des blocs, leur taille, la répartition et la puissance des symboles pilotes, ainsi que les caractéristiques du canal sont les mêmes que dans le cas général précédemment décrit. La supposition effectuée dans ce cas particulier permet d'obtenir une expression analytique de l'estimation de Gp#(d+1). Dans le cas où les symboles proviennent d'une constellation de type MDP-2, l'espression de Gp1(d+1) est la suivante: Dans le cas où les symboles proviennent d'une constellation de type MDP-4, l'expression de Gp#(d+1) est la suivante: où D#(k) est la valeur prise par le symbole pilote A#(k), # (k) # Sp et : Dans cette expression, la parenthèse représente le complexe conjugué ci(: de la # (k) ème composante du canal. La recombinaison des L branches permettant de tirer profit de la diversité en réception, du récepteur est réalisée lors du traitement itératif de la variable ##(k)(d). La décision est soit directement réalisée sur cette variable, si les données ne sont pas codées, soit par un décodeur si les données sont codées, après un nombre D d'itérations choisi de l'algorithme. De plus, dans le cas d'une modulation MDP-2, on peut utiliser la variante suivante: Le récepteur et le procédé de l' invention ont été simulés pour un canal multi-trajets avec un produit BdxTm égal à 10-5. Le bloc traité par le récepteur est un carré contenant 256 symboles dont 16 symboles pilotes. Les symboles pilotes sont équirépartis dans chaque bloc temps-fréquence comme illustré sur la figure 5 où ils figurent sous forme de carrés noirs, les carrés blancs représentant les symboles de données. Dans la simulation, le récepteur suppose que le canal est de spectre de puissance constant et de profil d'intensité multi-trajets constant. Sa fonction d'autocorrélation temps-fréquence s'écrit donc : Il est possible d'utiliser les vecteurs propres provenant d'un modèle de canal arbitraire. Par exemple, un modèle de canal à spectre de puissance Doppler et de profil d'intensité multi-trajets exponentiel est illustré par les figures 6 et 7. Ces figures représentent les vecteurs propres de la matrice de corrélation du canal à spectre Doppler classique et profil d'intensité multi-trajets exponentiel pour respectivement BdTm=10-5 (figure 6) et BdTm=10-3 (figure 7) . Pour illustrer les avantages de l'invention, les figures 8 et 9 comparent les performances d'un récepteur selon l'invention avec celles d'un récepteur classique de type MMSE. Ces figures donnent le taux d'erreur binaire TEB en fonction du rapport signal à bruit Es/I0. Dans les deux cas, les symboles pilotes sont supposés équirépartis comme illustré sur la figure 5. Pour la figure 8, le produit BdTm vaut 10-5 et il vaut 10-3 pour la figure 9. Les quatre courbes de ces figures correspondent respectivement à : 61, 71 : estimateur selon l'invention 62, 72 : estimateur à MLSE constant 63, 73 : estimateur à MLSE 64, 74 : courbe de limite théorique On voit que, dans les deux cas, l'invention conduit à de meilleures performances qu'avec les techniques classiques. A titre d'exemple, pour un TEB brut de 10-2, l'utilisation de l'invention garantit un gain en terme de rapport signal à bruit de 3 dB par rapport au plus favorable des récepteurs classiques dans le cas d'un canal dispersif avec BdXTm=10-3. Références [1] J. A.C. Bingham "Multicarrier modulation for data transmissions : an idea whose time has corne", IEEE Communication Magazine, 28 (5) May 1990. [2] "Digital Broadcasting systems for television, sound and data services", European Telecommunication Standard, prETS 300 744 (Draft, version 0. 0.3), April 1996. [3] "Radio Broadcasting Systems ; Digital Audio Broadcasting (DAB) to mobile, portable and fixed receivers", ETS 300 401, ETSI, European Télécommunication Institute, Valbonne, France, Feb. 1995. [4] EP-0 802 656, "Signal numérique à blocs de référence multiples pour l'estimation de canal, procédés d'estimation de canal et récepteurs correspondants". [5] A.P. Dempster, N. M. Laird and D. B. Rubin "Maximum Likelihood from incomplete data via the EM algorithm", Journal of the Royal Statistical Society, 39,1977. [6] G.K. Kaleh "Joint carrier phase estimation and symbole decoding of trellis codes", European Transactions on Télécommunications and Related Technologies, San Diego, Ca, Jan. 1990. [7] C. N. Georghiades and J.C. Han : "Séquence estimation in the présence of random parameters via the EM algorithm", IEEE Transactions on Communications, 45, n 3, mar. 1997. [8] J.G. Proakis : Digital Communications, McGraw-Hill, New York 1989. [9] L.R. Bahl, J. Cocke, F. Jelinek and J. Raviv : "Optimal decoding of linear codes for minimizing symbol error rate", IEEE Transactions on Information Theory, 20, Mar. 1974. REVENDICATIONS 1. Récepteur pour radiocommunications à multiplexage par répartition en fréquences orthogonales (OFDM) comprenant : i) une pluralité de L branches de diversité traitant des blocs de symboles numériques, chaque bloc comprenant des symboles de données et des symboles pilotes répartis dans un bloc bidimensionnel temps-fréquence à Nt intervalles de temps et à Nf intervalles de fréquence, chaque branche de diversité comprenant un capteur radioélectrique (301, 302, ... , 30L) , des moyens (321, 322, ... , 32L), (341, 342, 34L) délivrant un signal de sortie à N composantes constituant les composantes d'un vecteur R , où # désigne le rang de la branche de' diversité ( # allant de 0 à L-1), ii) un estimateur de canal (E) traitant les L signaux délivrés par les différentes branches de diversité et délivrant des estimations souples des symboles de données, iii) un organe de décision (42) recevant les estimations souples de symboles de données et délivrant une estimation (Â) des symboles de données, ce récepteur étant caractérisé en ce que : (a) l'estimateur de canal (E) traite un vecteur Ce à N composantes caractérisant le canal dans un bloc bidimensionnel temps-fréquence, ces moyens d'estimation étant aptes à définir une base de N vecteurs Bk qui sont les M vecteurs propres normalisés de la matrice de covariance temps- fréquence du canal, ces moyens décomposant chaque vecteur Cedans cette base, les N coefficients de cette décomposition étant notés Gk# avec k allant de 0 à N-l, les coefficients Gk# définissant, pour chaque branche de diversité, un vecteur Ge, qui est une représentation du canal vu en sortie de ladite branche de diversité, (b) l'estimateur de canal (E) traite un nombre fini (D) d'itérations selon un algorithme d'estimation-maximisation (EM) fondé sur le critère de probabilité maximum a posteriori (MAP), l'estimateur étant initialement mis en oeuvre en prenant en compte les symboles pilotes contenus dans le bloc bidimensonnel temps- fréquence considéré, ce qui conduit à une estimation d'ordre 0 (#0), l'estimateur prenant en compte, ensuite, les symboles estimés pour les autres itérations et ainsi de suite (3d) , l'estimateur délivrant finalement après une dernière itération D, les coefficients optimaux Gk#(D) (k de 0 à (N-l) et de 0 à L-l) définissant, pour chaque branche de diversité, le vecteur G# représentant le canal. 2. Récepteur selon la revendication 1, dans lequel l'estimateur de canal (E) est apte à calculer une estimation souple des symboles de données calculée à partir des probabilités P(A#(k)=A#{R# }#=0L-1, {G# (d)}#=0L-1) obtenues à partir de la recombinaison des contributions des L branches de diversité, ces contributions étant égales au produit des composantes R#(k)# du vecteur signal de sortie R du filtre adapté de chaque branche de diversité par le complexe conjugué C#(k)#* de l'estimation de la #(k) ème composante du canal obtenue après la dernière itération, (où 8(k) est une focntion d'indexation bidimensionnelle). 3. Récepteur selon la revendication 2, dans lequel ladite sortie souple est soit traitée directement par un organe de décision, si les données ne sont pas codées, soit par un décodeur, si les données sont codées, fournissant le symbole (Â) finalement délivré par le récepteur. 4. Récepteur selon la revendication 1, dans lequel, dans un bloc bidimensionnel, le nombre (Nt) d'intervalles de temps est égal au nombre (Nf) d'intervalles de fréquence. 5. Procédé de réception pour radiocommunications à multiplexage par répartition en fréquences orthogonales (OFDM) dans lequel : i) dans une pluralité de L branches de diversité on traite des blocs de symboles numériques, chaque bloc comprenant des symboles de données et des symboles pilotes répartis dans un bloc bidimensionnel temps-fréquence à Nt intervalles de temps et à Nf intervalles de fréquence, on reçoit le signal dans chaque branche de diversité par un capteur radioélectrique, on produit un signal de sortie à N composantes constituant les composantes d'un vecteur Re, où désigne le rang de la branche de diversité ( allant de 0 à L-l), on estime le canal radioélectrique emprunté par le signal reçu dans un bloc bidimensionnel temps-fréquence, ii) on recombine ~les signaux délivrés par les différentes branches de diversité, iii) on prend une décision à partir du signal de recombinaison et on délivre une estimation des symboles de données, ce procédé étant caractérisé en ce que : (a) on estime le canal dans chaque branche de diversité en traitant un vecteur Ce à N composantes caractérisant le canal dans un bloc bidimensionnel temps-fréquence, on définit une base de N vecteurs Bk qui sont les N vecteurs propres normalisés de la matrice de covariance temps-fréquence du canal, on décompose chaque vecteur Ce dans cette base, les N coefficients de cette décomposition étant notés Gk# avec k allant de 0 à N-l, les coefficients Gk# définissant, pour chaque branche de diversité, un vecteur Ge, qui est une représentation du canal vu en sortie de ladite branche de diversité, (b) pour effectuer l'estimation du canal on effectue un nombre fini (D) d'itérations selon un algorithme d'estimation-maximisation (E. M.) fondé sur le critère de probabilité maximum a posteriori (MAP), on initialise les itérations en prenant en compte les symboles pilotes contenus dans le bloc bidimensionnel temps- fréquence considéré, ce qui conduit à une estimation d'ordre 0, et l'on prend en compte, ensuite, les symboles de donnée pour les autres itérations et ainsi de suite, pour obtenir finalement, après une dernière itération D, les coefficients optimaux Gk#(D) (k de 0 à (N-1) et de 0 à L-1) définissant, pour chaque branche de diversité le vecteur Ge représentant le canal. 6. Procédé selon la revendication 5, dans lequel pour effectuer la recombinaison on calcule une variable ##(k)(D) en sommant les contributions des branches de diversité, ces contributions étant égales au produit des composantes R#(k)# du vecteur signal de sortie Re du filtre adapté de la branche de diversité considérée par le complexe conjugué C#(k)#* de la #(k)ième composante du canal obtenue après la dernière itération, (où #(k) est une fonction d'indexation bidimensionnelle). 7. Procédé selon la revendication 6, dans lequel la variable ##(k)(D) est traitée soit par un organe de décision si les données ne sont pas codées, soit par un décodeur si les données sont codées, fournissant les symboles estimé (##(k)) finalement délivrés par le récepteur. 8. Procédé selon la revendication 5, dans lequel on détermine le signe de la variable ##(k)(D), lequel signe constitue le symbole estimé (##(k)) que l'on délivre finalement. 9. Procédé selon la revendication 5, dans lequel, dans un bloc bidimensionnel, le nombre (Nt) d'intervalles de temps est égal au nombre (Nf) d'intervalles de fréquence.