FR 3038391 A1 20170106 FR 1556134 A 20150630 1 DOMAINE DE L'INVENTION La présente invention concerne un procédé de traitement de signaux de radionavigation émanant de satellites, et un dispositif adapté pour mettre en oeuvre un tel procédé. ETAT DE LA TECHNIQUE Dans le domaine des récepteurs GNSS (Global Navigation Satellite System), la tenue aux interférences est un point clé de la performance et de la continuité opérationnelle. Ces récepteurs mettent en oeuvre généralement une corrélation entre un signal de mesure et un signal de réplique, et des traitements de signal visant à réduire des interférences dans le signal de mesure avant qu'il soit corrélé avec le signal de réplique. Parmi les différents traitements pré-corrélatifs d'élimination d'interférences, il en existe un particulièrement simple à mettre en oeuvre tout en étant très efficace pour atténuer les interférences : l'excision en fréquence. L'excision en fréquence est un traitement appliqué à un signal spectral obtenu par application d'une transformée de Fourier à un signal temporel émanant d'un satellite et acquis par le récepteur. Le signal spectral est un signal d'amplitude en fonction de fréquences. Ce signal étant discret, il est formé par une pluralité de couples de valeurs d'amplitudes associées à des fréquences prédéfinies. Pour qualifier les valeurs de ces fréquences prédéfinies, on parle généralement de cases de fréquences. L'excision en fréquence consiste à annuler l'amplitude du signal spectral pour certaines cases de fréquences porteuses d'interférences, par exemple des cases de fréquences en lesquelles le signal spectral est d'amplitude supérieure à un seuil d'excision 25 prédéterminé. Cependant, un tel traitement d'excision en fréquence présente l'inconvénient de supprimer le signal utile dans la bande de la ou des interférences. Cela a pour effet de réduire les performances du récepteur (en acquisition et en poursuite) suite à la distorsion de la fonction d'autocorrélation. 30 D'autre part, en présence de brouillage autre que par un bruit large bande gaussien, un récepteur GNSS, du point de vue de la détection, n'est plus considéré comme optimal. Un récepteur est conçu pour détecter des signaux de type CDMA (Code Division Multiple Access) des satellites en présence de bruit thermique qui a les caractéristiques d'un bruit blanc gaussien. 3038391 2 Pour remédier à ces problèmes, il a été proposé d'ajuster l'amplitude du signal spectral sans pour autant l'annuler dans des plages fréquentielles du signal où des interférences sont présentes. Un tel traitement est connu sous le sigle FADP (pour - Amplitude Domain Processing » 5 dans le domaine Fréquentiel). On a illustré en figure 1 un traitement FADP connu appliqué à un signal spectral comprenant N cases de fréquences d'indices 1 à N. En une case de fréquence d'indice i, le signal spectral comprend une valeur de composante en phase Ii et une valeur de composante en quadrature Qi. 10 Le traitement FADP comprend les étapes suivantes : - une conversion du couple (Ii, Qi) en un couple de coordonnées polaires constitué d'un module et d'une phase ; - l'estimation d'une loi de répartition Fw, , - le calcul d'une fonction de pondération Gw(r), dépendant de la loi de 15 répartition estimée; - une multiplication du module par la fonction de pondération Gw(r), de sorte à produire un module filtré, - la production d'une paire d'amplitudes (lif, Qif) sur la base du résultat de la multiplication, et de la phase. 20 La loi de répartition utilisée pour ajuster l'amplitude du signal en la case fréquence d'indice i est déterminée de la façon suivante. K signaux spectraux ont été reçus avant le signal de référence à traiter. Chaque signal spectral comprend N échantillons, tout comme le signal à traiter. Les K valeurs d'amplitude en la case de fréquence d'indice i des K signaux spectraux précédents sont mémorisées. La loi de répartition pour le signal spectral de référence à appliquer en la fréquence d'indice i est estimée par une analyse restreinte aux K valeurs d'amplitude des K signaux spectraux précédents qui ont été mémorisées, en la même case fréquence d'indice i. Or ce procédé présente l'inconvénient de requérir des ressources numériques importantes du récepteur. Premièrement, les N traitements FADP mis en oeuvre pour les N cases de fréquences du signal de référence requièrent l'estimation de N lois de répartitions différentes, ce qui est particulièrement coûteux en charge de calcul lorsque N est grand. Deuxièmement, on constate qu'un historique de K signaux précédents est nécessaire pour procéder au calcul du signal de référence à traiter selon ce procédé : KN échantillons 3038391 3 préalables doivent en effet être mémorisés pour mettre en oeuvre un ajustement des N échantillons du signal de référence, ce qui nécessite de prévoir une mémoire tampon de taille importante dans le récepteur. Par ailleurs, l'utilisation d'un historique de K signaux acquis successivement au cours 5 d'un intervalle de temps passé de durée prolongée est inadaptée pour éliminer efficacement des interférences non stationnaires. EXPOSE DE L'INVENTION Un but de l'invention est de réduire des interférences dans un signal dans le domaine 10 fréquentiel moyennant des ressources matérielles en quantité réduite. A cet effet, il est proposé un procédé de traitement d'un signal temporel de radionavigation émanant d'un satellite, comprenant les étapes suivantes : - application d'une transformée de Fourier au signal temporel pour produire un signal spectral comprenant une pluralité d'échantillons associant chacun une amplitude 15 avec une fréquence, - estimation d'une loi de répartition d'amplitude du signal spectral, - pondération de l'amplitude d'un premier échantillon du signal spectral au moyen d'une fonction de pondération dépendant de la loi de répartition estimée, ladite fonction étant adaptée pour réduire une quantité d'interférences dans le premier 20 échantillon, et - pondération de l'amplitude d'au moins un deuxième échantillon du signal spectral au moyen d'une fonction de pondération qui dépend de la même loi de répartition estimée, pour réduire une quantité d'interférences dans le deuxième échantillon. Dans la mise en oeuvre de ce procédé, une seule et même loi de répartition a besoin 25 d'être estimée pour ajuster l'amplitude de plusieurs échantillons du signal spectral, voire tous les échantillons du signal spectral. La charge de calcul liée à l'établissement des lois de répartition est donc significativement diminuée. Les étapes du procédé peuvent être répétées pour plusieurs signaux spectraux. La loi 30 de répartition d'amplitude d'un des signaux spectraux, que l'on appelle par convention signal de référence, peut alors être estimée par une analyse d'échantillons restreinte aux échantillons dudit signal spectral de référence. Cette analyse diffère fondamentalement des procédés de l'art antérieur mentionnés en introduction par la provenance des échantillons analysés. 3038391 4 Les procédés de l'art antérieur calculent en effet non seulement des lois de répartitions différentes pour différents échantillons du signal de référence (et donc différentes fréquences), mais encore construisent chacune de ces lois de répartition au moyen d'un historique de valeurs passées se rapportant à une fréquence donnée. En 5 restreignant l'analyse de valeurs d'amplitudes au seul contenu du signal de référence pour construire une loi de répartition, aucune mémorisation de K signaux spectraux préalable n'est nécessaire, ce qui permet de diminuer considérablement la mémoire requise pour mettre en oeuvre l'ajustement des valeurs d'amplitude du signal de référence par rapport aux procédés décrits en introduction. 10 On pourrait penser de prime abord que le fait de s'abstenir de tout historique des valeurs d'amplitude pour construire la loi de répartition, et le fait d'utiliser cette loi de répartition pour l'ajustement de valeurs d'amplitudes en des fréquences différentes du signal spectral de référence, ne permet pas d'obtenir une loi de répartition efficace pour la réduire d'interférence dans le signal de référence. Or, contre toute attente, il a été 15 constaté que l'ajustement du signal de référence au moyen d'une loi de répartition établie sur la base de l'analyse échantillons restreinte aux valeurs d'amplitudes du seul signal de référence fournit des résultats satisfaisants pour l'ensemble des fréquences de ce signal de référence. Par ailleurs, restreindre l'analyse de valeurs d'amplitudes au seul contenu du signal 20 de référence pour construire la loi de répartition permet une élimination d'interférences non stationnaires. Le procédé de traitement ainsi proposé peut en outre être complété par les caractéristiques suivantes prises seules ou bien en combinaison lorsque cela est techniquement possible. 25 Le procédé peut comprendre une mémorisation de la fonction de pondération calculée, l'ajustement de l'amplitude du deuxième échantillon du signal spectral étant mis en oeuvre au moyen de la fonction de pondération mémorisée. Le procédé peut comprendre en outre un comptage d'occurrences de valeurs d'amplitudes du signal spectral au cours de l'estimation de la loi de répartition 30 d'amplitude du signal spectral. a w (k) , La fonction de pondération peut être de la forme G(k) = --ak(log(f k )) ou k est la valeur de l'amplitude à pondérer et f,' est la loi de répartition estimée. La transformée de Fourier peut être une transformée complexe de sorte que le signal spectral comprend une composante de signal en phase et une composante de signal en 3038391 5 quadrature ; un échantillon du signal spectral a alors pour amplitude un module du signal complexe calculé en coordonnées cartésiennes sur la base des deux composantes. Le procédé peut comprendre en outre un fenêtrage sélectionnant un bloc d'échantillons temporels dans le signal temporel, la transformée de Fourier étant appliquée 5 sélectivement au bloc sélectionné. L'étape de pondération étant répétée pour chaque échantillon du signal spectral, de sorte à produire une pluralité d'échantillons pondérés, le procédé peut comprendre en outre l'application d'une transformée de Fourier inverse à la pluralité d'échantillons pondérés, de sorte à produire un signal temporel de radionavigation corrigé. 10 Selon un autre aspect de l'invention, il est proposé un produit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé qui précède, lorsque ce produit programme est exécuté par un ordinateur. Selon un troisième aspect de l'invention, il est en outre proposé un dispositif de 15 traitement comprenant: - une entrée pour recevoir un signal temporel de radionavigation émanant d'un satellite, - un convertisseur configuré pour produire un signal spectral à partir du signal temporel par application d'un transformée de Fourier, le signal spectral comprenant 20 une pluralité d'échantillons associant chacun une amplitude avec une fréquence, - un module d'analyse spectrale configuré pour : o estimer une toi de répartition d'amplitude du signal spectral produit par le convertisseur, o pondérer l'amplitude d'un premier échantillon du signal spectral au moyen 25 d'une fonction de pondération qui dépend de la loi de répartition estimée, ladite fonction étant adaptée pour réduire une quantité d'interférences dans le premier échantillon, o pondérer l'amplitude d'au moins un deuxième échantillon du signal spectral reçu au moyen d'une fonction de pondération qui dépend de la même loi de 30 répartition, pour réduire une quantité d'interférences dans le deuxième échantillon. Le dispositif de traitement peut comprendre en outre un convertisseur inverse configuré pour produire un signal temporel corrigé par application d'une transformée de Fourier inverse aux échantillons dont l'amplitude a été pondérée par le module d'analyse 3038391 6 spectrale, et un corrélateur configure pour corréler le signal temporel corrigé avec au moins un signal de réplique prédéterminé. Enfin, il est également proposé un récepteur de signaux de radionavigation émanant de satellites, comprenant un dispositif selon le troisième aspect de l'invention, configure 5 pour traiter les signaux reçus par le récepteur. DESCRIPTION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en 10 regard des dessins annexés sur lesquels : - La figure 1, déjà discutée, est un diagramme d'étapes d'un procédé de traitement dans le domaine fréquentiel connu de l'état de la technique ; - La figure 2 représente schématiquement un récepteur de signaux de radionavigation selon un mode de réalisation de l'invention ; 15 - La figure 3 est un diagramme d'étapes d'un procédé de traitement d'un signal de radionavigation émanant de satellites, selon un mode de réalisation de l'invention ; - La figure 4 illustre deux signaux spectraux, l'un traité par un procédé d'excision en en fréquence connu de l'état de la technique, et l'autre traité via un procédé de traitement selon un mode de réalisation de l'invention. 20 Sur l'ensemble des figures, les éléments similaires portent des références identiques. DESCRIPTION DETAILLEE DE L'INVENTION En référence à la figure 2, un récepteur R de signaux de radionavigation comprend un capteur de signaux C et un dispositif de traitement de signaux 1. 25 Le capteur C est typiquement adapté pour acquérir des signaux de radionavigation temporels émis par un satellite S, par exemple de type GPS, GLONASS ou GNSS. Le capteur C est un capteur numérique, c'est-à-dire qu'il comprend des moyens de conversion analogique-numérique permettant de fournir un signal de radionavigation sous forme numérique au dispositif de traitement 1. 30 Le dispositif de traitement 1 comprend une entrée 2, une mémoire tampon 4, un premier convertisseur 6, un module d'analyse spectrale 8 et un deuxième convertisseur 10. La mémoire tampon 2 est adaptée pour mémoriser de façon temporaire des données numériques reçues via l'entrée 2 ou produites par le module 8, le convertisseur 6 et/ou le convertisseur 10. 35 La mémoire tampon 4 est une mémoire vive, par exemple de type RAM. 3038391 7 Le convertisseur 6 est configure pour convertir un signal temporel reçu par l'entrée 2 dans le domaine spectral par application d'une transformée de Fourier, et pour fournir le signal spectral ainsi converti au module d'analyse spectrale 8. Le module d'analyse spectrale 8 est configure pour traiter le signal spectral produit 5 par le convertisseur 6. Le module d'analyse 8 est prévu pour fournir des données traitées au convertisseur 10. Le convertisseur 10 est configure pour convertir un signal spectral en un signal temporel par application d'une transformée de Fourier inverse. 10 Le convertisseur 6 et/ou le convertisseur 10 et/ou le module 8 peuvent comprendre un ou plusieurs processeurs. Le dispositif de traitement 1 comprend en outre une mémoire 12 et un corrélateur 14. La mémoire 12 mémorise au moins un signal de réplique. La mémoire 12 est par exemple une mémoire morte d'un des types suivants : disque 15 dur, flash, EEPROM, etc. Le corrélateur 14 est agencé en sortie du convertisseur 10. Le corrélateur 14 est configure pour corréler un signal produit par le convertisseur inverse 10 avec un signal de réplique mémorisé par la mémoire 12. Le fonctionnement du récepteur R va maintenant être décrit. 20 Le capteur C acquiert un signal temporel émis par le satellite S. Le signal acquis est numérisé par le capteur C et fourni au dispositif de traitement 1 via son entrée 2. Le signal temporel numérique transporte plusieurs informations : - un signal utile émis par le satellite S, - un bruit thermique dans lequel le signal utile est noyé, 25 - des interférences indésirables, que l'on souhaite éliminer. En référence à la figure 3, le dispositif de traitement 1 met en oeuvre divers traitements sur le signal temporel acquis par le capteur C, dans le but d'atténuer les interférences présentes dans ce signal. Le signal numérique reçu via l'entrée 2 est mémorisé dans la mémoire tampon 4. 30 Dans une étape 104, le convertisseur 6 applique une transformée de Fourier (TFD) au signal mémorisé. La transformée de Fourier est par exemple une FFT (- Fast Fourier Transform »). La transformation 104 est de préférence mise en oeuvre sur un bloc d'échantillons du signal mémorisé dans la mémoire tampon 4 de taille prédéterminée, le bloc étant 35 préalablement sélectionné par un fenêtrage 102. 3038391 8 Le fenêtrage 102 a pour effet de limiter l'étalement spectral du signal spectral produit par la transformation de Fourier 104. Le fenêtrage 102 est par exemple du type Blackman-Harris. Le signal obtenu par l'application de la transformée de Fourier 104 est un signal 5 spectral comprenant N échantillons d'indices 1 à N. N est par exemple un multiple de 2. On a de préférence 256 REVENDICATIONS1. Procédé de traitement d'un signal temporel de radionavigation émanant d'un satellite, comprenant les étapes suivantes : - application (104) d'une transformée de Fourier au signal temporel pour produire un signal spectral comprenant une pluralité d'échantillons associant chacun une amplitude avec une fréquence, - estimation (108) d'une loi de répartition d'amplitude du signal spectral, - pondération (112) de l'amplitude d'un premier échantillon du signal spectral au moyen d'une fonction de pondération dépendant de la loi de répartition estimée, ladite fonction étant adaptée pour réduire une quantité d'interférences dans le premier échantillon, caractérisé en ce qu'il comprend en outre une étape de: - pondération (112) de l'amplitude d'au moins un deuxième échantillon du signal spectral au moyen d'une fonction de pondération qui dépend de la même loi de répartition estimée, pour réduire une quantité d'interférences dans le deuxième échantillon. 2. Procédé selon la revendication précédente, comprenant la pondération de l'amplitude de chaque échantillon du signal spectral au moyen d'une fonction de pondération qui dépend de la même loi de répartition estimée. 3. Procédé selon la revendication précédente, dans lequel les étapes d'estimation (108) et de pondération (112) sont répétées pour plusieurs signaux spectraux, et dans lequel la loi de répartition d'amplitude d'un des signaux spectraux, dit signal de référence, est estimée par une analyse d'échantillons restreinte aux échantillons dudit signal spectral de référence. 4. Procédé selon l'une des revendications précédentes, comprenant la mémorisation de la fonction de pondération calculée, l'ajustement de l'amplitude du deuxième échantillon du signal spectral étant mis en oeuvre au moyen de la fonction de pondération mémorisée. 5. Procédé selon l'une des revendications précédentes, dans lequel l'estimation (108) de la loi de répartition d'amplitude du signal spectral comprend un comptage d'occurrences de valeurs d'amplitudes du signal spectral. 3038391 13 6. Procédé selon l'une des revendications précédentes, dans lequel la fonction de pondération est de la forme G(k) = - -a (log (rw(k))) où k est la valeur de l'amplitude à ak k pondérer et f, est la loi de répartition estimée. 5 7. Procédé selon l'une des revendications précédentes, dans lequel : - la transformée de Fourier est une transformée complexe de sorte que le signal spectral comprend une composante de signal (I) en phase et une composante de signal (Q) en quadrature, 10 - un échantillon du signal spectral a pour amplitude un module du signal complexe calculé en coordonnées cartésiennes sur la base des deux composantes. 8. Procédé selon l'une des revendications précédentes, comprenant en outre un fenêtrage sélectionnant un bloc d'échantillons temporels dans le signal temporel, la transformée de 15 Fourier étant appliquée sélectivement au bloc sélectionné. 9. Procédé selon l'une des revendications précédentes, comprenant: - une répétition de l'étape de pondération pour chaque échantillon du signal spectral, de sorte à produire une pluralité d'échantillons pondérés, 20 - l'application d'une transformée de Fourier inverse à la pluralité d'échantillons pondérés, de sorte à produire un signal temporel de radionavigation corrigé. 10. Produit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une des revendications précédentes, 25 lorsque ce produit programme est exécuté par un ordinateur. 11. Dispositif de traitement comprenant: - une entrée (2) pour recevoir un signal temporel de radionavigation émanant d'un satellite, 30 - un convertisseur (6) configure pour produire un signal spectral à partir du signal temporel par application d'un transformée de Fourier, le signal spectral comprenant une pluralité d'échantillons associant chacun une amplitude avec une fréquence, - un module d'analyse spectrale (8) configure pour : 3038391 14 o estimer une loi de répartition d'amplitude du signal spectral produit par le convertisseur, o pondérer l'amplitude d'un premier échantillon du signal spectral au moyen d'une fonction de pondération qui dépend de la loi de répartition estimée, 5 ladite fonction étant adaptée pour réduire une quantité d'interférences dans le premier échantillon, caractérisé en ce que le module d'analyse spectrale est en outre configure pour : - pondérer l'amplitude d'au moins un deuxième échantillon du signal spectral reçu au moyen d'une fonction de pondération qui dépend de la même loi de répartition, 10 pour réduire une quantité d'interférences dans le deuxième échantillon. 12. Dispositif selon la revendication précédente, dans lequel le dispositif de traitement comprend en outre: - un convertisseur inverse (10) configure pour produire un signal temporel corrigé par 15 application d'une transformée de Fourier inverse aux échantillons dont l'amplitude a été pondérée par le module d'analyse spectrale (8), - un corrélateur (14) configure pour corréler le signal temporel corrigé avec au moins un signal de réplique prédéterminé. 20 13. Récepteur (R) de signaux de radionavigation émanant de satellites, comprenant un dispositif (1) selon l'une des revendications 11 à 12 configure pour traiter les signaux reçus par le récepteur (R).