FR 2969804 A1 20120629 FR 1061187 A 20101223 FILTRAGE PERFECTIONNE DANS LE DOMAINE TRANSFORME La présente invention concerne le filtrage des données numériques, notamment le filtrage des données audionumériques. L'usage du codage audio à réduction de débit nécessite des filtrages en post-traitement du décodage audio qui peuvent se manifester sous de diverses formes : - en personnalisant une écoute en changeant la couleur sonore de l'enregistrement, en particulier en atténuant ou amplifiant certaines bandes du spectre sonore, - en faisant intervenir le filtrage aussi en codage audio où des outils paramétriques permettent d'enrichir le signal audio reçu. Par exemple des techniques telles que le «Spectral Band Replication» (SBR) ou le «Parametric Stereo » (PS) (norme ISO/IEC 14496-3, norme MPEG-4 Audio) permettent respectivement de reconstruire les sons aigus à partir des basses fréquences et un son stéréo à partir d'un signal mono. Le traitement «MPEG Surround MPS» (norme ISO/IEC 23003-1, norme MPEG-D) étend l'approche du traitement PS à la reconstruction de plus de deux canaux sonores. Dans ces techniques, la reconstruction des parties manquantes du son est assurée par des opérations de recopie du signal et de filtrage. Par exemple le traitement SBR module vers les hautes fréquences des zones de basses fréquences et ajuste l'énergie du signal en fréquence. Cet ajustement permet, après décodage, d'obtenir un signal semblable au signal original (signal avant codage). Le traitement PS recrée à partir d'un signal mono deux signaux composites qui sont ajustés en énergie en fréquence afin, ici encore, de rendre le signal décodé ressemblant au signal original de référence. Le traitement MPS étend se principe à la génération de N signaux à partir de M voies sonores transmises (avec N?M). Les codeurs audio à réduction de débit utilisant des transformées, normalisés à MPEG, de type MP3, AAC ou USAC ou normalisés à l'ITU-T, comme le G.722.1, G.719, G.718 privilégient l'usage de transformées à échantillonnage critique. L'échantillonnage critique est, en codage à réduction de débit, une propriété importante. En effet, pour conserver une bonne efficacité de transmission, il convient de ne pas transmettre plus d'échantillons transformés qu'il n'y en avait dans le domaine temporel. Pour cette raison, dans les codeurs à réduction de débit actuels, seules des transformées à échantillonnage critiques sont employées. Il s'agit par exemple des transformées MDCT (pour « Modified Diserete Cosine Transform ») qui sont typiquement des transformées à coefficients réels. Ces transformées sont impropres à un filtrage sans artefact (car elles entraînent une distorsion de repliement dite « aliasing »). Pour mener des filtrages adéquats, deux familles de techniques peuvent être mises en oeuvre : 1. celles consistant à effectuer la transformation inverse, puis à appliquer un filtrage du type convolution ; 2. celles consistant à utiliser une transformation adaptée au filtrage (par exemple une transformée de Fourier à court-terme ou une autre transformation à valeur complexe, par exemple des filtres complexes de type PQMF pour « Pseudo Quadrature Mirror Filters »), qui n'est pas à échantillonnage critique, pour pouvoir réaliser cette opération de filtrage sans artefact : le filtrage consiste alors en une simple multiplication par coefficient transformé (égalisation). Cependant, il faut réaliser la transformation inverse de la transformation de codage, puis transposer dans le domaine à valeur complexe les échantillons, qui, après égalisation, sont rétablis, par transformation complexe inverse, dans le domaine temporel. Trois transformations sont donc nécessaires. Pour la mise en oeuvre de la première approche (point 1 ci-dessus), il faut réaliser un filtrage du type convolution après transformation inverse de codage, ce qui toutefois est couteux en opérations de calcul et peu polyvalent (peu de flexibilité dans les évolutions du filtre réalisé). La seconde approche (point 2) est bien plus polyvalente et il est aisé de changer les coefficients multiplicatifs (la fonction d'égalisation). En revanche, le nombre de transformations à appliquer entraîne une complexité importante. La présente invention vient améliorer la situation. Elle propose à cet effet un procédé de traitement d'un signal sous forme de blocs successifs d'échantillons, comportant un filtrage dans un domaine transformé de sous-bandes. Au sens de l'invention, le procédé comporte : - un traitement d'égalisation appliqué à un bloc courant dans le domaine transformé, et - un traitement d'ajustement du filtrage appliqué dans le domaine transformé à au moins un bloc adjacent au bloc courant. On entend par «bloc » toute succession d'échantillons, telle qu'une trame, ou encore une sous-trame dans certains types de formats de signal. Ainsi, l'invention propose un filtrage amélioré, dans le domaine transformé. Cette approche est avantageusement peu complexe car le traitement reste dans le domaine de la transformée initiale. Un avantage procuré consiste en une limitation des composantes de repliement audibles, tout en assurant fidèlement la caractéristique de filtrage initialement désirée. Dans une réalisation, le procédé comporte alors une étape préalable d'optimisation de paramètres de l'égalisation et de l'ajustement du filtrage, par estimation d'un repliement qu'induit l'égalisation. Le repliement est préférentiellement estimé dans un domaine obtenu à partir d'une transformée inverse du domaine des sous-bandes (par exemple dans le domaine temporel). Cette estimation dans le domaine direct permet une limitation plus efficace de la distorsion audible qu'induit le repliement (ou « aliasing »), et donc une optimisation plus fine des paramètres de l'ajustement de filtrage. Avantageusement, l'égalisation et l'ajustement de filtrage dans le domaine transformé comportent : - un traitement d'égalisation appliqué à un bloc courant, - un traitement d'ajustement de filtrage appliqué à au moins un bloc précédant dans le temps le bloc courant, et - un traitement d'ajustement de filtrage appliqué à au moins un bloc suivant dans le temps le bloc courant. Ainsi, dans cette réalisation, il est proposé de s'appuyer à la fois sur le bloc précédent, mais aussi sur le bloc qui suit immédiatement le bloc courant Dans une réalisation, l'égalisation et l'ajustement de filtrage comportent l'application d'un système matriciel comprenant : - une première matrice appliquée à un vecteur signal représentant le bloc courant, - une deuxième matrice appliquée à un vecteur signal représentant le bloc précédent, et - une troisième matrice appliquée à un vecteur signal représentant le bloc suivant. Avantageusement, la troisième matrice est la transposée de la deuxième matrice. Ainsi, l'invention propose en particulier des structures symétriques (par exemple par filtrage dans un banc de filtres modulé en cosinus à échantillonnage critique quelconque) qui permettent d'obtenir des fonctions simples à réaliser. Dans une réalisation, préalablement à l'application des matrices, les blocs sont transformés dans le domaine des sous-bandes par au moins une transformée modulée, par exemple de type MDCT. Dans une réalisation sophistiquée, la transformée peut être de type transformée modulée, à valeurs complexes (par exemple de type MCLT, ou encore de type PQMF). Plus généralement, la présente invention permet de mettre en oeuvre des structures pour corriger des filtres passe-bas, passe-bande ou autres, dans le domaine à valeurs réelles ou complexes, à l'aide de fonctions simples comme décrit ci-après. Dans une réalisation à transformée modulée et où le filtrage comporte une composante de coupure au-delà d'une sous-bande donnée : - la première matrice appliquée au vecteur signal du bloc courant comporte comme seuls éléments non nuls une succession d'éléments A, identiques, dans la diagonale de la 20 matrice, suivis d'un élément A-B pour la sous-bande donnée et d'un élément B pour la sous-bande qui suit la sous-bande donnée, et - la deuxième matrice appliquée au vecteur signal du bloc précédent comporte comme seuls éléments non nuls au moins deux éléments de valeur absolue identique et de signes opposés, disposés sur la diagonale de la matrice, respectivement pour la sous-bande donnée 25 et pour la sous-bande qui suit la sous-bande donnée. Avantageusement, dans le cas où une troisième matrice intervient, les deuxième et troisième matrices comportent un nombre d'éléments non nuls qui est fonction d'un degré choisi d'optimisation des paramètres de l'ajustement de filtrage, minimisant le repliement estimé.15 Ainsi, l'invention propose des structures efficaces en calcul, à nombre restreint de coefficients à ajouter. Mieux encore, il est possible de choisir le nombre de coefficients de matrice à gérer en fonction d'une complexité souhaitée, ou encore en fonction d'un compromis entre complexité et limitation du repliement. Dans un mode de réalisation permettant de réduire fortement (voire d'annuler pratiquement) la distorsion liée au repliement, pour un filtrage passe-bas, la première matrice s'exprime sous une forme : /1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1-ao 0 0 0 0 To = 0 0 0 0 ao 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ~0 0 0 0 0 0 0 0 OJ le coefficient 1-ao étant appliqué pour la sous-bande donnée, en la deuxième matrice s'exprime sous une forme : /0 0 0 0 0 0 0 0 0 0 0 0 a5 0 0 0 0 0 0 a3 0 -a4 0 0 0 0 -a3 -al -az 0 a5 0 0 -as 0 az a~ a3 0 0 0 0 a4 0 -a3 0 0 0 0 0 0 -a5 0 0 0 0 0 0 ~0 0 0 0 0 0 0 0 le coefficient -ai de la diagonale étant appliqué pour la sous-bande donnée, T1= la troisième matrice étant transposée de la deuxième matrice. Dans ces expressions, les coefficients ao, ai, az, a3, a4 et a5 sont des nombres réels positifs, le réel ai, au moins, étant non nul. En effet, dans le cas d'une combinaison linéaire de matrices comme décrit ci-après, le terme ao peut être nul pour une matrice donnée, ce qui peut être compensé par une autre matrice combinée à cette matrice donnée. Si tous ces coefficients sont non nuls pour une matrice donnée, une minimisation optimale de la distorsion de repliement pour un filtre passe-bas donne : ao=0,2548 al =0,1249 az= 0,0812 a3=0,0409 a4= 0,0192 a5=0,0132 Dans une réalisation où le filtrage à corriger résulte d'une combinaison linéaire de filtrages, le système matriciel de correction au sens de l'invention comporte au moins : - une combinaison linéaire correspondante de premières matrices appliquée au vecteur signal du bloc courant, - une combinaison linéaire de deuxièmes matrices appliquée au vecteur signal du bloc précédent, et - une combinaison linéaire de troisièmes matrices, transposées respectives des deuxièmes matrices, appliquée au vecteur signal du bloc suivant. Ainsi, l'approche au sens de l'invention peut être généralisée à des fonctions de filtrage et d'égalisation quelconques, en utilisant des coefficients d'ajustement de filtrage adaptés à partir d'une analyse de la distorsion à corriger. La présente invention vise aussi un programme informatique comportant des instructions pour la mise en oeuvre du procédé ci-avant lorsque ce programme est exécuté par un processeur. Un exemple d'organigramme de l'algorithme général d'un tel programme est décrit plus loin en référence à la figure 13. La présente invention vise aussi un dispositif de traitement d'un signal sous forme de blocs successifs d'échantillons, comportant des moyens de filtrage dans un domaine transformé de sous-bandes. Ces moyens appliquent en outre : - un traitement d'égalisation à un bloc courant dans le domaine transformé, et - un traitement d'ajustement de filtrage, dans le domaine transformé, à au moins un bloc adjacent au bloc courant. Un exemple de réalisation d'un tel dispositif est décrit plus loin en référence à la figure 14. D'ailleurs, d'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, et des dessins annexés sur lesquels : - la figure lA illustre schématiquement un premier traitement réalisant un filtrage S(z), puis réalisant une transformation directe suivie d'une transformation inverse, - la figure 1B illustre schématiquement un second traitement procédant à une transformée directe, suivie du traitement en sous-bandes désiré Sso(z), et enfin réalisant la transformée inverse, pour distinguer, avec la figure 1A, deux approches des systèmes polyphasés, - la figure 2 illustre schématiquement la multiplication d'un scalaire dans chaque sous-bande du domaine transformé pour représenter un filtrage quelconque, - la figure 3 illustre l'allure d'un filtrage linéaire (filtre passe-bas) appliqué sous forme matricielle dans le domaine transformé, - la figure 4 détaille l'allure fréquentielle du filtre de la figure 3, - la figure 5 représente la distorsion (en ordonnées), diminuée par optimisation du paramètre d'égalisation ao (en abscisses), dans une réalisation sans ajustement de filtrage, - la figure 6 représente les caractéristiques fréquentielles du filtre résultant de l'optimisation d'égalisation illustrée sur la figure 5, - la figure 7 représente les caractéristiques fréquentielles du filtre résultant de l'optimisation de l'égalisation et de l'ajustement de filtrage, - la figure 8 représente la réduction de distorsion observée du fait du repliement (ordonnées) en fonction du nombre de coefficients intervenant dans l'égalisation et l'ajustement de filtrage (abscisses), - la figure 9 illustre la fonction de filtrage, égalisation et ajustement de filtrage, réalisée à l'aide d'un jeu de coefficients ao, ai, az, a3, a4, as, dans le cas d'un filtre passe bande, - la figure 10 illustre le cas d'une transformée modulée à valeurs complexes (de type MCLT), - la figure 11 compare la réduction de distorsion observée du fait du repliement (ordonnées) en fonction du nombre de coefficients intervenant dans l'égalisation et l'ajustement de filtrage (abscisses), pour une transformée à valeurs réelles (de type MDCT, en trait plein) et pour une transformée à valeurs complexes (de type MCLT, en traits pointillés), pour un filtrage passe-bas, - la figure 12 compare la réduction de distorsion observée du fait du repliement (ordonnées) en fonction du nombre de coefficients intervenant dans l'égalisation et l'ajustement de filtrage (abscisses), pour une transformée à valeurs réelles (de type MDCT, en trait plein) et pour une transformée à valeurs complexes (de type MCLT, en traits pointillés), pour un filtrage passe-bande, - la figure 13 résume les étapes d'un procédé au sens de l'invention, dans un exemple de réalisation, et - la figure 14 illustre schématiquement un dispositif pour la mise en oeuvre de l'invention, à titre d'exemple de réalisation. On rappelle ci-après quelques principes des systèmes polyphasés. On cherche à mesurer en particulier ici l'effet d'une fonction de filtrage en sous-bandes (domaine transformé) dans le domaine temporel (après transformation inverse). A cette fin, on cherche à mettre en relation deux traitements : - un premier traitement (a) réalisant un filtrage S(z), puis réalisant une transformation directe suivie d'une transformation inverse, - un second traitement (b) procédant à une transformée directe, suivie du traitement en sous-bandes désiré Ssb(z), et enfin réalisant la transformée inverse. On en extrait ensuite une expression du filtrage S(z). Les deux traitements sont respectivement illustrés sur les figures lA et 1B. Pour identifier simplement ces traitements, on écrit les opérations réalisées dans le domaine polyphasé. Il s'agit d'une approche classique pour simplifier la résolution de systèmes incluant des transformées avec échantillonnage. Le banc de filtres d'analyse (ou la transformée directe) est exprimé par sa matrice polyphase à l'ordre M, E(z). Le banc de filtres de synthèse (ou transformée inverse) est exprimé par sa matrice polyphase à l'ordre M, R(z). M représente le nombre de coefficients transformés (c'est-à-dire le nombre de coefficients fréquentiels obtenus par la transformée). La décomposition polyphasée des transformées modulées, incluant la transformée MDCT, s'exprime par : M-1 Hk (z) = z-1 Ek,l (zM ) l=0 On écrit également les composantes polyphases des transformées, comme suit, en s'appuyant sur les réponses impulsionnelles des filtres d'analyse ha,k,n pour la sous-bande k et le coefficient n. On se restreint dans cet exemple, sans perte de généralité toutefois, à une transformée dont les réponses impulsionnelles ont une longueur 2M, comme la MDCT. Eo,o (z) Eo,1 (z) Eo,M-1 (z) E(z) = E1,0 ( z ) El,l ( z ) 4,11-1 ( z ) E0 11_1 ( z ) E1,M-1 ( z ) EM-1,M-1 ( z ) ha ,0 ,0 + z-1ha,0,M ha,0,1 + z 1ha,0,M+1 h a,O,M-1 + z-1ha,0,2M-1 ha ,1 ,0 + z-1ha,M ha,l,l + z 1ha,1,M+1 h a,1,M-1 + z-1ha,1,2M-1 ha,M-1,0 + z-1ha,M-1,M hM-1,1 + z 1hM-1,M+1 ha ,M -1,M -1 + z- ha ,M -1,2M-1 ha ,0,0 ha ,0,1 7ha,o,M-1 ha ,0,M ha ,O,M+l ha ,0,2M-1 ha ,1,0 ha ,l ha ,l -1 ha,M-1,0 ha,M-1,1 ha,M-1,M-1 ha,M-1,M ha,M-1,M+1 ha,M-1,2M-1 On introduit alors, pour une transformée modulée, les modulations ckn et le filtre prototype d'analyse ha,n, avec : =cos/~/ n+1 ~k+- ,avec 0 Le filtre hs n est ici un filtre prototype (appelé «fenêtre de synthèse ») contenant 2M échantillons, 15 certains pouvant être nuls (notamment ceux d'indices les plus faibles). Pour le cas de la transformée MDCT, la reconstruction est parfaite dans la mesure où les modulations et les filtres d'analyse et synthèse assurent les conditions suivantes : P(z)=R(z)E(z) _ (hs oC+ Z hs,1C0 ) (Coha,o + Z Cljla l = hs,OC Coha,o + z-1 [hs,OC Clha,l + hs,1COCoha,o 1+ Z~~hs,1COClha,1 impliquant alors les deux conditions : CCo = C'o C1 = 0 et hs,oC Clha,l + hs,lCoCoha,o = IM (pour un choix judicieux de ha,n et hs,n). IM étant la matrice identité de taille MxM. Ainsi : P(z)=R(z)E(z)=z-lIM Le transformée MDCT est donc à reconstruction parfaite (au prix d'un retard d'une trame, c'est-à-dire de M échantillons, dans le cas d'un signal comportant une succession de trames de M échantillons chacune). Un corollaire de la reconstruction parfaite de la transformée MDCT implique que : 15 R-«z)R(z)E(z)= z-IR-«z) On dégage ainsi une propriété utile pour ce qui suit : R-«z)= zE (z) On décrit maintenant l'influence, en termes de repliement (ou « aliasing »), d'une opération 20 d'égalisation de filtrage dans le domaine transformé. Cette influence est analysée dans le domaine temporel après transformée inverse. 5 On part donc d'un traitement en sous-bande Ssb (z) et on cherche à estimer la fonction de filtrage résultante S ( z ) . Les traitements étant exprimés dans le domaine polyphasé, les notations Sse (z) et S (z) visent des matrices de filtres. Il est à noter que S(z) ne représente pas nécessairement un filtre linéaire réalisable par une convolution. En repartant de l'identification des traitements (a) et (b) : Y(z)=R(z)Sse(z)E(z)X(z)=S(z)R(z)E(z)X(z) R(z)Sse(z)=S(z)R(z) et en multipliant par R-1(z) des deux côtés : 10 R(z)Sse(z)R-1(z)=S(z)R(z)R-1(z) il vient (avec R-1(z) = zE (z) , comme montré précédemment) R(z)Sse(z)zE(z)=S(z) S(z)=zR(z)Sse (z)E(z) En remplaçant R et E par leur expression, qui dépend du banc de filtres utilisé, on peut alors 15 calculer la matrice S(z). Les éléments de cette matrice sont indexés comme suit : s0,0 (z) s0,1 (2) s0,2 (z) S0,M-2 (2) so,M-1 (2) s1,0 (2) (2) s1,2 (2) s1,M-2 (2) ~(~~ _ s2,0 (z) s2,1 (2) S2,2 (2) S2,M-2 (2) S2,M-1 (2) SM-2,0 (2) SM-2,1 (2) SM-2,2 (2) SM-2,M-2 (2) SM-2,M-1 (2) Sm-1,0(z) Sm-1,1(z) Sm-1,2(z) SM-1,M-2 (z SM-1,M-1 ( 2 Cette matrice se décompose en : S(z) = Slin (z)+Salias (z), ou : _ Slin ( z ) est une matrice de filtrage prenant la forme circulante suivante (garantissant le fait qu'elle corresponde à un filtre linéaire réalisable sous forme d'une convolution Slin,O (Z) slin,l (Z) slin,2 (Z) Z 1Slin,M-1 (Z) Slin,O (Z) Slin,l (z) Slin,M-1 (Z ) Slin,M-2 (Z) Slin (Z) = Z 1Slin,M-2 (Z) Z 1Slin,M-1 (Z) Slin,O (Z) Slin,M-3 (Z) Z-1Slin,1 (z) Z 1Slin,2 (Z) Z-1Slin,M-1 (Z) Slin,O (Z) les coefficients du filtre linéaire s'écrivant : M-1 (Z) _ - l Slin ~j Slin,nZ-n _ - Z Slin,l ZM , - et Salias (z) est une matrice de filtrage quelconque qui représente des composantes correspondant à des repliements dus à des inversions dans le signal selon l'axe temporel et/ou fréquentiel et créant des composantes additionnelles associées. Pour l'estimation de la partie linéaire Si. (z) , on utilise préférentiellement une estimation au sens 10 des moindre carrés, en minimisant la puissance du terme Salias (z) . On cherche donc à observer la contribution principale de filtrage linéaire présente dans la matrice Ssb (z) . On cherche ainsi à minimiser une expression du type : >alias (z) r - S(Z)- Slin (z) r On peut calculer les termes de la matrice Slin (z) en estimant la moyenne des termes diagonaux de 15 la matrice S (z) , comme suit : 1 Slin,O ( z ) _ M M-1 1 Sk,k (Z ) k=0 1 et Slin,m ( z ) M M-1-m m-1 ~j Sk,k+m (Z) + ZI SM+k-m,k ( Z ) k=0 k=0 pour les indices 0 Ainsi, en observant les composantes de repliement créées par la matrice de filtrage en sous-bandes Ssb, il apparaît que l'on peut réduire ces composantes de repliement en ajoutant des composantes de filtrage 1-ao et ao anti-repliement à cette matrice. Au sens de l'invention, on va plus loin dans cette approche en tentant une correction de l'effet du repliement dans le domaine temporel. Ainsi, on s'intéresse non seulement au signal à traiter dans une trame donnée pour l'égalisation précédemment décrite, mais on s'intéresse en outre au signal de trame précédente et/ou suivante de cette trame donnée. Toujours dans cet exemple de réalisation d'une transformée de type MDCT, il est proposé une première matrice TI à appliquer à la trame précédente (sous forme d'un vecteur signal) et, du fait de la symétrie de ce type de transformée, une deuxième matrice T'i qui se déduit de la première matrice TI, à appliquer à la trame suivante. Ici, la deuxième matrice T'i correspond en particulier à la transposée de la première matrice TI. Ainsi, l'expression d'un tel traitement global s'exprime comme suit : Ssb (z) = 1 + To + z , où TI' est la matrice transposée de la matrice TI, les notations z-i et z visant respectivement la trame précédente et la trame suivante.30 On recherche, ici pour un filtrage passe-bas avec égalisation dans le domaine transformé, la forme que peut prendre cette matrice TI et qui permettrait une minimisation de l'aliasing (mesurée en «puissance de repliement» comme indiqué précédemment). Dans un premier temps, il est recherché une expression simple de la matrice TI ne comportant d'éléments non nuls que sur sa diagonale principale, pour limiter la complexité du traitement. Après estimation du repliement, il apparaît que cette matrice TI présente partout des éléments nuls sauf aux positions des coefficients 1-ao et ao de la matrice To. Les coefficients de la matrice TI occupant ces positions sont respectivement de type -ai et ai. Dans ce cas, bien entendu, la matrice transposée T'l est identique. Cette mise en oeuvre permet déjà de réduire avantageusement la distorsion de quelques dB par rapport à la simple égalisation basée uniquement sur la matrice To. On obtient en effet un niveau de distorsion abaissé à -32,71 dB, soit une amélioration de 8 dB vis-à-vis de la situation de l'art antérieur, avec des coefficients optimisés : ao=0,2845 et ai=0,1311. On réalise donc une fonction de filtrage sur les spectres consécutifs x_I,k, xo,k et xi,k suivant les composantes fréquentielles k, comme suit : y0,k - x0,k , si k i+1 On comprendra alors qu'une expression matricielle de l'ajustement de filtrage n'est pas nécessairement mise en oeuvre pour cette simple forme de réalisation. Il est possible alors de réduire encore la distorsion en enrichissant la forme de la matrice TI (en augmentant toutefois la complexité des calculs) et d'y ajouter d'autres termes non nuls au voisinage des coefficients -ai et ai. Dans l'exemple décrit ci-après, 14 coefficients non nuls sont obtenus dans la matrice TI et sa forme générale minimisant la distorsion est alors du type : T1 = 'o 0 0 0 0 0 0 0 0 0 a5 0 s 0 0 0 a3 0 -a4 0 0 -a3 -al az 0 0 -a5 0 az al a3 0 0 a4 0 -a3 0 0 0 0 -as 0 0 0 18 0 0 0 0 0 0 a5 0 s 0 0 0 0 0 0 0 ~0 0 0 0 0 0 0 0 OJ Le système matriciel de traitement global s'écrit alors : r0 0 0 0 0 0 0 0~ 0 0 0 0 0 a5 0 0 0 0 0 0 a3 0 -a4 0 0 Ssb (z) = z 0 0 -a3 -a, -a, 0 a5 0 0 -a5 0 a2 a, a3 0 0 0 0 a4 0 -a3 0 0 0 0 0 0 -a5 0 0 0 0 0 0 ~0 0 0 0 0 0 0 0 0, (.1 0 0 0 0 0 0 0` 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1-a, 0 0 0 0 + 0 0 0 0 ao 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, r0 0 0 0 0 0 0 0~ 0 0 0 0 0 -a5 0 0 0 0 0 0 -a3 0 a4 0 0 + z+' 0 0 a3 -a, az 0 -a5 0 0 as 0 -a2 a, -a3 0 0 0 0 -a4 0 a3 0 0 0 0 0 0 as 0 0 0 0 0 0 0 0 0 0 0 0 0 0, Selon une écriture simplifiée, on a : Ti, Ssb (z) = 1 +To + z , TI' étant la matrice transposée de la matrice TI, ce qui correspond à : - appliquer la matrice To, corrigée pour l'égalisation, à une trame courante dans le domaine transformé, - appliquer la matrice TI sur un signal transformé, calculé à partir du signal utilisant M échantillons antérieurs à celui utilisé pour la matrice To, donc finalement à la trame d'échantillons précédant une trame courante qui, elle, est traitée par la matrice To, - et à appliquer la matrice TI' sur un signal transformé, calculé à partir du signal utilisant M échantillons postérieurs à celui utilisé pour la matrice To, donc finalement à la trame d'échantillons qui succède la trame courante. On comprendra alors que les matrices TI, To et TI' sont respectivement appliquées à trois trames successives d'une transformée MDCT. La matrice TI a donc pour vocation de réduire le niveau de repliement introduit par la matrice To qui réalise une fonction d'égalisation. Dans l'expression des matrices ci-dessus, tous les éléments sont donc nuls, hormis : - quelques éléments de la diagonale principale de la matrice To et de la matrice TI , et - des éléments de la première et troisième diagonale supérieure de la matrice TI et de la première et troisième diagonale inférieure de la matrice TI . On rappelle que la matrice T'i est la transposée de la matrice TI. La matrice de filtrage, dans cet exemple de réalisation de l'invention, est donc parfaitement décrite par les éléments diagonaux suivants : - diag(T0,0) - diag(T1,0) - diag(T1,-1) - diag(T1,-3) - diag(T1,+1) - diag(T1,+3) On remarquera que les diagonales supérieures et inférieures sont identiques au signe prés : diag(T1,-1) = -diag(T1,+1) diag(T1,-3) = -diag(T1,+3)25 On obtient donc une forme plus compacte en écrivant uniquement les diagonales, comme suit : indice diag(T0,0) diag(TI,O) Diag(Tl,l) diag(T1,3) 0 1 0 0 0 M/4-2 1 0 0 as M/4-1 1 0 a~ -a4 M/4 1 ao -al a~ as M/4+1 ao a l a3 0 M/4+2 0 0 0 0 M-1 0 0 0 0 En minimisant l'erreur d'aliasing en fonction des variables libres ao à as, on obtient une puissance d'aliasing correspondant à -45,31 dB, soit un gain 20,6 dB par rapport à l'état de l'art. Dans cet exemple numérique, un bon choix des coefficients est le suivant : ao=0,2548 al =0,1249 az= 0,0812 a3=0,0409 a4= 0,0192 as=0,0132 On relèvera ici que les valeurs de ces coefficients a; décroissent avec l'indice i, ce qui signifie que pour réduire par exemple le nombre de coefficient a; à considérer (par exemple pour réduire la complexité du traitement), il suffit de : - fixer le coefficient as à 0 - et, éventuellement pour réduire encore la complexité, de fixer le coefficient a4 à 0, - etc. Les coefficients restant non nuls, ao, ai, az, etc., étant optionnellement ajustés pour réduire la distorsion. On peut calculer ainsi un jeu de coefficients ao, ..., ai pour chaque valeur de i allant de 1 à n (avec n=5 au maximum par exemple), de manière à choisir un jeu de coefficients ao, ..., ai adapté à une situation donnée (par exemple en termes de complexité de traitement dans un terminal ou autre équipement d'un réseau de télécommunication). On observe que le filtre linéaire résultant est toujours très proche de la fonction désirée, comme le montre la figure 7 illustrant un filtre passe-bas de caractéristiques très similaires. On comprendra en particulier que le contenu de la diagonale de la matrice To appliquée à la trame courante fixe finalement le gabarit fréquentiel désiré. On peut "dégrader" la représentation en augmentant le nombre de coefficients nuls, par exemple en imposant as à zéro. Dans ce cas, la solution de filtrage en sous-bande optimale introduit un niveau d'aliasing de -42,90 dB (au lieu -45,31 dB). On a représenté sur la figure 8 une réalisation dans laquelle, l'état de l'art préconisant l'usage de 16 coefficients « classiques », il est proposé d'ajouter ici entre 1 et 29 coefficients non nuls (avec alors le choix d'un jeu de six coefficients de base ao, ai, ..., as) pour obtenir une réduction de la puissance du repliement spectral allant de 4,47 à 20,6 dB. L'exemple de 29 coefficients non nuls ajoutés dans les matrices TI, To et TI' correspond à l'ajout : - du coefficient ao dans la matrice To, - des 14 coefficients non nuls présentés ci-avant dans la matrice TI, et - des 14 coefficients non nuls correspondant dans la matrice transposée ri. On montre ci-après que le cas d'un filtre passe-bas peut être transposé à toute forme de filtre, par exemple un filtre passe-bande. On décrit ci-après un mode de réalisation de l'invention dans le cadre de l'implémentation d'un filtre passe-bande. On prend alors pour effectuer l'égalisation un filtre de la forme : Tk=O pour 0k Ssb (Z) = TlZ-1 + To + Ti~z Seules certaines diagonales des matrices To et TI ne sont pas nulles. Leurs valeurs prennent la forme : indice diag(T0,0) diag(T1,0) Diag(T1,1) diag(T1,3) 0 0 0 0 0 M/8-3 0 0 0 0 M/8-2 0 0 0 -a5 M/8-1 0 0 -a; aa M/8 ao -a5 M/8+1 1-ao - al -a3 0 M/8+2 1 0 0 0 M/4-2 1 0 0 a5 M/4-1 1 0 a3 -a4 M/4 1-ao -al a, a5 M/4+1 ao al a~ 0 M/4+2 0 0 0 0 M-1 0 0 0 0 On a donc répliqué les fonctions anti-repliement de la structure passe-bas dans ce nouveau mode de réalisation. En cherchant les meilleures valeurs des coefficients an, on diminue significativement la quantité de repliement obtenu. En conservant les valeurs obtenues dans le mode de réalisation traitant d'un filtre passe-bas, on obtient un niveau de repliement réduit à -42,59 dB soit une amélioration de 20,91 dB en comparaison de la fonction de filtrage présentée plus haut (qui présente une distorsion de -21,68 dB). On présente ci-après un mode de réalisation appliqué à des bancs de filtres à valeur complexe. Un exemple de réalisation concerne les transformées MCLT (pour « Modulated Complex Lapped Transform ») décrites par exemple dans : H. Malvar, "A Modulated Complex Lapped Transform And Its Applications to Audio Processing". Proc. International Conference on Acoustics, Speech and Signal Processing, 1999. En référence à la figure 10, une transformée MCLT consiste en deux composantes : - une composante correspondant à la transformée MDCT du signal, identique aux modes de réalisation précédents, et - une composante correspondant à une transformée MDST du signal (pour «Modified Diserete Sine Transform ») : on remplace les termes en cosinus de la transformée MDCT par des termes en sinus pour obtenir son expression. On écrit alors une partie en cosinus de l'analyse comme suit (identique au cas exposé plus haut) : E, (z) =Cohao +z-iCIhai et une seconde partie en sinus comme suit : Es (z) = Sohao +z-iSihai Les transformations inverses s'écrivent de façon analogue : Re (z) = hs oCi + z-ihs Co iet Rs (z) = hs oSi + z-ihs iSo Les matrices Si contiennent les termes en sinus : 'z/ n+1 ~k+~ , 0 et on applique le type de traitement présenté dans les exemples de réalisation ci-avant : - pour paramétrer les traitements en sous-bandes Ssb - et pour estimer la composante linéaire de convolution et la quantité de repliement introduit. On peut appliquer simplement ainsi les résultats présentés pour le cas de la transformée MDCT dans les exemples de réalisation précédents. On utilise la matrice Ssb paramétrée comme s n = sm précédemment et la caractéristique de distorsion en fonction du nombre de coefficients peut être déduite de la même manière. Il convient de noter que pour le cas de la transformée MCLT, la matrice de filtrage Ssb a le même nombre de coefficients que pour la transformée MDCT. En revanche, du fait d'une double application à la fois sur la partie en cosinus et sur la partie en sinus, le nombre de coefficients à appliquer au total est double, comme illustré sur la figure 11 (pour un filtre passe-bas) et la figure 12 (pour un filtre passe-bande), lesquelles comparent le nombre de coefficients permettant de réduire l'effet d'aliasing pour une transformée MDCT (en trait plein) et pour une transformée MCLT (en traits pointillés). On décrit maintenant une mise en oeuvre dans le cas d'un filtrage quelconque. Les fonctions de filtrage présentées dans les modes de réalisation ci-avant peuvent se cumuler. Par exemple, le cas d'un filtre passe-bande présenté ci-avant peut être cumulé à celui d'un filtre passe-bas (sans restriction). Ainsi par exemple, la forme du filtre résultant peut être du type : - filtre passe-bas pour k allant de 0 à M/8, - et filtre passe-bande pour k allant de M/8 à M/4. On cumule simplement les contributions des deux filtres. Les deux tables correspondant aux contributions passe-bande et passe-bas sont présentées côte à côte ci-après. Passe Bande Passe Bas indice diag(T0,0) diag(TI,O) diag(Tl,l) Diag(T1,3) Diag(T0,0) diag(TI,O) diag(Tl,l) diag(T1,3) 0 0 0 0 0 1 0 0 0 M/8-3 0 0 0 0 1 0 0 0 M/8-2 0 0 0 -a5 1 0 0 a5 M/8-1 0 0 -a3 aa 1 0 a~ -a~ M/8 ao a, a~ -a5 I-ao -a~ a~ a5 M/8+1 I -ao -al -a3 0 ao a, a~ 0 M/8+2 1 0 0 0 0 0 0 0 M/4-2 1 0 0 a5 0 0 0 0 M/4-1 1 0 a~ a 0 0 0 0 M/4 I-ao a, a, a5 0 0 0 0 M/4+1 ao a, a~ 0 0 0 0 0 M/4+2 0 0 0 0 0 0 0 0 M-1 0 0 0 0 0 0 0 020 Les différentes matrices s'additionnent alors par linéarité et le résultat de la combinaison donne : indice diag(T0,0) Diag(T1,0) Diag(T1,1) diag(T1,3) 0 1 0 0 0 M/8-3 1 0 0 0 M/8-2 1 0 0 0 M/8-1 1 0 0 0 M/8 1 0 0 0 M/8+1 1 0 0 0 M/8+2 1 0 0 0 M/4-2 1 0 0 a5 M/4-1 1 0 a3 -a4 M/4 1-ap -al -a? 35 M/4+1 ao al a~ 0 M/4+2 0 0 0 0 M-1 0 0 0 0 On obtient bien une structure globale d'un filtre passe-bas. On relèvera les valeurs à zéro à la jointure des deux filtres combinés, ce qui signifie qu'une compensation de repliement n'est plus nécessaire ici, comme attendu, et la bande passante est de M/4 coefficients. Il est possible, aussi, de pondérer les deux filtres précédents en vue de réaliser une fonction de filtrage plus complexe. Par exemple, il peut être mis en oeuvre une fonction de filtrage : - appliquant un facteur g=10 sur les fréquences en-dessous de M/8, - conservant les fréquences entre M/8 et M/4 à leur valeur initiales, - et enfin coupant toute fréquence au-delà de M/4. Dans cet exemple de réalisation, on considère les éléments des deux matrices ci-dessous : Passe-Bande (pondération xl) Passe-Bas (pondération x10) indice diag(T0,0) diag(TI,O) diag(Tl,l) Diag(T1,3) Diag(T0,0) diag(TI,O) diag(Tl,l) diag(T1,3) 0 0 0 0 0 1 0 0 0 M/8-3 0 0 0 0 1 0 0 0 M/8-2 0 0 0 -a5 1 0 0 a5 M/8-1 0 0 -a3 aa 1 0 a~ -a4 M/8 ao al a~ -a5 I-ao -al az a5 M/8+1 I -ao -a l -a3 0 ao a l a~ 0 M/8+2 1 0 0 0 0 0 0 0 M/4-2 1 0 0 a5 0 0 0 0 M/4-1 1 0 a3 -a4 0 0 0 0 M/4 I-ao -al a, a5 0 0 0 0 M/4+1 ao a, a3 0 0 0 0 0 M/4+2 0 0 0 0 0 0 0 0 M-1 0 0 0 0 0 0 0 0 L'application de la pondération précitée mène aux nouvelles matrices : Passe Bande (pondération xl) Passe Bas (pondération x10) indice diag(T0,0) diag(TI,O) diag(Tl,l) diag(T1,3) Diag(T0,0) diag(TI,O) diag(Tl,l) diag(T1,3) 0 0 0 0 0 G 0 0 0 M/8-3 0 0 0 0 G 0 0 0 M/8-2 0 0 0 as G 0 0 (g-1).a5 M/8-1 0 0 a3 a G 0 (g-1).a, -(g-1).a~ M/8 ao ai az as g-(I-ao)+ao -(g-1).ai -(g-1 ).a2 (g-1 ).a.5 M/8+1 1-a. a, a3 0 (1-ao)+g-ao (g-1 ).al (g-1).a, 0 M/8+2 1 0 0 0 1 0 0 0 M/4-2 1 0 0 as 1 0 0 a5 M/4-1 1 0 a3 aa 1 0 a3 aa M/4 1-a' a, a a 1-a' a, a a M/4+1 ao a, a3 0 a~ ai a3 0 M/4+2 0 0 0 0 0 0 0 0 M-1 0 0 0 0 0 0 0 0 Les deux matrices peuvent être sommées comme dans la réalisation présentée précédemment. On peut reproduire ainsi toute sorte de gabarits d'égalisation en cumulant des fonctions passe-bande élémentaires (et/ou des fonctions passe-bas), pondérées, en les sommant. Ainsi, en termes généraux, partant d'une fonction de filtrage désirée traduite par une série de gains à appliquer à chaque bande d'une transformée, on répercute ces gains sur la diagonale principale de la matrice To. On peut ajuster ces gains afin de limiter les discontinuités d'une bande à l'autre (typiquement les valeurs des coefficients ao). Pour chaque différence de gains entre deux bandes, les gains respectifs sont pondérés par des gabarits de réduction de repliement définissant une modification des matrices TI et T'i destinées à pondérer les trames précédentes et suivantes dans le domaine transformé. On applique ensuite le filtrage. Ainsi, pour chaque trame de signal dans le domaine transformé MDCT, on multiplie trois trames consécutives (représentées sous forme de vecteur) respectivement par les matrices de filtrage : - TI, pour la trame précédant la trame courante à filtrer, - To, pour la trame courante à filtrer, et - T'1, pour la trame suivant la trame courante, et on somme les contributions des résultats. On obtient un spectre dans le domaine transformé MDCT, et on applique une transformation MDCT inverse pour obtenir un signal temporel présentant les caractéristiques de filtrage désirées. Cette démarche s'étend à une transformée à valeur complexe (par exemple de type MCLT) et plus généralement à toute transformée modulée, à valeur complexe ou non. On a résumé sur la figure 13 les principales étapes d'un exemple de réalisation d'un procédé au sens de l'invention. Dans une première étape 10, pour un filtrage donné, par exemple pour un passe-bande (un passe-bas ou un passe-haut étant considérés comme des cas particuliers de passe-bande), on détermine tous les jeux de coefficients : - ao, ai, - ao, ai, az, - ao, ai, az, as, - ao, ai, az, a3, .., an, par exemple avec n=5, permettant de construire les matrices To, TI, T'i qui minimisent la distorsion due au repliement. Cette étape 10 peut être réalisée préalablement, «hors-ligne ». Ensuite, pour une étape courante « en ligne » 11, où par exemple un critère de complexité doit être respecté, on choisit un jeu de coefficient ao, ai, ..., agi, avec i inférieur ou égal à n, tel que l'indice i respecte un compromis de complexité/qualité du filtrage, fixé par exemple en fonction de la capacité de calcul d'un terminal ou autre, ou encore en se fixant un niveau de qualité dépendant de la qualité de codage du son. En effet, un signal sonore codé est inévitablement distordu : il est donc inutile de réduire le repliement à des valeurs significatives en-dessous du niveau de bruit généré par le codage. Ainsi, la valeur de l'indice i peut être déterminée à cette étape 11, en fonction des conditions précitées, et, à l'étape suivante 12, on extrait par exemple d'une mémoire le jeu de coefficients ao, ai, ..., agi, correspondant à la valeur choisie pour l'indice i. Pour un filtrage donné F (étape 24), on détermine tous les filtrages élémentaires (F1, F2 ,..., Fk) constituant ce filtrage donné, par exemple un passe-bas FI (étape 13), un passe-haut F2 (formant donc un passe-bande, par soustraction, avec le filtrage élémentaire F1), un filtrage à partie complexe Fk (étape 20), ou autres. En particulier, le filtrage global F peut résulter d'une combinaison linéaire de filtrages élémentaires : F = 21F1 +...+2.kFk On applique donc aux matrices To , Tlk , T'lk correspondant à chaque filtrage élémentaire Fk un terme de pondération 2k (étapes 14 à 20) et on somme une à une les matrices pondérées pour obtenir les matrices To, T1, T'1 correspondant finalement au filtrage global F (étapes 21 à 23) et qui sont appliquées respectivement aux vecteurs représentant une trame donnée TRj, une trame précédente TR,_1 et une trame suivante TR;+1. Comme indiqué précédemment, la matrice T'i peut également être déduite de la matrice TI par transpositions. La présente invention vise aussi un programme informatique comportant des instructions pour la mise en oeuvre du procédé au sens de l'invention lorsque ce programme est exécuté par un processeur. On comprendra que la figure 13 peut correspondre à un organigramme de l'algorithme général d'un tel programme. La présente invention vise aussi un dispositif pour la mise en oeuvre du procédé et comportant alors des moyens de filtrage dans un domaine transformé de sous-bandes. En particulier, ces moyens appliquent: - un traitement d'égalisation à un bloc courant dans le domaine transformé, et - un traitement d'ajustement de filtrage, dans le domaine transformé, à au moins un bloc adjacent au bloc courant. On a représenté sur la figure 14 un exemple de réalisation d'un tel dispositif, avec, dans l'exemple représenté : - une entrée pour recevoir un signal S à traiter dans le domaine des sous-bandes, - une deuxième entrée optionnelle pour recevoir des conditions de détermination de l'indice i précité, de manière à choisir un jeu de coefficients ao, ai, ..., ai adaptés à un compromis de complexité/qualité du filtrage, - des moyens de traitement du signal tels que par exemple un processeur PROC et une mémoire de travail MEM, - et une sortie délivrant le signal filtré S'. La présente invention ne se limite pas aux formes de réalisation décrites ci-avant à titre d'exemple ; elle s'étend à d'autres variantes. Par exemple, on a présenté des modes de réalisation à trois spectres consécutifs issus des traitements de trames successives par les matrices TI, To, T'l. Néanmoins, le nombre de trames à traiter peut être plus grand s'il est souhaité réaliser des filtres à réponse impulsionnelle finie plus longs. En outre, dans un procédé de filtrage utilisant un filtre prototype d'analyse différent d'un filtre prototype de synthèse (formes de filtrage appelées «fenêtres d'analyse et de synthèse »), le traitement d'ajustement de filtrage par égalisation peut s'appliquer à un nombre de trames avant la trame courante, différent du nombre de trames après la trame courante. Par exemple, il est possible de ne traiter qu'une seule trame adjacente à la trame courante (antérieure ou postérieure). On obtient dans ce cas un filtre linéaire asymétrique. On a décrit ci-avant des matrices de traitement (en particulier les matrices TI et T'1) adaptées pour des transformées MDCT (notamment pour ce qui concerne la position des éléments non nuls de matrice). Toutefois, ces formes de matrices sont susceptibles de variantes pour d'autres types de transformée. Par exemple, la matrice T'i peut ne pas se présenter sous la forme de la transposée de la matrice TI pour un type de transformée différent de la transformée MDCT avec des filtres de Malvar comme réalisé ici. REVENDICATIONS1. Procédé de traitement d'un signal sous forme de blocs successifs d'échantillons, le procédé comportant un filtrage dans un domaine transformé de sous-bandes, caractérisé en ce qu'il comporte : un traitement d'égalisation (To) appliqué à un bloc courant dans le domaine transformé, et un traitement d'ajustement du filtrage (T1, T'1) appliqué dans le domaine transformé à au moins un bloc adjacent au bloc courant. 2. Procédé selon la revendication 1, caractérisé en ce qu'il comporte une étape préalable 10 d'optimisation de paramètres de l'égalisation et de l'ajustement du filtrage, par estimation d'un. repliement qu'induit l'égalisation. 3. Procédé selon la revendication 2, caractérisé en ce que le repliement est estimé dans un domaine obtenu à partir d'une transformée inverse du domaine des sous-bandes. 15 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'égalisation et l'ajustement de filtrage dans le domaine transformé comportent : un traitement d'égalisation (To) appliqué à un bloc courant, un traitement d'ajustement de filtrage (Tl) appliqué à au moins un bloc précédent dans le 20 temps le bloc courant, et un traitement d'ajustement de filtrage (T'1) appliqué à au moins un bloc suivant dans le temps le bloc courant. 5. Procédé selon la revendication 4, dans lequel les blocs courant, précédent et suivant sont 25 représentés par des vecteurs signaux, caractérisé en ce que l'égalisation et l'ajustement de filtrage comportent l'application d'un système matriciel comprenant : une première matrice (To) appliquée au vecteur signal du bloc courant, une deuxième matrice (Tl) appliquée au vecteur signal du bloc précédent, et une troisième matrice (T'1) appliquée au vecteur signal du bloc suivant. 30 6. Procédé selon la revendication 5, caractérisé en ce que la troisième matrice est la transposée de la deuxième matrice. 7. Procédé selon la revendication 6, caractérisé en ce que, préalablement à l'application des 35 matrices, lesdits blocs sont transformés dans le domaine des sous-bandes par au moins une transformée modulée. 31 8. Procédé selon la revendication 7, caractérisé en ce que la transformée modulée est de type MDCT. 9. Procédé selon l'une des revendications 7 et 8, caractérisé en ce que la transformée est de type transformée modulée à valeur complexe (MCLT). 10. Procédé selon l'une des revendications 6 à 8, dans lequel le filtrage comporte une composante de coupure au-delà d'une sous-bande donnée, caractérisé en ce que : - la première matrice (To) appliquée au vecteur signal du bloc courant comporte comme seuls éléments non nuls une succession d'éléments A, identiques, dans la diagonale de la matrice, suivis d'un élément A-B pour la sous-bande donnée et d'un élément B pour la sous-bande qui suit la sous-bande donnée, et - la deuxième matrice (TI) appliquée au vecteur signal du bloc précédent comporte comme seuls éléments non nuls au moins deux éléments de valeur absolue identique et de signes opposés, disposés sur la diagonale de la matrice, respectivement pour la sous-bande donnée et pour la sous-bande qui suit la sous-bande donnée. 11. Procédé selon l'une des revendications 5 à 10, prises en combinaison avec la revendication 2, caractérisé en ce que les deuxième et troisième matrices comportent un nombre d'éléments non nuls qui est fonction d'un degré choisi d'optimisation des paramètres de l'ajustement de filtrage, minimisant le repliement estimé. 12. Procédé selon l'une des revendications 10 et 11, caractérisé en ce que, pour un filtrage passe- bas, la première matrice s'exprime sous une forme : /1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1-ao 0 0 0 0 To = 0 0 0 0 ao 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ~0 0 0 0 0 0 0 0 OJ le coefficient 1-ao étant appliqué pour la sous-bande donnée,en ce que la deuxième matrice s'exprime sous une forme : /0 0 0 0 0 0 0 0 0 0 0 0 a5 0 0 0 0 0 0 a3 0 -a4 0 0 0 0 -a3 -al -az 0 a5 0 T1= 0 -a5 0 az al a3 0 0 0 0 a4 0 -a3 0 0 0 0 0 0 -a5 0 0 0 0 0 0 v0 0 0 0 0 0 0 0 le coefficient -ai de la diagonale étant appliqué pour la sous-bande donnée, et en ce que la troisième matrice s'exprime sous une forme : ~0 0 0 0 0 0 0 0 0 0 0 0 -a5 0 0 0 0 0 0 -a3 0 a4 0 0 0 0 a3 -al az 0 -a5 0 T'i = 0 as 0 -az al -a3 0 0 0 0 -a4 0 a3 0 0 0 0 0 0 as 0 0 0 0 0 0 v0 0 0 0 0 0 0 0 OJ le coefficient -ai de la diagonale étant appliqué pour la sous-bande donnée, les coefficients ao, ai, az, a3, a4 et a5 étant des nombres réels positifs, le réel ai, au moins, étant non nul. 13. Procédé selon la revendication 12, caractérisé en ce que les coefficients ao, ai, azi a3, a4 et a5 sont tels que : ao=0,2548 al =0,1249 az= 0,0812 a3=0,0409 a4= 0,0192 a5=0,0132 14. Procédé selon l'une des revendications 10 à 13, caractérisé en ce que, le filtrage comportant une combinaison linéaire de filtrages, le système matriciel comporte au moins : - une combinaison linéaire correspondante de premières matrices (To) appliquée au vecteur signal du bloc courant, - une combinaison linéaire de deuxièmes matrices (TI) appliquée au vecteur signal du bloc précédent, et - une combinaison linéaire de troisièmes matrices (ri), transposées respectives des deuxièmes matrices (TI), appliquée au vecteur signal du bloc suivant. 15. Programme informatique comportant des instructions pour la mise en oeuvre du procédé selon l'une des revendications 1 à 14 lorsque ce programme est exécuté par un processeur (PROC). 16. Dispositif de traitement d'un signal sous forme de blocs successifs d'échantillons, comportant des moyens de filtrage dans un domaine transformé de sous-bandes, caractérisé en ce que lesdits moyens appliquent en outre : - un traitement d'égalisation (To) à un bloc courant dans le domaine transformé, et - un traitement d'ajustement de filtrage (TI, T'1), dans le domaine transformé, à au moins un bloc adjacent au bloc courant.