FR 3012620 A1 20150501 FR 1302521 A 20131031 Procédé de détection d'interférences dans un signal de radionavigation par satellite par détection d'une déformation de la fonction de corrélation La présente invention concerne le domaine des systèmes de radio- navigation par satellite et plus précisément les récepteurs de signaux de radio-navigation par satellite. L'invention porte sur un procédé de détection d'interférences dans un signal de radio-navigation par satellite reçu par un tel récepteur. Lorsque le signal de radio-navigation est superposé à un signal interférant bande étroite, cela a pour effet d'engendrer des erreurs de mesure qui impactent les calculs de temps et de position effectués par le récepteur et in fine entrainent une erreur de positionnement. Il est donc important de pouvoir détecter la présence d'interférence dans le signal de radio-navigation reçu afin de pouvoir apporter des traitements correctifs ou d'exclure les mesures effectuées en présence d'interférence. Les solutions connues au problème précité sont basées sur la détection d'une mesure aberrante par rapport à une moyenne ou un écart attendu. Ces solutions présentent l'inconvénient de nécessiter un temps de détection trop important entre l'instant ou l'interférence apparait et l'instant où elle est effectivement détectée. L'invention propose un procédé de détection d'interférences basé sur l'évaluation de déformations de la fonction de corrélation du signal de radio- navigation. Le procédé proposé présente notamment l'avantage de permettre une détection plus rapide pour des signaux d'interférence bande étroite. A cet effet, l'invention a pour objet un procédé de détection d'interférences dans un signal de radio-navigation par satellite caractérisé en ce qu'il comporte les étapes suivantes : - Calculer une information représentative de la pente de la fonction de corrélation dudit signal avec un code d'étalement local, ladite fonction de corrélation ayant sensiblement la forme théorique d'un triangle isocèle, ladite pente étant estimée sur l'un quelconque des cotés de même longueur dudit triangle isocèle, - Comparer ladite information à un seuil de détection configuré au moins en fonction d'un rapport signal à bruit cible Selon un aspect particulier de l'invention, ladite information représentative de la pente est calculée au moins à partir d'une première mesure du maximum de la fonction de corrélation à une première position temporelle, d'une deuxième mesure de la fonction de corrélation à une deuxième position temporelle en avance de la première position temporelle d'une durée inférieure à la durée d'un créneau du code d'étalement et d'une troisième mesure de la fonction de corrélation à une troisième position temporelle en retard de la première position temporelle d'une durée inférieure à la durée d'un créneau du code d'étalement. Selon un aspect particulier de l'invention, le décalage temporel entre la première mesure et la deuxième mesure d'une part et entre la première mesure et la troisième mesure d'autre part sont sensiblement égaux en valeur absolue. Selon un aspect particulier de l'invention, ladite information représentative de la pente est calculée à partir du ratio entre l'énergie de la première mesure et la somme des énergies respectives des deuxième et 30 troisième mesures. Selon un aspect particulier de l'invention, ladite information représentative de la pente est calculée à partir d'un nombre supérieur à trois, de mesures de la fonction de corrélation prises selon différentes abscisses temporelles, en effectuant une régression linéaire. Selon un aspect particulier de l'invention, la dite information représentative de la pente est calculée au moins à partir dudit décalage 5 temporel et de la longueur d'un créneau du code d'étalement. Selon un aspect particulier de,l'invention, ledit seuil de détection est configuré en fonction du rapport des puissances du signal et de l'interférence. L'invention a également pour objet un de réception de signaux de 10 radio-navigation par satellite comprenant des moyens configurés pour mettre en oeuvre le procédé de détection d'interférences selon l'invention. D'autres caractéristiques et avantages de la présente invention apparaeont mieux à la lecture de la description qui suit en relation aux 15 dessins annexés qui représentent : - La figure 1, un synoptique d'un dispositif de réception de signaux de radio-navigation adapté pour mettre en oeuvre le procédé selon l'invention, - La figure 2a, un diagramme illustrant la forme théorique de la 20 fonction de corrélation du signal avec un code d'étalement local en l'absence d'interférences, La figure 2b, un diagramme illustrant l'impact d'interférences bande étroite sur la forme de la fonction de corrélation, La figure 3, un organigramme décrivant les étapes de mise en 25 oeuvre du procédé de détection d'interférences selon l'invention. La figure 1 représente un synoptique d'un dispositif de réception de signaux de radio-navigation adapté pour mettre en oeuvre le procédé selon l'invention. 30 Un tel dispositif 100 est basé sur un récepteur de signaux de radio- navigation usuel qui comporte notamment les fonctionnalités suivantes. Le signal reçu S est démodulé en phase par l'intermédiaire d'un démodulateur 101 et d'une porteuse locale 102. Le signal S démodulé est ensuite transmis à un corrélateur 103 pour effectuer une corrélation du signal S avec un code d'étalement généré localement 104. Avantageusement, plusieurs corrélations sont réalisées avec des retards différents à chaque fois. Typiquement, trois corrélations peuvent être réalisées, l'une à un instant temporel donné, l'une en avance de cet instant et la dernière en retard de ce même instant. Le résultat des corrélations est fourni à un module 105 de détection d'interférences selon l'invention qui exécute un procédé décrit plus loin dans 10 le texte. Le dispositif 100 comporte également un discriminateur de code 106 et un opérateur de contrôle numérique NCO 108 pour piloter la génération du code local, en particulier l'instant temporel auquel est réalisé la corrélation principale. 15 Le dispositif 100 peut comporter également un discriminateur de phase 107 et un opérateur de contrôle numérique NCO 109 pour piloter la génération de la porteuse locale 102, autrement dit la valeur du déphasage à imprimer au signal reçu S. Le dispositif 100 selon l'invention peut comporter toute autre 20 fonctionnalité habituellement implémentée dans un récepteur de signaux de radio-navigation par satellite. Les éléments décrits sur la figure 1 peuvent également être remplacés par toute autre fonction équivalente. La figure 2a représente, sur un diagramme, l'allure de la fonction de 25 corrélation du signal avec le code d'étalement généré localement, en fonction du retard appliqué. Sur la figure 2a sont représentés schématiquement trois points de corrélation. Le premier point P correspond à une corrélation du signal en phase avec le code d'étalement local, autrement dit qui donne la valeur maximale possible de la corrélation. Le deuxième point E correspond 30 à une corrélation du signal avec le code d'étalement local avancé d'un décalage temporel T et le troisième point L correspond à une corrélation du signal avec le code d'étalement local retardé d'un décalage temporel T. Ces trois points de corrélation sont typiquement utilisés dans un récepteur de signaux de radio-navigation par satellite pour estimer le décalage temporel entre le signal reçu et l'horloge associé au code généré localement. Les décalages temporels avance et retard sont inférieurs à la durée d'un créneau du code d'étalement, plus connu sous l'acronyme anglo-saxon « chip ». En l'absence d'interférence, la forme de la fonction de corrélation correspond sensiblement à un triangle isocèle en fonction du décalage temporel (avance ou retard) entre le code local et le signal reçu. La figure 2b représente les trois points de corrélation précités lorsqu'un signal interférent I bande étroite est superposé au signal de radionavigation par satellite. Dans ce cas, l'influence du signal interférent peut être telle que les 15 points de corrélation avance E et retard L sont augmentés et le point de corrélation en phase P est atténué. La forme de la fonction de corrélation est alors modifiée par rapport à sa forme théorique. Le principe à la base de l'invention consiste à détecter une 20 déformation de la fonction de corrélation par rapport à sa forme triangulaire théorique pour en déduire la présence d'interférences dans la bande du signal reçu. Le procédé selon l'invention comprend, comme illustré sur l'organigramme de la figure 3, deux étapes principales. 25 Dans une première étape 301, on calcule, à partir des résultats de corrélation du signal reçu avec le code d'étalement généré localement, une estimée de la pente de la fonction de corrélation théorique sur l'un des deux cotés du triangle (tel que représenté à la figure 2a). Pour calculer cette estimée -y, on peut utiliser la relation suivante : 30 Y= Rp .Rp R E RL.RL 3912620 Rp est le résultat de la corrélation principale, R .Rp* étant égal à l'énergie instantanée de ce résultat. RE est le résultat de la corrélation avec un retard positif de valeur prédéterminée -c, dite corrélation avance, RE.RE* étant égal à l'énergie 5 instantanée de ce résultat. RL est le résultat de la corrélation avec un retard négatif de valeur prédéterminée -T, dite corrélation retard, RL.RL* étant égal à l'énergie instantanée de ce résultat. Avantageusement, les valeurs absolues du retard positif et du retard 10 négatif peuvent être égales de sorte à ce que les corrélations avance et retard soient symétriques. Ainsi, la valeur de l'indicateur y diminue lorsque la puissance du signal interférent est élevée. Pour améliorer la précision de l'indicateur y, les trois résultats de corrélation peuvent être moyennés sur une durée K donnée. Y -- E RE (kT).RE* (kT) + RL (kT).RL* (kT) k=0 k=0 L'indicateur y peut encore être amélioré en intégrant un facteur de 20 normalisation qui dépend du retard et de la durée rc d'un créneau du code d'étalement : Rp(kT).R; (kT) k=0 2 iRp(kT).Rp (kT) Y°./ k =0 r, K RE (kT).RE* (kT) + RL (kT).RL* (kT) k=0 k =0 25 Dans une seconde étape 302, l'indicateur y est comparé à un seuil de détection d'interférences SDET configuré en fonction d'un rapport signal â bruit donné et d'un rapport des puissances entre le signal utile et le signal interférant. Si l'indicateur y est inférieur au seuil de détection SDET alors on conclut à la présence d'interférence. L'homme du métier saura, sans difficultés, déduire des indicateurs 5 décrits ci-dessus, d'autres indicateurs équivalents qui permettent d'estimer la pente de la fonctipn de corrélation ou une information caractéristique de cette pente. En particulier, un nombre supérieur à 3 de points de corrélation peuvent être utilisés, avec plusieurs retards différents. La pente de la 10 fonction de corrélation peut être estimée à partir de cette pluralité de résultats de corrélation en effectuant une régression linéaire. De façon équivalente, il est également possible d'implémenter l'indicateur inverse 1/-y et d'en déduire la présence d'interférence lorsque cet indicateur inverse est supérieur à un seuil de détection. 15 Dans les différentes variantes de réalisation du dispositif 100 selon l'invention, les modules de calcul peuvent être agencés selon différentes architectures, en particulier chaque étape du procédé selon l'invention peut être implémentée par un module distinct ou au contraire les deux étapes 20 peuvent être regroupées au sein d'un module de calcul unique. Chacun des modules de calcul que comporte le dispositif selon l'invention peut être réalisé sous forme logicielle et/ou matérielle. Chaque module peut notamment consister en un processeur et une mémoire. Le processeur peut être un processeur générique, un processeur spécifique, un 25 circuit intégré propre à une application (connu aussi sous le nom anglais d'ASIC pour « Application-Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array »). REVENDICATIONS1. Procédé de détection d'interférences dans un signal de radio-navigation par satellite caractérisé en ce qu'il comporte les étapes suivantes : - Calculer (301) une information représentative de laM pente de la fonction de corrélation dudit signal avec un code d'étalement local, ladite fonction de corrélation ayant sensiblement la forme théorique d'un triangle isocèle, ladite pente étant estimée sur l'un quelconque des cotés de même longueur dudit triangle isocèle, - Comparer (302) ladite information à un seuil de détection configuré au moins en fonction d'un rapport signal à bruit cible 2. Procédé de détection d'interférences selon la revendication 1 dans lequel ladite information représentative de la pente est calculée au moins à partir d'une première mesure du maximum de la fonction de corrélation à une première position temporelle, d'une deuxième mesure de la fonction de corrélation à une deuxième position temporelle en avance de la première position temporelle d'une durée inférieure à la durée d'un créneau du code d'étalement et d'une troisième mesure de la fonction de corrélation à une troisième position temporelle en retard de la première position temporelle d'une durée inférieure à la durée d'un créneau du code d'étalement. 3. Procédé de détection d'interférences selon la revendication 2 dans lequel le décalage temporel entre la première mesure et la deuxième mesure d'une part et entre la première mesure et la troisième mesure d'autre part sont sensiblement égaux en valeur absolue. 4. Procédé de détection d'interférences selon la revendication 3 dans lequel ladite information représentative de la pente est calculée à partir du ratioentre l'énergie de la première mesure et la somme des énergies respectives des deuxième et troisième mesures. 5. Procédé de détection d'interférences selon l'une des revendications 2 à 4 dans lequel ladite information représentative de la pente est calculée à partir d'un nombre supérieur à trois, de mesures de la fonction de corrélation prises selon différentes abscisses temporelles, en effectuant une régression linéaire. 6. Procédé de détection d'interférences selon l'une des revendications précédentes dans lequel la dite information représentative de la pente est calculée au moins à partir dudit décalage temporel et de la longueur d'un créneau du code d'étalement. 7. Procédé de détection d'interférences selon l'une des revendications précédentes dans lequel ledit seuil de détection est configuré en fonction du rapport des puissances du signal et de l'interférence. 8. Dispositif (100) de réception de signaux de radio-navigation par satellite comprenant des moyens (105) configurés pour mettre en oeuvre le procédé de détection d'interférences selon l'une des revendications 1 à 7.