WO2015044248A1

WO2015044248A1 - Method for non-supervised deinterleaving by n-dimensional enrichment

Info

Publication number: WO2015044248A1
Application number: PCT/EP2014/070450
Authority: WO
Inventors: Jean-François Grandin; Jean-Marie Lemoine
Original assignee: Thales
Priority date: 2013-09-27
Filing date: 2014-09-25
Publication date: 2015-04-02
Also published as: US20160209492A1; EP3049822A1; US10509102B2; FR3011338A1; EP3049822B1

Abstract

The invention relates to the field of listening systems, especially to a method for the non-supervised deinterleaving of pulse trains, comprising at least one step of N-dimensional enrichment, N being a whole number higher than 1.

Description

PROCÉDÉ DE DÉS ENTRELACEMENT NON S UPERVISÉ PAR ENRICHIS S EMENT N DIMENS IONS METHOD OF INTERLEAVING NON-UPERVISED BY N DIMENS ION ENRICHMENT

La présente invention concerne le domaine des systèmes d'écoute, et plus particulièrement le domaine du traitement de signal exploité par un tel système. La présente invention concerne plus particulièrement un procédé de désentrelacement non supervisé par enrichissement N dimensions. The present invention relates to the field of listening systems, and more particularly the field of signal processing operated by such a system. The present invention more particularly relates to a method of unsupervised deinterlacing by N-dimensional enrichment.

Le traitement considéré est celui du désentrelacement non supervisé des trains d'impulsions émises par des émetteurs électromagnétiques. En référence à la figure 1 , lorsqu'un capteur de Mesure de Renseignement Electronique (ou MRE), également connu sous le sigle anglo-saxon ESM ou ELINT pour respectivement Electronic Support Measures et ELectronic INTelligence, intercepte ces trains d'impulsions, il traite toutes les impulsions dont l'intensité est supérieure à son seuil de sensibilité. Les émissions étant simultanément présentes, ces trains d'impulsions sont imbriquées ou " entrelacés The treatment considered is that of the unsupervised deinterleaving of pulse trains emitted by electromagnetic emitters. With reference to FIG. 1, when an Electronic Intelligence Measurement (or MRE) sensor, also known by the abbreviation ESM or ELINT for Electronic Support Measures and ELectronic INTelligence respectively, intercepts these pulse trains, it processes all pulses whose intensity is greater than its sensitivity threshold. Emissions being simultaneously present, these pulse trains are interleaved or "interlaced"

Ces situations de mélange ou entrelacement de trains d'impulsions sont d'autant plus fréquentes que la sensibilité des capteurs MRE augmente induisant une augmentation de la densité des impulsions détectées, et que les émissions agiles étalent leurs impulsions sur le spectre, ce qui augmente la probabilité de trouver des impulsions de plusieurs radars à une fréquence donnée. These situations of mixing or interlacing of pulse trains are all the more frequent as the sensitivity of the MRE sensors increases inducing an increase in the density of the detected pulses, and the agile emissions spread their pulses on the spectrum, which increases the frequency of the pulses detected. probability of finding pulses of several radars at a given frequency.

Pour traiter ces situations l'extracteur doit en tout premier lieu isoler chaque émission d'intérêt pour pouvoir ensuite l'analyser et la caractériser. On parle de " Désentrelacement ". To deal with these situations, the extractor must first isolate each emission of interest in order to analyze and characterize it. We talk about "deinterlacing".

Les émissions radars interceptées sont grossièrement de trois types :Intercepted radar emissions are roughly of three types:

- TE1 : Cela correspond aux émissions de Haute et Moyenne fréquence de Récurrence (HFR/MFR) présentant des volumes de 30 000 à plus de 100 000 impulsions par seconde par émission. Ces émissions sont de caractéristiques simples, c'est-à-dire de Période de Répétition par Intervalle (PRI) et/ou Fréquence stable par train. Elles sont généralement traitées par des méthodes ad hoc, rapides et robustes. - TE2 : Cette classe correspond aux émissions Basses Fréquence de Récurrence (BFR) présentant des volumes de 1 000 impulsions par seconde par émission. Ces émissions peuvent présenter des caractéristiques complexes comme par exemple une période de répétition variable (ou " stagger " selon le terme anglo saxon) d'ordre élevé ou une agilité de fréquence impulsion à impulsion. - TE1: This corresponds to the High and Medium Frequency (HFR / MFR) emissions with volumes from 30,000 to more than 100,000 pulses per second per transmission. These emissions are simple features, ie Interval Repetition Period (IRP) and / or Stable Frequency per train. They are usually processed by ad hoc, fast and robust methods. - TE2: This class corresponds to Low Frequency Recurrence (BFR) transmissions with volumes of 1000 pulses per second per transmission. These emissions may have complex characteristics such as, for example, a high order variable repetition period (or "stagger") or impulse pulse frequency agility.

- TE3 : Ce type d'émission se réfère aux émissions BFR avec agilité impulsion à impulsion de PRI et de fréquence. Ou émissions à faible probabilité d'interception (ou LPI pour Low Probability of intercept selon la terminologie anglo saxonne) à faible nombre d'impulsions et fortement modulées. Diverses méthodes de traitement sont utilisables comme par exemple, la reconnaissance sur les paramètres intra-impulsions et DOA ou le tri par direction d'arrivée (ou DOA pour Direction Of Arrivai) et localisation. - TE3: This type of emission refers to BFR transmissions with pulse agility at pulse of PRI and frequency. Or emissions low probability of interception (or LPI for low probability of intercept according to English terminology) low pulse number and strongly modulated. Various processing methods can be used, such as, for example, recognition on intra-pulse and DOA parameters or direction of arrival (DOA) sorting and location.

Par la suite, on s'intéressera au désentrelacement des émissions de type TE2 et TE3. Les émissions de type TE1 sont traitées séparément par des algorithmes rapides. Thereafter, we will focus on the deinterleaving of TE2 and TE3 type emissions. TE1 emissions are handled separately by fast algorithms.

On distingue deux fonctions de désentrelacement complémentaires:There are two complementary deinterlacing functions:

- Le Désentrelacement Supervisé (DS), dans lequel on utilise les caractéristiques des signaux connues comme par exemple, les Fréquences- Supervised Deinterlacing (DS), in which the characteristics of the known signals are used, for example, the Frequencies

(F), les PRI, les Durées d'Impulsion (Dl) ou les paramètres intra impulsion (ou IP pour IntraPulse) pour reconnaître les impulsions provenant de ces signaux. (F), PRIs, Pulse Durations (Dl) or intra pulse parameters (or IP for IntraPulse) to recognize pulses from these signals.

- Le Désentrelacement Non Supervisé (DNS), dans lequel on n'a pas connaissance de ces caractéristiques. - Non-Supervised Deinterlacing (DNS), in which we are not aware of these characteristics.

Dans la plupart des extracteurs ESM ou ELINT ces deux fonctions coopèrent comme illustré figure 2. In most ESM or ELINT extractors these two functions cooperate as shown in Figure 2.

Les méthodes de DNS peuvent être segmentées entre : DNS methods can be segmented between:

- Celles qui utilisent des classifications (ou clustering) sur paramètres primaires de l'impulsion (DOA, F, Dl, IP) - Those that use classifications (or clustering) on primary parameters of the pulse (DOA, F, Dl, IP)

- Celles qui utilisent des classifications à partir des temps d'arrivée des impulsions (ou TOA pour Time Of Arrivai suivant l'expression anglo saxonne) correspondant au paramètre secondaire PRI. Dans la plupart des extracteurs les deux méthodes sont utilisées successivement comme par exemple selon le schéma figure 3. - Those that use classifications from the arrival times of the pulses (or TOA for Time Of Arrival according to the Anglo-Saxon expression) corresponding to the PRI secondary parameter. In most extractors, the two methods are used successively, for example according to the diagram in FIG.

De plus dans ces extracteurs le désentrelacement et le pistage sont réalisés de proche en proche en n'examinant pas la totalité des informations simultanément, Moreover in these extractors the deinterlacing and the tracking are carried out step by step by not examining all the information simultaneously,

Un grand volume de données nécessite d'être traité par un capteur MRE (ESM ou ELINT), avec des temps de réaction rapide. En général la charge de calcul nécessaire évolue plus rapidement que l'évolution en puissance des processeurs (loi de Moore) du fait de l'augmentation du nombre de données à traiter et de la complexité jusqu'ici non linéaire des traitements. De façon non limitative, les raisons suivantes peuvent expliquer l'augmentation drastique des volumes de données ! A large amount of data needs to be processed by an MRE sensor (ESM or ELINT), with fast reaction times. In general, the computing load required evolves more rapidly than the power evolution of the processors (Moore's law) because of the increase in the number of data to be processed and the hitherto non-linear complexity of the processing. In a non-limiting way, the following reasons may explain the drastic increase in data volumes!

- Les capteurs ELINT ont des sensibilités de plus en plus élevées pour la détection de certaines menaces. - ELINT sensors have increasing sensitivities for the detection of certain threats.

- Les émissions sont agiles ce qui nécessite, pour capter au maximum la forme d'onde : - Emissions are agile which requires, to capture the maximum waveform:

- d'avoir des bandes de réception importantes, - to have important reception bands,

- d'avoir des temps d'écoute longs. - to have long listening times.

L'augmentation de la densité temporelle des impulsions introduit également des effets de superposition d'impulsions simultanées conduisant à perdre une partie des impulsions d'une forme d'onde par mitage ou moirage. The increase in the temporal density of the pulses also introduces simultaneous pulse superposition effects leading to the loss of a portion of the pulses of a waveform by mime or moire.

La canalisation de certains récepteurs impose de balayer la bande de réception utile, ne permettant l'acquisition de la forme d'onde que partiellement. The channeling of certain receivers makes it necessary to scan the useful reception band, allowing the acquisition of the waveform only partially.

On conçoit que si on observe ces formes d'onde sur des temps longs, il va être possible de récupérer des indices statistiques et de les mettre en relation. It is conceivable that if we observe these waveforms over long times, it will be possible to retrieve statistical indices and put them in relation.

En résumé, on constate que : In summary, we find that:

- les volumes de données augmentent ; - the data volumes increase;

- il serait utile d'avoir des analyses sur des temps d'écoute longs ; - it would be useful to have analyzes on long listening times;

- Il faut traiter les paramètres dans une analyse conjointe. Les algorithmes actuels ne permettent pas de remplir ces différentes contraintes. - The parameters must be processed in a joint analysis. Current algorithms do not fulfill these different constraints.

Si la complexité de l'algorithme n'est que quadratique cela signifie en pratique que traiter 1 seconde de durée d'écoute coûte 100 fois plus cher, en temps de calcul, que traiter 100 ms de durée d'écoute. Ce qui devient 10 000 pour 10 secondes. Les algorithmes actuels ont souvent des complexités plus que quadratiques. Ceci interdit l'application de ces algorithmes sur des durées importantes bien que ce soit intéressant pour faire apparaître des discriminants statistiques. Donnons deux exemples d'indices statistiques en temps long non triviaux : certaines valeurs de fréquence (idem PRIs) sont souvent associées sur une même DOA dans une même écoute ; certains radars balaient l'espace régulièrement et présentent une Période de Rotation d'Antenne ou PRA. Ces indices peuvent être utilisés en temps long alors qu'ils ne peuvent pas l'être en temps court. If the complexity of the algorithm is only quadratic, it means in practice that processing 1 second of listening time is 100 times more expensive in computing time than processing 100 ms of listening time. Which becomes 10,000 for 10 seconds. Current algorithms often have more than quadratic complexities. This prohibits the application of these algorithms over long durations although it is interesting to reveal statistical discriminants. Let's give two examples of non-trivial long-term statistical indices: some frequency values (idem PRIs) are often associated on the same DOA in the same listening; some speed cameras scan the space regularly and have an Antenna Rotation Period or PRA. These indices can be used in a long time when they can not be used in a short time.

Dans ce contexte, il existe un besoin pour des implémentations efficaces minimisant la complexité de calcul. L'optimum de cette complexité est la complexité linéaire c'est à dire le fait que le nombre des calculs croît linéairement avec le débit des données entrantes. Avec une complexité linéaire il devient possible d'analyser de très importants volumes de données, donc de travailler sur des temps longs avec une grande densité impulsionnelle. Les principes classiquement utilisés en DNS reposent généralement sur la séparation en passages de lobes, la séparation en Direction d'Arrivée, la séparation par paramètres primaires (Fr, Ll, IP) ou la séparation par histogramme de différence de temps d'arrivée (ou HDTOA pour Histogram of Différence Time Of Arrivai selon la terminologie anglo saxonne). La séparation en passages de lobes est basée sur la décorrélation statistique des dates de lobes pour un environnement de menaces. Pour une densité moyenne de menaces (et ayant procédé à un premier filtrage grossier sur la DOA ou/et la fréquence) les éclairements des différentes menaces apparaissent rarement " entrelacés In this context, there is a need for efficient implementations that minimize computational complexity. The optimum of this complexity is the linear complexity, ie the fact that the number of computations increases linearly with the flow rate of the incoming data. With linear complexity it becomes possible to analyze very large volumes of data, thus working on long times with a high pulse density. The principles conventionally used in DNS are generally based on the separation in lobe passages, the separation in Direction of Arrival, the separation by primary parameters (Fr, L1, IP) or the separation by histogram of difference in arrival time (or HDTOA for Histogram of Difference Time Of Arrival according to the English terminology). The separation in lobe passages is based on the statistical decorrelation of lobe dates for a threat environment. For a medium density of threats (and having carried out a first rough filtering on the DOA and / or the frequency) the illuminations of the various threats appear rarely "intertwined

Cette technique est cependant plutôt adaptée aux détecteurs d'alerte (ou RWR pour Radar Warning Receiver) qui remplissent leur rôle d'alerte avec une sensibilité ne leur permettant pas en général d'intercepter les menaces sur diffus. Cette technique n'est évidemment pas adaptée à un système d'écoute qui doit intercepter les menaces sur diffus (ELINT) ou en environnement dense car dans ce cas le système voit toujours les émissions " entrelacées This technique, however, is rather adapted to the Radar Warning Receivers (RWR) which fulfill their alert role with a sensitivity that does not generally allow them to intercept the threats on diffuse. This technique is obviously not suitable for a listening system that must intercept the threats on diffuse (ELINT) or in dense environment because in this case the system always sees the transmissions "intertwined

Un paramètre fiable et facile à utiliser pour désentrelacer les impulsions est la direction d'arrivée (DOA). C'est une évidence de constater que c'est le seul paramètre qu'un radar moderne ne peut pas moduler ; c'est pourquoi la DOA quand elle est mesurée a un rôle primordial dans le procédé d'extraction. Mais dans de nombreux systèmes la DOA présente une qualité " moyenne " résolution de l'ordre de 10°. Il faut de plus faire attention quand la mesure de DOA est en défaut (mesure sur la polarisation croisée des aériens de réception) ou quand plusieurs émissions sont observées dans le même secteur angulaire. A reliable and easy-to-use parameter for deinterleaving pulses is the direction of arrival (DOA). It is obvious that this is the only parameter that a modern radar can not modulate; that is why the DOA when measured has a primary role in the extraction process. But in many systems the DOA has an "average" quality resolution of about 10 °. In addition, attention must be paid when the DOA measurement is in fault (measurement on the cross polarization of the receiving aerials) or when several emissions are observed in the same angular sector.

Une autre méthode consiste à utiliser les paramètres primaires pour séparer les impulsions. Outre le paramètre DOA, la plupart des algorithmes de désentrelacement " classiques " sont basés sur l'utilisation de tris simples sur paramètres mono-impulsions dont essentiellement fréquence, et durée d'impulsion. Ces extracteurs à base de tris sur paramètres mono impulsions ont été utilisés pendant de nombreuses années et peuvent encore être utilisés dans certaines situations simples. On notera l'existence d'algorithmes de tri plus élaborés que des tris simples en ligne (impulsion par impulsion) basés sur une approche statistique (par exemple découpage des modes d'un histogramme) mais encore réduits à l'utilisation des paramètres primaires. Another method is to use the primary parameters to separate the pulses. In addition to the DOA parameter, most "classical" deinterlacing algorithms are based on the use of simple sorts on single-pulse parameters including essentially frequency, and pulse duration. These extractors based on sorting on mono pulse parameters have been used for many years and can still be used in certain simple situations. Note the existence of more elaborate sorting algorithms than simple online sorting (impulse by pulse) based on a statistical approach (for example, cutting modes of a histogram) but still reduced to the use of primary parameters.

La généralisation de l'agilité impulsion à impulsion et de l'agilité par trains des formes d'onde rend aujourd'hui difficile voire caduque cette démarche. En effet dans un environnement dense il apparaît impossible dans certains secteurs angulaires (même avec un capteur de bonne résolution) de séparer deux émissions agiles inconnues par un simple tri fréquentiel. Soit les valeurs de fourchettes sont larges et on obtient un mélange des deux émissions. Soit les valeurs de fourchettes sont fines et on " pulvérise " chaque émission en une myriade de " groupements (ou clusters) d'impulsions " monofréquence. Pour résoudre l'agilité deux techniques sont apparues conjointement : The generalization of impulse impulse agility and train agility of waveforms makes today difficult or even obsolete this Steps. Indeed in a dense environment it appears impossible in certain angular sectors (even with a good resolution sensor) to separate two unknown agile emissions by simple frequency sorting. Either the range values are wide and you get a mix of both emissions. Either the range values are fine and each program is "sputtered" into a myriad of single-frequency "pulse" clusters. To solve the agility two techniques appeared jointly:

• d'une part l'utilisation du DTOA (pour Différence of Time Of Arrivai). • on the one hand the use of DTOA (for Difference of Time of Arrival).

• d'autre part l'utilisation de méthodes de recombinaison des groupes d'impulsions issus de tris monofréquence (+DOA ) en vue de résoudre le problème posé par les agilités. L'introduction de la différence de temps d'arrivée (ou DTOA) permet dans le cas d'un train d'impulsion à PRI fixe de tracker les impulsions correspondantes; il suffit alors de regarder quelles impulsions sont séparées par des " intervalles de temps égaux • on the other hand the use of recombination methods groups of pulses from tris single frequency (+ DOA) to solve the problem posed by agility. The introduction of the arrival time difference (or DTOA) makes it possible, in the case of a fixed PRI pulse train, to track the corresponding pulses; then just look at which pulses are separated by "equal time intervals

On retrouve les deux étapes précédentes : We find the two previous steps:

· Une étape de détection du groupement d'impulsions (ici un groupement autour d'une PRI) réalisée grâce à un histogramme des DTOA ou une Transformée de Fourier (TF) des TOA. · A step of detecting the grouping of pulses (here a grouping around a PRI) carried out by means of a histogram of the DTOA or a Fourier Transform (TF) of the TOA.

• Une étape de sélection des impulsions correspondant à une PRI, que nous appellerons " marquage " (ou " gating " selon le terme anglo saxon). • A step of selection of the pulses corresponding to a PRI, which we will call "marking" (or "gating" according to the Anglo-Saxon term).

On peut utiliser plusieurs types de transformation des DTOA abusivement tous appelés histogrammes de DTOA (ou HDTOA), pour ne citer que les plus remarquables : We can use several types of transformation of DTOA abusively called histograms of DTOA (or HDTOA), to quote only the most remarkable:

· Le HDTOA d'ordre k et ses variantes possibles (Histogramme complet, Histogramme progressif Alldiff, Histogramme progressif séquentiel) qui est une simple sommation des nombres de fois où un écart de temps est observé. Cet histogramme est celui le plus souvent cité ou utilisé. Il est équivalent à la fonction d'autocorrélation de la séquence des TOAs. Il présente le défaut de faire apparaître toutes les raies de PRIs ainsi que leurs combinaisons linéaires. · The HDTOA of order k and its possible variants (Complete Histogram, Progressive Histogram Alldiff, Sequential Progressive Histogram) which is a simple summation of the number of times a time difference is observed. This histogram is the one most often cited or used. It is equivalent to the autocorrelation function of the sequence of TOAs. It has the defect of showing all the lines of PRIs as well as their linear combinations.

• Le HDTOA dit " compressé " qui somme un argument complexe utilisé pour prendre en compte le caractère périodique des DTOA : Extraction de la Période de Répétition de Motif (ou PRM). Il présente des défauts en présence de mitage, induisant parfois des pertes de la PRM sur des périodes de répétition variable (ou " staggers " selon le terme anglo saxon). • The HDTOA says "compressed" which is a complex argument used to take into account the periodic character of the DTOA: Extraction of the Pattern Repetition Period (or PRM). It has defects in the presence of sprawl, sometimes inducing losses of the PRM on periods of variable repetition (or "staggers" according to the Anglo-Saxon term).

• Le HDTOA bidimensionnel qui étend la recherche des PRIs à la recherche de couples (ou transition) de PRIs significatifs. Il présente l'intérêt d'avoir un seuil d'interactions réduit par rapport aux précédents. Son défaut est le coût de calculs. • Two-dimensional HDTOA that extends the search for PRIs in search of significant PRI (or transition) couples. It has the advantage of having a threshold of interactions reduced compared to previous ones. Its default is the cost of calculations.

• L'histogramme des doublets qui permet de rechercher les PRM avec une grande efficacité. On montre qu'il correspond à une construction partielle de l'histogramme bidimensionnel et qu'il présente le même seuil d'interactions. Son intérêt est de réduire l'information contenue dans les TOAs aux seules raies PRM. • The doublet histogram which makes it possible to search the PRM with great efficiency. It is shown to correspond to a partial construction of the two-dimensional histogram and to have the same threshold of interactions. Its interest is to reduce the information contained in the TOAs to only PRM lines.

• L'Histogramme de Fourier. Les techniques fondées sur la construction d'histogrammes de DTOA ne sont pas robustes aux PRI jittées, surtout lorsqu'il y a un mélange. Le calcul de la transformée de Fourier d'un peigne de N Dirac centrés aux instants d'arrivée TOA des impulsions présente un maximum à la fréquence PRF correspondant à l'inverse de la PRI des impulsions. Le calcul de la TF à toutes les fréquences possibles est toutefois très coûteux. Différents algorithmes de ce type ont été décrits dans la littérature. • The Fourier Histogram. Techniques based on the construction of histograms of DTOA are not robust to PRI jittées, especially when there is a mixture. The calculation of the Fourier transform of a N Dirac comb centered at the arrival times TOA of the pulses has a maximum at the PRF frequency corresponding to the inverse of the PRI of the pulses. Calculating the TF at all possible frequencies, however, is very expensive. Different algorithms of this type have been described in the literature.

Ces histogrammes ont différentes propriétés. Tous ces histogrammes peuvent avoir une grande robustesse au mitage et au blanking. These histograms have different properties. All these histograms can have a great robustness to the mitage and the blanking.

Le marquage (ou " gating " selon le terme anglo saxon) à base de DTOA ayant isolé une PRI sur HDTOA ou une PRF sur périodogramme, on passe à une phase de recherche des impulsions correspondantes. Pour cela on extrait en plusieurs passes successives les trains d'impulsions correspondant à chaque PRM détectée. Le résultat est l'obtention de séquences d'impulsions à cette PRM. Le choix parmi les variantes présentées est étroitement lié à la sévérité des environnements, à la qualité des mesures primaires, et au coût de calcul admissible sur la ressource CPU allouée pour l'application. The marking (or "gating" according to the Anglo-Saxon term) based on DTOA having isolated a PRI on HDTOA or PRF on periodogram, we go to a search phase corresponding pulses. For this purpose, the pulse trains corresponding to each detected PRM are extracted in several successive passes. The result is the obtaining of pulse sequences for this PRM. The choice among the variants presented is closely related to the severity of the environments, the quality of the primary measurements, and the allowable calculation cost on the allocated CPU resource for the application.

L'extraction à base de DTOA peut être étendue à toute structure temporelle " contrastée " de la forme d'onde. On constate de plus que pour les formes d'ondes agiles et non agiles la structure temporelle (les PRIs) est plus discriminante, plus caractéristique et plus stable que la structure fréquentielle. En effet il est plus simple de faire évoluer les fréquences d'émission que les jeux de PRIs utilisés. DTOA-based extraction can be extended to any "contrasted" temporal structure of the waveform. It is further noted that for agile and non-agile waveforms the temporal structure (the PRIs) is more discriminating, more characteristic and more stable than the frequency structure. Indeed, it is easier to change the transmission frequencies than the sets of PRIs used.

Les techniques de désentrelacement basées sur le DTOA ont aussi leurs limites et voient ainsi leur efficacité diminuer avec l'apparition des Formes d'Onde à Agilité de PRI, agilité dans le sens d'une structure temporelle non stable. En effet une période de répétition variable (ou " Stagger ") de PRI est une émission dite agile en PRI mais reste très facile à traiter par un algorithme de désentrelacement basé sur le DTOA car le motif est stable dans le temps. A l'inverse les radars à balayage électronique peuvent émettre un signal qui vu d'un point de l'espace ne présente aucune structure temporelle apparente. Dans ce cas le désentrelacement ne peut se fonder sur le DTOA. De-interlacing techniques based on DTOA also have their limits and thus see their effectiveness decrease with the appearance of PRI Agility Wave Forms, agility in the sense of a non-stable temporal structure. Indeed a period of variable repetition (or "Stagger") of PRI is a so-called agile transmission in PRI but remains very easy to process by a deinterlacing algorithm based on the DTOA because the pattern is stable over time. On the contrary, the electronic scanning radars can emit a signal which from a point of space has no apparent temporal structure. In this case the deinterlacing can not be based on the DTOA.

Bien que l'environnement électromagnétique soit composé d'une grande quantité de formes d'ondes à motif de PRI stables, la prise en compte conjointe des paramètres de tris classiques sera très utile pour résoudre le cas des signaux peu nombreux à PRI instable. En outre les techniques à base de DTOA peuvent être " piégées " par le synchronisme fortuit d'émissions à PRI identiques et perçues simultanément par le capteur. Dans ce cas l'information de séparation doit être apportée par d'autres paramètres (DOA, Niveau, Fréquence, Di, lp...). Although the electromagnetic environment is composed of a large quantity of stable PRI pattern waveforms, the joint consideration of conventional sorting parameters will be very useful in resolving the case of unstable PRI low signal numbers. In addition, the DTOA-based techniques can be "trapped" by the coincidental synchronism of identical PRI transmissions simultaneously perceived by the sensor. In this case the separation information must be provided by other parameters (DOA, Level, Frequency, Di, lp ...).

Un but de l'invention est notamment de corriger tout ou partie des inconvénients précités en proposant une méthode de séparation des trains d'impulsions radar imbriqués fiable, rapide et nécessitant peu de ressource. An object of the invention is in particular to correct all or part of the aforementioned drawbacks by proposing a method of separation of nested radar pulse trains reliable, fast and requiring little resource.

A cet effet, l'invention a pour objet un procédé de désentrelacement non supervisé de trains d'impulsions comprenant : - une étape Etp10 de constitution, pour chaque paramètre ou groupe de paramètres, de N histogrammes, N étant un entier supérieur à 1 , à une ou plusieurs dimensions définissant des groupes d'impulsions par paramètre ou groupe de paramètres, lesdits histogrammes étant représentatifs des trains d'impulsions; For this purpose, the subject of the invention is a method of unsupervised deinterleaving of pulse trains comprising: a step Etp10 constituting, for each parameter or group of parameters, N histograms, N being an integer greater than 1, with one or more dimensions defining groups of pulses per parameter or group of parameters, said histograms being representative of pulse trains;

- une étape Etp20 d'extraction de modes à partir desdits histogrammes ; an etp20 step for extracting modes from said histograms;

- une étape Etp30 de constitution de groupes d'intérêt à partir desdits modes extraits ; a step Etp30 for forming groups of interest from said extracted modes;

- une étape Etp40 d'enrichissement de la description des impulsions à l'aide de labels traduisant leur appartenance à au moins un groupe d'intérêt ; an Etp40 step of enriching the description of the pulses by means of labels reflecting their membership in at least one interest group;

- une étape Etp50 de classement croisé non supervisé des modes afin de réorganiser les impulsions enrichies par similarité. an etp50 stage of unsupervised cross-classification of the modes in order to rearrange the enriched impulses by similarity.

Selon un mode de mise en œuvre, l'étape Etp20 comprend : According to one embodiment, the Etp20 step comprises:

- une étape Etp21 de seuillage des histogrammes ; an etp21 stage of thresholding of the histograms;

- une étape Etp22 de décomposition des histogrammes en modes ; an etp22 step of decomposing the histograms into modes;

- une étape Etp23 de caractérisation des modes ; a step Etp23 characterization modes;

Selon un mode de mise en œuvre, l'étape Etp50 de classement des modes consiste à créer progressivement au moins une classe d'impulsions- modes homogènes, l'étape comprenant : According to one embodiment, the mode classification step Etp50 consists in progressively creating at least one class of homogeneous mode pulses, the step comprising:

- à l'initialisation une recherche du mode ayant un maximum d'intérêt, ce mode formant la première classe, puis, pour chaque mode restant : at the initialization, a search of the mode having a maximum of interest, this mode forming the first class, then, for each remaining mode:

- une étape Etp51 de recherche, parmi les modes restants, du mode ayant un cardinal d'intersection avec la classe en cours le plus grand ; a step Etp51 of searching, among the remaining modes, of the mode having a cardinal intersection with the largest current class;

- une étape Etp52 de recherche de la partition du mode trouvé, permettant de maximiser un critère de regroupement avec la ou les classes en cours ; an Etp52 step of searching for the found mode partition, making it possible to maximize a grouping criterion with the current class or classes;

- une étape Etp53 de fusion du mode avec la ou les classes en cours. an etp53 step of merging the mode with the current class or classes.

Selon un mode de mise en œuvre, l'étape Etp50 de classement comprend : - une étape Etp56 de classification des impulsions en fonction de leurs liens de corrélation mutuelles ; According to one embodiment, the classification step Etp50 comprises: a step Etp56 of classification of the pulses according to their mutual correlation links;

- une étape Etp57 de classification des modes en fonction de leurs liens de corrélation mutuelles ; a mode classification step Etp57 as a function of their mutual correlation links;

- une étape Etp58 de regroupement des modes et impulsions en classes. a step Etp58 of grouping modes and pulses into classes.

De façon avantageuse, le procédé selon l'invention, grâce à sa faible complexité, permet d'analyser les signaux reçus sur des temps longs avec une grande densité impulsionnelle. Il permet également une analyse conjointe des différents paramètres ou indices statistiques. Advantageously, the method according to the invention, thanks to its low complexity, makes it possible to analyze the signals received over long times with a high pulse density. It also allows a joint analysis of different parameters or statistical indices.

D'autres particularités et avantages de la présente invention apparaîtront plus clairement à la lecture de la description ci-après, donnée à titre illustratif et non limitatif, et faite en référence aux dessins et annexes, dans lesquels : Other features and advantages of the present invention will appear more clearly on reading the following description, given by way of illustration and not limitation, and with reference to the drawings and annexes, in which:

La figure 1 illustre le désentrelacement ; Figure 1 illustrates the deinterlacing;

La figure 2 illustre la coopération entre les désentrelacements supervisé et non supervisé ; Figure 2 illustrates the cooperation between supervised and unsupervised deinterlacing;

La figure 3 représente un exemple de désentrelacement non supervisé utilisant des classifications RF-DOA suivi par des classifications HDTOA ; Figure 3 shows an example of unsupervised deinterleaving using RF-DOA classifications followed by HDTOA classifications;

La figure 4 représente un exemple de mise en œuvre du procédé selon l'invention ; FIG. 4 represents an exemplary implementation of the method according to the invention;

La figure 5 représente un exemple de représentation graphique d'un histogramme RF-DTOA ; Fig. 5 shows an exemplary graphical representation of an RF-DTOA histogram;

La figure 6 représente un exemple de mise en œuvre de l'étape d'extraction de modes ; FIG. 6 represents an exemplary implementation of the mode extraction step;

- Les figures 7a et 7b illustrent l'étape d'extraction des modes d'un histogramme ; FIGS. 7a and 7b illustrate the step of extracting the modes of a histogram;

La figure 8 représente différents choix de classification possible ; La figure 9 représente un exemple de structure pyramidale illustrant les liens de corrélation entre impulsions ; Figure 8 represents different possible classification choices; FIG. 9 represents an example of a pyramidal structure illustrating the correlation links between pulses;

La figure 10 représente un exemple de structure pyramidale illustrant les liens de corrélation entre modes. Figure 10 shows an example of a pyramidal structure illustrating the correlation links between modes.

Le procédé de désentrelacement non supervisé par enrichissement N dimensions selon l'invention repose sur l'utilisation conjointe des paramètres primaires et des paramètres secondaires comme par exemple la différence de temps d'arrivée (DTOA) pour séparer chaque émission d'intérêt. The method of non-supervised deinterleaving by N-size enrichment according to the invention is based on the joint use of the primary parameters and secondary parameters such as, for example, the difference in arrival time (DTOA) for separating each emission of interest.

En référence à la figure 4 le procédé de désentrelacement selon l'invention peut comprendre une étape Etp10 de constitution d'histogrammes, une étape Etp20 d'extraction de modes, une étape Etp30 de constitution de groupes d'intérêt, une étape Etp40 d'enrichissement de la description des impulsions et une étape Etp50 de classement non supervisé des modes. With reference to FIG. 4, the deinterleaving method according to the invention may comprise a step Etp10 for constituting histograms, a mode extraction step Etp20, an interest group formation step Etp30, a step Etp40 of enrichment of the description of the pulses and an Etp50 step of unsupervised classification of the modes.

Au cours de l'étape Etp10 des histogrammes sont constitués dans toutes les dimensions. Cette étape permet de réaliser une description multidimensionnelle conjointe (fréquence, DTOA, DOA, phase des trains temporels, paramètres intra-impulsion, Ll ....) d'un ensemble d'impulsions. Cette description conjointe permet en particulier de construire des histogrammes liant paramètres primaires, comme par exemple les fréquences, les DOA, les largeurs d'impulsion (Ll), ... et secondaires comme par exemple les PRI. On peut aussi ajouter des modes de regroupement comme la PRA, et plus généralement des indices statistiques construits en temps long comme des groupes de fréquence ou PRI cooccurentes, des groupes de point d'accumulation, des iso-mesures de type localisation. Pour comprendre le principe, le plus didactique de ces histogrammes est l'histogramme RF-DTOA associant un paramètre primaire c'est à dire mono-impulsion et un paramètre secondaire calculé sur 2 impulsions. A titre d'illustration, la figure 5 représente un tel histogramme dans un repère ayant pour abscisse les DTOA et pour ordonnée les fréquences. Dans cette représentation, la taille des points est relative au nombre d'impulsions se trouvant dans la cellule (Fréquence, DTOA). On remarque que cette représentation permet de distinguer aisément les différentes formes d'ondes (FO) entrelacées. Pour cela on s'appuie sur les alignements horizontaux ou verticaux de points observés, ce qui revient à s'intéresser aux projections sur les deux axes (ici fréquence et DTOA) et aux liens entre les modes en projection. During step Etp10 histograms are constituted in all dimensions. This step makes it possible to perform a joint multidimensional description (frequency, DTOA, DOA, phase of the time trains, intra-pulse parameters, L1, etc.) of a set of pulses. This joint description makes it possible in particular to construct histograms linking primary parameters, such as, for example, frequencies, DOAs, pulse widths (L1), etc. and secondary frequencies, for example PRIs. It is also possible to add grouping modes such as PRA, and more generally statistical indices built in long time as groups of frequency or PRI cooccurentes, groups of accumulation point, iso-measures of type localization. To understand the principle, the most didactic of these histograms is the histogram RF-DTOA associating a primary parameter that is to say mono-pulse and a secondary parameter calculated on 2 pulses. By way of illustration, FIG. 5 represents such a histogram in a reference frame having for abscissa the DTOAs and for ordinate the frequencies. In this representation, the size of the points is relative to the number of pulses found in the cell (Frequency, DTOA). Note that this representation makes it easy to distinguish the different interlaced (FO) waveforms. For this we rely on the horizontal or vertical alignment of observed points, which is to be interested in the projections on the two axes (here frequency and DTOA) and the links between the modes in projection.

Il est également aisé de reconnaître les structures temporelles ou fréquentielles les plus contrastées par rapport aux structures diffuses c'est-à- dire les structures comportant le plus d'impulsions. On pourra donc commencer le marquage en s'intéressant aux signaux aux caractéristiques les plus spécifiques donc facilement extractibles du mélange sans risque d'ambiguïté It is also easy to recognize the most contrasted temporal or frequency structures with respect to diffuse structures, that is, the structures with the most pulses. We can therefore start marking by looking at signals with the most specific characteristics so easily extractable from the mixture without risk of ambiguity

Bien entendu les histogrammes ne se limitent pas à ceux en deux dimensions (2D) et sur la représentation graphique on aurait pu ajouter Azimut (en z), Ll en(v)... On parle d'histogrammes à N Dimensions ou ND. Of course the histograms are not limited to those in two dimensions (2D) and on the graphical representation we could have added Azimut (in z), Ll in (v) ... We speak of histograms to N Dimensions or ND.

En référence à la figure 6, l'étape Etp20 d'extraction de modes peut comprendre une étape Etp 21 de seuillage des histogrammes, une étape Etp22 de découpe en modes et une étape Etp23 de caractérisation des modes. With reference to FIG. 6, the mode extraction step Etp20 may comprise a thresholding step Etp 21 of the histograms, a mode cutting step Etp22 and a mode characterization step Etp23.

Pour le seuillage des histogrammes, on peut distinguer deux cas, le cas à une dimension (1 D) et le cas à deux dimensions (2D) dans lequel on utilise systématiquement la DOA. For the thresholding of the histograms, one can distinguish two cases, the one-dimensional case (1 D) and the two-dimensional case (2D) in which the DOA is systematically used.

En référence aux figures 7a et 7b, les histogrammes sont seuillés de façon à éliminer les valeurs non significatives, c'est-à-dire inférieur à une valeur seuil 70 prédéterminée. With reference to FIGS. 7a and 7b, the histograms are thresholded so as to eliminate the insignificant values, that is to say less than a predetermined threshold value 70.

L'étape Etp 22 de découpe en modes consiste à sélectionner des séparatrices dans les " vallées " 71 de l'histogramme qui se trouvent entre deux pics 72 et pour lesquelles la différence de hauteur entre la vallée 71 et les pics 72 est significative. La signifiance statistique peut être déterminée selon un test du khi2, par exemple avec la formule : The mode etch step 22 consists in selecting separators in the "valleys" 71 of the histogram that are between two peaks 72 and for which the difference in height between the valley 71 and the peaks 72 is significant. The statistical significance can be determined by a chi square test, for example with the formula:

2(HV - HE)² _{^ 2} Où HV est égal à la hauteur de l'histogramme dans la " vallée " 71 et HE est la valeur moyenne des hauteurs dans la " vallée " 71 et pour le pic de plus faible hauteur. Pour une confiance de 95%, χ ₀₅₄ = 3.843 . 2 (HV - HE) ² _{^ 2} Where HV is equal to the height of the histogram in the "valley" 71 and HE is the average value of the heights in the "valley" 71 and for the lower peak height. For 95% confidence, χ ₀₅₄ = 3.843.

Chaque mode est ensuite caractérisé par sa moyenne et son écart type, par exemple par approximation gaussienne. Les modes sont étendues au-delà des séparatrices 73 en retenant les données à par exemple, plus ou moins 3 sigmas estimés. Dans la plupart des cas les modes extraits ne présentent pas de mélange. Quand deux pics sont proches on décide que les données appartiennent à la fois aux deux pics. Each mode is then characterized by its mean and its standard deviation, for example by Gaussian approximation. The modes are extended beyond the separators 73 by retaining the data at, for example, plus or minus 3 sigmas estimated. In most cases the extracted modes do not show a mixture. When two peaks are close we decide that the data belongs to both peaks.

Dans le cas où l'on utilise systématiquement l'Azimut pour chaque dimension (cas 2D) l'étape d'extraction des modes peut utiliser un algorithme de détection des maxima puis une découpe des modes autour de ces maxima par exemple par une méthode de segmentation par ligne de partage des eaux. In the case where Azimuth is systematically used for each dimension (2D case), the mode extraction step can use an algorithm for detecting maxima and then cutting modes around these maxima, for example by a method of segmentation by watershed.

Les modes de ces histogrammes définissent des groupes (ou clusters selon la terminologie anglo saxonne) d'impulsions par paramètre. Chaque groupe peut être décrit par la valeur moyenne du paramètre, l'écart type, l'effectif et par la liste des impulsions consécutives. On peut calculer ensuite le contraste de chaque mode et ces modes peuvent être classés par contraste décroissant. Ainsi, des groupes d'impulsions ou groupes d'intérêt sont formés au cours d'une étape Etp30, par exploration des modes toutes dimensions par ordre d'effectif décroissant. Les groupes formés sont caractérisés par leurs paramètres comme par exemple DOA, liste RF, PRM, Phase, paramètres intra impulsions.... The modes of these histograms define groups (or clusters according to the Anglo-Saxon terminology) of pulses per parameter. Each group can be described by the average value of the parameter, the standard deviation, the number and the list of consecutive pulses. The contrast of each mode can then be calculated and these modes can be classified by decreasing contrast. Thus, groups of pulses or groups of interest are formed during a step Etp30, by scanning all-size modes in decreasing order of strength. The formed groups are characterized by their parameters such as DOA, RF list, PRM, Phase, intra pulse parameters ....

Chaque impulsion est ensuite enrichie au cours d'une étape Etp40 d'enrichissement de la description des impulsions par son appartenance à un groupe (F, DOA) ou à plusieurs groupes (PRM, phase, DOA). Ainsi après l'enrichissement à N dimensions on dispose des modes dans les différentes dimensions et des impulsions à description enrichie. Each pulse is then enriched during a step Etp40 for enriching the description of the pulses by belonging to a group (F, DOA) or to several groups (PRM, phase, DOA). Thus after the N-dimensional enrichment modes are available in the different dimensions and pulses with enriched description.

Les impulsions enrichies peuvent se présenter sous la forme d'une table des impulsions enrichies. A titre illustratif l'annexe 1 présente un exemple d'une telle table. Les groupes d'impulsions formés dans les différentes dimensions ou modes sont non seulement caractérisés par leur valeur moyenne, écart type expérimental, effectif mais également par la liste des labels des impulsions appartenant au mode comme illustré annexe 2. The enriched pulses can be in the form of a table of enriched pulses. As an illustration, Appendix 1 presents an example of such a table. The groups of pulses formed in the different dimensions or modes are not only characterized by their mean value, experimental standard deviation, effective but also by the list of labels of pulses belonging to the mode as illustrated in Appendix 2.

Dans un autre mode de réalisation de l'invention, cette liste de labels peut être codée sous forme d'une table de valeurs binaires indiquant l'appartenance ou non de l'impulsion i au mode. In another embodiment of the invention, this list of labels may be coded in the form of a table of binary values indicating whether or not the pulse i belongs to the mode.

De façon avantageuse, la mémorisation des labels des impulsions sur les différents modes permet le calcul de liens de corrélation entre modes. Ceci permet d'éviter d'avoir à itérer sur des tris séquentiels d'impulsions. Advantageously, the memorization of the labels of the pulses on the different modes allows the calculation of correlation links between modes. This avoids having to iterate on sequential sorting of pulses.

La méthode repose sur la dualité individus/variables. Les impulsions sont ici les individus et les différents modes (DTOA, Fréquences, DOA, Ll ...) sont les variables. Les variables seront traitées en utilisant les individus qui sont liés à chaque valeur de la variable (mode) afin d'obtenir diverses informations. The method is based on the duality individuals / variables. The impulses are here the individuals and the different modes (DTOA, Frequencies, DOA, L1 ...) are the variables. The variables will be processed using the individuals that are linked to each value of the variable (mode) in order to obtain various information.

L'intérêt de l'enrichissement réside dans le fait de mesurer les propriétés des modes et leurs relations, c'est-à-dire leur appartenance à un même train, par des opérateurs d'intersection et d'union sur les listes d'impulsions appartenant à chaque mode significatif. The interest of the enrichment lies in the fact of measuring the properties of the modes and their relations, that is to say their belonging to the same train, by operators of intersection and union on the lists of impulses belonging to each significant mode.

On distingue le lien de corrélation intra dimension, par exemple entre deux fréquences ou deux PRM, et le lien de corrélation inter dimension, par exemple entre une fréquence et une PRM. The intra-dimensional correlation link, for example between two frequencies or two PRMs, and the inter-dimensional correlation link, for example between a frequency and a PRM, are distinguished.

Pour les paramètres bi-impulsions, le lien de corrélation entre deux modes peut être donnée par la distance ensembliste. Le lien de corrélation entre deux PRM peut être donnée par la formule: For bi-pulse parameters, the correlation link between two modes can be given by the set distance. The correlation link between two PRMs can be given by the formula:

di_{R R =} card (Rt u R₂ ) - card {R_X n R₂ ) Pour les paramètres mono-impulsion, le lien de corrélation entre deux modes peut être donnée par la distance euclidienne classique. Le lien de corrélation inter dimension entre une PRM détectée et une fréquence peut être réalisée en considérant la cardinalité de l'intersection entre les listes d'imp e considérées.

di _{RR =} card (Rt u R ₂ ) - card {R _X n R ₂ ) For single-pulse parameters, the correlation link between two modes can be given by the classical Euclidean distance. The inter-dimensional correlation link between a detected PRM and a frequency can be performed by considering the cardinality of the intersection between the considered tax lists.

Sur la base de ces liens de corrélation inter et intra-dimensions il est ainsi possible de regrouper les modes qui sont simultanément proches. On the basis of these inter and intra-dimensional correlation links it is thus possible to group the modes that are simultaneously close.

Une propriété importante de l'enrichissement est qu'on peut retrouver les modes multidimensionnels conjoints à partir des modes de faible dimension construits et du tableau des impulsions enrichies. Cette mémorisation est de plus linéaire c'est-à-dire qu'elle reste proportionnelle au nombre des impulsions. Un codage direct de la statistique multidimensionelle serait beaucoup plus coûteux en temps de calcul. An important property of enrichment is that one can find joint multidimensional modes from built-in low-dimensional modes and enriched-pulse array. This memorization is moreover linear, that is to say that it remains proportional to the number of pulses. A direct coding of the multidimensional statistics would be much more expensive in computing time.

En référence à l'annexe 4, l'étape Etp50 a pour but de réorganiser les impulsions enrichies ou tableaux impulsions-modes de façon à faire ressortir les différentes classes d'impulsions. Cette réorganisation utilise une classification croisée dont l'objectif principal est de fournir simultanément une partition des deux ensembles individus et variables. Cette classification peut être réalisée par plusieurs méthodes suivant le type de tableaux de données envisagés (tableau binaire, tableau de contingence...). With reference to Annex 4, Etp50 aims to reorganize the enriched pulses or pulse-mode tables so as to bring out the different classes of pulses. This reorganization uses a cross classification whose main objective is to simultaneously provide a partition of the two sets individuals and variables. This classification can be carried out by several methods depending on the type of data tables envisaged (binary table, contingency table ...).

Le but recherché est de former des classes ou blocs d'impulsions- modes homogènes. On veut que les impulsions d'une même classe se ressemblent et ne ressemblent pas à celles des autres classes. Par symétrie impulsions/modes on obtient également un critère sur les modes : les modes d'une même classe se ressemblent car leurs listes d'impulsions présentent de fortes intersections. The aim is to form classes or blocks of pulses-homogeneous modes. We want that the impulses of the same class are similar and do not resemble those of the other classes. By symmetry pulses / modes we also obtain a criterion on the modes: the modes of the same class are similar because their lists of pulses have strong intersections.

On défini donc les critères de similarité suivants : The following similarity criteria are defined:

Dans une même classe: In the same class:

- Les impulsions partagent un maximum de composantes ; - The impulses share a maximum of components;

- Les composantes partagent un maximum d'impulsions. - The components share a maximum of pulses.

Entre deux classes: - Les impulsions partagent un minimum de composantes ; Between two classes: - The impulses share a minimum of components;

- Les composantes partagent un minimum d'impulsions. - The components share a minimum of impulses.

De façon à accélérer le classement on dispose également de contraintes. A titre d'exemple, ces contraintes peuvent être : In order to accelerate the classification one also has constraints. For example, these constraints can be:

- Deux composantes DOA (idem pour deux composantes PRM) " disjointes " ne peuvent appartenir à une même classe ; - Two DOA components (idem for two PRM components) "disjointed" can not belong to the same class;

- Les fréquences d'une même " phase de PRM " ne doivent pas être séparées ; Etc.. - The frequencies of the same "PRM phase" must not be separated; Etc ..

Le traitement peut s'effectuer en utilisant seulement les critères de similarité ou en utilisant les critères de similarité et les contraintes. L'utilisation des critères de similarité permet à lui seul d'aboutir au résultat optimal, les contraintes ne sont qu'un accélérateur de calcul. Dans certains cas, le classement pourra être effectué en n'utilisant que les contraintes. Le critère de similarité entre impulsions ou lien de corrélation entre deux impulsions peut être donné par la formule : ∑ ^;.v.; Processing can be done using only similarity criteria or using similarity criteria and constraints. The use of similarity criteria alone leads to the optimal result, the constraints are just a computational accelerator. In some cases, the classification can be done using only the constraints. The criterion of similarity between pulses or correlation link between two pulses can be given by the formula: Σ ^; .v .;

Où xl représente un indicateur binaire indiquant si l'impulsion i appartient au mode j défini par : Where x1 represents a binary indicator indicating whether the pulse i belongs to the mode j defined by:

x{ =\ si l'impulsion i appartient au mode j x {= \ if the pulse i belongs to the mode j

xj =0 si l'impulsion i n'appartient pas au mode j xj = 0 if the impulse i does not belong to the mode j

Et où P représente le nombre total de modes. And where P represents the total number of modes.

De même, de façon symétrique, on peut définir un critère de similarité entre modes par la formule :

Similarly, symmetrically, a similarity criterion between modes can be defined by the formula:

Où xl représente un indicateur binaire indiquant si le mode j contient l'impulsion i défini par : Where x1 represents a binary indicator indicating whether the mode j contains the pulse i defined by:

J =1 si le mode j contient l'impulsion i J = 1 if mode j contains the pulse i

J / =0 si l'impulsion i n'appartient pas au mode j J / = 0 if the impulse i does not belong to the mode j

Et où N représente le nombre total d'impulsions. Calculer le critère de similarité revient à faire le produit scalaire entre deux impulsions ou deux modes. L'annexe 5a illustre par des exemples le critère de similarité entre les impulsions et l₂ ainsi que le critère pour les impulsions l₃ et l₄. C_{/1 /2} =4 signifie que les impulsions \ ^ et l₂ sont toutes les deux présentes dans 4 modes (cf. cases entourées). And where N is the total number of pulses. To calculate the similarity criterion is to make the scalar product between two pulses or two modes. Appendix 5a illustrates by examples the criterion of similarity between the pulses and l _{2 as} well as the criterion for the pulses l ₃ and l ₄ . C _{/ 1/2} = 4 means that pulses \ ^ and l ₂ are both present in 4 modes (see surrounded boxes).

De même, l'annexe 5b présente le critère de similarité entre les modes AOA1 et PRM1 ainsi que celui entre les modes AOA2 et F1 . Similarly, Appendix 5b presents the criterion of similarity between the modes AOA1 and PRM1 as well as that between the modes AOA2 and F1.

Pour réaliser un groupe de modes-impulsions homogène ou groupement d'intérêt, on va définir un critère de regroupement. Pour cela, on recherche une partition en classes où les individus de la classe sont les plus similaires les uns aux autres et les individus de classes différentes les plus dissimilaires. Une partition des impulsions peut être définie par une matrice carrée binaire X dont le terme général est défini par: To achieve a group of homogeneous pulse modes or group of interest, we will define a grouping criterion. For this, we search for a partition in classes where the individuals of the class are the most similar to each other and the individuals of the most dissimilar different classes. A partition of the pulses can be defined by a binary square matrix X whose general term is defined by:

il si i et i' sont dans la même classe it if i and i 'are in the same class

[0 sinon [0 otherwise

Plusieurs critères peuvent correspondre à cette définition. Parmi ceux- ci, se trouve le critère de Condorcet. Il peut s'écrire de la façon suivante : Several criteria may correspond to this definition. Among these is the Condorcet criterion. It can be written as follows:

^CW =∑ i=l∑ i'=l (^CA + ¾) ^C W = Σ i = lΣ i '= l ( ^C A + ¾)

Où : N représente le nombre d'impulsions, Where: N represents the number of pulses,

C_w représente le critère de similarité entre les impulsions l _n et Ι _η·C _w represents the similarity criterion between the pulses l _n and Ι _η ·

X_u, représente l'élément de rang i,i' de la matrice carré X précédemment définie. X _u , represents the element of rank i, i 'of the square matrix X previously defined.

Les tableaux annexe 6b et 6c présentent des exemples de partition en différentes classes d'une classification de 20 impulsions et 1 0 modes présentée annexe 6a. Appendix Tables 6b and 6c show examples of partitioning into different classes of a classification of 20 pulses and 10 modes presented in Annex 6a.

L'annexe 6b représente une classification des impulsions et des modes en 3 classes Cl-i , Cl₂ et Cl₃. Dans l'exemple les classes 3 et 2 se partagent la fréquence F2, et les classes Cl₂ et C\^ se partagent les fréquences F3 et F4. Pour cette partition, le critère est de 1 076. Ce classement est optimal c'est-à-dire qu'il présente un critère maximal. L'annexe 6c présente un classement alternatif en 5 classes différentes CI-1 1 , CI-I 2 , CI-I 3 , CI™ et CI 15- Cette partition n'est pas optimale, la valeur du critère pour cette partition est 912 donc inférieure à 1076. Un premier exemple de mode de mis en œuvre de l'étape Etp50 de classement non supervisé va être présenté. Ce classement des modes consiste principalement à créer progressivement au moins une classe d'impulsions-modes homogènes. Ce classement comprend à l'initialisation une recherche du mode ayant un maximum d'intérêt c'est-à-dire le mode ayant le plus grand cardinal. Ce mode va former la première classe, puis, pour chaque mode restantje classement comprend : Annex 6b shows a classification of pulses and modes in 3 classes Cl-i, Cl ₂ and Cl ₃ . In the example the classes 3 and 2 share the frequency F2, and the classes Cl ₂ and C \ ^ share the frequencies F3 and F4. For this partition, the criterion is 1 076. This classification is optimal, that is to say that it presents a maximum criterion. Annex 6c presents an alternative classification in 5 different classes CI-1 1, CI-I 2, CI-I 3, CI ™ and CI 15- This partition is not optimal, the criterion value for this partition is 912 therefore less than 1076. A first example of implementation mode Etp50 step of unsupervised classification will be presented. This classification of modes consists mainly of gradually creating at least one class of homogeneous mode pulses. This classification includes at initialization a mode search with a maximum of interest, that is to say the mode with the largest cardinal. This mode will form the first class, then for each remaining mode I rank includes:

L'initialisation de l'étape consiste à rechercher le mode qui est le plus intéressant à prendre en premier, c'est-à-dire le mode dont la contribution au critère est la plus forte. Ce mode formera la première classe C-i . The initialization of the step consists in finding the mode which is the most interesting to take first, that is to say the mode whose contribution to the criterion is the strongest. This mode will form the first class C-i.

On rappelle que la contribution au critère du mode q est donnée par la formule : Recall that the contribution to the mode q criterion is given by the formula:

Contrib(q) = ∑C_i{ - a ∑ min(C, , C_{{ ) Où: C_M* représente la similarité ou le lien de corrélation entre modes i et i', Contrib (q) = ΣC _{i {} - a Σ min (C,, C _{{ ) Where: C _M * represents the similarity or the correlation link between modes i and i ',

a représente un scalaire que l'on peut prendre égale à 1 a represents a scalar that can be equal to 1

On s'intéresse ensuite à l'intersection entre cette première classe et les modes restants, on s'intéresse plus particulièrement au cardinal de cette intersection. We are then interested in the intersection between this first class and the remaining modes, we are particularly interested in the cardinal of this intersection.

Suivant un premier mode de mise en œuvre du procédé on peut classer les modes restant par ordre décroissant des cardinaux d'intersection. Suivant un autre mode de mise en œuvre, on recherche parmi les modes restant, celui dont le cardinal d'intersection est le plus grand. Ce mode formera la classe C₂. According to a first mode of implementation of the method, the remaining modes can be classified in descending order of the intersection cardinals. According to another mode of implementation, one looks among the modes remaining, the one whose cardinal intersection is the largest. This mode will form the class C ₂ .

En référence à la figure 8, si on considère les classes Ci et C₂, et C_y=C₂), il existe quatre choix possible pour le regroupement de ces classes. L'intersection d nc₂ peut appartenir à la nouvelle classe C₂ (choix 1 ), à la nouvelle classe d (choix 2) ou former une nouvelle classe C₃ (choix 3). Une quatrième option consiste à considérer un nouvelle classe C'i regroupant les classes Ci et C₂ (choix 4). With reference to FIG. 8, if we consider the classes Ci and C ₂ , and C _y = C ₂ ), there are four possible choices for the grouping of these classes. The intersection d nc ₂ can belong to the new class C ₂ (choice 1), to the new class d (choice 2) or to form a new class C ₃ (choice 3). A fourth option is to consider a new class C'i grouping classes Ci and C ₂ (choice 4).

On considère donc ces quatre options et pour chacune d'elle on calcule le critère de regroupement correspondant et on retient le choix qui maximise ledit critère de regroupement. We therefore consider these four options and for each of them we calculate the corresponding grouping criterion and we retain the choice that maximizes the grouping criterion.

On continue ainsi, pour chaque mode restant, en considérant à chaque fois l'intersection entre la ou les classes (d_, C'_2i C₃ ) que l'on vient de former et le mode ayant le cardinal d'intersection de plus grand avec la ou les nouvelles classes. We thus continue, for each remaining mode, considering each time the intersection between the class or classes (d _, C ' _2i C ₃ ) that we have just formed and the mode having the intersection cardinal of greater with the new class or classes.

Si le choix 1 est retenu, les nouvelles classes formées sont: d = C-_\ - C₂ et C'₂= C₂ ; If choice 1 is retained, the new classes formed are: d = C- _\ - C ₂ and C ' ₂ = C ₂ ;

Si le choix 2 est retenu, les nouvelles classes deviennent : d = Ci et

If choice 2 is retained, the new classes become: d = Ci and

Si le choix 3 est retenu, les nouvelles classes deviennent : d = C-i , C'₂= C₂ et d= C-| flC₂ ; If choice 3 is retained, the new classes become: d = Ci, C ' ₂ = C ₂ and d = C- | flC ₂ ;

Si le choix 4 est retenu, la nouvelle classe devient : d = Ci UC₂. If choice 4 is retained, the new class becomes: d = Ci UC ₂ .

On ne s'intéresse qu'aux intersections non nulles. Dans le cas où le mode considéré aurait une intersection nulle, il formerait une nouvelle classe. We are only interested in non-zero intersections. In the case where the mode considered would have a zero intersection, it would form a new class.

Comme cité précédemment, on possède également des contraintes pour accélérer le calcul. Donc pour chaque mode considéré, on recherche le choix qui maximise le critère et qui respecte les contraintes. As mentioned above, there are also constraints to speed up the calculation. So for each mode considered, we look for the choice that maximizes the criterion and respects the constraints.

Il est à noter que l'on classe plusieurs données à la fois. En effet chaque classe Cx et Cy sont des classes d'impulsions. En référence aux annexes 8 à 15, l'étape Etp 50 de classement non supervisé va être illustrée à travers un exemple nullement limitatif. Dans cet exemple on va considérer 20 impulsions à \₂₀. It should be noted that several data are sorted at a time. Indeed each class Cx and Cy are classes of pulses. With reference to Annexes 8 to 15, the step Etp 50 unsupervised classification will be illustrated through a non-limiting example. In this example we will consider 20 pulses at \ ₂₀ .

L'annexe 6a correspond à l'information obtenue après le traitement d'enrichissement des impulsions. Ce tableau croise les impulsions et les différents modes (DOA,, F,, PRM,). Appendix 6a corresponds to the information obtained after the pulse enrichment treatment. This table crosses the pulses and the different modes (DOA ,, F ,, PRM,).

Après classement des modes par ordre croissant de leur contribution on obtient le tableau annexe 7. On commence le classement non supervisé avec le mode 10 (PRM3). After ranking the modes in ascending order of their contribution, we obtain the appendix table 7. We begin the unsupervised classification with mode 10 (PRM3).

Ce mode contient les impulsions ₂ à \₂₀ qui vont former la première classeThis mode contains pulses ₂ to \ ₂₀ which will form the first class

C-i . On considère ensuite le mode dont le cardinal d'intersection avec la classe C1 est le plus élevé à savoir le mode 9 (DOA3). Les modes 9 et 10 sont identiques, on est donc dans le cas de figure où C_x nC_y =C_x ^>C_y =C_x = C_y , les 4 choix sont donc identiques. Les modes 10 etThis . We then consider the mode whose cardinal intersection with the class C1 is the highest namely mode 9 (DOA3). The modes 9 and 10 are identical, so we are in the case where C _x nC _y = C _x ^> C _y = C _x = C _y , the 4 choices are identical. Modes 10 and

9 forment donc une seule classe Q= {/₁₂,/₁₃,/₁₄,/₁₅,/₁₆,/₁₇,/₁₈,/₁₉,/₂₀} (cf. annexe 8). 9 therefore form a single class Q = {/ ₁₂ , / ₁₃ , / ₁₄ , / ₁₅ , / ₁₆ , / ₁₇ , / ₁₈ , / ₁₉ , / ₂₀ } (see Annex 8).

On considère ensuite le mode 8 (F4) qui est le mode ayant la plus grande intersection avec la classe C1 { C₁nF₄ = 5 ) (cf. annexe 9). Then we consider mode 8 (F4) which is the mode with the greatest intersection with the class C1 {C ₁ nF ₄ = 5) (see appendix 9).

On teste les 4 choix suivants: We test the following 4 choices:

• C_v = C -F₄ , F₄, =F₄ (critère 243) • C _v = C -F ₄ , F ₄ , = F ₄ (criterion 243)

• C_v = C_x , F₄, = F₄-C_x (critère 303) • C _v = C _x , F ₄ , = F ₄ -C _x (criterion 303)

• C_r = C -F₄ , F₄, = F₄-C , C_z, = nF₄ (critère 263) • C _r = C -F ₄ , F ₄ , = F ₄ -C, C _z , = nF ₄ (criterion 263)

· C_r = C₁ F₄ (critère 235) · C _r = C ₁ F ₄ (criterion 235)

Le critère le plus élevé étant 303, on retient la solution consistant à séparer les modes en deux classes, à savoir The highest criterion being 303, we retain the solution of separating the modes into two classes, namely

Q ⁼ {^12 ' ^13 ' ^14 ' s ' ' ' ^19 ' -^20 } ®t C₂ = {/₁₀ , I_n } . On poursuit (cf. annexe 10 a) en considérant le mode 7 (F3) qui est le mode suivant ayant la plus grande intersection avec la nouvelle classe C1Q ⁼ {^ 12 '^ 13' ^ 14 's'''^19' - ^ ® T 20} C ₂ = {/ _10, I _n}. We continue (see appendix 10a) considering mode 7 (F3) which is the next mode having the greatest intersection with the new class C1

( Qn ₃ = 3 ) (Qn ₃ = 3)

Comme précédemment on teste les 4 choix suivant: As before we test the following 4 choices:

• C_v = C -F₃ , F_y =F₃ (critère 251 ) • C _v = C -F ₃ , F _y = F ₃ (criterion 251)

· C_r = C_x , F₃, = (critère 307) • C_y = C-F₃ , F_y = F₃- , C_z. = nF₃ (critère 267) · C _r = C _x , F ₃ , = (criterion 307) • C _y = CF ₃ , F _y = F ₃ -, C _z . = nF ₃ (criterion 267)

• C^ F, (critère 231) • C ^ F, (criterion 231)

A l'issue de cette étape (cf. annexe 10 b), en prenant en compte le critère le plus élevée, on aboutit à 3 nouvelles classes : j= {l₁₂,I₁₃,I₁₄,I₁₅,I₁₆,I₁₇,I_1&,I₁₉,I₂₀} , C₂={l_w,I } et C₃={/₈,/₉}. At the end of this step (see Annex 10b), taking into account the highest criterion, we obtain 3 new classes: j = {l ₁₂ , I ₁₃ , I ₁₄ , I ₁₅ , I ₁₆ , I ₁₇ , I ₁ , I ₁₉ , I ₂₀ }, C ₂ = {I _w , I} and C ₃ = {/ ₈ , / ₉ }.

On peut remarquer qu'en utilisant les contraintes en PRM ou en DOA (cf. annexe 11), à savoir "deux composantes DOA (ou deux composantes PRM) disjointes ne peuvent appartenir à une même classe aurait permis de réaliser plus vite cette étape pour le même résultat. La prise en compte de ces contraintes aurait interdit la formation d'une même classe avec des modes n'ayant pas la même DOA. Le même raisonnement peut être appliqué pour des modes n'ayant pas la même PRM. On poursuit (cf. annexe 12a) en considérant le mode 6It can be noticed that using PRM or DOA constraints (see Appendix 11), ie "two disassociated DOA components (or two PRM components) can not belong to the same class would have made it possible to achieve this step more quickly". The same result: Taking into account these constraints would have forbidden the formation of the same class with modes that do not have the same DOA The same reasoning can be applied for modes that do not have the same PRM. (see Annex 12a) considering mode 6

(C₂nPRM₂ = 2) et on teste les 4 choix suivants : (C ₂ nPRM ₂ = 2) and the following 4 choices are tested:

• C₂, = C₂-PRM₂ = 0 , PRM₂, =PRM₂ (critère 84) • C ₂ , = C ₂ -PRM ₂ = 0, PRM ₂ , = PRM ₂ (criterion 84)

• C₂, = C₂ , PRM₂, =PRM₂-C₂ (critère 68) • C ₂ , = C ₂ , PRM ₂ , = PRM ₂ -C ₂ (criterion 68)

• C₂, = C₂-PRM₂ = 0 , PRM₂, =PRM₂ - C₂ , C_z, = C₂nPRM₂ = C₂ (identique au choix précédent) • C ₂ , = C ₂ -PRM ₂ = 0, PRM ₂ , = PRM ₂ - C ₂ , C _z , = C ₂ nPRM ₂ = C ₂ (identical to the previous choice)

• C₂, = C₂ PRM₂ = PRM₂ (identique au premier choix)• C ₂ , = C ₂ PRM ₂ = PRM ₂ (identical to the first choice)

A ce stade on a 3 classes : Q= {/₁₂,/₁₃,/₁₄,/₁₅,/₁₆,/₁₇,/₁₈,/₁₉,/₂₀}, €₂={ΐ₆,Ι_Ί,Ι_&,Ι₉,Ι_ιο,Ι ) et C₃={/₈,/₉}. On continue en testant l'intersection du mode 6 avec la classe C3 {C₃ PRM₂ = 2). La résolution est identique, on n'a que deux choix. At this point we have 3 classes: Q = {/ ₁₂ , / ₁₃ , / ₁₄ , / ₁₅ , / ₁₆ , / ₁₇ , / ₁₈ , / ₁₉ , / ₂₀ }, € ₂ = {ΐ ₆ , Ι _Ί , Ι _& , Ι ₉ , Ι _ιο , Ι) and C ₃ = {/ ₈ , / ₉ }. We continue by testing the intersection of mode 6 with the class C3 {C ₃ PRM ₂ = 2). The resolution is the same, we only have two choices.

A l'issue de cette étape (cf. annexe 12b) on a donc 2 classes à savoir:

Les anciennes classes C₂ et C₃ qui étaient séparées se retrouvent maintenant réunies ensemble. At the end of this step (see appendix 12b) we thus have 2 classes namely:

The old classes C ₂ and C ₃ which were separated are now together.

On continue en considérant le mode 5 (DOA2). Comme DOA2=C2 la solution reste identique. On poursuit (cf. annexe 13) en considérant le mode 4 (F2). L'intersection entre le mode 4 et la classe 1 étant nulle {C₁nF₂ = 0), il n'y a aucun test à effectuer. C₂nF₂ = 2, on teste les 4 choix : We continue by considering mode 5 (DOA2). As DOA2 = C2 the solution remains the same. We continue (see appendix 13) considering mode 4 (F2). The intersection between mode 4 and class 1 being zero {C ₁ nF ₂ = 0), there is no test to perform. C ₂ nF ₂ = 2, we test the 4 choices:

. C₂, = C₂-F₂, F₂,=F₂ (critère 179) . C ₂ , = C ₂ -F ₂ , F ₂ , = F ₂ (criterion 179)

· C₂, = C₂, F₂, = F₂-C₂ (critères 207) · C ₂ , = C ₂ , F ₂ , = F ₂ -C ₂ (criteria 207)

• C₂, = C₂-F₂ , F₂, = F₂-C₂ , C_z, = C₂nF₂ (critère 191) • C ₂ , = C ₂ -F ₂ , F ₂ , = F ₂ -C ₂ , C _z , = C ₂ n F ₂ (criterion 191)

• C₂, = C_{2 2} (critère 123) es :

• C ₂ , = C _{2 2} (criterion 123) es:

On considère maintenant le mode 3 (F1). Ce mode n'intervient pas sur les classes Ci, C₂ et C₃ (Q F^O C₂nf_{1 =} 0 et C₃n ₁ =0). We now consider mode 3 (F1). This mode does not intervene on the classes Ci, C ₂ and C ₃ (QF ^ OC ₂ nf _{1 =} 0 and C ₃ n ₁ = 0).

On poursuit (cf. annexe 14) avec le mode 2 (PRM1). Ce mode n'intervient ni sur la classe Ci ni sur la classe C₃ {C₁r)PRM₁=0 et C₂(^PRM₁=Q). Le mode 2 a une intersection non nulle avec le classe C₃ {C₃ PRM₁ =3), on teste donc les 4 choix : We continue (see appendix 14) with mode 2 (PRM1). This mode does not intervene on the class Ci nor on the class C ₃ {C ₁ r) PRM ₁ = 0 and C ₂ (^ PRM ₁ = Q). Mode 2 has a non-zero intersection with the class C ₃ {C ₃ PRM ₁ = 3), so we test the 4 choices:

• C_y = C₃-PRM₁ = 0 , PRM_V =PRM₁ (critère 63) • C _y = C ₃ -PRM ₁ = 0, PRM _V = PRM ₁ (criterion 63)

• C_y = C₃ , PRM_V = PR ₁-C₃ (critère 51 ) • C _y = C ₃ , PRM _V = PR ₁ -C ₃ (criterion 51)

· Cy = C-PRM_X = 0 , PRM_r = PRM-C , C_z, =

= C₃ · Cy = C-PRM _X = 0, PRM _r = PRM-C, C _z , =

= C ₃

(Cas identique au deuxième choix) (Same case as second choice)

• C₃, = C₃ PRM_l (Cas identique au premier choix)

• C ₃ , = C ₃ PRM _l (Same case as first choice)

On remarque que cette solution aurait pu être directement accessible en n'utilisant que les contraintes sur les PRM et les DOA. Comme le mode 1 (DOA1) est identique à la classe C3, la solution reste inchangée. La solution finale (cf. annexe 15) est donc une partition en trois classes . C₁={l₁₂,I₁₃,I₁₄,I₁₅,I₁₆,I₁₇,I_li,I₁₉,I₂₀}, C₂={l₆,I₇,I_s,I₉,I ,I_n) etNote that this solution could have been directly accessible using only the constraints on PRM and DOA. Since mode 1 (DOA1) is identical to class C3, the solution remains unchanged. The final solution (see Appendix 15) is therefore a partition into three classes. C ₁ = {l _12, I _13, I _14, I _15, I _16, I _17, I _li, I _19, I _20}, C ₂ = {l _6, I _7, I _s, I _9, I, I _n ) and

^3^— ΙΛ ' ^2 ' ^3 ' ^4 ' ^5 }^■ En référence aux annexes 1 6 à 1 9 et aux figures 9 et 1 0 un nouvel exemple de mode de mise en œuvre de l'étape de classement non supervisé va être présenté. Ce classement repose sur une classification binaire croisée des modes. Cette méthode comprend principalement trois étapes, une étape Etp 56 de classification des vecteurs impulsion en fonction de leur lien de corrélation, une étape Etp 57 de classification des vecteurs mode en fonction de leur lien de corrélation, une étape Etp 58 de regroupement des modes et impulsions en classes. L'annexe 1 6 représente un tableau de n impulsions enrichies non classées, notées à l _N, comportant différents modes. ^ 3 ^- ΙΛ '^ 2' ^ 3 '^ 4' ^ 5} ^■ Referring to Appendices 16 to 19 and Figures 9 and 10, a new example of how to implement the unsupervised grading step will be presented. This ranking is based on a crossed binary classification of modes. This method mainly comprises three steps, a step Etp 56 of classification of the pulse vectors according to their correlation link, a step Etp 57 of classification of the mode vectors according to their correlation link, a step Etp 58 of grouping of the modes and pulses in classes. Annex 1 6 represents a table of n unclassified enriched pulses, denoted at l _N, comprising different modes.

L'annexe 1 7 représente un tableau des liens de corrélation entre vecteurs impulsions indiquant les liens de corrélation entre chaque impulsion. On peut organiser ces liens de corrélation de la plus forte corrélation à la plus faible corrélation. Par exemple, au niveau 1 , on peut rapprocher l₄ et l₆ ou h , l₂, et l₅ car ils ont le même lien de corrélation. De même, au niveau 0,4 les impulsions l₃, l₄ et l₆ se rejoignent. Après avoir fixé un seuil à une valeur S-i , on rapproche les impulsions dont le lien de corrélation est inférieure à un ce seuil de de corrélation S-i . On peut donc regrouper les impulsions {h , l₂, l₅}, {I 7 , I N} et {(I4, I6) , (13)}. Ainsi, dans le tableau annexe 1 6 lorsque l'on rapproche les vecteurs impulsions suivant les liens de corrélation, on obtient le tableau reclassé présenté annexe 1 8. Annex 17 represents a table of the correlation links between impulse vectors indicating the correlation links between each pulse. These correlation links can be organized from the highest correlation to the lowest correlation. For example, at level 1, we can compare l ₄ and l ₆ or h, l ₂ , and l ₅ because they have the same correlation link. Similarly, at level 0.4 pulses l ₃ , l ₄ and l ₆ meet. After setting a threshold to a value Si, we compare the pulses whose correlation link is less than this threshold Si correlation. The pulses {h, l ₂ , l ₅ }, {I 7, IN} and {(I4, I6), (13)} can thus be grouped together. Thus, in the appendix table 1 6 when we compare the impulse vectors according to the correlation links, we obtain the reclassified table presented appendix 1 8.

Comme le classement utilise une méthode de classification croisée, l'étape de classement utilise également les liens de corrélation entre modes. L'annexe 1 9 présente les liens de corrélation entre les différents modes. On rappelle que les liens de corrélation d entre deux impulsions ou entre deux modes peuvent être calculées à l'aide de la formule : Since the classification uses a cross-classification method, the ranking step also uses the correlation links between modes. Annex 1 9 presents the correlation links between the different modes. Recall that the correlation links d between two pulses or between two modes can be calculated using the formula:

di_{R R =} card {R_X U R₂ ) - card {R_X n R₂ ) Comme précédemment ces liens de corrélation entre modes peuvent être représentés sous la forme d'une structure pyramidale présentée figure 1 0. Cette structure fait apparaître des groupements de modes. Après avoir fixé un deuxième seuil de corrélation à une valeur S₂, on rapproche les modes dont le liens de corrélation est inférieur à un ce seuil de corrélation S₂. Dans notre exemple, on peut donc regrouper les modes {AOA _; PRM-i , F₂, IPo, F_N}, {AOA₂, F_{1 ;} IP₂, PRM_P} et {AOA_m, IP_{1 ;} PRM₂}. On obtient ainsi le tableau reclassé annexe 20. Dans ce tableau réorganisé on voit apparaître des regroupements d'intérêt (Cl-i , Cl₂, Cl₃) correspondant aux différentes classes. Les groupements d'intérêt peuvent se lire dans les blocs de 1 de la matrice réorganisée. di _{RR =} card {R _X UR ₂ ) - card {R _X n R ₂ ) As previously, these correlation links between modes can be represented in the form of a pyramidal structure presented in FIG. 1 0. This structure shows groups of modes. After fixing a second correlation threshold at a value S ₂ , the modes whose correlation links are less than a correlation threshold are compared. S ₂ . In our example, we can group modes {AOA _; PRM-1, F ₂ , IPo, F _N , AOA ₂ , F _1; IP ₂ , PRM _P } and {AOA _m , IP _1; PRM ₂ }. The table reclassified in Appendix 20 is thus obtained. In this reorganized table, clusters of interest (Cl-i, Cl ₂ , Cl ₃ ) corresponding to the various classes appear. Groups of interest can be read in the 1-block of the reorganized matrix.

La solution obtenue dépendra des seuils de corrélation choisis. En choisissant un seuil de 1 , on obtient une solution dans laquelle tous les modes et/ou toutes les impulsions sont regroupées, on a une confusion maximale. A l'inverse en choisissant un seuil de 0, on obtient un éclatement maximal des modes et des impulsions. The solution obtained will depend on the chosen correlation thresholds. By choosing a threshold of 1, we obtain a solution in which all the modes and / or all the pulses are grouped, we have a maximum confusion. Conversely, by choosing a threshold of 0, a maximum burst of the modes and the pulses is obtained.

De façon avantageuse, cette implémentation est beaucoup plus rapide que la méthode présentée précédemment car elle ne manipule que des valeurs binaires. Cette implémentation peut donc être projetée sur des matériels travaillant en binaire et des composants de calcul extrêmement élémentaires comme par exemple des réseaux logiques programmables ou FPGA pour " Field Programmable Gâte Arrays " selon l'acronyme anglo saxon qui vont traiter l'algorithme. De plus ces composants vont pouvoir paralléliser les calculs réduisant ainsi les temps de traitement. Advantageously, this implementation is much faster than the method presented above because it handles only binary values. This implementation can therefore be projected on hardware working in binary and extremely basic computing components such as programmable logic networks or FPGA for "Field Programmable Gaste Arrays" according to the acronym for the Anglo-Saxon that will process the algorithm. Moreover, these components will be able to parallelize the calculations thus reducing the processing times.

Dans l'algorithme de recherche du classement optimal on est amené dans différentes étapes à calculer les similarités ou les liens de corrélation de toutes les classes 2 à 2 soit pour un couple de classes C : ieC_{ ,i'eC_q, In the search algorithm of the optimal classification we are led in different steps to calculate the similarities or the correlation links of all classes 2 to 2, ie for a pair of classes C: ieC _{ , i'eC _q ,

Si les 2 classes comportent respectivement n et m impulsions, alors la complexité de ce calcul sera d'ordre O(nxm). Cette complexité peut être évitée. En effet : Q =∑∑-v/.v/ =∑∑.v,'.v/ If the two classes have respectively n and m pulses, then the complexity of this calculation will be of order O (nxm). This complexity can be avoided. Indeed : Q = ΣΣ-v / .v / = ΣΣ.v, '. V /

ί^' ί" j=l j=l i,ï ί ^' ί "j = lj = li, ï

On obtient alors un calcul d'ordre O(n+m) où n et m représentent le nombre d'impulsions des deux classes. On a ainsi transformé un calcul qui était quadratique en un calcul de complexité linéaire. We then obtain a computation of order O (n + m) where n and m represent the number of pulses of the two classes. We have thus transformed a calculation that was quadratic into a calculation of linear complexity.

Si on considère le critère de Condorcet utilisé précédemment, ce critère s'écrit de la façon suivante : c( ) =∑∑(cA + ¾.) If we consider the Condorcet criterion previously used, this criterion is written in the following way: c () = ΣΣ (cA + ¾.)

i=l i'=l i = l i '= l

Ce calcul est apparemment quadratique. Il peut être ramené en linéaire grâce au principe pré-établi : c(x)=∑∑c_iVx_ir +∑f c_irx_iV This calculation is apparently quadratic. It can be brought back to linear thanks to the pre-established principle: c (x) = ΣΣc _iV x _ir + Σfc _ir x _iV

i=l i'=l i=l i'=l i = l i '= l i = l i' = l

∑C_ir = R'.Rⁱ et C_i{ = Minic_u,C_Vi)-C_1{ il suit : C(x)=∑R^l.R^l -∑R^l.R^q +∑Μίη{& .R^l ,R^q .R^q ) ΣC _ir = R'.R ⁱ and C _{i {} = Minic _u , C _Vi ) -C _{1 {} it follows: C (x) = ΣR ^l .R ^l -ΣR ^l .R ^q + ΣΜίη {& .R ^l , R ^q .R ^q )

l 1,1*1 1,1*1 l 1,1 * 1 1,1 * 1

On obtient ainsi une expression calculable en O(n). This gives a calculable expression in O (n).

Dans la présente demande, des méthodes de traitement à complexité linéaire ont été décrites. Bien entendu, ces méthodes ne sont pas les seules à pouvoir être utilisées et le procédé de désentrelacement selon l'invention peut également être réalisé en utilisant des méthodes à complexité non linéaires. PRM principal PRM secondaire fréquence groupe f groupe PRM phase In the present application, linear complexity processing methods have been described. Of course, these methods are not the only ones that can be used and the deinterlacing process according to the invention can also be carried out using nonlinear complexity methods. PRM main PRM secondary frequency group f group PRM phase

imp 1 f1 1 1 P1 imp 1 f1 1 1 P1

imp i fi 1 1 Pi imp i fi 1 1 Pi

imp j fj 2 4 pj imp n fn 2 4 pn imp j fj 2 4 pj imp n fn 2 4 pn

Annexe 1 Annex 1

Annexe 2 Annex 2

Annexe 3 Annex 3

Classe 1 Classe 2 Classe 3Class 1 Class 2 Class 3

Impulsions enrichies Enriched impulses

DNS

DNS

Tableau impulsions-modes Tableau impulsions-modes réorganisé Table pulses-modes Table pulses-modes reorganized

Annexe 4 Annex 4

Annexe 5a Annex 5a

i=l i = l

Annexe 5b 11 I2 I3 I4 I5 I6 I7 I8 I9 110 11 1 112 113 114 115 11 6 11 7 118 119 I20Annex 5b 11 I2 I3 I4 I5 I6 I7 I8 I9 110 11 1 112 113 114 115 11 6 11 7 118 119 I20

F1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0F1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

F2 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0F2 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0

F3 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 1 1 1F3 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 1 1 1

F4 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0F4 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0

DOA1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0DOA1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

DOA2 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0DOA2 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0

DOA3 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1DOA3 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1

PRM1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0PRM1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

PRM2 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0PRM2 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0

PRM3 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 PRM3 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1

Annexe 6a Annex 6a

Critère Criterion

Annexe 7 Annex 7

Annexe 8 Annex 8

Critère 110 111 112 113 114 115 116 117 IIS 119 120

Criterion 110 111 112 113 114 115 116 117 IIS 119 120

Annexe 9

Annex 9

Critère Criterion

mode 10 fashion 10

Critère Criterion

mode 10

fashion 10

Critère 1S 19 112 113 114 ! 1S ne 17 1 18 119 20 Criterion 1S 19 112 113 114! 1S no 17 1 18 119 20

231 mode 1 0 0 0 0 0 0 0 0 0 0 0 DOA1 231 mode 1 0 0 0 0 0 0 0 0 0 0 0 DOA1

Annexe 10a Annex 10a

Critère 19 110 111 112 113 114 115 116 117 118 119 120 Criterion 19 110 111 112 113 114 115 116 117 118 119 120

mode 1 0 0 0 0 0 0 0 0 0 0 0 0 D0A1 mode 2 0 0 0 0 0 0 0 0 0 0 0 0 PR 1 mode 1 0 0 0 0 0 0 0 0 0 0 0 0 D0A1 mode 2 0 0 0 0 0 0 0 0 0 0 0 0 PR 1

Annexe 10b Critère 1S 19 112 113 114 15 116 117 ! 1S 19 !20 Annex 10b Criterion 1S 19 112 113 114 15 116 117! 1S 19! 20

307 mode 1 0 0 0 0 0 0 0 0 0 0 0 DOA1 mode 2 0 0 0 0 0 0 0 0 0 0 0 PRM 1 mode 3 0 0 0 0 0 0 0 0 0 0 Q F l 307 mode 1 0 0 0 0 0 0 0 0 0 0 0 DOA1 mode 2 0 0 0 0 0 0 0 0 0 0 0 PRM 1 mode 3 0 0 0 0 0 0 0 0 0 0 Q F l

Annexe 1 1 Annex 1 1

PRM3 moce 10 PRM3 MOCE 10

Annexe 12a Annex 12a

½» ½}

½ "½}

Annexe 12b Cr

Annex 12b Cr

mode 10 0 0 0 0 0 0 0 0 0 PP. M 3

mode 10 0 0 0 0 0 0 0 0 0 PP. M 3

Annexe 13

Annex 13

Critère

Criterion

mode 5 0 0 0 0 DOA2 mode 6 0 0 0 0 PR M2 mode 7 0 0 0 0 F3 mode 8 0 0 0 0 F4 mode 9 0 0 0 0 DOA3 mode 10 0 0 0 0 PR M3 mode 5 0 0 0 0 DOA2 mode 6 0 0 0 0 PR M2 mode 7 0 0 0 0 F3 mode 8 0 0 0 0 F4 mode 9 0 0 0 0 DOA3 mode 10 0 0 0 0 PR M3

Annexe 14 Annex 14

Annexe 16

Annex 16

Annexe 17

Annex 17

Annexe 18

Annex 18

Annexe 19

Annex 19

Claims

REVENDICATIONS 1 . Procédé de désentrelacement non supervisé de trains d'impulsions comprenant : CLAIMS 1. A method of unsupervised deinterleaving of pulse trains comprising:

- une étape (Etp10) de constitution, pour chaque paramètre ou groupe de paramètres, de N histogrammes, N étant un entier supérieur à 1 , à une ou plusieurs dimensions définissant des groupes d'impulsions par paramètre ou groupe de paramètres, lesdits histogrammes étant représentatifs des trains d'impulsions; a step (Etp10) of constitution, for each parameter or group of parameters, of N histograms, N being an integer greater than 1, with one or more dimensions defining groups of pulses per parameter or group of parameters, said histograms being representative of pulse trains;

- une étape (Etp20) d'extraction de modes à partir desdits histogrammes ; ledit procédé comprenant en outre : a step (Etp20) for extracting modes from said histograms; said method further comprising:

- une étape (Etp30) de constitution de groupes d'intérêt à partir desdits modes extraits ; a step (Etp30) for constituting groups of interest from said extracted modes;

- une étape (Etp40) d'enrichissement de la description des impulsions à l'aide de labels traduisant leur appartenance à au moins un groupe d'intérêt ; a step (Etp40) of enriching the description of the pulses using labels reflecting their membership in at least one interest group;

- une étape (Etp50) de classement croisé non supervisé des modes afin de réorganiser les impulsions enrichies par similarité. a step (Etp50) of unsupervised cross-classification of the modes in order to rearrange the enriched pulses by similarity.

2. Procédé selon la revendication précédente caractérisé en ce que l'étape (Etp20) comprend : 2. Method according to the preceding claim characterized in that the step (Etp20) comprises:

- une étape (Etp21 ) de seuillage des histogrammes ; a step (Etp21) for thresholding the histograms;

- une étape (Etp22) de décomposition des histogrammes en modes ; a step (Etp22) of decomposing the histograms into modes;

- une étape (Etp23) de caractérisation des modes ; a step (Etp23) of characterization of the modes;

3. Procédé selon une des revendications 1 à 2 caractérisé en ce que l'étape (Etp50) de classement des modes consiste à créer progressivement au moins une classe d'impulsions-modes homogènes, l'étape comprenant : - à l'initialisation une recherche du mode ayant un maximum d'intérêt, ce mode formant la première classe, puis, pour chaque mode restant : 3. Method according to one of claims 1 to 2 characterized in that the step (Etp50) mode classification is to gradually create at least one class of homogeneous mode pulses, the step comprising: at the initialization, a search of the mode having a maximum of interest, this mode forming the first class, then, for each remaining mode:

- une étape (Etp51 ) de recherche, parmi les modes restants, du mode ayant un cardinal d'intersection avec la classe en cours le plus grand ; a step (Etp51) of searching, among the remaining modes, of the mode having a cardinal intersection with the largest class in progress;

- une étape (Etp52) de recherche de la partition du mode trouvé, permettant de maximiser un critère de regroupement avec la ou les classes en cours ; a step (Etp52) of searching for the found mode partition, making it possible to maximize a grouping criterion with the current class or classes;

- une étape (Etp53) de fusion du mode avec la ou les classes en cours. a step (Etp53) of merging the mode with the current class or classes.

4. Procédé selon une des revendications 1 à 2 caractérisé en ce que l'étape (Etp50) de classement comprend : 4. Method according to one of claims 1 to 2 characterized in that the step (Etp50) classification comprises:

- une étape (Etp56) de classification des impulsions en fonction de leurs liens de corrélation mutuelles ; a step (Etp56) of classification of the pulses according to their mutual correlation links;

- une étape (Etp57) de classification des modes en fonction de leurs liens de corrélation mutuelles ; a step (Etp57) of classification of the modes according to their mutual correlation links;

- une étape (Etp58) de regroupement des modes et impulsions en classes. a step (Etp58) of grouping modes and pulses into classes.