Brevet FR2828327A1 — publié le 20030207

O est un vecteur nul. - 1 2828327 L' invention se rapporte au domaine de l'acoustique aérienne appliquée à des systèmes de sonorisation pour des salles en particulier. La présente invention se rapporte à un procédé de réduction d'écho mis en _uvre 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 _uvre 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. 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 1 5 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 même 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. I1 n'en est pas de même 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 2 2828327 général supérieure de 6 à 20 dB à la parole locale captée par les microphones, entrâîne l'effet Larsen si aucun traitement n'est appliqué. Si les techniques connues d'annulation d'écho ne permettent pas de traiter de manière s atisfaisante cet effet néfaste, une variation de gain importante do it é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 _uvre 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 l'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, même 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 perturtations 3 2828327 ne peuvent pas étre modélisées par un pro duit de convolution entre le signal provenant du haut parleur et une réponse impulsionnelle de filtrage. Qui plus est. le fonctionnement des transbucteurs 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 équi 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, à l'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 _uvre 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 4 2828327 microphone et le haut parleur. Un tel calcul est généralement suffis ant dans la plup art des situations puisque les microphones sont dans le champ direct du haut parleur. Cependant, cette solution nocessite 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 gênante pour un système individuel pour lequel cette disposition est imposée, il n'en est pas de même 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 pass ante et à une sensibilité plus grande aux obstacles situés à proximité. Enfin, une amplification du bruit incoLérent est à déplorer et n'est contrôlable qu'en modifiant la distance séparant les deux 1 5 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 (, 0) 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 l'espace qui est habituellement mesurée dans deux plans orthogonaux et représentée en coordonnées polaires (r,, 0) 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 (, 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 2828327 une efficacité certaine, mais en basses fréquences (f de champ lointain n'est pas vérifé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. A cet effet, l'invention a pour objet un procédé de réduction d'écho, mis en _uvre avec un dispo s it if de pris e 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 pro che du facteur de directivité Fd ( f) dont l'expression est la suivante: [id(f)= 1 1 (1) -W (f)D r (f)W(f) avec Dr(.f)=.t J W(f,,H)H(fr,)HH(fr,)sindd (2) o=o 4=o 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: C (f)W(f) = S(f) (3) et les dites pondérations w (f) satisfont une deuxTè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(f)= 0 (4) les formules (1) à (4) étant telles que: 6 2828327 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 et chaque capteur de l'antenne, calculées à la fréquence f, Dr(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) sont des pondérations spatiales permettant une atténuation plus importante des ondes provenant des directions des hautparleurs 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. Ains i, l' invent ion prend en co mpte 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 transtert 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 l'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 hautparleurs 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 hautparleurs et des capteurs. D'autres caractéristiques et avantages de l' invention apparfflront 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 _uvre 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 _uvre 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 _uvre d'un procédé selon l'invention. La figure S est un synoptique d'une mise en _uvre d'un procédé selon 1 5 l' invention. Le procédé de réduction d'écho selon l' invention est mis en _uvre 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 l'expression de la première contrainte (3) linéaire de gain désiré en champ lointain dans les directions correspondant à des sources utiles (lo cuteurs), déhinies par l' équat ion matricielle: C (f)W(f) = s(f) (3) é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. 8 2828327 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 transtert 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(f)W(f)= 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, O est un vecteur nul. Dans une troisième étape 3, le pro cédé selon l'invention consiste à maximiser le facteur de directivité (1) en champ proche. Le facteur de directivité en champ pro che F est défini par l'équation: d (f) 1 (1) 1 wH (f)D,,(fjW(f) r 27r r2N avec Dr(f)= | | W(f')H(fr,)HH(fr,) sin67dd J o é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 et chaque capteur de l'antenne, calculées à la fréquence f, Dr(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,ç,) 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 _uvre 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. (f) fixe la valeur, suivant l'équation: 9 2828327 W(f) Wtf)= () - 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 larg eur suffis ante du lo be principal de l' ant enne p our 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 _uvre 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: w(f)Ho(f)Ho (f)w(f)=go(f) (6) W (f)H (f)H, (f)W(f) gi (f) W (f)Hk (f)Hk (f)W(f) gk (f) système d'équations dans lequel, Hj, i=O,..., k représente le vecteur de propagation entre une source ponctuelle en champ lointain dans la direction angulaire (pj, O;) 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 l'expression de la matrice de directivité suivant la formule: Dr(f)= Dr(f)+ c((f)I + É"i(f)Hi(f)HiH (f) (7) [=0 o,u(f) est un coeffcient de Lagrange qui est ajusté pour satisfaire la contrainte de robustesse (5), et ai(f), i=O,..., 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: 2828327 W(f)=Dr (f)C(f)(cH(f)Dr (f)C(f)) S(f) (8) 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 _uvre d'un pro cé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 l'exemple retenu, par onze microphones 1 lm, 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 exprimoe 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 même 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 1 1 2828327 relation (2) les directions 0 et T, 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 _uvre 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 fgures 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 l'exemple, le générateur 14 génère une séquence pseudo-aléatoire de longueur maximale, dite MLS. Le signal ym(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 hesm(n) de chaque canal acoustique haut parleur / microphone m, m=1, 2, N avec N=11, est déterminée à partir du signal Ym (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 hesm(n) peut être déterminée par intercorrélation en utilisant en particulier une transformée de Hadamard. Un fenêtrage 18m temporel est pratiqué sur la réponse impulsionnelle estimée he5m(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ébuction du couplage s'appuie sur les éléments du couplage les plus stables (trajet acoustique direct et premières réflexions sur les obstacles fxes, par exemple les murs à proximité) et écarte les parties des éléments du couplage susceptibles de se modifer 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 fenêtre 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 fenêtres d'analyse que de contraintes de couplage. La fenêtre d'analyse peut être rectangulaire; cependant d'autres formes de fenêtres telles que Hamming, etc 12 2828327 peuvent être avantageusement utilisées. L'étendue de la fenêtre 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 l9m des zéros à la suite de la fenétre 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 Hes-direcm( ) 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 l'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 _uvre pour des fréquences inférieures ou égales à 1400 Hz. Les analyses fréquentielles ont été effectuées sur une fenêtre 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 R à 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 direct ivité en champ lo intain. 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 13 2828327 compte de deux contraintes sur le couplage au moyen de deux fenêtres d'analyse: la première est centrce 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 l'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 l'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 échogramm es 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 tra^nage en fm 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 même 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 transtert 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. I1 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-1, 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 14 2828327 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). 2828327 ANNEXE 1 Directivité Un capteur est caractérisé par sa directivité dans l'espace qui est habituellement mesurée dans deux plans orthogonaux et représentée en coordonnées polaires (r, p, O) 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 ( 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éfmition, 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: |M(O /)| = A()COSn (9) o A(c) 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: I I [ +] ( 1 O) avec, constante positive appelée constante de cardioicité. Le cas p=lcorrespond à une directivité cardioïde, p=à 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 d ou encore T d 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 même sensibilité et ont un même coefficient de pondération, s'écrit sous la forme suivante: 16 2828327 N-l in (X) = N,(x-nd) ( l l) En procédant à la transforrnée de Fourier spatiale, il est possible d'en déduire le spectre spatial de l'antenne Q(u/) N-l -j21rn_ 2 -j2nd- Q() = J in(X)e dx = e (12) n= Une autre interprétation de cette expression formelle peut être obtenue en considérant l'antenne comme un échantillonneur spatial tronqué par une fenêtre rectangulaire. En conséquence, la fonction d'excitation peut s'écrire sous la forme: 0o in (X) = N rectL(x) ú(x - nd) avec L =Nd et 1 S1 -L/2 () l O si x L / 2 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: Q(u/A) = TF{in (x)} = N TF{rectL(x)} * TF( (x - nd)| ( 1 4) Le premier terme conduit à un sinus cardinal (correspondant au spectre spatial de l'antenne continue): TE{re tL()} sin('zLu/) (15) Quand au second, il a pour expression: T1F{ (x - nd)l ( d) ( 1 6) n=-= n,=-m La transformée de Fourier d'un train d'impulsions équidistantes de d dans l'espace physique conduit également à un train d'impulsions équidistantes de ?/d 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: Q( /) 1 sin(u/) ( m2) 1 [)(À d)] (17) N rd u/À m=)() Le diagramme polaire correspond à la somme de sinus cardinaux décalés d'un intervalle de \/d. Le premier maximum centré en u=0 correspond au lobe principal de l'antenne de largueur \/L. L'allure du diagramme de la fonction TF{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 discrète est très pro che de celui de l'antenne continue. Un calcul direct de la fonction donnée par l'expression (12) normalisée (en remarquant qu'il s'agit d'une suite géométrique) conduit à un déphasage près à l'expression: Q(u/,)= sin (uN /) (18) Cette expression est très proche de celle obtenue pour une antenne continue dans le cas o d 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. Fact eur de direct ivit é 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 même pression efficace et de même fréquence provenant dans un angle 18 2828327 solide de 4. 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 même valeur en tous les points et o l'intensité acoustique est la méme dans toutes les directions. Champ lo intain/champ pro che Soit une antenne linéaire de longueur L, centrée en x = 0, 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 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 lié à la présence de la source S(xo,yo) a pour expression: s(xo,yo) = K | dx k= Àl, = 2 (19) -/2 |RS| 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: 1: 1-1 12 p=|RS|=|RO+OS| =|05|1+1 I: +2 1 12 (20) Dans le cas o la source se trouve suffisamment éloignée de l'origine au regard de la dimension de l'antenne (|OSI--PO>>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: |RS|A P! + RO 0S + | | l (21) 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: 19 2828327 RO.05 I Ro2 kP l 312 +2os2 | ( - kR02: (22) -L|2 -/2 En posant U=COS, il vient l'expression normalisée par rapport à l' axe du diagramme de directivité de l'anterme linéaire: t/2 jk xu- Q(k,u)= | in(x)e 2P Jdx (23) -12 o i', correspond au champ de pression qui satisfait à: IQ,O)I=I. 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 l'expression cidessus est donc l'analogue de la zone de Fraunhofer (approximation du premier ordre en x/po). En outre, la taille limitée de l'antenne peut être interprétée comme l'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. I1 n'est donc pas surprenant d'observer dans ce cas une grande similitude dans les résultats. En 1S l'occurrence, à l'instar de la théorie de la diffraction qui conduit à affrmer que la fgure 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 fgure 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 vo ire triple si l' antenne est tridimens ionnelle. La formule g énérale de la sortie de l'antenne est la suivante: Q(k, C0) = |||in (x,))e ' d Y (24) v Le deuxième terme de l'exponentielle sous l'intégrale de l'expression (23) caractérise la zone de Fresnel (approximation du second ordre en x/pO). 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=O, la contribution de l'intégrale de Fresnel peut s'exprimer à l'aide de la spirale de Cornu voir l'ouvrage de 2828327 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'anterme, l'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 lo intain peut d' ailleurs étre obtenue au moyen de cette spirale de Cornu lorsque u=O, 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 = r2/2i, l'ignorance du terme de second ordre implique une erreur inférieur à 1 dB. La valeur artitrairement retenue pour séparer champ lointain de champ pro che est alors donnée par: po > [2/; (25) I1 est important de noter que la limite entre champ lointain et champ proche est directement lice à 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 Fourier spatiale. Sous l'hypothèse (25), la fonction de directivité complexe Q(u/) 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 /,/2 _j2,U sin(LU/) Q(2)= l Ie dx ru/À (26) 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 21 2828327 Pour une antenne "endfire" à espacement constant entre les capteurs, le retard appliqué sur chaque capteur est de la forme: rn = (1 + '1) nd = ( d + rSup) (27) o est une constante positive déterminant la quantité d'" oversteering ". Cette procédure revient à appliquer des retards correspondant à une vitesse de propagation c/(l+) inférieure à la vitesse du son c. I1 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 apparâître des lobes d'aliasing à des fréquences plus basses conduisant à une condition sur d plus stricte: d 22 2828327 REVENDICATIONS 1. Procédé de réduction d'écho, mis en _uvre avec un dispositif de prise de son multi capteurs formant une antenne (10) et un dispositif (13) de restitution sonore, qui consiste à soumettre les signaux de sortie des capteurs (llm) à des pondérations complexes w(f), les dites pondérations w(f) étant calculées par une maximisation (3) sous contraintes aux basses fréquences et en champ pro che du facteur de directivité Fd ( f) dont l'expression est la suivante: [;d (f)= (1) -w (f)Dr(f)w(f) avec Dr(f)= | | W(f')H(fr,)HH(fr,)sinodd (2) e q=O le calcul étant tel que les dites pondérations w(f) satisfont (1) une première contrainte linéaire sur le module et la phase de la fonction de transtert 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: C (f)w(f) = s(f) (3) ledit procédé étant caractérisé en ce qu'en outre les dites pondérations w(f) sont calculées de façon à satisfaire (2) une deuxième contrainte déterminée à partir de mesures in situ des fonctions de transfert complexes de canaux acoustiques délimités par les entrées des hautparleurs du dispositif de restitution sonore et les sorties des capteurs composant le dispositif de prise de son, la formulation de cette deuxième contrainte à chaque fréquence f étant la suivante: MH(f)w(f)=0 (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 et chaque capteur de l'antenne, calculées à la fréquence f, Dr(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) sont des pondérations spatiales permettant une atténuation plus imp o rt ante des o ndes pro venant des direct io ns d es haut-parleurs p o ur rébuire le couplage direct avec ceux-ci, 23 2828327 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, dits canaux acoustiques de couplage, et, O est un vecteur nul. 2. Procédé de réduction d'écho selon la revendication 1, dans lequel les dites pondérations w(f) satisfont en outre une contrainte (4) non linéaire fixant une limite pour la réduction du bruit incohérent en sortie du dispositif de prise de son, la formulation de cette contrainte non linéaire à chaque fréquence f étant la suivante: W (f)W(f)= R (f) (5) formulation dans laquelle, Rmin(f) est la réduction de bruit incohérent minimale souhaitée, c'est-à-dire le rapport de l'énergie de bruit incohérent en sortie du dispositif de prise de son normalisé par le nombre de capteurs sur l'énergie moyenne de bruit incohérent d'un capteur seul. 3. Procédé de réduction d'écho selon l'une des revendications 1 et 2, dans lequel les dites pondérations w(f) satisfont en outre des contraintes (5) sur l'ouverture du lobe principal et sur la réduction des lobes d'aliasing et des lobes secondaires, ces dites contraintes s'exprimant sous la forme d'un système d'équations non linéaires: wH (f)Ho (f)Ho (f)w(f)=go(f) (6) W (f)H,(f)H,(f)W(f) g,(f) W (f)Hk (f)Hk (f)w(f) gk (f) système d'équations dans lequel, 24 2828327 Hi, i=O,..., k représente le vecteur de propagation entre une source ponctuelle en champ lointain dans la direction angulaire (pi, O;) et les capteurs de l'antenne et, gi le gain désiré dans cette direction. 4. Pro cédé de réduction d'écho selon l'une des revendications 1 à 3, dans lequel une pondération temporelle est appliquée sur au moins une partie des réponses impulsionnelles des canaux acoustiques de couplage avant application de la contrainte (4), pour privilégier l'influence de ces parties correspondant au couplage direct et / ou semi direct dans la détermination des vecteur des pondérations complexes w. 5. Procédé de réduction d'écho selon l'une des revendications 1 à 3, dans lequel les mesures des caractéristiques des canaux acoustiques de couplage sont effectuées au moyen de séquences de longueur maximale appliquées aux entrées de haut parleurs, du dispositif (13) de restitution sonore, situés au voisinage du dispositif (8) de prise de son à la mise en route du dit dispositif (8) de prise de son. 6. Dispositif (6) de réduction d'écho, destiné à fonctionner avec un dispositif (13) de restitution sonore pour la mise en _uvre d'un procédé selon l'une des revendications 1 à 5, le dit dispositif comprenant: - un dispositif (8) de prise de son multi-capteurs formant une antenne (10), - des moyens (21) de calcul d'une première contrainte, la formulation de cette première contrainte à chaque fréquence f étant la suivante: C (f Jw(f) = s(f) (3) - des moyens (14-21) de calcul d'une deuxième contrainte, la formulation de cette deuxième contrainte à chaque fréquence f étant la suivante: M H (f)w(f) = 0(4) - des moyens (22) de calcul de l'équation suivante: w (f) = Dr (f) C (f) (C (f) Dr (f) C (f)) s(f)(8) équation dans laquelle D (f)=Dr(f) et Dr (f) = J. W(f,,)H(f,r,,) HH (f,r,,)sindd(2) 2828327 les formules (2), (3), (4) et (8) étant telles que: 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, dits canaux acoustiques de couplage, et, O est un vecteur nul, 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 O et chaque capteur de l'antenne, calculées à la fréquence f, Dr(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) sont des pondérations spatiales permettant une atténnation plus importante des ondes provenant des directions des haut-parleurs pour réduire le couplage direct avec ceux-ci. 7. Dispositif (6) de réduction d'écho selon la revendication 6 comprenant en outre des moyens de calcul d'une contrainte non linéaire fixant une limite pour la réduction du bruit incohérent en sortie du dispositif de prise de son, la formulation de cette contrainte à chaque fréquence f étant la suivante: W (f)W(f)=R (f) (5) formulation dans laquelle, Rmin(f) est la réduction de bruit incohérent minimale souhaitée, c'est-à-dire le rapport de l'énergie de bruit incohérent en sortie du dispositif de prise de son normalisé par le nombre de capteurs sur l'énergie moyenne de bruit incohérent d'un capteur seul. 26 2828327 8. Dispositif (6) de réduction d'écho selon la revendication 6 comprenant en outre des mo yens de calcul de contraintes sur l'ouverture du lo be principal et sur la réduction des lobes d'aliasing et des lobes secondaires, ces dites contraintes s'exprimant sous la forme d'un système d'équations non linéaires: w (f)Ho (f)Ho (f)w(f) go (f) (6) W (f)H, (f)H (f)W(f) g (f) w (f)Hk(f)Hk(f)w(f) gk(f) système d'équations dans lequel, Hj, i=O,..., k représente le vecteur de propagation entre une source ponctuelle en champ lointain dans la direction angulaire (