Brevet FR2955442A1 — publié le 20110722

L'invention concerne un procédé de détermination de filtrage, ainsi qu'un programme d'ordinateur et un dispositif associés. Le champ sonore généré dans une salle en un point d'écoute résulte non seulement des sons directement émis par les haut-parleurs disposés dans la salle, mais aussi des sons réfléchis par les parois de la salle et le son perçu par un auditeur situé au point d'écoute peut de ce fait être très différent du son fourni par la source, tel que par exemple un lecteur CD ou DVD. Il a de ce fait été proposé des solutions de compensation des effets de la salle (en anglais "room compensation") qui visent à rapprocher autant que possible le son perçu par l'auditeur de celui reçu en entrée par le système de reproduction sonore. Une telle compensation est par exemple réalisée en jouant sur les filtres appliqués aux différents haut-parleurs. Le brevet US 6 195 435 propose dans ce cadre de configurer un ensemble de filtres associés chacun à un haut-parleur dans ce but. On connaît par ailleurs les documents US 2007/025 559 et US 2008/279 318 dans lesquels on définit séparément d'une part des filtrages à appliquer aux basses fréquences et d'autre part des filtrages à appliquer aux hautes fréquences. Dans ce contexte, l'invention propose un procédé de détermination d'un filtrage à appliquer à une pluralité de haut-parleurs, caractérisé par les étapes suivantes : a) détermination des réponses fréquentielles respectivement générées par les haut-parleurs de la pluralité en au moins un point ; b) détermination d'un profil d'énergie moyenne sur la base desdites réponses fréquentielles ; c) détermination du filtrage à appliquer sur la base du rapport entre un profil cible et le profil d'énergie moyenne. Du fait de sa faible variabilité dans l'espace, l'énergie moyenne associée aux réponses fréquentielles se révèle en effet une grandeur particulièrement appropriée pour déterminer un filtrage qui donne de bons résultats dans l'ensemble de la salle concernée (et non seulement aux points d'écoute considérés). Cet avantage se ressent en particulier à haute fréquence où d'autres grandeurs subissent des variations importantes dans l'espace d'écoute. Selon l'exemple de réalisation proposé et décrit plus loin, l'étape b) comprend la détermination d'un profil d'énergie moyenne associé audit point, incluant la détermination de la somme d'une partie au moins des modules des réponses fréquentielles générées audit point. Le profil d'énergie moyenne associé audit point est alors par exemple déterminé en divisant ladite somme par la racine carrée du nombre des modules supérieurs à une valeur de seuil afin de ne tenir compte que des sources prépondérantes, de même ordre de grandeur. Comme démontré plus loin dans la description, ces étapes permettent d'évaluer de manière simple le profil d'énergie moyenne. On peut prévoir par ailleurs que l'étape b) comprenne la détermination d'une pluralité de profils d'énergie moyenne locaux associés chacun à un point et que le profil d'énergie moyenne soit alors déterminé comme la moyenne des profils locaux. Le filtre déterminé à l'étape c) est par exemple en pratique un filtre à réponse impulsionnelle finie. L'invention propose également un procédé de détermination d'un ensemble de filtrages associés respectivement à une pluralité de haut-parleurs, caractérisé par les étapes suivantes : - détermination d'un filtrage haute-fréquence selon un procédé conforme à ce qui a été présenté ci-dessus ; - détermination d'une pluralité de filtrages basse-fréquence associés chacun à un haut-parleur et tels que la réponse fréquentielle générée en au moins un point par les haut-parleurs avec application du filtrage associé corresponde à une réponse cible sur une plage de fréquence déterminée ; - détermination du filtrage associé à un haut-parleur de manière à correspondre au filtrage basse-fréquence associé audit haut-parleur sur ladite plage de fréquence et au filtrage haute-fréquence en dehors de ladite plage de fréquence. On obtient ainsi pour les différents haut-parleurs des filtrages qui permettent à basse fréquence de compenser les effets de la salle, notamment par interaction entre les différents haut-parleurs ainsi filtrés, et à haute fréquence d'approcher le profil cible par une méthode particulièrement adaptée à ces hautes fréquences comme déjà indiqué. Ladite plage de fréquence peut comprendre en pratique les fréquences inférieures à une fréquence de séparation comprise entre 500 Hz et 2000 Hz, par exemple une fréquence de séparation de 1000 Hz, ou 1500 Hz dans l'exemple décrit ci-après. Ce type de valeur est particulièrement adapté pour la fréquence de séparation, notamment du fait qu'il corresponde à des longueurs d'onde de l'ordre de grandeur des dimensions la tête de l'auditeur (dimensions en deçà desquelles la variabilité déjà mentionnée devient particulièrement gênante). Selon une solution envisageable, particulièrement pratique, décrite ci-après, les filtrages basse-fréquence peuvent être déterminés par un processus d'optimisation de paramètres définissant ces filtrages, ledit processus comprenant une étape de génération des paramètres par tirage aléatoire. Par ailleurs, au moins un filtrage basse-fréquence peut par exemple comprendre la combinaison d'une pluralité de filtres à réponse impulsionnelle infinie. L'invention propose également un dispositif de détermination d'un filtrage à appliquer à une pluralité de haut-parleurs, caractérisé en ce qu'il comprend des moyens de détermination des réponses fréquentielles respectivement générées par les haut-parleurs de la pluralité en au moins un point, des moyens de détermination d'un profil d'énergie moyenne sur la base desdites réponses fréquentielles et des moyens de détermination du filtrage à appliquer sur la base du rapport entre un profil cible et le profil d'énergie moyenne. L'invention propose en outre un dispositif de détermination d'un ensemble de filtrages associés respectivement à une pluralité de haut-parleurs, caractérisé en ce qu'il comprend : - un dispositif de détermination d'un filtrage haute-fréquence conforme à ce qui vient d'être proposé ; - des moyens de détermination d'une pluralité de filtrages basse-fréquence associés chacun à un haut-parleur et tels que la réponse fréquentielle générée en au moins un point par les haut-parleurs avec application du filtrage associé corresponde à une réponse cible sur une plage de fréquence déterminée ; - des moyens de détermination du filtrage associé à un haut-parleur de manière à correspondre au filtrage basse-fréquence associé audit haut- parleur sur ladite plage de fréquence et au filtrage haute-fréquence en dehors de ladite plage de fréquence. Les avantages et caractéristiques optionnelles envisageables pour ces dispositifs sont analogues à ceux présentés ci-dessus en termes de procédés. L'invention propose enfin un programme d'ordinateur chargeable dans un système informatique, ledit programme contenant des instructions permettant la mise en oeuvre du procédé présenté ci-dessus, lorsque ce programme est chargé et exécuté par le système informatique. D'autres caractéristiques et avantages de l'invention apparaitront à la lecture de la description qui suit, faite en référence aux dessins annexés dans lesquels : - la figure 1 représente un exemple de procédé conforme aux enseignements de l'invention ; - la figure 2 représente un exemple de contexte dans lequel peut être mis en oeuvre le procédé de la figure 1 ; - la figure 3 présente un exemple de mise en oeuvre possible pour le détermination d'un ensemble de filtres par optimisation sur les basses fréquences ; - la figure 4 décrit un exemple de mise en oeuvre possible pour la génération d'un ensemble de filtres modifiés à partir d'un ensemble de filtres primaires ; - les figures 5 à 8 illustrent la mise en oeuvre d'une partie du procédé de la figure 4 sur un exemple de filtres primaires ; - la figure 9 illustre l'utilisation des filtres déterminés par le procédé de la figure 1 dans un système audio ; - la figure 10 représente un exemple de profil cible ; - la figure 11 représente un exemple de profil d'énergie moyenne ; - la figure 12 illustre la définition d'un profil cible corrigé ; - la figure 13 représente un exemple de filtre défini pour obtenir le profil cible corrigé sur la plage haute-fréquence ; - la figure 14 montre des filtres utilisables en pratique pour combiner des filtres adaptés respectivement à basse et haute fréquences. Le procédé représenté à la figure 1 débute par une étape E12 au cours de laquelle on mesure successivement la réponse fréquentielle générée par chaque haut-parleur d'une pluralité formée de Nc haut-parleurs HP1, HP2, HP3, HP4 en chaque point d'écoute d'une pluralité de N, points d'écoute PI, P2, P3 d'une salle représentée en figure 2. Les réponses fréquentielles mesurées sont celles obtenues sans application d'un filtrage avant application du signal électrique au haut-parleur concerné HP1, HP2, HP3, HP4. Chaque mesure mise en oeuvre pour un haut-parleur HP; à un point d'écoute Pi implique en pratique, pour une pluralité de fréquences réparties sur la bande de fréquence visée (par exemple par utilisation d'un bruit blanc sur toute ou partie de la bande de fréquences audibles 20 Hz û 20000 Hz), une mesure de l'amplitude et de la phase de l'onde acoustique au point concerné. La bande de fréquence étudiée est bornée par une fréquence inférieure Fmin (par exemple 20 Hz) et une fréquence supérieure Finax (par exemple 20 000 Hz). On peut donc à partir de telles mesures déterminer à l'étape E14 des réponses fréquentielles complexes ri,i(f) associées chacune à un haut-parleur HP; et à un point d'écoute Pi û l'amplitude de la mesure à la fréquence f au point Pi concerné correspondant au module de ri,i(f) et la phase correspondant à l'argument de ri,i(f). Ces étapes, de même que les étapes suivantes, sont par exemple mises en oeuvre au moyen d'un dispositif programmable, typiquement un système informatique (tel qu'un micro-ordinateur), doté d'une mémoire afin notamment de mémoriser les différentes valeurs traitées par les procédés décrits. Les étapes sont d'ailleurs par exemple mises en oeuvre du fait de l'exécution, par le dispositif programmable, d'instructions correspondantes d'un programme d'ordinateur également mémorisé dans le dispositif programmable. Le procédé se poursuit à l'étape E16 à laquelle on détermine un ensemble de Nc filtres (pi respectivement associés aux haut-parleurs HP; et tels que l'application de chaque filtre au haut-parleur correspondant permette d'optimiser les réponses acoustiques obtenues aux divers points d'écoute, selon un critère déterminé et sur une plage de fréquences inférieures à une fréquence de séparation FS (comprise par exemple entre 500 Hz et 2000 Hz, valant ici 1500 Hz). Les filtres ainsi déterminés permettent par exemple d'obtenir des réponses fréquentielles aux points d'écoute ayant un profil conforme approchant au mieux un profil cible (ou réponse souhaitée) T(f), éventuellement à un facteur (ou gain) près, sur la plage des fréquences étudiées et inférieures à la fréquence de séparation (fréquences comprises entre Fmin et Fs), soit ici entre 20 Hz et 1500 Hz. On note ici que le profil cible T(f) est de préférence défini sur l'ensemble de la plage de fréquences étudiée (ici 20 Hz û 20 000 Hz), comme visible en figure 10 pour un exemple de profil cible, mais que l'optimisation qui vient d'être mentionnée ne vise à atteindre ce profil cible que sur la partie de cette plage inférieure à la fréquence de séparation FS, l'autre partie du profil cible (ici 1500 Hz û 20 000 Hz) étant utilisée plus loin à l'étape E20. Un exemple de mise en oeuvre de l'étape E16 est décrit plus loin en référence à la figure 3. On procède ensuite à l'étape E18 à laquelle on détermine le profil d'énergie moyenne des réponses mesurées sur la plage des fréquences étudiées supérieures à la fréquence de séparation (ici 1500 Hz comme déjà indiqué). On entend par énergie moyenne la moyenne d'un paramètre représentatif de l'énergie (moyenne énergétique) sur au moins une longueur d'onde (moyenne dans l'espace à un instant donné), ou, ce qui est équivalent, sur au moins une période (en un point donné). Pour ce faire, on détermine par exemple d'abord le profil d'énergie moyenne en un point d'écoute Pi, ici estimé selon la formule : E;(f)= 1 Ni(f) i=1 >6 (f) où 6;(f) est un seuil qui permet de ne garder dans la somme que les réponses prépondérantes, de même ordre de grandeur (on peut prendre par 1 1 Nc exemple 6, (f) _ù.ù.L ri (f) , soit un seuil égal au dixième de la moyenne 10 Nc des modules des valeurs mesurées au point Pi), et Naf) est le nombre de telles réponses prépondérantes (soit N~(f) = card{i; r,i (f) > ai(f)}). En effet, la prise en compte des valeurs prépondérantes de même ordre de grandeur permet d'écrire pour chaque réponse : rit - ro cos(cot+0i), et le carré de l'énergie au point d'écoute concerné peut donc s'écrire : EN = ro2EN 1,... ,N ) avec Nc \2 2 a) 2n ENS =ù~o Lcos(cwt+0i) dt = Nc +2 Lcos(4i -4k) i=1 1 Aux fréquences considérées (ici supérieures à 1500 Hz), les valeurs de phase Oi varient de manière importante sur la zone d'écoute et la prise en compte d'une moyenne sur des dimensions au moins égales à celles de la tête d'un auditeur permettent d'écrire : (E)=r2(E(l,...,N))=Ncro2 (le terme en cosinus s'annulant dans la moyenne), ce qui implique que l'énergie moyenne varie avec la racine 5 carrée du nombre Nc de haut-parleurs, d'où la formule proposée ci-dessous pour le calcul de E1(f) . On remarque que l'on pourrait en variante ne pas limiter la sommation (utilisée pour le calcul de E;(f)) aux valeurs prépondérantes sans changement conséquent sur le résultat (puisque les valeurs non-10 prépondérantes ont une influence très limitée sur la somme). On déduit alors le profil d'énergie moyenne E(f) pour la salle en faisant la simple moyenne des profils d'énergie moyenne aux différents points Ni d'écoute Pi : E(f) = ùLE; (f) . Ni j=1 Un exemple d'un tel profil d'énergie moyenne est représenté en 15 figure 11 dans une salle où les niveaux sonores sont de l'ordre de 70 à 80 dB. On peut alors procéder à l'étape E20 où l'on choisit un filtre kit, par exemple de type à réponse impulsionnelle finie (ou FIR pour "Finite Impulse Response" selon une dénomination courante), adapté à transformer le profil d'énergie moyenne en un profil cible sur la plage des fréquences étudiées et 20 supérieures à la fréquence de séparation, soit la plage Fs û Fmax (ici 1500 Hz û 20 000 Hz). On remarque que le procédé d'optimisation présenté à la figure 3 en tant qu'exemple de mise en oeuvre de l'étape E16 permet d'obtenir des réponses Ri(f) conformes sur la plage Fmin û Fs (ici 20 Hz û 1500 Hz) au profil 25 cible T(f) à un facteur de gain près. Afin d'assurer une bonne continuité au niveau de la transition entre les deux plages de fréquences considérés ci-dessus, on propose ici d'utiliser comme profil cible sur la plage Fs û Fmax (ici 1500 Hz û 20 000 Hz) un profil cible corrigé F.T(f), où T(f) est le profil cible déjà mentionné et F est une 30 constante de gain qui maximise la corrélation sur la plage Fmin û Fs entre ce Ni profil corrigé F.T(f) et la réponse moyenne R(f) = 1 LR.(f) (obtenue par Ni j=1 optimisation sur cette dernière plage de fréquences), comme représenté dans l'exemple donné en figure 12. La corrélation maximum est donc en pratique obtenue pour la valeur de F qui minimise une distance d(F) entre F.T(f) et R(f), FS par exemple d(F) = J F.T(f) û R(f) df' (où on utilise une échelle logarithmique pour les fréquences f'= ln f ). Le filtre est donc (parmi la famille des filtres considérés, ici de type FIR) celui dont la fonction de transfert approxime au mieux sur la plage de fréquences considérée (Fs û Finax) le rapport F.T(f)/E(f). Pour plus de précisions sur ce type de filtres, on peut se référer par exemple à Digital Filters and Signal Processing, Third Edition, de Jackson, L.B., Kluwer Academic Publishers, Boston, 1996, pp. 301-307. Un exemple de filtre obtenu dans le cas de la figure 11 est représenté en figure 13. On procède alors à l'étape E22 de détermination des filtres Hi à appliquer respectivement aux haut-parleurs HPi. Chaque filtre Hi associé à un haut-parleur HPi donné est construit par combinaison du filtre (pi correspondant à ce haut-parleur HPi tel qu'il a été déterminé à l'étape E16, et du filtre déterminé à l'étape E20, de telle manière que le filtre Hi corresponde au filtre (pi sur la plage de fréquences Fmin û Fs et au filtre sur la plage de fréquences Fs ù Fmax. En pratique, on peut par exemple utiliser pour ce faire deux filtres de type FIR I, h représentés en figure 14, I étant un filtre passe-bas laissant passer les fréquences inférieures à Fs (ici de fréquence de coupure égale à Fs) et h étant un filtre passe-haut laissant passer les fréquences supérieures à Fs (ici de fréquence de coupure égale à Fs), et tels que l+h=l sur toute la plage de fréquences Fmin ù Fmax. On peut ainsi définir chaque filtre Hi sur toute la plage de fréquences Fmin ù Finax par : H; =1. (p; + h.. Le procédé proposé à la figure 3 pour mettre en oeuvre l'étape E16 débute par la détermination à l'étape E106 d'un ensemble de Nc filtres primaires cp;, chaque filtre primaire (pi étant associé à une voie, c'est-à-dire à un haut-parleur HP;. Chaque filtre primaire (pi est un filtre composite formé ici par l'application successive de N filtres à encoche, auxquels s'ajoutent un gain global G, une phase globale 0 et un retard global At. Les filtres utilisés dans le procédé de la figure 3 (et donc ceux produits par l'étape E16) sont ainsi des filtres à réponse impulsionnelle infinie (ou IIR pour "Infinite Impulse Response" selon la terminologie souvent utilisée). On choisit par exemple le gain G entre 0 et 1, le retard global entre 0 et 50 ms, la phase globale pouvant valoir 0° ou 180° (c'est-à-dire provoquer éventuellement une opposition de phase) dans un souci de simplicité (application d'un signe +1- au filtre) mais d'autres valeurs de phase pourraient être envisagées. Chaque filtre à encoche est défini par trois paramètres, à savoir sa fréquence centrale fo, son gain g (qui correspond au gain à la fréquence centrale fo) et sa largeur w. La fonction de transfert d'un tel filtre est : fol + 2igwfo f ù f 2 f02 + 2iwfo f ù f 2 Afin d'éviter des phénomènes de distorsion, on utilise ici des filtres à encoche dits à gain négatif, c'est-à-dire tel que leur gain g est inférieur ou égal à 1 (le terme "négatif' provenant du signe de la valeur de gain exprimée en dB). On prévoit par ailleurs de faire varier la largeur w sur une plage assez large afin de corriger non seulement les phénomènes existant sur une large plage de fréquences mais aussi les défauts intervenant sur une plage plus étroite. On prend par exemple des paramètres de largeur w compris entre 0,001 et 0,2. Chaque filtre primaire cp; est donc défini par 3N + 3 paramètres (trois paramètres pour chacun des N filtres à encoche et trois paramètres globaux). Les différents paramètres de chaque filtre primaire (pi sont déterminés arbitrairement à l'étape E106, par exemple pour chacun d'eux par tirage aléatoire sur la plage de valeurs envisagée (ou en variante en fixant chaque paramètre à une valeur prédéterminée, telle qu'une valeur médiane sur la plage envisagée). Une fois cet ensemble de filtres primaires déterminé, on procède à l'étape E108 à la génération d'une pluralité d'ensembles de filtres, chaque ensemble de filtres étant dérivé de l'ensemble de filtres primaires comme expliqué en détail plus bas en référence à la figure 4 et comprenant donc des filtres associés chacun à un haut-parleur HP; de la salle considérée, c'est-à-dire à une voie du système audio visé. On peut alors procéder à l'étape E110 à la détermination de l'ensemble de filtres (parmi les ensembles générés à l'étape E108) optimisant la qualité de reproduction sonore dans la salle. Pour ce faire, on procède tout d'abord au calcul, pour chaque ensemble de filtres généré, de la réponse obtenue en chaque point d'écoute Pj par application des filtres (pi de l'ensemble concerné au haut-parleur HP; associé : Nc R; (f) = L(P.(f ) , (f) . On évalue alors la qualité de ces réponses Rj(f) en les comparant à la réponse souhaitée T(f) en intégrant par exemple le rapport Rj(f)/T(f) en échelle linéaire (ou la différence en échelle logarithmique, typiquement exprimée en dB) sur la plage considérée F,,,in - Fs (à savoir ici entre 20 Hz et 1500 Hz), avec déduction de la moyenne de ce même rapport sur la plage afin de normaliser le résultat, ce qui permet d'obtenir pour chaque point d'écoute Pi une évaluation cj de la qualité de la réponse Rj(f) à ce point d'écoute, ici selon la formule : FS _ (f)-T(f)-R. +T df m,n 12 où Ri(f) et T(f) sont exprimés en dB et où Ri et T sont les moyennes sur la plage de fréquences considérée respectivement de Ri(f) et de T(f). L'évaluation globale de l'ensemble de filtres considéré (en considérant l'ensemble des points d'écoute Pi) est alors par exemple obtenue, en utilisant ici la moyenne des évaluations c m = Let /Ni et la moins bonne J évaluation c W = maxi cj , par : E =XEm+(1-X)Ew, où est une pondération adéquate, par exemple = 0,25. L'utilisation de la moins bonne évaluation (parmi les évaluations ci obtenues pour les points d'écoute Pi) permet d'éviter que la qualité recherchée soit obtenue par une forte amélioration en certains points d'écoute au détriment d'une mauvaise qualité en d'autres points d'écoute. On pourrait par ailleurs prévoir en variante que les évaluations aux différents points d'écoute soient affectées d'une pondération variable en fonction de l'importance relative des points d'écoute. L'ensemble de filtres déterminé à l'étape E110 comme celui optimisant la qualité est celui pour lequel l'évaluation globale c est minimale. On détermine alors à l'étape E112 si l'optimisation est terminée, par exemple du fait que la qualité obtenue pour l'ensemble déterminé à l'étape E110 est jugée suffisante (c'est-à-dire ici que l'évaluation globale c associée est inférieure à un seuil prédéterminé) ou qu'un nombre prédéterminé d'itérations des étapes E108 et E110 a été réalisé. S'il est déterminé à l'étape E112 que le processus d'optimisation n'est pas terminé, on procède à l'étape E114 à laquelle on remplace l'ensemble de filtres primaires utilisé lors de la précédente itération de l'étape E108 par l'ensemble de filtres déterminé à l'étape E110. On peut alors de nouveau procéder à l'étape E108 de génération d'une pluralité d'ensembles de filtres sur la base de ces nouveaux filtres primaires, puis à la détermination de l'ensemble parmi ces nouveaux ensembles générés optimisant un critère de qualité, ce qui constitue une nouvelle itération du processus d'optimisation. Si en revanche il est déterminé à l'étape E112 que le processus d'optimisation est terminé, le dernier ensemble de filtres déterminé à l'étape E110 est celui sélectionné comme le plus approprié pour compenser à basse fréquence les effets de la salle et qui est donc utilisé dans la suite du procédé représenté en figure 1 ; on procède donc alors à l'étape E18. La figure 4 propose un exemple de procédé de génération d'un ensemble de filtres à partir d'un ensemble de filtres primaires. Comme déjà annoncé, un tel procédé peut être utilisé dans le cadre de l'étape E108 précédemment décrite ; dans ce cadre, le procédé illustré à la figure 4 est mis en oeuvre pour générer chaque ensemble de la pluralité d'ensembles de filtres générés à l'étape E108. Le procédé de la figure 4 traite donc un ensemble de filtres primaires (pi dont un exemple est représenté en figure 5. Un tel filtre (pi est défini par plusieurs paramètres (ici 3N + 3 paramètres comme déjà indiqué), au sein desquels on peut définir des groupes de paramètres et des sous-groupes de paramètres utilisés au cours du traitement décrit ci-dessous. Dans le mode de réalisation décrit ici, les groupes sont des ensembles de trois paramètres qui correspondent soit aux paramètres globaux (G, cl), At) du filtre (pi, soit aux trois paramètres qui définissent un filtre à encoche. Les sous-groupes sont quant à eux ici formés chacun d'un seul paramètre. On pourrait en variante définir de différentes manières les groupes et les sous-groupes, certains groupes pouvant par exemple regrouper les paramètres de plusieurs filtres à encoche (ou d'autres filtres utilisés) et les sous-groupes pouvant correspondre à tout ou partie des ensembles formés par les éléments d'un groupe. Le procédé de la figure 4 débute à l'étape E302 par la détermination par tirage aléatoire d'un nombre n de groupes à modifier dans chaque filtre primaire, ce tirage aléatoire étant ici mis en oeuvre selon une demi-loi normale de probabilité (dont la variance est fixée à 10% du nombre de groupe de paramètres par exemple) sur une plage allant de 1 au nombre de groupes de paramètres (ici N+1). En pratique, on utilise par exemple une loi normale centrée sur zéro (espérance à zéro), de variance (N+1)/10 comme indiqué ci-dessus, et le nombre n est la valeur absolue, arrondie à l'entier supérieur (non-nul), du résultat de l'application de cette loi normale (les résultats strictement supérieurs à N+1, très rares en pratique du fait de la variance indiquée, sont par exemple considérés à la valeur N+1). L'algorithme de la figure 4 traite les Nc filtres primaires de l'ensemble de filtres primaires concernés (déterminés à l'étape E106 ou résultant d'une première phase d'optimisation comme indiqué à l'étape E114). On initialise donc une variable i à la valeur Nc à l'étape E304. On passe alors à l'étape E306 à laquelle on initialise à la valeur n déterminée à l'étape E302 une variable k. On procède alors à l'étape E308 au choix d'un groupe de paramètres du filtre primaire (pi (où i est la variable initialisée à l'étape E304) par tirage aléatoire (ici selon une loi uniforme de probabilité sur l'ensemble des groupes). On désigne par G le groupe ainsi choisi ; dans l'exemple de la figure 6, le groupe G correspond aux paramètres g2, w2, f2 d'un filtre à encoche. On procède ensuite à l'étape E310 au choix d'un sous-groupe du groupe choisi G, par tirage aléatoire ici selon une loi uniforme de probabilité sur l'ensemble des sous-groupes possibles. On désigne en figure 7 par SG le sous-groupe ainsi choisi, qui correspond dans l'exemple donné au paramètre de largeur w2 du filtre à encoche qui vient d'être mentionné. On remplace alors (étape E312) chaque paramètre du sous-groupe choisi par une valeur déterminée par tirage aléatoire dans le domaine de validité du paramètre (par exemple selon une loi uniforme de probabilité). Dans l'exemple illustré à la figure 8, le paramètre w2 précédemment mentionné est remplacé par une nouvelle valeur w'2 (le sous-groupe SG n'étant constitué que d'un paramètre). On procède alors à l'étape E314 à laquelle on décrémente la variable k. On teste à l'étape E316 si la variable k est nulle. Dans la négative, on retourne à l'étape E308 afin de modifier un nouveau groupe. Dans l'affirmative, les n groupes du filtre primaire concerné ont été modifiés et on passe à l'étape E318, à laquelle on décrémente la variable i afin de traiter un autre filtre primaire de l'ensemble. S'il résulte toutefois à l'étape E320 que la variable i atteint 0, tous les filtres primaires ont été traités et on met fin au procédé de la figure 4 (étape E322). Si en revanche la variable i n'est pas nulle, on retourne à l'étape E306 afin de traiter cet autre filtre primaire. On a représenté à la figure 9 un exemple de système audio auquel on applique l'ensemble des filtres H; déterminés par le processus de la figure 1. Une source audio S génère un signal sur plusieurs voies destinées chacune à un haut-parleur HP;. Sur chaque voie, on interpose, entre un amplificateur A; associé au haut-parleur HP; et la source S, un filtre H; associé au haut-parleur HP; dans l'ensemble de filtres obtenu par le processus de la figure 1. On obtient ainsi aux différents points d'écoute Pi la réponse acoustique souhaitée à basse fréquence comme à haute fréquence. Les modes de réalisation qui précèdent ne sont que des exemples possibles de mise en oeuvre de l'invention, qui ne s'y limite pas. REVENDICATIONS1. Procédé de détermination d'un filtrage à appliquer à une pluralité de haut-parleurs, caractérisé par les étapes suivantes : a) détermination des réponses fréquentielles respectivement générées par les haut-parleurs de la pluralité en au moins un point ; b) détermination d'un profil d'énergie moyenne sur la base desdites réponses fréquentielles ; c) détermination du filtrage à appliquer sur la base du rapport entre un profil cible et le profil d'énergie moyenne. 2. Procédé de détermination selon la revendication 1, caractérisé en ce que l'étape b) comprend la détermination d'un profil d'énergie moyenne associé audit point, incluant la détermination de la somme d'une partie au moins des modules des réponses fréquentielles générées audit point. 3. Procédé de détermination selon la revendication 2, caractérisé en ce que le profil d'énergie moyenne associé audit point est déterminé en divisant ladite somme par la racine carrée du nombre des modules supérieurs à une valeur de seuil. 4. Procédé de détermination selon l'une des revendications 1 à 3, caractérisé en ce que l'étape b) comprend la détermination d'une pluralité de profils d'énergie moyenne locaux associés chacun à un point et en ce que le profil d'énergie moyenne est déterminé comme la moyenne des profils locaux. 5. Procédé de détermination selon l'une des revendications 1 à 4, caractérisé en ce que le filtre déterminé à l'étape c) est un filtre à réponse impulsionnelle finie.30 6. Procédé de détermination d'un ensemble de filtrages associés respectivement à une pluralité de haut-parleurs, caractérisé par les étapes suivantes : - détermination d'un filtrage haute-fréquence selon un procédé conforme à l'une des revendications 1 à 5; - détermination d'une pluralité de filtrages basse-fréquence associés chacun à un haut-parleur et tels que la réponse fréquentielle générée en au moins un point par les haut-parleurs avec application du filtrage associé corresponde à une réponse cible sur une plage de fréquence déterminée ; - détermination du filtrage associé à un haut-parleur de manière à correspondre au filtrage basse-fréquence associé audit haut-parleur sur ladite plage de fréquence et au filtrage haute-fréquence en dehors de ladite plage de fréquence. 7. Procédé de détermination selon la revendication 6, caractérisé en ce que ladite plage de fréquence comprend les fréquences inférieures à une fréquence de séparation comprise entre 500 Hz et 2000 Hz. 8. Procédé de détermination selon la revendication 6 ou 7, caractérisé en ce que les filtrages basse-fréquence sont déterminés par un processus d'optimisation de paramètres définissant ces filtrages, ledit processus comprenant une étape de génération des paramètres par tirage aléatoire. 9. Procédé de détermination selon l'une des revendications 6 à 8, caractérisé en ce qu'au moins un filtrage basse-fréquence comprend la combinaison d'une pluralité de filtres à réponse impulsionnelle infinie. 10. Dispositif de détermination d'un filtrage à appliquer à une pluralité de haut-parleurs, caractérisé en ce qu'il comprend : - des moyens de détermination des réponses fréquentielles respectivement générées par les haut-parleurs de la pluralité en au moins un point ; - des moyens de détermination d'un profil d'énergie moyenne sur la base desdites réponses fréquentielles ; - des moyens de détermination du filtrage à appliquer sur la base du rapport entre un profil cible et le profil d'énergie moyenne. 11. Dispositif de détermination d'un ensemble de filtrages associés respectivement à une pluralité de haut-parleurs, caractérisé en ce qu'il comprend : - un dispositif de détermination d'un filtrage haute-fréquence conforme à la revendication 10 ; - des moyens de détermination d'une pluralité de filtrages basse- fréquence associés chacun à un haut-parleur et tels que la réponse fréquentielle générée en au moins un point par les haut-parleurs avec application du filtrage associé corresponde à une réponse cible sur une plage de fréquence déterminée ; - des moyens de détermination du filtrage associé à un haut-parleur de manière à correspondre au filtrage basse-fréquence associé audit haut-parleur sur ladite plage de fréquence et au filtrage haute-fréquence en dehors de ladite plage de fréquence. 12. Programme d'ordinateur chargeable dans un système informatique, ledit programme contenant des instructions permettant la mise en oeuvre du procédé selon l'une des revendications 1 à 9, lorsque ce programme est chargé et exécuté par le système informatique.