Brevet FR2821502A1 — publié le 20020830

d'enveloppe spectrale dans la bande non transmise par le codeur. La présente invention concerne un procédé d'estimation d'un canal de propagation à partir de ses statistiques. Elle trouve son application, par exemple dans les réseaux de téléphonie mobile de troisième génération communément appelée UMTS. Elle peut aussi s'appliquer aux récepteurs muiti-capteurs pour des fonctions de métrologie, d'augmentation de capacité ou encore de radio surveillance. Dans un système de transmission, notamment par ondes radio, un lo émetteur émet un signal dans un canal de transmission à destination d'un récepteur. Le signal émis subit des fluctuations d'amplitude et de phase dans le canal de transmission, le signal reçu par le récepteur est constitué de copies temporellement décalées et modifiées du signal émis. Les fluctuations du signai et les décalages génèrent ce que l'Homme du métier appelle I'interférence intersymbole. L'interférence provient notamment de la loi de modulation employée pour la transmission et également de la propagation de multi-trajets dans le canal. 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 à des retards variés au niveau des récepteurs. Ainsi la réponse impuisionnelle du canal représente l'ensemble des fluctuations auxquelles est soumis le signal émis. L'apparition prochaine des réseaux UMTS contraint les équipementiers à adapter les outils de métrologie d'augmentation de capacité ou de radio surveillance à ce nouveau standard. Les performances de ces outils reposent en partie sur une estimation de l'effet engendré par le milieu de propagation (le canal radio) sur les signaux émis par les différentes entités du réseau (station de base, mobile,...). Cette opération est appelée estimation du canal de propagation. Elle permet notamment de contrecarrer I'effet de la propagation pour améliorer la qualité du signal reçu (égalisation) - ou de faire ressortir des informations concernant le milieu de propagation (directions d'arrivée, retards des trajets) afin de mettre en _uvre des traitements spatiaux, tels que la goniométrie. Modélisation Sianal émis La composition du signal UMTS émis par une station de base vers les mobiles (liaison descendante) peut être modélisée par exemple selon le schéma donné à la figure 1. Le signal s(t) est constitué de plusieurs trames, chaque trame lO comportant un nombre de slots donné de durée fixée. Par exemple une trame de 10 ms comporte 15 slots. La figure 2 schématise un modulateur pour la liaison descendante (station de base vers les mobiles). Les symboles Q-PSK (Quadrature Phase Shift keying en anglo saxon), notés bq(O)'..., bq(Ns1) destinés à un utilisateur q sont d'abord multipliés par un facteur de puissance q, avec Ns le nombre de symboies envoyés à cet utilisateur. Chaque symbole est ensuite modulé par une séquence, appelée séquence d'étalement (désignée par le terme anglo- saxon spreading), notée cq de t1, et de taille Nq (facteur d'étalement). Ainsi, on forme à partir du symbole bq(l) la suite bq(l) cq(O),..., uq bq(l) cq(Nq1), 1 étant l'indice du symbole. Les séquences cq sont orthogonales, de sorte que: min(Nq,Np) 2, cq (n)cp (n) = Ny si p = q n=0 = 0 sinon Les suites ainsi formées sont alors multipliées terme à terme par une suite de symboles, (+1ti), appelée code d'embrouillage s (scrambling). Cette suite est périodique, de période trame, et est construite de façon à simuler un signal aléatoire. Le signal résultant à émettre s'écrit donc: yn) = Éuqbq (-g)cq (n[Nq s(77) avec s(n) embrouillage ou scrambling à un instant n donné o Ll désigne l'opérateur > et [] I'opérateur . Ceci peut encore s'écrire: Q y(7l) = >, /l qd q (n) q=o avec: dq (71) = Éq (LN)cq (n(Nq s(n) (2) Ce signal passe à travers un filtre d'émission puis est diffusé par la station de base à l'attention des mobiles. Signal reçu par un mobile Le signal est reçu sur un réseau de capteurs après avoir traversé le canal radio. Après échantillonnage, il peut se mettre sous la forme: ; x(n) = y(n - k)h(k) + b(n) (3) k=0 * o les caractères en gras désignent des vecteurs, et * L correspond à l'étalement du canal exprimé en nombre de chips, les indices n et k à des chips. * x(n) représente: Iorsque l'on échantillonne au rythme chip, x(n) est le vecteur des signaux reçu à l'instant n sur chaque capteur, Iorsque l'on suréchantillonne à un rythme correspondant à un chip/2, x(n) est le vecteur des signaux reçus aux instants n et n+chip/2. * h(k) est le canal multi-capteurs, et * b(n) est un bruit additif regroupant les interférences provenant des autres stations de base et le bruit thermique. Siqnal reçu par une station La figure 3 schématise une possibilité de modéliser la liaison montante (mobiles vers une station de base) détaillée après dans le cadre du procédé selon l'invention. lO Méthodes connues De nombreuses méthodes sont connues pour estimer la réponse d'un canal de propagation. Par exemple, une méthode ciassique procède par corrélation du signal reçu, avec des versions décalées d'une séquence d'apprentissage connue du récepteur, dont seules les étapes nécessaires à la compréhension de l'invention sont rappelées ci-après. Pour l'utilisateur 0, par exemple, I'estimation de la réponse impuisionnelle du canal est faite, classiquement, en corrélant des versions décalées du signal reçu x(n+k) (décalé de l'échantillon k), par la séquence d'apprentissage: h(k)=- s(n) cO(n[NO$bO(-) x(n+:) (4) o P représente le nombre de symboles d'apprentissage, et NoP le nombre de chips d'apprentissage, ou nombre de chips pilote En regroupant les équations (1), (2), (3) et (4), I'expression du canal estimé détaillée correspond à (5): h(k)=l h(l) V p l;0 dO(n+k-l)dO(n) + ZiqlIOh(l)N p Zo dq(+k-l)d I N..P-I * Nt, P '77o b('l + k)do (n) Les propriétés des séquences dq(n) font que (6) : I NuP-I * Nl)P nÉ-o d0(n+k-1)do(n)= I si k =1 v sinon I NoP-1 * I N')P nÉ-O dq(n+k-l)dO())- Ainsi, I'estimée de la réponse impuisionnelle du canal s'exprime en fonction notamment de la réponse du canal h (k) et de trois termes représentatifs des interférences existantes (7): h(k) = h(k) + h(l) 1, do(n+k-I)do(n) I=O,|=k NoP n=0 Q h(l) 1 NaP-1d ( k l)d* 1 NUP-1 +N p 1, b(n+k)dO(n) 10. o Q est le nombre d'utilisateurs du canal de propagation, I, k des indices correspondant à des chips, et n l'indice de l'instant considéré Trois types d'interférences venant perturber l'estimation sont ainsi identifiées: 15. Ie premier terme correspond aux auto-corrélations de la séquence do, Ie deuxième terme à l'interférence des autres utilisateurs, et Ie troisième terme représente la contribution du bruit extérieur et thermique. Cette technique est performante lorsque ia séquence d'apprentissage est longue et lorsque le canal de propagation n'évolue pas ou peu au cours du temps. Dans le cas d'une variation rapide du canal de propagation, il devient nécessaire de l'estimer sur des durées assez courtes. Une autre technique connue sous la désignation > permet d'améliorer la méthode précédente en s'affranchissant des auo-corrélations de la séquence d'apprentissage. Le brevet FR 2 762 164 divulgue un procédé pour estimer la réponse impuisionnelle d'un canal de transmission, qui utilise l'estimation de la matrice de covariance spatio temporelle de la réponse impuisionnelle du canal. Le procédé considère que le bruit d'estimation est blanc, de puissance B et indépendant du canal. Dans ce cas, la matrice est estimée par F=-BI avec la matrice de covariance spatio temporelle du canal estimé. Si B est inconnu, il peut être estimé par la plus petite valeur propre de la matrice estimée, il peut aussi être fixé à une valeur seuil. Une telle hypothèse convient parfaitement à un système, dans lequel chaque utilisateur a un canal de propagation qui lui est dédié, par exemple le système GSM. D'autre part, la technique proposée concerne les récepteurs mono-capteurs. Par contre une telle méthode ne convient plus pour des signaux incluant plusieurs utilisateurs, tel que l'UMTS, o le bruit n'est ni blanc, ni indépendant du canal, notamment en raison de l'interférence multi- utilisateurs. L'objet de la présente invention concerne un procédé pour estimer un canal de propagation à partir de ses statistiques, elles-mémes estimées en exprimant notamment le bruit à partir de la matrice de corrélation empirique des observations. L'invention concerne un procédé pour estimer la réponse impuisionnelle h d'un canal de propagation dans un système comportant au moins un ou plusieurs capteurs. Il est caractérisé en ce qu'il comporte au moins une étape d'estimation des statistiques du bruit additif résultant des interférences et du bruit thermique à partir des statistiques du signal reçu. L'étape d'estimation de la matrice de covariance du bruit est par exemple réalisée à partir de la matrice de covariance empirique des l o observations Rx et du nombre de chips pilotes d'une séquence d'apprentissage transmise avec le signal, la matrice du bruit s'exprimant sous ia forme -R x Le procédé est par exemple utilisé pour estimer la réponse impuisionnelle d'un canal de propagation dans le domaine de l'UMTS pour des liaisons montantes et/ou des liaisons descendantes entre une station de base et un ou plusieurs mobiles. L'invention concerne aussi un dispositif de réception adapté pour estimer la réponse impuisionnelle d'un canal de propagation. Il est caractérisé en ce qu'il comporte un ou plusieurs capteurs de réception d'un signal, un moyen pour échantillonner le signal reçu, un moyen adapté à estimer le bruit à partir des statistiques du canal. Le récepteur selon l'invention peut comporter un moyen adapté à estimer le bruit à partir de la matrice de covariance empirique des observations Rx et du nombre de chips pilote d'une séquence d'apprentissage transmise avec le signal, la matrice du bruit s'exprimant sous la forme -R NoP x Selon une variante de réalisation il comporte par exemple un moyen adapté à déterminer la réponse impuisionnelle du canal sous la forme Hm =[ MH", H", - N pRx]À H", Le récepteur est par exemple utilisé dans le domaine de l'UMTS. lo L'invention présente notamment comme avantage d'estimer un canal de propagation comportant plusieurs utilisateurs, pour lequel le bruit n'est ni blanc, ni indépendant. D'autres avantages et caractéristiques de l'invention apparatront mieux à la lecture de la description qui suit o sont proposés des exemples de mise en _uvre à titre illustratif, ceci en référence aux figures annexées qui représentent: La figure 1, un schéma représentant la composition du signal UMTS, La figure 2 la modélisation du signal pour une liaison descendante, La figure 3, un premier schéma de principe d'une première variante de mise en _uvre de l'invention, La figure 4 schématise un exemple de signal pour la liaison montante (entre une station mobile et une station de base), La figure 5, une deuxième variante de mise en _uvre de l'invention, La figure 6 des résultats comparatifs obtenus en utilisant différents procédés d'estimation. Afin de mieux faire comprendre l'objet de la présente invention, la description qui suit donnée à titre illustratif et nullement limitatif est appliquée pour estimer un canal de propagation dans le domaine de l'UMTS, notamment à partir de ses statistiques. Sous le terme " statistique >>, on désigne la covariance spatio temporelle de la réponse impuisionnelle du canal. Ce système fait appel à des signaux tels que décrits respectivement à la figure 1 pour une liaison descendante, c'est-à-dire le signal émis de la base vers des mobiles, et à la figure 5 pour une liaison montante correspondant à un signal émis d'un mobile vers la station. En appliquant au canal de propagation un modèle statistique lo réaliste, reposant sur l'hypothèse que le canal de propagation est assimilable à un filtre à réponse impulsionnelle finie, il est possible d'améliorer les techniques d'estimation du canal existantes, en tenant compte de ses statistiques. Modélisation du canal radio Le canal radio peut être modélisé comme une somme de trajets, chacun étant caractérisé par un temps d'arrivée, une atténuation et un déphasage. La modélisation considérée dans le cadre de l'invention utilise l'hypothèse que l'atténuation et le déphasage d'un trajet évoluent rapidement par rapport à l'évolution des temps d'arrivée des différents trajets. Sous ces hypothèses, le canal peut être représenté comme un processus aléatoire stationnaire, ce qui permet de dire que chaque réalisation du canal est régie par les mêmes statistiques (une réalisation désigne ici la valeur prise par le canal sur un slot donné i.e. Hm). Ainsi, pour un slot donné d'indice m, o m varie de 1 à M, par exemple, le vecteur des coefficients des canaux échantillonnées correspondant à chaque capteur du récepteur est noté (8): Hm=[hm(O)H hm(1)H... hm(L)H]H O l'exposant H correspond à la transposée hermitienne, L au nombre d'échantillons, m l'indice du slot du signal reçu. Nous supposons que Hm reste constant pendant toute la durée du slot, mais que pour mm', Hm=Hm De même, le canal est estimé à chaque slot m, et son estimée est notée H",. On utilise par exemple la méthode des corrélations successives décalées connue de i'Homme du métier. La réponse impuisionnelle du canal étant estimée pour l'ensemble des slots m du signal reçu, il est possible d'améliorer cette estimée H"" en mettant en _uvre des méthodes connues de l'Homme du métier. L 'exemple de mise en _uvre du procédé selon l'invention, donné 0 à titre illustratif, fait appel à une estimation de Wiener, connue de l'Homme du métier et dont les principes ne seront pas détaillés. La méthode d'estimation de Wiener est par exemple décrite dans le livre intitulé "Fonctions Aléatoires" ayant pour auteurs A. Blanc-Lapierre, et B. Pichinbono, aux Editions Masson 1981. Toute autre méthode d'estimation permettant d'améliorer l'estimée du canal peut étre utilisée sans sortir du cadre de l'invention. Estimation de Wiener du canal L'estimation de Wiener du canal consiste à améliorer l'estimée H,7 à l'aide d'une transformation linéaire de Hn. Cette transformation est choisie de manière à minimiser l'erreur quadratique moyenne entre le vrai canal et l'estimée améliorée. L'estimée améliorée de la réponse impuisionnelle du canal pour un slot m tient compte de la matrice de covariance du canal et de la matrice de covariance du canal estimé et s'écrit par exemple (9): Estimée améliorée de la réponse impulsionnelle du canal pour un slot m r Hn' o r représente la matrice de covariance du canal: = E[H,'HH] avec Hm le vecteur défini dans l'expression (8) qui reste constant ou pratiquement constant sur la durée du slot, et représente la matrice de covariance du canal estimé: = E[H",H,] L' idée de l' invention con siste notam m ent à estime r la covariance du bruit d'estimation afin de la retrancher de la matrice de covariance du canal estimé pour en déduire la matrice de covariance du canal L. A partir lo de ces deux valeurs il est alors possible d'en déduire l'estimée améliorée de la réponse impuisionnelle du canal de propagation. Le bruit d'estimation correspond à la différence entre le canal estimé H" et le canal réel Hn' Ce bruit résulte notamment des interférences et des autres sources de bruit affectant le signal lors de la procédure d'estimation du canal. Pour cela, le procédé utilise l'estimation des statistiques du canal comme il est décrit ci-après. Estimation de A, matrice de covariance du canal estimé Cette matrice peut être estimée de manière connue de l'Homme du métier, de la façon suivante, en supposant avoir accumulé pendant un nombre M de slots des estimées du canal H,n pour lesquels les temps d'arrivée n'ont pas ou peu évolué. A l'aide de ces estimées, m variant de 1 à M et, en utilisant un principe connu de l'Homme du métier (I'estimateur choisi est le meilleur non blaisé, par exemple), il est possible d'estimer la matrice contenant les statistiques du second ordre de 1'estimée du canal. La matrice est alors estimée par (10) l M M, Cet estimateur est consistant, dans la mesure o il converge vers la vraie matrice lorsque le nombre d'observations tend vers l'infini. Estimation de la covariance r du canal Elle s'exprime en tenant compte de la matrice estimée (10) M M, et lo de la matrice estimée de covariance du bruit (11) exprimée sous la forme 1 R - NoP o RX est ia matrice de corrélation empirique des observations, Soit L le no mbre d'échanti l lo ns du canal à esti mer et K le no mb re de capteurs, la matrice Rx est de taille KL*KL. On peut la considérer comme une matrice temporelle de taille L*L dont les éléments sont des matrices Àspatiales de taille K*K. Sur la ième ligne et la jème colonne, (i et j variant de 1 à L) de la matrice tem po rel le. on estim e la matrice spatiale (K* K) Rxjj par la som me s u r n des x(n+i)x(n+j)H ou de manière équivalente des x(n+i)x(n) H ou encore des x(n)x(n+j-i)H. La matrice de covariance du canal est estimée par (12): M, 't' "' N p R x Rx étant un estimateur consistant de la matrice, lorsque le nombre d'observations M tend vers l'infini, cette quantité converge vers E. Estimation de la réponse impuisionnelle du canal de propagation La réponse impuisionnelle du canal Hm est alors estimée à nouveau de la matrice de covariance estimée (12), de l'inverse de la matrice de covariance du canal estimé en appliquant par exemple la méthode d'estimation de Wiener ou tout autre méthode permettant d'améliorer une estimée à l'aide de ses statistiques. On conna^t par exemple une première valeur d'estimée H,' pour le canal de propagation Hm, par exemple celle qui a servi à estimer la covariance du canal estimé et on applique l'estimateur de Wiener, expression (9), pour en déduire la nouvelle valeur estimée Hn' du canal de propagation Hm (13). M Hm =[ M,HmHmH N pRX] Hm En résumé, les différentes étapes du procédé selon l'invention sont exposées à la figure 3 sous la forme d'un algoritUme comportant les étapes suivantes: Réception et échantillonnage du signal multi-capteurs, (a) À Stockage des vecteurs d'observation successifs par exemple dans une mémoire, (b) Estimation de la matrice Rx, ( c) Estimation du canal par corrélation du signal avec la séquence d'apprentissage, (d) À Stockage des estimées successives Hm du canal H"', (e) Estimation de la matrice A, (f) Estimation de la matrice, à partir des estimées ( c) et (f), et en utilisant l'expression de l'estimée du bruit sous la forme N p Rx (g) l o. Amélioration des estimées du canal en appliquant par exemple la méthode de Wiener aux estimoes (e), (f) et (g). Le calcul des estimées est réalisé par exemple au moyen d'un microprocesseur programmé en conséquence, ce dernier pouvant aussi être adapté pour exécuter les étapes du procédé selon l'invention. Le microprocesseur peut être en liaison avec un dispositif de traitement utilisant les valeurs estimées du canal pour réaliser de la démodulation ou de la gonlométrie telles que rappelées très brièvement ci après. Le procédé décrit précédemment dans le cadre des liaisons descendantes s'applique aussi sans sortir du cadre de l'invention pour des liaisons montantes, i.e des liaisons d'une station mobile vers une station de base. La figure 4 représente ie schéma du démodulateur dans le cas de la liaison montante. Le mobile émet un unique canal de contrôle, contenant une séquence d'apprentissage pour l'estimation du canal et des informations de signalisation, ainsi qu'un ou plusieurs canaux de données. Les canaux de donnces sont tous émis à la même puissance soit sur la voie > , soit sur la voie . Les symboles transmis sur chaque canal sont des symboles B-PSK (+1) (Binary Phase Shift Keying). Il est possible d'adapter le modèle choisi pour la liaison descendante à la liaison montante en considérant le canal de contrôle comme l'utilisateur auquel on s'intéresse, et les canaux de donnses comme lo les autres utilisateurs de la cellule. Les autres mobiles émettant vers la même station de base sont traités de la même façon que les autres cellules pour la liaison descendante: ils sont inclus dans le bruit global. Les étapes du procédé décrites dans le cadre de la liaison descendante s'appliquent avec ces hypothèses. La figure 5 représente un exemple de mise en _uvre d'étapes permettant une rébuction de rang des matrices utiilsées lors du procédé. Par exemple les étapes se déroulent de la manière suivante: Réception et échantillonnage du signal, (a) Estimation du canal de propagation, (d,e) 20. Estimation des paramètres du canal (temps d'arrivée, puissance des trajets,...) À Construction de la matrice de réduction, Réduction des observations et des estimées du canal, en exécutant les étapes précédemment décrites ( c) (f) (g) dans le cas normal en utilisant la matrice de réduction, Estimation de Wiener du canal réduit Il est possible d'améliorer à la fois l'efficacité et la complexité de ce procédé en recourant à la rébuction de rang de l'espace-canal. Cette tech n iq ue consiste à paramétre r l'es pace-can al à l 'aide d ' une base orthonormale constituée d'un petit nombre de vecteurs. En effet, généralement, dans le contexte d'une transmission par canal radio, en raison de la spécularité du canal (nombre de trajets dans le canal restreint), la matrice r est dite >, ce qui signifie que les vecteurs Hm se déduisent d'une combinaison linéaire d'un petit nombre de vecteurs u,,..., up formant une base de l'espace-canal. Si ces vecteurs étaient connus, il suffirait alors d'estimer la combinaison linéaire gm,, gm P permettant de reconstruire Hm par: Hni = I gn,pup p=, Regroupons gm,,..., gm,p dans un vecteur gm et u,,..., up dans une matrice U. On a alors: Hn, = Ugm Le vecteur gm est de taille (beaucoup) plus petite que Hm' il est donc plus simple de mettre en _uvre l'estimateur de Wiener de gm que de Hm. On peut en effet estimer gm par: gm =U Hm Dès lors, on peut estimer les matrices de covariance et de gm et de UHXm, o Xm représente le vecteur des observations [xm(0)H, xm(L)H]H reçues au slot m. On peut alors estimer la matrice Q de covariance de 9m 1 ^ Q=A NP et g m Q- g nl Dont on débuit une nouvelle estimée du canal: Hm = Ugm En pratique, les vecteurs u,,..., up ne sont pas connus, mais on peut les estimer de plusieurs fa,cons: Vecteurs propres associés aux plus grandes valeurs propres de la matrice estimée. Cette technique doit être mise en _uvre avec une méthode de détermination du nombre de vecteurs propres à retenir. Vecteurs canoniques associés aux plus grandes valeurs de la diagonale de E. Cette technique revient à ne conserver que les coefficients du canal dont la puissance dépasse un seuil fixé au préalable. Il n'est pas nécessaire d'estimer entièrement mais uniquement sa diagonale, ce qui o est beaucoup plus simple. Vecteurs associés au-filtre de mise en forme décalés des retards du canal. Cette technique ne peut être mise en _uvre qu'à condition d'avoir préalablement déterminé les retards des différents trajets du canal. Une, étape de no rm al isation est à p révoi r. L'étape de réduction de rang peut être mise en _uvre avant les différentes étapes du procédé selon l'invention ou encore après en fonction des cas d'application. Le procédé selon l'invention est utilisé par exemple pour démoduler des signaux ou encore pour effectuer de la goniométrie. Démodulation La figure 6 représente dans un diagramme taux d'erreur en fonction du rapport signal à bruit, les performances de différents estimateurs. La courbe (1) correspond à un estimateur classique, obtenu par la méthode de corrélation, la courbe (Il) une courbe obtenue en mettant en _uvre les étapes du procédé selon 1'invention et la courbe (111) les résultats obtenus en utilisant en plus une étape de rébuction de rang de la matrice. De nombreuses techniques de démodulation reposent sur la connaissance du canal de propagation. Le récepteur RAKE, par exemple, estime le symbole émis à partir de la phase de la quantité: I N"-' 1. bo(I)= >,is(lNO +n) cO(n)h(k) X(iNo +n+k) No,=0 l:=0 Lorsqu'on ne conna^t pas le canal de propagation, ce qui est toujours le cas en pratique, on recourt à une estimée de ce canal. Sur la figure 4, les performances de notre estimateur (et de sa version à rang réduit) sont comparées à celles de l'estimateur classique pour le récepteur RAKE. lo Le nouvel estimateur, et plus encore l'estimateur à rang rébuit, entranent un gain important en terme de performances. GonIométrie coopérative L'invention s'applique particulièrement bien pour améliorer les mesures de goniométrie. En effet, en goniométrie, on détermine les angles d'arrivée de différents trajets de propagation en estimant le vecteur h aux différents instants d'arrivée. L'estimation coopérative des temps d'arrivée et des canaux associés se fait au moyen d'une séquence d'apprentissage. Il est donc possible de l'améliorer en mettant en oeuvre le procédé selon I'invention. La précision de l'estimation des angles d'arrivée est directement liée à la qualité de l'estimation de h. Lorsque le canal estimé est bruité, I'angle mesuré se trouve dans un cone centré sur l'angle réel. Améliorer l'estimation de h permet de réduire l'écart entre l'angle mesuré et l'angle réel, ce qui améliore la résolution du traitement. REVENDiCATIONS 1 - Procédé pour estimer la réponse impulsionnelle h d'un canal de propagation dans un système comportant au moins un ou plusieurs capteurs caractérisé en ce qu'il comporte au moins une étape d'estimation des statistiques du bruit additif résultant des interférences et du bruit thermique à partir des statistiques du signal reçu. 2 - Procédé selon la revendication 1 caractérisé en ce qu'il comporte une étape d'estimation de la matrice de covariance du bruit à partir de la matrice de covariance empirique des observations Rx et du nombre de chips pilotes d'une séquence d'apprentissage transmise avec le signal, la matrice du bruit s'exprimant sous ia forme 1 R NoP x 3 - Procédé selon i'une des revendications 1 et 2 caractérisé en ce qu'il comporte une étape d'estimation de la matrice de covariance du canal estimée sous la forme = M >; H m H mH - N p R x 4- Procédé selon l'une des revendications 1 à 3 caractérisé en ce que l'estimation de la réponse impuisionnelle du canal de propagation utilise une méthode de Wiener et en ce que la réponse impuisionnelle du canal est égale à Hm =[ -HmHH' R]-iH M m=i NoP X m - Procédé selon l'une des revendications 1 à 4 caractérisé en ce qu'il comporte une étape o l'espace-canal est paramétré à l'aide d'une base orthonormale constituée d'un nombre de vecteurs donné ul u2,...up. 6 - Procédé selon la revendication 5 caractérisé en ce que les vecteurs u', u2,..., up correspondent aux vecteurs propres associés aux plus grandes valeurs propres de la matrice estimée F. 7 - Procédé selon la revendication 5 caractérisé en ce que les vecteurs u,, u2,..., up sont des vecteurs canoniques associés aux positions des plus grandes valeurs de la diagonale de la matrice E. lo 8 - Procédé selon la revendication 5 caractérisé en ce que les vecteurs u, u2,..., up sont formés à partir de versions échantillonnées du filtre de mise en forme, décalées des retards de propagation du canal, estimés au préalable, et normalisées. 9 - Procédé selon l'une des revendications 1 à 8 pour estimer la réponse impuisionnelle d'un canal de propagation dans le domaine de l'UMTS pour des liaisons montantes et/ou des iiaisons descendantes entre une station de base et un ou plusieurs mobiles. 10 - Dispositif d'émission et/ou de réception adapté pour estimer la. réponse impuisionnelle d'un canal de propagation caractérisé en ce qu'il comporte un ou plusieurs capteurs de réception d'un signal, un moyen pour échantillonner le signal reçu, un moyen adapté à estimer le bruit à partir des statistiques du canal. 11 - Récepteur selon la revendication 10 caractérisé en ce qu'il comporte un moyen adapté à estimer le bruit à partir de la matrice de covariance empirique des obseNations Rx et du nombre de chips pilote d'une séquence d'apprentissage transmise avec le signal, la matrice du bruit s'exprimant sous la forme l R. NoP 12 Récepteur selon la revendication 10 caractérisé en ce qu'il comporte un moyen adapté à déterminer la réponse impulsionnelle du canal sous la forme Hm =[ M>,H",H"I- N pRx]A Hn' 13 - Récepteur selon rune des revendications 11 et 12 applique dans le