FR 2831359 A1 20030425 FR 0113662 A 20011023 PROCEDE D'ESTIMATION D'UN CANAL DE TRANSMISSION ET DISPOSITIFS POUR SA MISE EN OEUVRE La présente invention concerne un procédé d'estimation d'un canal de transmission, ainsi que des dispositifs pour la mise en oeuvre du procédé. Elle se rapporte au domaine des transmissions numériques par ondes porteuses radiofréquences (transmissions radio-numériques). Elle trouve des applications, notamment, dans les récepteurs des systèmes de radiocommunications numériques avec les mobiles, par exemple les systèmes de radiocommunications professionnelles (systèmes PMR, de l'anglais "Professional Mobile Radio"). Dans ces systèmes, les données numériques sont transmises par modulation d'une onde porteuse radiofréquence. Dit autrement, un signal radio est émis sur le canal de transmission, ce signal étant modulé pour porter les informations numériques à transmettre. Par estimation du canal de transmission, on entend de façon classique l'estimation des conditions de propagation du signal radio à travers ce dernier, qui affectent le signal transmis. On cherche à mettre en oeuvre des techniques de modulation qui offrent une meilleure résistance vis-à-vis des perturbations subies par le signal radio lors de sa transmission à travers le canal de transmission. Pour l'essentiel, ces perturbations proviennent : - d'une part du phénomène d'évanouissement ("Fading", en anglais), qui est sélectif en fréquence dès lors qu'on dépasse la bande de cohérence (l'on parle dans ce premier cas de fading sélectif ( selective fading en anglais)), mais qui n'est pas sélectif en fréquence dès lors que la largeur du canal est inférieure à la bande de cohérence (l'on parle dans ce dernier cas de fading plat ( flat fading en anglais)). Ce phénomène d'évanouissement est du aux trajets multiples de propagation ("Multipaths", en anglais) qui engendrent des interférences inter-symboles (ISI, de l'anglais Intersymbol Interference, ) aussi connues sous le vocable de distorsion intersymboles ; - d'autre part du fait que l'amplitude et la phase du ou de chacun des trajets de propagation peuvent être statiques (en ce sens qu'elles ne varient pas au cours du temps) ou au contraire dynamiques (lorsque les conditions de propagation varient au cours du temps). Dans le cas dynamique, la fréquence de ce phénomène (appelé aussi fréquence du fading) et, d'une façon plus générale, le spectre fréquentiel du fading sont liés à la vitesse du mobile et à la fréquence porteuse du signal émis. Le modèle classique retenu pour le spectre de puissance du fading est décrit dans l'ouvrage"Microwave Mobile Communications", par William C. Jakes, Jr., Editions John Wiley & Sons, 1974, pp 19-25 ) et fait intervenir la fréquence Doppler fD donnée par : où V est la vitesse du mobile, c est la vitesse de la lumière, et fc est la fréquence de la porteuse radiofréquence. On cherche actuellement à mettre en oeuvre une modulation multiporteuse appelée OFDM (de l'anglais"Orthogonal Frequency Division Multiplexing"). Cette technique de modulation a été retenue pour la norme européenne sur les systèmes de diffusion audionumérique (systèmes DAB, de l'anglais"Digital Audio Broadcasting"). Elle consiste à répartir les données à transmettre sur un ensemble de sous-porteuses ("subcarriers"en anglais) émises en parallèle dans le signal radio. Il en résulte un effet d'évanouissement plat vis-à-vis de chaque sous-porteuse car la largeur de bande de chaque sous-porteuse est inférieure à la bande de cohérence. En outre, il en résulte une réduction de la sensibilité de la transmission vis-à-vis du phénomène de trajets multiples. Le signal à transmettre est construit sur un réseau temps-fréquence (aussi appelé"lattice"en anglais). Un tel réseau temps-fréquence comprend un ensemble de symboles, constituant un espace à deux dimensions qui est défini par un axe des fréquences (ou axe fréquentiel) et par un axe des temps (ou axe temporel). On rappelle qu'un symbole correspond à un nombre déterminé de bits d'information, par exemple huit bits, qui prend une valeur déterminée dans un alphabet ad-hoc. Par convention, l'axe des fréquences est représenté verticalement et l'axe des temps est représenté horizontalement. Chaque symbole est repéré par un indice m suivant l'axe des fréquence, et par un indice n suivant l'axe des temps. Par convention, un symbole dont la position suivant l'axe des fréquences est définie par l'indice m, et la position suivant l'axe des temps est définie par l'indice n est en général noté Sm, n. Enfin, l'espacement entre les symboles suivant l'axe des fréquence est noté yo. De même, l'espacement entre les symboles suivant l'axe des temps est noté T0. Si l'on note S (t) un signal construit sur un tel réseau de symboles, on peut décomposer le signal S (t) sous la forme : où le signe, désigne l'opération de sommation ; où les coefficients cm, nsont des coefficients correspondant à la valeur du symbole Sm, n ; et, où la fonction g (t) désigne l'impulsion de mise en forme pour la modulation Le signal à transmettre est structuré en trames qui sont transmises successivement à travers le canal de transmission. Chaque trame comprend un nombre M de sous-porteuses adjacentes à l'intérieur d'un canal de largeur spectrale déterminée, chacune de ces sous-porteuses étant divisée en N intervalles de temps, appelés temps symboles, qui sont transmis successivement à travers le canal de transmission. La durée d'un temps symbole correspond à la durée de transmission d'un symbole. Une trame du signal comprend donc MxN symboles. Le paramètre yo précité représente l'espacement entre deux sous-porteuses adjacentes, et le paramètre ro précité représente l'espacement entre deux symboles successifs sur une même sous-porteuse. Dans les systèmes utilisant une modulation de type OFDM, les impulsions de mise en forme pour la modulation sont choisies de façon que chaque symbole soit orthogonal avec tous les autres symboles. On dit alors que le réseau est orthogonal. Par définition, des symboles sont orthogonaux entre eux si leur produit scalaire est nul. Cette caractéristique permet de simplifier la démodulation. Les systèmes utilisant une modulation de type OFDM se subdivisent en deux catégories. D'une part, les systèmes utilisant un réseau temps-fréquence de densité 1 (dits dans la suite"systèmes de densité 1", par raccourci) pour lesquels le produit yo xro est égal à l'unité (10 x "Co = 1). Dans ces systèmes, les symboles modulés peuvent être des symboles complexes. Les coefficients cm, n précités sont alors des nombres complexes. On peut écrire Cm, n = am, n + i x bm, n, où am, n et bm, n sont des nombres réels. Ceci offre la possibilité d'employer à la fois une modulation en amplitude et une modulation en phase. En pratique, une garde doit toutefois être prise dans le domaine fréquentiel et/ou dans le domaine temporel entre deux symboles adjacents, suivant l'axe fréquentiel, respectivement suivant l'axe temporel consécutifs. Cette garde réduit sensiblement le débit maximal (exprimé en nombre de symboles par secondes, ou bauds) qui peut être écoulé à travers le canal de transmission. D'autre part, les systèmes utilisant un réseau temps-fréquence de densité 2 (dits dans la suite"systèmes de densité 2", par raccourci) pour lesquels le produit y pro est égal à (yoo=-)-Dans ces systèmes, le 2 2 débit maximal (exprimé en nombre de symboles par secondes, ou bauds) est deux fois plus élevé que dans les systèmes de densité 1. Par contre, dans les systèmes de densité 2, les symboles modulés doivent être monodimensionnels, c'est à dire qu'ils ont soit une valeur réelle (on parle alors de symboles réels), soit une valeur imaginaire pure (on parle alors de symboles imaginaires purs). On peut écrire cm, n =am, n pour les symboles réels ou Cmji =ixbm, n pour les symboles imaginaires purs, où am, n et bm, n sont des nombres réels. Plus précisément, si un symbole est réel, ses voisins immédiats, c'est-à-dire les symboles situés sur la même sous-porteuse dans les temps symboles immédiatement antérieur et immédiatement postérieur (par référence à l'ordre d'émission des symboles sur le canal de transmission, c'est- à-dire les symboles adjacents suivant la direction de l'axe des temps) et les symboles situés dans le même temps symbole sur les sous-porteuses placées sur les fréquences immédiatement supérieure et immédiatement inférieure (i. e., les symboles adjacents suivant la direction de l'axe des fréquences), sont imaginaires purs. A l'inverse, si un symbole est imaginaire pur, ses voisins immédiats, c'est à dire les symboles adjacents suivant la direction de l'axe des fréquences et les symboles adjacents suivant la direction de l'axe des temps, sont réels. Les systèmes de densité 2 ne nécessitent pas la présence d'une garde en fréquence ou en temps. Ils permettent donc de transporter un débit plus élevé que les systèmes de densité 1. Dans la suite, on considèrera uniquement le cas des systèmes de densité 2. L'invention s'applique en effet aux systèmes de ce type. Un exemple particulier d'une modulation de type OFDM dans un système de densité 2, est la modulation dite OFDM/IOTA (de l'anglais "OFDM/lsotropic Orthogonal Transform Algorithm"). La façon dont on peut définir un réseau temps-fréquence orthogonal avec une telle modulation est par exemple décrite dans l'article"Coded Orthogonal Frequency Division Multiplex", Bernard LE FLOCH et al., Proceedings of the IEEE, Vol. 83, No. 6, Juin 1995). Les coefficients cm, n sont alors des nombres soit réels, soit imaginaires purs, selon le placement du symbole Sm, n dans la trame. Ils sont donc toujours mono-dimensionnels. Ceci n'offre que la possibilité d'une modulation en amplitude. Néanmoins, il n'est pas nécessaire de garantir un temps de garde entre les symboles ou entre les sous-porteuses, ce qui présente l'avantage d'augmenter le débit de la transmission. Les symboles transmis sont donc pour moitié réels et pour moitié imaginaires purs. Ces symboles sont mis en forme par l'impulsion de modulation g (t), mentionnée plus haut, . Cette impulsion s'étend sur l'axe temporel, sur une durée correspondant à plusieurs symboles, et/ou sur l'axe fréquentiel, sur des fréquences correspondant à plusieurs sous-porteuses. A la réception d'un signal radio, on effectue une synchronisation temporelle et fréquentielle du signal reçu. Puis le signal reçu est corrélé avec le signal attendu, c'est-à-dire qu'on effectue une corrélation du signal reçu avec l'impulsion de modulation g (t). Cette corrélation peut être effectuée par différentes méthodes, par exemple en effectuant une multiplication par l'impulsion de modulation g (t) puis une FFT (de l'anglais"Fast Fourier Transform", qui signifie transformation de Fourier rapide). Il convient ensuite de procéder à l'estimation des conditions de propagation sur le canal de transmission, c'est-à-dire l'estimation du canal de transmission, aussi appelée estimation du fading car elle produit une valeur estimée de l'évanouissement (fading) du signal transmis à travers le canal de transmission). En effet, ces conditions de propagation sont à prendre en compte pour démoduler le signal reçu, et plus précisément pour estimer la valeur des symboles transmis. L'objet de la présente invention est de proposer une méthode d'estimation du canal convenant pour les systèmes de densité 2, c'est-à-dire les systèmes utilisant une modulation OFDM basée sur un réseau tempsfréquence de densité 2. Selon un premier aspect de l'invention, il est en effet proposé un procédé d'estimation d'un canal de transmission à partir d'un signal reçu après transmission à travers ledit canal de transmission, ledit signal étant un signal multiporteuses construit sur un réseau temps-fréquence défini par un axe des fréquences et un axe des temps, et comprenant des trames ayant MxN symboles répartis sur M sous-porteuses dont chacune est divisée en N temps symboles déterminés, chaque trame comprenant P symboles pilotes répartis en temps et en fréquence de manière à couvrir la trame selon une structure maillée, où les nombres M, N et P sont des nombres entiers non nuls, les symboles pilotes comprenant d'une part des symboles dits symboles pilotes réels, transmis comme des symboles ayant une valeur réelle, et d'autre part des symboles dits symboles pilotes imaginaires purs, transmis comme des symboles ayant une valeur imaginaire pure, le procédé comprenant les étapes consistant à : a) sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs zl du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique (en module en en phase) pour ces symboles pilotes ; b) déterminer des nombres complexes u et v et un nombre réel ^ en minimisant l'expression aux moindres carrés suivante : où le signe E désigne l'opérateur sommation, où Il xj) désigne l'opérateur valeur absolue de la variable réelle x ou le module de la variable complexe x, où Re (w) désigne l'opérateur partie réelle du nombre complexe w, où zu est un nombre réel, et où u et v sont des nombres complexes orthogonaux (c'est-à-dire que Re (u*. v) =0, où w* désigne le nombre complexe conjugué du nombre complexe w), tels que Il vil, c) déterminer une valeur estimée â de la valeur a de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en calculant : â = ./u. L'invention permet donc d'estimer des valeurs du fading dans un système de dimension 2. Les étapes du procédé sont renouvelées en sélectionnant d'autres paires ou groupes de symboles pilotes, de manière à produire suffisamment de valeurs estimées du fading pour permettre la poursuite du canal ("Channel Tracking"). Selon un second aspect de l'invention, il est aussi proposé un dispositif comprenant des moyens pour la mise en oeuvre de ce procédé. Le dispositif comprend : - des moyens pour sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives ck d'une part, et une ou plusieurs valeurs zl du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique (en module en en phase) pour ces symboles pilotes ; - des moyens pour déterminer les nombres complexes u et v et le nombre réel X minimisant l'expression aux moindres carrés suivante : où le signe désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où u et v sont des nombres complexes orthogonaux (c'est à dire tels que Re (u*. v) =0) tels que ##u## = ##v##, et où # est un nombre réel, - et des moyens pour déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en calculant : â =/u. Selon un premier mode de mise en oeuvre du procédé, ? L est égal à l'unité, u est égal à ss, et v est égal à - i. ss, où # désigne l'inverse de a, en sorte que l'étape b) consiste à déterminer des nombres réels Re (p) et Im (p) qui minimisent l'expression aux moindres carrés suivante : où le signe L désigne l'opérateur sommation, où Il déigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire de la variable complexe x ; et en sorte que l'étape c) consiste à déterminer la valeur estimée â de la valeur a de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en inversant le nombre complexe Re (p) +i-lm (p). Selon un troisième aspect de l'invention, il est aussi proposé un dispositif pour la mise en oeuvre du procédé selon ce premier mode de mise en oeuvre. Le dispositif comprend : - des moyens pour sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs zI du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique pour ces symboles pilotes ; - des moyens pour déterminer des nombres réels Re (ss) et Im (ss) minimisant l'expression aux moindres carrés suivante : où le signe E désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x, et où # désigne l'inverse de a ; et, - des moyens pour déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en inversant le nombre complexe Re (#) +i#lm(#). Selon un second mode de mise en oeuvre du procédé, qui constitue un mode préféré, Â. est égal à p, u est égal à e'', et v est égal à-i-e'', où p et (p sont des nombres réels qui désignent respectivement le module et la phase de a (a=p-e'), en sorte que l'étape b) et l'étape c) sont réalisées conjointement et consistent à déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, définie par a=p-e' où le signe désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x. Ce mode de mise en oeuvre est préféré car il permet d'obtenir directement la valeur de a, en ce sens qu'Il ne comprend pas d'étape finale d'inversion d'un nombre complexe. Il est donc plus rapide. Selon un quatrième aspect de l'invention, il est enfin proposé un dispositif pour la mise en oeuvre du procédé selon ce second mode de mise en oeuvre. Le dispositif comprend : - des moyens pour sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs zl du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique pour ces symboles pilotes ; et, - des moyens pour déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, définie par a=p-e' où le signe E désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x. D'autres caractéristiques et avantages de l'invention apparaîtront encore à la lecture de la description qui va suivre. Celle-ci est purement illustrative et doit être lue en regard des dessins annexés sur lesquels on a représenté : - à la figure 1 : un diagramme illustrant un réseau temps-fréquence sur lequel est construit le signal transmis sur le canal de transmission ; - à la figure 2 : un diagramme illustrant la structure d'une trame d'un signal multiporteuses selon un exemple de réalisation de l'invention ; -à la figure 3 : un diagramme illustrant les étapes d'un procédé de démodulation d'un signal radio reçu par un récepteur après transmission à travers un canal de transmission ; - aux figures 4a à 4c : des diagrammes d'étapes illustrant le procédé selon l'invention respectivement dans le cas général, selon un premier mode de mise en oeuvre et selon un second mode de mise en oeuvre ; - à la figure 5 : un schéma synoptique illustrant des dispositifs selon l'invention. Le diagramme de la figure 1 illustre un réseau temps-fréquence d'un système de densité 2, tel que par exemple un système utilisant une modulation OFDM/IOTA. Le réseau est défini par un axe fréquentiel (ici, l'axe vertical) et par un axe temporel (ici, l'axe horizontal). Le réseau comprend un ensemble de symboles représentés symboliquement ici par des petites flèches horizontales ou verticales. L'espacement entre les symboles suivant l'axe des fréquence est noté 10. De même, l'espacement entre les symboles suivant l'axe des temps est noté To. Selon des propriétés intrinsèques du réseau temps-fréquence, chaque symbole est orthogonal avec tous les autres symboles. A la figure 1, les symboles représentés par des flèches horizontales sont des symboles réels. Ceux représentés par des flèches verticales sont des symboles imaginaires purs. Comme il a été dit en introduction, si un symbole déterminé est réel, ses voisins immédiats, c'est à dire les symboles situés directement à sa droite ou à sa gauche suivant la direction de l'axe des temps et les symboles situés directement au dessus ou au dessous de lui suivant la direction de l'axe des fréquences, sont imaginaires purs. A l'inverse, si un symbole déterminé est imaginaire pur, ses voisins immédiats (définis de la même façon que ci-dessus) sont réels. Par exemple le symbole situé à l'intersection de l'axe temporel et de l'axe fréquentiel (appelé symbole utile) est un symbole réel. Un signal multiporteuse peut être construit sur un tel réseau tempsfréquence, en étant structuré en trames successives transmise sur le canal de transmission. Une trame est définie suivant l'axe des fréquence et suivant l'axe des temps, respectivement par une bande de fréquence B, et par une durée D. Elle comprend M sous-porteuses, où M est un nombre entier tel que B = M xyo. De plus, chaque sous-porteuse est divisée en N temps symboles, où N est un nombre entier tel que D =NxTo. La trame comprend donc MxN symboles. Le diagramme de la figure 2 représente la structure d'une trame d'un signal multiporteuse selon un exemple convenant pour la mise en oeuvre de l'invention. Par convention, on définit une double relation d'ordre pour repérer l'emplacement d'un symbole dans la trame suivant l'axe des fréquence d'une part, et suivant l'axe des temps d'autre part. Selon cette relation d'ordre, le symbole S1, 1 est le symbole qui est porté sur la première sous-porteuse (celle correspondant à l'indice m égal à l'unité (m=1)) et qui est transmis en premier sur cette sous-porteuse c'est-à-dire qui est situé dans le premier temps symbole (celui correspondant à l'indice n égal à l'unité (n=1)). Ce symbole S1, 1 est représenté en bas à gauche sur la figure. De même, le symbole SM, N est le symbole qui est porté sur la dernière sous-porteuse de la trame (celle correspondant à l'indice m égal à M (m=M)) et qui est transmis en dernier sur cette sous-porteuse c'est-à-dire qui est situé dans le dernier temps symbole (celui correspondant à l'indice n égal à N (n=N)). Ce symbole SM, N est représenté en haut à droite sur la figure. De façon générale, le symbole Sm, n est le symbole qui est porté sur la m-ième sous-porteuse de la trame (celle d'indice m) et qui est transmis sur cette sous-porteuse dans le n-ième temps symbole (celui d'indice n). Afin de permettre la poursuite du canal, la trame contient P symboles pilotes, où P est un nombre entier en principe très inférieur à MxN. On rappelle qu'un symbole pilote est un symbole dont l'emplacement dans la trame et dont la valeur sont connus du récepteur. Les symboles pilotes sont répartis en temps et en fréquence, de manière à couvrir la trame selon une structure maillée. Dans l'exemple représenté, le signal occupe une bande de fréquence B = 44 kHz (kilohertz) à l'intérieur d'un canal radio de 50 kHz de large. De plus l'espacement entre les sous-porteuses est yo =2 kHz. La trame comporte donc M = 22 sous-porteuses. En outre la durée de la trame est D = 20 ms (milliseconde). Le débit sur chaque sous-porteuse est de 4 kilosymboles/s (milliers de symboles par seconde), donc l'espacement temporel entre les symboles est To = 250 us. Dit autrement, la trame comporte N = 80 temps symboles. La trame comporte donc MxN = 1760 symboles. Sur la figure, les symboles pilotes sont représentés par des cellules grises, et les autres symboles, qui correspondent à de l'information utile, sont représentés par des cellules blanches. Parmi les 1760 symboles de la trame, il y a 206 symboles qui sont des symboles pilotes. Dit autrement, P=206. Certains des symboles pilotes, qui sont adjacents deux à deux suivant la direction de l'axe des fréquences et/ou suivant la direction de l'axe des temps, forment un bloc de symboles pilotes tels que 51 ou 53. Dans l'exemple, la trame comprend en effet des blocs de symboles pilotes. Un bloc de symboles pilotes est défini en ce sens qu'il s'agit d'un groupe de symboles pilotes, adjacents ou non suivant la direction de l'axe des fréquences et/ou suivant la direction de l'axe des temps, pour lesquels est satisfaite une double condition de stationnarité en temps et de stationnarité en fréquence du canal de transmission. Par convention, dans la suite, la position d'un bloc de symboles pilotes dans la trame est repérée par la position du symboles pilote de ce bloc qui est sur la porteuse de plus petit indice, et dans le temps symbole transmis en premier (sur les figures, c'est, pour chaque bloc, le symbole pilote qui est le plus bas et le plus à gauche). De même, la taille du bloc est définie par une dimension suivant l'axe des fréquences (ci-après"hauteur") exprimée en nombres de symboles, et par une dimension suivant l'axe des temps (ci-après "longueur"), exprimée en nombres de symboles. La taille du bloc est notée hxl, où h désigne la hauteur et 1 désigne la longueur du bloc. Cette convention est commode dans les cas où les blocs de symboles pilotes ont des dimensions régulières (formant pas exemples des lignes, ou des pavés de symboles pilotes c'est-à-dire des carrés ou des rectangles), comme c'est le cas dans l'exemple représenté. Néanmoins, il est bien entendu qu'un bloc de symboles pilotes peut avoir une structure irrégulière (par exemple trois symboles pilotes adjacents deux à deux mais non alignés). En outre, il est précisé que la notion de bloc de symboles pilotes selon l'invention ne correspond pas forcément à une notion d'adjacence mais plutôt à une notion de proximité à la fois suivant la direction de l'axe des fréquences et suivant la direction de l'axe des temps. En réalité, la définition d'un bloc de symboles pilotes est la suivante : les symboles pilotes d'un même bloc, qui peuvent être adjacents ou non, sont considérés comme des symboles satisfaisant une double condition de stationnarité en fréquence et de stationnarité en temps des conditions de propagation sur le canal de transmission. Ces deux conditions peuvent se traduire en termes d'espacement maximum des symboles pilotes, respectivement suivant la direction de l'axe fréquentiel et suivant la direction de l'axe temporel, ainsi qu'il va être explicité au paragraphe suivant. Il en résulte que les dimensions maximales d'un bloc de symboles pilotes selon l'invention dépendent des caractéristiques de la propagation, celles-ci devant donc être prises en compte lors du choix de la répartition des symboles pilotes dans la trame par le concepteur du système. On sait que les caractéristiques de la propagation à travers le canal de transmission sont définies par la fréquence maximum des variations de l'évanouissement (appelée"fréquence de fading", dans le jargon de l'homme du métier) et le retard maximum entre les trajets multiples. Dans un exemple, pour une propagation de type HT ("Hilly Terrain") qui est la plus contraignante, il est établi que la fréquence de fading est égale à 148,2 Hz (hertz) pour une vitesse maximum de déplacement des mobiles égale à 200 Km/h (kilomètre par heure) et pour une fréquence de porteuse égale à 400 MHz (Mégahertz), d'une part, et que le retard maximum entre les trajets multiples correspond à 7, 5 us (microseconde), soit un retard maximum entre le trajet le plus en avance et le trajet le plus en retard de 15 ils, d'autre part. Or, pour un débit de 4 kilosymboles/s par sous-porteuse, la trame doit comprendre des symboles pilotes (ou des blocs de pilotes) avec un espacement suivant la direction de l'axe temporel, appelé espacement temporel, qui doit être inférieur à l'inverse de la fréquence du fading, c'est-àdire qu'elle doit comprendre un symbole pilote tous les 27 symboles au maximum. Cet espacement maximum de 27 symboles suivant le direction de l'axe des temps correspond à un échantillonnage du canal de propagation (fading) effectué plus rapidement (même seulement légèrement plus rapidement) que l'occurrence des trous de fading (passage par zéro du fading sur l'axe temporel) successifs. Entre deux trous de fading successifs, la phase du fading a tourné de n (nombre Pl). Sur, par exemple, un dixième de cette période entre trous de fading, c'est-à-dire sur une période correspondant à 2,7 symboles successifs, le fading aura tourné de n/10. Dans la suite on considèrera pour des raisons pratiques un groupe ayant 2 symboles successifs suivant l'axe des temps. Entre ces deux symboles successifs, le fading aura tourné de n/27. Si le fading au milieu de cette période de 2 symboles a une certaine valeur Fm déterminée, le fading à la fin de cette période de 2 symboles aura une valeur Ff qui sera très proche de Fm xe''. Donc on a une erreur quadratique donnée par : soit un rapport signal sur bruit de : 24,71 dB. De même si le fading au début de cette période de 2 symboles a une certaine valeur Fd déterminée, on a au milieu de cette période une erreur quadratique donnée par : soit le même rapport signal à bruit de 24,71 dB. On peut donc considérer sans aucun inconvénient que le canal est stationnaire en temps sur une durée correspondant à deux temps symboles successifs, c'est-à-dire à deux symboles adjacents suivant la direction de l'axe des temps. De même, pour une propagation de type HT, qui présente un retard maximum entre trajets de 15 ils, et pour un espacement entre sous-porteuses de 2 kHz, la trame doit comprendre des symboles pilotes avec un espacement suivant l'axe des fréquences, appelé espacement fréquentiel, qui doit être inférieur à l'inverse du retard maximum entre les trajets multiples, soit un symbole pilote toutes les 33 sous-porteuses au maximum. Cet espacement de 33 sous-porteuses suivant la direction de l'axe des fréquences correspond à un échantillonnage du canal en fréquentiel plus fréquemment (même seulement légèrement plus fréquemment) que l'occurrence des trous de sélectivité fréquentielle (passage par zéro du niveau du signal reçu à certaines fréquences) successifs. Entre deux trous de sélectivité fréquentielle successifs la phase du fading a tourné de n. Sur, par exemple, un dixième de cet espace entre trous de sélectivité fréquentielle, c'est-à-dire sur une bande de fréquence correspondant à 3,3 sous-porteuses, le fading aura tourné de nit 0. On considèrera, pour des raisons pratiques, une bande correspondant à 3 sous-porteuses. Entre les sous-porteuses extrêmes de cette bande, le fading aura tourné de (n/33) x 2. Si le fading au milieu de la bande de fréquence comprenant ces 3 sous-porteuses a une certaine valeur F déterminée, le fading pour la sous-porteuse la plus élevée en fréquence de ce groupe de 3 sous-porteuses aura une valeur Ff qui sera très proche de Fxe''3). Donc on a une erreur quadratique donnée par : soit un rapport signal sur bruit de 20,43 dB. De même si le fading pour la sous-porteuse la moins élevée en fréquence de ce groupe de 3 sous-porteuses a une valeur Fd déterminée, on a en milieu de cette période une erreur quadratique donnée par : soit le même rapport signal à bruit de 20,43 dB. On peut donc considérer sans aucun inconvénient que le canal est stationnaire en fréquence sur une bande de fréquence correspondant trois sous-porteuses adjacentes, c'est-à-dire à trois symboles adjacents suivant la direction de l'axe des fréquences. En conséquence de ce qui précède, des symboles pilotes qui ne sont pas espacés de plus de deux symboles suivant l'axe des temps ni de plus de trois symboles suivant l'axe des fréquences peuvent être considérés comme satisfaisant une double condition de stationnarité en temps et en fréquence, des conditions de propagation sur le canal de transmission (i. e. du fading). On notera que les considérations sur la stationnarité du fading exposées ci-dessus, sont à évaluer en fonction du problème à traiter, c'est-à-dire notamment des caractéristiques de la propagation envisagée, de la vitesse du mobile, et de la fréquence porteuse. Dans le cas général, cela revient à dire que, dès lors que des symboles pilotes sont suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps, on peut considérer que l'évanouissement ("fading") du signal à travers le canal de transmission a une valeur complexe sensiblement identique pour ces symboles pilotes. En conclusion, on peut ainsi donner la définition d'un bloc de symboles pilotes au sens de la présente invention : un bloc de symboles pilotes selon l'invention est défini en ce sens qu'il s'agit d'un groupe de symboles pilotes, adjacents ou non suivant la direction de l'axe des temps et/ou suivant la direction de l'axe des fréquences, pour lesquels est satisfaite une double condition de stationnarité en temps et de stationnarité en fréquence des conditions de propagation sur le canal de transmission. On pourra alors faire l'hypothèse que les symboles d'un tel bloc ont été affectés d'un fading ayant une valeur identique (en module et en phase). Dans l'exemple représenté à la figure 2, la trame comprend au moins un bloc de six symboles pilotes de dimensions 3x2, c'est-à-dire dont les dimensions suivant l'axe des fréquences et suivant l'axe des temps correspondent respectivement à trois symboles (h=3) et à deux symboles (1=2). Selon l'exemple, la trame comporte plus exactement 32 blocs tels que 51, ayant six symboles pilotes chacun, dont les dimensions suivant l'axe des fréquences et suivant l'axe des temps correspondent respectivement à trois symboles et à deux symboles. Leurs emplacements respectifs dans la trame, repérés par l'emplacement du symbole pilote du bloc considéré qui est dans la sous-porteuse de fréquence la plus basse et dans le temps symbole transmis en premier (i. e., le symbole le plus en bas et le plus à gauche), sont les emplacements des symboles Sm, n (on rappelle que m et n sont des indices entiers qui repèrent la position du symbole suivant l'axe des fréquences et suivant l'axe des temps respectivement), avec m compris dans l'ensemble {1, 7, 14, 20} et avec n=1+11xj, où j est un nombre entier compris dans l'ensemble [0 ; 7]. La trame comporte en outre un premier bloc supplémentaire 52 de six symboles pilotes dont les dimensions suivant l'axe des fréquences et suivant l'axe des temps correspondent respectivement à trois symboles et à deux symboles. Elle comporte aussi un second bloc supplémentaire 53 de huit symboles pilotes dont les dimensions suivant l'axe des fréquences et suivant l'axe des temps correspondent respectivement à quatre symboles et à deux symboles. Les emplacements respectifs du bloc supplémentaire 52 et du bloc supplémentaire 53 dans la trame, repérés par l'emplacement du symbole pilote du bloc considéré qui est dans la sous-porteuse de fréquence la plus basse et dans le temps symbole transmis en premier (i. e., le symbole le plus en bas et le plus à gauche), sont les emplacements des symboles Sm, n avec le couple (m, n) compris dans l'ensemble de couples { (4, 1), (10, 1)}. Dit autrement, les blocs 52 et 53 sont positionnés sur les symboles S4, 1 et Sio. 1 respectivement. Les blocs de symboles pilotes supplémentaires 52 et 53, en combinaison avec les blocs 51 qui leurs sont adjacents, sont utilisés par le récepteur pour la synchronisation de trame. A la figure 3, on a représenté les étapes principales d'un procédé de démodulation mis en oeuvre par un récepteur d'un système de radiocommunications numériques. Dans un étape 31, le récepteur effectue une synchronisation temporelle et fréquentielle de ses circuits de traitement avec la structure de trame du signal reçu. Cette synchronisation est effectuée au moyen des blocs de symboles pilotes supplémentaires 52 et 53 en combinaison avec les blocs de symboles pilotes 51 qui leur sont adjacents suivant la direction de l'axe des fréquences, comme indiqué plus haut. La description détaillée de cette étape sortirait du cadre du présent exposé. Dans une étape 32, et pour chaque symbole transmis, noté Sp dans la suite (où l'indice p correspond à un couple d'indices m, n pour repérer la position du symbole dans la trame), le signal reçu est corrélé avec le signal attendu, c'est-à-dire qu'on effectue une corrélation du signal reçu avec l'impulsion de modulation g (t). Cette corrélation peut être effectuée par différentes méthodes, par exemple en effectuant une multiplication par l'impulsion de modulation g (t) puis une FFT. Le signal obtenu après cette corrélation, noté zp dans la suite et aux figures, peut s'écrire de la façon suivante : zp =ap-rp (7) où ap et rp sont des nombres complexes qui correspondent respectivement à la valeur du fading, et à la valeur du symbole utile affectée des interférences ayant dégradé le symbole au cours de la transmission à travers le canal de transmission. Du fait de la propriété d'orthogonalité des symboles, le nombre rp comprend en effet le symbole utile d'origine et en outre un terme d'interférences, qui provient des interférences dues à la transmission des symboles voisins. Ces interférences sont, par construction des systèmes OFDM de densité 2, orthogonales au symbole utile d'origine. Ainsi, si le symbole transmis Sp était réel (Le., si cp est un nombre réel égal à ap), on a alors : rp = ap +i. intp (8) où le terme i. intp représente les interférences et est un nombre imaginaire pur (c'est-à-dire que le nombre intp est un nombre réel). A l'inverse, si le symbole transmis était imaginaire pur (i. e., si cp est un nombre imaginaire pur, égal à i-bp), on a alors : rp = intp+ i. bp (9) où le terme intp représente les interférences et est un nombre réel. Dans une étape 33, on réalise ensuite l'estimation du fading pour chacun des symboles pilotes contenus dans la trame. C'est-à-dire qu'on estime les conditions de propagation à travers le canal de transmission pour les symboles pilotes, dont l'emplacement dans la trame et la valeur sont connus à l'avance du récepteur. Dans une étape 34, on procède alors à ce qu'on appelle la poursuite du canal. A cet effet, on effectue une ou plusieurs interpolations, à partir des valeurs estimées obtenues à l'étape 33, pour produire des valeurs estimées du fading pour les autres symboles de la trame (symboles correspondant à de l'information utile). On obtient ainsi une valeur estimée du fading, notée âp dans la suite et aux figures, pour chaque symbole de la trame. Enfin, dans une étape 35, on estime les symboles transmis (notamment les symboles autres que les symboles pilotes, puisque ce sont eux qui portent l'information utile), en effectuant pour chacun le calcul suivant : êp = zp/âp (10) L'invention concerne l'étape 33 ci-dessus, par laquelle on estime le fading pour les symboles pilotes de la trame. L'invention concerne en effet des procédés et des dispositifs d'estimation du canal. Le schéma de la figure 4a illustre les étapes d'un procédé selon un premier aspect de l'invention. Dans une étape 91, on sélectionne d'une part K valeurs du signal reçu, où K est un nombre entier supérieur ou égal à l'unité. Par signal reçu, on entend ici le signal radio reçu par le récepteur, considéré après les étapes de synchronisation 31 et de corrélation 32. Les K valeurs ainsi sélectionnées correspondent à un ou plusieurs symboles pilotes réels Sk. On note Zk ces K valeurs, où k est un indice entier compris entre 1 et K (1 s k s K) du signal reçu. De plus, on note Ck les valeurs respectives des symboles pilotes réels Sk transmis. D'autre part, on sélectionne également L valeurs du signal reçu (au sens indiqué ci-dessus). Ces valeurs sont notées ZI, et correspondent respectivement à un ou plusieurs symboles pilotes imaginaires purs notés SI, de valeurs respectives notées cl, où L est un nombre entier supérieur ou égal à l'unité et où 1 est un indice compris entre 1 et L (1 s I s L). Les symboles pilotes Sk et SI ne sont pas pris au hasard. Au contraire, il s'agit de symboles pilotes qui sont suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe sensiblement identique (en module en en phase) pour ces symboles pilotes. On note a cette valeur complexe. De plus, on note l'inverse de ce nombre complexe (c'est-à-dire P=Va). Dans une étape 92, on détermine ensuite des nombres complexes u et v et un nombre réel À qui minimisent l'expression aux moindres carrés suivante : où le signe L désigne l'opérateur sommation, où Il xii désigne l'opérateur valeur absolue de la variable réelle x ou le module de la variable complexe x, où dz est un nombre réel, où u et v sont des nombres complexes orthogonaux (c'est à dire que Re (u*. v) =0, où Re (x) désigne l'opérateur partie réelle du nombre complexe x, et où x* désigne le nombre complexe conjugué du nombre complexe x), tels quejuHH !. Par exemple, on pourra prendre v =-i-u, où i désigne la racine carrée - iV du nombre entier relatif-1, c'est-à-dire le nombre complexe e i 74. Dans une étape 93, on détermine une valeur estimée ci de la valeur a de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, c'est-à-dire pour les symboles pilotes Sk et SI sélectionnés à l'étape 91, en calculant : a=/u (12) Cette valeur estimée ci vaut pour les symboles pilotes Sk et SI. Bien entendu, les étapes 91 à 93 sont préférentiellement répétées de manière à produire des valeurs estimées ci de la valeur a de l'évanouissement du signal à travers le canal de transmission pour tous les symboles pilotes de la trame, ou du moins pour tous ceux de ces symboles pilotes qui sont pris en compte pour effectuer la poursuite du canal (étape 34 de la figure 3). L'expression aux moindres carrés (11) ci-dessus se comprendra mieux après la description suivante de deux modes de mise en oeuvre de l'invention, effectuée ci-dessous en regard des figures 4b et 4c, dont elle constitue la généralisation. Un premier mode de mise en oeuvre est décrit ci-dessous en regard du diagramme d'étapes de la figure 4b. Dans ce mode de mise en oeuvre, le procédé comprend une étape de sélection 41, qui est identique à l'étape 91 décrite plus haut, et des étapes 42 et 43, qui correspondent aux étapes respectivement 92 et 93 mentionnées plus haut. Dans un premier temps, supposons que K = L = 1. Dit autrement, supposons que, à l'étape 41, on a sélectionné une seule valeur Z1 du signal reçu correspondant à un seul et unique symbole pilote réel S1 de valeur ci d'une part, et une seule valeur Z2 du signal reçu correspondant à un seul et unique symbole pilote imaginaire pur 82 de valeur C2 d'autre part. Pour le symbole S1 de valeur ci transmis comme symbole réel, on peut noter : Re (zi. ss) = C1 (13) et, pour le symbole 82 de valeur C2 transmis comme symbole imaginaire pur on peut aussi noter : Im (Z2. ss) = C2 (14) où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x. On peut donc poser le système d'équations : Dans une étape de calcul 42, on résout le système d'équations cidessus, pour obtenir Re (p) et Im (p). On obtient ainsi le nombre complexe ss = Re (ss) + i. 1m (ss) qui est donné par : où Z1 et Z2 désignent respectivement le nombre complexe défini par le 1 2 conjugué de Z1 et le nombre complexe défini par le conjugué de Z2. Ce calcul se généralise aux autres cas, c'est-à-dire aux cas où (K, L) * (1, 1), en affectant les indices k aux symboles pilotes (sélectionnés à l'étape 41) qui ont été transmis comme symboles pilotes réels, et les indices 1 aux symboles pilotes (sélectionnés à l'étape 41) qui ont été transmis comme symboles imaginaires purs. L'étape 42 consistent alors à déterminer les nombres réels Re (p) et Im (ss) qui minimisent l'expression aux moindres carrés suivante : où le signe désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x. L'expression aux moindres carrés (17) correspond à l'expression aux moindres carrés (11) du cas général, en posant que # est égal à l'unité, u est égal à ss, et v est égal à - i. ss. Dit autrement, le premier mode de mise en oeuvre du procédé conforme au diagramme d'étapes de la figure 4b, se déduit du cas général en posant que X est égal à l'unité, u est égal à ss, et v est égal à-i. p. Les nombres Re (ss) et Im (p) cherchés sont ceux qui minimisent l'expression aux moindres carrés ci-dessous : Dans un exemple, l'étape 42 peut consister à résoudre le système de deux équations à deux inconnues suivant : On peut montrer que, en résolvant ce système d'équations, on obtient alors : d'une part, et d'autre part, où A=Ck-Re (Zk) + q. ! m (zt) ; k t oùB=Dk) +ER où B = LIZkl + LRe (z, f ; k 1 où C : =-Ck-Re (zk) +C). ! m (zj) ; k 1 où DRe (zk).) m (Zk)-DmM. Re (zl) ; et k 1 oùE=Re (zk) Dk) Le procédé permet donc de prendre en compte un nombre quelconque de symboles pilotes, et donc le cas échéant d'introduire de la diversité. Le mode de mise en oeuvre le plus simple, et donc le moins coûteux en temps de calcul, est toutefois celui où on ne sélectionne à l'étape 41 que deux symboles pilotes seulement (dont un symbole pilote réel et un symbole pilote imaginaire pur). Dans une étape 43, on détermine enfin une valeur estimée â de la valeur a de l'évanouissement du signal à travers le canal de transmission, en inversant le nombre complexe Re (p) +i#lm(#). On peut vérifier que cette étape 43 se déduit de la définition de l'étape 93 du cas général donnée plus haut, avec A égal à l'unité, et u égal à #. Cette valeur estimée â vaut pour les symboles pilotes Sk et SI. Bien entendu, les étapes 41 à 43 sont préférentiellement répétées de manière à produire des valeurs estimées â de la valeur a de l'évanouissement du signal à travers le canal de transmission pour tous les symboles pilotes de la trame, ou du moins pour tous ceux de ces symboles pilotes qui sont pris en compte pour effectuer la poursuite du canal (étape 34 de la figure 3). Dans un second mode de mise en oeuvre du procédé, illustré par le diagramme d'étapes de la figure 4c, on peut déterminer directement une valeur estimée â de l'évanouissement du signal à travers le canal de transmission. Ce mode de mise en oeuvre comprend une étape 61 qui est identique à l'étape 91 décrite ci-dessus, et une étape 62 qui est l'équivalent de l'étape 92 décrite ci-dessus. L'étape 62 permet néanmoins de s'affranchir de l'étape 93, qui n'a donc pas d'équivalent dans ce mode de mise en oeuvre. Selon cette variante en effet, on effectue l'étape 92 en posant que # est égal à p, u est égal à à e'', et v est égal à-i-e'', où p et (p sont des nombres réels qui désignent respectivement le module et la phase de â (â =p. e'). De la sorte l'étape 92 et l'étape 93 sont réalisées conjointement, dans une étape de calcul notée 62, et consistent à déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, définie par â =p-e' où le signe désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x. Dit autrement, on détermine à l'étape 62 une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes Sk et SI concernés, définie en coordonnées polaires par â = p-ei- Cette valeur estimée â vaut pour les symboles pilotes Sk et SI. Bien entendu, les étapes 61 et 62 sont répétées de manière à produire des valeurs estimées â de la valeur a de l'évanouissement du signal à travers le canal de transmission pour tous les symboles pilotes de la trame, ou du moins pour tous ceux de ces symboles pilotes qui sont pris en compte pour effectuer la poursuite du canal (étape 34 de la figure 3). Dans un exemple, l'étape 62 peut consister à résoudre le système d'équations suivant : Par calcul matriciel, on obtient alors une expression de tg (ça) d'une part, et de p sous la forme d'un fonction de (p d'autre part (p = f (p)). Dans le cas le plus simple où, à l'étape 91 on sélectionne la valeur Z1 du signal reçu correspondant à un seul et unique symbole pilote réel de valeur C1 d'une part, et la valeur Z2 du signal reçu correspondant à un seul et unique symbole pilote imaginaire pur de valeur C2 d'autre part, l'étape 62 comprend la résolution du système à deux équations suivant : Les symboles sélectionnés à l'étape 91 (dans le cas général en regard de la figure 4a), à l'étape 41 (dans le premier mode de mise en oeuvre décrit cidessus en regard de la figure 4b) ou à l'étape 61 (dans le second mode de mise en oeuvre décrit ci-dessus en regard de la figure 4c), peuvent appartenir à une même et unique trame. Néanmoins, ils peuvent aussi appartenir à deux trames transmises consécutivement à travers le canal de transmission, du moment que la condition de proximité à la fois suivant la direction de l'axe fréquentiel et suivant la direction de l'axe temporel, qui a été définie ci-dessus, est respectée pour ces symboles. Cette condition est respectée lorsque les symboles pilotes Sk et SI appartiennent à un bloc de symboles pilotes au sens défini plus haut, tel que les blocs 51 de l'exemple de trame représenté à la figure 2. De préférence, les symboles pilotes Sk et SI ainsi sélectionnés appartiennent à un groupe de symboles pilotes qui sont adjacents deux à deux dans la trame, suivant la direction de l'axe des fréquences et/ou suivant la direction de l'axe des temps. C'est le cas notamment des symboles pilotes appartenant aux blocs 51 (dont les dimensions sont égales à 3x2) de l'exemple de trame représenté à la figure 2. La figure 5 montre le schéma d'un dispositif selon l'invention, convenant pour la mise en oeuvre du procédé décrit ci-dessus en regard des figures 4a- 4c. Les moyens pour l'exécution de l'étape 93 de la figure 4a ou de l'étape 43 de la figure 4b étant représentés en traits discontinus à la figure 5, dans la mesure où ces moyens du dispositif n'existent pas en tant que tels dans un dispositif pour la mise en oeuvre du procédé selon le mode de mise en oeuvre conforme au schéma de la figure 4c. Le dispositif décrit est ici intégré dans un récepteur radio, qui comprend une antenne 71 pour recevoir le signal radio transmis à travers le canal de transmission constitué par l'air. Le signal capté par l'antenne 71 est amplifié par un amplificateur de réception 72, puis filtré par un filtre passe-bande 73 centré sur la bande de fréquence du canal radio. Le signal ainsi amplifié et filtré est ensuite fourni en entrée d'un module de synchronisation 74 assurant la synchronisation fréquentielle et temporelle du signal. Le signal en sortie du module 74 est transmis en entrée d'un module d'échantillonnage 75, qui reçoit une fréquence d'échantillonnage Fe délivrée par un oscillateur 76. Ce module 75 produit les valeurs zp précitées. Celles-ci sont stockées dans une mémoire 77 du récepteur. Le dispositif selon l'invention comprend spécifiquement un module de sélection 78, qui effectue l'étape de sélection 91,41 ou 61, à partir des valeurs zp stockées dans la mémoire 77. Il comporte en outre un module de calcul 79, qui effectue l'étape de calcul 92,42 ou 62 à partir des K valeurs Zk et des L valeurs ZI sélectionnées par le module de sélection 78. Ce module 79 délivre le nombre complexe p dans le cas du mode de mise en oeuvre du procédé conforme à la figure 4b, ou la valeur estimée â dans le cas du mode de mise en oeuvre du procédé conforme à la figure 4c. Pour la mise en oeuvre du procédé dans le cas général conforme à la figure 4a, ou selon le mode de mise en oeuvre conforme à la figure 4b, le dispositif comporte en outre un module d'inversion 80. Dans le cas général, le module 80 inverse le nombre complexe u. Selon le mode de mise en oeuvre conforme à la figure 4b, le module 80 inverse le nombre complexe ss fourni par le module de calcul 79, pour délivrer la valeur estimée â. Les modules 78,79 et 80 accèdent à la mémoire 77. Ils peuvent être réalisés sous la forme de modules logiciels, compris les moyens de démodulation du récepteur. REVENDICATIONS 1. Procédé d'estimation d'un canal de transmission à partir d'un signal reçu après transmission à travers ledit canal de transmission, ledit signal étant un signal multiporteuses construit sur un réseau temps-fréquence défini par un axe des fréquences et un axe des temps, et comprenant des trames ayant MxN symboles répartis sur M sous-porteuses dont chacune est divisée en N temps symboles déterminés, chaque trame comprenant P symboles pilotes répartis en temps et en fréquence de manière à couvrir la trame selon une structure maillée, où les nombres M, N et P sont des nombres entiers non nuls, les symboles pilotes comprenant d'une part des symboles dits symboles pilotes réels, transmis comme des symboles ayant une valeur réelle, et d'autre part des symboles dits symboles pilotes imaginaires purs, transmis comme des symboles ayant une valeur imaginaire pure, le procédé comprenant les étapes consistant à : a) sélectionner (91) une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs zl du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique (en module en en phase) pour ces symboles pilotes ; b) déterminer (92) des nombres complexes u et v et un nombre réel À en minimisant l'expression aux moindres carrés suivante : où le signe L désigne l'opérateur sommation, où il xi désigne l'opérateur valeur absolue de la variable réelle x ou le module de la variable complexe x, où X est un nombre réel, où u et v sont des nombres complexes orthogonaux (c'est à dire que Re (u*. v) =0, où Re (w) désigne l'opérateur partie réelle du nombre complexe w, et où w* désigne le nombre complexe conjugué du nombre complexe w), tels que ##u## = ##v##. c) déterminer (93) une valeur estimée â de la valeur a de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en calculant : â = /u. 2. Procédé d'estimation d'un canal de transmission selon la revendication 1, suivant lequel Â est égal à l'unité, u est égal à ss, et v est égal à-ioù désigne l'inverse de a, en sorte que l'étape b) consiste à déterminer (42) des nombres réels Re (p) et Im (ss) en minimisant l'expression aux moindres carrés suivante : où le signe désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire de la variable complexe x, ; et en sorte que l'étape c) consiste à déterminer (43) la valeur estimée a de la valeur a de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en inversant le nombre complexe Re (p) +i. lm (p). 3. Procédé selon la revendication 2, suivant lequel l'étape b) comprend la résolution du système de deux équations à deux inconnues suivant : 4. Procédé selon la revendication 2, suivant lequel l'étape a) comprend la sélection de la valeur Z1 du signal reçu correspondant à un seul et unique symbole pilote réel de valeur C1 d'une part, et de la valeur Z2 du signal reçu correspondant à un seul et unique symbole pilote imaginaire pur de valeur C2 d'autre part, et suivant lequel l'étape b) comprend la résolution du système d'équations à deux inconnues suivant : 5. Procédé selon la revendication 1, suivant lequel # est égal à p, u est égal à e'', et v est égal à -i#e-i.#, où p et (p sont des nombres réels qui désignent respectivement le module et la phase de â (â = p-e''), en sorte que l'étape b) et l'étape c) sont réalisées conjointement et consistent à déterminer (62) une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, définie par â=p#ei.# où p et (p minimisent l'expression aux moindres carrés suivante : où le signe L désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x. 6. Procédé selon la revendication 5, suivant lequel l'étape b) comprend la résolution du système de deux équations à deux inconnues suivant : 7. Procédé selon la revendication 5, suivant lequel l'étape a) comprend la sélection de la valeur Z1 du signal reçu correspondant à un seul et unique symbole pilote réel de valeur ci d'une part, et de la valeur Z2 du signal reçu correspondant à un seul et unique symbole pilote imaginaire pur de valeur C2 d'autre part, et suivant lequel l'étape b) comprend la résolution du système d'équations à deux inconnues suivant : 8. Procédé selon l'une quelconque des revendications 1 à 7, suivant lequel les symboles sélectionnés à l'étape a) appartiennent à une même et unique trame. 9. Procédé selon l'une quelconque des revendications 1 à 7, suivant lequel les symboles sélectionnés à l'étape a) appartiennent à deux trames transmises consécutivement à travers le canal de transmission. 10. Procédé selon l'une quelconque des revendications 1à9, suivant lequel les symboles sélectionnés à l'étape a) appartiennent à un groupe de symboles pilotes qui sont adjacents deux à deux dans la trame, suivant la direction de l'axe des fréquences et/ou suivant la direction de l'axe des temps. 11. Dispositif pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 1 à10, comprenant : - des moyens pour sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs ZI du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe et sensiblement identique pour ces symboles pilotes ; - des moyens pour déterminer des nombres complexes u et v et un nombre réel ^ minimisant l'expression aux moindres carrés suivante : où le signe désigne l'opérateur sommation, où il xi désigne l'opérateur valeur absolue de la variable réelle x ou le module de la variable complexe x, où À est un nombre réel, où u et v sont des nombres complexes orthogonaux (c'est à dire que Re (u*. v) =0, où Re (w) désigne l'opérateur partie réelle du nombre complexe w, et où w* désigne le nombre complexe conjugué du nombre complexe w), tels que lIulJ = Il vii. - des moyens pour déterminer une valeur estimée â de la valeur a de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en calculant : â = À 1 u. 12. Dispositif pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 2 à 4, comprenant - des moyens pour sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs zl du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique pour ces symboles pilotes ; - des moyens pour déterminer des nombres réels Re (ss) et Im (ss) minimisant l'expression aux moindres carrés suivante : où le signe L désigne l'opérateur sommation, où il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x, et où ss désigne l'inverse de a ; et, - des moyens pour déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, en inversant le nombre complexe Re (ss) +i-lm (P). 13. Dispositif pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications 5 à 7, comprenant : - des moyens pour sélectionner une ou plusieurs valeurs Zk du signal reçu correspondant à un ou plusieurs symboles pilotes réels de valeurs respectives Ck d'une part, et une ou plusieurs valeurs zl du signal reçu correspondant respectivement à un ou plusieurs symboles pilotes imaginaires purs de valeurs respectives cl d'autre part, ces symboles pilotes étant suffisamment rapprochés à la fois suivant l'axe des fréquences et suivant l'axe des temps pour qu'il puisse être considéré que l'évanouissement ("fading") du signal à travers le canal de transmission a eu une valeur complexe a sensiblement identique pour ces symboles pilotes ; et, - des moyens pour déterminer une valeur estimée â de l'évanouissement du signal à travers le canal de transmission pour les symboles pilotes concernés, définie par a=p-e'T où p et (p sont des nombres réels qui minimisent l'expression aux moindres carrés suivante : où le signe L désigne l'opérateur sommation, où Il x Il désigne l'opérateur valeur absolue de la variable réelle x où Re (x) désigne l'opérateur partie réelle de la variable complexe x, et où Im (x) désigne l'opérateur partie imaginaire pure de la variable complexe x.