Brevet FR3017716A1 — publié le 20150821

Procédé de traitement antibrouillage d'un signal radioélectrique, module et programme d'ordinateur associés. La présente invention concerne un procédé de traitement antibrouillage d'un signal radioélectrique émis par une source et obtenu par modulation d'une porteuse par un code d'étalement, le code d'étalement résultant de la multiplication d'une forme d'onde par chaque élément d'une séquence binaire propre à la source, ledit procédé de traitement antibrouillage comprenant les étapes suivantes : - réception, en M points d'entrée distincts dudit signal radioélectrique ; - détermination de la corrélation entre eux d'échantillons de signal issus de la réception, à L instants distincts, auxdits M points d'entrée, dudit signal, puis détermination de coefficients de pondération en fonction de ladite corrélation déterminée, M et L étant des entiers positifs, et un au moins de M et L étant strictement supérieur à 1 - obtention d'un signal traité (Sstap) par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée et pondérés par lesdits coefficients de pondération Un tel traitement antibrouillage est par exemple effectué dans des systèmes de positionnement par satellites. Pour rappel, dans de tels systèmes, plusieurs satellites transmettent des signaux radioélectriques. Un récepteur, positionné sur la position à localiser, capte les signaux satellites à l'aide de M capteurs, et les numérise. Puis le récepteur effectue un traitement antibrouillage sur ces signaux captés, par exemple à l'aide de la méthode spatio-temporelle STAP (en anglais « Spatio Temporal Adaptative Processing) selon laquelle le signal traité est une combinaison linéaire des signaux reçus sur les M capteurs et de leurs versions retardées, pondérés par des coefficients de pondérations calculés pour rejeter le brouillage en conservant les signaux satellite utiles. Une fois le traitement antibrouillage effectué, le récepteur effectue une démodulation, dans une boucle PLL (en anglais « Phase Lock Loop »), de la porteuse du signal reçu à l'aide d'une porteuse engendrée localement par un oscillateur piloté, et une démodulation du code d'étalement du signal reçu, dans une boucle DLL (en anglais « Delay Lock Loop »), à l'aide de séquences binaires pseudo-aléatoires engendrées localement par un générateur, afin de dés-étaler les signaux reçus. Puis le récepteur estime les distances qui le séparent des satellites à partir des temps de propagation des signaux satellites captés, et détermine alors la localisation recherchée par triangulation sur la base de ces distances. Le traitement numérique STAP nécessite le calcul préalable d'une matrice d'intercorrélation, nommée R', pour caractériser les brouilleurs à éliminer. Ce calcul est consommateur de ressources matérielles, notamment du fait des nombreux multiplicateurs complexes (les signaux issus de la numérisation sont en effet sous forme complexe). En outre, la méthode d'inversion de puissance, ou méthode de Capon, souvent utilisée pour calculer les coefficients de pondération suppose une puissance de signal utile constante dans la bande passante avant numérisation, ce qui n'est pas le cas avec des signaux de positionnement par satellite à bande étroite et qui conduit à une dégradation des performances. Il est donc nécessaire d'accroître les performances des traitements antibrouillage de l'art antérieur. A cet effet, suivant un premier aspect, l'invention propose un procédé de traitement anti-brouillage du type précité caractérisé en ce que ladite corrélation entre eux desdits échantillons est déterminée après filtrage desdits échantillons par un filtre ayant une réponse impulsionnelle fonction de ladite forme d'onde. Un tel procédé permet de faire prendre en compte, lors du calcul des intercorrélations entre les signaux, les caractéristiques des signaux sources émis. Les performances du traitement antibrouillage se trouvent améliorées et le nombre d'opérations nécessaires au calcul de la matrice d'intercorrélation peut aussi être réduit par sous-échantillonnage. Dans des modes de réalisation, le procédé de traitement antibrouillage suivant l'invention comporte en outre une ou plusieurs des caractéristiques suivantes ledit filtre est associé à une réponse impulsionnelle sensiblement égale à la forme d'onde ; ladite obtention du signal traité est effectuée par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée, filtrés par ledit filtre et pondérés par lesdits coefficients de pondération ; ledit filtre est un filtre à horizon glissant ; ledit filtre est un filtre comprenant un additionneur, un soustracteur, une ligne à retard provoquant un retard égal à la durée de la forme d'onde et un élément à retard provoquant le retard d'une durée séparant deux échantillons du signal sur une boucle de rétroaction de la sortie du soustracteur vers l'entrée de l'additionneur ; ladite corrélation entre eux desdits échantillons est déterminée après un sous-échantillonnage desdits échantillons filtrés ; la détermination de la corrélation entre eux desdits échantillons issus de la réception, à L instants distincts, auxdits M points d'entrée, comprend la détermination d'une matrice d'intercorrélation et la multiplication de l'inverse de ladite matrice d'intercorrélation par une fonction vectorielle de contrainte, ladite fonction vectorielle de contrainte étant fonction de la densité spectrale de puissance de ladite forme d'onde ; ladite fonction vectorielle de contrainte est en outre fonction de la réponse fréquentielle d'un filtre analogique réalisé sur chaque voie entre la réception au point d'accès et la numérisation ; - le signal radioélectrique est de type modulation BOC comportant une première composante centrée sur une fréquence f -carrier fsub-carrier et une deuxième composante centrée sur une fréquence f -carrier + fsub-carrier, et selon lequel ledit procédé de traitement antibrouillage est réalisé de manière indépendante d'une part sur la première composante, et d'autre part sur la deuxième composante ; - le signal radioélectrique source est émis par un satellite d'un système de positionnement par satellite, ledit procédé comportant une étape de détermination de position en fonction d'au moins ledit signal traité obtenu. Suivant un deuxième aspect, la présente invention propose un programme d'ordinateur à installer dans un module de traitement antibrouillage, ledit programme comportant des instructions pour mettre en oeuvre les étapes d'un procédé suivant le premier aspect de l'invention lors d'une exécution du programme par des moyens de traitement dudit module. Suivant un troisième aspect, la présente invention propose un module de traitement antibrouillage d'un signal radioélectrique émis par une source et obtenu par modulation d'une porteuse par un code d'étalement, le code d'étalement résultant de la multiplication d'une forme d'onde par chaque élément d'une séquence binaire propre à la source, ledit module comprenant : - M points d'entrée distincts adaptés pour recevoir ledit signal radioélectrique ; ledit module de traitement antibrouillage étant adapté pour déterminer la corrélation entre eux d'échantillons de signal issus de la réception, à L instants distincts, auxdits M points d'entrée, dudit signal, puis pour déterminer des coefficients de pondération en fonction de ladite corrélation déterminée, M et L étant des entiers positifs, et un au moins de M et L étant strictement supérieur à 1 ; ledit module de traitement antibrouillage étant en outre adapté pour obtenir un signal traité par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée et pondérés par lesdits coefficients de pondération ; ledit module de traitement antibrouillage étant caractérisé en ce qu'il comprend un filtre ayant une réponse impulsionnelle fonction de ladite forme d'onde et en ce qu'il est adapté pour filtrer lesdits échantillons par ledit filtre, ladite corrélation déterminée par le module de traitement antibrouillage étant la corrélation entre eux desdits échantillons filtrés. Ces caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple, et faite en référence aux dessins annexés, sur lesquels : la figure 1 représente une vue d'un récepteur de système de positionnement par satellites implémentant un traitement antibrouillage dans un mode de réalisation de l'invention ; la figure 2 est une vue détaillée de la fonction réalisée par un filtre FIR n°m de la figure 1 ; - la figure 3 est une vue détaillée d'une partie de la vue représentée en figure 1 dans un mode de réalisation, illustrant le calcul de la matrice d'intercorrélation ; - la figure 4 illustre plus en détail le calcul de la matrice d'intercorrélation représenté en figure 3 ; - la figure 5 est une vue d'un filtre Fadaptm adapté à la forme d'onde du signal source dans un mode de réalisation ; la figure 6 est une vue d'un procédé dans un mode de mise en oeuvre de l'invention ; la figure 7 illustre le traitement réalisé par un sous-bloc Cal, du bloc 5 de traitement par canal satellite représenté en figure 1 ; la figure 8 est une vue sur l'axe des temps d'un code d'étalement, d'une porteuse, d'une sous-porteuse et d'un signal émis ; la figure 9 est une vue de la fonction d'autocorrélation en fonction du temps dans le cas de modulation BOC ou BPSK et de la densité spectrale de puissance correspondante. La figure 1 représente une vue d'un récepteur 1 de système GNSS de positionnement par satellites (en anglais « Global Navigation Satellite System »). Le récepteur 1 implémente un traitement antibrouillage dans un mode de réalisation de l'invention. Les s satellites Sat,, i = 1 à s, du système GNSS émettent sur les mêmes fréquences. Chaque satellite Sat, du système génère un signal radiofréquence à l'aide d'une porteuse sinusoïdale (fréquence de 1,2 à 1, 6GHz), éventuellement d'une sous-porteuse dans le cas d'une modulation BOC ou à double décalage de porteuse, modulée par un code d'étalement binaire connu propre à chaque satellite (fréquence de 10 à 1 MHz), puis émet le signal radiofréquence ainsi généré. Les codes d'étalement sont constitués d'un grand nombre de divisions temporelles élémentaires, appelées « chips » de code qui ont une durée égale à Tc. Le code d'étalement est obtenu par multiplication d'une forme d'onde, nommée FO, de durée Tc par successivement chaque élément +1 ou -1 d'une séquence binaire propre au satellite. Dans le cas d'une modulation BPSK (« Binary Phase Shift Keying »), la forme d'onde FO est rectangulaire. La fonction d'autocorrélation du signal émis est donc triangulaire et la densité spectrale de puissance du signal émis est un sinus cardinal au carré. Une modulation BOC (en anglais « Binary Offset Carrier ») utilise une porteuse, une sous-porteuse carrée et un code d'étalement. La forme d'onde FO est en signe de sinus (créneaux). La densité spectrale de puissance du signal émis comporte deux lobes principaux et la fonction d'autocorrélation correspondante comporte des pics multiples. A titre d'illustration, sur la figure 8, les courbes suivantes en fonction du temps t sont représentées : - courbe Cl : code d'étalement de période élémentaire Tc, - courbe C2 : sous-porteuse, de période À -sub-carrier, - courbe C3 : porteuse, de période À -carrier, - courbe C4 : signal émis suite à la modulation BOC selon ces code, porteuse et sous-porteuse. La figure 9 représente, dans sa partie supérieure, la fonction d'autocorrélation y en fonction du temps t dans le cas d'une modulation BOC (traits pleins) ou BPSK (tirets) et, dans sa partie inférieure, la densité spectrale de puissance DSP correspondant à la modulation BOC (traits pleins) ou BPSK (tirets) en fonction de la fréquence f (tarder étant l'inverse de À -carrier et fsub-carrier étant l'inverse de Àsub-camer), pour des code, porteuse et ou sous-porteuse (BOC) correspondant à ceux de la figure 8. Le récepteur 1 comporte des moyens de réception 2 des signaux satellites, par exemple dans le cas considéré une antenne à diagramme de rayonnement contrôlé CRPA (en anglais « Controled Radiation Pattern Array ») comportant M capteurs Am, m = 1 à M, avec M nombre entier supérieur ou égal à 1. Le récepteur comporte en outre un module 3 de prétraitement comportant M voies, un module 4 de traitement antibrouillage et un module 5 de détermination de la localisation du récepteur. Chacun des M. capteurs est adapté pour capter les signaux radiofréquences présents dans la bande de fréquences de fonctionnement des capteurs, puis pour délivrer un signal radiofréquence analogique tel que capté Sm(t), m = 1 à M, au module de prétraitement 3. Le module 3 de prétraitement comprend M voies parallèles. Chaque voie m du module 3 de prétraitement comporte au moins un filtre analogique FRFm, un convertisseur analogique-numérique CANm réalisant un échantillonnage et une numérisation du signal Sm(t), à la fréquence Fe, et un module Fm de passage en bande de base dans lequel le signal reçu numérisé réel est multiplié par un signal numérique complexe exp( 2jn Fe/4 tk ) suivant l'état de l'art, et délivre en sortie de la voie n° m les échantillons Sm(k) d'un signal numérique complexe (Sm(k) = Sm(tk) = Sm(TO+k/Fe)), k entier, m = 1 à M, au module de traitement antibrouillage 4. Le module de traitement antibrouillage 4 comporte un bloc principal 10 de rejet des brouillages et un bloc secondaire 11 de détermination de coefficients de pondération utilisés par le bloc principal 10. Le bloc principal 10 de rejet des brouillages est adapté pour effectuer le traitement STAP. Il comporte M voies parallèles, opérant à la fréquence Le signal Sm(k) y est reçu en entrée de la mème voie, puis est fourni à l'entrée d'un filtre numérique à réponse impulsionnelle finie FIR n° m qui est illustré plus en détail à la figure 2. Il comporte L-1 éléments à retards unitaires. La sortie de chaque élément à retard est multipliée par un coefficient de pondération, de même que le signal Sm(k). Ainsi le FIR rem délivre en sortie un signal S',(k) = E w - b,-'o à L-1 miS' m(k -) . F, Les signaux S'm(k), m = 1 à M, sont ensuite sommés dans le bloc principal 10 : Sstap(k) = m`sem(k) , soit Sstap(k) = Y '...d m.1 à Ael Ei=oà/..114,ffms ,(k - ) = F, b I , Emoà44,1=0 à L-1 Wrn,i'Sm - -Fe) m étant l'indice des capteurs, I étant l'indice des retards. Le signal résultant du traitement antibrouillage Ssta(k) est délivré en sortie du bloc principal 10 et est fourni en entrée du module 5 de détermination de la localisation du récepteur. Le bloc secondaire 11 de détermination de coefficients de pondération est adapté pour déterminer les valeurs actualisées des coefficients de pondérations Win' utilise par le bloc principal 10. Dans le mode de réalisation considéré, il comprend un sous-bloc 12 de filtres adaptés à la modulation comportant M voies de traitement parallèle n° 1 à M, un sous-bloc 13 d'intercorrélation et un sous-bloc 14 de calcul des coefficients de pondération. Dans le mode de réalisation considéré, dans lè sous-bloc 12 de filtres adaptés à la modulation, le signal Sm(k) est délivré, en entrée de la voie n° m, à un fi ltre Fadaptm, dont la réponse impulsionnelle est adaptée, par exemple sensiblement égale, à la forme d'onde FO du signal émis correspondant au signal reçu. Ce filtre Fadaptm permet de tenir compte de la forme d'onde FO dans le calcul des valeurs actualisées des coefficients de pondération Wmt. Dans un mode de réalisation particulier, le filtre adapté à la forme d'onde Fadaptn, est un filtre à horizon glissant, opérant à la fréquence Fe. Dans l'exemple représenté en figure 5, il reçoit en entrée la partie réelle ou imaginaire (I ou Q) du signal complexe Sm(k) et comporte : une ligne à retard comportant un nombre d'éléments unitaires de retard (fonction de transfert égale à z-1 provoquant le retard d'une durée 1/Fe séparant deux échantillons du signal) en série de nature à introduire un retard global égal à la durée Tc d'un chip de code - un soustracteur, un additionneur et un élément à retard unitaire sur une boucle de rétroaction de la sortie du soustracteur vers l'entrée de l'additionneur. Le sous-bloc 13 d'intercorrélation est adapté pour calculer la corrélation entre les signaux issus des M différents capteurs, et/ou leurs versions retardées de 0 à L-1 échantillons, par exemple via le calcul de la matrice d'intercorrélation La matrice R, peut être représentée sous forme d'une matrice à coefficients matriciels [R0] [Ri*T] [Ri+i*T1 [R1÷2*Ti [RL-i*r - - - - [ Rr I [ Ro [Ri*T T R2* [R1,1] [Ri+21 [ [ [ Ro ] R1 * R1+2*1 R 1+1*T [ RI*T avec, [ R2 j [ R1 j [ ROI [ [ R1+2] [R1+1] [R1j [ RO _ L.M x L.M pour I = 0 à L-1, R1= r1,1,1 rmlf rnm,/ rm,A4,1 où = fS'/ t . Sn( t -1/Fd* dt MxM .. rem, (la notation S* indique le conjugué de S ). La durée d'intégration dt (par exemple d'une durée de 100 échantillons) détermine la fréquence de rafraîchissement des valeurs des coefficients de pondération. Rxx a ses coefficients matriciels identiques sur chaque diagonale (matrice de Toeplitz). Rxx est symétrique, au sens hermitien (i.e. RT = R*) , Ro aussi : Ro = Ro*T Le sous-bloc 14 de calcul des coefficients de pondération est adapté pour calculer les coefficients de pondération actualisés Wm,, en fonction de la matrice Rxx actualisée et fournir ces coefficients de pondérations au bloc principal 10. Le module 5 de détermination de la localisation du récepteur comporte un bloc par canal satellite Cali, i = 1 à s, adapté pour démoduler le signal émis par le satellite Set, à l'aide d'une boucle PLL et d'une boucle DLL. Le décalage temporel entre le signal émis par le satellite Sat, à une date connue et le signal reçu par le récepteur à une date à déterminer correspond à la durée recherchée de propagation du signal. Dans le bloc Cal' une réplique du code émis est générée localement. La date de réception du signal est mesurée en mettant en phase le signal reçu et le signal local. Le critère de mise en phase correspond à la maximisation de la fonction de corrélation des deux signaux, des hypothèses de décalage différentes étant considérées pour chaque calcul de corrélation. Ces opérations faisant partie de l'état de l'art sont illustrées succinctement en figure 7. Dans le mode de réalisation particulier illustré en figure 7, les opérations au sein du bloc Cal,m sont réalisées en matériel, celles hors du bloc et au sein du bloc Cali sont réalisées en logiciel. Le module 5 de détermination de la localisation du récepteur est en outre adapté pour déterminer la localisation du récepteur 1 en fonction des durées de propagations des signaux déterminées pour les différents satellites, par triangulation. Les étapes d'un procédé mis en oeuvre selon l'invention sont décrites ci-dessous en référence à la figure 6. Dans une étape 101, les M capteurs Am, m = 1 à M, captent les signaux émis par les satellites sati, i = 1 s et les délivrent au bloc de prétraitement 3. Dans une étape de pré-traitement 102, chacun des signaux analogiques Sm(t), m = 1 à M, transmis par le capteur Am, est filtré à l'aide du filtre FRFm, numérisé à la fréquence d'échantillonnage Fe par le CAN,', puis est transposé en bande de base dans le module Fm. Les échantillons de signal numérique Sm(k), m = 1 à M, résultant de ce pré-traitement sont fournis en parallèle au module de traitement antibrouillage 4. Dans une étape 103 de traitement anti-brouillage, un traitement STAP est mis en oeuvre dans le bloc principal 10 de rejet des brouillages qui effectue, dans une sous-étape 103_1 s'effectuant à la fréquence Fe, une combinaison linéaire des échantillons des signaux Sm(k) et des échantillons précédents, ces échantillons étant pondérés par des coefficients de pondération respectifs. Les valeurs considérées pour ces coefficients de pondérations prennent les dernières valeurs fournies par le bloc secondaire 11 de détermination de coefficients de pondération. Le signal Sstap(k) résultant de cette combinaison linéaire est délivré au module 5 de détermination de la localisation du récepteur. Dans une sous-étape 103_2 réalisée en parallèle à une fréquence F,,= 1/T,,, les coefficients de pondération Wm,, m = 1 à M, I = 0 à L-1, sont actualisés. A cette fin, dans le mode de réalisation considéré, sur chaque voie de traitement parallèle m, les signaux Sm(k) sont filtrés par le filtre Fadaptm du bloc 12, puis, dans le cas considéré en référence à la figure 3, ils sont sous-échantillonnés. La matrice d'intercorrélation R,0, est déterminée par le sous-bloc 13 d'intercorrélation sur de ces signaux filtrés et ici en outre sous-échantillonnés, puis les coefficients de pondération Wm,i, sont déduits par le sous-bloc 14 de calcul des coefficients à partir de la matrice d'intercorrélation Les coefficients de pondération Wm,, ainsi actualisés sont fournis au bloc principal 10 de rejet des brouillages pour qu'ils soient pris en compte dans la suite du traitement STAP à la place des anciennes valeurs précédemment calculées pour ces coefficients. L'étape 103_2 de calcul des coefficients de pondération est réalisée selon une période T, qui dans un mode de réalisation correspond au meilleur compromis entre charge de calcul, réactivité à l'environnement de brouillage et performance de réjection. Dans une étape 104, dans le module 5 de détermination de la localisation du récepteur, un traitement par canal satellite est réalisé, puis la localisation du récepteur 1 est déterminée. La matrice R' définie ci-dessus peut être calculée de différentes manières, en matériel ou en logiciel, par corrélation entre les voies ou par une méthode utilisant la racine carrée triangulaire. Dans un mode de réalisation permettant d'optimiser les performances du traitement antibrouillage, R' est calculée en matériel en temps réel par intercorrélation entre les M voies retardées de 0 à L-1 retards. La figure 3 illustre d'un mode de réalisation particulier dans lequel les M signaux issus des capteurs fournis en entrée du sous-bloc secondaire 11, échantillonnés à la fréquence Fe, sont représentés par le vecteur Z(k), de dimension M, avec Z(k) = [S1(k), S2(k),..., SM(k)]. Le vecteur de signaux Z(k) est fourni au sous-bloc 12 de filtrage adapté à la modulation, puis fourni à l'entrée du sous-bloc 13 d'intercorrélation. Dans le sous-bloc 13 d'intercorrélation, le vecteur de signaux Z(k) est fourni à une ligne à retard comprenant L-1 éléments à retard unitaire 15,, (i = 1 à L-1) disposés en série. Le signal Z(k-I+1) résultant de chaque (1-1)Ième retard unitaire (I = 1 à L) est fourni en entrée d'un sous-échantillonneur, élisant un échantillon sur n échantillons successifs. Le sous-échantillonnage du signal Z(k-I+1) délivre le signal sous-échantillonné nommé Zn(k1+1), I = 1 à L (Ze(k-I+1) = [ Sni(k-1+1), Se2(k-I+1),..., Sem(k-I+1)] où Snm(k) représente l'échantillon k issu du sous-échantillonnage du signal Sem, m = 1 à M) En aval du sous-échantillonnage, la fréquence de travail est donc Fe/n. Chaque vecteur de signaux Ze(k-1+1), I = 1 à L, est fourni en entrée d'un multiplicateur complexe dans lequel chacune de ses coordonnées Snm(k-I+1), m = 1 à M, est multipliée par chaque coordonnée Snm(k) de Ze(k). Le résultat de cette multiplication est une matrice pl, de taille MxM avec pi[u,v]= Pu,v,i = Snu(k-I+1)xSev(k), u= 1 à M etv = 1 à M. La figure 4 illustre plus en détail l'opérateur de corrélation spatiale donnant lieu à l'obtention des termes pu,',, de la matrice pl' u = 1 à M et v = 1 à M. En référence à la figure 3, dans le sous-bloc 13 d'intercorrélation, la matrice pl, ainsi définie est délivrée à un intégrateur c qui est adapté pour réaliser sur chaque terme de la matrice une intégration de durée T. L'intégrateur délivre ainsi le résultat de l'intégration de la matrice pl, égale à la matrice [R1_1], sous-élément de la matrice d'intercorrélation R. Dans un mode de réalisation, la méthode de Capon est utilisée pour déterminer les coefficients de pondération Wei m = 1 à M, I = 0 à L-1, sur la base de l'équation , suivante : W = R 1C , ou CTRxx-1C - W est le vecteur des coefficients de pondération : W = Wo W1 ] et = [ W1,1 W2,I --- ], I = 0 à L-1 ; - C est le vecteur de contrainte : C = Ctemporel Cspatial (® produit de Kronecker), avec ct1 Cté L-1)12 est le vecteur de contrainte temporelle, de dimension Lx1 ; CtL CSi est le vecteur de contrainte spatiale, de dimension Lx1 ; Ctemporel = Cspatial - CSm_ cti.csi => C = Ctemporel ® Cspatial de dimension (L.M)x1 ; cti.csm Les termes du vecteur de contrainte spatiale sont obtenus à partir des déphasages Om des capteurs Am, m = 1 à M, par rapport à la direction satellite choisie : C spatial - Les déphasages em sont calculés à partir d'une direction privilégiée, par exemple la verticale en repère antenne, suivant l'état de l'art. On peut aussi prendre comme vecteur de contrainte spatiale le vecteur [ 1 , 0 , , 0 ]T pour ne privilégier aucune direction. Le vecteur Ctemporel = [ Ctl ct2 ]T est fonction de la densité spectrale de puissance de la forme d'onde FO des signaux satellites. Une telle disposition permet de tenir compte de la densité spectrale des signaux satellites utiles contenus dans les signaux numérisés dans le calcul des pondérations W et ainsi avoir un performance de réjection de brouillage améliorée. Par exemple, dans un mode de réalisation : Ctl = J Hfiltre RF(f) - DSPsignal satellite(f) . exp( 2jit .( I - lo ) . f /Fe ) . df, I = 1 à Loù Hfiltre RF(f) est la réponse fréquentielle des filtres analogiques FRFm, m = 1 à M, supposée identique sur toutes les filtres. DSPsignal satellite(f) est la densité spectrale de puissance des signaux satellites utiles, (par exemple DSPsignal signal satellite) = Tc . sinc( dans le cas d'un signal BPSK) ; lo = (L+1)/2. Il n'est pas nécessaire d'avoir une modélisation fine de la fonction de transfert des filtres analogiques FIRFm,. Un gabarit valant 1 ou 0 selon que l'on est dans la bande passante du filtre ou pas est suffisant : H filtre RF(f) = 1 si f appartient à la bande passante du filtre, H filtre RF(f) = 0 sinon. L'ordre du traitement STAP est défini comme le nombre de retards unitaires (i.e. L- 1). Le nombre L de coefficients de pondération par voie est donc égal à l'ordre augmenté de 1. Un ordre égal à 0 (i.e. L=1) signifie un traitement sans décalage temporel, c'est-à-dire un traitement purement spatial (SAP). Dans un mode de réalisation, on choisit un ordre de valeur paire de manière à avoir un nombre de coefficients impair pour la symétrie du vecteur contrainte temporelle Ctemporel- Cette symétrie évite d'introduire des temps de propagation parasites lors du traitement STAP. On notera que lorsque l'antenne n'a qu'un seul capteur (M = 1), le traitement STAP selon l'invention se réduit à un traitement TAP (Temporal Adaptive Processing) qui permet d'éliminer des brouilleurs à bande étroite (excision de fréquence). La méthode de Capon permet en fait de réaliser une inversion de puissance dans le domaine fréquentiel, plus efficace et robuste qu'une technique de « blanking » fréquentiel : C = Ctemporel Cspatial = Ctemporel, avec Cspatiai = Nul et ct1 Ctemporel = Ct(L-1)/2 CtL Dans le mode de réalisation décrit ci-dessus, le filtrage adapté à la modulation (filtres Fadaptm du sous-bloc 12) était réalisé dans le bloc secondaire 11 de détermination de coefficients de pondération. Dans un autre mode de réalisation, le filtrage adapté est réalisé à la place en amont du module de traitement antibrouillage 4 : le sous-bloc 12 de filtres adaptés à la modulation est alors disposé entre le module 3 de prétraitement et le module de traitement antibrouillage 4. Le traitement STAP est alors réalisé sur les signaux ainsi filtrés. Un sous-échantillonnage est mis en oeuvre dans des modes de réalisation dans le bloc secondaire 11 de détermination de coefficients de pondération comme illustré en figure 3. Dans des modes de réalisation, un sous-échantillonnage synchronisé sur les chips de code local grâce au NCO code est réalisé sur le signal Sstap(k) dans chaque bloc Cal,, permettant de bénéficier du gain d'étalement de la forme d'onde en limitant le nombre de calculs dans ce bloc. Dans un mode de réalisation où la modulation BOC est utilisée, le traitement des signaux reçus pourra être dupliqué : la composante correspondant au lobe gauche, respectivement la composante correspondant au lobe droit, de signal reçu sera fourni en entrée d'un module de traitement antibrouillage 4, après passage en bande de base autour de la porteuse gauche (f carrier fsub-carrier), respectivement la porteuse droite (f carrier fsub-carrier), avec utilisation d'un filtre adapté à la modulation BPSK associée à chaque lobe assimilé à un signal BPSK. Le traitement antibrouillage est donc réalisée de façon indépendante d'une part relativement au lobe droit, d'autre part relativement au lobe gauche. Chacune des composantes résultant du traitement antibrouillage sera fournie en entrée du module 5 de détermination de la localisation du récepteur. Dans le mode de réalisation décrit en référence aux figures, il a été mis en oeuvre un filtrage adapté à la forme d'onde par le bloc 12 de filtres, un sous-échantillonnage, et la prise en compte dans le vecteur de contrainte de la densité spectrale de puissance de la forme d'onde. Dans d'autres modes de réalisation de l'invention, seule l'une ou certaines de ces dispositions sont mises en oeuvre. Certains des traitements détaillés ci-dessus sont mis en oeuvre, dans des modes de réalisation, par le module de traitement anti brouillage 4 suite à l'exécution sur des moyens de calcul du module de traitement antibrouillage 4, d'instructions logicielles stockées dans un mémoire du module 4. La présente invention a été décrite ci-dessus en référence à un récepteur de système de positionnement par satellite. Elle peut bien sûr être mise en oeuvre dans tout autre type de récepteur de signaux de navigation effectuant un traitement antibrouillage, par exemple dans le cas de signaux de navigation émis par des balises au sol. REVENDICATIONS1.- Procédé de traitement antibrouillage d'un signal radioélectrique émis par une source et obtenu par modulation d'une porteuse par un code d'étalement, le code d'étalement résultant de la multiplication d'une forme d'onde par chaque élément d'une séquence binaire propre à la source, ledit procédé de traitement antibrouillage comprenant les étapes suivantes : - réception, en M points d'entrée distincts (Ai, i = 1 à M) dudit signal radioélectrique - détermination de la corrélation entre des échantillons de signal issus de la réception, à L instants distincts, auxdits M points d'entrée, dudit signal, puis détermination de coefficients de pondération en fonction de ladite corrélation déterminée, M et L étant des entiers positifs, et un au moins de M et L étant strictement supérieur à 1 ; - obtention d'un signal traité (Sstap) par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée et pondérés par lesdits coefficients de pondération ; ledit procédé étant caractérisé en ce que ladite corrélation entre des échantillons est déterminée après filtrage desdits échantillons par un filtre (Fadaptm, m = 1 à M) ayant une réponse impulsionnelle fonction de ladite forme d'onde. 2.- Procédé de traitement antibrouillage selon la revendication 1, selon lequel ledit filtre (Fadaptm, m = 1 à M) est associé à une réponse impulsionnelle sensiblement égale à la forme d'onde. 3.- Procédé de traitement antibrouillage selon la revendication 1 ou 2, selon lequel ladite obtention du signal traité (Ssap) est effectuée par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée, filtrés par ledit filtre (Fadaptm, m = 1 à M) et pondérés par lesdits coefficients de pondération. 4.- Procédé de traitement antibrouillage selon l'une des revendications 1 à 3, selon lequel ledit filtre (Fadaptim m = 1 à M) est un filtre à horizon glissant. 5.- Procédé de traitement antibrouillage selon la revendication 4, selon lequel ledit filtre (Fadaptm, m = 1 à M) est un filtre comprenant un additionneur, un soustracteur, une ligne àretard provoquant un retard égal à la durée de la forme d'onde (Tc) et un élément à retard provoquant le retard d'une durée séparant deux échantillons du signal sur une boucle de rétroaction de la sortie du soustracteur vers l'entrée de l'additionneur. 6.- Procédé de traitement antibrouillage selon l'une quelconque des revendications précédentes, selon lequel ladite corrélation entre des échantillons est déterminée après un sous-échantillonnage desdits échantillons filtrés. 7.- Procédé de traitement antibrouillage selon l'une quelconque des revendications 10 précédentes, selon lequel la détermination de ladite corrélation entre des échantillons issus de la réception, à L instants distincts, auxdits M points d'entrée, comprend la détermination d'une matrice d'intercorrélation et la multiplication de l'inverse de ladite matrice d'intercorrélation par une fonction vectorielle de contrainte, ladite fonction vectorielle de contrainte étant fonction de la densité spectrale de puissance de ladite 15 forme d'onde. 8.- Procédé de traitement antibrouillage selon la revendication 7, selon lequel ladite fonction vectorielle de contrainte est en outre fonction de la réponse fréquentielle d'un filtre analogique réalisé sur chaque voie entre la réception au point d'accès et la 20 numérisation. 9.- Procédé de traitement antibrouillage selon l'une quelconque des revendications précédentes, selon lequel le signal radioélectrique est de type modulation BOC comportant une première composante centrée sur une fréquence f -carrier fsub-carrier et une 25 deuxième composante centrée sur une fréquence f -carrier + tub -carrier, et selon lequel ledit procédé de traitement antibrouillage est réalisé de manière indépendante d'une part sur la première composante, et d'autre part sur la deuxième composante. 10.- Procédé de traitement antibrouillage selon l'une quelconque des revendications 30 précédentes, selon lequel le signal radioélectrique source est émis par un satellite d'un système de positionnement par satellite, ledit procédé comportant une étape de détermination de position en fonction d'au moins ledit signal traité obtenu. 11.- Programme d'ordinateur à installer dans un module de traitement antibrouillage, ledit 35 programme comportant des instructions pour mettre en oeuvre les étapes d'un procédésuivant l'une des revendications 1 à 10 lors d'une exécution du programme par des moyens de traitement dudit module. 12.-Module de traitement antibrouillage d'un signal radioélectrique émis par une source et obtenu par modulation d'une porteuse par un code d'étalement, le code d'étalement résultant de la multiplication d'une forme d'onde par chaque élément d'une séquence binaire propre à la source, ledit module comprenant : - M points d'entrée distincts (Ai, i = 1 à M) adaptés pour recevoir ledit signal radioélectrique ; ledit module de traitement antibrouillage étant adapté pour déterminer la corrélation entre des échantillons de signal issus de la réception, à L instants distincts, auxdits M points d'entrée, dudit signal, puis pour déterminer des coefficients de pondération en fonction de ladite corrélation déterminée, M et L étant des entiers positifs, et un au moins de M et L étant strictement supérieur à 1 ; ledit module de traitement antibrouillage étant en outre adapté pour obtenir un signal traité (Sstap) par combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée et pondérés par lesdits coefficients de pondération ; ledit module de traitement antibrouillage étant caractérisé en ce qu'il comprend un filtre (Fadaptm, m = 1 à M) ayant une réponse impulsionnelle fonction de ladite forme d'onde et en ce qu'il est adapté pour filtrer lesdits échantillons par ledit filtre, ladite corrélation déterminée par le module de traitement antibrouillage étant la corrélation entre lesdits échantillons filtrés. 13.- Module de traitement antibrouillage selon la revendication 12, dans lequel ledit filtre (Fadaptm, m = 1 à M) est associé à une réponse impulsionnelle sensiblement égale à la forme d'onde. 14.- Module de traitement antibrouillage selon la revendication 12 ou 13, adapté pour obtenir ledit signal traité (Sstap) en effectuant une combinaison linéaire desdits échantillons de signal issus de la réception à L instants distincts, auxdits M points d'entrée, filtrés par ledit filtre (Fadaptm, m = 1 à M) et pondérés par lesdits coefficients de pondération. 15.- Module de traitement antibrouillage selon l'une des revendications 12 à 14, dans lequel ledit filtre (Fadaptim m = 1 à M) est un filtre à horizon glissant. 16.- Module de traitement antibrouillage selon la revendication 15, dans lequel ledit filtre (Fadaptm, m= 1 à M) comprend un additionneur, un soustracteur, une ligne à retard provoquant un retard égal à la durée de la forme d'onde (Te) et un élément à retard provoquant le retard d'une durée séparant deux échantillons du signal sur une boucle de rétroaction de la sortie du soustracteur vers l'entrée de l'additionneur. 17.- Module de traitement antibrouillage selon l'une quelconque des revendicationsl2 à 16, comprenant un sous-échantillonneur adapté pour sous-échantillonner lesdits échantillons filtrés, ladite corrélation déterminée par le module de traitement antibrouillage étant la corrélation entre lesdits échantillons filtrés et sous-échantillonnés. 18.- Module de traitement antibrouillage selon l'une quelconque des revendications12 à 17, adapté pour, lors de la détermination de la corrélation entre lesdits échantillons issus de la réception, à L instants distincts, auxdits M points d'entrée, déterminer une matrice d'intercorrélation et pour multiplier l'inverse de ladite matrice d'intercorrélation par une fonction vectorielle de contrainte, ladite fonction vectorielle de contrainte étant fonction de la densité spectrale de puissance de ladite forme d'onde. 19.- Module de traitement antibrouillage selon la revendication 18, dans lequel ladite fonction vectorielle de contrainte est en outre fonction de la réponse fréquentielle d'un filtre analogique réalisé sur chaque voie entre la réception au point d'accès et la numérisation. 20.- Module de traitement antibrouillage selon l'une quelconque des revendicationsl2 à 19, dans lequel le signal radioélectrique est de type modulation BOC comportant une première composante centrée sur une fréquence f -carrier fsub-carrier et une deuxième composante centrée sur une fréquence f -carrier + tub-carrier, ledit module de traitement antibrouillage étant adapté pour déterminer des premiers coefficients de pondération en fonction de ladite corrélation des échantillons déterminée pour la première composante exclusivement et pour déterminer des seconds coefficients de pondération en fonction de ladite corrélation des échantillons déterminée pour la deuxième composante exclusivement, ledit module de traitement antibrouillage étant adapté pour obtenir une première composante de signal traité (Sstap) par combinaison linéaire desdits échantillons de la première composante exclusivement pondérés par lesdits premiers coefficients depondération et pour obtenir une seconde composante de signal traité (Sstap) par combinaison linéaire desdits échantillons de la seconde composante exclusivement pondérés par lesdits seconds coefficients de pondération. 21.- Module de traitement antibrouillage selon l'une quelconque des revendicationsl2 à 20, pour un système de positionnement par satellite, dans lequel le signal radioélectrique source est émis par un satellite du système de positionnement par satellite, ledit module de traitement étant adapté en outre pour déterminer une position en fonction d'au moins ledit signal traité obtenu.10