FR2955442A1

FR2955442A1 - Method for determining filtration to be applied to set of loudspeakers in room in listening station, involves determining filtration to be applied to set of loudspeakers based on ratio between target profile and average energy profile

Info

Publication number: FR2955442A1
Application number: FR1050398A
Authority: FR
Inventors: Sebastien Lasserre
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2010-01-21
Filing date: 2010-01-21
Publication date: 2011-07-22
Anticipated expiration: 2030-01-21
Also published as: FR2955442B1

Abstract

The method involves determining frequency responses (E14) respectively generated by a set of loudspeakers at one point. An average energy profile is determined (E18) based on the frequency responses by determining the sum of frequency response modulus generated at the point. The average energy profile associated with the point is determined by dividing the sum by the square root of the number of modules greater than a threshold value. The filtration to be applied to the set of loudspeakers is determined (E22) based on the ratio between a target profile and the average energy profile. Independent claims are also included for the following: (1) a device for determining filtration to be applied to a set of loudspeakers (2) a computer program comprising a set of instructions for performing a filtration determining method.

Description

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. The invention relates to a filter determination method, as well as to a computer program and an associated device. The sound field generated in a room at a listening point results not only from the sounds directly emitted by the speakers arranged in the room, but also from the sounds reflected from the walls of the room and the sound perceived by a listener located in the room. The listening point can therefore be very different from the sound provided by the source, such as for example a CD or DVD player. It has therefore been proposed solutions for compensation of the effects of the room (in English "room compensation") which aim to bring as close as possible the sound perceived by the listener that received by the input sound reproduction system. Such compensation is for example made by playing on the filters applied to the different speakers. US Pat. No. 6,195,435 proposes in this context to configure a set of filters each associated with a loudspeaker for this purpose.

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. Documents US 2007/025 559 and US 2008/279 318 are also known, in which filterings to be applied at low frequencies and filterings to be applied at high frequencies are defined separately. In this context, the invention proposes a method for determining a filtering to be applied to a plurality of loudspeakers, characterized by the following steps: a) determining the frequency responses respectively generated by the loudspeakers of the plurality in minus one point; b) determining an average energy profile based on said frequency responses; c) determining the filtering to be applied based on the ratio between a target profile and the average energy profile.

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. Because of its low variability in space, the average energy associated with frequency responses is indeed a particularly appropriate quantity to determine a filtering that gives good results in the whole of the room concerned (and not only at the points considered). This advantage is felt particularly at high frequency where other quantities undergo significant variations in the listening space. According to the exemplary embodiment proposed and described below, step b) comprises the determination of a mean energy profile associated with said point, including the determination of the sum of at least a portion of the modules of the frequency responses generated. audit point. The average energy profile associated with said point is then for example determined by dividing said sum by the square root of the number of modules greater than a threshold value in order to take into account only the dominant sources of the same order of magnitude. As demonstrated later in the description, these steps make it possible to evaluate in a simple manner the average energy profile. It can further be provided that step b) comprises determining a plurality of local average energy profiles each associated with a point and that the average energy profile is then determined as the average of the local profiles. The filter determined in step c) is, for example, in practice a finite impulse response filter. The invention also proposes a method for determining a set of filterings associated respectively with a plurality of loudspeakers, characterized by the following steps: - determination of a high-frequency filtering according to a method according to what has been presented above; determination of a plurality of low-frequency filters each associated with a loudspeaker and such that the frequency response generated at at least one point by the loudspeakers with application of the associated filtering corresponds to a target response over a frequency range determined; determining the filtering associated with a loudspeaker so as to correspond to the low-frequency filtering associated with said loudspeaker over said frequency range and the high-frequency filtering outside said frequency range.

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é. For the various loudspeakers, this results in low-frequency filtering which makes it possible to compensate for the effects of the room, in particular by interaction between the different loudspeakers thus filtered, and at high frequency to approach the target profile by a particular method. adapted to these high frequencies as already indicated.

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. Said frequency range may comprise in practice the frequencies below a separation frequency of between 500 Hz and 2000 Hz, for example a separation frequency of 1000 Hz, or 1500 Hz in the example described below. This type of value is particularly suitable for the frequency of separation, especially since it corresponds to wavelengths of the order of magnitude of the dimensions of the listener's head (dimensions below which the aforementioned variability becomes particularly troublesome). According to one conceivable solution, particularly practical, described below, the low-frequency filtering can be determined by a parameter optimization process defining these filterings, said process comprising a step of generating the parameters by random draw. Furthermore, at least one low-frequency filtering may for example comprise the combination of a plurality of infinite impulse response filters. The invention also proposes a device for determining a filtering to be applied to a plurality of loudspeakers, characterized in that it comprises means for determining the frequency responses respectively generated by the loudspeakers of the plurality in at least a point, means for determining an average energy profile based on said frequency responses and filtering determining means to be applied based on the ratio of a target profile to the average energy profile. The invention furthermore proposes a device for determining a set of filterings associated respectively with a plurality of loudspeakers, characterized in that it comprises: a device for determining a high-frequency filtering according to which has just been proposed; means for determining a plurality of low-frequency filters each associated with a loudspeaker and such that the frequency response generated at at least one point by the loudspeakers with application of the associated filtering corresponds to a target response on a determined frequency range; means for determining the filtering associated with a loudspeaker so as to correspond to the low-frequency filtering associated with said loudspeaker over said frequency range and to the high-frequency filtering outside said frequency range. The optional advantages and characteristics that can be envisaged for these devices are similar to those presented above in terms of methods.

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). The invention finally proposes a computer program loadable in a computer system, said program containing instructions for implementing the method presented above, when this program is loaded and executed by the computer system. Other features and advantages of the invention will appear on reading the description which follows, made with reference to the accompanying drawings in which: - Figure 1 shows an example of a method according to the teachings of the invention; FIG. 2 represents an exemplary context in which the method of FIG. 1 can be implemented; FIG. 3 presents an exemplary implementation possible for the determination of a set of filters by optimization on the low frequencies; FIG. 4 describes an exemplary possible implementation for generating a set of modified filters from a set of primary filters; - Figures 5 to 8 illustrate the implementation of part of the method of Figure 4 on an example of primary filters; FIG. 9 illustrates the use of the filters determined by the method of FIG. 1 in an audio system; FIG. 10 represents an example of a target profile; FIG. 11 represents an example of average energy profile; FIG. 12 illustrates the definition of a corrected target profile; FIG. 13 represents an example of a filter defined to obtain the corrected target profile over the high-frequency range; FIG. 14 shows filters that can be used in practice to combine filters adapted respectively to low and high frequencies. The method represented in FIG. 1 starts with a step E12 during which the frequency response generated by each loudspeaker of a plurality formed of Nc loudspeakers HP1, HP2, HP3, HP4 is measured successively at each point of contact. listening to a plurality of N, listening points P1, P2, P3 of a room shown in FIG. 2. The frequency responses measured are those obtained without applying a filtering before applying the electrical signal to the speaker concerned HP1 , HP2, HP3, HP4. Each measurement implemented for an HP speaker; at a listening point Pi in practice involves, for a plurality of frequencies distributed over the target frequency band (for example by using white noise on all or part of the audible frequency band 20 Hz - 20000 Hz), a measurement of the amplitude and the phase of the acoustic wave at the point concerned. The frequency band studied is bounded by a lower frequency Fmin (for example 20 Hz) and a higher frequency Finax (for example 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). From such measurements it is therefore possible to determine, in step E14, complex frequency responses ri, i (f) each associated with an HP loudspeaker; and at a listening point Pi, the amplitude of the measurement at the frequency f at the relevant point Pi corresponding to the modulus of ri, i (f) and the phase corresponding to the argument of ri, i (f). These steps, as well as the following steps, are for example implemented by means of a programmable device, typically a computer system (such as a microcomputer), provided with a memory in particular for memorizing the different values. processed by the methods described. The steps are also for example implemented because of the execution, by the programmable device, corresponding instructions of a computer program also stored in the programmable device. The method is continued in step E16 at which a set of Nc filters (pi respectively associated with HP loudspeakers) is determined, and such that the application of each filter to the corresponding loudspeaker makes it possible to optimize the acoustic responses obtained. at the various listening points, according to a determined criterion and over a frequency range below a separation frequency FS (for example between 500 Hz and 2000 Hz, worth 1500 Hz here).

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. The filters thus determined make it possible, for example, to obtain frequency responses to the listening points having a conformal profile that best approximates a target profile (or desired response) T (f), possibly at a factor (or gain), on the frequency range studied and below the separation frequency (frequencies between Fmin and Fs), here between 20 Hz and 1500 Hz. It is noted here that the target profile T (f) is preferably defined over the entire frequency range studied (here 20 Hz - 20 000 Hz), as can be seen in figure 10 for an example of a target profile, but that the optimization just mentioned aims to reach this target profile only on the part of this target profile. range lower than the separation frequency FS, the other part of the target profile (here 1500 Hz - 20 000 Hz) being used later in step 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é). An exemplary implementation of step E16 is described below with reference to FIG. 3. Next, step E18 is used to determine the average energy profile of the responses measured over the range of frequencies studied above. at the separation frequency (here 1500 Hz as already indicated). Mean energy means the average of a parameter representative of the energy (energy average) over at least one wavelength (average in space at a given moment), or, which is equivalent, on at least one a period (at a given point).

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<i<k<Nc où le dernier terme tient compte des interférences entre sources (i.e. haut-parleurs). To do this, for example, the average energy profile is first determined at a listening point Pi, here estimated according to the formula: E (f) = 1 Ni (f) i = 1> 6 (f) where 6, (f) is a threshold which allows to keep in the sum only the preponderant answers, of the same order of magnitude (we can take by 1 1 Nc example 6, (f) _ù.ù.Li (f) , a threshold equal to one-tenth of the average of 10 Nc of the modules of the values measured at the point Pi), and Naf) is the number of such preponderant responses (ie N ~ (f) = card {i; r, i (f) > ai (f)}). Indeed, the taking into account of the dominant values of the same order of magnitude makes it possible to write for each answer: rit-ro cos (cot + 0i), and the square of the energy at the listening point concerned can therefore be write: EN = ro2EN 1, ..., N) with Nc \ 2 2 a) 2n ENS = ù ~ o Lcos (cwt + 0i) dt = Nc + 2 Lcos (4i -4k) i = 1 1 <i < k <Nc where the last term takes into account interferences between sources (ie loudspeakers).

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 : At the frequencies considered (here above 1500 Hz), the phase values O i vary considerably over the listening area and the taking into account of an average over dimensions at least equal to those of the head of a listener allow to write:

(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) . (E) = r2 (E (l, ..., N)) = Ncro2 (the cosine term being averaged), which implies that the mean energy varies with the square root of the Nc number speakers, hence the formula proposed below for the calculation of 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 Note that one could alternatively not limit the summation (used for the calculation of E; (f)) to preponderant values without

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. consequent change in the result (since the non-dominant values have a very limited influence on the sum). We then deduce the average energy profile E (f) for the room by simply averaging the average energy profiles at the different listening points Pi: E (f) = ùLE; (f). Ni j = 1 An example of such a mean energy profile is shown in FIG. 11 in a room where the sound levels are in the range of 70 to 80 dB. It is then possible to proceed to step E20 where a filter kit, for example of finite impulse response type (or FIR for "Finite Impulse Response" according to a current denomination), adapted to transform the average energy profile, is chosen. in a target profile over the frequency range studied and above the separation frequency, ie the range Fs-Fmax (here 1500 Hz - 20 000 Hz). It should be noted that the optimization method presented in FIG. 3 as an example of implementation of step E16 makes it possible to obtain responses Ri (f) conforming over the range Fmin-Fs (here 20 Hz - 1500). Hz) to the target profile T (f) at a gain factor close to. In order to ensure good continuity in the transition between the two frequency ranges considered above, it is proposed here to use as a target profile over the range Fs - Fmax (in this case 1500 Hz - 20 000 Hz) a target profile. corrected FT (f), where T (f) is the target profile already mentioned and F is a gain constant which maximizes the correlation over the range Fmin-Fs between this corrected profile FT (f) and the average response R ( f) = 1 LR. (f) (obtained by Ni j = 1 optimization on this last frequency range), as represented in the example given in FIG. 12. The maximum correlation is thus obtained in practice for the value of F which minimizes a distance d (F) between FT (f) and R (f), FS for example d (F) = J FT (f) - R (f) df '(where a logarithmic scale is used for the frequencies f' = ln f). The filter is (among the family of the considered filters, here FIR type) the one whose transfer function best approximates over the frequency range considered (Fs-Finax) the ratio F.T (f) / E (f). For more details on this type of filters, reference may be made, for example, to Digital Filters and Signal Processing, Third Edition, by Jackson, L.B., Kluwer Academic Publishers, Boston, 1996, pp. 301-307. An example of a filter obtained in the case of FIG. 11 is shown in FIG.

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. This is followed by the step E22 for determining the filters Hi to be applied respectively to the speakers HPi. Each filter Hi associated with a given loudspeaker HPi is constructed by combining the filter (pi corresponding to this speaker HPi as determined in step E16, and the filter determined in step E20, of in such a way that the filter Hi corresponds to the filter (p1 over the frequency range Fmin-Fs and to the filter over the frequency range Fs-Fmax In practice, for example, it is possible to use two filters of the FIR I, h type represented in FIG. 14, I being a low-pass filter allowing the frequencies lower than Fs to be passed (in this case having a cut-off frequency equal to Fs) and h being a high-pass filter allowing the frequencies above Fs to be passed (here the cut-off frequency equal to Fs), and such that l + h = 1 over the entire frequency range Fmin ù Fmax, so that each filter Hi can be defined over the entire frequency range Fmin ù Finax by: H; = 1. (p; + h .. The method proposed in Figure 3 to implement the E16 step debuted te by the determination in step E106 of a set of Nc primary filters cp ;, each primary filter (pi being associated with a channel, that is to say an HP speaker ;. Each primary filter (pi is a composite filter formed here by the successive application of N notch filters, to which are added a global gain G, a global phase 0 and a global delay At. The filters used in the method of FIG. 3 (and thus those produced by the step E16) are thus infinite impulse response filters (or IIR for "Infinite Impulse Response" according to the terminology often used), for example the gain G between 0 and 1, the delay between 0 and 50 ms, the global phase being able to be 0 ° or 180 ° (that is to say possibly to cause a phase opposition) for the sake of simplicity (application of a + 1 sign to the filter) but other phase values could be considered: Each notch filter is defined by three parameters, namely its center frequency fo, its gain g (which corresponds to the gain at the center frequency fo) and its width 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). The transfer function of such a filter is: fol + 2igwfo f ù f 2 f02 + 2iwfo f ù f 2 In order to avoid distortion phenomena, negative gain gain filters are used here, that is, that is, their gain g is less than or equal to 1 (the term "negative" from the sign of the gain value expressed in dB), and it is also expected to vary the width w over a fairly wide range in order to correct not only the phenomena existing over a wide range of frequencies but also the defects occurring on a narrower range, for example, parameters of width w between 0.001 and 0.2, each primary filter cp is therefore defined by 3N + 3 parameters (three parameters for each of the N notch filters and three global parameters).

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. The different parameters of each primary filter (pi are arbitrarily determined in step E106, for example for each of them by random drawing over the envisaged range of values (or alternatively by setting each parameter to a predetermined value, such as a median value over the range envisaged.) Once this set of primary filters has been determined, step E108 is carried out at the generation of a plurality of sets of filters, each set of filters being derived from the set of filters. as explained in detail below with reference to Figure 4 and thus comprising filters each associated with a speaker HP, the room in question, that is to say, a channel of the audio system referred to. proceed to step E110 to the determination of the set of filters (among the sets generated in step E108) optimizing the quality of sound reproduction in the room.

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) . To do this, it is first calculated, for each set of filters generated, the response obtained at each point of listening Pj by applying the filters (pi of the set concerned to the HP speaker; 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. The quality of these responses Rj (f) is then evaluated by comparing them with the desired response T (f) by integrating, for example, the ratio Rj (f) / T (f) into a linear scale (or the difference in logarithmic scale, typically expressed in dB) over the considered range F ,,, in - Fs (namely here between 20 Hz and 1500 Hz), with deduction of the average of this same ratio in the range in order to normalize the result, which makes it possible to obtain for each listening point Pi an evaluation cj of the quality of the response Rj (f) at this listening point, here according to the formula: FS _ (f) -T (f) -R. + T df m, n 12 where Ri (f) and T (f) are expressed in dB and where Ri and T are the averages over the considered frequency range of Ri (f) and T (f), respectively. The overall evaluation of the set of filters considered (considering the set of listening points Pi) is then obtained for example, by using here the average of the evaluations cm = Let / Ni and the least good J evaluation c W = maxi cj, by: E = XEm + (1-X) Ew, where is an adequate weighting, for example = 0.25. The use of the worst evaluation (among the evaluations obtained for the listening points Pi) makes it possible to prevent the quality sought from being obtained by a strong improvement in certain listening points at the expense of a poor quality in other listening points. Alternatively, the evaluations at the different listening points could be varied according to the relative importance of the listening points. The set of filters determined in step E110 as that optimizing the quality is that for which the overall evaluation c is minimal. It is then determined in step E112 whether the optimization is complete, for example because the quality obtained for the set determined in step E110 is considered sufficient (that is to say here that the overall evaluation c is less than a predetermined threshold) or a predetermined number of iterations of steps E108 and E110 has been performed. If it is determined in step E112 that the optimization process is not completed, proceed to step E114 to which the set of primary filters used during the previous iteration of step E108 is substituted. by the set of filters determined in step E110. It is then again possible to proceed to step E108 of generating a plurality of filter sets on the basis of these new primary filters, and then to determine the set of these new generated sets optimizing a quality criterion, which constitutes a new iteration of the optimization process. If, on the other hand, it is determined in step E112 that the optimization process is completed, the last set of filters determined in step E110 is the one selected as the most appropriate for compensating the effects of the room at low frequency and which is therefore used in the following process shown in Figure 1; we then proceed to step E18. Figure 4 provides an example of a method of generating a set of filters from a set of primary filters. As already announced, such a method can be used in the context of step E108 previously described; in this context, the method illustrated in FIG. 4 is used to generate each set of the plurality of filter sets generated in step E108. The method of FIG. 4 therefore deals with a set of primary filters (pi, an example of which is represented in FIG. 5). Such a filter (pi is defined by several parameters (here 3N + 3 parameters as already indicated), within which one can define parameter groups and subgroups of parameters used during the processing described below In the embodiment described here, the groups are sets of three parameters that correspond to either the global parameters (G, cl), At) of the filter (pi, or at the three parameters that define a notch filter), the subgroups are each formed of a single parameter, alternatively the groups and subgroups can be defined in different ways. , some groups may for example group the parameters of several notch filters (or other filters used) and subgroups that may correspond to all or part of the sets formed by the ele The process of FIG. 4 starts at step E302 by randomly determining a number n of groups to be modified in each primary filter, this random selection being implemented here in a half way. normal probability distribution (whose variance is fixed at 10% of the number of parameter groups for example) over a range from 1 to the number of groups of parameters (here N + 1). In practice, we use for example a normal distribution centered on zero (expectation at zero), variance (N + 1) / 10 as indicated above, and the number n is the absolute value, rounded up to the nearest integer ( non-zero), the result of the application of this normal law (the results strictly greater than N + 1, very rare in practice because of the variance indicated, are for example considered to the value N + 1). The algorithm of FIG. 4 processes the Nc primary filters of the set of primary filters concerned (determined in step E106 or resulting from a first optimization phase as indicated in step E114). Therefore, a variable i is initialized to the value Nc in step E304. We then go to step E306 to which is initialized to the value n determined in step E302 a variable k. We then proceed to step E308 at the choice of a group of parameters of the primary filter (pi (where i is the variable initialized in step E304) by random draw (here according to a uniform law of probability over all groups) G is the group thus selected, in the example of Figure 6, the group G corresponds to the parameters g2, w2, f2 of a notch filter, then proceed to step E310 at the choice of a subgroup of the chosen group G, randomly drawn here according to a uniform law of probability over the set of possible subgroups, and designated in Figure 7 by SG the subgroup thus chosen, which corresponds in the example given to the width parameter w2 of the notch filter which has just been mentioned, then replacing (step E312) each parameter of the subgroup chosen by a value determined by random draw in the range of validity of the parameter (for example according to a uniform law of probability.) In the illustrative example in Figure 8, the previously mentioned parameter w2 is replaced by a new W'2 value (SG subgroup being constituted as a parameter). We then proceed to step E314 at which we decrement the variable k. It is tested in step E316 if the variable k is zero. If not, return to step E308 to modify a new group.

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. If so, the n groups of the primary filter concerned have been modified and proceed to step E318, at which the variable i is decremented in order to treat another primary filter of the set. If, however, it results in step E320 that the variable i reaches 0, all the primary filters have been processed and the process of FIG. 4 is terminated (step E322). If, on the other hand, the variable i is not zero, we return to step E306 in order to process this other primary filter. FIG. 9 shows an example of an audio system to which all the filters H are applied; determined by the process of FIG. 1. An audio source S generates a multi-channel signal each intended for an HP speaker. On each channel, one interposes between an amplifier A; associated with the HP speaker; and the source S, a filter H; associated with the HP speaker; in the set of filters obtained by the process of FIG. 1. This gives the different listening points Pi the desired acoustic response at low frequency as well as at high frequency. The preceding embodiments are only possible examples of implementation of the invention, which is not limited thereto.

Claims

REVENDICATIONS1. A method of determining a filter to be applied to a plurality of loudspeakers, characterized by the steps of: a) determining the frequency responses respectively generated by the plurality of loudspeakers in at least one point; b) determining an average energy profile based on said frequency responses; c) determining the filtering to be applied based on the ratio between a target profile and the average energy profile.

2. Determination method according to claim 1, characterized in that step b) comprises determining a mean energy profile associated with said point, including determining the sum of at least part of the response modules. frequencies generated at that point.

3. Determination method according to claim 2, characterized in that the average energy profile associated with said point is determined by dividing said sum by the square root of the number of modules greater than a threshold value.

4. Determination method according to one of claims 1 to 3, characterized in that step b) comprises determining a plurality of local average energy profiles each associated with a point and in that the profile of mean energy is determined as the average of the local profiles.

5. Determination method according to one of claims 1 to 4, characterized in that the filter determined in step c) is a finite impulse response filter.

6. A method of determining a set of filterings associated respectively with a plurality of loudspeakers, characterized by the following steps: - determining a high-frequency filtering according to a method according to one of claims 1 to 5; determination of a plurality of low-frequency filters each associated with a loudspeaker and such that the frequency response generated at at least one point by the loudspeakers with application of the associated filtering corresponds to a target response over a frequency range determined; determining the filtering associated with a loudspeaker so as to correspond to the low-frequency filtering associated with said loudspeaker over said frequency range and the high-frequency filtering outside said frequency range.

7. Determination method according to claim 6, characterized in that said frequency range comprises the frequencies below a separation frequency between 500 Hz and 2000 Hz.

8. Determination method according to claim 6 or 7, characterized in that the low-frequency filtering is determined by a parameter optimization process defining these filtering, said process comprising a step of generating the parameters by random draw.

9. Determination method according to one of claims 6 to 8, characterized in that at least a low-frequency filtering comprises the combination of a plurality of infinite impulse response filters.

10. Device for determining a filtering to be applied to a plurality of loudspeakers, characterized in that it comprises: means for determining the frequency responses respectively generated by the loudspeakers of the plurality in at least one point ; means for determining an average energy profile on the basis of said frequency responses; means for determining the filtering to be applied on the basis of the ratio between a target profile and the average energy profile.

11. Device for determining a set of filterings associated respectively with a plurality of loudspeakers, characterized in that it comprises: a device for determining a high-frequency filtering device according to claim 10; means for determining a plurality of low-frequency filterings each associated with a loudspeaker and such that the frequency response generated at at least one point by the loudspeakers with application of the associated filtering corresponds to a target response on a determined frequency range; means for determining the filtering associated with a loudspeaker so as to correspond to the low-frequency filtering associated with said loudspeaker over said frequency range and the high-frequency filtering outside said frequency range.

12. A computer program loadable in a computer system, said program containing instructions for implementing the method according to one of claims 1 to 9, when the program is loaded and executed by the computer system.