FR 2782227 A1 20000211 FR 9810224 A 19980804 Méthode d'estimation radio de la vitesse d'un mobile La présente invention concerne une méthode d'estimation de la vitesse d'une station mobile par rapport à une station distante au moyen d'un signal de transmission radio électrique empruntant un canal de transmission qui relie ces deux stations. Le domaine de l'invention est donc celui des radiocommunications mobiles qui permettent d'établir une liaison radioélectrique entre un terminal portable ou10 portatif et un équipement central qui, la plupart du temps, déssert plusieurs terminaux. Cet équipement central fait souvent partie d'une infrastructure telle qu'un réseau et, généralement fixe, il sera référencé en tant que station distante. Par contre, le terminal peut être amené à se15 déplacer alors qu'il est en cours de fonctionnement, si bien qu'on le dénommera station mobile. Le déplacement relatif de la station mobile par rapport à la station distante est caractérisé par une vitesse qui constitue une information importante à plus d'un20 titre. En premier lieu, il s'agit d'un paramètre de fonctionnement de la station mobile et du véhicule dans lequel elle est disposée. Par ailleurs, on mentionnera le domaine plus particulier des réseaux de radiotéléphonie cellulaire.25 Dans ces réseaux, la couverture d'une station mobile est assurée par une pluralité de stations distantes chacune affectée à un secteur géographique identifié dit cellule. Lorsque la station mobile se trouve à la limite de plusieurs cellules, se pose la question de savoir à quelle station30 distante la connecter. On comprend bien que la vitesse de la station mobile est un élément fondamental dans le choix de la station distante appropriée. En effet, si cette vitesse est sensiblement nulle, il est peu probable que la station mobile quitte la cellule dans laquelle elle est inscrite alors que si cette vitesse est élevée, il est beaucoup plus probable que cette station soit en train de rejoindre une nouvelle cellule. A titre complémentaire, on citera les réseaux microcellulaires dans lesquels chaque emplacement est couvert simultanément par une microcellule et par une cellule parapluie qui se superpose à une pluralité de microcellules adjacentes. On pressent ici que des stations mobiles animées d'une vitesse nulle ou faible seront préférentiellement raccordées à des microcellules tandis que10 les stations mobiles se déplaçant rapidement devraient être connectées à la cellule parapluie pour éviter de trop nombreux et trop fréquents changements de microcellules. Il apparaît ainsi que la vitesse d'une station mobile est un paramètre d'une grande utilité lorsqu'il s'agit de sélectionner la station fixe à laquelle la raccorder, procédure connue sous le vocable anglo-saxon de "handover". La présente invention a ainsi pour objet une méthode d'estimation de la vitesse d'une station mobile reliée à une station distante par un canal de transmission qui véhicule20 un signal de transmission radioélectrique. Selon l'invention, cette méthode comprend la succession des étapes suivantes: * réception et filtrage du signal de transmission empruntant le canal de transmission pour obtenir un signal de mesure reflétant les variations temporelles de ce canal, * détermination d'une fonction de corrélation de ce signal de mesure paramétrée par la vitesse de la station mobile, * recherche de la valeur de la dérivée seconde de cette fonction de corrélation à l'origine, * estimation de la dérivée première du signal de mesure et calcul de la variance de cette estimation, * identification des résultats obtenus lors des deux étapes précédentes. Avantageusement, le filtrage du signal de transmission consiste à supprimer le bruit additif du canal de transmission. Selon une caractéristique additionnelle de l'invention, le signal de mesure est la réponse impulsionnelle dynamique du canal de transmission. Dans ce cas, la méthode comprend de préférence une étape pour produire la réponse statique r du canal de transmission; le canal de transmission présentant une dispersion temporelle déterminée, le signal de transmission s étant une séquence d'apprentissage, une matrice de mesure M étant établie à partir de la séquence d'apprentissage10 compte tenu de la dispersion temporelle, cette étape pour produire la réponse statique r du canal se fait selon l'expression suivante: s = M.r + n o n représente le bruit additif du canal. De plus, on prévoit une étape supplémentaire pour rechercher le vecteur propre r0 associé à la plus grande valeur propre de la covariance de cette réponse statique r; la dérivée première du signal de mesure est égale à un coefficient a0 déterminé de sorte que le signal de20 transmission s soit maintenant défini par l'expression suivante: s = M.r + a0TM.r0 + n o T est une matrice représentant le temps. L'invention apparaîtra maintenant avec plus de détails dans le cadre de la description qui suit d'un exemple de réalisation donné à titre illustratif en référence à l'unique figure annexée qui représente sous forme symbolique les principales étapes de cet exemple de réalisation. En notant X(t) une fonction complexe du temps, on adoptera les notations suivantes: * X'(t): dérivée première de la fonction X(t), * X"(t): dérivée seconde de la fonction X(t), * IX(t)J: module de la fonction X(t), * X*(t): fonction conjuguée complexe de la fonction X(t), * E[X(t) X*(t)]: variance de la fonction X(t). On rappelle que la fonction de corrélation temporelle R(T) de X(t) s'écrit comme suit: R(x) = E[X(t)X*(t - T)] En dérivant cette expression par rapport à t, on obtient: R'(T) = - E[X(t)X*'(t - t)] En effectuant le changement de variable t' = t - t, l'expression précédente devient: R'(T) = E[X(t' + t) X*'(t')] En procédant à une nouvelle dérivation de cette expression par rapport à t, il vient que: R"(t) = - E[X'(t' + t) X*'(t')] Pour la valeur T = 0, l'expression ci-dessus prend donc la valeur suivante: R"(0) = - E[X'(t) X*'(t)], ou encore R"(0) = - E[IX'(t)12] (1) Ainsi, la dérivée seconde R"(0) de la fonction de corrélation à l'origine a même valeur que la variance de la dérivée première de la fonction X(t).20 L'invention propose de réaliser une estimation de la vitesse de la station mobile à partir de ce résultat. On fait ici l'hypothèse que la station mobile reçoit un signal de transmission en provenance de la station distante par un canal de transmission. Il faut cependant25 remarquer que le problème est entièrement symétrique et que la solution serait la même si c'était la station distante qui recevait un signal en provenance de la station mobile. La première opération consiste à sélectionner un modèle de corrélation du signal reçu. Un tel modèle qui a pour paramètre la vitesse de la station mobile par rapport à la station distante peut-être déterminé de manière empirique. Il peut également correspondre à un modèle théorique et, dans le cadre du présent mode de réalisation, on retiendra le modèle dit "modèle de Jakes" qui est bien35 adapté à la situation. Selon ce modèle, pour un signal d'énergie normalisée, en notant J0 la fonction de Bessel d'ordre 0, la fonction de corrélation s'écrit: R(T) = J0(-x) o * v est la vitesse, * X la longueur d'onde du signal. La dérivée première de cette fonction de corrélation s'écrit alors: La dérivée seconde s'écrit, quant à elle: 4 t2v2 2 zv% R" (t) =2 ( O) Or la dérivée seconde de la fonction de Bessel d'ordre 0 à l'origine vaut - 1/2: J"0(0) = - 1/2 Il vient que: 2 i2v2 R"(0) = 2 (2) On obtient ici une relation (2) qui explicite la dépendance entre la dérivée seconde à l'origine de la fonction de corrélation et la vitesse. Dans le cas du modèle de Jakes cette relation est particulièrement simple mais,20 quel que soit le modèle retenu, il existe une relation biunivoque entre ces deux grandeurs. Il convient ensuite de définir un signal de mesure X(t) qui reflète les variations temporelles du canal de transmission.25 Naturellement ce signal de mesure peut-être le signal de transmission lui-même. Dans ce cas, il est préférable de procéder à une opération de filtrage de manière à éliminer le bruit additif du canal de transmission.30 On remarque que si le signal de transmission est obtenu par l'émission d'une impulsion, alors ce signal, après le filtrage mentionné ci-dessus, correspond à la réponse impulsionnelle du canal de transmission. Ainsi, selon une caractéristique avantageuse de l'invention, le signal de mesure est la réponse impulsionnelle du canal. En effet, cette réponse est généralement employée à d'autres fins, notamment pour 5 l'égalisation du signal de transmission, si bien qu'elle est déjà disponible. A titre d'exemple, on indique maintenant un moyen d'estimer cette réponse impulsionnelle dynamique qui dépend du temps par opposition à une réponse impulsionnelle10 statique estimée abstraction faite des variations temporelles du canal de transmission. Il sera maintenant fait référence au système cellulaire numérique de radiocommunications GSM car ce système a le mérite d'être bien connu de l'homme du métier.15 Cette présentation est 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. Ce système fait appel à des séquences d'apprentissage TS formées de 26 symboles notés a0 à a25 prenant la valeur +1 ou -1. Ces symboles en provenance de l'émetteur sont connus du récepteur et l'on englobera donc sous le terme "séquence d'apprentissage" toute séquence de bits qui sont connus à priori de ce récepteur par quelque moyen que ce soit.25 En référence à la figure, le signal de transmission est donc la séquence s de symboles reçus par le récepteur. Cette séquence s, qui correspond à la séquence d'apprentissage TS émise par l'émetteur, est formée elle aussi de 26 symboles.30 En matière de rappel, les techniques d'estimation font appel à une matrice de mesure M construite à partir de la séquence d'apprentissage TS de longueur n. Cette matrice comprend (n-d) lignes et (d+1) colonnes, d représentant la dispersion temporelle du canal. L'élément figurant à la ième35 ligne et à la jième colonne est le (d+i- j)ième symbole de la séquence d'apprentissage: a4 a3 a2 al ao a5 a4 a3 a2 al a6 a5 a4 a3 a2 M = a7.... *................ a25.À......... a21 La séquence d'apprentissage est choisie telle que la matrice MtM soit inversible o l'opérateur.t représente la transposition.5 Classiquement, la dispersion temporelle d du canal valant 4 dans le cas du GSM, l'estimation de la réponse impulsionnelle statique prend la forme d'un vecteur à 5 composantes. Si l'on retient la technique des moindres carrés, ce vecteur r vaut (MtM)-1 Mt.s.10 On construit maintenant une matrice de lissage L par lissage des différentes réponses r obtenues pour les séquences d'apprentissage successivement émises, ceci afin d'obtenir une estimation de la covariance associée à cette réponse statique. On entend ici lissage dans un sens très15 général, c'est-à-dire toute opération permettant de lisser ou de moyenner la réponse statique. Un premier exemple de lissage consiste à effectuer la moyenne de la matrice rrh sur une période supposée comprendre m séquences d'apprentissage, l'opérateur.h20 représentant la transformation hermitienne: L(rrh) = 1 rr m Un second exemple de lissage consiste à actualiser, à la ième séquence d'apprentissage reçue, la matrice de lissage obtenue à la (i-1)ième séquence d'apprentissage au25 moyen d'un coefficient multiplicatif B, ce facteur étant généralement connu sous le nom de facteur d'oubli de lissage et étant compris entre 0 et 1: Li(rrh) = Xririh + (1-X)Li-l(rrh) L'initialisation peut se faire par tous moyens, notamment au moyen de la première estimation r obtenue ou bien par une moyenne obtenue comme ci-dessus pour un faible nombre de séquences d'apprentissage. 5 On recherche alors le vecteur propre r0 associé à la plus grande valeur propre de cette matrice L(rrh). En notant u0 = M.r0, on postule que le signal reçu se présente sous la forme suivante: s = M.r + a0T.u0 + n (3) o n figure le bruit additif dans le canal de transmission et T est une matrice qui figure le temps auquel son reçus les différents symboles, ce temps étant exprimé en durées symbole: -10,5 0 0........ ...... 0..DTD: 0 -9,5 0................ 0 O O - 8,5 0.*.......0. T =......................................... ................,..DTD: 0 9,5 0 ..................... 0 0 10,5..DTD: En fait, T est une matrice diagonale de dimension 22 dont l'élément qui figure à la ième ligne et à la ième colonne représente le temps qui correspond au (d+i)ième symbole de la séquence d'apprentissage, soit a(d+il), l'origine du temps étant arbitrairement fixée entre les20 quinzième et seizième symboles. Il apparaît que la dérivée temporelle de la réponse impulsionnelle dynamique a pour valeur le coefficient O0. En notant I la matrice identité, on introduit un opérateur de transformation A: A = I - M(MtM)-lMt L'équation (3) se traduit alors par l'expression suivante: A.s = c0AT.u0 + A.n En notant u'0 = T.u0, la résolution de cette équation A au sens des moindres carrés donne la valeur estimée a0 du coefficient a0: A usIohA.s a0 = h u 0 A.u'o L'opération suivante consiste à calculer la variance C de la dérivée de la réponse impulsionnelle dynamique: C = E[ X] En reprenant les équations (1) et (2): R"(0) = - E[IX'(t)12] 2 ir2V2 R" (0) =- 2 En reprenant d'autre part l'équation (3): A X(t) = r + tao0r0 Lorsque r0 est normalisé: E[IX'(t)12] = E[a 0 Il vient que: 272v2 C = On obtient ainsi la vitesse de la station mobile: V = - - Il apparaît ainsi que, connaissant un modèle de la fonction de corrélation d'un signal reçu par un canal de transmission, l'on peut procéder à une estimation de la vitesse relative des stations émettant et recevant ce signal. L'invention ne se limite pas au mode de réalisation décrit. En particulier, tout moyen peut-être remplacé par un moyen équivalent. REVENDICATIONS 1) Méthode d'estimation de la vitesse (v) d'une station mobile reliée par un canal de transmission radio à une station distante, cette méthode comprenant la succession des5 étapes suivantes: * réception et filtrage d'un signal de transmission (s) empruntant ledit canal de transmission pour obtenir un signal de mesure reflétant les variations temporelles dudit canal de transmission,10 * détermination d'une fonction de corrélation (R(T)) dudit signal de mesure paramétrée par la vitesse (v) de la station mobile, * recherche de la valeur (R"(0)) de la dérivée seconde de cette fonction de corrélation à l'origine, * estimation de la dérivée première (ao0) dudit signal de mesure et calcul de la variance (C) de cette estimation, * identification des résultats obtenus lors des deux étapes précédentes (R"(0) = C). 2) Méthode selon la revendication 1, caractérisée en ce que le filtrage dudit signal de transmission (s) consiste à supprimer le bruit additif dudit canal de transmission. 3) Méthode selon la revendication 1, caractérisée en ce que ledit signal de mesure est la réponse impulsionnelle dynamique dudit canal de transmission.25 4) Méthode selon la revendication 3 caractérisée en ce que, ledit canal de transmission présentant une dispersion temporelle déterminée, ledit signal de transmission (s) étant une séquence d'apprentissage (TS), une matrice de mesure (M) étant établie à partir de ladite séquence d'apprentissage (TS) compte tenu de ladite dispersion temporelle, elle comprend une étape pour produire la réponse statique (r) dudit canal de transmission selon l'expression: s = M.r + n o n représente le bruit additif du canal. ) Méthode selon la revendication 4 caractérisée en ce que, comprenant une étape pour rechercher le vecteur propre (ro) associé à la plus grande valeur propre de la covariance (L) de ladite réponse statique (r), la dérivée première dudit signal de mesure est égale à un coefficient do déterminé de sorte que ledit signal de transmission (s) soit maintenant défini par l'expression: s = M.r + aoTM.ro + n o T est une matrice représentant le temps.