Brevet FR3011338A1 — publié le 20150403

PROCÉDÉ DE DÉSENTRELACEMENT NON SUPERVISÉ PAR ENRICHISSEMENT N DIMENSIONS La présente invention concerne le domaine des systèmes d'écoute, et 5 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. Le traitement considéré est celui du désentrelacement non supervisé 10 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 15 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 ". Ces situations de mélange ou entrelacement de trains d'impulsions sont d'autant plus fréquentes que la sensibilité des capteurs MRE augmente 20 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. Pour traiter ces situations l'extracteur doit en tout premier lieu isoler 25 chaque émission d'intérêt pour pouvoir ensuite l'analyser et la caractériser. On parle de " Désentrelacement ». Les émissions radars interceptées sont grossièrement de trois types : - TE1 : Cela correspond aux émissions de Haute et Moyenne 30 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. - 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. Par la suite, on s'intéressera au désentrelacement des émissions de 15 type TE2 et TE3. Les émissions de type TE1 sont traitées séparément par des algorithmes rapides. On distingue deux fonctions de désentrelacement complémentaires: - Le Désentrelacement Supervisé (DS), dans lequel on utilise les 20 caractéristiques des signaux connues comme par exemple, les Fréquences (F), les PRI, les Durées d'Impulsion (Dl) ou les paramètres intra impulsion (ou IP pour IntraPulse) pour reconnaltre les impulsions provenant de ces signaux. - Le Désentrelacement Non Supervisé (DNS), dans lequel on n'a 25 pas connaissance de ces caractéristiques. Dans la plupart des extracteurs ESM ou ELINT ces deux fonctions coopèrent comme illustré figure 2. Les méthodes de DNS peuvent être segmentées entre : 30 - Celles qui utilisent des classifications (ou clustering) sur paramètres primaires de l'impulsion (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. 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 5 simultanément, 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 10 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 ! - Les capteurs ELINT ont des sensibilités de plus en plus élevées 15 pour la détection de certaines menaces. - Les émissions sont agiles ce qui nécessite, pour capter au maximum la forme d'onde : - d'avoir des bandes de réception importantes, - d'avoir des temps d'écoute longs. 20 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. La canalisation de certains récepteurs impose de balayer la bande de 25 réception utile, ne permettant l'acquisition de la forme d'onde que partiellement. 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. 30 En résumé, on constate que : - les volumes de données augmentent ; il serait utile d'avoir des analyses sur des temps d'écoute longs ; Il faut traiter les paramètres dans une analyse conjointe. Les algorithmes actuels ne permettent pas de remplir ces différentes contraintes. Si la complexité de l'algorithme n'est que quadratique cela signifie en 5 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 apparaltre des 10 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 15 qu'ils ne peuvent pas l'être en temps court. 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 croc 20 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. 25 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, LI, IP) ou la séparation par histogramme de différence de temps d'arrivée (ou HDTOA pour Histogram of Difference Time Of Arrivai selon la terminologie anglo saxonne). 30 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 ". 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 ". 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. 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. La généralisation de l'agilité impulsion à impulsion et de l'agilité par 35 trains des formes d'onde rend aujourd'hui difficile voire caduque cette démarche. En effet dans un environnement dense il apparat 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 : - d'une part l'utilisation du DTOA (pour Difference of Time Of Arrivai). - 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 retrouve les deux étapes précédentes : - 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. - Une étape de sélection des impulsions correspondant à une PRI, que nous appellerons " marquage " (ou « gating " selon le terme anglo saxon). 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 : - 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. - 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). - 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. - L'histogramme des doublets qui permet de rechercher les PRM avec une grande efficacité. On montre qu'il correspond à une construction 15 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. - L'Histogramme de Fourier. Les techniques fondées sur la construction d'histogrammes de DTOA ne sont pas robustes aux PRI jittées, 20 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 25 algorithmes de ce type ont été décrits dans la littérature. Ces histogrammes ont différentes propriétés. Tous ces histogrammes peuvent avoir une grande robustesse au mitage et au 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 30 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. 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. 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. 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, Ip...). Un but de l'invention est notamment de corriger tout ou partie des 30 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. A cet effet, l'invention a pour objet un procédé de désentrelacement non supervisé de trains d'impulsions caractérisé en ce qu'il comprend une 35 étape d'enrichissement N dimensions, N étant un entier supérieur à 1. Selon un mode de mise en oeuvre, le procédé comprend: - une étape Etp10 de constitution d'histogrammes à plusieurs dimensions représentatifs des trains d'impulsions; 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 ; - une étape Etp50 de classement croisé non supervisé des modes. Selon un mode de mise en oeuvre, l'étape Etp20 comprend : - une étape Etp21 de seuillage des histogrammes ; io une étape Etp22 de décomposition des histogrammes en modes ; une étape Etp23 de caractérisation des modes ; Selon un mode de mise en oeuvre, l'étape Etp50 de classement des modes consiste à créer progressivement au moins une classe d'impulsions-modes homogènes, l'étape comprenant : 15 - à 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 : - une étape Etp51 de recherche, parmi les modes restants, du mode ayant un cardinal d'intersection avec la classe en cours le plus grand ; - une étape Etp52 de recherche de la partition du mode trouvé, 20 permettant de maximiser un critère de regroupement avec la ou les classes en cours ; - une étape Etp53 de fusion du mode avec la ou les classes en cours. Selon un mode de mise en oeuvre, l'étape Etp50 de classement comprend : 25 - une étape Etp56 de classification des impulsions en fonction de leurs liens de corrélation mutuelles ; - une étape Etp57 de classification des modes en fonction de leurs liens de corrélation mutuelles ; - une étape Etp58 de regroupement des modes et impulsions en classes. De façon avantageuse, le procédé selon l'invention, grace à sa faible 5 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. D'autres particularités et avantages de la présente invention 10 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 : - La figure 1 illustre le désentrelacement ; 15 - La figure 2 illustre la coopération entre les désentrelacements supervisé et non supervisé ; - La figure 3 représente un exemple de désentrelacement non supervisé utilisant des classifications RF-DOA suivi par des classifications HDTOA ; 20 - La figure 4 représente un exemple de mise en oeuvre du procédé selon l'invention ; - La figure 5 représente un exemple de représentation graphique d'un histogramme RF-DTOA ; - La figure 6 représente un exemple de mise en oeuvre de l'étape 25 d'extraction de modes ; - Les figures 7a et 7b illustrent l'étape d'extraction des modes d'un histogramme ; - La figure 8 représente différents choix de classification possible ; - La figure 9 représente un exemple de structure pyramidale 30 illustrant les liens de corrélation entre impulsions ; - La figure 10 représente un exemple de structure pyramidale illustrant les liens de corrélation entre modes. Le procédé de désentrelacement non supervisé par enrichissement N dimensions selon l'invention repose sur l'utilisation conjointe des paramètres 5 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. 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, 10 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. Au cours de l'étape Etp10 des histogrammes sont constitués dans 15 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, LI....) d'un ensemble d'impulsions. Cette description conjointe permet en particulier de construire des histogrammes liant paramètres primaires, comme par exemple les 20 fréquences, les DOA, les largeurs d'impulsion (LI),... 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. 25 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 30 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 35 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. Il est également aisé de reconnaUe 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 "ité Bien entendu les histogrammes ne se limitent pas à ceux en deux 10 dimensions (2D) et sur la représentation graphique on aurait pu ajouter Azimut (en z), LI en(v)... On parle d'histogrammes à N Dimensions ou 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 15 Etp22 de découpe en modes et une étape Etp23 de caractérisation des modes. Pour le seuillage des histogrammes, on peut distinguer deux cas, le cas à une dimension (1D) et le cas à deux dimensions (2D) dans lequel on 20 utilise systématiquement la DOA. 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. 25 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 : 30 = 2(HV - HE )2 2 X HE 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%, É05,1 = 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 5 présentent pas de mélange. Quand deux pics sont proches on décide que les données appartiennent à la fois aux deux pics. 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 10 maxima par exemple par une méthode de segmentation par ligne de partage des eaux. Les modes de ces histogrammes définissent des groupes (ou clusters selon la terminologie anglo saxonne) d'impulsions par paramètre. Chaque 15 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 20 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.... Chaque impulsion est ensuite enrichie au cours d'une étape Etp40 25 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. Les impulsions enrichies peuvent se présenter sous la forme d'une 30 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 35 appartenant au mode comme illustré annexe 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. 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. La méthode repose sur la dualité individus/variables. Les impulsions 10 sont ici les individus et les différents modes (DTOA, Fréquences, DOA, LI ...) 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. 15 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. On distingue le lien de corrélation intra dimension, par exemple entre 20 deux fréquences ou deux PRM, et le lien de corrélation inter dimension, par exemple entre une fréquence et une PRM. 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: 25 d(RpR2)= card(R,o R2)- card (R,n R2 ) card(R, U R2) 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 30 fréquence peut être réalisée en considérant la cardinalité de l'intersection entre les listes d'impulsions de la PRM et de la Fréquence considérées. card(RF n Rj,pRA,I) 42,F, R j,pRm) 1 card(R,,F) 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. 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. 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...). 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. On défini donc les critères de similarité suivants : Dans une même classe: - Les impulsions partagent un maximum de composantes ; - Les composantes partagent un maximum d'impulsions. Entre deux classes: - Les impulsions partagent un minimum de composantes ; - Les composantes partagent un minimum d'impulsions. De façon à accélérer le classement on dispose également de contraintes. A titre d'exemple, ces contraintes peuvent être : Deux composantes DOA (idem pour deux composantes PRM) « disjointes " ne peuvent appartenir à une même classe ; Les fréquences d'une même « phase de PRM » ne doivent pas être séparées ; 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 10 deux impulsions peut être donné par la formule : Cn =IXj Xj, j=1 Où xi représente un indicateur binaire indiquant si l'impulsion i appartient au mode j défini par : xi =1 si l'impulsion i appartient au mode j 15 =0 si l'impulsion i n'appartient pas au mode j Et où P représente le nombre total de modes. De même, de façon symétrique, on peut définir un critère de similarité entre modes par la formule : 20 C"' =IX' Xi i=1 Où x, représente un indicateur binaire indiquant si le mode j contient l'impulsion i défini par : x =1 si le mode j contient l'impulsion i xi =0 si l'impulsion i n'appartient pas au mode j 25 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 12 ainsi que le critère pour les 30 impulsions 13 et 14. C11,12=4 signifie que les impulsions et 12 sont toutes les deux présentes dans 4 modes (cf. cases entourées). De même, l'annexe 5b présente le critère de similarité entre les modes A0A1 et PRM1 ainsi que celui entre les modes AOA2 et 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: 1 si i et i' sont dans la même classe Xll, 0 sinon 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 : N N C(X)=EZ(C i=1 i'-=1 Où: N représente le nombre d'impulsions, Cu, représente le critère de similarité entre les impulsions In et In' représente l'élément de rang i,i' de la matrice carré X précédemment définie. Les tableaux annexe 6b et 6c présentent des exemples de partition en 20 différentes classes d'une classification de 20 impulsions et 10 modes présentée annexe 6a. L'annexe 6b représente une classification des impulsions et des modes en 3 classes C11, Cl2 et C13. Dans l'exemple les classes 3 et 2 se partagent la fréquence F2, et les classes Cl2 et CI1 se partagent les 25 fréquences F3 et F4. Pour cette partition, le critère est de 1076. 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 Cl11, CI12, C113, Cl14 et CI15. Cette partition n'est pas optimale, la valeur du critère pour cette partition est 912 donc inférieure à 1076. 30 Un premier exemple de mode de mis en oeuvre 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 restant,le classement comprend : - une étape Etp51 de recherche, parmi les modes restants, du mode ayant un cardinal d'intersection avec la classe en cours le plus grand ; - une étape Etp52 de recherche de la partition du mode trouvé, permettant de maximiser un critère de regroupement avec la ou les classes 10 en cours ; - une étape Etp53 de fusion du mode avec la ou les classes en cours. 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 15 critère est la plus forte. Ce mode formera la première classe C1. On rappelle que la contribution au critère du mode q est donnée par la formule : Contrib(q)= C', -a Emire ', ieCq,YeCq ieCq,i'eCq Où: représente la similarité ou le lien de corrélation entre 20 modes i et i', a représente un scalaire que l'on peut prendre égale à 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. 25 Suivant un premier mode de mise en oeuvre du procédé on peut classer les modes restant par ordre décroissant des cardinaux d'intersection. Suivant un autre mode de mise en oeuvre, on recherche parmi les modes restant, celui dont le cardinal d'intersection est le plus grand. Ce mode formera la classe C2. 30 En référence à la figure 8, si on considère les classes C1 et C2, (Cx=C1 et Cy=C2), il existe quatre choix possible pour le regroupement de ces classes. L'intersection c1nc2 peut appartenir à la nouvelle classe C'2 (choix 1), à la nouvelle classe C'i (choix 2) ou former une nouvelle classe C'3 (choix 3). Une quatrième option consiste à considérer un nouvelle classe C'i regroupant les classes Ci et C2 (choix 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 5 maximise ledit critère de regroupement. On continue ainsi, pour chaque mode restant, en considérant à chaque fois l'intersection entre la ou les classes (C't C3 ) que l'on vient de former et le mode ayant le cardinal d'intersection de plus grand avec la ou 10 les nouvelles classes. Si le choix 1 est retenu, les nouvelles classes formées sont: C'1= C1' C2 et C'2= C2 ; Si le choix 2 est retenu, les nouvelles classes deviennent : C'1= C1 et C'2= C2- C1; 15 Si le choix 3 est retenu, les nouvelles classes deviennent : C'1= C1, C'2= C2 et C'3= Ci nC2; Si le choix 4 est retenu, la nouvelle classe devient : C'1= Ci UC2. On ne s'intéresse qu'aux intersections non nulles. Dans le cas où le 20 mode considéré aurait une intersection nulle, il formerait une nouvelle classe. 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. Il est à noter que l'on classe plusieurs données à la fois. En effet 25 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 li à 120. 30 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,). Après classement des modes par ordre croissant de leur contribution on obtient le tableau annexe 7. 35 On commence le classement non supervisé avec le mode 10 (PRM3). Ce mode contient les impulsions 112 à 120 qui vont former la première classe C1. On considère ensuite le mode dont le cardinal d'intersection avec la classe Cl 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'nCy =CzuCi, =Cx= Cy , les 4 choix sont donc identiques. Les modes 10 et 9 forment donc une seule classe C1= {112,113,114,115,116,117,118, /19 ,120 } (cf. annexe 8). On considère ensuite le mode 8 (F4) qui est le mode ayant la plus grande intersection avec la classe Cl (C1nF4 =5) (cf. annexe 9). On teste les 4 choix suivants: - Cl, = C1-F4, F4, =F4 (critère 243) - Cr = F4-C1 (critère 303) - Cl, = F4, = F4 -C1 , Cz, = C1nF4 (critère 263) - C1, = CmF4 (critère 235) Le critère le plus élevé étant 303, on retient la solution consistant à séparer les modes en deux classes, à savoir C1= 1/12, /13, /14, /15 , /16, /17, /18,4, /20 / et C2= {110,111}- 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 Cl (C1nF3 =3) Comme précédemment on teste les 4 choix suivant: - Cl, = F3, =F3 (critère 251) - C1, = C1 , F3. = F3-C1 (critère 307) - C1, =Cl-F3F3, = F3-C1, Cz, = C1nF3 (critère 267) - Cl, = CluF3 (critère 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 : Ci= 1/12, /i3, /14, /1' /16 /17 , /18 c2={110,I11} et C3= 1/8, /91. On peut remarquer qu'en utilisant les contraintes en PRM ou en DOA (cf. annexe 11), à savoir « deux composantes DOA (ou deux composantes 35 PRM) disjointes ne peuvent appartenir à une même classe », aurait permis de 301 13 3 8 21 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. 5 On poursuit (cf. annexe 12a) en considérant le mode 6 (C2nPRM2 = 2) et on teste les 4 choix suivants : - C2, = C2 -PRM2 = , PRM2, =PRM2 (critère 84) - C2, = C2 , PRM2, =PRM2 - C2 (critère 68) 10 - C2, = C2 -PRM2 = J , PRM2, =PRM2 - C2 , Cz, = C2nPRM2 = C2 (identique au choix précédent) - C2, = C2 UPRM2 = PRM2 (identique au premier choix) A ce stade on a 3 classes : C1= {[12 ,1139114,115111691179118,119,120}, C2 = {Io 110 111} et C3=1/8,/91. On continue en testant l'intersection 15 du mode 6 avec la classe C3 (C3nPRM2 = 2). La résolution est identique, on n'a que deux choix. A l'issue de cette étape (cf. annexe 12b) on a donc 2 classes à savoir: ci= 1112,1-13,114,115,116,117,118,119,1201 et c2= {/8,4,4,4}. Les anciennes classes C2 et C3 qui étaient séparées se retrouvent maintenant réunies 20 ensemble. On continue en considérant le mode 5 (DOA2). Comme DOA2=C2 la solution reste identique. 25 On poursuit (cf. annexe 13) en considérant le mode 4 (F2). L'intersection entre le mode 4 et la classe 1 étant nulle (C1nF2= 0), il n'y a aucun test à effectuer. C2nF2 = 2 , on teste les 4 choix : - C2, = C2 -F2 , F2, =F2 (critère 179) - C2, = C2 , F2, = F2 -C2 (critères 207) 30 - C2, = C2 -F2 , F2, = F2 -C2 , Cz, = C2 nF2 (critère 191) - C2, = C2 UF2 (critère 123) A l'issue de cette étape, une nouvelle classe C3. On a les classes : cl= {'12,113,114,115,'16,'17,'18,'19'12o}' C2=1/6,4/8,4,4,41 et C3= {4/4, 15 On considère maintenant le mode 3 (F1). Ce mode n'intervient pas sur les classes Ci, 02 et C3 (C1nF1 = 0 C2nF1= 0 et C3nF1 = 0). On poursuit (cf. annexe 14) avec le mode 2 (PRM1). Ce mode 5 n'intervient ni sur la classe Ci ni sur la classe C3 (C1nPRM1 = 0 et C2nPRM1 = 0). Le mode 2 a une intersection non nulle avec le classe C3 (C3nPRM1 = 3), on teste donc les 4 choix : - C3, = C3-PRMI = QS , PRM1, =PRM1 (critère 63) - C3, = C3 , PRM1. = PRM1-C3 (critère 51) 10 - Cy = C3-PRM1 = 0 , PRM1, = PRM,-C3, Cz, = C3nPRM1 = C3 (Cas identique au deuxième choix) - C3, = C3L.)PRIld1 (Cas identique au premier choix) On retient le premier choix, on a donc 3 classes : C1= 1/12, 113 ,114 , 115 , 116 , 118 ,119 120 , C2= {16,17,18,19,110,111} et 15 C3-= tt1,/2,/3,/4, /51. On remarque que cette solution aurait pu être directement accessible en n'utilisant que les contraintes sur les PRM et les DOA. 20 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 : Ci= {112,113,114,115,116,117,118,119,120}, C2 =116,418,19,110,1j et C3= {/1,12,13,14"5}- 25 En référence aux annexes 16 à 19 et aux figures 9 et 10 un nouvel exemple de mode de mise en oeuvre 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 30 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 16 représente un tableau de n impulsions enrichies non classées, notées à IN, comportant différents modes. L'annexe 17 représente un tableau des liens de corrélation entre 5 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 14 et 16 ou 12, et 15 car ils ont le même lien de corrélation. De même, au niveau 0,4 les impulsions 13, 14 et 16 se rejoignent. Après avoir fixé un seuil à une valeur 10 St on rapproche les impulsions dont le lien de corrélation est inférieure à un ce seuil de de corrélation S1. On peut donc regrouper les impulsions {li, 12, 15}, {17, IN} et {(14,16), (13)}. Ainsi, dans le tableau annexe 16 lorsque l'on rapproche les vecteurs impulsions suivant les liens de corrélation, on obtient le tableau reclassé présenté annexe 18. 15 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 19 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 : 20 d(Ri, R2) ,card(Ri v R2 )- card (Ri n R2 ) card(Ri u R2 ) 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 10. Cette structure fait apparaitre des groupements de modes. Après avoir fixé un deuxième seuil de corrélation à une valeur S2, on rapproche les 25 modes dont le liens de corrélation est inférieur à un ce seuil de corrélation S2. Dans notre exemple, on peut donc regrouper les modes {AOA1, PRM1, F2, IPo, FN}, {AOA2, F1, 1P2, PRMp} et {AOAm, IP1, PRM2}. On obtient ainsi le tableau reclassé annexe 20. Dans ce tableau réorganisé on voit apparaitre des regroupements d'intérêt (Cli, C12, C13) correspondant aux différentes 30 classes. Les groupements d'intérêt peuvent se lire dans les blocs de 1 de la matrice réorganisée. 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. 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 Gate 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. 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 : ieC1 ,YE Cg, Si les 2 classes comportent respectivement n et m impulsions, alors la 20 complexité de ce calcul sera d'ordre O(nxm). Cette complexité peut être évitée. En effet : P P puis j=1 et finalement = , = ryr =Rq.Rq' 1E C7,Z~ECq J=1 25 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. Si on considère le critère de Condorcet utilisé précédemment, ce critère s'écrit de la façon suivante : n n C (X) E E (c',x', i =1 ;'=1 Ce calcul est apparemment quadratique. Il peut être ramené en 5 linéaire grâce au principe pré-établi : n n n n C (X) X ±EIC 1=1 '=1 1=1 '=1 E Cil, = hi.hq et C11. ieC EC il suit : C(X)---- hl n(1?> 1 . ) 1 1 ,q *I 1 ,e1 On obtient ainsi une expression calculable en 0(n). 10 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 15 linéaires. PRM principal PRM secondaire fréquence groupe f groupe PRM phase imp 1 f1 1 1 p1 imp i fi 1 1 pi imp j fj 2 4 pi imp n fn 2 4 pn Annexe 1 Groupe RF Effectif moyenne ect Liste des impulsions du groupe G1 1 2 18 21 ............... ...... ... 392 G2 186 201 209 312 ......... 415 NA NA Gn 39 I 53 I ........... 178 NA ... ...... ... I NA Annexe 2 Groupe RF Effectif moyenne ect 11 12 13 14 15 16 17 18 ................... In G1 1 0 1 0 0 0 1 0 0 G2 0 0 0 1 0 0 0 1 1 Gn 0 1 0 0 1 0 0 0 0 Annexe 3 Classe 1 Classe 2 Classe 3 Impulsions enrichies Modes DOA Modes Freq Modes Li { Modes PRM moc....see.ILM. 11 12 13 14 15 16 17 ... In DOM 0 0 1 1 0 1 0 0 DOA2 1 1 00 1 00 0 DOA 3 0 0 0 0 0 0 1 1 DOAn ... ... ... ... ... ... ... ... Fréq.1 1 1 0 0 1 00 0 Fréq.2 0 0 0 1 0 1 0 0 Fréq.3 0 0 1 0000 0 Fn... ... ... ... ... ... ... LII 0 0 0 0 0 0 1 1 L12 1 1 00 1 00 0 Li3 00 0 1 0 1 Cl 0 Lin... ... ... ... ... ... ... ... PRMI 0 0 1 1 0 1 0 0 PRM2 0 0 0 0 0 0 1 1 PRM3 1 1 0 0 1 00 0 PRMn ... ... ... ... ... --- ... 7jeceb...L......,TIP 13 '14 16 11 12 15 In DOM 1 1 1 0 0 0 0 Fréq.2 0 1 1 0 0 0 0 Fréq.3 1 0000000 Li3 0 1 1 0 0 0 0 PRM1 1 1 1 0 0 0 0 DOA2 0001 1 1 00 Fréq.1 0001 1 1 00 Lit 0 O 01 1 1 00 PRM3 0001 1 1 00 DOA3 000000 1 1 Fréq.4 0 0 0 0 0 0 1 1 Lit 0000001 1 PRM2 0000001 1 0000000 Tableau impulsions-modes Tableau impulsions-modes réorganisé Annexe 4 11 12 13 14 15 16 17 ... IN A0A1 0 0 nn 0 1 0 0 A0A2 .: :: (1)p ô 1 0 0 0 AOAm 0 0 0 0 0 0 1 1 F1 el ) CL ) 0 0 1 0 0 0 F2 0 0 0 1 0 1 0 0 Fn 0 0 1 0 0 0 0 0 IP1 0 0 0 0 0 0 1 1 IP2 .)0)0 0 1 0 0 0 Ipo 0 0 0 1 0 1 0 0 PRM 0 0CD(00 1 0 0 PRM 0 0 0 0 0 0 1 1 PRM ei)e00 0 1 0 0 0 Annexe 5a 11 12 13 14 15 16 17 ... IN AOA 1 - e - .: -t : :1 : 0 0 AOA 2 exp 0 0II 0 0 0 AOAm 0 0 0 0 0 0 1 1 F1 010 0 0 1 0 0 0 F2 0 0 0 1 0 1 0 0 Fn 0 0 1 0 0 0 0 0 IP1 0 0 0 0 0 0 1 1 IP2 1 1 0 0 1 O. 0 0 Ipo 0 0 0 1 0 1 0 0 PRM1 0 0 CID(1) 0 (0 0 0 PRM2 0 0 0 0 0 Q 1 1 PRMp _ 1 1 0_ 0 1 0 0 0 C3,4 2 Ci ----- xi - xi Cit =ICI( cA0A2,F1 =__ 3 Annexe 5b 11 12 13 14 15 16 17 18 19 110 111 112 113 114 115 116 117 118 119 120 F1 1 1 00000000 0 0 0 0 0 0 0 0 0 0 F2 00 1 1 1 1 1 000 0 0 0 0 0 0 0 0 0 0 F3 0000000 1 1 0 0 0 0 0 0 0 1 1 1 1 F4 000000000 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 0 DOA200000 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 DOA30000000000 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 0 PRM200000 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 PRM30000000000 0 1 1 1 1 1 1 1 1 1 Annexe 6a ) 11 12 13 14 15 16 CI3 17 18 19 1:0 11: 12 0 0 0 0 0 0 0 0 0 0 0 0 Critère 1076 mode 1 mode 2 mode 3 mode 4 mode 5 mode s mode 7 mode. F. mode 9 mode 10 0 0 1 1 0 U 0 0 0 U 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 113 114 115 116 117 118 119 120 0 0 0 0 0 0 0 b 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 U U U U U U U 0 3 D G 0 11 0 0 1 1 1 1 ï 1 0 0 T 2 1 I U D0A1 PRM1 F1 F2 DOA2 PRM2 Cl2 Annexe 6b 0 1 0 0 0 0 1-3 F4 DOA3 PRM3 CI13 Annexe 6c C112 CI15 11 r 12 13 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 17 18 19 110 111 112 113 0 0 0 0 0 0 0 0 0 0 0 0 0 0 II 0 U 0 0 0 114 115 116 117 118 119 120 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 U (} (1 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 F3 F4 DOA3 PRM3 00A1 PRIV11 1-1 F2 00A7 PRM2 Critère 912 mode 1 mode 2 mode. mode 4 mode 5 mode 6 mode 7 mode 8 mode 9 mode 10 11 12 13 14 15 16 17 18 19 110 111 112 113 114 115 116 117 118 119 120 Mode 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 DOA1 Mode 2 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 PRM1 Mode 3 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 F1 Mode 4 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 F2 Mode 5 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 DOA2 Mode 6 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 PRM2 Mode7 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 1 1 1 F3 Mode8 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0 F4 Mode9 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 DOA3 Mode10 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 PRM3 Critere 203 mode 1 mode 2 mode 3 mode 4 mode 5 mode 6 mode 7 mode 8 mode 9 mode 10 Annexe 7 11 12 13 14 15 16 17 18 19 110 111 1 1 1 1 1 Ci Ci 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0 U U 0 0 U 0 U U 0 U 0 0 0 0 0 0 0 0 0 0 0 112 113 114 115 116 117 118 119 120 0 0 CI 0 0 0 0 Cl 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 DOA1 PRM1 F1 F2 DOA2 PRM2 F3 :4 JOA3 r>Rm3 Annexe 8 w 0 1-1 1-1 w w 00 NJ n (b w Il C Il C O 3 3 3 3 3 3 3 3 3 3 o o o o o o o o o o Q il CL il 11 il . CL lb lb fb II 0 0 0 lb lb fb O tO CO Cl Ln A tO I-' 3 3 3 3 3 3 3 3 3 3 0 0 0 0 0 0 0 0 0 0 CL 0_ 11 Il 0_ 0_ 0_ 0_ lb lb lb lb fb fb fb M 0 0 P N Ob V 01 Ul A UT N P 0 3 3 3 3 3 3 3 3 O o o 0 o o o o o man_b 0 0 0 ro ro 0 3 3 3 3 3 3 3 3 3 3 o o o o o o o o o 0 CL CL CL CL CL 0_ CL CL 0_ 0_ 0 lb 0 0 0 0 ro ro /b - LO CO v 01 Ln A (.0 1,) UD CO CID Ln 4a W hi o O Q Q Q il Ne O000000r->- O0000ao O000000.i; 0 0 0 0 CD CD 0 Fr; O000loo^-"" o o o co O00000ieir' N o O 0 0 0 0 0 O o o ce; o 000 o 000 ; o 000 O000 ; O000tn; O0000; O000 o o 11 Q T T -0 Q 11 'Tl DO 0 A LU D) 0 1,1 W W N 1. r r r ti t) 71 0 et W > W O 0 0 0 0 0 0 0 0 0 0 i5 -g 0 m T lj 0 IV N 0 CP Ci 0 Ci 1)1 el Fe Fe MMM000 bb Me Ci di 0 Fe 1.1. if 0 0 0 1.1 là M O O O blif'14000 !-- 1,- CD X CO 1-, - C> 0 CD CD 0 0 in. 0 0 - 0 1 0 0 0 I:: I I o o o o 0 OIS, 0 e- 010 0 ^-" fi-n 0Q c) 1,4 ID Q f> )> ba w (:) 1-1 1-1 w w 00 Q TIr' 'V m MI 0 A lu xi 0 NT -ce '11 W -13 70 o 0 0i-, 0 s. O 0 I-% D 1-1 Critère Cr = F3.=F3-C1 Critère 251 Cr = Cl-F3 P;1 = Critère 267 I C- = F3 1 L 1 .1 18 19 mode 1 mode 2 O mode 3 O 0 mode 4 0 mode 5 mode 6 mode 7 mode 8 mode 9 mode 10 O 0 18 19 mode 1 mode 2 mode 3 mode 4 mode 5 mode 6 mode 7 mode 8 u 'J mode 9 mode 10 O 112 113 114 115 116 117 118 119 120 o o o 0 o o o o o o o o o o o o o o o o o 0 0 o o o o O 0 00 0 0 o 0 o o o o o o o o o o o o o o o o o o o O0 O 0 01111 1111100 00 111111111 11111 1 1 1 1 112 113 114 115 116 117 118 119 120 o o o o o o o o a o o a o o o o O o o o o o o o o O o o o o o o o o o o o o o o o o o o o o o o o o O o o 0 o 1 1 1 1 1 0 C) 1 1 1 1 1 1 1 1 1 1 DOA1 PRM1 F1 F2 DOA2 PRM2 F3 F4 DOA3 PRM3 DOA1 PRM1 F1 F2 DOA2 PRM2 F3 F4 DOA3 PRM3 18 19 112 113 114 115 116 117 118 119 120 mode 1 O 0 O O O 0 O O 0 0 O DOA1 mode 2 o O O 0 O O O O 0 O O PRM1 mode 3 o 0 O 0 0 0 O O O O O F1 mode 4 o O O 0 0 O 0 0 O O 0 F2 mode 5 1 0 0 0 0 o 0 0 0 0 DOA2 mode 6 1 O 9 9 0 0 O PRM2 mode 7 0 0 0 0 T 1 1 1 1 IFS mode 8 0 111114 0 O 0 0 I F4 mode 9 0 0 1 1 1 1 1 1 1 1 1 1 DOA3 i mode 10 0 0 r 1. 1. 1. 1. 1 1 j PRM3 Critère 231 C = C11_)F3 18 19 mode 1 O 0 mode 2 mode 3 mode 4 O 0 mode 5 mode 6 mode 7 mode 8 mode 9 mode 10 120 0 o 112 113 114 115 116 117 118 119 o o o o o o o 0 o o o o 0 o o o o o o o o o o o o o o o o o o o Annexe 10a DOA1 PRM1 F1 F2 DOA2 PRM2 F3 F4 DOA3 PRM3 110000000-0 11000 0000-0 11000001111 o o 1 1 1 1 1 0 0 9 0 O 0 1 1 1 1 1 1 1 / 1 O 0111122111 Critère mode 1 mode 2 mode 3 mode 4 mode 5 mode 6 mode 7 mode 8 18 19 110 111 112 113 114 115 116 117 118 119 120 0 0 0 0 0 0 0 0 0 0 0 0 0 DOA1 0 0 0 0 0 0 0 0 0 0 0 0 0 PRM1 0 0 0 0 0 0 0 0 0 0 0 0 0 F1 O 1 0 0 0 0 0 0 0 0 0 F2 1 1 1 1100000000 0 DOA2 0 0 0 0 0 0 0 0 0 PRM2 F3 F4 0 mode 9 0 0 0 0 1 1 1 1 1 DOA3 mode 10 0 0 0 o C3={4444 I c2 110,111 IV Annexe 10b Critère 18 19 112 113 114 115 116 117 118 119 120 307 mode 1 0 000 0 0 0 0 0 0 0 DOA1 mode 2 0 000 0 0 0 0 0 0 0 PRM1 mode 3 0 00 0 0 0 0 0 0 0 CI F1 mode 4 0 Q 0 0 0 0 0 0 0 0 O F2 1 mode 5 1 1 00 0 0 0 0 0 0 0 DOA21 mode 6 0 0 0 0 0 0 0 0 0 PRM2 mode 7 F3 mode 8 0 0 Fit_ mode 9 0 0 DOA3 mo-e 1 3 Annexe 11 mode 1 mode 2 mode 3 mode 4 mode 5 mode 6 mode 7 mode 8 mode 9 mode 10 mode 1 mode 2 mode 3 mode 4 mode 5 mode 6 mode 7 mode 8 mode 9 mode 10 Critère IE 17 18 19 1:0 111 68 mode 1 0 0 0 0 0 0 DOA1 mode 2 O 0 0 0 0 0 PRM1 mode 3 0 0 F1 mode 4 1 0 c 0 0 F2 mode 5 1 1 1 1 1 DOA2 mode 6 PRM2 mode 7 0 4 1 1 0 0 F3 mode 8 0 0 0 0 1 1 F4 mode 9 mode 10 0 0 0 0 O 0 DOA3 0 0 0 0 O 0 PRM3 16 17 18 19 110 111 DOA1 O 0 0 o PRM1 O 0 0 0 o o F1 O 0 0 0 F2 DOA2 PRM2 »e3 F4 DOA3 PRM3 O 0 0 0 0 0 O00000 16 17 18 19 110 111 68 moce 1 0 0 0 0 0 OMM moce 2 ODO 0 0 0 PRM1 moce 3 0 0 0 0 F1 moce 4 1 1 0 0 F2 moce 5 1 1 1 1 1 1 DOA2 moce 6 PRM2 moce 7 t1 0 1 1 0 0 F3 moce 8 0 D 0 0 1 1 F4 moce 9 0 D 0 0 0 0 DOA3 moce 10 0 0 0 0 0 0 PRM3 16 17 18 19 110 111 Cr tère O00000D0A1 O00000PRM1 F1 DOA2 O 00A3 - PRM3 PRM2 F3 F4 } 0 0 0 0 0 Annexe 12a 16 17 18 19 110 111 112 113 114 I.5 116 117 118 119 20 model 0 0 0 0 0 0 0 0 0 0 0 0 G 0 0 DOA1 mode 2 0 0 0 0 0 0 0 0 0 0 C 0 C 0 0 PRM1 mode 3 On 0 o 1-1 n 0 0 0 0 C 0 C D 0 F1 mode 4 1 0 0 0 0 0 0 0 0 0 0 0 0 0 F2 mode 5 1 1 0 0 0 0 0 0 0 D 0 DOA2 0 PRM2 F3 F4 DOA3 PRM3 0 o 0 O 0 0 0 0 0 O 0 0 0 0 0 mode 6 mode 7 mode 8 mode 9 mode 10 Annexe 12b n .3 3 3 3 3 3 3 3 3 3 0 0 0 o 0 o 0 o 0 o CL 0 I-. 0 O. CL a a_ D- a aga 0 ID CD CD il Dl 0 CO éD US CO --/ Ch ln A Lu 1,, 1-, 0 0 0 0 0 0 0 0 0 0 CD X CD 0 0 0 0 0 o 0 0 0 0 0 0 0 0 0 0 0 O o I 0 0 zoo0c7, 0 0 z000re 0 0 3--- i, 0 0 U -n T Q T 71 "Ci U » 0 .1o. W 70 0 Ki > > > N Ni .^^--^^^ n ID 3 3 3 3 3 3 3 3 3 3 O o o o o 0 o o o o . CL CL . - 0 CD M CD 0 CD 0 0 P IO CO V CIO Ln A Ln N P 0 3 3 3 3 3 3 3 3 3 3 O o o 0 o o 0 o o D_ CL D D. D_ D CI D_ D_ C_ aamaammmmm 10 ta CO ---1 01 In 0 0 0 0 0 O 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 n O 0 9 0 0 0 r 0 r r N 0 0 0 0 0 0 r O r r Cr.)C.A) NJ n 0- 3 3 3 3 3 3 3 3 3 3 o o o o o 0 o 0 o a a a 0. cL a a Q. 0. - D D di Co D à D e-' l0 Co V G1 lA la Lu Al 0 13 ID m,-,:10-n-n-p xi 0 e. uà motvi-s,30 toLe 1,1u w O 1-1 1-1 w w 00 0 0 7' -o 4". W -n 7' rsi O N I"-, O O o -p O -n -n O O -II il 17 7" 0 'n A0Nr 710 17' 0 0 00 XI 0 > O 0 0 0 P 0 0 0 0 P 0 0 213 0 0 0 0 03 0 0 1p 0 0 O 0 O 0 0 Critère mode 1 mode 2 mode 3 1 1 mode 4 mode 5 0 mode 6 0 0 mode 7 0 0 mode 8 0 0 mode 9 mode 10 13 14 15 DOA1 PRM1 0 0 0 F1 F2 0000D0A2 0 0 0 PRM2 0 0 0 F3 0 0 0 F4 00000D0A3 00000PRM3 13 14 15 34 Critère 11 12 11 12 51 mode 1 DOA1 mode 2 PRM1 mode 3 1 1 0 0 O F1 mode 4 0 0 1 1 1 F2 mode 5 00000D0A2 mode 6 00000PRM2 mode 7 00000F3 mode 8 00000F4 mode 9 00000D0A3 mode 10 00000PRM3 Annexe 14 Critère 1076 mode 1 mode 2 mode 3 mode 4 mode 5 0 0 0 0 mode 6 0 0 0 0 mode 7 0 0 0 0 mode 8 0000 mode 9 00000 mode 10 0 0 0 0 17 18 19 0000000 00000000 112 110 111 0 0 0 0000 113 000 000 000 114 115 116 117 000000000A' 0000 0000 118 119 120 0 0 0 0 0 PRM1 0 0 0 0 OF1 0 F2 OD0A2 PRM2 F3 F4 DOA3 PRM3 Il 12 13 14 15 16 Annexe 15 11 12 13 14 15 16 17 ... IN AOA 1 0 0 1 1 0 1 0 0 AOA 2 1 1 0 0 1 0 0 0 ..... AOAm 0 0 0 0 0 0 1 1 F1 1 1 0 0 1 0 0 0 F2 0 0 0 1 0 1 0 0 Fn 0 0 1 0 0 0 0 0 IP1 0 0 0 0 0 0 1 1 IP2 1 1 0 0 1 0 0 0 ..... Ipo 0 0 0 1 0 1 0 0 PRM1 0 0 1 1 0 1 0 0 PRM2 0 0 0 0 0 0 1 1 ..... PRMp 1 1 0 0 1 0 0 0 Annexe 16 11 12 13 14 15 16 17 ... IN 11 1 1 0 0 1 0 0 0 12 1 0 0 1 0 0 0 13 1 0,4 0 0,4 0 0 14 1 0 1 0 0 15 1 0 0 0 16 1 0 17 1 1 ... IN 1 Annexe 17 13 14 16 11 12 15 17 IN ... AOA 1 1 1 1 0 0 0 0 0 AOA 2 0 0 0 1 1 1 0 0 ..... AOAm 0 0 0 0 0 0 1 1 F1 0 0 0 1 1 1 0 0 F2 0 1 1 0 0 0 0 0 Fn 1 0 0 0 0 0 0 0 IP1 0 0 0 0 0 0 1 1 IP2 0 0 0 1 1 1 0 0 ..... Ipo 0 1 1 0 0 0 0 0 PRM1 1 1 1 0 0 0 0 0 PRM2 0 0 0 0 0 0 1 1 ..... PRMp 0 0 0 1 1 1 0 0 Annexe 18 AOA1 AOA2 AOAm F1 F2 FN IP1 IP2 I PO PRM1 PRM2 PRMP AOA 1 1 0 0 0 0,66 0,33 0 0 0,66 1 0 0 AOA2 1 0 1 0 0 0 1 0 0 0 1 AOAm 1 0 0 0 1 0 0 0 1 0 F1 1 0 0 0 1 0 0 0 1 F2 1 0 0 1 0 0,66 0 0 Fn 1 0 0 0 0,33 0 0 I PI 1 0 0 0 1 0 I P2 1 0 0 0 1 Ipo 1 1 0 0 PRM1 1 0 0 PRM2 1 0 PRMp 1 Annexe 19 11 13 14(16 111 12 15 17 IN AOA 1 1 1 1 0 0 0 0 0 PRM1 1 1 1 0 0 0 0 0 F2 0 1 1 0 0 0 0 0 Ipo 0 1 1 0 0 0 0 0 Fn 1 0 0 0 0 0 0 0 -A0A2 0 0 0 1 1 1 0 --0-----'6 0 F1 0 0 0 1 1 1 1P2 0 0 0 1 1 1 0 0 PRMp 0 0 0 1 1 1 0 0 AOAm 0 0 0 0 0 0 1 f 1P1 0 0 0 0 0 0 1 1 PRM2 0 0 0 0 0 0 1 .-'---- 1 CI3 Annexe 20 Cl2 REVENDICATIONS1. Procédé de désentrelacement non supervisé de trains 5 d'impulsionscaractérisé en ce qu'il comprend : - une étape Etp10 de constitution de N, N étant un entier supérieur à 1, histogrammes à une ou plusieurs dimensions représentatifs des trains d'impulsions; 10 - une étape Etp20 d'extraction de modes à partir desdits histogrammes ; - une étape Etp30 de constitution de groupes d'intérêt à partir desdits modes extraits ; - une étape Etp40 d'enrichissement de la description des impulsions 15 par son appartenance à au moins un groupe d'intérêt ; - une étape Etp50 de classement croisé non supervisé des modes afin de réorganiser les impulsions enrichies. 2. Procédé selon la revendication précédente caractérisé en ce que l'étape 20 Etp20 comprend : une étape Etp21 de seuillage des histogrammes ; une étape Etp22 de décomposition des histogrammes en modes ; une étape Etp23 de caractérisation des modes ; 25 3. Procédé selon une des revendications 0 à 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 : 30 - une étape Etp51 de recherche, parmi les modes restants, du mode ayant un cardinal d'intersection avec la classe en cours le plus grand ;- 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 ; - une étape Etp53 de fusion du mode avec la ou les classes en cours. 5 4. Procédé selon une des revendications 0 à 2 caractérisé en ce que l'étape Etp50 de classement comprend : - une étape Etp56 de classification des impulsions en fonction de leurs liens de corrélation mutuelles ; 10 - une étape Etp57 de classification des modes en fonction de leurs liens de corrélation mutuelles ; - une étape Etp58 de regroupement des modes et impulsions en classes. 15