FR 2999724 A1 20140620 FR 1203408 A 20121214 2999724 PROCEDE POUR DETERMINER LA POSITION D'UN EMETTEUR SANS UTILISER DE VOIE DE REFERENCE L'objet de l'invention concerne un procédé pour déterminer la position d'un émetteur au sol et son système associé sans utiliser de voie de référence. L'invention concerne notamment le domaine de la localisation passive d'émetteurs fixes ou mobiles. L'art antérieur décrit différentes méthodes pour déterminer la position d'émetteurs au sol, et pour localiser ces émetteurs à partir de porteurs défilant et d'un capteur embarqué produisant des vecteurs directeurs associés. L'un des problèmes résultant des procédés de localisation connus exploitant la signature spatiale d'un réseau antennaire est que le procédé repose sur une estimation ambiguë du vecteur directeur. Cette estimation est en effet indéterminée, à un facteur complexe près. Pour s'affranchir de cette indétermination, les procédés de localisation reposant sur la mesure de vecteur directeur Ak, aux instants k=1, K, mettent en oeuvre une première étape de normalisation consistant à choisir une voie de référence et à élaborer des vecteurs Bk construits à partir de cette référence. Une manière de procéder est, par exemple, décrite dans le brevet US 7907089 du Demandeur. Le procédé détermine la position d'un émetteur à partir des composantes des vecteurs al jusqu'à ai( Ces vecteurs ak ont la particularité de dépendre de l'instant tk et surtout de la position (XMe yl115 Zrn) de l'émetteur. En particulier pour un réseau composé de N=2 capteurs espacés d'une distance de d dans l'axe du porteur le vecteur ak: vérifie (1) . ( ak = exp j27-1- -d cos(0(t k , x11, , y., z.))cos(A(t- k , = a(tk,xm,Ym,zrn) z.))j D'après la figure 1, l'incidence (9(tk,xm,ym,zn7)4(tk,xm,Ym,zm) ) peut être directement calculée à partir de la position (xk,Yk,zk) du porteur à l'instant tk et de la position (xm,ym,zm) de l'émetteur. Dans le procédé initial un vecteur bk est constitué à partir des composantes du vecteur ak. En particulier, le vecteur bk de dimension (N-1)x1 est construit en choisissant un capteur de référence en n=i : 2 2999724 bkm = a k',(1) I ak',(i) = b(tk, xm, ym, zm) où akm(n) est la ilème composante de akm ak',(i -1)1 ak',(i) ak', (i + 1)/ ak. (i) akm (N) 1 ak,n(i) Sachant que les vecteurs directeurs akm sont estimés avec une certaine erreur ekm telle que akm= a(tk,xm,Ym,zm)+ ekm, 1 est possible de déduire qu'il en est de même pour le vecteur transformé bkm : bkm = V(tk,Xm,ym,Zm)+Wkm où wkm = et ekm(n) est la iième composante de ekm. En particulier, pour un réseau composé de N=2 capteurs espacés d'une distance d dans l'axe du porteur, la signature théorique a vérifie pour un modèle de front d'onde plan parfait, dans l'hypothèse d'absence de couplage et amplitudes sur les différentes voies indépendantes de l'incidence : exp j2 -d 1-t- cos(0(t k , x ',, y', ,z ',))cos(A(t- k , x y ,,,, , z. a(tk,xm,ym,zm) m = , et le vecteur d'observation théorique vk est donné en fixant i=1 par : vk = [expP2n- -d cos(61(t k , x ,n, y', ))cos(Mtk , X , yll, z', bj] Le procédé de localisation, tel que le procédé d'évolution de localisation d'évolution temporelle des vecteurs d'onde ou VOAE (Vector Of Arrivai Evolution) initial décrit ci-dessus, repose alors sur l'exploitation des vecteurs Bk ou sur des variantes opérant sur une fonction de ces vecteurs Bk. Afin d'éliminer l'indétermination liée à la mesure des signatures spatiales, les procédés de localisation tels que VOAE mettent en oeuvre une étape de normalisation pour éliminer cette indétermination. Pour cela, on choisit une voie particulière servant de voie de référence, ce qui entraine une limitation dans les résultats obtenus, notamment, du fait de l'asymétrie introduite par ce choix. Si cette voie de référence ne capte pas de signal 3 2999724 ou un signal beaucoup plus faible que les autres voies, la localisation sera très dégradée, voire impossible. Ceci se produit, par exemple, lorsque la voie est défaillante, du fait d'une antenne et/ou une tête de réception déficiente, ou lors de l'utilisation d'un réseau antennaire à diversité de polarisation, si la voie choisie est dans une polarisation très différente de celle de l'émission reçue. L'un des problèmes non résolus par les procédés selon l'art antérieur est le cas où le signal reçu a un niveau d'amplitude très faible sur l'antenne de référence pour l'émission considérée. De manière plus générale, à l'heure actuelle il existe un besoin de disposer d'un procédé permettant de déterminer la position d'un émetteur au sol en s'affranchissant de paramètres inconnus. Dans le procédé selon l'invention, le facteur d'ambigüité va être estimé, à chaque mesure, sans faire appel à une voie de référence. L'invention concerne un procédé pour déterminer la position P d'un ou plusieurs émetteurs dans un réseau de capteurs en mouvement, le réseau comprenant au moins N aériens adaptés pour mesurer des signatures spatiales ou vecteur d'observation AK le procédé est caractérisé en ce qu'il comporte au moins les étapes suivantes : E0: Exprimer le vecteur d'observation complexe Ak de dimension N supérieur à 1, reçu sur le réseau de capteurs pour un instant k, en utilisant un coefficient scalaire complexe ak inconnu et un vecteur modèle complexe dépendant de la position P Ak Cek-'Àk(P)+ Wk avec k = 1, K les indices des instants de K mesures de vecteur d'observation, ;1k (P) un vecteur modèle complexe dépendant du paramètre inconnu P, Wk un paramètre de bruit additif, E1 : définir le critère d'estimation conjointe des valeurs des paramètres de nuisance ak et de la valeur de la position P, en utilisant la méthode de maximum de vraisemblance en présence de paramètres de nuisance et en minimisant le critère quadratique correspondant, (a p) 9 MC = arg min {c(a, P P,a 4 2999724 E2: estimer le paramètre de nuisance a pour une position donnée en minimisant le critère «a, P) , où a = [ai OEk, E4: utiliser les valeurs obtenues pour le paramètre de nuisance a dans le critère: «a, Met estimer la valeur de la position P en utilisant les K mesures Ak et les mesures attendues corrigées A'k, calculées de la façon suivante :4k (17)11 F k-1 *Ak F 1c-1 *Àk(P) ) avec Fk-1 la matrice de covariance du vecteur d'observation Ak. Selon une variante de réalisation, pour déterminer le paramètre de nuisance pour une position donnée, le procédé utilise la méthode des moindres carrés et exécute les étapes suivantes : estimer le paramètre de nuisance : arg min {c(a, P)} arg min {min (c(a,P))} P,a a en développant le critère C(OE,P), et en utilisant, par exemple, la méthode des moindres carrés, on trouve : (P) Ak (P) C (a , P) = Ak 2 FI alc 2 jlk (P)H Fkl;ik(P) 2Réelle {ai A rk 1 Ak (P) k=1 en minimisant ce critère C(a,P) par rapport au paramètre de nuisance a, pour une position P donnée, on obtient une valeur estimée pour le paramètre de nuisance otK pour un instant k donné, A , puis 4(P) 1k 1 4 a k Àk (e n-1 (P) on détermine les valeurs estimées pour chacune des valeurs k variant de 1 à K où Fk-1 est la matrice de covariance du vecteur de mesure Ak. Il est possible d'estimer la valeur de la position P en utilisant un critère des moindres carrés portant sur les mesures Ak et sur les mesures attendues corrigées A'k. 2999724 D'autres caractéristiques et avantages du dispositif 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 : - La figure 1, un exemple de localisation d'un émetteur au sol par un avion équipé d'un réseau de capteurs adaptés, - La figure 2, une illustration des angles d'incidence d'un émetteur dans un réseau de cinq antennes, et - La figure 3, un exemple d'enchaînement d'étapes mises en oeuvre par le procédé selon l'invention. Un exemple de mise en oeuvre du procédé selon l'invention est donné en relation avec les figures 1, 2 et 3. Le système permettant la mise en oeuvre du procédé comporte, par exemple, comme représenté sur la figure 1, un avion 2 équipé d'un réseau de N capteurs et d'un processeur adapté à exécuter les étapes du procédé détaillé ci-après. L'émetteur au sol, dont on souhaite déterminer la position P porte la référence 1. L'émetteur 1 est à la position (x0,yo,z0) et, le porteur 2 à l'instant tk, est à la position (xk, yk, zk) et voit l'émetteur sous l'incidence (61(tk,x0,y0,z0),A(tk,x0,yo,z0)). Les angles 6(t, xo, yo, zo) et zl(t, xo, Yo, zo) évoluent au cours du temps et dépendent de la position de l'émetteur ainsi que de la trajectoire du porteur 2. Une manière de représenter les angles est donnée à la figure 2, par rapport à un réseau de N capteurs ou antennes pouvant être fixées par le porteur. En présence de M émetteurs, l'avion reçoit à l'instant t en sortie des N capteurs du réseau, un ensemble de signatures spatiales ou vecteurs de mesures du signal reçu qui peuvent s'écrire selon le modèle suivant : Ak ak.Àk(P)+ Wk (1) où - P la position de l'émetteur, - k = 1, K sont les indices des instants des K mesures de signatures spatiales Ak est le vecteur d'observation complexe (signature spatiale), de dimension N>1 où N est le nombre d'aériens du réseau antennaire, - OEk est un coefficient scalaire complexe inconnu, dont la dépendance éventuelle avec la position P de l'émetteur 1 est inconnue a priori, 6 2999724 - k(P) est un vecteur modèle complexe, dépendant du paramètre inconnu P (position de l'émetteur) que l'on cherche à estimer. Cette étape de modélisation correspond à l'étape E0 du procédé selon l'invention. On fait l'hypothèse que le vecteur d'observation est bruité par un vecteur bruit W additif gaussien (complexe), centré, circulaire et de matrice de covariance r: W N(0,['). Wk correspond au bruit pour un instant k de mesure. Sous cette hypothèse, l'estimateur maximum de vraisemblance ou MV de la position P est égal à l'estimateur des moindres carrés ou MC, et est efficace au sens de l'atteinte de la borne de Cramer-Rao. En dehors de la position P à déterminer, les paramètres inconnus OEk, qui n'ont pas d'intérêt particulier pour l'estimation (en l'absence de modèle les reliant à P), sont des « paramètres de nuisance » empêchant une estimation directe de P. Le procédé selon l'invention va, dans une première étape El, définir le critère d'estimation conjointe des valeurs de l'ensemble des paramètres inconnus de nuisance a = [ai ak, OEK] et de la valeur de la position P en utilisant la méthode du Maximum de Vraisemblance MV. Sous les hypothèses précédentes, dans cet exemple au sens du critère de moindre carré MC, ceci revient à minimiser le critère quadratique (norme hermitienne) : Ak OC k-74k(P)) 2 (2) C (a , P) = E rÏl k=1 ce qui conduit à (a,P),,, = arg min {c(a, P)} où P,a (ce,P)e est l'estimation de a et de P selon le critère de moindre carré, 11.111- désigne la norme : X 12= XH.F.X Pour estimer la position P au sens du Maximum de Vraisemblance MV (et du moindre carré MC), le procédé va exécuter les étapes suivantes: E2, on estime le paramètre de nuisance a= [ai ak, . OEK] en minimisant le critère quadratique C(a,P) en fonction de P, pour une position donnée de P, dans cet exemple en utilisant le critère du moindre carré, 7 2999724 - E3, on injecte la valeur estimée du paramètre de nuisance a(P) (a pour une position P donnée) ainsi obtenue dans l'expression (2) de C, qui ne dépend plus alors que de P, d'où le critère C'(P), - E4, on minimise l'expression du critère C'(P) en utilisant, par exemple, la méthode des moindres carrés. Etape On estime le paramètre de nuisance pour une position donnée arg min {c(a,P)} = arg min {min (c(a, P))} (4) P,a a En développant le critère C(oc,P), et en utilisant, par exemple, la méthode des moindres carrés, on trouve : C(ot,P)- 2 k 2 :4-k(P)H FA (P) (P)-2Réelle{a AkH r k (P)}) (5) k=1 en minimisant ce critère C(oc,P), par rapport à a, pour une position P donnée, on obtient une valeur estimée pour le paramètre de nuisance OEK pour un instant k donné, A (6) a k (P) ';4k (P) H Fk 1 .Ak ;1k (en-151k (P) On détermine les valeurs estimées pour chacune des valeurs k variant de 1 à K, Fk-1 est la matrice de covariance du vecteur de mesure Ak Etape E3 En injectant les valeurs des paramètres de nuisance dans le critère de l'équation (5) on obtient un critère dépendant uniquement de la position P : 2 c(p)= k=1 A k ( re, \H r-, -1 A (7) k -zik k (P)II F -;1k (P) On somme sur les K paramètres de nuisance désigne la norme : r2 = .X 8 2999724 Etape E4 Le procédé consiste alors à appliquer, par exemple, un critère des moindres carrés portant sur les mesures Ak et sur les mesures attendues corrigées A'k, calculées de la façon suivante : r :4k (P)H F 1:1 *Ak 14ke (P) = Àk (p)H F k-1 51k (P) On en déduit une valeur estimée pour la position P. Dans cet exemple, elle correspond à une estimation au sens des moindres carrés F'mc, = arg min C' (P) (9) avec C' (P) = Ak-Àlk (P)) rk- I 2 k=1 (10) A partir des valeurs de position P obtenues, il est ensuite possible d'appliquer des techniques connues de l'Homme du métier pour localiser les émetteurs une fois la position connue. Le fait d'estimer le facteur d'ambigüité au lieu de l'éliminer en choisissant une voie de référence, permet de disposer d'un procédé symétrique vis-à-vis des différentes voies de réception du système et offre au procédé une indépendance par rapport au niveau de la qualité de réception du signal reçu sur une voie particulière. Le procédé est robuste aux dysfonctionnements d'antenne et/ou de tête de réception radio fréquence. Il offre de plus la possibilité de mettre en oeuvre une localisation de type VOAE sur des réseaux à diversité de polarisation. iik(P) (8) 9 REVENDICATIONS1 - Procédé pour déterminer la position P d'un ou plusieurs émetteurs (1) dans un réseau de capteurs (2) en mouvement comprenant au moins N aériens adaptés pour mesurer des signatures spatiales ou vecteur d'observation Ak caractérisé en ce qu'il comporte au moins les étapes suivantes : E0: Exprimer le vecteur d'observation complexe Ak de dimension N supérieur à 1, reçu sur le réseau de capteurs pour un instant k en utilisant un coefficient scalaire complexe ak inconnu et un vecteur modèle complexe dépendant de la position P Ak ak.Àk(P)+Wk avec k = 1, K les indices des instants des K mesures de vecteur d'observation, Ak(P)unvecteur modèle complexe dépendant du paramètre inconnu P, Wk un paramètre de bruit additif, El: définir le critère d'estimation conjointe des valeurs des paramètres de nuisance ak et de la valeur de la position P, en utilisant la méthode de maximum de vraisemblance en présence de paramètres de nuisance et en minimisant le critère quadratique correspondant, G p) \ 9 MC = argminMa, P)} P,a E2: estimer le paramètre de nuisance a pour une position donnée en minimisant le critère C(a, P), où a = E4: utiliser les valeurs obtenues pour le paramètre de nuisance a dans le critère P)et estimer la valeur de la position P en utilisant les K mesures Ak et les mesures attendues corrigées A'k, calculées de la façon suivante : (A(P)H rk 1 *Ak Fk-1'71k(P) ) avec Fk-1 la matrice de covariance du vecteur d'observation Ak. 2 - Procédé selon la revendication 1 caractérisé en ce que pour déterminer le paramètre de nuisance pour une position donnée, on utilise la méthode des moindres carrés et on exécute les étapes suivantes : On estime le paramètre de nuisance k' (P) Ak(P) 10 2999724 arg min {C (a, P)} = arg min {min («a, pb} En développant le critère C(OE,P), et en utilisant la méthode des moindres carrés, on trouve : K f rki 2 + lak 7ik (p)H n-1-À C(a,P)= Ak k( P) 2Réelle{akAkH Fic Àk (P)}) k=1 en minimisant ce critère C(P) par rapport au paramètre de nuisance a, pour une position P donnée, on obtient une valeur estimée pour le paramètre de nuisance ak pour un instant k donné, ;4; (P)II F k 1 *Ak 1( k (p) ;4k (P)HFk-1 (P) on détermine les valeurs estimées pour chacune des valeurs k variant de 1 à K où rk-1 est la matrice de covariance du vecteur de mesure Ak. 3 - Procédé selon l'une des revendications 1 ou 2 caractérisé en ce que l'on estime la valeur de la position P en utilisant un critère des moindres carrés portant sur les mesures Ak et sur les mesures attendues corrigées kk. , puis