FR 3031257 A1 20160701 FR 1403062 A 20141231 -I^ 3031257 1 Procédé de désentrelacement de signaux radars La présente invention concerne un procédé de désentrelacement de signaux radars, le procédé comprenant : - la réception de signaux électromagnétiques par un récepteur et l'extraction des impulsions des signaux reçus, et - la formation de trains d'impulsions regroupant au moins trois impulsions espacées d'un même intervalle de répétition des impulsions, chaque train d'impulsions étant défini par l'intervalle de répétition des impulsions. La présente invention concerne également un dispositif de désentrelacement associé. L'un des enjeux de la guerre électronique consiste à intercepter les émissions radioélectriques en provenance de systèmes de détection tels que des émetteurs radar. La présence de nombreux émetteurs fait que les signaux interceptés sont entrelacés, c'est-à-dire que les signaux émis par un émetteur radar d'intérêt sont brouillés par d'autres signaux émanant du milieu électromagnétique ambiant. Il convient donc de désentrelacer les signaux interceptés pour séparer les différents signaux émis par les différents émetteurs radars. Or, des signaux émis par un même émetteur radar peuvent avoir des caractéristiques, définissant une forme d'onde, variables, notamment en termes d'intervalle de répétition des impulsions des signaux ou de fréquence porteuse des impulsions des signaux, ce qui rend le processus de désentrelacement complexe. A la grande richesse de formes d'ondes du monde électromagnétique correspond une grande variété de traitements de désentrelacement visant à extraire des impulsions d'une même forme d'onde du milieu électromagnétique ambiant. Plus particulièrement, le domaine technique, objet de ce procédé, concerne le désentrelacement des formes d'ondes dont l'intervalle de répétition des impulsions est moyen ou court (jusqu'à quelques centaines de micro-secondes). Ces formes d'ondes sont généralement constituées de plusieurs trains d'impulsions. Il est connu d'utiliser des extracteurs de signaux radars mettant en oeuvre un processus de désentrelacement de signaux radars en deux étapes. La première étape consiste en la formation de trains d'impulsions à partir de l'ensemble des impulsions interceptées. La deuxième étape consiste à regrouper les trains d'impulsions formés pour obtenir des signaux radars désentrelacés. La première étape de formation des trains d'impulsions exploite les informations statistiques des signaux interceptés telles que les fréquences des porteuses des 3031257 impulsions des signaux, les intervalles de répétition des impulsions et les directions d'arrivée des impulsions. La deuxième étape de regroupement des trains d'impulsions regroupe les trains d'impulsions formés suivant leur proximité pour former des signaux désentrelacés. 5 Néanmoins, les algorithmes existants de regroupement de trains d'impulsions n'offrent pas le même niveau de maturité que ceux de formation des trains d'impulsions. En particulier, l'étape de regroupement des trains d'impulsions est généralement abordée comme un problème de « clustering » ou de « classification de données », où chaque train est comparé à un autre selon un unique critère de distance. Pourtant, les formes 10 d'ondes formées par les trains d'impulsions sont d'une grande diversité, certaines familles de formes d'ondes pouvant avoir des caractéristiques totalement antagonistes entre elles. Par conséquent, un critère de distance unique peut mener à la reconstruction de signaux erronés. Le problème technique porte sur le regroupement de trains d'impulsions issus d'un 15 même signal radar en environnement électromagnétique dense où plusieurs formes d'ondes distinctes peuvent apparaître simultanément, la difficulté consistant à ne pas effectuer de mauvais regroupements de trains d'impulsions. Il existe donc un besoin pour un procédé de désentrelacement de signaux radars permettant de regrouper des trains d'impulsions avec une meilleure fiabilité, en limitant les 20 risques d'obtention d'un signal désentrelacé erroné, tout en étant de mise en oeuvre rapide. A cet effet, l'invention a pour objet un procédé de désentrelacement de signaux radars du type précité, dans lequel le procédé comporte en outre : - le regroupement des trains d'impulsions ayant un, même intervalle de répétitions 25 des impulsions suivant une loi de regroupement prédéfinie pour former des paliers d'impulsions, et - l'association des paliers d'impulsions suivant au moins une loi d'association prédéfinie pour obtenir des signaux radars désentrelacés formés de la concaténation des trains d'impulsions des paliers d'impulsions associés. 30 Suivant des modes particuliers de mise en oeuvre, le procédé de désentrelacement comprend une ou plusieurs des caractéristiques suivantes, prises isolément ou suivant toutes les combinaisons techniquement possibles : - chaque train d'impulsions est également défini par au moins un élément choisi dans un groupe constitué de : le temps d'arrivée de la première impulsion du train d'impulsions, le temps d'arrivée de la dernière impulsion du train d'impulsions, la 3031257 3 fréquence des impulsions du train d'impulsions, la durée des impulsions du train d'impulsions et la direction d'arrivée des impulsions du train d'impulsions. - le procédé comporte avant l'étape de regroupement, une étape de classification des trains d'impulsions selon leur fréquence porteuse pour obtenir deux classes de trains 5 d'impulsions : une classe regroupant les trains d'impulsions de fréquence porteuse fixe et l'autre classe regroupant les trains d'impulsions de fréquence porteuse variable, l'étape de regroupement étant mise en oeuvre pour chacune des deux classes de trains d'impulsions et permettant d'obtenir des paliers d'impulsions mono-fréquence à partir de la classe des trains d'impulsions mono-fréquence et des paliers d'impulsions agiles en 10 fréquence à partir de la classe des trains d'impulsions agiles en fréquence. - l'étape d'association comprend une phase de regroupement des paliers d'impulsions d'intervalles de répétition des impulsions différents et qui s'enchaînent dans le temps pour obtenir des groupes de paliers d'impulsions en commutation. - l'étape d'association comprend une phase de regroupement des paliers 15 d'impulsions de fréquence porteuse fixe, ayant des intervalles de répétitions des impulsions identiques et se superposant dans le temps pour obtenir des groupes de paliers d'impulsions en recouvrement. - chacune des lois de regroupement et d'association est mise en oeuvre par au moins un algorithme permettant d'obtenir des groupes à partir d'éléments, les éléments 20 désignant des trains d'impulsions lors de l'étape de regroupement et des paliers d'impulsions lors de l'étape d'association, les groupes désignant des palier d'impulsions lors de l'étape de regroupement et des groupes de paliers d'impulsions lors de l'étape d'association, l'algorithme comprenant : o le choix d'un élément de référence parmi un ensemble d'éléments, 25 o la suppression de l'élément de référence de l'ensemble d'éléments et l'ajout, dans un ensemble de groupes, d'un groupe de référence comportant l'élément de référence, o la sélection dans l'ensemble d'éléments, des éléments compatibles avec le groupe de référence selon un ensemble de critères pour obtenir un 30 ensemble d'éléments candidats, o l'évaluation de la distance entre le groupe de référence et chaque élément de l'ensemble d'éléments candidats, o l'annexion de l'élément de l'ensemble d'éléments candidats minimisant une distance au groupe de référence et la suppression de l'élément annexé de 35 l'ensemble d'éléments, 3031257 4 o la répétition des phases de sélection, d'évaluation et d'annexion tant que l'ensemble des éléments candidats comporte des éléments, et o la répétition de l'ensemble des phases précédentes tant que l'ensemble d'éléments comporte des éléments. 5 - l'élément de référence est l'élément de l'ensemble d'éléments dont le temps d'arrivée de la première impulsion est le plus petit. - l'ensemble de critères évalue la compatibilité des éléments de l'ensemble d'éléments avec le groupe de référence en fonction d'une ou de plusieurs caractéristiques, les caractéristiques étant choisies dans un groupe comprenant : la 10 direction d'arrivée des éléments, la superposition temporelle des éléments, la fréquence porteuse des éléments, la largeur d'impulsion des éléments, l'intervalle de répétition des impulsions des éléments, la phase des éléments et le nombre d'impulsions des éléments. - les critères sont choisis d'après des statistiques sur les caractéristiques des formes d'ondes radars d'une base de données. 15 - pour la loi de regroupement, la distance est l'écart en temps séparant la dernière impulsion du groupe de référence de la première impulsion de l'élément candidat de l'ensemble d'éléments candidats, et pour la loi d'association, la distance est un taux de recouvrement entre le groupe de référence et l'élément candidat de l'ensemble d'éléments candidats ou un score permettant de sélectionner parmi l'ensemble 20 d'éléments candidats l'élément partageant le plus de caractéristiques communes avec le groupe de référence. - le procédé comporte avant l'étape de regroupement, une étape de rejet des trains d'impulsions incohérents en termes d'intervalle de répétition des impulsions. - le procédé comporte avant l'étape de regroupement, une étape de rejet des 25 trains d'impulsions dont l'intervalle de répétition des impulsions est supérieur à un seuil d'intervalle de répétition des impulsions et dont le nombre d'impulsions est inférieur à un seuil de nombre d'impulsions. L'invention concerne aussi un dispositif de désentrelacement de signaux radars comprenant : 30 - un récepteur propre à recevoir des signaux électromagnétiques, - une unité de traitement numérique du signal propre à extraire les impulsions des signaux reçus par le récepteur, et - un support lisible d'informations, sur lequel est mémorisé un programme d'ordinateur comprenant des instructions de programme, le programme d'ordinateur étant 35 chargeable sur une unité de traitement de données et adapté pour entraîner la mise en 3031257 5 oeuvre d'un procédé selon l'une quelconque des revendications précédente lorsque le programme d'ordinateur est mis en oeuvre sur l'unité de traitement des données. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description qui suit de modes de réalisation de l'invention, donnée à titre d'exemple 5 uniquement et en référence aux dessins qui sont : - figure 1, une représentation schématique d'un dispositif de désentrelacement de signaux radars permettant la mise en oeuvre d'un procédé de désentrelacement selon l'invention, - figure 2, un organigramme d'un exemple de mise en oeuvre d'un procédé de 10 désentrelacement selon l'invention, - figure 3, une représentation schématique de trois trains d'impulsions, - figure 4, une représentation schématique d'un train d'impulsions et des grandeurs définissant ce train d'impulsions, - figure 5, une représentation schématique d'un signal et d'un palier d'impulsions 15 résultant du regroupement de trois trains d'impulsions de ce signal, - figure 6, une représentation schématique illustrant le regroupement de deux trains d'impulsions lors d'une étape de regroupement du procédé de désentrelacement selon l'invention, - figure 7, un organigramme du fonctionnement d'un algorithme mis en oeuvre au 20 cours du procédé de désentrelacement selon l'invention, - figure 8, une représentation schématique d'une phase d'élimination d'éléments incompatibles avec un groupe de référence, - figure 9, une représentation schématique d'une phase de réunion d'un élément compatible avec un groupe de référence lorsque la distance entre l'élément et le groupe 25 de référence est minimale par rapport aux autres éléments compatibles, - figure 10, une représentation schématique d'un groupe de paliers d'impulsions formés à partir de paliers d'impulsions d'une forme d'onde de type FMICW, - figure 11, une représentation schématique illustrant l'association de trois paliers d'impulsions pour former des groupes de paliers d'impulsions lors d'une autre étape 30 d'association du procédé de désentrelacement selon l'invention, - figure 12, une représentation schématique illustrant l'association de huit paliers d'impulsions pour former deux groupes de paliers d'impulsions lors d'une étape d'association du procédé de désentrelacement selon l'invention, et 3031257 6 - figure 13, une représentation schématique du processus de désentrelacement de signaux radars depuis la réception des signaux par un récepteur jusqu'à l'obtention de signaux désentrelacés. Un dispositif 10 de désentrelacement de signaux radars est représenté sur la 5 figure 1. Le dispositif 10 de désentrelacement de signaux radars est propre à mettre en oeuvre un procédé de désentrelacement de signaux radars. Les données en entrée du dispositif 10 de désentrelacement sont des mesures des caractéristiques des impulsions reçues par le dispositif 10. Les caractéristiques sont, par exemple, la fréquence porteuse des impulsions, la largeur des impulsions, la 10 puissance des impulsions, la direction d'arrivée des impulsions ou les périodes de répétition des impulsions. De telles mesures caractéristiques d'impulsions sont issues, à titre d'exemple : de la réception de signaux radar par un récepteur, puis de la numérisation des signaux et de l'extraction des impulsions des signaux par une unité de 15 traitement numérique du signal, de la réception de signaux analogiques créés en laboratoire simulant la réception de signaux radars, puis de la numérisation des signaux et de l'extraction des impulsions des signaux par une unité de traitement numérique du signal, 20 de la génération, via un logiciel, de données numériques simulant la réception et la numérisation de signaux radar, et de l'extraction des impulsions des signaux par une unité de traitement numérique du signal, ou de la génération, via un logiciel, d'impulsions qui sont ensuite enregistrées sur une interface de réception. 25 Comme illustré sur la figure 1, le dispositif 10 comprend un récepteur 12 d'ondes électromagnétiques, un calculateur 14 et un support lisible d'informations 16 en interaction avec le calculateur 14. Le récepteur 12 est propre à recevoir des signaux électromagnétiques émanant, par exemple, de systèmes de détection radioélectriques tels que des radars. Les signaux 30 électromagnétiques sont, par exemple, issus d'émetteurs radars ou sont des signaux analogiques crées en laboratoire et simulant des signaux radars. Le récepteur 12 est connecté au calculateur 14 et est propre à envoyer les signaux reçus par ce récepteur 12 au calculateur 14. Le récepteur 12 est, par exemple, une antenne. A titre d'exemple, le récepteur 12 35 est une antenne élémentaire, une antenne en réseaux, une antenne à réflecteur, une 3031257 7 antenne à polarisation circulaire, une antenne à guide d'onde, une antenne active, une antenne raccourcie, une antenne large bande, une antenne à patch, une antenne cadre ou une antenne boucle ou un système constitué d'une ou plusieurs des antennes précédentes. 5 Le calculateur 14 est propre à recevoir des signaux en provenance par exemple du récepteur 12 ou des données numériques en provenance de logiciels de simulation. Le calculateur 14 est un ordinateur comportant un processeur 18 et, optionnellement, une interface homme-machine 20 et une unité d'affichage 22. Le calculateur 14 comprend, en outre, une unité de traitement numérique du 10 signal 23 et optionnellement une interface 32. Le processeur 18 comprend une unité de traitement de données 24, des mémoires 26 et un lecteur 28 de support d'informations. Le lecteur 28 de support d'informations est propre à recevoir et à lire le support lisible d'informations 16. Le lecteur 28 de support d'informations est relié à l'unité de 15 traitement de données 24.Le support lisible d'informations 16 est un support lisible par le lecteur 28 de support d'informations. Le support lisible d'informations 16 est un médium adapté à mémoriser des instructions électroniques et capable d'être couplé à un bus d'un système informatique. A titre d'exemple, le support lisible d'informations 16 est une disquette ou disque souple (de la dénomination anglaise floppy disc), un disque optique, 20 un CD-ROM, un disque magnéto-optique, une mémoire ROM, une mémoire RAM, une mémoire EPROM (acronyme de l'anglais Erasable Programmable Read-Only Memory), une mémoire EEPROM (acronyme de l'anglais Electrically-Erasable Programmable ReadOnly Memory), une carte magnétique ou une carte optique. Sur le support lisible d'informations 16 est mémorisé un produit programme 25 d'ordinateur comprenant des instructions de programme. Le programme d'ordinateur est chargeable sur l'unité de traitement de données 24 et est adapté pour entraîner la mise en oeuvre d'un procédé de désentrelacement de signaux radars selon l'invention. L'interface homme-machine 20 est, par exemple, un clavier. L'unité d'affichage 22 est, par exemple, un écran. 30 L'unité de traitement numérique du signal 23 est configurée pour traiter numériquement les signaux reçus par le calculateur 14. Plus précisément, l'unité de traitement numérique du signal 23 est configurée pour numériser les signaux reçus, extraire les impulsions des signaux et mesurer les caractéristiques de chaque impulsion extraite. 3031257 8 En variante, l'unité de traitement numérique du signal 23 est également configurée pour générer des données numériques simulant la réception et la numérisation de signaux radars ou pour générer directement des impulsions. L'interface 32 permet, d'une part de mémoriser des impulsions résultant du 5 traitement numérique du signal effectué par l'unité de traitement numérique du signal 23 et d'autre part de recevoir des impulsions préalablement mémorisées ou de recevoir des impulsions générées par un logiciel de simulation. Le fonctionnement du dispositif de désentrelacement de signaux radars est maintenant décrit en référence à la figure 2 qui est un organigramme d'un exemple de 10 mise en oeuvre d'un procédé de désentrelacement de signaux radars selon l'invention. Le procédé de désentrelacement comprend une étape de réception 100 d'une pluralité de signaux électromagnétiques par le récepteur 12, de numérisation des signaux reçus et d'extraction des impulsions ..., lm électromagnétiques des signaux numérisés par traitement numérique du signal. La numérisation des signaux et l'extraction des 15 impulsions I,, ..., lm sont effectués par l'unité de traitement numérique du signal 23 et sont éventuellement enregistrées dans l'interface 32. De telles impulsions I,, ..., Im électromagnétiques émanent ou sont représentatives, notamment, de signaux émis par des systèmes de détection radioélectriques tels que des radars. En variante, les impulsions 1' sont déjà enregistrées dans l'interface 32. 20 Les étapes du procédé de désentrelacement décrites dans ce qui suit sont mises en oeuvre par le calculateur 14 en interaction avec le support lisible d'informations 16. Le procédé de désentrelacement comprend ensuite une étape 110 de formation de trains d'impulsions T1, Tn à partir des impulsions lm reçues par le récepteur 12 ou enregistrées dans l'interface 32. 25 Lors de l'étape de formation 110, au moins trois impulsions 12, 13 sont regroupées pour former un train d'impulsions T. Le critère de formation des trains d'impulsions T1, Tri consiste à regrouper dans un même train d'impulsions TX, les impulsions lm ayant un même intervalle de répétition des impulsions PRI. Il est entendu par le terme « intervalle de répétition des impulsions », la période séparant deux 30 impulsions consécutives. En d'autres termes, les impulsions ..., I, regroupées dans un même train d'impulsions TX sont telles que deux impulsions consécutives quelconques du train d'impulsions T, sont espacées du même intervalle de répétition des impulsions PRI que deux autres impulsions consécutives quelconques du train d'impulsions T. Les trains d'impulsions T1, ..., T, formés constituent un ensemble ET de trains 35 d'impulsions T1, Tn. Les trains d'impulsions T1, Tn formés sont définis par l'intervalle 3031257 9 de répétition des impulsions PRI séparant deux impulsions consécutives du train d'impulsions T1, ..., Tn. Il est illustré sur la figure 3, trois trains d'impulsions T1, T2 et T3 ayant des intervalles de répétition des impulsions, PRI,, PRI2 et PRIS, éventuellement différents les 5 uns des autres. Comme illustré sur la figure 4, un train d'impulsions Tx est également défini par le temps d'arrivée de la première impulsion du train d'impulsions Tx, appelé TOAdéb et le temps d'arrivée de la dernière impulsion 1m du train d'impulsions Tx, appelé TOAtin- En outre, chaque train d'impulsions Tx est éventuellement défini par au moins un 10 élément choisi dans un groupe constitué de : la fréquence des impulsions ..., lm du train d'impulsions Tx, la durée des impulsions I,, ..., lm du train d'impulsions Tx et la direction d'arrivée des impulsions ..., lm du train d'impulsions Tx. Le procédé de désentrelacement comprend ensuite une étape 120 de rejet des trains d'impulsions T1, ..., T, incohérents. 15 L'étape de rejet 120 vérifie la cohérence en intervalle de répétition des impulsions PRI des trains d'impulsions T1, ..., T, formés lors de l'étape de formation 110. En effet, certains trains d'impulsions Tx peuvent avoir été formés par le mélange d'au moins deux intervalles de répétition des impulsions PRI proches mais néanmoins différents. Un train d'impulsions Tx incohérent en termes d'intervalle de répétition des 20 impulsions PRI est, par exemple, détecté grâce à un test statistique se basant sur les intervalles séparant toutes les impulsions du train Tx. Un test du khi2 peut par exemple être utilisé. Les trains d'impulsions Tx jugés incohérents en termes d'intervalle de répétition des impulsions PRI sont éliminés de l'ensemble ET de trains d'impulsions T1, - - -, Tn- 25 Au cours de cette étape de rejet 120, optionnellement, la cohérence des trains d'impulsions T1, ..., Tri en terme de fréquence porteuse des impulsions est également vérifiée. Il est entendu par le terme « fréquence porteuse d'une impulsion », la fréquence de la porteuse de l'impulsion, la porteuse étant une onde modulée par un signal d'entrée. En effet, certains trains d'impulsions Tx peuvent avoir été formés par le mélange d'au 30 moins deux fréquences porteuses proches mais néanmoins différentes. Un train d'impulsions Tx incohérent en termes de fréquence porteuse est, par exemple, détecté grâce à un test statistique se basant sur les fréquences des impulsions du train Tx. Un test du khi2 peut par exemple être utilisé. Les trains d'impulsions Tx jugés incohérents en terme de fréquence porteuse des impulsions sont éliminés de l'ensemble ET de trains 35 d'impulsions T1, Tn. 3031257 10 Dans le cas où les intervalles de répétition des impulsions PRI, respectivement les fréquences porteuses des impulsions, sont modélisés comme des variables gaussiennes et indépendantes entre elles, la cohérence en intervalle de répétition des impulsions PRI, respectivement en fréquence porteuse des impulsions, est évaluée par un test statistique 5 du khi deux. Il est entendu par l'expression « test du khi deux », abrégé en « test du x2 », ou en « test du khi2 », un test statistique permettant de tester l'adéquation d'une série de données à une famille de lois de probabilités ou de tester l'indépendance entre deux variables aléatoires. Le procédé de désentrelacement comprend ensuite une étape 130 de rejet des 10 trains d'impulsions Tx n'appartenant pas à la catégorie des formes d'ondes de type haute ou moyenne fréquence de récurrence. Une forme d'onde regroupe des trains d'impulsions T1, Tn ayant des caractéristiques communes en terme, notamment, d'intervalle de répétition des impulsions PRI, de fréquence porteuse des impulsions et de nombre d'impulsions. Par exemple, l'ensemble des trains d'impulsions T1, ..., Tri ayant des 15 intervalles de répétition des impulsions PRI courts et un grand nombre d'impulsions constitue une forme d'onde de type haute ou moyenne fréquence de récurrence. A l'inverse, l'ensemble ET des trains d'impulsions T1, ..., T,, ayant des intervalles de répétition des impulsions PRI élevé et un nombre d'impulsions faible constitue une autre forme d'onde de type basse fréquence de récurrence. 20 L'étape de rejet 130 comprend la définition de seuils caractéristiques d'une forme d'onde : un seuil d'intervalle de répétition des impulsions SpRi et un seuil de nombre d'impulsions S'impulsions. Le seuil d'intervalle de répétition des impulsions SpRi et le seuil de nombre d'impulsions Simpuisions sont définis à partir d'une base de données représentative d'un 25 certain nombre de formes d'ondes radars d'intérêt. Dans une variante de réalisation, cette base de données de formes d'ondes est stockée dans la mémoire 26 du processeur 18. Dans une autre variante de réalisation, seuls les seuils déterminés à partir de la base de données sont stockés dans la mémoire 26 du processeur 18. Le seuil SpRi est choisi de façon à être supérieur aux valeurs majoritairement 30 observées dans la base de données de formes d'ondes. Le seuil d'intervalle de répétition des impulsions SpRI est, par exemple, compris au sens large entre 1 microsecondes (ps) et 1 millisecondes (ms) pour tester l'appartenance aux formes d'ondes de type haute ou moyenne fréquence de récurrence. Le seuil S'impulsions est choisi de façon à être inférieur aux valeurs majoritairement 35 observées dans la base de données de formes d'ondes. Le seuil de nombre d'impulsions 3031257 11 S'impulsions est, par exemple, compris au sens large entre 1 impulsion et 100 impulsions pour tester l'appartenance aux formes d'ondes de type haute ou moyenne fréquence de récurrence. Dans ce cas, de tels seuils d'intervalle de répétition des impulsions SARI et de nombre d'impulsions Simpuls,ons permettent d'exclure les trains d'impulsions T1, Tn de 5 type basse fréquence de récurrence. L'étape de rejet 130 comprend ensuite le rejet des trains d'impulsions Ti, Tn dont l'intervalle de répétition des impulsions PRI est supérieur au seuil d'intervalle de répétition des impulsions SARI. L'étape de rejet 130 comprend également le rejet des trains d'impulsions T1, ..., Tn dont le nombre d'impulsions est inférieur au seuil de nombre 10 d'impulsions S'impulsions. Le seuil d'intervalle de répétition des impulsions SARI et le seuil de nombre d'impulsions S'impulsions forment ainsi un gabarit de sorte que tout train d'impulsions TX sortant de ce gabarit est rejeté. Ainsi, les trains d'impulsions TX rejetés sont retirés de l'ensemble ET de trains d'impulsions T1, Tn et sont traités par la suite par un autre 15 procédé de désentrelacement spécifique. Le procédé de désentrelacennent comprend ensuite une étape 140 de classification des trains d'impulsions T1, Tn. L'étape de classification 140 consiste en un tri des trains d'impulsions T1, ..., T, de l'ensemble ET de trains d'impulsions T1, ..., I-, suivant la fréquence porteuse des 20 impulsions de chaque train d'impulsions T1, Tn. Au cours de cette étape de classification 140, les trains d'impulsions T1, Tn sont répartis en deux classes Ci, C2 de trains d'impulsions T1, Tn. La première classe C1 regroupe les trains d'impulsions T1, de fréquence porteuse fixe, c'est-à-dire les trains d'impulsions T1, ..., T, formés d'impulsions ayant la même fréquence porteuse à 25 une incertitude de mesure près. L'incertitude de mesure est, par exemple, égale à plus ou moins 5 pourcents (%) de la fréquence porteuse. Les trains d'impulsions T1, Tn de fréquence porteuse fixe sont également appelés trains d'impulsions T1, Tn mono- fréquence. La deuxième classe C2 de trains d'impulsions T1, ..., T,, regroupe les trains 30 d'impulsions T1, ..., Tr, de fréquence porteuse variable, c'est-à-dire les trains d'impulsions T1, ..., Tr, formés d'impulsions de fréquences porteuses différentes. Les trains d'impulsions T1, ..., 1-0 de fréquence porteuse variable sont également appelés trains d'impulsions T1, ..., T0 agiles en fréquence. La fréquence porteuse des impulsions des trains d'impulsions T1, ..., T, agiles en fréquence est généralement aléatoire ou pseudo 35 aléatoire. 3031257 12 Le procédé de désentrelacement comprend, ensuite, une étape 150 de regroupement des trains d'impulsions T1, T, de l'ensemble ET de trains d'impulsions T1, .- - , Tn- Lors de l'étape de regroupement 150, les trains d'impulsions T1, Tn ayant un 5 même intervalle de répétition des impulsions PRI sont regroupés suivant au moins une loi de regroupement prédéfinie pour former des paliers d'impulsions P1, - - Pp. Il est entendu par l'expression « palier d'impulsions », un ensemble d'impulsions appartenant à au moins un train d'impulsions formé lors de l'étape de regroupement 150. Il est illustré sur la figure 5, un exemple de signal reçu R par le récepteur 12 et un 10 palier d'impulsions P1 formé par le regroupement de trains d'impulsions T4, T5 et T6 de ce signal. De même, il est illustré sur la figure 6, un palier d'impulsions P2 issu du regroupement des trains d'impulsions T6 et T7. Plus précisément, l'étape 150 de regroupement des trains d'impulsions T1, Tn est mise en oeuvre, d'une part, pour les trains d'impulsions T1, Tn de fréquence 15 porteuse fixe lors d'une première sous-étape 160 pour obtenir des paliers d'impulsions P1, Pp de fréquence porteuse fixe. La loi de regroupement utilisée est appelée « loi de regroupement fixe ». L'étape 150 de regroupement des trains d'impulsions T1, Tn est mise en oeuvre, d'autre part, pour les trains d'impulsions T1, ..., Tn de fréquence porteuse variable lors d'une deuxième sous-étape 170 pour obtenir des paliers d'impulsions 20 P1, Pp de fréquence porteuse variable. La loi de regroupement utilisée est appelée «'loi de regroupement variable ». Les trains d'impulsions T1, Tn regroupés lors de l'étape de regroupement 150 n'ont pas été mis ensemble pour former un seul train d'impulsions T1, Tn lors de l'étape de formation 110 car de tels trains d'impulsions T1, Tn présentaient, par 25 exemple, des mitages M, c'est-à-dire des impulsions manquantes, du fait de raisons techniques. Ces raisons techniques dépendent notamment de limitations d'écoute du récepteur 12, de la qualité des mesures par l'unité de traitement numérique du signal 23, de la proximité du récepteur 12 par rapport à l'émetteur, du positionnement du récepteur 12 par rapport à l'émetteur et de perturbations électromagnétiques. 30 Les mitages sont, par exemple, représentés sur les trains d'impulsion T6 et T7. Les impulsions manquantes sont reconstituées lors de la formation du palier d'impulsions P2. De manière générale, une loi de regroupement consiste en un algorithme dont la structure générale est décrite en référence à la figure 7. Initialement, l'algorithme comprend une phase 150A de choix d'un élément de 35 référence eréf parmi un ensemble d'éléments E. 3031257 13 L'algorithme comprend, ensuite, une phase 150B de suppression de l'élément de référence eréf de l'ensemble d'éléments E et d'ajout, dans un ensemble de groupes Eg, d'un groupe de référence gréf comportant l'élément de référence eréf- L'algorithme comprend, ensuite, une phase 150C de sélection dans l'ensemble 5 d'éléments E, des éléments compatibles avec le groupe de référence gréf selon un ensemble de critères C. Les éléments compatibles avec le groupe de référence gréf forment un ensemble d'éléments candidats Ec. Les critères C sont, par exemple, choisis d'après des statistiques sur les caractéristiques de formes d'ondes radars de la base de données de formes d'ondes. 10 Par exemple, la figure 8 illustre un groupe G3 déjà formé et un groupe G4 en cours de formation qui est donc à ce stade le groupe de référence gréf de l'algorithme. Les éléments el, e2, e3 et e4 sont les éléments de l'ensemble d'éléments E de l'algorithme. Comme illustré par une croix sur cette figure 8, les éléments el et e3 sont incompatibles avec le groupe G4 en formation. A l'inverse, comme illustré par une coche sur cette figure 15 8, les éléments e2 et e4 sont compatibles pour être réunis avec le groupe G4 en formation. L'ensemble des éléments candidats Ec de l'algorithme comprend donc deux éléments : e2 et e4. L'algorithme comprend, ensuite, une phase 150D d'évaluation de la distance entre le groupe de référence gréf et chaque élément de l'ensemble d'éléments candidats Ec. 20 Puis, l'algorithme comprend une phase 150E d'annexion d'un élément de l'ensemble d'éléments candidats E, au groùpe de référence gréf. L'élément annexé est retiré de l'ensemble d'éléments E. L'élément de l'ensemble d'éléments candidats E, réuni avec le groupe de référence gréf est l'élément ayant la plus petite distance d. Il est entendu par l'expression « annexion d'un élément à un groupe », le rattachement de l'élément au 25 groupe : à l'issue de l'annexion le groupe final est formé du groupe initial et de l'élément annexé qui est donc réuni au groupe initial. Par exemple, la figure 9 illustre la détermination de l'élément à réunir au groupe de référence gréf, c'est-à-dire au groupe G4 en formation, parmi l'ensemble d'éléments candidats Ec. La figure 9 comprend les mêmes groupes et éléments que ceux de la 30 figure 8, mais représentés de manière différente. Ainsi, le groupe G3 déjà formé comprend trois éléments représentés par trois ronds. Le groupe G4 en cours de formation comprend déjà quatre éléments représentés par quatre ronds. Comme visible sur la figure 9, la distance d2 séparant l'élément e2 du groupe G4 est inférieure à la distance d4 séparant l'élément e4 du groupe G4. Par conséquent, l'élément de l'ensemble d'éléments 3031257 14 candidats En qui minimise la distance avec le groupe de référence gréf est l'élément e4. L'élément e4 sera donc réuni avec le groupe G4 en formation. Les phases de sélection d'éléments compatibles 150C, d'évaluation de distance 150D et d'annexion 150E sont, ensuite, répétées avec les éléments restant de 5 l'ensemble d'éléments E, tant que l'ensemble des éléments candidats En calculé à la phase de sélection d'éléments compatibles 150C, comporte des éléments. Enfin, les phases de choix 150A, de suppression et d'ajout 150B, de sélection d'éléments compatibles 150C, d'évaluation de distance 150D, et d'annexion 150E sont répétées tant que l'ensemble d'éléments E comporte des éléments. 10 La loi de regroupement fixe appliquée lors de la première sous-étape 160 est décrite de manière spécifique dans ce qui suit. Chaque élément de l'algorithme est un train d'impulsions T1, T,, de la première classe C1 des trains d'impulsions T1, ..., T, de fréquence porteuse fixe et chaque groupe de l'algorithme est un palier d'impulsions P1, Pp de fréquence porteuse fixe. Le groupe 15 de référence gréf de l'algorithme est un palier d'impulsions P1, Pp. En outre, l'ensemble d'éléments E initial de l'algorithme est l'ensemble des trains d'impulsions T1, Tr, de la première classe C1 des trains d'impulsions T1, Tn de fréquence porteuse fixe. Lors de la phase de choix 150A, l'élément choisi parmi l'ensemble d'éléments E est l'élément dont le temps d'arrivée de la première impulsion est le plus petit. En d'autres 20 termes, l'élément choisi est l'élément de l'ensemble d'éléments E arrivant le premier sur le récepteur 12 ou encore l'élément enregistré sur l'interface 32 avec le temps d'arrivée de la première impulsion TOA le plus petit. L'ensemble de critères C comprend les critères suivants : un critère de direction d'arrivée, un critère temporel, un critère fréquentiel, un critère de largeur d'impulsion, un 25 critère d'intervalle de répétition des impulsions et un critère de phase. Seuls les éléments de l'ensemble d'éléments E vérifiant l'intégralité des critères de l'ensemble de critères C sont compatibles avec le groupe de référence gréf et sont ajoutés dans l'ensemble d'éléments candidats En. En variante, l'ensemble de critères C ne comprend pas l'intégralité des critères 30 énoncés ou comprend des critères différents de ceux énoncés. Le critère de direction d'arrivée permet de tester la compatibilité de chaque élément de l'ensemble d'éléments E avec le groupe de référence gréf en fonction de leur direction d'arrivée DOA. Le critère de direction d'arrivée stipule que pour être compatibles, le groupe de référence gréf et l'élément doivent être formés d'impulsions ayant la même 35 direction d'arrivée DOA. 3031257 15 Le critère de direction d'arrivée est vérifié pour un élément de l'ensemble d'éléments E lorsque l'élément vérifie un test de comparaison des directions d'arrivées DOA de l'élément et du groupe de référence gréf. Le test de comparaison se base sur un modèle statistique de la direction d'arrivée DOA et est adapté pour prendre en compte la 5 proportion de valeurs aberrantes. Pour vérifier le test de comparaison, l'élément de l'ensemble d'éléments E doit vérifier l'équation eqi suivante : 2 Ir (1 1 (D0Ai. - DOA2)2 5 2 * * * - + *erf-1(P(1/011/0)) 1- 2r 2 ni. n2 où DOA1 est la direction d'arrivée du groupe de référence gréf, DOA2 est la direction d'arrivée de l'élément de l'ensemble d'éléments E à tester, 62 est la variance de la direction d'arrivée des impulsions 10 r est la proportion de mesures de direction d'arrivée dégradées, n1 est le nombre d'impulsions présentes dans le groupe de référence gréf, n2 est le nombre d'impulsions présentes dans l'élément de l'ensemble d'éléments E à tester, errl est la réciproque de la fonction d'erreur, la fonction d'erreur étant donnée par 15 l'équation erf(x) = 712 g e-t2 dt, et P(1101110) est la probabilité de détection d'égalité entre DOA1 et DOA2, H, désignant l'hypothèse que les valeurs DOA1 et DOA2 soient égales, et 1)(1101H0) désignant la probabilité de faire le choix H, sachant que l'on est dans le cas H,. En variante, un autre test de comparaison consiste à calculer, puis à comparer, la 20 direction d'arrivée DOA moyenne des impulsions du groupe de référence gréf par rapport à la direction d'arrivée DOA moyenne des impulsions de chaque élément de l'ensemble d'éléments E. En variante, encore un autre test de comparaison consiste à calculer, puis à comparer, la direction d'arrivée DOA médiane des impulsions du groupe de référence gréf 25 par rapport à la direction d'arrivée DOA médiane des impulsions de chaque élément de l'ensemble d'éléments E. L'emploi de la médiane et non pas de la moyenne permet d'éviter la prise en compte de données aberrantes fortement présentes dans les mesures de directions d'arrivée DOA. Le critère temporel permet de tester la compatibilité temporelle de chaque élément 30 de l'ensemble d'éléments E avec le groupe de référence gréf. Le critère temporel stipule qu'un élément de l'ensemble d'éléments E qui se superpose en temps avec le groupe de référence gréf est incompatible avec le groupe de référence gréf. Le critère temporel stipule, en outre, qu'un écart en temps trop important entre un élément de l'ensemble 3031257 16 d'éléments E et le groupe de référence gréf entraîne également une incompatibilité de l'élément avec le groupe de référence gréf. En effet, un écart en temps important entre le groupe de référence gréf et l'élément peut signifier qu'il n'y a pas eu d'impulsions émises. Dans ce cas, le groupe de référence gréf et l'élément, même s'ils appartiennent au même 5 émetteur, sont distincts. L'écart en temps maximal séparant le groupe de référence gréf d'un élément compatible de l'ensemble d'élément E est, par exemple, un multiple de l'intervalle de répétition des impulsions PRI moyen du groupe de référence gréf. L'écart en temps maximal est, par exemple, égal à vingt fois la valeur moyenne de l'intervalle de répétition 10 des impulsions PRI du groupe de référence gréf. La valeur du multiple est fixée de manière empirique ou à partir de la base de données de formes d'ondes. En variante, d'autres valeurs de multiples ou d'autres valeurs de seuils ne dépendant pas de l'intervalle de répétition des impulsions PRI peuvent encore être envisagées. Le critère fréquentiel permet de tester la compatibilité fréquentielle des éléments 15 de l'ensemble d'éléments E avec le groupe de référence gréf. Le critère fréquentiel stipule que pour être compatibles deux éléments doivent être formés d'impulsions ayant la même fréquence porteuse f. Le critère fréquentiel est donc vérifié pour un élément de l'ensemble d'éléments E lorsque l'élément vérifie un test de comparaison des fréquences porteuses des impulsions de l'élément avec les fréquences porteuses des impulsions du groupe de 20 référence gréf. Le test de comparaison consiste à comparer la fréquence porteuse moyenne des impulsions du groupe de référence gréf par rapport à la fréquence porteuse moyenne des impulsions de chaque élément de l'ensemble d'éléments E. Un tel test de comparaison repose sur le fait qu'une mesure de fréquence suit un modèle gaussien de variance 25 connue. Néanmoins, dans l'hypothèse où les mesures de fréquence porteuses sont obtenues en s'aidant de fenêtres fréquentielles, les possibilités d'exclusion de certaines valeurs fréquentielles sont à prévoir. Dans ce cas, la distribution des fréquences porteuses mesurées n'est plus gaussienne, car partiellement tronquée, et dès lors la comparaison des moyennes est biaisée. 30 En variante, notamment dans le cas où la distribution des fréquences n'est plus gaussienne, un test du x2 est utilisé. Pour vérifier le critère fréquentiel, l'élément de l'ensemble d'éléments E doit vérifier l'équation eq2 suivante : nl n2 2 2 N2 f1,i f2,i * Xni+n2-1,P(Ho 'Ho) i=o i=o 3031257 17 où hi est la fréquence porteuse de l'impulsion d'ordre i du groupe de référence gréf, f2,i est la fréquence porteuse de l'impulsion d'ordre i de l'élément de l'ensemble d'éléments E à comparer, 5 n1 est le nombre d'impulsions du groupe de référence gréf, n2 est le nombre d'impulsions de l'élément de l'ensemble d'éléments E à comparer, cri est la variance de la distribution des mesures de fréquence porteuse des impulsions, 10 ep est le quantile xp pour la probabilité p de la distribution du khi2 à y degrés de liberté, c'est à dire que si une variable aléatoire X suit une loi du khi2 à y degrés de liberté, la probabilité d'avoir )( xp est p, et P(1-101110) est la probabilité de détection d'égalité entre ri/10A2,i et E720f2i, H, désignant l'hypothèse que les deux valeurs de fréquence soient égales, et P(1101110) 15 désignant la probabilité de faire le choix Ho sachant que l'on est dans le cas Ho. Le critère de largeur d'impulsion permet de tester la compatibilité en largeur d'impulsion LI des éléments de l'ensemble d'éléments E avec le groupe de référence gréf. Le critère de largeur d'impulsion stipule que pour être compatibles le groupe de référence gréf et l'élément doivent être formés d'impulsions ayant la même largeur d'impulsion. 20 Comme illustrée sur la figure 4, la largeur d'impulsion LI d'une impulsion est la durée de l'impulsion. Le critère de largeur d'impulsion est vérifié pour un élément de l'ensemble d'éléments E lorsque les largeurs d'impulsions LI des impulsions de l'élément sont du même ordre de grandeur que les largeurs d'impulsions LI des impulsions du groupe de 25 référence gréf. Il est entendu par le terme « même ordre de grandeur », un écart relatif moindre qu'un seuil déterminé à partir de la base de données de formes d'ondes, ou valant par exemple 50%. Le critère d'intervalle de répétition des impulsions PRI permet de tester la compatibilité de l'intervalle de répétition des impulsions PRI de chaque élément de 30 l'ensemble d'éléments E avec le groupe de référence gréf. Le critère d'intervalle de répétition des impulsions PRI stipule que les impulsions de chaque élément d'un même groupe doivent avoir le même intervalle de répétition des impulsions PRI. Le critère d'intervalle de répétition des impulsions PRI est vérifié pour un élément de l'ensemble d'éléments E lorsque l'élément vérifie un test de comparaison de l'intervalle de répétition 3031257 18 des impulsions PRI de l'élément avec l'intervalle de répétition des impulsions PRI du groupe de référence gréf. Le test de comparaison consiste en premier lieu à estimer, pour chaque élément de l'ensemble d'éléments E, l'intervalle de répétition des impulsions PRI théorique issu de 5 la réunion des impulsions de l'élément avec les impulsions du groupe de référence gréf- L'intervalle de répétition des impulsions PRI théorique est estimé, par exemple, par la méthode de Fourier ou par la méthode des moindres carrés. Le test de comparaison consiste, ensuite, à effectuer un test d'adéquation pour vérifier, par exemple, si la somme des erreurs quadratiques est supérieure ou inférieure à un seuil. Le test d'adéquation est 10 par exemple un test du x2. Ce test de comparaison présente l'avantage de déterminer, en outre, l'intervalle de répétition des impulsions PRI du groupe obtenu lorsque l'élément de l'ensemble d'éléments E est réuni avec le groupe de référence gréf. Cependant, un tel test de comparaison est relativement lent et consommateur de ressources. En variante est présentée un test de comparaison moins long et moins 15 consommateur de ressources. Pour vérifier le critère d'intervalle de répétition des impulsions PRI, l'élément de l'ensemble d'éléments E doit vérifier l'équation eq3 suivante : (PRI1 - PRI2)2 5 4* oloa * (erf-1(Pd))2 où PR/1 est l'intervalle de répétition des impulsions du groupe de référence gréf, PR/2 est l'intervalle de répétition des impulsions de l'élément de l'ensemble d'éléments E à comparer, 20 ''TOA est la variance de la distribution des mesures du temps d'arrivée des impulsions TOA, erf-1 est la réciproque de la fonction d'erreur, et Pd est est la probabilité de détection d'égalité de deux valeurs de PRI. En variante, si les précisions des intervalles de répétition des impulsions de 25 chaque élément ont été estimées, encore un autre test est proposé. Pour vérifier le critère d'intervalle de répétition des impulsions PRI, l'élément de l'ensemble d'éléments E doit vérifier l'équation eq4 suivante : (PRI1 - PRI2)2 5 2 * (var(PRI1) + var(PRI2)) * (erf-1(Pd))2 où var(PRI1) est la variance de l'intervalle de répétition des impulsions du groupe de référence gréf, 30 var(PRI2) est la variance de l'intervalle de répétition des impulsions de l'élément de l'ensemble d'éléments E à comparer, et les autres notations sont celles de l'équation eq3. 3031257 19 Le critère d'intervalle de répétition des impulsions tient également compte du fait que certains paliers peuvent avoir des impulsions manquantes et que par conséquent l'intervalle de répétition des impulsions PRI mesuré est susceptible d'être un multiple de l'intervalle de répétition des impulsions PRI réel. 5 Ainsi, le critère d'intervalle de répétition des impulsions PRI compare également les multiples des intervalles de répétition des impulsions de chaque élément de l'ensemble d'éléments E avec l'intervalle de répétition des impulsions PRI du groupe de référence gréf, et inversement. Par exemple, le test de comparaison compare l'intervalle de répétition des impulsions PRI du groupe de référence gréf avec respectivement deux 10 fois, trois fois, quatre fois, l'intervalle de répétition des impulsions PRI de chaque élément de l'ensemble d'éléments E ou encore l'intervalle de répétition des impulsions PRI de chaque élément de l'ensemble d'éléments E avec respectivement deux fois, trois fois, quatre fois, l'intervalle de répétition des impulsions PRI du groupe de référence gréf. Le nombre de multiples à comparer dépend des conditions d'écoute et du pourcentage de 15 mitage M des éléments. Néanmoins, trois est une limite de multiple raisonnable. Il est, en effet, peu probable qu'un élément mité conserve une impulsion sur quatre à intervalle régulier. Le critère de phase permet de tester la compatibilité en phase des temps d'arrivée des impulsions de chaque élément de l'ensemble d'éléments E avec le groupe de 20 référence gréf. Les phases en temps d'arrivée des impulsions de deux éléments ont été préalablement calculées par rapport à un temps de référence tref. Cependant, si le temps de référence tref est éloigné de la position temporelle des éléments à comparer, cela peut s'avérer problématique pour la comparaison. D'un point de vue mathématique, la 25 précision de la phase en temps d'arrivée des impulsions est obtenue par la méthode des moindres carrés ou encore par la méthode de Fourier. La précision de la phase en temps d'arrivée des impulsions vérifie l'équation eq5 suivante : E Icî 2 var ésp = n E ki2 ' 'TOA - (Eki)2 où rp est l'estimée de la phase en fonction du temps d'arrivée des impulsions, t,11- K est le rang ou l'ordre de l'impulsion i par rapport au temps zéro tel que ki = [PRI 30 avec ti le temps d'arrivée de l'impulsion i et PRI l'intervalle de répétition des impulsions et [ ] désignant la partie entière, n est le nombre d'impulsions du groupe obtenu, et 3031257 20 672,0A est la variance de la distribution des mesures du temps d'arrivée des impulsions TOA. De l'équation eq6, il est possible de montrer que la précision de la phase est meilleure, c'est-à-dire, plus petite, quand les ordres des impulsions sont faibles, c'est-à- 5 dire quand le temps de référence tref est proche de la position temporelle de l'élément. Le critère de phase suppose que les deux éléments à comparer vérifient en premier lieu le critère d'intervalle de répétition des impulsions PRI. Si tel n'est pas le cas, le groupe de référence gréf et l'élément sont incompatibles en phase. Le critère de phase est vérifié pour un élément de l'ensemble d'éléments E lorsque l'élément vérifie un test de 10 comparaison des phases de l'élément et du groupe de référence gréf. Le test de comparaison comporte en premier lieu la détermination d'un nouveau temps de référence tref, qui sera positionné au centre du la position temporelle des deux éléments à comparer. Le test de comparaison comprend ensuite le calcul des phases du groupe de référence gréf et de l'élément dans ce nouveau temps de référence tref. Enfin, 15 l'élément de l'ensemble d'éléments E est compatible avec le groupe de référence gréf lorsqu'il vérifie l'équation eq6 suivante : (çPioa,1 4oa,2)2 5 2(varCeto.,i) + varFPioa,2)) * (erf-1(P(H011-10)))2 OÙ (pitoaa = ((Ptoa,i tref) * mod(PRIi) et ePtoa,2 = (ePtoa,2 tref) * mod(PRI2), toafin,i+toadeb,2 OÙ tref = est la moyenne entre le temps d'arrivée de la dernière 2 20 impulsion du groupe de référence gréf et le temps d'arrivée de la première impulsion de l'élément de l'ensemble d'éléments E à comparer, (x)mod(y) indique que l'on calcule un modulo en retranchant y autant de fois que nécessaire à x pour que (x)mod(y) soit compris entre 0 et y, ePtoa,1 est la phase en temps d'arrivée des impulsions TOA du groupe de référence 25 gréf, ePtoa,2 est la phase en temps d'arrivée des impulsions TOA de l'élément de l'ensemble d'éléments E à comparer, var désigne la variance, où el-ri est la réciproque de la fonction d'erreur, 3031257 21 et où P(1101Ho) est la probabilité de détection d'égalité entre les phases, H, désignant l'hypothèse que les phases soient égales, et P(1-101H0) désignant la probabilité de faire le choix H, sachant que l'on est dans le cas H,. Encore en variante, notamment lorsque les données de variance ne sont pas 5 accessibles, l'élément de l'ensemble d'éléments E est déclaré compatible avec le groupe de référence gréf lorsqu'il vérifie l'équation eq7 suivante : Pioaa (Ptoa,21 5 2 * Utoa * erf 1 (Kilo I /10)) où les notations sont identiques aux notations des équations eq5et eq6. Seuls les éléments de l'ensemble d'éléments E vérifiant l'intégralité des critères de l'ensemble de critères C sont compatibles avec le groupe de référence gréf. L'éventualité 10 d'avoir au moins deux éléments compatibles avec le groupe de référence gréf est peu fréquente. Cependant, si tel est le cas, la distance d est l'écart en temps séparant la dernière impulsion lm du groupe de référence gréf de la première impulsion If de l'élément candidat de l'ensemble d'élément candidat Ec. L'élément réunit avec le groupe de référence gréf est alors l'élément qui minimise la distance d, c'est-à-dire l'élément pour 15 lequel l'écart en temps entre la première impulsion de l'élément et la dernière impulsion I, du groupe de référence gréf est minimum. Lorsqu'un élément de l'ensemble d'éléments candidats E, est annexé au groupe de référence gréf, les données caractéristiques du groupe obtenu sont mises à jour. Ainsi, le temps d'arrivée de la première impulsion du groupe obtenu est le 20 minimum du temps d'arrivée de la première impulsion du groupe de référence gréf et de l'élément annexé. Le temps d'arrivée de la dernière impulsion du groupe obtenu est le maximum des temps d'arrivée de la dernière impulsion du groupe de référence gréf et de l'élément annexé. La direction d'arrivée DOA du groupe obtenu est la moyenne des directions 25 d'arrivée du groupe de référence gréf et de l'élément annexé. La fréquence porteuse du groupe obtenu est donnée par l'équation eq5 suivante : fgroupe - N1*f1+N2*f2 Ni +N2 OÙ fgroupe désigne la fréquence porteuse du groupe obtenu, N1 désigne le nombre d'impulsions du groupe de référence gréf, 30 f1 désigne la fréquence porteuse du groupe de référence gréf, f2 désigne la fréquence porteuse de l'élément annexé, et N2 désigne le nombre d'impulsions de l'élément annexé. La somme des fréquences porteuses du groupe obtenu est l'addition de la somme des fréquences porteuses du groupe de référence gréf et de la somme des fréquences 3031257 22 porteuses de l'élément annexé. La somme des fréquences porteuses au carré du groupe obtenu est l'addition de la somme des fréquences porteuses au carré du groupe de référence gréf et de la somme des fréquences porteuses au carré de l'élément annexé. Le nombre d'impulsions du groupe obtenu est la somme du nombre d'impulsions 5 du groupe de référence gréf et du nombre d'impulsions de l'élément annexé. La largeur d'impulsion LI du groupe obtenu est donnée par l'équation eq9 suivante : Llg''pe N1 + N2 OÙ Llgroupe désigne la longueur d'impulsions LI du groupe obtenu, N1 désigne le nombre d'impulsions du groupe de référence gréf, 10 Lli désigne la longueur d'impulsions LI du groupe de référence gréf, LI2 désigne la longueur d'impulsions LI de l'élément annexé, et N2 désigne le nombre d'impulsions de l'élément annexé. L'intervalle de répétition des impulsions PRI du groupe obtenu est calculé au moyen d'une régression linéaire par critère des moindres carrés sur l'ensemble des 15 impulsions du nouveau groupe. Cette régression linéaire est donnée par l'équation echo suivante : N1 * -I- N2 *LI2 PRI = Ek7 OEki)2 (kiti - t. Ereik )) avec k, l'ordre des impulsions, t1 le temps d'arrivée des impulsions et n le nombre d'impulsions du groupe obtenu. En variante, une approximation de l'intervalle de répétition des impulsions PRI est 20 calculée en supposant les estimées de l'intervalle de répétition des impulsions PRI du groupe de référence gréf et de l'intervalle de répétition des impulsions PRI de l'élément annexé indépendantes entre elles. L'approximation consiste, en premier lieu, à corriger les intervalles de répétition des impulsions PRI des éléments du groupe et les variances d'intervalle de répétition des impulsions PRI des éléments du groupe lorsque l'intervalle 25 de répétition des impulsions PRI de l'un des éléments du groupe est un multiple de l'intervalle de répétition des impulsions PRI de l'autre élément du groupe. Par exemple, si l'intervalle de répétition des impulsions PRI du groupe de référence gréf est égal à deux fois l'intervalle de répétition des impulsions PRI de l'élément annexé, l'intervalle de répétition des impulsions PRI du groupe de référence gréf est divisée par deux et la 30 variance d'intervalle de répétition des impulsions PRI du groupe de référence gréf est 3031257 23 divisée par quatre. Ensuite, l'intervalle de répétition des impulsions PRI du groupe obtenu est calculé à partir de l'équation echi suivante : var (PRI2) * PRI1 + var (PRIi) * PRI2 PRI = var (PRIi) + var(PRI2) où PRI est l'estimée de l'intervalle de répétition des impulsions du groupe obtenu, PR/1 est l'estimée de l'intervalle de répétition des impulsions du groupe de 5 référence g réf, var (PRIi) est la variance de PRII, PRI2 est l'estimée de l'intervalle de répétition des impulsions de l'élément annexé, et var (PRI2) est la variance de PR/2. 10 En outre, la variance de l'estimée de l'intervalle de répétition des impulsions PRI obtenue avec l'équation eqiivérifie l'équation eq12 suivante : var (PRIi) * var (PRI2) var(PRI) = var (PRIi) + var(PRI2) où var(PRI) est la variance de l'estimée de l'intervalle de répétition des impulsions du groupe obtenu et les autres notations sont identiques à celles de l'équation eqii. En variante, dans le cas où les variances des intervalles de répétition des 15 impulsions PRI n'ont pas été estimées pour chaque élément du groupe obtenu, un autre calcul d'intervalle de répétition des impulsions PRI consiste à corriger les intervalles de répétition des impulsions multiples comme dans la variante précédente et à calculer l'estimée de l'intervalle de répétition des impulsions PRI du groupe obtenu en utilisant l'équation eq13 suivante : PRI = (T0Afin - TOAdeb + 2) 20 où TOAdeb est le temps d'arrivée de la première impulsion du groupe obtenu, TOAfin est le temps d'arrivée de la dernière impulsion du groupe obtenu, et les autres notations sont identiques aux notations de l'équation eqii. La variance de l'estimée de l'intervalle de répétition des impulsions PRI obtenue avec l'équation equvérifie l'équation eq14 suivante : 2 (TOAfin - TOAdeb1)2 PRI + 2) 25 où TQA est la variance de la mesure du temps d'arrivée des impulsions, et les autres notations sont identiques aux notations de l'équation eq13. TOAfin - TOAdeb var(PRI) = 2 3031257 24 La phase en temps d'arrivée TOA des impulsions du groupe obtenu est calculée suite aux étapes suivantes. En premier lieu, un nouveau temps dè référence tref est pris centré entre le premier élément et le deuxième élément formant le groupe. Puis, les phases de chacun des éléments formant le groupe sont calculées avec ce nouveau temps de référence tref. Ensuite, la phase du groupe obtenu par rapport à ce temps de référence tref vérifie l'équation eq16 suivante : ePTOA = + N2 où (pj,0A est la phase en temps d'arrivée des impulsions du groupe obtenu calculée par rapport à un temps de référence tref, (p-0,4,1 est la phase du groupe de référence gréf recalculée par rapport à un temps 10 de référence tref, N1 est le nombre d'impulsions du groupe de référence gréf, et N2 est le nombre d'impulsions de l'élément annexé. Enfin, la phase du groupe obtenu calculée par rapport à un temps d'origine vérifie l'équation eq16 suivante : (Proa = ((P.0,4 + tref) * mod (Fru) 15 où (proA est la phase en temps d'arrivée des impulsions du groupe obtenu calculée par rapport au temps d'origine, (x)mod(y) désigne le modulo comme dans l'équation eq6, tref est le temps de référence par rapport auquel la phase cproA a été calculée, er est l'estimée de l'intervalle de répétition des impulsions du groupe obtenu, et 20 les autres notations sont identiques aux notations de l'équation eq16, Le temps d'origine est le temps de réception de la première impulsion sur le récepteur 12 ou le plus petit temps des impulsions enregistrées sur l'interface 32. La loi de regroupement variable appliquée lors de la deuxième sous-étape 170 est décrite de manière spécifique dans ce qui suit. 25 Chaque élément de l'algorithme est un train d'impulsions T1, ..., Tri de la deuxième classe C2 des trains d'impulsions T1, ..., Tr, de fréquence porteuse variable et chaque groupe de l'algorithme est un palier d'impulsions P1, Pp de fréquence porteuse variable. En outre, l'ensemble d'éléments E de l'algorithme est l'ensemble des trains d'impulsions T1, Tn de la deuxième classe C2 des trains d'impulsions T1, Tn de 30 fréquence porteuse variable. Dans ce qui suit, seules les différences de la deuxième sous-étape 170 par rapport à la première sous-étape 160 sont mises en évidence. A[1* (1/7'0A,1+ N2* (P;'0A,1 3031257 25 Le critère fréquentiel stipule que la différence entre la fréquence porteuse maximale des impulsions d'un élément et la fréquence porteuse minimale des impulsions de l'élément doit être inférieure à un seuil. Ce seuil est lié aux limites technologiques de balayage fréquentiel des émetteurs radar d'intérêt. Ce seuil est choisi d'après les 5 statistiques sur les formes d'ondes de la base 'de données utilisée. Ce seuil est, par exemple, compris au sens large entre 100 Mégahertz (MHz) et 10 Gigahertz (GHz). Le procédé de désentrelacement comprend, ensuite une étape 180 d'association en recouvrement des paliers d'impulsions P1, Pp selon une loi d'association. Lors de l'étape d'association en recouvrement 180, les paliers d'impulsions 10 P1, Pp ayant un même intervalle de répétition des impulsions PRI et qui se superposent dans le temps sont regroupés pour former des groupes de paliers d'impulsions G1,..., G9 en recouvrement. Plus précisément, l'étape 180 d'association en recouvrement des paliers d'impulsions Pl, Pp est mise en oeuvre, uniquement pour les paliers d'impulsions 15 P1, ..., Pp de fréquence porteuse fixe obtenus lors de l'étape de regroupement 150. L'étape d'association en recouvrement 180 permet, notamment, de regrouper des formes d'ondes modulée en fréquence appelée FMICW (abréviation de l'anglais Frequency Modulated Interrupted Continuous Wave). Les ondes FMICW sont formées d'impulsions espacées d'un même intervalle de répétition des impulsions PRI et dont la 20 fréquence porteuse augmente continument dans le temps pendant une durée donnée. Le profil fréquentiel d'une onde FMICW est illustré en fonction du temps sur la figure 10. Les impulsions de même fréquence et ayant la même période de répétition des impulsions PRI forment un train d'impulsions ou un palier d'impulsions. Comme visible sur la figure 10, chaque impulsion appartient à un palier 25 d'impulsions P3. Les impulsions d'un même palier d'impulsions P3 ont la même fréquence. Les paliers, d'impulsions P3 de la figure 10 n'ont pas été regroupés pour former un seul palier d'impulsions lors de l'étape de regroupement 150 car le critère temporel de l'étape de regroupement 150 élimine les éléments se superposant dans le temps. Or, les paliers d'impulsions ou trains d'impulsions d'une onde FMICW se superposent dans le temps. 30 L'étape d'association en recouvrement 180 permet donc, dans le cas des ondes FMICW, d'assembler des paliers d'impulsions se superposant dans le temps de sorte à reconstruire un signal FMICW. Ainsi, sur cette figure 10, l'assemblage des paliers d'impulsions P3 permet de former un groupe de paliers d'impulsions G1. 3031257 26 La loi d'association utilisée pour associer en recouvrement les paliers d'impulsions P1, Pp de fréquence porteuse fixe est par la suite appelée « loi d'association en recouvrement ». La loi d'association en recouvrement est mise en oeuvre par un algorithme de 5 même structure générale que la structure générale de l'algorithme de la loi de recouvrement fixe décrit dans l'étape de regroupement 150. Chaque élément de l'algorithme est un palier d'impulsions P1, Pp de fréquence porteuse fixe et chaque groupe de l'algorithme est un groupe de paliers d'impulsions G1, ..., Gg en recouvrement. Le groupe de référence gréf de l'algorithme est un groupe de 10 paliers d'impulsions G1,..., Gg. En outre, l'ensemble d'éléments E de l'algorithme est l'ensemble des paliers d'impulsions P1, Pp de fréquence porteuse fixe. Dans ce qui suit, les caractéristiques de la loi d'association en recouvrement identiques aux caractéristiques de la loi de regroupement fixe ne sont pas décrites à nouveau. Seules les différences de la loi d'association en recouvrement par rapport à la 15 loi de regroupement fixe sont décrites dans ce qui suit. Le critère temporel stipule que pour être compatibles, deux éléments doivent se superposer dans le temps. Lorsque le groupe de référence gréf est un groupe de paliers d'impulsions, le critère fréquentiel est identique au critère fréquentiel de la deuxième sous-étape 170 de l'étape 20 de regroupement 150. Dans ce cas, le critère de phase n'est pas appliqué. En effet, il n'existe aucune contrainte sur les phases en temps d'arrivée des impulsions au niveau de groupes de paliers. Lorsque le groupe de référence gréf est un palier d'impulsions P1, Pp, le critère de phase est identique au critère de phase de la première sous-étape 160 de l'étape de 25 regroupement 150. Dans ce cas, si le critère de phase est vérifié et si les critères de direction d'arrivée, temporel, de largeur d'impulsion et d'intervalle de répétition des impulsions sont également vérifiés, alors le critère fréquentiel s'évalue selon deux cas de figures. Dans le premier cas de figure, les fréquences porteuses des impulsions des deux 30 paliers d'impulsions P1, Pp à comparer sont identiques et le nombre théorique d'impulsions des deux paliers réunis est supérieur ou égal au nombre réel d'impulsions des deux paliers d'impulsions P1, Pp réunis. Dans ce premier cas de figure, les deux paliers sont compatibles. Dans le deuxième cas de figure, les fréquences porteuses des impulsions des 35 deux paliers d'impulsions P1, Pp à comparer sont différentes ou le nombre théorique 3031257 27 d'impulsions des deux paliers réunis est inférieur au nombre réel d'impulsions des deux paliers d'impulsions P1, Pp réunis. Dans ce deuxième cas de figure, les deux paliers sont incompatibles. Entre plusieurs éléments compatibles pour être réunis avec un groupe de 5 référence gréf, l'étape d'association en recouvrement 180 donne la priorité à l'élément dont la phase est déjà présente dans le groupe de référence gréf. Cette première règle de priorité permet de réunir en premier les paliers d'impulsions P1, Pp éclatés et de même intervalle de répétition des impulsions PRI. En outre, si un groupe correspondant à un FMICW est déjà en formation, une deuxième règle de priorité s'applique. Selon cette 10 deuxième règle de priorité, les informations des paliers déjà présents dans le groupe correspondant à un FMICW sont consolidées, avant de rajouter un nouveau palier d'impulsions P1, ..., Pp au groupe. Une fois les règles de priorité appliquées, la préférence est donnée à l'élément compatible qui a le meilleur taux de recouvrement avec le groupe de référence gref. Le 15 taux de recouvrement est donné par l'équation eqi, suivante : ATrecouvement tauxrecouvrement Min(ATi, AT2) où tauxrecouvrement est le taux de recouvrement entre un premier élément et un deuxième élément, ce taux de recouvrement est compris entre 0 et 100% LTrecouvement est l'écart temporel entre le temps de réception de la première impulsion de l'élément compatible et le temps de réception de la dernière impulsion du 20 groupe de référence gréf, LiTi est la durée du groupe de référence gréf, AT2 est la durée de l'élément compatible, et Min(AT1,AT2) est le minimum de ATiet AT2. La distance d que doit minimiser un élément de l'ensemble d'éléments candidats 25 pour être réuni au groupe de référence gréf est donnée par l'équation eq18 suivante : d = a + (1 - taux recouvrement) où a prend la valeur zéro si les phases sont égales, et prend la valeur 1 si les phases ne sont pas égales, les autres notations étant identiques aux notations de l'équation eq17. Le choix d'une telle distance d permet de réunir au groupe de référence gréf, 30 l'élément qui fait croître le plus rapidement un taux d'explication. Le taux d'explication d'un groupe est le rapport du nombre d'impulsions réelles du groupe sur le nombre d'impulsions théoriques du groupe. 3031257 28 Lorsqu'un élément est annexé au groupe de référence gréf, les caractéristiques du groupe obtenu en terme de direction d'arrivée DOA des impulsions, de temps d'arrivée de la première impulsion, de temps d'arrivée de la dernière impulsion, de fréquence des impulsions, de nombre d'impulsion, d'intervalle de répétition des impulsions PRI, de 5 largeur d'impulsion LI, de phase en temps d'arrivée des impulsions sont mises à jour en procédant de la même manière que lors de la première sous-étape 160 de l'étape de regroupement 150. Le procédé de désentrelacement comprend, ensuite, une étape 190 d'association en commutation des paliers d'impulsion. 10 Lors de l'étape d'association en commutation 190, les paliers d'impulsions P1, Pp ayant un intervalle de répétition des impulsions PRI différent et s'enchaînant dans le temps sont regroupés suivant au moins une loi d'association prédéfinie pour former des groupes de paliers d'impulsions G1, ..., Gg en commutation. Par exemple, comme illustré sur la figure 12, les paliers d'impulsions P7, P3, Pg, et 15 P10 sont reliés par des flèches en traits fins pour former un groupe de paliers d'impulsions G1, ..., Gg en commutation. De même, les paliers d'impulsions 11 P tP - , P 12, P 13 et P14 sont reliés par des flèches en traits gras pour former un autre groupe de paliers d'impulsions G1,..., Gg en commutation. Comme visible sur cette figure 12, chaque groupe de paliers G1, ... Gg est formé de paliers qui ne se superposent pas dans le temps et qui ont des 20 intervalles de répétition des impulsions PRI différents. Plus précisément, l'étape 190 d'association en commutation des paliers d'impulsions P1, Pp est mise en oeuvre, d'une part, lors d'une première sous-étape 200, pour les paliers d'impulsions P1, Pp de fréquence porteuse fixe n'ayant pas été regroupés lors de l'étape d'association en recouvrement 180. A l'issue de cette première 25 sous-étape 200, des groupes de paliers d'impulsions G1, ..., Gg en recouvrement de fréquence porteuse fixe sont obtenus. Dans ce cas, la loi d'association mise en oeuvre est appelée « loi d'association en commutation fixe ». L'étape 190 d'association en commutation des paliers d'impulsions P1, Pp est mise en oeuvre, d'autre part, pour les paliers d'impulsions P1, Pp de fréquence 30 porteuse variable lors d'une deuxième sous-étape 210 pour obtenir des groupes de paliers d'impulsions G1, ..., Gg en commutation de fréquence porteuse, variable. Dans ce cas, la loi d'association mise en oeuvre est appelée « loi d'association en commutation variable ». 3031257 29 Chacune des lois d'association en commutation est mise en oeuvre par un algorithme dont la structure générale est identique à la structure générale de l'algorithme décrit lors de l'étape de regroupement 150. Chaque élément de l'algorithme est un palier d'impulsions P1, Pp et chaque 5 groupe de l'algorithme est un groupe de paliers d'impulsions G1, ..., G9. Le groupe de référence gréf est un groupe de paliers d'impulsions G1,..., G9. Dans ce qui suit, les caractéristiques de la loi d'association en commutation fixe identiques aux caractéristiques de la loi de regroupement fixe ne sont pas décrites à nouveau. Seules les différences de la loi d'association en commutation fixe par rapport à 10 la loi de regroupement fixe sont décrites dans ce qui suit. Chaque élément de l'algorithme est un palier d'impulsions P1, Pp de fréquence porteuse variable n'ayant pas été réuni lors de l'étape d'association en recouvrement 180. Chaque groupe de l'algorithme est un groupe de paliers d'impulsions G1,..., G9 en commutation de fréquence porteuse fixe. En outre, l'ensemble d'éléments E de 15 l'algorithme est l'ensemble des paliers d'impulsions P1, Pp de fréquence porteuse fixe n'ayant pas été réunis lors de l'étape de regroupement 150. Le critère de largeur d'impulsion permet de déterminer les paliers ayant des largeurs d'impulsions LI proches ou des facteurs de forme proches. Le facteur de forme d'un palier est défini comme le rapport entre la largeur d'impulsion LI du palier et 20 l'intervalle de répétition des impulsions PRI. En effet, l'observation de formes d'ondes spécifiques, notamment des formes d'ondes de type haute fréquence de récurrence et moyenne fréquence de récurrence, montre une constante dans les valeurs des largeurs d'impulsions d'un même signal et/ou dans les valeurs des facteurs de forme des paliers d'un même signal. Les mesures de largeurs d'impulsions LI étant relativement peu fiables, 25 le critère de largeurs d'impulsion consiste en un test de comparaison des ordres de grandeur des largeurs d'impulsion de deux éléments ou des facteurs de forme de deux éléments. Le critère d'intervalle de répétition des impulsions PRI permet de déterminer les paliers ayant des intervalles de répétition des impulsions proches. En effet, l'observation 30 de différentes formes d'ondes permet de définir une limite supérieure dans les écarts d'intervalles de répétition des impulsions des paliers d'un même signal. En effet, pour un même signal, les intervalles de répétition des impulsions des différents paliers formant le signal ne doivent pas être espacés de plus d'une certaine valeur seuil. Par conséquent, deux paliers sont compatibles lorsque les intervalles de répétition des impulsions PRI de 35 chacun des deux paliers vérifient l'équation eq19suivante : 3031257 30 Max(PRI1,PRI2)-Min(PRI1,PRI2) REVENDICATIONS1.- Procédé de désentrelacement de signaux radars, le procédé comprenant : - la réception (100) de signaux électromagnétiques par un récepteur (12) et l'extraction des impulsions (li, Im) des signaux reçus, - la formation (110) de trains d'impulsions (T1, Tn) regroupant au moins trois impulsions espacées d'un même intervalle de répétition des impulsions (PRI), chaque train d'impulsions (T1, ..., Te) étant défini par l'intervalle de répétition des impulsions (PRI), caractérisé en ce que le procédé comporte en outre : - le regroupement (150) des trains d'impulsions (T1, Tn) ayant un même intervalle de répétitions des impulsions (PRI) suivant une loi de regroupement prédéfinie pour former des paliers d'impulsions (P1, Pp), et - l'association (180, 190) des paliers d'impulsions (P1, Pp) suivant au moins une loi d'association prédéfinie pour obtenir des signaux radars désentrelacés formés de la concaténation des trains d'impulsions (T1, ..., Tn) des paliers d'impulsions (P1, ..., Pp) associés. 2.- Procédé selon la revendication 1, dans lequel chaque train d'impulsions (T1, Tn) est également défini par au moins un élément choisi dans un groupe constitué de : le temps d'arrivée (TOA) de la première impulsion (11) du train d'impulsions (T1, ..., Ta), le temps d'arrivée (TOA) de la dernière impulsion (lm) du train d'impulsions (T1, ..., Te), la fréquence des impulsions (I1, ..., Im) du train d'impulsions (T1, ..., Te), la durée des impulsions (Il , lm) du train d'impulsions (T1, Tn) et la direction d'arrivée des impulsions (Il, ..., lm) du train d'impulsions (T1, .- Ta). 3.- Procédé selon la revendication 1 ou 2, dans lequel le procédé comporte avant l'étape de regroupement (150), une étape de classification (140) des trains d'impulsions (T1, ..., Te) selon leur fréquence porteuse pour obtenir deux classes (C1, C2) de trains d'impulsions (T1, Tn) : une classe (C1) regroupant les trains d'impulsions (T1, ..., Te) de fréquence porteuse fixe et l'autre classe (C2) regroupant les trains d'impulsions (T1, ..., Te) de fréquence porteuse variable, l'étape de regroupement (150) étant mise en oeuvre pour chacune des deux classes (C1, C2) de trains d'impulsions (T1, ..., Tn) et permettant d'obtenir des paliers d'impulsions (P1, ..., Pp) mono-fréquence à partir de la classe (C1) des trains d'impulsions (T1, Tn) mono- 3031257 36 fréquence et des paliers d'impulsions (P1, Pp) agiles en fréquence à partir de la classe (C2) des trains d'impulsions (T1, ..., Tr,) agiles en fréquence. 4.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape d'association (180, 190) comprend une phase de regroupement des paliers d'impulsions (P1, Pp) d'intervalles de répétition des impulsions (PRI) différents et qui s'enchaînent dans le temps pour obtenir des groupes de paliers d'impulsions (G1,..., G9) en commutation. 10 5.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape d'association (180, 190) comprend une phase de regroupement des paliers d'impulsions (P1, Pp) de fréquence porteuse fixe, ayant des intervalles de répétitions des impulsions (PRI) identiques et se superposant dans le temps pour obtenir des groupes de paliers d'impulsions (G1,..., G9) en recouvrement. 15 6.- Procédé selon l'une quelconque des revendications précédentes, dans lequel chacune des lois de regroupement et d'association est mise en oeuvre par au moins un algorithme permettant d'obtenir des groupes à partir d'éléments, les éléments désignant des trains d'impulsions (T1, Tn) lors de l'étape de regroupement et des paliers 20 d'impulsions (P1, Pp) lors de l'étape d'association, les groupes désignant des palier d'impulsions (P1, Pp) lors de l'étape de regroupement et des groupes de paliers d'impulsions (G1,..., G9) lors de l'étape d'association, l'algorithme comprenant : o le choix (150A) d'un élément de référence (eréf) parmi un ensemble d'éléments (E), 25 o la suppression (150B) de l'élément de référence (eréf) de l'ensemble d'éléments (E) et l'ajout, dans un ensemble de groupes (E9), d'un groupe de référence (gréf) comportant l'élément de référence (eréf), o la sélection (150C) dans l'ensemble d'éléments (E), des éléments compatibles avec le groupe de référence (gréf) selon un ensemble de 30 critères (C) pour obtenir un ensemble d'éléments candidats (Es), o l'évaluation (150D) de la distance entre le groupe de référence (gréf) et chaque élément de l'ensemble d'éléments candidats (Es), o l'annexion (150E) de l'élément de l'ensemble d'éléments candidats (Es) minimisant une distance (d) au groupe de référence (gréf) et la suppression 35 de l'élément annexé de l'ensemble d'éléments (E), 3031257 37 la répétition des phases de sélection (150C), d'évaluation (150D) et d'annexion (150E) tant que l'ensemble des éléments candidats (Es) comporte des éléments, et la répétition de l'ensemble des phases précédentes (150A, 150B, 150C, 5 150D, 150E) tant que l'ensemble d'éléments (E) comporte des éléments. 7.- Procédé selon la revendication 6, dans lequel l'élément de référence (eréf) est l'élément de l'ensemble d'éléments (E) dont le temps d'arrivée (TOA) de la première impulsion (11) est le plus petit. 10 8.- Procédé selon l'une quelconque des revendication 6 ou 7, dans lequel l'ensemble de critères (C) évalue la compatibilité des éléments de l'ensemble d'éléments (E) avec le groupe de référence (gréf) en fonction d'une ou de plusieurs caractéristiques, les caractéristiques étant choisies dans un groupe comprenant : la direction d'arrivée 15 (DOA) des éléments, la superposition temporelle des éléments, la fréquence porteuse (f) des éléments, la largeur d'impulsion (LI) des éléments, l'intervalle de répétition des impulsions (PRI) des éléments, la phase (1) des éléments et le nombre d'impulsions des éléments. 20 9.- Procédé selon l'une quelconque des revendications 6 à 8, dans lequel les critères (C) sont choisis d'après des statistiques sur les caractéristiques des formes d'ondes radars d'une base de données. 10.- Procédé selon l'une quelconque des revendications 6 à 9, dans lequel pour la 25 loi de regroupement, la distance (d) est l'écart en temps séparant la dernière impulsion (lu) du groupe de référence (gréf) de la première impulsion (If) de l'élément candidat de l'ensemble d'éléments candidats (Es), et dans lequel pour la loi d'association, la distance (d) est un taux de recouvrement entre le groupe de référence (gréf) et l'élément candidat de l'ensemble d'éléments candidats (Ec) ou un score permettant de sélectionner parmi 30 l'ensemble d'éléments candidats (Ec) l'élément partageant le plus de caractéristiques communes avec le groupe de référence (gréf). 11.- Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comporte avant l'étape de regroupement (150), une étape de rejet (120) des 35 trains d'impulsions incohérents en termes d'intervalle de répétition des impulsions (PRI). 3031257 38 12.- Procédé selon l'une quelconque des revendications précédentes, dans lequel le procédé comporte avant l'étape de regroupement (150), une étape de rejet (130) des trains d'impulsions dont l'intervalle de répétition des impulsions (PRI) est supérieur à un 5 seuil d'intervalle de répétition des impulsions (SARI) et dont le nombre d'impulsions (Ii, ..., lm) est inférieur à un seuil de nombre d'impulsions (Simpuision)- 13.- Dispositif de désentrelacement de signaux radars comprenant : - un récepteur (12) propre à recevoir des signaux électromagnétiques, 10 - une unité de traitement numérique du signal (23) propres à extraire les impulsions (li, 1,) des signaux reçus par le récepteur (12), et - un support lisible d'informations (26), sur lequel est mémorisé un programme d'ordinateur comprenant des instructions de programme, le programme d'ordinateur étant chargeable sur une unité de traitement de données (24) et adapté pour entraîner la mise 15 en oeuvre d'un procédé selon l'une quelconque des revendications précédente lorsque le programme d'ordinateur est mis en oeuvre sur l'unité de traitement des données (24).