Brevet FR2729804A1 — publié le 19960726

ANNULEUR D'éCHO ACOUSTIQUE A FILTRE ADAPTATIF ET PASSADE DANS LE DOMAINE FREQUENTIEL La présente invention concerne les annuleurs d'écho acoustique destinés à des terminaux téléphoniques utilisés dans des conditions telles que le signal de sortie vers une ligne téléphonique est affecté d'un écho acoustique et souvent de bruits dus à l'environnement. Elle trouve une application particulièrement importante dans les installations susceptibles de comporter des terminaux mains-libres portables et dans les installations de télé-conférence. On utilise souvent, dans les annuleurs d'écho acoustique, des filtres adaptatifs transverses travaillant dans le domaine temporel. Leur réalisation se heurte à de nombreuses difficultés. Le nombre de coefficients nécessaire pour modéliser la réponse impulsionnelle d'un trajet d'écho acoustique est très élevé. Les signaux de parole ont un degré élevé d'auto-corrélation et ils sont non-stationnaires. Enfin le signal reçu par le capteur d'entrée (microphone en général) comporte non seulement l'écho acoustique, mais aussi le bruit ambiant et la parole du locuteur proche, en cas de parole simultanée. Pour tenir compte de la longueur de la réponse impulsionnelle et du degré élevé de corrélation, on a déjà proposé d'utiliser des filtres mettant en oeuvre un algorithme travaillant dans le domaine fréquentiel. Mais des problèmes de résistance au bruit ambiant et de prise en compte de la situation de paroles simultanées du locuteur proche et du locuteur lointain subsistent. Des tentatives ont été faites pour améliorer les caractéristiques des filtres transverses adaptatifs travaillant dans le domaine temporel ou dans le domaine fréquentiel. Dans le domaine temporel, on a en particulier proposé d'appliquer l'algorithme du gradient stochastique, dit LMS, par blocs d'échantillons plutôt qu'échantillon par échantillon. Le gradient est alors estimé bloc à bloc et fournit une indica tion de l'inter-corrélation entre les signaux d'entrée et d'erreur sur le bloc correspondant. On a également proposé des annuleurs d'écho adaptatifs opérant dans le domaine fréquentiel basés sur des techniques de convolution rapide de type dit OLS (recouvrement et conservation) ou OLA (recouvrement et addition). La figure 1 montre le schéma général d'un filtre adaptatif utilisable dans un annuleur d'écho de ce genre.Le signal numérique entrant x(n), où n est l'ordre de l'échantillon est appliqué à un haut parleur 12 après conversion numérique-analogique CNÂ et amplification par un amplificateur 10. L'écho acoustique arrive, par un trajet 13, au détecteur 14, constitué généralement par un microphone. Ce microphone reçoit également le bruit ambiant, et, à certains moments, un signal de parole du locuteur proche. Le signal, amplifié en 16, est mis sous forme numérique z(n) par un convertisseur analogique-numérique 22. L'annuleur d'écho du terminal proche est interposé entre la ligne de réception LR recevant le signal entrant et la ligne d'émission LE fournissant le signal sortant. Du fait que le filtre opère sur des blocs, il comporte, sur la ligne LE, un convertisseur d'entrée série-parallèle 23 et un convertisseur de sortie parallèle-série 24 entre lesquels est interposé un additionneur 26 qui reçoit sur une entrée additive les blocs d'échantillons d'entrée, et sur une entrée soustractive les blocs de signaux de correction, élaborés par un filtre adaptatif 28. L'annuleur d'écho montré en figure 1 travaille dans le domaine fréquentiel et utilise une procédure de partition avec recouvrement. Cela conduit à calculer une transformée de Fourier (ou plus généralement une transformée de passage du domaine temporel au domaine fréquentiel), ayant une taille égale à deux fois la longueur de la réponse impulsionnelle estimée N. Le filtre 28 reçoit le signal entrant x(n) par l'intermédiaire d'un convertisseur série-parallèle 32 permettant de constituer des blocs, d'un circuit de partition en blocs, avec recouvrement, et d'un circuit 36 de calcul de transformée de Fourier discrète, dite TFD ou DFT. Symétriquement, la sortie du filtre 28 est appliquée à un circuit 38 de calcul de transformée de Fourier inverse, ou DFT-I et, par l'intermédiaire d'un circuit de recouvrement 40, récupérant le dernier bloc à l'additionneur 26. Les moyens d'adaptation des coefficients du filtre fréquentiel 28 ont une constitution générale connue. Ils comportent un multiplieur 42 qui reçoit dune part les échantillons de sortie de l'additionneur 26 par l'intermédiaire de moyens de partition et d'un circuit de calcul de TFD 46 et d'autre part la conjuguée 48 de la TFD élaborée par le circuit 36. La sortie du multiplieur 42 attaque le circuit 50 d'adaptation des coefficients du filtre 28 ; ce circuit peut utiliser un des algorithmes habituels, tels que lalgo- rithme fréquentiel du gradient stochastique, dit FLMS. Dans le mode de réalisation montré en figure 1, une contrainte supplémentaire sur le gradient est apportée par un circuit 52. Ce circuit comporte une TFD inverse 54 et une TFD directe 56, nécessaires de façon à forcer à zéro les derniers N échantillons. de la réponse impulsionnelle, les transformées étant calculées sur 2N, c c'est-à-dire deux fois la longueur de la réponse impulsionnelle estimée. La contrainte sur le gradient peut être mémorisée sous forme d'une matrice de contrainte dans un circuit 58 interposé entre les transformées inverse et directe. Au lieu d'une transformée de Fourier, on peut utiliser une autre transformée dans le domaine fréquentiel, par exemple la transformée en cosinus ou la transformée de Hartley rapide. Diverses réalisations de filtres acoustiques ayant la constitution généralement montrée en figure 1 ont déjà été décrites. On peut par exemple se reporter à l'arti- cle de Borrallo et al, "On the implementation of a partitioned block frequency domain adaptive filtre (PBFDAF) for long acoustic echo cancellation, Signal processing, volume 27, No. 3, pages 301-315, juin 1992. Lorsque la réponse impulsionnelle est courte, ce qui est notamment le cas d'un téléphone mains-libres destiné à être utilisé dans un véhicule, une technique à recouvrement de blocs peut être utilisée sans qu'une partition des blocs soit nécessaire pour réduire la complexité. Elle met par exemple en oeuvre 1algorithme dit FLMSO (contraint) ou UFLMSO (non contraint). Il est difficile de donner aux annuleurs d'écho fréquentiels du genre décrit ci-dessus une résistance suffisante au bruit. Les algorithmes habituels sont par ailleurs très sensibles à la valeur du pas d'incrémentation. Aux faibles rapports signal un bruit, il est nécessaire de réduire cette valeur, ce qui ralentit l'adaptation ; même dans ce cas, des débuts de divergence peuvent s'amorcer. La présente invention vise notamment à fournir un annuleur d'écho fréquentiel amélioré de ce point de vue. Dans ce but l'invention propose d'utiliser une technique de soustraction spectrale, dont l'emploi a déjà été envisagé dans d'autres domaines tels que la reconnaissance de paroles, en profitant de ce que le passage dans le domaine fréquentiel est déjà disponible dans les annuleurs déchois du genre cidessus définis. En conséquence, l'invention propose un annuleur d'écho acoustique disposé à être placé entre une interface acoustique mains libres et un réseau de communication, comprenant - un filtre adaptatif recevant les coefficients d'une transformée directe discrète dans le domaine fréquentiel, calculés sur des blocs successifs d'un signal d'entrée, et attaquant, par l'intermédiaire d'un circuit de transformée inverse, l'entrée soustractive d'un additionneur de formation d'un signal de sortie allant vers le réseau ; et - un circuit d'adaptation des coefficients du filtre, caractérisé en ce que ledit circuit d'adaptation reçoit les coefficients des transformées successives des blocs d'échantillons du signal de sortie par l'intermédiaire d'un réducteur de bruit.Ce réducteur de bruit est avantageusement un soustracteur spectral dont l'entrée soustractive reçoit les coefficients représentatifs du spectre estimé du bruit. On peut utiliser divers modèles de bruit, suivant l'environnement. Le modèle peut être estimé a priori ou plus généralement déterminé au cours de périodes de silence (pause de parole de la part des locuteurs proche et éloigné). La constitution des moyens de soustraction spectrale de bruit et les algorithmes mis en oeuvre peuvent être ceux décrits dans l'article de Lockwood et Boudy J. "Experiments with a Nonlinear Specral Subtractor (NSS). Hidden Markov Models and the projection, for robust speech recognition in cars-, Speech Communication, Vol. 11, Nos 2-3, pp 215-228, June 1992, ou dans le document EP-A-0 534 837. L'algorithme d'adaptation du poids des coefficients du filtre travaille dans le domaine fréquentiel peut alors être o (k+l) = H (k) + 2F.g.F-1.m(k).XH(k).SPS [E (k)] où k est lordre du bloc courant, F désigne la matrice de transformée de Fourier discrète, ou DFT g est la matrice de contrainte sur le gradient m est la matrice diagonale des pas dincrémentation normalisés en puissance d'entrée H est le vecteur des coefficients du filtre adaptatif 28 SPS désigne l'opérateur de soustraction spectrale XH est le vecteur des échantillons d'entrée x(n) E est le vecteur des échantillons de sortie e(n), en lab- sence de double parole. Quel que soit l'algorithme d'adaptation utilisé, la soustraction spectrale n'apporte pas de complication notable, puisque les transformées directe et inverse sont de toute façon nécessaires. Le processus d'adaptation d'un filtre adaptatif, dans le domaine fréquentiel ou dans le domaine temporel, est perturbé par la présence simultanée d'une activité vocale aux terminaux proche et éloigné (situation dite de double parole). Pour éviter ces perturbations, il est courant dinter- rompre ou de ralentir le processus d'adaptation des coefficients du filtre en cas de double parole. Mais la plupart des détecteurs de double parole existants ont une réponse lente, sont souvent incapables de différencier la présence de double parole d'une modification du canal d'écho augmentant le niveau d'écho et/ou ne fonctionnent pas dans un environnement très bruité du fait qu'ils sont fondés sur une observation d'énergie. L'invention propose en conséquence un détecteur de double parole utilisable dans un annuleur d'écho du genre ci-dessus, utilisantun critère fiable d'identification de présence de parole du locuteur proche. Ce détecteur effectue une mesure de distorsion spectrale en comparant des paramètres spectraux extraits du signal de sortie de l'interface (dûs au locuteur éloigné) aux paramètres spectraux correspondants extraits du signal d'entrée (dû au bruit, à l'écho et au locuteur proche en cas de double parole). Un accroissement brusque de la discordance entre les signaux indique la présence d'une interférence de forte énergie, due à la parole proche ou au bruit.En conséquence, il est proposé un annuleur d'écho comprenant des moyens d'extraction de paramètres spectraux et des moyens d'évaluation de distorsion par comparaison des dits paramètres et de commande d'arrêt d'adaptation en cas de distorsion élevée. L'évaluation de distorsion peut utiliser différents critères. Un premier critère possible est la mesure de distance au sens d'Itakura-Saito entre deux modèles auto-régressifs (constituant l'extraction de paramètres) calculés sur les signaux d'entrée et de sortie. Une autre approche consiste à utiliser des distances cepstrales, utilisant des paramètres cepstraux déjà utilisés pour la reconnaissanec de parole, comme indiqué par exemple dans l'article de Alexandre et al shoot Cepstral Analysis: A Unified View - Application to Speech Processing in Car Noise Environnments" dans Speech Communication 12, pp 277288, 1993. L'invention sera mieux comprise à la lecture de la description qui suit d'un mode particulier de réalisation, donné à titre d'exemple non limitatif. La description se réfère aux dessins qui l'accompagnent, dans lesquels - la figure 1, déjà mentionnée, montre la constitution de principe d'un annuleur d'écho dans le domaine fréquentiel - la figure 2 similaire à une fraction de la figure 1, est un synoptique de principe montrant une modification possible de l'annuleur d'écho de la figure 1 suivant un aspect de l'invention Sur la figure 2, les éléments correspondants à ceux de la figure 1 sont désignés par le même numéro de référence. On supposera par la suite que l'annuleur est destiné à une installation téléphonique mains-libres, destinée à fonctionner dans la bande allant de 300 Hz à 3,4 Khz, avec échantillonnage à 8 Khz ; mais la même constitution pourrait être utilisée dans un terminal mains-libres prévu pour opérer en bande large de 16 Khz. Les blocs peuvent être par exemple de 256 échantillons, ce qui correspond à une durée de 32ms à la fréquence d'échantillonnage courante de 8 Khz. L'annuleur d'écho de la figure 2 ne comporte pas de moyens imposants une contrainte sur le gradient. Mais il comporte un soustracteur spectral 60 qui reçoit d'une part les coefficients de la transformée directe calculée sur un bloc, d'autre part les coefficients correspondants représentatifs du bruit, provenant d'un estimateur 62. Cet estimateur peut être remis à jour pendant des périodes de non activité vocale, par des moyens non représentés, tels que ceux décrits dans les documents mentionnés plus hauts. Le soustracteur 60 fournit ainsi, au circuit 50 d'adaptation des coefficients du filtre, une matrice de coefficient de transformation directe corrigée de l'effet du bruit. L'absence d'activité vocale peut être déterminée par différents moyens connus, fondés par exemple sur une mesure d'énergie dans un domaine fréquentiel déterminé, d'une part sur le signal arrivant au haut parleur 12, d'autre part sur le signal de sortie du microphone. Cependant, l'inhibition de l'estimateur 62 peut également être provoquée d'une part par mesure de l'énergie du signal attaquant le haut parleur, par détection d'une situation de double parole. Cette situation peut être notamment décelée par un bloc 64 de mesure de distorsion spectrale. Ce bloc comporte deux extracteurs de paramètres spectraux 66 et 68, recevant l'un le signal z(n) de sortie du convertisseur analogique, numérique 22, l'autre le signal montrant x(n). Le bloc comporte enfin un circuit 69 de mesure de similarité entre les paramètres représentatifs des signaux x(n) et z(n). Ces paramètres peuvent être notamment ceux mentionnés plus haut. En cas de détection d'un écart représentatif d'une forte distorsion, le circuit 68 applique un signal de blocage ou de ralentissement d'adaptation de la matrice diagonale m qui comporte les pas d'adaptation du filtre, sur une entrée 70 du circuit 50. L'annuleur d'écho peut facilement être prévu pour fournir un signal de parole débruitée en mettant simplement en oeuvre une transformée inverse supplémentaire 72 et un circuit de conversion parallèle-série 74. REVENDICATIONS 1. Annuleur d'écho acoustique disposé à être placé entre une interface acoustique mains libres et un réseau de communication, comprenant - un filtre adaptatif recevant les coefficients d'une transformée directe discrète dans le domaine fréquentiel, calculés sur des blocs successifs d'un signal d'entrée, et attaquant, par l'intermédiaire d'un circuit de transformée inverse, l'entrée soustractive d'un additionneur de formation d'un signal de sortie allant vers le réseau, et - un circuit d'adaptation des coefficients du filtre, caractérisé en ce que ledit circuit d'adaptation reçoit les coefficients des transformées successives des blocs d'échantillons du signal de sortie par l'intermédiaire d'un réducteur de bruit. 2. Annuleur d'écho selon la revendication 1, caractérisé en ce que le réducteur de bruit est un soustracteur spectral dont l'entrée soustractive reçoit les coefficients représentatifs du spectre estimé du bruit. 3. Annuleur d'écho selon la revendication 2, caractérisé en ce que ledit circuit d'adaptation des coefficients met en oeuvre l'algorithme : e e H(k+l)= H(k)+2F.g.F~. m(k).XH(k).SPS[E(k)] OÙ: k est l'ordre du bloc courant F désigne la matrice de transformée de Fourier discrète, ou DFT ; g est la matrice de contrainte sur le gradient m est la matrice diagonale des pas d'incrémentation normalisés en puissance d'entrée ; H est le vecteur des coefficients du filtre adaptatif 28 SPS désigne l'opérateur de soustraction spectrale XH est le vecteur des échantillons d'entrée x(n) E est le vecteur des échantillons de sortie e(n), en lab- sence de double parole. 4. Annuleur d'écho selon la revendication 1, 2 ou 3, caractérisé par des moyens pour contrôler l'adaptation en cas d'interférence d'énergie élevée due au bruit ou à un état de double parole. 5. Annuleur d'écho selon la revendication 4, caractérisé en ce qu'il comprend des moyens d'extraction de paramètres spectraux et des moyens d'évaluation de distorsion par comparaison des dits paramètres des signaux d'entrée (locuteur proche) et de sortie (locuteur éloigné) de l'interface. 6. Annuleur d'écho selon la revendication 5, caractérisé en ce que les dits moyens de comparaison mettent en oeuvre une mesure de la distance en sens dItakura-Saito entre deux modèles auto-régressifs, constituant l'extraction de paramètres calculés sur les signaux d'entrée et de sortie. 7. Annuleur d'écho selon la revendication 5, caractérisé en ce que les dits moyens de comparaison utilisent des paramètres cepstraux des signaux d'entrée et de sortie. 8. Annuleur d'écho selon l'une quelconque des revendications 1 à 6, caractérisé en ce qu'il comprend, de plus, un circuit de fourniture de signal de parole de locuteur proche débruité, le dit circuit étant prévu pour effectuer une transformée inverse (72), à partir du domaine fréquentiel et une conversion parallèle-série (74).