FR 2819357 A1 20020712 FR 0100540 A 20010109 La présente invention concerne un procédé et un récepteur permettant l'identification autodidacte des paramètres caractéristiques d'un canal de propagation à partir notamment des statistiques d'ordre 2 du signal reçu par un capteur. Les paramètres identifiés sont par exemple le retard et l'atténuation des trajets du canal. L'invention trouve son application notamment dans le domaine des radiocommunications mobiles ou des liaisons ionosphériques en HF. Canal spéculaire Les radiocommunications mobiles ou ionosphériques HF sont sujettes au phénomène de trajets multiples. En zone urbaine, le signal émis est réfléchi et diffracté sur des obstacles fixes ou mobiles présents dans l'environnement. Dans le cas des transmissions HF, les réflexions se font sur les différentes couches de la ionosphère. Dans le cas de canaux de propagation, pouvant être qualifiés de spéculaire, i.e. la transmission s'effectue selon un nombre limité de trajets discrets ou ponctuels caractérisés par un retard et une atténuation complexe. En considérant un temps d'observation compatible avec la durée de stationnarité du canal, une réception multicapteurs et une propagation spéculaire, le canal multi-capteurs a pour expression: d d c(t) = E a-t t-) k=l o k est l'indice d'un trajet, ak le vecteur dont les composantes sont les atténuations complexes du trajet k pour les différents canaux,:k le retard associé à ce kième trajet et d le nombre de trajets dans un canal. 2 2819357 En outre, si chaque trajet k est incident sur le réseau suivant un cône de diffusion spatiale réduit, I'expression (1) prend la forme suivante d C(t) =I3Dka(ok)(t -rTc) (2) k=l o a(Ok) est le vecteur d'arrivée pour un angle associé au kième trajet et 1k l'atténuation complexe du trajet. Dans le premier cas (Eq. 1), le canal n'est pas paramétré par la direction d'arrivée des trajets, chaque trajet de propagation est paramétré par un temps d'arrivée X et un Ainsi, la calibration n'est plus nécessaire et les algorithmes ne sont plus limités par la dispersion spatiale ou les trajets cohérents. Dans le deuxième cas (Eq. 2), les trajets sont paramétrés par leurs directions d'arrivée, ce qui suppose que la variété de l'antenne est connue et entraîne donc en général la mise en ceuvre d'une calibration pour estimer les valeurs de Ok. Identification autodidacte Dans les systèmes actifs ou entraînés, les paramètres du canal sont calculés au cours d'une phase d'apprentissage durant laquelle I'émetteur transmet une séquence connue du récepteur. Si le canal de propagation fluctue au cours du temps, notamment à cause des mouvements des stations mobiles ou des couches ionosphériques, la séquence doit être émise périodiquement afin de mettre à jour la valeur des paramètres. Si de tels systèmes sont efficaces, l'émission régulière d'une séquence d'apprentissage peut entraîner une réduction importante du débit efficace. Par exemple, dans la norme STANAG 4285 pour les transmissions HF ionosphériques coopératives la moitié des symboles transmis sont des symboles d'apprentissage. 3 2819357 L'art antérieur décrit aussi différentes méthodes et systèmes autodidactes, ou systèmes aveugles, dans lesquels, l'estimation des paramètres est effectuée à partir des statistiques du signal reçu, sans connaissance a priori de la séquence d'apprentissage. Techniques d'ordre deux Les techniques proposées utilisent, par exemple, uniquement les statistiques d'ordre deux (matrice de covariance spatiotemporelle) du signal reçu. Par rapport à des algorithmes aux ordres supérieurs, les algorithmes à l'ordre deux ont de meilleures propriétés de convergence, i.e. pour un nombre de symboles donné, la variance des estimateurs à l'ordre deux est en général inférieure à celle des estimateurs aux ordres supérieurs. En autre, elles présentent moins de problèmes locaux des problèmes d'optimisation que les techniques aux ordres supérieurs. De nombreuses méthodes et algorithmes ont été développés, tels que ceux décrits dans le document ayant pour titre "Multichannel Blind Identification: from subspace to maximum likelihood methods." de L.Tong et S.Perreau; Proceedings of the IEEE, 86(10): 1951- 1967, October 1998. L'un des inconvénients de ces algorithmes est qu'ils nécessitent la connaissance de la longueur du canal de transmission global et que la plupart ne sont pas du tout robustes à une erreur d'estimation de cette longueur. Or, lorsque la transmission est à bande limitée, la longueur du canal n'est définie qu'approximativement et ces algorithmes ne peuvent être utilisés. Etat de l'art en identification paramétrique autodidacte Les principaux travaux d'identification autodidacte des paramètres caractéristiques de la propagation sont des algorithmes découplés, i.e. estimant indépendamment les temps d'arrivées et les vecteurs d'antennes ou les directions d'arrivée. On peut citer les algorithmes de goniométrie pour l'estimation des directions d'arrivée, les algorithmes de séparation de sources (cohérentes) pour 4 2819357 l'estimation des vecteurs d'antennes. Il existe aussi de nombreux travaux concernant l'estimation des temps d'arrivée. Les méthodes d'estimation conjointes conduisent à améliorer la précision et la résolution des estimateurs permettant par exemple, l'estimation des paramètres même lorsque les retards (ou les angles) sont très proches. Les principaux travaux d'estimation conjointe autodidacte des paramètres (0,c) sont présentés, par exemple, dans les références suivantes: "Identification spatio-temporelle de canaux de propagation à trajets multiples", de J.Gouffraud, PhD thesis, Ecole Normale Supérieure de Cachan, 1997, ou * "Improved blind channel identification using a parametric approach" de M.C.Vanderveen et A.Paulraj, IEEE Communications Letters, pages 226-228, August 1998. L'idée est que les critères utilisés pour l'estimation autodidacte de la réponse impulsionnelle peuvent être minimisés directement en fonction des angles et des retards, par des critères de type sous- espace, tels que décrits dans la référence "Subspace methods for the blind identication of multichannel FIR filters" de E. Moulines, P. Duhamel, J-F.Cardoso et S. Mayrargue, IEEE Trans.on signal Processing, 43(2): 516-525, February 1995. Ces algorithmes nécessitent la connaissance a priori du filtre d'émission/réception et une calibration de l'antenne. D'autre part, un exemple de méthode d'estimation conjointe autodidacte des paramètres (a,') est par exemple décrit dans la référence "Methods for blind equalization and resolution of overlapping echoes of unknown shape, de A.Swindlehurst et J.Gunther, IEEE Trans.on Signal Processing, 47(5): 1245-1254, May 1999. Les auteurs travaillent dans le domaine fréquentiel et proposent un algorithme itératif de type IQML pour approximer la solution du 2819357 maximum de vraisemblance ainsi qu'un algorithme d'initialisation explicite basé sur l'algorithme ESPRIT. Ces algorithmes ne nécessitent pas la connaissance a priori du filtre d'émission, mais présentent l'inconvénient de nécessiter une transformation de Fourier avant traitement. L'invention concerne un procédé autodidacte ou partiellement autodidacte permettant de déterminer au moins les paramètres spatio- temporels caractéristiques d'un canal de propagation dans un système comportant au moins un capteur de réception de signaux y(t) caractérisé en ce qu'il comporte au moins une étape o la structure de type spéculaire du canal est prise en compte et une étape de détermination conjointe des paramètres spatio-temporel, tels que les vecteurs d'antenne (a) et/ou les vecteurs temporels (c) à partir des statistiques d'ordre 2 des signaux reçus sur un ou plusieurs capteurs. Selon un mode de mise en oeuvre, il comporte une étape o le signal reçu y(t) est suréchantillonné, par exemple avec une période d'échantillonnage correspondant à T/p, o T est la période symbole. Selon un autre mode de réalisation, le procédé comporte une étape de réception sur au moins deux capteurs et une étape de suréchantillonnage. Le procédé selon l'invention trouve son application, par exemple, dans le cadre des liaisons normalisées à des fins d'égalisation, de localisation et de filtrage spatial ou dans le cas de la surveillance du spectre d'un canal de propagation, à des fins de localisation à partir d'une ou de plusieurs stations HF. L'invention présente notamment les avantages suivants: il n'est pas nécessaire de transmettre une séquence d'apprentissage pouvant être la cause d'une diminution du débit 6 2819357 efficace de transmission, ou de connaître la longueur du canal de transmission globale, elle s'avère robuste à la surestimation de l'ordre du canal de transmission. D'autres avantages et caractéristiques de l'invention apparaîtront à la lecture de la description se référant à la figure unique donnée à titre illustratif et nullement limitatif qui représente le schéma d'un récepteur de radio-communication selon l'invention. Avant de décrire les étapes du procédé selon l'invention quelques rappels, hypothèses et définitions sont données. Hypothèses et modèle Le modèle de signal mono-capteur Le signal émis est de la forme (t)= Eg(t - IT)s(l) (3) LEZ o g(t) est le filtre d'émission/réception et {s(I)} sont les symboles d'information qui ont été émis à la période-symbole T. Le système suppose que la fréquence porteuse et la période symbole ont été estimées de façon préalable. Le signal reçu en bande de base s'écrit y(t) = 3 h(t - lT)s(l) + w(t) (4) lEZ o {s(l)} sont les symboles d'information qui ont été émis à la période- symbole T, h(t) est la réponse impulsionnelle du canal de 7 2819357 transmission, et w(t) est un bruit additif. Le signal reçu échantillonné à la période-symbole suit le modèle discret: y(k) = y(kT) = h (k -)s(l) + w(k) (5) lEZ o h(l) = h(lT) (en supposant que la phase d'échantillonnage soit égale àO). le système multi-voies La figure 1 représente un système de réception multi-voies comportant un nombre q de capteurs de réception et un moyen permettant de suréchantillonner le signal s(k) correspondant au trajet d'indice k. Sans sortir du cadre de l'invention, il est bien entendu que le procédé concerne aussi des récepteurs comportant plusieurs capteurs sans moyen permettant de suréchantillonner ou encore un seul capteur et un moyen de suréchantillonner le signal reçu. Dans l'exemple de réalisation donné ci-après, le dispositif considéré est un système multi-voies combinant une diversité spatiale et une diversité temporelle, i.e. le signal est reçu sur un réseau de q capteurs Ci ayant chacun une réponse h(i)(t) et il est ensuite suréchantillonné par rapport à la période symbole T, au rythme T/p comme représenté sur la figure 1. Le suréchantillonnage est réalisé à un moyen d'un dispositif adapté. La réponse impulsionnelle globale du canal de transmission et le vecteur des observations multi-capteurs sont définis respectivement par les vecteurs colonnes qxl suivants: h(t) = [h(ó)(t),..., h(q)(t)]T (6) y(t) = [y (1)(t) ,..., y(q)(t)] (7) L'indice inférieur indique le numéro du capteur et l'exposant T représente la transposée d'une matrice. Les caractères gras indiquent une dimension spatiale. Les échantillons obtenus par échantillonnage au rythme T/p sont vectorisés en rassemblant les signaux 8 2819357 échantillonnés sur la même phase d'échantillonnage à l'intérieur d'une période-symbole. Plus précisément, à chacune des p différentes phases d'échantillonnage au rythme symbole est associée une voie, créant ainsi sur chaque capteur un canal multi-voies à p sorties. On obtient ainsi un canal multi-voies à p.q sorties. La réponse impulsionnelle globale du canal de transmission et le vecteur des observations multi-voies sont définis par les vecteurs colonnes p.q xl suivants: h(k) = [h(1) (kT)... h(P)(kT)]T (8) y(k) = [y(') (kT)... y(P)(kT)]T (9) L'indice supérieur indique la phase d'échantillonnage, h (i) (kT) = h(kT +(i - 1)T/p), y (i) (kT) = y(kT + (i - 1)T/p). Hypothèses Le procédé selon l'invention utilise certaines hypothèses notamment les suivantes. Source de symboles H 1 - Le processus { s(k)} est par exemple un bruit blanc complexe qui peut aussi être circulaire au second-ordre de variance unité, E [s(k)] = 0, E [s(k)s(j)H] = 8(k - j), E [s(k)s(j)] = 0 (10) Bruit H 2- Le bruit w(k) est un processus stationnaire blanc spatialement et temporellement, Gaussien, circulaire, de variance 2, E [w(k)w(k)H] = (2 Ipq (11) avec Ipq la matrice identité Filtre d'émission/réception H 3 - Le filtre d'émission/réception g(t) a un support fini, g(t) = 0 pour t e [0, LgT[. o LgT correspond à la longueur du filtre et T la période d'échantillonnage. H 4 - Le filtre d'émission/réception g(t) G(f) est approximativement à bande limitée, de bande Bg 9 2819357 G(f) 0 Vf ó [-Bq,g1 (12) et G(f) ne s'annule pas dans une certaine bande j, s ( f) #0 v! [-_l.z (13) Canal de propagation H 5 - Le signal est bande étroite pour le réseau de capteurs, i.e. on suppose que le temps de propagation du signal d'un capteur au suivant est très inférieur à l'inverse de la bande du signal. H 6 - Le canal est constant sur l'intervalle d'observation. Le décalage Doppler et les résidus de porteuse sont supposés négligeables. H 7- L'écart temporel entre le premier et le dernier trajet est fini, et on en connaît une borne supérieure: A'rmax d étant le nombre de trajets et Tk et ak le retard et le vecteur d'antennes de dimension q x I associés au kième trajet la réponse du canal s'exprime de la façon suivante. d h(t) = Saig9(t- Tk) (14) k=l En échantillonnant à T/p, d hT(z) = h(kT/p) z-k= EakgT(rk;z) (15) O gT(r;z) - T pr =,r z 9 (n -r)z (16) n P est la transformée en z de filtre d'émission/réception décalé de r et échantillonné à T/p. Par ailleurs, h(z) correspond à la réponse en z au rythme symbole du système combiné, h(z) h(kT)z-k (17) k=-oo Le signal reçu suit le modèle discret équivalent mono-source multi- sorties: y(k) = [h(z)]s(k) + w(k) (18) il est intéressant de noter que p hT (z) = 11(i) (zp)Z-i (19) i=1 o jh(O(z) = %Looh(i)(k)zk sont les composantes polyphases du canal suréchantillonné. D'après les hypothèses H3 et H7, la réponse impulsionnelle du canal de transmission global multi-capteurs est à longueur finie. Donc, les réponses discrètes h(z) et hT/p(z) sont de degré fini. Nous notons L = deg(_) i.e. au moins une des voies a son degré égal à L, les autres pouvant être de degré inférieur, h(z) = h(O) + h(1)z1 +... + h(L)z'L De même, nous notons LT/p = deg(hT/p). On a pL Afin de simplifier la présentation, nous supposons que LT/p + 1 = p(L+1). Par ailleurs, les algorithmes à l'ordre deux sont basés sur la condition de diversité du canal de transmission global suivante: H 8 Le vecteur polynômial h(z) est irréductible, i.e. ces composantes n'ont pas de racines communes. Quelques notations et outils Soient {b(1) (z),..., b(r)(z)} r polynômes scalaires, b(z) = [b (1) (z),..., b(r)(z)]T est un vecteur polynomial (un vecteur dont les coordonnées sont des polynômes) de dimension r. On Il appelle ordre ou degré du vecteur, deg(b(z)), le degré maximum des polynômes {b()(z),..., b(r)(z)}. Si Lb est le degré de b(z), alors Lb b(z) = b(k)z-k (20) k=O o b(k) = [b(1) (k),..., b(r) (k)]'. Le vecteur associé à b(z), i.e. contenant l'ensemble de ses coefficients est désigné à l'aide d'une flèche, = [b(0)... b(Lb)']T (21) est de dimension r(Lb + 1) x 1. Donc, en particulier h correspond au vecteur de dimension q.p(L + 1) xl associé à h(z) et hT/p(z) i.e. contenant l'ensemble de leurs coefficients, h= L[h(O)T. h(L)TIT (22) [hT(0)T...hT(LT)T]T (23) p P P De plus, soit K un Entier correspondant au nombre d'observations multidimensionnelles, K(b) est la matrice de Sylvester (i.e. Toeplitz par bloc) de dimension r(K + 1) x (K + Lb + 1) associée au polynôme b(z), et définie par b(0)... b(Lb) -b-) =0 b(o0) b(Lb)". :O-.. ". (24h o.. o b(0)... b(Lb) 12 2819357 Enfin, soit B(z) une matrice polynomiale de dimension r x D et de degrée Lb, B(z) = [b(z),..., bD(z)] o b1(z) est un vecteur polynomial de dimension r. Nous notons B = [bl,..., bD] la matrice associée de dimension r(Lb + 1) x D et DK(B) la matrice de Sylvester par bloc de dimension r(K + 1) x (K + Lb + 1)D définie par VDK(B) = [T(bl),..., z(bD) (25) Xmin(A) désigne la plus petite valeur propre de la matrice A et Vpmin (A) le vecteur propre associé à Xmin(A). Procédé selon l'invention Le procédé selon l'invention repose notamment sur la prise en compte de la structure de canal de type spéculaire pour déterminer les paramètres spécifiques d'un canal de propagation ou de transmission d'un signal notamment dans un système de communication tel que décrit à la figure 1. Les paramètres déterminés conjointement sont par exemple le vecteur temporel (X) des vecteurs retards ('k) et le vecteur d'antenne (a) des vecteurs (ak). Une première variante de mise en oeuvre du procédé, ou "procédé a priori", est particulièrement bien adaptée lorsque le filtre d'émission/réception est connu. Une seconde variante de mise en oeuvre du procédé ou "procédé sans a priori" correspond mieux au cas o les caractéristiques du filtre d'émission/réception ne sont pas connus. Procédé avec a priori Le procédé suppose que le filtre d'émission g(t) est connu. 13 2819357 Les différentes étapes du procédé ou algorithme avec à priori selon l'invention utilisent le fait que l'équation (15) s'écrit sous forme "vectorielle" de la façon suivante h = GL,,(-)a (26) h est le vecteur de dimension q.p.(L + 1) contenant l'ensemble des coefficients de la réponse hT/p(z), a est le vecteur de dimension q.d contenant l'ensemble des vecteurs d'antennes a = [a... aT]T et GL,d(Tr) est une matrice de dimension q.p(L + 1) x q.d contenant les versions retardées échantillonnées à T/p du filtre d'émission: g(O - -ri)... g(O - -d) g(0- g(0 - -rd) GLd(T) = Pd) q (27) k g((L + 1)T - -)g((L + 1)T- -rd) q est la matrice identité de dimension q. Nous notons t le vecteur de dimension d contenant l'ensemble des retards = [...] Les différentes étapes du procédé sont par exemple les suivantes 1. Choisir la longueur du canal global de propagation tenant compte de l'émission et de la réception, L > Lg + Armax et la valeur du nombre d'observations K > L avec A'rma la plus grande valeur possible du retard relatif entre deux trajets. 2. Calculer la matrice QL(R) du critère sous-espace par exemple de la manière suivante: - Calculer l'estimée empirique de la matrice de covariance: 14- K R N il et N-K le nombre d'observations indépendantes. - Construire la matrice L en prenant par exemple les vecteurs propres associés aux pq(K + 1) - (K + L + 1) plus petites valeurs propres de R GL peut aussi être obtenue à partir de la décomposition en éléments propres de la matrice des données ^^H Y = [:(K). y(N)].) Puis former la matrice de projection tIL = GLG. - Former la matrice du critère sous-espace QL(R) = DL(:lL)DL(IL)H. 3. Estimer le nombre de trajets d par exemple en mettant en oeuvre la méthode décrite par exemple dans la référence "Detection of signals by information theoric criteria" de M. Wax et T. Kailath, IEEE Trans.on Acoust.Speech and Sig. Proc., 33(2): 387-392, April 1985. 4. Estimer les retards = arg tain JLd (r) avec min(GL d(T)HQL (R)G^A(T)) (31) J,,(0) L min (GLd(T)HGLd(T)) Pour ce faire, appliquer le processus itératif décrit ci-dessous - Poser d = 1. Obtenir une estimation TX de 'rrn en minimisant 4L,(-1). - Poser d = 2. Obtenir une estimation 12 de r02 en minimisant JL,(Ir, - 2] ) en fonction de r2 et en conservant 'r égal à l'estimée 'f obtenue à la première étape. Puis, obtenir une nouvelle estimation de ['roo02] en mettant en oeuvre une minimisation multidimensionnelle de JLd([-1,T2]) en utilisant les valeurs [If,f2] que l'on vient d'obtenir comme valeurs initiales. - Itérer ce processus jusqu'à ce que d = d, pour tous les chemins possibles de manière à obtenir le vecteur des retards. 5. Estimer la valeur du vecteur d'antenne par = VPmin(GLd(T)HQL(R)GL,d()). 6. Former l'estimée de la réponse impulsionnelle: h = GL dU)â Procédé sans a priori o0 La deuxième variante du procédé correspond à un algorithme d'estimation paramétrique totalement ou en majorité autodidacte. L'algorithme associé tient compte, comrnme dans la première variante de réalisation du caractère spéculaire du canal de propagation afin de déterminer ses paramètres, notamment le vecteur d'antenne et le vecteur temporel caractéristiques. Contrairement à la première variante, il ne nécessite pas la connaissance a priori du filtre d'émission. Il suppose simplement que le filtre d'émission est à bande limitée et que sa bande est approximativement connue. Les différentes étapes du procédé ou de l'algorithme sont par exemple les suivantes: 1. Estimer Lg la longueur du filtre d'émission,et choisir L > Lg + A-Tmax etf > L. 2. Appliquer la méthode sous-espace décrite au point 2 du procédé avec a priori: - Estimer la matrice de covariance R par exemple en mettant en 2s oeuvre les étapes décrites précédemment correspondant à l'équation(Eq. 28). - Calculer la matrice de projection sur l'espace bruit fIL à partir des vecteurs propres associés aux pq(K + 1) - (K + L + 1) valeurs propres nulles de la matrice. R. - Former la matrice QL(R) = DL(IIL)DL(rlL)H. 16 2819357 - Obtenir h le vecteur propre associé à la plus petite valeur propre de la matrice QL(R), et former hT (z). P 3. Former le nouveau critère paramétrique: - Choisir la valeur de la longueur du canal R > LT = p(L + 1) - 1 P - et former la matrice TR(hr). P - Calculer la matrice úiL, contenant les vecteurs propres associés P aux q(R + 1) - (R + LT +1) p P plus petites valeurs propre de - TR(h)f. P - Former la matrice DLT(QLT). p Y 4. Estimer le nombre de trajets d, par exemple à l'aide de la méthode précédemment référencée pour l'étape 3 de la première variante de réalisation. 5. Déterminer une valeur minimale de la bande du filtre d'émission:, généralement 3 = 1/T. 6. Choisir v (t), un filtre continu à bande limitée de bande A. Former le filtre v(t), v(t) = vi(t) O 17 2819357 en appliquant le processus itératif proposé pour la méthode sous- espace paramétrique. (La matrice VL,a(r) est définie de la même façon que la matrice GL,dQ(r) mais à partir du filtre v(t).) 8. Estimer les vecteurs d'antennes: â = Vpmin(VL,,(-T)//Lz (fL1)J'LT (LT)/VLd()) (34) P p p' p Applications possibles Le procédé selon l'invention trouve son application, par exemple dans les domaines mentionnés ci-après, donnés à titre illustratif et nullement limitatifs. Réseaux cellullaires Dans le contexte des transmissions cellulaires, la norme de transmission est connue. Le filtre d'émission normalisé est donc connu et il est préférable d'appliquer le procédé avec a priori, utilisant la connaissance a priori du filtre d'émission/réception. Les applications potentielles sont les suivantes: Egalisation La propagation joue le rôle d'une convolution temporelle. A l'exception des communications à très bas-débit (quelques centaines de bits/seconde), cette dispersion temporelle induite par la propagation crée de l'interférence entre symboles (IES). Afin de pouvoir récupérer au mieux les symboles transmis, I'effet dispersif du canal doit être compensé par une étape d'égalisation telle que décrite dans l'ouvrage ayant pour titre "Digital Communications", pour auteurs J.G.Proakis, publié aux éditions McGraw Hill, 1989. Si le filtre d'émission/réception est connu, à partir de la connaissance des paramètres du canal de propagation, il est possible de reconstruire la réponse globale puis de mettre en oeuvre une étape d'égalisation. L'intérêt d'une telle approche est qu'en général en estimant les paramètres de la propagation plutôt que la réponse globale, le nombre de paramètres à estimer est réduit, ce qui permet en général de 18 2819357 réduire la variance des estimateur. On peut ainsi améliorer significativement la qualité de l'identification et les performances de l'égalisation. De plus, certains paramètres caractéristiques du canal de propagation (les retards et les angles d'arrivée) évoluent peu d'une trame à l'autre. La modélisation paramétrique peut permettre de suivre de façon plus efficace les variations des paramètres du canal et d'améliorer encore la qualité de l'identification, par exemple tel que le décrit la référence "Multi-channel MLSE equalizer with parametric FIR channel identification" de J.T. Chen et A. Paulraj, In Proc. VTC, pages 710-713, 1997. Localisation La localisation de l'émetteur, problème classique des applications militaires, est devenue récemment d'un grand intérêt pour des applications civiles en communications mobiles. En effet, la localisation des mobiles dans un réseau cellulaire permet de rendre de nombreux services tels que la localisation des appels d'urgence, la gestion de la taxation, la détection des fraudes, I'aide à la planification cellulaire, I'amélioration de la mise en oeuvre du " handover ", telle que décrite dans la référence " Position location using wireless communications on highways of the future de T.S.Rappaport, J.H.Reed et B.D.Woerner, IEEE Commun. Mag., October 1996. Il existe plusieurs techniques de localisation. La plus classique est la triangularisation. Dans ce cas, la localisation est e.ectuée à l'aide de deux ou trois récepteurs (stations de base) distincts estimant chacun les angles ou/et les temps d'arrivée du signal émis telle que celle décrite dans la référence "Mobile Localization in a GSM Network "de H.EI Nahas, PhD thesis ENST, june 1999. Ce système présente certains inconvénients qui ont poussé au développement d'algorithmes de localisation à partir d'une station unique. Ces techniques de localisation nécessitent toujours la connaissance a priori des caractéristiques de la propagation. Un procédé de 19 2819357 localisation est proposé par exemple dans la référence" Procédés de localisation spatio-temporel de radiomobiles en milieu urbain " de M. Chenu-Toumier, A. Ferréol et J-J. Monot, Technical report, TCC, 1997. Filtrage d'antenne - SDMA La connaissance des caractéristiques spatiales des différents trajets peut être utilisée afin de focaliser l'émission dans la direction principale d'arrivée et ainsi minimiser l'effet des trajets multiples au niveau des récepteurs (uplink beamformer design, downlink transmission design). Surveillance du spectre Dans le contexte de l'analyse radio-électrique du spectre, le récepteur n'a aucune information a priori sur l'émetteur. Le système de réception scrute le spectre, intercepte un signal transmis parun émetteur inconnu, et l'analyse afin d'en extraire certaines informations qui permettront d'écouter l'émission ou de localiser l'émetteur. Dans ce contexte, il est peu probable que l'on connaisse le filtre d'émission, et le procédé sans a priori pour estimer les paramètres du canal de propagation (.; a) est plus adapté. Dans ce cas, le filtre d'émission n'étant pas connu, il n'est pas possible d'égaliser le signal reçu, mais la connaissance des paramètres (.; a) permet par exemple la mise en ouvre des applications données ci-dessous. Localisation La principale application est la localisation. En HF, la localisation à station unique utilise classiquement une étape initiale de prévisions ionosphériques. Connaissant l'angle d'incidence et le profil de la densité électronique (par prévision ionosphérique), une méthode de tracé de rayons utilisant la loi de Descartes permet de reconstituer le trajet de l'onde reçue sur le réseau de capteurs selon l'une des références suivantes: 2819357 * "Comparison of the.xing accuracy of single-station location of long range transmitters, de H.C.H"ring, IEEE Proceedings, 37(3): 173- 176, June 1990. À "loniospheric modelling in support of single station location of long range transmitters de L.F.McNamara, Journal of Atmospheric and terrestrial Physics, 50(9): 781-795, 1988. En outre, la connaissance des temps d'arrivée permet de lever certaines ambiguïtés dues à des multibonds. Autres La connaissance des vecteurs d'antennes permet de faire du filtrage spatial. La connaissance des retards permet de recombiner les trajets de propagation en phase. L'utilisation d'une des deux connaissances ou leur utilisation conjointe peut constituer une première étape à la mise en oeuvre d'un système d'analyse du signal. 21 2819357 REVENDICATIONS 1. Procédé autodidacte ou partiellement autodidacte pour déterminer les paramètres spatio-temporels caractéristiques d'un canal de s propagation dans un système comportant au moins un capteur de réception recevant un signal y(t) caractérisé en ce qu'il comporte au moins une étape o la structure de type spéculaire du canal est prise en compte et une étape de détermination conjointe des paramètres, tels que les vecteurs d'antennes (a) et/ou les vecteurs temporels (t) à partir des o statistiques d'ordre 2 des signaux reçus. 2. Procédé selon la revendication 1 caractérisé en ce qu'il comporte une étape de suréchantillonnage du signal reçu. 3. Procédé selon la revendication 1 caractérisé en ce qu'il comporte une étape o les signaux sont reçus sur au moins deux capteurs et une étape de suréchantillonnage du signal reçu. 4. Procédé selon l'une des revendications 2 et 3 caractérisé en ce que la période d'échantillonnage correspond à T/p avec T la période symbole. 5. Procédé selon l'une des revendications 1 à 4 caractérisé en ce qu'il comporte au moins les étapes suivantes: À choisir la longueur du canal global de propagation L > Lg ±"Lmg et la valeur du nombre d'observation K L. Avec Atmax la plus grande valeur possible du retard relatif entre deux trajets, déterminer la matrice QL(R) du critère sous-espace, * estimer le nombre de trajets d, 22 2819357 À estimer les retards = arg ninJLî,('r) () min(L.^t LRG1() JLd(T) = j (GL J(T)HGîd(T) * estimer la valeur du vecteur d'antenne par = VPmi(G ()qL() GL( former l'estimée de la réponse impulsionnelle: h = GLad(-)â 6. Procédé selon l'une des revendications 1 à 4 caractérisé en ce qu'il comporte au moins les étapes suivantes: Estimer Lg la longueur du filtre d'émission, et choisir LÄLg+Azr et K>L. Appliquer la méthode sous-espace décrite au point 2 du procédé avec a priori - Estimer la matrice de covariance R. - Calculer la matrice de projection sur l'espace bruit -L à partir des vecteurs propres associés aux pq(K+1) - (K + L+1) valeurs propres nulles de la matrice R. - Former la matrice - QL(R) = DL(fL)DL(iL)H. - Obtenir h le vecteur propre associé à la plus petite valeur propre de la matrice QL(R), et former h1 (z) p Former le nouveau critère paramétrique: - Choisir la valeur de la longueur du canal R> LT = p(L + 1) - 1 p et former la matrice TR(hT) P p - Calculer la matrice LT, contenant les vecteurs propres associés aux q(R + 1) - (R + Lr + 1) p 23 2819357 plus petites valeurs propres de TR - Former lamatrice DLt(.Lz) > P s. Estimer le nombre de trajets d, À Déterminer une valeur minimale de la bande du filtre d'émission: s. * Choisir v(t), un filtre continu à bande limitée çle bande z. Former le filtre v(t),{ t=it i1er filtre v(t), { v(t) = v(t) 0 1L,a()= " " - Lx.5,rH La( Estimer les vecteurs d'antennes: a =VPmin (VLad(')VDL (CL,)D, (L)V,(+)) 7. Application du procédé selon l'une des revendications 1 à 6 pour la surveillance du spectre d'un canal de propagation à des fins de localisation à partir d'une ou de plusieurs stations en HF. 8. Application du procédé selon l'une des revendications 1 à 6 dans le cadre de liaisons de communication normalisées à l'égalisation ou la localisation ou encore le filtrage spatial.