FR2885225A1 - Detecteur sequentiel markovien - Google Patents
Detecteur sequentiel markovien Download PDFInfo
- Publication number
- FR2885225A1 FR2885225A1 FR0504390A FR0504390A FR2885225A1 FR 2885225 A1 FR2885225 A1 FR 2885225A1 FR 0504390 A FR0504390 A FR 0504390A FR 0504390 A FR0504390 A FR 0504390A FR 2885225 A1 FR2885225 A1 FR 2885225A1
- Authority
- FR
- France
- Prior art keywords
- criterion
- detection
- signal
- data
- state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/66—Sonar tracking systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
- Peptides Or Proteins (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Medicines Containing Material From Animals Or Micro-Organisms (AREA)
- Investigating Or Analysing Materials By The Use Of Chemical Reactions (AREA)
- Transmission And Conversion Of Sensor Element Output (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
La présente invention concerne le problème général de la détection en sonar ou en radar Elle s'intéresse en particulier à la détection initiale d'un signal d'intérêt dans un flot de données correspondant au signal reçu par un tel équipement. Cette détection initiale est transmise aux équipements chargés d'effectuer le pistage des objets d'intérêt ayant fait l'objet d'une détection.Le procédé selon l'invention considère qu'un signal d'intérêt est caractérisé par un modèle d'état x(t) à valeurs dans un ensemble fini {x1, ..., xN}, auquel correspond un ensemble fini {omega1, ..., omegaN} de N valeurs d'une variable d'observation omega, caractéristique de ce signal.Pour détecter la présence d'un signal d'intérêt le procédé calcule simultanément de deux façons différentes la valeur d'un critère de détection An(t):- comme une fonction du produit de la probabilité Pn de l'état xn, et du rapport des probabilités conditionnelles p1 et p0 de la donnée y(t, omegan) contenant le signal, en l'absence et en présence de signal,- comme une fonction du produit d'une des N valeurs du critère calculées à l'instant précédent t - 1, de la probabilité de transition Pn,m entre l'état xm pour lequel cette valeur a été calculée à l'instant t-1 et l'état xn considéré, et du rapport des probabilités conditionnelles p1 et p0 de la donnée y(t, omegan) en l'absence et en présence de signal,la plus grande des deux valeurs calculées est attribuée au critère de An(t) qui est alors comparé à un seuil de détection.
Description
Détecteur séquentiel markovien.
DOMAINE DE L'INVENTION La présente invention concerne le problème général de la détection en sonar ou en radar Elle s'intéresse en particulier à la détection initiale d'un signal d'intérêt dans un flot de données correspondant au signal reçu par un tel équipement, encore appelé senseur.
CONTEXTE DE L'INVENTION - ART ANTERIEUR La fonction générale de Détection et de Pistage Automatique (ADT selon la dénomination anglo-saxonne) est une fonction couramment utilisée dans les sonars et les radars. Cette fonction a pour objet de traiter le flot de données reçu par le système considéré de façon à détecter la présence de signaux d'intérêt et à analyser l'évolution du contact à l'origine de ce signal. Cette évolution se produit généralement dans le domaine spatial et/ou dans le domaine fréquentiel. Les données analysées peuvent être représentées par des fonctions de la forme y(t, co) dépendant d'une variable de temps t et d'une variable d'observation w qui peut être par exemple la fréquence du signal reçu ou encore la direction d'arrivée (gisement) du signal.
De manière générale la fonction ADT opère sur un contact en deux phases successives.
Tout d'abord, la fonction ADT détermine pour chaque instant t considéré et 25 pour chaque valeur de la variable d'observation w, si un signal émanant d'un objet d'intérêt est présent ou non. C'est la phase de détection.
Ensuite, lorsqu'une décision positive a été prise pour un instant tdet, la fonction ADT procède à l'évaluation de la valeur w(t) de la variable d'observation w, et éventuellement de l'état x(t), de l'objet détecté, pour les instants t ultérieurs (t tdet). C'est la phase de pistage. La phase de détection est également appelée phase d'initialisation de piste.
La phase d'initialisation et la phase de pistage constituent toutes deux des processus d'association de données au cours du temps. Mais, alors que le pistage travaille seulement sur un nombre restreint de données aux alentours des pistes, l'initialisation de piste est réalisée sur la totalité des données présentes. C'est pourquoi l'initialisation de piste (la détection) est réalisée par des traitements simples à mettre en oeuvre, alors que le pistage est généralement mis en oeuvre au moyen d'algorithmes de traitement puissants, nécessitant une grosse charge de calcul, tel que des algorithmes basés sur l'exploitation de modèles d'états.
Ceci constitue en soi une situation peu satisfaisante car la qualité de la détection souffre bien évidemment de la relative simplicité des algorithmes d'initialisation de piste.
Dans la plupart des sonars et des radars existants, la phase d'initialisation est basée sur un algorithme mettant en oeuvre de façon plus ou moins directe un principe connu de détection: la détection du type "P parmi N". L'entier N représente ici la durée de la fenêtre de test de détection le long d'un axe des temps discrétisé. Selon ce principe connu, la détection est effectuée en prenant en compte, pendant la durée de la fenêtre, le nombre d'instants où survient un événement répondant à un critère simple déterminé. La détection peut donc, par exemple, consister dans le calcul du nombre n d'instants pendant l'intervalle de temps de la fenêtre, où au moins un échantillon de signal, inclus dans la fenêtre, excède un seuil de détection r1, et dans la comparaison de ce nombre n à un nombre P (P N). Ainsi, un signal est détecté quand il y a de tels événements au moins P fois dans la fenêtre (n P).
Dans certaines applications connues, ce principe général est affiné par ajout de la condition supplémentaire que la valeur moyenne des P plus grands événements doit être supérieure à un deuxième seuil de détection r2, lui-même supérieur au premier seuil r1.
De manière connue, la probabilité de fausse alarme des procédés mettant en oeuvre ce principe connu, dépend de ri, de r2 et de la taille de la fenêtre de test dans le domaine des données d'observation w. En particulier, l'augmentation de la taille de la fenêtre de test M) dans le domaine des données d'observation conduit à l'augmentation de la probabilité de fausse alarme, mais permet de mieux s'affranchir des pertes de détection liées à la dérive au cours du temps du signal dans le domaine d'observation.
Dans la pratique, ces données contradictoires conduisent à rechercher des compromis plus ou moins satisfaisants. Ainsi, pour assurer une probabilité de détection convenable sans trop favoriser l'apparition de fausses alarmes, Aco est généralement choisi de façon à être, au plus, égal à la résolution du senseur et la durée N du test est choisie suffisamment petite pour que la dérive du signal pendant la durée de la mesure ne puisse pas excéder M). De plus, la durée N de la fenêtre ne doit pas excéder beaucoup la durée du signal le plus court à détecter, sinon la probabilité de détection de ce type de signaux, pour une probabilité de fausse alarme fixée, diminue.
Ce compromis conduit en particulier à choisir une durée N limitée à quelques unités, contrainte qui empêche de tirer pleinement parti de l'accroissement potentiel de la probabilité de détection sur les signaux de longue durée. Ainsi, quelle que soit la variante du principe de détection "P parmi N" utilisée celle-ci conduit toujours à limiter la durée N de la fenêtre de test utilisée, N étant généralement inférieur à 10..
La mise en oeuvre d'une fonction ADT peut également être réalisée en utilisant un procédé connu consistant à former avant détection des pistes d'essai, appelées pistes candidates, pour lesquelles on effectue des tests de validation. Lorsqu'une piste candidate passe positivement son test de validation, une détection est signalée et la piste devient une piste validée, gérée ultérieurement par le pistage.
Dans cette approche, connue sous l'appellation anglo-saxonne "Track before detect" (TBD), l'initialisation de piste consiste donc à former et tester des pistes candidates sur les données reçues. Les données reçues sont intégrées le long de ces pistes sélectionnées conformément à un modèle déterminé de dynamique, attribué au type d'objet que l'on souhaite détecter.
L'utilisation d'un tel modèle permet d'intégrer ces données sur une durée de test plus longue que dans le cas de l'emploi du critère "P parmi N". Cette approche est donc particulièrement adaptée pour traiter des signaux ayant un faible rapport signal sur bruit. En revanche, comme pour la mise en oeuvre du critère "P parmi N", le traitement est ici le plus souvent réalisé sur des blocs de données de durée fixe. On doit alors soit faire un compromis entre la détection des signaux courts et celle des signaux longs en fixant une durée intermédiaire, soit faire en parallèle plusieurs traitements ayant des durées d'intégration différentes, ce qui multiplie par autant le coût de calcul.
PRESENTATION DE L'INVENTION Un des buts de l'invention est de proposer un procédé de détection permettant en particulier d'effectuer une initialisation de piste de type "Track before Detect" dont la durée d'intégration n'est pas fixe mais adaptée pour tirer parti de la durée du signal d'intérêt reçu, de façon à obtenir une probabilité de détection élevée tout en maintenant une probabilité de fausse alarme faible.
A cet effet l'invention a pour objet un procédé continu de détection d'un signal s(t) présent dans un flot de données y(t, w) dépendant d'une variable de temps discrète t et d'une variable d'observation w, le signal étant caractérisé par un modèle d'état x(t) à valeurs dans un ensemble fini {xi, ..., xN}, auquel correspond un ensemble fini {col, ..., wN} de N valeurs de la variable d'observation (o.
Le procédé selon l'invention élabore pour chaque instant t et pour chaque état xn un critère de détection An(t), dont la valeur à chaque instant t, pour un 25 état xn, est calculée simultanément de deux façons différentes: - comme une fonction du produit de la probabilité Pn de l'état xn, et du rapport des probabilités conditionnelles pi et po de la donnée y(t, wn) en l'absence et en présence de signal, - comme une fonction du produit d'une des N valeurs du critère calculées à l'instant précédent t 1, de la probabilité de transition Pn,m entre l'état xm pour lequel cette valeur a été calculée à l'instant t 1 et l'état xn considéré, et du rapport des probabilités conditionnelles pi et po de la donnée y(t, (On) en l'absence et en présence de signal.
Les résultats des deux calculs sont comparés et la plus grande des deux valeurs calculées est attribuée au critère de An(t).
Selon l'invention, le procédé comporte au moins: - une étape A d'initialisation où les N valeurs du critère de détection An à l'instant to 1 précédant le premier instant des données à traiter to sont fixées à des valeurs prédéterminées, - une étape B de calcul des N valeurs du critère An à l'instant courant 10 t, - une étape C de détection par comparaison à un seuil des N valeurs du critère déterminées à la phase B, le dépassement du seuil de détection par une valeur du critère déterminant la présence d'un signal dans l'état xn correspondant à cette valeur du critère, à l'instant courant.
Les étapes B et C sont réalisées de manière itérative jusqu'à ce que le critère de détection An dépasse le seuil de détection fixé à l'instant tdet. Une piste est alors créée avec l'état x à cet instant. Les valeurs du critère de détection Am(tdet) pour les états xm dans un voisinage de xn, qui sont influencées par le maximum atteint à l'état xn, sont en outre mises à zéro afin d'empêcher que ces valeurs proches du seuil de détection puissent causer l'initialisation d'une nouvelle piste à l'instant suivant tdet +1.
Selon l'art connu, le pistage, qui consiste à mettre à jour l'estimation de l'état x(t) de chaque piste déjà initialisée, est effectué à chaque instant avant la mise en oeuvre du procédé d'initialisation de piste de la présente invention. Les données associées aux pistes par le pistage sont remplacées par la valeur moyenne du bruit, afin que ces données ne puissent favoriser l'initialisation de nouvelles pistes. Le traitement d'initialisation de piste de l'invention peut être couplé selon ce principe avec un traitement de pistage quelconque.
Selon l'invention, le critère de détection An est défini par les relations générales: 6 p1(Y(to,wn)) \K1 po(Y(to,wn))Pn ' et
A
(t) = max F K2 p1(y(t, wn)) max {P Am (t -1)} , F/K Ap1(yt, w" )) po (Y(t, on)) 1<m<N ne po (Y(t, wn)) i où to représente un instant initial, où p1 et po représentent les lois de probabilité conditionnelles des données y(t, w) en l'absence et en présence de signal, où F est une fonction croissante et K1 et K2 sont deux constantes. Le procédé selon l'invention met avantageusement en oeuvre un algorithme itératif basé sur l'exploitation d'un modèle de Markov caché (HMM selon la terminologie anglo-saxonne) auquel il associe une opération de réinitialisation automatique du critère An conditionnée par la comparaison entre la probabilité initiale Pn et la quantité: K. max{Pn m É Am(t -1)}. 1< m<N
Selon une variante de mise en oeuvre du procédé, le facteur K1 ou le le facteur K2 sont régulés de façon à maintenir un taux de réinitialisation du calcul récursif de An, voisin d'une valeur spécifiée R. Le taux de réinitialisation est déterminé par la condition F/K1. p1(Y(t,wn))pnl> F"K2 P1(Y(t,wn)) max{Pn,m'Am(t 1)}I, po(Y(t,(on)) J po(Y(t,wn)) 15m<N l Selon une variante de mise en oeuvre du procédé R est choisi égal à 0,5.
Le procédé selon l'invention peut avantageusement être utilisé pour l'initialisation de pistes sonar sur des raies spectrales à partir de données 25 d'observation y(t, w) pour lesquelles co est une donnée fréquentielle.
Le procédé selon l'invention peut également avantageusement être utilisé pour l'initialisation de pistes sonar sur des raies spectrales à partir de données d'observation y(t, w) pour lesquelles w est un couple de données (fréquence, gisement).
An (to) = F
DESCRIPTION DES FIGURES
D'autres caractéristiques et avantages apparaîtront au travers de la description qui suit et de la figure 1 annexée qui représente un 5 organigramme de principe du procédé selon l'invention.
DESCRIPTION DETAILLEE
Nous présentons dans la suite de la description le principe de détection utilisant un modèle de Markov caché, sur lequel est basé le procédé de détection selon l'invention. Puis le principe de détection selon l'invention lui-même.
La détection au moyen d'un modèle de Markov caché est basée sur les 15 considérations suivantes: Un objet d'intérêt que l'on essaie de détecter au travers d'un signal s(t) le long d'un chemin (Ds(t) dans l'espace d'observation a une évolution x(t) dans l'espace des états au cours du temps qui n'est a priori pas connue. Cette évolution peut être modélisée par un modèle de processus markovien pouvant prendre à chaque instant t d'observation, une valeur comprise dans un espace d'états E comportant N états possibles {xi, ..., xN}. A chaque état xn correspond une valeur con de la variable d'observation w. Cette observation peut par exemple porter sur une variable unique telle que l'ensemble des composantes fréquentielles du signal reçu ou encore sur plusieurs variables associées comme la fréquence et la direction d'observation (gisement).
On appelle piste candidate la suite d'états x(t) que peut prendre l'objet au cours de son évolution pendant un temps compris entre un instant initial to et un instant to+At. On note cette suite X(to,At) = (x(t))t-to, K,to+At Les N probabilités Pn = P[x(to) = xr] pour chacun des états xn de l'espace des états de correspondre à l'état initial x(to) ainsi que les N2 probabilités Pn,m = P[x(t) = xn 1 x(t-1) = xm] associées aux différentes transitions possibles entre états, constituent des paramètres connus associés au modèle de Markov caché choisi.
Dès lors, la probabilité a priori pour que le comportement de l'objet suive une 35 piste candidate quelconque X(to, At), est égale au produit de la probabilité de l'état x(to) par les probabilités des transitions entre les états successifs x(t-1) et x(t) pour to+1 <- t <- to+At.
On note wx(t) la valeur de la variable d'observation correspondant à l'état x(t). 5 La suite de données Yx (to, At) _ (y(t, wx (t)))t=to,K, to +At correspondant à la piste candidate X(to, At) est modélisée par un processus aléatoire indépendant et on note po et pi ses densités de probabilité respectivement en présence de bruit uniquement et en présence de signal et de bruit: po(y(t,(ox(t)))=P[y(t,wx(t))I Ho] [1] Pi(y(t,wx(t)) )=P[y(t,wx(t))I x(t)] [2] où Ho représente l'hypothèse d'une absence de signal correspondant à un 15 objet d'intérêt dans le flot de données traité.
Il résulte du caractère markovien du processus aléatoire d'état x(t) et de l'indépendance du processus aléatoire de données y(t,wx(t)) que la probabilité conjointe de X(to, At) et Y(to, At) a pour expression: P[X(t0, At),Yx(t0,At)]= P[y(to +At,wx(to +At)) X(t o, At)]. P[Yx (t o, At 1) IX(to,At)].
P[x(to + At) l x(to + At -1)] É P[X(to, At -1)]. [3] Dans un modèle de Markov caché, on fait de plus l'hypothèse que la dépendance d'une donnée y(to + At, wx(t)) par rapport au processus d'état est sa dépendance par rapport à l'état x(to + At) seulement. Par suite la probabilité conjointe de X(to, At) et Y(to, At) prend la forme récursive suivante: P[X(t0,At), Yx(to,At),]= p1(y(to +At,wx(t0 +At))).P[x(t0 +At) Ix(to +At 1)].
P[X(to,At 1),Yx(to,At 1)] . [4] La probabilité de la suite Yx(to, At) dans l'hypothèse Ho étant égale à t=to+At P[Yx (to, At) I H0]= fl po (y(t, (ox (t))), le rapport de vraisemblance Ax,y(to, t=to At) de l'association de l'observation Yx(to, At) et du chemin modèle X(to, At) est alors déterminé par les relations suivantes: p1Y( o, . P[x(t0)] [5] ox({0))) Y(t0,0)=
AX
Po(y(o, x(0))) P1(Y(t0 +At,cox(t0 +At))) nXY(tO,ot = po(Y(t0 +At, wx(t0 + At))) [6] P[x(to +At)Ix(to +At 1)]É AX,y(to,At 1).
On rappelle ici qu'à t = to+At, x(t) peut a priori prendre n'importe quelle valeur 1 o xn de l'espace des états.
II est par suite possible de déterminer pour tous les chemins X(to, At) pour lesquels x(t) prend la valeur xn à l'instant t = to+At, le chemin pour lequel le critère de détection prend sa valeur maximum, notée A(xn, to, At). Cette détermination peut par exemple être faite en utilisant un algorithme itératif tel que l'algorithme de Viterbi par exemple. A(xn, to, At) sera alors calculé au moyen des relations suivantes: A x,t 0) = p1(Y(t0,n))P ( n 0, P0(Y(t0, n)) n A(x,t At) p1(y(to+At,n)) max {P A(x t, At 1)} A(xn 0, P0(Y(t0 +Ot,wn))1_m_N n m m 0 Ainsi, au bout du temps d'observation At on connaît l'ensemble des chemins qui maximisent, pour chaque valeur possible de x(to+At), le critère de détection Ax,y(to, At). La détection peut alors être achevée en ne retenant de ces chemins que ceux pour lesquels la valeur A(xn, to, At) est supérieure à un seuil donné.
Chaque détection peut ensuite être utilisée pour initialiser une piste qui sera traitée par ailleurs.
[7] [8] Comme on peut donc le constater, l'utilisation d'un modèle de Markov caché pour rendre compte de l'évolution des paramètres d'un objet d'intérêt permet de réaliser une détection de cet objet par analyse du rapport de vraisemblance de seulement N hypothèses d'évolution des paramètres de l'objet à chaque instant de la mesure, ces hypothèses d'évolution se traduisant par N hypothèses sur l'évolution o (t) de la variable d'observation du signal.
En comparant à un seuil les valeurs de A(xn, to, At) pour 1 n N, il est possible d'exécuter un test de détection portant sur l'ensemble des données y(t,co) présentes pendant l'intervalle de temps At. Si le test est positif à l'état xn, alors le flot de données analysé pendant le temps At est considéré comme incluant le signal d'un objet à détecter ayant l'état xn à l'instant to + At Avantageusement, parmi les tests de détection connus qui opèrent sur des données reçues dans un intervalle de temps fixé, ce test de détection est celui qui permet de maximiser la probabilité de détection d'un objet d'intérêt pour une probabilité de fausse alarme donnée. En revanche comme tous les tests portant sur une fenêtre de durée finie, l'utilisation d'un modèle de Markov caché associé à un critère de détection tel que celui défini par les relations 7 et 8, peut être amélioré si l'on adapte la durée de la fenêtre d'observation pour tenir compte de la durée du signal d'intérêt. Pour ce faire, il est connu de réaliser le traitement décrit pour toutes les valeurs possibles de tO et At, ou du moins en pratique, pour un certain nombre de valeurs choisies judicieusement. Ce choix résulte généralement d'un compromis à trouver entre le coût de calcul et la performance de détection.
Afin d'éviter le recours à un compromis sur le choix de la taille de la fenêtre d'observation, le procédé selon l'invention met en oeuvre, à partir d'un modèle de Markov caché, un critère de détection modifié. Ce nouveau critère de détection est calculé de manière continue depuis un instant d'origine to jusqu'à ce qu'une détection se produise. Comme précédemment cette détection est caractérisée par le dépassement par le critère de détection, d'un seuil donné. Le procédé selon l'invention combine ainsi le principe de la détection séquentielle avec celui de la détection à partir d'un modèle markovien caché.
Selon l'invention, la valeur maximale pour xn, à l'instant t, du critère de 5 détection n'est plus déterminé à partir des relations [7] et [8] de l'algorithme de Viterbi, mais, par exemple, par les relations suivantes: p1(Y(to,wn)) A (to) = P n p0(Y(t0,wn)) An (t) = p1 (Y(t, wll max K É max { Pn,m Am(t -1)}, Pn [10] p0(Y(t+ wnl 1<m<N Dans la relation [10], K correspond à un facteur constant inférieur ou égal à 1 Comme on peut le constater, la relation [9] correspond à un cas particulier d'expression de la relation [10] pour lequel An(to-1) = 0 pour toutes les valeurs de n.
Comme pour une détection séquentielle classique, le calcul du critère An ne requiert pas de fixer une durée At de valeur donnée.
Selon l'invention, le calcul du critère de détection An selon la relation [10] fait intervenir un facteur A = max K max { Pn m É Am(t -1)}, Pn. 1<m<N
Le facteur A de la relation [10] diffère du facteur correspondant contenu dans la relation [8]. Ce facteur A selon l'invention permet avantageusement de pouvoir réaliser le test des données d'observation de manière continue. L'intervalle de temps pendant lequel le test est réalisé n'est donc plus ici un intervalle de durée fixe At mais un intervalle permettant de prendre en compte la durée réelle du signal.
Le facteur A permet également avantageusement de prendre en compte le fait qu'un signal d'intérêt peut n'être présent que pendant une fraction de la fenêtre d'observation et permet de ne pas intégrer au test des données [9] correspondant au signal d'intérêt, des données présentes pendant la durée de la fenêtre d'observation et correspondant seulement à du bruit.
Comme on peut le constater le critère An(t) selon la relation [10], prise comme exemple non limitatif, est lié au critère A(xn, to, At) de l'algorithme de Viterbi par la relation suivante: An (t) = Kt-t A(xn, to, t to) [Il] où to représente l'instant de la dernière réinitialisation de An(t). Selon l'invention, le facteur K est choisi de façon à ce que la relation suivante soit vérifiée: P[K1maN{Pn,m'Am(t 1)}>Pn I Ho] = 2 m où Ho est l'hypothèse qu'il n'y a pas de signal d'intérêt dans le flot de données observé.
Cette relation impose que les probabilités qu'une piste candidate X(to, t to) soit abandonnée ou continue à être explorée soient égales en l'absence de 20 signal.
Le critère de détection selon l'invention, An(t), peut être calculé à partir des relations récursives [9] et [10]. Ce calcul peut être effectué de manière similaire à l'algorithme de Viterbi, mis à part le processus de réinitialisation.
Le critère de détection An(t) tel que décrit dans ce qui précède, est avantageusement mis en oeuvre sur les données y(t, a') par le procédé de détection selon l'invention, selon un schéma général de fonctionnement semblable à celui illustré par la figure 1.
Durant la phase d'initialisation A, on calcule les valeurs du critère pour les N états selon la relation [10] en prenant An(to-1) = 0 pour toutes les valeurs de n.
Durant la phase B qui suit, le critère de détection est calculé pour les N états aux instants successifs selon la relation [10]. Ce calcul est effectué de [12] manière récursive jusqu'à ce que, pour un instant donné t, An(t) dépasse le seuil de détection pour au moins une valeur de n. L'apparition d'un tel dépassement cause la création d'une piste validée, dont l'état à cet instant est celui où est apparu le dépassement du seuil de détection, et la remise à zéro de la valeur du critère pour ce même état ainsi que pour les états voisins. Le calcul récursif de la phase B peut alors reprendre.
Les relations [9] et [10] donnent à titre d'exemple, une expression particulière du critère de détection mis en oeuvre par le procédé selon l'invention. Plus généralement, le critère An(t), peut être défini par les relations suivantes: An(to 1) = 0, [13] Pl Mt' 'n p 15 An(t)=max F[KZ Po(Y(twn 1mm {Pn,m'nm(t 1)}l, FI K1'Po(Y(t,(0n n [14] où F est une fonction croissante, KI et K2 sont deux constantes, éventuellement régulées en fonction des valeurs du critère de détection pour les N états. On peut ainsi, par exemple, fixer KI et K2 de façon à maintenir le taux de réinitialisation du critère proche d'une valeur fixée R, égale par exemple à 0,5.
Les relations [9] et [10] sont un cas particulier des relations [13] et [14] où F est la fonction identité, KI étant égal à let K2 à K. Parmi les autres choix possibles de la fonction F, il y a notamment les fonctions logarithmes qui rendent le traitement totalement additif.
Le synoptique de la figure 1 résume de manière schématique le 30 procédé de détection selon l'invention.
Ce procédé comporte principalement une phase d'initialisation A, une phase B de calcul de la valeur d'un critère de détection par implémentation d'un algorithme de calcul récursif, et une phase C de test du critère calculé par rapport à un seuil de détection.
Durant la phase A, le critère de détectionAn(t) est initialisé pour un instant to-5 1 donné, en posant An(to-1) = 0 pour toutes les valeurs de n.
Ensuite, durant la phase B, pour chaque instant t to, le calcul du critère de détection An(t) est effectué de manière itérative pour chaque état xn, par deux opérations de calcul 11 et 12 menées simultanément et correspondant à la relation [14]. Les résultats sont comparés et la valeur maximale pour chaque état n est conservée A l'issue de chaque itération, durant la phase C, les valeurs A,(t), ..., AN(t) conservées sont comparées au seuil de détection.
Si le test est négatif pour tous les états xn la boucle d'itération se poursuit à l'instant suivant t+1 par un nouveau passage par la phase B. Si en revanche le test est positif pour un état xn, une détection est signalée donnant lieu à la création d'une piste validée à l'état xn pour l'instant t. Les valeurs du critère pour l'état xn et les états voisins font alors l'objet d'une remise à zéro 13. Puis, la boucle d'itération se poursuit à l'instant suivant (phase B).
Le procédé de détection selon l'invention peut avantageusement prendre place dans un grand nombre d'applications. On peut en particulier citer à titre d'exemple la détection sonar d'objets dont l'évolution est observée au travers de la variation au cours du temps de la fréquence f du signal reçu et de la pente de variation f' de cette fréquence, encore appelée "frequency slope" selon la dénomination anglo-saxone.
Dans une telle application la variable d'état x(t) est définie par la paire {f(t), f'(t)}. Le procédé selon l'invention est alors appliqué aux données observées qui représentent par exemple l'amplitude quadratique des composantes spectrales d'un signal, obtenues par transformée de Fourier rapide (FFT) opérant sur une fenêtre glissante et ont pour expression 2885225 15 y(t, f) = 1 S(t,f) 12. La variable d'observation co est dans ce cas la fréquence f des raies spectrales obtenues par FFT.
Dans le cas où la composante de bruit du signal A(t,f), dont y(t,f) représente le module au carré, est gaussienne, de moyenne nulle et de variance unité, et que la composante de A(t,f) constituée par le signal à détecter est d'amplitude constante, de rapport signal sur bruit ro, les densités po(y) et pi(y) suivent des lois de type x2 à deux degrés de liberté, respectivement centrées et non centrées, ayant pour expression: po (y) = exp( Y) [15] pi(y)=exp( y r0).'o(2 roi) [16] avec 10(z 1 f 0 ezcos(0)dO.
Soit TFT la longueur de la fenêtre temporelle glissante de la transformation de Fourier, kt et kt les coefficients tels que le pas d'échantillonnage temporel des données y (t, f) est égal TFT/kt et le pas d'échantillonnage fréquentiel des états est égal à TFT-'/kf. On fixe alors le pas de pente de fréquence des états égal à TFT-2.kt/kt, rapport du pas fréquentiel sur le pas temporel.
L'espace d'états E est alors défini comme un réseau fini, avec des mailles de dimensions (TFT1/kt, TFT-2.kt/kf). Si on n'a pas de connaissance a priori sur la fréquence ou la pente des
signaux que l'on cherche à détecter, on pose Pn = 1/N. Les probabilités de transition d'un état à l'autre peuvent être par exemple définies de la façon suivante: Si fn+f'm TFT < 1 [17] -1 fn fm 2kfTFT 2 kt 2kfTFT Alors Pn,m = h /kfT2 [18] FT Ifin_f,ml kt i Sinon Pn,m = 0, Où h représente une fonction décroissante telle que: h(0) + 2Eh(i) =1. [19] i=1 Une telle expression de la probabilité de transition traduit le fait que la probabilité pour passer de l'état (fm, fm) à l'état (fn, fn) n'est pas nulle si et seulement si la différence entre l'écart de fréquence fn- fm et la valeur moyenne (fn+fm)/2 des deux pentes multipliée par le pas temporel TFT/kt est inférieur à la moitié du pas fréquentiel TFT-'/kf. Dans ce cas, la probabilité de transition est une fonction décroissante de I fn-fm L. Pour un état (fm, fm) donné et un écart de pente (fn-fm) donné, il y a seulement une fréquence fn qui satisfait la relation [17]. La relation [19] est donc équivalente à la condition EPn,m = 1, qui exprime le fait que la somme n des probabilités conditionnelles P[x(t)=xn 1 x(t-1) = xm]des N états possibles, connaissant l'état précédent xm, est égale 1.
En pratique, dans ce cas particulier, la fonction h qui détermine les probabilités de transition est fixée de façon arbitraire car on ne dispose que rarement d'un modèle statistique de l'évolution de la fréquence des signaux à détecter. On peut cependant rappeler que plus large est le sommet de la fonction h en 0, meilleur est la capacité du traitement à faire face à la fluctuation rapide de la variation de fréquence, mais plus basse est la performance du traitement sur la détection des signaux ayant une pente de fréquence f' constante, voire nulle.
Claims (7)
1. Procédé continu de détection d'un signal s(t) présent dans un flot de données y(t, w) dépendant d'une variable de temps discrète t et d'une variable d'observation w, le signal étant caractérisé par un modèle d'état x(t) à valeurs dans un ensemble fini {xi, ..., xN}, auquel correspond un ensemble fini {col, ..., CON) de N valeurs de la variable d'observation co, ce procédé élaborant pour chaque instant t et pour chaque état xn un critère de détection An(t), caractérisé en ce qu'à chaque instant t la valeur du critère An(t) pour un état xn est calculée simultanément de deux façons différentes: o - comme une fonction du produit de la probabilité Pn de l'état xn, et du rapport des probabilités conditionnelles pi et po de la donnée y(t, con) en l'absence et en présence de signal, - comme une fonction du produit d'une des N valeurs du critère calculées à l'instant précédent t 1, de la probabilité de transition Pn,m entre l'état xm pour lequel cette valeur a été calculée à l'instant t 1 et l'état xn considéré, et du rapport des probabilités conditionnelles pi et po de la donnée y(t, con) en l'absence et en présence de signal, la plus grande des deux valeurs calculées étant attribuée au critère de An(t).
2. Procédé selon la revendication 1 comportant au moins: - une étape A d'initialisation où les N valeurs du critère de détection An à l'instant to 1 précédant le premier instant des données à traiter to sont fixées à des valeurs prédéterminées, - une étape B de calcul des N valeurs du critère An à l'instant courant t suivant les deux façons et de détermination de la valeur maximale, - une étape C de détection par comparaison à un seuil des N valeurs du critère déterminées à la phase B, le dépassement du seuil de détection par une valeur du critère déterminant la présence d'un signal dans l'état xn correspondant à cette valeur du critère, à l'instant courant, les étapes B et C étant réalisées de manière itérative.
3. Procédé selon l'une des revendications 1 ou 2, dans lequel le critère de détection An(t) est défini par les relations suivantes: An(to 1) = 0, An(t)=max F/K2. P1Y(t,wn max{Pn,m.Am(t 1)} ,F/K1. P1(Y(twn p Po(Y(t,wn) ) 1<m<N i \ Po(Y(t,wn)) n où to est le premier instant des données à traiter, con est le lieu dans l'espace d'observation correspondant à l'état xn, où pi et po représentent les lois de probabilité conditionnelles des données y(t, w) en l'absence et en présence de signal, où F est une fonction croissante, et où KI et K2 sont deux facteurs choisis arbitrairement.
4. Procédé selon la revendication 3 pour lequel le facteur KI ou le facteur K2 sont régulés de façon à maintenir un taux de réinitialisation du calcul récursif de An, déterminée par la condition FK1. P1(Y(t'wnpnl>FK2. P1(Y(t wn. max{Pn,m'Am(t 1)}J, Po(Y(t,w,, Po(Y(t, (0n)) 1<m<N voisin d'une valeur spécifiée R.
5. Procédé selon la revendication 4 pour lequel R est égal à 0,5.
6. Application du procédé selon l'une des revendications précédentes à l'initialisation de pistes sonar sur des raies spectrales à partir de données d'observation y(t, w) pour lesquelles co est une donnée 30 fréquentielle.
7. Application du procédé selon l'une des revendications précédentes à l'initialisation de pistes sonar sur des raies spectrales à partir de données d'observation y(t, w) pour lesquelles w est un couple de données (fréquence, gisement).
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0504390A FR2885225B1 (fr) | 2005-04-29 | 2005-04-29 | Detecteur sequentiel markovien |
AT06754748T ATE423982T1 (de) | 2005-04-29 | 2006-04-18 | Nachweis von markov-sequenzen |
DE602006005350T DE602006005350D1 (de) | 2005-04-29 | 2006-04-18 | Nachweis von markov-sequenzen |
EP06754748A EP1877826B1 (fr) | 2005-04-29 | 2006-04-18 | Détecteur séquentiel markovien |
CA2611978A CA2611978C (fr) | 2005-04-29 | 2006-04-18 | Detecteur sequentiel markovien |
US11/912,924 US7649806B2 (en) | 2005-04-29 | 2006-04-18 | Markov sequential detector |
PCT/EP2006/061630 WO2006117290A1 (fr) | 2005-04-29 | 2006-04-18 | Détecteur séquentiel markovien |
NO20075449A NO339653B1 (no) | 2005-04-29 | 2007-10-29 | Markov-sekvensdetektor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0504390A FR2885225B1 (fr) | 2005-04-29 | 2005-04-29 | Detecteur sequentiel markovien |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2885225A1 true FR2885225A1 (fr) | 2006-11-03 |
FR2885225B1 FR2885225B1 (fr) | 2007-06-22 |
Family
ID=35395782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0504390A Expired - Fee Related FR2885225B1 (fr) | 2005-04-29 | 2005-04-29 | Detecteur sequentiel markovien |
Country Status (8)
Country | Link |
---|---|
US (1) | US7649806B2 (fr) |
EP (1) | EP1877826B1 (fr) |
AT (1) | ATE423982T1 (fr) |
CA (1) | CA2611978C (fr) |
DE (1) | DE602006005350D1 (fr) |
FR (1) | FR2885225B1 (fr) |
NO (1) | NO339653B1 (fr) |
WO (1) | WO2006117290A1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7649806B2 (en) * | 2005-04-29 | 2010-01-19 | Thales | Markov sequential detector |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9659063B2 (en) * | 2010-12-17 | 2017-05-23 | Software Ag | Systems and/or methods for event stream deviation detection |
US9389215B2 (en) | 2011-09-23 | 2016-07-12 | Mastinc | Multi-modal fluid condition sensor platform and system thereof |
US9020766B2 (en) * | 2011-09-23 | 2015-04-28 | Mastinc. | Multi-modal fluid condition sensor platform and system therefor |
US9568461B2 (en) * | 2012-12-31 | 2017-02-14 | Mastinc | Multi-modal fluid condition sensor platform and system therefor |
US9792259B2 (en) | 2015-12-17 | 2017-10-17 | Software Ag | Systems and/or methods for interactive exploration of dependencies in streaming data |
CN107202989B (zh) * | 2017-05-08 | 2020-01-14 | 电子科技大学 | 一种适用于被动拖曳线列阵声呐的复杂弱目标检测和跟踪方法 |
CN108460773B (zh) * | 2018-02-28 | 2022-03-18 | 哈尔滨工程大学 | 一种基于偏移场水平集的声纳图像分割方法 |
CN109901154B (zh) * | 2019-03-29 | 2023-03-03 | 中国人民解放军海军航空大学 | 基于递推rtht-tbd的自适应调整方法 |
GB202001226D0 (en) | 2020-01-29 | 2020-03-11 | Givaudan Sa | Process |
WO2024096938A2 (fr) | 2022-07-22 | 2024-05-10 | Numerica Corporation | Détection multi-cible à l'aide d'une dispersion convexe antérieure |
CN115407299B (zh) * | 2022-09-15 | 2024-06-18 | 中国人民解放军国防科技大学 | 基于伯努利滤波器算法的弱小目标检测方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5900835A (en) * | 1998-07-09 | 1999-05-04 | The United States Of America As Represented By The Secretary Of The Navy | Coherent hidden markov model |
US6038526A (en) * | 1998-06-24 | 2000-03-14 | The United States Of America As Represented By The Secretary Of The Navy | Method for detecting weak signals in a non-gaussian and non-stationary background |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5832430A (en) * | 1994-12-29 | 1998-11-03 | Lucent Technologies, Inc. | Devices and methods for speech recognition of vocabulary words with simultaneous detection and verification |
US5675706A (en) * | 1995-03-31 | 1997-10-07 | Lucent Technologies Inc. | Vocabulary independent discriminative utterance verification for non-keyword rejection in subword based speech recognition |
FR2738918B1 (fr) * | 1995-09-19 | 1997-12-05 | Thomson Csf | Procede d'autofocalisation pour sonar a antenne synthetique |
FI101429B (fi) * | 1995-09-29 | 1998-06-15 | Instrumentarium Oy | Törmäyslevenemän korjaus kaasujen ei-dispersiivisessä absorptiomittauk sessa |
US6038528A (en) * | 1996-07-17 | 2000-03-14 | T-Netix, Inc. | Robust speech processing with affine transform replicated data |
FR2769372B1 (fr) * | 1997-10-07 | 1999-12-03 | Thomson Marconi Sonar Sas | Procede de correction des effets des mouvements parasites de l'antenne dans un sonar a antenne synthetique |
FR2818486B1 (fr) * | 2000-12-19 | 2004-07-23 | Thomson Marconi Sonar Sas | Antenne sonar synthetique interferometrique |
FR2821163B1 (fr) * | 2001-02-16 | 2003-05-09 | Thomson Marconi Sonar Sas | Procede de detection d'objets mobiles au moyen de bouees sous-marines passives |
FR2843461B1 (fr) * | 2002-08-06 | 2004-09-10 | Thales Sa | Antenne sonar synthetique |
FR2885225B1 (fr) * | 2005-04-29 | 2007-06-22 | Thales Sa | Detecteur sequentiel markovien |
-
2005
- 2005-04-29 FR FR0504390A patent/FR2885225B1/fr not_active Expired - Fee Related
-
2006
- 2006-04-18 US US11/912,924 patent/US7649806B2/en not_active Expired - Fee Related
- 2006-04-18 EP EP06754748A patent/EP1877826B1/fr active Active
- 2006-04-18 WO PCT/EP2006/061630 patent/WO2006117290A1/fr not_active Application Discontinuation
- 2006-04-18 AT AT06754748T patent/ATE423982T1/de not_active IP Right Cessation
- 2006-04-18 CA CA2611978A patent/CA2611978C/fr active Active
- 2006-04-18 DE DE602006005350T patent/DE602006005350D1/de active Active
-
2007
- 2007-10-29 NO NO20075449A patent/NO339653B1/no not_active IP Right Cessation
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6038526A (en) * | 1998-06-24 | 2000-03-14 | The United States Of America As Represented By The Secretary Of The Navy | Method for detecting weak signals in a non-gaussian and non-stationary background |
US5900835A (en) * | 1998-07-09 | 1999-05-04 | The United States Of America As Represented By The Secretary Of The Navy | Coherent hidden markov model |
Non-Patent Citations (2)
Title |
---|
ANGELOVA D ET AL: "Monte carlo-based filter for targettracking with feature measurement", INFORMATION FUSION, 2002. PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON JULY 8-11, 2002, PISCATAWAY, NJ, USA,IEEE, vol. 2, 8 July 2002 (2002-07-08), pages 1499 - 1505, XP010594371, ISBN: 0-9721844-1-4 * |
JOHNSTON L A ET AL: "Performance analysis of a track before detect dynamic programming algorithm", ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, 2000. ICASSP '00. PROCEEDINGS. 2000 IEEE INTERNATIONAL CONFERENCE ON 5-9 JUNE 2000, PISCATAWAY, NJ, USA,IEEE, vol. 1, 5 June 2000 (2000-06-05), pages 49 - 52, XP010507265, ISBN: 0-7803-6293-4 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7649806B2 (en) * | 2005-04-29 | 2010-01-19 | Thales | Markov sequential detector |
Also Published As
Publication number | Publication date |
---|---|
US20080186225A1 (en) | 2008-08-07 |
US7649806B2 (en) | 2010-01-19 |
WO2006117290A1 (fr) | 2006-11-09 |
CA2611978A1 (fr) | 2006-11-09 |
EP1877826A1 (fr) | 2008-01-16 |
NO20075449L (no) | 2008-01-29 |
ATE423982T1 (de) | 2009-03-15 |
FR2885225B1 (fr) | 2007-06-22 |
EP1877826B1 (fr) | 2009-02-25 |
CA2611978C (fr) | 2014-06-03 |
NO339653B1 (no) | 2017-01-16 |
DE602006005350D1 (de) | 2009-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1877826B1 (fr) | Détecteur séquentiel markovien | |
EP3155608B1 (fr) | Procedé de suivi d'une partition musicale et modélisation associée | |
EP2696344B1 (fr) | Procede et systeme pour detecter des evenements sonores dans un environnement donne | |
EP2100161B1 (fr) | Procede de traitement radar passif multivoies d'un signal d'opportunite en fm | |
EP2309499B1 (fr) | Procédé de filtrage optimisé des bruits non stationnaires captés par un dispositif audio multi-microphone, notamment un dispositif téléphonique "mains libres" pour véhicule automobile | |
FR2944903A1 (fr) | Systeme et methode pour detecter des evenements audio anormaux | |
EP3440472B1 (fr) | Procede de detection de defauts non francs dans un cable par fusion de donnees | |
EP3555645B1 (fr) | Procede de calcul d'un reflectogramme pour l'analyse de defauts dans une ligne de transmission | |
WO2003048711A2 (fr) | System de detection de parole dans un signal audio en environnement bruite | |
EP1792278B1 (fr) | Procede de detection et de pistage de cibles ponctuelles, dans un systeme de surveillance optronique | |
EP2196821A1 (fr) | Procédé de filtrage cinématique temporel multidimensionel de plots radar, de tour d'antenne à tour d'antenne | |
WO2021023478A1 (fr) | Procede et systeme de surveillance d'un reseau de cables, par analyse en composantes principales | |
EP2452204B1 (fr) | Traitement de donnees multi-cibles pour radars passifs multi-statiques et multi-canaux | |
EP3627510A1 (fr) | Filtrage d'un signal sonore acquis par un systeme de reconnaissance vocale | |
EP3709049B1 (fr) | Systeme de traitement radar et procede de debruitage associe | |
EP3324361A1 (fr) | Procede de detection et de pistage de cibles | |
FR2913115A1 (fr) | Methode automatique de fusion de pistes et estimation de peignes d'harmoniques. | |
FR2955950A1 (fr) | Procede et dispositif permettant de determiner un parametre central associe a des coefficients d'identification obtenus par un algorithme de burg | |
EP0514252B1 (fr) | Procédé pour détecter des sources en mouvement et estimer leur trajectoire | |
FR3128564A1 (fr) | Procédé de suivi d’au moins un élément dans une succession temporelle de signaux physiques | |
FR2952190A1 (fr) | Procede de determination d'un parametre associe a la vibration d'un objet par vibrometrie laser | |
FR2853793A1 (fr) | Procede et dispositif de detection d'un changement de plan dans une sequence video, programme d'ordinateur et procede d'indexation correspondants | |
FR2733600A1 (fr) | Procede et dispositif de traitement de signal pour lever l'ambiguite vitesse d'un radar doppler |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |
Effective date: 20121228 |