Brevet FR3051611A1 — publié le 20171124

Procédé de traitement d'un signal formé d'une séquence d'impulsions La présente invention concerne un procédé de traitement d'un signal formé d'une séquence d'impulsions, comportant au moins un motif répétitif formé d'au moins une impulsion, ledit motif étant répété dans le signal avec une période de répétition de motif. L'invention se situe dans le domaine du traitement des signaux, en particulier des signaux numériques, et trouve une application particulière dans le domaine de la caractérisation de séquences d'impulsions radar. En effet, un signal périodique ou cyclo-périodique, formé de motifs d'impulsions qui se répètent, est en général caractérisé par des paramètres dont la période de répétition de motif, ou PRM, et les intervalles de répétition des impulsions, ou PRI du signal. En pratique, la caractérisation d'un signal radar intercepté dans un créneau temporel d'écoute permet notamment d'identifier le mode d'émission du signal intercepté ou de pister des séquences d'impulsions provenant d'une même émission. De manière connue, les intervalles de répétition des impulsions sont utilisés pour caractériser une telle séquence d'impulsions. Cependant, les valeurs de PRI sont faussées en cas de perte aléatoire d'impulsions du signal, phénomène également connu sous le nom de mitage, qui se produit en pratique. L'invention vise à améliorer la caractérisation des périodes de répétition temporelle des impulsions dans une séquence d'impulsions, notamment dans le cas de pertes aléatoires. A cet effet, selon un premier aspect, l'invention propose, un procédé de traitement d'un signal formé d'une séquence d'impulsions, comportant au moins un motif répétitif formé d'au moins une impulsion, ledit motif étant répété dans le signal avec une période de répétition de motif. Le procédé comporte les étapes suivantes : - estimation de la période de répétition de motif dudit signal, - calcul d'une séquence de phases en fonction d'une date d'arrivée de chaque impulsion par rapport à une date d'arrivée de référence choisie et de la période de répétition de motif estimée, - à partir de ladite séquence de phases calculées, estimation d'au moins une valeur de phase et d'un écart-type associé, ladite valeur de phase étant associée à un moment de phase représentatif du motif répétitif ; - obtention et exploitation de paramètres de caractérisation du signal en utilisant les valeurs de phase estimées. 3051611 2 Avantageusement, le procédé de l'invention permet de caractériser les répétitions temporelles caractéristiques d'une séquence d'impulsions de manière robuste aux pertes aléatoires d'impulsions. Le procédé proposé s'applique pour les séquences d'impulsions à trains 5 d'impulsions imbriqués ou successifs, la période de répétition des impulsions d'un train d'impulsions pouvant être ou non affectée d'un bruit aléatoire ou pseudo-aléatoire induisant une fluctuation de la période de répétition par rapport à une valeur moyenne. Le procédé selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prises indépendamment ou dans toutes leurs combinaisons techniquement 10 acceptables. Il comprend, après l'étape de calcul d'une séquence de phases, une étape de construction d'un histogramme de phases, et d'extraction d'un ou plusieurs pic(s) de l'histogramme, chaque pic étant associé à un moment de phase de la séquence d'impulsions. 15 L'étape d'estimation d'une phase associée à un moment de la séquence d'impulsions, ledit moment étant associé à un pic de l'histogramme de phases, comprend une extraction d'un sous-ensemble de dates d'arrivée autour dudit pic, et une estimation précise de la valeur de phase associée audit moment, et d'une variance associée à la valeur de phase estimée. 20 Lorsque la séquence d'impulsions comprend une pluralité de trains d'impulsions imbriqués, chaque train d'impulsions ayant un moment de phase associé, l'estimation précise de la phase associée au moment comprend une sommation vectorielle de vecteurs unités, chaque vecteur ayant une orientation égale à une valeur de phase calculée, avec la période de répétition de motif dudit signal préalablement estimée, à 25 partir d'une date d'arrivée dudit sous-ensemble, ladite sommation vectorielle fournissant un vecteur résultant. La valeur de phase associée audit moment estimée est égale à l'argument dudit vecteur résultant. Lorsque la séquence d'impulsions comprend une pluralité de trains d'impulsions 30 successifs, chaque train d'impulsions ayant un moment de phase associé, l'estimation précise de la phase associée à un moment traité comprend en outre une estimation d'une période de répétition de motifs correspondant audit moment à partir du sous-ensemble de dates d'arrivée considérées. L'estimation de la période de répétition de motif du signal comporte une première 35 estimation grossière d'une première période de répétition de motif, et une deuxième 3051611 3 estimation raffinée d'une deuxième période d'estimation d'une première période de répétition de motif. Lors de l'étape d'obtention d'une séquence de phases, pour chaque date d'arrivée d'impulsion tn, on applique un repliement de période égale à la période de répétition de 5 motif estimée selon la formule suivante : ( tn - to et'),, = tn. j Où Lx est la partie entière inférieure du nombre réel x, tO est la date d'arrivée de référence et X la période de répétition de motif estimée. Selon un autre aspect, l'invention concerne un dispositif de traitement d'un signal 10 formé d'une séquence d'impulsions successives, apte à mettre en oeuvre un procédé tel que brièvement décrit ci-dessus. Ce dispositif comporte au moins un processeur de calcul apte à mettre en oeuvre : - une unité d'estimation de la période de répétition de motif dudit signal, - une unité de calcul, en fonction d'une date d'arrivée de chaque impulsion par 15 rapport à une date d'arrivée de référence choisie, et de la période de répétition de motif estimée, d'une séquence de phases, -une unité d'estimation, à partir de ladite séquence de phases calculées, d'au moins une valeur de phase et d'un écart-type associé, ladite valeur de phase étant associée à un moment de phase représentatif du motif répétitif ; 20 - une unité d'obtention et d'exploitation de paramètres de caractérisation du signal numérique en utilisant les valeurs de phase estimées. Selon un mode de réalisation le dispositif comprend en outre une unité de construction d'un histogramme de phases, et d'extraction d'un ou plusieurs pic(s) de l'histogramme, chaque pic étant associé à un moment de phase de la séquence 25 d'impulsions. Selon un autre mode de réalisation, l'invention concerne un programme d'ordinateur comportant des instructions logicielles qui, lorsqu'elles sont mises en oeuvre par un dispositif programmable, mettent en oeuvre un procédé de traitement d'un signal formé d'une séquence d'impulsions tel que brièvement décrit ci-dessus. 30 Selon un autre aspect, l'invention concerne un support d'enregistrement d'informations, comportant des instructions logicielles pour l'exécution d'un procédé de traitement d'un signal formé d'une séquence d'impulsions tel que brièvement décrit ci-dessus, lorsque ces instructions sont exécutées par un dispositif programmable. 3051611 4 D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles : 5 - la figure 1 est une vue schématique d'une séquence d'impulsions à trains d'impulsions imbriqués ; - la figure 2 est un schéma synoptique des principales étapes d'un procédé selon un mode de réalisation de l'invention ; - la figure 3 illustre schématiquement plusieurs périodes de répétition d'un signal 10 numérique et un histogramme de phases associé ; - la figure 4 est un schéma représentant les blocs fonctionnels d'un dispositif programmable apte à mettre en oeuvre l'invention. La figure 1 illustre schématiquement un signal numérique S comportant une 15 séquence d'impulsions 10, composée d'impulsions I,, 12, ...ln reçues successivement dans le temps. Chaque impulsion lk est illustrée schématiquement sur un axe temporel. La séquence d'impulsions 10 est une séquence à répétition de motif, et comprend, dans l'exemple simplifié illustré, quatre trains binaires périodiques imbriqués. Bien entendu, dans la pratique, des signaux numériques réels traités comprennent 20 un nombre P quelconque de signaux imbriqués. Dans l'exemple de la figure 1, le motif 12 est composé de quatre impulsions successives, à 14, et de période de répétition de motif X. Les différences entre dates d'arrivée sur l'axe temporel, connues sous l'abréviation DTOA (« difference of time of arrivai ») entre deux impulsions successives permettent de 25 calculer les valeurs d'intervalle entre impulsions successives, notés respectivement PRI1 à PRI4. Chacun des trains d'impulsions imbriqués dans la séquence d'impulsions 10 est constitué d'impulsions se répétant avec une période de répétition X, égale à la période de répétition de motif de la séquence d'impulsions 10. 30 En d'autres termes, chaque train d'impulsions Tp est constitué des impulsions ~Ik2+pJ' La date d'arrivée de référence to est choisie arbitrairement, par exemple égale à la date d'arrivée de la première impulsion de la séquence d'impulsions considérée. La perte aléatoire d'une sous-partie des impulsions induit une erreur sur le calcul 35 des valeurs de PRIk d'une séquence d'impulsions reçue. 3051611 5 Il est proposé de calculer des valeurs de phase pour chaque impulsion, représentatives de la différence P,, pour chaque impulsion, entre sa date d'arrivée et une date de référence périodisée par rapport à la période de répétition de motif. En particulier, P1 à P4 sont illustrées sur la figure 1 pour l'exemple considéré. 5 Les valeurs de phase sont calculées comme suit. A partir d'une séquence de dates d'arrivée {tk}, on calcule des phases, avec un repliement de période X donnée, par rapport à une date d'arrivée de référence to par la formule suivante : 0, 9)(tic =2 2 tk (EQ 1) 10 Où Lx] est la partie entière inférieure (ou par défaut) du nombre réel x. Les valeurs de phase ek sont supérieures ou égales à 0 et strictement inférieures à X. Les valeurs de phase ainsi calculées ont une meilleure robustesse aux pertes aléatoires d'impulsions que les valeurs DTOA mesurées entre impulsions successives, 15 lorsque les valeurs respectives de date d'arrivée de référence to et de période X sont connues ou préalablement déterminées. Les principales étapes d'un mode de réalisation du procédé de traitement de signal selon l'invention seront décrites ci-après en référence à la figure 2. Ce procédé s'applique dans le cas d'une séquence d'impulsions formée d'une 20 pluralité de P trains d'impulsions imbriqués, avec P un entier supérieur à 1, chaque train d'impulsions formant un moment de phase m(i) de la séquence d'impulsions. Une telle séquence d'impulsions est également appelée « stagger à P moments ». Par exemple, la figure 1 illustre un signal stagger à 4 moments. Le procédé de l'invention s'applique également au cas d'une séquence 25 d'impulsions formée d'une pluralité de P trains d'impulsions qui se succèdent temporellement, chacun des trains d'impulsions ayant une période de répétition de motif associée, égale à la période de répétition des impulsions dudit train d'impulsions. Il n'est pas nécessaire de connaître le type de séquence à l'avance, le procédé permet de le déduire. 30 Dans un cas d'application, la séquence d'impulsions est en plus affectée, volontairement à sa création par le concepteur du radar, par un bruit appelé jitter en anglais, la période de répétition des impulsions fluctuant de manière aléatoire ou pseudoaléatoire autour d'une valeur moyenne. 3051611 6 Chaque moment m(i) a une valeur de phase associée, les impulsions du train d'impulsion T, associé à m(i) ayant une période de répétition 2' qui est la même que la période de répétition associée à la séquence d'impulsions. Dans le mode de réalisation décrit, on dispose des dates d'arrivée : {tk k 5 n} 5 des impulsions d'une séquence d'impulsions à traiter, dans un référentiel temporel donné. Ces dates d'arrivée sont par exemple mémorisées dans une mémoire d'un dispositif de calcul apte à mettre en oeuvre l'invention. Lors d'une première étape 20 d'estimation de la période de répétition de motif, mise dans oeuvre dans le mode de réalisation décrit, une valeur de PRM de la séquence 10 d'impulsions, notée , est calculée. Ce calcul comprend deux sous-étapes, une première sous-étape 22 de première estimation, également appelée estimation grossière, permettant d'obtenir une première valeur estimée 2g, de la PRM, suivie d'une deuxième sous-étape 24 de deuxième estimation, appelée également estimation fine, permettant d'obtenir une deuxième valeur 15 estimée, , de la PRM. La première sous-étape 22 met en oeuvre toute méthode connue d'estimation de la période de répétition de motif. Par exemple, des méthodes d'estimation de motif sont décrites dans le document « New techniques for the deinterleaving of repetitive sequences », de H. K. Mardia et al, publié dans IEE Proceedings F, vol 136 N°4, pp. 149- 20 154, en août 1989.Cette estimation 22 est réalisée par exemple par les méthodes d'histogrammes des différences de dates d'arrivée (ou DTOA). En alternative, toute autre méthode d'estimation de la période de répétition de motif, utilisant une transformation de Fourier ou une transformation en ondelettes peut être utilisé. Cependant, comme dans cette première sous-étape 22 la précision de l'estimation n'est pas recherchée, il est 25 préférable d'utiliser une méthode simple et peu coûteuse en calcul, plutôt qu'une méthode basée sur des transformations qui ont un coût calculatoire important. Lors de la sous-étape 24 de deuxième estimation de la PRM, on calcule, pour chaque date d'arrivée mémorisée, un indice entier associé en fonction de la première valeur 2g de la PRM : 30 k(i) =[ , i =1,....,n (EQ 2) g Ensuite on effectue le calcul suivant : k(i)- k (n) -1) t (EQ 3) k (n) - (k (n) 2 - 1) j 3051611 7 La variance associée à cette estimation est donnée par : n3 -n 12 62 Var(A) (EQ 4) n étant le nombre d'impulsions de la séquence d'impulsions, et cr étant l'écart- 5 type de la mesure de temps t. L'écart-type o- est celui de l'instrument de mesure des dates d'arrivée. Il est connu pour chaque instrument, pour les conditions de mesures réalisées en interception. Les calculs indiqués ci-dessus s'appliquent notamment pour une modélisation des dates d'arrivée selon un modèle linéaire en fonction des valeurs de k(i), autrement dit 10 pour une suite d'impulsions périodique de période : = to +2. k(i)+ c (EQ 5) Les calculs ci-dessus s'appliquent pour des séquences d'impulsions à P trains imbriqués, éventuellement bruitées par un bruit non cumulatif. Dans le cas d'un bruitage par bruit non-cumulatif, chaque date d'arrivée peut 15 s'écrire de la manière suivante : = (i -1)2 + 0m(i) + (EQ 6) Où y, est un bruit centré d'écart-type Ifs, et 0,nwla valeur phase du moment m(i). Selon une variante, la séquence d'impulsions reçue est affectée par un bruit (jitter) cumulatif : 20 t = 2+ tfri + + (EQ 7) Où y, est un bruit blanc gaussien centré d'écart-type cs, et 0 la phase du train d'impulsions considéré de la séquence. Dans cette variante, on applique le calcul suivant pour la deuxième estimation de la PRM : 25 = tn-t, (EQ 8) n -1 La variance associée à cette estimation est donnée par : Var(2)= cr2 (EQ 9) n-1 n étant le nombre d'impulsions de la séquence et o" étant l'écart-type de la mesure de temps t. Cet écart-type cumule quadratiquement l'écart type de mesure du 30 temps t et l'écart-type ai, de la loi de jitter. 3051611 8 De retour à la figure 2, l'étape 20 d'estimation de la PRM est suivie d'une étape 30 de calcul d'une séquence de phases, par calcul selon la formule (EQ1), prenant en compte la date d'arrivée de référence to choisie et la valeur de la PRM estimée j, obtenue à l'étape 20. 5 On obtient donc une séquence de phases : {Bk ,1 k n}, Bk E [0, > [ L'étape d'obtention d'une séquence de phases 30 est suivie d'une étape 40 de construction d'histogramme de phases. La figure 3 illustre de manière schématique une séquence d'impulsions mitée, comprenant 4 périodes successives, notées Per-1 à Per-4, superposées dans la 10 représentation illustrée de manière à illustrer le repliement, et l'histogramme de phases H associé. Il apparaît de cette illustration schématique que le motif Mg de la séquence d'impulsions, qui est un motif dit stagger 8 dans cet exemple, est reconstituable à partir de l'histogramme de phases. 15 Bien entendu, il s'agit d'un exemple schématique. En pratique, les séquences d'impulsions traitées comprennent des impulsions plus rapprochées temporellement et/ou un bien plus grand nombre de trains d'impulsions imbriqués. De retour à la figure 2, l'étape 40 de construction d'histogramme est suivie d'une étape 50 de recherche de pics dans l'histogramme de phases. 20 On applique toute méthode connue à cet effet. Dans un mode de réalisation, les pics retenus dans l'histogramme sont les raies de l'histogramme dont la hauteur dépasse strictement un seuil choisi. Lorsqu'il n'existe qu'un très faible nombre de raies dans l'histogramme à proportion de la taille de l'histogramme, par exemple moins de dix pour cent, le seuil choisi vaut zéro et les pics 25 retenus sont toutes les raies de l'histogramme (ce cas correspond souvent à la présence d'un seul signal radar). A l'inverse lorsqu'il y a un grand nombre de raies (ce cas correspond souvent à la présence de plusieurs signaux radar), le seuil choisi est la moyenne des hauteurs des raies de l'histogramme, ou alpha fois cette moyenne, avec alpha supérieur ou égal à 1. 30 Dans une variante de traitement, les pics retenus dans l'histogramme sont des groupes de raies voisines de l'histogramme dont la hauteur dépasse le seuil calculé comme précédemment et entourant un maximum local de hauteur dans l'histogramme. Chaque pic d'histogramme correspond à une valeur de phase d'un moment de la séquence d'impulsions. 3051611 9 L'étape 50 est suivie des étapes 60 et 70 d'estimation de la valeur de phase pour chaque moment de phase, qui sont itérées pour chaque pic trouvé à l'étape 50. Pour un moment courant m(p) traité, on effectue une extraction (étape 60) des impulsions isophases, c'est-à-dire des impulsions dont la phase calculée est proche de la 5 phase associée au moment courant m(p). En effet, en pratique, à cause de la quantification du temps mesuré, les valeurs de phases calculées sont susceptibles d'être légèrement différentes pour les impulsions d'un même moment de phase sur des intervalles de répétition de motif successifs. De plus, les valeurs de phases peuvent être affectées par un jitter, avec bruit non 10 cumulatif ou cumulatif. Dans un mode de réalisation, à l'étape 60 on sélectionne les impulsions situées autour du pic d'histogramme traité, dans une fenêtre de taille déterminée par l'écart type du jitter, par exemple une taille de fenêtre proche de trois fois l'écart-type de jitter. L'étape 60 d'extraction des impulsions isophases est suivie d'une étape 70 15 d'estimation précise de la valeur de phase associée au moment m(p). Lors de l'étape 70, pour une séquence d'impulsions composée d'une pluralité P de trains d'impulsions imbriqués, on utilise la période de répétition de motif 51: , issue de l'étape 20, pour chaque moment m(p) traité :Âni(p) = . Si la séquence d'impulsions initiale est formée d'une pluralité de trains 20 d'impulsions qui se succèdent dans le temps, chaque moment m(p) étant associé à un train d'impulsions distinct, lors de l'étape 70 on effectue d'abord l'estimation précise de la période de répétition de motif /Îm(p) associée au moment m(p), ainsi que de la variance associée var(/ m(p)) . Les calculs effectués sont analogues à ceux de l'étape 24, sur l'ensemble des 25 impulsions extraites à l'étape 60. Ensuite, pour tout type de séquence d'impulsions, on applique les calculs suivants pour estimer la phase An(p),, associée au moment m(p), en utilisant la date d'arrivée de référence ta. Le calcul est effectué sur les impulsions sélectionnées à l'étape 60, qui sont au 30 nombre de N(p) impulsions. On calcule la somme vectorielle suivante, permettant d'obtenir un vecteur résultant S(p) N(p) (EQ 10) S(p) = EexpVetkp),i; (p),o k=1 3051611 10 Où tkp est la date d'arrivée de la kième impulsion du moment m(p), et les valeurs 0(tkp)Â sont calculées selon l'équation (EQ 1) avec la valeur de période de répétition m(pm. de motif 2n1(p) estimée. L'estimation de la phase ,ûrn(p),0 associée au moment m(p) est donnée par : 5 = arg[S(p)] (EQ 11) Où arg[S(p)] est l'argument du vecteur S(p). La variance associée à cette estimation est donnée par : var(A',(,),,0 ) = -21n ( 1 S(p)1 (EQ 12) N(p) Où S(p) est la norme de la somme vectorielle S(p) et In(x) est le logarithme 10 népérien de x. De retour à la figure 2, l'étape 70 d'estimation précise de la phase du moment m(p) est suivie d'une étape 80 de caractérisation de la séquence d'impulsions initiale par l'enregistrement d'un résumé comprenant des paramètres caractéristiques de cette séquence : 15 - Dans le cas où la séquence d'impulsions traitée est constituée de trains d'impulsions imbriqués, le résumé comporte : o Le nombre de moments P ; o Les valeurs calculées de la période de répétition de motif , et la variance associée ; 20 o Pour chaque moment, la valeur de phase estimée `"m(p),,,, , et la variance associée. - Dans le cas où la séquence d'impulsions traitée est constituée de trains d'impulsion successifs, le résumé comporte : o Le nombre de moments P, 25 o Pour chaque moment : - la valeur calculée de la période de répétition de motif Âm(p) , la variance associée, - la valeur de phase estimée fini(p),to , et la variance associée. Ces résumés permettant de caractériser la séquence d'impulsions traitée sont 30 mémorisés, et sont exploitables à une étape d'exploitation 90 pour diverses applications, 3051611 11 par exemple la comparaison avec des caractéristiques signal pré-enregistrées permettant d'identifier le mode d'émission du signal intercepté. Bien évidemment, les résumés sont exploitables pour toute autre application de traitement sur les impulsions radar. 5 Le procédé de l'invention est mis en oeuvre par un dispositif programmable, qui peut être soit intégré dans un dispositif de réception radar, soit indépendant et utilisé pour un traitement ultérieur. Un dispositif programmable 100 apte à mettre en oeuvre l'invention, typiquement un ordinateur, comprend une unité centrale de traitement 110, ou CPU, apte à exécuter 10 des instructions de programme informatique lorsque le dispositif 100 est mis sous tension. Le dispositif 100 comporte également des moyens de stockage d'informations 120, par exemple des registres ou des mémoires, aptes à stocker des instructions de code exécutable permettant la mise en oeuvre de programmes comportant des instructions de code aptes à mettre en oeuvre les procédés selon l'invention. 15 Dans un mode de réalisation, les instructions logicielles pour l'exécution du procédé de traitement d'un signal formé d'une séquence d'impulsions tel que brièvement décrit ci-dessus sont stockées sur un support d'enregistrement d'informations, par exemple une mémoire externe, une mémoire de stockage de masse, un disque externe. Lorsque le support d'enregistrement d'informations est connecté au dispositif 20 programmable, ce dispositf est apte à exécuter ces instructions logicielles, comme expliqué ci-dessus. En option, le dispositif programmable 100 comprend un écran 130 et un moyen 140 de saisie des commandes d'un opérateur, par exemple un clavier, optionnellement un moyen supplémentaire de pointage 150, tel une souris, permettant de sélectionner des 25 éléments graphiques affichés sur l'écran 130. Les divers blocs fonctionnels 110 à 150 du dispositif 100 décrits ci-dessus sont connectés via un bus de communication 160. En variante non représentée, le dispositif programmable 100 est réalisé sous forme de composants logiques programmables, tel qu'un ou plusieurs FPGA (de l'anglais 30 Field-Programmable Gate Array), ou encore sous forme de circuits intégrés dédiés, de type ASIC (de l'anglais Application-Specific Integrated Circuit). Avantageusement, le procédé proposé permet de traiter tout type de séquence d'impulsions, qu'il s'agisse de trains d'impulsions imbriqués ou successifs. Avantageusement, le procédé proposé permet d'obtenir des résumés fiables y 35 compris dans le cas de pertes aléatoires d'impulsions, et en plus, le niveau de fiabilité de 3051611 12 chaque estimation est quantifiable grâce à la variance d'estimation associée qui est également calculée et mémorisée. REVENDICATIONS1.- Procédé de traitement d'un signal formé d'une séquence d'impulsions, comportant au moins un motif répétitif formé d'au moins une impulsion, ledit motif étant répété dans le signal avec une période de répétition de motif, caractérisé en ce qu'il comporte les étapes suivantes : - estimation (20) de la période de répétition de motif dudit signal, - calcul (30) d'une séquence de phases en fonction d'une date d'arrivée de chaque impulsion par rapport à une date d'arrivée de référence choisie et de la période de répétition de motif estimée, - à partir de ladite séquence de phases calculées, estimation (60, 70) d'au moins une valeur de phase et d'un écart-type associé, ladite valeur de phase étant associée à un moment de phase représentatif du motif répétitif ; - obtention (80) et exploitation (90) de paramètres de caractérisation du signal en utilisant les valeurs de phase estimées. 2.- Procédé selon la revendication 1, comprenant, après l'étape de calcul (30) d'une séquence de phases, une étape de construction (40) d'un histogramme de phases, et d'extraction d'un ou plusieurs pic(s) de l'histogramme, chaque pic étant associé à un moment de phase de la séquence d'impulsions. 3.- Procédé selon la revendication 2, dans lequel l'étape d'estimation d'une phase associée à un moment de la séquence d'impulsions, ledit moment étant associé à un pic de l'histogramme de phases, comprend une extraction d'un sous-ensemble de dates d'arrivée autour dudit pic, et une estimation précise de la valeur de phase associée audit moment, et d'une variance associée à la valeur de phase estimée. 4.- Procédé selon la revendication 3, dans lequel, lorsque la séquence d'impulsions comprend une pluralité de trains d'impulsions imbriqués, chaque train d'impulsions ayant un moment de phase associé, l'estimation précise de la phase associée au moment comprend une sommation vectorielle de vecteurs unités, chaque vecteur ayant une orientation égale à une valeur de phase calculée, avec la période de répétition de motif dudit signal préalablement estimée, à partir d'une date d'arrivée dudit sous-ensemble, ladite sommation vectorielle fournissant un vecteur résultant. 3051611 14 5.- Procédé selon la revendication 4, dans lequel la valeur de phase associée audit moment estimée est égale à l'argument dudit vecteur résultant. 6.- Procédé selon la revendication 3, dans lequel, lorsque la séquence 5 d'impulsions comprend une pluralité de trains d'impulsions successifs, chaque train d'impulsions ayant un moment de phase associé, l'estimation précise de la phase associée à un moment traité comprend en outre une estimation d'une période de répétition de motifs correspondant audit moment à partir du sous-ensemble de dates d'arrivée considérées. 10 7.-Procédé selon l'une quelconque des revendications 1 à 6, dans lequel l'estimation de la période de répétition de motif du signal comporte une première estimation grossière d'une première période de répétition de motif, et une deuxième estimation raffinée d'une deuxième période d'estimation d'une première période de 15 répétition de motif. 8.- Procédé selon l'une quelconque des revendications 1 à 7, dans lequel lors de l'étape d'obtention d'une séquence de phases, pour chaque date d'arrivée d'impulsion tn, on applique un repliement de période égale à la période de répétition de motif estimée 20 selon la formule suivante : 0(t,,),1,to = (tn tn [ jy Où Lx] est la partie entière inférieure du nombre réel x, tO est la date d'arrivée de référence et X la période de répétition de motif estimée. 25 9.- Dispositif de traitement d'un signal formé d'une séquence d'impulsions, comportant au moins un motif répétitif formé d'au moins une impulsion, ledit motif étant répété dans le signal avec une période de répétition de motif, caractérisé en ce qu'il comporte au moins un processeur de calcul apte à mettre en oeuvre : - une unité d'estimation de la période de répétition de motif dudit signal, 30 - une unité de calcul, en fonction d'une date d'arrivée de chaque impulsion par rapport à une date d'arrivée de référence choisie, et de la période de répétition de motif estimée, d'une séquence de phases, 3051611 15 -une unité d'estimation, à partir de ladite séquence de phases calculées, d'au moins une valeur de phase et d'un écart-type associé, ladite valeur de phase étant associée à un moment de phase représentatif du motif répétitif ; - une unité d'obtention et d'exploitation de paramètres de caractérisation du signal 5 numérique en utilisant les valeurs de phase estimées. 10.- Dispositif selon la revendication 9 comprenant en outre une unité de construction d'un histogramme de phases, et d'extraction d'un ou plusieurs pic(s) de l'histogramme, chaque pic étant associé à un moment de phase de la séquence 10 d'impulsions. 15