Brevet FR3025030A1 — publié le 20160226

PROCEDE ET DISPOSITIF POUR ESTIMER LA TRAJECTOIRE D'UN OBJET EN MOUVEMENT L'invention concerne un procédé pour estimer la trajectoire d'un 5 objet par mesures d'angles et de la connaissance de la trajectoire d'un système équipé d'un moyen de mesures d'angle. Elle peut être utilisée pour déterminer la trajectoire d'objets qui communiquent par voie radioélectrique avec d'autres dispositifs dans un système de réseau de radiocommunications. Le procédé sert notamment à déterminer les 10 paramètres de la trajectoire tels que les paramètres d'accélération, etc. Le procédé trouve son utilisation dans le domaine de la radiogoniométrie. La localisation d'un objet en déplacement avec une vitesse similaire à celle d'un système de détection rencontre certaines difficultés techniques. La technique connue de localisation par densité ne conduit pas à 15 des résultats satisfaisant du fait notamment des vitesses similaires entre l'objet à détecter et le système de détection. Il est connu d'utiliser des capteurs électrostatiques pour déterminer la trajectoire d'objets en déplacement. Il est aussi connu d'employer des capteurs optroniques pour déterminer la trajectoire avec laquelle un objet se déplace. 20 Certaines méthodes connues de l'art antérieur, par exemple celle décrite dans le brevet FR 2833084, s'appliquent pour les cas où l'objet dont on cherche à déterminer la trajectoire se déplace selon un mouvement rectiligne uniforme. De l'avis du demandeur, les techniques décrites dans l'art 25 antérieur sont mieux adaptées à des cas particuliers et ne fournissent pas de méthode de trajectographie pour des applications plus générales, pour des trajectoires non rectilignes: l'art antérieur ne propose pas de méthode permettant de détecter le type de trajectoire suivie par l'objet dont on veut obtenir une estimation de trajectoire. 30 Le procédé selon l'invention repose notamment sur l'utilisation d'une quantité de mesures d'angles limitée pour estimer la trajectoire 3025030 2 quelconque d'un objet. Les mesures sont notamment effectuées à l'aide d'un goniomètre. Dans la suite de la description, on rappelle tout d'abord la notion complète d'observabilité déterministe : « Un système sans bruit d'état ni de 5 mesure est considéré comme complètement observable sur un intervalle temporel [fini tfin] si les mesures disponibles durant cette période permettent d'obtenir l'état initial de l'objet observé de façon unique. ». Ensuite, on introduit la notion de potentiel d'observabilité d'une plateforme qui correspond à la capacité à détecter des mouvements d'objets. Ce potentiel 10 est directement lié au mouvement du système équipé du senseur. Enfin on appelle complexité d'une trajectoire le nombre minimum de dimension nécessaire afin de décrire une trajectoire donnée. Par exemple, un point peut être simplement décrit par un triplet (x,y,z) soit une complexité de 1. Aussi, un mouvement rectiligne uniforme est complètement caractérisé par un point 15 initial (xo,yo,zo) et une vitesse initiale (vxo,vyo, vzo) soit une complexité de 2. Il existe bien des façons de compresser l'information pour décrire une trajectoire. De fait le choix des fonctions de compression (par exemple les polynômes d'interpolation) a un impact sur la complexité. On appelera transformation compressive toute opération permettant de compresser la 20 donnée que représentent les points d'une trajectoire. Il faut noter que les transformations compressives sont, par essence, non bijectives. On entendra donc par le terme « transformation inverse » la décompression d'une trajectoire à partir des paramètres compressés. L'invention concerne un procédé pour estimer une trajectoire d'un 25 premier objet au moyen d'une plateforme en déplacement équipée de moyens de mesure d'au moins un angle entre le premier objet et la plateforme selon une direction donnée comportant au moins les étapes suivantes : - Déterminer un ensemble de points correspondant à la trajectoire de la 30 plateforme pour un intervalle de temps donné [tl,..tN], 3025030 3 - Acquérir un ensemble de valeurs d'angles sur le même intervalle de temps [ti,..tN]' - Définir la complexité Cg de la trajectoire de la plateforme en calculant l'erreur existant entre les points de la trajectoire partitionnés en plusieurs groupes de points et des points de la trajectoire obtenus par extrapolation après compression de la trajectoire, - Déterminer pour le même intervalle de temps [tl,..tN], la complexité cp de la trajectoire du premier objet à partir d'une valeur de complexité estimée Cg de la trajectoire de la plateforme et de l'erreur calculée entre des valeurs d'angles 0 = [8,,..., 9N] mesurées et partitionnés en K groupes et des valeurs d'angles obtenues par extrapolation d'une trajectoire estimée de l'objet, - Comparer la valeur de la complexité Cg de la plateforme et si la valeur de la complexité de la plateforme est supérieure à la valeur de la complexité Cp de l'objet, alors estimer la trajectoire de l'objet, sinon, transmettre un signal ou une information pour ajuster la trajectoire de la plateforme. Le procédé peut comporter les étapes suivantes : a) estimer la complexité Cg associée à la trajectoire de la plateforme en distribuant, pour une valeur de complexité définie, les points de la trajectoire de façon aléatoire en K groupes, appliquer une transformation compressive et son inverse sur un groupe de points pour obtenir un deuxième groupe de points pour chaque complexité, et extrapoler les points, puis déterminer l'erreur existant sig entre les points définissant la trajectoire de la plateforme et les points obtenus par extrapolation, réitérer ces étapes en travaillant sur un ensemble de groupes de points par tirage aléatoire et sans recouvrement, puis définir une deuxième valeur d'erreur 52g, définir la valeur d'erreur minimale et la transformation compressive associée, b) estimer la complexité Cp de la trajectoire de l'objet en mouvement à partir 30 d'une valeur de complexité estimée de la plateforme, en construisant K groupes d'angles par tirage aléatoire sans répétition, pour chaque groupe k, 3025030 4 en appliquant une transformation compressive inverse afin de déterminer une première trajectoire pour l'objet et en extrapolant la valeur de la première trajectoire pour les indices temporels propres au groupe k, en reconstruisant les valeurs d'angles observées par la plateforme si l'objet avait suivie cette 5 deuxième trajectoire extrapolée, et en calculant la valeur de l'erreur pour obtenir un premier score 51p, puis en réitérant ces étapes pour un ensemble de groupes d'angles obtenus par un découpage dans l'ordre des indices des angles et sans recouvrement et pour chaque valeur de complexité choisie, afin d'obtenir un deuxième score, définir la valeur minimale, 10 c) comparer la valeur de la complexité Cg de la trajectoire de la plateforme et la valeur de la complexité Cp de la trajectoire de l'objet et si Cg est supérieure à Cp alors exécuter un algorithme d'estimation de la trajectoire de l'objet, sinon modifier la complexité de la trajectoire de la plateforme. Les valeurs d'erreur peuvent être déterminées en utilisant une 15 méthode de calcul d'erreurs quadratiques. On mesure des valeurs d'angles par goniométrie ou radiogoniométrie, pour un objet équipé d'un émetteur. Dans certains cas d'utilisation, on utilise une signature de l'objet afin de déterminer les angles. 20 Selon un mode de mise en oeuvre du procédé, l'objet est en mouvement par rapport à la plateforme. L'invention concerne aussi un système pour estimer une trajectoire d'un premier objet en mouvement par rapport à une plateforme en déplacement équipé de moyens de mesure d'au moins un angle entre le 25 premier objet et la plateforme selon une direction donnée caractérisé en ce que la plateforme comporte des moyens adaptés à exécuter les étapes du procédé selon l'invention. D'autres caractéristiques et avantages du procédé selon l'invention apparaîtront mieux à la lecture de la description qui suit d'un exemple de réalisation donné à titre illustratif et nullement limitatif annexé des figures qui représentent : 3025030 5 - La figure 1 un exemple de système pour la mise en oeuvre du procédé selon l'invention, - La figure 2, une représentation dans un repère x, y, z, des paramètres utilisés pour la résolution du problème de trajectographie passive, 5 - La figure 3, un premier schéma bloc représentant les étapes globales du procédé selon l'invention, - La figure 4, une représentation des étapes pour déterminer la complexité de la trajectoire d'une plateforme équipée d'un dispositif de mesures d'angle et d'une plateforme de navigation, 10 - La figure 5, une représentation des étapes permettant de déterminer la complexité d'un objet en mouvement. La description qui suit est donnée à titre d'exemple illustratif et nullement limitatif pour estimer la trajectoire d'un objet en mouvement équipé d'un émetteur radioélectrique grâce à une plateforme équipée d'un 15 radiogoniomètre, la plateforme effectuant une ou plusieurs manoeuvres entrainant un mouvement dont la complexité telle que définie précédemment devient supérieure à celle de l'objet. Le procédé permet par ailleurs d'indiquer à la plateforme si sa manoeuvre est suffisante pour estimer la trajectoire de l'objet et de proposer une idée de manoeuvre optimale afin 20 d'améliorer la précision de localisation de l'objet. Le procédé peut être utilisé dans le cas où l'on estime que la trajectoire d'un objet peut être représentée par un point. La figure 1 représente un exemple de système pour la mise en oeuvre du procédé selon l'invention. Le système 1 comporte un premier objet 25 10 dont on cherche à estimer la trajectoire, comprenant un émetteur de signaux et un système 11 ou plateforme équipée d'une centrale de navigation 12 permettant de connaître au cours du temps la trajectoire de la plateforme et la position d'un goniomètre 13 ou d'un radiogoniomètre permettant la mesure angulaire. La plateforme comporte un capteur 14 30 recevant le signal en provenance au moins de l'objet observé 10. Les informations de la centrale de navigation 12 et les informations de 3025030 6 goniométrie fournies par le goniomètre 13 sont redirigées vers un module 15 de la plateforme 11. Le module 15 comprend un processeur 16 qui est adapté à traiter les mesures angulaires et à exécuter les étapes mises en oeuvre par le procédé selon l'invention. La plateforme 11 est aussi munie 5 d'un espace mémoire 17 de stockage des données. Le procédé peut aussi être utilisé pour estimer la trajectoire d'objets qui n'ont pas d'émetteur, mais qui peuvent être identifiés, suivis, au moyen d'une signature qui leur est propre et qui permet de faire une mesure d'angle par rapport à un observateur. 10 Le procédé selon une variante de mise en oeuvre comporte en résumé, par exemple, une première phase d'initialisation au cours de laquelle on va déterminer les caractéristiques intrinsèques au radiogoniomètre, les caractéristiques du problème à observer, par exemple, la plage de vitesses des objets à observer, la plage de distance à couvrir 15 autour de la plateforme, etc.). Une deuxième phase active comprend par exemple les étapes suivantes : - On effectue l'acquisition de mesures angulaires via le radiogoniomètre de la plateforme pendant une durée T donnée, en parallèle on collecte 20 les données de navigation de la plateforme sur la même durée T, - A partir des données de navigation, on caractérise le mouvement de la plateforme équipée du radiogoniomètre. Ceci permet d'évaluer la capacité Cg de la plateforme à détecter des mouvements d'objets plus ou moins complexes, ou potentiel d'observabilité, 25 - On estime la complexité Cp de la trajectoire de l'objet en tenant compte des valeurs d'angle mesurées, - On compare la valeur de la complexité Cg de la trajectoire de la plateforme et la valeur de la complexité Cp de la trajectoire de l'objet, - Si la valeur Cg est supérieure à Cp alors, on estime la trajectoire de 30 l'objet en mouvement, 3025030 7 - Si la valeur Cg est inférieure ou égale à Cp, alors il est possible de proposer un ajustement de la manoeuvre de la plateforme afin d'obtenir un potentiel d'observabilité suffisant pour estimer la trajectoire de l'objet, puis on estime la trajectoire de l'objet. 5 La figure 2, illustre le problème de la trajectographie passive dans un repère à trois dimensions, les paramètres utilisés, avec capteurs multiples et trajectoire curviligne, non planaire. Le problème de la trajectographie posée dans la présente invention peut être énoncé comme suit : on considère un objet dont on 10 souhaite connaître la trajectoire avec une capacité d'émission radio p et une plateforme munie d'un radiogoniomètre g ou récepteur. La notion de trajectoire à deux dimensions est définie comme un ensemble de fonctions continues t E [0;l] ---> C OÙ C désigne le corps bien connu des complexes. 15 Les trajectoires de l'objet émetteur et de la plateforme récepteur sont définies dans le cas continu: {p : t -> p(t) g :t -> g (t) Les fonctions sont échantillonnées sur N échantillons pour obtenir les valeurs de la trajectoire à temps discret : 20 [P1,--,PN] G Où: p' = [x' , y' ]T avec [ l'opérateur transposé. Lorsque P ou G sera représenté en gras et sans indice, l'ensemble des positions aux instants de mesures des angles est considéré ici sous la forme 25 d'une matrice de dimension 2 x N . On définit une fonction de mesure d'angle entre p' et g,, par rapport à une direction fixée, par exemple le nord magnétique, On = 4pn - g' 3025030 8 L'ensemble des mesures d'angles est noté 0 = [0'...,ON] Un bruit ri additif gaussien vient s'ajouter à la mesure d'azimuts, la séquence bruitée de mesures est la suivante : Ô=O+ri 5 Le procédé selon l'invention va permettre d'estimer la trajectoire de l'objet en mouvement P connaissant la trajectoire de la plateforme G et les valeurs d'angles Ô : 10 La figure 3 représente selon un schéma bloc les étapes principales du procédé selon l'invention. Les étapes sont exécutées au niveau du module 15 de la plateforme 11. Le procédé prend en compte les positions du goniomètre G = les directions d'arrivées 0 = [8,,..., ON] des signaux reçus sur le goniomètre et les instants [tI,...,tN] de 15 mesures= [8j,...,0N], valeurs qui sont par exemple mémorisées au niveau de la plateforme, 30. Le procédé va caractériser, 31, le mouvement du goniomètre de la plateforme par une méthode de validation croisée de type I décrite ci-après, et déterminer la complexité Cg du goniomètre, 32. Puis le procédé va caractériser, 33, le mouvement de l'objet dont on cherche à 20 connaître la trajectoire en tenant compte du mouvement du goniomètre, de la complexité du goniomètre ou de la plateforme par une méthode de validation croisée 2. On compare ensuite la valeur de la complexité Cg de la trajectoire de l'objet à celle de la trajectoire du goniomètre C9. Si le goniomètre a une complexité de trajectoire supérieure à celle de l'objet à observer, alors, 34, 25 35, le procédé va déterminer les valeurs de positions de l'objet observé pour différents instants. Dans le cas contraire, le procédé va transmettre à l'opérateur un signal ou message l'avertissant que la complexité de la trajectoire de la plateforme n'est pas suffisante pour permettre de déterminer 3025030 9 la trajectoire de l'objet à observer. L'opérateur peut alors ordonner le changement de la trajectoire de la plateforme. La figure 4 schématise les étapes du procédé mises en oeuvre pour la Validation croisée de type 1 permettant l'estimation de la complexité 5 Cg de la trajectoire du goniomètre. Le procédé prend en entrée les positions mesurées du goniomètre aux instants [ti,...,tN]pour la plateforme. Soit G = [g'...,gN] l'ensemble des positions du goniomètre, Les valeurs des complexités possibles pour la trajectoire de la plateforme 10 équipée du goniomètre sont choisies dans un ensemble de complexité de la trajectoire défini par un utilisateur et l'application souhaitée, Cgmax]. Pour chaque transformation compressive T, on va exécuter les étapes décrites ci-après. Les transformations compressives peuvent être choisies, par exemple, dans un dictionnaire. On peut compresser dans l'espace de 15 Fourier : cet espace est adapté aux trajectoires circulaires. De manière générale la transformation compressive sera choisie en fonction du type supposé de trajectoire (circulaire, continue par morceaux, fractal...). a. boucle sur la complexité 401- Pour chaque complexité Cg allant de 1 à Cg.. 20 i. Construire K groupes de taille Cg à partir de l'ensemble des positions G =[gt,...,gis] de la façon suivante : 402, 403 - choisir aléatoirement des points dans l'ensemble G et les distribuer dans les K groupes, G = [g j(l), , g j(N)1 Gk = [glik (1), ,glik (c.01 La distribution peut se faire avec recouvrement de points entre les groupes : plusieurs groupes 25 peuvent contenir le même point. Le choix aléatoire peut être une permutation aléatoire. L'étape avec partitionnement aléatoire en K groupes avec répétition/recouvrement entre groupe conduit à b- k ii. 404, Pour tous les groupes formés, k variant de 1 à K : 3025030 10 ii.1. Effectuer la transformation compressive à partir des Cg points du kième groupe, la transformation compressera la trajectoire selon Cg dimensions : Gk -> Wk qui correspond à la version compressée de la trajectoire 5 définie par les points contenus dans le k-ème groupe Gk , 405, ii.2. Effectuer la transformation compressive inverse pour obtenir un ensemble de points Gk , 406, ii.3.Extrapoler les points de Gk, 407, à tous les indices temporels disponibles 10 ii.4.Calculer, 408, l'erreur entre les points ainsi extrapolés et les points des autres groupes qui n'appartiennent pas à Gk , le calcul pouvant être un calcul d'erreur quadratique, On obtient un premier score sig (somme des erreurs quadratiques), 409, pour chaque groupe et pour chaque valeur de complexité. 15 Partitionner 410, G =[gL(1)' g L(N)1, 411 aléatoirement la trajectoire G =[gi,...,gN] en K groupes contenant[ -N points sans recouvrement, K [on note âkfl ces groupes, Gk = gLk(i), iii. 412, Pour tous les groupes formés, k variant de 1 à K 413, Sur chaque complémentaire des groupes, points qui n'appartiennent 20 pas au groupe considéré et notés Gk, iv.1.Effectuer la transformation compressive puis son inverse pour extrapoler, 414, 415 , Gk vers Wk puis Wk vers Gk, iv.2.Calculer l'erreur quadratique 416 entre chacun des groupes complémentaires et le groupe Gk' transformé. 3025030 11 417, On obtient le score s2g (somme des erreurs quadratiques) pour tous les groupes et pour toutes les hypothèses de valeurs de complexité, b. 418, Sortir le score final sfg pour la plateforme lié à toutes les complexités en sommant les deux scores sfg et s2g, 5 c. Choisir le meilleur score smg (le score minimal) et sortir le couple (Transformée, Cg complexité de la trajectoire de la plateforme) qui minimise les smg. Cette valeur de complexité estimée Cg va être utilisée pour déterminer la trajectoire de l'objet en mouvement. 10 Lors du calcul d'erreur toute méthode autre que la méthode d'erreur quadratique peut être utilisée, sans sortir du cadre de l'invention, ceci pour toutes les étapes de calcul d'erreur lors de la validation croisée de type 1 ou lors de la validation croisée de type 2. La figure 5 représente les étapes de la validation croisée de type 2 15 pour l'estimation de la complexité de la trajectoire de l'objet les entrées du bloc sont les suivantes, 500 : - Les angles 0 4,...,ON]estimés par le radiogoniomètre aux instants k'..., tN - Les positions G = ki,...,gAdu goniomètre aux instants , 20 - La complexité estimée du goniomètreÔg ou une complexité maximum telle que Cquekonque Cg valeur estimée lors de l'étape de validation croisée de type 1j - Le type de transformation compressive. Pour chaque complexité C, allant de 1 à Cs, , 501 [ 25 a.Construire K groupes d'angles de taille . Pour construire ces groupes, K on tire aléatoirement les angles dans l'ensemble 0 = [01,...,0N] sans répétition, 502, b.Pour chaque groupe k, k=1 à K , 503 3025030 12 i. 504, Effectuer l'estimation de la trajectoire de l'objet, en supposant un type de trajectoire de complexité Cp, à l'aide des angles qui ne sont pas dans le k-ième groupe (complémentaire du groupe k) et des positions du goniomètre à ces instants, 5 On obtient après transformation compressive inverse une trajectoire estimée PkCP pour l'objet et pour un groupe k, pour une valeur de complexité, ii. 505, On extrapole la trajectoire PkCP estimée aux indices temporels propres au k-ième groupe d'angle. --^^^ 10 On obtient pour une valeur de complexité la trajectoire complète Picc° , iii. 506, A partir de pkci, et de la trajectoire du goniomètre on reconstruit les angles qu'auraient estimé le goniomètre si l'objet observé était passé par tous les points de la trajectoire estimée complète licp , iv. 507, On calcule l'erreur quadratique entre les angles du groupe k et les 15 angles déterminés à l'aide de Cp aux mêmes indices temporels, On obtient un score pour le couple (k, Cp ) correspondant à cette scoremsi(k,Cphypoll. ) erreur, . c. 508, Une fois que l'on a traité tous les groupes à Cp fixé, on obtient le score s1p en effectuant la moyenne des scores sur les couples (k, Cp ); puis 20 le procédé va, d. 509, Construire K groupes de taille [-1 V] à partir de [81,...,ON] . Pour construire ces groupes, on découpe l'ensemble des angles mesurés [e1,...,ON] dans l'ordre des indices et sans recouvrement en K groupes de taille , 25 e. 510 Pour chaque groupe k, k=1 à K,: y. 511, Effectuer l'estimation de la trajectoire à l'aide des angles qui ne sont pas dans le k-ième groupe (complémentaire du groupe k) et des positions du goniomètre, 3025030 On obtient après transformation trajectoire i3 P, 13 vi. 512, On extrapole la trajectoire compressive inverse de l'estimé une ième groupe d'angle, hkcâux indices temporels propres au k- 5 On obtient &P complète, vii. 513, A partir de bkCP et de la trajectoire du goniomètre, on reconstruit les angles qu'aurait estimé le goniomètre si l'objet observé était passé par tous les points de 15kcp , viii. 514, On calcule l'erreur quadratique entre les angles du groupe k et les 10 angles déterminés à l'aide de aux mêmes indices temporels,. scoreRms2(k,Cphypo,'s) On obtient alors un autre type de score pour le couple (k, Cp), f. Une fois que l'on a traité tous les groupes à Cp fixé, on obtient le score s2p en effectuant la moyenne des scores sur les couples (k, Cp ), 515, après avoir 15 traité tous les Cp on dispose de Cp scores s2p, 2. Le score global est obtenu en sommant les scores s1p et s2p, 3. Cp est obtenu en minimisant la fonction (s1p et s2p ) qui dépend des Cp testés et du score associé, 516, Connaissant les valeurs de complexité pour la plateforme et pour 20 l'objet, le procédé peut ensuite dérouler les étapes présentées à la figure 3 afin de déterminer la trajectoire de l'objet ou encore d'avertir un opérateur afin de modifier la trajectoire de la plateforme. Les étapes décrites dans le cadre de mesures d'angle d'azimut s'appliquent aussi avec des angles d'élévation pour une trajectoire en 3 25 dimensions ou d'autres mesures angulaires pouvant être utilisées. Le procédé a défini la notion d'observabilité qui est une grandeur binaire : lorsqu'elle est vérifiée, il est possible d'estimer la trajectoire; lorsqu'elle n'est pas vérifiée, il n'est pas possible de converger vers une solution unique ou un espace contenant cette solution. Cette grandeur est 30 définie, par exemple, de la façon suivante : 3025030 14 = o si cg cg 11 si Cg > Cp Puis en fonction des conditions d'estimation : - Si les conditions sont bonnes: on estime la trajectoire via une méthode asymptotiquement non biaisée connue de l'homme du métier (variable 5 instrumentale VI, variable instrumentale pondérée VIP, variable instrumentale pondérée compensée du biais ou en anglo-saxon Bias Compensated Weightened Instrumental Variable BCWIV ou ML méthode de vraisemblance maximum), - Si les conditions sont mauvaises on donne des conseils de manoeuvre 10 élaborés pour augmenter la complexité de la trajectoire du dispositif de détection, un goniomètre par exemple. L'estimabilité est, quant à elle, une grandeur graduelle qui fournit une indication sur les performances d'estimation optimales que l'utilisateur peut espérer avoir. Il est possible de choisir la Borne de Cramer-Rao (BCR) 15 connue de l'homme du métier ou tout autre outil mathématique équivalent. Elle se définit à partir de la matrice d'information de Fisher J de l'estimateur : e = tr(I-1) où l'opérateur tr(.)désigne la trace d'une matrice soit la somme des éléments diagonaux. Pour estimer la valeur de la complexité de l'objet observé pendant 20 l'étape de validation croisée de type 2, le procédé peut mettre en oeuvre un estimateur pseudo-linéaire. On définit à partir des notations vectorielles scalaires pour les points de trajectoires une notation complexe: p n = [xn .n prplexe = xn iyn x {pcomplexe _ 25 Si bien que le vecteur x = s'écrit G"'"'" }et Y Z{pCOMPleXe G complexe } où J?{ } et 3{ } désignent respectivement la partie réelle et imaginaire. Il est alors possible de déterminer un espace orthogonal sur lequel la projection de ces coordonnées est nulle : 3025030 ( - sinK 0 ... 15 0 sin(02) 0 , 0 0 -sin(ON) =0 (x y) + cos(00) 0 0 0 + cosK ' - 0 o o COSK y)( diag(sin(0))1 + diag(cose = 0 où l'opérateur diag désigne la matrice carré dont les éléments diagonaux sous les éléments du vecteur en paramètre. 5 On peut réécrire cette équation après application de la transformation compressive T(x y). (cox ffl). (7' - diag(sin(0))1 0 (oex oy 0 T diag(cos(0))) Soit encore (cox co y ( Tdiag(sin(01 + Tdiag(cosW , 0 10 Que l'on peut mettre sous la forme : wy)3± = ° Ceci fournit déjà un espace contenant la solution qui est celui de du noyau Ol . Cet espace donnera des indications sur la solution dans le cas où la complexité est insuffisante pour l'estimation. Une information 15 supplémentaire sur la vitesse ou distance entre le capteur et l'objet fixerait un espace pour les trajectoires compatibles. Dans le cas où l'observabilité est assurée, l'estimateur pseudo-linéaire ou EPL permet d'effectuer l'estimation. L'EPL peut être effectué avec tout type de transformation compressive linéaire suggérée par la forme de la 20 trajectoire et les tests effectués. 3025030 16 De plus il est possible de définir assez simplement à partir des équations ci-dessus la matrice d'information de Fisher correspondante : J = (01)diag(p)-1Ecliag(p) 1(Olf où : 5 - p = P - G sont les distances entre l'objet et le capteur - E est la matrice de covariance du bruit de mesure sur les angles Cette matrice peut bien évidement être calculée avec toutes les approximations possibles : 10 - Limitation de la complexité à min(Cp,Cg) - Distance entre le capteur et l'estimée au lieu de la vraie distance - Matrice de covariance du bruit de mesure estimée Les trajectoires estimées par le procédé selon l'invention sont 15 alors transmises pour être, soit affichées à un opérateur, soit injectées dans un système de niveau supérieur (système d'information ou autre). Elles contribuent ainsi à la connaissance de la situation radioélectrique entourant le porteur. Sans sortir du cadre de l'invention, le procédé peut s'appliquer 20 pour estimer la trajectoire d'un objet reconnaissable grâce à sa signature et non équipé d'émetteur pour la mise en oeuvre de l'invention. Toute méthode de calcul d'erreur entre valeurs peut être utilisée. Il en est de même pour les transformations compressives permettant de diminuer le nombre de données mises en oeuvre dans les calculs qui seront choisies, par exemple en 25 fonction du type de trajectoire. La description a utilisé à titre d'exemple un estimateur pseudo linéaire EPL, mais le procédé pourrait mettre en oeuvre tout type d'estimateur adapté. L'invention présente l'avantage d'être utilisée pour différentes plateformes en mouvement équipée d'un radiogoniomètre ou un dispositif 30 similaire. Elle peut être utilisée dans tous types de plateformes mouvantes 3025030 17 munies d'un radiogoniomètre. Elle permet d'estimer le degré de complexité de la trajectoire d'un objet de manière automatique, qu'il s'agisse d'un point, d'une droite ou d'un mouvement plus complexe. 5 REVENDICATIONS1 - Procédé pour estimer une trajectoire d'un premier objet (10) au moyen d'une plateforme (11) en déplacement équipée de moyens de mesure (13) d'au moins un angle entre le premier objet (10) et la plateforme selon une direction donnée comportant au moins les étapes suivantes : - Déterminer un ensemble de points correspondant à la trajectoire de la plateforme pour un intervalle de temps donné - Acquérir un ensemble de valeurs d'angles sur le même intervalle de temps [tl,..tN], - Définir la complexité Cg de la trajectoire de la plateforme (11) en calculant l'erreur existant entre les points de la trajectoire partitionnés en plusieurs groupes de points et des points de la trajectoire obtenus par extrapolation après compression de la trajectoire, - Déterminer pour le même intervalle de temps [tl,..tN], la complexité cp de la trajectoire du premier objet (10) à partir d'une valeur de complexité estimée Cg de la trajectoire de la plateforme et de l'erreur calculée entre des valeurs d'angles 0= [8,,...,ON] mesurées et partitionnés en K groupes et des valeurs d'angles obtenues par extrapolation d'une trajectoire estimée de l'objet, - Comparer la valeur de la complexité Cg de la plateforme et si la valeur de la complexité de la plateforme est supérieure à la valeur de la complexité cp de l'objet, alors estimer la trajectoire de l'objet, sinon, transmettre un signal ou une information pour ajuster la trajectoire de la plateforme. 2 - Procédé selon la revendication 1 caractérisé en qu'il comporte au moins les étapes suivantes : a) estimer la complexité Cg associée à la trajectoire de la plateforme en distribuant, pour une valeur de complexité définie, les points de la trajectoire de façon aléatoire en K groupes, appliquer une transformation compressive 3025030 19 et son inverse sur un groupe de points pour obtenir un deuxième groupe de points pour chaque complexité, et extrapoler les points, puis déterminer l'erreur existant sig entre les points définissant la trajectoire de la plateforme et les points obtenus par extrapolation, réitérer ces étapes en travaillant sur 5 un ensemble de groupes de points par tirage aléatoire et sans recouvrement, puis définir une deuxième valeur d'erreur 52g, définir la valeur d'erreur minimale et la transformation compressive associée, b) estimer la complexité Cp de la trajectoire de l'objet en mouvement à partir d'une valeur de complexité estimée de la plateforme, en construisant K 10 groupes d'angles par tirage aléatoire sans répétition, pour chaque groupe k, en appliquant une transformation compressive afin de déterminer une première trajectoire pour l'objet et en extrapolant la valeur de la première trajectoire pour les indices temporels propres au groupe k, en reconstruisant les valeurs d'angles observées par la plateforme si l'objet avait suivie cette 15 deuxième trajectoire extrapolée, et en calculant la valeur de l'erreur pour obtenir un premier score 51p, puis en réitérant ces étapes pour un ensemble de groupes d'angles obtenus par un découpage dans l'ordre des indices des angles et sans recouvrement et pour chaque valeur de complexité choisie, afin d'obtenir un deuxième score, définir la valeur minimale, 20 c) comparer la valeur de la complexité Cg de la trajectoire de la plateforme et la valeur de la complexité Cp de la trajectoire de l'objet et si Cg et supérieure à cp alors exécuter un algorithme d'estimation de la trajectoire de l'objet, sinon modifier la complexité de la trajectoire de la plateforme. 25 3 - Procédé selon la revendication 1 caractérisé en ce que l'on calcule les valeurs d'erreur par une méthode de calcul d'erreurs quadratiques. 4 - Procédé selon l'une des revendications précédentes caractérisé en ce que l'on mesure les angles par goniométrie ou radiogoniométrie. 30 3025030 20 5 - Procédé selon l'une des revendications 1 à 3 caractérisé en ce que l'on utilise une signature de l'objet afin de déterminer les angles. 6 -Procédé selon l'une des revendications 1 à 5 caractérisé en ce que l'objet 5 est en mouvement par rapport à la plateforme. 7 - Système pour estimer une trajectoire d'un premier objet (10) en mouvement par rapport à une plateforme (11) en déplacement équipé de moyens de mesure (13) d'au moins un angle entre le premier objet (10) et la 10 plateforme selon une direction donnée caractérisé en ce que la plateforme (11) comporte des moyens adaptés à exécuter les étapes du procédé selon l'une des revendications 1 à 5. 8 - Système selon la revendication 7 caractérisé en ce que la plateforme 15 (11) est équipée d'un goniomètre.