FR 2646512 A1 19901102 FR 8905677 A 19890428 PROCEDE ET DISPOSITIF D'ANALYSE SPECTRALE A HAUTE RESOLUTION POUR SIGNAUX DE SYSTEME DE DETECTION ELECTROMAGNETIQUE NON STATIONNAIRES La présente invention se rapporte à un procédé et à un dispositif d'analyse spectrale à haute résolution pour signaux de système de détection électromagnétique non stationnaires o ces signaux subissent une première analyse spectrale glissante classique sur une première durée d'analyse correspondant à la durée de validité du modèle d'évolution choisi pour lesdits signaux non stationnaires. Par "système de détection électromagnétique", on entendra dans toute la description et les revendications tout système permettant la détection d'une cible par analyse, en particulier en fréquence, des signaux électromagnétiques réfléchis ou émis par la cible. De tels systèmes comprennent notamment les radars ainsi que les barrières électromagnétiques dans lesquelles un émetteur émet un faisceau large vers un récepteur situé à distance, les cibles étant détectées par les perturbations qu'elles induisent, par leur passage dans le faisceau, dans la transmission entre récepteur et émetteur. Dans le traitement de certains signaux issus de tels systèmes de détection électromagnétique, les performances de détection et de dénombrement de cibles en situation multicibles dépendent de la qualité de l'analyse spectrale. Ceci est particulièrement vrai pour les barrières électromagnétiques o le seul critère de détection est obtenu par analyse spectrale. Un des facteurs essentiels pour l'obtention de bonnes performances est alors le pouvoir de résolution fréquentielle. Depuis un certain nombre d'années, on a mis au point et perfectionné des méthodes d'analyse spectrale à haute résolution pour des signaux radars stationnaires. En effet, les méthodes classiques utilisant la transformée de Fourier sont intrinsèquement limitées en résolution par l'inverse du temps 2646512. d'analyse T du signal. Dans certaines conditions de stationnarité des signaux, on a donc proposé des méthodes d'analyse spectrale telles que celles basées sur une modélisation autorégressive du processus étudié. Cependant l'hypothèse de stationnarité conduit dans la pratique à une limitation du temps d'analyse et donc à une limitation de la résolution due aux non-stationnarités du signal. L'invention a pour objet de remédier à ces inconvénients en prévoyant un traitement simple adapté aux signaux non stationnaires. Un autre objet de l'invention est un procédé et un - dispositif d'analyse spectrale à haute résolution consistant à rendre les signaux sensiblement stationnaires avant de leur appliquer une analyse spectrale classique. Selon l'invention, il est donc prévu un procédé d'analyse spectrale à haute résolution pour signaux de système de détection électromagnétique non stationnaires, dans lequel on réalise une première analyse spectrale glissante de type connu desdits signaux sur une première durée d'analyse correspondant 2r à la durée de validité du modèle d'évolution choisi pour lesdits signaux non stationnaires, ledit procédé étant caractérisé en ce qu'il consiste à: - effectuer un pistage Doppler, en fonction du temps, de la fréquence instantanée obtenue par ladite première analyse; - estimer pour une piste choisie Ies paramètres de la loi d'évolution de la fréquence Doppler retenue; - déterminer à partir desdits paramètres un signal synthétique de fréquence instantanée opposée à la fréquence instantanée estimée théorique correspondant auxdits paramètres; - - stationnariser lesdits signaux fournis par ledit système par démodulation à l'aide dudit signal synthétique; -effectuer une seconde analyse spectrale de type connu dudit signal stationnarisé, sur une durée d'analyse longue pour obtenir une résolution améliorée; et -effectuer un dénombrement des cibles dans le signal analysé selon ladite seconde analyse par comparaison à un seuil déterminé. Selon un autre aspect de l'invention il est également prévu un dispositif d'analyse spectrale à haute définition mettant en oeuvre le procédé ci-dessus. L'invention sera mieux comprise et d'autres caractéristiques et avantages apparaîtront à l'aide de la description ci-après et des dessins joints o - la figure I est un schéma par blocs du dispositif d'analyse spectrale selon l'invention; - la figure 2 est un diagramme temps-fréquence de deux motifs correspondant à des signaux non séparés par transformée de Fourier classique; - la figure 3 représente le spectre d'un des motifs à signaux non séparés; et - la figure 4 représente le spectre d'un des motifs après traitement selon l'invention. Comme on l'a déjà expliqué ci-dessus, lorsque les 2G signaux fournis par un système de détection électromagnétique sont fortement non stationnaires, c'est-à-dire que la fréquence Doppler évolue rapidement au cours du temps, la résolution en fréquence des procédés d'analyse spectrale classiques est 1/ t dans le cas de la transformée de Fourier. Or t, durée d'analyse, est limitée par le temps pendant lequel l'hypothèse faite sur la trajectoire de la cible est vraie. t étant fortement limitée du fait de la non-stationnarité du signal, la résolution est également fortement limitée. Pour simplifier et préciser l'explication, on va supposer, sans que cela soit en rien limitatif de l'invention, que l'on est en présence de deux cibles multiples comprenant chacune cinq cibles proches se déplaçant ensemble. Si l'on suppose aussi que la fréquence Doppler des deux groupes de cibles évolue rapidement, l'analyse spectrale classique, par exemple par transformée de Fourier ne permet pas de séparer les cibles composant chacune des deux cibles multiples. La figure 1 représente le schéma du dispositif selon l'invention. Le système de détection électromagnétique i fournit des signaux s(t) à un dispositif d'analyse spectrale classique 2 à durée d'analyse courte At. Par analyse glissante, on dispose de spectres successifs au cours du temps qui sont représentés sur le diagramme de la figure 2 o on retrouve deux motifs formés des cinq signaux non séparés. Ces spectres successifs 1O sont envoyés à un système de pistage Doppler multicible 3 qui permet d'analyser séparément chaque piste correspondant à un motif. Un tel système de pistage est décrit par exemple dans le brevet français N 80 25039. Une unité de gestion 4 permet de sélectionner une des pistes et de commander un dispositif d'estimation des paramètres de la cible (multicible en fait) choisie qui, à partir de la loi d'évolution de la fréquence Doppler au cours du temps fournie par le système 3, permet de modéliser cette loi d'évolution et de caleuler les paramètres correspondants de la cible. Le modèle utilisé peut être un modèle simple mais général de modulation linéaire de fréquence du type fi(t) = fi + (aXit o fi(t) est la fréquence instantanée du signal, fi est une valeur donnée initiale et a. est la pente de modulation. Dans les zones o le modèle linéaire est valide, le dispositif 5 estime Ies valeurs fi et a,' par exemple par la méthode des moindres carrés. On peut aussi prendre un modèle plus complet et spécifique du système de détection considéré: dans le cas d'une barrière électromagnétique par exemple, la loi d'évolution de la fréquence Doppler peut s'exprimer simplement en fonction de paramètres de la cible tels que vitesse de la cible, angle de la trajectoire par rapport à la ligne joignant émetteur et récepteur, position du point de coupure de cette ligne par la trajectoire, altitude de la cible ou distance de l'origine de la trajectoire audit point de coupure. Quel que soit le modèle 26465 1 2 choisi, si on appelle e les paramètres à déterminer et fth (t, ô) la loi d'évolution théorique (modèle) de la fréquence Doppler, le dispositif 5 estime les paramètres en minimisant l'expression t22 [f(t) - fth (t,)12 dt t1 té 1 mesurée o f(t) est la loi d'évolution fournie par le système de pistage mesurée 3 entre les instants t1 et t2 Ces paramètres e sont envoyés à un système 6 de restitution d'un signal synthétique qui calcule un signal s'(t) de fréquence instantanée opposée à la fréquence théorique du modèle estimé. Par exemple dans le cas du modèle simple, le signal synthétique est de la forme - 2 S - J( i t2/ s'(t) = e- 211j(fit t /2) ou, dans le cas plus général, de la forme s'(t) = e-j (p(t) avec t - Q(t) = 2r f (u. 5) du t0 Grâce à ce signal synthétique ainsi calculé, on peut alors stationnariser le signal s(t) reçu par démodulation par s'(t), c'est-à-dire par simple multiplication. Pour cela, on stocke le signal s(t) reçu dans une mémoire 7 et on le restitue en même temps que le signal synthétique s' (t) à un multiplicateur ou mélangeur 8. Celui-ci fournit donc un signal stationnarisé s"(t) qui est envoyé à un dispositif d'analyse spectrale classique 9 (par exemple par transformée de Fourier) qui peut effectuer l'analyse sur une durée d'analyse T longue, donc avec une haute résolution en fréquence. Les signaux de sortie de l'analyse spectrale sont envoyés à un dispositif à seuil 10 qui permet le dénombrement des cibles constitutives du motif (multicible) analysé. Ce dénombrement et l'affichage du résultat sont effectués par un circuit 11. A titre d'exemple illustratif, la figure 3 montre le spectre d'un motif fourni par le dispositif d'analyse 2: on ne peut pas séparer les cibles constituant le motif. Après traitement selon l'invention, le spectre fourni par le dispositif d'analyse spectrale 9 est représenté pour le même motif sur la figure 4. On voit distinctement les pics correspondant aux cinq cibles constituant le motif, ce qui permet avec le dispositif à seuil 10 de bien dénombrer ces cibles. Ainsi, le procédé selon l'invention permet une amélioration sensible de la résolution en fréquence pour le traitement de signaux non stationnaires et il apparaît bien adapté aux contraintes de traitement en temps réel du signal radar notamment. Le dispositif à haute résolution prévu selon l'invention est tout à fait simple et robuste. Bien entendu, l'exemple de réalisation décrit n'est nullement limitatif de l'invention. REVENDICATIONS 1. Procédé d'analyse spectrale à haute résolution pour signaux de système de détection électromagnétique non stationnaires, dans lequel on réalise une première analyse spectrale glissante de type connu desdits signaux sur une première durée d'analyse (Lt) correspondant à la durée de validité du modèle d'évolution choisi pour lesdits signaux non stationnaires, ledit procédé étant caractérisé en ce qu'il consiste à: - effectuer un pistage Doppler, en fonction du temps, de la fréquence instantanée obtenue par ladite première analyse -estimer pour une piste choisie les paramètres (E) de la loi d'évolution de la fréquence Doppler retenue; - déterminer à partir desdits paramètres un signal synthétique de fréquence instantanée opposée à la fréquence instantanée estimée théorique correspondant auxdits paramètres; - stationnariser lesdits signaux fournis par ledit systènie par démodulation à l'aide dudit signal synthétique; - effectuer une seconde analyse spectrale de type connu dudit signal stationnarisé, sur une durée d'analyse (T) longue pour obtenir une résolution améliorée; et effectuer un dénombrement des cibles dans le signal analysé selon ladite seconde analyse par comparaison à un seuil déterminé. 2. Procédé selon la revendication 1, caractérisé en ce que dans les étapes consistant à estimer les paramètres et à déterminer un signal synthétique, la loi d'évolution de la fréquence Doppler est une loi de modulation linéaire de fréquence de la forme fi(t) '= fi + ai.t pour la piste de rang i, o fi est une valeur initiale de fréquence et ai la pente de modulation, ces deux valeurs fi et ai étant estimées à partir desdits paramètres (E). 3. Dispositif d'analyse spectrale à haute résolution pour signaux de système de détection électromagnétique non. stationnaires, mettant en oeuvre le procédé selon l'une des revendications 1 ou 2, caractérisé en ce qu'il comprend: - un premier dispositif d'analyse spectrale (2) glissante à durée d'analyse (At) courte recevant lesdits signaux non stationnaires; - un système de pistage Doppler (3) multicibles recevant les signaux dudit premier dispositif (2) et fournissant au moins une courbe de variation fréquence instantanée-temps correspondant à une piste vraisemblable; - un dispositif (5) d'estimation des paramètres (0) d'une cible correspondant à ladite courbe; - une unité de gestion (4) pour sélectionner une des pistes fournies par ledit système de pistage (3) et commander ledit dispositif d'estimation (5); - des moyens (6) de restitution d'un signal synthétique de fréquence instantanée opposée à partir desdits paramètres - des moyens de mémorisation (7) pour stocker lesdits signaux 29 non stationnaires fournis par le système de détection (1); - des moyens (8) de démodulation desdits signaux stockés dans lesdits moyens de mémorisation (7) par ledit signal synthétique pour obtenir un signal stationnarisé; - un second dispositif d'analyse spectrale (9) pour analyser ledit signal stationnarisé; et -des moyens de comparaison à seuil (10) pour extraire des signaux fournis par ledit second dispositif (9) des informations sur le nombre de cibles ayant contribué à fournir les signaux correspondant à ladite piste sélectionnée. 4. Dispositif selon la revendication 3, caractérisé en ce que lesdits premier et second dispositifs d'analyse spectrale (2, 9) comprennent chacun un calculateur de transformée de Fourier. 26 4 6 5 1 2 5. Dispositif selon l'une des revendications 3 ou 4, caractérisé en ce que lesdits moyens de démodulation (8) comprennent un multiplicateur de signaux.