Brevet FR2674710A1 — publié le 19921002

PROCEDE ET SYSTEME DE TRAITEMENT DES PREECHOS D'UN SIGNAL AUDIO-NUMERIQUE CODE PAR TRANSFORMEE FREQUENTIELLE L'invention est relative à un procédé et à un sys- tème de traitement des prééchos d'un signal audio-numérique codé par transformée fréquentielle. Les processus de stockage, de transmission ou de traitement des signaux sonores numérisés, de plus en plus utilisés de nos jours, se heurtent au problème du débit numérique nécessaire toujours plus important selon les applications envisagées. Des techniques très élaborées, telles que le codage de ces signaux par transformée fréquentielle, telle que transformée en cosinus discrète modifiée, TCDM, sont utilisées afin de diminuer le débit d'information résultant, tout en conservant sensiblement la qualité du signal d'origine Ainsi qu'on l'a représenté de manière illustrative en figure la, le codage par transformée20 procède par découpage du signal temporel, succession d'échantillons x(n), en blocs comportant N échantillons. Sur chacun des blocs précités, après application d'une fenêtre de pondération spatio-temporelle d'amplitude, telle que la fenêtre du bloc m sur les échantillons O à N-1 ainsi25 que représenté en figure la, ces fenêtres étant réitérées avec un recouvrement sur N/2 échantillons, une transformée temps/fréquence est appliquée au bloc d'échantillons ainsi pondérés, pour obtenir un ensemble de coefficients y(k) lesquels sont alors codés puis transmis avec la réduction30 de débit recherchée. Au décodage, une transformation inverse est appliquée aux coefficients reçus y'(k), ce qui permet de régénérer les échantillons d'origine x'(n) en temps réel, au temps de traitement codage-décodage, temps de35 transmission près On notera toutefois que plusieurs types de transformées fréquentielles peuvent être utilisées, la transformée TCDM permettant cependant d'obtenir la réduction de débit la plus importante Pour une description plus détaillée de ce type de traitement, on pourra se reporter par exemple aux publications de Johnston "Transfer Coding of Audio Singals Using Perceptual Noise Criteria" gi IEEE Journal on selected Areas in Communications, Vol 6, No 2, Février 1988, p 314-323 et de Princen et Bradley "Adaptive Transform Coding incorporating Time Domain Aliasing Cancellation" Speech communication, Décembre 1987. Dans les systèmes de codage par transformée 1 i O fréquentielle, le bruit de codage résultant de la quantification des coefficients est réparti uniformément sur toute la durée du bloc d'échantillons. Lorsque ce bloc contient une non stationnarité, telle qu'une attaque brusque, ainsi que représenté en l 5 figure lb, le spectre fréquentiel du signal, c'est-à-dire la distribution de l'amplitude des coefficients y(k) obtenus après transformée fréquentielle, en fonction de la fréquence f est à peu près plat. En outre, les algorithmes de codage précités 2:0 réalisent généralement une mise en forme spectrale du bruit (cf publication de Johnston précitée) et, même relati- vement au bloc contenant une transition, telle que men- tionnée précédemment, le spectre du bruit est à peu près plat et d'un niveau constant pour toute la durée du bloc. 255 En conséquence, pour la partie précédant l'attaque, le spectre du bruit est bien supérieur à celui du signal. Ceci se traduit dans le domaine temporel, ainsi que représenté en figure lc et ld relatives à un signal original x(n) d'une attaque vibraphone respectivement à un 3 signal correspondant décodé puis soumis à transformée inverse x'(n), par l'apparition d'une dégradation très importante appelée préécho. En outre, lorsque le système de codage utilise la transformée TCDM, les deux blocs précédant la transition 355 sont, du fait du recouvrement de N/2 échantillons entre fenêtres successives, affectés par le phénomène de prééchos. Différentes méthodes de traitement ont été propo- sées pour réduire ou éliminer le phénomène de prééchos. Parmi celles-ci on peut citer la méthode préconisée par Schroeder décrite dans l'article intitulé "Stereo Coding with CD quality" International Conference on Consumer Electronics, Chicago 1987 et par Sugiyama décrite dans l'article "Adaptive Transform Coding with an adaptive Block Size", Proceedings of ICASSP 90, Albuquerque p 1093 -10 1096 Dans la méthode précitée, un facteur d'échelle est appliqué au bloc d'échantillons avant la transformation, de façon à réduire la différence de niveau avant et après la transition engendrée par l'attaque L'application des facteurs d'échelle inverses au décodage permet de réduire15 le niveau de bruit dans les zones de faible énergie et donc l'amplitude, d'o la diminution correspondante du phénomène de préécho Cependant une telle méthode ne peut être valablement utilisée lorsque la transformée fréquentielle est une transformée TCDM, une reconstitution parfaite du20 signal, en l'absence de codage, ne pouvant ête effectuée, du fait de la présence même des facteurs d'échelle. La présente invention a pour objet la mise en oeuvre d'un procédé et d'un système de traitement des prééchos d'un signal audio-numérique codé par transformée25 TCDM dans lesquels les prééchos sont sensiblement suppri- més. Un autre objet de la présente invention est la mise en oeuvre d'un procédé et d'un système de traitement des prééchos d'un signal audio- numérique codé par transformée30 TCDM dans lesquels aucune modification des niveaux de signal, avant ou après attaque, n'étant effectuée, une reconstitution quasi-parfaite du signal après décodage est ainsi obtenue. Le procédé de traitement des prééchos d'un signal audio-numérique codé par transformée fréquentielle, objet 4 de la présente invention, est remarquable en ce qu'il consiste, après codage de ce signal x(n) par tranformée fréquentielle en coefficients y(k), sur des fenêtres de pondération formées sur des blocs d'échantillons comportant5 N échantillons x(n) successifs, à effectuer au décodage des coefficients de transformée y(k), un filtrage optimal de façon à réduire, lors de la présence d'une transition ou variation d'amplitude du signal, le bruit additif résultant lors du codage de la quantification des coefficients de10 transformée. Le système de traitement des prééchos de filtrage d'un signal audio-numérique codé par transformée fréquentielle, objet de la présente invention, ce système comportant, d'une part, disposés en cascade au niveau du15 codeur un module de lecture des échantillons x(n) de signal audio-numérique constituant ces derniers en blocs de N échantillons, un module de calcul de la transformée TCDM délivrant, à partir desdits blocs d'échantillons, des coefficients de transformée y(k) correspondants et un20 module de codage par quantification de ces coefficients en vue de leur transmission par l'intermédiaire d'un canal de transmission, et, d'autre part, disposés en cascade au niveau du décodeur, en réception du canal de transmission, un module de décodage des coefficients y(k) et un module de25 calcul de la transformée inverse TCDM-1, ces modules constituant le décodeur et délivrant un signal décodé bruité x'(n), est remarquable en ce qu'il comporte en outre, d'une part, au niveau du codeur, un module de détection de tran- sitions repérées par le rang N-q de la transition dans le30 bloc courant considéré, un module d'extraction, à partir des blocs d'échantillons x(n), des paramètres d'énergie x(n) de l'échantillon x(n) ou d'un groupe d'échantillons successifs du signal d'origine, un module de codage, à partir des paramètres d'énergie x 2 x(n) de l'échantillon35 x(n) ou du groupe d'échantillons du bloc d'échantillons, permettant la transmission codée des paramètres d'énergie 2 x(n), et d'un signal indicateur de la présence d'une transition dans le bloc d'échantillons considéré, et, d'autre part, au niveau du décodeur, un module de reconstitution de paramètres d'énergie o 2 x(n) transmis, un5 module d'extraction des paramètres de filtrage al, a 2 à partir des échantillons x(n) du bloc précédent b-i et d'évaluation des énergies du bruit de codage e 2 b(n) et du bruit o 2 w(n) pour l'échantillon de rang N considéré, un module de filtrage de Kalman recevant en entrée de filtrage10 le signal décodé bruité x'(n) et en entrée de commande lesdits paramètres de filtrage al, a 2 et l'évaluation des énergies du bruit de codage o 2 w(n), ledit module de filtra- ge délivrant en sortie le signal décodé filtré xr(n). Le procédé et le système objets de la présente invention trouvent application à la transmission de signaux audionumériques, de qualité haute fidélité, à la transmis- sion, au codage, au décodage de signaux audio-numériques, à toutes les applications de la transmission de signaux audio-numériques, telles que le réseau numérique à20 intégration de services RNIS. Le procédé et le système objet de l'invention seront mieux compris à la lecture de la description et à l'observation des dessins ci-après dans lesquels, outre les figures la à lc relatives à l'art antérieur,25 la figure 2 représente un schéma illustratif de la mise en oeuvre du procédé objet de l'inven- tion, la figure 3 a représente un mode de réalisation du système selon l'invention au niveau du codeur des échantillons par transformée TCDM, la figure 3 b représente un mode de réalisation du système selon l'invention au niveau du décodeur des échantillons codés transmis par transformée inverse TCDM-1, la figure 4 a représente un détail de réalisation des circuits d'extraction des paramètres de filtrage optimal, filtrage de KALMAN, au niveau du décodeur tel que représenté en figure 3 b, la figure 4 b représente un détail de réalisation des circuits d'estimation des paramètres d'éner- gie des échantillons audio-numériques x(n) au niveau du codeur, tel que représenté en figure 3 a, et la figure 5 représente un chronogramme d'échan- tillons représentatifs d'un signal d'attaque de guitare échantillonné à 32 k Hz,d'un signal correspondant codé par transformée TCDM sans traitement des prééchos avec un débit de trans- mission de 64 Kbits/s et d'un signal correspon- dant codé avec traitement des prééchos au moyen du procédé et du système objets de la présente invention. Une description plus détaillée de la mise en oeuvre du procédé de traitement des prééchos d'un signal audio- numérique, codé de manière non limitative par transformée20 TCDM, conformément à l'objet de la présente invention sera donnée en liaison avec la figure 2 D'autres types de transformée fréquentielle peuvent cependant être utilisés, sans sortir du cadre de l'objet de la présente invention. De manière classique, le procédé comprend des étapes de codage du signal audio-numérique x(n), ces échan- tillons étant, par un processus de lecture des échantillons conformés en bloc, notés x(n), ces blocs comportant par exemple N échantillons successifs, chaque bloc étant soumis à une fenêtre de pondération formée sur ces blocs, ainsi que représenté par exemple en figure la Les blocs précités sont ensuite soumis à un traitement par transformée TCDM, cette transformation délivrant les coefficients de trans- formée y(k) Les coefficients de transformée précités sont ensuite soumis à un processus de codage noté H, ce proces-35 sus de codage délivrant des coefficients de transformée codés Cy(k) destinés à être transmis vers et en vue d'un 7 processus de décodage. Ainsi qu'on l'observera sur la figure 2, le procédé objet de la présente invention consiste à effectuer, au décodage des coefficients de transformée y(k), ce processus5 de décodage étant noté H-1 au niveau du décodage proprement dit, et plus particulièrement suite à un processus de transformation inverse noté TCDM-1 appliqué au coefficient de transformée y(k) pour obtenir les échantillons x'(n) de rang N du signal décodé bruité, un filtrage optimal de ces10 échantillons du signal décodé bruité de façon à réduire lors de la présence d'une transition d'amplitude du signal audio-numérique le bruit additif résultant lors du codage de la quantification des coefficients de transformée. De manière générale, on notera que la présence des prééchos au niveau du signal audio-numérique codé est en fait due à la dégradation du signal précédant la transition par introduction d'un bruit additif b(n) résultant de la quantification des coefficients de transformée y(k). Ainsi, le traitement et la réduction de ce phéno- mène, conformément à l'objet du procédé selon la présente invention, est assimilable à un problème de débruitage lequel peut être conduit au moyen d'un filtrage optimal du signal audio-numérique décodé. Ainsi qu'on l'a représenté en figure 2, le procédé objet de la présente invention implique de manière générale au niveau du codage un processus d'identification de la transition Une telle identification peut être effectuée de manière non limitative par la détection d'une telle transition dans le bloc courant considéré, par la création30 par exemple d'une variable t représentative de l'existence ou de la non existence d'une telle transition dans le bloc précité En outre, et de manière avantageuse non limitative, l'identification de la transition dans le bloc courant considéré pourra être réalisée par la35 discrimination de la position de la transition dans le bloc précité, c'est-à-dire en fait de l'abscisse M de l'échan- 8 tillon représentatif de la première occurence de cette transition On notera enfin que, outre l'identification de la transition dans le bloc considéré, le procédé objet de la présente invention implique également l'évaluation d'un5 paramètre d'énergie de chacun des échantillons x(n), ce paramètre d'énergie pour chacun des échantillons x(n) du signal audio-numérique précité étant noté o 2 x(n). Bien entendu, l'ensemble des paramètres précités, paramètres d'identification de la transition, est alors transmis au décodage sous forme du paramètre t de présence ou d'absence de transition précité et du paramètre noté C,M lesquels représentent les paramètres codés de l'abscisse dans le bloc de la transition ainsi que d'un paramètre Cî 2 x(n), lequel représente le paramètre d'énergie codé pour15 chaque échantillon x(n) du signal audio-numérique, l'ensemble des paramètres précités désignés par la référence Cp étant alors transmis vers le processus de décodage conjointement avec les coefficients de transformée codés notés Cy(k).20 Une mise en oeuvre particulièrement avantageuse du procédé objet de l'invention peut être réalisée compte tenu du fait que le bruit additif b(n) est non stationnaire En effet, l'énergie du bruit précité est plus forte pour n proche de la transition dans le bloc considéré que pour les25 échantillons situés en début du bloc. Conformément à la remarque précitée, une mise en oeuvre particulièrement avantageuse du procédé objet de la présente invention consiste à réaliser le filtrage optimal précédemment mentionné sous la forme d'un filtrage de30 KALMAN Le filtrage de KALMAN effectué au décodage permet alors, conformément au procédé de la présente invention, d'obtenir la réduction du phénomène des prééchos de manière particulièrement satisfaisante, ainsi qu'il sera décrit ci-après dans la description.35 D'une manière générale, on rappellera que le fil- trage de KALMAN suppose une modélisation de l'expression des signaux audio-numériques sous forme d'une équation d'état et d'une équation d'observation. Le modèle d'état est associé au signal original x(n) tandis que l'équation d'observation décrit le signal décodé x'(n) bruité, c'est-à- dire affecté des prééchos. L'équation d'état peut être définie sous forme auto-régressive d'ordre p La valeur de p doit être choisie de façon à assurer un compromis entre une modélisation suffisamment exacte de la suite des échantillons du signal10 audio-numérique x(n) et un volume et une charge de calcul limité pour assurer un traitement compatible en temps réel. Pour une description plus détaillée de ce type de défini- tion autorégressive on pourra se reporter à la publication de Jayant et Noll intitulée "Digital Coding of Waweforms"15 1984 Prentice Hall Signal Processing series. Les équations matricielles du filtrage sont données par les relations: X(n) = A X(n-l) + W(n) x'(n)= C X(n) + b(n) avec X(n) = x(n) C=l 1,0,0 0 l A= al a 2 ap x(n-l) 1 O O O x(n-2) O 1 O O i x(n-p) x'(n) et b(n) étant des scalaires Les équations générales du filtre s'écrivent: * estimation: Xr(n)= A Xr(n-l)+ K(n)( x'(n) C A Xr(n-l)) * calcul du gain: K(n) = Pl(n) CT lC Pl(n) CT + b 2 (n)l-l * matrice de l'erreur de prédiction: Pl(n)=A P(n-1) AT+Q(n) avec Q(n) = 2 w(n) 0 O 0 0 0 0 O * matrice de l'erreur d'estimation: P(n)=Pi(n)-K(n) C Pi(n) A titre d'exemple non limitatif, on considèrera comme adapté pour la mise en oeuvre du procédé objet de la présente invention une valeur p = 2. Dans ces conditions, l'équation d'état définie sous forme régressive d'ordre 2 s'écrit: x(n) = al x(n-1) + a 2 x(n-2) + w(n) ( 1) Dans la relation précitée, al et a 2 désignent les coefficients de prédiction du filtre et w(n) désigne l'am- plitude du bruit d'innovation pour l'échantillon x(n) de rang N considéré.10 En outre, l'équation ou relation d'observation s'écrit sous forme scalaire: x'(n) = x(n) + b(n) ( 2) Dans la relation précitée, x(n) désigne bien entendu le signal audio-numérique original, c'est-à-dire l'échantillon d'ordre N courant considéré, b(n) désigne l'amplitude du bruit de codage ou préécho pour ce même échantillon et x'(n) désigne l'échantillon de rang N du signal décodé bruité, c'est-à-dire l'échantillon obtenu au décodage, ainsi que représenté en figure 2 après transfor-20 mation inverse TCDM-1. Les équations du filtrage de KALMAN sont connues de l'état de la technique et compte tenu de la modélisation précédemment décrite, l'échantillon de rang N du signal décodé filtré, c'est-à-dire obtenu après filtrage optimal25 par filtrage de KALMAN est donné par la relation: xr(n) = al ( 1-Kl) xr(n-1) + a 2 ( 1-Kl) xr(n-2) + Kl x'(n) ( 3) xr(n-1)= xr(n-1) + K 2 ( x'(n)-(al xr(n-l)+a 2 xr(n-2))) Dans la relation précitée, K 1, K 2 désignent la valeur de gain du filtre vérifiant la relation: K 1 = Pl( 1)/(Pl( 1) + c 2 b(n)) ( 4) K 2 = Pl( 3)/(Pl( 1) + a 2 b(n)) ( 5) relations dans lesquelles les coefficients Pl( 1) et P( 3) sont donnés par la matrice de l'erreur de prédiction et par la matrice de l'erreur d'estimation du filtrage de KALMAN.35 On rappellera pour mémoire que la matrice de l'erreur de prédiction est donnée par les relations: 11 Pl( 1) = a 12 P( 1) + ala 2 (P( 2) + P( 3)) + a 22 P( 4) + ow 2 (n)( 6) P 1 ( 2) = al P(l) + a 2 P( 3) ( 7) Pl( 3) = al P(l) + a 2 P( 2) ( 8) P 1 ( 4) = P( 1) ( 9) 5 En outre, la matrice de l'erreur d'estimation est donnée par les relations: P( 1) = (l-Kl) Pl(l) ( 10) P( 2) = ( 1-K 1) Pl( 2) ( 11) P( 3) = Pl( 3) K 2 Pi(l) ( 12)10 P( 4) = Pl( 4) K 2 Pl( 2) ( 13) Les équations précitées sont calculées pour chaque valeur de N pour N compris, dans le bloc courant considéré, entre zéro et M-i, o M-1 est le numéro de l'échantillon immédiatement antérieur à la transition correspondant à15 l'échantillon d'ordre M, échantillon par exemple pour lequel une variation d'amplitude est supérieure à un rapport 4 par exemple. On notera que dans le cas d'un codage par transfor- mée TCDM les équations précitées doivent être appliquées aux deux blocs en recouvrement dans lesquels la transition est présente, c'est-à-dire les deux blocs en recouvrement correspondants. On notera que dans les relations précitées, 2 b(n) désigne l'énergie du bruit de codage et a 2 w(n) l'énergie du bruit d'innovation w(n) pour l'échantillon x(n) de rang n considéré. On notera que les équations précitées sont établies dans l'hypothèse o b(n) est considéré comme un bruit blanc indépendant de x(n), ces conditions d'application étant en30 général vérifiées En effet, lors de l'attaque, c'est-à-dire de l'apparition de la transition, le spectre à court terme résultant de la transformée par TCDM est sensiblement plat, il en est de même pour celui du bruit. D'autre part, on remarquera que ce bruit dépend largement35 de la portion la plus énergétique du signal, c'est-à-dire celle située après la transition, et par conséquent le 12 bruit b(n) est généralement décorrellé du signal x(n) précédant l'attaque proprement dite. Ainsi, conformément au procédé objet de la présente invention telle que représenté en figure 2, dans le cas o le filtrage optimal est un filtrage de KALMAN, ce procédé consiste lors du processus de codage à déterminer, d'une part, l'énergie 2 x(n) de l'échantillon x(n) du signal original et, d'autre part, l'existence d'une transition au niveau d'au moins un échantillon x(N-q) de rang N-q d'un10 bloc courant déterminé, avec: N q = M ( 14) Le procédé selon l'invention consiste alors, lors du processus de codage, à transmettre au niveau décodage la valeur des paramètres d'énergie précités e 2 x(n) et d'exis- tence de la transition ainsi que décrit précédemment.15 Puis, le procédé objet de la présente invention consiste lors du processus de décodage à estimer les para- mètres de filtrage, c'est-à-dire les coefficients de prédiction al, a 2 du filtrage à partir des échantillons du bloc b-l précédant, conformément aux relations précédemment20 données dans la description De même, les énergies de bruit de codage 2 b(n) et du bruit c 2 w(n) pour l'échantillon de rang N considéré sont évaluées ce qui permet d'assurer, au niveau du décodage, le filtrage adaptatif de type filtrage de KALMAN du signal décodé bruité x'(n) pour reconstituer25 le signal décodé filtré dans lequel les prééchos sont sensiblement supprimés. De manière symbolique et conformément au procédé objet de la présente invention le signal S obtenu après décodage est donné par la relation logique30 S = t xr(n) + t x'(n) ( 15) On comprendra bien sûr que dans la relation préci- tée, t prend la valeur -0 ou 1 en fonction de l'absence respectivement de la présence d'une transition pour le bloc courant considéré. On comprendra que d'une manière générale, afin 13 d'assurer une mise en oeuvre efficace du procédé objet de l'invention, certains paramètres nécessaires au filtrage peuvent être estimés au niveau du décodage alors que les paramètres destinés à suivre la non stationnarité du signal5 audio-numérique original x(n) doivent être évalués au niveau du codage puis transmis afin d'assurer le décodage et le filtrage optimal précédemment mentionné. Une description plus détaillée d'un système de traitement des prééchos d'un signal audio-numérique codé par transformée fréquentielle conformément au procédé objet de la présente invention sera donné en liaison avec les figures 3 a représentant plus particulièrement un codeur et 3 b représentant au contraire un décodeur constitutif du système selon l'invention. Bien entendu, les figures 3 a et 3 b sont représen- tées sous forme de schémas blocs, les différents schémas blocs permettant de réaliser les étapes du procédé objet de la présente invention tel que réalisé en figure 2. En conséquence, ainsi qu'on l'a représenté en figure 3 a par exemple, le système comporte, d'une part, disposé en cascade au niveau du codeur 1, un module 10 de lecture des échantillons x(n) de signal audio-numérique constituant ces derniers en bloc de N échantillons, un module 11 de calcul de la transformée TCDM délivrant à25 partir des blocs d'échantillons les coefficients de transformée y(k) correspondants et un module 12 de codage par quantification de ces coefficients en vue de leur transmission par l'intermédiaire d'un canal de transmission numérique noté CT sur la figure 3 a Le canal de transmis-30 sion numérique CT ne sera pas décrit en détail car il correspond à un canal de transmission de type classique. En outre, ainsi qu'on l'a représenté en figure 3 b, le système selon l'invention comporte, d'autre part, dis- posé en cascade au niveau du décodeur 2 et en réception du35 canal de transmission CT précité un module 20 de décodage des coefficients y(k) et un module 21 de calcul de la 14 transformée inverse TCDM-1, ces modules constituant le décodeur et délivrant un signal décodé bruité x'(n). Conformément à l'objet de la présente invention, le système comporte en outre, au niveau du codeur 1, un module 13 de détection des transitions repérées par le rang N-q de la transition dans le bloc courant considéré et un module 14 d'extraction à partir des blocs d'échantillons x(n) des paramètres d'énergie o 2 x(n) de l'échantillon x(n) courant correspondant ou d'un groupe d'échantillons successifs du10 signal d'origine. Enfin, un module 15 de codage est prévu à partir des paramètres d'énergie o 2 x(n) de l'échantillon x(n) ou du groupe d'échantillons du bloc d'échantillons, ce module de codage permettant d'engendrer des codes correspondants15 Cîy 2 x(n) et C,M afin de permettre la transmission codée des paramètres précités et du signal indicateur t de la présence d'une transition dans le bloc courant par l'intermédiaire du canal de transmission CT. En ce qui concerne le module 13 de détection des transitions précédemment mentionné, on indiquera simplement que celui-ci engendre l'indicateur t lequel prend par exemple la valeur 1 dans le cas d'un bloc comportant une transition, c'est-à-dire susceptible d'engendrer des prééchos, cet indicateur étant simplement codé sur un bit25 et transmis au décodeur Le module 13 délivre également la position de la transition dans le bloc considéré, c'est-à-dire l'abscisse du point M précité L'adresse de la transition, M est alors délivrée au module 15 de codage précédemment mentionné et codé par exemple sur log 2 (N) bits30 o N est de la taille du bloc de transformée pour transmis- sion sous la forme du code C,M précédemment mentionné. On notera que la détection des attaques peut être réalisée par examen de la forme d'onde temporelle, par exemple comparaison de la valeur de l'amplitude de35 l'échantillon considéré à une valeur de seuil dans le rapport 4 précédemment mentionné, le module 13 de détection de transition recevant bien entendu les échantillons x(n) du bloc considéré. D'autres solutions peuvent être envisagées par exemple comparaison du spectre à court terme du bloc considéré par rapport à celui calculé sur le bloc précédent. Ainsi qu'on l'a en outre représenté en figure 3 b, le système selon l'invention comprend au niveau décodeur 2 un module 22 de reconstitution des paramètres d'énergie10 transmis, ce module recevant à partir du canal de transmis- sion CT les codes d'adresse de transition C,M et d'énergie Cî 2 x(n) et délivre les paramètres correspondants d'adresse M et d'énergie e 2 x(n) On notera que le module 22 permet de reconstituer les paramètres indispensables au filtrage et15 qui sont inconnus au niveau du décodeur. En outre, ainsi que représenté en figure 3 b, un module 23 d'extraction des paramètres de filtrage al, a 2 est prévu avant de réaliser l'évaluation de ces paramètres à partir des échantillons x(n) du blocd'échantillons précédent et d'évaluation des énergies du bruit de codage 2 b(n) et du bruit a 2 w(n) pour l'échantillon de rang n considéré On notera d'une manière générale que l'extrac- tion des paramètres de filtrage précités peut être effec- tuée soit à partir du signal décodé bruité x'(n) soit des25 coefficients de transformée quantifiés y'(k) délivrés par le module 20 de décodage. Enfin un module 24 de filtrage de KALMAN est prévu lequel reçoit en entrée de filtrage le signal décodé bruité x'(n) délivré par le module 21 de transformation inverse30 TCDM-1, ce même module de filtrage 24 recevant en entrée de commande les paramètres de filtrage al, a 2 et d'évaluation des énergies du bruit de codage e 2 b(n) et du bruit W 2 w(n) délivré par le module 23 d'extraction des paramètres préci- tés En outre, le module-24 de filtrage reçoit les paramè-35 tres d'énergie reconstitués délivrés par le module 22 et d'adresse de la transition M Le module 24 de filtrage 16 délivre en sortie le signal décodé filtré xr(n). De manière non limitative et selon un aspect avantageux du système objet de la présente invention tel que représenté en figure 3 b, le décodeur comporte en sortie un interrupteur commandé 25 lequel reçoit, d'une part, sur une première entrée de commutation le signal décodé filtré xr(n) et sur une deuxième entrée de commutation le signal x'(n) délivré par le décodeur et plus particulièrement par le module de transformation inverse 21 En entrée de com-10 mande de commutation, l'interrupteur commandé 25 reçoit l'indicateur de présence de transition t Le signal S délivré correspond à celui indiqué en liaison avec la figure 2. Une description plus détaillée du module 23 d'ex- traction des paramètres de filtrage al, a 2 et d'évaluation des énergies de bruit de codage 2 b(n) et du bruit 2 w(n) sera donnée en liaison avec la figure 4. Ainsi que représenté sur la figure précitée, le module 23 d'extraction des paramètres de filtrage précité comporte en parallèle sur un port d'entrée des échantillons x'(n) du signal décodé bruité une première voie A de calcul des paramètres de filtrage al, a 2 et du bruit o 2 w(n) pour l'échantillon de rang N considéré et une deuxième voie B de calcul du bruit de codage a 2 b(n) Les paramètres de filtra-25 ge al, a 2 du bruit a 2 w(n) et du bruit de codage a 2 b(n) sont délivrés au module de filtrage de KALMAN 24. La structure du module 23 d'extraction des paramètres de filtrage telle que représentée en figure 4 a peut être justifiée de la façon ci-après. L'équation d'état donnée par la relation ( 1) nécessite la connaissance des coefficients de prédiction ai, c'est-à-dire al, a 2, lorsque cette équation d'état est définie sous forme auto-régressive d'ordre 2. Ces coefficients doivent par contre être calculés sur le signal non bruité En supposant qu'avant l'appari- tion de la transition d'attaque le signal est stationnaire, 17 alors les coefficients ai, a 2 peuvent être estimés au niveau du décodeur à partir des échantillons reconstitués du bloc précédent puique ces échantillons ne sont pas affectés par le phénomène de préécho.5 Ainsi, la première voie A de calcul des paramètres de filtrage comporte une mémoire d'échantillon 230 recevant les échantillons bruités x'(n) délivrés par le module 21 de transformation inverse A chaque bloc à l'issue du décoda- ge, les R derniers échantillons sont conservés en mémoire10 au niveau de la mémoire 230 Pour une taille de bloc N = 1024 par exemple, on peut prendre R = 256 Ces R échantillons servent d'entrée au dispositif de calcul des coefficients de prédiction La première voie A comporte en conséquence un circuit 231 de calcul des coefficients de15 prédiction ai ou plus précisément al, a 2 et de gain de prédiction Gp Les R échantillons mis en mémoire au niveau de la mémoire 230 servent d'entrée au dispositif de calcul des coefficients de prédiction si le bloc suivant est détecté comme étant un bloc à préécho, c'est-à-dire pour20 t=_. Le circuit 231 de calcul peut par exemple être constitué par un circuit de calcul numérique permettant le calcul des coefficients précités au moyen d'un algorithme de Levinson L'algorithme de Levinson précité ne sera pas25 décrit car cet algorithme est connu de l'état de la tech- nique Les coefficients de prédiction al, a 2 ou plus géné- ralement ai constituent l'une des entrées du module de filtrage 24 Ces coefficients sont supposés constants sur la durée du préécho.30 En outre, la première voie A comprend en cascade avec le circuit 231 de calcul des coefficients de prédic- tion et de gain de prédiction Gp un circuit 232 de calcul de la variance du bruit d'innovation W 2 w(n) Ce circuit reçoit d'une part la valeur reconstituée du paramètre35 d'énergie de l'échantillon 2 x'(n) transmis délivré par le module 22 correspondant et le gain de prédiction Gp La 18 variance du bruit a 2 w(n) vérifie la relation donnée par le modèle auto-régressif: Ow 2 (n) = ex 2 (n)/Gp ( 16) Dans cette relation c 2 x(n) est l'énergie du signal original audio-numérique non bruité, Gp représentant le gain de prédiction associé au modèle auto-régressif. On précise que le calcul du gain Gp peut être effectué au niveau du circuit 231 à partir des coefficients de prédiction al et a 2 ou plus généralement ai Une méthode10 de calcul utilisable est décrite par Jayant et Noll dans l'ouvrage Digital Coding of Waveforms 1984, Prentice Hall Signal Processing Series page 263. Ainsi qu'on l'a en outre représenté en figure 4 a, la deuxième voie B comporte successivement connecté en cascade un circuit 233 de calcul du paramètre d'énergie 2 x'(n) de l'échantillon x'(n) bruité Le circuit 233 de calcul est suivi d'un circuit 235 de calcul du paramètre de bruit de codage o 2 b(n) des échantillons. L'énergie du bruit de codage W 2 b(n) peut être déduite des énergies des signaux non bruités d'origine c 2 x(n) et du signal bruité e 2 x'(n) du signal bruité x'(n). En effet, dans l'hypothèse d'une décorrélation totale entre le signal original et le bruit de codage b(n) on peut écrire:25 2 x'(n) = 2 x(n) + c 2 b(n) ( 17) C 2 b (n) =2 x' (n),2 x(n) ( 18) On notera que les énergies des signaux audio- numériques d'origine et bruités peuvent être évaluées de diverses manières.30 Par exemple, l'énergie du signal bruit 2 x'(n) peut être calculée au niveau du décodeur à l'aide d'une équation récursive de la forme: e 2 x'(n) = a o 2 x'(n-1) + ( 1-a) x'2 (n) ( 19) Dans cette équation a est une constante comprise entre O et 1 et peut par exemple être prise égale à 0,95. Dans ce but, au circuit de calcul 233 est associé un circuit de retard 234 lequel restitue la valeur o 2 x' (n- 1) relative à l'échantillon d'ordre n-l du signal bruité x'(n-l) Bien entendu, le circuit 233 de calcul du paramè- tre d'énergie a 2 x'(n) permet alors de déterminer la valeur actuelle de l'énergie de l'échantillon bruité 2 x'(n) selon la relation 19. La valeur de l'énergie de l'échantillon non bruité o 2 x(n) est, on le rappelle, calculée au niveau du codeur, transmise au décodeur et délivrée par le circuit corres-10 pondant 22 On peut considérer c 2 x(n) constant sur toute la durée du préécho ce qui permet de ne transmettre ainsi qu'une seule valeur correspondante par bloc d'échantillons. Cependant, ainsi qu'on le remarquera à l'observa- tion de la figure lc I par exemple, le signal précédant l'attaque, c'est-à-dire l'apparition de la transition, peut contenir des variations énergétiques relativement impor- tantes dont le dispositif de filtrage précédemment décrit doit être informé. Dans ce but, une description plus détaillée du module 14 d'extraction des paramètres d'énergie au niveau du codeur sera donnée en liaison avec la figure 4 b. Selon la figure précitée, le module 14 comprend un module 140 de découpe de chaque bloc b de N échantillons en sous blocs d'un nombre déterminé d'échantillons successifs.25 Ainsi, le module 140 effectue cette découpe en sous bloc de dimension L La valeur du paramètre L peut être choisie de façon à ce que les variations locales d'énergie puissent être isolées dans un sous bloc A titre d'exemple non limitatif pour une fréquence d'échantillonnage égale à 32 K Hz la valeur de L peut être prise égale à L = 128. Le traitement des prééchos par filtrage est alors effectué sur R = ElM/Ll sous blocs de L échantillons, o la notation El l signifie partie entière On notera que les M-r L derniers échantillons précédant la transition peuvent35 ne pas être filtrés En effet, du fait des propriétés de masquage temporel et pour une valeur de L suffisamment faible, le bruit qui les affecte n'est pas audible Une description des effets du masquage temporel plus détaillée pourra être trouvée dans l'ouvrage publié par Zwicker intitulé "Psycho-acoustique, l'oreille récepteur d'infor- 5 mation" publié à Paris aux Editions MASSON en 1981. Dans chacun des M/L sous-blocs, le paramètre d'énergie e 2 x(n) peut être calculé au moyen de la rela- tion: (j+ 1)L-1 02 x(j) = (l/L) > x 2 ( 1) ( 20) l=j L Dans cette relation j représente le numéro du sous bloc parmi l'ensemble des M/L sous blocs ainsi constitué, O est codé et transmis au décodeur par l'intermédiaire du module de codage 15. Selon le mode de réalisation de la figure 4 b, le codage des paramètres d'énergie précités peut par exemple être réalisé conformément au mode de réalisation de cette figure par quantification des rapports successifs. Dans ce but, le module de codage comporte avanta- geusement successivement un diviseur 150 recevant les paramètres d'énergie 02 x(j) précédemment décrit, un quan- tificateur uniforme 151 et un codeur de Huffmann 152 La30 sortie du quantificateur uniforme 151 est en outre rebou- clée sur une entrée du diviseur 150 par l'intermédiaire d'un circuit de quantification inverse 153 suivi d'un circuit de retard d'un échantillon 154. On notera que le diviseur 150 délivre ainsi au quantificateur uniforme 151 les valeurs des rapports successifs e(j) donnés par la relation: e(j) = o 2 x(j)/c 2 x(j-1) ( 21) Ces rapports peuvent être quantifiés par le quantificateur uniforme précité constitué par un quantificateur à faible nombre de niveaux Les mots de code en sortie du 5 quantificateur uniforme 151 peuvent alors subir un codage par entropie par exemple, un codage de Huffmann au moyen du codeur 152. D'autres techniques telles que la quantification vectorielle peuvent cependant être utilisées Le module de codage 15 délivre alors le code Cî 2 x(n) correspondant pour transmission par le canal de transmission CT précédemment décrit. Au niveau du décodeur, en particulier au niveau du module 23 d'extraction des paramètres de filtrage al, a 2 tel que représenté en figure 4 a, l'énergie du bruit de codage e 2 b(n) est donnée par la relation: 2 b(n) = 2 x' (n) 2 x' (j) ( 22) Dans cette relation les indices j,n sont liés par l'équation: j L Dans la relation ( 22) précitée, on notera en particulier que e 2 x'(j) est la version quantifiée de W 2 x(J) au sous bloc j reçu au décodeur et e 2 x'(n) est estimé par l'équation ( 19) précédemment décrite. Afin de disposer d'une estimation de bruit lissé délivré par le module 23, notamment à la jonction des sous blocs en N = j L une équation récursive du même type que l'équation 19 peut être utilisée, la valeur de bruit lissé vérifiant la relation: a 2 b(n) = p î 2 b(n-1) + ( 1-g)(a 2 x'(n) o 2 x'(j)) ( 23) Dans cette relation O est une constante comprise entre zéro et un ayant pour valeur par exemple la valeur p = 0,85. La première valeur 2 x( 0) permettant l'initialisa- tion de la forme récursive donnée par la relation 24 doit être codée par rapport à une variable commune au décodeur 1 et au décodeur 2. A titre d'exemple non limitatif, il peut s'agir de l'énergie du bloc précédant estimé par la somme des éner- gies des coefficients de transformée quantifiés y'(k) Dans une telle hypothèse le rapport e(O) vérifie la relation:5 e(O) = e 2 x(O)/p 2 x ( 24) Dans la relation précitée p 2 = 2 y,2 (k) ( 25) Dans la relation précédente, les termes y'2 (k) sont les coefficients quantifiés du bloc précédent En effet, d'après la relation de Parseval, P 2, = W 2 au bruit de codage près Toutefois dans le cas d'un codeur utilisant la transformée TCDM la relation de Parseval n'est pas vérifiée mais du fait des particularités de cette transformée le terme p 2, peut être approximé par:15 P 2 x = 4y 2 (k) ( 26) Le procédé et le système de la présente invention permettent d'obtenir des résultats particulièrement remar- quables car ils permettent d'obtenir une forte réduction des prééchos. 20 Sur la figure 5 on a représenté respectivement en I le signal d'origine d'une attaque de guitare échantillonné à 32 K Hz sur une durée de 32 msec L'axe des amplitudes est gradué en valeur relative d'amplitude ce qui permet de constater, au niveau des échantillons immédiatement anté-25 rieurs à la transition proprement dite, une forte variation d'amplitude par rapport aux échantillons du même bloc en début de celui-ci. Au point II de cette même figure 5 on a représenté au contraire le signal codé sans traitement des prééchos, le codage étant effectué au moyen d'une transformée TCDM et la transmission étant effectuée avec un débit de 64 Kb/s. Ainsi qu'on le notera à l'observation du point II de la figure 5 précité, le phénomène des prééchos est très apparent, dans la mesure o le bruit de préécho, sur les35 mêmes échantillons immédiatement antérieurs à la transition elle-même, présente une amplitude crête à crête sensible- 23 ment 4 à 5 fois supérieure à l'amplitude du niveau de bruit des échantillons présents en début de bloc correspondant. Au point III de cette même figure 5 on a représenté le signal codé avec traitement des prééchos et l'on cons- tate que le bruit de préécho a pratiquement disparu pour donner une signature de signal pratiquement semblable à celle du point I de cette même figure 5 Bien entendu sur les points I, II et III les axes sont gradués selon les mêmes unités en valeur relative de façon à donner une bonne10 appréciation de la réduction du niveau des prééchos ainsi obtenu. Du fait du caractère adaptatif du procédé et du système objets de la présente invention la réduction du niveau des prééchos est particulièrement performante pour15 tout type de signal L'amélioration de la qualité qui en résulte est particulièrement sensible surtout lorsque le codeur utilise une transformée par TCDM sans que pour cela il y ait un accroissement sensible du retard du système de codage On comprendra en effet que le traitement, tant au20 niveau du codage qu'au niveau du décodage, est effectué en temps réel la transmission des échantillons débruités étant ainsi effectuée en temps réel au temps de transmission et de traitement d'un bloc près. REVENDICATIONS 1 Procédé de traitement des prééchos d'un signal audio-numérique codé par transformée fréquentielle, caractérisé en ce que ledit procédé consiste, après codage 5 dudit signal par transformée fréquentielle en coefficient y(k), sur des fenêtres de pondération formées sur des blocs d'échantillons x(n) comportant N échantillons successifs, à effectuer, au décodage des coefficients de transformée y(k), un filtrage optimal, de façon à réduire, lors de la10 présence d'une transition d'amplitude du signal, le bruit additif résultant lors du codage de la quantification des coefficients de transformée. 2 Procédé selon la revendication 1, caractérisé en ce que le filtrage effectué au décodage est un filtrage de KALMAN. 3 Procédé selon la revendication 2, caractérisé en ce que le processus de filtrage de KALMAN étant défini à partir d'une équation d'état définie sous forme régressive d'ordre 2 de la forme:20 x(n) = al x(n-l)+a 2 x(n-2)+w(n) o al et a 2 désignent les coefficients de prédiction du filtre et w(n) le bruit d'innovation pour l'échantillon x(n) de rang N considéré et d'une relation d'observation de la forme:25 x'(n)=x(n)+b(n) o b(n) désigne l'amplitude du bruit de codage ou préécho pour ce même échantillon x(n), x'(n) désignant l'échantillon de rang N du signal décodé bruité et l'échantillon de rang N du signal décodé filtré étant donné par la relation: xr(n)=al (l-K 1) xr(n-1)+a 2 ( 1-Kl) xr(n-2)+Kl x'(n) o K 1 désigne la valeur de gain du filtre vérifiant les relations: K 1 = Pl(l)/(Pl( 1)+ 2 b(n)) K 2 = Pl( 3)/(Pl( 1)+e 2 b(n)) les coefficients Pl(l) et P 1 ( 3) étant donnés par la matrice de l'erreur de prédiction et par la matrice de l'erreur d'estimation du filtrage de KALMAN, c 2 b(n) désignant l'énergie du bruit de codage et c 2 w(n) désignant l'énergie du bruit d'innovation w(n) pour l'échantillon x(n) de rang n considéré, ledit procédé consiste:5 lors du processus de codage, à déterminer, d'une part, l'énergie c 2 x(n) de l'échantillon x(n) du signal original et, d'autre part, l'existence d'une transition au niveau d'au moins un échantillon x(N-q) de rang N-q d'un bloc courant déterminé, puis à trans- mettre au niveau décodage la valeur des paramè- tres d'énergie o 2 x(n) et d'existence de la transition puis, lors du processus de décodage, à estimer les paramètres de filtrage, coeffi- cients de prédiction al, a 2 du filtrage à partir des échantillons du bloc précédent, à évaluer les énergies du bruit de codage a 2 b(n) et du bruit e 2 w(n) pour l'échantillon de rang n considéré, ce qui permet d'assurer, au niveau du décodage, le filtrage adaptatif du signal décodé bruité x'(n) pour reconstituer le signal décodé filtré dans lequel les prééchos sont sensi- blement supprimés. 4 Système de traitement des prééchos d'un signal audio-numérique codé par transformée fréquentielle, ce système comportant, d'une part, disposés en cascade au niveau du codeur ( 1) un module ( 10) de lecture des échantillons x(n) de signal audio-numérique constituant ces derniers en blocs de N échantillons, un module ( 11) de calcul de la transformée fréquentielle délivrant à partir desdits blocs d'échantillons, des coefficients de transformée y(k) correspondants et un module ( 12) de codage par quantification de ces coefficients en vue de leur35 transmission par l'intermédiaire d'un canal de transmission numérique et, d'autre part, disposés en cascade au niveau 26 du décodeur ( 2), en réception du canal de transmission, un module ( 20) de décodage des coefficients y(k) et un module ( 21) de calcul de la transformée fréquentielle inverse, ces modules constituant le décodeur et délivrant un signal décodé bruité x' (n), caractérisé en ce que ledit système comporte en outre, d'une part au niveau du codeur ( 1), un module ( 13) de détection des transitions repérées par le rang N-q de la transition dans le bloc courant considéré, un module ( 14) d'extraction, à partir des d'échantillons x(n), des paramètres d'énergie v 2 x(n) de l'échantillon x(n) ou d'un groupe d'échantillons successifs du signal d'origine, un module ( 15) de codage, à partir des paramè- tres d'énergie u 2 x(n) de l'échantillon x(n) ou du groupe d'échantillons du bloc d'échantillons permettant la transmission codée des paramètres d'énergie a 2 x(n) et d'un signal indicateur de la présence d'une transition dans le bloc d'échan- tillons considéré et, d'autre part, au niveau du décodeur ( 2), un module ( 22) de reconstitution de paramètres d'énergie e 2 x(n) transmis, un module ( 23) d'extraction des paramètres de filtrage al, a 2, à partir des échantillons x(n) du bloc précédent et d'évaluation des énergies du bruit d'innovation c 2 w(n) et du bruit de codage a-2 b(n) pour l'échantillon de rang n considéré, un module ( 24) de filtrage de KALMAN recevant en entrée de filtrage le signal décodé bruité x'(n) et en entrée de commande lesdits paramètres de filtrage al, a 2 et d'évaluation des énergies du bruit de codage e 2 b(n) et du bruit d'innovation a 2 w(n), ledit module ( 24) de filtrage délivrant en sortie le signal décodé filtré xr(n). Système selon la revendication 4, caractérisé en ce qu'il comprend en outre: un commutateur deux voies ( 25) dont une première voie commutée reçoit ledit signal décodé filtré xr(n) et dont une deuxième voie commutée reçoit ledit signal décodé bruité x'(n) délivré par le module ( 21) de transformation inverse, la voie de commande de commutation étant pilotée par ledit signal indicateur (t) de la présence d'une transition dans le bloc d'échantillons considéré.10 6 Système selon la revendication 4 ou 5, caractérisé en ce que ledit module ( 14) d'extraction des paramètres d'énergie 2 x(n) des échantillons x(n) du signal comporte disposés en cascade: un module ( 140) de découpe de chaque bloc de N échantillons en sous-blocs d'un nombre déterminé d'échantillons successifs, un module ( 141) de calcul de l'énergie o'2 x(j) d'un échantillon ou d'un groupe d'échantillons du signal d'origine. 7 Système selon l'une des revendications 4 à 6, caractérisé en ce que le module de codage comporte: un quantificateur uniforme ( 151) délivrant, à partir des valeurs de paramètre, des mots de code, un circuit ( 152) de codage de Huffmann. 8 Système selon l'une des revendications 4 à 7, caractérisé en ce que ledit module ( 23) d'extraction des paramètres de filtrage al, a 2 et d'évaluation des énergies du bruit de codage e 2 b(n) et du bruit e 2 w(n) comporte, en30 parallèle sur un port d'entrée des échantillons x'(n) du signal décodé bruité, une première voie (A) de calcul des paramètres de filtrage al,a 2 et de la variance du bruit d'innovation e 2 w(n) pour l'échantillon de rang n considéré, une deuxième voie (B) de calcul de l'énergie du bruit de codage e 2 b(n), les paramètres de filtrage al,a 2, du bruit d'innovation c 2 w(n) et du bruit de codage o 2 b(n) étant délivrés au module de filtrage de KALMAN. 9 Système selon la revendication 8, caractérisé en ce que ladite première voie (A) comporte successivement connectés en cascade: une mémoire d'échantillons ( 230), recevant les échantillons bruités x'(n) délivrés par le module ( 21) de transformation inverse, un circuit ( 231) de calcul des coefficients de prédiction (ai, al, a 2) et de gain de prédiction Gp, un circuit ( 232) de calcul de la variance du 2 bruit d'innovation d'échantillon 2 w(n), ce circuit recevant, d'une part, la valeur reconstituée du paramètre d'énergie de l'échantillon c 2 x'(n) transmis délivrée par le module ( 22) correspondant et, d'autre part, le gain de prédiction Gp, le bruit d'échantillon vérifiant la relation: a 2 w(n) = a 2 x,(n) Gp Système selon la revendication 8, caractérisé en ce que ladite deuxième voie (B) comporte successivement connectés en cascade: un circuit ( 233) de calcul du paramètre d'énergie e 2 x'(n) de l'échantillon x'(n) bruité, ce paramètre d'énergie vérifiant une relation récursive de la forme: 2 x'(n) = a a 2 x'(n-1)+( 1-a) x'2 (n), o a est une constante de valeur comprise entre zéro et un, un circuit ( 235) de calcul du paramètre de bruit de codage e 2 b(n) des échantillons, ce paramètre de bruit vérifiant une relation de la forme: 02 b(n) = e 2 x'(n) e 2 x(n).