FR 2848715 A1 20040618 FR 0215618 A 20021211 L'invention concerne un procédé de correction multi références des déformations spectrales de la voix introduites par un réseau de communication. Elle concerne également un système pour la mise en oeuvre du procédé. La présente invention a pour but l'amélioration de la qualité de la parole transmise sur les réseaux de communication, en offrant des moyens pour corriger les déformations spectrales du signal de parole, 10 déformations provoquées par différents maillons de la chaîne de transmission du réseau. La description qui en est donnée dans la suite fait explicitement référence à la transmission de la parole sur lignes téléphoniques "classiques" (c'est-à-dire 15 filaires), mais elle s'applique également à tout type de réseau de communication (fixe, mobile ou autre) introduisant des déformations spectrales dans le signal, les paramètres pris comme référence pour qualifier le réseau devant être modifiés en fonction du 20 réseau. On va rappeler dans ce qui suit les différentes déformations rencontrées dans le cas du Réseau Téléphonique Commuté (RTC). 1.1. Les dégradations du timbre de la voix sur le 25 Réseau RTC: La figure 1 représente un schéma d'une liaison RTC. La parole émise par un locuteur est transmise par un terminal émetteur 10, transportée par la ligne d'abonné 20, subit une conversion 30 analogiquenumérique (loi 30 A), est transmise par le réseau numérique 40, subit une conversion numérique (loi A)- analogique 50, est transmise par la liaison d'abonné 60, passe par le terminal récepteur 70 pour être enfin reçue par le destinataire. Chaque locuteur est relié par une ligne analogique (paire torsadée) au central téléphonique le plus proche. Il s'agit d'une transmission analogique en bande de base référence 1 et 3 sur la figure 1. La 5 liaison entre les centraux, emprunte un réseau entièrement numérique 40. Le spectre de la voix est affecté par deux types de distorsions lors de la transmission analogique du signal en bande de base. Le premier type de distorsion est le filtrage 10 passe-bande des terminaux et des points d'accès à la partie numérique du réseau. Les caractéristiques typiques de ce filtrage sont décrites par 1'UIT-T sous le nom de "système de référence intermédiaire" (SRI) [UIT-T, Recommandation P.48, 1988]. Ces caractéristiques fréquentielles, issues de mesures réalisées dans les années 70, tendent cependant à devenir obsolètes. C'est pourquoi 1'UIT-T préconise depuis 1996 d'utiliser un SRI "modifié" [UIT-T, Recommandation P.830, 1996], dont la caractéristique 20 nominale est représentée sur la figure 2 pour la partie émission, et sur la figure 3 pour la partie réception. Entre 200 et 3400 Hz, la tolérance est de 2,5 dB; en dessous de 200 Hz, la décroissance de la caractéristique du système global doit être d'au moins 25 15 dB par octave. On appelle respectivement, selon la terminologie de l'UIT-T, "système émetteur" et "système récepteur" les parties émission et réception du SRI. La seconde distorsion affectant le spectre de la voix est l'atténuation des lignes d'abonné. Dans un 30 modèle simple de la ligne analogique locale [donné dans une Note Technique CNET NT/LAA/ELR/289 par Cadoret, 1983], on considère que celle-ci introduit un affaiblissement du signal dont la valeur en dB dépend de sa longueur et est proportionnelle à la racine carrée de la fréquence. L'affaiblissement est de 3 dB à 800 Hz pour une ligne moyenne (environ 2 km), de 9,5 dB à 800 Hz pour les lignes les plus longues (jusqu'à 10 km). Selon ce modèle, l'affaiblissement d'une ligne, représenté sur la figure 4, a pour expression AdB (f) = AdB (800HZ). 800 (0.1) A ces distorsions s'ajoute le filtrage antirepliement du codeur MIC (réf 30). Celui-ci est 10 typiquement un filtre passe-bande 200-3400 Hz avec une réponse presque plate sur la bande passante et une forte atténuation en dehors de la bande, selon le gabarit de la figure 5 par exemple [National Semiconductor, aot 1994 Documentation technique 15 " TP3054,TP3057]. Au final, la voix subit une distorsion spectrale telle que représentée sur la figure 6 pour les différentes combinaisons de trois types de ligne analogique en émission et en réception (soit 6 20 distorsions), sous l'hypothèse d'équipements respectant la caractéristique nominale du SRI modifié. La voix apparaît ainsi étouffée si une des lignes analogiques est longue et souffre dans tous les cas d'un manque de "présence" d à l'affaiblissement des composantes basse 25 fréquence. 1.2. Les dégradations du timbre de la voix sur le Réseau Numérique à Intégration de Services (RNIS) et le réseau mobile GSM. Dans le RNIS et le réseau GSM, le signal est 30 numérisé dès le terminal. Les seules parties analogiques sont les transducteurs en émission et en réception associés à leurs chaînes d'amplification et de conditionnement respectives. L'UIT-T a défini des gabarits d'efficacité en fréquence à l'émission représentés sur la figure 7, et à la réception représenté sur la figure 8, valables à la fois pour les téléphones numériques filaires [UIT-T, Recommandation 5 P.310, mai 2000] et les terminaux numériques mobiles ou sans fil [UIT-T, Recommandation P.313, septembre 1999]. Par ailleurs, pour les réseaux GSM, il est reconnu que le codage et le décodage modifient légèrement l'enveloppe spectrale du signal. Cette altération est 10 représentée sur la figure 9 pour un bruit rose codé puis décodé en mode EFR (Enhanced Full Rate). L'effet de ces filtrages sur le timbre est principalement un affaiblissement des composantes basse fréquence, moins marqué cependant que dans le cas du 15 RTC. L'invention concerne la correction de ces distorsions spectrales par un traitement centralisé, c'est-à-dire un dispositif installé dans la partie numérique du réseau, comme indiqué sur la figure 10 20 pour le RTC. L'objectif d'une correction du timbre de la voix est que le timbre de la voix en réception soit le plus proche possible de celui de la voix émise par le locuteur, que l'on appellera voix originale. 25 2. tat de la technique Une compensation des distorsions spectrales introduites dans le signal de parole par les divers éléments de la liaison téléphonique est permise à ce 30 jour par des dispositifs à base d'égalisation. Celle-ci peut être fixe ou s'adapter en fonction des conditions de transmission. 2.1 L'égalisation fixe Des dispositifs d'égalisation centralisée ont été proposés dans les brevets US 5333195 [Duane O. Bowker] et US 5471527 [Helena S. Ho]. Ces égaliseurs sont des filtres fixes qui restaurent le niveau des basses 5 fréquences atténuées par l'émetteur. Bowker propose par exemple un gain de 10 à 15 dB sur la bande 100-300 Hz. Ces méthodes présentent deux inconvénients: * L'égaliseur ne compense que le filtrage de l'émetteur, de sorte qu'à la réception, les composantes 10 basse- fréquence restent fortement affaiblies par le filtrage SRI de réception. * Cette égalisation fixe compense des conditions de transmission (ligne et système d'émission) moyennes. Si les conditions réelles sont trop différentes (par 15 exemple si les lignes analogiques sont longues) le dispositif ne corrige pas suffisamment le timbre, voire l'altère plus que la liaison sans égalisation. 2.2 L'égalisation adaptative. Le dispositif décrit dans le brevet US 5915235 20 [Andrew P De Jaco] vise à corriger la réponse fréquentielle non idéale d'un transducteur de téléphone mobile. L'égaliseur est décrit comme étant placé entre le convertisseur analogique-numérique et le codeur CELP, mais peut être aussi bien dans le terminal que 25 dans le réseau. Le principe de l'égalisation est de rapprocher le spectre du signal reçu d'un spectre idéal. Deux méthodes sont proposées. La première méthode (illustrée par la figure 4 du brevet précité de De Jaco), consiste à calculer des 30 coefficients d'autocorrélation à long terme RLT RLT (ni) = G RLT(n- 1i) + (1-a)R(ni), (0.2) avec RLT(ni) i'è" coefficient d'autocorrélation à long terme à la nième trame, R(n,i) ieème coefficient d'autocorrélation spécifique à la nième trame, et ax constante de lissage fixée par exemple à 0,995. De ces coefficients sont dérivés les coefficients LPC à long terme, qui sont les coefficients d'un filtre blanchisseur. la sortie de ce filtre, le signal est 5 filtré par un filtre fixe qui lui imprime les caractéristiques spectrales à long terme idéales, ie celles qu'il aurait à la sortie d'un transducteur ayant la réponse fréquentielle idéale. Ces deux filtres sont complétés par un gain multiplicatif égal au rapport 10 entre les énergies à long terme de l'entrée du blanchisseur et de la sortie du deuxième filtre. La deuxième méthode, illustrée par la figure 5 du brevet précité de De Jaco, consiste à diviser le signal en sous-bandes, et, pour chaque sousbande, appliquer 15 un gain multiplicatif de manière à atteindre une énergie cible, ce gain étant défini comme le rapport entre l'énergie cible de la sous-bande et l'énergie à long terme (obtenue par un lissage de l'énergie instantanée) du signal dans cette sous-bande. Ces deux méthodes présentent l'inconvénient de ne corriger que la réponse non idéale du système d'émission, et pas celle du système de réception. Le dispositif du brevet US 5905969 [Chafik Mokbel] a pour objet de compenser le filtrage du système 25 d'émission et de la ligne d'abonné pour améliorer la reconnaissance centralisée de la parole et/ou la qualité de la parole transmise. Comme le présente la figure 3a de Mokbel, le spectre du signal est divisé en 24 sous-bandes, et chaque énergie de sous-bande est 30 multipliée par un gain adaptatif. L'adaptation du gain est réalisée selon l'algorithme du gradient stochastique, par minimisation de l'erreur quadratique, l'erreur étant définie comme la différence entre l'énergie de sous-bande et une énergie de référence définie pour chaque sous-bande. L'énergie de référence est modulée à chaque trame par l'énergie de la trame courante, de manière à respecter les variations naturelles de niveau à court terme du signal de parole. 5 La convergence de l'algorithme permet d'obtenir en sortie les 24 signaux de sous-bande égalisés. Si l'application visée est l'amélioration de la qualité vocale, le signal de parole égalisé est obtenu par transformée de Fourier inverse des énergies de 10 sous-bande égalisées. Le brevet de Mokbel, ne mentionne pas de résultats en termes d'amélioration de la qualité vocale, et reconnaît que la méthode est sousoptimale, en ce qu'elle réalise une convolution circulaire. Par 15 ailleurs, il est douteux qu'un signal de parole puisse être reconstruit correctement par transformée de Fourier inverse d'énergies de bandes distribuées selon l'échelle MEL. Enfin, le dispositif décrit ne corrige pas le filtrage du système de réception et de la ligne 20 analogique de réception. La compensation de l'effet de ligne est réalisée dans la méthode " Mokbel ", de soustraction cepstrale, dans le but d'améliorer la robustesse de la reconnaissance de la parole. On montre que le cepstre 25 du canal de transmission peut être estimé par le cepstre moyen du signal reçu, celuici étant préalablement blanchi par un filtre de préaccentuation. Cette méthode permet une nette amélioration des performances des systèmes de reconnaissance, mais est 30 considérée comme une méthode "off-line", 2 à 4 s étant nécessaires pour estimer le cepstre moyen. 2.3 Un autre état de la technique combine une préégalisation fixe et une égalisation adaptée et a fait l'objet du dépôt d'une demande de brevet FR 2822999 par le Déposant. Le dispositif décrit vise à corriger le timbre de la parole en combinant deux filtres. Un filtre fixe, appelé pré-égaliseur, compense les distorsions d'une liaison téléphonique moyenne, définie 5 comme étant constituée de deux lignes d'abonné moyennes et de systèmes d'émission et de réception respectant les réponses fréquentielles nominales définies dans [UIT-T, Recommandation P.48,App.I,1988]. Sa réponse fréquentielle, sur la bande [Fc-3150 Hz], est l'inverse 10 de la réponse globale de la partie analogique de cette liaison moyenne, Fc étant la fréquence basse limite d'égalisation. Cette pré-égalisation est complétée par un égaliseur adapté, qui adapte la correction de manière 15 plus précise aux conditions réelles de transmission. La réponse fréquentielle de l'égaliseur adapté est donnée par: I EQ(ft) I YJD=tdref (f) (0.3) EQ(f)[ S_ RX(f).L _RX(f)j 7(f) ' avec L RX la réponse fréquentielle de la ligne de 20 réception, SRX la réponse fréquentielle du système de réception et Yx(f) le spectre à long terme de la sortie x du pré-égaliseur. Le spectre à long terme est défini comme la moyenne temporelle des spectres à court terme des trames 25 successives de signal; Yref(f), appelé spectre de référence, est le spectre moyen de la parole défini par 'UIT [UIT-T/P.50/App. I, 1998], pris comme approximation du spectre à long terme original du locuteur. Du fait de cette approximation, la réponse 30 fréquentielle *de l'égaliseur adapté est très irrégulière et seule sa forme générale est pertinente. C'est pourquoi elle doit être lissée. L'égaliseur adapté étant réalisé sous la forme d'un filtre temporel RIF, ce lissage dans le domaine fréquentiel est obtenu par un étroit fenêtrage (symétrique) de la réponse impulsionnelle. Cette méthode permet de restaurer un timbre proche 5 de celui du signal original sur la bande d'égalisation [Fc-3150 Hz], mais: pour certains locuteurs, l'approximation de leur spectre à long terme original par le spectre de référence est très grossière, de sorte que l'égaliseur 10 introduit une distorsion perceptible; - le fort lissage de la réponse fréquentielle de l'égaliseur, rendu nécessaire par l'erreur d'approximation, interdit de corriger des distorsions spectrales fines. L'invention a pour but de remédier aux inconvénients de l'état de la technique. Elle a pour objet un procédé et un système pour améliorer la 20 correction du timbre en réduisant l'erreur d'approximation du spectre à long terme original des locuteurs. A cet effet, il est proposé de classer les locuteurs selon leur spectre à long terme et 25 d'approcher celui-ci non plus par un spectre de référence unique mais par un spectre de référence par classe. Le procédé proposé permet de réaliser un traitement d'égalisation apte à déterminer la classe du locuteur et à égaliser suivant le spectre de référence 30 de la classe. Cette réduction de l'erreur d'approximation permet de lisser moins fortement la réponse fréquentielle de l'égaliseur adapté, le rendant apte à corriger des distorsions spectrales plus fines. La présente invention a plus particulièrement pour objet un procédé de correction des déformations spectrales de la voix, introduites par un réseau de communication, comprenant une opération d'égalisation 5 sur une bande de fréquence [Fl-F21, adaptée à la distorsion réelle de la chaîne de transmission, cette opération étant réalisée au moyen d'un filtre numérique ayant une réponse fréquentielle fonction du rapport entre un spectre de référence et un spectre 10 correspondant au spectre à long terme du signal voix des locuteurs, principalement caractérisé en ce qu'il comprend: * Préalablement à l'opération d'égalisation du signal voix d'un locuteur en communication: - La constitution de classes de locuteurs avec une référence de voix par classe, * Puis, pour un locuteur donné en communication - Le classement de ce locuteur c'est-à-dire son affectation à une classe à partir de critères de 20 classement prédéfinis pour lui faire correspondre une référence de voix qui est la plus proche de la sienne, - L'égalisation du signal numérisé de la voix du locuteur opérée avec comme spectre de référence, la référence de voix de la classe à laquelle ledit 25 locuteur a été affecté. Selon une autre caractéristique, la constitution de classes de locuteurs comprend: - Le choix d'un corpus de N locuteurs enregistrés dans des conditions non dégradées et la détermination 30 de leur spectre de fréquence à long terme, - La classification des locuteurs du corpus selon leur cepstre partiel c'est à dire le cepstre calculé à partir du spectre à long terme restreint à la bande d'égalisation [Fl-F2] et en appliquant un critère de classification prédéfini sur ces cepstres pour obtenir K classes, - Le calcul du spectre de référence associé à chaque classe de manière à obtenir une référence de voix correspondant à chacune des classes. Selon une autre caractéristique, le spectre de référence sur la bande de fréquences d'égalisation [FlF2], associé à chaque classe, est calculé par transformée de Fourier du centre de la classe défini 10 par son cesptre partiel. Selon une autre caractéristique, la classification d'un locuteur comprend: - L'utilisation du pitch moyen du signal voix et du cepstre partiel de ce signal comme paramètres de 15 classement, - L'application d'une fonction discriminante à ces paramètres pour classer ledit locuteur. Selon l'invention le procédé comprend en outre une étape de pré-égalisation du signal numérique par un 20 filtre fixe ayant une réponse fréquentielle dans la bande de fréquence [Fl-F21, correspondant à l'inverse d'une déformation spectrale de référence introduite par la liaison téléphonique. Selon une autre caractéristique, l'égalisation du 25 signal numérisé de la voix d'un locuteur comprend: - la détection d'une activité vocale sur la ligne pour déclencher un enchaînement de traitements comprenant le calcul du spectre à long terme, la classification du locuteur, le calcul du module de la 30 réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [Fl-F21 et le calcul des coefficients du filtre numérique différenciés selon la classe du locuteur, à partir de ce module, - la commande du filtre avec les coefficients obtenus, - le filtrage du signal sortant du pré-égaliseur par ledit filtre. Selon une autre caractéristique, le calcul du module [EQ] de la réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [F1F2] est réalisé par la mise en ouvre de la relation suivante: 1 EQ(f) J= Yie.(f) (0.3) SQ() _RX(f).L _ RX(f)l f) dans laquelle yref(f)est le spectre de référence de la classe à laquelle appartient ledit locuteur, et dans laquelle LRX est la réponse fréquentielle 15 de la ligne de réception, SRX la réponse fréquentielle du système de réception etx(f) le spectre à long terme du signal x d'entrée du filtre. Selon une variante, le calcul du module de la réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [F1-F2] est réalisé par la mise en oeuvre de la relation suivante: Ce = Crqef -Cx Cs- Cl r. (0.13) dans laquelle &eq, Cx, CI, et CPLr sont les cepstres partiels respectifs de l'égaliseur adapté, du signal d'entrée x du filtre égaliseur, du système de réception et de la ligne de réception, CPref étant le cepstre 30 partiel de référence, centre de la classe du locuteur. Le module [EQ] restreint à la bande F1-F2 est alors calculé par Transformée de Fourier discrète de CPEQ. L'invention a également pour objet un système de correction des déformations spectrales de la voix, introduites par un réseau de communication, comprenant 5 des moyens d'égalisation adaptée dans une bande de fréquence [Fl-F2] qui comportent un filtre numérique dont la réponse fréquentielle est fonction du rapport entre un spectre de référence et un spectre correspondant au spectre à long terme d'un signal voix, 10 principalement caractérisé en ce que ces moyens comprennent en outre: - des moyens de traitements du signal pour le calcul des coefficients du filtre numérique munis: * d'un bloc de traitement du signal pour 15 calculer le module de la réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [Fl-F2] selon la relation suivante 20} EQ(f) Yr=f(f) (0.3) ISRX(f).L RX(f) ,f dans laquelle Yref (f)est le spectre de référence, qui peut être différent d'un locuteur à un autre et qui correspond à une référence de classe prédéterminée à laquelle appartient ledit locuteur, 25 et dans laquelle LRX est la réponse fréquentielle de la ligne de réception, SRX la réponse fréquentielle du système de réception et yx(f) le spectre à long terme du signal x d'entrée du filtre; * d'un deuxième bloc de traitements pour le calcul de la réponse impulsionnelle à partir du module de réponse fréquentielle ainsi calculé, afin de déterminer les coefficients du filtre différenciés selon la classe du locuteur. Selon une autre caractéristique, le premier bloc de 5 traitement comprend des moyens de calcul du cepstre partiel du filtre égaliseur selon la relation: Cep, = Crpef - cxp CeP lx CIP rx S(0.13) dans laquelle Cq, Cx, C et Cd,, sont les cepstres partiels respectifs de l'égaliseur adapté, du signal 10 d'entrée x du filtre égaliseur, du système de réception et de la ligne de réception, CPref étant le cepstre partiel de référence, centre de la classe du locuteur, le module de [EQ] restreint à la bande F1-F2 est alors calculé par Transformée de Fourier Discrète de CFeq. Selon une autre caractéristique, le premier bloc de traitements comprend un sous-ensemble pour calculer les coefficients du cepstre partiel d'un locuteur en communication et un deuxième sous-ensemble pour opérer le classement de ce locuteur, ce deuxième sous-ensemble 20 comprenant un bloc de calcul du pitch F0, un bloc d'estimation du pitch moyen à partir du pitch calculé F0o, et un bloc de classement appliquant une fonction discriminante sur le vecteur x ayant pour composantes le pitch moyen et les coefficients du cepstre partiel 25 pour classer ledit locuteur. Selon l'invention, le système comprend en outre un pré-égaliseur, le signal égalisé à partir de spectres de référence différenciés selon la classe du locuteur étant le signal x de sortie du pré-égaliseur. 30 D'autres particularités et avantages de l'invention apparaîtront clairement dans la description suivante, qui est donnée à titre d'exemple illustratif et non limitatif et qui est faite en regard des figures annexées qui représentent: - La figure 1, une liaison téléphonique schématisée pour un réseau téléphonique commuté (RTC), - La figure 2, la courbe de réponse fréquentielle en émission du système de référence intermédiaire SRI modifié, - La figure 3, la courbe de réponse fréquentielle en réception du système de référence intermédiaire SRI 10 modifié, -La figure 4, la réponse en fréquence des lignes d'abonnés selon leur longueur, - La figure 5,le gabarit du filtre anti-repliement du codeur MIC, - La figure 6, les distorsions spectrales subies par la parole sur le réseau téléphonique commuté avec SRI moyen et différentes combinaisons de lignes analogiques, - La figure 7, le gabarit à l'émission pour les 20 terminaux numériques, - La figure 8, le gabarit à la réception pour les terminaux numériques, - La figure 9, la distorsion spectrale introduite par le codage-décodage GSM en mode EFR (Enhanced Full 25 Rate) - La figure 10, le schéma d'un réseau de communication avec un système de correction des distorsions de la parole, - La figure 11, les étapes de calcul du cepstre 30 partiel, - La figure 12, la classification des cepstres partiels selon le critère de la variance, - Les figures 13a et 13b, les spectres à long terme correspondant aux centres des classes de locuteurs respectivement pour hommes et femmes, - La figure 14, les caractéristiques fréquentielles 5 des filtrages appliqués au corpus pour définir le corpus d'apprentissage, - La figure 15, la réponse fréquentielle du préégaliseur pour différentes fréquences Fc, - La figure 16, le schéma de réalisation du système 10 de correction par égalisation différenciée par classe de locuteurs, - La figure 17, une variante d'exécution du système selon la figure 16. Dans toute la suite les mêmes références portées sur les dessins correspondent aux mêmes éléments. La description qui va suivre va d'abord présenter l'étape préalable de classification d'un corpus de locuteurs selon leur spectre à long terme. Cette étape 20 permet de définir K classes et une référence par classe. Un enchaînement de traitements permet de traiter le signal de parole (dès la détection d'une activité vocale par le système) de chaque locuteur pour d'une 25 part classer les locuteurs c'est à dire les affecter à une classe selon des critères prédéterminés et pour d'autre part corriger la voix en utilisant la référence de la classe du locuteur. Etape préalable de classification des locuteurs. * Choix du corpus de définition des classes. Le spectre de référence étant une approximation du spectre à long terme original des locuteurs, la définition des classes de locuteurs et de leurs spectres de référence respectifs nécessite de disposer d'un corpus de locuteurs enregistrés dans des conditions non dégradées. En particulier, le spectre à long terme d'un locuteur mesuré sur cet enregistrement 5 doit pouvoir être considéré comme son spectre original, i.e. celui de sa voix à l'extrémité d'émission d'une liaison téléphonique. Définition de l'individu: le cepstre partiel. Le traitement proposé permet de disposer, dans 10 chaque classe, d'un spectre de référence le plus proche possible du spectre à long terme de chaque membre de la classe. Cependant, seule la partie du spectre comprise dans la bande d'égalisation Fl-F2 est prise en compte dans le traitement d'égalisation adaptée. Les classes 15 sont donc constituées selon le spectre à long terme restreint à cette bande. D'autre part, la comparaison entre deux spectres est effectuée à un faible niveau de résolution spectrale, de manière à ne refléter que l'enveloppe 20 spectrale. C'est pourquoi, on se place de préférence, dans l'espace des premiers coefficients cepstraux d'ordre supérieur à O (le coefficient d'ordre O représentant l'énergie), le choix du nombre de coefficients dépendant de la résolution spectrale 25 souhaitée. On détermine donc dans le traitement, le 'cepstre partiel à long terme", que l'on note Cp, comme la représentation cepstrale du spectre à long terme restreint à une bande de fréquence. Si l'on note kl et 30 k2 les indices de fréquence correspondant respectivement aux fréquences Fl et F2, et y le spectre à long terme de la parole, le cepstre partiel est défini par la relation: CP =TFD-' (lOlog(y (k1...k2)oy(k2 -1... k +1))) (0.4) o - désigne l'opération de concaténation. La (TFD) Transformée de Fourier Discrète Inverse est calculée par exemple par IFFT après interpolation des échantillons du spectre tronqué de manière à 5 atteindre un nombre d'échantillons puissance de 2. Par exemple, en choisissant la bande d'égalisation 187-3187 Hz, correspondant aux indices fréquentiels 5 à 101 pour une représentation du spectre (symétrisé) sur 256 points (de C à 255), l'interpolation se fait simplement 10 en intercalant une raie fréquentielle (interpolée linéairement) toutes les trois raies dans le spectre restreint à 187-3187 Hz. Les étapes du calcul du cepstre partiel sont représentées sur la figure 11. Pour que les coefficients cepstraux reflètent l'enveloppe spectrale mais pas l'influence de la structure harmonique du spectre de la parole sur les spectres à long terme, on ne conserve pas les coefficients d'ordre élevé. Les locuteurs à classer 20 sont donc représentés par les coefficients d'ordres 1 à L de leur cepstre partiel à long terme, L valant typiquement 20. * La Classification. Les classes sont constituées par exemple de manière 25 non-supervisée, selon une classification hiérarchique ascendante. Celle-ci consiste à créer, à partir de N individus disjoints, une hiérarchie de partitions selon le processus suivant: à chaque étape, on agrège les deux 30 éléments les plus proches, un élément étant soit un individu non agrégé, soit un agrégat d'individus constitué lors d'une précédente étape. La proximité entre deux éléments est déterminée par une mesure de dissimilarité que l'on appelle distance. Le processus se poursuit jusqu'à l'agrégation de toute la population. La hiérarchie de partitions ainsi créée peut se représenter sous la forme d'un arbre comme celui de la figure 12, contenant N-1 partitions 5 imbriquées. Chaque coupure de l'arbre fournit une partition, d'autant plus fine que l'on coupe bas. Dans ce type de classification, on choisit comme mesure de distance entre deux éléments, la variation d'inertie intra-classes résultant de leur agrégation. 10 Une partition est en effet d'autant meilleure que les classes créées sont homogènes, c'est-à-dire que l'inertie intra-classes est faible. Dans le cas d'un nuage de points xi de masses respectives mi, répartis en classes q de centres de gravité respectifs gq, 15 l'inertie intra-classes est définie par intra = Y, mi|xi -gj|l (0.5) q ieq L'inertie intra-classes, nulle à l'étape initiale de l'algorithme de calcul, s'accroît inévitablement à chaque agrégation. On utilise de préférence le principe connu de l'agrégation selon la variance. Selon ce principe on recherche, à chaque étape de l'algorithme mis en oeuvre, les deux éléments dont l'agrégation produit l'augmentation d'inertie intra-classes la plus faible. La partition ainsi obtenue est améliorée par une procédure d'agrégation autour des centres mobiles, qui permet de réduire la variance intraclasses. Le spectre de référence, sur la bande Fl-F2, associé à chaque classe est calculé par transformée de 30 Fourier du centre de la classe. * Exemple de classification Le traitement décrit ci-dessus est appliqué à un corpus de 63 locuteurs. L'arbre de classification du corpus est représenté sur la figure 12. Dans cette représentation, la hauteur d'un segment horizontal agrégeant deux éléments est choisie proportionnelle à leur distance, ce qui permet de visualiser la proximité 5 des éléments regroupés dans une même classe. Cette représentation facilite le choix du niveau de coupure de l'arbre, et donc des classes retenues. La coupure doit être faite au-dessus des agrégations de niveau faible, qui regroupent des individus proches, et en 10 dessous des agrégations de niveau élevé, qui associent des groupes d'individus bien distincts. De cette manière, on obtient de façon nette quatre classes (K= 4). Ces classes sont assez homogènes du point de vue du sexe des locuteurs, et une coupure de 15 l'arbre en deux classes fait apparaître à peu près une classe hommes et une classe femmes. La consolidation de cette partition par une procédure d'agrégation autour des centres mobiles aboutit à quatre classes de cardinaux 11, 18, 18 et 16, 20 plus homogènes que précédemment du point de vue du sexe : seuls un homme et deux femmes sont affectés à des classes ne correspondant par à leur sexe. Les spectres restreints à la bande 187-3187 Hz correspondant aux centres de ces classes sont 25 représentés sur les figures 13a et 13b pour les classes hommes et femmes ainsi que pour leurs sous-classes respectives. Ces spectres, résultats de la classification, sont utilisés comme référence multiple par l'égaliseur adapté. * Utilisation de critères de classement des locuteurs. Les classes de locuteurs étant définies, le traitement prévoit l'utilisation de paramètres et des critères pour affecter un locuteur à l'une ou l'autre des classes. Cette affectation n'est pas réalisée simplement selon la proximité du cepstre partiel avec un des 5 centres de classes, puisque ce cepstre est dévié par la partie de la liaison téléphonique en amont de l'égaliseur. Il est proposé avantageusement d'utiliser des critères de classement robustes à cette déviation. 10 Cette robustesse est assurée à la fois par le choix des paramètres de classement et par celui du corpus d'apprentissage des critères de classement. * De préférence on utilise les paramètres de classement: pitch moyen et cepstre partiel. Les classes précédemment définies sont homogènes du point de vue du sexe. Le pitch moyen étant à la fois assez discriminant pour un classement homme / femme et insensible aux distorsions spectrales induites par une liaison téléphonique, il est donc utilisé comme 20 paramètre de classement, conjointement avec le cepstre partiel. * Choix du corpus d'apprentissage des critères de classement. On applique à ces paramètres une technique de 25 discrimination par exemple la technique usuelle de l'analyse linéaire discriminante. D'autres techniques connues peuvent être utilisées telles qu'une technique non linéaire utilisant un réseau de neurones. Si l'on dispose de N individus décrits par des vecteurs de dimension p et répartis a priori en K classes, l'analyse linéaire discriminante consiste: - dans un premier temps, à chercher les K-l fonctions linéaires indépendantes qui séparent au mieux les K classes. Il s'agit de déterminer quelles sont les combinaisons linéaires des p composantes des vecteurs qui minimisent la variance intra-classes et maximisent la variance interclasses. - dans un deuxième temps, à déterminer la classe d'un nouvel individu par application des fonctions linéaires discriminantes au vecteur le représentant. Dans le cas présent, les vecteurs représentatifs des individus ont pour composantes le pitch et les 10 coefficients 1 à L (typiquement, L = 20) du cepstre partiel. La robustesse des fonctions discriminantes à la déviation des coefficients cepstraux est assurée à la fois par la présence du pitch dans les paramètres et par le choix du corpus d'apprentissage. Celui-ci est 15 composé d'individus dont la voix originale a subi une grande diversité de filtrages représentatifs des distorsions occasionnées par les liaisons téléphoniques. Plus précisément, à partir d'un corpus de voix 20 originales (non dégradées) de N locuteurs, on définit un corpus de N vecteurs de composantes [F;CP(l); ... ; C(L)], avec F. le pitch moyen et Ca le cepstre partiel. La construction du corpus d'apprentissage desdites fonctions consiste à définir un ensemble de M biais 25 cepstraux qui s'ajouteront chacun à chaque cepstre partiel représentatif d'un locuteur du corpus original, ce qui permet d'obtenir un nouveau corpus de NM individus. Ces biais dans le domaine du cepstre partiel 30 correspondent à une large gamme de distorsions spectrales sur la bande Fl-F2, proches de celles pouvant résulter de la liaison téléphonique. A titre d'exemple, on propose l'ensemble de réponses fréquentielles représentées sur la figure 14 pour la bande 187-3187 Hz: chaque réponse fréquentielle correspond à un chemin de gauche à droite dans le treillis. L'amplitude de leurs variations sur cette bande n'excède pas 20 dB, à l'instar des 5 caractéristiques extrémales des systèmes d'émission et lignes. A partir de ces 81 caractéristiques fréquentielles sont calculés les 81 biais correspondants dans le domaine du cepstre partiel, selon le traitement décrit 10 pour la mise en oeuvre de la relation (0.4). Par addition de ces biais au corpus de 63 locuteurs précédemment utilisé, on obtient un corpus d'apprentissage comptant 5103 individus représentatifs de diverses conditions (locuteur, filtrage de la 15 liaison). Dans le cas de classement par analyse linéaire discriminante: * Application des critères de classement. Soit (ak)l La figure 16 illustre le système de correction apte à mettre en oeuvre le procédé. La figure 17 illustre ce système selon une variante de réalisation comme cela va être détaillé dans la suite. Ces variantes portent sur le mode de calcul du module de la réponse fréquentielle 20 de l'égaliseur adapté restreinte à la bande F1-F2. Le pré-égaliseur 200 est un filtre fixe, dont la réponse fréquentielle, sur la bande F1-F2, est l'inverse de la réponse globale de la partie analogique d'une liaison moyenne telle que définie précédemment 25 [UITT/P.830, 1996]. La raideur de la réponse fréquentielle de ce filtre implique une réponse impulsionnelle longue; c'est pourquoi, de manière à limiter le retard introduit par le traitement, le pré-égaliseur est réalisé typiquement 30 sous forme d'un filtre RII, d'ordre 20 par exemple. La figure 15 représente les réponses fréquentielles typiques du préégaliseur pour trois valeurs de Fl. La dispersion des retards de groupe est inférieure à 2 ms, de sorte que la distorsion de phase résultante n'est pas perceptible. La chaîne de traitement 400 qui suit permet le classement du locuteur et l'égalisation adaptée différenciée. Cette chaîne comprend deux blocs de traitements 400A et 400B. Le bloc 400A permet de 10 calculer le module de la réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation EQ dB [F1-F2]. Le deuxième bloc 400B permet de calculer la réponse impulsionnelle du filtre égaliseur afin d'obtenir les 15 coefficients eq(n) du filtre différenciés selon la classe du locuteur. Un détecteur de trames d'activité vocale 401 permet de déclencher les différents traitements. Le bloc de traitement 410 permet le classement du 20 locuteur. Le bloc de traitement 420 permet de calculer le spectre à long terme suivi du calcul du cepstre partiel de ce locuteur. La sortie de ces deux blocs est appliquée sur 25 l'opérateur 428a ou 428b. La sortie de cet opérateur fournit le module de la réponse fréquentielle de l'égaliseur adapté en dB restreinte à la bande d'égalisation F1-F2 via le bloc 429 pour 428a, via le bloc 440 pour 428b. Les blocs de traitement 430 à 435 permettent de calculer les coefficients eq(n) du filtre. La sortie x(n) du pré-égaliseur est analysée par trames successives d'une durée typique de 32 ms, avec un recouvrement inter-trames de 50 % typiquement. On ouvre pour cela une fenêtre d'analyse représentée par les blocs 402 et 403. L'opération d'égalisation adaptée est mise en oeuvre par un filtre RIF 300 dont les coefficients sont 5 calculés à chaque trame d'activité vocale par la chaîne de traitement illustrée sur les figures 16 et 17. Le calcul de ces coefficients correspond au calcul de la réponse impulsionnelle du filtre à partir du module de la réponse fréquentielle. Le spectre à long terme de x(n), Yx, est d'abord calculé (à partir de l'instant initial de fonctionnement) sur une fenêtre temporelle croissant de 0 à une durée T d'activité vocale (typiquement 4 s), puis ajusté récursivement à chaque trame d'activité 15 vocale, ce qui se traduit par la formule générique suivante: yr(f,n) =a(n)X(f,n)2 +(1-a(n))y,(f,n-1), (0.10) o yx (f,n) est le spectre à long terme de x à la nième trame d'activité vocale, X(f,n) la transformée de 20 Fourier de la nième trame d'activité vocale, et a(n) est défini par la relation (0.11). En notant N le nombre de trames dans la durée T, a(n)= (0. 11) (nmin(n,N) (0.11) Ce calcul est opéré par les blocs 421, 422, 423. On calcule ensuite, à partir de ce spectre à long terme le cepstre partiel Cp, selon la relation (0.4) mise en oeuvre par les blocs de traitement 424, 425, 426. Le pitch moyen F0 est estimé par le bloc de traitement 412 à chaque trame voisée selon la formule: Fo (m) =a(m)Fo(m)+ (1 -a(m))Fo(m-1) (0.12) o F0(m) est le pitch de la mième trame voisée et est calculé par le bloc 411 selon une méthode appropriée de l'état de l'art (par exemple méthode de 5 l'autocorrélation, avec détermination du voisement par comparaison de l'autocorrélation normalisée à un seuil [UIT-T/G.729, 1996]). Ainsi, à chaque trame d'activité vocale, on dispose d'un nouveau vecteur x de composantes le pitch 10 moyen et les coefficients 1 à L du cepstre partiel, auquel on applique la fonction discriminante a définie à partir du corpus d'apprentissage. Ce traitement est mis en oeuvre par le bloc 413. Le locuteur est alors affecté à la classe q de score discriminant minimal. Le module en dB de la réponse fréquentielle de l'égaliseur adapté restreinte à la bande F1-F2, noté IEQIdB[Fl-F2], est calculé selon l'une des deux méthodes suivantes: La première méthode (figure 16) consiste à calculer 20 IEQIF1-F2 selon l'équation (0.3), o YrefD est le spectre de référence de la classe du locuteur (transformée de Fourier du centre de la classe). Cette méthode de calcul est mise en oeuvre dans cette variante représentée sur la figure 16 avec les opérateurs 414a, 25 428a, 427 et 429. La deuxième méthode (figure 17) consiste à transcrire l'équation (0.3) dans le domaine du cepstre partiel, puisque l'on dispose du cepstre partiel de la sortie x du pré-égaliseur, nécessaire au classement du 30 locuteur. Ainsi, l'équation (0.3) devient: cp- =cP - cP - cP -cp (1 eq ref x s _rx l_rx' (0.13) O Ceq, Cx, s et Ca sont les cepstres partiels respectifs de l'égaliseur adapté, de la sortie x du pré-égaliseur, du système de réception et de la ligne de réception, CPref étant le cepstre partiel de 5 référence, centre de la classe du locuteur. Les cepstres partiels sont calculés comme indiqué précédemment, en sélectionnant la bande de fréquences F1-F2. Ce calcul est effectué uniquement pour les coefficients 1 à 20, les coefficients suivants étant 10 inutiles car représentatifs d'une finesse spectrale qui sera éliminée par la suite. Les 20 coefficients du cepstre partiel de l'égaliseur adapté sont obtenus par les opérateurs 414b et 428b selon la relation (0.13). Le bloc de traitement 441 complète ces 20 coefficients par des zéros, les symétrise et calcule, à partir du vecteur ainsi formé, le module en dB de la réponse fréquentielle de l'égaliseur adapté restreinte à la bande F1F2 en mettant en ouvre la relation 20 suivante: EQdBIFF =TFD- (C). (0.14) Cette réponse est décimée d'un facteur 34 par 25 l'opérateur 442. Pour les deux variantes qui viennent d'être décrites, les valeurs de IEQI hors de la bande F1-F2 sont calculées par extrapolation linéaire de la valeur en dB de IEQI F1-F2, notée EQdB par la suite, par le bloc 30 430 et de la manière suivante: Pour chaque indice de fréquence k, l'approximation linéaire de EQdB s'exprime par: EQdB (k) = a, + a2k (0.15) Les coefficients ai et a2 sont choisis de manière à minimiser l'erreur quadratique de l'approximation sur l'intervalle F1-F2, définie par k2 e = EQdB (k)-EQd (k) (0.16) k=k, Les coefficients al et a2 sont donc définis par: k2-kl+1 t k l EQdB (k) ak=k (0.17) =al k=kl k =k] katk 9, k È kEQdB (k) k=kl k=k, k=k, Les valeurs de IEQI, en dB, hors de la bande F1-F2, sont alors calculées à partir de la formule (0.15). La caractéristique en fréquence ainsi obtenue doit 10 être lissée. Le filtrage étant réalisé dans le domaine temporel, le moyen permettant ce lissage est de multiplier par une fenêtre étroite la réponse impulsionnelle correspondante. La réponse impulsionnelle est obtenue par une 15 opération IFFT appliquée sur IEQI réalisée par les blocs 431 et 432 suivie d'une symétrisation réalisée par le bloc de traitement 433, de manière à obtenir un filtre causal à phase linéaire. La réponse impulsionnelle résultante est multipliée, opérateur 20 435, par une fenêtre temporelle 434. La fenêtre utilisée est typiquement une fenêtre de Hamming de longueur 31 centrée sur le pic de la réponse impulsionnelle et est appliquée sur la réponse impulsionnelle au moyen de l'opérateur 435. 25 REVENDICATIONS 1. Procédé de correction des déformations spectrales de la voix, introduites par un réseau de communication, comprenant une opération d'égalisation 5 sur une bande de fréquence [F1- F2], adaptée à la distorsion réelle de la chaîne de transmission, cette opération étant réalisée au moyen d'un filtre numérique ayant une réponse fréquentielle fonction du rapport entre un spectre de référence et un spectre 10 correspondant au spectre à long terme du signal voix des locuteurs, caractérisé en ce qu'il comprend: * Préalablement à l'opération d'égalisation du signal voix d'un locuteur en communication: - La constitution de classes de locuteurs avec une 15 référence de voix par classe, * Puis, pour un locuteur en communication donné - Le classement de ce locuteur c'est-à-dire son affectation à une classe à partir de critères de classement prédéfinis pour lui faire correspondre une 20 référence de voix qui est la plus proche de la sienne, - L'égalisation du signal numérisé de la voix du locuteur opérée avec comme spectre de référence, la référence de voix de la classe à laquelle ledit locuteur a été affecté. 2. Procédé de correction des déformations spectrales de la voix selon la revendication 1, caractérisé en ce que: * La constitution de classes de locuteurs comprend: - Le choix d'un corpus de N locuteurs enregistrés dans des conditions non dégradées et la détermination de leur spectre de fréquence à long terme, - La classification des locuteurs du corpus selon 5 leur cepstre partiel c'est à dire le cepstre calculé à partir du spectre à long terme restreint à la bande d'égalisation [F1-F2] et en appliquant un critère de classification prédéfini sur ces cepstres pour obtenir K classes, - Le calcul du spectre de référence associé à chaque classe de manière à obtenir une référence de voix correspondant à chacune des classes. 3. Procédé de correction des déformations 15 spectrales de la voix selon la revendication 2, caractérisé en ce que le spectre de référence sur la bande de fréquences d'égalisation [F1-F2], associé à chaque classe, est calculé par transformée de Fourier du centre de la classe défini par son cesptre partiel. 20 4. Procédé de correction des déformations spectrales de la voix selon la revendication 1, caractérisé en ce que: * La classification d'un locuteur comprend - L'utilisation du pitch moyen du signal voix et du cepstre partiel de ce signal comme paramètres de classement, - L'application d'une fonction discriminante à ces paramètres pour classer ledit locuteur. 30 5. Procédé de correction des déformations spectrales de la voix selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre une étape de préégalisation du signal numérique par un filtre fixe ayant une réponse fréquentielle dans la bande de fréquence [Fl-F2], correspondant à l'inverse d'une déformation spectrale de référence introduite par la liaison téléphonique. 6. Procédé de correction des déformations spectrales de la voix selon l'une quelconque des revendications précédentes, caractérisé en ce que l'égalisation du signal numérisé de la voix d'un 10 locuteur comprend: la détection d'une activité vocale sur la ligne pour déclencher un enchaînement de traitements comprenant le calcul du spectre à long terme, la classification du locuteur, le calcul du module de la 15 réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [Fl-F2] et le calcul des coefficients du filtre numérique différenciés selon la classe du locuteur, à partir de ce module, - la commande du filtre avec les coefficients 20 obtenus, - le filtrage du signal sortant du pré-égaliseur par ledit filtre. 7. Procédé de correction des déformations 25 spectrales de la voix selon la revendication 6, caractérisé en ce que le calcul du module de la réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [Fl-F2] est réalisé par la mise en oeuvre de la relation suivante 1 Yref () (0.3) _ EQ()IS_RX(f).LRX(f)I y̋ (f) dans laquelle yref(f)est le spectre de référence de la classe à laquelle appartient ledit locuteur, et dans laquelle L RX est la réponse fréquentielle de la ligne de réception, SRX la réponse fréquentielle du système de réception et yx(f) le spectre à long terme du signal x d'entrée du filtre. 8. Procédé de correction des déformations spectrales de la voix selon la revendication 6, caractérisé en ce que le calcul du module [EQ] de la 10 réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [F1-F2] est réalisé par la mise en oeuvre de la relation suivante: CeP = Cr'ef - cx, Cs'-rx -ci,,., (0.13) dans laquelle Ceq, CP, CPsr et Cdzx sont les cepstres partiels respectifs de l'égaliseur adapté, du signal d'entrée x du filtre égaliseur, du système de réception et de la ligne de réception, CPref étant le cepstre partiel de référence, centre de la 20 classe du locuteur; le module [EQ] restreint à la bande F1-F2 étant calculé par Transformée de Fourier Discrète de CPeq.. 9. Système de correction des déformations 25 spectrales de la voix, introduites par un réseau de communication, comprenant des moyens d'égalisation adaptée dans une bande de fréquence [F1-F2] qui comportent un filtre numérique (300) dont la réponse fréquentielle est fonction du rapport entre un spectre 30 de référence et un spectre correspondant au spectre à long terme d'un signal voix, caractérisé en ce que ces moyens comprennent en outre: - des moyens de traitement du signal (400) pour le calcul des coefficients du filtre numérique munis: * d'un premier bloc de traitements (400A) du signal pour calculer le module de la réponse fréquentielle du filtre égaliseur restreinte à la bande d'égalisation [F1-F2] selon la relation suivante: I EQ(f) 1= SgX(f) (0.3) iS RX(f).L _ RX( ( dans laquelle Yref (f)est le spectre de référence, qui peut être différent d'un locuteur à un autre et 10 qui correspond à une référence de classe prédéterminée à laquelle appartient ledit locuteur, et dans laquelle L RX est la réponse fréquentielle de la ligne de réception, SRX la réponse fréquentielle du système de réception et YX(f) le 15 spectre à long terme du signal x d'entrée du filtre; * d'un deuxième bloc de traitements (400B) pour le calcul de la réponse impulsionnelle à partir du module de réponse fréquentielle ainsi 20 calculé, afin de déterminer les coefficients du filtre différenciés selon la classe du locuteur 10. Système de correction des déformations spectrales de la voix selon la revendication 9, 25 caractérisé en ce que le premier bloc de traitements (400A) comprend des moyens (414b, 428b) de calcul du cepstre partiel du filtre égaliseur selon la relation: Ce^ = Cref -Cx C_- lx Clx 1 (0.13) dans laquelle Ceq, Cx, Cs_ et CaP sont les cepstres 30 partiels respectifs de l'égaliseur adapté, du signal d'entrée x du filtre égaliseur, du système de réception et de la ligne de réception, CPref étant le cepstre partiel de référence, centre de la classe du locuteur le module [EQI restreint à la bande F1-F2 étant calculé par Transformée de Fourier Discrète de CPeq.. 11. Système de correction des déformations spectrales de la voix selon la revendication 9 ou 10, caractérisé en ce que le premier bloc de traitements comprend un sous-ensemble (420) pour calculer les coefficients du cepstre partiel d'un locuteur en 10 communication et un deuxième sous-ensemble (410) pour opérer le classement de ce locuteur, ce deuxième sous ensemble comprenant un bloc (411) de calcul du pitch Fo, un bloc (412) d'estimation du pitch moyen à partir du pitch calculé Fo, et un bloc (413) de classement 15 appliquant une fonction discriminante sur le vecteur x ayant pour composantes le pitch moyen et les coefficients du cepstre partiel pour classer ledit locuteur. 12. Système de correction des déformations spectrales de la voix selon l'une quelconque des revendications 9 à 11, caractérisé en ce qu'il comprend un pré-égaliseur (200) et en ce que le signal égalisé à partir de spectres de référence différenciés selon la 25 classe du locuteur est le signal x de sortie du préégaliseur.