Brevet FR3031225A1 — publié le 20160701

Procédé de séparation amélioré et produit programme d'ordinateur La présente invention a pour domaine celui des procédés de séparation d'une pluralité de contributions dans un signal acoustique de mélange, et, en particulier, la séparation d'une contribution vocale, d'une contribution musicale de fond sonore, dans un signal acoustique de mélange. Une bande son d'une chanson comporte une contribution vocale (les paroles chantées par un ou plusieurs chanteurs) et une contribution musicale (la musique d'accompagnement jouée par un ou plusieurs instruments). Une bande son d'un film comporte une contribution vocale (les dialogues entre acteurs) superposée à une contribution musicale (les effets spéciaux sonores et/ou une musique de fond). Il est connu des algorithmes de séparation permettant de séparer la contribution vocale, de la contribution musicale, dans une bande son originale. Par exemple l'article de Jean-Louis Durrieu et al. "An iterative approach to monaural musical mixture de-soloing," in International Conference on Acoustics, Speech, and Signal Processing (ICASSP),Taipei, Taiwan, April 2009, pp. 105 - 108 divulgue un algorithme de séparation du type algorithme de séparation de sources sous-déterminée fondé sur une décomposition en matrices non-négatives, permettant de séparer la contribution vocale de la contribution de fond sonore. Cependant, les algorithmes de séparation connus ne permettent pas de prendre correctement en compte le phénomène de réverbération affectant les composantes du mélange. Dans le cas particulier d'une composante vocale, celle-ci résulte de la superposition de la voix sèche, ou pure dans ce qui suit, correspondant à l'enregistrement du son émis par le chanteur et qui s'est propagé directement vers le microphone d'enregistrement, et de la réverbération, correspondant à l'enregistrement du son émis par le chanteur mais qui s'est propagé indirectement vers le microphone d'enregistrement, c'est-à-dire par réflexion, éventuellement multiples, sur les parois de la salle d'enregistrement. La réverbération, constituée des échos de la voix pure à un instant donné, s'étale sur un intervalle de temps pouvant être important (par exemple trois secondes). Dit autrement, à un instant donné, la contribution vocale résulte de la superposition de la voix pure à cet instant et des différents échos de la voix pure à des instants précédents. Or, les algorithmes de séparation connus ne prennent pas en compte les effets à long terme de la réverbération affectant une composante du mélange. Le document de 3031225 2 Ngoc Q. K. Duong, Emmanuel Vincent, et Remi Gribonval, "Underdetermined reverberant audio source separation using a full-rank spatial covariance model," IEEE Transactions on Audio, Speech, and Language Processing, vol. 18, no. 7, pp. 1830 - 1840, Sept 2010 s'intéresse aux effets instantanés de diffusions spatiales de la réverbération, mais ne modélise pas les effets de 5 mémoire, c'est-à-dire la prise en compte du temps de latence entre l'enregistrement d'un son et l'enregistrement des échos associé à ce son. Ainsi, le type d'algorithme proposé par ce document ne s'applique qu'à des signaux multicanaux et ne permet pas une extraction correcte des effets de réverbération, que l'on peut trouver dans la musique. Dans le cas d'une composante vocale, la réverbération qui affecte cette composante est répartie dans les différentes composantes 10 obtenues à l'issue de la séparation. La composante vocale séparée perd de sa richesse et la composante musicale d'accompagnement n'est pas de bonne qualité. Il est à noter que la réverbération peut avoir pour cause les conditions dans lesquelles est réalisée la prise de son, mais peut également être ajoutée artificiellement au cours de la post-production de la bande son, essentiellement pour des raisons 15 esthétiques. Il y a donc un besoin pour un procédé permettant de séparer des contributions dans un mélange, ces contributions intégrant une réverbération du signal sonore pure correspondant. Plus particulièrement, il y a un besoin pour séparer une contribution vocale pure affectée par de la réverbération, d'une contribution musicale de fond sonore, 20 dans un signal sonore. L'invention a donc pour but de pallier ce problème. L'invention a donc pour objet un procédé de séparation, dans un signal acoustique de mélange, d'une contribution spécifique pure affectée par de la réverbération et d'une contribution de fond sonore, caractérisé en ce qu'il consiste à séparer la contribution 25 spécifique pure x(t) et la contribution de fond sonore z(t) en utilisant un spectrogramme de modélisation du signal acoustique de mélange V correspondant à la somme d'un spectrogramme d'une contribution spécifique réverbérée grev'Y et d'un spectrogramme de la contribution de fond sonore 17z, le spectrogramme de la contribution spécifique réverbérée dépendant du spectrogramme de la contribution pure Vx selon le modèle : T Vrev.y f,t =Ifyx t-r-F1Rf,t 2=1 30 où R est une matrice de réverbération, f est un pas de fréquence, t est un pas de temps, et r un entier entre 1 et T ; et en minimisant une fonction de coût entre le spectrogramme du signal de mélange et le spectrogramme de modélisation du signal de mélange. 3031225 3 Selon d'autres modes de réalisation, le procédé de séparation comporte une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou selon toutes les combinaisons techniquement possibles : la fonction de coût utilise une divergence entre le spectrogramme du signal 5 de mélange et le spectrogramme de modélisation du signal de mélange, notamment la divergence d'ITAKURA-SAITO. la contribution spécifique étant une contribution vocale, le spectrogramme de la contribution pure 17x est modélisé par : gx = (WF01-1F0)0(WKIIK) où WFO est une matrice d'atomes harmoniques, HFO est une matrice d'activation des 10 atomes harmoniques de la matrice WFO, WK est une matrice d'atomes de filtrage, HK est une matrice d'activation des atomes de filtrage de la matrice WK, et où O est un opérateur correspondant au produit terme à terme entre matrices. la minimisation de la fonction de coût met en oeuvre des règles de mise à jour multiplicatives du type : w70((wKHK) (R *t (17 grevO)3-1))) HFO HFO WFO ((wKHK) (R *t grev0I3-2)) WIT ((WFOHFO) O (R *t grev0)3-1))) 1/1/1 ((WFOHFO) (R *t grev03-2)) WK WK 0 ((WFOHFO) O (R *t (vo grevO)3-1))) Hic ((WFOHFO) O (R *t grev0 )3-2)) Hic 15 où O est un opérateur correspondant au produit terme à terme entre matrices (ou vecteur) ; .00 est un opérateur correspondant à l'exponentiation terme à terme d'une matrice par un scalaire ; (DT est la transposée d'une matrice ; et 1T est un vecteur T x 1 dont tous les éléments sont égaux à 1. la séparation de la contribution spécifique pure x(t) et la contribution de 20 fond sonore z(t) en utilisant un spectrogramme de modélisation du signal acoustique de mélange V constituant une seconde partie du procédé, celui-ci comporte une première partie consistant à séparer, dans le signal acoustique de mélange, une contribution spécifique et la contribution du fond sonore, sans tenir compte de la réverbération, le spectrogramme de la contribution spécifique étant utilisé comme valeur initiale du 25 spectrogramme de la contribution spécifique réverbérée lors de la minimisation de la fonction de coût dans une seconde partie du procédé. HK HK 0 3031225 4 la première partie comporte la minimisation d'une fonction de coût similaire à celle de la seconde partie. pour la minimisation de la fonction de coût, la première partie met en oeuvre des règles de mise à jour multiplicatives du type : WFO ((wKIIK) O (V O Vw-1)) HK HK W; ((WFOHFO) O (17W-2)) WK wK ((WFOHFO) O (V O 9(15'-1)))) H ((WFOHFO) O (f7w-2))Hii le procédé comporte, dans la première partie, à la suite de la minimisation de la fonction de coût, l'application d'un algorithme de suivi du maximum de puissance dans le spectrogramme de la contribution spécifique, ledit algorithme étant de préférence du type algorithme de Viterbi, puis la mise à zéro de tous les termes du spectrogramme 10 qui sont trop éloignés du maximum de puissance trouvé. le spectrogramme de la contribution de fond sonore qz est modélisé par une factorisation en matrices non-négatives : f7z = (WRHR) où WR est une matrice de modèles spectraux élémentaires et HR est une matrice d'activation des modèles élémentaires de la matrice WR, et en ce que la minimisation de 15 la fonction de coût met en oeuvre des règles de mise à jour multiplicatives du type : HR HR O WR ((7 C) - 2 ) ((V 017(15'-1)H17?' WR WR O ((VW-2)1117; où O est un opérateur correspondant au produit terme à terme entre matrices (ou vecteur) ; .00 est un opérateur correspondant à l'exponentiation terme à terme d'une matrice par un scalaire ; (. )T est la transposée d'une matrice ; et 1T est un vecteur T x 1 dont tous les éléments sont égaux à 1. 20 L'invention a également pour objet un produit programme d'ordinateur permettant la mise en oeuvre du procédé précédent. 5 HFO H -FO -n 1/1/70((WKHK) O (17(15'-2)) Wx ((WFOHFO) O (V 09(15'-1)) WR ((V O gw-1-) 3031225 5 L'invention sera mieux comprise à la lecture de la description qui va suivre d'un mode de réalisation particulier, donné uniquement à titre d'exemple illustratif et non limitatif, et faite en se référant aux dessins annexés sur lesquels : la figure 1 est une représentation sous forme de blocs des différentes étapes 5 du procédé de séparation selon l'invention ; et, les figures 2 et 3 correspondent à des graphes qui résultent de tests permettant de comparer, selon des critères normatifs connus, les résultats de la mise en oeuvre du procédé de la figure 1. En se référant à la figure 1, le procédé de séparation 100 utilise un signal 10 acoustique temporel de mélange w(t), pour délivrer un signal acoustique vocal y(t) et un signal acoustique musical z(t). Les signaux sont tous des signaux acoustiques, de sorte que le qualificatif d'acoustique sera omis dans ce qui suit. Ces signaux sont des signaux temporels. Ils dépendent du temps t. 15 Le signal acoustique de mélange est une bande son source, ou tout au moins un extrait d'une bande son. Le signal acoustique de mélange w(t) comprend une première contribution dite spécifique et une seconde contribution dite d'accompagnement. Dans la présente description, la première contribution est une contribution vocale 20 et correspond à des paroles chantées par un chanteur. La seconde contribution est une contribution musicale et correspond à l'accompagnement musical du chanteur. Le signal acoustique vocal y(t) correspond à la seule contribution vocale, isolée du reste du signal de mélange w(t), et le signal acoustique musical z(t) correspond à la 25 seule contribution musicale, isolée du reste du signal de mélange w(t). Dans le présent mode de réalisation, on considère que seule la contribution vocale est réverbérée. La réverbération est modélisée de la manière suivante : y(t) = r(t) * x(t) où x(t) est le signal vocal pur, c'est-à-dire le signal sonore généré par le chanteur est qui 30 s'est propagé directement vers le microphone d'enregistrement ; et où r(t) est une réponse impulsionnelle, qui est une distribution donnant l'amplitude des échos pour chaque instant d'arrivé de l'écho correspondant sur le microphone d'enregistrement, et où * correspond au produit de convolution. Le signal vocal pur x(t) est le signal en champ libre et la réponse impulsionnelle 35 r(t) est caractéristique de l'environnement acoustique de l'enregistrement. 3031225 6 Dans le domaine temps fréquence, pour des spectrogrammes non-négatifs, ce modèle de réverbération peut être approximé, tel que proposé dans le document de Rita Singh, Bhiksha Raj, et Paris Smaragdis, "Latent-variable decomposition based dereverberation of monaural and multi-channel signais," in IEEE International Conference 5 on Audio and Speech Signal Processing, Dallas, Texas, USA, March 2010, par : T ures,y =lux R V f,t f,t-r-F1f,r r=1 où Vres'Y est le spectrogramme du signal y(t), considéré comme affecté par de la réverbération, Vx est le spectrogramme du signal x(t), R est une matrice de réverbération correspondant au spectrogramme de la réponse impulsionnelle r(t), et T la dimension temporelle de H. 10 La première étape 110 du procédé 100 consiste à échantillonner le signal de mélange w(t) et à calculer un spectrogramme V du signal de mélange w(t). Ce spectrogramme est défini comme la valeur absolue (ou bien le carré de la valeur absolue) de la transformée de Fourier à court terme du signal w(t) échantillonné. Pour chaque pas d'échantillonnage temporel, le spectrogramme comporte une 15 trame en fréquence, indiquant pour chaque pas d'échantillonnage en fréquence, la puissance instantanée du signal. Le spectrogramme V est donc une matrice F x U, de nombres réels positifs U représente le nombre total de trames qui subdivisent la durée du signal du mélange w(t). F est le nombre total de pas d'échantillonnage en fréquence, qui vaut en 20 général entre 200 et 2000. Le procédé 100 comporte ensuite une première partie dans laquelle le signal vocal est considéré comme un signal vocal pure, sans réverbération. Dans cette première partie, le spectrogramme de modélisation du signal de mélange est la somme du spectrogramme du signal vocal 931, et du spectrogramme du 25 signal musical gz. f7Y est le spectrogramme du signal y(t), considéré comme non affecté par de la réverbération. Cette modélisation est finalement la modélisation usuelle dans le cadre des méthodes de décomposition par factorisation en matrices non-négatives. Il est à noter que â se réfère à une quantité qui est une estimation de la quantité a. Ainsi, dans les étapes de la première partir du procédé 100, on cherche à estimer 30 les deux spectrogrammes de sortie dont la somme approxime au mieux le spectrogramme du mélange : V V =93 +17z 3031225 7 La modélisation du signal vocal est fondée sur un modèle de production de la voix du type source / filtre, tel que proposé dans le document de Jean-Louis Durrieu et al. "An iterative approach to monaural musical mixture de-soloing," in International Conference on Acoustics, Speech, and Signal Processing (ICASSP),Taipei, Taiwan, April 2009, pp. 105 - 5 108 : 9Y = (WRIDHR0)0(WKHK) Le premier terme de cette modélisation est la source de la voix, qui correspond à l'excitation des cordes vocales : WFO est une matrice d'atomes harmoniques, qui est prédéfinie et est spécifique au chanteur ; HFO est une matrice d'activation indiquant à chaque instant les atomes harmoniques de la matrice WFO qui sont activés. 10 Le second terme de cette modélisation est le filtre de la voix, qui correspond au filtrage effectué par le conduit vocal : WK est une matrice d'atomes de filtrage ; HK est une matrice d'activation indiquant à chaque instant les atomes de filtrage de la matrice WK qui sont activés. L'opérateur G correspond à la multiplication matricielle terme à terme de deux 15 matrices (aussi dénommé produit d'Hadamard). La modélisation du signal musicale est fondée sur un modèle générique de factorisation par matrices non-négatives : 9z = (WRHR) Des colonnes de WR peuvent être vues comme des modèles spectraux élémentaires et HR comme une matrice d'activation de ces modèles élémentaires au fil du 20 temps. La première partie du procédé consiste alors à estimer les matrices HFO, WK, HK, WR et HR. Afin d'estimer les paramètres de ces matrices, une fonction de coût C, fondée sur une divergence d par élément, est utilisée : C = D(1719Y+17z) = Ef,td(vft I vft + eit) 25 Dans le mode de réalisation actuellement envisagé, la divergence d'Itakura-Saito, bien connue de l'homme du métier, est utilisée. Elle s'écrit : a a d(alb) = Tp- log Tp- 1 A l'étape 120, la fonction de coût C est ainsi minimisée de manière à déterminer la valeur optimale de chaque paramètre de chaque matrice. Cette minimisation est effectuée par itérations, avec des règles de mise à jour multiplicatives qui sont successivement 30 appliquées à chacun des paramètres des matrices HFO, WK, HK, WR et HR. 3031225 8 Ces règles de mise à jour sont par exemple élaborées en considérant le gradient (c'est-à-dire la dérivée partielle) de la fonction de coût C par rapport à chaque paramètre. Plus précisément, le gradient de la fonction de coût par rapport au paramètre considéré est écrit sous la forme d'une différence entre deux termes positifs, et la règle de mise à 5 jour correspondante est une multiplication du paramètre considéré par le rapport de ces deux termes. Cela permet notamment que les paramètres restent non négatifs à chaque mise à jour et deviennent constants lorsque le gradient de la fonction de coût par rapport au paramètre considéré tend vers zéro. 10 De cette manière, les paramètres évoluent vers un minimum local. Les règles de mise à jour sont ainsi les suivantes : W70 ((WKFIK) O (V O -1/W-1)) HR HR O WR ((7 C) - 2) ((V 017(15'-1)H17?' WR WR O ((VW-2)1117; où O est un opérateur correspondant au produit terme à terme entre matrices (ou vecteur) ; .00 est un opérateur correspondant à l'exponentiation terme à terme d'une matrice par un scalaire ; (.)T est la transposée d'une matrice. 15 Pour cette première partie, tous les paramètres sont initialisés avec des valeurs non-négatives choisies de manière aléatoire. Puis, à l'étape 130, la matrice HFO est contrainte en utilisant un algorithme de suivi tel que l'algorithme de « tracking » de Viterbi afin de sélectionner, pour chaque pas temporel, le pas en fréquence dans lequel on retrouve un maximum de puissance, sans 20 être trop éloigné en fréquence des maxima de puissance sélectionnés pour les pas temporels précédents. Puis, à l'étape 140, les coefficients de la matrice HFO qui sont à une distance en fréquence supérieure à une distance de référence sont fixés à 0. Une matrice 1-11F0 est obtenue. HFO HK HK W; ((WFOHFO) O (17W-2)) ((WFOHFO) O (V O 9(15'-1)))) Hx WK wK ((WFOHFO) O (f7w-2))1-iii H -FO -n 1/1/70((WKHK) O (17(15'-2)) Wx ((WFOHFO) O (V 09(15'-1)) WR ((V O 90fl-1-) 3031225 9 Dans la seconde partie du procédé 100, le signal vocal est considéré comme affecté par de la réverbération. La modélisation du signal vocal considéré comme réverbéré, grev.y, en fonction du signal vocal pure gx s'écrit alors : [grev.311t = [,73, *t R = x 17 f,t-r-FtR f,r 2=1 5 où *t dénote un opérateur de convolution ligne par ligne tel qu'explicité dans le membre de droite de l'équation ci-dessus. La matrice de réverbération R comporte T pas de temps (de même durée qu'un pas d'échantillonnage du signal de mélange), et F pas d'échantillonnage en fréquence. T est prédéterminé par l'utilisateur et vaut généralement entre 20 et 200, par exemple 100. 10 De plus, comme ci-dessus, le spectrogramme f7x du signal pure est modélisé par : gx = (WF01-1F0)0(WKIIK) La seconde partie du procédé consiste alors à estimer les matrices HFO, WK, HK, W125 HR et R qui permettent d'approximer le spectrogramme du mélange V : grev = grev,y gz Afin d'estimer les paramètres de ces matrices, une fonction de coût C, fondée sur une divergence d par élément, est utilisée : C = D(17117rev'Y +17z) =Ef,td(Vftlgfrtv'Y + Vit) 15 Dans le mode de réalisation actuellement envisagé, la divergence d'Itakura-Saito, bien connue de l'homme du métier, est utilisée. Elle s'écrit : a a d(alb) = Tp- log Tp- 1 Avantageusement, la fonction de coût de la seconde partie est similaire à celle utilisée dans la première partie. A l'étape 220, la fonction de coût C est alors minimisée de manière à déterminer la 20 valeur optimale de chaque paramètre de chaque matrice, en particulier les paramètres de la matrice de réverbération. Cette minimisation est effectuée par itérations avec des règles de mise à jour multiplicatives, qui sont successivement appliquées à chacun des paramètres des matrices. Pour les matrices de la composante vocale intégrant une réverbération, on a : 3031225 WK WK 0 ((WFOHFO) O (R *t grevC))3-2)) Hic où *t désigne l'opérateur de convolution ligne par ligne tel que défini ci-dessus. Pour la composante musicale, on a, comme dans la première partie du procédé : 10 Wjo Weil() O (R VrevC))3-1))) HFO HFO W70 ((WKHK) O (R *t grevC))3-2)) WIT ((WFOHFO) O (R *t grevC))3-1))) WIT ((WFOHFO) (R *t grevC))3-2)) ((WFOHFO) O (R *t grevC))3-1))) HK HK 0 wiT((v grevC))3-1) HR HR 0 WRerevC))3-2) ((wRHR) (V VrevC))3-1)Hii WR WR erevC))3-2)H7R' En ce qui concerne la matrice HFO, les itérations partent de la matrice Hi Fo déterminée dans la première partie du procédé. Il est à noter que, puisque les règles de 5 mise à jour sont multiplicatives, les coefficients de la matrice HFO fixés initialement à 0 resteront à 0 au cours de la minimisation de la fonction de coût dans la seconde partie du procédé. Lorsque la distance entre le spectrogramme de mélange V et le spectrogramme estimé gr" + gz est inférieure à un seuil prédéterminé ou lorsqu'un nombre d'itérations 10 limite fixé à l'avance est atteint, le procédé sort de la boucle d'itération et les valeurs des matrices R, HFO, WK, HK, WR et HR sont les valeurs finales. A l'étape 230, des traitements adaptés classiques (en particulier un traitement du type filtrage de Wiener) sont appliqués sur les spectrogrammes précédents pour obtenir notamment les spectrogrammes d'intérêt gx, gz. Puis, à l'étape 240, une transformation 15 inverse de celle de l'étape 110 est réalisée sur ces spectrogrammes pour obtenir les signaux de sorties, signal vocal pure x(t) et signal musical z(t). Dans les modes de réalisation décrits ici en détail, ces signaux acoustiques sont des signaux monophoniques. En variante, ces signaux sont stéréophoniques. Plus généralement encore, ils sont multicanaux. L'homme du métier sait comment adapter à 20 des signaux stéréophoniques ou multicanaux les traitements présentés pour le cas de signaux monophoniques. Le mode de réalisation préféré est relatif à une composante spécifique ou d'intérêt qui est une composante vocale. Cependant, la modélisation de la réverbération d'une 3031225 11 composante est générale et s'applique à tout type de composante. En particulier, la composante de fond sonore peut également être affectée par une réverbération. De plus, n'importe quel type de modélisations non-négatives des spectrogrammes des sons non réverbérés peut également être utilisées, en lieu et place de celles utilisées 5 ci-dessus. Par ailleurs, dans le mode de réalisation présenté ci-dessus, le mélange comporte deux composantes. La généralisation à un nombre quelconque de composantes est directe. Des tests comparatifs ont été menés afin de comparer les résultats de la mise en 10 oeuvre du présent procédé : - le premier procédé est une séparation, fondée sur une méthode de type NMF, sans inclure de modélisation sur la réverbération ; - le second procédé est une séparation selon le procédé décrit ci-dessus, c'est-à-dire incluant une modélisation de la réverbération du signal vocal ; et, 15 - le troisième procédé est une limité mathématique théorique. Afin de quantifier les résultats obtenus pour les différents procédés, des indicateurs standards du domaine de la séparation de sources ont été calculés. Ces indicateurs sont le rapport signal sur distorsion SDR (selon l'acronyme anglais « Signal to Distorsion Ratio »), et qui correspond à un test quantitatif ; le rapport signal sur artefact 20 SAR (selon l'acronyme « Signal to Artefact Ratio »), et qui correspond aux artefacts dans les composantes séparées ; et le rapport signal sur interférence SIR (selon l'acronyme anglais « Signal to Interference Ratio »), et qui correspond aux interférences résiduelles entre les composantes séparées. Les résultats sont présentés sur les figures 2 pour le signal vocal et la figure 3 25 pour le signal musical. Le procédé selon l'invention améliore donc les résultats obtenus, quelle que soit la manière de les analyser. REVENDICATIONS1.- Procédé de séparation (100), dans un signal acoustique de mélange (w(0), d'une contribution spécifique pure, affectée par de la réverbération, et d'une contribution de fond sonore, caractérisé en ce qu'il consiste à séparer la contribution spécifique pure x(t) et la contribution de fond sonore z(t) en utilisant un spectrogramme de modélisation du signal acoustique de mélange V correspondant à la somme d'un spectrogramme d'une contribution spécifique réverbérée firev'Y et d'un spectrogramme de la contribution de fond sonore 17z, le spectrogramme de la contribution spécifique réverbérée dépendant du spectrogramme de la contribution spécifique pure 17x selon le modèle : T 9fr,te12.31 =II2fx,t-r-F1Rf,t T=1 où R est une matrice de réverbération, f est un pas de fréquence, t est un pas de temps, et T un entier entre 1 et T ; et en minimisant une fonction de coût (C) entre un spectrogramme du signal de mélange et le spectrogramme de modélisation du signal de mélange. 2.- Procédé selon la revendication 1, caractérisé en ce que la fonction de coût (C) utilise une divergence (d) entre le spectrogramme du signal de mélange et le spectrogramme de modélisation du signal de mélange, notamment la divergence d'ITAKURA-SAITO. 3.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, la contribution spécifique pure étant une contribution vocale, le spectrogramme de la contribution spécifique pure gx est modélisé par : f7x = (WFTIFIF0)0(WKI/K) où WFO est une matrice d'atomes harmoniques, HFO est une matrice d'activation des atomes harmoniques de la matrice WFO, FO, WK est une matrice d'atomes de filtrage, 1/K est une matrice d'activation des atomes de filtrage de la matrice WK, et où O est un opérateur correspondant au produit terme à terme entre matrices. 4.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que, la séparation de la contribution spécifique pure x(t) et la contribution de fond sonore z(t) en utilisant un spectrogramme de modélisation du signal 3031225 13 acoustique de mélange 17 constituant une seconde partie du procédé, celui-ci comporte une première partie consistant à séparer, dans le signal acoustique de mélange (w(t)), une contribution spécifique et la contribution du fond sonore, sans tenir compte de la réverbération, le spectrogramme de la contribution spécifique étant utilisé comme valeur 5 initiale du spectrogramme de la contribution spécifique réverbérée lors de la minimisation de la fonction de coût dans une seconde partie du procédé. 5.- Procédé selon la revendication 4, caractérisé en ce que la première partie comporte la minimisation d'une fonction de coût similaire à celle de la seconde partie. 10 6.- Procédé selon la revendication 5, caractérisé en ce qu'il comporte, dans la première partie du procédé, à la suite de la minimisation de la fonction de coût, l'application d'un algorithme de suivi du maximum de puissance dans le spectrogramme de la contribution spécifique, le dit algorithme étant de préférence du type algorithme de 15 Viterbi, puis la mise à zéro de tous les termes du spectrogramme qui sont trop éloignés du maximum de puissance trouvé. 7.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le spectrogramme de la contribution de fond sonore qz est 20 modélisé par une factorisation en matrices non-négatives : = (WRHR) où WR est une matrice de modèles spectraux élémentaires et HR est une matrice d'activation des modèles élémentaires de la matrice WR 8.- Produit programme d'ordinateur, caractérisé en ce qu'il comporte des instructions propres à être stockées dans la mémoire d'un calculateur et exécutées par le processeur dudit calculateur pour mettre en oeuvre un procédé de séparation conforme à l'une quelconque des revendications précédentes.