FR 3006802 A1 20141212 FR 1301306 A 20130607 PROCEDE ET SYSTÈME DE DIFFUSION DIRECTIVE DE SIGNAUX SONORES L'invention porte sur un procédé de diffusion directive de signaux sonores, et notamment de messages vocaux. L'invention porte 5 également sur un système pour la mise en oeuvre d'un tel procédé. Il existe des situations où il est souhaitable de diffuser une pluralité de signaux sonores dans des directions distinctes avec une faible interférence mutuelle. Par exemple, dans une gare, il peut être nécessaire d'adresser des annonces différentes aux passagers en attente sur les 10 différents quais. Or, les ondes acoustiques audibles sont très peu directives. L'article de F. Joseph Pompei « Spotlight of sound », Sound&Communication Vol. 55, No. 7 (juillet 2009) décrit un système émettant de façon très directionnelle des ultrasons modulés à fréquence acoustique, la démodulation étant assurée par les propriétés non-linéaires du 15 milieu de propagation (l'air). En effet, les ultrasons sont beaucoup plus directifs que les sons audibles. Cette solution est cependant complexe et couteuse à mettre en oeuvre, car elle nécessite l'utilisation de transducteurs ultrasonores à haute puissance, qui ne sont pas des équipements conventionnels. En outre, la diffusions d'ultrasons très intenses est 20 susceptible de déranger les animaux domestiques et a été soupçonnée d'être nuisible à la santé humaine. Le document EP 1 224 037 divulgue l'utilisation d'une « barre sonore » constituée d'un nombre important de transducteurs acoustiques (haut-parleurs), bénéficiant d'une directivité proportionnelle à la largeur de la 25 barre, permettant d'orienter la direction d'émission sonore par application d'une loi de retard linéaire sur les différents haut-parleurs. Plusieurs signaux sonores peuvent être diffusés dans plusieurs directions par simple superposition des entrées des transducteurs. Cette solution est principalement adaptée à la diffusion de musique. Elle nécessite un nombre 30 élevé de haut-parleurs, bien supérieur à celui des signaux sonores distincts à émettre, ce qui induit un coût élevé et un encombrement important. L'invention vise à surmonter les inconvénients précités de l'art antérieur et à procurer un procédé et un système simples et économiques pour la diffusion directive de signaux sonores. Le procédé et le système de 35 l'invention conviennent tout spécialement à la diffusion de messages vocaux, mais ils peuvent convenir également à d'autres applications, telles que la diffusion de musique. Comme dans le cas du document EP 1 224 037, l'invention utilise une pluralité de transducteurs acoustiques, par exemple regroupés au sein d'une « barre sonore ». Elle s'en différencie par un filtrage préalable particulier des signaux d'entrée réalisant, de façon exacte ou 5 approchée, une inversion spatio-temporelle - ou, de manière équivalente, spatio-fréquentielle - du système d'équations exprimant les phénomènes d'interférence entre les sons émis par les différents transducteurs dans différentes directions de propagation. Ce filtrage permet d'obtenir une bonne directivité et une bonne discrimination des différents signaux sonores avec un 10 plus petit nombre de transducteurs acoustiques, et donc à un moindre cout. Conformément à l'invention, un tel but est atteint par un procédé de diffusion directive d'une pluralité de signaux sonores au moyen d'une pluralité de transducteurs acoustiques, comprenant les étapes suivantes : 15 a) acquérir ou générer N>1 signaux audio; b) filtrer chaque dit signal audio au moyen d'un banc respectif de M â\J filtres, M étant le nombre desdits transducteurs acoustiques et chaque filtre de chaque dit banc étant associé à un transducteur acoustique respectif ; 20 c) pour chaque dit transducteur acoustique, générer un signal de pilotage respectif en combinant linéairement les signaux de sortie des N filtres qui lui sont associés ; et d) fournir à chaque dit transducteur acoustique le signal de pilotage correspondant. 25 Avantageusement, lesdits filtres peuvent être configurés de telle sorte que les sons émis par lesdits transducteurs acoustiques en réponse auxdits signaux de pilotage interfèrent pour générer des signaux sonores correspondant auxdits signaux audio et se propageant dans des directions respectives. 30 Plus particulièrement, lesdits N signaux audio peuvent comprendre : - NS signaux représentatifs de signaux sonores à diffuser selon des directions de propagation distinctes ; et - No signaux identiquement nuls correspondant à des zones de silence intercalées entre lesdites directions de propagation de signaux sonores ou bordant lesdites directions ; où Ns>0, No>0 et Ns+No=N. Encore plus particulièrement, on peut avoir No=Ns-1 et donc M=N=2Ns-1, ce qui signifie que le nombre de transducteurs acoustiques peut être égal à deux fois le nombre de signaux sonores à diffuser, moins un. Cela permet d'avoir une zone de silence interposée entre les lobes du diagramme d'émission sonore, minimisant ainsi les interférences entre signaux sonores. En variante, on peut avoir No=Ns-+1 et donc M=N=2Ns+1, ce qui signifie que le nombre de transducteurs acoustiques peut être égal à deux fois le nombre de signaux sonores à diffuser, moins un. Cette deuxième variante permet d'assurer que même les deux lobes extrêmes du diagramme d'émission sonore sont bordés par deux zones de silence (on considère le cas où l'émission est directive seulement dans un plan horizontal, ce qui est le plus fréquent en pratique, même si l'invention peut également s'appliquer à une directivité selon deux directions orthogonales). Selon une troisième variante, No=0 et Ns=M. Dans ce cas, le nombre de transducteur est égal à celui de signaux sonores à diffuser, mais il n'y a pas de zones de silence entre deux directions d'émission distinctes. Le niveau d'interférence entre les signaux sonores est donc plus important que dans les cas considérés plus haut. Selon différentes caractéristiques avantageuses de l'invention, prises isolément ou en combinaison : - Ladite étape a) peut comprendre également une opération de filtrage passe-haut desdits signaux audio permettant de compenser des pôles à basse fréquence des fonctions de transfert des filtres desdits bancs. - Les signaux audio peuvent être des signaux numériques, lesdites opérations b) et c) sont mises en oeuvre sous forme numérique et 30 ladite étape d) comprend une opération de conversion desdits signaux de pilotage en forme analogique. - Selon un premier mode de réalisation, ladite étape b) peut être mise en oeuvre dans le domaine temporel au moyen de filtres numériques à réponse impulsionnelle finie. Dans ce cas, le procédé peut comprendre également une étape préliminaire de détermination des coefficients desdits filtres numériques à réponse impulsionnelle finie comportant les opérations suivantes : I. acquérir ou générer M signaux audio dits d'apprentissage ; Il. calculer, au moyen d'un modèle mathématique approprié, 10 les N signaux sonores qui se propageraient dans lesdites directions de propagations distinctes si lesdits M signaux audio d'apprentissage étaient appliqués aux entrées desdits transducteurs acoustiques ; et III. déterminer des coefficients desdits filtres tels que, en prenant comme entrées desdits filtres les N signaux sonores ainsi déterminés 15 et en combinant linéairement leurs sorties, on obtient des signaux de pilotage desdits transducteurs acoustiques approchant lesdits M signaux audio d'apprentissage. Ladite opération III peut être mise en oeuvre par une méthode d'optimisation au sens des moindres carrés. Plus particulièrement, elle peut 20 être mise en oeuvre par une méthode d'optimisation au sens des moindres carrés avec un paramètre de régularisation dont la valeur est ajustée de manière itérative de façon à assurer que les amplitudes de crête des signaux de sortie desdits filtres numériques à réponse impulsionnelle finie n'excédent pas des seuils prédéfinis. 25 Lesdits signaux audio d'apprentissage peuvent être des signaux aléatoires présentant une distribution de probabilité gaussienne. En variante, ils peuvent être les signaux sonores à diffuser eux-mêmes, ou des portions desdits signaux. - Selon un deuxième mode de réalisation, ladite étape b) 30 peut être mise en oeuvre au moyen de filtres numériques opérant dans un domaine fréquentiel. Dans ce cas, le procédé peut comprendre également une étape préliminaire de détermination des fonctions de transfert desdits filtres numériques opérant dans un domaine fréquentiel comportant les opérations suivantes : i. choisir N direction de propagation distinctes pour N 5 signaux sonores à diffuser ; ii. déterminer un système d'équations exprimant les signaux sonores se propageant dans lesdites directions de propagations distinctes sous la forme de combinaisons linéaires respectives des sons émis par lesdits transducteurs acoustiques, avec les retards relatifs de propagation 10 correspondants ; iii. convertir ledit système dans le domaine fréquentiel ; et iv. résoudre, de façon exacte ou approchée, le système converti dans le domaine fréquentiel. - Lesdits signaux sonores à diffuser peuvent être des 15 messages vocaux. Un autre objet de l'invention est un système de diffusion directive d'une pluralité de signaux sonores comprenant : - M>1 transducteurs acoustiques ; et - un dispositif de traitement de signaux pour recevoir 20 en entrée ou générer N signaux audio, avec 1 N (il y a plus de transducteurs que de signaux acoustiques indépendants), la matrice F n'est pas carrée et ne peut donc pas être inversée. Dans ces conditions, le système (3) ne peut être résolu que de façon approchée, par exemple au sens des moindres carrés, en calculant une pseudo-inverse de F. La figure 2 illustre un schéma conceptuel de mise en oeuvre d'un procédé selon un mode de réalisation de l'invention basé sur la résolution du système (1) dans le domaine fréquentiel, conformément à l'équation (4) ci-dessus. On prend en entrée les signaux - considérés ici analogiques - yi(t), y2(t) et y3(t) qui correspondent aux signaux sonores à diffuser. Ces signaux passent d'abord par un filtre analogique passe-bande FPB, qui a une double fonction : - un filtrage passe-haut, qui élimine les basses fréquences pouvant donner lieu à des divergences en raison des pôles précités à f=0 ; et un filtrage passe-bas anti-repliement du spectre. Les signaux ainsi filtrés sont ensuite échantillonnés et convertis au format numérique par un échantillonneur-bloqueur-convertisseur EAN. En variante, les signaux pourraient être directement acquis, 5 voire générés, par exemple au moyen d'un synthétiseur vocal, au format numérique, auquel cas le bloc EAN n'est pas nécessaire et le filtrage passe-bas peut être réalisé numériquement. Même dans ce cas, il peut être utile de ré-échantillonner les signaux numériques pour alléger le traitement qui suit en réduisant le nombre d'échantillons qui constituent le message, au moins si 10 une haute fidélité n'est pas exigée (comme dans le cas des messages vocaux, par exemple). Les signaux filtrés et numérisés sont convertis dans le domaine des fréquences par un bloc TFD réalisant une transformée de Fourier discrète, par exemple au moyen de l'algorithme bien connu de 15 transformé de Fourier rapide (d'autres transformations, telle la transformée en cosinus discrets, peuvent être utilisées). Les signaux exprimés dans le domaine fréquentiel, Y1(f), Y2(f) et Y3(f), passent ensuite par des bancs de filtres respectifs : BF1 pour le signal Y1(f), BF2 pour le signal Y2(f) et BF3 pour le signal Y3(f). Le banc de filtres BF1 correspond à la première ligne de la 20 matrice F-1 et comprend trois filtres, dont chacun est associé à une ligne de cette matrice : - filtre FF11, correspondant à F-1(1,1) : fonction de transfert z-2 FF11(z) = (1-z_2)(1-z_1) - filtre FFR12, correspondant à F-1(1,2) : fonction de 25 transfert FF12 (z) = -z- 1(1+z- 1) (1-z-2)(1-z-1) - filtre I1R13, correspondant à F-1(1,3) : fonction de transfert -1 FF13(z) = (1-z-2)(1-z-1) ' et ainsi de suite pour les autres bancs de filtres. Les sorties des « premiers » filtres de chaque banc (c'est-à-dire des filtres FF11, FF21, FF31), dont les fonctions de transfert correspondent aux éléments de la première colonne de la matrice F-1, sont sommées entre elles (plus généralement, sont combinées linéairement). Le 5 résultat de cette somme est transformé à nouveau dans le domaine du temps (bloc de transformée de Fourier discrète inverse TFD-1) et converti au format analogique (convertisseur NA) pour donner le signal de pilotage du premier haut-parleur HP1, xi(t). Il en va de même pour les sortes des « deuxièmes » filtres FF12, FF22, FF32, qui permettent d'obtenir le signal de pilotage du 10 deuxième haut-parleur HP2, x2(t), et pour celles des « troisièmes » filtres FF13, FF23, FF33, qui permettent d'obtenir le signal de pilotage du troisième haut-parleur HP3, x3(t). Les sons émis par les haut-parleurs se propagent et interfèrent comme illustré sur la figure 1, de telle sorte qu'après une certaine 15 distance de propagation (conditions de « champ lointain ») on obtient un premier signal sonore reproduisant le signal audio yi(t), un deuxième signal sonore reproduisant le signal audio y2(t) et un troisième signal sonore reproduisant le signal audio y3(t). Selon un mode de réalisation particulier, le signal audio y2(t) 20 peut être identiquement nul pour générer une zone de silence séparant les « vrais » signaux yi(t) et y3(t) de façon à empêcher toute interférence ; dans ce cas, le deuxième banc de filtres BF2 peut être omis, tout comme les blocs de filtrage, conversion analogique-numérique et transformée de Fourier situés en amont. Réciproquement, il est possible d'utiliser des signaux audio y1(t) et 25 y3(t) identiquement nuls pour délimiter un « vrai » signal sonore y2(t) se propageant de manière directionnelle. La généralisation à des cas prévoyant un nombre plus élevé de signaux audio, signaux acoustiques et/ou haut-parleurs ne pose pas de difficultés particulières. Seule contrainte, le nombre de signaux sonores indépendants à générer (y compris les éventuelles zones 30 de silence) ne peut pas excéder le nombre de haut-parleurs. Les directions de propagation des signaux sonores peuvent être modifiées simplement en agissant sur les fonctions de transfert des filtres numériques. Une limite de ce mode de réalisation réside dans le fait que 5 les signaux audio comprennent en général un grand nombre d'échantillons (80.000 pour une annonce de seulement 10s, avec une bande passante limitée à 4 kHz et donc une fréquence d'échantillonnage de 8 kHz), ce qui nécessite un traitement multi-blocs. Un tel traitement est également nécessaire en présence de contraintes de latence : dans ce cas, on ne peut 10 généralement pas attendre la fin du message pour calculer sa transformée de Fourier, il est nécessaire de le découper en blocs à traiter au fil de l'eau. Cela ajoute de la complexité au traitement, et pose des problèmes de qualité au niveau des transitions entre blocs. Une autre limite est constituée par les problèmes de stabilité 15 précités, qui ne sont pas toujours parfaitement résolus par le filtrage passe-haut (passe-bande) en amont. Pour surmonter ces inconvénients, un mode de réalisation préféré de l'invention utilise des filtres de type à réponse impulsionnelle finie, opérant dans le domaine temporel. Cela est illustré sur la figure 3A, 20 généralement semblable à la figure 2 mais où on remarque l'absence des blocs de transformée de Fourier/transformée de Fourier inverse. La longueur de ces filtres (c'est-à-dire le nombre de coefficients constituant leur réponse impulsionnelle) doit être suffisante pour permettre une approximation acceptable des fonctions de transfert théoriques 25 (à réponse impulsionnelle infinie), tout en minimisant la complexité computationnelle. Pour des signaux vocaux, ayant une bande passante comprise entre quelques centaines de Hz et quelques kHz, une longueur comprise entre 100 et 1000 et de préférence entre 100 et 500 est généralement adaptée. 30 La figure 3B illustre une méthode de détermination des coefficients des filtres à réponse impulsionnelle finie permettant de trouver un compromis entre deux exigences pouvant être contradictoire : celle d'approcher aussi bien que possible les fonctions de transfert théoriques données par la matrice F-1 et celle, tout aussi importante, d'éviter des amplitudes transitoires trop fortes, qui pourraient saturer ou même endommager les transducteurs acoustiques (stabilité). Cette méthode se base sur une technique d'apprentissage reposant sur une optimisation au sens des moindres carrés. Avantageusement, l'apprentissage est réalisé à partir de signaux aléatoires ou pseudo-aléatoires présentant des caractéristiques (bande spectrale, distribution de probabilité) représentatives des signaux sonores à diffuser ; il peut s'agir par exemple des échantillons réalisés suivant un tirage aléatoire gaussien ; en variante, il peut s'agir des signaux audio « réels » à diffuser, ou de fragments de ces signaux. L'expérience a montré que, pour des filtres ayant une réponse impulsionnelle de CF=129 coefficients, des signaux d'apprentissage de NE=2000 échantillons sont suffisants. La première étape (I) du procédé de la figure 3B est la génération ou l'acquisition des signaux d'apprentissage xi'(t), x2'(t), x3'(t) qui, comme expliqué plus haut, peuvent être aléatoires ou pseudo-aléatoires. Ensuite (étape Ila) un modèle mathématique MMP - correspondant par exemple au système (1) - est utilisé pour calculer les signaux sonores yi'(t), y2'(t), y3'(t) qui seraient obtenus si les signaux d'apprentissage étaient utilisés pour piloter les haut-parleurs. Dans un souci de simplification, on considère que la période d'échantillonnage des signaux acoustiques est égale à At, de telle sorte que le système (1) peut se réécrire de la façon suivante : ly'l [i] = xi [i] + x2 [i + 1] + x3 [i + 2] y'2 [i] = xi [i + 1] + x2 [i + 1] + x3 [i + 1] (5a) y'3 [i] = xi [i + 2] + x2 [i + 1] + x3[i] où i est un indice entier qui prend des valeurs comprises entre 1 et NE-2. Ainsi, les vecteurs y'1, y'2 et y'3 sont des vecteurs colonnes à NE-2 éléments. 3006 802 13 L'étape Ilb consiste à représenter le filtrage FIR de la figure 3a (qui, sur le plan mathématique, est constitué par un ensemble de produits de convolution) sous forme matricielle. Pour ce faire, il est nécessaire d'arranger les éléments des vecteurs y'i, y'2 et y3 dans une matrice Y donnée 5 par : Yi(1) Yi (2) Yi(NE - CF - 1)\ Yi(CF) Y1 (CF + 1) Yi(NE - 2) Y2(1) Y2(2) Y2 (NE - CF - 1) Y= (5b) Y2(CF) Y2 (CF + 1) y2 (NE - 2) Y3(1) Y3 (2) y3 (NE - CF - 1) \y3 (CF) y3 (CF + 1) y3 (NE. - 2) Comme on peut le voir chaque colonne de la matrice Y est formée par une concaténation de trois vecteurs constitués de CF éléments 10 successifs de y'i, de y'2 et de y3. L'indice du premier élément de chaque vecteur est décalé d'une unité à chaque changement de colonne. Cette matrice a donc 3.CF=387 lignes et NE-CF-1=1870 colonnes. Dans l'étape Ilb on définit également trois vecteurs d'entrée (vecteurs ligne) de dimension NE-CF-1=1870, construits à partir des signaux 15 d'entrée x2 et x3 en excluant les premier (CF+1)12 et les derniers [(CF+1)12+1] éléments : Xi=[xi((CF+1)/2), , xi(N-(CF+1)/2-1)] X21x2((CF+1)/2), , x2(N-(CF+1)/2-1)] (5c) X3=k3((CF+1)/2), , x3(N-(CF+1)/2-1)] 20 Le filtrage de la figure 3A peut alors s'écrire, sous forme matricielle : X1(t)=FIR1 Y (6a) X2(t)=FIR2 Y (6b) X3(t)=FIR3 Y (6c) 25 où FIR1, vecteur ligne de 3-CF éléments, est la concaténation selon un vecteur ligne unique des trois vecteurs de CF éléments chacun contenant les coefficients de filtres FIR11, FIR21, FIR31 et ainsi de suite. Dans la suite on s'intéressera à la seule détermination de FIR1, car le problème se pose exactement dans les mêmes termes pour FIR2 et FIR3. L'équation (6a), dont l'inconnue est le vecteur FIR11, est 5 surdéterminée, car il y a plus d'équations (NE-CF-1=1870) que d'inconnues (CF=129) ; elle n'admet donc pas de solution exacte. Cependant, il est possible de trouver une solution au sens des moindres carrés, et plus particulièrement des moindres carrés régularisés. Cette solution est donnée par : 10 FIR1=X1(t) Yi (7a) avec Yi=y1"(yyT+ ki 1)1 (7b) où « T » est l'opérateur de transposition, I la matrice unité et k1 un coefficient de régularisation. Comme indiqué sur la figure 3B, l'étape Illa comprend le 15 calcul de la matrice Yi avec une valeur prédéterminée du paramètre k1 par application de l'équation (7b) et l'étape Illb comprend le calcul du vecteur FIR1 par application de l'équation (7a). L'étape Illc comprend la vérification d'une condition de régularisation ; cette vérification peut se faire en appliquant les filtres FIR11, FIR12 et FIR13 (dont les coefficients forment le vecteur 20 FIR1) à au moins un signal de contrôle, qui peut être le signal d'apprentissage xi'(t) lui-même, et en vérifiant que les amplitudes de crête des signaux de sortie de ces filtres n'excèdent pas un seuil prédéfini. Si la condition n'est pas satisfaite, le paramètre k1 est incrémenté d'une valeur positive (étape Ille) et les étapes Illa, Illb et Illc sont répétées. Si la condition est satisfaite, par 25 contre, les coefficients des filtres FIR11, FIR12 et FIR13 sont mémorisés. Ces opérations sont répétées pour les autres bancs de filtres. Ensuite, il est possible de tester la fidélité obtenue en appliquant les signaux yi'(t), y2'(t), y3'(t), calculés à l'étape Ilb, à l'entrée de ces filtres et en comparant les signaux de sortie aux signaux d'apprentissage 30 xi'(t), x2'(t), x3'(t). Si l'écart (par exemple, l'erreur quadratique) entre ces signaux dépasse un seuil prédéterminé, cela signifie que le procédé de la figure 3B doit être itéré en augmentant la longueur des filtres. La figure 4 représente, très schématiquement, un système selon un mode de réalisation particulier de l'invention. Ce système comprend 5 essentiellement un dispositif de traitement de signaux DTS et une barre sonore BS. Le dispositif DTS comprend un processeur généraliste ou spécialisé, une mémoire stockant un programme pour la mise en oeuvre d'un filtrage selon l'invention, des ports d'entrée et de sortie. Il peut également comprendre des convertisseurs analogique/numérique en entrée et/ou des 10 convertisseurs numérique/analogique en sortie. Dans les cas les plus simples, il peut être basé sur un ordinateur de bureau ou portable programmé de manière opportune, le cas échéant relié à des circuits d'acquisition audio (en entrée) et/ou d'interface avec des chaines audio. Dans le mode de réalisation illustré, le dispositif DTS présente 15 deux entrées pour des signaux audio numériques yi(t) et y2(t) respectifs, représentatifs de messages vocaux (annonces) devant être diffusés dans des directions distinctes. Le dispositif génère en interne un troisième « signal » identiquement nul destiné à être « diffusé » entre les deux messages vocaux pour créer une zone tampon de silence évitant toute interférence entre ces 20 derniers. Le dispositif DTS présente en outre trois sorties pour des signaux analogiques xi(t), x2(t), x3(t) qui sont fournis en entrée à la barre sonore BS, le cas échéant par l'intermédiaire de chaines d'amplification de puissance (non représentées). La barre sonore BS comprend trois haut-parleurs alignés, 25 montés sur un support commun. Il est important que les haut-parleurs présentent une réponse aussi linéaire que possible et que leurs couplages mutuels soient minimisés (utilisation d'un support rigide et de préférence de sous-caissons individuels). De préférence, les haut-parleurs devront avoir des propriétés sensiblement identiques, ou alors un procédé d'égalisation devra 30 être mis en oeuvre, le cas échéant en ajustant de manière opportune les fonctions de transfert des filtres numériques mis en oeuvre par le dispositif DTS. Dans les exemples de réalisation considérés ici on a considéré, dans un souci de simplicité, que les directions angulaires de 5 propagation des signaux sonores gisent dans un même plan. Il ne s'agit pas là d'une limitation essentielle, et l'invention peut être aisément généralisée au cas de directions de propagation orientables en deux dimensions angulaires (par exemple site et gisement). Dans ce cas, les transducteurs acoustiques ne doivent pas être alignés, mais être agencés selon un motif bi- ou 10 tridimensionnel. REVENDICATIONS1. Procédé de diffusion directive d'une pluralité de signaux sonores au moyen d'une pluralité de transducteurs acoustiques, comprenant les étapes suivantes : a) acquérir ou générer N>1 signaux audio (yi, y21 y3); b) filtrer chaque dit signal audio au moyen d'un banc respectif (BF1, BF2, BF3) de M â\J filtres (FIR11...), M étant le nombre desdits transducteurs acoustiques et chaque filtre de chaque dit banc étant associé à un transducteur acoustique respectif ; c) pour chaque dit transducteur acoustique (HP1, HP2, HP3), générer un signal de pilotage respectif (xi, x2, x3) en combinant linéairement les signaux de sortie des N filtres qui lui sont associés ; et d) fournir à chaque dit transducteur acoustique le signal de pilotage correspondant. 2. Procédé selon la revendication 1 dans lequel lesdits filtres sont configurés de telle sorte que les sons émis par lesdits transducteurs acoustiques en réponse auxdits signaux de pilotage interfèrent pour générer des signaux sonores correspondant auxdits signaux audio et se propageant dans des directions respectives. 3. Procédé selon l'une des revendications précédentes dans lequel lesdits N signaux audio comprennent : - NS signaux représentatifs de signaux sonores à diffuser 25 selon des directions de propagation distinctes ; et - No signaux identiquement nuls correspondant à des zones de silence intercalées entre lesdites directions de propagation de signaux sonores ou bordant lesdites directions ; où Ns>0, No>0 et Ns+No=N. 30 4. Procédé selon la revendication 3 dans lequel la valeur No et choisie parmi Ne-1 et Ne+1. 5. Procédé selon l'une des revendications précédentes dans 5 lequel ladite étape a) comprend également une opération de filtrage passe-haut desdits signaux audio permettant de compenser des pôles à basse fréquence des fonctions de transfert des filtres desdits bancs. 6. Procédé selon l'une des revendications précédentes dans 10 lequel les signaux audio sont des signaux numériques, lesdites opérations b) et c) sont mises en oeuvre sous forme numérique et ladite étape d) comprend une opération de conversion desdits signaux de pilotage en forme analogique. 7. Procédé selon la revendication 6 dans lequel ladite étape 15 b) est mise en oeuvre dans un domaine temporel au moyen de filtres numériques à réponse impulsionnelle finie (FIR11...). 8. Procédé selon la revendication 7 comprenant également une étape préliminaire de détermination des coefficients desdits filtres 20 numériques à réponse impulsionnelle finie comportant les opérations suivantes : I. acquérir ou générer M signaux audio dits d'apprentissage (x'i, x'2, x'3) II. calculer, au moyen d'un modèle mathématique approprié 25 (MMP), les N signaux sonores (yi, Y2, y3) qui se propageraient dans lesdites directions de propagations distinctes si lesdits M signaux audio d'apprentissage étaient appliqués aux entrées desdits transducteurs acoustiques ; et III. déterminer des coefficients desdits filtres tels que, en 30 prenant comme entrées desdits filtres les N signaux sonores ainsi déterminés et en combinant linéairement leurs sorties, on obtient des signaux de pilotagedesdits transducteurs acoustiques approchant lesdits M signaux audio d'apprentissage. 9. Procédé selon la revendication 8 dans lequel ladite 5 opération III est mise en oeuvre par une méthode d'optimisation au sens des moindres carrés. 10. Procédé selon la revendication 9 dans lequel ladite opération III est mise en oeuvre par une méthode d'optimisation au sens des 10 moindres carrés avec un paramètre de régularisation (k1) dont la valeur est ajustée de manière itérative de façon à assurer que les amplitudes de crête des signaux de sortie desdits filtres numériques à réponse impulsionnelle finie n'excédent pas des seuils prédéfinis. 15 11. Procédé selon l'une des revendications 8 à 10 dans lequel lesdits signaux audio d'apprentissage sont des signaux aléatoires présentant une distribution de probabilité gaussienne. 12. Procédé selon l'une des revendications 8 à 10 dans lequel 20 lesdits signaux audio d'apprentissage sont les signaux sonores à diffuser, ou des portions desdits signaux. 13. Procédé selon la revendication 6 dans lequel ladite étape b) est mise en oeuvre dans un domaine fréquentiel au moyen de filtres 25 numériques opérant dans un domaine fréquentiel (FF11...). 14. Procédé selon la revendication 13 comportant également une étape préliminaire de détermination des fonctions de transfert desdits filtres numériques opérant dans un domaine fréquentiel comportant les 30 opérations suivantes :i. choisir N direction de propagation distinctes pour N signaux sonores à diffuser ; ii. déterminer un système d'équations exprimant les signaux sonores se propageant dans lesdites directions de propagations distinctes sous la forme de combinaisons linéaires respectives des sons émis par lesdits transducteurs acoustiques, avec les retards relatifs de propagation correspondants ; iii. convertir ledit système dans le domaine fréquentiel ; et iv. résoudre, de façon exacte ou approchée, le système 10 converti dans le domaine fréquentiel. 15. Procédé selon l'une des revendications précédentes dans lequel lesdits signaux sonores à diffuser sont des messages vocaux. 15 16. Système de diffusion directive d'une pluralité de signaux sonores comprenant : M>1 transducteurs acoustiques (HP1, HP2, HP3) ; et un dispositif de traitement de signaux (DTS) pour 20 recevoir en entrée ou générer N signaux audio (yi, y2), avec 1 17. Système selon la revendication 16 dans lequel lesdits 30 filtres sont configurés de telle sorte que les sons émis par lesdits transducteurs acoustiques en réponse auxdits signaux de pilotage interfèrentpour générer des signaux sonores correspondant auxdits signaux audio et se propageant dans des directions respectives.