FR 2952491 A1 20110513 FR 0905346 A 20091106 PROCEDE ET SYSTEME POUR LE BROUILLAGE DE SIGNAUX GNSS La présente invention concerne un procédé et un système pour le brouillage de signaux de type GNSS (« Global Navigation Satellite System »). Ces signaux sont utilisés dans les systèmes de navigation par satellite et permettent à un récepteur mobile de déterminer ses coordonnées géographiques. Un des systèmes de navigation par satellite connu est le système américain GPS (« Global Positioning System ») mais l'invention s'applique à tout type de systèmes GNSS, notamment le système européen en cours de développement GALILEO. Un brouilleur est un dispositif actif qui génère un signal de brouillage dont les caractéristiques permettent de perturber la réception d'un signal GNSS par un récepteur mobile. Une application possible de ce type de dispositif est de dérégler les moyens de navigation d'un véhicule. En particulier le brouillage des signaux GNSS permet de supprimer la connaissance, par le véhicule visé, de ses coordonnées géographiques. Un brouilleur se caractérise le plus souvent par une antenne active dont la puissance est réglable et constitue un paramètre important. La puissance émise doit être suffisamment importante pour perturber le signal visé même si le récepteur associé est à une distance importante du brouilleur. A contrario, la puissance émise ne doit pas être trop importante pour ne pas impacter l'autonomie du dispositif de brouillage, ni générer trop d'effets collatéraux. Un brouilleur évolué peut également générer un signal de brouillage adapté à une forme d'onde donnée de façon à prendre en compte au mieux les caractéristiques intrinsèques du signal que l'on souhaite perturber. Considérons un espace géographique donné dans lequel se déplacent un premier groupe de récepteurs GNSS dont on souhaite perturber la réception d'un signal de positionnement par satellite tout en préservant un second groupe de récepteurs GNSS que l'on veut éviter de brouiller. De plus la puissance totale du ou des brouilleurs utilisés doit être optimale, car il ne faut pas de déperdition d'énergie pour économiser la puissance du ou des brouilleurs dans le but de préserver leur autonomie. Les solutions de l'art antérieur ne permettent pas de déterminer de façon optimale le niveau de puissance d'un ensemble de brouilleurs utilisés dans le but de perturber un premier groupe de récepteurs GNSS tout en préservant un second groupe de récepteurs situés dans la même zone géographique. Au contraire, le signal de brouillage généré interfère sur tous les signaux qui utilisent des fréquences proches du signal à perturber. Les solutions actuelles conduisent donc à appliquer un brouillage vers une zone géographique donnée et à exclure de cette zone les récepteurs que l'on ne veut pas perturber. De plus les solutions connues entrainent un surdimensionnement de la puissance de brouillage qui n'est pas déterminée de façon optimale et adaptée au récepteur visé. La présente invention propose une solution permettant de supprimer les limitations de l'art antérieur. L'effet du brouillage ainsi que ses contraintes sont modélisés selon une approche linéaire afin de déterminer une répartition optimale des puissances d'un ensemble de brouilleurs utilisés pour perturber un premier groupe de récepteurs tout en épargnant un second groupe de récepteurs, les deux groupes de récepteurs étant situés dans une zone géographique donnée. De plus la puissance globale générée par l'ensemble de brouilleurs est minimisée de façon à ne pas gaspiller de ressources inutilement. Ainsi l'autonomie globale du système de brouillage est améliorée. A cet effet l'invention a pour objet un procédé de brouillage, par un nombre K au moins égal à un, de brouilleurs, des signaux de type GNSS 2952491 reçus par un premier groupe, comportant un nombre N au moins égal à un, de récepteurs mobiles i de positions connues au sein d'une zone géographique donnée, lesdits récepteurs mobiles étant adaptés à recevoir au moins un signal GNSS émis selon une forme d'onde donnée pour supporter 5 un service I, lesdits brouilleurs k étant adaptés à générer au moins un signal de brouillage s de puissance Jk,s au plus égale à une puissance maximale Jk max ,ledit signal de brouillage étant adapté à brouiller ledit signal GNSS supportant le service I, ledit procédé étant caractérisé en ce qu'il comporte au moins les étapes suivantes : 10 o déterminer, pour chacun desdits récepteurs mobiles i et chacun des signaux GNSS reçus par ledit récepteur i, une marge de fonctionnement ANi,i = G` Cl ù Nrh•FR en fonction du gain G, r7.LR procuré par l'antenne dudit récepteur i sur le signal GNSS qu'il reçoit, de la puissance CI du signal GNSS reçu devant l'antenne 15 dudit récepteur i, LR un terme modélisant les pertes de puissance internes audit récepteur i, FR le facteur de bruit dudit récepteur i, Nth la puissance de bruit thermique qui affecte le signal GNSS reçu par ledit récepteur i et 'n la valeur minimum du rapport signal à bruit plus interférence en dessous duquel ledit récepteur i ne reçoit plus 20 ledit signal GNSS, o déterminer, par une étape de calcul, la puissance optimale Jk,S affectée à chaque signal de brouillage s produit par chacun desdits brouilleurs k et qui respecte les contraintes suivantes : K s ^ la puissance de brouillage totale J =EEJk,s est la plus k=1 s=1 25 faible possible, avec S le nombre de signaux de brouillage différents générés par ledit brouilleur k, ^ la contrainte de brouillage suivante est respectée pour chaque récepteur i et chaque service GNSS I utilisé par ledit récepteur i: K S d 1EGEk,i.Pk,i.SSCIk,s l . jk,s / Ak,i k=1 s=1 avec GEk,; le gain d'antenne du brouilleur k dans la direction dudit récepteur i, Pk,; le niveau de réjection du signal de brouillage par l'antenne dudit brouilleur k, Ak,; un terme modélisant les pertes de propagation entre le brouilleur k et le récepteur i, SSC rks,l le coefficient de corrélation spectral entre le signal de brouillage et le signal supportant le service I et reçu par ledit récepteur i. Dans une variante de réalisation de l'invention, un second groupe, comportant un nombre M au moins égal à 1, de récepteurs mobiles de positions connues sont situés dans ladite zone géographique donnée et les puissances Jk,s desdits signaux de brouillage émis par lesdits brouilleurs sont déterminées de façon à ce qu'ils n'affectent pas la réception correcte des services GNSS utilisés par lesdits récepteurs mobiles dudit second groupe, le procédé selon l'invention comportant les étapes supplémentaires suivantes : o déterminer ledit seuil de fonctionnement ON;,s pour chacun des récepteurs mobiles dudit second groupe, o déterminer, par une étape de calcul, la puissance optimale Jk,s affectée à chaque signal de brouillage s produit par chacun desdits brouilleurs k et qui respecte en plus les contraintes de non brouillage suivantes: ^ les inégalités suivantes sont respectées pour chaque récepteur i dudit second groupe et chaque service GNSS I utilisé par ledit récepteur i : K E E GEk,i .Pk,i .S SCIk , ,l ..%k,s / Ak,i i=N+1,...,N+M k=1 s=1 Dans une variante de réalisation de l'invention, l'étape de calcul permettant de déterminer lesdites puissances optimales de brouillage Jk,s est réalisée, si une solution existe, par l'exécution de l'algorithme du Simplexe appliqué aux contraintes suivantes : Max( Ct.x) A.x = b x >=0, où Max(Ct.x) désigne la contrainte selon laquelle le produit scalaire des vecteurs C transposé et x doit être maximum, avec x = ~E~ un vecteur de dimension égale à K.S + (N1+ M1+ K) avec J= J! un vecteur de dimension égale à K.S dont les composantes sont les puissances JI=Jk,s des dits brouilleurs, un vecteur de dimension N1+M1+K dont les composantes sont des variables spécifiques à l'algorithme du simplexe, N1 est le nombre de contraintes de brouillage à respecter, 0 +1 à (N+M).S et I variant de 1 à K.S, si au moins une condition de non 5 brouillage existe, sinon Ab est nulle, INI est la matrice identité de dimension N1xN1, IMl est la matrice identité de dimension M1xM1, IK est la matrice identité de dimension KxK, 10 sinon, ù Da b= Db avec Jmax un vecteur de dimension N1 dont les composantes sont les seuils de fonctionnement AN;,s de chaque récepteur GNSS dudit premier groupe, 15 Db= un vecteur de dimension M1 dont les composantes sont les seuils de fonctionnement ON; de chaque récepteur GNSS dudit second groupe, avec q;,i = 1 pour j variant de (i-1).S+1 à i.S et qi= 0 J = max un vecteur de dimension K dont les composantes sont les Jkmax puissances maximales de chacun des brouilleurs. L'invention a également pour objet un système de brouillage comportant une pluralité de brouilleurs adaptés à générer au moins un signal de brouillage permettant de perturber le signal reçu par au moins un récepteur mobile de type GNSS et un organe centralisé arrangés en réseau par des moyens de communication, ledit système étant caractérisé en ce que ledit organe centralisé comporte des moyens pour exécuter les étapes du procédé selon l'invention. D'autres caractéristiques apparaîtront à la lecture de la description détaillée donnée à titre d'exemple et non limitative qui suit faite en regard d'un dessin annexé qui représente : La figure 1, un schéma du système de brouillage selon l'invention mis en oeuvre sur une zone géographique dans laquelle se déplace une pluralité de récepteurs GNSS. Le procédé de brouillage selon l'invention est basé sur une modélisation linéaire globale du problème relatif au brouillage d'un nombre N de récepteurs GNSS avec un réseau comportant un nombre K de brouilleurs. L'introduction d'une modélisation linéaire permet ensuite de résoudre le problème en appliquant l'algorithme du simplexe, connu de l'Homme du métier. En outre, le procédé selon l'invention permet de prendre en compte la présence éventuelle d'un nombre M de récepteurs GNSS que le brouillage ne doit pas affecter. Plus précisément, l'invention consiste à mettre en oeuvre une modélisation linéaire de l'effet cumulé des interférences, de formaliser le problème de minimisation de la puissance totale de brouillage émise puis de déterminer une solution optimale, si elle existe, qui permette de rendre inopérant les N récepteurs GNSS à brouiller tout en assurant le fonctionnement des M autres récepteurs GNSS à préserver. Modèle linéaire d'interférence d'un récepteur GNSS Le comportement d'un récepteur peut de façon générale être défini par un seuil de fonctionnement ri et un rapport signal à bruit SNR, en anglais « Signal to Noise Ratio ». En première approximation, le récepteur ne fonctionne que si, et seulement si, le rapport signal à bruit SNR est supérieur au seuil r.. En présence d'interférence volontaire ou involontaire, on définit le rapport signal à bruit plus interférence SNIR, en anglais « Signal to Noise & Interference Ratio ». En supposant que le récepteur GNSS est idéal, et en présence d'une source d'interférence, par exemple un brouilleur, le rapport signal à bruit plus interférence en sortie de la chaîne de traitement interne au récepteur GNSS peut être exprimé par la relation suivante : GS .CS /L R SINR = N F + SSC C (1), th' R I,s. Où: GS est le gain procuré par l'antenne, ou le réseau antennaire, du récepteur sur le signal utile reçu par ce même récepteur, Cs est la puissance du signal utile devant l'antenne du récepteur, Lr est un terme qui modélise les pertes internes au récepteur, Nth est la puissance de bruit thermique, Fr est le facteur de bruit du récepteur, C, est la puissance du signal interférant devant le récepteur, SSC1,s est le coefficient de corrélation spectral entre le signal interférant I et le signal utile S. Sa valeur est comprise entre 0 et 1 et est donnée par la relation suivante : SSC,,s = ft's(f).1(f)df (2) Où: fis( f) est la densité spectrale normalisée de puissance du signal utile S (D,( f) est la densité spectrale normalisée de puissance du signal interférant Tous les termes de la relation (1) sont supposés connus, soit par mesure ou caractérisation du récepteur GNSS, soit en prenant l'hypothèse la plus pessimiste. En présence d'un bruit thermique stationnaire, le seuil de io fonctionnement du récepteur GNSS peut se traduire par une marge de tolérance sur la puissance du signal interférant tel que l'illustrent les relations suivantes : SINR > i1 p SSC, s .CI En combinant les relations (1) et (4) on obtient une nouvelle formule du rapport signal à bruit plus interférences : 15 GS.Cs /LR (5) SINR = Nth.FR + GEk .Pk .SSCI,S .Jk / Ak La relation (3) donne la puissance minimale de brouillage requise pour que le récepteur GNSS ne fonctionne plus, ou la puissance maximale pour qu'il ne 5 soit pas perturbé. Modèle linéaire du brouillage multiple d'un récepteur de service GNSS Dans le cas où un nombre K de signaux de brouillage est considéré, le niveau d'interférence résultant en sortie de la chaine de traitement du 10 récepteur est égal à la somme des niveaux d'interférence individuels, la relation (5) devient alors : Gs .Cs / LR SINR = K (6) Nth.FR +EGEk.Pk.SSC,,s.Jk /Ak k=1 La condition de fonctionnement du récepteur, modélisée par la relation (3) dans le cas d'un seul brouilleur, devient alors : K GEk .Pk .SSCI S .Jk l Ak AN (8) k=1 K La puissance globale de brouillage mise en oeuvre est égale à 4 =EJk k=1 Modélisation du problème d'optimisation Le contexte dans lequel se positionne l'invention est le suivant. Il concerne un système de brouillage fixe constitué d'un nombre K de 15 20 brouilleurs de puissances maximales données. Les brouilleurs sont de positions et d'orientations connues. Ce système de brouillage a pour objectif de perturber un nombre N de récepteurs GNSS tout en préservant le fonctionnement d'un autre groupe de récepteurs GNSS contenant M unités. Les positions et les caractéristiques de tous les récepteurs GNSS sont supposées connues. Le problème technique que vise à résoudre l'invention est formulé de la façon suivante. Existe-t-il une solution de brouillage, c'est-à-dire une configuration de brouillage qui respecte les contraintes de puissance maximum, le brouillage des N récepteurs cible et l'absence de brouillage collatéral des M récepteurs à préserver ? Si plusieurs solutions existent, quelle est celle qui minimise la puissance globale émise par le système de brouillage ? Afin de répondre au problème technique précédent, les puissances Jk des brouilleurs doivent respecter les contraintes suivantes. La puissance globale J doit être minimisée, c'est-à-dire qu'elle doit être la plus faible possible tout en respectant les autres contraintes liées au problème ce qui se traduit par la contrainte suivante : K min 1 Jk (9) \k=1 Pour chaque récepteur GNSS indexé par l'indice i et que l'on souhaite 20 brouiller, l'inégalité suivante doit être respectée : K E GEk i .Pk i .SSC1, .Jk / Ak,i >_ ONi i =1,..., N (10) k=1 De façon analogue, pour chaque récepteur GNSS d'indice i variant de N+1 à 25 N+M, que l'on souhaite préserver l'inégalité suivante doit être respectée : K EGEki•Pki.SSCIks.Jk/Aki Mise sous forme standard pour résolution par l'algorithme du simplexe Le problème d'optimisation traduit par les relations (9), (10) et (11) peut se réécrire sous la forme matricielle standard de l'algorithme du simplexe qui est une méthode de programmation linéaire connue de l'Homme du métier. Ainsi les relations précédentes s'écrivent sous la forme suivante : Max( Ct.x) (12) Sous les contraintes : A.x = b (13) x >=O (14) Où Max(Ct.x) représente le maximum du produit scalaire entre le transposé du vecteur C et le vecteur x. b est un vecteur, A est une matrice. Le lien entre les relations (9),(10),(11) et (12),(13),(14) se fait en assignant aux variables des relations (12),(13),(14) les valeurs suivantes. x = [J] est un vecteur de dimension égale à K + (N+M+K) avec J= Ji un vecteur de dimension égale à K dont les composantes sont les 10 puissances Jk de chaque brouilleur, E_ ei un vecteur de dimension N+M+K dont les composantes sont des variables spécifiques à l'algorithme du simplexe. C = [[ù1] [0] est un vecteur de dimension K +N+M+K avec 15 [-l] un vecteur de dimension K dont toutes les composantes sont égales à -1 et [0] un vecteur de dimension N+M+K dont toutes les composantes sont égales à O. La contrainte (12) est ainsi équivalente à la contrainte (9). 20 On pose ensuite A _Aa_ Ab IK qui est une matrice de IN 0 0 o IM o 0 IK 0 dimension (N+M+K) x (N+M+2K) avec : et ai,k ù GEk,i.Pk,i.SSCk i / Ak i pour I variant de 1 à N et et A,k = GEk,i.Pk,i.SSCk,i /Ak,i pour I variant de N+1 à N+M et k variant de 1 à K, IN est la matrice identité de dimension NxN, lm est la matrice identité de dimension MxM, IK est la matrice identité de dimension KxK. Enfin le vecteur b est un vecteur de dimension N+M+K donné ùDa 10 b= Db Jmax avec un vecteur de dimension N dont les composantes sont les seuils de tolérance AN; de chaque récepteur GNSS à brouiller, un vecteur de dimension M dont les composantes sont les 15 seuils de tolérance AN; de chaque récepteur GNSS à préserver, Jkm un vecteur de dimension K dont les composantes sont les puissances maximales de chacun des brouilleurs. A a a ... ai,k k variant de 1 à K, Ab = A,k " Avec les définitions précédentes, les relations (13) et (14) sont équivalentes aux relations (10) et (11) ainsi qu'à la contrainte de puissance maximale de chaque brouilleur. Une fois les relations (13) et (14) établies sous la contrainte (12), l'algorithme du simplexe, connu de l'Homme du métier peut s'appliquer afin de déterminer si le problème posé possède des solutions et lorsque c'est le cas, de déterminer une solution de brouillage de puissance globale minimale. 1 o L'algorithme du simplexe est une méthode créée par George Dantzig dans les années 50 pour résoudre les problèmes de programmation linéaire. Etant donné un ensemble d'inégalités linéaires sur n variables réelles, l'algorithme permet de trouver la solution optimale pour une fonction de coût donnée qui est elle aussi linéaire. 15 Modélisation du problème d'optimisation pour plusieurs services Dans une variante de réalisation de l'invention, les récepteurs GNSS utilisent différents services de navigation qui diffèrent par leur performance, leur qualité et leur accessibilité et qui se caractérisent par un signal de 20 puissance, de forme d'onde et de fréquence donnée. On considère par la suite un nombre S de services différents. Les brouilleurs ont également la capacité de générer un signal de brouillage adapté à chacun desdits services. On parle dans ce cas de forme d'onde de brouillage. Chaque brouilleur peut être paramétré pour émettre 25 une ou plusieurs formes d'onde de brouillage avec des puissances respectives différentes. La configuration de brouillage de chaque brouilleur représente la liste des puissances allouées à chaque forme d'onde de brouillage. La modélisation du problème précédemment décrite concerne le cas ou les 30 brouilleurs n'ont la capacité de générer qu'une seule forme d'onde de brouillage adaptée à perturber un seul service utilisé par le récepteur GNSS. 5 Dans le cas où plusieurs services sont considérés, les relations (9),(10) et (11) deviennent , en intégrant la dimension du service en supplément : (K S I ElJk,s k,s ~k 1 s / Pour chaque récepteur GNSS indexé par l'indice i et que l'on souhaite brouiller et pour chaque service indicé par la variable I et utilisé par le récepteur, l'inégalité suivante doit être respectée : (9') K S E GEk i .Pk i .SSCIk,s .Jk s' Ak i k=1 s=1 10 i =1,..., N (10') De façon analogue, pour chaque récepteur GNSS d'indice i variant de N+1 à N+M, que l'on souhaite préserver et pour chaque service I utilisé par le récepteur, l'inégalité suivante doit être respectée : K S E GEk,i .Pk,i .SSCIk,s .Jk,s / Ak,i i,1 i=N+1,...,N+M (11') k=1 s=1 Enfin la somme des puissances allouées à chaque forme d'onde de brouillage ne doit pas dépasser la puissance maximale du brouilleur ce qui 20 se traduit par la relation suivante : s 15 EJk,s - Jk max s=~k =1..K (15) La modélisation du problème selon les relations (9'), (10'), (11') et (15) conduit donc à introduire un nombre égal à N1+M1+K contraintes dont N1 contraintes de brouillage avec N1 Mise sous forme standard pour résolution par l'algorithme du simplexe dans le cas de plusieurs services La résolution du problème d'optimisation dans le cas où plusieurs services sont considérés se réécrit sous la même forme que pour le cas d'un unique service, c'est-à-dire par l'intermédiaire des relations (12),(13) et (14) avec les différences suivantes dans la définition des variables. Le vecteur J devient J = J1 avec I=k,s k variant de 1 à K et s variant de 1 à S. Ce vecteur est de dimension KxS. 1 o Le vecteur E est de dimension N 1 +M 1 +K Le vecteur C comporte KxS composantes égales à -1 et N1+M1+K composantes égales à O. La matrice A est de dimension (N 1 +M 1 +K)x(K.S +N 1 +M 1 +K), Aa - 'Ni --- - Ab 0 Q 0 15 La composante au de la matrice Aa est égale à ai I = GEk i.Pk i.SSCI,i / Ak,i avec i variant de 1 à N.S et I variant de 1 à K.S. La composante (i;,, de la matrice Ab est égale à Al = GEk,i.Pk,i.SSCI,i /Ak,i avec i variant de N.S +1 à (N+M).S et I variant de 1 à K.S. La matrice Q est égale à Q = qi,i avec q;= 1 pour j variant de (i- 20 1).S+1 à i.S et q;,i = 0 sinon. IN1 est la matrice identité de dimension N1xN1, lMi est la matrice identité de dimension M1xM1, IK est la matrice identité de dimension KxK. A= b est un vecteur identique à celui définie pour le cas d'un unique service mis à part que le vecteur Da est de dimension N1 et le vecteur Db est de dimension M1. L'algorithme du simplexe s'applique ensuite de la même façon que pour le cas où un unique service est considéré. La figure 1 illustre un exemple de mise en oeuvre du procédé selon l'invention. Un ensemble de véhicules mobiles, dans l'exemple décrit des aéronefs, sont équipés de récepteurs GNSS par l'intermédiaire duquel ils reçoivent les signaux, d'un ou plusieurs services, provenant d'une pluralité de satellites 101,102. Le procédé selon l'invention a pour objectif de dénier les services GNSS d'un premier groupe de récepteurs mobiles 103,104,105,c'est à dire de brouiller la réception des signaux supportant ces services et provenant des satellites 101,102, tout en préservant les services GNSS d'un second groupe de récepteurs mobiles 106,107,108, c'est à dire en garantissant une réception correcte des signaux supportant ces services et provenant desdits satellites 101,103 et à destination. Un système de brouillage 100 selon l'invention est composé d'une pluralité de brouilleurs 109,110,111, disposés en réseau au moyen de liens de communications 112 qui leur permet de communiquer entre eux et vers un organe centralisé 113. Cet organe 113 comporte au moins une fonction d'observation qui a pour objectif de déterminer ou de recevoir depuis un dispositif externe au système les positions des différents véhicules mobiles 103,104,105,106,107,108 situés dans une zone géographique donnée. A partir de ces positions, des caractéristiques des récepteurs GNSS embarqués dans lesdits véhicules mobiles et des caractéristiques des brouilleurs 109,110,111, l'organe centralisé 113 exécute le procédé selon l'invention afin de déterminer la répartition des puissances de brouillage sur l'ensemble des brouilleurs 109,110,111 disponibles, ainsi que l'allocation de puissance affectée à chaque forme d'onde de brouillage exploitée par chacun desdits brouilleurs. Une fois la liste des puissances déterminées, l'organe centralisé 113 communique cette liste aux brouilleurs 109,110,111, par l'intermédiaire des moyens de communication 112. Si aucune solution n'est trouvée à l'issu de la mise en oeuvre du procédé selon l'invention, alors l'organe centralisé 113 communique aux brouilleurs 109,110,111 une valeur par défaut de la puissance qu'ils doivent utiliser. Dans une variante de réalisation de l'invention, si aucune solution n'est trouvée ou si les besoins spécifiés ne nécessitent pas d'épargner certains récepteurs de l'effet du brouillage, les contraintes qui consistent à io préserver un nombre M de récepteurs GNSS peuvent être relâchées, ce qui revient à supprimer de la formulation du problème les inéquations (11) ou (11'). L'algorithme du simplexe est exécuté à partir des nouvelles contraintes relâchées afin de déterminer une liste de puissances qui résolvent le problème posé. La suppression des inéquations (11) et (11') entraine alors 15 une valeur nulle du nombre M ou M1 de contraintes de non brouillage. REVENDICATIONS1. Procédé de brouillage, par un nombre K au moins égal à un, de brouilleurs, des signaux de type GNSS reçus par un premier groupe, comportant un nombre N au moins égal à un, de récepteurs mobiles i de positions connues au sein d'une zone géographique donnée, lesdits récepteurs mobiles étant adaptés à recevoir au moins un signal GNSS émis selon une forme d'onde donnée pour supporter un service I, lesdits 1 o brouilleurs k étant adaptés à générer au moins un signal de brouillage s de puissance Jk,s au plus égale à une puissance maximale Jk max ,ledit signal de brouillage étant adapté à brouiller ledit signal GNSS supportant le service I, ledit procédé étant caractérisé en ce qu'il comporte au moins les étapes suivantes : 15 o déterminer, pour chacun desdits récepteurs mobiles i et chacun des signaux GNSS reçus par ledit récepteur i, une marge de fonctionnement ~Ni,c = 1. ` ù Nth•FR , en fonction du gain G, procuré par l'antenne dudit récepteur i sur le signal GNSS qu'il reçoit, de la puissance C, du signal GNSS reçu devant l'antenne 20 dudit récepteur i, LR un terme modélisant les pertes de puissance internes audit récepteur i, FR le facteur de bruit dudit récepteur i, Nth la puissance de bruit thermique qui affecte le signal GNSS reçu par ledit récepteur i et rt la valeur minimum du rapport signal à bruit plus interférence en dessous duquel ledit récepteur i ne reçoit plus 25 ledit signal GNSS, o déterminer, par une étape de calcul, la puissance optimale Jk,s affectée à chaque signal de brouillage s produit par chacun desdits brouilleurs k et qui respecte les contraintes suivantes :K S ^ la puissance de brouillage totale J =EEJk,s est la plus k=1 s=1 faible possible, avec S le nombre de signaux de brouillage différents générés par ledit brouilleur k, ^ la contrainte de brouillage suivante est respectée pour chaque récepteur i et chaque service GNSS I utilisé par ledit récepteur i: K S GEk,i.Pk,i.SSCIk s l.Jks / =1,...,N k=1 s=1 avec GEk,; le gain d'antenne du brouilleur k dans la direction dudit récepteur i, Pk,; le niveau de réjection du signal de brouillage par l'antenne dudit brouilleur k, Ak,; un terme modélisant les pertes de propagation entre le brouilleur k et le récepteur i, SSC Ik ,l le coefficient de corrélation spectral entre le signal de brouillage et le signal supportant le service I et reçu par ledit récepteur i. 15 2. Procédé selon la revendication 1 caractérisé en ce qu'un second groupe, comportant un nombre M au moins égal à 1, de récepteurs mobiles de positions connues sont situés dans ladite zone géographique donnée et que les puissances Jk,s desdits signaux de brouillage émis par lesdits 20 brouilleurs sont déterminées de façon à ce qu'ils n'affectent pas la réception correcte des services GNSS utilisés par lesdits récepteurs mobiles dudit second groupe, ledit procédé comportant les étapes supplémentaires suivantes : o déterminer ledit seuil de fonctionnement AN;,s pour chacun des 25 récepteurs mobiles dudit second groupe, o déterminer, par une étape de calcul, la puissance optimale Jk,s affectée à chaque signal de brouillage s produit par chacun desdits 10 brouilleurs k et qui respecte en plus les contraintes de non brouillage suivantes: ^ les inégalités suivantes sont respectées pour chaque récepteur i dudit second groupe et chaque service GNSS I utilisé par ledit récepteur i : K S Ek,i .Pk,i .S.SCIk s l "Jk,s / Ak,i ~ A i,l i = N + 1,..., N + M k=1 s=1 3. Procédé selon l'une des revendications précédentes caractérisé en ce que ladite étape de calcul permettant de déterminer lesdites puissances 1 o optimales de brouillage Jk,s est réalisée, si une solution existe, par l'exécution de l'algorithme du Simplexe appliqué aux contraintes suivantes : Max( Ct.x) A.x = b 15 x >=0, où Max(Ct.x) désigne la contrainte selon laquelle le produit scalaire des vecteurs C transposé et x doit être maximum, avec x = E un vecteur de dimension égale à K.S + (N1+ M1+ K) avec J = Jl un vecteur de dimension égale à K.S dont les composantes 20 sont les puissances JI=Jk,s des dits brouilleurs, E= ei un vecteur de dimension N1+M1+K dont les composantes sont des variables spécifiques à l'algorithme du simplexe, N1 est le nombre de contraintes de brouillage à respecter, 0 4. Système de brouillage (100) comportant une pluralité de brouilleurs (109,110,111) adaptés à générer au moins un signal de brouillage permettant de perturber le signal reçu par au moins un récepteur mobile de type GNSS et un organe centralisé (113) arrangés en réseau par des moyens de communication (112), ledit système étant caractérisé en ce que ledit organe centralisé (113) comporte des moyens pour exécuter les étapes du procédé selon l'une des revendications 1 à 3.