FR 3051959 A1 20171201 FR 1751073 A 20170209 1 Procédé et dispositif pour estimer un signal déréverbéré. DOMAINE DE L'INVENTION La présente invention est relative aux procédés et 5 dispositifs pour estimer un signal déréverbéré. ARRIERE PLAN DE L'INVENTION Lorsqu'un signal acoustique originel est émis dans un milieu réverbérant puis capté par un microphone, le microphone capte un signal réverbéré dépendant du milieu 10 réverbérant. Dans la suite, on entend par « signal acoustique anéchoïque » le signal acoustique originel non réverbéré par un milieu. Un signal acoustique anéchoïque peut parfois être directement enregistré par un microphone, par 15 exemple lorsque le signal acoustique originel est émis dans une chambre anéchoïque. Toutefois, dans des conditions courantes d'enregistrement un microphone enregistre un signal acoustique réverbéré qui est un signal constitué du signal 20 acoustique originel reçu directement, mais également des réflexions du signal acoustique originel sur les éléments réverbérant du milieu, par exemple les murs d'une salle. Une forte réverbération acoustique du milieu peut être particulièrement gênante puisqu'elle dégrade la 25 qualité du son enregistré et réduit l'intelligibilité de la parole et sa reconnaissance par des machines. Pour résoudre ce problème, on connait des procédés et dispositifs qui permettent de reconstruire l'amplitude d'un signal déréverbéré à partir d'un signal acoustique 30 réverbéré par un milieu. Dans la présente demande, on entend par « signal déréverbéré » une estimation du signal acoustique originel, ou signal anéchoïque, obtenue par traitement analogique ou 3051959 2 digital d'un signal acoustique réverbéré enregistré par un microphone. A titre d'exemple, le document US2016035367 décrit un procédé de déreverbération qui permet de reconstruire un signal déréverbéré à partir d'un signal acoustique réverbéré par un milieu en calculant l'amplitude du signal déréverbéré dans plusieurs bandes de fréquences. Il existe un besoin pour améliorer encore les performances de tels procédés en estimant de manière plus précise les caractéristiques du signal déréverbéré à partir d'un signal acoustique réverbéré enregistré par un microphone. Un autre procédé est décrit dans le document « Restoration of instantaneous amplitude and phase of speech signal in noisy reverberant environments » de Yang Liu et al. paru dans les rapport de la 23ème conférence européenne sur le traitement des signaux (« 23rd European Signal Processing Conference »). Ce document décrit une méthode supervisée pour apprendre à un filtre de Kalman à reconstruire la phase et l'amplitude d'un signal déréverbéré en utilisant une base d'entraînement constituée de pair de signaux réverbéré et anéchoïque. Une telle base est toutefois compliquée à rassembler et les résultats obtenus dépendent de manière importante de la qualité de la base d'apprentissage et de l'adéquation entre les types de réverbération présents dans les signaux de la base d'entraînement et les réverbérations apparaissant dans les applications réelles. En outre, le procédé de déréverbération au moyen d'un filtre de Kalman décrit dans ce document permet uniquement de tenir compte de modulations linéaires de l'amplitude et de la phase, c'est-à-dire dans lesquelles les dérivées temporelles de l'amplitude et de la phase déréverbérées sont constantes au 3051959 3 cours du temps. La présente invention vient ainsi améliorer cette situation. OBJETS ET RESUME DE L'INVENTION 5 A cette fin un premier objet de l'invention est un procédé pour estimer une phase instantanée de signal acoustique déréverbéré. Le procédé comprend les étapes suivantes : (a) on mesure un signal acoustique réverbéré par 10 une propagation dans un milieu, (b) on estime au moins une transformée de Fourier à court terme du signal acoustique réverbéré avec au moins une fonction fenêtre, (c) on calcule au moins une fréquence instantanée 15 de signal déréverbéré à partir de ladite transformée de Fourier à court terme et d'un facteur d'influence du milieu, ledit facture d'influence étant fonction d'un temps de réverbération dudit milieu, et (d) on détermine au moins une phase instantanée de 20 signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps. Dans des modes de réalisation préférés de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : 25 - pour calculer au moins une fréquence instantanée de signal déréverbéré à partir de ladite transformée de Fourier à court terme, on estime pour chaque bande fréquentielle k d'une pluralité de N bandes fréquentielles une fréquence 30 instantanée lissée du signal réverbéré dans ladite bande fréquentielle k ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, 3051959 4 on calcule une fréquence instantanée de signal déréverbéré dans ladite bande fréquentielle k à partir de ladite fréquence instantanée lissée du signal acoustique réverbéré, du taux de variation au cours du 5 temps de ladite fréquence instantanée lissée du signal réverbéré et du facteur d'influence du milieu, et on détermine une phase instantanée de signal déréverbéré dans ladite bande fréquentielle k en intégrant la fréquence instantanée de signal déréverbéré dans la 10 bande fréquentielle k au cours du temps ; - le facteur d'influence du milieu est donné par : 1 min(t, Th) R(t) = 28 + 1 - e2Smin(t,Th) où ô et Th sont respectivement un facteur d'amortissement et une durée d'une décroissance exponentielle p(t)=e-st1[0,Thi de la réponse impulsionnelle du milieu, et le facteur 15 d'amortissement ô est calculé à partir d'un temps de réverbération mesuré dans le milieu, notamment un temps de réverbération RT60, par exemple tel que ô = 3.1og(10)/RT60 ; - pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k de la 20 pluralité de N bandes fréquentielles, on applique un algorithme de vocodeur réassigné; - pour calculer ladite au moins une fréquence instantanée de signal déréverbéré, on détermine un facteur de correction en multipliant le taux de variation au cours 25 du temps de la fréquence instantanée lissée du signal réverbéré par le facteur d'influence du milieu, en particulier on additionne ledit facteur de correction avec ladite fréquence instantanée lissée du signal acoustique réverbéré; 30 - pour calculer au moins une fréquence instantanée de signal déréverbéré à partir de ladite transformée de 3051959 5 Fourier à court terme, on calcule une pluralité de termes quadratiques de ladite au moins une transformée de Fourier à court terme pour chaque bande fréquentielle k d'une pluralité de N 5 bandes fréquentielles et chaque période de temps m d'une pluralité de périodes de temps, et on détermine, pour chaque bande fréquentielle k et chaque instant temporel m, une fréquence instantanée du signal déréverbéré et un taux de variation au cours du 10 temps de ladite fréquence instantanée du signal déréverbéré, en calculant une dérivée première et une dérivée seconde d'un paramètre dual solution d'un système linéaire dont les coefficients sont fonctions de ladite pluralité de termes quadratiques et du facteur d'influence 15 du milieu, ladite fréquence instantanée du signal déréverbéré étant une partie imaginaire de la dérivée première du paramètre dual et ledit taux de variation au cours du temps étant une partie imaginaire de la dérivée seconde du paramètre dual, 20 en particulier on inverse une matrice construite à partir de ladite pluralité de termes quadratiques et du facteur d'influence du milieu pour résoudre ledit système linéaire; - on estime au moins cinq transformées de Fourier à 25 court terme du signal acoustique réverbéré avec respectivement une première fonction fenêtre, une deuxième fonction fenêtre étant une dérivée première de la première fonction fenêtre, une troisième fonction fenêtre étant une dérivée seconde de la première fonction fenêtre, une 30 quatrième fonction fenêtre étant un produit de la première de la première fonction fenêtre avec une fonction linéairement croissante du temps et une cinquième fonction fenêtre étant une dérivée première de la quatrième fonction 3051959 6 fenêtre, et on calcule ladite pluralité de termes quadratiques à partir desdites au moins cinq transformées de Fourier à court terme; 5 - pour chaque bande fréquentielle k et chaque instant temporel m, on détermine une amplitude instantanée du signal déréverbéré à partir de ladite pluralité de termes quadratiques, et des dérivées premières et secondes du paramètre dual pour chaque bande fréquentielle k et 10 chaque instant temporel m; - pour déterminer au moins une phase instantanée de signal déréverbéré pour une bande de fréquence k, on détermine une bande de fréquence précédente k' permettant de minimiser une différence entre les fréquences centrales 15 h des fonctions fenêtre gi(t) et une fréquence estimée dans la bande de fréquence k, et on intègre une fréquence instantanée de signal déréverbéré et un taux de variation de ladite fréquence instantanée de signal déréverbéré pour ladite bande de fréquence précédente k'. 20 L'invention a également pour objet un dispositif pour estimer une phase instantanée de signal acoustique déréverbéré, comprenant : des moyens de mesure pour capter au moins un signal acoustique réverbéré par une propagation dans un milieu, 25 des moyens pour estimer au moins une fréquence instantanée lissée du signal réverbéré ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, des moyens pour calculer au moins une fréquence 30 instantanée de signal déréverbéré à partir de ladite fréquence instantanée du signal acoustique réverbéré, du taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré et d'un facteur 3051959 7 d'influence du milieu, ledit facture d'influence étant fonction d'un temps de réverbération dudit milieu, des moyens pour déterminer au moins une phase instantanée de signal déréverbéré en intégrant la fréquence 5 instantanée de signal déréverbéré au cours du temps. BREVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'un de ses modes de réalisation, donné à titre 10 d'exemple non limitatif, en regard des dessins joints. Sur les dessins : - la figure 1 est une vue schématique illustrant la réverbération du son dans une pièce lorsqu'un sujet parle de façon que ses paroles soient captées par un 15 dispositif selon une forme de réalisation de l'invention, - la figure 2 est un schéma de principe du dispositif de la figure 1, et - la figure 3 est un organigramme d'un procédé de reconstruction de signal déréverbéré selon un mode de 20 réalisation de l'invention, mettant en particulier en oeuvre un procédé d'estimation d'une phase instantanée de signal déréverbéré selon une forme de réalisation de l'invention. DESCRIPTION PLUS DETAILLEE Sur les différentes figures, les mêmes références 25 désignent des éléments identiques ou similaires. L'invention a pour but d'estimer une phase instantanée de signal acoustique déréverbéré à partir d'une mesure d'un signal acoustique réverbéré par une propagation dans un milieu 7, par exemple une pièce dans un bâtiment 30 telle que schématisée sur la figure 1. L'invention permet ainsi de traiter les signaux acoustiques captés par un dispositif électronique 1 pourvu d'un microphone 2. Le dispositif électronique 1 peut être 3051959 8 par exemple un téléphone dans l'exemple représenté, ou un ordinateur ou autre. Lorsqu'un son est émis dans le milieu 7, par exemple par une personne 3, ce son se propage jusqu'au 5 microphone 2 selon divers trajets 4, soit directs, soit après réflexion sur une ou plusieurs parois 5, 6 du milieu 7. Comme représenté sur la figure 2, le dispositif électronique 1 peut comporter par exemple une unité 10 centrale électronique 8 telle qu'un processeur ou autre, reliée au microphone 2 et à divers autre éléments, incluant par exemple un haut-parleur 9, un clavier 10, un écran 11. L'unité centrale électronique 8 peut communiquer avec un réseau externe 12, par exemple un réseau téléphonique. 15 L'invention permet au dispositif électronique 1 d'estimer une phase instantanée de signal acoustique déréverbéré. Dans une première application de grand intérêt, la phase instantanée de signal déréverbéré peut être utilisée 20 pour reconstruire un signal déréverbéré à partir d'un signal acoustique réverbéré. Pour cela, on mesure tout d'abord un signal acoustique réverbéré par une propagation dans le milieu. Puis, on détermine un spectre d'amplitude de signal 25 déréverbéré pour une pluralité de N bandes fréquentielles à partir du signal acoustique réverbéré. De nombreux procédés permettant de déterminer un spectre d'amplitude de signal déréverbéré à partir d'un signal acoustique réverbéré sont connus de l'art antérieur. 30 Ces procédés consistent par exemple à estimer un spectre de réverbération à partir du signal acoustique réverbéré puis à soustraire ledit spectre de réverbération du signal acoustique réverbéré. 3051959 9 On connait ainsi des procédés de détermination d'un spectre d'amplitude de signal déréverbéré utilisant : - la prédiction à long-terme tel que décrit dans l'article "Suppression of late reverberation effect on 5 speech signal using long-term multiple-step linear prediction" de K. Kinoshita, M. Delcroix, T. Nakatani, et M. Miyoshi paru dans IEEE Transactions on Audio, Speech, and Language Processing, vol. 17, no. 4, pp. 534-545, Mai 2009, 10 une modélisation stochastique de la réponse impulsionelle du milieu tel que décrit dans "A new method based on spectral subtraction for speech dereverberation" de K. Lebart et J.M. Boucher, paru dans ACUSTICA, vol. 87, no. 3, pp. 359-366, 2001, ou encore 15 - des réseaux de neurones profonds tel que décrit dans "Speech dereverberation for enhancement and recognition using dynamic features constrained deep neural networks and feature adaptation" de X. Xiao, S. Zhao, D.H. Ha Nguyen, X. Zhong, D.L. Jones, E.S. Chang, et H. Li, paru 20 dans EURASIP Journal on Advances in Signal Processing, vol. 2016, no. 1, pp. 1-18, 2016. Dans ces procédés de l'art antérieur, on reconstitue ensuite un signal déréverbéré à partir du spectre d'amplitude de signal déréverbéré obtenu et de la 25 phase du signal réverbéré. Il existe toutefois un besoin pour améliorer encore la qualité et l'intelligibilité du signal déréverbéré obtenu par ce procédé. A cette fin, selon l'invention, on détermine par 30 ailleurs, à partir du signal acoustique réverbéré une phase instantanée de signal déréverbéré pour chaque bande fréquentielle k de la pluralité de N bandes fréquentielles au moyen d'un procédé tel que décrit ci-après. 3051959 10 Puis, on reconstruit un signal déréverbéré partir du spectre d'amplitude de signal déréverbéré et de la phase estimée à l'aide de la méthode selon l'invention. On obtient de cette manière un signal déréverbéré 5 reconstruit de qualité nettement supérieure. La phase instantanée de signal déréverbéré déterminée par le procédé selon l'invention peut également avoir d'autres utilisations que la reconstruction du signal déréverbéré et peut par exemple être utilisée pour 10 améliorer la qualité et la précision d'un algorithme de localisation de source sonore tel que connu par ailleurs dans la littérature. Il est connu que le milieu réverbérant peut être modélisé par un modèle stochastique en définissant une 15 réponse impulsionnelle h(t) de la forme : h(t) = b(t)p(t) (1) dans laquelle b(t)-N(0,0-2) est un bruit blanc avec une distribution gaussienne centrée de variance 0-2 et p(C)= e-e110,71] est une décroissance exponentielle de la réponse 20 impulsionnelle du milieu où ô et Th sont respectivement un facteur d'amortissement et une durée de la réponse impulsionnelle du milieu. Un tel modèle stochastique est par exemple exposé dans la thèse de J.D. Polack, « La transmission de 25 l'énergie sonore dans les salles », soutenue à l'Université du Maine en 1988. Le facteur d'amortissement ô et la durée de la réponse impulsionnelle Th peuvent être déterminés à partir d'un temps de réverbération mesuré dans le milieu. 30 Un temps de réverbération couramment utilisé est le temps de réverbération à 60dB dénoté R160. Le temps de réverbération à 60dB correspond au temps requis pour que la 3051959 11 courbe de décroissance d'énergie (Energy Decay Curve) EDC diminue de 60 dB. Le temps de réverbération à 60dB peut par exemple être défini par la méthode d'intégration inverse de Manfred R. 5 Schroeder (New Method of Measuring Reverberation Time, The Journal of the Acoustical Society of Ameri.ca, 37(3) :409, 1965) par la courbe de décroissance d'énergie (Energy Decay Curve) EDC(n)=ENk_hnh(k)2 où h est la réponse impulsionnelle d'un milieu de longueur Nh et n est un indice temporel, par 10 exemple un nombre d'échantillons obtenus par un échantillonnage de pas temporel constant, n étant compris entre 1 et Nh. RT60 est alors le temps à l'indice temporel n requis pour que EDC(n) diminue de 60 dB. Des valeurs typiques du temps de réverbération RT60 15 sont par exemple des valeurs entre 0.4 s et 2 s. Bien que le temps de réverbération RT60 soit le plus couramment utilisé, on peut également utiliser un autre temps de réverbération caractéristique du milieu 7. On peut ensuite calculer le facteur d'amortissement du 20 milieu ô à partir du temps de réverbération RT60 par la formule ô = 3.1og(10)/RT60 . La durée de la réponse impulsionnelle Th peut également être définie à partir du temps de temps de réverbération, par exemple comme Th =a.RT60 où a peut 25 notamment être supérieur à 1, par exemple égale à 1.5. Toutefois, le facteur d'amortissement du milieu ô et la durée de la réponse impulsionnelle Th peuvent également être calculés par d'autres méthodes connues de l'art antérieur. 30 A partir du modèle statistique donné par l'équation (1), le signal acoustique réverbéré peut être relié au signal acoustique anéchoïque par l'équation de 3051959 12 convolution : y(t) = (h * s) (t) ( 2 ) où y(t) est le signal acoustique réverbéré et s(t) est le signal acoustique anéchoïque. 5 La phase instantanée du signal réverbéré peut par ailleurs être exprimée en fonction de la transformée de Hilbert du signal réverbéré comme : (P ev = arctan (9 (3) 3/(0 où cpr'(t) est la phase instantanée du signal 10 réverbéré et .9(t) est la transformée de Hilbert du signal réverbéré. on peut par ailleurs relier la fréquence instantanée du signal réverbéré à la phase instantanée du signal réverbéré par l'expression : 15 frev (0 = 1 dePrev(t) (t) (4) 2n dt Dans un premier mode de réalisation de l'invention, on peut estimer tout d'abord le taux de variation au cours du temps de la fréquence instantanée lissée du signal réverbéré. On peut alors déterminer la fréquence 20 instantanée du signal anéchoïque en fonction de l'espérance mathématique de la fréquence instantanée du signal réverbéré à partir des équations (1) à (4) comme : f E [frev (01 f + min(t,Th) 28 1-e2Smin(t,Th) (5) où f(t) est la fréquence instantanée du signal 25 anéchoïque estimée à l'instant t, E[fre'(t)] est l'espérance mathématique de la fréquence instantanée du signal réverbéré à l'instant t et f est le taux de variation au cours du temps de la fréquence instantanée du signal réverbéré. 30 L'espérance mathématique de la fréquence 3051959 13 instantanée du signal réverbéré à l'instant t ne peut pas être mesurée mais peut être approximée par un lissage temporelle de la fréquence instantanée du signal réverbéré mesurée. 5 On peut ainsi estimer une fréquence instantanée d'un signal déréverbéré en fonction d'une fréquence instantanée du signal réverbéré à partir des équations (1) (5) comme : (6) lf(t) = frev(t) + f ( -28 + -e2S( i(t,T h)) - 1 .. 10 où f(t) est la fréquence instantanée du signal déréverbéré estimée à l'instant t, f",'(0 est une fréquence instantanée lissée du signal réverbéré à l'instant t maintenant on lisse directement la stft et f est le taux de variation au cours du temps de la fréquence instantanée 15 lissée du signal réverbéré. L'équation (6) permet d'estimer une fréquence instantanée du signal déréverbéré en fonction de la fréquence instantanée lissée du signal réverbéré, du taux de variation au cours du temps de la fréquence instantanée et d'un facteur d'influence du milieu 20 R est donné par 1 min(t,T h) R (t) = 28 + 1-e28min(t,Th) (7 ) On peut ainsi réécrire l'équation (6) comme : f (0 = frev(t) + f R(t) (8 ) Une phase instantanée du signal déréverbéré 0(0 25 peut par la suite être déterminée par intégration temporelle comme : 0(0 = 2if fot - AT) cit + (P(0) ( 9) où (P(0) est une phase originelle du signal déréverbéré. 30 La fréquence et la phase du signal déréverbéré 3051959 14 estimées au moyen des équations (6) à (9) sont donc des estimations de la fréquence et de la phase du signal acoustique originel, ou signal anéchoïque. Les tests réalisés par les inventeurs indiquent que 5 ces estimations sont particulièrement bonnes puisqu'elles conduisent à un signal déréverbéré de qualité nettement supérieure à l'état de l'art. Un tel procédé peut encore être amélioré en déterminant directement à la fois la fréquence instantanée 10 du signal déréverbéré et le taux de variation de la fréquence instantanée du signal déréverbéré. Ceci permet notamment d'estimer de manière plus précise à la fois la phase et l'amplitude du signal déréverbéré. 15 Pour cela, on calcule plusieurs transformées de Fourier discrètes à court terme du signal réverbéré y(t) pour plusieurs fonctions fenêtre associées. Plus précisément, on définit une première fonction fenêtre gk(t) pour chaque bande fréquentielle k d'une 20 pluralité de N bandes fréquentielles, kE [0,N-1] et pour tout temps t, tE I La fonction fenêtre gk(t) est une fonction de réponse complexe d'un filtre passe bande analogique centré sur une fréquence fk. On définit ensuite en outre, une deuxième, une troisième, une quatrième et une 25 cinquième fonction fenêtre à partir de la première fonction fenêtre de la manière suivante : la deuxième fonction fenêtre gk(t) est une dérivée première de la première fonction fenêtre, la troisième fonction fenêtre gk(t) est une dérivée 30 première de la première fonction fenêtre, la quatrième fonction fenêtre g'k(t)=t.gk(t) est un produit de la première fonction fenêtre par la fonction 3051959 15 temps et la cinquième fonction fenêtre b'k(t) est une dérivée première de la quatrième fonction fenêtre. On calcule alors cinq transformées de Fourier à 5 court terme du signal acoustique réverbéré respectivement pour chacune desdites cinq fonctions fenêtre: Yg[m, k] = (g k * y)(tm) (10) 17,[m, k] = (g k * y)(tm) (11) yq [M, k] = (g k * y)(tm) (12) 10 Yg,[m, k] = (g' k * y)(tm) (13) yq, [M, k] = (g. ' k * y)(tm) (14) pour chaque bande fréquentielle k de la pluralité de bandes fréquentielles et chaque période de temps m (de manière équivalente tm) d'une pluralité de périodes de temps, où 15 trn=rnii et R est un facteur d'échantillonnage ou nombre fs. d'échantillons par période de temps et f", est une fréquence d'échantillonnage. A partir de la forme de la réponse impulsionnelle donnée en (1) et de la relation entre le signal acoustique 20 réverbéré et le signal acoustique anéchoïque donné par l'équation (2), on peut en déduire des relations entre les termes quadratiques des transformées de Fourrier discrète à court terme du signal acoustique anéchoïque et du signal acoustique réverbéré comme : 1.5'912 = -12 E [261'912 + 29Î(Y*gYj] S*S-9 2 = 1 -E [26Y* Y-9 + Y*99 Y- 9 + 1Y- 12] 9 a 1 Sg*Sg, = 2 E[28Y9*Y9, + r9179, + Yg* l'A Ci 3051959 16 1Sg,12 = 1_2 E [26117912 + 29Î(Y9*,rg,)] 1 r Sg*,Slq = --2 E[26179,Yà + Y*,.g,Yg + Yg*,Y:g1 où chaque terme est défini pour chaque bande fréquentielle k de la pluralité de bandes fréquentielles et chaque période de temps m d'une pluralité de périodes de temps 5 mais où les dépendances en k et m ont été masquées pour alléger la notation (c'est-à-dire que par exemple isgi2 dans 1 1 l'équation ci-dessus doit être vu en détail comme 1S9[m, k]12 ) Ici également, l'espérance mathématique des termes 10 peut être approximée par un lissage temporel et on peut obtenir les estimations : Is9l2 1 = 28117,12 + 29î(y9y9)) 7 ( (15) =6 f2(5'Yg* Yg + Yg* Yg + I Yg 12 ) (16) 1 SA, = 7 (28Y*9179, + Y*:gYg, + Y* Yg,) 1 (17) 15 Isg,12 = 071 (261y9,12 + 2wt(y9,,,rg,)) (18) 1 sg*,sg = 7 (2 8179*,17à + Y*,.g,Y'g + Y9, Yÿ) (19) Ici également, on peut définir un facteur d' influence du milieu R donné par 1 R 75 A partir de ces termes quadratiques et en réalisant un 20 développement de Taylor au second ordre du signal 3051959 17 anéchoïque s(t), on peut alors établir un système linéaire vérifié par les dérivées premières et seconde d'un paramètre dual Bk(t) =4(0+ Liiii(0 représentant le signal déréverbéré en notation exponentielle : s(t) = ek(t) = exp (E1/ REVENDICATIONS1. Procédé pour estimer une phase instantanée de 5 signal acoustique déréverbéré, le procédé comprenant les étapes suivantes : (a) on mesure un signal acoustique réverbéré par une propagation dans un milieu, (b) on estime au moins une transformée de Fourier à 10 court terme du signal acoustique réverbéré avec au moins une fonction fenêtre, (c) on calcule au moins une fréquence instantanée de signal déréverbéré à partir de ladite transformée de Fourier à court terme et d'un facteur d'influence du 15 milieu, ledit facture d'influence étant fonction d'un temps de réverbération dudit milieu, (d) on détermine au moins une phase instantanée de signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps. 20 2. Procédé selon la revendication 1, dans lequel pour calculer au moins une fréquence instantanée de signal déréverbéré à partir de ladite transformée de Fourier à court terme, on estime pour chaque bande fréquentielle k d'une 25 pluralité de N bandes fréquentielles une fréquence instantanée lissée du signal réverbéré dans ladite bande fréquentielle k ainsi qu'un taux de variation au cours du temps de ladite fréquence instantanée lissée du signal réverbéré, 30 on calcule une fréquence instantanée de signal déréverbéré dans ladite bande fréquentielle k à partir de ladite fréquence instantanée lissée du signal acoustique réverbéré, du taux de variation au cours du temps de ladite 3051959 30 fréquence instantanée lissée du signal réverbéré et du facteur d'influence du milieu et dans lequel on détermine une phase instantanée de signal déréverbéré dans ladite bande fréquentielle k en 5 intégrant la fréquence instantanée de signal déréverbéré dans la bande fréquentielle k au cours du temps. 3. Procédé selon la revendication 2, dans lequel, le facteur d'influence du milieu est donné par : 1 min(t,Th) R(t) = 28 + 1 - e où ô et Th sont respectivement un facteur d'amortissement 10 et une durée d'une décroissance exponentielle pM=e-st1[0,Thi de la réponse impulsionnelle du milieu, et dans lequel le facteur d'amortissement ô est calculé à partir d'un temps de réverbération mesuré dans le milieu, notamment un temps de réverbération RT60, par exemple tel que ô = 3.1og(10)/RT60 . 15 4. Procédé selon la revendication 2 ou 3, dans lequel, pour estimer une fréquence instantanée lissée du signal réverbéré pour chaque bande fréquentielle k de la pluralité de N bandes fréquentielles, on applique un algorithme de vocodeur réassigné. 20 5. Procédé selon l'une quelconque des revendications 2 à 4, dans lequel, pour calculer ladite au moins une fréquence instantanée de signal déréverbéré, on détermine un facteur de correction en multipliant le taux de variation au cours du temps de la fréquence instantanée 25 lissée du signal réverbéré par le facteur d'influence du milieu, en particulier dans lequel on additionne ledit facteur de correction avec ladite fréquence instantanée lissée du signal acoustique réverbéré. 30 6. Procédé selon la revendication 1, dans lequel pour calculer au moins une fréquence instantanée de signal 2 Smin(t,Th ) 3051959 31 déréverbéré à partir de ladite transformée de Fourier à court terme, on calcule une pluralité de termes quadratiques de ladite au moins une transformée de Fourier à court terme pour chaque bande fréquentielle k d'une pluralité de N bandes fréquentielles et chaque période de temps m d'une pluralité de périodes de temps, et on détermine, pour chaque bande fréquentielle k et chaque instant temporel m, une fréquence instantanée du signal déréverbéré et un taux de variation au cours du temps de ladite fréquence instantanée du signal déréverbéré, en calculant une dérivée première et une dérivée seconde d'un paramètre dual solution d'un système linéaire dont les coefficients sont fonctions de ladite pluralité de termes quadratiques et du facteur d'influence du milieu, ladite fréquence instantanée du signal déréverbéré étant une partie imaginaire de la dérivée première du paramètre dual et ledit taux de variation au cours du temps étant une partie imaginaire de la dérivée seconde du paramètre dual, en particulier on inverse une matrice construite à partir de ladite pluralité de termes quadratiques et du facteur d'influence du milieu pour résoudre ledit système linéaire. 7. Procédé selon la revendication 6, dans lequel on estime au moins cinq transformées de Fourier à court terme du signal acoustique réverbéré avec respectivement une première fonction fenêtre, une deuxième fonction fenêtre étant une dérivée première de la première fonction fenêtre, une troisième fonction fenêtre étant une dérivée seconde de la première fonction fenêtre, une quatrième fonction fenêtre étant un produit de la première de la première fonction fenêtre avec une fonction linéairement croissante 3051959 32 du temps et une cinquième fonction fenêtre étant une dérivée première de la quatrième fonction fenêtre, et dans lequel on calcule ladite pluralité de termes quadratiques à partir desdites au moins cinq 5 transformées de Fourier à court terme. 8. Procédé selon l'une quelconque des revendications 6 et 7, dans lequel pour chaque bande fréquentielle k et chaque instant temporel m, on détermine une amplitude instantanée du signal déréverbéré à partir de 10 ladite pluralité de termes quadratiques, et des dérivées premières et secondes du paramètre dual pour chaque bande fréquentielle k et chaque instant temporel m. 9. Procédé selon l'une quelconque des revendications 6 à 8, dans lequel pour déterminer au moins 15 une phase instantanée de signal déréverbéré pour une bande de fréquence k, on détermine une bande de fréquence précédente k' permettant de minimiser une différence entre les fréquences centrales h des fonctions fenêtre gi(t) et une fréquence estimée dans la bande de fréquence k, et on 20 intègre une fréquence instantanée de signal déréverbéré et un taux de variation de ladite fréquence instantanée de signal déréverbéré pour ladite bande de fréquence précédente k'. 10. Dispositif pour estimer une phase instantanée 25 de signal acoustique déréverbéré, comprenant : des moyens de mesure (2) pour capter au moins un signal acoustique réverbéré par une propagation dans un milieu, des moyens (8) pour estimer au moins une 30 transformée de Fourier à court terme du signal acoustique réverbéré avec au moins une fonction fenêtre, des moyens (8) pour calculer au moins une fréquence instantanée de signal déréverbéré à partir de ladite 3051959 33 transformée de Fourier à court terme et d'un facteur d'influence du milieu, ledit facture d'influence étant fonction d'un temps de réverbération dudit milieu, des moyens (8) pour déterminer au moins une phase 5 instantanée de signal déréverbéré en intégrant la fréquence instantanée de signal déréverbéré au cours du temps.