Brevet FR3018015A1 — publié le 20150828

PROCEDE ET SYSTEME D'EGALISATION ACOUSTIQUE AUTOMATISE 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 « Infinite 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 oeuvre 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 Jose 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 RI1, RI2, RIN 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 Cec à partir d'un certain nombre de points choisis par un utilisateur ; - Traduction en échelle fréquentielle logarithmique de ladite courbe de signature sonore cible Cec ; - Comparaison de ladite réponse moyennée M et de ladite réponse cible Cec, en calculant la différence entre ladite réponse moyennée M et ladite réponse de signature sonore cible Cec ; - Analyse de la courbe Cdiff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; 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 Cecen 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 Cec, 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 oeuvre 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 oeuvre particulier, ledit procédé met en outre en oeuvre 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 RI1, RI2, ..., 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 Cec à partir d'un certain nombre de points définis par un utilisateur; - traduire en échelle fréquentielle logarithmique ladite courbe cible Cec - comparer ladite réponse moyennée M et ladite réponse cible Cec, en calculant la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; - analyser la courbe Cdiff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; 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 Cecen 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 Cec 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 RI1, RI2, RIN 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 15 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 Cec est interpolée, puis est traduite en échelle fréquentielle logarithmique. 20 Ensuite, on compare la réponse moyennée M et la réponse de signature sonore cible Cec, en calculant la différence entre la réponse moyennée M et la réponse cible Cec. La courbe Cdiff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible Cec est analysée. 25 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 Cecen 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 3 0 d'optimisation des paramètres des filtres afin d'améliorer la performance du système. Les réponses fréquentielles 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 Cec 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 Cec, en calculant la différence entre la réponse moyennée M et la réponse cible Cec. 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 Cec. Des itérations successives sont réalisées. Il 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 3 0 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. ii) 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 = mean(Ce,(nl: nf)) où ni et nf sont respectivement le premier et le dernier des points de fréquence de la bande de fréquence logarithmique d'égalisation. iii) L'algorithme d'optimisation consiste à trouver l'offset optimal, Offset, qui minimise l'erreur entre M et Cec (=Shape+Offset), définie comme suit : nf emean = 1 nf - ni + 1 IM(fk) Ceak) k=ni Avec Cec=Offset+Shape Ceci est réalisé avec l'algorithme d'optimisation, qui calcule de façon itérative l'erreur emean, et recherche l'offset optimal dans un intervalle de +/- 100 dB autour de la valeur initiale. En plus de la minimisation de l'erreur emean, une contrainte est ajoutée au problème d'optimisation, afin de limiter le gain des pics à Gmax en dB. Il est défini comme suit : max I M (fk) - Ce, (fk ) I Dans un mode de réalisation, l'objectif est de trouver les paramètres optimaux (fcopt, Gopt, Qopt) d'un filtre et l'offset optimal de la courbe cible Offsetopt Les limites des paramètres sont déterminées comme ceci : If \ lc 1 max 1 , FreqRange(1) 0 QRange(1) 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 Cec, 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 oeuvre 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 oeuvre 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 RI1, RI2, 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 Cec à partir d'un certain nombre de points définis par un utilisateur; - traduire en échelle fréquentielle logarithmique ladite courbe cible Cec - comparer ladite réponse moyennée M et ladite réponse cible Cec, en calculant la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; - analyser la courbe Cdiff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; 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 Cecen 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. REVENDICATIONS1. Procédé d'égalisation acoustique automatisé, caractérisé en ce qu'il comporte les étapes suivantes : - Mesure de N réponses impulsionnelles RI1, RI2, RIN 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 Cec à partir d'un certain nombre de points choisis par un utilisateur ; - Traduction en échelle fréquentielle logarithmique de ladite courbe de signature sonore cible Cec ; - Comparaison de ladite réponse moyennée M et de ladite réponse cible Cec, en calculant la différence entre ladite réponse moyennée M et ladite réponse de signature sonore cible Cec ; - Analyse de la courbe Cdiff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; 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 Cecen 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. 2. Procédé d'égalisation acoustique automatisé selon la 3 0 revendication 1, caractérisé en ce que l'étape d'interpolation de la courbe de signature sonore cible est réalisée au moyen de la méthode d'Hermite. 3. Procédé d'égalisation acoustique automatisé selon la revendication 1 ou 2, caractérisé en ce qu'il comporte en outre une étape d'optimisation automatique de l'offset de la réponse cible Cec, répétée à chaque itération. 4. Procédé d'égalisation acoustique automatisé selon la revendication 1, 2 ou 3, caractérisé en ce qu'il comporte en outre une étape de lissage des N réponses fréquentielles. 5. Procédé d'égalisation acoustique automatisé selon l'une au moins des revendications précédentes, caractérisé en ce qu'il met en oeuvre 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. 6. Procédé d'égalisation acoustique automatisé selon l'une au moins des revendications précédentes, caractérisé en ce qu'il met en outre en oeuvre un algorithme d'optimisation globale pour minimiser l'erreur. 7. Système d'égalisation acoustique automatisé, caractérisé en ce qu'il comporte des moyens pour : - mesurer N réponses impulsionnelles RI1, RI2, ..., 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 Cec à partir d'un certain nombre de points définis par un utilisateur; - traduire en échelle fréquentielle logarithmique ladite courbe cible Cec ,- comparer ladite réponse moyennée M et ladite réponse cible Cec, en calculant la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; - analyser la courbe Cdiff résultant de la différence entre ladite réponse moyennée M et ladite réponse cible Cec ; 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 Cecen 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.