Brevet FR2762167A1 — publié le 19981016

Sondage à pondération spatiale d'une voie de transmission La présente invention concerne un procédé de sondage d'une voie de transmission formée de plusieurs canaux. Il s'agit donc d'une technique connue sous le nom de diversité de réception, technique selon laquelle un récepteur comporte plusieurs antennes chacune associée à un canal de transmission distinct. En d'autres termes, l'invention propose une méthode d'estimation de la réponse impulsionnelle de ces canaux de transmission. En effet, dans un système de transmission, notamment par ondes radio, le récepteur reçoit pour chaque canal de transmission un signal émis par un émetteur. Le signal émis subit des fluctuations d'amplitude et de phase dans le canal de transmission, si bien que le signal reçu par le récepteur ne lui est pas identique. Les fluctuations du signal sont essentiellement dues à ce que l'homme de métier appelle l'interférence intersymbole. Cette interférence peut provenir de la loi de modulation employée pour la transmission et elle est également due à la propagation multi-trajets dans le canal. Il se trouve que le signal reçu est généralement issu d'un grand nombre de réflexions dans le canal, les différents trajets empruntés par le signal émis conduisant ainsi à des retards variés au niveau du récepteur. La réponse impulsionnelle du canal représente l'ensemble de ces fluctuations, auxquelles est soumis le signal émis. Il s'agit donc là de la caractéristique fondamentale représentant les transmissions entre l'émetteur et le récepteur. La réponse impulsionnelle du canal est utilisée notamment par un égaliseur qui a précisément pour fonction de corriger l'interférence intersymbole dans le récepteur. Une méthode classique pour réaliser une estimation de cette réponse impulsionnelle consiste à disposer dans le signal émis une séquence d'apprentissage formée de symboles connus. Cette séquence est choisie en fonction de la loi de modulation et de la dispersion du canal, dispersion devant s'entendre ici comme le retard d'un symbole émis empruntant le trajet le plus long du canal par rapport à ce même symbole empruntant le trajet le plus court. La dispersion est couramment exprimée comme un multiple de la durée qui sépare deux symboles émis successifs, soit un nombre de "durée symbole". A titre d'exemple, on citera deux techniques connues d'estimation de la réponse impulsionnelle d'un canal de transmission. La première technique fait appel à des séquences d'apprentissage particulières dites séquences CAZAC, pour l'expression anglo-saxonne "Constant Amplitude Zero Autocorrelation". De telles séquences sont décrites dans l'article de A. MILEWSKI : "Periodic sequences with optimal properties for channel estimation and fast start-up equalization", IBM Journal of Research and Development, vol.27, N05, Sept.83, pages 426-431. Le système de radiocommunication cellulaire numérique GSM fait appel à des séquences d'apprentissage TS formées de 26 symboles notés aO à a25 prenant la valeur +1 ou -1. Ces séquences possèdent les propriétés suivantes En notant d la dispersion du canal qui vaut 4 dans le cas du GSM, l'estimation de la réponse impulsionnelle prend la forme d'un vecteur X à 5 composants notés xO à x4. La séquence de symboles reçus S correspondant à la séquence d'apprentissage TS est formée elle aussi de 26 symboles notés sO à s2s. On suppose naturellement ici que l'émetteur et le récepteur sont parfaitement synchronisés et dans ce cas l'estimation de la réponse impulsionnelle X est donnée par l'expression suivante pourOSk La technique CAZAC présente l'avantage d'une grande simplicité de mise en oeuvre. Cependant, on remarque que chaque composante de la réponse impulsionnelle est établie à partir de seulement 16 symboles reçus. Etant donné que la séquence d'apprentissage comprend 26 symboles et que la dispersion du canal vaut 4, il y a de l'information dans le signal reçu qui n'est pas prise en compte et cela conduit à une réduction des performances par rapport à l'idéal théorique. La deuxième technique connue fait appel au critère des moindres carrés. Elle est décrite notamment dans les demandes de brevet FR 2 696 604 et EP O 564 849. En matière de rappel, cette technique fait appel à une matrice de mesure A construite à partir de la séquence d'apprentissage TS de longueur n. Cette matrice comprend (n-d) lignes et (d+l) colonnes, d représentant toujours la dispersion du canal. L'élément figurant à la ième ligne et à la jième colonne est le (d+i-j)ième symbole de la séquence d'apprentissage La séquence d'apprentissage est choisie telle que la matrice AtA soit inversible où l'opérateur .t représente la transposition. C'est naturellement le cas pour les séquences CAZAC mais c'est également le cas pour d'autres séquences. Dans la séquence de symboles reçus, on ne prend pas en compte les quatre premiers sg à s3 car ceux-ci dépendent également de symbole inconnus émis avant la séquence d'apprentissage, étant donné que la dispersion du canal vaut 4. Par un abus de langage on définira donc dorénavant le signal reçu comme un vecteur S ayant pour composantes les symboles reçus, s4, sg, s6, ... I 25 Dès lors, l'estimation de la réponse impulsionnelle prend la forme suivante X = (At A)-1 At . S Cette technique des moindres carrés est un peu plus complexe que la précédente mais il faut noter que la matrice (At A)-l At est calculée une seule fois. On remarque ici que chacune des composantes de l'estimation de la réponse impulsionnelle X est établie à partir de 22 symboles reçus et non pas de 16 comme dans le cas de la technique CAZAC. On doit donc s'attendre à une amélioration des performances. Cependant, quelle que soit la technique utilisée, les réponses impulsionnelles de chacun des canaux de la voie de transmission sont considérées comme indépendantes. La présente invention a ainsi pour premier objet un procédé de sondage d'une voie de transmission qui prend en compte le fait que les différentes antennes sont spatialement liées. Selon l'invention, le procédé de sondage d'une voie de transmission s'applique lorsque cette voie est formée de plusieurs canaux et il nécessite une estimation de la réponse impulsionnelle C1, C2, ..., Cn de ces canaux. Le procédé comprend les étapes suivantes - acquisition d'une statistique spatiale de cette voie de transmission, - établissement d'une réponse impulsionnelle corrigée au moins par pondération de ces estimations de réponses impulsionnelles au moyen de cette statistique spatiale et d'une estimation du bruit additif de ces canaux. Avantageusement, cette statistique spatiale correspond à une estimation de la corrélation des canaux de transmission pris deux à deux. Suivant un mode de réalisation préférentiel, cette estimation de la corrélation des canaux de transmission prend la forme d'une matrice de corrélation spatiale ainsi définie que l'élément figurant à la ième ligne et à la jième colonne soit obtenu par lissage du produit de la transposée hermitienne de la réponse impulsionnelle estimée du ième canal et de la réponse impulsionnelle estimée du jième canal. Selon une caractéristique additionnelle, un signal reçu S par un premier canal correspondant à une séquence d'apprentissage émise, en notant A la matrice de mesure associée à cette séquence d'apprentissage, l'estimation du bruit additif N01 de ce canal est obtenue par normalisation de l'énergie du vecteur (S - AC1) On peut prévoir que cette normalisation soit suivie d'une étape de moyennage. De plus, une matrice de bruit N étant formée à partir de l'estimation du bruit additif Nol, N02, . .., N0n des canaux, une matrice de pondération spatiale G' étant ainsi définie à partir de la matrice de corrélation spatiale G et de la matrice de bruit G' = G (G + les réponses impulsionnelles corrigées C'1, C'2, ... , C'n sont obtenues par l'expression suivante Le procédé de sondage d'une voie de transmission est donc mis en oeuvre à partir de l'estimation de la réponse impulsionnelle des différents canaux considérés comme autant de canaux indépendants. Or les erreurs d'estimation sont inévitables. De manière générale, la détermination de la réponse impulsionnelle d'un canal de transmission unique est un problème qui ne peut être résolu de façon exacte en présence de bruit additif. De plus les techniques antérieures font l'hypothèse implicite que cette réponse impulsionnelle peut prendre une forme quelconque. La présente invention a ainsi pour deuxième objet une méthode d'estimation de la réponse impulsionnelle d'un canal de transmission qui présente une meilleure résistance au bruit additif ou, autrement dit, qui mène à une erreur réduite comparée à l'erreur d'estimation des techniques connues. Cette méthode sera avantageusement appliquée pour l'un au moins des canaux formant la voie de transmission avant l'établissement des réponses impulsionnelles corrigées par pondération au moyen de la statistique spatiale et d'une estimation du bruit additif de ces canaux. Selon l'invention, cette méthode nécessite un signal reçu par un canal, ce signal reçu correspondant à une séquence d'apprentissage émise. La méthode comprend les étapes suivantes - acquisition d'une statistique temporelle de ce canal de transmission, - établissement de l'estimation de la réponse impulsionnelle de ce canal, cette estimation étant pondérée par la statistique temporelle du canal au moyen du signal reçu. La statistique temporelle du canal représente une valeur de la réponse impulsionnelle antérieure à l'acquisition du signal reçu. La susdite pondération introduit le fait que la réponse impulsionnelle afférente au signal reçu a une valeur probablement plus proche de cette valeur antérieure qu'une valeur qui en serait très éloignée. Ainsi, statistiquement, l'erreur d'estimation est diminuée. Avantageusement, cette statistique correspond à une estimation de la covariance de cette réponse impulsionnelle. Selon une première variante du procédé, celui-ci comprend les étapes suivantes - lissage de la réponse impulsionnelle et orthonormalisation au moyen d'une matrice de transformation W pour obtenir l'estimation de la covariance qui prend alors la forme d'une matrice L', - recherche des vecteurs propres vi' et valeurs propres ki' associées de cette matrice L', - estimation de la réponse impulsionnelle instantanée du canal à partir du signal reçu et application de cette matrice de transformation W pour former un vecteur X', l'établissement de l'estimation pondérée Xp étant ainsi réalisée où Ng est un nombre réel strictement positif représentant le bruit additif du canal. On peut ici prévoir que le bruit additif soit choisi égal à la plus petite des valeurs propres ki' Lorsqu'un sous-ensemble de ces valeurs propres ki' présente une contribution inférieure à un seuil prédéterminé, on peut également prévoir que chacune de ces valeurs propres soit forcée à la valeur du bruit additif. La complexité s'en trouve réduite d'autant. Selon une seconde variante du procédé, l'estimation de la covariance prenant la forme d'une matrice R, en notant A la matrice de mesure associée à la séquence d'apprentissage, l'établissement de l'estimation pondérée est ainsi réalisé Xp = (AtA + NOR'1-1 At.S où Ng est un réel strictement positif représentant le bruit additif du canal. Par ailleurs, il est possible d'effectuer une étape d'orthonormalisation de la matrice R au moyen d'une matrice de transformation W pour obtenir une nouvelle matrice R', l'estimation pondérée prenant alors la nouvelle forme suivante Xp = W(I + NoR'-l)-l A't.S où la matrice A' est égale au produit de la matrice de transformation W et de la matrice de mesure A. Avantageusement, l'expression (I + NOR'-1)-1 est calculée au moyen du lemme d'inversion matricielle. La présente invention apparaîtra maintenant de manière plus détaillée dans le cadre de la description qui suit où sont proposés des exemples de mise en oeuvre à titre illustratif, ceci en référence aux figures annexées qui représentent - la figure 1, un diagramme identifiant les principales étapes d'un mode de mise en oeuvre du procédé de sondage d'une voie de transmission selon l'invention, - la figure 2, une première variante d'estimation de la réponse impulsionnelle d'un canal selon l'invention, - la figure 3, une deuxième variante d'estimation de la réponse impulsionnelle d'un canal selon l'invention. L'invention sera présentée dans son application au GSM car ce système a le mérite d'être bien connu de l'homme du métier. Il s'agit donc là d'une présentation adoptée dans un souci de clarté mais il ne faut y voir en aucun cas une limitation de l'invention à ce seul système. En référence à la figure 1, le procédé de sondage d'une voie de transmission s'applique lorsque cette voie comporte au moins deux canaux de transmission, n canaux dans le cas général. Les canaux correspondent chacun à une antenne distincte. On considère donc que chaque canal a fait l'objet d'une estimation de sa réponse impulsionnelle respective C1, C2, ..., C n selon l'une quelconque des techniques disponibles. En premier lieu le procédé prévoit l'acquisition d'une statistique spatiale de cette voie de transmission. Par statistique spatiale on entend un ensemble de données reflétant le comportement de cette voie sur une période prédéterminée que l'on conviendra d'appeler période de corrélation. En effet, les différentes antennes étant fixes, les signaux reçus sur ces antennes présentent un certain degré de corrélation. L'invention entend précisément tirer parti de ce fait pour améliorer la qualité d'estimation de la réponse impulsionnelle de l'un des canaux au moins. A titre d'exemple, une telle statistique peut être obtenue au moyen d'une matrice de corrélation spatiale G où l'opérateur .h représente la transposition hermitienne. Cette matrice G carrée de dimension (n, n) peut ainsi être représentée de manière générique par l'élément figurant à la ième ligne et à la jième colonne g jj = V(CiCj ) Cet élément gij est obtenu par lissage du produit h au moyen des réponses impulsionnelles estimées Ci, Cj des ième et jième canaux obtenues durant la période de corrélation. Ce lissage est une estimation de la corrélation de ces deux canaux. On entend ici lissage dans un sens très général, c'est-à-dire toute opération permettant de lisser ou de moyenner le produit CihCj sur la période de corrélation. Un premier exemple de lissage consiste à faire la moyenne de ce produit sur cette période supposée comprendre e estimations successives de chacune des réponses impulsionnelles C1, C2, ..., Cn Un deuxième exemple de lissage consiste à actualiser, à la pième estimation reçue pour chacun des ième et jième canaux, l'expression de lissage Vp1(ChC;) obtenue à la (p - l)ième estimation au moyen d'un coefficient multiplicatif a, ce coefficient compris entre 0 et 1 étant souvent connu sous le nom de facteur d'oubli de lissage Vp(CihCj) = &alpha;CihCj + (1-&alpha;)Vp-1(CihCj) L'initialisation peut se faire par tous moyens, notamment au moyen de la première estimation obtenue ou bien par une moyenne obtenue comme dans le premier exemple pour les premières estimations reçues. Le procédé de sondage propose ensuite de réaliser une estimation du bruit additif Nol, N02, ..., NOn présent dans chacun des canaux au moyen des estimations des réponses impulsionnelles respectives C1, C2, ..., C n de ces canaux. Il sera proposé différentes solutions pour mener à bien cette estimation du bruit et on traitera le cas d'un seul canal, le premier par exemple, étant donné que les mêmes solutions s'appliquent à chacun des canaux. Une solution simple consiste à affecter N01 d'une valeur prédéterminée qui reflète un seuil en dessous duquel il est peu probable que le bruit additif puisse descendre. Cette valeur pourrait être déterminée par une mesure de rapport signal à bruit, ou par les performances du récepteur, ceci à titre d'exemple. Par ailleurs, l'estimation du bruit N01 peut être obtenue au moyen de l'estimation de la réponse impulsionnelle C1 du premier canal et du signal reçu correspondant S. Cette estimation du bruit a donc pu être réalisée avant la mise en oeuvre du procédé de sondage à pondération spatiale selon quelque technique que ce soit. Toutefois, si ce n'était pas le cas, on propose ici une méthode appropriée lorsque l'estimation de la réponse impulsionnelle C1 est acquise selon la technique des moindres carrés. Il s'agit donc d'estimer directement le bruit additif à partir du signal reçu S et de la matrice de mesure A. En effet, en notant N1 le vecteur bruit affectant le signal reçu, il vient que S = AC1 + N1 Compte tenu du fait que les vecteurs S et N1 ont 22 composantes, l'estimation du bruit N01 prend la forme suivante NOl = (22) (S - ACl)h (S - AC1) 22 Naturellement cette estimation du bruit additif N01 peut être moyennée ou lissée. A partir des estimations du bruit additif des différents canaux Nol, N02, ..., NOn, on construit une matrice de bruit N dont l'élément figurant à la ième ligne et à la jième colonne est noté nij. Lorsque i est égal à j, cet élément nij prend la valeur du bruit estimé Noi. Lorsque i et j sont différents, l'élément nij est nul. Par ailleurs, toujours dans un souci de simplification, il est possible de calculer la valeur moyenne Ng des bruits estimés Nol, N02, ..., NOn et de forcer chacun des éléments nij de la diagonale de cette matrice N à cette valeur moyenne Ng. En notant I la matrice identité, il s'ensuit que la matrice de bruit N prend la forme suivante N = Ng I La matrice de corrélation G et la matrice de bruit N permettent de définir une nouvelle matrice, la matrice de pondération spatiale G' G' = G(G + N)-1 Les réponses impulsionnelles corrigées C'1, C'2, C'n au moyen de cette matrice de pondération spatiale G' sont ainsi définies où l'opérateur .t représente toujours la transposition. L'une au moins des réponses impulsionnelles corrigées C'1 est utilisée à la place de l'estimation de la réponse impulsionnelle C1 dans le récepteur. A titre optionnel, l'invention permet également d'améliorer l'estimation de la réponse impulsionnelle C1, C2, ..., Cn de chacun des canaux, ceci avant la mise en oeuvre du procédé de sondage à pondération spatiale d'une voie de transmission, c'est-à-dire avant l'établissement de la ou des réponses impulsionnelles corrigées C'1, C'2, . On prévoit ainsi l'acquisition d'une statistique temporelle pour l'un au moins de ces canaux, le premier par exemple. Par statistique temporelle, on entend un ensemble de données reflétant le comportement de ce canal considéré indépendamment des autres canaux sur une période prédéterminée que l'on conviendra d'appeler période d'analyse. Il s'agit donc d'une représentation du comportement moyen du canal pendant la période d'analyse. Cette statistique peut être établie par quelque moyen que ce soit et en quelque lieu que ce soit. En effet, l'établissement de la statistique peut prendre place dans un équipement quelconque du réseau de radiocommunication. Ce qui importe c'est que le récepteur puisse acquérir cette statistique. A titre d'exemple, une telle statistique temporelle peut être obtenue de la manière suivante. Au moyen du signal reçu S pendant la période d'analyse, on calcule une estimation X de la réponse impulsionnelle selon une méthode connue. Si l'on retient la technique des moindres carrés cette estimation X vaut X = (At A)'1 At.S (1) On rappelle ici que l'émetteur et le récepteur sont supposés synchronisés à mieux que un demi-symbole prés, auquel cas le signal de réception est formé par le vecteur S dont les composantes sont les symboles reçus s4 à s25 synchrones des symboles a4 à a25 de la séquence d'apprentissage TS. Si une telle synchronisation n'était pas acquise, plusieurs solutions sont disponibles pour l'acquérir et on en citera deux à titre d'exemple. La première solution consiste à décaler le signal reçu en avance ou en retard de j périodes symboles de sorte que Sjt = (sq-j, sg-jr s6~j, ... s25-j) On calcule alors l'estimation Xj pour chacun des vecteurs Sj et on retient la valeur jM pour laquelle Xjh Xj est maximal. Cette valeur jM donne la synchronisation escomptée et il suffit de remplacer le vecteur S dans l'équation (1) par le vecteur SjM. La deuxième solution consiste à augmenter artificiellement la dispersion d du canal d'une quantité 2q prédéterminée. On peut alors définir une matrice de mesure modifiée Am comprenant (n-d-2q) lignes et (d+2q+l) colonnes. En reprenant les valeurs 26 et 4 respectivement pour n et Il faut alors réduire le nombre de composantes du signal reçu S de cette même quantité 2q et on peut convenir de retenir le vecteur Sm modifié Smt = (S'4+q, s'5+q, ... s 25-q) On obtient ainsi une estimation modifiée Xm Xm = (Amt Am) Amt.sm Cette estimation modifiée Xm comporte d+2q+l composantes Xmt = (x~q, ... x0, xl, ..., x4, ... X4+q) On recherche alors la valeur jM de j comprise entre -q et +q qui maximise l'expression suivante où l'opérateur .* représente la conjugaison complexe. Cette valeur jM détermine l'estimation X de la réponse impulsionnelle pour une dispersion d=4 xt = (xjM, xjM+l ... XjM+4) La synchronisation s'en déduit immédiatement en appliquant le décalage jM au signal reçu S. On peut dès lors appliquer à nouveau l'équation (1). On construit maintenant une matrice de lissage L par lissage des différentes estimations X obtenues pendant la période d'analyse pour obtenir une estimation de la covariance associée à cette réponse impulsionnelle. On entend ici lissage dans un sens très général, c'est-à-dire toute opération permettant de lisser ou de moyenner la réponse impulsionnelle sur la période d'analyse. On obtient ainsi une représentation statistique du comportement du canal de transmission. Ce lissage peut être obtenu notamment par l'une des deux méthodes proposées ci-dessus, l'expression la plus simple de cette matrice de lissage L étant la suivante où m correspond au nombre de séquences d'apprentissage sur lequel est calculé le lissage. On admet ici que cette matrice de lissage peut être approchée par l'équation suivante L ~ (AtA)-1 Ng + R (2) où Ng représente toujours le bruit présent dans le canal de transmission ou bruit additif et où R est une matrice que l'on a coutume d'appeler statistique à priori du canal car elle représente le comportement du canal abstraction faite du bruit. On admet également que la matrice de mesure A est bien conditionnée, c'est-à-dire que les valeurs propres de la matrice AtA sont très proches les unes des autres. Dans ce cas, il est intéressant de procéder à l'orthonormalisation des vecteurs constitués par les lignes de la matrice de mesure A, mais il ne faut pas voir là une limitation de l'invention. Pour ce faire, on emploie une matrice de transformation W telle que A' = AW et A'tA' = I où I représente la matrice identité. En notant L' la matrice ainsi définie L = WL'Wt, il vient que l'équation (2) peut maintenant s'écrire L' ~ NoI + R' (3) Selon une première variante représentée à la figure 2, on remarque que les vecteurs propres vi' de L' et vi de R' sont identiques tandis que les valeurs propres ki' de L' et ki de R' sont décalées de Ng. Soit en prenant toujours la même valeur de 4 pour la dispersion du canal, pour tout i compris entre 0 et 4 Vil = Vi x i Xi + NO Il apparaît ainsi que la détermination des vecteurs propres et valeurs propres de R' et celle de L' sont identiques sous réserve que Ng soit connu. L'étape d'estimation du bruit sera décrite plus loin pour rendre l'exposé plus clair, bien que cette étape précède celle qui va maintenant être explicitée. Il convient donc maintenant de rechercher les couples valeur-propre/vecteur-propre pour l'une ou l'autre des matrices L' ou R'. Cette étape ne sera pas plus détaillée car bien connue de l'homme du métier. Par ailleurs, il va sans dire que l'on peut annuler les valeurs propres dont la contribution est jugée non significative. Par exemple, si ces valeurs propres sont classées par ordre décroissant, on supprime les dernières qui sont telles que leur somme soit inférieure à un seuil prédéterminé. L'étape suivante consiste à estimer la réponse impulsionnelle instantanée X établie selon l'une quelconque des techniques connues à partir du signal reçu correspondant à la dernière séquence d'apprentissage reçue. En notant X = WX', cette dernière estimation est pondérée par la méthode suivante pour obtenir une pondération temporelle Xp de la réponse impulsionnelle instantanée Pour obtenir la pondération Xp il faut donc estimer le bruit additif Ng. Cette estimation du bruit peut être effectuée notamment selon l'une des méthodes exposées plus haut. Une autre possibilité consiste à considérer que la dernière valeur propre, (la plus faible) de la matrice de lissage L est égale à No #4' = Ng ou X4 = 0. Quelle que soit la méthode retenue, la pondération temporelle Xp de l'estimation de la réponse impulsionnelle instantanée peut alors être réalisée comme mentionné cidessus. Selon une seconde variante représentée à la figure 3, on établit l'estimation pondérée Xp directement comme suit Xp = (AtA + N0R-î)-lAt.S ou bien en reprenant la matrice de transformation W définie ci-dessus Xp = W (I + N0R'-1)-1 A't.S (4) Conformément à l'équation (3) R' = L' - NoI Il convient là aussi d'estimer le bruit additif Ng. Une solution avantageuse pour obtenir la pondération temporelle Xp consiste à adopter la méthode suivante. On divise la matrice R' par N0 : R' B = No Il s'ensuit que I + NOR'-1 = I + B-1 On utilise le lemme d'inversion matricielle pour calculer la matrice de pondération P = (I + B-l)-l. Ainsi, en notant ei les vecteurs canoniques, on procède à l'itération suivante - initialisation P=B - pour i variant de O à d (4 dans le cas présent) P étant connu, il reste à établir la pondération Xp selon l'équation (4). On remarquera que la matrice de pondération P n'est pas nécessairement calculée à l'occasion de chaque nouvelle séquence d'apprentissage émise. Elle peut être calculée à un rythme plus lent car elle varie sensiblement au même rythme que R' et donc plus lentement que le signal reçu S. On remarquera également que l'estimation pondérée est réalisée sans avoir recours à la réponse impulsionnelle instantanée. Elle est produite directement à partir du signal reçu S. Quelle que soit la variante retenue, c' est donc l'estimation pondérée Xp qui sera avantageusement utilisée comme estimation de la réponse impulsionnelle C1 pour la mise en oeuvre du procédé de sondage REVENDICATIONS 1) Procédé de sondage d'une voie de transmission formée de plusieurs canaux, ce procédé nécessitant une estimation de la réponse impulsionnelle (C1, C2, ... ,Cn) de ces canaux, caractérisé en ce qu'il comprend les étapes suivantes - acquisition d'une statistique spatiale de cette voie de transmission, - établissement d'une réponse impulsionnelle corrigée (C'l, C'2, ... ,C'n) au moins par pondération desdites estimations de réponses impulsionnelles (C1, C2, ... ,Cn) au moyen de ladite statistique spatiale et d'une estimation du bruit additif (Nol, N02, ..., NOn) desdits canaux. 2) Procédé selon la revendication 1, caractérisé en ce que ladite statistique spatiale correspond à une estimation de la corrélation desdits canaux de transmission pris deux à deux. 3) Procédé selon la revendication 2, caractérisé en ce que cette estimation de la corrélation des canaux de transmission prend la forme d'une matrice de corrélation spatiale (G) ainsi définie que l'élément (gij) ) figurant à la ième ligne et à la jième colonne soit obtenu par lissage du produit (C.C-) de la transposée hermitienne de la réponse impulsionnelle estimée (Ci) du ième canal et de la réponse impulsionnelle estimée (Cj) du jième canal. 4) Procédé selon la revendication 3 caractérisé en ce que, un signal reçu S par un canal correspondant à une séquence d'apprentissage émise, en notant A la matrice de mesure associée à ladite séquence d'apprentissage, l'estimation du bruit additif (N01) de ce canal est obtenue par normalisation de l'énergie du vecteur (S - AC1) 5) Procédé selon la revendication 4, caractérisé en ce que ladite normalisation est suivie d'une étape de moyennage. 6) Procédé selon l'une quelconque des revendications 3 à 5 caractérisé en ce que, une matrice de bruit (N) étant formée à partir de l'estimation du bruit additif (Nol, N02, ..., NOn) des canaux, une matrice de pondération spatiale (G') étant ainsi définie à partir de ladite matrice de corrélation spatiale (G) et de ladite matrice de bruit G' = G (G + lesdites réponses impulsionnelles corrigées (C'l, C'2, C'n) sont obtenues par l'expression suivante 7) Procédé selon l'une quelconque des revendications 1 à 6 caractérisé en ce que, avant l'établissement de ladite réponse impulsionnelle corrigée (C'1) d'un canal, le signal reçu (S) par ce canal correspondant à une séquence d'apprentissage émise, il comprend les étapes suivantes - acquisition d'une statistique temporelle de ce canal de transmission, - établissement de l'estimation de la réponse impulsionnelle (C1) dudit canal, cette estimation (Xp) étant pondérée par la dite statistique temporelle du canal au moyen dudit signal reçu (S). 8) Procédé selon la revendication 7, caractérisé en ce que ladite statistique temporelle correspond à une estimation de la covariance de ladite réponse impulsionnelle. 9) Procédé selon la revendication 8 caractérise en ce qu'il comprend les étapes suivantes - lissage de ladite réponse impulsionnelle et orthonormalisation au moyen d'une matrice de transformation W pour obtenir ladite estimation de la covariance qui prend alors la forme d'une matrice L', - recherche des vecteurs propres (vi') et valeurs propres (khi') associées de cette matrice L', - estimation de la réponse impulsionnelle instantanée du canal à partir dudit signal reçu (S) et application de cette matrice de transformation W pour former un vecteur X', l'établissement de ladite estimation pondérée (Xp) étant ainsi réalisée où Ng est un nombre réel strictement positif représentant le bruit additif dudit canal. 10) Procédé selon la revendication 9, caractérisé en ce que ledit bruit additif (Ng) est choisi égal à la plus petite desdites valeurs propres (i') 11)Procédé selon l'une quelconque des revendications 9 ou 10, caractérisé en ce qu'un sous-ensemble desdites valeurs propres (khi') présentant une contribution inférieure à un seuil prédéterminé, chacune de ces valeurs propres est forcée à la valeur dudit bruit additif (No). 12) Procédé selon la revendication 8 caractérisé en ce que ladite estimation de la covariance prenant la forme d'une matrice R, en notant A la matrice de mesure associée à ladite séquence d'apprentissage, l'établissement de ladite estimation pondérée (Xp) est ainsi réalisé Xp = (AtA + NOR-1'1 At.S où Ng est un réel strictement positif représentant le bruit additif dudit canal. 13) Procédé selon la revendication 12, caractérisé en ce qu'il comprend une étape d'orthonormalisation de ladite matrice R au moyen d'une matrice de transformation W pour obtenir une nouvelle matrice R', cette estimation pondérée prenant alors la nouvelle forme suivante Xp = W(I + NOR'-1-1 A't.S où la matrice A' est égale au produit de cette matrice de transformation W et de ladite matrice de mesure A. 14) Procédé selon la revendication 13, caractérisé en ce que l'expression (I+N0R'-1)-1 est calculée au moyen du lemme d'inversion matricielle.