FR 3038390 A1 20170106 FR 1501416 A 20150703 1 Procédé de localisation d'une source de brouillage de signaux d'un système de navigation par satellites et système associé L'invention porte sur un procédé de localisation d'une source de 5 brouillage de signaux d'un système de navigation par satellites et sur un système associé. On entend par source de brouillage, une source de brouillage intentionnelle ou une source d'interférence. Un système de navigation par satellites, ou GNSS comme 10 acronyme de "Global Navigation Satellites System" en langue anglaise, utilise une constellation de satellites qui tournent autour de la Terre sur des orbites très précisément déterminées. Ainsi, il est possible de connaître à tout instant la position d'un satellite quelconque. Les orbites des satellites sont choisies de façon à ce qu'à toute heure, 6 à 12 satellites soient visibles 15 en tout point de la terre. Chaque satellite émet plusieurs signaux radioélectriques de type de modulation et de fréquences déterminés. Au sol ou sur un véhicule terrien, maritime ou aérien, un récepteur reçoit les signaux émis par des satellites visibles. Un récepteur de système de navigation par satellites embarqué 20 mesure la durée de propagation requise pour qu'une marque horaire transmise par un satellite lui parvienne. Les marques horaires sont codées sur des ondes porteuses par la technique de modulation de phase. Chaque satellite transmet ainsi un ensemble de codes pseudo-aléatoires qui lui sont propres. Une réplique de la séquence du code est générée par le récepteur 25 et le décalage que doit subir la réplique afin de coïncider avec le code reçu correspond à la durée de propagation du signal pour parcourir la distance satellite-récepteur. Cette durée multipliée par la vitesse de la lumière dans le milieu traversé donne une mesure de distance appelée pseudo-distance. A partir des mesures des pseudo-distances le séparant de chaque satellite 30 visible, et de la connaissance de la position des satellites, le récepteur déduit sa position précise en latitude, longitude, et en altitude dans un repère 3038390 2 terrestre par une résolution numérique voisine de la triangulation. A partir des mesures de phase (Doppler) des porteuses, et de la connaissance précise de la vitesse apparente des satellites, le récepteur calcule précisément la vitesse. Il peut également en déduire la date et l'heure précise dans le repère temporel du système de navigation par satellites. La réception des signaux satellites et la précision des mesures reste très sensible, malgré l'élargissement des codes d'étalement et l'accroissement des puissances d'émission, à la présence de sources de brouillages et d'interférences. Les services utilisant à la base les systèmes GNSS de localisation et de transfert de temps sont multiples et occupent une place importante dans l'économie et le fonctionnement infrastructures des états. A ce titre, la protection des services GNSS vis-à-vis des sources de brouillage fait partie des impératifs des états qui devront déployer des moyens de détection et de localisation sur tout leur territoire, de façon à assurer la conformité de l'environnement de réception des signaux avec les standards d'utilisation. L'invention proposée s'inscrit dans le cadre d'un système de surveillance étendu de la qualité de l'environnement de réception des signaux, incluant l'ensemble des constellations GNSS, dédié à la détection et à la localisation des sources d'interférence ou de brouillage susceptibles de dégrader la qualité des mesures de position, vitesse et temps fournies par les systèmes GNSS. L'invention s'intéresse en particulier à la capacité d'utiliser les plateformes de transports ferroviaires, aériens, et maritimes, voire les transports routiers pour relaxer le maillage d'un réseau de surveillance fixe, et concentrer la surveillance autour des voies de transports de biens et de personnes. Face aux premières démonstrations de vulnérabilités des 30 récepteurs et à la présence de brouilleurs personnels à bas coût, les moyens classiques d'analyse du spectre des signaux GNSS sont progressivement 3038390 3 déployés à proximité des infrastructures critiques comme les aéroports, les ports, les centre de communications, les centrales nucléaires, ... Cependant de tels moyens (maillage du territoire national par un réseau de stations d'analyse spectrale permettant la détection de brouilleurs 5 forts) sont relativement peu efficaces dans le cas de brouilleurs lointains à large bande, puisque le niveau de brouillage suffisant à la perturbation des signaux GNSS est bien inférieur au niveau du bruit thermique standard. Ceci rend très difficile de détecter de faibles sources de brouillage, pourtant susceptibles de perturber le fonctionnement des récepteurs, ainsi 10 que de prédire leurs impacts sur le fonctionnement des récepteurs utilisateurs situés au voisinage de l'antenne de l'analyseur. De telles solutions nécessitent un nombre important de stations pour couvrir le territoire, donc un coût élevé. En outre, la zone de détection est limitée autour de chaque station. La localisation d'une source de 15 brouillage forte est faite par croisement des niveaux mesurés sur plusieurs stations. Un but de l'invention est de détecter et localiser des sources de brouillage GNSS sur de larges zones, par exemple autour de voies de communication, à partir de récepteurs mobiles. 20 Aussi, il est proposé, selon un aspect de l'invention, un procédé de localisation d'une source de brouillage de signaux d'un système de navigation par satellites à partir de puissances détectées reçues par au moins un récepteur de signaux du système de navigation par satellites embarqué à bord d'un porteur, comprenant une estimation de la puissance 25 du bruit en sortie de chaque récepteur selon l'angle de gisement par rapport au porteur et la distance par rapport au porteur, pour chaque satellite du système de navigation, utilisant une antenne à ouverture synthétique réalisant une intégration cohérente dans la direction de chaque angle de gisement visé en utilisant le déplacement connu du porteur. Un tel procédé permet d'obtenir une cartographie des sources de brouillage à proximité de voies de communication susceptibles de perturber 3038390 4 la navigation des mobiles, en profitant du déplacement du récepteur pour réaliser une localisation spatiale de ces sources. Selon un mode de mise en oeuvre, ladite estimation de la puissance du bruit en sortie de chaque récepteur utilise une pluralité 5 d'intégrateurs cohérents et une pluralité d'intégrateurs non cohérents, de durées adaptées au temps de passage d'une source dans le faisceau de l'antenne pour différentes distances. Ainsi, il est possible de réaliser une trajectographie des sources sans avoir à déployer une antenne à traitement spatial complexe, sur la base 10 d'une antenne de réception GNSS standard. Dans un mode de mise en oeuvre, ladite estimation comprend une correction la phase de porteuse dudit signal correspondant au déplacement du porteur projeté orthogonalement dans la direction de chaque angle de gisement visé, en phase d'acquisition ou en phase de poursuite. 15 Ainsi, le déplacement du porteur est mis à profit pour réaliser une poursuite en angle et distance des cibles en exploitant leur vitesse de défilement. Selon un mode de mise en oeuvre, ladite correction comprend une détermination d'un vecteur vitesse de déplacement du porteur. 20 Ainsi, la connaissance de la vitesse de déplacement du porteur permet à la fois de reconstituer une antenne synthétique pour viser simultanément dans différentes directions permettant de surveiller l'espace environnant le porteur en déplacement ainsi que de réaliser une trajectographie des cibles une fois détectées. 25 Dans un mode de mise en oeuvre, ladite détermination d'un vecteur vitesse de déplacement du porteur comprend une synchronisation en temps pour relever la trajectoire de référence du porteur et appliquer des corrections de déplacement nécessaires pour réaliser une antenne synthétique, aux instant de réception des signaux du système de navigation 30 par satellites. 3038390 5 Ainsi, les signaux sont maintenus à phase stationnaire sur une durée compatible de la longueur apparente de l'antenne synthétique et de chacune de ses directions de visée. 5 Selon un mode de mise en oeuvre, ladite détermination d'un vecteur vitesse de déplacement du porteur utilise des données fournies par un récepteur de signaux d'un système de navigation par satellites et/ou des données fournies par une centrale inertielle. Ainsi, le système de localisation des sources de brouillage est 10 autonome et peut facilement s'adapter à différents types de porteurs. Il est également proposé, selon un autre aspect de l'invention, un système de localisation d'une source de brouillage de signaux d'un système de navigation par satellites à partir de puissances détectées reçues par au moins un récepteur de signaux du système de navigation par satellites 15 embarqué à bord d'un porteur, comprenant un calculateur configuré pour effectuer une estimation de la puissance du bruit en sortie de chaque récepteur selon l'angle de gisement par rapport au porteur et la distance par rapport au porteur, pour chaque satellite du système de navigation, par antenne à ouverture synthétique réalisant une intégration cohérente dans la 20 direction de chaque angle de gisement (visé) en utilisant le déplacement connu du porteur. Dans un mode de réalisation, le système de localisation des sources de brouillage comprend des moyens de poursuite des pistes selon les angles en gisement et trajectographie par azimétrie de ces pistes. 25 On entend par azimétrie la méthode, principalement appliquée dans les domaines radar et sonar, qui permet grâce à l'évolution des mesures d'azimut et à la connaissance de la trajectoire de l'antenne de réception aux mêmes instants, de remonter à la localisation de la source du signal, sous certaines hypothèses de déplacements simples (fixe, 30 mouvement rectiligne uniforme). 3038390 6 L'invention sera mieux comprise à l'étude de quelques modes de réalisation décrits à titre d'exemples nullement limitatifs et illustrés par les dessins annexés sur lesquels : - la figure 1 illustre schématiquement un mode de réalisation d'un 5 système selon un aspect de l'invention; et - la figure 2 illustre schématiquement le fonctionnement du système, selon un aspect de l'invention. Sur l'ensemble de figures, les éléments ayants les mêmes références sont similaires. 10 La figure 1 représente schématiquement une implémentation de la présente invention. Le traitement spatial pour la formation des voies directives est basé sur celui de l'antenne synthétique en réception, comme décrit dans le brevet français FR2962812 131. 15 Le traitement d'antenne synthétique permet de réaliser des intégrations cohérentes des signaux reçus, spatialement, par exploitation du déplacement du porteur, supposé connu. Le principe de base consiste à compenser le mouvement du porteur projeté sur les directions attendues du signal supposé stationnaire, de façon à conserver une phase stationnaire 20 pendant la durée d'intégration. Il est donc possible de définir les angles de visée de l'antenne correspondant aux directions dans lesquelles on veut réaliser la détection des signaux. La fonction du traitement spatial est de réaliser une couverture 25 complète des angles en gisement (3600) par rapport à l'axe de déplacement du porteur, grâce à des faisceaux de directivité de l'antenne synthétique se recoupant à -3dB les uns des autres. Ces faisceaux sont réalisés en parallèle et en continu. 3038390 7 Le récepteur GNSS étant accordé sur les fréquences des signaux des satellites, on ne s'intéresse qu'aux sources d'interférence susceptibles d'être reçues au voisinage de ces fréquences de réceptions. Cependant, les Doppler relatifs dépendant des satellites, le 5 traitement spatial est mené de façon adapté à chaque satellite poursuivi. Selon les capacités de calcul du récepteur GNSS, il est possible d'envisager de cumuler les sorties de détection des traitements spatiaux sur l'ensemble ou seulement sur une partie des satellites poursuivis. L'estimation de la distance de la source de brouillage est réalisée 10 par adaptation de la durée d'intégration non cohérente de la chaîne de détection dans chacune des voies de visée de l'antenne synthétique, en gisement, en fonction de la vitesse de défilement de la cible dans la voie, pour différentes hypothèses de distances. La durée d'intégration est déterminée par la vitesse du porteur, la 15 largeur et la direction en gisement de la voie, et bien sûr des hypothèses de distance de la source de brouillage. A chaque voie est donc associée une batterie d'intégrateurs non cohérents spécifique à la direction en gisement de la voie. L'intégration non cohérente (qui permet l'adaptation en distance) 20 est lui-même précédé d'une étape d'intégration cohérente, sur une durée élémentaire constante et identique pour toutes les voies, qui permet de reconstituer l'étalement apparent (taille) de l'antenne, sur la durée de l'intégration cohérente. Ce traitement d'intégration cohérente permet ainsi de fournir une 25 estimation du champ de puissance du signal en gisement, réactualisée à chaque période d'intégration cohérente. L'intégration non cohérente permet de cumuler le champ de puissance sur des durées correspondant au temps de passage de la cible sous diverses hypothèses d'angle et de distance, pour la vitesse courante du porteur, ceci sous l'hypothèse que la trajectoire 30 reste rectiligne pendant la durée d'intégration non cohérente, ou encore sous condition de réaliser l'accumulation non cohérente des puissances selon des 3038390 8 directions stables en azimut (donc après compensation de la variation de direction de déplacement du porteur). La figure 1 représente le principe d'application de l'antenne à ouverture synthétique pour la détection et localisation des sources de 5 brouillage et/ou d'interférence. L'intégration cohérente est réalisée en tenant compte du déplacement du porteur. Il est effectué une détermination 1 des valeurs attendues en retard et Doppler pour réception de signaux GNSS. Cette détermination des 10 paramètres de poursuite en fréquence et distance réalise l'entretien de la valeur du Doppler de porteuse correspondant à la vitesse du satellite, calculée à partir de la position du porteur au temps courant, associée à un "top" physique ("1PPS" en entrée). Cet entretien est réalisé pour chacun des satellites du système GNSS, à partir de la connaissance des éphémérides 15 qui permet de calculer leurs positions précises à chaque instant (fournit par 1PPS de référence) et des positions et vitesses du porteur qui peuvent être fournies par toute référence fournissant une position datée (récepteur GPS par exemple, hybridé ou non avec une inertie). Ce calcul étant spécifique à chaque satellite, il existe autant de calculs que de satellites. 20 Un calcul 2 de la correction de la phase de porteuse supplémentaire permet de compenser l'évolution de la phase liée au déplacement du porteur dans la direction de la voie de visée ; cette correction s'ajoute à la correction classique de Doppler du satellite pour reconstituer la direction des voies de directivité (partie supérieure gauche du 25 schéma : NCO porteuse). Cette évolution de phase correspond à la projection du déplacement du porteur (défini par une vitesse calculée 3 du porteur) sur la direction visée par l'antenne (il existe autant de corrections de phase que de directions visées, ceci pour chacun des satellites poursuivis). 30 Cette correction 2 de phase est calculée par l'expression suivante : 2z.f0 .(t, - tk )2, . 3038390 9 dans laquelle : f0 représente la fréquence du canal Doppler considéré (on procède de façon identique sur l'ensemble des K canaux doppler), en Hz; 5 - tk représente l'instant de début d'intégration cohérente 4 qui fournit la date d'origine de la phase, en s; - ti représente l'instant courant, en s; cÏ7 représente le vecteur unitaire de la direction du satellite vue du porteur (vecteur unitaire), en coordonnées cartésiennes ; 10 - Vp représente le vecteur vitesse du porteur, en m/s; et - C représente la vitesse de la lumière, en m/s. Cette phase est ensuite introduite et corrigée 4 par un oscillateur à commande numérique NCO de phase de porteuse qui permet de générer et entretenir un signal de correction. 15 Pour chaque satellite pris en compte, il existe autant d'oscillateurs à commande numérique 4 de phase de porteuse que de directions de voie en gisement. La phase de porteuse d'un satellite est ainsi corrigée de la phase correspondant à la direction de visée, par les oscillateurs à commande 20 numérique 4 et multiplicateurs 5. Cette correction de phase correspond à la projection de la vitesse de déplacement du porteur dans la direction du signal reçu, ceci pendant une durée T d'intégration cohérente 6 relative à une position particulière du retard du signal local en cours de test 13. La mise en oeuvre d'une corrélation entre 25 le signal interférent (brouilleur) et le code local synchronisé sur chacun des signaux reçus des satellites permet de réaliser une estimation précise de l'effet du signal perturbateur sur la chaine de réception du récepteur GNSS, dont la sensibilité aux sources de brouillage dépend à la fois du spectre du signal de brouillage et de l'étalement spectral de ce signal après la 30 convolution spectrale réalisée par la corrélation avec les codes locaux. 3038390 10 L'application des compensations de phase dans chaque direction visée est synchronisée avec l'application de l'intégration cohérente 6. La durée T de l'intégration cohérente est déterminée par le nombre de voies gisement que l'on souhaite former. Par exemple, la 5 génération de 16 voies en gisement, nécessite une durée d'intégration cohérente de 20 ms pour un véhicule se déplaçant à 35m/s, i.e. 126 km/h. La synchronisation du code local (pour chacune des directions de voie) est réalisée par commande du retard de code 13 à appliquer sur la génération 14 d'un signal du code local GNSS (modulation type BPSK) par 10 un oscillateur à commande numérique NCO code 14, et calculée à partir de l'étape 1 de détermination des valeurs attendues en retard et doppler des signaux GNSS. Cette estimation est réalisée indépendamment de celle du récepteur GNSS qui pourrait avoir été affectée par la source de brouillage. La génération 14 du code local avec le retard attendu est faite par 15 un générateur à commande numérique de code. Le signal reçu est multiplié par un multiplicateur 15 par le signal généré 14 par le générateur à commande numérique de code. L'intégration non cohérente 12 est réalisée sur chaque voie en gisement, après intégration cohérente et détection quadratique 11. La 20 puissance reçue peut être alors estimée par sommation des sorties de voie en gisement, sous l'hypothèse d'un déplacement rectiligne pendant la durée d'intégration non cohérente 12. Pour être plus robuste aux changements de trajectoire du porteur, une variante de l'intégration non cohérente 12, il est possible de réaliser une 25 interpolation de la puissance en sortie des voies de gisement, de façon à reconstituer une répartition de la puissance reçue stable en azimut, ceci en utilisant le cap du déplacement disponible à chaque période d'intégration cohérente. Cela nécessite au préalable une transformation de l'échelle des mesures en gisement (obtenues après intégration cohérente et quadration 6 30 et 11) en une échelle de mesure en azimut, ceci en utilisant le cap connu du porteur. 3038390 11 Cette détection est menée en parallèle dans toutes les directions de voies. Cette poursuite en parallèle de voies dont la directivité se recoupe à -3dB du lobe principal permet de réaliser une interpolation angulaire pour une localisation précise de l'angle d'incidence du signal, par recherche des 5 maxima locaux selon les directions en gisement ou en azimut. Après la détermination 1 des valeurs attendues en retard et Doppler pour acquisition de signaux GNSS, une commande en Doppler 7 permet de piloter la valeur de la fréquence porteuse de référence selon la valeur attendue du Doppler du signal reçu depuis chaque satellite. 10 La génération 8 d'un signal à la fréquence attendue du canal Doppler est faite par un oscillateur à commande numérique de porteuse. Le signal reçu S(t) est multiplié par un multiplicateur 9 par le signal généré 8 par l'oscillateur à commande numérique de porteuse. Une synchronisation 10 en temps de mesures permet de 15 synchroniser l'estimation 3 de trajectoire de référence du porteur sur la base de temps du récepteur par la fourniture d'une horloge commune. Le signal de sortie de la synchronisation 10 permet de déterminer la référence de trajectoire 3 du porteur, i.e. la position, la vitesse et la direction de déplacement du porteur. 20 En sortie de l'intégration non cohérente 12, on dispose de la puissance du signal reçue dans chacune des directions de voies en gisement et pour différentes hypothèses de défilement de la cible dans la voie, adaptée aux différentes durées d'intégration non cohérentes. Sur la figure 2, est représentée schématiquement le procédé de 25 localisation de sources de brouillage utilisant les sorties angle-distance du traitement de la figure 1. A partir des puissances détectées en gisement et distance par satellite k du système GNSS, la figure 2 illustre comment est déterminée la présence d'un brouilleur et sa localisation. 3038390 12 La figure 1 illustre le principe de calcul de la puissance ou énergie reçue selon en gisement et distance des durées d'intégration non cohérentes Nj qui fournit une grille d'énergies reçues Ek(di, Nj) ceci par satellite k, le code généré pour la corrélation étant spécifique à chaque satellite. 5 Sur la figure 2 est représenté en entrée les puissances détectées en gisement et distance par satellite k, correspondant à la sortie de la figure 1. Un cumul 20 des énergies Ek(di, Nj) sur l'ensemble des satellites k, en utilisant une table des énergies gisement et distance par satellite k. 10 Ensuite, sont effectuées des extractions 21 des maxima locaux en gisement et distance en utilisant une table des énergies gisement et distance globale. Ensuite, est effectuée une initialisation 22 de pistes de poursuite correspondant aux maxima locaux issus d'une même source de brouillage 15 identifiée par une détection continue et cohérente en direction sur plusieurs récurrences successives; un principe d'initialisation de pistes consiste à identifier l'existence de N détections de maxima locaux dans la même direction en gisement (par exemple, sur 3 récurrences successives dans un intervalle dépendant du rapport signal à bruit du maximum local). 20 Puis, les pistes sont entretenues et poursuivies 23; une fois la piste initialisée, la proximité des nouveaux maxima locaux (des récurrences successives) avec chacune des pistes existantes est évaluée et le nouvel événement est associé ou non, à chaque récurrence t, selon un critère de distance dépendant de l'âge (depuis le moment d'initialisation de la piste) de 25 la piste et du rapport signal à bruit du maximum local testé. Les maxima retenus sont alors intégrés dans un filtre récursif de type Alpha-Beta qui permet de filtrer le bruit de mesure en gisement associé à la piste et de prédire la position en gisement pour la récurrence suivante (t+1) afin de renouveler le test d'association des futurs évènements. Les pistes en 30 gisement sont ensuite converties en azimut de sorte à pouvoir visualiser 3038390 13 l'évolution de leur défilement de façon indépendante des variations de trajectoire du porteur. Enfin la position de l'émetteur ou des émetteurs est déterminée 24 à partir de la séquence de mesures filtrées constituant chacune des pistes 5 23; pour cela on suppose que la source est fixe et que la trajectoire du porteur est connue sur la durée d'observation des pistes. Plusieurs méthodes d'azimétrie angle-distance, du type de celles appliquées en radar ou sonar, peuvent être utilisées à cette fin. Les étapes du procédé décrit ci-dessus peuvent être effectuées 10 par un ou plusieurs processeurs programmables exécutant un programme informatique pour exécuter les fonctions de l'invention en opérant sur des données d'entrée et générant des données de sortie. Un programme informatique peut être écrit dans n'importe quelle 15 forme de langage de programmation, y compris les langages compilés ou interprétés, et le programme d'ordinateur peut être déployé dans n'importe quelle forme, y compris en tant que programme autonome ou comme un sous-programme, élément ou autre unité appropriée pour une utilisation dans un environnement informatique. Un programme d'ordinateur peut être 20 déployée pour être exécuté sur un ordinateur ou sur plusieurs ordinateurs à un seul site ou répartis sur plusieurs sites et reliées entre elles par un réseau de communication. Le mode de réalisation préféré de la présente invention a été 25 décrit. Diverses modifications peuvent être apportées sans s'écarter de l'esprit et de la portée de l'invention. Par conséquent, d'autres mises en oeuvre sont dans la portée des revendications suivantes. 30 REVENDICATIONS1. Procédé de localisation d'une source de brouillage de signaux d'un système de navigation par satellites à partir de puissances détectées reçues par au moins un récepteur de signaux du système de navigation par satellites embarqué à bord d'un porteur, comprenant une estimation (12) de la puissance du bruit en sortie de chaque récepteur selon l'angle de gisement par rapport au porteur et la distance par rapport au porteur, pour chaque satellite (k) du système de navigation, utilisant une antenne à ouverture synthétique réalisant une intégration cohérente (6) dans la direction de chaque angle de gisement en utilisant le déplacement connu du porteur. 2. Procédé selon la revendication 1, dans lequel ladite estimation (12) de la puissance du bruit en sortie de chaque récepteur utilise une pluralité d'intégrateurs cohérents et une pluralité d'intégrateurs non cohérents, de durées adaptées au temps de passage d'une source dans le faisceau de l'antenne pour différentes distances. 3. Procédé selon la revendication 2, dans lequel ladite estimation (12) comprend une correction (4) de la phase de porteuse dudit signal correspondant au déplacement du porteur projeté orthogonalement dans la direction de chaque angle de gisement, en phase d'acquisition ou en phase de poursuite. 4. Procédé selon la revendication 3, dans lequel ladite correction (4) comprend une détermination (10) d'un vecteur vitesse de déplacement du porteur. 3038390 15 5. Procédé selon la revendication 4, dans lequel ladite détermination (10) d'un vecteur vitesse de déplacement du porteur 5 comprend une synchronisation en temps pour relever la trajectoire de référence du porteur et appliquer des corrections de déplacement nécessaires pour réaliser une antenne synthétique, aux instant de réception des signaux du système de navigation par satellites. 10 6. Procédé selon la revendication 4 ou 5, dans lequel ladite détermination (10) d'un vecteur vitesse de déplacement du porteur utilise des données fournies par un récepteur de signaux d'un système de navigation par satellites et/ou des données fournies par une centrale inertielle. 15 7. Système de localisation d'une source de brouillage de signaux d'un système de navigation par satellites à partir de puissances détectées reçues par au moins un récepteur de signaux du système de navigation par satellites embarqué à bord d'un porteur, comprenant un 20 calculateur configuré pour effectuer une estimation (12) de la puissance du bruit en sortie de chaque récepteur selon l'angle de gisement par rapport au porteur et la distance par rapport au porteur, pour chaque satellite (k) du système de navigation, par antenne à ouverture synthétique réalisant une intégration cohérente (6) dans la direction de 25 chaque angle de gisement en utilisant le déplacement connu du porteur. 8. Système de localisation des sources de brouillage selon la revendication 7, comprenant des moyens de poursuite des pistes selon les angles en gisement et trajectographie par azimétrie de ces pistes.