Brevet FR3008800A1 — publié le 20150123

DISPOSITIF DE DETECTION DE SIGNAUX ELECTROMAGNETIQUES La présente invention concerne un dispositif de détection de signaux électromagnétiques. Le domaine de l'invention est celui de la séparation et de la localisation de sources électromagnétiques à partir d'antennes de type réseau. Le domaine d'exploitation couvre potentiellement l'ensemble du 5 domaine des récepteurs de signaux électromagnétiques. Plus particulièrement, l'invention concerne le domaine des intercepteurs électromagnétiques de type actif ou passif, dans lesquels on cherche à la fois une large couverture angulaire instantanée et une capacité de séparation angulaire importante, en particulier pour discriminer les signaux reçus de 10 faible puissance et à impulsions de longue durée en environnement dense. Dans ce contexte, un problème technique à résoudre est la minimisation du nombre de voies de réception dans un dispositif d'antenne de type réseau, ayant pour objectif de réduire à la fois la complexité, le volume de calcul à 15 effectuer et finalement le coût. Le problème est habituellement résolu par des techniques de formation de faisceaux par le calcul ou par des techniques d'interférométrie. Un pouvoir séparateur angulaire important nécessite un faisceau d'antenne étroit et donc une grande antenne. Ainsi, plus le pouvoir séparateur angulaire 20 exigé par une application est important, plus la dimension L de l'antenne doit être importante. En effet l'ouverture de l'antenne, notée AO, une fois focalisée est, pour une longueur d'onde À, typiquement de la forme AO = À/L en radians, donc inversement proportionnelle à la dimension L. La grande dimension de l'antenne conduit alors utiliser un grand nombre 25 d'éléments rayonnant constituant cette antenne. Ces éléments rayonnants doivent être espacés entre eux d'une fraction de longueur d'onde, typiquement À/2 pour une couverture angulaire totale de 90°. Cette condition est en effet nécessaire pour éviter les ambiguïtés de direction angulaire. Ainsi pour une antenne d'ouverture 1/N, le nombre d'éléments rayonnant 30 sera en général typiquement de 2N pour une antenne linéaire et de 4N2 pour une antenne surfacique. A chaque élément rayonnant correspond un récepteur, ce qui conduit à une grande complexité et à un coût élevé. Il existe une solution permettant de diminuer le nombre de récepteurs en faisant appel à des antennes dites « lacunaires » ou encore à des réseaux d'interférométrie ambigus. En ce qui concerne les antennes lacunaires, les dégradations du diagramme 5 de rayonnement augmentent en fonction du taux de lacunarité, ce qui en pratique ne permet pas de diminuer le nombre de récepteurs de façon significative. En ce qui concerne les antennes interférométriques, elles permettent d'obtenir de bonnes précisions de localisation avec un faible nombre 10 d'éléments rayonnants, la valeur de la précision étant inversement proportionnelle à la dimension de l'antenne. Cependant elles nécessitent de lever les ambiguïtés de mesure de direction angulaire à l'aide d'au moins deux bases d'interférométrie de dimensions différentes. Cette méthode de levée d'ambiguïté peut être mise en échec lorsque la densité de signaux à 15 traiter est importante ou lorsque le rapport signal sur bruit est faible. En présence de signaux de longue durée, par exemple de type LPI (« Low Probability of lntercept »), le risque de superposition temporelle de signaux augmente, ce qui complique encore l'analyse des signaux reçus. Il peut alors devenir impossible de séparer les signaux par leur direction ou 20 d'estimer cette direction d'arrivée. Un but de l'invention est de pallier les inconvénients précités, à cet effet l'invention a pour objet un dispositif de détection de signaux électromagnétiques comportant une antenne de réception réseau à N 25 éléments rayonnants et comportant M voies de réception en aval de l'antenne de réception, M étant inférieur à N, ledit dispositif comportant par ailleurs au moins : - des moyens de commutation des M voies de réception sur les éléments rayonnants selon des cycles de séquences successives, M 30 éléments rayonnants étant connectés auxdites voies de réceptions à chaque séquence, un même élément rayonnant dit de référence étant connecté auxdites voies de réception pour toutes les séquences, un cycle étant achevé lorsque tous les éléments rayonnants ont été connectés au moins une fois à une desdites voies de réception ; des moyens de traitement effectuant pour chaque séquence l'estimation des corrélations spatiale deux à deux du signal reçu sur la voie de référence et des signaux reçus sur les autres M-1 voies de réception puis estimant la densité spectrale de puissance spatiale dans N directions d'arrivée à partir d'une somme cohérente des N termes de corrélations ainsi obtenus. Le réseau d'éléments rayonnant étant linéique, l'élément de référence est par exemple le premier éléments, les M premiers éléments étant connectés aux voies de réception dans la première séquence, puis les M-1 éléments rayonnant étant connectés aux voies de réceptions dans la deuxième séquence, et ainsi de suite. Dans un mode de réalisation possible, les moyens de traitement estiment les corrélations spatiales après avoir effectué l'élimination des signaux de courte durée stationnaire de durée inférieure à une durée t donnée, l'élimination des signaux de courte durée stationnaire étant effectuée en multipliant le signal reçu à un t, x(t), par un signal de la forme x(t)x(t-T)/x(t)*2, où t est un retard choisi pour supprimer les signaux de durée inférieure à la durée T. Dans chaque séquence, ladite corrélation spatiale est par exemple estimée en multipliant directement le signal reçu sur ledit élément rayonnant de 20 référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne. Dans un autre mode de réalisation possible, dans chaque séquence, ladite corrélation spatiale est estimée en multipliant le signal reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément 25 rayonnant de l'antenne, après séparation fréquentielle desdits signaux reçus. Dans chaque séquence, ladite corrélation spatiale est par exemple estimée en multipliant le signal reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne, l'estimation de corrélation étant effectuée en sortie de bancs de filtres de 30 différentes largeurs et de différentes fréquences centrales. Les filtres sont par exemple réalisés par transformées de Fourrier numérique d'ordres différents ou par filtres polyphasés d'ordres différents. Le dispositif est par exemple un radar. Dans ce cas, dans chaque séquence, ladite corrélation spatiale est par exemple estimée en multipliant le signal 35 reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne, après séparation temporelle desdits signaux reçus. La séparation temporelle est par exemple obtenue par filtrage adapté et échantillonnage des impulsions reçues par ledit radar. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite en regard de dessins annexés qui représentent : - la figure 1, par un schéma synoptique, un réseau d'antenne linéique comportant N éléments rayonnant ; - la figure 2, un diagramme d'antenne relatif à un dispositif du type de celui de lafigure 1 ; - la figure 3, un exemple de dispositif d'interférométrie ; les figures 4a et 4b, un diagramme d'antenne obtenu respectivement par un interféromètre ambigu et par un interféromètre non ambigu ; la figure 5, une illustration du principe de l'invention ; - la figure 6, par un synoptique, un premier exemple de réalisation d'un dispositif selon l'invention ; - la figure 7, un exemple de réalisation des moyens de commutations utilisés dans un dispositif selon l'invention ; la figure 8, un autre exemple de réalisation d'un dispositif selon l'invention. La figure 1 représente, par schéma synoptique, un réseau d'antenne linéique comportant N élément rayonnants. L'invention s'applique à une antenne surfacique, de dimension 2, ou à une antenne linéique, de dimension 1, comme celle illustrée par la figure 1. A titre d'exemple, l'invention sera décrite par la suite pour une antenne linéique. Le réseau 10 comporte N éléments rayonnants 1 espacés entre eux d'une demi-longueur d'onde. Dans un récepteur à formation de faisceaux par le calcul, chaque élément rayonnant 1 est connecté à un récepteur effectuant l'amplification, le filtrage et la transposition en fréquence du signal reçu avant codage. Ainsi en réception pour effectuer la transposition de fréquence, chaque élément rayonnant 1 est connecté à une entrée d'un mélangeur hyperfréquence 2, l'autre entrée du mélangeur recevant une fréquence intermédiaire. La sortie d'un mélangeur 2 est connectée à l'entrée d'un amplificateur faible bruit 3. Les amplificateurs 3 sont connectés en sortie à l'entrée d'un codeur 4 à N voies. Ce dernier assure la conversion analogique-numérique des signaux issus des amplificateurs. Les signaux reçus et codés sont traités simultanément par une unité de traitement 5 pour former les faisceaux. L'opération de formation de faisceaux par le calcul consiste à sommer de façon cohérente l'ensemble des signaux x reçus pour chaque direction de pointage du faisceau. Elle peut s'écrire selon la relation (1) suivante, S représentant l'amplitude du signal, fonction de la direction de visée et de la 10 période d'échantillonnage : N-1 sin(ek S(0k,nT,.)=Exi(nTr)exp i=o où : - ek est l'angle de visée dans la direction k8 (8k=k0) ; 15 - Tr est la période d'échantillonnage des signaux reçus, à la réception au niveau des codeurs 4 ; - n est l'indice de l'instant d'échantillonnage auquel est effectuée la formation de faisceau ; i est l'indice de l'élément rayonnant ; 20 - d la distance entre deux éléments rayonnants successifs ; la longueur d'onde du signal ; D'autre part, les directions de pointage sont régulièrement espacées et leur nombre est égal au nombre N d'éléments rayonnants du réseau. Dans ces 25 conditions, la formation de faisceau correspond simplement à une transformée de Fourier et peut s'écrire selon la relation (2) ci-dessous : N-1 (2) S(0 k,n7;) = xi (nTr)exp N i=0 30 où la fréquence fk et le temps t s'expriment de la façon suivante : sin ek J k = t = i.d (1) L'indice k correspond à l'indice de la direction d'arrivée ek dans laquelle on souhaite former un faisceau. Quand sine = AkiNd, on obtient un maximum dans la direction ek. Le dispositif illustré par la figure 1 permet ainsi de former N faisceaux directifs, d'ouverture typique égale à 2/N radians, sans ambiguïté et avec des niveaux de lobes secondaires ajustables par pondération numérique. Ce dispositif permet donc de séparer angulairement les sources d'émission et de les localiser avec une grande efficacité. La figure 2 présente par une courbe 21 un exemple de diagramme d'antenne obtenu sans pondération pour un dispositif du type de la figure 1, pour un pointage particulier, les ordonnées représentant le gain en dB et les abscisses la position angulaire. Malheureusement un tel dispositif peut s'avérer très complexe, en particulier 15 lorsque le nombre de voies est grand du fait d'une fréquence d'échantillonnage élevée et/ou d'une ou plusieurs des exigences suivantes : - résolution angulaire fine ; - couverture angulaire importante ; 20 - fréquence RF élevée. Par exemple, dans un dispositif passif d'interception des signaux radar, ces conditions se trouvent réunies, ce qui rend très compliqué la mise en oeuvre d'une telle solution. C'est notamment une raison pour laquelle on utilise 25 généralement des bases interférométriques, conformément à la figure 3, afin de localiser angulairement les différentes sources d'émission dans un vaste domaine angulaire. La figure 3 présente donc un dispositif formé de bases interférométriques, 30 utilisé pour localiser angulairement les différentes sources d'émission dans un domaine angulaire important. Dans un tel dispositif, au moins une mesure très précise mais ambigüe est réalisée par une base de grande dimension constituées par les éléments d'antenne les plus éloignés, et au moins une mesure non ambigüe mais peu 35 précise est réalisée par une base constituée par les deux éléments d'antennes les plus proches. Le dispositif de la figure 3 comporte quatre éléments rayonnant 31, 32, 33, 34, sachant que trois éléments au moins pourraient suffire pour effectuer une mesure de direction d'arrivée de signaux par interférométrie. En considérant la distance élémentaire d, formant un pas, et en considérant par ailleurs la position du premier élément 31 comme position d'origine, le deuxième élément est alors à la position d, le troisième éléments à une position kd et le quatrième élément placé à une position Nd. Les éléments rayonnants sont connectés chacun à un mélangeur 2, recevant sur une autre entrée une fréquence FI issue d'un oscillateur local, suivi d'un amplificateur faible bruit 3 de façon analogue au dispositif de la figure 1. Les signaux issus des amplificateurs 3 sont ensuite convertis numériquement par un convertisseur analogique-numérique 4 à plusieurs voies, de même type par exemple que celui de la figure 1. Les signaux ainsi numérisés sont pris en compte par des moyens de traitement 35 pour la mesure de direction d'arrivée par interférométrie. La combinaison des mesures de phases obtenues sur les différentes bases, constituées des éléments rayonnant 31, 32, 33, 34, permet de lever l'ambiguïté de la mesure. Typiquement, en considérant le schéma synoptique de la figure 3, il est 20 possible pour une longueur d'onde donnée À, d'estimer les différences de phase suivantes, respectivement entre le premier élément et les éléments aux positions d, kd et Nd : (1) (D2 = À, 25 -k = 27rdk sin(0) 274c1N sin(0) ci) cDN = 11.- puis d'en déduire après levée d'ambiguïté conforme aux méthodes connues de l'homme de l'art, d'en déduire la direction d'arrivée du signal 8. Le dispositif de la figure 3 qui ne comprend qu'un faible nombre de voies de 30 réception est beaucoup plus simple que celui de la figure 1. Cependant, s'il permet de localiser les différentes directions d'arrivée de signaux issus de sources d'émission avec une bonne précision, il ne permet généralement pas 27-cc/ sin(0) à lui seul de séparer ces différentes sources d'émission, et peut être mis en échec si le rapport signal à bruit n'est pas suffisant. Les figures 4a et 4b représentent respectivement un diagramme d'antenne 5 obtenu respectivement par un interféromètre ambigu et par un interféromètre non ambigu. Plus particulièrement la courbe 41 de la figure 4a représente le diagramme d'antenne pour un interféromètre constitué des deux bases 31, 34 les plus éloignées, étant distancées de Nd, dans un système d'axes où les ordonnées représentent le gain d'antenne et les abscisses l'angle O. La 10 courbe 42 de la figure 4b représente le diagramme d'antenne pour un interféromètre constitué des deux bases 31, 32 les plus rapprochées, distancées de d, dans le même système d'axes. Le principe de l'invention est illustré par la figure 5. Selon l'invention, à partir 15 d'un dispositif du type de celui de la figure 1, à N éléments rayonnants, on réduit le nombre de récepteur de N à M, M étant inférieur à N, en commutant les M récepteurs d'une certaines façon sur les différents éléments rayonnants 1, au cours du temps. Dans la suite de la description, on prendra à titre d'exemple M = 4 et N = 16. Les M récepteurs forment ce que l'on 20 appellera par la suite le dispositif de réception. A condition de connecter chaque élément rayonnant au dispositif de réception au moins une fois pendant le temps d'observation, il est alors possible de séparer les signaux d'arrivée stationnaire pendant le temps d'observation selon N faisceaux directifs conformément au résultat obtenu 25 par une formation de faisceau par le calcul conventionnelle. La figure 5 illustre un exemple de séquences de commutation 51, 52, 53, 54, 55 incluses dans le temps d'observation permettant connecter au moins une fois chaque élément rayonnant au dispositif de réception. Les différentes séquences illustrées représentent la position des N récepteurs, 16 dans 30 l'exemple de la figure 5, les récepteurs repérés par un point noir étant ceux qui sont connectés au dispositif de réception. Les éléments sont par exemples espacés entre eux d'une distance A/2. Pour toutes les séquences, le premier élément rayonnant 50 est en connecté au dispositif de réception. Ce premier élément 50 constitue ainsi l'élément de 35 référence. Il n'est pas nécessaire de prendre le premier élément 50 comme élément de référence, d'autres éléments du réseau pourraient constituer un élément de référence, pourvu que cet élément soit connecté au dispositif de connexion à chaque séquence. Dans la première séquence 51, les M premiers éléments sont connectés au 5 dispositif de commutation. Dans la deuxième séquence 52, le premier élément 50 est connecté et les M-1 éléments suivant 521 sont connectés. Dans la troisième séquence 53, le premier élément 50 est connecté et les M1 éléments suivant 531 sont connectés. Dans la quatrième séquence 54, le premier élément 50 est connecté et les M-1 éléments suivant 541 sont 10 connectés. En fin, dans la dernière séquence, le premier élément 50 est connecté et les M-1 éléments suivant 551 sont connectés. Ces derniers M-1 éléments 551 sont les M-1 éléments du réseau à N éléments, N étant égal à 16 et M étant égal à 4. Le mode de connexion des éléments au dispositif de réception illustré par la 15 figure 5 a l'avantage d'être particulièrement simple à mettre en oeuvre. D'autres modes de connections sont bien sûr possibles pourvu qu'ils permettent que chaque élément au moins soit connecté au cours de l'ensemble des séquences 51,52, 53, 54, 55. 20 A chaque séquence on calcule l'espérance mathématique du produit du signal reçu sur l'élément d'antenne de référence par le complexe conjugué de chacun des signaux reçus des M-1 autres éléments d'antennes connectés au dispositif de réception, soit pour i compris entre 1 et M : 25 Rx(i,nT,.)= Etx1(nTr ) +1(nTr )] (3) Où: x1(nTr) est le signal reçu sur la voie 1 à la date nTr, Tr étant la période d'échantillonnage précédemment définie ; xi,;(nTr) est le signal reçu sur la voie i+1 à la date nTr, x*i.i(nTr) étant son 30 conjugué ; E[x] représente l'espérance mathématique sur nTr échantillons temporels à partir d'une position pTr : E[x(p.nT,.)]=1Ëx(pT,.+1.7;) (4) n 1=0 Rx(i,nTr) représente la corrélation entre le signal reçu de la voie 1 et le complexe conjugué du signal reçu de la voie i+1 à la date nTr. En supposant que le signal est stationnaire, à bande étroite vis-à-vis de la porteuse, et que l'amplitude du signal reçu et sa direction d'arrivée sont constantes pendant le temps d'observation, on peut considérer que la corrélation entre les signaux reçus de deux éléments rayonnants quelconques est invariante au cours du temps, au bruit de mesure près. Ainsi, quels que soient les instants temporels nTr et pTr : Rx(i,nT,.)= Rx(i,p.nT,.)= Rx(i) (5) Et pour un signal reçu d'amplitude A, de longueur d'onde À et de direction d'arrivée 8, la corrélation Rx(i) s'exprime comme suit : R x(i) = A2 exp- (6) On peut alors estimer la densité spectrale S(8k) de puissance spatiale dans un filtre spatial de direction 8k selon la relation suivante : N-1 2.j.z.i.k s(ek)=ERx(i)exp N (7) i=0 k ayant été défini relativement à la relation (2). Il est à remarquer que, même si l'amplitude du signal A n'est pas rigoureusement constante au cours du temps, la formation de faisceau reste possible, avec cependant une dégradation de la résolution spatiale. La figure 6 présente par un synoptique un exemple de réalisation d'un dispositif selon l'invention dans le cas où le dispositif de réception comporte quatre récepteurs, soit M = 4, dispositif auquel peut s'appliquer les séquences de la figure 5. Le dispositif comporte un réseau 10 de N éléments rayonnant 1, N étant égal à 16. Les éléments rayonnant sont reliés au dispositif de réception via une matrice de commutation 61 permettant de réaliser les différentes séquences de commutation, par exemples les séquences 51, 52, 53, 54, 55 décrites précédemment. La sortie de la matrice de commutation est reliée aux quatre récepteurs, chaque récepteur comportant au moins un mélangeur 2, transposant le signal reçu à une 5 fréquence intermédiaire FI, et un amplificateur 3. Les signaux en sortie d'un récepteur sont numérisés par un convertisseur analogique-numérique 62 à M voies. Les signaux de réception numérisés sont traités par des moyens de traitement 63 effectuant les calculs définis aux relations (3) à (7) définies précédemment. Les moyens de traitement peuvent être réalisés par un 10 circuit de type FPGA ou encore par un processeur de traitement du signal. En cas d'application des séquences se la figure 5, les moyens de traitement calculent les corrélations R'(i) comme indiqué ci-dessous pour chaque séquence, où ils effectuent une corrélation des voies de réception deux à deux avec la voie de référence. 15 1ère séquence 51 : Réception des éléments rayonnants des positions 1, 2, 3 et 4. Calcul de Rx(i) sur l'horizon temporel [0, nTr] : Rx(0) = E(xi(nTr)x* 1 (nTr )) 20 R,(1)= E(x1(nTr)x*2(nT')) R, (2) = E(xl(nT')x*3(nTr)) R, (3) = E(xi(nT,)x* 4 (nTr )) 2ème séquence 52 : Réception des éléments rayonnants des positions 1, 5, 6 et 7. 25 Calcul de Rx(i) sur l'horizon temporel [nTr, 2nTr] : R' (4) = E(x1(2nTr)x5* (2nT,.)) R,(5) = E(x1(2nTr)x; (2nT,.)) R, (6) = E(xi (2nTr)x; (2nT,)) 3ème séquence 53 : 30 Réception des éléments rayonnants des positions 1, 8, 9 et 10. Calcul de Rx(i) sur l'horizon temporel [3nTr, 4nTr] : R'(7)= E(xi(3nT,.)x; (3nTr )) Rx (8) = E(xi (3nT,.)x9* (3nTr )) Rx (9) = E(xl (3nT,, )xio* (3n Tr)) 4è" séquence 54 : Réception des éléments rayonnants des positions 1, 11, 12 et 13. Calcul de Rx(i) sur l'horizon temporel [4nTr, 5nTr] : Rx(10) = E(xi(4nTr)xii* (4nTr)) Rx(11)= E(x1(4n7)x12* (4nTr)) Rx (12) = E(xi (4nT,.)x13* (4nTr )) sème séquence 55 : Réception des éléments rayonnants des positions 1, 14, 15 et 16. Calcul de Rx(i) sur l'horizon temporel [5nTr, 6nTr] : Rx (13) = E(x1(5nT,.)x14* (5nTr )) Rx (14) = E(x1(5nT,.)x15* (5nTr )) Rx(15)= E(x1(5nT,.)x16* (5nTr )) Finalement, l'estimation de la densité spectrale de puissance spatiale dans 20 différentes directions ek permet de déterminer la direction d'arrivée de l'onde reçue, comme l'exprime la relation suivante : 15 s(ek) = ex(t) exp 16 (8) r=o Le calcul peut également être effectué séquentiellement, en estimant à 25 chaque séquence les sommes partielles et en les cumulant de séquence à séquence avant d'aboutir à la somme totale. La figure 7 illustre un exemple de réalisation de la matrice de commutation 61. La matrice de commutation comporte N entrées et M sorties, plus particulièrement 16 entrées et 4 sorties dans l'exemple de la figure 7. Elle est réalisée à l'aide de blocs commutateurs 701, 702, 703, 704 à diodes PIN, du type SPXT. Chaque bloc comporte trois commutateurs élémentaires. La première entrée 71 est reliée par une ligne directe 72 à une première sortie 73, permettant à l'élément rayonnant de référence d'être en permanence connecté au dispositif de réception. Les trois entrées suivantes 711 sont reliées chacune à un commutateur d'un premier bloc 701. Les trois entrées suivantes 712 sont reliées à ces mêmes commutateurs si bien que la 2ème entrée est commutée avec la Sème entrée et ainsi de suite. Les trois entrées suivantes 713 sont reliées chacune à un commutateur d'un deuxième bloc 702. Les trois entrées suivantes 714 sont reliées à ces mêmes commutateurs si bien que la Sème entrée est commutée avec la 12ème entrée et ainsi de suite. Les sorties du premier bloc 701 sont reliées chacune à un commutateur d'un troisième bloc 703. Les sortie du deuxième bloc 702 sont reliées à ces mêmes commutateurs, si bien que la première sortie du premier bloc 701 est commutée avec la première sortie du deuxième bloc 702 et ainsi de suite. Les sorties du troisième bloc 703 sont reliées chacune à un commutateur d'un quatrième bloc 704. Les trois dernières entrées 715 de la matrice sont connectées à ces mêmes commutateurs si bien que la 14è" entrée est commutée avec la première sortie du troisième bloc 703 et ainsi de suite. Les sorties du 4ème bloc 704 sont connectées chacune à un récepteur 2, 3 du dispositif de réception. L'architecture de la matrice de commutation de la figure 7 permet notamment de mettre en oeuvre les séquences de commutation illustrées par la figure 5. La première séquence est ainsi réalisée en commutant le premier bloc 701 sur les 2è", 3ème et 4ème entrées, puis en commutant le troisième bloc 703 sur les sortie du premier bloc 701 et enfin en commutant le quatrième bloc 704 sur les sortie du troisième bloc 703. L'invention s'applique notamment à la détection et à la localisation de 35 signaux radar. Dans ce cadre, elle peut s'appliquer à un intercepteur actif de type radar ou à un intercepteur passif, par exemple de type RESM (Radar Electronic Support Measure), adapté à la détection de signaux à faible rapport signal à bruit et à impulsions de longue durée, typiquement supérieure à 100 ps. Dans le cas d'un intercepteur actif de type radar, les caractéristiques temporelles et fréquentielles de signaux reçus sont connues. C'est le cas classiquement, en radar, où les signaux reçus font l'objet d'un filtrage adapté à la forme d'onde émise, ce filtrage permettant à la fois d'optimiser le rapport signal à bruit et de séparer ces signaux dans différentes cases distance et/ou Doppler avant détection. Dans ce cas, la corrélation spatiale est estimée après séparation temporelle des signaux reçus. A l'issue de chaque période de traitement Tr, on dispose donc d'une matrice distance-Doppler dans laquelle se répartissent l'ensemble des signaux reçus. Si on considère, comme précédemment une antenne à N éléments rayonnant, à chaque élément rayonnant connecté au récepteur correspond à chaque période de traitement une matrice distance-Doppler. Le traitement de formation de faisceau décrit précédemment s'applique en parallèle sur chacune des sorties des cases distance-Doppler. Dans le cas d'un intercepteur passif de type RESM, les caractéristiques temporelles et fréquentielles des signaux reçus sont inconnues. Un premier tri grossier des signaux est par exemple effectué par filtrage fréquentiel analogique, les récepteurs étant à bande limitée. Après numérisation du signal, un second tri peut être effectué pour limiter à l'entrée du traitement la présence de signaux de courte durée, correspondant à des signaux non stationnaires. Ceci peut être réalisé en multipliant le signal reçu x(t) à l'instant t par un signal de pondération Pond, défini par la relation suivante : Pond=x(t)x(t-r)/x(t)*2 (9) OÙ t est un retard choisi pour supprimer les signaux de durée inférieure à la durée T. Un troisième tri, plus fin, visant à séparer les signaux et à en optimiser le rapport signal à bruit est effectué par filtrage numérique, en fonction de leur fréquence centrale et de leur largeur de bande. Cette opération s'effectue par exemple à l'aide de filtres FFT ou de filtres 5 polyphasés de plusieurs largeurs différentes AF1. Par exemple, pour un récepteur de bande passante AF = 100 MHz, on peut disposer d'une batterie de filtres de 25 MHz, d'une autre batterie de filtres de 12,5 MHz et d'une troisième batterie de filtres de 1,56 MHz pour le troisième tri, ces filtres pouvant être cascadés ou non. 10 Ainsi, si la période d'échantillonnage du signal est égale à Te = 1/2 OF, on peut former l'ensemble de filtres suivant : - 4 filtres de bande AF/4 = 25 MHz avec une période de répétition de 4Te = 40 ns ; - 8 filtres de bande AF/8 = 12,5 MHz avec une période de répétition de 15 8Te = 80 ns ; - 64 filtres de bande AF/64 = 1,56 MHz avec une période de répétition de 64Te = 640 ns. La figure 8 illustre un exemple de réalisation fonctionnel dans le cas d'une 20 application à un intercepteur passif. Le dispositif est représenté en sorties 80 des M récepteurs. Les signaux de sorties des M récepteurs sont issus des N éléments rayonnants, via la matrice de commutation 61. Ils sont numérisés par le convertisseur analogique-numérique à M voies 62. Chaque voie de sortie de ce convertisseur est reliée à un banc de filtres 81. Ainsi, à chaque 25 élément rayonnant d'ordre i, connecté au dispositif de réception à un instant T, correspond un ensemble de filtres 811 de fréquence central Fm et de largeur de bande AF1. Un banc de filtre 81 comporte par exemple l'ensemble de filtres décrit ci-dessus. Les termes de corrélation entre les signaux reçus des éléments rayonnants 30 d'ordre i et d'ordre j, soit entre les voies i et j, sont obtenus en calculant la valeur moyenne du produit 83 de la sortie de chaque filtre Smd,i,k connecté à la voie i à l'instant kT avec le conjugué de la sortie du filtre de même fréquence centrale et de même bande passante filtres connecté à la voie j à l'instant kT. Cela se traduit par la relation suivante : 35 où : R,J,,,,j(pT)= E(S * p+R = (10) k=p représente la sortie du filtre de fréquence centrale Fm, de largeur de bande AF1, connecté à la voie i à l'instant kT ; S*rn,l,j,k représente le conjugué de la sortie du filtre de fréquence centrale Fm, de largeur de bande AF1, connecté à la voie j à l'instant kT , - RT représente l'intervalle de temps sur lequel est effectuée la moyenne ; - pT représente le temps auquel est effectuée la moyenne ; (pT) représente la corrélation entre les signaux reçus des voies d'indice i et d'indice j pour les filtres de fréquence centrale Fm et de bande AF1 , à l'instant pT. Les signaux étant supposés à bande étroite et constants en amplitude et en 15 direction d'arrivée pendant le temps d'observation, les termes de corrélation sont indépendants du temps, soit : R,J,n,J(pT)= (11) L'ensemble des termes de corrélation est ainsi estimé 82 selon un ensemble de séquences successives pendant lesquelles le dispositif de réception à M 20 voies est connecté successivement aux différents éléments rayonnants constituant l'antenne. L'ensemble des résultats est alors sommé de façon cohérente 83 pour chaque direction de l'espace et pour chaque filtre. Si on pose Ri,m,, , la corrélation des signaux provenant des éléments 25 rayonnants j et j+i pour les filtres de fréquence centrale Fm et de bande AF1, l'estimation de la densité spectrale de puissance spatiale 84 en sortie des filtres de fréquence centrale Fm et de bande AF1 s'écrit pour la direction Ok : N bijik S(9 k) = Rimiexp N (12) i=0 Le résultat de cette estimation est ensuite utilisé pour détecter puis localiser angulairement les signaux reçus 85, selon les principes classiques en traitement radar. Les moyens de traitement réalisent notamment les bancs de filtre 81, les corrélations 82, les intégrations cohérentes 83 et l'estimation de la densité spectrale dans les différentes directions liées aux éléments d'antenne. Typiquement, dans l'exemple d'une antenne à 16 voies selon le schéma de la figure 6 qui utilise un récepteur à 4 voies simultanées, 5 séquences sont nécessaires pour calculer tous les termes de corrélation permettant la 10 formation de faisceaux, comme le montrent les estimations de ces termes Rx(i) dans les séquences 51, 52, 53, 54, 55 décrites précédemment. Si par ailleurs on s'intéresse à l'analyse de signaux de durée supérieure à 100 ps, on peut allouer 20 ps par séquence, c'est-à-dire estimer chaque terme de corrélation pendant un temps total RT=20 ps. 15 Si par exemple, on réalise la séparation angulaire dans chaque filtre de bande AF/4 = 25 Mhz, on dispose du signal en sortie de chaque filtre d'un échantillon toutes les 40 ns et la moyenne des termes produits peut être effectuée sur 500 échantillons successifs. En plus de la séparation angulaire, ce procédé permet une intégration 20 temporelle cohérente du signal reçu pendant un temps voisin de la durée totale de ce signal, ce qui est très intéressant pour détecter les signaux à faible puissance crête (de type LPI) et notamment de modulation quelconque. Avantageusement, l'intégration des différents termes de corrélation peut être adaptée pour détecter des signaux à impulsions de 25 longue durée, par exemple supérieure à 100 ps, de faible puissance crête, par exemple inférieure à quelques watts, et de modulation inconnue quelconque. REVENDICATIONS1. Dispositif de détection de signaux électromagnétiques comportant une antenne de réception réseau (10) à N éléments rayonnants (1) , caractérisé en ce que, comportant M voies de réception (2, 3) en aval de l'antenne de 5 réception, M étant inférieur à N, ledit dispositif comporte au moins : des moyens de commutation (61) des M voies de réception sur les éléments rayonnants selon des cycles de séquences successives (51, 52, 53, 54, 55), M éléments rayonnants étant connectés auxdites voies de réceptions à chaque séquence, un même élément rayonnant 10 (50) dit de référence étant connecté auxdites voies de réception pour toutes les séquences, un cycle étant achevé lorsque tous les éléments rayonnants ont été connectés au moins une fois à une desdites voies de réception ; des moyens de traitement (63, 82, 83, 84) effectuant pour chaque 15 séquence l'estimation des corrélations spatiale deux à deux du signal reçu sur la voie de référence et des signaux reçus sur les autres M-1 voies de réception puis estimant la densité spectrale de puissance spatiale dans N directions d'arrivée à partir d'une somme cohérente des N termes de corrélations ainsi obtenus. 20 2. Dispositif selon la revendication 1, caractérisé en ce que le réseau d'éléments rayonnant étant linéique, l'élément de référence (50) est le premier éléments, les M premiers éléments étant connectés aux voies de réception dans la première séquence, puis les M-1 éléments rayonnant étant 25 connectés aux voies de réceptions dans la deuxième séquence, et ainsi de suite. 3. Dispositif selon la revendication précédente, caractérisé en ce que les moyens de traitement estiment les corrélations spatiales après avoir effectué 30 l'élimination des signaux de courte durée stationnaire de durée inférieure à une durée T donnée, l'élimination des signaux de courte durée stationnaire étant effectuée en multipliant le signal reçu à un t, x(t), par un signal de la forme x(t)x(t-r)/x(t)*2, où r est un retard choisi pour supprimer les signaux de durée inférieure à la durée r. 4. Dispositif selon l'une quelconque des revendications précédentes, caractérisé en ce que, dans chaque séquence, les dits moyens de traitement effectuent l'estimation de ladite corrélation spatiale en multipliant directement le signal reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne. 5. Dispositif selon l'une quelconque des revendications 1 à 3, caractérisé en ce que, dans chaque séquence, les dits moyens de traitement effectuent l'estimation de ladite corrélation spatiale en multipliant le signal reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne, après séparation fréquentielle desdits signaux reçus. 6. Dispositif selon la revendication 5, caractérisé en ce que, dans chaque séquence, les dits moyens de traitement effectuent l'estimation de ladite corrélation spatiale en multipliant le signal reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne, l'estimation de corrélation étant effectuée en sortie de bancs de filtres (81) de différentes largeurs et de différentes fréquences centrales. 7. Dispositif selon la revendication 6, caractérisé en ce que les filtres sont réalisés par transformées de Fourrier numérique d'ordres différents. 8. Dispositif selon la revendication 6, caractérisé en ce que les filtres sont réalisés par filtres polyphasés d'ordres différents. 9. Dispositif selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que ledit dispositif est un radar. 10. Dispositif selon la revendication 9, caractérisé en ce que, dans chaque séquence, les dits moyens de traitement effectuent l'estimation de ladite corrélation spatiale en multipliant le signal reçu sur ledit élément rayonnant de référence par le signal reçu conjugué sur un autre élément rayonnant de l'antenne, après séparation temporelle desdits signaux reçus. 11. Dispositif selon la revendication 10, caractérisé en ce que les dits moyens de traitement obtiennent la séparation temporelle par filtrage adapté et échantillonnage des impulsions reçues par ledit radar.5