Domaine de l'invention La présente invention se rapporte au domaine du traitement du signal sonore.

La présente invention se rapporte plus particulièrement à un procédé et à un système d'égalisation acoustique automatisé.

Un cas d'usage de la présente invention est le suivant : dans le domaine de l'automobile : on place une pluralité de microphones dans un véhicule et on diffuse des séquences sonores pré-calibrées dans les haut-parleurs du véhicule. Un système compare les signaux sonores émis et les signaux sonores reçus et enregistrés. On en déduit « la signature acoustique » de l'habitacle du véhicule. L'utilisateur définit ensuite une courbe de signature acoustique cible, qui est différente de la signature acoustique native du véhicule. Un second algorithme calcule des coefficients de filtres numériques de manière à ce que, lorsqu'on applique ces filtres avant diffusion de signaux sonores dans les haut- parleurs du véhicule, la signature acoustique du véhicule devienne la courbe de signature sonore cible, et non la signature acoustique native du véhicule.

Dans le cadre de la présente invention, on utilise des filtres « IIR » ou « Infinité Impulse Response » en terminologie anglo-saxonne.

Plus particulièrement, on utilise dans le cadre de la présente invention des filtres dits « Biquad », d'ordre deux.

Le procédé et le système selon la présente invention concernent l'égalisation de l'amplitude de la réponse fréquentielle de l'habitacle.

Etat de la technique

On connaît dans l'état de la technique la demande de brevet français N° FR 2 967 848 (Centre Scientifique et Technique du Bâtiment), qui se rapporte à un système de correction de spectre destiné notamment à une salle de spectacle. Cette demande de brevet de l'art antérieur décrit un système électroacoustique comportant une pluralité de cellules. Dans ces cellules sont prévus : un dispositif d'égalisation, un émetteur, un récepteur, un circuit d'amplification pour amplifier les signaux issus du récepteur vers ledit émetteur et un organe de calcul qui va agir, entre autres, sur le dispositif d'égalisation. Cette solution technique de l'art antérieur se propose de mesurer la réponse émetteur, salle de spectacle observée au niveau du récepteur en utilisant, notamment, un signal de mesure fourni par un générateur de bruit ou toute autre méthode de mesure permettant d'observer la réponse du système en boucle ouverte. Le dispositif d'égalisation œuvre pour que ladite réponse se rapproche le plus d'une réponse désirée. Cette solution technique de l'art antérieur a pour application la modification des propriétés acoustiques de salles de spectacle.

L'art antérieur connaît également, par le brevet américain N° US 6 721 428 B1 (Texas Instruments), un égaliseur automatique de haut-parleurs. Ce brevet américain de l'art antérieur se rapporte plus particulièrement à un procédé pour générer des filtres numériques pour égaliser un haut-parleur. Des premières données numériques sont fournies, pour un intervalle de tolérance pour une courbe de réponse cible de signal sonore en fonction de la fréquence pour le haut-parleur. Des secondes données numériques sont générées, pour une courbe de réponse réelle du signal sonore en fonction de la fréquence pour le haut-parleur. Les premières données numériques sont comparées avec les secondes données numériques, et il est déterminé si la courbe de réponse réelle est dans l'intervalle de tolérance. Si la courbe de réponse réelle n'est pas dans l'intervalle de tolérance, des filtres audio numériques sont générés de façon itérative, et les filtres audio numériques sont appliqués aux secondes données numériques pour générer des troisièmes données numériques pour une courbe de réponse compensée. La fréquence, le gain et la bande passante des filtres audio numériques sont optimisés de façon automatique jusqu'à ce que la courbe de réponse compensée soit dans l'intervalle de tolérance ou une limite prédéterminée du nombre de filtres audio numériques ait été atteinte, celle des deux ayant lieu en premier. On connaît également dans l'état de la technique la publication scientifique « Filter Design Method for Loudspeaker Equalization Based on IIR Parametric Filters » de German Ramos et José J. Lopez. Exposé de l'invention

La présente invention propose un procédé permettant de fournir une égalisation d'un signal en déterminant des paramètres de filtres permettant de réduire l'écart entre l'amplitude d'une réponse fréquentielle représentant la signature acoustique d'un ensemble de haut-parleurs dans leur environnement et une courbe de signature sonore cible.

A cet effet, la présente invention concerne, dans son acception la plus générale, un procédé d'égalisation acoustique automatisé, caractérisé en ce qu'il comporte les étapes suivantes :

• Mesure de N réponses impulsionnelles R , Rl ₂ , RI N après émission d'un signal sonore précalibré reçu par N microphones ;

• Calcul des N réponses fréquentielles correspondantes par Transformée de Fourier rapide ;

· Etablissement d'une moyenne M des N réponses fréquentielles ;

• Traduction en échelle fréquentielle logarithmique de ladite moyenne M des N réponses fréquentielles ;

• Interpolation d'une courbe de signature sonore cible C _ec à partir d'un certain nombre de points choisis par un utilisateur ;

· Traduction en échelle fréquentielle logarithmique de ladite courbe de signature sonore cible C _ec ;

• Comparaison de ladite réponse moyennée M et de ladite réponse cible C _ec, en calculant la différence entre ladite réponse moyennée M et ladite réponse de signature sonore cible C _ec ; · Analyse de la courbe Cditt résultant de la différence entre ladite réponse moyennée M et ladite réponse cible C _{ec ;} et

• Détermination de paramètres de filtres pour la réduction de la différence entre ladite réponse moyennée M et ladite réponse cible C _ec en traitant tout d'abord les maxima locaux par ordre décroissant suivant leur gain, puis les minima locaux, et en réalisant des itérations successives ;

ledit procédé comportant en outre une étape d'optimisation des paramètres des filtres afin d'améliorer la performance du système.

Ainsi, le procédé selon la présente invention permet d'obtenir une égalisation acoustique automatisée grâce à un calcul précis et optimisé de paramètres de filtres.

De préférence, l'étape d'interpolation de la courbe de signature sonore cible est réalisée au moyen de la méthode d'Hermite.

Avantageusement, ledit procédé comporte en outre une étape d'optimisation automatique de l'offset de la réponse cible C _ec , répétée à chaque itération.

Avantageusement, ledit procédé comporte en outre une étape de lissage des N réponses fréquentielles.

De préférence, ledit procédé met en œuvre des filtres correspondants aux types suivants : « peak », « notch » « high-shelf » et « low-shelf » en fonction de la forme des maxima locaux et des minima locaux. Selon un mode de mise en œuvre particulier, ledit procédé met en outre en œuvre un algorithme d'optimisation globale pour minimiser l'erreur.

La présente invention se rapporte également à un système d'égalisation acoustique automatisé, caractérisé en ce qu'il comporte des moyens pour :

· mesurer N réponses impulsionnelles Rl-i , Rl ₂ , RIN après émission d'un signal sonore précalibré reçu par N microphones ;

• calculer les N réponses fréquentielles correspondantes par Transformée de Fourier rapide ;

• établir une moyenne M des N réponses fréquentielles ; • traduire en échelle fréquentielle logarithmique ladite moyenne M des N réponses fréquentielles ;

• interpoler une courbe cible C _ec à partir d'un certain nombre de points définis par un utilisateur;

· traduire en échelle fréquentielle logarithmique ladite courbe cible

• comparer ladite réponse moyennée M et ladite réponse cible C _ec, en calculant la différence entre ladite réponse moyennée M et ladite réponse cible C _ec ;

· analyser la courbe C _d itt résultant de la différence entre ladite réponse moyennée M et ladite réponse cible C _{ec ;} et

• déterminer des paramètres de filtres pour la réduction de la différence entre ladite réponse moyennée M et ladite réponse cible C _ec en traitant tout d'abord les maxima locaux par ordre décroissant suivant leur gain, puis les minima locaux, et en réalisant des itérations successives ;

ledit procédé comportant en outre des moyens pour optimiser des paramètres des filtres afin d'améliorer la performance du système.

Brève description des dessins

On comprendra mieux l'invention à l'aide de la description, faite ci-après à titre purement explicatif, d'un mode de réalisation de l'invention, en référence aux Figures dans lesquelles :

• la Figure 1 illustre les différentes étapes du procédé selon la présente invention ;

• la Figure 2 représente la courbe de signature sonore cible C _ec au sens de la présente invention, les réponses fréquentielles dérivées des N mesures de réponses impulsionnelles, ainsi que la moyenne M ; et

• la Figure 3 illustre la détection et le classement des maxima locaux (« pics ») et des minima locaux (« creux »). Description détaillée des modes de réalisation de l'invention La Figure 1 illustre les différentes étapes du procédé selon la présente invention.

Le procédé d'égalisation acoustique automatisé conforme à la présente invention comporte les étapes suivantes :

Dans un premier temps, on mesure N réponses impulsionnelles Rl-i , Rl ₂ ,

RI _N après émission d'un signal sonore précalibré reçu par N microphones.

Ensuite on calcule les N réponses fréquentielles correspondantes par Transformée de Fourier rapide.

Ensuite, on établit une moyenne M des N réponses fréquentielles calculées.

Une traduction en échelle fréquentielle logarithmique de ladite moyenne M des N réponses fréquentielles -est effectuée.

Une courbe de signature sonore cible C _ec est interpolée, puis est traduite en échelle fréquentielle logarithmique.

Ensuite, on compare la réponse moyennée M et la réponse de signature sonore cible C _ec, en calculant la différence entre la réponse moyennée M et la réponse cible C _ec .

La courbe C _d itt résultant de la différence entre ladite réponse moyennée M et ladite réponse cible C _ec est analysée.

Enfin, on détermine des paramètres de filtres pour la réduction de la différence entre ladite réponse moyennée M et ladite réponse cible C _ec en traitant tout d'abord les maxima locaux par ordre décroissant suivant leur gain, puis les minima locaux, et en réalisant des itérations successives.

Le procédé selon la présente invention comporte en outre une étape d'optimisation des paramètres des filtres afin d'améliorer la performance du système. Les réponses frequentielles peuvent être moyennées de façon « standard » (c'est-à-dire avec des poids identiques), ou bien avec des poids différents. La Figure 2 représente la courbe de signature sonore cible C _ec au sens de la présente invention, les réponses fréquentielles dérivées des N mesures de réponses impulsionnelles, ainsi que la moyenne M. Dans le cadre de la présente invention, on procède à une comparaison de la réponse moyennée M et de la réponse cible C _ec , en calculant la différence entre la réponse moyennée M et la réponse cible C _ec .

La Figure 3 illustre la détection et le classement des maxima locaux (« pics ») et des minima locaux (« creux »). Conformément à la présente invention, on traite tout d'abord les maxima locaux (pics) par ordre décroissant suivant leur gain, puis on traite les minima locaux (creux) par ordre croissant. Ceci permet de déterminer des paramètres de filtres pour la réduction de la différence entre la réponse moyennée M et la réponse cible C _ec . Des itérations successives sont réalisées. II a été démontré dans des études scientifiques qu'il est préférable d'égaliser d'abord les pics, puis ensuite les creux. En effet, l'oreille humaine est plus sensible aux pics qu'aux creux.

Dans un mode de réalisation, l'optimisation de l'offset de la courbe cible est réalisée comme suit :

La courbe cible et la réponse fréquentielle moyenne sont recalculées sur une échelle logarithmique afin d'approximer la résolution non uniforme du système auditif. Ceci est réalisé par une fonction de lissage qui rééchantillonne la réponse fréquentielle sur une échelle logarithmique avec par exemple une résolution fréquentielle de 1 /48 octave.

i) La bande de fréquence d'optimisation « FreqRange » est appliquée comme un vecteur de poids FreqWeight qui vaut 0 en dehors de la bande de fréquence et 1 à l'intérieur de la bande de fréquence. La valeur initiale de l'offset « Offset » (en dB) est calculée comme la valeur moyenne de la réponse fréquentielle moyenne dans la bande de fréquence d'égalisation :

Offset =

où n, et sont respectivement le premier et le dernier des points de fréquence de la bande de fréquence logarithmique d'égalisation.

L'algorithme d'optimisation consiste à trouver l'offset optimal, Offset, qui minimise l'erreur entre M et C _ec {=Shape+ Offset), définie comme suit : emean = {fk) ^~ C _ec (f _k \ Avec C _ec = Offset+ Shape

Ceci est réalisé avec l'algorithme d'optimisation, qui calcule de façon itérative l'erreur e _mea et recherche l'offset optimal dans un intervalle de +/- 100 dB autour de la valeur initiale. En plus de la minimisation de l'erreur e _mea n, une contrainte est ajoutée au problème d'optimisation, afin de limiter le gain des pics à G _max en dB. Il est défini comme suit :

max \M{f _k ) - C _ec {f _k ) \ ≤

tii≤k≤tif

Les paramètres et le nombre des filtres sont optimisés au moyen d'un algorithme. Les paramètres f, Q et G (respectivement fréquence centrale, facteur de qualité et gain des filtres biquads) sont optimisés à partir d'intervalles de valeurs qui peuvent être prédéfinis par un utilisateur, et les plages de valeurs de Q et G peuvent dépendre de la fréquence. Ainsi, par exemple dans les hautes fréquences, les filtres de faible gain sont plus facilement éliminés car ils ne sont pas perceptibles.

Dans un mode de réalisation, l'objectif est de trouver les paramètres optimaux (fc _op t, G _op t, Qopt) d'un filtre et l'offset optimal de la courbe cible Offset _op t. Les limites des paramètres sont déterminées comme ceci : J

max j—f-, FreqRange(l) I < f _c < min ( _c 2^, FreqRange l)

—G x 0.9 < G _opt ≤ G x 1.1 si G≤ 0 et - G x 1.1 < G _opt ≤ G x 0.9 si G > 0

QRangeÇl)≤ Q _opt ≤ Q Range (2)

TargetGain— 100 < TargetGain _opt ≤ TargetGain + 100

où f _c et G sont respectivement la fréquence centrale et le gain d'un filtre biquad modélisant le nième pic, et QRange est la plage de valeurs admissibles du facteur de qualité Q.

Un processus de post-optimisation est réalisé. Ce processus de postoptimisation consiste à reclasser les filtres par fréquence croissante et à réoptimiser les coefficients. Si un filtre est annulé lors de ce processus, un nouveau pic/creux est recherché afin d'avoir en sortie le nombre maximum de filtres. Le processus d'optimisation est mis en œuvre jusqu'à ce que le nombre maximum de filtres soit atteint. Dans un mode de réalisation, l'étape d'interpolation de la courbe cible est réalisée au moyen de la méthode d'Hermite.

Dans un mode de réalisation, le procédé selon la présente invention comporte en outre une étape d'optimisation automatique de l'offset de la réponse cible C _ec , répétée à chaque itération.

Dans un mode de réalisation, le procédé selon la présente invention comporte en outre une étape de lissage des N réponses fréquentielles. Dans le cadre du procédé selon la présente invention, on met en œuvre des filtres correspondants aux types suivants : « peak », « notch » « high- shelf » et « low-shelf » en fonction de la forme des maxima locaux (pics) et des minima locaux (creux).

Dans certains cas, il est préférable de choisir un filtre de type « peak ». Dans d'autres cas, il est préférable de choisir un filtre de type « shelf ». La sélection du filtre est effectuée en fonction du dépassement ou non d'un certain seuil par le facteur de qualité.

Dans un mode de réalisation, le procédé selon la présente invention met en outre en œuvre un algorithme d'optimisation globale pour minimiser l'erreur.

La présente invention se rapporte également à un système d'égalisation acoustique automatisé, comportant des moyens pour :

• mesurer N réponses impulsionnelles Rl-i , Rl ₂ , RIN après émission d'un signal sonore précalibré reçu par N microphones ;

• calculer les N réponses fréquentielles correspondantes par Transformée de Fourier rapide ;

• établir une moyenne M des N réponses fréquentielles ;

• traduire en échelle fréquentielle logarithmique ladite moyenne M des N réponses fréquentielles ;

• interpoler une courbe cible C _ec à partir d'un certain nombre de points définis par un utilisateur;

• traduire en échelle fréquentielle logarithmique ladite courbe cible

· comparer ladite réponse moyennée M et ladite réponse cible C _ec, en calculant la différence entre ladite réponse moyennée M et ladite réponse cible C _ec ;

• analyser la courbe C _d iff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible C _{ec ;} et · déterminer des paramètres de filtres pour la réduction de la différence entre ladite réponse moyennée M et ladite réponse cible C _ec en traitant tout d'abord les maxima locaux par ordre décroissant suivant leur gain, puis les minima locaux, et en réalisant des itérations successives ;

ledit système comportant en outre des moyens pour optimiser des paramètres des filtres afin d'améliorer la performance du système.

L'invention est décrite dans ce qui précède à titre d'exemple. Il est entendu que l'homme du métier est à même de réaliser différentes variantes de l'invention sans pour autant sortir du cadre du brevet.