Brevet FR2897999A1 — publié le 20070831

PROCEDE ET DISPOSITIF D'ESTIMATION DE LA FONCTION DE TRANSFERT DU CANAL DE TRANSMISSION POUR DEMODULATEUR COFDM Domaine de l'invention La présente invention concerne un procédé et un dispositif d'estimation de la fonction de transfert du canal de transmission pour un démodulateur dit COFDM ("Coded Orthogonal Frequency Division Multiplex" ou multiplexage par division de fréquence orthogonale codée). Exposé de l'art antérieur Une modulation COFDM peut, par exemple, être mise en oeuvre pour la transmission hertzienne de données vidéo numériques selon le standard DVB (acronyme anglais pour Digital Video Broadcasting). Un tel standard prévoit également la transmission hertzienne de données vidéo numériques vers des récepteurs mobiles (standard DVB-H). La figure 1 est destinée à illustrer le principe d'une modulation COFDM. Des paquets de données à émettre sont mis sous la forme de N coefficients complexes associés à N fréquences (ou porteuses) respectives. Le nombre N des fréquences est égal par exemple à 1705 pour le mode dit "2K" et à 6817 pour le mode dit "8K", en transmission herzienne de télévision numérique. Chaque coefficient complexe correspond à un vecteur qui est illustré en figure 1 comme partant d'un axe des fréquences en un point indiquant la fréquence associée au coefficient. L'ensemble de ces N coefficients est traité par transformée de Fourier rapide inverse (IFFT), ce qui produit un "symbole" constitué d'une somme de porteuses modulées, chaque porteuse ayant une amplitude et une phase déterminées par le coefficient complexe associé. Le symbole ainsi généré est émis. De manière classique, en transmission herzienne, la largeur du canal d'émission est de 6, 7 ou 8 MHz et chaque porteuse est séparée de la suivante par un écart de fréquence Af = 1/Tu. Tu est la durée d'émission d'un symbole et est appelée durée utile. La durée utile est de l'ordre de 224 ps en mode 2K et 896 ps en mode 8K, pour une bande passante de 8MHz. A la réception, un récepteur fait subir au symbole le traitement inverse, c'est-à-dire principalement une transformée de Fourier rapide (FFT) pour reconstituer les coefficients complexes de départ. Comme cela est représenté en figure 1, certains vecteurs P1, P2, P3... régulièrement répartis ont une valeur constante connue. Ces vecteurs, ou les porteuses correspondantes, sont appelés pilotes. Ils servent à refléter les dis-torsions subies par le signal transmis et par l'information qu'ils donnent sur la réponse du canal de transmission, ils permettent de corriger les vecteurs inconnus situés entre les pilotes. La figure 2 représente schématiquement la place de pilotes dans les symboles. Les symboles sont groupés en trames de 68 symboles, de façon classique en transmission hertzienne de télévision numérique (norme ETSI EN 300 744, V1.4.1). En figure 2, chaque ligne représente un symbole et chaque case représente la position d'une porteuse. Les porteuses sont définies comme allant d'une position 0 à une position N-1. Les pilotes, représentés par des disques noirs, sont de deux types. D'une part, il y a, dans chaque symbole, des pilotes continus Pc. Les pilotes continus correspondent à des fréquences particulières réparties dans le canal. Dans la norme ETSI évoquée ci-dessus, il y a en 45 en mode 2K et 177 en mode 8K. Les pilotes continus sont présents dans tous les symboles et occupent toujours la même position en fréquence. En figure 2, seuls les pilotes continus correspondants aux positions 0 et N-1 ont été représentés. D'autre part, il y a, dans chaque symbole, des pilotes Pr dits "répartis" ("scattered pilots" en anglais), qui sont disposés toutes les 12 porteuses, et décalés de trois positions entre deux symboles successifs. Ainsi, tous les quatre symboles, on retrouve la même disposition des pilotes répartis Pr. Les autres porteuses, représentées par des disques blancs, véhiculent le signal utile proprement dit et sont appelées par la suite porteuses utiles. Comme cela apparaît en figure 2, pour des premières positions en fréquence, les porteuses associées correspondent toujours à des porteuses utiles quel que soit le symbole considéré, tandis que pour des secondes positions en fréquence, les porteuses associées correspondent à des porteuses utiles ou à des pilotes répartis selon le symbole considéré. Les premières positions sont appelées par la suite positions de porteuses utiles. Les secondes positions sont appelées par la suite positions de pilotes répartis. Lors de la transmission de symboles, chaque symbole est précédé d'un intervalle de garde qui est généralement une recopie d'une partie de la fin du symbole correspondant. Les intervalles de garde sont souvent définis par une fraction de la durée utile Tu. Des valeurs classiques de l'intervalle de garde sont Tu/32, Tu/16, Tu/8 ou Tu/4. Les intervalles de garde servent à éviter des distorsions de modulation inter-symboles provoquées par un écho lors de la transmission. Si la réponse impulsionnelle du canal est inférieure à la longueur de l'intervalle de garde et si le canal n'évolue pas ou peu pendant la transmission d'un symbole, on peut considérer que le canal de transmission est équivalent à N canaux multiplicatifs parallèles. De ce fait, si on appelle Xn k le coefficient complexe associé à la porteuse de position k du symbole d'indice n et Yn k le coefficient complexe obtenu après l'application de la transformée de Fourrier rapide au symbole reçu, on peut alors écrire : Yn k Hn kXn k + Bn k Equation (1) où Hn k est la fonction de transfert du canal de transmission pour la porteuse de position k du symbole d'indice n, et Bn k est le bruit dû au canal de transmission. Chaque coefficient complexe émis est donc multiplié par la fonction de transfert du canal correspondante qui dépend seulement du symbole et de la porteuse considérés. Au niveau du récepteur, on réalise une estimation du coefficient complexe émis associé à chaque porteuse. Notons Xn k l'estimation du coefficient Xn k . Pour déterminer Xn k il est tout d'abord nécessaire d'estimer la fonction de transfert Hn k du canal pour la porteuse de position k du symbole d'indice n. On appelle Hn k l'estimation de la fonction de transfert Hn k du canal de transmission pour la porteuse de position k du symbole d'indice n, ou estimation de canal de la porteuse de position k du symbole d'indice n. On détermine alors l'estimation Xn k de la façon suivante : Xn,k = s(Hn,k)Yn,k Equation (2) où s est une fonction, par exemple la fonction inverse, l'équation (2) correspondant alors à l'équation suivante : Xn k Yn,k Hn,k Xn k + Bn'k Equation (3) H n,k H n,k H n,k Il est donc nécessaire de déterminer une estimation du canal Hn,k la plus juste possible. Pour un pilote, une façon simple de déterminer l'estimation de canal consiste à diviser le coefficient complexe reçu par le coefficient complexe émis. En effet, le coefficient complexe émis est connu pour les pilotes. Une telle estimation, appelée Hn k, est obtenue de la façon suivante : Hn,k Xn,k ' Hn,k + n'k Equation (4) n,k Xn,k Une telle estimation de canal est appelée estimation de canal bruitée car elle comporte généralement une composante de bruit importante. Après avoir déterminé, pour un symbole donné, les estimations de canal bruitées des pilotes continus et répartis, il est nécessaire de déterminer les estimations de canal des autres porteuses du symbole qui ne correspondent pas à des pilotes continus ou répartis, c'est-à-dire les porteuses utiles. Si la fonction de transfert du canal ne variait sensiblement pas en fonction de la fréquence, les estimations de canal des porteuses utiles d'un symbole pourraient être déterminées en réalisant la moyenne des estimations de canal bruitées des pilotes du symbole. Toutefois, de façon générale, la fonction de transfert du canal varie en fonction de la fréquence et du temps et un procédé adapté de détermination des estimations de canal des porteuses utiles doit être mis en oeuvre. Un exemple de procédé consiste à utiliser un filtre de Wiener à deux dimensions. Dans ce cas, l'estimation de canal d'une porteuse utile d'un symbole donné correspond à une combinaison des estimations de canal de porteuses du même symbole et de plu-sieurs symboles voisins du symbole considéré. Un tel procédé est complexe à mettre en oeuvre. Un procédé plus simple consiste à séparer l'opération de filtrage à deux dimensions en deux opérations de filtrage à une dimension, une première opération de filtrage réalisée, pour une position de porteuse donnée, sur plusieurs symboles, et une seconde opération de filtrage réalisée, pour un symbole donné, sur plusieurs porteuses du symbole. Plus précisément, le procédé de détermination des 30 estimations de canal des porteuses d'un symbole donné comprend alors les deux étapes successives suivantes : (a) une étape d'interpolation temporelle qui consiste, pour chaque porteuse utile du symbole donné à une position de pilote réparti, à déterminer l'estimation de canal de la 35 porteuse par interpolation à partir des estimations de canal bruitées de pilotes répartis d'autres symboles à la même position en fréquence ; et (b) une étape d'interpolation fréquentielle qui consiste, pour le symbole donné, et pour chaque porteuse utile à une position de porteuse utile, à déterminer l'estimation de canal de la porteuse par interpolation à partir des estimations de canal des porteuses du symbole déterminées à l'étape précédente. Les estimations de canal bruitées des pilotes continus et répartis du symbole donné sont généralement déterminées au fur et à mesure du traitement du symbole donné à l'étape (a). Un procédé d'interpolation temporelle est donc mis en oeuvre à l'étape (a), pour un symbole donné, réalisée indépendamment pour chaque position de porteuse en fréquence. Pour le symbole en cours de décodage ou symbole courant, le procédé d'interpolation temporelle consiste à réaliser une combinaison linéaire d'estimations de canal bruitées de porteuses à la même position en fréquence de plusieurs symboles reçus avant le symbole courant, ou symboles passés, et d'un ou de plusieurs symboles reçus après le symbole courant, ou symboles futurs. La figure 3 est une vue de détail d'une partie de la figure 2. On appelle A, B, C, D, des positions successives en fréquence de pilotes répartis et E la position d'un pilote continu. Pour un procédé d'interpolation temporelle pour lequel la détermination de l'estimation de canal d'un pilote utile à une position de pilote réparti utilise une seule estimation de canal bruitée d'un symbole futur, la détermination des estimations de canal de toutes les porteuses utiles à des positions de pilotes répartis d'un symbole courant donné nécessite au moins que les trois symboles futurs suivants soient connus. En d'autres termes, si le dernier symbole reçu par le démodulateur est le symbole d'indice n, le procédé d'interpolation temporelle est mis en oeuvre pour le symbole d'indice n-3 lorsqu'une seule estimation de canal bruitée d'un symbole futur est utilisée, et, de façon plus générale, le symbole d'indice n- 4s+1, où s est un entier, lorsque les estimations bruitées de s symboles futurs sont utilisées. Un tel procédé d'interpolation temporelle est mis en oeuvre de la façon suivante : pour la porteuse à la position E qui correspond à un pilote continu et pour la porteuse à la position B qui correspond à un pilote réparti, l'estimation de canal est, par exemple, égale à l'estimation de canal bruitée ou peut correspondre au résultat d'un filtrage temporel d'estimations de canal bruitée de pilotes continus ou répartis ; et pour les porteuses utiles aux positions A, C et D, l'estimation de canal est déterminée à partir de l'estimation de canal bruitée d'un pilote réparti d'un symbole futur (c'est-à-dire pour les positions A, C et D, de l'estimation de canal bruitée de même position correspondant respectivement au symbole d'indice n, n-2 et n-1) et des estimations de canal bruitées de pilotes répartis de même position de symboles passés. On appelle Un k un vecteur complexe de dimension M+1, dont chaque composante correspond à une estimation de canal bruitée. Le vecteur Un k est donné, pour une porteuse correspondant à un pilote du symbole d'indice n à la position k, par la relation suivante : Hn-4M,k Hn-4(M-1),k Un,k = Equation (5) Hn-4,k H n,k où M est un entier, par exemple égal à 5. De façon classique, les estimations de canal Hn-1 k, Hn_2 k et Hn_3 k de pilotes utiles sont obtenues par combinaison linéaire des composantes du vecteur Un k et sont données par les équations suivantes : Equation (7) Hn-1,k = O Un,k Equation (6) T Hn-2,k = 62 Un,k Hn-3,k = 03Un,k Equation (8) où 01, 02 et 03 sont des vecteurs à coefficients complexes de dimension M+1, appelés par la suite filtres d'interpolation, et où AT, i étant un entier compris entre 1 et 3, correspond à la transposée de O. La figure 4 représente une courbe d'évolution C en traits pointillés de l'amplitude de la fonction de transfert du canal de transmission pour la position en fréquence A. Les disques noirs Hn-24,A, Hn-20,A, Hn-16,A, Hn-12,A, Hn-8,A Hn-4 A et Hn A correspondent aux estimations de canal bruitées de pilotes répartis qui ont été obtenues, par exemple, à partir de l'équation (4). Les disques blancs Hn-1,A, Hn-2,A et Hn-3,A correspondent aux estimations de canal de pilotes utiles qui sont à déterminer. Chaque estimation de canal Hn-1 A , Hn-2 A et Hn-3 A est alors égale à une combinaison linéaire spécifique des six estimations de canal bruitées Hn-20,A, Hn-16,A Hn-12,A Hn-8,A, Fin-4,A et Hn,A. En théorie, les propriétés statistiques du canal de transmission sont sensiblement indépendantes de la fréquence. Ceci signifie que même si l'amplitude de la fonction de transfert du canal peut être différente pour chaque porteuse d'un symbole donné, la valeur absolue du taux de variation de la fonction de transfert est indépendante de la fréquence. Pour un symbole donné, les filtres d'interpolation 01, 02, 03 ne dépendent donc pas de la position en fréquence des porteuses. Les filtres d'interpolation peuvent être déterminés seulement pour les pilotes continus et être utilisés pour la détermination des estimations de canal des autres porteuses du symbole. Par contre, les propriétés statistiques du canal de transmission peuvent varier dans le temps, en particulier lorsque le récepteur est mobile. Il est donc nécessaire de déterminer de nouvelles valeurs des filtres d'interpolation 01, 02 et 03, par exemple, pour chaque symbole. Les filtres d'interpolation 01, 02, 03 peuvent être 35 déterminés par un procédé itératif. De ce fait, si on attribue le rang p à la dernière valeur déterminée du filtre 6i, i étant compris entre 1 et 3, une nouvelle valeur du filtre Oi correspond au rang p+1 et est fournie par la relation suivante : ei(p+1)=f(ei(p)) Equation (9) où f est une fonction d'itération. De nouvelles valeurs du filtre Oi peuvent être déterminées lors de la mise en oeuvre du procédé d'interpolation temporelle pour chaque pilote continu d'un symbole, et ce, par exemple, pour chaque symbole. Une contrainte de la mise en oeuvre d'un procédé itératif est que la fonction d'itération f doit permettre une convergence rapide de la valeur du filtre O. Un inconvénient du procédé d'interpolation temporelle précédemment décrit qui effectue une combinaison linéaire d'estimations de canal bruitées passées et future est que la précision obtenue peut être insuffisante. Ceci peut être le cas lorsque l'amplitude de la fonction de transfert du canal de transmission varie de façon importante dans le temps, ce qui se produit, par exemple, lorsque le récepteur est mobile. Résumé de l'invention La présente invention vise à pallier tout ou partie des inconvénients des procédés de détermination d'estimations de la fonction de transfert du canal de transmission connus. La présente invention vise plus particulièrement un procédé d'interpolation temporelle qui permet de déterminer des estimations de la fonction de transfert du canal de transmission avec une précision améliorée. Selon un autre objet de la présente invention, la présente invention vise un procédé d'interpolation temporelle utilisant des filtres d'interpolation déterminés par un procédé itératif qui permet de faire converger rapidement les composantes des filtres d'interpolation vers des valeurs représentatives des propriétés du canal de transmission. Pour atteindre tout ou partie de ces objets ainsi que d'autres, la présente invention prévoit un procédé de démodu- lation COFDM de symboles successifs d'un signal reçu depuis un canal de transmission, chaque symbole comportant des premières porteuses modulées en phase et/ou en amplitude véhiculant des données qui dépendent du symbole et des secondes porteuses de référence dites pilotes dont les positions en fréquence varient au moins en partie d'un symbole au symbole suivant, au moins une porteuse donnée des secondes porteuses correspondant, pour chaque symbole, au même pilote, le procédé comprenant, pour chaque symbole, une étape de détermination d'une première estimation de la fonction de transfert du canal pour chaque porteuse d'un ensemble de porteuses parmi les premières porteuses dudit symbole telles que, pour les positions en fréquence des porteuses considérées, des symboles différents dudit symbole comprennent des pilotes, correspondant à une combinaison linéaire de secondes estimations de la fonction de transfert du canal déterminées pour des pilotes à la fréquence de ladite porteuse. Les coefficients de la combinaison linéaire sont déterminés de façon itérative, une nouvelle valeur d'au moins un des coefficients étant égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme comprenant le produit entre un pas d'itération et un terme d'erreur égal à la différence entre la deuxième estimation déterminée pour ladite porteuse donnée d'un symbole qui dépend du coefficient et la première estimation déterminée pour ladite porteuse donnée dudit symbole, le pas d'itération étant déterminé de façon itérative. Selon un exemple de réalisation de la présente invention, chaque seconde estimation déterminée pour une porteuse d'un symbole est égale au rapport entre un coefficient reçu associé à la porteuse et un coefficient théorique prédéterminé associé à la porteuse. Selon un exemple de réalisation de la présente invention, la première estimation est déterminée à partir d'au moins la seconde estimation déterminée pour un pilote à la même fréquence que ladite porteuse d'un symbole reçu après ledit symbole. Selon un exemple de réalisation de la présente invention, la nouvelle valeur du coefficient est égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme égal au produit du pas d'itération, du terme d'erreur et du conjugué de la seconde estimation déterminée pour la porteuse donnée d'un symbole qui dépend du coefficient. Selon un exemple de réalisation de la présente invention, le procédé comprend les étapes, pour au moins la porteuse donnée de chaque symbole, consistant à déterminer un vecteur dont au moins certaines composantes correspondent à des secondes estimations déterminées pour la porteuse donnée de symboles reçus avant ledit symbole ; et à déterminer une nouvelle valeur du pas d'itération égale à la somme de la dernière valeur du pas d'itération déterminée et d'un terme égal au produit d'un terme constant et de la partie réelle du produit de la transposée du vecteur, d'une matrice de pondération et du conjugué du terme d'erreur. La présente invention prévoit également un démo- dulateur COFDM destiné à recevoir des symboles successifs d'un signal reçu depuis un canal de transmission, chaque symbole comportant des premières porteuses modulées en phase et/ou en amplitude véhiculant des données qui dépendent du symbole et des secondes porteuses de référence dites pilotes dont les positions en fréquence varient au moins en partie d'un symbole au symbole suivant, au moins une porteuse donnée des secondes porteuses correspondant, pour chaque symbole, au même pilote, le démo- dulateur comprenant un circuit de détermination, pour chaque symbole, d'une première estimation de la fonction de transfert du canal pour chaque porteuse d'un ensemble de porteuses parmi les premières porteuses dudit symbole telles que, pour les positions en fréquence des porteuses considérées, des symboles différents dudit symbole comprennent des pilotes, correspondant à une combinaison linéaire de secondes estimations de la fonc- tion de transfert du canal déterminées pour des pilotes à la fréquence de ladite porteuse. Le démodulateur comprend un circuit de détermination des coefficients de la combinaison linéaire de façon itérative, une nouvelle valeur d'au moins un des coefficients étant égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme comprenant le produit entre un pas d'itération et un terme d'erreur égal à la différence entre la deuxième estimation déterminée pour ladite porteuse donnée d'un symbole qui dépend du coefficient et la première estimation déterminée pour ladite porteuse donnée dudit symbole, le démodulateur comprenant en outre un circuit de détermination du pas d'itération de façon itérative. Selon un exemple de réalisation de la présente invention, le démodulateur comprend un circuit adapté à fournir chaque seconde estimation pour une porteuse d'un symbole comme étant égale au rapport entre un coefficient reçu associé à la porteuse et un coefficient théorique prédéterminé associé à la porteuse. Selon un exemple de réalisation de la présente invention, le circuit de détermination des coefficients de la combinaison linéaire est adapté à fournir la nouvelle valeur du coefficient comme étant égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme égal au produit du pas d'itération, du terme d'erreur et du conjugué de la seconde estimation déterminée pour la porteuse donnée d'un symbole qui dépend du coefficient. Selon un exemple de réalisation de la présente invention, le circuit de détermination du pas d'itération comprend un circuit adapté, pour au moins la porteuse donnée de chaque symbole, à déterminer un vecteur dont au moins certaines composantes correspondent à des secondes estimations déterminées pour la porteuse donnée de symboles reçus avant ledit symbole ; et un circuit adapté à déterminer une nouvelle valeur du pas d'itération égale à la somme de la dernière valeur du pas d'itération déterminée et d'un terme égal au produit d'un terme constant et de la partie réelle du produit de la transposée du vecteur, d'une matrice de pondération et du conjugué du terme d'erreur. Selon un exemple de réalisation de la présente invention, le circuit de détermination du pas d'itération comprend un circuit adapté à déterminer une nouvelle valeur de la matrice de pondération égale à la somme de la dernière valeur de la matrice de pondération déterminée, d'un premier terme égal au produit du terme d'erreur et du conjugué du vecteur et d'un second terme égal au produit de l'opposée du pas d'itération, du conjugué du vecteur, de la transposée du vecteur et de la dernière valeur de la matrice de pondération déterminée. Brève description des dessins Ces objets, caractéristiques et avantages, ainsi que d'autres de la présente invention seront exposés en détail dans la description suivante d'exemples de réalisation particuliers faite à titre non-limitatif en relation avec les figures jointes parmi lesquelles : la figure 1, précédemment décrite, représente des porteuses modulées en phase et en amplitude dans un système de 20 transmission COFDM ; la figure 2, précédemment décrite, représente schématiquement la position de pilotes dans des symboles ; les figures 3 et 4, précédemment décrites, illustrent le principe de détermination d'estimations de canal ; 25 la figure 5 représente un exemple de démodulateur COFDM selon la présente invention ; la figure 6 illustre les étapes de détermination des filtres d'interpolation mises en oeuvre par un premier exemple de procédé d'interpolation temporelle selon l'invention ; 30 la figure 7 illustre les étapes du premier exemple de procédé d'interpolation temporelle selon l'invention ; et la figure 8 illustre les étapes de détermination des filtres d'interpolation mises en oeuvre par un second exemple de procédé d'interpolation temporelle selon l'invention. Description détaillée Par souci de clarté, de mêmes éléments ont été désignés par de mêmes références aux différentes figures. La figure 5 représente un exemple de démodulateur COFDM d'un récepteur selon la présente invention. Une entrée E du démodulateur reçoit un signal IF de fréquence intermédiaire permettant un échantillonnage, par exemple à 36 MHz. Le signal IF correspond au signal reçu après divers changements ou transpositions de fréquence. L'entrée E est couplée à un convertisseur analogique-numérique 10 (ADC) qui numérise le signal d'entrée IF. Le convertisseur analogique-numérique 10 attaque un module 12 de changement de fréquence. Le module 12 fournit un signal sensiblement en bande de base, le spectre du signal en sortie du module 12 étant centré sur une fréquence sensiblement égale à zéro. Le module 12 est couplé à un module 14 permettant d'une part un réglage fin de la fréquence centrale du spectre du signal et, d'autre part, de fournir des échantillons temporels à des instants appropriés au traitement ultérieur. En sortie du module 14, le spectre du signal est centré sur une fréquence égale à 0 et le nombre et la position temporelle des échantillons sont adaptés à la transformation par transformée de Fourier qui a lieu dans le module suivant. Le module 14 est commandé par des liaisons 15 et 15' reliant le module 14 à un module 16 de traitement des pilotes continus et répartis. La sortie du module 14 attaque un module 20 de transformée de Fourier rapide (FFT) qui fournit les fréquences correspondant à un symbole. Le module 20 est piloté par un module 22 qui fournit, par l'intermédiaire d'une liaison 24, un signal de réglage de la fenêtre d'analyse de la transformée de Fourier. La sortie du module 20 est couplée au module 16 qui réalise l'extraction et le traitement des pilotes continus et répartis. Le module 16 fournit sur les liaisons 15 et 15' les signaux destinés à corriger la fréquence centrale du spectre et la fréquence d'échantillonnage du signal. La sortie du module 20 attaque un module 28 dans lequel le signal est corrigé à l'aide d'une estimation de la réponse fréquentielle du canal, c'est-à-dire l'estimation de la fonction de transfert du canal. L'estimation de la fonction de transfert de canal est réalisée dans le module 16 à l'aide des pilotes. Cette estimation est fournie par le module 16 sur une liaison 30, dont une branche 30a est couplée au module 28. En sortie du module 28, le signal comprend les porteuses véhiculant les données. L'estimation de la fonction de transfert du canal, fournie par le module 16, alimente, par l'intermédiaire de la liaison 30 et d'une branche 30b de la liaison 30, un module 32 de transformée de Fourier inverse (IFFT), pour déterminer la réponse impulsionnelle du canal. Le module 32 fournit la réponse impulsionnelle du canal au module 22, pour ajuster dynamiquement le positionnement de la fenêtre d'analyse FFT. Le traitement des porteuses utiles est assuré dans un circuit 40 de traitement et de fourniture de données. Le circuit 40 a une structure classique et peut comporter, comme cela est représenté en figure 5, un module 42 de désentrelacement de symboles, un module 44 dit de "démappage", un module 46 dedésentrelacement de bits, et un module 48 (FEC) de correction d'erreur. La sortie du module 48 constitue la sortie S du circuit 40 et du démodulateur et fournit des données correspondant aux données émises. Pour chaque symbole, le procédé d'interpolation temporelle selon l'invention est mis en oeuvre par le module 16. Le procédé comporte des étapes de détermination de nouvelles valeurs des filtres d'interpolation 6i, i étant compris entre 1 et 3, chaque fois que la porteuse traitée par le module 16 correspond à un pilote continu. La figure 6 illustre les étapes de détermination des 35 filtres d'interpolation 6i(p), i étant compris entre 1 et 3, mises en oeuvre par un premier exemple de procédé d'interpolation temporelle selon l'invention. A l'étape 100, on réalise une initialisation des vecteurs, matrices et paramètres utilisés par le procédé itératif. A titre d'exemple, chaque filtre 0i(0), i étant compris entre 1 et 3, est initialisé de la façon suivante : 0i(0)=(1/6,..., 1/6)T Equation (10) Le procédé itératif utilise un vecteur G de dimension M+1 qui est modifié au cours du procédé en même temps que le filtre 01(p). L'indice temporel p est donc également utilisé pour le vecteur G. Le vecteur G est initialisé de la façon suivante : G(0) = (0,...,0)T Equation (11) Le procédé utilise des paramètres scalaires max et a qui sont respectivement mis, par exemple, à la valeur 0,2 et max/300. Le procédé utilise un paramètre scalaire p qui est modifié au cours du procédé en même temps que le vecteur G. L'indice temporel p est donc également utilisé pour le paramètre p. Le paramètre p est initialisé de la façon suivante : P(0) = Pmax Equation (12) L'indice i est initialement mis à 1. Le procédé se poursuit à l'étape 102 lorsque le module 16 a reçu suffisamment de symboles pour déterminer le vecteur Un k, au moins 4M+1 symboles dans le premier exemple de réalisation. Toutefois, le procédé peut également être mis en oeuvre lorsque moins de 4M+1 symboles sont reçus. Dans ce cas, le vecteur Un k peut être complété avec des valeurs initiales déterminées. A l'étape 102, on détermine une estimation de canal bruitée Hn,k du pilote continu du symbole d'indice n à la position k à partir de l'équation (4). Le procédé se poursuit à l'étape 104. A l'étape 104, on détermine le vecteur Un k à partir de l'équation (5). Le procédé se poursuit à l'étape 106. A l'étape 106, on détermine l'estimation de canal Hn-ik du pilote continu du symbole d'indice n-i à la position k à partir de l'une des équations (6), (7) ou (8) selon la valeur de l'indice i. Le procédé se poursuit à l'étape 108. A l'étape 108, on détermine une erreur d'estimation en k selon l'équation suivante : en,k = Hnùi,k ù Hnùi,k Equation (13) le vecteur Hnùik ayant été déterminé à un cycle précédent du procédé itératif. Le procédé se poursuit à l'étape 110. A l'étape 110, on détermine une nouvelle valeur 0i(p+l) du filtre d'interpolation 0i selon l'équation suivante : 0i(p+1) = Oi(p)+ (p)Un ken,k Equation (14) où Un k correspond au conjugué transposé du vecteur Un k. Le paramètre (p) correspond donc au pas du procédé itératif. Le procédé se poursuit à l'étape 112. A l'étape 112, on détermine si l'indice i est égal à 1. Dans la négative, le procédé se poursuit à l'étape 114. A l'étape 114, on détermine si l'indice i est égal à 3. Dans la négative, le procédé se poursuit à l'étape 116. A l'étape 116, l'indice i est augmenté d'une unité et le procédé recommence à l'étape 106. Si à l'étape 114, l'indice i est égal à 3, le procédé se poursuit à l'étape 118. A l'étape 118, on attend le prochain pilote continu du symbole d'indice n ou le premier pilote continu du symbole d'indice n+1. L'indice i est mis à 1. Le procédé se poursuit à l'étape 102. Si à l'étape 112, i est égal à 1, le procédé se poursuit à l'étape 120. A l'étape 120, le paramètre old est modifié de la façon suivante : Nold =N(P) Equation (15) Le procédé se poursuit à l'étape 122. A l'étape 122, on détermine une nouvelle valeur p(p+1) du paramètre p selon la relation suivante : (p + 1) =I-t(p) + ai(Un kG(p)en,k ) Equation (16) où en,k correspond au conjugué de enk , 9î(Un kG(p)en,k) correspond à la partie réelle de Un kG(p)en k . Le procédé se poursuit à l'étape 124. A l'étape 124, on compare p(p+1) à max. Si (p+l) est supérieur à max, le procédé se poursuit à l'étape 126. A l'étape 126, (p+1) est mis à max. Le procédé se poursuit à l'étape 128. Si à l'étape 124, p(p+1) est inférieur à max, le procédé se poursuit à l'étape 128. A l'étape 128, on détermine si p(p+1) est inférieur à O. Dans l'affirmative, le procédé se poursuit à l'étape 130. A l'étape 130, (p+1) est mis à old• Le procédé se poursuit à l'étape 132. Si à l'étape 128, p(p+1) est supérieur à 0, le procédé se poursuit à l'étape 132. A l'étape 132, on détermine une nouvelle valeur G(p+l) du vecteur G de la façon suivante : G(p+1)=G(p)- (p+l)Un,kUnkG(p)+Un,ken,k Equation (17) Le procédé se poursuit à l'étape 114. Le fait de modifier, aux étapes 120 à 132, le pas p utilisé pour la détermination des filtres d'interpolation 6i, i étant compris entre 1 et 3, permet de faire converger beaucoup plus rapidement le procédé de détermination des filtres d'interpolation Oi que par un procédé dans lequel le pas serait fixe. Le procédé itératif selon l'invention est alors relativement peu sensible aux valeurs initiales utilisées pour les filtres d'interpolation O. La figure 7 représente les étapes du premier exemple de procédé d'interpolation temporelle mise en oeuvre par le module 16. Le procédé va être décrit en relation à la figure 3. On suppose que les symboles jusqu'au symbole d'indice n-1 ont déjà été reçus par le module 16. A l'étape 150, le module 16 reçoit le symbole d'indice n et en extrait les pilotes continus et répartis. La détermination des estimations de canal va donc être réalisée pour le symbole d'indice n-3. On poursuit alors en considérant la porteuse à la position A à l'étape 152. A l'étape 152, comme la porteuse à la position A du symbole d'indice n correspond à un pilote réparti, on détermine l'estimation de canal bruitée Hn A à partir de l'équation (4). La valeur Hn A est alors mémorisée. En outre, pour le symbole d'indice n-3, on détermine l'estimation de canal H1-3A en utilisant 03 selon la relation suivante qui se déduit de l'équation (8) : Hn-3,A = 03TUn,A On poursuit alors en considérant la porteuse à la position B à l'étape 154. A l'étape 154, l'estimation de canal Hn-3 B de la porteuse à la position B est égale à l'estimation de canal 15 bruité Hn-3 B qui a été déterminée auparavant. On poursuit alors en considérant la porteuse à la position C à l'étape 156. A l'étape 156, on détermine l'estimation de canal Hn-3 C en utilisant 01 selon la relation suivante qui se déduit de l'équation (6) : Hn-3,C = 0T 1 Un-2,C On poursuit alors en considérant la l'étape 158. A l'étape 158, on la relation suivante qui se déduit 30 porteuse à la position E est égale à l'estimation de canal bruitée Hn-3 E• De nouvelles valeurs des filtres d'interpolation 01, 02 et 03 sont déterminées par le procédé décrit précédemment en relation à la figure 6. Le présent procédé peut être poursuivi pour toutes les porteuses du symbole n-3 correspondant Equation (18) 20 Equation (19) porteuse à la position D à détermine l'estimation de canal Hn-3 D en utilisant 02 selon 25 de l'équation (7) : Hn-3,D - 02TUn-1,D On poursuit alors en considérant la l'étape 160. A l'étape 160, l'estimation Equation (20) porteuse à la position E à Hn-3,E de canal de la à des pilotes continus, des pilotes répartis ou des porteuses utiles à des positions de pilotes répartis. En effet, on retrouve toujours la succession des étapes 152, 154, 156 et 158 et occasionnellement l'étape 160. Lorsque toutes les porteuses du symbole d'indice n-3 ont été considérées, le procédé continue à l'étape 162. A l'étape 162, le module 16 reçoit le symbole d'indice n+1 et en extrait les pilotes continus et répartis. La détermination des estimations de canal va donc être réalisée pour le symbole d'indice n-2. On poursuit alors en considérant la porteuse à la position A à l'étape 164. A l'étape 164, on détermine l'estimation de canal Hn_2 A en utilisant 02 selon la relation suivante qui se déduit de l'équation (7) fin-2,A - 02 Un,A Equation (21) On poursuit alors en considérant la porteuse à la position B à l'étape 166. A l'étape 166, comme la porteuse à la position B du symbole d'indice n+1 correspond à un pilote réparti, on détermine l'estimation de canal bruitée Hn+1,B à partir de l'équation (4). La valeur Hn+1,B est alors mémorisée. En outre, pour le symbole d'indice n-2, on détermine l'estimation de canal Hn_2 B en utilisant 03 selon la relation suivante qui se déduit de l'équation (8) Hn-2,B - 03 Un+l,B Equation (22) On poursuit alors en considérant la porteuse à la position C à l'étape 168. A l'étape 168, l'estimation de canal Hn_2 C de la porteuse à la position C est égale à l'estimation de canal bruité Hn_2 C qui a été déterminée auparavant. On poursuit alors en considérant la porteuse à la position D à l'étape 170. A l'étape 170, on détermine l'estimation de canal Hn_2 D en utilisant 01 selon la relation suivante qui se déduit de l'équation (6) Hn-2,D = 0T 1 Un-1,DEquation (23) On poursuit alors en considérant la porteuse à la position E à l'étape 172. A l'étape 172, l'estimation de canal Hn_2 E de la porteuse à la position E est égale à l'estimation de canal bruitée Hn_2,E. De nouvelles valeurs des filtres d'interpolation 01, 02 et 03 sont déterminées par le procédé décrit précédemment en relation à la figure 6. Le présent procédé peut être poursuivi pour toutes les porteuses du symbole n-2 correspondant à des pilotes continus, des pilotes répartis ou des porteuses utiles à des positions de pilotes répartis. En effet, on retrouve toujours la succession des étapes 164, 166, 168 et 170 et occasionnellement l'étape 172. Lorsque toutes les porteuses du symbole d'indice n-2 ont été considérées, le procédé continue à l'étape 174. A l'étape 174, le module 16 reçoit le symbole d'indice n+2 et en extrait les pilotes continus et répartis. La détermination des estimations de canal va donc être réalisée pour le symbole d'indice n-1. On poursuit alors en considérant la porteuse à la position A à l'étape 176. A l'étape 176, on détermine l'estimation de canal Hn-1 A en utilisant 01 selon la relation suivante qui se déduit de l'équation (6) : Hn-1,A =A1 Un,A Equation (24) On poursuit alors en considérant la porteuse à la position B à l'étape 178. A l'étape 178, on détermine l'estimation de canal lin-1B B en utilisant 02 selon la relation suivante qui se déduit de l'équation (7) : Hn-1,B = 02 Un+1,B Equation (25) On poursuit alors en considérant la porteuse à la position C à l'étape 180. A l'étape 180, comme la porteuse à la position C du symbole d'indice n+2 correspond à un pilote réparti, on détermine l'estimation de canal bruitée Hn+2 C à partir de l'équation (4). La valeur Hn+2 C est alors mémorisée. En outre, pour le symbole d'indice n-1, on détermine l'estimation de canal Hn-1 C en utilisant 03 selon la relation suivante qui se déduit de l'équation (8) : Hn-1,C - 03 Un+2,C Equation (26) On poursuit alors en considérant la porteuse à la position D à l'étape 182. A l'étape 182, l'estimation de canal Hn_1 D de la porteuse à la position D est égale à l'estimation de canal bruité Hn-1 D qui a été déterminée auparavant. On poursuit alors en considérant la porteuse à la position E à l'étape 184. A l'étape 184, l'estimation de canal fin-1,E de la porteuse à la position E est égale à l'estimation de canal bruitée Hn1 E• De nouvelles valeurs des filtres d'interpolation 01, 02 et 03 sont déterminées par le procédé décrit précédemment en relation à la figure 6. Le présent procédé peut être poursuivi pour toutes les porteuses du symbole n-1 correspondant à des pilotes continus, des pilotes répartis ou des porteuses utiles à des positions de pilotes répartis. En effet, on retrouve toujours la succession des étapes 176, 178, 180 et 182 et occasionnellement l'étape 184. Lorsque toutes les porteuses du symbole d'indice n-1 ont été considérées, le procédé continue à l'étape 186. A l'étape 186, le module 16 reçoit le symbole d'indice n+3 et en extrait les pilotes continus et répartis. La détermination des estimations de canal va donc être réalisée pour le symbole d'indice n. On poursuit alors en considérant la porteuse à la position A à l'étape 188. A l'étape 188, l'estimation de canal Hn A à la porteuse de position A est égale à l'estimation de canal bruité Hn,A qui a été déterminée auparavant. On poursuit alors en considérant la porteuse à la position B à l'étape 190. A l'étape 190, on détermine l'estimation de canal fin ,B en utilisant 01 selon la relation suivante qui se déduit de l'équation (6) : Hn,B '0T Un+l,B Equation (27) On poursuit alors en considérant la porteuse à la position C à l'étape 192. A l'étape 192, on détermine l'estimation de canal Hn C en utilisant 02 selon la relation suivante qui se déduit de l'équation (7) : Hn,C -02Un+2,C Equation (28) On poursuit alors en considérant la porteuse à la position D à l'étape 194. A l'étape 194, comme la porteuse à la position D du symbole d'indice n+3 correspond à un pilote réparti, on détermine l'estimation de canal bruitée Hn+3 D à partir de l'équation (4). La valeur Hn+3 D est alors mémorisée. En outre, pour le symbole d'indice n, on détermine l'estimation de canal Hn,D en utilisant 03 selon la relation suivante qui se déduit de l'équation (8) : Hn,D ù 02 Un+3,D Equation (29) On poursuit alors en considérant la porteuse à la position E à l'étape 196. A l'étape 196, l'estimation de canal Hn,E de la porteuse à la position E est égale à l'estimation de canal bruitée Hn E . De nouvelles valeurs des filtres d'interpolation 01, 02 et 03 sont déterminées par le procédé décrit précédemment en relation à la figure 6. Le présent procédé peut être poursuivi pour toutes les porteuses correspondant à des pilotes continus, à des pilotes répartis ou à des porteuses utiles à des positions de pilotes répartis du symbole n. En effet, on retrouve toujours la succession des étapes 188, 190, 192 et 194 et occasionnellement l'étape 196. Lorsque toutes les porteuses du symbole d'indice n ont été considérées, le procédé reprend à l'étape 150 pour le traitement du symbole d'indice n+1 pour lequel la position relative des pilotes répartis est identique à celle du symbole d'indice n-3. Après la mise en oeuvre du premier exemple de procédé d'interpolation temporelle du canal selon l'invention, on obtient pour un symbole d'indice n les estimations de canal pour toutes les porteuses correspondant à des pilotes continus, des pilotes répartis ou des porteuses utiles à des positions de pilotes répartis, c'est-à-dire une estimation de canal toutes les trois porteuses. Les estimations de canal pour les autres porteuses utiles du symbole d'indice n sont obtenues à l'étape d'interpolation fréquentielle suivante par combinaison linéaire des estimations de canal de porteuses du symbole d'indice n fournies par le procédé d'interpolation temporelle en utilisant un filtre d'interpolation à réponse impulsionnelle finie (ou filtre FIR de l'anglais Finite Impulse Response). Pour un symbole d'indice n et une porteuse à la position k, k peut être écrit sous la forme k = 3r+j, où r est un entier positif, éventuellement nul, et j est un entier compris entre 0 et 2. La porteuse à la position r correspond à un pilote continu, à un pilote réparti ou à une porteuse utile à une position de pilote réparti. De ce fait, l'estimation de canal de la porteuse à la position r a été fournie par le procédé d'interpolation temporelle. On détermine un vecteur Zn,r de la façon suivante : Hn, 3r-3(P-1) Hn, 3r-3(P-2) Zn,r = Hn, 3r-3 Hn, 3r Hn, 3r+3 Equation (30) Hn, 3r+3(P-1) Hn, 3r+3P où P est un entier qui dans le présent exemple de réalisation est pris égal à 6. On appelle Hn k l'estimation de canal de la porteuse du symbole d'indice n à la position k fournie à l'étape d'interpolation fréquentielle pour la différencier d'une estimation de canal fournie par le procédé d'estimation temporelle. L'estimation de canal Hn k est obtenue selon la relation suivante : Hn,k =(TZn,r Equation (31) où (D1, i étant compris entre 0 et 2, est un vecteur à coefficients complexes de dimension 2P qui représente la réponse impulsionnelle du filtre utilisé pour réaliser l'interpolation en fréquence. Il s'agit, par exemple, d'un filtre de Wiener. Les filtres d'interpolation fréquentielle (D1, i étant compris entre 0 et 2, sont déterminés en fonction du type de démodulation COFDM mise en oeuvre (par exemple, du type 2K, 4K, 8K) et de la durée de la réponse impulsionnelle choisie (par exemple, 5 }'s, 20 ps ou 40 }'s) . Pour k inférieur à 3P-3, les estimations de canal utilisées pour la détermination de Znr, pour lesquelles la position en fréquence seraient inférieures à 0, sont remplacées par Hn p. De façon analogue, Pour k supérieur à N-1-3P, les estimations de canal utilisées pour la détermination de Zn,r, pour lesquelles la position en fréquence seraient supérieure à N-1, sont remplacées par Hn N-1 A la fin de l'étape d'interpolation fréquentielle, on dispose d'une estimation de canal pour chaque porteuse du symbole d'indice n. On notera que, pendant l'étape d'inter- polation fréquentielle, les estimations de canal des porteuses utiles à des positions de pilotes répartis, qui avaient été précédemment fournies par le procédé d'interpolation temporelle, peuvent être à nouveau déterminées à l'étape d'interpolation fréquentielle. Selon un second exemple de procédé d'interpolation temporelle, on utilise pour la détermination de l'estimation de canal d'une porteuse utile à une position de pilote réparti non seulement des estimations de canal bruitées futures et passées mais également des estimations de canal de porteuses utiles passées fournies à des étapes d'interpolation fréquentielle précédentes. Ceci permet d'améliorer la précision des estimations de canal fournies par le procédé d'interpolation temporelle et donc les performances de la démodulation COFDM. A titre d'exemple, en plus des cinq estimations de canal bruitées passées et d'une estimation de canal bruitée future utilisées par le premier exemple de procédé d'interpolation temporelle, le second exemple de procédé d'interpolation temporelle utilise les deux estimations de canal de porteuses utiles qui précèdent le symbole considéré. Par rapport au premier exemple de procédé d'interpolation temporelle, les filtres d'interpolation 01, 02 et 03 sont, dans le second exemple de procédé, sont des vecteurs de dimension 8. La figure 8 illustre les étapes mises en oeuvre par le second exemple de procédé d'interpolation temporelle pour la détermination des filtres d'interpolation 0i(p), i étant compris entre 1 et 3. De telles étapes ne sont mises en oeuvre que lorsque la porteuse traitée par le module 16 correspond à un pilote continu. Seules les différences par rapport aux étapes décrites précédemment pour le premier exemple de procédé d'interpolation temporelle et illustrées à la figure 6 vont être décrites. A l'étape 200, par rapport à l'étape 100, les filtres d'interpolation 0i(0), pour i compris entre 1 et 3, sont initialisés de la façon suivante : 0i(0)=(1/8,..., 1/8)T Equation (32) En outre, le vecteur G est, dans le second exemple de procédé d'interpolation temporelle, un vecteur de dimension 8 initialisé de la façon suivante : G(0)=(0,..., 0)T Equation (33) L'étape 202 est identique à l'étape 102. A l'étape 204, par rapport à l'étape 104, au lieu du vecteur Un k, on utilise le vecteur Un ki qui est déterminé selon la relation suivante : Un,k,i = Hn-4M,k Hn-4(M-1),k Hn-8, k Hn-3-i,k Equation (34) Hn-l+i,k H n,k Les étapes 206 à 232 sont identiques aux étapes 106 à 132 du premier exemple de procédé d'interpolation temporelle. A la différence du premier exemple de procédé d'interpolation temporelle, dans le second exemple de procédé d'interpolation temporelle, après l'étape 216, le procédé se poursuit à l'étape 204. En effet, un nouveau vecteur Un ki doit être déterminé lorsque l'indice i est incrémenté. Le second exemple de procédé d'interpolation tempo- relie comporte par ailleurs les mêmes étapes que le premier exemple de procédé d'interpolation temporelle décrit précédemment en relation avec la figure 7. Toutefois, à la place des équations (6), (7) et (8), on utilise respectivement les équations suivantes : Hn-1,k =81 Vn,k Equation (35) 14n-2,k = 02 V'n,k Equation (36) 14n-3,k =e3 V"n,k Equation (37) où Vn,k, V'n,k et V"n,k sont des vecteurs de dimensions 8 définis de la façon suivante : Vn,k = V'n,k = Hn-4,k Hn-3,k Hn-2,k Hn,k Les vecteurs Hn-8,k Hn-5,k Hn-4,k Hn-3,k Hn,k Vn,k , V'n,k V" n,k = chacun deux composantes correspondant fournies à des étapes d'interpolation Les équations (17) à (28) Hn-20,k Hn-16,k Hn-12,k Hn-8,k Hn-20,k Hn-16,k Hn-12,k Hn-20,k Hn-16,k Hn-12,k Hn-8,k Hn-6,k Equations (38) 14n-5,k Hn-4,k Hn,k et V"n,k comportent donc à des estimations de canal fréquentielle précédentes. utilisées pour le premier 15 20 exemple de procédé d'interpolation temporelle selon l'invention sont respectivement donc suivantes : Hn-3,A = e3V"n,A Hn-3,C = 0T 1 Vn-2,C T Hn-3,D = 02 Vin-1,D Hn-2,A = 0T 2 V'n,A Hn-2,B = 03V"n+l,B Hn-2,D = 6T 1 Vn-1,D Hn-1,A = 01 Vn,A Hn-1,B = 02TVin+l,B Hnù1,C = 03V"n+2,C T Hn,B = 01 Vn+l,B Hn,C = 02TVin+2,C Hn,D =03 v n+3,D diverses variantes et modifications qui apparaîtront à l'homme de l'art. En particulier, la présente invention a été décrite remplacées par les équations Bien entendu, la présente Equation (39) Equation (40) Equation (41) Equation (42) Equation (43) Equation (44) Equation (45) Equation (46) Equation (47) Equation (48) Equation (49) Equation (50) invention est susceptible de pour une démodulation COFDM pour laquelle les pilotes répartis sont disposés toutes les 12 porteuses et décalés de 3 porteuses d'un symbole au symbole suivant de sorte que tous les 4 symboles, on retrouve la même disposition des symboles répartis. Il est clair que la disposition des pilotes répartis peut être différente. En particulier, le décalage des pilotes répartis d'un symbole à l'autre peut être différent. Le nombre de filtres d'interpolation, 6i, à utiliser par le procédé d'interpolation temporelle est alors adapté en conséquence. REVENDICATIONS 1. Procédé de démodulation COFDM de symboles successifs d'un signal reçu depuis un canal de transmission, chaque symbole comportant des premières porteuses modulées en phase et/ou en amplitude véhiculant des données qui dépendent du symbole et des secondes porteuses de référence (Pr, Pc) dites pilotes dont les positions en fréquence (A, B, C, D, E) varient au moins en partie d'un symbole au symbole suivant, au moins une porteuse donnée des secondes porteuses correspondant, pour chaque symbole, au même pilote, le procédé comprenant, pour chaque symbole, une étape de détermination d'une première estimation (Hn-1,A, lin-2,A, lin-3,A) de la fonction de transfert du canal pour chaque porteuse d'un ensemble de porteuses parmi les premières porteuses dudit symbole telles que, pour les positions en fréquence des porteuses considérées, des symboles différents dudit symbole comprennent des pilotes, correspondant à une combinaison linéaire de secondes estimations (Hn-24,A, Hn-20,A, Hn-16,A, Hn-12,A, Hn-8,A, Hn-4,A, Hn,A ) de la fonction de transfert du canal déterminées pour des pilotes à la fréquence de ladite porteuse, caractérisé en ce que les coefficients de la combinaison linéaire sont déterminés de façon itérative, une nouvelle valeur d'au moins un des coefficients étant égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme comprenant le produit entre un pas d'itération et un terme d'erreur égal à la différence entre la deuxième estimation déterminée pour ladite porteuse donnée d'un symbole qui dépend du coefficient et la première estimation déterminée pour ladite porteuse donnée dudit symbole, le pas d'itération étant déterminé de façon itérative. 2. Procédé selon la revendication 1, dans lequel chaque seconde estimation (Hn-24,A, Hn-20,A, Hn-16,A, Hn-12,A Hn-8,A , Hn-4,A , Hn,A) déterminée pour une porteuse d'un symbole est égale au rapport entre un coefficient reçu associé à la porteuse et un coefficient théorique prédéterminé associé à la porteuse. 3. Procédé selon la revendication 1, dans lequel la première estimation (Hn-1,A, fin-2,A, fin-3,A) est déterminée à partir d'au moins la seconde estimation (Hn-24,A, Hn-20,A Hn-16,A, Hn-12,A, Hn-8,A, Hn-4,A, Hn,A) déterminée pour un pilote à la même fréquence que ladite porteuse d'un symbole reçu après ledit symbole. 4. Procédé selon la revendication 2, dans lequel la nouvelle valeur du coefficient est égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme égal au produit du pas d'itération, du terme d'erreur et du conjugué de la seconde estimation déterminée pour la porteuse donnée d'un symbole qui dépend du coefficient. 5. Procédé selon la revendication 4, comprenant les étapes suivantes, pour au moins la porteuse donnée de chaque 15 symbole : déterminer un vecteur dont au moins certaines composantes correspondent à des secondes estimations (Hn-24,A Hn-20,A, Hn-16,A, Hn-12,A, Hn-8,A, Hn-4,A, Hn,A) déterminées pour la porteuse donnée de symboles reçus avant ledit 20 symbole ; et déterminer une nouvelle valeur du pas d'itération égale à la somme de la dernière valeur du pas d'itération déterminée et d'un terme égal au produit d'un terme constant et de la partie réelle du produit de la transposée du vecteur, 25 d'une matrice de pondération et du conjugué du terme d'erreur. 6. Démodulateur COFDM destiné à recevoir des symboles successifs d'un signal reçu depuis un canal de transmission, chaque symbole comportant des premières porteuses modulées en phase et/ou en amplitude véhiculant des données qui dépendent du 30 symbole et des secondes porteuses de référence (Pr, Pc) dites pilotes dont les positions en fréquence (A, B, C, D, E) varient au moins en partie d'un symbole au symbole suivant, au moins une porteuse donnée des secondes porteuses correspondant, pour chaque symbole, au même pilote, le démodulateur comprenant un 35 circuit de détermination, pour chaque symbole, d'une première estimation (Hn-1,A, Hn-2,A, Hn-3,A) de la fonction de transfert du canal pour chaque porteuse d'un ensemble de porteuses parmi les premières porteuses dudit symbole telles que, pour les positions en fréquence des porteuses considérées, des symboles différents dudit symbole comprennent des pilotes, correspondant à une combinaison linéaire de secondes estimations (Hn-24,A, Hn-20,A, Hn-16,A, Hn-12,A, Hn-8,A, Hn-4,A, Hn,A de la fonction de transfert du canal déterminées pour des pilotes à la fréquence de ladite porteuse, caractérisé en ce qu'il comprend un circuit de détermination des coefficients de la combinaison linéaire de façon itérative, une nouvelle valeur d'au moins un des coefficients étant égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme comprenant le produit entre un pas d'itération et un terme d'erreur égal à la différence entre la deuxième estimation déterminée pour ladite porteuse donnée d'un symbole qui dépend du coefficient et la première estimation déterminée pour ladite porteuse donnée dudit symbole, le démodulateur comprenant en outre un circuit de détermination du pas d'itération de façon itérative. porteuse d'un symbole comme étant égale au rapport entre un 25 coefficient reçu associé à la porteuse et un coefficient théorique prédéterminé associé à la porteuse. 8. Démodulateur selon la revendication 6, dans lequel le circuit de détermination des coefficients de la combinaison linéaire est adapté à fournir la nouvelle valeur du coefficient 30 comme étant égale à la somme de la dernière valeur déterminée dudit coefficient et d'un terme égal au produit du pas d'itération, du terme d'erreur et du conjugué de la seconde estimation déterminée pour la porteuse donnée d'un symbole qui dépend du coefficient. 7. Démodulateur selon la revendication 6, comprenant circuit adapté à fournir chaque seconde estimation (Hn-24 A Hn-20,A Hn-16,A , Hn-12,A, Hn-8,A, Hn-4,A, Hn,A) pour une un 9. Démodulateur selon la revendication 6, dans lequel le circuit de détermination du pas d'itération comprend : un circuit adapté, pour au moins la porteuse donnée de chaque symbole, à déterminer un vecteur dont au moins certaines composantes correspondent à des secondes estimations (Hn-24 A, Hn-20,A, Hn-16,A, Hn-12,A, Hn-8,A, Hn-4,A, Hn,A) déterminées pour la porteuse donnée de symboles reçus avant ledit symbole ; et un circuit adapté à déterminer une nouvelle valeur du pas d'itération égale à la somme de la dernière valeur du pas d'itération déterminée et d'un terme égal au produit d'un terme constant et de la partie réelle du produit de la transposée du vecteur, d'une matrice de pondération et du conjugué du terme d'erreur. 10. Démodulateur selon la revendication 9, dans lequel le circuit de détermination du pas d'itération comprend un circuit adapté à déterminer une nouvelle valeur de la matrice de pondération égale à la somme de la dernière valeur de la matrice de pondération déterminée, d'un premier terme égal au produit du terme d'erreur et du conjugué du vecteur et d'un second terme égal au produit de l'opposée du pas d'itération, du conjugué du vecteur, de la transposée du vecteur et de la dernière valeur de la matrice de pondération déterminée.