FR 3093819 A1 20200918 FR 1902597 A 20190314 Description Titre de l'invention : Procédé de contrôle d'intégrité d'un signal de radionavigation par satellite Domaine technique de l'invention [0001] L'invention concerne le domaine de systèmes de traitement de signaux de radiona- vigation par satellites ou pseudolites, et concerne plus précisément les récepteurs GNSS comprenant une fonction de traitement anti-interférence. [0002] L'invention porte plus spécifiquement sur un procédé et une fonction de contrôle d'intégrité au regard d'un risque de faux accrochage pour un récepteur GNSS comprenant une fonction de traitement anti-interférence. Etat de la technique [0003] Les systèmes de positionnement par satellites (désignés GNSS pour « Global Na- vigation Satellite System » en anglais) mettent en oeuvre plusieurs satellites transmettant des signaux radioélectriques et un récepteur GNSS placé à la position à localiser estimant les distances dites « pseudo-distances », qui le séparent des satellites à partir des temps de propagation des signaux satellites reçus multipliés par la vitesse de la lumière dans le milieu traversé. A partir des mesures des pseudo-distances séparant le récepteur des satellites visibles et de la connaissance des positions de ces satellites, le récepteur déduit sa position précise en latitude, longitude, et en altitude dans un repère terrestre par une résolution numérique de type triangulation. Plus les positions des satellites sont connues du récepteur avec précision et plus les mesures des pseudo-distances faites par le récepteur sont précises, plus la localisation obtenue est précise. Les positions des satellites sont déterminées à partir d'un réseau de stations au sol indépendant des récepteurs de positionnement. Les positions des satellites sont communiquées aux récepteurs de positionnement par les satellites eux-mêmes par les données modulant les signaux émis [0004] Les pseudo-distances sont déduites par les récepteurs de positionnement des retards apparents présentés par les signaux reçus. Ces signaux sont datés à l'émission par les horloges des satellites, qui sont synchronisées sur le temps système, et à la réception par l'horloge du récepteur, pouvant présenter un biais par rapport au temps du système. Les distances déduites des temps de propagation ainsi mesurés sont appelées pseudodistances car elles sont entachées d'une erreur commune pouvant être importante, due au biais de l'horloge du récepteur, biais qui peut être éliminé si les signaux d'au moins quatre satellites sont reçus. [0005] Les signaux radioélectriques émis par des satellites parcourant de grandes distances étant émis avec des puissances limitées, ils parviennent avec de très faibles puissances 2 aux récepteurs, noyés dans un bruit radioélectrique dû à l'environnement physique. Pour faciliter la réception des signaux émis par les satellites et les rendre moins sensibles aux parasites (ou interférences), notamment à bande étroite, les largeurs de bande des signaux émis sont étalées au moyen de la technique de la bande étalée, également désignée « technique d'étalement de spectre ». Ainsi, les signaux émis par les satellites sont constitués d'une porteuse de signal modulée par un code d'étalement. Plus spécifiquement, les signaux émis par les satellites sont formés par modulation de la porteuse du signal avec un code d'étalement formé par une séquence binaire pseudoaléatoire. La modulation de la porteuse par un code d'étalement étale le spectre, ce qui accroît la résistance du système aux interférences. Et, en outre, cela permet de dissocier les satellites (en utilisant un code d'étalement différent par satellite). En d'autres termes, l'identification du satellite émetteur d'un signal radioélectrique peut se faire à partir du code d'étalement qui lui est réservé. Dans la suite de la présente invention, on pourra utiliser par simplification le terme « code » pour désigner le terme « code d'étalement ». [0006] Dans le récepteur, le signal satellite est reçu via au moins une antenne (une antenne pouvant comprendre un ou plusieurs capteurs) et les informations contenues dans un signal radioélectrique du satellite sont extraites par deux démodulations successives, afin d'en déduire une mesure de temps de propagation, de Doppler ou de phase de porteuse, et éventuellement de données transmises : une première démodulation à l'aide d'une porteuse engendrée localement par un oscillateur piloté par une boucle de poursuite en fréquence, dite FLL (sigle issu des termes anglais Frequency Lock Loop), ou en phase, dite PLL (sigle issu des termes anglais Phase Lock Loop), également désignées par « boucle de phase de porteuse », ou « boucle de porteuse » permettant de transposer le signal reçu en bande de base et une deuxième démodulation à l'aide d'un code d'étalement engendré localement par un générateur de code piloté par une boucle de poursuite en retard dite DLL (sigle issu des termes anglais Delay Lock Loop), également désignée par «boucle de code » permettant de démoduler la suite de données binaires présente dans le signal reçu. En général, il y a au niveau du récepteur un canal par satellite (désigné « canal satellite »), avec un code d'étalement par canal propre au satellite poursuivi, ainsi qu'une boucle de code et une boucle de porteuse par canal, comme cela est expliqué par exemple dans la demande de brevet FR2892202. [0007] La démodulation consiste à asservir un signal généré localement, de manière à former une image du signal reçu en provenance du satellite considéré caractérisé par un code d'étalement propre et une porteuse, c'est-à-dire en répliquant localement le code, et en recherchant le maximum de corrélation entre ce signal reçu et le signal répliqué (nommé « signal local »). L'asservissement est réalisé par la boucle de porteuse, qui pilote la phase de la porteuse locale, et par la boucle de code qui pilote la 3 phase du code local. La boucle de porteuse mesure un écart de phase de porteuse entre le signal local et le signal reçu grâce à la corrélation avec un signal local. La boucle de code mesure un écart de phase de code entre le signal local et le signal reçu grâce à la corrélation avec des signaux locaux modulés par des codes dérivés (avance, retard ou ponctuel). On appelle cette phase de démodulation la « phase de poursuite ». Le ternie « poursuite » représente la capacité à se synchroniser avec les signaux satellites en utilisant par rétroaction les écarts de phase entre un signal local et un signal reçu. [0008] La démodulation comprend, préalablement à la phase de poursuite, une phase d'acquisition qui permet d'initialiser le fonctionnement des boucles de poursuite. En effet, ni la position (ou phase) du code reçu, ni la fréquence Doppler ne sont connus a priori, or les boucles de poursuite ne peuvent fonctionner que si la position du code et la fréquence Doppler sont proches de celles du signal issu du satellite considéré. Le principe mis en oeuvre lors de la phase d'acquisition est d'essayer plusieurs hypothèses sur la position du code d'étalement, ainsi que sur la valeur du Doppler et de rechercher un pic de corrélation entre le signal local et le signal reçu, avec un pas d'échantillonnage suffisamment fin pour ne pas manquer le pic principal de la fonction d'auto-correlation du code d'étalement. [0009] Par fonction d'auto-corrélation, on entend une corrélation du code d'étalement par le même code décalé dans le temps et filtré. La fonction donne un résultat de corrélation en fonction du décalage. Le décalage désigné T correspond à une avance du code reçu par rapport au code local. Les valeurs de la fonction d' autocorrélation peuvent être complexes, avec une partie réelle et une partie imaginaire, notamment à cause du filtrage appliqué au signal reçu. [0010] Une fois qu'un pic a été trouvé, la recherche de la phase du code et du Doppler est affinée en diminuant le pas d'échantillonnage, autour du premier pic détecté. Quand la précision obtenue est jugée suffisante, on passe alors en phase de poursuite : on ferme les boucles pour les faire converger vers le maximum de corrélation. Lorsque le pic de corrélation est déterminé et que les boucles sont fermées, on dit que les boucles sont accrochées », on considère alors que le code d'étalement engendré localement est en phase avec le code du satellite recherché contenu dans le signal reçu. Cette recherche de pic de corrélation est également nommée « accrochage ». On précise que le terme pic » désigne un extremum local, qui peut être un maximum local lorsqu'il concerne des valeurs absolues. [0011] Le temps de propagation d'un satellite est obtenu à chaque instant par la différence entre la date d'émission du signal reçu, calculée à partir de la phase du code local, image du code reçu, et la date de réception du signal reçu, donnée par l'horloge du récepteur. Et le positionnement du récepteur est obtenu en mesurant la pseudo-distance entre un satellite et le récepteur à partir du temps de propagation du signal entre ce satellite et le récepteur. [0012] Un exemple de calcul de corrélation réalisé lors de la phase d'acquisition est par exemple décrit dans la demande de brevet FR2898998. Un code d'étalement est typiquement constitué d'un grand nombre de divisions temporelles élémentaires, appelées « chips » de code. Dans la demande de brevet FR2898998, Da, désigne la durée moyenne de la division temporelle élémentaire du code (pendant chaque chip, le code conserve une valeur constante), Cl (0 désigne le code local généré par un récepteur d'un système de positionnement par satellite, et S(t) le signal de radionavigation reçu par le récepteur, ce signal analogique étant échantillonné à une fréquence F, supérieure à 2/13,1'. Un calcul de corrélation entre un signal reçu et un code généré localement est réalisé en plusieurs étapes : - réaliser un échantillonnage du code Cl synchrone de l'échantillonnage du signal reçu, - multiplier deux à deux les échantillons du code Cl et du signal reçu, correspondant à une même date d'échantillonnage, sur la durée Dcad, pour déterminer N,,h résultats de produits d'échantillons, Noel, étant le produit de la fréquence F, et de la durée Dcad, ; - sommer les résultats de produits d'échantillons. [0013] Le résultat du calcul de corrélation peut ainsi s'exprimer sous la forme suivante : [0014] [Math 1]:Corrélationl - NeuhS t- ) citt \où i est un indice = [f d'échantillonnage, ti est la i-ème date d'échantillonnage du signal reçu S, S(t;) est la valeur du signal S à la date ti et Cl (ti) est la valeur du code à la date t. 100151 Malgré la technique d'étalement de spectre, les systèmes de radionavigation par satellite peuvent être perturbés par des sources interférentes, intentionnelles ou non, par exemple des sources émettant un signal sur une fréquence proche de celle du signal de radionavigation ou présentant des harmoniques autour de la fréquence du signal de radionavigation, ce qui peut compromettre la fiabilité du système de radionavigation par satellite. 10016] Des techniques existent afin de caractériser et de supprimer ces sources interférentes. 10017] Par exemple, il est connu d'appliquer un traitement aux signaux reçus, et ce, avant leur démodulation, qui consiste à couper le signal reçu lorsqu'un tel signal interférent est détecté, ce qui entraine une perturbation de la réception des signaux utiles. Ce système est connu sous le terme anglais de « blanking ». Selon le nombre de signaux interférents reçus, la proportion de coupures ou « blanking » est plus ou moins grande et à la limite on ne reçoit plus assez de signal et il n'y a alors aucun moyen de récupérer un signal utile. [0018] D'autres techniques existent, telles que celles explicitées par exemple dans les demandes de brevet FR2981168 ou FR3017716, qui utilisent une technique de filtre linéaire appliquée aux signaux reçus par le récepteur, et ce, avant leur démodulation. Les coefficients du filtre linéaire sont calculés de manière à ce que la réponse en fréquence atténue l'énergie de l'interférence dans les bandes de fréquence où elle est présente. 11 peut s'agir d'un filtre STAP des termes anglais « Spatio Temporal Adaptive Processing » utilisé typiquement lorsque le récepteur comprend plusieurs antennes ou une antenne à plusieurs capteurs ou d'un filtre TAP des termes anglais « Temporal Adaptive Processing » utilisé lorsque le récepteur comprend une antenne à un seul capteur. Une autre technique alternative au TAP est un filtre FDAF des termes anglais « Frequency Domain Adaptive Filtering ». [0019] Si l'énergie d'une interférence se situe au centre du spectre du signal satellite, il en résulte un trou au milieu du spectre dont la conséquence est l'apparition de pics secondaires à côté du pic principal de la fonction d'auto-corrélation. [0020] Or comme indiqué plus avant, lors de la phase d'acquisition, un récepteur ou un canal satellite du récepteur recherche un pic de corrélation pour s'accrocher dessus en verrouillant les boucles de code et de porteuse. Les boucles restent accrochées sur un pic local de la fonction de corrélation lors de la phase de poursuite. Lorsque l'on a plusieurs pics dans la fonction de corrélation, il y a un risque que le canal satellite s'accroche sur un pic secondaire, pouvant conduire à des erreurs de positionnement significatives voire inacceptables, jusqu'à quelques dizaines voire quelques centaines de mètres. Selon l'invention, cet évènement sera nommé « faux accrochage ». Un faux accrochage peut dégrader le niveau d'intégrité du système de navigation par satellite. L'invention vise à surmonter l'inconvénient précité. [0021] L'invention vise à disposer d'une fonction de contrôle d'intégrité d'un récepteur GNSS comprenant un traitement anti-interférence, plus précisément une fonction apte à détecter un risque de faux accrochage, afin d'éviter un risque de mesures erronées non détectées. [0022] De manière avantageuse, si un risque de faux accrochage a été détecté, l'invention vise à déterminer si le risque est avéré et donc si un faux accrochage a effectivement eu lieu. [0023] De manière avantageuse, l'invention vise à corriger le défaut d'accrochage avéré. Exposé de l'invention [0024] L'invention permettant de résoudre les problème précités concerne un procédé de contrôle d'intégrité du traitement d'un signal de radionavigation émis par un satellite d'un système de positionnement par satellites, ledit signal étant reçu par un récepteur comportant des moyens de réception et des moyens de traitement, lesdits moyens de traitement comprenant un filtre linéaire anti-interférence, ledit procédé de contrôle d'intégrité étant caractérisé en ce qu'il comporte au moins une première phase de 6 détection d'un risque de faux accrochage comprenant les étapes suivantes : - une étape de récupération d'une fonction d'auto-corrélation théorique nominale du signal reçu non traité par le filtre linéaire anti-interférence ; - une étape de détermination d'une fonction d'auto-corrélation théorique moyenne dudit signal reçu et traité par le filtre linéaire anti-interférence sur une période d'intégration définie ; - une étape de détermination du nombre de maximum locaux du module ou du module au carré de la fonction d'auto-corrélation théorique moyenne, un risque de faux accrochage étant détecté si ledit nombre de maximum locaux est supérieur ou égal à deux. [0025] Le signal reçu non filtré a une fonction d'auto-corrélation théorique nominale. Le signal reçu et filtré a une fonction d'auto-corrélation théorique qui est modifiée du fait du filtrage. L'invention nécessite de connaître cette fonction d'auto-corrélation théorique modifiée. [0026] Ainsi l'invention consiste à calculer la fonction de corrélation théorique moyenne du signal reçu et traité par le filtre linéaire anti-interférence sur une période d'intégration définie. Cette fonction d'autocorrélation théorique moyenne renseigne sur la présence de pics secondaires (donc d'au moins deux pics, le pic principal et au moins un pic secondaire) et donc un risque de faux accrochage. Plus largement, l'invention consiste donc en la recherche des pics de la fonction d'autocorrélation théorique moyenne. [0027] La fonction d'autocorrélation d'un signal de radionavigation émis par un satellite n'est généralement pas réelle mais complexe (avec une partie réelle et une partie imaginaire), notamment lorsqu'un filtrage anti-interférence a été appliqué au signal émis C'est pourquoi, on recherche le (ou les) maximum du module ou du module au carré de la fonction d'autocorrélation théorique moyenne. Le cas où la fonction d'autocorrélation théorique moyenne est réelle est un cas particulier. [0028] Un maximum local se définit comme étant un maximum sur un intervalle non nul. En pratique, comme on travaille généralement sur des valeurs discrètes espacées de la période d'échantillonnage Te, l'intervalle peut avoir une largeur égale à 2 fois la période d'échantillonnage Te (qui peut correspondre à l'espacement minimal entre les corrélateurs avance et retard comme cela sera expliqué plus après). [0029] Selon un mode de réalisation, le procédé comprend une étape de détermination d'un filtre anti-interférence moyen sur la période d'intégration, préalablement à l'étape de détermination de la fonction d'auto-con-élation théorique moyenne, ladite fonction d'auto-corrélation théorique moyenne étant déterminée à partir de la fonction d'auto-corrélation théorique nominale et dudit filtre anti-interférence moyen. [0030] Selon un mode de réalisation, le filtre anti-interférence moyen sur la période d'intégration est déterminé à partir des différentes valeurs des coefficients du filtre 7 anti-interférence calculés sur ladite période d'intégration. [0031] Selon un mode de réalisation, le procédé comprend en outre : - une étape de définition d'un premier seuil, ledit premier seuil étant fonction de l'écart entre le rapport signal sur bruit maximum et un seuil de poursuite ; - le nombre de maximum locaux déterminés correspondant au nombre de maximum locaux pour lesquels la valeur du module ou du module au carré de la fonction d'auto-corrélation moyenne finale est supérieure ou égale audit premier seuil. [0032] En particulier, le premier seuil est défini par la relation mathématique Math40 définie plus après et reproduite ici : max( II max S1 - 1010 (A0) II 2 ) où x représente l'écart entre le rapport signal sur bruit maximum et un seuil de poursuite et AO la fonction d'auto-corrélation théorique nominale. [0033] Selon un mode de réalisation, le procédé comprend en outre une étape de dé- clenchement d'alarme si le nombre de maximum locaux déterminés est supérieur ou égal à deux. [0034] Selon un mode de réalisation dans lequel les moyens de traitement du signal reçu comprennent un canal de traitement par satellite, chaque canal comportant une boucle de code apte d'une part à générer localement un code d'étalement caractéristique du satellite associé audit canal de manière à former un signal local et apte d'autre part à asservir ledit signal local avec ledit signal reçu en utilisant une voie de corrélation entre ledit signal reçu et ledit signal local, la période d'intégration est déterminée en fonction de la période de ladite boucle de code, également dénommée « période de boucle ». [0035] Selon un mode de réalisation particulier, la période d'intégration est sensiblement égale à la période de boucle. [0036] Selon un mode de réalisation dans lequel les moyens de traitement du signal reçu comprennent un canal de traitement par satellite, chaque canal comportant une boucle de code apte à générer localement des codes d'étalement caractéristiques du satellite associé audit canal, lesdits codes comprenant au moins un code ponctuel, un code retard et/ou un code avance formant chacun un signal local, chaque signal local ponctuel, retard et/ou avance étant corrélé avec le signal reçu en utilisant respectivement une voie de con-élation ponctuelle, une voie de corrélation retard et/ou une voie de con-élation avance, le procédé comporte en outre une deuxième phase de détection d'un faux accrochage, ladite deuxième phase comprenant : - une étape de calcul de vraisemblance entre la fonction d'auto-corrélation théorique moyenne et une fonction d'auto-con-élation mesurée pour chacune de l'ensemble des voies de corrélation ponctuelle, retard et/ou avance ; 8 - une étape de comparaison du résultat de calcul de vraisemblance avec un deuxième seuil, un faux accrochage étant détecté par un résultat du calcul dc vraisemblance inférieur audit deuxième seuil. [0037] Selon un mode de réalisation particulier, le résultat de l'étape de calcul dc vrai- semblance utilise un critère nominal de vraisemblance défini par la relation mathématique Math41 définie plus après et reproduite ici : (, = 1A )'A Ines(Cii)/ ti'CNC ( i 0 II A (g i) II ) 0 II A mes(qi) Il ) ci est la position du corrélateur n°i et NC le nombre de voies de corrélations. 100381 Selon un mode de réalisation, le procédé comporte en outre une troisième phase de recherche du maximum local principal, ladite troisième phase effectuant une recherche du maximum local principal si le résultat du calcul de vraisemblance est inférieur au deuxième seuil, et comprenant les étapes suivantes : - une j-ème étape de calcul de vraisemblance entre la fonction d'auto-corrélation théorique moyenne translatée d'une j-ème valeur de translation par rapport à la position des voies de corrélation, et la fonction d'auto-corrélation mesurée pour chacune de l'ensemble des voies de corrélation, la j-ème valeur de translation correspondant à un j-ème maximum local ; - une étape de comparaison du résultat du j-ème calcul de vraisemblance avec un troisième seuil, un résultat dudit j-ème calcul dc vraisemblance supérieur audit troisième seuil permettant de déterminer si la j-ètne valeur de translation correspond au maximum local principal. [0039] Le calcul est réalisé pour chaque j-ème maximum local déterminé lors de la première phase de détection d'un risque de faux accrochage. Ainsi j est au moins égal à deux. On fixe de préférence un troisième seuil de telle sorte qu'on soit assuré de déterminer le maximum local principal. Si on trouve au moins deux maximum locaux au-dessus du troisième seuil on peut relever le troisième seuil ou reproduire le calcul dc vraisemblance en considérant un quatrième seuil supérieur au troisième seuil. [0040] Selon un mode de réalisation particulier, le résultat de l'étape de calcul dc vrai- semblance utilise un j-ème critère de vraisemblance défini par la relation mathématique Math42 définie plus après et reproduite ici : 1;1 A (q + cliyA','(q in Critère j = (17c II A II 2 ) ;N/C0 II A in e5(ci ?) ) où di est la j-ème valeur de translation, qi est la position de la i-ème voie de corrélation et NC le nombre de voies de corrélations. [0041] L'invention concerne également un récepteur d'un système de positionnement par sa- Critère0 = où 9 teintes comportant des moyens de réception aptes à capter un signal émis par lesdits satellites et des moyens de traitement aptes à traiter le signal reçu, lesdits moyens de traitement comprenant un filtre linéaire anti-interférence, ledit récepteur comprenant une première fonction apte à implémenter la première phase de détermination d'un risque de faux accrochage du procédé selon l'invention. [0042] Selon un mode de réalisation, le récepteur comprend en outre une deuxième fonction couplée à la première fonction et apte à implémenter la deuxième phase de détection d'un faux accrochage du procédé selon l'invention. [0043] Selon un mode de réalisation, le récepteur comprend en outre une troisième fonction couplée à la première et la deuxième fonction et apte à implémenter la troisième phase de recherche du maximum local principal du procédé selon l'invention. [0044] Ainsi, l'invention permet de disposer d'un procédé et d'une fonction de contrôle d'intégrité pour un récepteur d'un système de traitement de signaux de radionavigation par satellites ou pseudolites, et plus précisément pour un récepteur comprenant un traitement anti-interférence. L'invention permet de détecter un risque de faux accrochage, comme défini plus avant, afin d'éviter un risque de mesures erronées non détectées. L'invention permet a minima de transformer une non intégrité en non continuité, à savoir que le récepteur préfère ne pas fournir de mesure, plutôt que de risquer de fournir une mesure erronée sans prévenir. Par exemple, dans ce cas, un pilote d'aéronef pourra choisir de dérouter l'aéronef ou de temporiser l'atterrissage le temps de retrouver des mesures de position garanties intègres par le récepteur, en évitant ainsi un risque d'accident. [0045] En outre, grâce à la deuxième phase de détection d'un faux accrochage, l'invention permet de vérifier si le risque de faux accrochage est avéré ou non. Cela permet d'éviter une non-continuité de la mesure en réduisant le risque de fausse alarme. Par exemple, dans le cas d'un atterrissage d'aéronef, cela permet d'éviter un déroutage coûteux, voire dangereux. [0046] Enfin, en cas de faux accrochage avéré, la troisième phase de recherche du maximum local principal permet de retrouver rapidement une mesure intègre en évitant de refaire une phase d'acquisition, généralement chronophage, réduisant ainsi la durée de non continuité de la mesure. [0047] D'une manière générale, le procédé selon l'invention comprend des opérations qui peuvent être exécutées de manière rapide et qui sont simples de mise en oeuvre. Elles n'augmentent pas significativement la charge du (ou des) processeur(s) en place. L'invention peut notamment être mise en oeuvre sous la forme de lignes de codes supplémentaires ajoutées à un (ou des) processeur(s) et/ou circuit(s) imprimé(s) existant(s), ou plus généralement à tout support adapté pour implémenter des lignes de codes dans le récepteur. Ainsi, l'invention ne nécessite pas de revoir toute la conception du récepteur. Brève description des figures [0048] D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la des- cription qui suit, donnée à titre illustratif et non limitatif, faite en regard des figures annexées parmi lesquelles : [0049] Ifie.11 représente un exemple de récepteur GNSS selon l'invention ; [0050] [fig.21 est une vue détaillée de la fonction réalisée par un filtre linéaire anti- interférence FIR n°m de la figure 1 ; [0051] Ifie.31 est une vue détaillée d'une partie de la vue représentée en figure 1 illustrant un mode de calcul (le la matrice d'intercorrélation Rxx ; [0052] [fig.4] illustre plus en détail le calcul de la matrice d'intercorrélation Rxx représenté en figurc 3 ; [0053] [fig.51 illustre un module de détermination de la localisation d'un récepteur GNSS ; [0054] [fig.6A] représente un premier exemple de fonction théorique nominale d'auto-corrélation. [0055] [fig.613J représente un second exemple de fonction théorique nominale d'auto-corrélation. [0056] [fig.7A] représente deux fonctions d'auto-corrélation théorique et mesurée, illustrant un bon accrochage. [0057] [fig.713J illustre deux fonctions d'auto-corrélation théorique (en trait plein) et théorique décalée (en pointillés), le décalage illustrant un faux accrochage. Description détaillée de l'invention 100581 La description détaillée de l'invention est basée sur un exemple de procédé de traitement asti-interférence d'un signal radioélectrique émis par un satellite, ledit procédé comprenant les étapes suivantes : - réception, en M points d'entrée (antennes et/ou capteurs) distincts dudit signal radioélectrique ; - détermination de la corrélation entre eux d'échantillons de sienal issus de la réception, à L instants distincts et auxdits M points d'entrée, dudit signal, puis détermination de coefficients de pondération W en fonction de ladite corrélation déterminée, M et L étant des entiers positifs, M étant strictement supérieur ou égal à 1, et L étant strictement supérieur à 1 ; - obtention d'un signal traité asti-interférence (Straité) par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée et pondérés par lesdits coefficients de pondération W. 100591 Le traitement anti-interférence est effectué au niveau d'un récepteur. 100601 Un exemple de traitement anti-interférence est un filtre STAP (en anglais « Spatio 11 Temporal Adaptative Processing ») et dans ce cas SLuaté est S,,,p; ou un filtre TAP (en anglais « Temporal Adaptative Processing ») si M est égal à 1 (une seule antenne à un seul capteur) et dans ce cas Sihthé est Stap; ou encore un filtre FDAF (en anglais « Frequency Domain Adaptive Filtering ») et dans ce cas ,S,''té est Sfd,r. [0061] Pour rappel, dans un système de positionnement par satellites, plusieurs satellites transmettent des signaux radioélectriques (RF). Un récepteur, positionné sur la position à localiser, reçoit les signaux satellites à l'aide de M points d'entrée (antennes ou capteurs), et les numérise. Puis le récepteur peut effectuer un traitement antiinterférence sur ces signaux reçus, par exemple à l'aide d'une parmi les méthodes STAP, TAP ou FDAF selon laquelle le signal traité est une combinaison linéaire des signaux reçus sur les M capteurs et de leurs L versions retardées, pondérés par des coefficients de pondérations W calculés pour rejeter l'interférence tout en conservant les signaux satellite utiles. Le filtre linéaire est appliqué avec des coefficients W remis à jour périodiquement selon une période de rafraîchissement TW en fonction de l'interférence rencontrée et caractérisée par exemple par sa matrice d' intercorrélation calculée sur le signal reçu. Ainsi, le traitement de filtrage nécessite le calcul préalable d'une matrice d'intercorrélation, nommée Rn, pour caractériser les interférences à éliminer. [0062] En aval du traitement anti-interférence, le récepteur effectue un traitement en deux phases successives du signal traité anti-interférence, une phase d'acquisition et une phase de poursuite. La phase d'acquisition a pour but de synchroniser un code et une porteuse générés dans le récepteur avec le code et la porteuse du signal reçu du satellite et prétraité comme illustré plus après. Le critère de synchronisation correspond à la maximisation de la corrélation des deux signaux, c'est-à-dire à la recherche d'un pic dans des résultats de corrélation entre le signal reçu et le signal local, qui comme indiqué plus avant est également nommée « accrochage ». [0063] Le but de l'invention est de détecter un risque de faux accrochage. Une détection de faux accrochage constitue un moyen de contrôle d'intégrité du traitement du signal reçu dans le récepteur. [0064] Ensuite, la phase de poursuite a pour but de maintenir la meilleure synchronisation du code local et de la phase de la porteuse locale avec le code et la phase de la porteuse du signal reçu, de manière à produire une mesure de la position du code et de la phase de la porteuse pour le calcul de positionnement du récepteur. Cette phase de poursuite consiste à fermer les boucles de poursuite de code (DLL) et de phase de porteuse (PLL). Puis le récepteur estime les distances qui le séparent des satellites à partir des temps de propagation des signaux satellites démodulés, et détermine alors la localisation recherchée par triangulation sur la base de ces distances. [0065] La figure 1 représente un exemple de récepteur 1 d'un système de navigation par 12 satellite utilisant des signaux provenant de S satellites SAT, (s = 1 à S) et comprenant une fonction de contrôle d'intégrité selon l'invention. Le récepteur 1 illustré implémente un traitement anti-interférence. Chaque satellite SAT, du système génère un signal radiofréquence à l'aide d'une porteuse sinusoïdale (fréquence de 1,2 à 1,6 GHz), modulée par un code d'étalement binaire connu propre à chaque satellite (fréquence de 1 à 10 MHz), puis émet le signal radiofréquence RF ainsi généré. Les codes d'étalement sont constitués d'un grand nombre de divisions temporelles élémentaires, appelées « chips » de code qui ont une durée égale à Tc. [0066] Le récepteur 1 comporte des moyens de réception 2 des signaux satellites, par exemple une antenne à diagramme de rayonnement contrôlé CRPA (en anglais « Controled Radiation Pattern Array ») comportant M capteurs A', m variant entre 1 et M, avec M nombre entier supérieur ou égal à 1. Dans la figure 1, on a représenté 3 capteurs (M est donc égal à 3). [0067] Le récepteur 1 comporte en outre un module 3 de prétraitement comportant M voies, un module 4 de traitement anti-interférence et un module 5 de détermination de la localisation du récepteur qui comprend des moyens pour la démodulation du signal reçu, prétraité et traité anti-interférence. [0068] Chacun des M capteurs est adapté pour capter les signaux radiofréquences RF présents dans sa bande de fréquences de fonctionnement, puis pour délivrer un signal radiofréquence analogique tel que S,'(t), m = 1 à M, au module 3 de prétraitement. [0069] Le module 3 de prétraitement comprend M voies parallèles correspondant chacune à un capteur A. Chaque voie m du module 3 de prétraitement comporte au moins un filtre analogique FRF, , un convertisseur analogique-numérique CAN,' réalisant un échantillonnage et une numérisation du signal S,'(t), à la fréquence d'échantillonnage Fe (la période d'échantillonnage Te étant l'inverse de la fréquence d'échantillonnage Fe), et un module Fm de passage en bande de base dans lequel le signal reçu numérisé réel est multiplié par un signal numérique complexe de type exp( 2jot Fe/4 tk ) et délivre en sortie de la voie n°m, au module de traitement anti-interférence 4, les échantillons S,'(k) d'un signal numérique complexe, où k entier indique l'indice d'un instant tk d'échantillonnage, m est l'indice de capteur (m = 1 à M) et tk est égal à : [0070] [Math.2] th = où t0 est le début de l'échantillonnage. [0071] Le module de traitement anti-interférence 4 comporte un bloc principal 40 de rejet des interférences et un bloc secondaire 41 de détermination de coefficients de pondération utilisés par le bloc principal 40. [0072] Le bloc principal 40 de rejet des interférences est adapté pour effectuer le traitement 13 de filtrage linéaire STAP. Il comporte M voies parallèles, opérant à la fréquence Fe. Le signal Sm(k) y est reçu en entrée de la même voie m, puis est fourni à l'entrée d'un filtre numérique à réponse impulsionnelle finie (dénommé « FIR n°m ») qui est illustré plus en détail à la figure 2. Le FIR n°m comporte L-1 éléments à retards unitaires complexes réalisés par l'opérateur Z-1. La sortie complexe de chaque élément à retard d'indice I est multipliée par un coefficient de pondération complexe Le signal en sortie du filtre FIR n°m est égal à la somme des sorties complexes des L-1 éléments multipliées par les coefficients de pondération complexes. [0073] Ainsi le FIR n°m délivre en sortie un signal : [0074] [Math.3] S ( k ) = o 1 W ( k - ) - ". 100751 Les signaux S'm(k) sont ensuite sommés pour m compris entre 1 et M dans le bloc principal 40. Cela équivaut à la formule : 100761 IMath.41 S ( k ) = (-I) Stap I=OiL-1Wrel,/:"Sm(k Où m est l'indice des capteurs (compris entre 1 et M), 1 est l'indice des retards (compris entre 0 et L-1), k est l'indice de l'instant d'échantillonnage. [0077] Le bloc secondaire 41 de détermination de coefficients de pondération est adapté pour déterminer les valeurs actualisées (selon la période de rafraîchissement Tw) des coefficients de pondérations utilisés par le bloc principal 40. [0078] Dans le mode de réalisation considéré, le bloc secondaire 41 comprend un sous-bloc 42 d'intercorrélation et un sous-bloc 43 de calcul des coefficients de pondération. [0079] Le sous-bloc 42 d'intercorrélation est adapté pour calculer la corrélation entre les signaux issus des M différents capteurs, et/ou leurs versions retardées de 0 à L-1 échantillons, par exemple via le calcul de la matrice d'intercorrélation 100801 La matrice R'' peut être représentée sous forme d'une matrice à coefficients ma- triciels : 100811 14 [Math.5] R;kx= [iflT] [R7,2] [Ro] [RIT] [R>211 [R1] [Ro} [RIT] [Ri [Ri] [Ro] [R1,2] [RI d [Ri] L.M u LAI 100821 R' a ses coefficients matriciels identiques sur chaque diagonale (matrice de Toeplitz). En outre, Rxx est symétrique, au sens hermitien (i.e. RT est égal à W), de même que Ro. 100831 Pour 1 variant de 0 à L-1, on a: 100841 IMath.61 1,1., 1,1,1 1, m,/ "' r m, 1,1 r rzi, rtI I rM. 1. 1 r 1, ° - - r sr, ! m, r 2,4 I R1 M 100851 Les coefficients rm.o.j de la matrice I R11 sont obtenus par la formule : 100861 IMath.71 I rit Te r114' rm n J m, n,l(t) dt = S'(k. -1) où S* est le conjugué de S. p est le produit du signal S(t) avec une version décalée du signal comme illustré plus après en relation avec les figures 3 et 4, m et n correspondent à des capteurs (entre 1 et M), 1 est un retard du signal (entre 0 et L-1), Te est la fréquence échantillonnage de numérisation du signal S, k est l'indice d'instant d'échantillonnage, Tw est la période de rafraîchissement des coefficients W, i est un indice de l'intervalle de renouvellement des coefficients W. [0087] La figure 3 illustre un mode de réalisation particulier dans lequel les M signaux issus des capteurs fournis en entrée du sous-bloc secondaire 41, échantillonnés à la fréquence Fe, sont représentés par le vecteur Z(k), de dimension M, avec : [0088] [Math.8] Z(k) = [Si(k), S Sm(k)] [0089] Le vecteur de signaux Z(k) est fourni à l'entrée du sous-bloc 42 d'intercorrélation. Dans le sous-bloc 42 d'intercorrélation, le vecteur de signaux Z(k) est fourni à une ligne à retard comprenant L-1 éléments Z-1 disposés en série. Le vecteur de signaux Z(k-1) résultant de chaque 1 retard unitaire (1 variant entre 0 et L-1) est fourni en entrée d'un ensemble de multiplicateurs complexe dans lequel chacune de ses coordonnées complexe Sm(k-1) (m variant entre 1 et M), est multipliée par le conjugué de chaque coordonnée complexe S'(k) de Z(k), où m et n correspondent à des capteurs et sont compris entre 1 et M. Le résultat de ces multiplications est une matrice pi de taille MxM dont les éléments pna, sont obtenus par la formule suivante : [0090] [Math.9] ,i(k _Te) = Avec n et m compris entre 1 et M. [0091] La figure 4 illustre plus en détail l'opérateur de corrélation spatiale permettant l'obtention des termes de la matrice p,. [0092] En référence à la figure 3, dans le sous-bloc 42 d'intercorrélation, la matrice p, ainsi définie est délivrée à un intégrateur TNT qui est adapté pour réaliser sur chaque terme de la matrice une intégration sur un intervalle de temps de durée Tw. L'intégrateur délivre ainsi le résultat de l'intégration de la matrice p, qui permet d'obtenir la matrice [R], sous-élément de la matrice d'intercorrélation R. par la formule suivante : [0093] [Math.10] (t)dt [0094] A partir de la matrice Rxx, on calcule ensuite les coefficients de pondération Wm,i à la fin de chaque période de rafraîchissement T. 100951 Le sous-bloc 43 de calcul des coefficients de pondération est adapté pour calculer les coefficients de pondération actualisés W',., en fonction de la matrice R, actualisée et fournir ces coefficients de pondérations au bloc principal 40. 100961 Dans un mode de réalisation, la méthode de Capon est utilisée pour déterminer les coefficients de pondération W sur la base de l'équation suivante : 100971 IMath.11l RY 1xC PI'- C RxxC 16 Où C est le vecteur contrainte et W est le vecteur des coefficients de pondération, cc dernier pouvant être exprimé sous forme de vecteur de la manière suivante : [0098] [Math.12] - [Se [0099] Et où Ww, peut également être exprimé sous forme de vecteur de la manière suivante : [0100] [Math.13] V1(n [0101] Le vecteur contrainte C est un vecteur déterminé pour chaque satellite (il est déterminé de manière à conserver la puissance des signaux utiles en provenance des satellites). Il peut s'exprimer par le produit de Kronecker suivant : [0102] [Math.14] C Temporel (2) Cspaüal [0103] [Math.15] Cti -1)/2 Ctewpore C tL [0104] [Math.16] [ C S ] 1 C spath Cs 1f 2%1 x 1 101051 IMath.171 Lx1 Ectpcsi [ctpcsi ctL.CSd CtL.CSO L.Mxl [0106] Les coefficients du vecteur contrainte temporel Ciemwiei peuvent être obtenus à partir de la fonction de transfert HERE du filtre analogique F. soit : [0107] 17 [Math.18] + Fe / 2 -(2/Tr.(1- ici-r.- if) Ce/ = 1_, ,)H F RF( f ),E' Lof où 10 est égal à (L+1)/2 et où f est la fréquence. [0108] Il n'est pas nécessaire d'avoir une modélisation fine de la fonction de transfert des filtres analogiques FRF',. Un gabarit valant 1 ou 0 selon que l'on est dans la bande passante du filtre ou pas est suffisant : Rut,. RF(f) est égal à 1 si f appartient à la bande passante du filtre, H filtre RF(f) est égal à 0 sinon. 101091 Les coefficients du vecteur contrainte spatial Cspaual sont obtenus à partir des dé- phasages Off, entre les capteurs pour la direction satellite choisie, soit : 101101 IMath.191 101111 On notera que lorsque l'antenne n'a qu'un seul capteur (M est égal à 1), le traitement STAP selon l'invention se réduit à un traitement TAP (Temporal Adaptive Processing). Toutes les formules données plus avant sont donc les mêmes, avec M égal à 1. En particulier dans ce cas, le vecteur contrainte spatiale Csi,'fiai est égal à 1 et ainsi le vecteur contrainte C est égal au vecteur contrainte temporel Ctemporel. 101121 Le signal résultant du traitement antibrouillage SsLap(k) est délivré en sortie du bloc principal 40 et est fourni en entrée du module 5 de détermination de la localisation du récepteur. 101131 Le module 5 de détermination de la localisation du récepteur comporte généralement plusieurs canaux de traitement CAL' s = 1 à S, un canal correspondant à un satellite. Chaque canal de traitement CAL; est adapté pour démoduler le signal émis par le satellite SAT, à l'aide d'une boucle PLL ou FLL et d'une boucle DLL. Le décalage temporel entre le signal émis par le satellite SAT, à une date à déterminer et cc signal reçu par le récepteur à une date connue correspond à la durée recherchée de propagation du signal. Dans chaque canal CAL' une réplique du code d'étalement émis est générée localement. Une mise en phase entre le signal reçu et le signal local est réalisée. Le critère de mise en phase correspond à la maximisation de la corrélation des deux signaux, des hypothèses de décalage différentes étant considérées pour chaque calcul de corrélation. [0114] Le module 5 est en outre adapté pour déterminer la localisation du récepteur 1 en fonction des durées de propagations des signaux déterminées pour les différents satellites, par triangulation. 18 [0115] Ces opérations sont illustrées succinctement en figure 5 qui illustre un canal de traitement CAL,. Dans le mode de réalisation illustré, le canal de traitement CAL, comprend un module matériel 51 et un module logiciel 52. [0116] Le module matériel 51 reçoit le signal complexe S,,,,p(k) délivré en sortie du bloc principal 40. Il est configure pour générer le code local et la porteuse locale, générant des signaux locaux, et pour les corréler ensuite avec le signal qu'il reçoit. Il peut être cadencé à une fréquence relativement élevée, de quelques MHz, par exemple de l'ordre de 50 MHz. [0117] Le module matériel 51 comprend un intégrateur numérique de phase de porteuse locale 512 générant à partir d'un premier signal de commande approprié une phase de la porteuse locale 9,,,,,,,,,.,Ipermettant la génération d'une porteuse locale complexe à partir d'un générateur de porteuse 514. Le module matériel 51 permet en premier lieu une démodulation de la porteuse, par une multiplication du signal numérisé par la porteuse locale complexe, au moyen d'un multiplicateur M. [0118] Le module matériel 51 comprend également un intégrateur numérique de phase de code 511 générant à partir d'un deuxième signal de commande approprié une phase de code local pcorkinc', permettant la génération d'un code d'étalement local propre au satellite auquel le canal considéré est dédié à partir d'un générateur de codes locaux 513. Le module matériel 51 permet en second lieu une démodulation du code, par une multiplication du signal démodulé en porteuse issu du multiplicateur M précité par un code local géré par le générateur de codes locaux 513. [0119] Le générateur de codes locaux 513 peut générer une pluralité de codes locaux : le code local dit «ponctuel », ainsi qu'un ou plusieurs codes locaux avancés et/ou retardés par rapport au code local ponctuel. Dans l'exemple illustré par la figure 5, le générateur de codes locaux 513 génère un code local ponctuel, un code local dit retard », correspondant au code local ponctuel retardé d'un délai d par exemple égal à 0,5 chip, et un code local dit « avance », correspondant au code local ponctuel avancé du délai d. Ainsi, la démodulation du code évoquée précédemment peut être réalisée en parallèle sur différentes voies correspondant aux différents codes locaux générés, au moyen d'une pluralité de multiplicateurs de corrélation Mv disposés en parallèle. [0120] Les signaux complexes en sortie des multiplicateurs de corrélation Mv peuvent alors être intégrés de manière cohérente par un intégrateur de corrélation 519, produisant des échantillons de signaux numérisés formant les voies de corrélation. Dans l'exemple illustré par la figure 5, trois voies de corrélation sont ainsi issues de l'intégrateur de corrélation 519: une voie ponctuelle, une voie retard et une voie avance, formalisables par des signaux complexes, respectivement : Zp , ZR et ZA. [0121] Ces signaux complexes de corrélation peuvent alors être exploités par le module logiciel 52. Ainsi, l'intégrateur de corrélation 519 peut être cadencé à une fréquence 19 plus basse, correspondant à la fréquence de travail du module logiciel 52, typiquement de l'ordre de quelques dizaines de Hz, par exemple 50Hz. La période correspondante est nommée période de boucle De la sorte, l'intégrateur de corrélation 519 réalise une intégration toutes les 20 ms (Tb'aek.). [0122] Le module logiciel 52 permet de générer les signaux de commande de l'intégrateur numérique de phase de porteuse locale 512 et de l'intégrateur numérique de phase de code 511, respectivement via l'asservissement d'une boucle de code DLL et d'une boucle de phase de porteuse PLL, de manière à ce que la porteuse locale soit en phase avec la porteuse reçue, et à ce que le code local de la voie ponctuelle soit en phase avec le code reçu. [0123] La boucle de code DLL comprend un discriminateur de code 521 recevant en entrée les différents signaux complexes issus des voies de corrélation formées par les sorties de l'intégrateur de corrélation 519, et restituant en sortie un signal Deode représentatif de la différence T entre la phase du code reçu pc'dc_rocii et la phase du code local pcode_incal. Ledit signal Dcode est fourni à un correcteur de code 523, ce dernier générant des commandes de code en vitesse, à la fréquence de travail du module logiciel 52, c'est-à-dire par exemple toutes les 20 ms. Les commandes de code en vitesse sont appliquées à l'entrée d'un amplificateur de commandes de code 525 restituant les deuxièmes signaux de commande en terme de nombre de chips par période d'échantillonnage à l'attention de l'intégrateur numérique de phase de code 511. [0124] D'une manière similaire, la boucle de phase de porteuse PLL comprend un discri- minateur de porteuse 522 recevant en entrée le signal issu de la voie de corrélation ponctuelle formée par la sortie correspondante de l'intégrateur de corrélation 519, et restituant en sortie un signal D,0',,, représentatif de la différence de phase Ap entre la phase de porteuse du signal reçu et la phase de porteuse locale Ledit signal D,,,,,''' est fourni à un correcteur de porteuse 524, ce dernier générant des commandes de porteuse en vitesse, à la fréquence de travail du module logiciel 52. Les commandes de porteuse en vitesse sont appliquées à l'entrée d'un amplificateur de commandes de porteuse 526 restituant les premiers signaux de commande en terme de cycles par période d'échantillonnage à l'attention de l'intégrateur numérique de phase de porteuse locale 512. [0125] Le récepteur 1 comprend également une première fonction permettant de mettre en oeuvre une première phase de détermination d'un risque de faux accrochage du procédé selon l'invention. [0126] On a vu que les calculs de corrélation sont réalisés sur la période Tso'cie des boucles de code et de porteuse et sur la base du signal S',,(k) délivré en sortie du bloc principal 40, qui du fait de son traitement au niveau du filtre anti-interférence, est renouvelé toutes les périodes de rafraîchissement Tw des coefficients W du filtre. Ainsi, à chaque période de rafraîchissement T, le changement des coefficients W du filtre est susceptible de changer les caractéristiques fréquentielles du signal S,,,'(k). Ces changements ont lieu à une fréquence beaucoup plus rapide que le temps de réponse des boucles de poursuite, ou période des boucles(TeuuPour donner un ordre de grandeur, la période de rafraîchissement T' peut varier entre 1 microseconde (1MHz) et 1 milliseconde (1kHz) alors que la période T'uck. des boucles de poursuite peut varier entre 1 milliseconde (1 kHz) et 200 millisecondes (5Hz). [0127] Selon le mode de réalisation présenté, la première phase de détermination d'un risque de faux accrochage du procédé comprend une première étape de détermination du filtre moyen sur une période d'intégration [0128] La première fonction peut comprendre une première sous-fonction 44 adaptée pour mettre en oeuvre la première étape de détermination du filtre moyen sur une période d'intégration T',,,'.'. [0129] La période d'intégration T,,,' est égale à la période sur laquelle les filtres sont moyennes. La période d'intégration est définie par rapport à la période de boucle , voire est sensiblement égale à la période de boucle La période d'intégration T ,',_,'' peut être confondue avec la période de boucle ou les deux périodes peuvent avoir des valeurs différentes. [0130] La période d'intégration T,,,' peut varier typiquement entre 1 et 200 millisecondes, par exemple être égale à 20 millisecondes. [0131] Dans le mode illustré en figure 1, la deuxième sous-fonction 44 (ou sous-fonction de moyennisation) est mise en oeuvre dans le module 4 de traitement anti-interférence et plus précisément dans le bloc secondaire 41 dudit module. La deuxième sous-fonction 44 est couplée au sous-bloc 43 de calcul des coefficients de pondération de manière à recevoir en entrée les coefficients W (et donc pouvoir les moyenner). [0132] La première étape de détermination du filtre moyen peut comprendre la déter- mination de la fonction de transfert moyenne du filtre anti-interférence. Le calcul est facilité en partant de la réponse impulsionnelle du filtre. [0133] Dans le cas d'un filtre numérique de type FIR, à réponse impulsionnelle finie, les co- efficients W donnent directement la réponse impulsionnelle RI. Il suffit de calculer la moyenne des coefficients sur la période d'intégration T, pour obtenir la réponse im- pulsionnelle moyenne. Les coefficients moyens peuvent être déterminés à l'aide de la formule suivante : [0134] [Math.20] 14 o } e in 1,17(u = 1 [0135] La réponse impulsionnelle moyenne RI,',' est donnée par : 21 [0136] [Math.21] RI = W in y e 7t moyen 101371 Dans le cas d'un filtre numérique de type IIR, à réponse impulsionnelle infinie, il est nécessaire de calculer la réponse impulsionnelle pour chaque jeu de coefficients du filtre (numérateur N et dénominateur D), puis d'en faire la moyenne, soit 101381 IMath.221 Ln -1 Id ijk - 1) o 1 S sortie (k) = r, _o S' entrée (k i DPSso = [0139] [Math.23] ll N = [ Nb N1 N2 ri [0140] [Math.24] D = [DO Di D2 ... DL; 101411 On calcule la réponse impulsionnelle moyenne RI,,,oy,' à partir de la réponse impul- sionnelle instantanée RI(u) (et non plus à partir des coefficients W), soit : [0142] [Math.25] moyen PT if RI(u)RTmoyen u=1 [0143] Pour calculer S(k) on utilise l'équation du filtre, avec en entrée : [0144] [Math.26] Sen/rée = [ 1 0 ... 0 0] 101451 Et pour conditions initiales : 101461 IMath.271 (q ) =0 pour q REVENDICATIONSProcédé dc contrôle d'intégrité du traitement d'un signal de radionavigation émis par un satellite (SATD d'un système de positionnement par satellites, ledit signal étant reçu par un récepteur (1) comportant des moyens de réception (2) et des moyens dc traitement (3, 4, 5), lesdits moyens de traitement comprenant un filtre linéaire anti-interférence (4), ledit procédé de contrôle d'intégrité étant caractérisé en ce qu'il comporte au moins une première phase de détection d'un risque de faux accrochage comprenant les étapes suivantes : - une étape de récupération d'une fonction d'auto-corrélation théorique nominale (AD) du signal reçu non traité par le filtre linéaire antiinterférence ; - une étape de détermination d'une fonction d'auto-corrélation théorique moyenne (A) dudit signal reçu et traité par le filtre linéaire antiinterférence sur une période d'intégration (T'','') définie ; - une étape de détermination du nombre de maximum locaux du module ou du module au carré de la fonction d'auto-corrélation théorique moyenne, un risque de faux accrochage étant détecté si ledit nombre de maximum locaux est supérieur ou égal à deux. Procédé selon la revendication 1, caractérisé en ce qu'il comprend une étape de détermination d'un filtre anti-interférence moyen sur la période d'intégration (T'',,,,), préalablement à l'étape de détermination de la fonction d'auto-corrélation théorique moyenne (A), et en ce que ladite fonction d'auto-corrélation théorique moyenne (A) est déterminée à partir de la fonction d'auto-corrélation théorique nominale (A0) et dudit filtre anti-interférence moyen. Procédé selon la revendication 2, caractérisé en ce que le filtre antiinterférence moyen sur la période d'intégration (T,,,0 1) est déterminé à partir des différentes valeurs des coefficients du filtre anti-interférence calculés sur ladite période d'intégration. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'il comprend en outre : - une étape de définition d'un premier seuil (S1), ledit premier seuil étant fonction de l'écart entre le rapport signal sur bruit maximum et un seuil de poursuite ; et en ce que le nombre de maximum locaux déterminés correspond au nombre (le maximum locaux pour lesquels la valeur du module ou du 31 [Revendication 5] [Revendication 6] [Revendication 7] [Revendication 8] module au carré de la fonction d'auto-corrélation moyenne finale est supérieure ou égale audit premier seuil (S1). Procédé selon la revendication 4, caractérisé en ce que ledit premier seuil (S1) est défini par défini la formule Math40 : Max ( II IllaX (Al,) II 2)où x représente l'écart entre le S1 = 10 ) rapport signal sur bruit maximum et un seuil de poursuite et Ao la fonction d'auto-corrélation théorique nominale. Procédé selon l'une des revendications 1 à 5, caractérise en ce qu'il comprend en outre une étape de déclenchement d'alarme si le nombre de maximum locaux détermines est supérieur ou égal à deux. Procédé selon l'une des revendications 1 à 6, caractérise en ce que les moyens de traitement du signal reçu comprennent un canal de traitement par satellite, chaque canal comporte une boucle de code apte d'une part à générer localement un code d'étalement caractéristique du satellite (SAT) associe audit canal de manière à former un signal local et apte d'autre part à asservir ledit signal local avec ledit signal reçu en utilisant une voie de corrélation entre ledit signal reçu et ledit signal local, et en ce que la période d'intégration (T,,,,') est fonction de la période (Tho',,,) de ladite boucle de code. Procédé selon l'une des revendications précédentes, caractérisé en ce que les moyens de traitement du signal reçu comprennent un canal de traitement par satellite, chaque canal comporte une boucle de code apte à générer localement des codes d'étalement caractéristiques du satellite (SATs) associe audit canal, lesdits codes comprenant au moins un code ponctuel, un code retard et/ou un code avance formant chacun un signal local, chaque signal local ponctuel, retard et/ou avance étant corrélé avec le signal reçu en utilisant respectivement une voie de corrélation ponctuelle, une voie de corrélation retard et/ou une voie de con-élation avance, et en ce que ledit procédé comporte en outre une deuxième phase de détection d'un faux accrochage, ladite deuxième phase comprenant : - une étape de calcul de vraisemblance entre la fonction d'auto-corrélation théorique moyenne (A) et une fonction d'auto-corrélation mesurée (Aimes) pour chacune de l'ensemble des voies de corrélation ponctuelles, retard et/ou avance ; - une étape de comparaison du résultat de calcul de vraisemblance avec 32 [Revendication 9] un deuxième seuil (S2), un faux accrochage étant détecté par un résultat du calcul de vraisemblance inférieur audit deuxième seuil. Procédé selon la revendication 8, caractérisé en ce que le résultat de l'étape de calcul de vraisemblance utilise un critère nominal (Critère0) de vraisemblance défini selon la formule Math41 : Critère° - NC ( - ( qj Aines ( cli) )2 Ir Li=o A(qi) Il ) iNf 0 V Ames 2) où qi est la position de la voie de corrélation n°i et NC le nombre de voies de corrélations. [Revendication 10] Procédé selon l'une des revendications 8 ou 9, caractérisé en cc qu'il comporte en outre une troisième phase de recherche du maximum local principal, ladite troisième phase effectuant une recherche du maximum local principal si le résultat du calcul de vraisemblance est inférieur au deuxième seuil (S2), ladite troisième phase comprenant les étapes suivantes : - une j-ème étape de calcul de vraisemblance entre la fonction d'auto-corrélation théorique moyenne (A) translatée d'une j-ème valeur de translation (d) par rapport à la position (qi) des voies de corrélation, et la fonction d'auto-corrélation mesurée (A'') pour chacune de l'ensemble des voies de corrélation, la j-ème valeur de translation (d) correspondant à un j-ème maximum local ; - une étape de comparaison du résultat du j-ème calcul de vraisemblance avec un troisième seuil (S3), un résultat du j-ème calcul de vraisemblance supérieur audit troisième seuil permettant de déterminer si la j-ème valeur de translation (d) correspond au maximum local principal. [Revendication 11] Procédé selon la revendication 10, caractérisé en ce que le résultat de la j-ème étape de calcul de vraisemblance utilise un j-ème critère (Critère j) de vraisemblance défini selon la formule Math42 : NU ' - + di) .A.'114) ) Critère j = N(' NU = II A(Cli + di) II 2 )/ =o II Aualc(Cli) II 2 ) où d est la j-ème valeur de translation, qi est la position de la i-ème voie de corrélation et NC le nombre de voies de corrélations. [Revendication 12] Récepteur (1) d'un système de positionnement par satellites comportant des moyens de réception (2) aptes à capter un signal émis par lesdits satellites et des moyens de traitement (3, 4, 5) couplés auxdits moyens de réception et aptes à traiter un signal reçu, lesdits moyens de traitement 33 [Revendication 13] comprenant des moyens de filtrage anti-interférence (4), ledit récepteur comprenant une première fonction (44, 53) apte à implémenter la première phase de détermination d'un risque de faux accrochage du procédé selon l'une des revendications 1 à 7. [Revendication 14] Récepteur (1) selon la revendication 12 comprenant en outre une deuxième fonction (54s) couplée à la première fonction (44, 53) et apte à implémenter la deuxième phase de détection d'un faux accrochage du procédé selon l'une des revendications 8 ou 9. Récepteur (1) selon la revendication 13 comprenant en outre une troisième fonction (553 couplée à la première fonction (44, 53) et à la deuxième fonction (54s) et apte à implémenter la troisième phase de recherche du maximum local principal du procédé selon l'une des revendications 10 ou 11.