FR 2882156 A1 20060818 FR 0501523 A 20050215 PROCEDE DE LOCALISATION D'UN OBJET MOBILE A PARTIR DE MESURES DE DISTANCE ET DE VITESSE BISTATIQUES DE CET OBJET A L'AIDE D'UN SYSTEME MULTISTATIQUE La présente invention se rapporte à un procédé de localisation d'un objet mobile à partir de mesures de distance et de vitesse bistatiques de cet objet, à l'aide d'un système multistatique. La technique de localisation d'un objet mobile à l'aide de systèmes multistatiques fait classiquement appel à l'exploitation conjointe de mesures élémentaires bistatiques de la distance entre chacun des émetteurs et cet objet et entre cet objet et respectivement chacun des récepteurs. Cette exploitation consiste à rechercher l'intersection des lieux d'ambiguïté correspondant à chaque mesure élémentaire, en l'occurrence des hyperboles ou des ellipses, selon la référence de temps utilisée, pour localiser de façon unique dans l'espace la position de cet objet. Une telle technique est utilisée par exemple dans les systèmes de localisation des transpondeurs secondaires des aéronefs et utilisés en contrôle de trafic aérien. Cette technique connue présente les limitations suivantes: - la géométrie quasi plane de l'implantation des récepteurs (si ce sont des récepteurs au sol) ne permet qu'une localisation dans un plan horizontal, sans 20 fournir l'altitude de l'objet, si les signaux reçus de différentes cibles présentes simultanément dans un même espace n'ont pas un support temporel disjoint, l'association des mesures élémentaires issues d'une même cible devient quasiment impossible, car la combinatoire des cas possibles devenant beaucoup trop importante pour une exploitation pratique. La présente invention a pour objet un procédé de localisation d'un objet mobile, en particulier d'un objet volant, à partir de mesures de distance et de vitesse bistatiques de cet objet à l'aide d' un système multistatique. Par localisation , on entend selon la présente invention la position 3D de cet objet et la détermination de son vrai vecteur vitesse en coordonnées 3D également. Ce procédé doit permettre de localiser cet objet le plus précisément possible dans l'espace, c'est-àdire en mesurant ses coordonnées dans un repère à trois dimensions, et ce, même lorsque cet objet se situe dans un environnement pouvant comporter beaucoup d'autres objets mobiles. Le procédé conforme à l'invention est un procédé de localisation d'un objet mobile à partir de mesures de distance et de vitesse bistatiques de cet objet à l'aide d'un système multistatique, et il est caractérisé en ce que l'on associe au moins deux couples de distances bistatiques tels que pour chaque couple les ellipsoïdes qui les déterminent aient un foyer commun, que l'on résout pour chaque couple l'intersection de ces deux ellipses à différentes altitudes, que l'on optimise par rapport à la variable d'altitude z les différentes intersections, qu'ayant trouvé la meilleure valeur de z, on calcule les coordonnées cartésiennes en trois dimensions de cet objet, et que, connaissant les coordonnées de l'objet et au moins trois vitesses bistatiques, on déduit la direction et le module du vrai vecteur vitesse 3D de l'objet mobile. Selon une caractéristique avantageuse de l'invention, on prédit la position de l'objet mobile à l'étape suivante du calcul, et donc son altitude z prédite, et on affine l'optimisation du calcul de z autour du z prédit en ne gardant dans tous lesdits couples de distances bistatiques que les n plus petites, avec n=3 de préférence. La présente invention sera mieux comprise à la lecture de la description détaillée d'un mode de réalisation, pris à titre d'exemple non limitatif et illustré par le dessin annexé, sur lequel: - la figure 1 est un schéma simplifié de trois ellipses résultant d'une coupe dans un plan horizontal des ellipsoïdes relatifs aux distances bistatiques d'un objet mobile, telles qu'utilisées par le procédé de l'invention, et - la figure 2 est un schéma simplifié explicitant la façon de déterminer la direction et le module d'un objet volant à partir de la connaissance des ellipsoïdes de la figure 1. La présente invention est décrite ci-dessous en référence à la localisation dans l'espace d'un aéronef (appelé ici: avion), mais il est bien entendu qu'elle n'est pas 30 limitée à cette seule application, et qu'elle peut être mise en oeuvre pour la localisation de tout objet mobile, en particulier un objet volant, mais non exclusivement. Cette localisation est mise en oeuvre à l'aide d'un système multistatique à plusieurs émetteurs et récepteurs. L'invention consiste à exploiter conjointement les mesures de distances bistatiques et de vitesses Doppler bistatiques des avions à localiser, ces mesures étant effectuées sur un ensemble de plusieurs (au moins trois) couples émetteur-récepteur . Ces couples sont déterminés de la façon suivante. On a représenté en figure 1 l'intersection d'un plan horizontal (qui est ici le plan du dessin) et de trois ellipsoïdes (El, E2, E3) particuliers relatifs à un même avion à localiser, dans un système multistatique (non représenté en détailsur la figure). Les traces, dans le plan du dessin, des ellipsoïdes, sont trois ellipses ell, e12 et e13. Ces ellipsoïdes sont tels que leurs foyers sont, de façon connue en soi, les émetteurs et les récepteurs du système bistatique de localisation. On obtient ainsi un système de trois équations (définies par les équations respectives des ellipsoïdes) à trois inconnues, qui a, théoriquement, une solution permettant de calculer la position 3D de l'avion en question, mais il n'existe pas de méthode directe pour le résoudre (il s'agit en fait de la résolution d'une équation du 8è degré). Selon l'invention, les ellipsoïdes sont tels que, pris deux à deux, ils ont un foyer commun. Dans l'exemple de la figure 1, on a pris trois tels couples d'ellipsoïdes (El, E2 et E3), mais on pourrait en prendre davantage. Dans cet exemple, le système multistatique comprend deux émetteurs Em1 et Em2 et deux récepteurs R1 et R2. Les ellipsoïdes E1 et E2 ont comme foyer commun R1, tandis que El et E3 ont comme foyer commun Eml. Grâce à cette caractéristique de l'invention, la résolution de ce système d'équations peut se faire plus facilement qu'avec ladite méthode directe. En outre, la connaissance des vecteurs vitesse bistatiques pour chacun des trois systèmes bistatiques permet de calculer le vecteur vitesse 3D de l'avion. Pour la mise en oeuvre de l'invention, on commence par rechercher les intersections des couples d'ellipsoïdes précités, à foyer commun pour chacun de ces couples, et ce, dans un plan horizontal. On procède alors de la façon suivante: on appelle (x, y, z) les coordonnées cartésiennes en 3D de l'avion à localiser, (xei, yei, zei) celles de l'émetteur Em1 et (xri, yri, zri) celles du récepteur R1. Soient: - da: la distance entre Eml et l'avion + la distance entre l'avion et R1, - db: la distance entre R1 et l'avion + la distance entre l'avion et Em2, - d: la distance entre Eml et l'avion + la distance entre l'avion et R2. On obtient, de façon classique, l'équation des trois ellipsoïdes précités, respectivement El, E2 et E3: da =xel)2+(Y Yei)2+(z zei)2 + \/(x xri)2+ (Y Yrl)2+(z zri)2 db =.J(x xri)2 +(y Yrl)2 +(z zri)2 + . J(x xe2)2 +(Y ye2)2 +(z ze2)2 d,, = . /(x xei)2 +(y )2ei)2 +(z zel)2 + (x xr2)2 +(y Yr2) 2 +(z zr2)2 Dans une phase dite d'initialisation, c'est-à-dire pour la première mesure effective des distances bistatiques de l'avion, et pour une altitude z donnée de cet avion, on obtient les coordonnées des points d'intersection des couples d'ellipsoïdes de la façon suivante, pour trois plans horizontaux coupant ces ellipsoïdes et ayant les altitudes respectives zi, z2, z3. Dans le présent exemple, on se contente de trois valeurs de z, qui sont suffisantes pour la précision recherchée en phase d'initialisation. De préférence, dans cette phase d'initialisation, ces trois valeurs de z sont prises très différentes les unes des autres, par exemple 0, 500m et 1000m d'altitude. On obtient alors les valeurs suivantes de ces points d'intersection, pour les différents couples d'ellipsoïdes précités: z=zl E1n E2 Eln E3 E2n E3 (xl, yl) (x'l, y'l) (x"l, ylll) z=z2 (x2, y2) (x'2, y'2) (x"2, y"2) z=z3 (x3, y3) (x'3, y'3) (x"3, y"3) Chaque coordonnée (xi, yi), (x'i, y'i), (x"i, y"i) est donc fonction de z et on peut représenter toutes ces coordonnées par un polynôme du second degré (3 valeurs disponibles par coordonnée) : x=aXZ2+ bXz+cX x'=a'Xz2+b'Xz+c'X x"=a"Xz2+b"Xz+c" X y=ayz2+byz+cy y'=a'yz2+b'yz+c'y y"=a"yz2+b"yz+c"y On calcule ensuite les équations des différentes distances Di, en fonction de z, entre les points d'intersection M (x, y), M'(x', y') et M"(x", y"), puis on dérive, selon z, les équations des différentes distances Di (z) et on les résout pour obtenir les minima de ces distances Di (z). Ces minima correspondent à des distances minimales entre les points M, M' et M" pris deux à deux. Selon une caractéristique avantageuse de l'invention, on retient parmi ces distances minimales la plus petite. Soit Zavion cette plus petite distance; elle correspond à un couple de points, par exemple M(x, y) et M'(x', y'), qui résultent du calcul des polynômes du second degré mentionnés ci-dessus, dans lesquels on prend z= Zavion. On en déduit la position 3D de l'avion: - selon l'axe des x, c'est le barycentre des deux valeurs calculées de x, par exemple (x+x')/2, selon l'axe des y, c'est le barycentre des deux valeurs calculées de y, par exemple (y+y')/2, - selon l'axe des z, c'est bien entendu Zavion. Enfin, on calcule le vrai vecteur vitesse 3D de l'avion à partir de la connaissance de la position 3D et du module de la vitesse bistatique de l'avion dans au moins trois des systèmes bistatiques précités. Ce calcul est illustré par le schéma de la figure 2. Dans ce schéma, qui, pour des raisons de clarté du dessin, ne représente que deux des trois systèmes bistatiques choisis dans un plan coupant deux ellipsoïdes, et qui est le plan du dessin. Dans ce plan, on a tracé les ellipses el l et e12 relatives aux émetteurs et récepteurs Eml, Em2, R1, R2, l'avion à localiser étant situé à l'une des deux intersections de ces ellipses dont le foyer commun est R1. De façon connue en soi, on calcule les valeurs des modules v1 et v2 des vitesses Doppler bistatiques de l'avion pour les ellipses el1 et e12, et on reporte sur le schéma les segments correspondants de telle façon que le segment v1 soit la bissectrice de l'angle formé par les segments (avion-Eml), (avion-R1) et v2 la bissectrice de l'angle formé par les segments (avion-RI), (avion-Em2). On trace ensuite la droite DI(qui est la trace, seule visible dans le plan du dessin, du plan PLI perpendiculaire au segment v1 et passant par l'extrémité de celui-ci), puis la droite D2 (qui est la trace, seule visible dans le plan du dessin, du plan PL2 perpendiculaire au segment v2 et passant par l'extrémité de celui-ci). Les droites Dl et D2 se coupent en un point P. Ces deux plans PLI et PL2, qui passent respectivement par Dl et D2, sont les lieux géométriques des positions potentielles de l'extrémité du vrai vecteur vitesse de l'avion (et ayant pour origine cet avion) et se coupent selon une droite D1,2 passant par le point P et perpendiculaire au plan du dessin. On procède de même pour un autre couple de vitesses bistatiques (par exemple v2 et une troisième vitesse v3), et on obtient un premier plan PL2 passant par D2 et un autre plan PL3 (qui n'est pas dans le plan du dessin). Ces deux plans PL2 et PL3 se coupent selon une autre droite D2,3 (non visible sur le dessin) qui, ellemême, coupe D1,2 en un point Pl (non visible sur le dessin) L'extrémité du vrai vecteur vitesse est ainsi ce point P1. Après cette phase d'initialisation, lorsque l'avion a été localisé et son vecteur vitesse déterminé, on entretient le processus de poursuite de l'avion pour être à même de déterminer le plus précisément et le plus facilement possible à chaque instant la position et la vitesse de l'avion. En pratique, un système multistatique comporte généralement un grand nombre d'émetteurs et de récepteurs, et plus ce nombre est élevé, plus il existe de systèmes bistatiques intrinsèques, et donc un nombre encore plus élevé de distances bistatiques calculées. Il en résulte une complexité encore plus rapidement croissante du procédé de localisation de l'invention. Pour résoudre ce problème, une autre caractéristique de l'invention permet de simplifier le procédé de calcul d'entretien grâce à la prédiction de la position de l'avion. Cette prédiction se fait de la façon suivante. A partir de la connaissance de la position 3D de l'avion et de son vecteur vitesse à l'instant t, on calcule la position de l'avion à l'instant t+At (la valeur de At dépendant en particulier de la vitesse et de la trajectoire de l'avion et de la précision des mesures, ce afin de ne pas prendre At trop élevé pour éviter une prédiction trop approximative, ou trop petit pour ne pas faire de calculs intermédiaires inutiles). On calcule ainsi, comme exposé ci-dessus, l'altitude de l'avion à l'instant t+At, à savoir Z prédit. Connaissant z prédit, on refait tous les calculs précédents (depuis la première étape de calcul des points d'intersection de couples d'ellipsoïdes jusqu'au calcul du vecteur vitesse 3D), en prenant à la première étape une fourchettes de valeurs de z (zl, z2, z3) beaucoup moins large que celle mentionnée ci- dessus en exemple, et maintenant centrée autour de z prédit, par exemple 0,8 z prédit, Z prédit et 1,2 z prédit On simplifie les calculs en prenant parmi toutes les valeurs de distances bistatiques disponibles de préférence les trois les plus faibles (présentant le meilleur rapport 10 signal/bruit). Parmi les applications possibles de l'invention, on peut citer son utilisation, dans un système multistatique de détection et de localisation de cibles à faible coût, destiné à la détection dans une zone de surveillance, par exemple au voisinage de la zone d'approche d'un aéroport, d'intrus à faible surface équivalente radar, ou destiné à la protection d'un site sensible. Le dispositif permettant d'effectuer les mesures de distances bistatiques et de vitesse Doppler bistatiques est basé sur l'exploitation d'un ou de plusieurs récepteurs, associés, simultanément ou non, à : - un ensemble d'émetteurs coopératifs (émettant une forme d'onde déterministe et maîtrisée), - un ensemble d'émetteurs classiques de diffusion à modulation de fréquence, - un ensemble d'émetteurs de diffusion à formes d'onde OFDM (tels que des émetteurs DAB ou DVBT), - un ensemble d'émetteurs de diffusion dont les signaux émis 25 comportent des raies à position déterministe, par exemple des émetteurs de télévision analogique. En conclusion, l'invention permet, y compris en situation de détection d'objets volants multiples, d'obtenir une mesure pratiquement instantanée de la position des objets à localiser en 3D (en coordonnées dans un plan horizontal et en altitude) et de leur vecteur vitesse en module et en direction (en 3D). La précision obtenue grâce au procédé de l'invention est directement liée à la précision de mesure des distances bistatiques. Dans le cas de configurations typiques de distribution géométrique (c'està-dire des configurations pour lesquelles les émetteurs et récepteurs sont disposés aux sommets de figures géométriques simples, telles que des rectangles), une précision de distance bistatique élémentaire de quelques mètres permet d'atteindre des précisions dans un plan horizontal de quelques mètres, et des précisions en altitude de l'ordre de grandeur d'une dizaine de mètres. REVENDICATIONS 1. Procédé de localisation d'un objet mobile à partir de mesures de distance et de vitesse bistatiques de cet objet à l'aide d'un système multistatique, caractérisé en ce que l'on associe au moins deux couples de distances bistatiques tels que pour chaque couple les ellipsoïdes qui les déterminent aient un foyer commun, que l'on résout pour chaque couple l'intersection de ces deux ellipses à différentes altitudes, que l'on optimise par rapport à la variable d'altitude z les différentes intersections, qu'ayant trouvé la meilleure valeur de z, on calcule les coordonnées cartésiennes en trois dimensions de cet objet, que, connaissant les coordonnées de l'objet et au moins trois vitesses bistatiques, on déduit la direction et le module du vrai vecteur vitesse 3D de l'objet mobile. 2. Procédé selon la revendication 1, caractérisé en ce que l'on prédit la position de l'objet mobile à l'étape suivante du calcul, et donc son altitude z prédite, et que l'on affine l'optimisation du calcul de z autour du z prédit en ne gardant dans tous lesdits couples de distances bistatiques que les n plus petites. 3. Procédé selon la revendication 2, caractérisé en ce que l'on prend n=3.