La présente invention se rapporte à un procédé de réduction d'écho mis en oeuvre avec un dispositif de prise de son multi-capteurs formant une antenne et un dispositif de restitution sonore. L'invention se rapporte en outre à un dispositif de réduction d'écho, destiné à fonctionner avec un dispositif de restitution sonore pour la mise en oeuvre d'un procédé selon l'invention. Le domaine fréquentiel auquel se rapporte l'invention est celui de la parole et plus particulièrement celui des signaux de parole basses fréquences.

Description de l'art antérieur L'invention trouve une application très avantageuse en ce qu'elle permet, par la réduction d'écho dans un dispositif de prise de son d'un système de téléconférence, de s'affranchir d'un microphone déporté pour permettre des communications dites mains- libres.

L'écho acoustique est l'obstacle majeur au bon fonctionnement des terminaux de communication mains-libres. Cet écho acoustique résulte de la capture par les capteurs du dispositif de prise de son, d'une partie du signal émis par les transducteurs du dispositif de restitution sonore. Les capteurs et les transducteurs utilisés en acoustique aérienne sont respectivement des microphones et des hauts parleurs.

Dans un terminal de communication mains-libres, les microphones sont en couplage acoustique, et éventuellement en couplage mécanique avec le système de restitution du son soit de manière directe lorsque les microphones et le système de restitution sont dans la mme coque, soit de manière indirecte comme dans le cas d'un terminal de visioconférence posé sur un moniteur de télévision. L'écho résultant se décompose en deux parties : -la réponse de la salle, -le couplage direct haut parleur/microphone.

La réponse de la salle est en général traitée efficacement par des techniques classiques de contrôle d'écho telles que l'annulation d'écho ou la variation de gain, en particulier du fait du niveau assez faible de cette partie de l'écho. Il n'en est pas de mme pour le couplage direct qui résulte à la fois du trajet acoustique aérien mais également des propagations de vibrations à travers la coque du terminal de

communication ainsi que d'éventuelles résonances de cavités ou d'éléments mécaniques dans le cas de coffrets de moniteurs de télévision. Cette partie de l'écho, en général supérieure de 6 à 20 dB à la parole locale captée par les microphones, entraîne 1'effet Larsen si aucun traitement n'est appliqué.

Si les techniques connues d'annulation d'écho ne permettent pas de traiter de manière satisfaisante cet effet néfaste, une variation de gain importante doit tre implantée pour éviter le départ en Larsen. Cette variation de gain doit intervenir d'une part pendant la phase de convergence du dispositif d'annulation d'écho et d'autre part à la suite d'une variation du chemin d'écho ou pendant les moments de double parole.

Du fait de cette contrainte, le terminal de communication est rendu peu interactif et l'utilisateur peut considérer que la fonction mains libres n'est pas performante.

Pour pallier cet inconvénient, il existe plusieurs techniques qui portent sur la réduction du couplage direct.

Une première technique s'attache à réduire le couplage mécanique. La solution habituellement mise en oeuvre consiste à faire usage de matériaux antivibratoires tels que des mousses ou du caoutchouc qui désolidarisent le microphone de la coque. Ces matériaux ont pour effet de supprimer ou de diminuer fortement les vibrations de la coque engendrées par le haut parleur et leur propagation vers le microphone.

La solution du découplage mécanique diminue effectivement le couplage par les vibrations, mais elle s'avère une solution industriellement chère. Par ailleurs, elle ne réduit pas le couplage acoustique qui peut se révéler important si le microphone est placé près du haut parleur, ce qui est aujourd'hui la tendance du fait de l'attrait pour des terminaux de communication de taille réduite. En outre, lorsque le système de restitution fait partie intégrante d'un ensemble comprenant le dispositif de prise de son, comme dans le cas d'un moniteur de télévision incluant des haut-parleurs et des microphones ou d'un moniteur sur lequel est placé un dispositif incluant des microphones, il n'est pas envisageable de désolidariser les hauts parleurs de la coque de 1'ensemble.

Une seconde technique consiste à utiliser un filtre numérique de compensation dont la réponse est le résultat du calcul de l'inverse de la réponse impulsionnelle du couplage mécanique et acoustique entre le microphone et le haut parleur.

En théorie, à la sortie du filtre de compensation le signal en provenance du haut parleur dû à l'écho est annulé. En pratique, cette technique ne donne pas satisfaction si le couplage change dans ses caractéristiques, mme d'une façon minime, suite par exemple à un montage/démontage du terminal de communication ou suite à un

changement des caractéristiques des microphones. Par ailleurs cette technique n'est pas adaptée si la nature des perturbations n'est pas linéaire, c'est-à-dire si les perturbations ne peuvent pas tre modélisées par un produit de convolution entre le signal provenant du haut parleur et une réponse impulsionnelle de filtrage. Qui plus est, le fonctionnement des transducteurs est rarement linéaire, car ils sont généralement sujets à des distorsions et/ou à des saturations, qui sont des exemples types de fonctions non linéaires. Enfin, lors d'une production en grande série, le couplage est forcément différent d'un terminal de communication à un autre, du fait par exemple de l'utilisation de haut-parleurs ou de microphones légèrement différents entre eux.

Finalement, ces limitations rendent cette technique de compensation par filtre peu efficace.

Il existe d'autres systèmes de pondération de microphones basés sur une equi- répartition des microphones autour d'un haut-parleur associée à une pondération simple des microphones ou à un déphasage des microphones. De tels systèmes font l'objet du brevet français 93 02054. Ces systèmes sont destinés aux terminaux de communication de groupe qui possèdent une symétrie géométrique bien particulière due à leur prise de son omnidirectionnelle. L'efficacité de ces systèmes est très sensible aux disparités des microphones ou des voies d'acquisition. Ces systèmes nécessitent par conséquent un tri des composants et un calibrage précis des microphones qui peut tre mis en défaut au cours du temps. En outre, l'ajout de cadran de visualisation, de clavier, de caméra pivotante peut rompre les propriétés de symétrie des hauts parleurs vis-à-vis des microphones. Par ailleurs, dans le cas où le nombre de microphones N dépasse deux ces dispositifs ne sont pas optimaux car N degrés de liberté sont disponibles pour satisfaire deux contraintes, la fixation d'un gain dans une direction donnée et l'annulation du couplage direct. Enfin, comme le signal capté correspond, vis-à-vis de la source utile, à 1'ensemble des contributions de chaque microphone et comme ces contributions sont déphasées entre elles et sont éventuellement d'amplitudes différentes, le signal résultant a un spectre dégradé en hautes fréquences.

Un autre système visant des applications plutôt individuelles, connu par exemple du brevet français 98 14321, offre l'avantage de prendre en compte les spécificités du terminal final, c'est-à-dire du terminal de communication après son intégration. Ce système met en oeuvre deux microphones positionnés à des distances différentes du haut-parleur. Les contributions des deux microphones sont combinées par pondération afin d'annuler l'onde directe de couplage. La pondération différente des deux signaux provenant des deux microphones permet de capter les sources utiles

situées dans un champ lointain. Le calcul du filtre de voie appliqué au deuxième microphone peut tre obtenu en inversant la réponse d'une mesure de couplage entre le microphone et le haut parleur. Un tel calcul est généralement suffisant dans la plupart des situations puisque les microphones sont dans le champ direct du haut parleur.

Cependant, cette solution nécessite l'existence d'une différence d'amplitude suffisante de l'onde de couplage entre les deux microphones. Pour obtenir une telle différence, il est nécessaire de positionner les microphones suffisamment près du haut parleur ce qui entraîne des conditions très défavorables en terme de couplage. Alors que cette condition n'est pas gnante pour un système individuel pour lequel cette disposition est imposée, il n'en est pas de mme pour d'autres systèmes tels que des terminaux de visioconférence de type"Set-top-box", c'est-à-dire des terminaux dont la forme permet de les poser sur un moniteur de télévision. Suivant ce principe, plus les microphones sont éloignés du haut parleur plus la distance d'espacement inter capteurs devra tre grande. Ceci conduit à une limitation de la bande passante et à une sensibilité plus grande aux obstacles situés à proximité. Enfin, une amplification du bruit incohérent est à déplorer et n'est contrôlable qu'en modifiant la distance séparant les deux microphones.

Finalement, la directivité de ces dispositifs n'est que partiellement contrôlable et est peu différente de la directivité propre des microphones.

Parallèlement à ces solutions, il existe des techniques multi capteurs qui cherchent à se prémunir du couplage acoustique en rendant le système globalement directif. En effet, lorsque la source sonore est située dans un environnement bruyant ou réverbérant, la directivité d'un capteur seul peut ne pas tre suffisante pour extraire le signal du bruit.

Pour remédier à ce problème, une méthode consiste à associer N capteurs pour former une antenne acoustique. Cette association consiste à additionner en cohérence de phase les signaux de sortie des différents capteurs pour une direction (, S) donnée.

Cette addition est une des techniques regroupées sous l'appellation formation de voies.

Un capteur, ou un transducteur, est caractérisé par sa directivité dans 1'espace qui est habituellement mesurée dans deux plans orthogonaux et représentée en coordonnées polaires (r, ?, ) sous la forme de deux diagrammes. La directivité d'un capteur donne une"image"du niveau de signal capté par ce capteur pour une source sonore ponctuelle située dans une direction (o, 0) et à une distance r du centre du capteur. La directivité de l'antenne obtenue à l'issue de la formation de voie a des performances meilleures que la directivité d'un seul capteur. Les N microphones, séparés entre eux

d'une distance d, captent la pression en une région localisée de l'espace et opèrent un échantillonnage spatial du champ sonore. Cette méthode offre en hautes fréquences une efficacité certaine, mais en basses fréquences (f < 1000 Hz), la difficulté persiste pour plusieurs raisons : -un manque de directivité du fait de la taille réduite de l'antenne, -une moins bonne robustesse du fait de la sensibilité importante aux disparités entre capteurs en basses fréquences, du fait de 1'effet éventuel de l'écran de télévision, -un diagramme directionnel de l'antenne pour une source située à 1'emplacement des enceintes différent du diagramme directionnel à 3 m car l'hypothèse de champ lointain n'est pas vérifiée dans le premier cas.

Ainsi un but de l'invention est de réduire l'écho produit par le couplage acoustique entre le dispositif de prise de son et le dispositif de restitution sonore sans les inconvénients des dispositifs et méthodes connus.

Résumé de l'invention A cet effet, l'invention a pour objet un procédé de réduction d'écho, mis en oeuvre avec un dispositif de prise de son multi-capteurs formant une antenne et un dispositif de restitution sonore, qui consiste à soumettre les signaux de sortie des capteurs à des pondérations complexes w (f), les dites pondérations w (f) étant calculées par une maximisation sous contraintes aux basses fréquences et en champ proche du facteur de directivité Fd (f) dont 1'expression est la suivante : le calcul étant tel que les dites pondérations w (f) satisfont une première contrainte linéaire (3) sur le module et la phase de la fonction de transfert du dispositif de prise de son dans des directions données, la formulation de cette première contrainte à chaque fréquence f étant la suivante : CH(~)w(~)=s(~) et les dites pondérations w (f) satisfont une deuxième contrainte (4) déterminée à partir de mesures"in situ"des fonctions de transfert complexes de canaux acoustiques délimités par les entrées des haut-parleurs du dispositif de restitution

sonore et les sorties des capteurs composant le dispositif de prise de son, cette contrainte ayant pour expression : MH (f) W (w) = O (4) les formules (1) à (4) étant telles que : H est le vecteur de propagation dont les éléments sont les valeurs complexes des fonctions de transfert en champ libre entre une source ponctuelle située à la distance r du centre de l'antenne dans la direction définie par les angles polaires (p et A et chaque capteur de l'antenne, calculées à la fréquence f, D, (f) est la matrice de directivité qui caractérise les propriétés de sélectivité spatiale de l'antenne à la distance r, W (f, 0, 0) sont des pondérations spatiales permettant une atténuation plus importante des ondes provenant des directions des haut-parleurs pour réduire le couplage direct avec ceux-ci, w (f) est un vecteur des pondérations complexes des signaux de sortie des capteurs à la fréquence f, C (f) est une matrice dite"de contraintes"contenant les vecteurs de propagation théoriques calculés à partir d'un modèle de propagation de type champ libre ou mesurés en condition de champ libre, s (f) sont les gains complexes désirés dans des directions données, M (f) est une matrice de contrainte de couplage contenant les fonction de transfert complexes mesurées"in situ"des canaux acoustiques délimités par les entrées des haut-parleurs du dispositif de restitution sonore et les sorties des capteurs composant le dispositif de prise de son et, 0 est un vecteur nul.

Ainsi, l'invention prend en compte des contraintes sur le couplage direct lors du calcul de la formation de voies. Le calcul de la formation de voies n'est pas effectué uniquement sur la base de modèles théoriques de propagation ou de mesures effectuées en chambre sourde mais aussi sur la base d'estimées des fonctions de transfert obtenues "in situ". Ainsi, l'invention permet de réduire le couplage direct et/ou semi-direct tout en imposant une directivité souhaitée et en contrôlant l'amplification maximale du bruit incohérent.

L'invention s'applique en particulier à des dispositifs de prise de son pour des visioconférences de groupe.

L'invention a pour avantage de prendre en compte automatiquement 1'environnement acoustique proche (forme du moniteur de télévision et du boîtier

contenant l'antenne, position de celui-ci au regard du moniteur, obstacle proche ou présence d'un mur à proximité), la réponse des haut-parleurs et la disparité des capteurs. En outre, l'optimisation étant effectuée en aveugle, le procédé selon l'invention ne nécessite pas de faire des hypothèses préalables sur la connaissance de la position des haut-parleurs, sur les caractéristiques électro-acoustiques des haut-parleurs et des capteurs.

Brève description des dessins D'autres caractéristiques et avantages de l'invention apparaîtront lors de la description qui suit faite en regard de figures annexées données à titre d'exemples non limitatifs.

La figure 1 est un organigramme d'une première mise en oeuvre d'un procédé selon l'invention.

La figure 2 est un repère en coordonnées polaires.

La figure 3 est un organigramme d'une seconde mise en oeuvre d'un procédé selon l'invention.

La figure 4 est un schéma d'une réalisation particulière d'un dispositif de réduction d'écho fonctionnant avec un système acoustique pour la mise en oeuvre d'un procédé selon l'invention.

La figure 5 est un synoptique d'une mise en oeuvre d'un procédé selon l'invention.

Description détaillée d'un mode de réalisation Le procédé de réduction d'écho selon l'invention est mis en oeuvre avec un dispositif de prise de son multi-capteurs formant une antenne et un dispositif de restitution sonore. Le procédé selon l'invention consiste à soumettre les signaux de sortie des capteurs du dispositif de prise de son à des pondérations complexes w (f), les dites pondérations étant calculées par une maximisation sous contraintes aux basses fréquences et en champ proche du facteur de directivité Fd (f). Le procédé se déroule suivant plusieurs étapes décrites ci-après et illustrées par la figure 1.

Dans une première étape 1, le procédé selon l'invention consiste à déterminer 1'expression de la première contrainte (3) linéaire de gain désiré en champ lointain dans les directions correspondant à des sources utiles (locuteurs), définies par l'équation matricielle : CH(~)w(~)=s(~) équation dans laquelle,

C (f) est la matrice dite"de contraintes"contenant les vecteurs de propagation théoriques calculés à partir d'un modèle de propagation de type champ libre ou mesurés en condition de champ libre, w (f) est le vecteur des pondérations complexes des signaux de sortie des capteurs à la fréquence f, s (f) sont les gains complexes désirés dans des directions données.

Dans une deuxième étape 2, le procédé selon l'invention consiste à déterminer l'expression de la deuxième contrainte (4) déterminée à partir de mesures"in situ"des fonctions de transfert complexes de canaux acoustiques délimités par les entrées des haut-parleurs du dispositif de restitution sonore et les sorties des capteurs composant le dispositif de prise de son, cette contrainte étant définie par l'équation matricielle : MH (H) w (/) =0 (4) équation dans laquelle, M (f) est une matrice de contraintes de couplage contenant les fonction de transfert complexes mesurées in situ des canaux acoustiques délimités par les entrées des haut-parleurs du dispositif de restitution sonore et les sorties des capteurs composant le dispositif de prise de son et, 0 est un vecteur nul.

Dans une troisième étape 3, le procédé selon l'invention consiste à maximiser le facteur de directivité (1) en champ proche.

Le facteur de directivité en champ proche Fd est défini par l'équation : équations dans lesquelles, H est le vecteur de propagation dont les éléments sont les valeurs complexes des fonctions de transfert en champ libre entre une source ponctuelle S, représentée sur la figure 2, située à la distance r du centre de l'antenne dans la direction définie par les angles polaires (p et 9 et chaque capteur de l'antenne, calculées à la fréquence f, D (f) est la"matrice de directivité"qui caractérise les propriétés de sélectivité spatiale de l'antenne à la distance r,

W (f, ç, O) sont des pondérations spatiales permettant une atténuation plus importante des ondes provenant des directions des haut-parleurs pour réduire le couplage direct avec ceux-ci.

Lors de la maximisation du facteur de directivité, les contraintes suivantes peuvent en outre tre prises en compte à titre facultatif lors de la mise en oeuvre du procédé, un organigramme du procédé correspondant est illustré par la figure 3 : -une troisième contrainte (4) non linéaire de robustesse au bruit incohérent dont le facteur de réduction minimal R, i (f) fixe la valeur, suivant l'équation : -des quatrièmes contraintes (5) non-linéaires sur l'ouverture du lobe principal et/ou sur les gains en amplitude dans des directions données et sur la réduction des lobes d'aliasing et des lobes secondaires. Ces contraintes permettent en premier lieu de garantir une largeur suffisante du lobe principal de l'antenne pour couvrir correctement la zone angulaire où peuvent se trouver les directions des sources utiles, elles complètent en ce sens la contrainte linéaire correspondant à l'équation (3). En second lieu, elles permettent de maintenir actifs des microphones à des fréquences pour lesquelles l'espacement est trop grand pour respecter les conditions de Shannon, en limitant le niveau des lobes d'aliasing en résultant ainsi que des lobes secondaires. Ces contraintes permettent ainsi d'aboutir à des facteurs de directivité plus élevés notamment dans les bandes fréquentielles où certains capteurs passent de l'état actif à l'état inactif. Elles autorisent ainsi la mise en oeuvre d'un nombre plus limité de capteurs, donc une optimisation du coût du dispositif. Elles s'expriment sous la forme d'un système d'équations non linéaires : wH(~)H0(~)H0H(~)w(~) = g0(~) , (f) HI (f) H" (f) w (f) = gi (f) wH(~)Hk(~)HkH(~)w(~) = gk(~) système d'équations dans lequel, Hi, i=0, ..., k représente le vecteur de propagation entre une source ponctuelle en champ lointain dans la direction angulaire ((pi, Oi) et les microphones de l'antenne et, gi le gain désiré dans cette direction. Ainsi, go représente le gain désiré, par exemple-3 dB, dans la direction correspondant à l'ouverture du lobe principal, par exemple 30 degrés par rapport au centre de celui-ci.

Ces différentes contraintes facultatives sont prises en compte en modifiant 1'expression de la matrice de directivité suivant la formule : où ll (f) est un coefficient de Lagrange qui est ajusté pour satisfaire la contrainte de robustesse (5), et ai (f), i=0,..., k est un jeu de coefficients de Lagrange ajustés pour satisfaire les contraintes non-linéaires sur l'ouverture du lobe principal et sur la réduction des lobes d'aliasing et des lobes secondaires (6).

La solution optimale du problème est donnée par la formule : w (f) =Dr (f) C (~)(CH(~)Dr@-1(~)C(~))-1s(~) Une telle optimisation permet de réduire la sensibilité de l'antenne pour des sources proches sans modifier de manière significative les performances de l'antenne en champ lointain.

Le schéma de la figure 4 est un exemple d'un dispositif 6 de réduction d'écho fonctionnant avec un système 7 acoustique pour la mise en oeuvre d'un procédé selon l'invention.

Le dispositif 6 de réduction d'écho se présente sous la forme d'une boîte de type "Set-top-box". La boîte comprend un dispositif 8 de prise de son aérien et un module 9 électronique de calcul. Le dispositif 8 de prise de son aérien consiste en une antenne 10 formée, suivant 1'exemple retenu, par onze microphones 11 m, vus par transparence sur le schéma, répartis suivant trois plans parallèles, m=l, 2,... N avec N=11. Le module 9 électronique de calcul comprend des moyens de calcul pour calculer les différentes contraintes et, en particulier, la première contrainte exprimée par la relation (3) et la deuxième contrainte exprimée par la relation (4). Les moyens de calcul consistent typiquement en des circuits intégrés de traitement du signal qui effectuent les calculs en temps réel. Selon une variante, ils peuvent consister en un processeur d'usage général suffisamment rapide tel que ceux actuellement intégrés dans les terminaux de téléconférence évolués.

Le système 7 acoustique comprend un moniteur 12 de télévision et un dispositif 13 de restitution sonore constitué de deux hauts-parleurs disposés de part et d'autre du moniteur de télévision.

En basses fréquences, du fait d'indices de superdirectivité élevés, l'antenne 10, bien que pointant vers un locuteur situé en champ lointain, possède une sensibilité accrue pour les sources situées en champ proche. Les haut-parleurs constituent des sources proches de mme que les sources images des hauts-parleurs lorsque le système

7 est disposé près d'un mur. Ces réflexions/diffusions sont fortement énergétiques du fait de leur précocité. Cet effet est d'autant plus notable que dans ces bandes de fréquences, les ondes acoustiques ne subissent qu'une faible atténuation par absorption des surfaces réfléchissantes.

Pour réduire la sensibilité de l'antenne 10 pour des sources proches sans modifier de manière conséquente les performances de l'antenne 10 en champ lointain, le dispositif 6 de réduction d'écho maximise sous contraintes le facteur de directivité en champ proche défini par la relation (1) selon le procédé selon l'invention. Dans la relation (2) les directions 0 et 71, qui correspondent aux directions des hauts parleurs du système 7, sont pondérées de façon plus importante par la fonction de pondération spatiale w de telle sorte que le couplage acoustique direct soit minimisé.

La mise en oeuvre d'un procédé selon l'invention avec le dispositif de réduction d'écho et le dispositif de restitution sonore est détaillé en regard des figures 4 et 5.

Après mise sous tension du moniteur 12 de télévision et de la boîte 6, un signal de mesure x (n) généré par un générateur 14 est converti en analogique par un convertisseur 15 numérique analogique. La sortie x (t) du convertisseur 15 est transmise séparément sur chacun des haut-parleurs 13 du moniteur 12 de télévision.

Selon 1'exemple, le générateur 14 génère une séquence pseudo-aléatoire de longueur maximale, dite MLS. Le signal Vn, (t) issu de chaque microphone llm constituant l'antenne acoustique est numérisé par un convertisseur 16m analogique numérique sous la forme d'un signal numériqueym (n) La réponse impulsionnelle estimée hestm (n) de chaque canal acoustique haut- parleur/microphone m, m=l, 2,... N avec N=11, est déterminée à partir du signal y,,, (n) par un calculateur 17m effectuant un calcul du domaine de l'état de l'art. Lorsque le signal de mesure x (t) est une séquence MLS, la réponse impulsionnelle estimée hestm (n) peut tre déterminée par intercorrélation en utilisant en particulier une transformée de Hadamard.

Un fentrage 18m temporel est pratiqué sur la réponse impulsionnelle estimée hestm (n) de chaque canal acoustique haut-parleur/microphone, afin de ne tenir compte que du trajet direct et d'un petit nombre de réflexions précoces désignés respectivement dans la suite par les termes couplage direct et couplage semi-direct. Ainsi, cette réduction du couplage s'appuie sur les éléments du couplage les plus stables (trajet acoustique direct et premières réflexions sur les obstacles fixes, par exemple les murs à proximité) et écarte les parties des éléments du couplage susceptibles de se modifier rapidement dans le temps (par exemple en raison des mouvements et déplacements des

locuteurs). A cet fin, le maximum à court terme de la réponse impulsionnelle du microphone central est recherché et le numéro no de l'échantillon correspondant est relevé. Pour chaque réponse impulsionnelle, il est procédé à une transformée de Fourier rapide (FFT) sur une zone temporelle délimitée par une fentre d'analyse centrée sur l'échantillon no.

Dans le cas où plusieurs contraintes de couplage sont prises en compte, il y a autant de fentres d'analyse que de contraintes de couplage. La fentre d'analyse peut tre rectangulaire ; cependant d'autres formes de fentres telles que Hamming, etc peuvent tre avantageusement utilisées. L'étendue de la fentre d'analyse est choisie courte, typiquement de l'ordre de 64 à 256 échantillons, soit de 4 à 16 millisecondes pour une fréquence d'échantillonnage fe = 16 kHz. Il est donc nécessaire d'ajouter 19m des zéros à la suite de la fentre d'analyse pour améliorer la résolution fréquentielle avant de procéder à la FFT 20m. Le résultat de la FFT est la fonction de transfert complexe Hest directm (0Z) correspondant à chaque canal acoustique haut-parleur/ microphone m de l'antenne, après pondération temporelle privilégiant le couplage direct et/ou semi-direct.

Les coefficients w (f) de pondération de la formation de voie sont calculés 21, 22 à partir de ces différentes fonctions de transfert suivant la méthode d'optimisation du facteur de directivité en champ proche décrite en regard de la figure 1. La matrice M (f) de la relation (4) contient les fonctions de transfert complexes de couplage obtenues précédemment, ces fonctions correspondant aux vecteurs de propagation mesurés. La solution optimale est calculée suivant 1'expression (8) dans laquelle la matrice C (f) contient à la fois les vecteurs de propagation théoriques de type champ libre, c'est-à-dire les contraintes en champ lointain pour les sources utiles, et les vecteurs de propagation mesurés. Le vecteur de contrainte s (f) contient typiquement des retards purs représentant la contrainte (3) sur les locuteurs utiles, et des zéros pour la contrainte sur le couplage suivant l'équation (4).

Les données numériques suivantes données à titre d'exemple sont extraites de mesures effectuées sur des enregistrements réalisés dans une salle au moyen d'un dispositif 6 de réduction d'écho fonctionnant avec un système 7 acoustique comme représentés sur la figure 4. Le procédé selon l'invention a été mis en oeuvre pour des fréquences inférieures ou égales à 1400 Hz. Les analyses fréquentielles ont été effectuées sur une fentre de 256 points étendue à 1024 points par l'ajout de zéros comme décrit précédemment. Les filtres de voies calculés sont de type FIR à 152

coefficients, ce qui permet d'obtenir une approximation satisfaisante du gabarit de filtrage souhaité.

Une première série de mesures prise comme référence se rapporte à la réponse impulsionnelle obtenue en sortie de la formation de voie correspondant à la maximisation sans contrainte du facteur de directivité en champ lointain. Une deuxième série de mesures se rapporte à la réponse impulsionnelle obtenue en sortie de la formation de voie calculée après prise en compte d'une seule contrainte sur le couplage direct. Une troisième série de mesures se rapporte à la réponse impulsionnelle obtenue en sortie de la formation de voie calculée après prise en compte de deux contraintes sur le couplage au moyen de deux fentres d'analyse : la première est centrée sur l'onde directe et la seconde est décalée par rapport à la première de telle sorte qu'elle soit centrée sur une réflexion précoce.

La comparaison de la réponse impulsionnelle obtenue soit par application de deux contraintes, soit par application d'une seule contrainte sur le couplage direct, avec la réponse impulsionnelle obtenue au moyen de la maximisation sans contrainte du facteur de directivité en champ lointain, fait ressortir la réduction de 1'effet du couplage. Cette réduction est principalement mesurable sur le début de la réponse impulsionnelle, bien qu'une atténuation de la réponse impulsionnelle persiste au-dela.

En effet, le procédé selon l'invention permet non seulement d'atténuer l'effet du couplage direct mais en outre il est efficace pour atténuer 1'effet du couplage global.

Une analyse en fréquence indique que la réduction d'écho obtenue au moyen des deux contraintes est de l'ordre de 6 à 8 dB pour les fréquences inférieures à 250 Hz.

L'examen des échogrammes correspondant à la première série de mesures et à la troisième série de mesures confirme une atténuation importante au niveau de l'onde directe pouvant aller jusqu'à 20 dB. Le tramage en fin de réponse impulsionnelle est également réduit de façon importante.

L'émission d'une séquence de longueur maximale MLS lors de la mise sous tension du système, pour mesurer les réponses impulsionnelles de couplage, peut fournir une estimation robuste mme si simultanément les locuteurs se mettent à parler ; la séquence MLS étant décorrélée du signal de parole. Néanmoins, afin d'éviter que la variance introduite de ce fait dans l'estimée des fonctions de transfert ne dégrade les performances du procédé, des précautions peuvent tre prises telles que le rejet de la mesure si une parole locale a été détectée au niveau de la prise de son, le moyennage de plusieurs mesures, etc. L'estimation des fonctions de transfert est peu sensible au mouvement des locuteurs tant que ceux-ci ne se trouvent pas à proximité directe du

système puisque seul le couplage direct ou proche est analysé. En outre, le traitement ayant lieu uniquement en basses fréquences, il y a peu de chance de modifier de manière importante les fonctions de transfert lors du déplacement au cours de la communication d'un obstacle de dimensions réduites.

Il est à noter que le nombre de contraintes possibles, y compris les contraintes sur le gain en champ lointain, dépend du nombre de degrés de liberté disponibles. Le nombre de contraintes exactes (en module et en phase) est en théorie inférieur à N-l, le nombre de contraintes partielles (sur le module ou sur la phase exclusivement), est en théorie inférieur à 2N-1, N étant le nombre de microphones. Il est nécessaire de limiter le nombre de contraintes sous peine de dégrader sensiblement les propriétés directionnelles de l'antenne en champ lointain. Néanmoins, si le nombre de contraintes s'avère supérieur, le problème d'optimisation ne peut plus tre résolu de façon exacte et il convient de minimiser l'écart quadratique entre le vecteur C (f) w (f) et le vecteur de contrainte désiré s (f).

ANNEXE 1 Directivité Un capteur est caractérisé par sa directivité dans 1'espace qui est habituellement mesurée dans deux plans orthogonaux et représentée en coordonnées polaires (r, (p, 0) sous la forme de deux diagrammes. Lorsque le capteur est de révolution autour d'un axe, la directivité mesurée dans un plan passant par cet axe est alors dite directivité en site. La directivité mesurée dans le plan orthogonal est dite directivité en gisement.

Habituellement, seule la directivité en site dans le plan passant par le maximum de la directivité en gisement est mesurée. La directivité donne une image du niveau de signal capté par le capteur pour une source sonore ponctuelle située dans une direction ((p, 0) et à une distance r du centre du capteur.

Les microphones utilisés en acoustique aérienne présentent généralement une symétrie de révolution autour de leur axe de référence. Par définition, le microphone sera dit : -omnidirectionnel si sa sensibilité est indépendante de la direction d'incidence de l'onde acoustique, -bidirectionnel d'ordre n si sa sensibilité peut se mettre sous la forme : où A (m) est une variable fonction de la fréquence et n un entier positif, -unidirectionnel d'ordre n si sa sensibilité peut s'écrire sous la forme : avec ß constante positive appelée constante de cardioïcité. Le cas (3=1 correspond a une directivité cardioïde, p=/3a une directivité supercardioïde et ß=3 à une directivité hypercardioïde.

Les antennes sont généralement constituées de microphones discrets. Les N microphones séparés entre eux d'une distance d captent la pression en un endroit localisé de l'espace et opèrent un échantillonnage spatial du champ sonore. Le comportement de ces antennes discrètes peut s'analyser par analogie avec la numérisation des signaux échantillonnés sur une période T : T r d ou encore T # # pour le cas tridimensionnel La fonction d'excitation qui traduit le champ de pression reçu par les microphones, pour une antenne acoustique dont les capteurs ont la mme sensibilité et ont un mme coefficient de pondération, s'écrit sous la forme suivante :

En procédant à la transformée de Fourier spatiale, il est possible d'en déduire le spectre spatial de l'antenne Q (u/2) Une autre interprétation de cette expression formelle peut tre obtenue en considérant l'antenne comme un échantillonneur spatial tronqué par une fentre rectangulaire. En conséquence, la fonction d'excitation peut s'écrire sous la forme : L correspond à la longueur"effective"de l'antenne qui diffère de la longueur physique Lp = (N-l) d.

Par suite, le spectre spatial de l'antenne peut s'écrire sous la forme d'un produit de convolution : Le premier terme conduit à un sinus cardinal (correspondant au spectre spatial <BR> <BR> <BR> de l'antenne continue) :<BR> sin(#Lu/#)<BR> TF{rectL(x)}= (15)<BR> <BR> <BR> <BR> <BR> <BR> #Lu/# Quand au second, il a pour expression : La transformée de Fourier d'un train d'impulsions équidistantes de d dans 1'espace physique conduit également à un train d'impulsions équidistantes de kId dans le domaine des fréquences spatiales.

La convolution de l'équation (15) avec l'équation (16) permet d'obtenir le diagramme directionnel Q (u) de l'antenne :

Le diagramme polaire correspond à la somme de sinus cardinaux décalés d'un intervalle de k/d. Le premier maximum centré en u=0 correspond au lobe principal de l'antenne de largueur B/L. L'allure du diagramme de la fonction TF I rectl (x)) permet de tirer les conclusions suivantes : -l'antenne sera d'autant plus directive que sa longueur L sera importante, -par définition, la gamme de fréquence spatiale appelée « partie visible » est délimitée par-l<u<l. Si la distance intercapteur devient trop grande au regard de la longueur d'onde, un deuxième lobe apparaît dans la plage visible : c'est l'aliasing spatial. Ce phénomène est tout à fait analogue au repliement de spectre dans le domaine temps-fréquence lorsque le théorème de Shannon n'est pas respecté, -Si d est suffisamment petit, le diagramme de directivité de l'antenne discrète est très proche de celui de l'antenne continue.

Un calcul direct de la fonction donnée par 1expression (12) normalisée (en remarquant qu'il s'agit d'une suite géométrique) conduit à un déphasage près à 1'expression : Cette expression est très proche de celle obtenue pour une antenne continue dans le cas où d<<X/2.

La discrétisation de l'antenne a pour conséquence de périodiser le diagramme de directivité de l'antenne continue et impose donc une limitation en hautes fréquences du fait de la présence d'aliasing spatial.

Facteur de directivité Pour rendre compte des propriétés directionnelles d'un microphone, une méthode consiste à définir son facteur de directivité, lequel caractérise la faculté qu'a le microphone à privilégier une onde perpendiculaire à l'axe de référence (la direction visée) et se dirigeant vers la face avant du microphone, par rapport à l'ensemble des ondes pouvant provenir de toutes les directions de l'espace. Le facteur de directivité s'exprime comme le rapport du carré de la tension produite en sortie du microphone pour une onde arrivant dans l'axe par la moyenne des carrés des tensions produites par les ondes de mme pression efficace et de mme fréquence provenant dans un angle

solide de 4 TE. Elle peut donc s'exprimer également comme le rapport du carré de l'efficacité en champ libre par l'efficacité en champ diffus. Le champ diffus correspond à l'état d'un champ dans lequel l'énergie volumique moyenne a la mme valeur en tous les points et où l'intensité acoustique est la mme dans toutes les directions.

Champ lointain/champ proche Soit une antenne linéaire de longueur L, centrée en x = O, disposée le long de l'axe des abscisses x. Soit un point S de coordonnées (xo, yo) tel que la direction OS fait un angle (p avec l'axe des abscisses. Le signal reçu par l'antenne placée dans un champ acoustique créé par une source placée au point S correspond à l'ensemble des signaux captés par chaque récepteur élémentaire R (x, y) de l'antenne. Le signal électrique s (xo, yo) capté en sortie d'antenne et hé à la présence de la source S (xo, yo) a pour expression : où K est une constante indépendante de x et i (x) est le champ de pression reçu par chaque récepteur élémentaire.

La distance p séparant la source S (xo, yo) de chaque récepteur élémentaire R (x, y) de l'antenne s'écrit : Dans le cas où la source se trouve suffisamment éloignée de l'origine au regard de la dimension de l'antenne (###=#0#L), la dernière expression peut se développer sous la forme d'une série binomiale en ne gardant que les termes prédominants : Pour le cas de sources en champ lointain (po>>L), le terme d'amplitude présent au dénominateur de l'équation (19) peut tre remplacé par le terme po. Il vient :

En posant u=cossp, il vient 1'expression normalisée par rapport à l'axe du diagramme de directivité de l'antenne linéaire : ou 1, correspond au champ de pression qui satisfait à : ##(k,0)#=1.

Le premier terme de l'exponentielle sous l'intégrale représente l'intégrale de Fourier spatiale. Celle-ci est associée au diagramme de directivité de l'antenne en champ lointain. Le deuxième terme est nécessaire pour rendre compte du comportement directionnel de l'antenne lorsque la source est en champ proche. Ce terme correctif permet de tenir compte de la courbure du front d'onde. Par analogie avec la théorie de la diffraction de Kirchhoff-Fresnel, le premier terme de 1'expression ci-dessus est donc l'analogue de la zone de Fraunhofer (approximation du premier ordre en x/ro). En outre, la taille limitée de l'antenne peut tre interprétée comme 1'effet d'un écran virtuel parfaitement anéchoïque masquant une antenne de longueur infinie et les deux hypothèses de Kirchhoff sont en conséquence respectées. Il n'est donc pas surprenant d'observer dans ce cas une grande similitude dans les résultats. En l'occurrence, à l'instar de la théorie de la diffraction qui conduit à affirmer que la figure de diffraction est la transformée de Fourier double de la fonction d'ouverture, voir l'ouvrage de M. Bruneau"Introduction aux théories de l'acoustique", éditions Université du Maine, Le Mans, 1983 p. 370, la figure de directivité de l'antenne est la transformée de Fourier spatiale de la fonction d'excitation. Cette transformée de Fourier est monodimensionnelle si l'antenne est rectiligne, double si l'antenne est planaire et voire triple si l'antenne est tridimensionnelle. La formule générale de la sortie de l'antenne est la suivante : Le deuxième terme de l'exponentielle sous l'intégrale de 1'expression (23) caractérise la zone de Fresnel (approximation du second ordre en x/p0). Tout comme le rayonnement d'un piston oscillant, si la fonction d'excitation (la pression acoustique) reste constante le long de l'antenne, pour u=0, la contribution de l'intégrale de Fresnel peut s'exprimer à l'aide de la spirale de Cornu voir l'ouvrage de

M. Bruneau"Introduction aux théories de l'acoustique", éditions Université du Maine, Le Mans, 1983 p. 301.

Si la distance po est suffisamment grande au regard de l'ouverture L de l'antenne, 1'expression (23) montre que le terme quadratique de la phase a peu d'influence dans l'intégrale. La distance po minimale à partir de laquelle il est envisageable de considérer uniquement le cas champ lointain peut d'ailleurs tre obtenue au moyen de cette spirale de Cornu lorsque u=0, voir l'ouvrage de B. Steinberg,"Principles of aperture and array system design", Wiley-interscience publications, John Wiley & Sons, New York, 1976. Ce dernier a montré que pour p0 = L2/2R, l'ignorance du terme de second ordre implique une erreur inférieur à 1 dB. La valeur arbitrairement retenue pour séparer champ lointain de champ proche est alors donnée par : Po > (25) Il est important de noter que la limite entre champ lointain et champ proche est directement liée à la taille de l'antenne et évolue comme le carré de sa longueur.

A titre d'exemple, pour une antenne de longueur L=40 cm et po 3m, la condition (25) implique f<6300 Hz. En conséquence, compte tenu de la bande de fréquences de la parole retenue pour la mise en oeuvre du procédé selon l'invention, un locuteur situé à une telle distance est considéré comme étant en champ lointain au regard du comportement directionnel d'une antenne de 40 cm. Le diagramme de directivité d'une telle antenne peut donc raisonnablement tre obtenu par l'intégrale de Fourier spatiale.

Sous l'hypothèse (25), la fonction de directivité complexe Q (u/2) de l'antenne en champ lointain obtenue au moyen de l'équation normalisée devient pour le cas d'une excitation homogène I le long de l'antenne Antennes superdirectives Une antenne superdirective est une antenne pour laquelle il est effectué une remise en phase des capteurs supérieure à celle correspondant au retard lié au trajet acoustique intercapteurs. Une telle procédure revient à pointer l'antenne vers une région invisible puisque le maximum du lobe principal se trouve désormais en dehors de l'intervalle-1 < u < 1.

Pour une antenne"endfire"à espacement constant entre les capteurs, le retard appliqué sur chaque capteur est de la forme : où 77 est une constante positive déterminant la quantité d' « oversteering ». Cette procédure revient à appliquer des retards correspondant à une vitesse de propagationc/(1+17) inférieure à la vitesse du son c. Il apparaît clairement qu'un tel pointage permet de rendre le lobe principal plus étroit au détriment de lobes secondaires de niveau relatif plus élevé. En outre, l'oversteering peut également faire apparaître des lobes d'aliasing à des fréquences plus basses conduisant à une condition sur d plus stricte : <BR> <BR> d<q#min<BR> <BR> <BR> 2+ruz Puisqu'une partie de l'énergie utile se trouve dans la partie invisible, la sensibilité est réduite et la réduction du bruit incohérent s'en trouve détériorée. Enfin, en pratique, les performances théoriques optimales sont difficilement réalisables du fait de la grande sensibilité d'une telle synthèse aux petites imperfections du système (disparités entre capteurs, etc...). Les erreurs sur l'amplitude ou la phase des signaux sont généralement non corrélées d'un capteur à l'autre et affectent la formation de voie de manière similaire à un bruit incohérent. L'estimation de la réduction de bruit incohérent constitue donc un bon indicateur de degré de faisabilité d'une telle synthèse.