FR 2857101 A1 20050107 FR 0307979 A 20030701 PROCEDE DE REJECTION D'INTERFERENCES PERTURBANT LA RECEPTION D'UN SIGNAL DE TRANSMISSION ET DISPOSITIF L'invention concerne un procédé de réjection d'interférences perturbant la réception d'un signal issu d'un système de transmission. Le domaine de l'invention est celui de la réjection d'interférences perturbant la réception de signaux de transmission, notamment des signaux de type GPS (Global Positionning System). Dans le cas d'interférences provenant d'une direction donnée, une solution consiste à disposer d'un réseau de capteurs et à pondérer les différentes voies des capteurs de manière à rejeter partiellement ou en totalité le signal perturbateur provenant de cette direction. Il s'agit alors d'une o pondération spatiale. II existe actuellement différents algorithmes adaptatifs de calcul des poids à attribuer aux différentes voies pour diminuer l'impact des interférences; généralement on effectue une réjection. On connaît les algorithmes itératifs tels que les algorithmes 15 d'Appelbaum, du gradient stochastique, des moindres carrés, ..., basés sur une minimisation de l'erreur quadratique moyenne ou sur la maximisation du rapport signal à bruit. Une autre solution consiste à appliquer l'algorithme d'inversion de puissance aussi appelé algorithme de Capon basé sur la minimisation de puissance en sortie du réseau de capteurs. Cet algorithme peut être aussi bien utilisé en instantané ( snapshot en anglais) que dans une implémentation itérative. On va décrire les principales étapes de l'algorithme de Capon en instantané. On a schématiquement représenté figure 1 un réseau de capteurs 25 Ci recevant un signal perturbateur modélisé par une onde plane incidente de longueur d'onde X. En sortie du réseau de K capteurs le signal Sout est de la forme: Saut (t) = WappiSI + Wapp2S2 + ..., soit K Sou, (t) = E Wappi.Si (t)) i=1 Wappi étant le gain (ou pondération) W; appliqué au capteur Ci et Si étant le signal temporel provenant du capteur Ci. 2857101 2 La puissance de sortie P du type E[ Sout.S*out l], E étant l'espérance intégrée sur un temps long (tendant vers l'infini) s'exprime sous la forme: P = WappH. Rss. Wapp, Avec Wapp=(Wapp1, Wapp2, ... WappK), WappH représentant l'hermitien (c'est-à-dire le transposé conjugué) du vecteur pondération Wapp, et Rss étant la matrice de corrélation entre les signaux Si des différents capteurs. Rss = r11 r21 rKI ru r22 r, K rKK avec r,k (t) = S, (t).Sk (t) Une solution triviale permettant de minimiser la puissance P est: Wapp=O. Pour éviter cette solution triviale, on impose aux coefficients de pondération une contrainte C. Par exemple, une contrainte possible est telle que: Wapp.0 = 1 La solution sous la contrainte C est donnée par: Wcal CHRSS et on a: Wapp=Wcal* Les réjections obtenues par les différents algorithmes sont illustrées figure 2: les algorithmes itératifs ont pour résultat d'abaisser le niveau de l'interférence lnt quasiment au niveau du bruit thermique (cas a de la figure 2) et l'algorithme de Capon abaisse le niveau de l'interférence sous le niveau du bruit thermique à une valeur symétrique de celle de départ (cas b de la figure 2). La limite de sensibilité de ces algorithmes est fixée par rapport au bruit thermique. Ces algorithmes présentent alors des difficultés de réjection pour des interférences dont la puissance est faible, mais déjà suffisante pour 30 dégrader les performances des récepteurs. RSS'C 2857101 3 Un but important de l'invention est donc de proposer un procédé permettant de rejeter les interférences notamment de faible puissance, perturbant la réception d'un signal de transmission. Pour atteindre ce but, l'invention propose un procédé de réjection d'interférences perturbant la réception au moyen d'un système de réception, d'un signal issu d'un système de transmission, le procédé comportant un calcul de pondération des signaux Si(t) i variant de 1 à K, issus du système de réception, basé sur la matrice RSS de corrélation des signaux Si(t) entre eux, la matrice RSS comprenant sur sa diagonale les termes d'autocorrélation des signaux S;(t). Le procédé est principalement caractérisé en ce que le calcul de pondération comporte une étape consistant à modifier les termes d'autocorrélation de la diagonale de la matrice RSS. La puissance du bruit n'intervenant que sur les termes d'autocorrélation (si on intègre sur un temps infini), c'est-à-dire sur la diagonale de la matrice, on peut en modifiant ces termes contrôler la sensibilité des algorithmes de réjection (qui dépend du bruit) et/ou adapter le niveau de réjection des interférences à l'application voulue. Selon une caractéristique de l'invention, la modification des 20 termes de la diagonale comporte une étape pour rendre les termes de la diagonale égaux à une constante, par exemple égale à 1. On obtient ainsi la même contribution de la puissance du bruit sur les différentes voies i. K peut être le nombre de capteurs du système de réception et Si(t) le signal temporel issu du capteur i. Selon une caractéristique de l'invention, Si(t) peut en outre être lui-même déterminé en fonction de M intervalles temporels, chaque capteur étant échantillonné dans le temps. Le système de réception peut ne comporter qu'un capteur échantillonné dans le temps: K est le nombre d'intervalles de temps et Si(t) le signal temporel issu du capteur et correspondant à l'intervalle i. L'invention concerne également un récepteur d'un signal issu d'un système de transmission, susceptible d'être perturbé par des interférences, comportant au moins un capteur et relié à ce (ou chaque) capteur, un dispositif comportant au moins un élément apte à pondérer le signal issu du 2857101 4 capteur correspondant et, relié à l'(aux) élément(s) de pondération, un processeur apte à mettre en oeuvre le procédé décrit. D'autres caractéristiques et avantages de l'invention apparaîtront 5 à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels: la figure 1 déjà décrite représente schématiquement un réseau de capteurs recevant un signal perturbateur modélisé par une onde plane incidente et de signal de sortie Sout, la figure 2 déjà décrite illustre schématiquement les réjections obtenues par les algorithmes itératifs (cas a) et par l'algorithme de Capon (cas b), la figure 3 illustre schématiquement les réjections obtenues par les algorithmes itératifs (cas a) et par l'algorithme de Capon (cas b) à partir du 15 procédé selon l'invention, la figure 4 représente schématiquement un récepteur selon l'invention en cas de réjection spatiale, la figure 5 représente schématiquement un récepteur selon l'invention en cas de réjection temporelle, la figure 6 représente schématiquement un récepteur selon l'invention en cas de réjection spatio-temporelle. On va tout d'abord décrire l'invention sur la base de l'algorithme de Capon instantané et en considérant que le signal de transmission est un 25 signal radio-satellitaire. En négligeant le signal radio-satellitaire très faible et en supposant unitaire le gain des différentes voies I (gain du capteur i inclus), le signal Si(t) s'exprime sous la forme: S,(t) =Aie +b;(t) avec: w la fréquence du signal perturbateur cp sa phase a, sa longueur d'onde A;, l'amplitude du signal perturbateur vu sur le capteur i d; la différence de marche entre le capteur i et le capteur k 35 b;(t) le bruit thermique reçu sur le capteur i. Par conséquent, on peut écrire l'intercorrélation du signal entre le capteur i et le capteur k: rik (0 = Si (t)ÉSik (t) = AiAke,dkd' +bi (t). Ake'rwr+q)+ ) +bk (t) A,e-i(mt+rp+X) +1); (t) bk (t) L'autocorrélation sur le capteur i s'exprime donc sous la forme: rii(t)=Si(t)ÉSH(t)=Ai2+ bi(t).Aie'(wr+v+a)+(t).Aie +11b(t) L'espérance E des rik s'écrit alors: dk -d; dk -d; dk -d; E[rik]=A;Ake' ?t +E[bibk]=A;Ake3 X +E[bi].E[bk]=A;A. ej E[rii]=Ai2 +ELbi21=Ai2 +Pi-bruit car les bruits sont blancs, indépendants et de moyenne nulle. La 15 puissance du bruit n'intervient que sur les termes d'autocorrélation, c'est-à-dire sur la diagonale de la matrice RSSÉ En intégrant sur un temps long qui tend vers l'infini, la matrice RSS s'écrit alors sous la forme: RSS,- _ d k -d; avec Cik = e' 1 Jusqu'à présent on a supposé unitaire le gain sur chaque voie. En 25 réalité chaque voie i présente un certain gain Gi et rii s'écrit: ni = Gi2(Ai2+Pi-bruit La puissance du bruit n'intervient donc pas équitablement sur chaque voie. Pour obtenir la même contribution de la puissance du bruit sur 30 les différentes voies, on peut normaliser la matrice de la façon suivante. On pose r;; = g;2 On peut calculer tous les g;, i=1 à K avant normalisé rik rik = gi.gk La matrice Rss normalisée s'écrit alors: normalisée Rss = normalisé normalisé 1 r21... rK1 normalisé r1, r1 K normalisé 1 On peut multiplier cette matrice par Al2+P1_bruit pour faciliter l'étape suivante de mise en forme de la matrice. On obtient alors: normalisé normalisé 1 r21... r1 1 normalisé r12 1 normalisé rI K 1 Les termes de la diagonale de la matrice sont maintenant normalisés (c'est-à-dire égaux entre eux) à 1 ou à Al2+Pl_bruit: la puissance du bruit a ainsi la même contribution sur les différentes voies. La normalisation est une étape de mise en forme de la matrice qui, appliquée à un algorithme de pondération permet d'améliorer les performances de réjection; cependant, la mise en forme présentée dans le paragraphe suivant peut être effectuée sans normalisation. Pour adapter la réjection à l'application souhaitée, on modifie les termes d'autocorrélation, par exemple par la soustraction d'une matrice diagonale. Dans le cas de la matrice de corrélation non normalisée, on obtient alors: normalisée Rss = (A + Pl bruit) 2857101 7 AI2 +Pl-bruit A1A2C21... AIAKCKI R A2A1C12 A2 2 +P2-bruit AKAICIK A K 2 + PK-bruit 0... 0 PK Dans le cas de la matrice de corrélation normalisée à 1, on obtient Pl 0.. 0 0 P2 alors: normalisé normalisé r21... Tm. 1 0 0 0 1 normalisé r12 1 -P R finale = R ss r1 K normalisé 1 0 0 1 avec P Pour 0 Pour P P peut être déterminé selon plusieurs méthodes, par exemple en fonction de la puissance estimée du bruit. On peut également souhaiter diminuer la réjection dans le cas par exemple où l'on ne dispose pas d'un temps d'intégration assez long: on augmente alors le niveau de la diagonale par exemple par une addition. On peut aussi modifier la réjection en fonction de critères autres que la puissance estimée du bruit ou que le temps d'intégration. On peut par exemple la modifier en fonction de la puissance du signal GPS reçu. On peut également modifier la réjection en fonction de critères déterminés au début du traitement de pondération et restant inchangés pendant ce traitement; selon une variante de l'invention, les critères peuvent aussi être déterminés adaptativement pendant le traitement. 2857101 8 On peut aussi modifier les termes d'autocorrélation de la matrice normalisée ou non par une multiplication de ces termes par un vecteur. Une fois la matrice mise en forme, on applique alors l'algorithme de Capon et on obtient pour le vecteur pondération Wcal: finale 1 R ss ss Les réjections ainsi obtenues sont illustrées figure 3 (cas b). On a montré un exemple d'application à l'algorithme de Capon du 10 procédé selon l'invention. Le procédé selon l'invention s'applique également aux algorithmes itératifs. Ceux-ci comportent une étape de calcul de la corrélation des signaux entre les voies, une étape de calcul de la pondération à appliquer au moyen de l'algorithme choisi et une étape d'application de la pondération aux différentes voies pour la réjection des interférences. Le poids affecté à chaque voie est mis à jour en tenant compte du jeu de poids précédents et à chaque itération, la réjection de l'interférence augmente puis se stabilise. Pour ces algorithmes aussi, le calcul de la pondération est basé sur la corrélation des signaux entre eux, par exemple sous la forme d'une matrice. De même que le procédé selon l'invention s'applique à l'algorithme de Capon par une mise en forme de la matrice de corrélation, de même il s'applique aux algorithmes itératifs par une modification des éléments de corrélation en diminuant (ou augmentant) la valeur des termes d'autocorrélation. Les réjections ainsi obtenues sont illustrées figure 3 (cas a). On a considéré jusqu'à présent une pondération spatiale des signaux issus d'un réseau de capteurs. Le procédé selon l'invention s'applique également à une pondération temporelle des signaux issus d'un capteur échantillonné dans le temps. Les K signaux temporels S,..., Si, ..., SK proviennent du même capteur C et correspondent à K différents intervalles de temps T comme illustré figure 5. Les intervalles de temps sont aussi désignés par décalages temporels. La suite du traitement est alors identique à celui présenté pour le traitement spatial, la pondération étant effectuée sur les différents signaux Wcal c H R finale 1 c ss 2857101 9 décalés dans le temps par le vecteur Wapp déterminé comme précédemment par R finale 1 C _ ss Wcal CHR finale 1C sset Wapp=Wcal* et le signal final se présente sous la forme: K Sou, (t) = E (Wappi.Ji (t)) i=1 Plus généralement, le procédé selon l'invention s'applique à une 10 pondération spatio-temporelle de signaux issus d'un réseau de capteurs C1, Ci, ...CK, chaque capteur étant échantillonné dans le temps comme illustré figure 6. Le signal S1(t) du capteur i s'exprime alors sous la forme: Si(t) = Si1(t) ÉWappii + S12(t).Wappi2 +.. .+ Sij(t)ÉWappij + .. .+ S1M(t).WapPiM, sij(t) étant le signal issu du capteur i et correspondant au décalage temporel j, wij étant la pondération du capteur i et correspondant au décalage temporel j, M étant le nombre d'intervalles temporels du capteur i. Chaque capteur a le même nombre M de décalages temporels T. Le signal Soit s'exprime alors sous la forme: K M sous (t = W aPP ü s (t)) K étant le nombre de capteurs. La matrice Rss de corrélation entre les signaux de chaque capteur et de chaque décalage temporel s'écrit comme suit: i=1 j=1 2857101 10 x11111 x11112 "' rll/IM r11/21 x11/22 "' rll/2M "' rll/K1 rl1/K2 r12/11 r12/12 rIMnM r1M/11 r21/11 x22/11 x21/21 x22/22 R = SS r2M/11 r2M/2M rKl/11 rK2n1 rKM/l l... ... ... rKM/KM _ avec ruh,p (t) = su (t).sRP (t) , i et n correspondant respectivement aux capteurs i ou n et variant de 1 à K, j et p correspondant respectivement aux décalages temporels j ou p et variant de 1 à M. Les termes de corrélation sont formés du produit des signaux de chaque capteur et de chaque décalage temporel, soit pour K capteurs et M 10 décalages temporels, un total de (KM)2 termes. Dans ce cas également les termes d'autocorrélation sont modifiés en fonction de l'application souhaitée. Le signal de transmission peut être un signal radio-satellitaire comme par exemple un signal GPS, un signal GALILEO ou un signal provenant de tout autre système de positionnement par satellites. Le signal de transmission peut être issu de tout système de transmission radio- électrique ou autre, optique par exemple. L'invention concerne également un récepteur R apte à mettre en oeuvre le procédé décrit. On a représenté figure 4 un récepteur R dans le cas d'une réjection spatiale. Il comprend K capteur C1, ... CK et relié à chaque capteur C;, un dispositif comportant au moins un élément Wappi apte à "' rK1/K1 rK2/K2 ri 1/KM 2857101 11 pondérer le signal Si(t) issu du capteur C; correspondant et, relié à aux élément(s) Wappi de pondération, un processeur P apte à mettre en oeuvre le procédé décrit. 2857101 12 REVENDICATIONS 1. Procédé de réjection d'interférences perturbant la réception au moyen d'un système de réception, d'un signal issu d'un système de transmission, le procédé comportant un calcul de pondération des signaux Si(t) i variant de 1 à K, issus du système de réception, basé sur la matrice RSS de corrélation des signaux Si(t) entre eux, la matrice RSS comprenant sur sa diagonale les termes d'autocorrélation des signaux Si(t), caractérisé en ce que le calcul de pondération comporte une étape consistant à modifier les termes d'autocorrélation de la diagonale de la matrice RSS. 2. Procédé selon la revendication précédente, caractérisé en ce que la modification des termes de la diagonale comporte une étape de normalisation de ces termes qui consiste à rendre les termes de la diagonale de RSS égaux à une constante. 3. Procédé selon la revendication précédente, caractérisé en ce que la constante égalel. 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que la modification comporte une étape de soustraction de la matrice RSS par une matrice diagonale prédéterminée ou une étape de multiplication de la matrice RSS par un vecteur prédéterminé. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que le signal étant en outre perturbé par du bruit thermique, les termes d'autocorrélation de la diagonale de la matrice RSS sont modifiés en fonction de la puissance estimée du bruit thermique. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce que la matrice de corrélation étant calculée sur un temps d'intégration déterminé, les termes d'autocorrélation de la diagonale de la matrice RSS sont modifiés en fonction de ce temps d'intégration. 2857101 1 3 7. Procédé selon l'une des revendications précédentes, caractérisé en ce que le signal reçu du système de transmission présentant une puissance déterminée, les termes d'autocorrélation de la diagonale de la matrice Rss sont modifiés en fonction de la puissance de ce signal. 8. Procédé selon l'une des revendications précédentes, caractérisé en ce que la modification est adaptative. 9. Procédé selon l'une des revendications précédentes, caractérisé en ce que le système de réception comprenant un capteur échantillonné dans le temps, K est le nombre d'intervalles de temps et Si(t) le signal temporel issu du capteur et correspondant à l'intervalle de temps i. 10. Procédé selon l'une des revendications précédentes 1 à 8, 15 caractérisé en ce que le système de réception comprenant K capteurs, Si(t) est le signal temporel issu du capteur i. 11. Procédé selon la revendication précédente, caractérisé en ce que chaque capteur i est échantillonné dans le temps, M étant le nombre 20 d'intervalles temporels du capteur i et le signal Si(t) est de la forme: S1(t) = Si1(t).Wappi1 + Si2(t).Wappi2 +...+ Sij(t)ÉWappij + ...+ SIM(t).WappiM, sij(t) étant le signal issu du capteur i et correspondant à un 25 intervalle temporel j, Wappij étant la pondération appliquée au capteur i et correspondant à l'intervalle temporel j. 12. Procédé selon l'une des revendications précédentes, caractérisé en ce que le signal issu du système de transmission est un signal 30 radiosatellitaire ou radio-électrique ou optique. 13. Procédé selon l'une des revendications précédentes, caractérisé en ce que la matrice modifiée est utilisée dans des applications de traitement des interférences basé sur la minimisation de la puissance des interférences. 2857101 14 14. Procédé selon la revendication précédente, caractérisé en ce que le traitement est basé sur un algorithme de type Capon ou de type itératif. 15. Récepteur (R) d'un signal issu d'un système de transmission et susceptible d'être perturbé par des interférences, comportant au moins un capteur (C) et relié à ce (ou chaque) capteur (C), un dispositif comportant au moins un élément apte à pondérer le signal issu du capteur correspondant et, relié à l'(aux) élément(s) de pondération, un processeur (P) apte à mettre en oeuvre le procédé selon l'une quelconque des revendications précédentes. 16. Récepteur (R) selon la revendication précédente, caractérisé en ce qu'il comporte K capteurs (CI, ...C;, ...CK). 17. Récepteur (R) selon l'une des revendications 15 ou 16, caractérisé en ce qu'il comporte pour chaque capteur, M éléments (T) de décalages temporels.