FR 3003960 A1 20141003 FR 1300712 A 20130327 Procédé de détection de signaux destinés à leurrer un récepteur de signaux d'un système de navigation par satellites et récepteur associé La présente invention porte sur un procédé et système de détection de signaux leurrés d'un système de navigation par satellites de signaux reçus par un récepteur muni d'un réseau d'antennes à rayonnement contrôlé. La radionavigation par satellites permet d'obtenir la position d'un récepteur par une résolution voisine de la triangulation, utilisant les pseudo-distances mesurées à partir des signaux envoyés par les satellites. Les signaux civils émis par les satellites sont prévisibles donc reproductibles. Cela constitue une menace pour les utilisateurs car il est de ce fait possible de générer des signaux pirates ou leurrés ressemblant aux signaux des satellites. Lorsqu'un récepteur s'accroche dessus il mesurerait alors une position erronée, ce qui peut avoir de graves conséquences s'il s'agit d'un avion. Il est connu plusieurs procédés pour détecter une tentative de leurre de signaux satellitaires ou leurrage. Un procédé consiste à détecter une similitude sur les profils de rapports signaux à bruit mesurés sur les signaux satellites. Ce procédé ne fonctionne que dans des environnements avec des obstacles générant des atténuations, masquages et multi-trajets intempestifs, comme en ville. En outre, si les signaux de leurrage sont modulés suivant des profils temporels différents des puissances de chaque signal satellite généré, alors la détection ne verra rien. Un autre procédé consiste à utiliser la différence d'effet Doppler en fonction des directions d'arrivée des signaux satellitaires. Si les effets Doppler ne sont pas cohérents avec la vitesse estimée par le récepteur alors il y détection de leurrage. Cependant il est facile pour le dispositif de leurrage de générer des signaux donnant des mesures d'effet Doppler cohérents dans le récepteur. Ce dernier mesure alors une position et une vitesse erronées sans voir d'incohérence sur les mesures. L'utilisation d'une antenne à plusieurs capteurs permet de mesurer les directions d'arrivées des signaux grâce aux déphasages et de vérifier qu'elles sont cohérentes par apport aux positions des satellites et à l'attitude de l'antenne supposée connue. Selon l'état de l'art ce dernier procédé utilise un traitement du signal séparé entre les capteurs et ne bénéficie pas de la redondance des signaux pour améliorer la robustesse de la démodulation. Un but de l'invention est de pallier cet inconvénient. Il est proposé, selon un aspect de l'invention, un procédé de détection de signaux destinés à leurrer un récepteur de signaux d'un système de navigation par satellites, le récepteur étant muni d'une antenne à capteurs à rayonnement contrôlé, comprenant les étapes suivantes consistant à, pour chaque satellite : - calculer la position du satellite dans un premier repère TGL centré sur le récepteur à partir de la position du récepteur et de la position du satellite; - déterminer l'attitude de l'antenne dans le premier repère; - calculer des premiers angles de site et d'azimut du satellite dans un deuxième repère lié à l'antenne, à partir de la position du satellite dans le premier repère et de l'attitude de l'antenne dans le premier repère; - déterminer des deuxièmes angles de site et d'azimut du satellite dans le deuxième repère, par recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne; - calculer la valeur d'une fonction de vraisemblance entre les premiers angles de site et d'azimut du satellite et les deuxièmes angles de site et d'azimut du satellite; et - détecter un risque de réception de signaux leurrés lorsque ladite valeur de la fonction de vraisemblance est inférieure à un seuil. Un tel procédé permet donc de nettement améliorer la détection de signaux émis dans le but de leurrer le récepteur de signaux d'un système de navigation par satellites, et ainsi éviter de graves conséquences à l'utilisateur du récepteur, notamment dans le cas d'un aéronef. Dans un mode de mise en oeuvre, ladite recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne est effectuée à partir d'un discriminateur d'angle de site et d'un discriminateur d'angle d'azimut. Ainsi, le procédé permet d'améliorer la robustesse au brouillage grâce à la redondance des capteurs de l'antenne. Selon un mode de mise en oeuvre, l'étape de détermination de l'attitude de l'antenne dans le premier repère utilise des données fournies par une référence inertielle. L'attitude est ainsi non leurrable. Dans un mode de mise en oeuvre, l'étape de détermination de l'attitude de l'antenne dans le premier repère utilise ladite détermination des deuxièmes angles de site et d'azimut des satellites dans le deuxième repère et procède par une recherche itérative des angles d'attitude de l'antenne donnant le minimum de la fonction de vraisemblance. Ainsi, il n'y a pas besoin de référence inertielle. Selon un mode de mise en oeuvre, l'étape de détermination de l'attitude de l'antenne dans le premier repère procède par une recherche itérative des angles d'attitude de l'antenne donnant le maximum d'énergie après les sommes pondérées en sortie des corrélateurs des canaux de traitement respectivement associés aux satellites. Ainsi, on améliore la robustesse au brouillage. Dans un mode de mise en oeuvre, ladite maximisation utilise une méthode de résolution par les moindres carrés. Cette mise en oeuvre est simple. Selon un mode de mise en oeuvre, ladite maximisation utilise une méthode de résolution par filtrage de Kalman. Une telle mise en oeuvre permet d'augmenter la précision et la robustesse au brouillage. Dans un mode de mise en oeuvre, on utilise des discriminateurs à base de cosinus directeurs représentatifs de la position des satellites dans le deuxième repère lié à l'antenne. Dans un mode de mise en oeuvre, on utilise une méthode à modèle d'observation de la relation entre les cosinus directeurs des directions des satellites et les angles d'Euler des directions des satellites.35 Il est également proposé, selon un autre aspect de l'invention, un récepteur de signaux d'un système de navigation par satellites, muni d'une antenne à capteurs à rayonnement contrôlé, comprenant, pour chaque satellite : - des moyens de calcul de la position du satellite dans un premier repère TGL centré sur le récepteur à partir de la position du récepteur et de la position du satellite; - des moyens de détermination de l'attitude de l'antenne dans le premier repère - des moyens de détermination des premiers angles de site et d'azimut du satellite dans un deuxième repère lié à l'antenne, à partir de la position du satellite dans le premier repère et de l'attitude de l'antenne dans le premier repère; - des moyens de détermination des deuxièmes angles de site et d'azimut du satellite dans le deuxième repère, par recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne; - des moyens de calcul de la valeur d'une fonction de vraisemblance entre les premiers angles de site et d'azimut du satellite et les deuxièmes angles de site et d'azimut du satellite; et - des moyens de détection d'un risque de réception de signaux leurrés lorsque ladite valeur de la fonction de vraisemblance est inférieure à un seuil. L'invention sera mieux comprise à l'étude de quelques modes de réalisation décrits à titre d'exemples nullement limitatifs et illustrés par les dessins annexés sur lesquels : - la figure 1 illustre schématiquement un procédé de détection de leurrage, selon un aspect de l'invention; - la figure 2 illustre les angles d'Euler; - la figure 3 illustre une antenne à quatre capteurs; - la figure 4 illustre une démodulation avec boucles de poursuite; - la figure 5 illustre la détermination des angles de site et d'azimut, selon un aspect de l'invention; - la figure 6 illustre un contrôle de cohérence sur les directions des satellites, selon un aspect de l'invention; - les figures 7, 8, 9 et 10 illustrent la détermination de l'attitude de l'antenne à partir des directions satellites mesurées; - la figure 11 illustre un synoptique du procédé selon un aspect de l'invention; et - la figure 12 illustre un récepteur selon un aspect de l'invention. Sur la figure 1, est illustré un procédé de détection de signaux destinés à leurrer un récepteur de signaux d'un système de navigation par satellites, selon un aspect de l'invention, le récepteur étant muni d'une antenne à capteurs à rayonnement contrôlé. Le procédé comprend, pour chaque satellite du système de navigation par satellite les étapes qui suivent : On calcule, dans une étape 1, la position du satellite dans un premier repère TGL (Trièdre Géographique Local, repère Nord, Est, Vertical centré sur la position du récepteur) centré sur le récepteur à partir de la position du récepteur et de la position du satellite. Les positions du récepteur et du satellite sont connues grâce aux éphémérides ou aux almanachs, contenus dans les messages de navigation transmis par les signaux émis par les satellites. On détermine, dans une étape 2, l'attitude, ou orientation définie par trois angles, de l'antenne dans le premier repère. Cela peut être fait de plusieurs manières, comme décrit ultérieurement dans la présente demande de brevet. On détermine, lors d'une étape 3, des premiers angles de site et d'azimut du satellite dans un deuxième repère lié à l'antenne, à partir de la position du satellite dans le premier repère et de l'attitude de l'antenne dans le premier repère. On détermine, dans un étape 4, des deuxièmes angles de site et d'azimut du satellite dans le deuxième repère, par recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne. Cette étape peut être réaliser de différentes manières, comme décrit ultérieurement dans la présente demande. On calcule, dans un étape 5, la valeur d'une fonction de vraisemblance entre les premiers angles de site et d'azimut du satellite et les deuxièmes angles de site et d'azimut du satellite. Enfin, on détecte, lors d'une étape 6, un risque de réception de signaux leurrés lorsque la valeur de la fonction de vraisemblance est inférieure à un seuil. Le vecteur directeur entre l'antenne du récepteur et le satellite d'indice n du système de navigation par satellites, dans le premier repère TGL, est égal par définition à : def X n TGL lp 14n TGL = Y n TGL Z n TGL Le dénominateur p représente la distance séparant l'antenne du récepteur et le satellite n, indépendante du repère : def = x n2 + yn2 z n2 = n TGL2 ± Y n TGL 2 Z n TGL2 Dans l'ensemble de la présente demande, la présence d'un indice def au-dessus d'un signe d'égalité dans une expression mathématique signifie une égalité par définition et non par conséquence. Les positions des satellites (Xn TGL TGLI Zn TGL) en repère TGL sont calculées grâce aux éphémérides. L'étape 2 de détermination de l'attitude de l'antenne utilise traditionnellement une référence inertielle d'attitude et de cap, d'acronyme AHRS pour "Attitude and Heading Reference System" en langue anglaise, hybridée avec d'autres capteurs (anémométrie, GPS, radar Doppler,...) ou par une référence inertielle autonome (IRS - lnertial Reference System). Il est aussi possible de déterminer l'attitude de l'antenne du récepteur grâce aux signaux émis par les satellites eux-mêmes, qui arrivent avec des déphasages différents sur les capteurs de l'antenne, en fonction de la direction incidente. Plusieurs variantes sont décrites dans la suite de la description. L'étape 3 de calcul des directions des satellites dans le deuxième repère lié à l'antenne dit repère antenne ou repère capteur, peut être réalisée comme suit : On représente, comme sur la figure 2, l'attitude de l'antenne par rapport au premier repère TGL grâce aux angles d'Euler : Cap (v), Assiette (0) et Inclinaison (9). Ces angles permettent de construire la matrice de changement de repère M(yr, 0, (p) du premier repère TGL au deuxième repère lié à l'antenne ou aux capteurs d'antenne (idem) : def --1 0 0 - 0 - sin 0 cos yi - sin sin cos 1,11 0- M (Vi O, ço) = 0 cos ço sin ÇO - sin ço cos yo cos 0 1 0 0 0 0 0 0 0 cos 1 sin 0 - X capteur - X TGL Ycapteur =-M-((e,t9e g") - YTGL _Z capteur _ Z TGL Le vecteur directeur entre l'antenne du récepteur et le satellite d'indice n du système de navigation par satellites, dans le deuxième repère lié à l'antenne, est égal par définition à: def X n capteur 'p 171> n capteur = Y n capteur Z n capteur _ Le dénominateur p représente la distance séparant l'antenne du récepteur et le satellite n, indépendante du repère : 2 Vxn 2 ± yn 2 + zn 2 = .\iy 2 v 2 _L 7 n capteur n capteur n capteur On peut ainsi exprimer le vecteur directeur du satellite d'indice n en fonctions des angles d'Euler : !Ti capteurs 91(9) = 1111 (g (°) - il. n TGL def Un n capteu' (V , O, (0) = 25 On calcule alors, dans l'étape 2, dans le deuxième repère lié à l'antenne, les angles de site yn et d'azimut an à partir des angles d'attitude de l'antenne (ou des capteurs) dans le deuxième repère : def cos(yn).cos(crn ) def ux,7 çor n capteurs (VI, 6), (t)) sin (y, ).cos(an ) = uy n (y t, 0 , ça) uz (v, e, ço) sin(o) _ L'angle de site en repère antenne ou capteur s'écrit : an(v,e,cp ) = sin-1( uzn(v,0,(p ) ) L'angle d'azimut en repère antenne ou capteur s'écrit : yn(1g,8,(p ) = tan-1( uxn(v,e,q) ) uYn(Ilf,0,(1) ) ) Le vecteur directeur calculé à partir des angles de site et d'azimut variables s'écrit alors : - cos(y). cos (CF) def- - ux(cF,r) def capteurs (a, = sin (y). cos() = uy(o-, sin(a)uz(cr, y) On peut également prendre comme paramètres de direction les cosinus directeurs sur les axes X et Y du repère antenne ou repère capteur. Le cosinus directeur sur l'axe Z se déduit grâce à la norme égale à 1. Pour éviter toute ambiguïté sur le signe de Z il faut supposer que les satellites se trouvent toujours au-dessus du plan de l'antenne, donc ignorer les satellites en-dessous. Cela est légitime dans la mesure où les satellites bas en repère capteur ne sont pas exploitables pour la détermination des directions d'arrivée à cause de gains d'antenne inadaptés à un angle de site bas. def fi .'J1_a2- ,62 def UX(a, fi) el> capteurs (a, 16) = = uy(a, ,6) _uz (a, ,6) _ L'angle de site en repère antenne ou capteur s'écrit : a(ocim = tan-i (Vcr2 ± /32 , V]. _ cx2 4., /32 ) L'angle d'azimut en repère capteur s'écrit : y(a,f3) = tan-1( a , p ) _ de! X m capteur En notant Fm capteurs = Y ni capteur le vecteur position du capteur m en repère Z m capteur _ capteurs, connu par exemple grâce à des mesures de centre de phase en usine, dans un exemple à quatre capteurs d'antenne, on a : il capteurs - r -i .; 0 r F3 capteurs = _ 7. _ 0 2 capteurs = 0 -r 4 capteurs = 0 0 0 -r 0 _ 0 comme illustré sur la figure 3. Les déphasages supposés des signaux reçus entre les capteurs d'antenne sont calculés à partir des angles de site et d'azimut variables. La différence de chemin optique entre le centre du deuxième repère antenne et le capteur d'indice m pour le signal reçu du satellite n s'exprime comme suit : def e m,n (g f 9 9, g)) = u,, ri capteurs (g f 56 e 1 g)) - Fni capteurs (en m) Le déphasage correspondant s'exprime comme suit : de! 8,,,, 0 I f 5 t 9 50 = n capteurs ( e 3 t 9, g') - 7'm capteurs - (2z/2) (en rad) X, représentant la longueur d'onde de la porteuse On peut alors calculer le déphasage pour des angles de site et d'azimut variables cr, y. On calcule alors la différence de chemin optique : de! m (a, Y) = 'il-capteurs (cf, Y) - in, capteurs (en m) et le déphasage correspondant : de! m (Cr = 11 capteurs (Ci. e 7) n capteurs - (27r I 2) (en rad) L'utilisation de l'indice n car la position du satellite n'intervient plus à partir du moment où on a les angles de site et d'azimut en deuxième repère capteur. On calcule le déphasage pour des cosinus directeurs variables. On commence par le calcul de différence de chemin optique : def em (a, /3) = iicapteurs (a, /3) - Fm capteurs (en m) et le déphasage correspondant : de! 8 m (C fi) = 17 capteurs (a, /3) -m capteurs - (271- I 2) (en rad) L'utilisation de l'indice n disparaît car la position du satellite n'intervient plus à partir du moment où on utilise les angles de site et d'azimut variables. La démodulation des signaux, numérisés à la fréquence d'échantillonnage Fe, est réalisée de manière matérielle, dans un composant numérique "câblé", ASIC ou FPGA. Elle demande une fréquence de calcul très élevée, difficile à réaliser de manière logicielle. La démodulation est propre à chaque signal satellite, caractérisé par son code d'étalement et son Doppler. Dans la présente demande, on démodule tous les signaux reçus par les capteurs en parallèle dans chaque canal satellite matérielle. La démodulation procède par multiplication du signal reçu par une porteuse locale complexe et par un code local réel puis intégration du produit par intervalles de temps consécutifs (c'est la corrélation). On a en fait deux voies de corrélation par signal reçu : une voie ponctuelle avec un code local en phase avec le code reçu et une voie delta avec un code "delta" ou différence entre code avancé et code retardé, selon l'état de l'art. La figure 4 représente un exemple de démodulation avec les boucles de poursuite. Les sommes pondérées des signaux démodulés des différents capteurs sont réalisées en logiciel à basse fréquence (typiquement Fw = 50 Hz). Les coefficients de pondération complexes ont pour rôle de mettre en phase les signaux des capteurs en sortie des voies de corrélation dans le but d'avoir le maximum d'amplitude après la somme. On utilise les directions des satellites calculées à partir des angles de site et d'azimut pour déterminer les coefficients de pondération. Les résultats en sortie des corrélateurs avant pondération s'écrivent : _Z pm Mx! ZAM Mx! zpl def zp2 Z P capteurs = zp3 def Z A capteurs = zm z,2 zà3 La pondération sur la voie ponctuelle s'écrit alors : def Zp = CnT(VM,(p ) - ZP capteurs (produit scalaire) def Zp = C1,n(Veee )- ZP1 C2,n(VMe ). ZP2 CM,n(11110,(1) )2PM La pondération sur la voie Delta (discriminateur de code) s'écrit alors : def = Cn-r(V,O-() ) - Z4 capteurs (produit scalaire) def Zn, = Ci,n(WM,Y ). Zel C2,(V,0,(1) ). ...+ Cm,n(V,0,(1) )2Am Les coefficients de pondération complexes s'écrivent alors : def C n (V , 0,ço)= Mx! def C m,n(g,0,y))= Acapteur ',(u n,yn).exp i( (5,n, (yt, 0, ço) + Bcapteu,. ,n (0 n,yn)) Acapteur m (47 n rn) représentant l'amplitude relative due au capteur m dans la direction du satellite n calculée à partir de la connaissance à priori des diagrammes d'antenne, y compris les atténuations dues aux voies RF analogiques. Bcapteur m(er n rt1 ) représentant le déphasage dû au capteur m dans la direction du satellite n calculé à partir de la connaissance à priori des diagrammes d'antenne, y compris les déphasages dus aux voies RF analogiques. Il est possible d'identifier les biais des capteurs en usine. Mais cela augmente le coût de production de l'équipement (moyens de calibration et temps passé) et ne permet pas de compenser le vieillissement du capteur. Il est aussi possible d'identifier les déphasages des capteurs grâce aux signaux des satellites par un traitement spécifique (cf document FR 2972809) Les boucles de poursuite (code et porteuse) sont réalisées de manière logicielle car elles fonctionnent à basse fréquence (typiquement de l'ordre Fw = 50 Hz), elles utilisent des fonctions mathématiques plus complexes, et présentent différents modes (recherche, transition, poursuite, poursuite dégradées, ...) On utilise un discriminateur de phase : Dp (k) = argument( Zp ) Dy (k) = arc tangente( lp , Qp ) Zp = lp + j.Qp I représentant la partie réelle, et Q représentant la partie imaginaire Dip(k)=O Zp = pele p représentant le module, et 0 représentant l'argument K représentant l'indice de la période de recalage (typiquement des périodes 5 de 20 ms pour Fw = 50 Hz) On utilise également un discriminateur de code : Dcode (k) = partie réelle( Z,. Zp. ) Zip* = conjugué( Zp* ) = lp - j-Qp 10 Dcode (k) = lp . IA + j.( lp . là ) On utilise un correcteur de boucle de phase (PLL pour "Phase Lock Loop" en langue anglaise) : 15 Aporteuse(k) = Aporteuse(k-1 ) + T.Kporteuse 3. D(p(k) . Iporteuse/(27t) (en m/s2) (en m/s) Vporteuse(k) = Vporteuse(k-1) + TA () .p porteuse.k + T K orteuse 2. D(P(k) - ?porteuse/(21t) Cporteuse(k) = Vporteuse(k) + Kporteuse 1. D(P (k) - Xporteuse/(27c) (en m/s) 20 On utilise un correcteur de boucle de code (DLL pour "Delay Lock Loop" en langue anglaise) : Vcode(k) = Vcode (k-1 ) TAcode (k) + T. K-ode 2. -code(k) - L chip (en m/s) Ccode (k) = Vcode (k) kcode 1. Dcode(k) - Lchip (en m/s) 25 Concernant les oscillateurs à commande numérique ou NCO pour " numerically controlled oscillators" en langue anglaise : CNCO porteuse (k) = Cporteuse(k) I Iporteuse (en Hz) 30 CNco code (k) = ( Cporteuse(k) Ccode(k) ) Lchip (en Hz) On détermine alors, dans l'étape 4, les deuxièmes angles de site et d'azimut comme suit. Il s'agit de déterminer les directions d'arrivée des signaux poursuivis. On cherche, satellite par satellite, les deuxièmes angles de site et d'azimut (angles d'arrivée dans le deuxième repère antenne ou capteurs) qui maximisent la puissance après la somme pondérée. Ce traitement est effectué de manière logicielle, en même temps que les boucles de poursuite, comme illustré sur la figure 5. Les coefficients de pondération complexes à optimiser s'écrivent selon le critère qui suit : de! C(o-, y) = Cm (o-, y) M x 1 de! C (u, y) = exp i( ,n(o , y) + B capteur ni(o' , y) ) On cherche les valeurs des angles de site a et d'azimut y qui maximisent le critère C(cr, y) : de! j (erY 11 Z(I5 ) 112 de! -p capteurs IF(a ) Il CT(a ,y ) . 7 112 de! 111(a ,y ) = [ CT(a ,y ) . Zpcapteurs ] [ CT(a ,y ) . ZP capteurs ] Cette maximisation peut être faite de plusieurs manières différentes. Il est possible de faire cette maximisation par une recherche de maximum dans l'ensemble des valeurs possibles de a et a comprises respectivement dans [0,7c] et [0,2n], avec un pas d'échantillonnage suffisamment fin pour la précision demandée. Il est également possible d'effectuer une recherche de maximum par la méthode itérative du gradient, en s'approchant du maximum à chaque itération. [Gr i+11.[cr'1+ .Ôracel (cri,yi) 7i-Fi qi def ÔraclIP (a-, y) = i représente l'indice de l'itération [ acT Re [ (o, y).Z p capteurs] [CT (o-, y).Zp capteurs]* Do- [acr Re [ (o-, y).Zp capteurs]. [(CT (a-, y).Z p capteurs]* Dy Ôrace- (o y) = [Re( Zs (cr, y) Ôrad'Y (a-, y) = Re( Z (o-, Z (c r, y)* Z (o-, y)* Avec : def Z (0-, y) = y).z p capteurs def S = ao- y).Z p capteurs def Z (o-, y) = (0", y).Zp capteurs a Y Notation : Re(z) représente la partie réelle de z On construit des discriminateurs pour estimer l'erreur de site et d'azimut et converger plus rapidement vers les valeurs vraies (ob,y0) donnant le maximum de puissance après pondération. On utilise des discriminateurs de site et d'azimut : def Ds(a,y) = Re[ Zs(G,7)-Z(G,7)*1/11 Z 112I Ps = Re( Zs I Z*)I Ps def DA(G,y) = Re[ ZA(G,Y) - Z(,7)*I /II Z 112 / PA = Re( ZA Z*) PA On utilise la somme pondérée "ponctuelle" : def z (0-, y) = CT vl , 'Z P capteurs Et les sommes pondérées avec des pondérations complexes "delta" : def T def acT zS (a-, y) = C (cr -kr (0-, r) - P capteurs def Z A(0", y) = C AT (a , y).Z P capteurs def acT c AT (o-, y) = Pour calculer la sensibilité des discriminateurs, on pose : def n(cY,Y) = CT (e Ar(0-0, y.) def acr a As(a,y) =-a cs (a r).c(cro , ro) = ( )(0c,13) def aCT AA(0,11) = -acr (0-,A.c(0-o,yo) = -aay( ) (00) Le signal reçu s'exprime: Z P capteurs = vrai 5 yvrai).T + bruit P représentant la puissance du signal reçue, et T représentant le temps d'intégration cohérente -, Le vecteur C(vrai Ivraiv 1cs est représentatif des amplitudes complexes reçues sur chaque capteur. Les angles (Gvrai,ïvrai) correspondent à la direction vraie donc au maximum de 111(cy,7 ). D'où : Z (o-, y) = CT (o-, y) .0 (a- 'a' yvrai).T acT Z s (o-, y) = (o-, y) .0 (a- vra, , y'ai).T acs aCT Z A(CF, y) = (o-, y) .0 (a- 'ai , yvra,).T a y Dans l'hypothèse où on connaît la v 1 en posant (ao, yo) = v n , -vrai, ivraii, -vrai, ivrai, -n a : Z (o-, y) = II(a- , y).T .,17) Z s (o-, y) = às (cr, y).T Ji; Z A(0-, y) =(cl, y).T D'où : Ds(a,y) = Re[ Asn] /II n112/ Ps DA(a,y) = Re[ / Il n 112/ p. Pour normaliser la sensibilité des discriminateurs à 1 on pose : def a ps = Re[ As.n ]/ lin 112 )(y) def PA = Re[ Il n112) (a0,70) ay Ainsi : Ds(a,y) (a- 'Geai) + bruit DA(a,y) (y- "(vrai) + bruit 20 Les discriminateurs permettent bien d'estimer les erreurs de site et d'azimut. Pour calculer ces pondérations "delta" on peut utiliser la méthode des différences finies : def ac c( + do-, - C(a, y) 25 C( cl, -a 0,(.0 cy do- def ac ay dy Pour da et da petits, le résultat est très proche de la dérivée. Cette méthode qui évite de recoder un calcul de dérivée dans le récepteur. 30 'vrai, Comme on ne connaît pas les vraies valeurs la de (cs,y) on prend à la place les valeurs estimée à l'étape précédente (abyi) 15 La méthode du gradient converge vers un maximum local qui est le maximum absolu si les valeurs initiales ao et yo en sont proches. Pour cela on peut commencer par une recherche globale avec un pas d'échantillonnage grossier, par exemple un pas de l'ordre de 30°. Du choix du gain de recalage [..ti dépend la vitesse de convergence. On peut prendre par exemple un gain constant égal à 1. Le critère de convergence peut être la norme de l'incrément pi.ÔraclED(crs, ) Une fois que l'algorithme, appliqué au même vecteur Zp capteur(k) de la période k, a convergé, on peut utiliser les valeurs de a et y ainsi obtenue pour réinitialiser la recherche à l'étape k+1 suivante. Si l'antenne est sur un porteur avec une dynamique d'attitude normale les nouvelles valeurs de site et d'azimut ne doivent pas être trop éloignées des précédentes, ce qui autorise une itération par période k: [Crk+11=[Crk1±[PSk-DS(Crk 5 rk)1 Îlc-F1 rk IlAkDA(ek,rk) On peut prendre par exemple comme valeurs de gain de recalage Ls=11A = 1 / ('r :constante de temps). On peut aussi utiliser un filtre de Kalman pour optimiser les gains de recalage et 1-tA Il est possible de faire la même recherche en utilisant comme paramètres les cosinus directeurs a et [3. On contrôle la cohérence sur les directions des satellites. Une fois que le récepteur a déterminé les directions d'arrivée des signaux des satellite en deuxième repère antenne, représentés par les angles de site et d'azimut pour chaque satellite, il peut vérifier si ces directions son cohérentes des angles d'attitude de l'antenne déterminés par ailleurs. On peut utiliser, à ces fins, comme critère la somme des carrés des écarts entre les angles de site an et d'azimut yn déterminés sur le signal démodulé en sortie des corrélateurs et ceux calculés à partir des angles d'attitude et de la direction des satellites. Néanmoins du fait de discontinuités sur les angles, on préfère effectuer la comparaison directement sur les vecteurs directeurs : def ' c1) = E',__,,. (v e (0)- 0- '2 +0 rn(i,19,(0)-rk,n11 Cela est illustré sur la figure 6. Si le critère dépasse un seuil prédéterminé alors on détecte un leurrage et le récepteur informe l'utilisateur du risque encouru, au moyen d'une alarme, par exemple visuelle ou sonore. Le récepteur peut aussi décider de relancer une acquisition en espérant retrouver le signal non leurré, jusqu'à obtention de mesures cohérentes. Ensuite, dans les étapes 4 et 5, on fait un contrôle de cohérence, comme illustré sur la figure 7. Une fois que le récepteur a déterminé, dans l'étape 4, les directions d'arrivée des signaux satellite en repère capteur, représentés par les angles de site et d'azimut pour chaque satellite, il peut vérifier si ces directions son cohérentes des angles d'attitude de l'antenne déterminé par ailleurs dans l'étape 3. On peut utiliser pour cela, une fonction de vraisemblance entre les premiers et deuxièmes angles de site et d'azimut, dans l'étape 5. Par exemple, comme critère, on peut prendre la somme des carré des écarts entre les angles de site an et d'azimut oc,-, déterminés sur le signal démodulé en sortie des corrélateurs et ceux calculés à partir des angles d'attitude et de la direction des satellites. Néanmoins du fait de discontinuité sur les angles, on peut effectuer la comparaison directement sur les vecteurs directeurs : de! 2 e = lin capteurs (Vi g)) - capteurs (Cr k ,te 5 k ,')11 Si le critère dépasse un seuil prédéterminé, alors on détecte un leurrage et le récepteur informe l'utilisateur du risque encouru, au moye d'une alarme, par exemple visuelle ou sonore. Le récepteur peut aussi décider de relancer une acquisition en espérant retrouver le signal non leurré, jusqu'à obtention de mesures cohérentes. Concernant l'étape 2 de détermination de l'attitude de l'antenne dans le premier repère TGL, celle-ci peut être fournie par une centrale inertielle ou référence d'attitude. Lorsque le récepteur ne dispose pas d'une référence d'attitude, il est possible de déterminer les angles d'attitude à partir des directions satellites en repère capteur déterminés à partir des signaux reçus, dans l'étape 4. Dans ce cas, en situation de leurrage, le critère de cohérence sera diminué par rapport au cas où on dispose d'une référence externe, mais du fait de la multiplicité des satellites on détectera encore l'incohérence. Dans l'étape 4, on détermine l'attitude de l'antenne comme suit. La relation entre les angles d'Euler (v,0,y) et les directions d'arrivée n'est pas linéaire. On ne peut pas utiliser un algorithme de résolution directe, par exemple la méthode des moindres carrés. Il est possible d'effectuer une recherche dans le domaine des angles d'Euler : pour chaque hypothèse (ie,0,(p) considérée, on calcule les directions supposées d'arrivée des signaux en deuxième repère ou repère capteur, puis le critère ce cohérence. On retient les valeurs des angles donnant le critère minimal. def '2 (1)(V) 199 = capteurs 9, (°) capteurs ,11 rk,n)il Néanmoins, une recherche complète du domaine dans lequel les angles d'Euler prennent leurs valeurs serait très longue et coûteuse en calcul : Le cap N prend ses valeurs dans [0,2n] ou [0°,3601 L'assiette 0 prend ses valeurs dans [-rc,n] ou [-180°,1801 L'inclinaison cp prend ses valeurs dans [-n,n] ou [-180°,1801 Si on échantillonne de 10 degrés en 10° cela fait 36 x 36 x 36 = 46656 hypothèses à traiter. On préfère faire une recherche itérative en se dirigeant vers le minimum à chaque itération. On peut utilise la méthode de Newton dans un espace à plusieurs dimensions : Ôradc1)(K , 0, , çoi) -411)W9ei,q)i)*1 2 Ôrad(1)(,ço,)11 (V 6 (0) Ôrado(v,0,ço) = def ae -(v,e,V) aço avec : a(13 Lorsque le minimum du critère de vraisemblance est nul, la méthode converge très rapidement vers ce minimum, comme illustré sur la figure 8. C'est le cas idéal dans lequel il existe une solution exacte (résidu nuls), c'est à dire le cas dans lequel les d m,n mesuré sont sans bruit de mesure et où les biais d'antennes 50 capteur m sont parfaitement tabulés. Dans la réalité, il est peu probable que la somme des carrés des résidus soit nulle. Il faut donc s'attendre à ce que la méthode ne converge plus difficilement, comme illustré sur la figure 9. Pour éviter la divergence, on peut avantageusement utiliser une méthode moins rapide mais plus robuste : e, - a . Ôrad(13(v'01,q),) _so, comme illustré sur la figure 10. Pour calculer le gradient on peut utiliser la méthode des différences finies comme illustré sur la figure 11 en utilisant : de! (13(yi + dyi3O,ço)-(13(e,0,ço) Ôrad(13(v,(9,(0)= liM dyi d .-40 (13(v,0 + dO, y)) - (13(v,0,ço) dO (13(w,0, q)+ dço)-cte(yi3O,ço) dço Ôradolè(yf,0,w)L- 22 pour dy, de et dy petits ENyf + d y f,0,0- (1)(v, 0, y)) dyf 1:13(v , 0 + d ,y)) - (DM 0 , ça) de 1:13(v , 0 , + d ça) - (I)(v ,ço) dço En variante, l'étape 2 de détermination de l'attitude de l'antenne dans le premier repère ou repère antenne ou capteurs. Au lieu de déterminer les angles d'attitude de l'antenne à partir des directions d'arrivée estimées des signaux en repère capteur de l'étape 4, qui peuvent être fausses en cas de leurrage, on peut préférer chercher les angles d'attitude donnant le maximum d'énergie après les sommes pondérées en sortie des corrélateurs. On cherche alors à maximiser la somme des énergies sur tous les signaux satellites démodulés : def 111(q),0,11) = In=1,..,N II Zp,(ip,O,iiJ) II2 def j ((PM = En=1,..,N II CnT((),O,v ) . Zp capteurs n 112 def If(q),0,111 ) = En=1,..,N E CnT((p,e,y ) . Zp capteurs n ]-E CT( cs ,y ) - ZP capteurs n On peut utiliser la méthode du gradient : Re([acT n (v, ayi ReaDC: (vt, ae Red acnT Dço 6 ))-Z P capteurs g)) P capteurs n].[C nT (y', O, ÇO).Z p capteur s ) n].[C nT (V, ÇO).Z p capteursnr ) Ôrad`l (y f,0,ço)= n=1, .,N , q)).Z capteursn1.[C nT (V, O, ça) .Z capteurs nn C tn(yf , O,y)) C (V , O, q)) Cn(V, 0,y))= Czn(V,O,ÇO) M ,n(y13 09 V)- M x 120 cni,'(v,19,0=expi(gm,'(v,e9,0+ ;6Pcctp'u r m(e n(V38 3 v),r,,(y.',e,s9))) 8 m,n(Vf 919 3 0 = 17 n capteurs (y.', t 950 - Fm capteurs - (2z/2) (en rad) _ cos(yn ). cos (o-n ) sin (yn ).cos(cin ) - sin(0n) _ capteurs(v,e,s0)= L'angle de site en repère capteur s'écrit alors: Gri(1V,0,(p) = sin1( uzn(WM,Y) ) L'angle d'azimut en repère capteur s'écrit alors : Ya(W,e,y) = tan-1( uxn(w,e,y) , uYa(W,O,(p) ) Il faut alors appliquer trois pondérations "delta" par satellite ce qui est coûteux en nombre de calculs. On peut alors utiliser des discriminateurs. On prend comme paramètres de directions des satellites les cosinus directeurs (ot,p) en repère capteurs ou deuxième repère. On calcule des discriminateurs en sortie des corrélateurs pour mesurer les erreurs commise sur les directions des satellites. def Da n(,3) = Re[ Zan(4) . Zn(oc,5)* l III Z 112/ Pa = Re( Za / Z* ) / Pa def Db n(a,13) = Re[ Zb a(a,13) . Zn(a,13)* l / II Z Ir / Pb = Re( Zb / Z* ) / Pb Les calculs sont effectués avec les mêmes formules que celles utilisées avec les angles de site et d'azimut (szy,y) Somme pondérée ponctuelle : def p n(a, fi) = CT (a, 13).Z p n capteurs Sommes pondérées avec les pondérations complexes "delta" : def T Za n (a, fi) = C a (a, f3).Z p n capteurs def T Z b n (a, fi) = (a, fl)zp n capteurs def CT C a (a, fl) = a a (a, fi) def aCT C (a, fi) = (a, f3) afi 10 C ,(a, 16) C 2(a, ig) C(a, fl) = C 3(a, fi) _C (a, 16) Mx! Cm (a, fi) = exP i( g ',(a, B capteur ni( cr(a, fi), Ka, 16)) gm (a, /3) = il-capteurs (a, 16) - im capteurs - (27r I 2) 15 capteurs (a, fi) = fi ux(a, 16) uy(ot, fi) - a2 - fl2 uz (a, fi) _ cs(a,13 ) = tan-1( Va2 +,62 V1- a2 + fl2 ) 20 Les pentes Pa et Pb sont calculées de la manière suivante : def CT (a, ,6).0 (a o , A ) F1(00) = def ac T moc,p)= aa (a,fl).C(ao,flo) = 2--( def acr Ab(CC,5) = a u (a,fl).C(ao,flo) = 2- (11) (43) afl def a Re[ Aa.II] )(ao,P / II n(P,o,Ro) 112 aoe def = 4.(Re[Ab.n] ) (ao,P / II n(a0,130) 112 pour des valeurs (ao, [30) correspondant aux dernière estimées de (a,13) Si les valeurs de (a, (3) sont proches des vraies valeurs (m-, R : ,raie r vrai, - Da n(a,p) a- ( a - avrai n ) + bruit Dbn(a,13) p - Pvrai ) + bruit On utilise également un modèle d'observation. La relation entre les cosinus directeurs (an,f3n) et les angles d'Euler (v,e,y) n'est pas linéaire. uxi,(y1,0,ço) cet, (V, (0) Ûn capteurs (V f >t 9, g') = uy(v,0,ço) fin (v,(9,(0) 15 _uz,,(v,0,ço)_ 111-an (v, 0, yo) 2 - fin (o)2 Ps PA 10 1 n capteurs (V1 6 gel) = M (V 0,Ço). n TGL M (1,ti , 0, ço) = 1 0 0 cos 0 0 -. 9 cos yi sin 0- 0 cos y° sin y) -sinço cossinO 0 1 0 - sin yt cos 0 0 0 cos 6, 0 0 1 On linéarise cette relation autour d'un point de fonctionnement (Nr,O,y) pour établir un modèle d'observation linéaire : def dan = ap,(v+dv,0+d0,ço+dço) - def n = J3n(yi-i-dv,0+d0,v+dço) - fin(v,0,ço) 20 25 da =auxn (v,19,q)).dv + aux (V1,0,v).cle9 + aux" (v, 0,y0).dq) n a ao Dy) dfin = auy" (v, 0,q)).dyi + auy" (v,O, ).dO + auy" t9,(0).4 a DO D Ça Pour le calcul des dérivées partielles on peut utiliser la méthode des différences finies : aux n (V1 e50 = dV19 9, (°) eiXii(yIle5(°) )1 CIV a yf En utilisant d'autre notations : dan= ana . d + an,2 .dt9 + an,3 dço On= b,1 .dyt + bn,2 .dt9 + bn,3 .dço Si on considère tous les satellites visibles on a 2.N équations pour 3 inconnues. dot, dfl, a1,1 a1,2 a1,3 dço daN dfiN b1,2 b1,2 b1,3 dO aN,1 aN,2 aN,3 dv _hN,1 bN,2 bN,3 soit : da, =H dço dO avec H= a1,1 a1,2 a1,3 dfl, b1,2 b1,2 b1,3 daN dv aN,1 aN,2 aN,3 dfiN _bN,1 bN,2 bN,3 _ iv vrai , dO = 0 - °vrai On cherche à remonter au erreurs sur les angles d'Euler (dy = w w - , dy) = - (pvrai) à partir des erreurs sur les cosinus directeurs ( dan = an(w,0,9) rai n , dpn = 13n(w,0,y) - Pevrai ) Les erreurs sur les cosinus directeurs sont observables grâce aux discriminateurs : dan = an(lif,e,(P) - ai) n = Da n( an(111,9,(P) 13n(111,0,(P) ) + bruit dtn = Pn(w,e,(P) - 130 n = Db n( an(110,(P) 511()11,0e) ) + bruit Aussi, on peut utiliser une résolution par les moindres carrés : d ç°estimé dal dfil avec H* = (HT WT WH)_1 H T WT pour pseudo d0''mé = H* doeN dvf inverse. La matrice W permet de pondérer les mesures en fonction du rapport signal sur bruit estimé sur chaque signal satellite. 1/0%21 o O 1/0-1,21 w= O 0 1/o-,_ 1 "an="bn= C I IVO estimé n Après la résolution on corrige les angles d'Euler estimés a priori : _ _ _ _ _ q'estimé q'estimé précédent dq'estimé °estimé °estimé précédent + d °estimé _V estimé _ Y' estimé précédent _ _d V estimé _ _ En variante des moindres carrés, on peut utiliser une résolution par filtre de Kalman (étendu), comme suit : Etat : (01 (i) estimé i Réel :Xi = (9, Estimé :Xi,i = °estimé i à l'instant i _ Vf estimé i Initialisation : Poio = o OOoo 0- Xo/o q) estimé initial °estimé initial estimé initial 0 o . Propagation : = Vecteur d'état propagé +Q Matrice de covariance propagée O o Qr= O ge 0 Bruit d'état 0 0 qw_ Le bruit d'état est choisi en fonction de la dynamique sur l'attitude du porteur. Recalage : a1,2 a1,3 dfil b1,2 b1,2 b1,3 Y = Hi= daN aN,1 _ bN 1 aN ,2 aN,3 bN,2 bN,3 " a 1 o Crbl 2 29 Ri= o. a2N 0 0 .. 0 0-2 bN _ L'innovation Yi = Zi est déjà contenue dans les discriminateurs, il n'est pas nécessaire de la recalculer. K = (H iT + Ri)-1 Pi, iT Gain de recalage = +Ki.Yi Vecteur d'état recalé Pui = (Id -Kil Matrice de covariance recalée Le synoptique de la figure 11 illustre ce qui précède Sur la figure 12 est illustré un récepteur 10 de signaux d'un système de navigation par satellites, muni d'une antenne 11 à capteurs à rayonnement contrôlé, comprenant, pour chaque satellite : - des moyens de calcul 12 de la position du satellite dans un premier repère TGL centré sur le récepteur à partir de la position du récepteur et de la position du satellite; - des moyens de détermination 13 de l'attitude de l'antenne dans le premier repère - des moyens de détermination 14 des premiers angles de site et d'azimut du satellite dans un deuxième repère lié à l'antenne, à partir de la position du satellite dans le premier repère et de l'attitude de l'antenne dans le premier repère; - des moyens de détermination 15 des deuxièmes angles de site et d'azimut du satellite dans le deuxième repère, par recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne; - des moyens de calcul 16 de la valeur d'une fonction de vraisemblance entre les premiers angles de site et d'azimut du satellite et les deuxièmes angles de site et d'azimut du satellite; et - des moyens de détection 17 d'un risque de réception de signaux leurrés lorsque ladite valeur de la fonction de vraisemblance est inférieure à un seuil. Un tel récepteur met en oeuvre le procédé précédemment décrit. Dans la présente demande, les notations utilisées ont pour définitions : expi( ) = exp( j. ) = dg' (j2 = -1) Définitions : Fr: Fréquence résiduelle des signaux numérisés en bande de base (Fr petit devant Fe) Dop : Doppler sur le signal satellite n (Pm,n Déphasage du signal satellite n sur le capteur m C'( t ) : code d'étalement du satellite n Ret n : retard du signal satellite n R(T) : Fonction d'autocorrélation du code R(t) = 1 - 'Tl si ITI 1000) et quel que soit k R(t) 1«/ L [kL, (k+1)1_] Cn( U ) . Cp( U +'r) du 0 si p # n Signaux reçus : Le signal reçu sur l'antenne est la somme des signaux reçus des satellites visibles plus du bruit Sreçu m (t) = E n=1,..,N Sreçu m,n (t) + bruit Signal reçu du satellite n sur le capteur m : Sreçu m,n (t) = A.exPi( (Pporteuse reçu m,n(t)) - Cn( (Pcode reçu n(t) ) Sreçu m,n (t) = A.expi( 27c.(Fr+Dop n).t + 0,',n ) . Cn( Fcode-(t+Ret n) ) Signal local dans le canal n: Slucal n (t) = expi( Yporteuse locale n(t) ) - Cn( (Pcode local n(t) ) Slocal n (t) = expi( -27r.(Fr+Dop n).t - en). Cn( Fcode.(t+Ret n) Corrélations : ZnIn = 111- I [kT, (k+1)1] Sreçu m(t) - Slocal n(t) dt T : temps d'intégration cohérent (typiquement T = 20 ms) Les codes satellites étant décorrélés (R,,p('r) a 0) on peut écrire : Zn,,n(k)= 1/T I [kT, (k+1)T1 Sreçu m,n(t) - Slocal n(t) dt + bruit Zm,n(k)= lir I [kT, (k+1)1] expi( (Pport reçue m,n(t) (Pport locale n(t) )- Cn( (Pcode reçu n(t) )- Cn( (Pcode local n(t) ) dt Zn,,,,(k) = A. expi( On',n - On ) . 111_ J [kL, (k+1)L] Cn( u + t). Cn( U ) du + bruit = A. expi( Oro - en ) . R(T) + bruit On suppose que le code local du canal n est en phase avec le code reçu du satellite n, grâce à la boucle de code, réalisée suivant l'état de l'art. D'où 'G = 0 et R(T) = 1 Zrn,n(k) = A. expi( Onn - en ) + bruit REVENDICATIONS1. Procédé de détection de signaux destinés à leurrer un récepteur de signaux d'un système de navigation par satellites, le récepteur étant muni d'une antenne à capteurs à rayonnement contrôlé, comprenant les étapes suivantes consistant à, pour chaque satellite : - calculer (1) la position du satellite dans un premier repère TGL centré sur le récepteur à partir de la position du récepteur et de la position du lo satellite; - déterminer (2) l'attitude de l'antenne dans le premier repère; - calculer (3) des premiers angles de site et d'azimut du satellite dans un deuxième repère lié à l'antenne, à partir de la position du satellite dans le premier repère et de l'attitude de l'antenne dans le premier repère; 15 - déterminer (4) des deuxièmes angles de site et d'azimut du satellite dans le deuxième repère, par recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne; - calculer (5) la valeur d'une fonction de vraisemblance entre les premiers angles de site et d'azimut du satellite et les deuxièmes angles de 20 site et d'azimut du satellite; et - détecter (6) un risque de réception de signaux leurrés lorsque ladite valeur de la fonction de vraisemblance est inférieure à un seuil. 2. Procédé selon la revendication 1, dans lequel ladite recherche 25 itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne est effectuée à partir d'un discriminateur d'angle de site et d'un discriminateur d'angle d'azimut. 3. Procédé selon la revendication 1 ou 2, dans lequel l'étape de 30 détermination (2) de l'attitude de l'antenne dans le premier repère utilise des données fournies par une référence inertielle. 4. Procédé selon la revendication 1 ou 2, dans lequel l'étape de détermination (2) de l'attitude de l'antenne dans le premier repère utilise 35 ladite détermination (4) des deuxièmes angles de site et d'azimut dessatellites dans le deuxième repère et procède par une recherche itérative des angles d'attitude de l'antenne donnant le minimum de la fonction de vraisemblance. 5. Procédé selon la revendication 1 ou 2, dans lequel l'étape de détermination (2) de l'attitude de l'antenne dans le premier repère procède par une recherche itérative des angles d'attitude de l'antenne donnant le maximum d'énergie après les sommes pondérées en sortie des corrélateurs des canaux de traitement respectivement associés aux satellites. 6. Procédé selon la revendication 5, dans lequel ladite maximisation utilise une méthode de résolution par les moindres carrés. 7. Procédé selon la revendication 5, dans lequel ladite maximisation utilise une méthode de résolution par filtrage de Kalman. 8. Procédé selon la revendication 6 ou 7, dans lequel on utilise des discriminateurs de cosinus directeurs représentatifs de la position des satellites dans le deuxième repère lié à l'antenne. 9. Procédé selon la revendication 8, dans lequel on utilise une méthode à modèle d'observation de la relation entre les cosinus directeurs des directions des satellites et les angles d'Euler des directions des satellites. 10. Récepteur (10) de signaux d'un système de navigation par satellites, muni d'une antenne (11) à capteurs à rayonnement contrôlé, comprenant, pour chaque satellite : - des moyens de calcul (12) de la position du satellite dans un premier repère TGL centré sur le récepteur à partir de la position du récepteur et de la position du satellite; - des moyens de détermination (13) de l'attitude de l'antenne dans le premier repère - des moyens de détermination (14) des premiers angles de site et d'azimut du satellite dans un deuxième repère lié à l'antenne, à partir de la position du satellite dans le premier repère et de l'attitude de l'antenne dansle premier repère; - des moyens de détermination (15) des deuxièmes angles de site et d'azimut du satellite dans le deuxième repère, par recherche itérative d'un maximum d'une somme complexe pondérée des signaux démodulés reçus par l'antenne; - des moyens de calcul (16) de la valeur d'une fonction de vraisemblance entre les premiers angles de site et d'azimut du satellite et les deuxièmes angles de site et d'azimut du satellite; et - des moyens de détection (17) d'un risque de réception de signaux leurrés lorsque ladite valeur de la fonction de vraisemblance est inférieure à un seuil.