Brevet FR2704111A1 — publié le 19941021

PROCEDE DE DETECTION ENERGETIQUE DE SIGNAUX NOYES DANS DU BRUIT La présente invention se rapporte à un procédé de détection énergétique de signaux noyés dans du bruit. Lorsqu'on dispose d'un modèle du signal, des outils de détection de ce signal sont largement disponibles dans la littérature, les méthodes les plus connues étant basées sur la notion de filtre adapté et, plus généralement, sur la théorie de la décision en traitement du signal (P. Y. ARQUES, Collection Technique et Scientifique des Télécommunications, MASSON). Ces techniques président à l'élaboration des récepteurs cohérents et non cohérents en communications numériques (Principle of Coherent Communication A.J.VITERBI, MacGraw-Hill). Par contre, la présente invention se place dans le cas o l'on ne dispose pas de modèle susceptible de permettre l'application directe de la théorie de la détection. On suppose être en présence d'un bruit de fond, et de temps en temps se produit une "anomalie", qui, suivant le contexte, peut en fait représenter un signal qu'il serait souhaitable de détecter. Des exemples de détection d'un signal dit "utile" dans un bruit se trouvent largement dans la littérature concernant la détection de la parole. En effet, le signal de parole, par sa grande variabilité, ne se prête pas à une modélisation efficace et un des moyens les plus naturels pour le détecter consiste à effectuer un seuillage énergétique. Ainsi, beaucoup de recherches actuelles portent soit sur l'amplitude instantanée par référence à un seuil déterminé expérimentalement (La discrimination parole-bruit et ses applications V. PETIT, F. DUMONT Revue Technique THOMSON-CSF - Vol. 12 - N 4 - Dec. 1980), soit sur un seuillage empirique de l'énergie ("Suppression of Acoustic Noise in Speech Using Spectral Substraction", S.F. BOLL, IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol.ASSP-27, N .2, April 1979), soit sur l'énergie du signal total sur une tranche temporelle de durée T, en seuillant, toujours expérimentalement, cette énergie à l'aide d'histogrammes locaux, par exemple ("Problème de détection des frontières de mots en présence de bruits additifs",.P.WACRENIER, Mémoire de D.E.A. de l'université de PARIS- SUD, Centre d'ORSAY). D'autres techniques sont présentées dans "A Study of Endpoint Detection Algorithms in Adverse Conditions: Incidence on a DTW and HMM Recognizer", J.C. JUNQUA, B. REAVES, B. MAK EUROSPEECH 1991. Dans toutes ces approches, une grande part est faite à l'heuristique, et peu d'outils théoriques puissants sont utilisés. Il faut rappeler les travaux présentés dans "Evaluation of Linear and Non-Linear Spectral Substraction Methods for Enhancing Noisy Speech", A. LE FLOC'H, R. SALAMI, B. MOUY and J-P. ADOUL, Proceedings of "Speech Processing in Adverse Conditions", ESCA WORKSHOP, CANNES- MANDELIEU, 10-13 November 1992, o toutes les énergies supérieures à un certain seuil expérimental sont considérées comme révélatrices de la présence de signal utile et o toutes les énergies inférieures à ce seuil sont considérées comme des énergies dues au bruit seul lorsque la distance usuelle (valeur absolue de leur différence) les séparant est inférieure à un seuil, expérimental lui aussi. Cependant dans ce document de Le Floc'h et al, les auteurs travaillent sur la notion de distances entre énergies, mais la distance utilisée est une simple valeur absolue de la différence des énergies et leurs travaux font grandement appel à l'heuristique. La présente invention a pour objet un procédé de détection énergétique de signaux, utiles noyés dans du bruit, procédé qui mette essentiellement en oeuvre des moyens rigoureux sans pratiquement faire appel à l'heuristique, et qui soit optimisé, c'est-à-dire qui permette de détecter pratiquement tous les signaux utiles noyés dans du bruit, même intense, avec le plus faible taux possible de fausses détections. Le procédé conforme à l'invention consiste, à partir d'un ensemble d'échantillons d'un signal bruité groupés en trames successives, à effectuer une pré-classification en comparant les énergies des trames successives les unes par rapport aux autres au sens d'une distance qui est la valeur absolue de la différence des logarithmes des deux énergies, de manière à classer dans une première classe "bruit seul" les trames qui présentent une forte probabilité d'appartenir à cette classe, on détecte ensuite, pour les autres trames, celles présentant une énergie suffisamment élevée par rapport à une énergie de référence calculée à partir des énergies des trames "bruit seul" de manière que ces trames détectées présentent une forte probabilité d'appartenance à une deuxième classe "bruit + signal utile". Le procédé de l'invention suppose que lorsque le signal utile est présent, l'énergie du signal observé appartient à une certaine classe notée C1, et que lorsque le signal utile est absent, l'énergie observée appartient à une classe notée C2. Une des caractéristiques nouvelles de la présente invention est de pouvoir mettre en évidence de telles énergies de la classe C2 (donc des énergies de bruit seul), utilisées alors de manière à optimiser la détection d'énergies de la classe C1 (donc d'énergies révélant la présence de signal utile), par un procédé optimisé. On considère une distance entre énergies, mais la distance dont il est question dans l'invention, entre deux énergies U et V, n'est pas la distance usuelle IU-VI mais ILog(U/V)[, ce qui revient à considérer que deux énergies U et V sont proches l'une de l'autre lorsqu'on a 1/s U-V EN(1Ncs2,2N(0-s2 + 0_2)2 +2Nc7X4). Désignons par r = rs2 / cx,2 le rapport signal à bruit. On peut encore écrire: U-VEN(Nrco-x2,2Nco-4[(r+ 1)2 + 1]J). Le résultat dépend de ax2 et de r, ce qui montre qu'un seuillage de la distance [U-VI n'est pas valide lorsqu'on ne connalt pas ax2. Par contre, si on considère le rapport U/V, on démontre que la densité de probabilité de U/V ne dépend plus que de r, et est donc indépendante de ax2. Ce résultat remarquable valide l'utilisation d'un seuil sur U/V lorsqu'on ne connaît que r. De façon résumée, selon le procédé de l'invention, on observe L*N échantillons u(n) d'un signal. Chaque ensemble i = {u(N + k) /k {0,..., N - 1}}, lorsque i varie de 0 à L- 1, est appelé trame et est associé à une énergie E(Ti) notée Ui = E(Ti), ce qui permet de définir E = {Ui /i e {0,...,L - Il}}. Lorsque le signal utile est absent, les échantillons u(iN+k) sont exactement égaux aux échantillons du bruit noté x(iN+k) (u(iN+k) = x(iN+k)). Lorsque le signal utile (noté s(iN+k)) est présent, les échantillons u(iN+k) sont exactement égaux à u(iN+k) = s(iN+k) + x(iN+k). On met en évidence par un premier procédé, décrit ci-après (procédé dit de pré- classification) un sous-ensemble A d'éléments de E qui sont vraisemblablement des énergies de la classe C2. Il est alors possible de calculer un modèle autorégressif du bruit x(n) qui viendra blanchir les trames que l'on traitera par la suite, ou un spectre moyen du bruit x(n) qui peut servir à débruiter les trames suivantes (ni le blanchiment ni le débraitage ne sont impératifs mais sont utilisés suivant le contexte particulier traité). On utilise ensuite un second procédé (procédé dit de détection) décrit ci-dessous, qui détectera au mieux parmi les éléments de E (blanchis ou non, débruités ou non) les énergies de la classe C1. Soient alors N nouveaux échantillons, réunis sous forme d'une trame associée à une nouvelle énergie. On peut soit utiliser cette nouvelle énergie pour venir réactualiser l'ensemble A en utilisant le procédé de pré-classification, soit décider au sens d'un aspect particulier du procédé si cette nouvelle énergie appartient ou non à Ci, après un éventuel débruitage ou un éventuel blanchiment. Ce processus est réitéré pour chaque trame de N échantillons acquis. Le procédé de l'invention est caractérisé par l'utilisation d'outils théoriques nouveaux de traitement du signal et de statistiques. Ainsi, il fait appel à un modèle des lois statistiques que suivent les énergies des signaux, celui des Variables Aléatoires Gaussiennes Positives (VAGP) décrit ci-après. On utilise alors une propriété originale concernant le rapport de deux VAGP. On va maintenant définir les Variables Aléatoires Gaussiennes "Positives" (VAGP) utilisées par l'invention. Une variable aléatoire X sera dite positive lorsque Pr{ X Pr { X Dès que m/a est suffisamment grand, X peut être considérée comme positive. Lorsque X est gaussienne, on désigne par F(x) la fonction de répartition de la variable gaussienne normale et on a: Pr { X sorte que l'on peut encore écrire X e N(m,m2/a2). Modèles des énergies: exemples de variables gaussiennes "positives" Signal à énergie déterministe Soient les échantillons x(0),...x(N-l) d'un signal quelconque, dont l'énergie est déterministe et constante, ou approximable par une énergie déterministe ou constante (comme précisé ci-dessous). On a donc U = YO équiréparti entre [0,2Xr]. Pour N suffisamment grand, on a: (l/N) X0 n Pour N assez grand, U peut être assimilé à NA2/2 et donc à une énergie constante. On va maintenant examiner le cas de l'énergie d'un Processus Gaussien quelconque. Considérons un processus x(n), stationnaire du second ordre, mais gaussien, de variance crx2. On démontre le résultat suivant: U = yo Tr(Cx,N) = Ncx2. Donc U e N( Ncx2, 2Tr(Cx,N2)) Un calcul simple amène à Tr(Cx,N2) = 10 corrélation le permet. I1 existe des cas particuliers intéressants décrits ci- dessous, permettant d'accéder à cette fonction d'autocorrélation. Cas de l'énergie d'un Processus Blanc Gaussien. On considère le cas d'un processus blanc gaussien x(n) o n est compris entre 0 et N-1. Les échantillons sont indépendants et sont tous de même variance Cx2 = E[x(n)2]. On a alors Cx N = cyx2 IN, o IN est la matrice identité de dimension NxN. On en déduit: Tr(Cx,N2) = Ncx4 de sorte que: U = Io Le paramètre a est a = (N/2)1/2 Cas de l'énergie d'un Processus Gaussien Bande Etroite. On suppose que le signal numérique x(n) est issu de l'échantillonnage du processus x(t), lui-même issu du filtrage d'un bruit blanc gaussien b(t) par un filtre passe-bande h(t) de fonction de transfert: H(f) = U[-f0.B/2,-f0+B/2](f) + U[f0-B/2,f0+B/2](f), o U désigne la fonction caractéristique de l'intervalle en indice et f0 la fréquence centrale du filtre. La fonction de corrélation Fx(r) de x(t) vaut Fx(t) =Fx(0)cos(2if0x)sinc(nrBr) o sinc(x) = sin(x)/x. La fonction de corrélation de x(n) est alors: Fx(k) =]Fx(O)cos(2xkfoTe).sinc(7rkBTe). Si gfoBTe(k) = cos(2irkf0Te)sinc(nkBTe), on a: Tr(CxN2) = Fx(0)20o On a: U E N( Nax2, 2cx410 Ces relations restent valables même si f0 = 0. Cas de l'énergie d'un processus gaussien quelconque "sous- échantillonné". Ce modèle est plus pratique que théorique. Si la fonction de corrélation est inconnue, on sait cependant que: lim k --> + oo Fx(k) = 0. Donc, pour k assez grand tel que k > k0, la fonction de corrélation tend vers 0. Aussi, au lieu de traiter la suite d'échantillons x(0)...x(N-l), peut-on traiter la sous-suite x(0), x(k0), x(2k0),..., et l'énergie associée à cette suite reste une variable aléatoire positive gaussienne, à condition qu'il reste dans cette sous- suite suffisamment de points pour pouvoir appliquer les approximations dues au théorème central-limite. Cette façon de procéder peut permettre dans certains cas difficiles d'appliquer les règles de décision qui sont décrites ci-dessous. Résultat théorique fondamental. Si X = X1/X2 o X1 et X2 sont toutes deux gaussiennes, indépendantes, telles que: X1 N( ml; 122) et X2 e N(m2; c22). On pose m = ml/m2, al = ml/cl, c2 = m2/c2. - Lorsque cal et a2 sont suffisamment grands pour pouvoir supposer que X1 et X2 sont des variables aléatoires gaussiennes positives, la densité de probabilité fX(x) de X = X1/X2 peut alors être approchée par: a12 ca22 (x-m)2 2(a 2x2 + a 2m2) fx (x) = (2)-1/2 ala2mal2 x + a22m e U(m) (al2x2 +a22m2)3/2 o U(x) est la fonction indicatrice de R+: U(x) = 1 six x-m 2 (xmaa)F[(yIp) Sih(x,ml al,a2)= al a2 +2m2 p(x, mlal, a2)= F[h(x, y2a,))] (al 2X + a22m)1/ o F désigne la fonction de répartition de la variable normale, et o P(x, m al, a2) = Pr{X de ces paramètres fixées par une connaissance a priori ou par l'heuristique. On va maintenant exposer l'étape de pré-classification du procédé de l'invention. On suppose que C1 =N(m1, c12) représente les énergies observables en présence de signal utile et que C2 = N(m2, c22) représente les énergies observables en absence de signal utile. On pose m = ml/m2, al = ml/crl et a2 = m2/c2 et on suppose al1 et a que les éléments de C1 et de C2 soient des VAGP. E={U{,...,U,} est l'ensemble des énergies dont on dispose. Chacune de ces énergies Ui vaut Ui = 0 des énergies de la classe C2. Cette étape fait appel aux notions présentées ci- dessous. Notion de compatibilité entre énergies: Soit (U, V) e (C1UC2)X(C1UC2) et X = U/V. On définit les hypothèses suivantes: H1: (U,V) e (C1XC1)U(C2UC2) et H2: (U,V) e (C1XC2)U(C2UC1). Si on a: 1/s à-dire que H1 est considérée conmme vraie. On dira que U et V sont compatibles. Cette décision sera notée D1. Mais si on a: X s > on décide que U et V n'appartiennent pas à la même classe, c'est-à- dire que H2 est considérée comme vraie. On dira que U et V sont incompatibles. Cette décision sera notée D2. Si I = [1/s,s], la règle s'exprime selon: x e I On cherche- à optimiser cette règle de décision qui permnnet d'associer les réalisations de variables aléatoires entre elles. Pour ce faire, on calcule le seuil s optimal. Ce calcul est différent suivant que l'on connaît ou non la probabilité p. Lorsque p est connue, on applique directement le critère du maximum de vraisemblance. Lorsque p est inconnue, et comme les hypothèses sont réduites au nombre de 2, on utilise le critère de Neyman-Pearson. Critère du Maximum de vraisemblance: On démontre que la probabilité de décision correcte est: PC = p2[2P(s, lial, al) - 1]+(1-p)2[2P(s, lia2, a2) - 1] +2p(1- p)[2- P(s, I/mlal, a2)- P(s,mlal, a2)] Le seuil optimal s vérifie RPc = O. ôs ôP9c O p2f(s ll, al) +(1 _ p)2 f(s, l[a2, a2) = ôs p(l -p)[f (s, 1/mial, a2) + f(s,mal, a2)] Cette équation se résout sur calculateur, une fois fixées les valeurs m, p, al et a2. Critère de Neyman-Pearson: Lorsque p est inconnu, on utilise une approche du type Neyman- Pearson. On dira qu'il y a détection si on prend la décision D1, c'est-àdire si on décide que les deux variables aléatoires sont de la même classe. On définit alors les probabilités de non détection Pnd et de fausse alanne Pfa par: Pnd-Pr{D2[Hi} (probabilité de décider de l'incompatibilité, alors que les variables sont de la même classe) et Pfa= Pr {D1 [ H2} (probabilité de décider de la compatibilité alors que les variables sont incompatibles). Le critère de Neyman-Pearson consiste à minimiser Pnd lorsque Pfa est fixée (ou inversement). Ce type de critère s'applique lorsqu'une erreur est beaucoup plus grave que l'autre. Comme il est question ici de savoir si les variables aléatoires observées appartienent ou non à une même classe, il est évident que l'on cherchera à n'avoir, dans les réalisations retenues comme étant des réalisations de variables d'une même classe, que peu d'erreurs. C'est donc la Pfa que l'on va fixer de manière à n'avoir que très peu de fausses alarmes. Pfa = 1+P(l/ s,mIal, a2) - P(s,mlal, a2) et _ 1_p2[2P(sliai, al) - 1] +(i-P)2 [2P(s, Ila2, a2) - Pnd+1- 112 p p)2 de sorte que lorsque a1 # a2, Pnd dépend de p, qui est inconnu et est inaccessible. Dans le cas o acl = a2 = ac, alors Pnd = 2.P (s,lia,a) - 1 et est donc accessible. On peut fixer Pnd dans ce cas. Ayant l'expression de Pfa (ou de Pnd), on fixe cette probabilité, cequi permet d'obtenir leseuil s correspondant. Compatibilités entre plusieurs énergies. Lorsque le seuil a été calculé suivant un des deux procédés précités, il est intéressant de généraliser cette notion de compatibilité entre plusieurs énergies. Soient alors U1,..., UN, N énergies, on dira que ces énergies sont compatibles entre elles, si et seulement si, V i et j, Ui et Uj sont compatibles au sens évoqué ci-dessus, autrement dit, si toutes ces énergies sont compatibles deux à deux. Pour la mise en oeuvre du procédé, on fait les hypothèses suivantes: - les énergies de la classe C2 sont plus faibles statistiquement que celles de la classe C1; - la trame présentant l'énergie la plus faible est une trame de la classe C2. Soit Tio, cette trame. Le calcul se déroule alors comme suit: L'ensemble A est initialisé: A = {Tio}; POUR i décrivant {E(T1),..., E(Tn)} - {E(Ti0)} FAIRE Si E(Ti) est compatible avec chaque élément de A: A = A U {E(Ti)}. FIN POUR Le procédé de confirmation de bruit fournit un certain nombre de trames qui peuvent être considérées comme du bruit avec une très forte probabilité. On calcule, à partir de la donnée des échantillons temporels, un modèle autorégressif du bruit. Si x(n) désigne les échantillons de bruit, on modélise x(n) selon: x(n)= l_ méthodes sont disponibles (Burg, Levinson-Durbin, Kalman, Fast Kalman... ). De façon avantageuse, on met en oeuvre des procédés du type Kalmnan et Fast Kalman: "Le Filtrage Adaptatif Transverse", O. MACCHI, M. BELLANGER, Traitement du signal, Vol. 5, N 3, 1988 et "Analyse des signaux et filtrage numérique adaptatif', M. BELLANGER, Collection CNET-ENST, MASSON, qui présentent de très bonnes performances temps réel. Disposant d'un modèle autorégressif du bruit, il est alors aisé de blanchir ce bruit, ce qui permet de travailler sur un bruit blanc gaussien aisément manipulable. Soit u(n) = s(n) + x(n) le signal total, composé du signal utile s(n) et du bruit x(n). Soit le filtre H(z) = 1- -l vient H(z)U(z) = H(z)S(z) + H(z)X(z). Or H(z)X(z) = E(z) > H(z)U(z) = H(z)S(z) + E(z). Le filtre réjecteur H(z) blanchit le signal, de sorte que le signal en sortie de ce filtre est un signal utile (filtré donc déformé), additionné d'un bruit généralement blanc et gaussien. Travailler sur un bruit blanchi permet de se rapprocher d'hypothèses idéales, notamment lors de l'application du procédé de détection. Cependant ce blanchiment n'est pas obligatoire et on peut appliquer directement le procédé de détection sans passer par cet intermédiaire. Etant donné qu'après mise en oeuvre du procédé de l'invention, on dispose d'un certain nombre de trames confirmées comme étant des trames de bruit, on peut aussi calculer un spectre moyen de ce bruit, de manière à l1 implanter un filtrage spectral du type soustraction spectrale ou filtrage de WIENER, type largement décrit dans la littérature: "Suppression of Acoustic Noise in Speech Using Spectral Substraction" S. F. BOLL, IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-27, N 2, April 1979; "Enhancement and Bandwidth Compression of Noisy Speech", J. S. LIM, A. V. OPPENHEIM, Proceedings of the IEEE, Vol. 67, N 12, Dec. 1979, et "Noise Reduction For Speech Enhancement In Cars: Non-Linear Spectral Subtraction, Kalman Filtering", P. LOCKWOOD, C. BAILLARGEAT, J. M. GILLOT, J. BOUDY, G. FAUCON, EUROSPEECH 91. Cet aspect peut être intéressant dans ces certaines applications, voir par exemple: "Procédé de détection de la parole", D. PASTOR, demande de brevet français N 92 12582, déposée le 21.10.92. Détection selon le procédé de l'invention. Disposant d'un ensemble A dont les composantes sont vraisemblablement des énergies de la classe C2 (après un éventuel blanchiment), on cherche à détecter grâce à ces références les énergies de la classe C1. Si V représente la valeur moyenne des énergies de l'ensemble A, cette variable est aussi une VAGP. Si A = {V1,..., VM}, on a, avec les notations déjà utilisées, V i e { 1,..., M}, Vi e N(m2, cv22). Eo =(1/M)l On pose m = ml/m2, cl = ml/a et 02 = m2/c2. On passe alors'à la règle de décision optimale. Application du critère du maximum de vraisemblance (on connaît la probabilité p de décision correcte): soit p = Pr {U e C1}. La règle de décision optimale est alors: pf(x,m[cul, MI/2cz2) >(1-p)f(x, la22, MI/2a2) D = D2 Application du critère de Neyman-Pearson: Lorsqu'on ne connaît pas la valeur de p. on peut: - soit la fixer arbitrairement par une approche heuristique, - soit la fixer à p = 0,5, ce qui place dans le cas le plus défavorable, - soit utiliser le critère de Neyman-Pearson ou le critère de la médiane qui consiste à avoir: probabilité de fausse alarme = probabilité de non détection. Si l'on veut appliquer le critère Neyman-Pearson ou le critère de la médiane, la règle de détection sera de la forme: f(xm|al, M112 a2) / f (x, l1ax2, M1/2a2) >A tD = D, f(x,mlxa,M122a2)/f(x, ila2, M12a2)>Z <>D = D2 Le seuil X est fixé pour atteindre une valeur a priori de la probabilité de fausse alarme (ou la probabilité de décision correcte). Cette probabilité PFA de fausse alarme vaut: [f(X,mla",M/2 a2)] PFA =Pr1 f(XMiaî,M1/2a2) >À X C2 f Il n'y a a priori pas de calcul théorique simple de cette expression, donc pas de moyen théorique d'évaluer le seuil X. Le calcul de X peut cependant se faire sous simulation, suivant le cas particulier traité. La règle de décision simplifiée exposée ci-dessous est plus pratique à utiliser dans ce cas. Règle de décision simplifiée. Cette règle est: x>s<> U C1, x>s≤>U eC2 Cas du critère du maximum de vraisemblance: la probabilité Pc de décision correcte est: Pc= pl- P(s, m|alMîI2a2)] +(1- p)P(s, lia2, MI/2a2) Le seuil optimal est obtenu pour: C = pf(sm lalM2a2)-(î-p)f(sla2,M a2)=0 ôs Cas du critère de Neyman-Pearson: Lorsque la probabilité p est inconnue, on peut: - soit la fixer arbitrairement par une approche heuristique, - soit la fixer à p = 0,5, ce qui place dans le cas le plus défavorable, - soit utiliser le critère de Neyman-Pearson ou le critère de la médiane qui consiste à avoir Probabilité de fausse alarme = Probabilité de nondétection. Pour appliquer le critère de Neyman-Pearson ou le critère de la médiane on défminit les probabilités de non-détection et de fausse alarme: Pnd = {x sH2} On a: Pnd= P(s, lca2, Ml/2 a2)etPfa= 1- P(s, mil a,,M 1/2a2) On fixe alors Pfa ou Pnd, pour déterminer la valeur du seuil. Le critère de la médiane amène à: Pfa = Pnd > P(s,1|a2, M1 a2) = 1- P(s, mial, M1I2a2) Mise en oeuvre. Lorsque la règle de décision a été définie à l'aide des outils théoriques rappelés ci-dessus, et disposant d'une énergie E0 de "référence" de bruit, on effectue la détection sur E(T1),..., E(Tn), avec: E(Ti) =1O o les ui (n) sont les N échantillons constituant la trame Ti. Parmi les trames dont on dispose au départ, l'algorithme de préclassification a mis en évidence un ensemble A de trames presque sûrement de la classe "bruit". La moyenne des énergies des trames de l'ensemble A permet d'obtenir une valeur de référence Eo qui va être utilisée par l'algorithme de détection de manière à classifier les énergies des trames autres que celles de l'ensemble A ainsi que les nouvelles trames acquises ultérieurement. POUR E(Ti) décrivant {E(T1),..., E(Tn)} FAIRE X = E(Ti)/E0 Cas de la détection optimale: SI pf(X, mla, Ml/2a2) >(1- p)f(X, 1la2, M/2 a2), détection sur la trame Ti. Cas d'une détection à seuil: SIX > s détection sur la trame Ti FIN POUR POUR chaque nouvelle trame T représentée par les échantillons u(0),,u (N-l), FAIRE E = no Cas de la détection optimale. SI pf(X,mlal,Ml/2a2)>(1-p)f(X,1la2,Ml/2a2), détection sur la trame Ti Cas d'une détection à seuil. SIX > s Détection sur la trame Ti (s'il n'y a pas de détection, la trame acquise peut être considérée comme du bruit et éventuellement venir réactualiser A ainsi que la valeur Eo de référence). FIN POUR Exemples d'application. Il est possible de traiter un grand nombre d'exemples permettant de mettre en évidence l'intérêt du procédé de l'invention. En fait, il existe autant d'exemples que de couples demodèles que l'on peut former à partir de ceux décrits ci-dessus (voir les exemples de VAGP donnés ci-dessus): - détection d'un bruit blanc gaussien dans un autre bruit blanc gaussien; - détection d'un bruit blanc gaussien dans un bruit gaussien bande étroite; - détection d'une énergie déterministe dans un bruit gaussien bande étroite... Détection d'un signal d'énergie bornée dans un bruit gaussien bande étroite: Hypothèse 1: nous supposons ne pas connaître le signal utile dans sa forme, mais nous ferons l'hypothèse suivante: pour toute réalisation s(0),..., s(N-1) de s(n), l'énergie S définie par: S = (1/N)YO NPs2À Hypothèse 2: Le signal utile est perturbé par un bruit additif noté x(n), que l'on suppose gaussien et bande étroite. On suppose que le processus x(n) traité est obtenu par filtrage bande étroite d'un bruit blanc gaussien. La fonction de corrélation d'un tel processus est alors: Fx(k) = Fx(0)cos(27rkfoTe)sinc(rkBTe) Si on considère N échantillons x(n) de ce bruit, on a alors V = (1/N)o_ C- =0 - (l 1/2,1'/21 (Z0_ E[s(n)x(n)]/(E[s(n)2]E[x(n)2])1/2 lorsque les processus sont ergodiques. Soit u(n) = s(n) + x(n) le signal total, et U = 0_ On approxime U par: U = 0o ps2 on aura: U > Npus2 + 0_ Hypothèse 4: Comme nous supposons que le signal présente une énergie moyenne bornée, nous supposerons qu'un procédé capable de détecter une énergie ps2, sera capable de détecter tout signal d'énergie supérieure. Compte tenu des hypothèses précédentes, on définit la classe C1 comme étant la classe des énergies lorsque le signal utile est présent. Selon l'hypothèse 3, U > Nps2 + 20 2ax40o Donc C1 = N( NUPs2 + Ncx2, 2ax40o le rapport signal à bruit. C2 est la classe des; énergies correspondant au bruit seul. D'après l'hypothèse 2, si les échantillons de bruit sont x(0),...,x(M-1), il vient: V = (1/M)o0 cx2 = m2/a2 = M/[2Eo On peut alors mettre en oeuvre les étapes du procédé de l'invention exposées ci-dessus. Détection de code PN On considère une modulation BPSK étalée par un code PN de longueur L très grande devant 1. La durée d'émission d'un élément binaire dn est Tb. La durée d'émission d'un élément binaire du code PN est T'. Sur l'intervalle de [nTb,(n+1)Tb], le signal émis est: m(t) = (2Eb/Tb)1/2 dn o0 - K désigne le nombre d'échantillons du code PN vus sur cet intervalle, et. - (p la phase aléatoire équirépartie sur [0,27r] Ce signal émis est noyé dans un bruit de fond qui est b(t), que l'on suppose blanc et gaussien. On cherche alors à détecter le signal s(t) à partir du signal reçu r(t) = m(t) + b(t), en supposant que l'on ne connaît pas le code PN, donc ni les valeurs ck ni la durée L, ni le temps Tb, ni la fréquence coo. Soit alors la variable aléatoire: u(n) = (2/T)1/2 (nT)T r(t)cos(cnt)dt, o: - T est une durée d'intégration assez grande pour que les échantillons du code PN vus sur cet intervalle soient suffisamment nombreux et décorrélés, tout en restant assez faibles pour qu'on n'atteigne pas la prériodicité L du code PN. Si K est le nombre d'éléments binaires du code PN vu sur cet intervalle, on suppose ddnc: L >> 1, K >1. T vérifie en outre ooT >> 1 - co est une fréquence qui sert à essayer de récupérer la porteuse, telle que oT" >> 1 On pose: s(n) =(2/T)/2 f (nl (t)coTs(ot)d et: x(n) =(21T)121 ( n+1).T)cot)dt u(n) = s(n) + x(n), T, (k + rp) co 0 t)dt s(n) = (2Eb/Tb)112 dn YO variance: as2 = (Tb/T)(Eb/2K)sinc2(n(co - o)/K). En pratique, on suppose que les s(n) sont indépendants, de sorte que la suite des échantillons s(n) constitue un processus discret blanc gaussien. De même, la suite des échantillons x(n) constitue un processus blanc gaussien de moyenne nulle et de variance as2 = Ob2. Détecter le code PN, c'est détecter s(n), donc détecter un bruit blanc gaussien noyé dans un autre bruit blanc gaussien. Soit alors la variable U = TO U e N(N(as2+ ax2); 2N(cs2+ cyx2)2). Le paramètre aE de cette variable est al = (N/2)112 Soit alors la variable V = 70 V e N(Mcrx2; 2Max4). Le paramètre aE de cette variable est o2 = (M/2)112. On peut donc modéliser comme dans le cas d'une détection entre les deux classes: C1 = N(N(as2+ ax2); 2N((as2+ ax2)2) et C2 = N(Mcax2; 2Max4). On a alors: m = (N/M)(1+r), cal = (N/2)112, OE2 = (M/2)112 On notera que si N = M => al = a2 = (N/2)112 Le procédé décrit ci-dessus est donc applicable à ce problème. REVENDICATIONS 1. Procédé de détection énergétique de signaux utiles noyés dans du bruit, caractérisé par le fait qu'il consiste, à partir d'un ensemble d'échantillons d'un signal bruité groupés en trames successives, à effectuer une pré- classification en comparant les énergies des trames successives les unes par rapport aux autres au sens d'une distance qui est la valeur absolue de la différence des logarithmes des deux énergies, de manière à classer dans une première classe "bruit seul" les trames qui présentent une forte probabilité d'appartenir à cette classe, on détecte ensuite, pour les autres trames, celles présentant une énergie suffisamment élevée par rapport à une énergie de référence calculée à partir des énergies des trames "bruit seul" de manière que ces trames détectées présentent une forte probabilité d'appartenance à une deuxième classe "bruit + signal utile". 2. Procédé selon la revendication 1, caractérisé par le fait que le calcul d'un seuil optimal fait par application du critère du maximum de vraisemblance lorsque l'on connaît la probabilité de décision correcte. 3. Procédé selon la revendication 1, caractérisé par le fait que le calcul du seuil optimal est fait par application du critère de Neyman- Pearson lorsque l'on ne connaît pas la probabilité de décision correcte. 4. Procédé selon l'une des revendications précédentes, caractérisé par le fait que l'appartenance à la deuxième classe est déterminée par application du critère du maximum de vraisemblance lorsque l'on connaît la probabilité de décision correcte. 5. Procédé selon l'une des revendications 1 à 3, caractérisé par le fait que l'appartenance à la deuxième classe est déterminée par application du critère de Neyman-Pearson lorsque l'on ne connaît pas la probabilité de décision correcte. 6. Procédé selon l'une des revendications précédentes, caractérisé par le fait qu'après avoir traité une première trame d'échantillons, on calcule un modèle autorégressif du bruit que l'on utilise pour blanchir les trames suivantes, à partir d'un ensemble de trames constituées, avec une forte probabilité, uniquement de bruit. 7. Procédé selon l'une des revendications 1 à 5, caractérisé par le fait qu'après avoir mis en évidence un ensemble de trames constituées, avec une forte probabilité, de bruit seul, on calcule un spectre moyen du bruit et qu'on applique un débruitage aux trames suivantes. 8. Procédé selon l'une des revendications précédentes, caractérisé par le fait que qu'on se sert des résultats acquis pour une trame d'échantillons pour réactualiser les résultats de la trame précédente.