FR2998973A1 - Procede de determination des caracteristiques d un signal radar en presence d interferences - Google Patents

Procede de determination des caracteristiques d un signal radar en presence d interferences Download PDF

Info

Publication number
FR2998973A1
FR2998973A1 FR1203239A FR1203239A FR2998973A1 FR 2998973 A1 FR2998973 A1 FR 2998973A1 FR 1203239 A FR1203239 A FR 1203239A FR 1203239 A FR1203239 A FR 1203239A FR 2998973 A1 FR2998973 A1 FR 2998973A1
Authority
FR
France
Prior art keywords
pulses
signal
agility
repetition period
aprim
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
FR1203239A
Other languages
English (en)
Other versions
FR2998973B1 (fr
Inventor
Jean Pierre Dubost
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Thales SA
Original Assignee
Thales SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Thales SA filed Critical Thales SA
Priority to FR1203239A priority Critical patent/FR2998973B1/fr
Publication of FR2998973A1 publication Critical patent/FR2998973A1/fr
Application granted granted Critical
Publication of FR2998973B1 publication Critical patent/FR2998973B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/38—Jamming means, e.g. producing false echoes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

Procédé de détermination des caractéristiques d'un signal comprenant une pluralité d'impulsions, consistant à déterminer au moins la période de répétition des impulsions (PRI) et comportant les étapes: - Calculer (302) la période de répétition moyenne PRIm des impulsions en présence d'interférences, - Calculer (303) l'écart type ΔPRIm de la période de répétition moyenne PRIm des impulsions en présence d'interférences, - Définir (304) un domaine d'agilité [PRImin ; PRImax] de la période de répétition des impulsions (PRI) de la façon suivante PRImin = (PRIm + k.ΔPRIm)/√α et PRImax = √α.(PRIm - k.ΔPRIm), avec α un facteur d'excursion strictement positif et inférieur ou égal à 2 et k un nombre strictement positif, - Calculer (305) les différences de temps d'arrivée (DTOA) entre deux fronts avants ou deux fronts arrières de deux impulsions successivement reçues, - - Eliminer (306) les valeurs des différences de temps d'arrivée qui sont hors du domaine d'agilité [PRImin ; PRImax].

Description

Procédé de détermination des caractéristiques d'un signal radar en présence d'interférences Le domaine technique de l'invention concerne les systèmes de 5 surveillance de l'environnement électromagnétique et plus précisément les systèmes de détection et d'identification d'un signal radar. L'invention a pour objet un procédé de détermination des caractéristiques d'un signal radar en présence d'interférences causées par des systèmes environnants rayonnants qui provoquent une gêne pour la 10 capture et l'analyse des impulsions dont le signal est composé. L'invention s'applique avantageusement à un système d'identification de radar embarqué sur une plateforme multi-systèmes qui comprend par ailleurs d'autres systèmes rayonnants tels un système de communications, d'aide à la navigation ou un dispositif de brouillage. 15 L'invention s'applique également pour un système d'identification de radar localisé dans un environnement électromagnétique soumis à des interférences provenant par exemple d'infrastructures de communications telles une antenne relai et pour lequel un indicateur de la présence d'interférences est disponible. 20 Les systèmes d'identification de radar ont pout but d'analyser les signaux émis par un radar situé dans leur portée pour caractériser le mode d'émission du radar. Un signal radar est le plus souvent formé d'une succession d'impulsions émises périodiquement. La caractérisation d'un 25 mode d'émission radar passe par l'estimation de la fréquence du signal, de la période de répétition des impulsions et de la durée temporelle d'une impulsion, encore appelée longueur d'impulsion. Ces paramètres peuvent varier au cours du temps. L'estimation du paramètre de période d'impulsion nécessite la mesure de l'instant de réception de chaque impulsion qui est 30 obtenue indirectement en mesurant la différence de temps entre la réception de deux impulsions successives. L'estimation du paramètre de longueur d'Impulsion nécessite la mesure des instants de début et de fin d'impulsions. Lorsque l'environnement électromagnétique est soumis à des 5 interférences provenant d'autres systèmes situés à portée du système d'identification de radar, le signal radar reçu est dégradé de sorte que certaines impulsions sont totalement ou partiellement atténuées ou que leur détection est rendue impossible car non fiable. Dans un tel cas, l'estimation des caractéristiques du signal radar et donc l'identification du mode radar 10 deviennent moins fiables voire impossibles. Le problème de la vulnérabilité aux interférences est notamment présent dans le cas de plateformes multi-systèmes qui embarquent, dans un espace réduit, à la fois un système d'identification radar et d'autres systèmes rayonnants qui sont source de gênes. 15 Un problème similaire existe lorsque les interférences proviennent d'infrastructures indépendantes du système d'identification de radar. Dans le cas d'une plateforme multi-systèmes, une solution connue pour limiter la présence d'interférences, consiste à faire coopérer les 20 différents systèmes entre eux afin d'interdire aux systèmes interférents d'émettre lorsque le système d'identification cherche à analyser un signal radar, autrement dit à imposer des silences d'émission aux systèmes interférents. Une telle solution présente de nombreuses limitations. Tout d'abord, 25 lorsque le nombre de systèmes interférents est important et que le nombre de radars à identifier l'est également, il devient difficile de faire coopérer l'ensemble des systèmes de sorte à imposer des intervalles temporels de durée suffisante pendant lesquels aucun signal interférent n'est émis. En outre, les systèmes d'identification de radar peuvent être couplés à 30 un brouilleur qui a pour but, lorsqu'un radar est détecté et identifié comme fonctionnant selon un mode donné, de neutraliser ce radar en émettant un signal de brouillage qui est caractérisé par une forte puissance et donc un niveau d'interférences important. Même lorsqu'un signal radar est brouillé, il peut être souhaitable de continuer à identifier les caractéristiques de ce signal car celui-ci peut évoluer vers un mode différent du mode initial. Il y a donc un intérêt à poursuivre en permanence l'analyse du signal radar. Ce cas de figure est notamment adapté aux brouillages de type non synchrone. L'invention apporte une solution au problème de la caractérisation d'un signal radar en présence d'interférences qui ne nécessite pas d'imposer 10 aux systèmes environnants d'interrompre l'émission d'un signal potentiellement interférant. L'invention permet d'estimer les paramètres de période de répétition des impulsions et de longueur des impulsions d'un signal radar. Elle s'applique à la caractérisation de tout type de signaux radar composés d'une 15 pluralité d'impulsions successivement émises. L'invention s'applique en particulier pour des radars à balayage ou des radars pointés. L'invention a pour objet un procédé de détermination, en présence d'interférences, des caractéristiques d'un signal comprenant une pluralité 20 d'impulsions émises, ledit procédé consistant à déterminer au moins la période de répétition des impulsions, ledit procédé étant caractérisé en ce qu'il comporte les étapes suivantes, exécutées pour une portion dudit signal de durée donnée TE : - Calculer la période de répétition moyenne PRIm des impulsions en 25 présence d'interférences, - Calculer l'écart type APRIm de la période de répétition moyenne PRIm des impulsions en présence d'interférences, Définir un domaine d'agilité [PRIn,m ; PRImax] de la période de répétition des impulsions de la façon suivante PRImm = (PRIm + 30 k.APRIm)Na et PRImax = -Va.(PRIm - k.APRIm), avec a un facteur d'excursion strictement positif et inférieur ou égal à 2 et k un nombre strictement positif, Calculer, pour l'ensemble des impulsions comprises dans ladite portion de signal et dont le front avant ou le front arrière est non tronqué, les différences de temps d'arrivée entre deux fronts avants ou deux fronts arrières de deux impulsions successivement reçues, Eliminer les valeurs des différences de temps d'arrivée qui sont hors du domaine d'agilité [PRImin ; PRImax], la ou les valeurs restantes donnant une estimée de la période de répétition des impulsions. L'invention a aussi pour objet un procédé de détermination, en présence d'interférences, des caractéristiques d'un signal comprenant une pluralité d'impulsions émises, ledit procédé consistant à déterminer au moins la période de répétition des impulsions, ledit procédé étant caractérisé en ce qu'il comporte les étapes suivantes, exécutées pour une portion dudit signal de durée donnée TE : Calculer une estimée de la période de répétition moyenne PRIm des impulsions en présence d'interférences, Calculer une estimée de l'écart type APRIm de la période de répétition moyenne PRIm des impulsions en présence d'interférences, Si la variance APRIm2 ou l'écart type APRIm ou l'écart type normalisé par la moyenne APRIm/PRIm est inférieur à un seuil de précision Pmax prédéterminé, alors exécuter les étapes suivantes : Définir un domaine d'agilité [PRImin ; PRImax] de la période de répétition des impulsions, de la façon suivante PRImm = (PRIm + k.APRIm)hla et PRImax = '\ia.(PRIm - k.APRIm), avec a un facteur d'excursion strictement positif et inférieur ou égal à 2 et k un nombre strictement positif, Calculer, pour l'ensemble des impulsions comprises dans ladite portion de signal et dont le front avant ou le front arrière est non tronqué, les différences de temps d'arrivée entre deux fronts avants ou deux fronts arrières de deux impulsions successivement reçues, Eliminer les valeurs des différences de temps d'arrivée qui sont hors du domaine d'agilité [PRImin ; PRImax], la ou les valeurs restantes donnant une estimée de la période de répétition des impulsions.
Selon un aspect particulier de l'invention, ledit seuil de précision Prnax est déterminé en fonction de l'excursion PRImax/PRImin minimale souhaitée du domaine d'agilité. Selon un autre aspect particulier de l'invention, la période de répétition moyenne PRI' des impulsions est calculée à l'aide de la relation suivante : PRIm = (1-Ti).TE/n, avec Ti un indicateur du taux de présence d'interférences pendant la durée TE et n le nombre d'impulsions, comprises dans ladite portion, dont le front avant est non tronqué ou le nombre d'impulsions, comprises dans ladite portion, dont le front arrière est non tronqué ou une combinaison du nombre nav d'impulsions dont le front avant est non tronqué et du nombre nar d'impulsions dont le front arrière est non tronqué sous la forme n=(a.nav+[3.nar)/(a+13), où a et i3 sont deux entiers positifs ou nuls. Selon un autre aspect particulier de l'invention, l'écart type APRI' de la période de répétition moyenne PRI, des impulsions est calculé à l'aide de la relation suivante : APRI' = (Ti/n)in.PRIm, Dans une variante de réalisation de l'invention, l'écart type APRIm de la période de répétition moyenne PRI, des impulsions est pondéré par le coefficient : 2 /22 lorsque n est une combinaison du nombre a+fl d'impulsions dont le front avant est non tronqué et du nombre na, d'impulsions dont le front arrière est non tronqué sous la forme 30 n=(a.nav+13.nar)/(a-Ff3), où a et 13 sont deux entiers positifs ou nuls. nav Selon un autre aspect particulier de l'invention, le nombre k est pris égal à 1. Selon un autre aspect particulier de l'invention, ledit facteur correctif a est pris égal à 2.
Dans une variante de réalisation de l'invention, les différences de temps d'arrivée sont en outre calculées pour les impulsions dont le front avant ou le front arrière est tronqué et sont dans ce cas assorties d'une fenêtre d'incertitude de plus ou moins la longueur d'impulsion moyenne Llm. Dans une variante de réalisation de l'invention, lorsqu'aucune valeur 10 n'est conservée après l'étape d'élimination des valeurs des différences de temps d'arrivée hors domaine d'agilité [PRImin ; PRImax], une estimée de la période de répétition des impulsions est obtenue en recherchant un sous-multiple, situé dans ledit domaine d'agilité [PRImin PRImaxi, d'au moins deux mesures successives de différences de temps d'arrivée. 15 Dans une variante de réalisation de l'invention, ladite recherche d'un sous-multiple est effectuée en exécutant au moins les étapes suivantes : Définir deux variables X,Y et les initialiser à deux mesures successives PPRI1,PPRI2 de différences de temps d'arrivée, Calculer l'opération modulo des deux variables Mod(X,Y), 20 Si le résultat de l'opération modulo est situé dans ledit domaine d'agilité D1=[(PRIm + APRIm)hia ; .Nia.(PRIm - APRIn,)], l'estimée de la période de répétition des impulsions est égale à ce résultat, Si le résultat de l'opération modulo est situé hors du domaine d'agilité, réaffecter la variable X à la seconde mesure PPRI2 et la 25 variable Y au résultat de l'opération modulo Mod(PPRIi, PPRI2) entre la première mesure PPRIi et la seconde mesure PPRI2, Si la valeur de X est différente de la valeur de Y, alors réitérer les étapes précédentes, Si la valeur de X est égale à la valeur de Y, alors calculer 30 l'ensemble des sous-multiples de la variable X, Exécuter l'étape de test suivante : L Si un seul desdits sous-multiples est situé dans ledit domaine d'agilité D1, l'estimée de la période de répétition des impulsions est égale audit seul sous-multiple, ii. Si plusieurs desdits sous-multiples sont situés dans ledit domaine d'agilité D1, calculer le résultat de l'opération modulo Mod(X,PPRI3) entre la variable X et une troisième mesure successive de différences de temps d'arrivée PPRI3, l'estimée de la période de répétition des impulsions étant égale à ce résultat.
Dans une variante de réalisation de l'invention, ledit domaine d'agilité D1 est remplacé par un second domaine d'agilité D2= RPRIm - APRIm)hla ; -Va.(PRIm + APRI,)]. Dans une variante de réalisation de l'invention, l'étape de test est modifiée de la façon suivante: Si un seul desdits sous-multiples est situé dans le premier domaine d'agilité D1=[(PRIm + APRIm)ka ; -Va.(PRIm - APRIm)], et qu'aucun autre desdits sous-mutliples n'est situé dans un second domaine d'agilité D2= RPRIm - APRIm)hla ; -Va.(PRIm + APRIm)] alors l'estimée de la période de répétition des impulsions (est égale audit seul sous-multiple, Si plusieurs desdits sous-multiples sont situés dans ledit domaine d'agilité D1 et qu'aucun autre desdits sous-mutliples n'est situé dans un second domaine d'agilité D2= RPRIm - APRIm)Not ; -Va.(PRIm + APRIm)], calculer le résultat de l'opération modulo Mod(X,PPRI3) entre la variable X et une troisième mesure successive de différences de temps d'arrivée PPRI3, l'estimée de la période de répétition des impulsions étant égale à ce résultat. Si un seul desdits sous-multiples est situé dans le second domaine d'agilité D2= RPRIm - APRIm)hla ; -Va.(PRIm + APRIm)] alors l'estimée de la période de répétition des impulsions est égale audit seul sous-multiple, - Si plusieurs desdits sous-multiples sont situés dans le second domaine d'agilité D2= RPRIm - APRIm)ka ; \la.(PRIm + APRIm)], calculer le résultat de l'opération modulo Mod(X,PPRI3) entre la variable X et une troisième mesure successive de différences de temps d'arrivée PPRI3, l'estimée de la période de répétition des impulsions étant égale à ce résultat. Dans une variante de réalisation de l'invention, ledit procédé consiste à déterminer en outre la longueur d'impulsion et comporte les étapes suivantes : Si ladite portion de signal comporte des impulsions non tronquées, utiliser les longueurs de ces impulsions comme estimée de la longueur d'impulsion, Si ladite portion de signal ne comporte que des impulsions tronquées, calculer une estimée de la valeur moyenne Llm de la longueur d'impulsion. Dans une variante de réalisation de l'invention, ledit procédé consiste à déterminer en outre la longueur d'impulsion et comporte les étapes suivantes : - Si ladite portion de signal comporte des impulsions non tronquées, utiliser les longueurs de ces impulsions comme estimée de la longueur d'impulsion, - Calculer un indicateur de précision ALIm de la valeur moyenne Llm de la longueur d'impulsion, - Si ledit indicateur est inférieur à un seuil de précision prédéterminé, alors exécuter l'étape suivante : Si ladite portion de signal ne comporte que des impulsions tronquées, calculer une estimée de la valeur moyenne Llm de la longueur d'impulsion. Ledit indicateur est soit un écart type ALlm, soit un écart type 30 normalisé par la moyenne ALIm/Llm soit une variance ALIm2.
Selon un aspect particulier de l'invention, ledit indicateur est un écart type ALln, calculé à l'aide de l'une des deux relations suivantes : ALI. =.ii--qELI(i)1 OU 4L1=\/21-''3 avec Zu(i)une mesure de la n somme des longueurs des impulsions reçues comprises dans ladite portion 5 et n le nombre d'impulsions, comprises dans ladite portion, dont le front avant est non tronqué ou le nombre d'impulsions, comprises dans ladite portion, dont le front arrière est non tronqué ou une combinaison du nombre nav d'impulsions dont le front avant est non tronqué et du nombre har d'impulsions dont le front arrière est non tronqué sous la forme 10 n=(a.nav+13.nar)/(a-Fr3), où a et 13 sont deux entiers positifs ou nuls. Selon un aspect particulier de l'invention, la valeur moyenne LI,, de la longueur d'impulsion est calculée en divisant la somme des longueurs des impulsions comprises dans ladite portion par le nombre d'impulsions, comprises dans ladite portion, dont le front avant est non tronqué ou par le 15 nombre d'impulsions, comprises dans ladite portion, dont le front arrière est non tronqué ou par une combinaison du nombre nav d'impulsions dont le front avant est non tronqué et du nombre nar d'impulsions dont le front arrière est non tronqué sous la forme n=(a.nav-FrInar)/(a+13), où a et p sont deux entiers positifs ou nuls. 20 Selon un aspect particulier de l'invention, un indicateur d'interférences synchrone dudit signal est utilisé pour déterminer si une impulsion est tronquée ou non. Dans une variante de réalisation de l'invention, ledit procédé comporte une étape préalable d'interception et de recueillement de ladite portion de 25 signal pendant une durée donnée TE. L'invention a également pour objet un système d'identification d'un signal, notamment un signal radar, comprenant des moyens adaptés pour mettre en oeuvre le procédé de détermination des caractéristiques dudit signal en présence d'interférences selon l'invention. 30 D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent : La figure 1, un diagramme temporel représentant un exemple de signal radar caractérisable par l'invention, - La figure 2, un diagramme temporel illustrant un exemple d'influence d'interférences sur la réception du signal radar décrit à la figure 1, - La figure 3a, un organigramme illustrant les étapes de mise en oeuvre du procédé selon l'invention selon le meilleur mode de réalisation, - La figure 3b, un organigramme illustrant les étapes supplémentaires de mise en oeuvre dans une variante de réalisation de l'invention, La figure 4, un organigramme illustrant les étapes de mise en oeuvre d'une méthode alternative applicable lorsqu'une connaissance a priori des modes du signal radar possibles est disponible, - La figure 5, un synoptique d'un système d'identification d'un radar selon l'invention. La figure 1 représente, sur un diagramme temporel, un exemple de signal radar constitué d'une succession d'impulsions 101-105. Un tel signal est caractérisé d'une part par la longueur d'une impulsion LI, autrement dit sa durée temporelle, et d'autre part par la période de répétition des impulsions PRI, c'est-à-dire la différence de temps entre l'émission ou la réception de deux impulsions successives. Le procédé selon l'invention a pour but de caractériser ces deux paramètres LI,PRI qui peuvent être constants ou variables au cours du temps selon que la forme d'onde ou mode employé par le radar est fixe ou variable au cours du temps.
La figure 2 représente, sur un diagramme temporel similaire à celui de la figure 1, un exemple de signal radar reçu en présence d'interférences. Les rayonnements émis par des dispositifs dans l'environnement proche du détecteur de signal radar provoquent une altération du train d'impulsions principalement sous deux formes. Une impulsion émise 201 peut être totalement altérée par l'interférence et être rendue invisible pour le détecteur. C'est le cas lorsque la durée d'une plage d'interférence 210 coïncide totalement avec l'arrivée de cette impulsion 201. Une impulsion émise 202,203 peut également être altérée partiellement ce qui engendre une modification de sa longueur d'impulsion LI et indirectement de la période de répétition d'impulsions mesurée entre l'impulsion partiellement altérée 202,203 et l'impulsion suivante. On distingue deux cas de troncature d'impulsion. Le premier cas correspond à une plage d'interférence 220 qui coïncide avec le début d'une impulsion 202. On dit alors que le front avant de l'impulsion 202 est tronqué. Le second cas correspond à une plage d'interférence 230 qui coïncide avec la fin d'une impulsion 203. On dit que cette impulsion 203 est tronquée sur son front arrière. Dans la suite de la description, l'hypothèse est prise que le dispositif selon l'invention reçoit au moins une impulsion dont le front avant ou le front arrière n'est pas tronqué. Lorsque l'invention est exécutée par un dispositif d'interception de signal radar embarqué dans une plateforme multi-systèmes, alors une information sur la présence ou l'absence d'interférences est disponible en permanence pour indiquer au dispositif d'interception qu'une impulsion de signal reçue est valide ou au contraire est une impulsion tronquée front avant ou front arrière. L'information sur la présence d'interférences peut être matérialisée sous la forme d'un signal binaire synchrone du signal radar reçu, tel que représenté sur le bas de la figure 2. Lorsque ce signal est à 0, il indique l'absence d'interférences et permet de valider qu'une impulsion reçue à ce moment est une impulsion non tronquée. Lorsque ce signal 210,220,230 présente une amplitude différente de zéro, cela indique la présence d'interférences. Si une impulsion 202 est reçue immédiatement après la fin d'une période 220 d'interférences, on en déduit que le front avant de l'impulsion a probablement été tronqué. Si une impulsion 203 est reçue immédiatement avant le début d'une période 230 d'interférences, on en déduit que le front arrière de l'impulsion a probablement été tronqué. L'indicateur d'interférences peut être fourni par la plateforme elle-même lorsqu'il consiste à indiquer l'émission d'au moins un signal gênant par un dispositif de la plateforme. L'indicateur d'interférences peut également être obtenu par le biais d'une analyse du spectre électromagnétique dans la bande de fréquences visée par le détecteur de signal radar ou plus généralement par toute méthode apte à fournir une indication sur la présence, et éventuellement la puissance, de signaux interférents. L'indicateur d'interférences permet de calculer un taux moyen d'interférences Ti sur la durée d'acquisition égal à la durée totale occupée 15 par au moins un signal interférent divisée par la durée totale de l'acquisition. La figure 3 représente un organigramme du procédé de caractérisation d'un signal radar selon l'invention. Dans une première étape 301, le signal radar à caractériser est 20 intercepté et les impulsions de signal sont recueillies. A l'aide de l'indicateur d'interférences, il est possible de distinguer les impulsions reçues valides d'une part et les impulsions reçues tronquées ou ayant une forte probabilité d'être tronquées d'autre part. De plus il est possible de distinguer, parmi les impulsions reçues tronquées, celles dont le front avant ou le front arrière 25 n'est pas altéré. Dans une variante de réalisation de l'invention, la première étape 301 peut être omise si la caractérisation du signal radar ne se fait pas instantanément. Dans un tel cas, les impulsions de signal sont recueillies antérieurement à la mise en oeuvre du procédé et sont assorties d'un 30 indicateur d'interférences. Les impulsions mesurées et l'indicateur d'interférences sont alors fournis comme point d'entrée du procédé selon l'invention. Dans une étape suivante 302, la période moyenne de répétition des impulsions PRIm est calculée en prenant en compte les mesures brutes des 5 impulsions reçues ainsi que le taux de présence d'interférences, c'est-à-dire la proportion de temps disponible pendant lequel aucune interférence ne vient gêner la capture du signal radar. La période moyenne PRIm est égale au temps d'éclairement TE que divise le nombre N d'impulsions émises pendant ce temps. Le temps d'éclairement TE correspond à la durée 10 d'écoute du signal radar. Le nombre N d'impulsions émises pendant cette durée est égale à N = n/(1-Ti) avec n le nombre d'impulsions dont le front avant est non tronqué et qui sont reçues pendant la durée TE et Ti le taux de présence d'interférences pendant la durée TE. Ainsi la période moyenne de répétition d'impulsions peut être obtenue à l'aide de la relation : PRIm = (1- 15 Ti).TE/n. Alternativement, il est possible de remplacer n par le nombre d'impulsions dont le front arrière est non tronqué ou encore par une combinaison, par exemple une moyenne, du nombre d'impulsions reçues dont le front avant est non tronqué d'une part et du nombre d'impulsions reçues dont le front arrière est non tronqué d'autre. Une telle combinaison 20 peut s'exprimer sous la forme n'= (a.nav+13.nar)/(a-FE3), avec nav le nombre d'impulsions non tronquées front avant reçues, nar le nombre d'impulsions non tronquées front arrière reçues et a,I3 deux coefficients de pondération positifs ou nuls. 25 Dans une étape suivante 303, l'écart type APRI, de la période moyenne PRIm est calculé afin d'évaluer la précision de la moyenne précédemment estimée. Ce calcul peut être effectué par le biais de la relation suivante : APRIm = (Ti/n)1/2.PRIm. Dans le cas où n désigne une combinaison du nombre d'impulsions 30 reçues dont le front avant est non tronqué d'une part et du nombre d'impulsions reçues dont le front arrière est non tronqué d'autre part, le calcul d'écart type précité peut être pondéré d'un facteur qui dépend des p coefficients a,r3 : APRI',= 2 n2 F afin de tenir compte de l'indépendance des évènements respectifs associés à la réception d'un front avant d'une impulsion ou d'un front arrière lorsque la longueur d'une impulsion est grande 5 par rapport à la période d'apparition des interférences. Dans une étape suivante 304, on détermine un domaine d'agilité [PRImin,PRImax] de la période de répétition des impulsions, autrement dit un domaine dans lequel il est probable que la valeur réelle de la période de 10 répétition se situe. La période de répétition des impulsions du signal reçu varie statistiquement dans l'intervalle suivant PRIm - APRIm < PRI < PRIm + APRIm. Si la période de répétition est fixe pendant la durée d'acquisition TE, alors le domaine d'agilité peut être défini par : 15 PRImin= PRIm - APRIm PRImax= PRIm + APRIm Lorsque la période de répétition varie pendant la durée d'acquisition TE, il est nécessaire de définir un domaine d'agilité qui permette de suivre ces variations tout en filtrant les valeurs qui ne sont pas pertinentes. 20 Si la période de répétition des impulsions varie dans un rapport égal à un facteur a pour un même émetteur radar, alors on peut définir une seconde plage de variation ou domaine d'agilité de la période de répétition, autour de sa valeur moyenne, égale à [PRIm/./a; -Va.PRI,]. Si la période de répétition est estimée à sa valeur basse PRImin alors le domaine d'agilité associé est 25 égal à [(PRIm - APRIm)hia ; -Nia.(PRIm - APRIm)]. Si au contraire, la période de répétition est estimée à sa valeur haute PRImax alors le domaine d'agilité associé est égal à [(PRIm + APRI')hla ; + APRIm)]. Au final, le domaine d'agilité retenu est l'intersection des domaines d'agilité associés respectivement à la valeur basse et à la valeur haute de la 30 période de répétition moyenne, soit : [(PRIm + APRIm)/'ici ; -Va.(PRIn, - APRIm)]. En choisissant le domaine d'agilité de la sorte, la variation de la période de répétition est limitée à un facteur a, il est ainsi possible d'éliminer les valeurs non pertinentes tout en permettant un suivi des variations éventuelles de la valeur réelle de la période de répétition.
Avantageusement, le facteur a est pris au plus égal à 2 afin d'éliminer, du domaine d'agilité, les multiples de la valeur réelle de la période de répétition des impulsions. Dans une variante de réalisation de l'invention, il est également possible de définir le domaine d'agilité de la façon suivante : [(PRIm + k.APRIm)ka ; - k.APRIm)], avec k un entier strictement positif. Dans une étape suivante 305, on mesure, pour l'ensemble des impulsions non tronquées front avant ou front arrière reçues, les différences de temps d'arrivée DTOA entre deux impulsions successivement reçues. Les différences de temps d'arrivée peuvent être mesurées indifféremment entre deux fronts avant successifs ou deux fronts arrière successifs. Les résultats sont stockés, par exemple sous forme d'un histogramme qui contient alors d'une part les mesures de la période de répétition réelle des impulsions, c'est-à-dire la différence de temps d'arrivée entre deux impulsions successivement émises, mais aussi les mesures des multiples de la période de répétition réelle qui sont obtenues lorsqu'une ou plusieurs impulsions successives sont masquées à cause des interférences. Dans une étape suivante 306, le domaine d'agilité est appliqué à l'histogramme des différences de temps d'arrivée afin de filtrer les mesures 25 hors du domaine. De cette façon, on élimine les mesures des multiples de la période de répétition réelle. A partir des mesures DTOA filtrées, on en déduit la période de répétition réelle PRI du mode de propagation du signal radar capturé. Cette période peut être unique ou peut présenter plusieurs valeurs si le mode de 30 propagation a varié au cours de la capture du signal.
Dans une variante de réalisation de l'étape 305 du procédé selon l'invention, il est possible d'intégrer, dans l'histogramme des différences de temps d'arrivée DTOA, également les mesures effectuées sur des impulsions tronquées front avant ou front arrière. Pour de telles impulsions, la mesure 5 de l'instant de réception d'un front avant ou arrière de l'impulsion est incertaine car une partie de l'impulsion est tronquée du fait des interférences. Dans ce cas, la mesure DTOA réalisée entre deux impulsions dont au moins l'une est tronquée front avant ou front arrière, est assortie d'une fenêtre d'incertitude égale à plus ou moins la longueur d'impulsion moyenne 10 estimée. De cette façon on aboutit à l'issue du procédé à une estimée de la ou des périodes de répétition d'impulsion assorties d'une fenêtre d'incertitude. Les fenêtres d'incertitude peuvent être comparées à une bibliothèque contenant l'ensemble des valeurs possibles que peut prendre la période de répétition d'impulsion d'un signal radar afin d'identifier la plus 15 plausible. Dans une autre variante de réalisation du procédé selon l'invention, lorsque l'étape de filtrage 306 aboutit à ce qu'aucune mesure DTOA n'est conservée car l'ensemble des mesures sont hors du domaine d'agilité, alors 20 on applique la méthode suivante illustrée par l'organigramme de la figure 3b. Lorsqu'aucune mesure DTOA n'est présente dans le domaine d'agilité, cela signifie que les mesures effectuées correspondent toutes à des multiples de la période de répétition d'impulsion réelle. Ce cas de figure peut se produire notamment lorsque le taux d'interférences est important. 25 Dans un tel cas, on cherche alors à identifier une valeur de période de répétition qui correspond à un sous-multiple des mesures DTOA et qui soit situé dans la fenêtre du domaine d'agilité. La recherche peut se faire sur deux mesures DTOA successives car on prend l'hypothèse que la période de répétition d'impulsion reste fixe pendant une durée suffisante pour que la 30 probabilité d'un changement de valeur de la période entre deux mesures successives soit faible. Typiquement, la période de répétition reste fixe sur un horizon d'une dizaine à une centaine d'impulsions. C'est notamment le cas pour des radars fonctionnant en mode cohérent. La recherche d'un sous-multiple de deux mesures DTOA successives appartenant au domaine d'agilité peut se faire par une méthode dite du 5 PGCD, pour Plus Grand Diviseur Commun, dont l'organigramme est représenté à la figure 3b. Dans une étape préalable on définit deux variables X et Y qui sont initialisées respectivement aux valeurs de deux mesures successives de différence de temps d'arrivée X=PPRIi et Y=PPRI2. Comme indiqué ci- 10 dessus, ces mesures peuvent être réalisées indifféremment entre deux fronts avants non tronqués successivement reçus ou deux fronts arrières non tronqués successivement reçus. On calcule 360 l'opération modulo des deux variables X et Y, autrement dit on calcule l'opération Mod(X,Y)= X - Ent[XM.Y, où EntO 15 désigne la partie entière d'un nombre réel, c'est-à-dire l'entier qui lui est immédiatement inférieur ou égal. On compare 361 le résultat de l'opération Mod(X,Y) avec la fenêtre du domaine d'agilité [PRImin,PRImax]. Si le résultat est situé dans cette fenêtre alors il s'agit de la période de répétition d'impulsion PRI= Mod(X,Y). 20 Dans le cas contraire, on réaffecte les variables X et Y en posant X= PPRI2 et Y = Mod(PPRI1,PPRI2) et, si ces deux variables ne sont pas égales en valeur on revient à l'étape 360 de calcul du modulo des variables X et Y. On réitère les étapes 360,361,362 jusqu'à obtenir un résultat de l'opération modulo qui soit compris dans la fenêtre du domaine d'agilité. 25 Dans le cas où, à l'issue de l'étape 362, les deux variables X et Y sont égales en valeur, alors cela signifie qu'il n'a pas été possible de trouver un sous-multiple situé dans le domaine d'agilité. Dans ce cas, la variable X=Y correspond à un multiple de la période de répétition d'impulsion recherchée. Dans une étape 363 supplémentaire, on calcule alors tous les sous-30 multiples de la variable X.
On teste 364 si un seul de ces sous-multiples est situé dans un domaine d'agilité prédéterminé auquel cas il s'agit de la période de répétition d'impulsions PRI recherchée. Dans le cas contraire, il existe plusieurs sous-multiples situés dans le domaine d'agilité et il faut alors lever l'ambigüité restante en exploitant une troisième mesure PPRI3 de différence de temps d'arrivée. La période de répétition d'impulsions PRI recherchée est alors obtenue en calculant 365 l'opération modulo entre la variable X et cette troisième mesure PPRI3. Le domaine d'agilité prédéterminé D' peut être égal au domaine 10 d'agilité D1= [(PRIm + APRIm)hla ; .Va.(PRIm - APRIm)] précédemment introduit à l'étape 304. Dans une variante, si le rapport P= APRIm/PRIm est trop élevé, la fenêtre de filtrage définie à l'étape 304 risque d'être trop étroite et donc de rejeter des sous-multiples ambigus de la variable X. Dans ce cas, et pour 15 éliminer un tel risque, on retient comme domaine d'agilité, non pas l'intersection des plages [(PRIm - APRIm)/'/a; .Nlat.(PRIm - APRIm)] et [(PRIm + APRIm)hla ; -Va.(PRIm + APRIm)] comme expliqué ci-dessus pour l'étape 304 mais l'union de ces deux plages : [(PRIm - APRIm)/Va; -Va.(PRIm - APRIm)] l) [(PRIm + APRIm)/'la ; -Va.(PRIm + APRIm)] ce qui conduit à choisir un second 20 domaine d'agilité D2=[(PRIm - APRIm)/-Ja ; -\la.(PRIm + APRIm)], Dans une autre variante, afin de diminuer encore d'avantage le risque d'aboutir à une valeur de période de répétition d'impulsions erronée, l'étape 364 de test est modifiée de la façon suivante. Si un seul sous-multiple de la 25 variable X est situé dans le premier domaine d'agilité D1= [(PRIm + APRIm)/Va; -Va.(PRIm - APRIm)] et qu'aucun autre sous-multiple de la variable X n'est situé dans le second domaine d'agilité D2=[(PRI, - APRIm)/Va ; -\la.(PRIm + APRIm)], alors on conserve ce sous-multiple comme étant la période de répétition d'impulsions PRI recherchée. 30 Dans le cas contraire, on retient uniquement le ou les sous-multiples compris dans le second domaine d'agilité D2 et on applique, au besoin, la dernière étape 365 pour ne retenir in fine qu'une seule valeur parmi celles situées dans le second domaine d'agilité D2. En outre, si un seul desdits sous-multiples est situé dans le second domaine d'agilité D2 alors l'estimée de la période de répétition des impulsions est égale audit seul sous-multiple.
Enfin, si plusieurs desdits sous-multiples sont situés dans le second domaine d'agilité D2, on applique également la dernière étape 365. Cette variante de réalisation de l'étape de test 364 permet d'écarter le risque que la période de répétition d'impulsions réelle soit en dehors du domaine d'agilité D1 et qu'un sous-multiple de la solution issue de la méthode dite du PGCD fournisse à la fois un sous-multiple dans le domaine d'agilité D1(ce sous-multiple serait alors une fausse valeur) et à la fois un sous-multiple dans le domaine d'agilité D2 (ce sous-multiple serait alors la vraie valeur).
Dans une autre variante de réalisation, le procédé selon l'invention consiste en outre à caractériser la longueur des impulsions LI émises par le radar. La longueur d'impulsion LI peut être constante ou variable lorsque le mode d'émission du radar change au cours du temps. Si des impulsions non tronquées sont reçues, alors on mesure leurs 20 longueurs d'impulsions LI qui vont permettre de caractériser directement le signal radar. Si aucune impulsion non tronquée n'est disponible, alors, dans une nouvelle étape 310, on détermine une valeur moyenne de la longueur d'impulsion Llm. 25 Cette valeur moyenne Llm peut être estimée en calculant la somme des longueurs des impulsions reçues et en divisant le résultat par le nombre n d'impulsions non tronquées front avant ou front arrière reçues c'est-à-dire les impulsions dont le front avant ou arrière est reçu en l'absence d'interférences mais dont la longueur est tronquée du fait de l'apparition 30 d'interférences pendant la réception de l'impulsion. Dans une variante de réalisation, comme déjà explicité ci-dessus, le nombre n peut être remplacé par toute combinaison du nombre d'impulsions non tronquées front avant et du nombre d'impulsions non tronquées front arrière. Cette estimation découle du fait que la variable définie par la somme des longueurs des impulsions reçues suit une statistique identique à celle de la variable définie par le nombre de fronts avant ou de nombre de fronts arrière d'une impulsion reçue, c'est-à-dire une loi binomiale. En moyenne, la somme des longueurs des impulsions reçues est donc donnée par la relation suivante : ELI = (1-Ti).N.Llm = n. Llm On en déduit directement une estimée de la valeur moyenne Llm des longueurs d'impulsions. ZH(j) (i) LIMOY - - (1- Ti).N n OU (1) (a+ ME LM Limoy = anA, +tJflAr Dans une étape supplémentaire 311, il est également possible de 15 déterminer l'écart type ALIm associé à la valeur moyenne de la longueur d'impulsion. Comme la somme des longueurs des impulsions reçues suit une statistique de loi binomiale, on peut démontrer que sa variance est donnée par la relation suivante : 2 20 Var E LI(J)1= EL/(J)1 _N(1 -T,)LI} =[N(1-T,)T,ILI2 (2) _ A partir de la relation (2), on en déduit la variance de la longueur d'impulsion moyenne : 2 E LI(j) - LI N(1-7;) 1LI2 N(1-T,) Var (3) La relation (3) peut encore s'écrire sous la forme : Varkl moyi= EspiLl moy- LI12}- No _Ti) LI2 , où Esp{} désigne l'espérance mathématique En remplaçant la longueur d'impulsion LI par sa valeur moyenne 5 estimée, on obtient la relation (4) : ELI(j)- 2 i Var[LIAloyl= Esp{[LImoy - LI12} = N(1 - (4) L'écart type de la valeur moyenne peut alors être obtenu à partir de la relation (5) suivante : z11,1,7,= .ear[Llmoy]- Euw _r,[,LI(,)] (5) 10 avec LI(j) une mesure de la longueur d'impulsion reçue d'indice j. Le calcul développé ci-dessus présente une méthode permettant de déterminer l'écart type de la valeur moyenne de la longueur d'impulsion. Des variantes de calcul existent, en particulier la relation (5') suivante donne une 15 estimation légèrement différente de celle obtenue via la relation (5) : E Lm] (5') _ J 4L1,7, =11124 n' Le facteur multiplicatif 2 a pour objectif de majorer l'erreur quadratique d'estimation de la longueur d'impulsion moyenne dans le cas où les deux 20 sources d'erreurs, à savoir le nombre d'impulsions d'une part et la durée du signal reçu d'autre part, sont décorrélées. C'est le cas notamment lorsque les longueurs des impulsions sont importantes par rapport à la période d'apparition des interférences.
Le calcul de l'écart type de la valeur moyenne de la longueur d'impulsion développé ci-dessus est utilisé pour assortir la mesure de la valeur moyenne d'une précision. Cette précision permet, par exemple, de définir un intervalle d'incertitude autour de la valeur moyenne, et d'utiliser cet intervalle pour le comparer à des valeurs présumées de la longueur d'impulsion réelle. Un exemple d'intervalle de précision consiste à définir l'intervalle [Llm - ALIm ; Llm + ALIm J. Le calcul de l'écart type de la valeur moyenne de la longueur d'impulsion sert également à décider si la précision de la valeur moyenne est suffisante et si ce n'est pas le cas, à envisager une méthode d'estimation alternative. A cet effet, l'écart type peut être remplacé par la variance ou l'écart type normalisé ou tout autre indicateur statistique de la précision de la moyenne.
Le procédé selon l'invention tel que décrit ci-dessus à l'appui de la figure 3, s'applique préférentiellement lorsque la précision des valeurs moyennes des paramètres de la forme d'onde du signal radar est suffisante. La précision de la période de répétition moyenne et de la longueur d'impulsion moyenne est inversement proportionnelle à l'écart type calculé sur ces deux moyennes. Dans le cas de la période de répétition des impulsions, l'écart type sur la valeur moyenne APRIm détermine la largeur du domaine d'agilité [PRImin,PRImax]. Si on pose P= APRIm/PRIm, l'écart type normalisé par la moyenne, le rapport entre la borne supérieure et la borne inférieure du domaine d'agilité est alors égal à PRImax/PRImin = a[(1-P)/(1+P)]. Plus l'écart type diminue, plus la précision sur la mesure de la moyenne augmente et plus la largeur du domaine d'agilité est augmentée. Au contraire, plus l'écart type augmente, plus la précision sur la mesure de la moyenne décroit et plus le domaine d'agilité est restreint.
Lorsque le mode d'émission du signal radar évolue au cours du temps et que l'on souhaite suivre des variations de la période de répétition, alors il n'est pas souhaitable de restreindre trop le domaine d'agilité. Il est alors possible de fixer une valeur seuil Pmax de l'écart type normalisé P au delà duquel la précision de la valeur moyenne devient peu fiable et le domaine d'agilité devient trop étroit. Dans ce cas le procédé décrit à la figure 3 n'est plus 5 performant. Ainsi, dans une variante de l'étape 303 de calcul de l'écart type de la valeur moyenne de la période de répétition des impulsions, il est possible d'ajouter une étape de calcul du ratio P et de comparaison de ce ratio à la valeur seuil Pmax. Si P> Pmax alors la suite des étapes 304-306 n'est pas appliquée et peut être remplacée par une méthode alternative, un 10 exemple d'une telle méthode est donné plus loin dans la description. La valeur du seuil Pmax est déterminée de sorte à fixer l'excursion du domaine d'agilité qui correspond à la plage de variation de la période de répétition d'impulsion que l'on souhaite suivre. Autrement dit la valeur du seuil Pmax est calculée en fonction de l'excursion minimale souhaitée du 15 domaine d'agilité. Par exemple, le rapport PRImax/PRImin peut être choisi égal à 1,33, ce qui équivaut, lorsque le facteur correctif a est égal à 2, à un seuil Pmax égal à 20% de la valeur moyenne. Pour déterminer si la précision de la valeur moyenne de la période de 20 répétition des impulsions ou de la longueur moyenne des impulsions est suffisante, le test de comparaison P> Pmax décrit ci-dessus peut être appliqué indifféremment à l'écart type, à l'écart type normalisé à la variance ou à tout autre indicateur équivalent de la précision de la valeur moyenne. 25 On décrit à présent une méthode alternative à celle décrite ci-dessus qui peut avantageusement s'appliquer lorsque la précision sur la valeur moyenne de la période de répétition des impulsions est mauvaise, autrement dit lorsque l'écart type sur cette moyenne, ou l'écart type normalisé, dépasse un seuil maximum. De façon similaire, cette méthode s'applique également 30 lorsque l'écart type associé à la moyenne de la longueur d'impulsion est également trop élevé rendant la mesure de moyenne peu fiable.
Préférentiellement on estime que la moyenne de la longueur d'impulsion est peu précise lorsque l'écart type normalisé par la moyenne est supérieur à 30%.
On prend comme hypothèse que le radar que l'on cherche à caractériser a déjà pu être analysé précédemment ou en tout cas qu'une connaissance a priori des différents modes d'émission possibles que peut mettre en oeuvre ce radar est disponible. Les paramètres de période de répétition d'impulsions et de longueur 10 d'impulsions de chaque mode sont stockés dans une base de données 400 ou une mémoire que comprend le système selon l'invention. Pour chacun des modes disponibles en mémoire, les étapes suivantes sont effectuées. Dans une première étape 401, on calcule le nombre N d'impulsions 15 supposées reçues si le mode supposé est activé par le radar. Ce nombre N est égal à la durée d'éclairement TE que divise la période de répétition des impulsions du mode supposé, période qui est disponible en mémoire du système selon l'invention. Dans une deuxième étape 402, on détermine la plage de variations 20 dans laquelle le nombre d'impulsions reçues n, dont le front avant est non tronqué, est censé varier, compte tenu du taux d'interférences Tl. Cette plage est centrée sur la valeur moyenne (1-Ti)N du nombre d'impulsions reçues avec une variation plus ou moins égale à l'écart type de cette moyenne [(1-Ti).Ti.N]1"2. 25 Ainsi les bornes inférieure et supérieure de cette plage sont respectivement égales à: nmin= (1-Ti)N - [(1-Ti).Ti.Nr2 et nmax= (1-Ti)N + [(1-Ti).Ti.N]112 Dans une variante de réalisation de la deuxième étape 403, la plage de variations peut être choisie centrée sur la valeur moyenne du nombre 30 d'impulsions reçues mais avec une variation plus ou moins égale à un multiple k de l'écart type de cette moyenne, ce qui donne les bornes inférieure et supérieure de la plage choisie respectivement égales à: nmin= (1-Ti)N - k[(1-Ti).Ti.N]1/2 et nmax= (1-Ti)N + k[(1-TO.Ti.N]112 k est un nombre positif différent de 0.
Dans une troisième étape 403, on compare le nombre d'impulsions n reçues avec un front avant non tronqué à la plage de variations déterminée à l'étape 402 précédente. Les étapes 401,402,403 sont réitérées pour chaque mode disponible 10 en mémoire et le mode le plus probable est retenu, c'est-à-dire celui qui vérifie le critère nmm n < nmax. Dans une variante de réalisation de l'invention, le nombre n d'impulsions reçues avec un front avant non tronqué peut être remplacé par 15 le nombre d'impulsions reçues avec un front arrière non tronqué. Dans une autre variante de réalisation de l'invention, il est également possible de remplacer le nombre n d'impulsions reçues avec l'un des deux front tronqués, par une combinaison, par exemple une moyenne, du nombre d'impulsions reçues dont le front avant est non tronqué d'une part et du 20 nombre d'impulsions reçues dont le front arrière est non tronqué d'autre part. Autrement dit, le nombre n d'impulsions reçues avec l'un des deux front tronqués peut être remplacé par le nombre n'= (a.nav+13.nar)/(a+13), avec nav le nombre d'impulsions non tronquées front avant reçues, na, le nombre d'impulsions non tronquées front arrière reçues et a,13 deux coefficients de 25 pondération positifs ou nuls. L'utilisation d'une combinaison du nombre d'impulsions dont le front avant est non tronqué et du nombre d'impulsions dont le front arrière est non tronqué permet de renforcer les contrôles de la validité de la période de répétition d'impulsions dans le cas où la longueur d'impulsion est grande par 30 rapport à la période de répétition d'impulsion. En effet, dans un tel cas, les évènements associés à l'apparition d'un front avant ou d'un front arrière peuvent être considérés indépendants, les sources d'erreurs impactant les deux mesures sont donc décorrélées. Lorsque le contrôle a été effectué pour deux valeurs quelconques du couple a,f3 il est vérifié pour toutes les autres valeurs du couple.
Si alternativement ou additionnellement, les longueurs d'impulsions de chaque mode supposé sont connues, on applique, alternativement ou en supplément, les étapes suivantes. Dans une première étape 404, on détermine la plage de variations dans laquelle la somme des longueurs des impulsions reçues est censée varier, compte tenu du taux d'interférences Tl. Cette plage est centrée sur la valeur moyenne (1-Ti)N.LI, où LI est la longueur d'impulsion supposée du mode radar, avec une variation plus ou moins égale à l'écart type de cette moyenne [(1-Ti).Ti.N]112.LI.
Ainsi les bornes inférieure et supérieure de cette plage sont respectivement égales à: Smin= ({1-Ti)N - [(1-TO.Ti.N]1/2}11 et Smax= {(1-Ti)N + [(1-TO.Ti.N11/2}11. Dans une variante de réalisation de l'invention les bornes inférieure et supérieure de la plage peuvent également être choisies de la façon suivante, 20 avec k un nombre positif différent de 0. Smin= ({1-Ti)N -k [(1-Ti).Ti.N]1/2111 et S.= {(1-Ti)N + k[(1-Ti).Ti.N]1/2111. Dans une étape suivante 405, on compare la somme des longueurs des impulsions reçues, à la plage de variations déterminée à l'étape 404 25 précédente. Les étapes 401,404,405 sont réitérées pour chaque mode disponible en mémoire et le mode le plus probable est retenu, c'est-à-dire celui qui vérifie le critère suivant : la somme des longueurs des impulsions reçues est comprise dans la plage de variations [Smin, Smax].
L'estimation du nombre d'impulsions reçues et de la somme des longueurs d'impulsions reçues peut se faire indépendamment ou conjointement afin d'améliorer la fiabilité des résultats.
La figure 5 représente schématiquement, sur un synoptique, un exemple de système d'identification de radar 500 selon l'invention. Un tel système comporte des moyens de réception d'un signal radar constitués d'une antenne 501 et de moyens 502 de recueil des impulsions incidentes dudit signal.
Un tel système 500 comporte en outre un moyen 503 de regroupement des impulsions appartenant à un même émetteur. Un tel moyen 503 est en charge de séparer les signaux provenant d'émetteurs différents. Le système 500 selon l'invention comporte en outre un moyen 504 pour exécuter le procédé, selon l'invention, de détermination des caractéristiques de période de répétition des impulsions et/ou de longueur des impulsions et pour fournir ces caractéristiques à un moyen 505 d'identification de la forme d'onde mise en oeuvre par le signal radar, par exemple en comparant les caractéristiques estimées à une base de données contenant l'ensemble des caractéristiques possibles ou présumées. Un moyen 507 interne ou externe au système selon l'invention est utilisé pour fournir un indicateur d'interférences au moyen 504 pour exécuter le procédé selon l'invention. Le système 500 selon l'invention peut également comporter un moyen 25 506 de rafraîchissement de la piste dédiée au suivi de l'émetteur. Les moyens 503, 504, 505, 506 peuvent être implémentés à partir d'éléments logiciel. 30

Claims (21)

  1. REVENDICATIONS1. Procédé de détermination, en présence d'interférences, des caractéristiques d'un signal comprenant une pluralité d'impulsions émises (101-105), ledit procédé consistant à déterminer au moins la période de répétition des impulsions (PRI), ledit procédé étant caractérisé en ce qu'il comporte les étapes suivantes, exécutées pour une portion dudit signal de durée donnée TE: Calculer (302) la période de répétition moyenne PRIm des impulsions en présence d'interférences, Calculer (303) l'écart type APRIm de la période de répétition moyenne PRIm des impulsions en présence d'interférences, Définir (304) un domaine d'agilité [PRImin ; PRImax] de la période de répétition des impulsions (PRI) de la façon suivante PRImin = (PRIm + kAPRIm)hla et PRImax = - k.APRIm), avec a un facteur d'excursion strictement positif et inférieur ou égal à 2 et k un nombre strictement positif, Calculer (305), pour l'ensemble des impulsions comprises dans ladite portion de signal et dont le front avant ou le front arrière est non tronqué, les différences de temps d'arrivée (DTOA) entre deux fronts avants ou deux fronts arrières de deux impulsions successivement reçues, Eliminer (306) les valeurs des différences de temps d'arrivée qui sont hors du domaine d'agilité [PRImin ; PRImax], la ou les valeurs restantes donnant une estimée de la période de répétition des impulsions (PRI).
  2. 2. Procédé de détermination, en présence d'interférences, des caractéristiques d'un signal comprenant une pluralité d'impulsions émises (101-105), ledit procédé consistant à déterminer au moins la période derépétition des impulsions (PRI), ledit procédé étant caractérisé en ce qu'il comporte les étapes suivantes, exécutées pour une portion dudit signal de durée donnée TE: Calculer (302) une estimée de la période de répétition moyenne PRI, des impulsions en présence d'interférences, Calculer (303) une estimée de l'écart type APRI, de la période de répétition moyenne PRI, des impulsions en présence d'interférences, - Si la variance APRI'2 ou l'écart type APRI, ou l'écart type normalisé par la moyenne APRI'/PRI, est inférieur à un seuil de précision Pmax prédéterminé, alors exécuter les étapes (304,305,306) suivantes : - Définir (304) un domaine d'agilité [PRImin ; PRI.] de la période de répétition des impulsions (PRI), de la façon suivante PRImin = (PRI' + k.APRI,)/-Va et PRI. = - k.APRI'), avec a un facteur d'excursion strictement positif et inférieur ou égal à 2 et k un nombre strictement positif, Calculer (305), pour l'ensemble des impulsions comprises dans ladite portion de signal et dont le front avant ou le front arrière est non tronqué, les différences de temps d'arrivée (DTOA) entre deux fronts avants ou deux fronts arrières de deux impulsions successivement reçues, - Eliminer (306) les valeurs des différences de temps d'arrivée qui sont hors du domaine d'agilité [PRImin ; PRImax], la ou les valeurs restantes donnant une estimée de la période de répétition des impulsions (PRI).
  3. 3. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 2 dans lequel ledit seuil de précision Pmax est déterminé en fonction de l'excursion PRImax/PRImin minimale souhaitée du domaine d'agilité.
  4. 4. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes dans lequel la période de répétition moyenne PRIm des impulsions est calculée (302) à l'aide de la relation suivante : PRIm = (1-Ti).TE/n, avec Ti un indicateur du taux de présence d'interférences pendant la durée TE et n le nombre d'impulsions, comprises dans ladite portion, dont le front avant est non tronqué ou le nombre d'impulsions, comprises dans ladite portion, dont le front arrière est non tronqué ou une combinaison du nombre nav d'impulsions dont le front avant est non tronqué et du nombre nar d'impulsions dont le front arrière est non tronqué sous la forme n=(a.nav+6.nar)/(a+6), où a et sont deux entiers positifs ou nuls.
  5. 5. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 4 dans lequel l'écart type APRIm de la période de répétition moyenne PRIm des impulsions est calculé (303) à l'aide de la relation suivante : APRIm = (Tiin)1/2.pRim,
  6. 6. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 5 dans lequel l'écart type APRIm de la période de répétition moyenne PRIm des impulsions est pondéré par le 11,, coefficient : 2 a2lorsque n est une combinaison du nombre a+/3 d'impulsions dont le front avant est non tronqué et du nombre d'impulsions dont le front arrière est non tronqué sous la forme n=(a.nav+6.nar)/(a+6), où a et 6 sont deux entiers positifs ou nuls.
  7. 7. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes dans lequel le nombre k est pris égal à 1. nav nar30
  8. 8. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes dans lequel ledit facteur correctif a est pris égal à 2.
  9. 9. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes dans lequel les différences de temps d'arrivée (DTOA) sont en outre calculées (305) pour les impulsions dont le front avant ou le front arrière est tronqué et sont dans ce cas assorties d'une fenêtre d'incertitude de plus ou moins la longueur d'impulsion moyenne Llm.
  10. 10. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications 1 à 8 dans lequel, lorsqu'aucune valeur n'est conservée après l'étape d'élimination (306) des valeurs des différences de temps d'arrivée hors domaine d'agilité [PRImin ; PRImax], une estimée de la période de répétition des impulsions (PRI) est obtenue en recherchant un sous-multiple, situé dans ledit domaine d'agilité [PRI,,in ; PRI.], d'au moins deux mesures successives de différences de temps d'arrivée.
  11. 11. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 10 dans lequel ladite recherche d'un sous-multiple est effectuée en exécutant au moins les étapes suivantes : Définir deux variables X,Y et les initialiser à deux mesures successives PPRI1,PPRI2 de différences de temps d'arrivée, Calculer (360) l'opération modulo des deux variables Mod(X,Y), Si (361) le résultat de l'opération modulo est situé dans ledit domaine d'agilité D1=[(PRIm + APRIm)ka ; - APRIm)], l'estimée de la période de répétition des impulsions (PRI) est égale à ce résultat,Si (361) le résultat de l'opération modulo est situé hors du domaine d'agilité, réaffecter (362) la variable X à la seconde mesure PPRI2 et la variable Y au résultat de l'opération modulo Mod(PPRIi, PPRI2) entre la première mesure PPRIi et la seconde mesure PPRI2, Si la valeur de X est différente de la valeur de Y, alors réitérer les étapes (360,361,362) précédentes, Si la valeur de X est égale à la valeur de Y, alors calculer (363) l'ensemble des sous-multiples de la variable X, Exécuter l'étape de test (364) suivante : L Si (364) un seul desdits sous-multiples est situé dans ledit domaine d'agilité D1, l'estimée de la période de répétition des impulsions (PRI) est égale audit seul sous-multiple, ii. Si (364) plusieurs desdits sous-multiples sont situés dans ledit domaine d'agilité D1, calculer (365) le résultat de l'opération modulo Mod(X,PPRI3) entre la variable X et une troisième mesure successive de différences de temps d'arrivée PPRI3, l'estimée de la période de répétition des impulsions (PRI) étant égale à ce résultat.
  12. 12. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 11 dans lequel ledit domaine d'agilité D1 est remplacé par un second domaine d'agilité D2= [(PRI,, - APRIm)hla ; -\/a.(PRIn, + APRIm)].
  13. 13. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 11 dans lequel l'étape de test (364) est modifiée de la façon suivante: Si (364) un seul desdits sous-multiples est situé dans le premier domaine d'agilité D1=[(PRIm + APRIm)hla ; -Va.(PRIm - APRIm)], et qu'aucun autre desdits sous-mutliples n'est situé dans un seconddomaine d'agilité D2= RPRIm - APRIm)hla ; + APRIm)] alors l'estimée de la période de répétition des impulsions (PRI) est égale audit seul sous-multiple, Si (364) plusieurs desdits sous-multiples sont situés dans ledit domaine d'agilité D1 et qu'aucun autre desdits sous-mutliples n'est situé dans un second domaine d'agilité D2= RPRIm - APRIm)/Va; -Va.(PRIn, + APRIm)], calculer (365) le résultat de l'opération modulo Mod(X,PPRI3) entre la variable X et une troisième mesure successive de différences de temps d'arrivée PPRI3, l'estimée de la période de répétition des impulsions (PRI) étant égale à ce résultat. Si (364) un seul desdits sous-multiples est situé dans le second domaine d'agilité D2= RPRIm - APRIm)/Va; -\la.(PRIm + APRIm)] alors l'estimée de la période de répétition des impulsions (PRI) est égale audit seul sous-multiple, Si (364) plusieurs desdits sous-multiples sont situés dans le second domaine d'agilité D2= RPRIm - APRIm)hla ; -Nla.(PRIm + APRIm)], calculer (365) le résultat de l'opération modulo Mod(X,PPRI3) entre la variable X et une troisième mesure successive de différences de temps d'arrivée PPRI3, l'estimée de la période de répétition des impulsions (PRI) étant égale à ce résultat.
  14. 14. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes consistant à déterminer en outre la longueur d'impulsion (LI), ledit procédé comportant les étapes suivantes : Si ladite portion de signal comporte des impulsions non tronquées, utiliser les longueurs de ces impulsions comme estimée de la longueur d'impulsion (LI), 2 9989 73 34 - Si ladite portion de signal ne comporte que des impulsions tronquées, calculer (310) une estimée de la valeur moyenne Llm de la longueur d'impulsion. 5
  15. 15. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications 1 à 13 consistant à déterminer en outre la longueur d'impulsion (LI), ledit procédé comportant les étapes suivantes : Si ladite portion de signal comporte des impulsions non tronquées, 10 utiliser les longueurs de ces impulsions comme estimée de la longueur d'impulsion (LI), Calculer un indicateur de précision ALIm de la valeur moyenne Llm de la longueur d'impulsion, Si ledit indicateur est inférieur à un seuil de précision 15 prédéterminé, alors exécuter l'étape suivante (310) : Si ladite portion de signal ne comporte que des impulsions tronquées, calculer (310) une estimée de la valeur moyenne Llm de la longueur d'impulsion. 20
  16. 16. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 15 dans lequel ledit indicateur est un écart type ALIm, un écart type normalisé par la moyenne ALImiLlm ou une variance ALI,2. 25
  17. 17. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon la revendication 16 dans lequel ledit indicateur est un écart type ALIm calculé à l'aide de l'une des deux relations suivantes : 41,4, = JTf[LJo] OU LILI', _\1211,3 EL/(i)] avec ELI(j)une mesure de la n j n somme des longueurs des impulsions reçues comprises dans ladite 30 portion et n le nombre d'impulsions, comprises dans ladite portion, dont lefront avant est non tronqué ou le nombre d'impulsions, comprises dans ladite portion, dont le front arrière est non tronqué ou une combinaison du nombre na, d'impulsions dont le front avant est non tronqué et du nombre nar d'impulsions dont le front arrière est non tronqué sous la forme n=(a.na,-F13.nar)/(a+8), où a et 13 sont deux entiers positifs ou nuls.
  18. 18.Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications 14 à 17 dans lequel la valeur moyenne Llm de la longueur d'impulsion est calculée (310) en divisant la somme des longueurs des impulsions comprises dans ladite portion par le nombre d'impulsions, comprises dans ladite portion, dont le front avant est non tronqué ou par le nombre d'impulsions, comprises dans ladite portion, dont le front arrière est non tronqué ou par une combinaison du nombre na, d'impulsions dont le front avant est non tronqué et du nombre nar d'impulsions dont le front arrière est non tronqué sous la forme n=(a.na,-1-8.nar)/(a-1-13), où a et rs sont deux entiers positifs ou nuls.
  19. 19. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes dans lequel un indicateur d'interférences synchrone dudit signal est utilisé pour déterminer si une impulsion est tronquée ou non.
  20. 20. Procédé de détermination des caractéristiques d'un signal en présence d'interférences selon l'une des revendications précédentes dans lequel ledit procédé comporte une étape préalable d'interception et de recueillement (301) de ladite portion de signal pendant une durée donnée TE.
  21. 21.Système (500) d'identification d'un signal, notamment un signal radar, comprenant des moyens (504) adaptés pour mettre en oeuvre le procédéde détermination des caractéristiques dudit signal en présence d'interférences selon l'une des revendications précédentes.
FR1203239A 2012-11-30 2012-11-30 Procede de determination des caracteristiques d un signal radar en presence d interferences Active FR2998973B1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
FR1203239A FR2998973B1 (fr) 2012-11-30 2012-11-30 Procede de determination des caracteristiques d un signal radar en presence d interferences

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR1203239A FR2998973B1 (fr) 2012-11-30 2012-11-30 Procede de determination des caracteristiques d un signal radar en presence d interferences

Publications (2)

Publication Number Publication Date
FR2998973A1 true FR2998973A1 (fr) 2014-06-06
FR2998973B1 FR2998973B1 (fr) 2015-01-23

Family

ID=49448198

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1203239A Active FR2998973B1 (fr) 2012-11-30 2012-11-30 Procede de determination des caracteristiques d un signal radar en presence d interferences

Country Status (1)

Country Link
FR (1) FR2998973B1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2527630A (en) * 2014-11-21 2015-12-30 Dizic Co Ltd Ranging method, ranging device, location device and location method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5323161A (en) * 1986-08-29 1994-06-21 The Marconi Company Limited Method of assessing radar target pulses
US20040201517A1 (en) * 2003-04-08 2004-10-14 Globespanvirata, Inc. Radar detector having a multi-period peridocity validator and method therefor

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5323161A (en) * 1986-08-29 1994-06-21 The Marconi Company Limited Method of assessing radar target pulses
US20040201517A1 (en) * 2003-04-08 2004-10-14 Globespanvirata, Inc. Radar detector having a multi-period peridocity validator and method therefor

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KAZIMIERZ BANASIAK ET AL: "Radar signals monitoring", INTERNATIONAL RADAR SYMPOSIUM. IRS 2006, IEEE, PISCATAWAY, NJ, USA, 1 May 2006 (2006-05-01), pages 1 - 4, XP031138439, ISBN: 978-83-7207-621-2 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2527630A (en) * 2014-11-21 2015-12-30 Dizic Co Ltd Ranging method, ranging device, location device and location method
GB2527630B (en) * 2014-11-21 2017-03-01 Dizic Co Ltd Ranging method, ranging device, location device and location method

Also Published As

Publication number Publication date
FR2998973B1 (fr) 2015-01-23

Similar Documents

Publication Publication Date Title
EP2100161B1 (fr) Procede de traitement radar passif multivoies d&#39;un signal d&#39;opportunite en fm
EP2318854B1 (fr) Procede de protection d&#39;un utilisateur de recepteur de radionavigation vis-a-vis de mesures de pseudo-distances aberrantes
EP2176680B1 (fr) Procedes et dispositifs pour determiner la reponse impulsionnelle de canaux de propagation impliquant des emetteurs, des reflecteurs et des capteurs fixes ou mobiles
EP2869085B1 (fr) Procédé de détection d&#39;interférences dans un signal de radio-navigation par satellite par detection d&#39;une déformation de la fonction de corrélation
EP3238094A1 (fr) Procédé non linéaire d&#39;estimation d&#39;un mélange de signaux
EP3709058A1 (fr) Procédé de contrôle d&#39;intégrité d&#39;un signal de radionavigation par satellite et recepteur apte à implementer ce procédé
EP3412012B1 (fr) Procédé d&#39;estimation de paramètres de signaux contenus dans une bande de fréquences
EP2168337B1 (fr) Méthode de détection de signaux cyclostationnaires
EP0933885B1 (fr) Procédé de détection cyclique en diversité de polarisation de signaux radioélectriques numériques cyclostationnaires
EP2168336B1 (fr) Méthode de détection de présence de signaux étalés spectralement
EP2438461B1 (fr) Procédé de détection d&#39;un message émis par un interrogateur ou un répondeur en mode s
EP3036840B1 (fr) Procédé de détection d&#39;un signal électromagnétique par un réseau antennaire et dispositif mettant en oeuvre le procédé
FR2998973A1 (fr) Procede de determination des caracteristiques d un signal radar en presence d interferences
EP1085341B1 (fr) Procédé de traitement de signaux de navigation à trajets multiples dans un récepteur présentant une pluralité d&#39;antennes
EP1742376B1 (fr) Procédé de fusion d&#39;informations issues de récepteurs radioélectriques et dispositif de réception radioélectrique
EP2851703A1 (fr) Procédé réalisant conjointement la synchronisation, l&#39;identification, la mesure, l&#39;estimation du filtre de propagation et la localisation d&#39;émetteurs utiles et interferants
EP3036841B1 (fr) Procédé de détection d&#39;un signal électromagnétique par un réseau antennaire a diversité de diagramme et dispositif mettant en oeuvre le procédé
EP2869086B1 (fr) Procédé de détection d&#39;interférences dans un signal de radio-navigation par satellite basé sur la surveillance d&#39;un coefficient de corrélation temporel
EP2453261A1 (fr) Procédé de détection de la deformation d&#39;un signal GNSS
FR2998974A1 (fr) Procede de determination des caracteristiques parmi un ensemble des caracteristiques presumees d un signal radar en presence d interferences
EP3236280B1 (fr) Calcul automatique d&#39;une dimension d&#39;une plateforme mobile
EP3400457B1 (fr) Procédé d&#39;élimination d&#39;un signal en provenance d&#39;un radar de bord
WO2015144649A1 (fr) Procédé de détection d&#39;un signal cible dans un signal de mesure d&#39;un instrument embarqué dans un engin spatial et système de mesure
FR2955950A1 (fr) Procede et dispositif permettant de determiner un parametre central associe a des coefficients d&#39;identification obtenus par un algorithme de burg
EP1949548A1 (fr) Procede de detection de trajets en transmission impulsionnelle et dispositif correspondant

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 4

PLFP Fee payment

Year of fee payment: 5

PLFP Fee payment

Year of fee payment: 6

PLFP Fee payment

Year of fee payment: 7

PLFP Fee payment

Year of fee payment: 8

PLFP Fee payment

Year of fee payment: 9

PLFP Fee payment

Year of fee payment: 10

PLFP Fee payment

Year of fee payment: 11

PLFP Fee payment

Year of fee payment: 12