FR 2715784 A1 19950804 FR 9401160 A 19940202 PROCEDE ET DISPOSITIF D'ANALYSE D'UN SIGNAL DE RETOUR ET ANNULEUR D'ECHO ADAPTATIF EN COMPORTANT APPLICATION La présente invention concerne un procédé et un dispositif pour détecter, dans un signal de retour, la présence de composantes autres que des composantes d'écho d'un signal direct. Elle trouve une application particulière mais non exclusive dans le domaine des annuleurs d'écho adaptatifs. Les annuleurs d'écho sont utilisés dans des liaisons bidirectionnelles. Leur but est d'annuler, ou du moins d'atténuer sensiblement, dans le signal circulant dans l'un des deux sens (le signal de retour) les composantes d'écho du signal circulant dans l'autre sens (le signal direct) apparaissant au niveau du site de réception du signal direct. Les échos peuvent être d'origine acoustique, par exemple dans l'application aux postes téléphoniques mains libres, ou d'origine électrique. L'annuleur d'écho comporte habituellement un filtre qui modélise les trajets d'écho au niveau du site de réception du signal direct. Ce filtre reçoit le signal direct, et sa sortie, qui est une estimation des composantes d'écho, est soustraite du signal de retour. Dans un annuleur d'écho adaptatif, les coefficients du filtre de modélisation sont adaptés en temps réel pour minimiser l'énergie du signal de retour résiduel. Lorsque le signal de retour contient des composantes utiles et non seulement des composantes d'écho du signal direct, il convient de figer la procédure d'adaptation des coefficients pour éviter de dégrader les composantes utiles et pour ne pas provoquer de divergences ou d'oscillations indésirables dans les valeurs des coefficients. Or, la détection dans le signal de retour de la présence de composantes utiles autres que les composantes d'écho du signal direct se révèle être un problème délicat. On peut certes envisager des calculs de corrélations entre les signaux directs et de retour sur la base desquels pourrait s'effectuer la détection. Mais une telle méthode implique des volumes de calculs considérables, actuellement incompatibles avec les applications pratiques. Un but de la présente invention est d'apporter une solution au problème ci-dessus, qui soit réaliste en termes de volume de calculs. L'invention propose ainsi un procédé pour détecter, dans un signal de retour, la présence de composantes autres que des composantes d'écho d'un signal direct, caractérisé en ce qu'on soumet le signal direct à un filtre de prédiction linéaire à réponse impulsionnelle finie pour déterminer un premier signal résiduel d'énergie minimale, on soumet le signal de retour à un filtre à réponse impulsionnelle finie dont les coefficients sont les mêmes que ceux du filtre de prédiction linéaire du signal direct pour déterminer un second signal résiduel, on calcule le rapport des énergies contenues dans les premier et second signaux résiduels, et on détermine si le signal de retour contient des composantes autres que les composantes d'écho du signal direct en comparant le rapport calculé à un seuil de détection. La détection repose sur une adaptation des techniques de prédiction linéaire habituellement utilisées pour résoudre des problèmes de compression de l'information. Lorsque le signal de retour ne contient essentiellement que des composantes d'écho du signal direct, le signal direct et le signal de retour sont spectralement voisins, de sorte que le rapport des énergies des signaux résiduels est sensiblement constant à une valeur dépendant de l'amplitude de la réponse des trajets d'écho. Lorsque le signal de retour contient en outre des composantes utiles, le rapport des énergies s'écarte sensiblement de la valeur précédente, dans la mesure o les coefficients communs des filtres, qui sont choisis pour minimiser l'énergie du premier signal résiduel, s'éloignent de ceux qui minimiseraient l'énergie du second signal résiduel. La comparaison du rapport des énergies calculé à un seuil de détection permet donc de discriminer entre les deux situations. Le recours à une technique de prédiction linéaire implique un volume de calcul très inférieur à la détermination directe de corrélations entre le signal direct et le signal de retour. L'ordre de la prédiction linéaire est typiquement compris entre 5 et 10. De préférence, la valeur du seuil de détection est adaptée de façon récursive en fonction des valeurs calculées du rapport d'énergie, l'adaptation du seuil de détection étant interrompue lorsqu'on détermine que le signal de retour contient des composantes autres que les composantes d'écho du signal direct et/ou lorsqu'on détecte la présence de composantes utiles dans le signal direct. Cette disposition permet au processus de détection de s'adapter aux variations éventuelles de la réponse des trajets d'écho. Un dispositif de détection pour la mise en oeuvre du procédé selon l'invention comprend un filtre de prédiction linéaire à réponse impulsionnelle finie recevant le signal direct et produisant un premier signal résiduel d'énergie minimale, un second filtre à réponse impulsionnelle finie ayant les mêmes coefficients que le filtre de prédiction linéaire du signal direct, recevant le signal de retour et produisant un second signal résiduel, des moyens de calcul du rapport des énergies contenues dans les premier et second signaux résiduels, et des moyens de comparaison du rapport d'énergies calculé à un seuil de détection pour déterminer si le signal de retour contient des composantes autres que des composantes d'écho du signal direct. Un autre aspect de la présente invention concerne un annuleur d'écho adaptatif, pour atténuer, dans un signal de retour, les composantes d'écho d'un signal direct, comprenant un filtre adaptatif de modélisation des trajets d'écho auxquels est adressé le signal direct et dont la sortie est soustraite du signal de retour, caractérisé en ce qu'il comprend en outre un dispositif de détection du type défini ci-dessus, pour détecter la présence dans le signal de retour de composantes autres que les composantes d'écho du signal direct, les coefficients du filtre adaptatif de modélisation étant figés lorsque ledit dispositif de détection révèle présence dans le signal de retour de composantes autres que les composantes d'écho du signal direct. D'autres particularités et avantages de la présente invention apparaîtront dans la description ci-après d'un exemple de réalisation non limitatif, lue conjointement aux dessins annexés, dans lesquels: - la figure 1 est un schéma synoptique d'un annuleur d'écho adaptatif incorporant un dispositif de détection selon l'invention; et - la figure 2 est un organigramme illustrant les étapes d'un procédé de détection selon l'invention. L'invention est décrite ci-après dans son application particulière, mais non exclusive, à un annuleur d'écho d'un poste téléphonique de type mains libres. Ce poste émet et reçoit des signaux acoustiques numérisés à une fréquence d'échantillonnage de 8 000 Hz par exemple. Le signal direct Sd(n) est adressé au haut-parleur 5 du poste par l'intermédiaire d'un convertisseur numérique-analogique 6 et d'un amplificateur 7. Le signal de sortie analogique du microphone 8 du poste est amplifié par l'amplificateur 9 puis numérisé par le convertisseur 10 pour former le signal numérique de retour Sr(n). L'annuleur d'écho comprend un filtre de modélisation 12 à réponse impulsionnelle finie qui reçoit les échantillons successifs du signal direct Sd(n), et produit en sortie une estimation A(n) de l'écho du signal direct. Chaque échantillon de l'estimation A(n) est soustrait de l'échantillon correspondant du signal de retour Sr(n) par le soustracteur 13 pour former le signal de retour résiduel Rr(n) dépourvu d'écho. Les coefficients du filtre 12 modélisent la réponse impulsionnelle des trajets d'écho entre le haut-parleur 5 et le microphone 8 (symbolisés par la flèche F sur la figure 1). Le filtre 12 comporte typiquement 2 000 coefficients. Les coefficients du filtre 12 sont mis à jour de manière adaptative pendant la communication, comme indiqué par le bloc fonctionnel 16 sur la figure 1. L'adaptation consiste à minimiser l'énergie contenue dans le signal de retour résiduel Rr(n). Elle peut se faire de façon classique par un algorithme de moindres carrés tel que l'algorithme du gradient. L'annuleur d'écho comporte un détecteur d'activité vocale 15 recevant le signal direct Sd(n). Le détecteur d'activité vocale 15 sert à détecter la présence de composantes utiles autres que du bruit dans le signal direct. Les détecteurs d'activité vocale sont bien connus dans le domaine du traitement des signaux acoustiques. On peut en trouver des exemples dans les articles "A voice activity detector based on cepstral analysis" de J. A. Haigh et al (Proc. Eurospeech 1993 - Berlin, Vol. 2, pages 1103- 1106), "Speech enhancement using a soft-decision noise suppresssion filter" de R.J. Mc Aulay et al (IEEE Trans. ASSP, Vol. 28, N' 2, pages 137-147, avril 1980), "Modification of piecewise LPC" de J. Roberts (MITRE Working paper WP-21752, Mai 1978), "Evaluation of linear and non- linear spectral subtraction methods for enhancing noisy speech" de A. Le Floc'h et al (ESCA, pages 131-134, Novembre 1992), "Voice activity detection using a periodicity measure" de R. Tucker (IEEE Proc., Vol. 139, N' 4, Août 1992), et "Speech/non-speech detection for voice response systems" de L. Mauuary et al (Proc. Eurospeech 1993 - Berlin, Vol. 2, pages 1097-1100). Le détecteur d'activité vocale 15 est utilisé pour commander les moyens d'adaptation 16, de façon que les coefficients du filtre 12 restent figés lorsqu'aucune composante utile n'est détectée dans le signal direct Sd(n). L'adaptation des coefficients n'a pas lieu non plus lorsqu'on détecte que le signal de retour Sr(n) contient des composantes autres que les composantes d'écho du signal direct. A cette fin, l'invention prévoit de munir l'annuleur d'écho d'un dispositif de détection 20 permettant de déterminer si le signal de retour contient ou non des composantes autres que les composantes d'écho du signal direct. Le dispositif de détection 20 comprend un filtre de prédiction linéaire 21 recevant le signal direct Sd(n). Les échantillons de sortie en du filtre 21 sont donnés par: p en = Sd (n) +_ aiSd(n-.i) i-i o p désigne l'ordre de la prédiction linéaire, et ai (0 A la fin de chaque bloc de longueur N. un diviseur 27 calcule le rapport d'énergies p = f/E. Un comparateur 28 compare la valeur calculée du rapport p à un seuil de détection Ps. Lorsque p > Ps, il est déterminé que le signal de retour Sr(n) contient des composantes autres que les composantes d'écho du signal direct Sd(n). Le comparateur 28 adresse alors aux moyens 16 d'adaptation des coefficients du filtre de modélisation 12 un signal X pour que les coefficients du filtre 12 restent figés lorsque le dispositif 20 a détecté la présence dans le signal de retour de composantes autres que des composantes d'écho du signal direct. Lorsque p Dans la version préférée de l'invention, la valeur du seuil de détection p. est adaptée de façon récursive en fonction des valeurs calculées du rapport p, comme schématisé par le bloc de calcul 29 sur la figure 1. Mais l'adaptation du seuil de détection Ps est interrompue lorsqu'on détermine que le signal de retour Sr(n) contient des composantes autres que les composantes d'écho du signal direct Sd(n). Il faut en effet éviter que la valeur du seuil p. augmente de manière inopportune lorsque le signal de retour contient des composantes utiles, ce qui dégraderait les performances de détection du dispositif 20. De préférence, l'adaptation du seuil de détection Ps est également interrompue lorsque le détecteur d'activité vocale révèle la présence de composantes utiles dans le signal direct S,(n). Dans ce cas, le seuil de détection p. n'est modifiable que lorsqu'il n'y a des composantes utiles ni dans le signal direct ni dans le signal de retour, c'est-à- dire pendant les périodes de silence. Le calcul du seuil ps comporte une estimation récursive de la valeur moyenne puy des rapports d'énergies p calculés pendant les périodes d'adaptation du seuil p, à l'aide d'une fenêtre d'oubli exponentielle: Pmoy = XPmoy+(1-X)p (1) o x désigne un coefficient d'oubli compris entre 0 et 1. La valeur maximum Pmax des rapports p calculés pendant les périodes d'adaptation du seuil ps est également déterminée, en comparant la valeur courante de p à la valeur p. mémorisée. La valeur du seuil ps est alors figée lorsqu'on détecte des composantes utiles dans l'un des deux signaux: Ps = Pmoy + (Pmx - Pmoy) (2) o c désigne un second coefficient compris entre 0 et 1. En prenant une valeur p. supérieure à Pmoy, on minimise les risques de fausse détection. Diverses expériences réalisées par les demandeurs ont révélé que les valeurs optimales des coefficients a et x sont comprises entre 0,8 et 0,9. Ces expériences ont montré d'excellentes performances de détection avec un modèle de prédiction linéaire d'ordre p = 10 appliqué à des blocs d'échantillons non recouvrants de N = 256 échantillons (soit 32 ms). La méthode a donné des résultats satisfaisants pour des blocs d'analyse variant de N = 64 à N = 256 échantillons avec des ordres de prédiction linéaire p allant de 5 à 10. Des ordres supérieurs à 10 alourdissent la charge de calcul sans apporter d'améliorations significatives. Matériellement, l'annuleur d'écho selon l'invention sera le plus souvent réalisé en programmant de façon appropriée un ou plusieurs processeurs de traitement de signal (DSP). La programmation des parties autres que le dispositif de détection 20 est classique et ne sera pas détaillée davantage. La partie correspondant au dispositif peut être réalisée conformément à l'organigramme représenté sur la figure 2. Au début du processus, les compteurs E, r et n sont initialisés à 0 (étape 40). Ensuite, à chaque réception d'un couple d'échantillons Sd(n), Sr(n), les signaux résiduels en et gn sont calculés (étape 41), puis les compteurs E, r sont incrémentés respectivement de en2 et de ga2 (étape 42) . Si l'indice n reste inférieur à N-1 (test 43), on l'incrémente d'une unité (étape 44) pour revenir à l'étape 41 dès réception du couple d'échantillons suivant. Quand n devient égal à N-l, le bloc courant est terminé, et on passe à l'étape 45 de minimisation de l'énergie résiduelle E et de détermination des coefficients a1, ....,. ap des filtres 21, 24. Le rapport p=F/E est alors calculé (étape 46), puis comparé, en 47, au seuil p.. Si p > ps, on décide que le signal de retour Sr(n) contient des composantes autres que des composantes d'écho du signal direct, et on fige les coefficients du filtre de modélisation 12 (étape 48), avant de revenir à l'étape d'initialisation 40 pour le traitement du bloc suivant. Si p Pm REVENDICATIONS 1. Procédé pour détecter, dans un signal de retour (Sr(n)), la présence de composantes autres que des composantes d'écho d'un signal direct (Sd(n)), caractérisé en ce qu'on soumet le signal direct à un filtre de prédiction linéaire à réponse impulsionnelle finie (21) pour déterminer un premier signal résiduel (en) d'énergie minimale, on soumet le signal de retour à un filtre à réponse impulsionnelle finie (24) dont les coefficients sont les mêmes que ceux du filtre de prédiction linéaire du signal direct pour déterminer un second signal résiduel (gn), on calcule le rapport (p) des énergies (T,E) contenues dans les premier et second signaux résiduels, et on détermine si le signal de retour contient des composantes autres que les composantes d'écho du signal direct en comparant le rapport calculé à un seuil de détection (Ps) . 2. Procédé selon la revendication 1, caractérisé en ce que la valeur du seuil de détection (ps) est adaptée de façon récursive en fonction des valeurs calculées dudit rapport d'énergies (p), l'adaptation du seuil de détection (p.) étant interrompue lorsqu'on détermine que le signal de retour (S.(n)) contient des composantes autres que les composantes d'écho du signal direct (Sd(n)). 3. Procédé selon la revendication 2, caractérisé en ce qu'on interrompt l'adaptation du seuil de détection (p.) lorsqu'on détecte la présence de composantes utiles dans le signal direct (S,(n)). 4. Procédé selon la revendication 2 ou 3, caractérisé en ce que, pendant les périodes d'adaptation du seuil de détection (p.), on calcule une moyenne Pmy des rapports d'énergies p successivement calculés par une formule récursive de la forme pmoy = X Pmoy + (1-X)p o X désigne un coefficient compris entre 0 et 1, on détermine une valeur maximale Pm 5. Procédé selon la revendication 4, caractérisé en ce que les coefficients X et a sont compris entre 0,8 et 0,9. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que l'ordre (p) du filtre de prédiction linéaire (21) est compris entre 5 et 10. 7. Dispositif de détection de la présence de composantes autres que des composantes d'écho d'un signal direct (Sd(n)) dans un signal de retour (Sr(n)), caractérisé en ce qu'il comprend un filtre de prédiction linéaire à réponse impulsionnelle finie (21) recevant le signal direct (S,(n)) et produisant un premier signal résiduel (en) d'énergie minimale, un second filtre à réponse impulsionnelle finie (24) ayant les mêmes coefficients que le filtre de prédiction linéaire du signal direct, recevant le signal de retour (Sr(n)) et produisant un second signal résiduel (gn), des moyens (22,25,27) de calcul du rapport (p) des énergies (r,E) contenues dans les premier et second signaux résiduels, et des moyens (28) de comparaison du rapport d'énergies calculé à un seuil de détection (p.) pour déterminer si le signal de retour contient de composantes autres que les composantes d'écho du signal direct. 8. Dispositif selon la revendication 7, caractérisé par des moyens (29) d'adaptation récursive de la valeur du seuil de détection (ps) en fonction des valeurs calculées du rapport d'énergie (p), dont le fonctionnement est interrompu lorsque les moyens de comparaison (28) révèlent que le signal de retour (Sr(n)) contient des composantes autres que les composantes d'écho du signal direct (Sd(n)). 9. Dispositif selon la revendication 8, caractérisé en ce que le fonctionnement des moyens d'adaptation (29) est interrompu lorsque des composantes utiles sont détectées dans le signal direct (Sd(n)). 10. Dispositif selon la revendication 8 ou 9, caractérisé en ce que les moyens d'adaptation (29) sont agencés pour déterminer la valeur du seuil de détection (ps) selon une formule du type ps = Pmy + a (Pmax - Pmoy) o Pmoy est une moyenne des rapports d'énergies p calculés pendant les périodes de fonctionnement des moyens d'adaptation, obtenue par une formule récursive de la forme Pmoy = x Pmoy + (1-) p Pmax est la valeur maximale des rapports d'énergies p calculés pendant les périodes de fonctionnement des moyens d'adaptation, et c et X sont deux coefficients compris entre 0 et 1. 11. Dispositif selon la revendication 10, caractérisé en ce que les coefficients a et A sont compris entre 0,8 et 0,9. 12. Dispositif selon l'une quelconque des revendications 7 à 11, caractérisé en ce que l'ordre (p) du filtre de prédiction linéaire (21) est compris entre 5 et 10. 13. Annuleur d'écho adaptatif, pour atténuer, dans un signal de retour (Sr(n)), les composantes d'écho d'un signal direct (Sd(n)), comprenant un filtre adaptatif (12) de modélisation des trajets d'écho auquel est adressé le signal direct et dont la sortie est soustraite du signal de retour, caractérisé en ce qu'il comprend en outre un dispositif de détection (20) selon l'une quelconque des revendication 7 à 12 pour détecter la présence dans le signal de retour (Sr(n)) de composantes autres que les composantes d'écho du signal direct (S,(n)), les coefficients du filtre adaptatif de modélisation (12) étant figés lorsque ledit dispositif de détection révèle la présence dans le signal de retour de composantes autres que les composantes d'écho du signal direct.