Brevet FR2729246A1 — publié le 19960712

PROCEDE DZ CODAGE DE PAROLE A ANALYSE PAR SYNTHESL La présente invention concerne le codage de la parole utilisant analyse par synthèse. La demanderesse a notamment décrit de tels codeurs de parole qu'elle a développés dans ses demandes de brevet européen 0 195 487, 0 347 307 et 0 469 997. Dans un codeur de parole à analyse par synthèse, on effectue une prédiction linéaire du signal de parole pour obtenir les coefficients d'un filtre de synthèse à court terme modélisant la fonction de transfert du conduit vocal. Ces coefficients sont transmis au décodeur, ainsi que des paramètres caractérisant une excitation à appliquer au filtre de synthèse à court terme. Dans la plupart des codeurs actuels, on recherche en outre les corrélations à plus long terme du signal de parole pour caractériser un filtre de synthèse à long terme rendant compte de la hauteur tonale de la parole.Lorsque le signal est voisé, l'excitation comporte en effet une composante prédictible pouvant être représentée par l'excitation passée, retardée de TP échantillons du signal de parole et affectée dun gain gp Le filtre de synthèse à long terme, également reconstitué au décodeur, a alors une fonction de transfert de la forme 1/B(z) avec B(z)=l-gp. TP La partie restante, non prédictible, de l'excitation est appelée excitation stochastique. Dans les codeurs dits CELP ("Code Excited Linear Prediction1), lexcitation stochastique est constituée par un vecteur recherché dans un dictionnaire prédéterminé.Dans les codeurs dits MPLPC ("Multi-Pulse Linear Prediction Coing1), lexcitation stochastique comporte un certain nombre d' impulsions dont les positions sont recherchées par le codeur. En général, les codeurs CELP sont préférés pour les bas débits de transmission, mais ils sont plus complexes à mettre en oeuvre que les codeurs MPLPC. Pour déterminer le retard de prédiction à long terme, on utilise une analyse en boucle ouverte, une analyse en boucle fermée ou une combinaison des deux. Analyse en boucle ouverte est peu exigente en volume de calculs, mais sa précision est limitée. A 11 inverse l'analyse en boucle fermée requiert beaucoup de calculs, mais elle est plus fiable car elle contribue directement à minimiser l'écart pondéré perceptuellement entre le signal de parole et le signal synthétique. Dans certains cas, une analyse en boucle ouverte est d'abord effectuée pour limiter l'intervalle dans lequel analyseur en boucle fermée recherchera le retard de prédiction. Cet intervalle de recherche doit néanmoins rester relativement large car on doit tenir compte de ce que le retard peut varier rapidement. L'invention vise notamment à trouver un bon compromis entre la qualité de la modélisation de la partie à long terme de l'excitation et la complexité de la recherche du retard correspondant dans un codeur de parole. L'invention propose ainsi un procédé de codage à analyse par synthèse d'un signal de parole numérisé en trames successives divisées en nst sous-trames, comprenant les étapes suivantes : analyse par prédiction linéaire du signal de parole pour déterminer des paramètres d'un filtre de synthèse à court terme ; analyse en boucle ouverte du signal de parole pour détecter les trames voisées du signal et pour déterminer, pour chaque trame voisée, un degré de voisement du signal et un intervalle de recherche d'un retard de prédiction à long terme ; analyse prédictive en boucle fermée du signal de parole pour sélectionner, pour certaines au moins des sous-trames des trames voisées, un retard de prédiction à long terme contenu dans l'intervalle de recherche et constituant un paramètre d'un filtre de synthèse à long terme ; et détermination d'une excitation stochastique pour chaque sous-trame, de façon à minimiser un écart pondéré perceptuellement entre le signal de parole et l'excitation stochastique filtrée par les filtres de synthèse à long terme et à court terme. A l'étape d'analyse en boucle ouverte, on détermine l'intervalle de recherche relatif à chaque trame voisée de façon qu'il contienne un nombre de retards dépendant du degré de voisement de ladite trame. Ainsi, le nombre de retards qui sont à tester en boucle fermée est adaptable au mode de voisement de la trame. En général, la largeur de l'intervalle de recherche sera plus faible pour les trames les plus voisées afin de tenir compte de leur plus grande stabilité harmonique. Pour ces trames très voisées, on peut gagner un ou plusieurs bits sur la quantification différentielle du retard dans l'intervalle de recherche, et réaffecter ce ou ces bits gagnés à des paramètres perceptuellement importants comme le gain de prédiction à long terme, ce qui améliore la qualité de restitution de la parole. D'autres particularités et avantages de l'invention apparaîtront dans la description ci-après d'exemples de réalisation préférés, mais non limitatifs, en référence aux dessins annexés, dans lesquels - la figure 1 est un schéma synoptique d'une station de radiocommunication incorporant un codeur de parole mettant en oeuvre l'invention - la figure 2 est un schéma synoptique d'une station de radiocommunication apte à recevoir un signal produit par celle de la figure 1 - les figures 3 à 6 sont des organigrammes illustrant un processus d'analyse LTP en boucle ouverte appliqué dans le codeur de parole de la figure 1 - la figure 7 est un organigramme illustrant un processus de détermination de la réponse impulsionnelle du filtre de synthèse pondéré appliqué dans le codeur de parole de la figure 1 - les figures 8 à 11 sont des organigrammes illustrant un processus de recherche de l'excitation stochastique appliqué dans le codeur de parole de la figure 1. Un codeur de parole mettant en oeuvre l'invention est applicable dans divers types de systèmes de transmission et/ ou de stockage de parole faisant appel à une technique de compression numérique. Dans l'exemple de la figure 1, le codeur de parole-16-fait partie d'une station mobile de radiocommunication. Le signal de parole S est un signal numérique échantillonné à une fréquence typiquement égale à 8kHz. Le signal S est issu d'un convertisseur analogique-numérique 18 recevant le signal de sortie amplifié et filtré d'un microphone 20. Le convertisseur 18 met le signal de parole S sous forme de trames successives elles-mêmes subdivisées en nst sous-trames de 1st échantillons. Une trame de 20 ms comporte typiquement nst=4 sous-trames de lst=40 échantillons de 16 bits à 8kHz.En amont du codeur 16, le signal de parole S peut également être soumis à des traitements classiques de mise en forme tels qu'un filtrage de Hamming. Le codeur de parole 16 délivre une séquence binaire de débit sensiblement plus faible que celui du signal de parole S, et adresse cette séquence à un codeur canal 22 dont la fonction est dintro- duire des bits de redondance dans le signal afin de permettre une détection et/ou une correction d'éventuelles erreurs de transmission. Le signal de sortie du codeur canal 22 est ensuite modulé sur une fréquence porteuse par le modulateur 24, et le signal modulé est émis sur l'interface air. Le codeur de parole 16 est un codeur à analyse par synthèse. Le codeur 16 détermine d'une part des paramètres caractérisant un filtre de synthèse à court terme modélisant le conduit vocal du locuteur, et d'autre part une séquence d'excitation qui, appliquée au filtre de synthèse à court terme, fournit un signal synthétique constituant une estimation du signal de parole S selon un critère de pondération perceptuelle. Le filtre de synthèse à court terme a une fonction de transfert de la forme 1/A(z), avec Les coefficients ai sont déterminés par un module 26 d'analyse par prédiction linéaire à court terme du signal de parole S. Les ai sont lercoefficients de prédiction linéaire du signal de parole S. L'ordre q de la prédiction linéaire est typiquement de l'ordre de 10. Les méthodes applicables par le module 26 pour la prédiction linéaire à court terme sont bien connues dans le domaine du codage de la parole. Le module 26 met par exemple en oeuvre l'algorithme de Durbin-Levinson (voir J. Makhoul : ZLinear Prediction : A tutorial review", Proc. IEEE, Vol.63, N04, Avril 1975, p. 561-580). Les coefficients ai obtenus sont fournis à un module 28 qui les convertit en paramètres de raies spectrales (LSP). La représentation des coefficients de prédiction ai par des paramètres LSP est fréquemnent utilisée dans des codeurs de parole à analyse par synthèse. Les paramètres LSP sont les q nombres cos(2sfi) rangés en ordre décroissant, les q fréquences de raies spectrales (LSF) normalisées fi(lsi5q) étant telles que les nombres complexes exp(2xjfi), avec i=1,3,...,q-l,q+1 et fi+1=0,5, soient les racines du polynôme Q(z) défini par Q(z)=A(z)+z (q+l)A(z 1) et que les nombres complexes exp(2Xjf1), avec i=0,2,4,...,q et fi=0, soient les racines du polynôme Q (z) défini par Q"(z)=A(z)-z-(q+l) A(z-1) Les paramètres LSP peuvent être obtenus par le module de conversion 28 par la méthode classique des polynômes de Chebyshev (voir P. Kabal et R.P.Ramachandran : 1The computation of line spectral frequencies using Chebyshev polynomials", IEEE Trans. ASSP, vol.34, NO 6, 1986, pages 14191426). Ce sont des valeurs de quantification des paramètres LSP, obtenues par un module de quantification 30, qui sont transmises au décodeur pour que celui-ci retrouve les coefficients ai du filtre de synthèse à court terme.Les coefficients ai peuvent être retrouvés simplement, étant donné que: Q(z)=(1+z-1) n (1-2cos(2sfi)z-1+z-) i=1,3,...,q-1 9*(z)=(1-z-1) IL (l-2cos(2z.)z1+z2) et A(z) = [O(z > +0* (z) ]/2 Pour éviter des variations brusques dans la fonction de transfert du filtre de synthèse à court terme, les paramètres LSP font l'objet d'une interpolation avant qu'on en dé duise les coefficients de prédiction ai Cette interpolation est effectuée sur les premières sous-trames de chaque trame du signal. Par exemple, si LSPt et LSPt~l désignent respectivement un paramètre LSP calculé pour la trame t et pour la trame précédente t-1, on prend : LSPt(0)=0,5.LSPt 1+0,5.LSPt, LSPt(l)=0,25.LSPt-1+0,75.LSPt et LSPt(2)=...=LSTt(nst-l)=LSPt pour les sous-trames 0,1,2,...,nst-l de la trame t.Les coefficients ai du filtre 1/A(z) sont alors déterminés, soustrame par sous-trame à partir des paramètres LSP interpolés. Les paramètres LSP non quantifiés sont fournis par le module 28 à un module 32 de calcul des coefficients d'un filtre de pondération perceptuelle 34. Le filtre de pondération perceptuelle 34 a de préférence une fonction de transfert de la forme W(z)=A(z/71)/A(z/y2) où T1 et t2 sont des coefficients tels que Y 2 > (par exemple rl=o,9 et 82=0,6). Les coefficients du filtre de pondération perceptuelle sont calculés par le module 32 pour chaque sous-trame après interpolation des paramètres LSP reçus du module 28. Le filtre de pondération perceptuelle 34 reçoit le signal de parole S et délivre un signal SW pondéré perceptuellement qui est analysé par des modules 36, 38, 40 pour déterminer la séquence d'excitation. La séquence d'excitation du filtre à court terme se compose d'une excitation prédictible par un filtre de synthèse à long terme modélisant la hauteur tonale (pitch) de la parole, et d'une excitation stochastique non prédictible, ou séquence d'innovation. Le module 36 effectue une prédiction à long terme (LTP) en boucle ouverte, c'est-à-dire qu'il ne contribue pas directement à la minimisation de l'erreur pondérée. Dans le cas représenté, le filtre de pondération 34 intervient en amont du module d'analyse en boucle ouverte, mais il pourrait en être autrement : le module 36 pourrait opérer directement sur le signal de parole S ou encore sur le signal S débarrassé de ses corrélations à court terme par un filtre de fonction de transfert A(z). En revanche, les modules 38 et 40 fonctionnent en boucle fermée, c'est-à-dire qu'ils contribuent directement à la minimisation de l'erreur pondérée perceptuellement. Le filtre de synthèse à long terme a une fonction de transfert de la forme l/B(z) avec B(z)=1-gp,z TP où gpoù désigne un gain de prédiction à long terme et TP désigne un retard de prédiction à long terme. Le retard de prédiction à long terme peut typiquement prendre N=256 valeurs comprises entre rmin et rmax échantillons. Une résolution fractionnaire est prévue pour les plus petites valeurs de retard de façon à éviter les écarts trop perceptibles en termes de fréquence de voisement. On utilise par exemple une résolution 1/6 entre rmin=21 et 33+5/6, une résolution 1/3 entre 34 et 47+2/3, une résolution 1/2 entre 48 et 88+1/2, et une résolution entière entre 89 et rmax=142. Chaque retard possible est ainsi quantifié par un index entier compris entre 0 et N-1=255. Le retard de prédiction à long terme est déterminé en deux étapes. Dans la première étape, le module 36 dana- lyse LTP en boucle ouverte détecte les trames voisées du signal de parole et détermine, pour chaque trame voisée, un degré de voisement MV et un intervalle de recherche du retard de prédiction à long terme. Le degré de voisement MV d'une trame voisée peut prendre trois valeurs : 1 pour les trames faiblement voisées, 2 pour les trames modérément voisées, et 3 pour les trames très voisées. Dans les notations utilisées ci-après, on prend un degré de voisement MV=0 pour les trames non voisées.L'intervalle de recherche est défini par une valeur centrale représentée par son index de quantification ZP et par une largeur dans le domaine des index de quantification, dépendant du degré de voisement MV. Pour les trames faiblement ou modérément voisées (Mv=l ou 2) la largeur de l'intervalle de recherche est de N1 index, c'est-à-dire que l'index du retard de prédiction à long terme sera recherché entre ZP-16 et ZP+15 si N1=32. Pour les trames très voisées (MV=3), la largeur de l'intervalle de recherche est de N3 index, c'est-à-dire que l'index du retard de prédiction à long terme sera recherché entre ZP-8 et ZP+7 si N3=16. Une fois que le degré de voisement MV d'une trame a été déterminé par le module 36, le module 30 opère la quantification des paramètres LSP qui ont auparavant été déterminés pour cette trame. Cette quantification est par exemple vectorielle, c'est-à-dire qu'elle consiste à sélectionner, dans une ou plusieurs tables de quantification prédéterminées, un jeu de paramètres quan-tifiés LSPQ qui présente une distance minimale avec le jeu de paramètres LSP fourni par le module 28. De façon connue, les tables de quantification diffèrent suivant le degré de voisement MV fourni au module de quantification 30 par l'analyseur en boucle ouverte 36. Un ensemble de tables de quantification pour un degré de voisement MV est déterminé, lors d'essais préalables, de façon à être statistiquement représentatif de trames ayant ce degré MV. Ces ensembles sont stockés à la fois dans les codeurs et dans les décodeurs mettant en oeuvre l'invention. Le module 30 délivre le jeu de paramètres quantifiés LSPQ ainsi que son index Q dans les tables des quantification applicables. Le codeur de parole 16 comprend en outre un module 42 de calcul de la réponse impulsionnelle du filtre composé du filtre de synthèse à court terme et du filtre de pondération perceptuelle. Ce filtre composé a pour fonction de transfert W(z)/A(z). Pour le calcul de sa réponse impulsionnelle' h= (h(O), h (1), ..., h(lst-l)) sur la durée d'une sous-trame, le module 42 prend pour le filtre de pondération perceptuelle W(z) celui correspondant aux paramètres LSP interpolés mais non quantifiés, c'est-à-dire celui dont les coefficients ont été calculés par le module 32, et pour le filtre de synthèse 1/A(z) celui correspondant aux paramètres LSP quantifiés et interpolés, c'est-à-dire celui qui sera effectivement reconstitué par le décodeur. Dans la deuxième étape de la détermination du retard TP de prédiction à long terme, le module 38 d'analyse LTP en boucle fermée détermine le retard TP pour chaque sous-trame des trames voisées (MV=1, 2 ou 3). Ce retard TP est caractérisé par une valeur différentielle DP dans le domaine des index de quantification, codée sur 5 bits si MV=1 ou 2 (N1=32), et sur 4 bits si MV=3 (N3=16). L'index du retard TP vaut ZP+DP.De façon connue, l'analyse LTP en boucle fermée consiste à déterminer, dans l'intervalle de recherche des retards T de prédiction à long terme, le retard TP qui maximise, pour chaque sous-trame d'une trame voisée, la corrélation normalisée où x(i) désigne le signal de parole pondéré SW de la soustrame auquel on a soustrait la mémoire du filtre de synthèse pondéré (c'est-à-dire la réponse à un signal nul, due à ses états initiaux, du filtre dont la réponse impulsionnelle h a été calculée par le module 42), et yT(i) désigne le produit de convolution u(j-T) désignant la composante prédictible de la séquence d'excitation retardée de T échantillons, estimée par la technique bien connue du répertoire adaptatif ("adoptive codebook"). Pour les retards T inférieurs à la longueur d'une sous-trame, les valeurs manquantes de u(j-T) peuvent être extrapolées à partir des valeurs antérieures. Les retards fractionnaires sont pris en compte en suréchantillonnant le signal u(j-T) dans le répertoire adaptatif. Un suréchantillonnage d'un facteur m est obtenu au moyen de filtres polyphasés interpolateurs. Le gain gp de prédiction à long terme pourrait être déterminé par le module 38 pour chaque sous-trame, en appliquant la formule connue Toutefois, dans une version préférée de l'invention, le gain gp est calculé par le module d'analyse stochastique 40. L'excitation stochastique déterminée pour chaque sous-trame par le module 40 est de type multi-impulsionnelle. Une séquence d'innovation de lst échantillons comprend np impulsions de positions p(n) et d'amplitude g(n). Autrement dit, les impulsions ont une amplitude de 1 et sont associées à des gains respectifs g(n). Etant donné que le retard LTP n'est pas déterminé pour les sous-trames des trames non voisées, on peut prendre un nombre d'impulsions supérieur pour l'excitation stochastique relative à ces sous-trames, par exemple np=5 si MV=1, 2 ou 3 et np=6 si Mu=0. Les positions et les gains calculés par le module 40 d'analyse stochastique sont quantifiés par un module 44. Un module d'ordonnancement des bits 46 reçoit les différents paramètres qui seront utiles au décodeur, et constitue la séquence binaire transmise au codeur canal 22. Ces paramètres sont - - l'index Q des paramètres LSP quantifiés pour chaque trame - le degré MV de voisement de chaque trame - l'index ZP du centre de l'intervalle de recherche des retards LTP pour chaque trame voisée - l'index différentiel DP du retard LTP pour chaque sous-trame d'une trame voisée, et le gain associé - les positions p(n) et les gains g(n) des impulsions de l'excitation stochastique pour chaque sous-trame. Certains de ces paramètres peuvent avoir une importance particulière dans la qualité de restitution de la parole ou une sensibilité particulière aux erreurs de transmission. On prévoit ainsi dans le codeur un module 48 qui reçoit les différents paramètres et qui ajoute à certains d'entre eux des bits de redondance permettant de détecter et/ou de corriger d'éventuelles erreurs de transmission. Par exemple, le degré de voisement MV codé sur deux bits étant un paramètre critique, on souhaite qu'il parvienne au décodeur avec aussi peu d'erreurs que possible. Pour cette raison, des bits de redondance sont ajoutés à ce paramètre par le module 48. On peut par exemple ajouter un bit de parité aux deux bits codant MV et répéter une fois les trois bits ainsi obtenus. Cet exemple de redondance permet de détecter toutes les erreurs simples ou doubles et de corriger toutes les erreurs simples et 75% des erreurs doubles. L'allocation du débit binaire par trame de 20 ms est par exemple celle indiquée dans le tableau I. Dans l'exemple considéré ici, le codeur canal 22 est celui utilisé dans le système paneuropéen de radiocommuni- cation avec les mobiles (GSM). Ce codeur canal, décrit en détail dans la Recommandation GSM 05.03, a été mis au point pour un codeur de parole à 13 kbit/s de type RPE-LTP qui produit également 260 bits par trame de 20 ms. La sensibilité de chacun des 260 bits a été déterminée à partir de tests d'écoute. Les bits issus du codeur source ont été regroupés en trois catégories. La première de ces catégories IA regroupe 50 bits qui sot.codés convolutionnellement sur la base d'un polynôme générateur donnant une redondance d'un demi avec une longueur de contrainte égale à 5.Trois bits de parité sont calculés et ajoutés aux 50 bits de la catégorie IA avant le codage convolutionnel. La seconde catégorie (IB) compte 132 bits qui sont protégés à un taux d'un demi par le même polynôme que la catégorie précédente. La troisième catégorie (II) contient 78 bits non protégés. Après application du code convolutionnel, les bits (456 par trame) sont soumis à un entrelacement. Le module d'ordonnancement 46 du nouveau codeur source mettant en oeuvre l'invention distribue les bits dans les trois catégories en fonction de l'importance subjective de ces bits. paramètres MV=0 MV=1 ou 2 MV=3 quantifiés LSP - 34 34 34 MV + redondance 6 6 6 ZP 8 X 8 DP 20 16 20 20 20 24 positions 80 72 72 impulsions gains 140 100 100 impulsions Total 260 260 0 260 TABLEAU I Une station mobile de radiocommunication apte à recevoir le signal de parole traité par le codeur source 16 est représentée schématiquement sur la figure 2.Le signal radio reçu est d'abord traité par un démodulateur 50 puis par un décodeur canal 52 qui effectuent les opérations duales de celles du modulateur 24 et du codeur canal 22. Le décodeur canal 52 fournit au décodeur de parole 54 une séquence binaire qui, en l'absence d'erreurs de transmission ou lorsque les éventuelles erreurs ont été corrigées par le décodeur canal 52, correspond à la séquence binaire qu'a délivrée le module d'ordonnancement 46 au niveau du codeur 16. Le décodeur 54 comprend un module 56 qui reçoit cette séquence binaire et qui identifie les paramètres relatifs aux différentes trames et sous-trames. Le module 56 effectue en outre quelques contrôles sur les paramètres reçus.En particulier, le module 56 examine les bits de redondance introduits par le module 48 du codeur1 pour détecter et/ou corriger les erreurs affectant les paramètres associés à ces bits de redondance. Pour chaque trame de parole à synthétiser, un module 58 du décodeur reçoit le degré de voisement MV et l'index de Q de quantification des paramètres LSP. Le module 58 retrouve les paramètres LSP quantifiés dans les tables correspondant à la valeur de MV, et, après interpolation, les convertit en coefficients ai pour le filtre de synthèse à court terme 60. Pour chaque sous-trame de parole à synthétiser, un générateur d'impulsions 62 reçoit les positions p(n) des np impulsions de l'excitation stochastique. Le générateur 62 délivre des impulsions d'amplitude unitaire qui sont chacune multipliées en 64 par le gain associé g(n). La sortie de l'amplificateur 64 est adressée au filtre de synthèse à long terme 66. Ce filtre 66 a une structure à répertoire adaptatif. Les échantillons u de sortie du filtre 66 sont mémorisés dans le répertoire adaptatif 68 de façon à être disponibles pour les sous-trames ultérieures. Le retard TP relatif à une soustrame, calculé à partir des index de quantification ZP et DP, est fourni au répertoire adaptatif 68 pour produire le signal u convenablement retardé. L'amplificateur 70 multiplie le signal ainsi retardé par le gain gp de prédiction à long terme.Le filtre à long terme 66 comprend enfin un addi tionneur 72 qui ajoute les sorties des amplificateurs 64 et 70 pour fournir la séquence d'excitation u. Lorsque l'analyse LTP n'pas été effectuée au codeur, par exemple si Mu=0, un gain de prédiction gp nul est imposé à l'amplificateur 70 pour les sous-trames correspondantes. La séquence d'excitation est adressée au filtre de synthèse à court terme 60, et le signal résultant peut encore, de façon connue, être soumis à un post-filtre 74 dont les coefficients dépendent des paramètres de synthèse reçus, pour former le signal de parole synthétique S'. Le signal de sortie S' du décodeur 54 est ensuite converti en analogique par le convertisseur 76 avant d'être amplifié pour commander un haut-parleur 78. On va maintenant décrire, en référence au figures 3 à 6, le processus analyse LTP en boucle ouverte mis en oeuvre par le module 36 du codeur suivant un premier aspect de l'invention. Dans une première étape 90, le module 36 calcule et mémorise, pour chaque sous-trame st=0,l nst-l de la trame courante, les autocorrélations Cst(k) et les énergies retardées Gst(k) du signal de parole pondéré SW pour les retards entiers k compris entre rmin et rmax Les énergies par sous-trame Rosit sont également calculées A l'étape 90, le module 36 détermine en outre, pour chaque sous-trame st, le retard entier Kst qui maximise l'estimation en boucle ouverte Pst (k) du gain de prédiction à long terme sur la sous-trame st, en excluant les retards k pour lesquels l'autocorrélation Cst (k) est négative ou plus petite qu'une petite fraction de l'énergie R0St de la soustrame.L'estimation PSt(k) exprimée en décibels s'écrit Pst (k) = 20.log10[RO,st/ROst-Cst)/G,t(k))] Maximiser Pst (k) revient donc à maximiser l'expression Xst (k) =CSt2 (k) /Gst (k) comme indiqué sur la figure 6. Le retard entier Kst est le retard de base en résolution entière pour la sous-trame st. L'étape 90 est suivie par une comparaison 92 entre une première estimation en boucle ouverte du gain de prédiction global sur la trame courante et un seuil prédéterminé SO typiquement compris entre 1 et 2 décibels (par exemple 50 =1,5 dB).La première estimation du gain de prédiction global est égale à où RO est l'énergie totale de la trame (RO = R00+ R01+... + R0nst-1), et X,t(K,t)= Cst(Kst)/Gst(Kstl désigne le maximum déterminé à l'étape 90 relativement à la sous-trame st. Comme l'indique la figure 6, la comparaison 92 peut être effectuée sans avoir à calculer le logarithme. Si la comparaison 92 montre une première estimation du gain de prédiction inférieure au seuil S0, on considère que le signal de parole contient trop peu de corrélations à long terme pour être voisé, et le degré de voisement MV de la trame courante est pris égal à 0 à l'étape 94, ce qui termine dans ce cas les opérations effectuées par le module 36 sur cette trame. Si au contraire le seuil S0 est dépassé à l'étape 92, la trame courante est détectée comme voisée et le degré MV sera égal à 1, 2 ou 3. Le module 36 calcule alors, pour chaque sous-trame st, une liste I5t contenant des retards candidats pour constituer le centre ZP de linter- valle de recherche des retards de prédiction à long terme. Les opérations effectuées par le module 36 pour chaque sous-trame st (st initialisé à 0 à l'étape 96) d'une trame voisée commencent par la détermination 98 d'un seuil de sélection SE st en décibels égal à une fraction déterminée de l'estimation P St (Kst) du gain de prédiction en décibels sur la sous-trame, maximisée à l'étape 90 (=0,75 typiquement). Pour chaque sous-trame st d'une trame voisée, le module 36 détermine le retard de base rbf en résolution entière pour la suite du traitement. Ce retard de base pourrait être pris égal à l'entier K St obtenu à l'étape 90. Le fait de rechercher le retard de base en résolution fractionnaire autour de Kst permet toutefois de gagner en précision.L'étape 100 consiste ainsi, à rechercher, autour du retard entier K St obtenu à l'étape 90, le retard fractionnaire qui maximise l'expression Cst2/G,t. Cette recherche peut être effectuée à la résolution maximale des retards fractionnaires (1/6 dans l'exemple décrit ici) même si le retard entier K St n'es t pas dans le domaine où cette résolu-tion maximale s'applique.On détermine par exemple le nombre Ast qui maximise CSt2(Kst+6/6)/Gst{Kst+6/6) pour -6 Pour les valeurs fractionnaires T du retard, les autocorrélations Cst (T) et les énergies retardées G St (T) sont obtenues par interpolation à partir des valeurs mémorisées à l'étape 90 pour les retards entiers. Bien entendu, le retard de base relatif à une sous-trame pourrait également être déterminé en résolution fractionnaire dès l'étape 90 et pris en compte dans la première estimation du gain de prédiction global sur la trame. Une fois que le retard de base rbf a été déterminé pour une sous-trame, on procède à un examen 101 des sousmultiples de ce retard afin de retenir ceux pour lesquels le gain de prédiction est relativement important (figure 4), puis des multiples du plus petit sous-multiple retenu (figure 5). A l'étape 102, l'adresse j dans la liste Zist et l'index m du sous-multiple sont initialisés à 0 et 1, respectivement. Une comparaison 104 est effectuée entre le sous-multiple rbf/m et le retard minimal rmin. Le sous-multiple rbf/m est à examiner s'il est supérieur à rmin. On prend alors pour l'entier i la valeur de l'index du retard quantifié ri le plus proche de rbf/m (étape 106), puis on compare, en 108, la valeur estimée du gain de prédiction Pst(ri) associée au retard quantifié ri pour la sous-trame considérée au seuil de sélection SEst calculé à l'étape 98 PSt(ri) = 20.10g10[ROst/(Rost-Cst(ri)/G5t(ri)]] avec, pour les retards fractionnaires une interpolation des valeurs Cst et Gst calculées à l'étape 90 pour les retards entiers.Si Pst(ri) L'examen des sous-multiples du retard de base est terminé lorsque la comparaison 104 montre rbf/m rmax, on effectue le test 118 pour déterminer si l'index mO du plus petit sous-multiple est un multiple entier de n. Dans l'affirmative, le retard n.rbf/mO a déjà été examiné lors de l'examen des sousmultiples de rbf, et on passe directement à l'étape 120 d'incrémentation de l'index n avant d'effectuer de nouveau la comparaison 116 pour le multiple suivant.Si le test 118 montre que mO n'est pas un multiple entier de n, le multiple n.rbf/m0 est à examiner. On prend alors pour l'entier i la valeur de l'index du retard quantifié ri le plus proche de n.rbf/m0 (étape 122), puis on compare, en 124, la valeur estimée du gain de prédiction Pst(ri) au seuil de sélection Sert. Si Pst(ri) Pst(ri) 2 SESt, le retard ri est retenu et on exécute l'étape 126 avant d'incrémenter l'index n à l'étape 120. A l'étape 126, on mémorise l'index i à l'adresse j dans la liste ISt, puis on incrémente d'une unité l'adresse j. L'examen des multiples du plus petit sous-multiple est terminé lorsque la comparaison 116 montre que n.rbf/m0rmax. A ce moment, la liste Ist contient j index de retards candidats. Si on souhaite limiter à jmax la longueur maximale de la liste ISt pour les étapes suivantes, on peut prendre la longueur jst de cette liste égale à min(j,jmax) (étape 128) puis, à l'étape 130, ordonner la liste Ist dans l'ordre des gains Cst(rIst(j))/Gst(rIst (j)) décroissants pour 05j Une fois que les sous-multiples et les multiples ont été examinés et que la liste ISt a ainsi été obtenue (figure 3), le module d'analyse 36 calcule une quantité Ymax déterminant une seconde estimation en boucle ouverte du gain de prédiction à long terme sur l'ensemble de la trame, ainsi que des index ZP, ZPO et ZP1 dans une phase 132 dont le déroulement est détaillé sur la figure 6. Cette phase 132 consiste à tester des intervalles de recherche de longueur N1 pour déterminer celui qui maximise une deuxième estimation du gain de prédiction global sur la trame. Les intervalles testés sont ceux dont les centres sont les retards candidats contenus dans la liste ISt calculée lors de la phase 101. La phase 132 commence par une étape 136 où l'adresse j dans la liste 1St est initialisée à 0. A l'étape 138, on vérifie si l'index Ist(j) a déjà été rencontré en testant un intervalle précé dent centré sur ISt,(j') avec st' Lorsque le test 138 est négatif, on teste l'intervalle centré sur ISt(j) en commençant par l'étape 148 où on détermine, pour chaque sous-trame st', l'index zist, du retard optimal qui maximise sur cet intervalle l'estimation en boucle ouverte Pst(ri) du gain de prédiction à long terme, c'est-àdire qui maximise la quantité Yst.(i)=Cst.2(ri)/Gst,friJ où ri désigne le retard quantifié d'index i pour ISt(j)-Nl/2 #i#Ist(j)+N1/2 et Osi Ensuite, en 150, la quantité Y déterminant la deuxième estimation du gain de prédiction global pour l'intervalle centré sur ISt(j) est calculée selon puis comparée à Ymax, où Ymax représente la valeur à maximiser. Cette valeur Ymax est par exemple initialisée à O en même temps que l'index st à l'étape 96. Si Y#Ymax, on passe directement à l'étape 140 d'incrémentation de l'index j. Si la comparaison 150 montre que Y > Ymax, on exécute l'étape 152 avant d'incrémenter l'adresse j à l'étape 140. A cette étape 152, l'index ZP est pris égal à Ist(j) et les index ZPO et ZP1 sont respectivement pris égaux au plus petit et au plus grand des index ist, déterminés à l'étape 148. A la fin de la phase 132 relative à une sous-trame st, l'index st est incrémenté d'une unité (étape 154) puis comparé, à l'étape 156, au nombre nst de sous-trames par trame. Si st ZP. A l'étape 158, le module 36 détermine le degré de voisement MV, sur la base de la seconde estimation en boucle ouverte du gain exprimée en décibels : Gp=20.loglO(RO/RO- Ymax). On fait appel à deux autres seuils S1 et S2. Si GpSSl, le degré de voisement Mv est pris égal à 1 pour la trame courante. Le seuil S1 est typiquement compris entre 3 et 5 dB; par exemple S1=4 dB. Si Sl S2, on examine la dispersion des retards optimaux pour les différentes sous-trames de la trame courante.Si ZPl-ZP S2). Sinon, si ZP1-ZP2N3/2 ou ZP-ZPO > N3/2, le degré de voisement est pris égal à 2 (si Gp > S2). L'index ZP du centre de l'intervalle de recherche du retard de prédiction pour une trame voisée peut être compris entre 0 et N-1=255, et l'index différentiel DP déterminé pour le module 38 peut aller de -16 à +15 si MV=1 ou 2, et de -8 à +7 si MV=3 (cas N1=32, N3=16). L'index ZP+DP du retard TP finalement déterminé peut donc dans certains cas être plus petit que O ou plus grand que 255. Ceci permet à l'analyse LTP en boucle fermée de porter également sur quelques retards TP plus petits-que rmin où plus grands que rmax. On améliore ainsi la qualité subjective de la restitution des voix dites pathologiques et des signaux non vocaux (fréquences vocales DTMF ou fréquences de signalisation utilisées par le réseau téléphonique commuté).Une autre possibilité est de prendre pour l'intervalle de recherche les 32 premiers ou derniers index de quantification des retards si ZP 240 avec MV=1 ou 2, et les 16 premiers ou derniers index si ZP 248 avec MV=3. Le fait de réduire lintervalle de recherche des retards pour les trames très voisées (typiquement 16 valeurs pour MV=3 au lieu de 32 pour MV=1 ou 2) permet de diminuer la complexité de 1' analyse LTP en boucle fermée effectuée par le module 38 en réduisant le nombre de convolutions yT(i) à calculer suivant la formule (1). Un autre avantage est qu'un bit de codage de l'index différentiel DP est économisé. Le débit de sortie étant constant, ce bit peut être réalloué au codage d'autres paramètres. On peut en particulier allouer ce bit supplémentaire à la quantification du gain de prédiction à long terme gp calculé par le module 40.En effet, une meilleure précision sur le gain gp grâce à un bit de quantification supplémentaire est appréciable car ce paramètre est perceptuellement important pour les sous-trames très voisées (MV=3). Une autre possibilité est de prévoir un bit de parité pour le retard TP et/ou le gain gpv permettant de détecter d'éventuelles erreurs affectant ces paramètres. I1 est possible d'apporter quelques modifications au processus d'analyse LTP en boucle ouverte décrit ci-dessus en référence aux figures 3 à 6. Suivant une première variante de ce processus, les premières optimisations effectuées à l'étape 90 relativement aux différentes sous-trames sont remplacées par une seule optimisation portant sur l'ensemble de la trame. Outre les paramètres CSt C5tAk) et Gst (k) calculés pour chaque sous-trame st, on calcule également les autocorrélations C(k) et les énergies retardées G(k) pour l'ensemble de la trame On détermine alors le retard de base en résolution entière K qui maximise X(k)=C2(k)/G(k) pour rmin s k s rmax. La première estimation du gain comparée à SO à l'étape 92 est alors P(K)=20.log10[RO/(RO-X(K)]]. On détermine ensuite, autour de K, un seul retard de base en résolution fractionnaire rbf et l'examen 101 des sous-multiples et des multiples est effectué une seule fois et produit une seule liste I au lieu de nst listes Ist. La phase 132 est ensuite effectuée une seule fois pour cette liste I, en ne distinguant les soustrames qu'aux étapes 148, 150 et 152. Cette variante de réalisation a pour avantage de réduire la complexité de l'analyse en boucle ouverte. Suivant une seconde variante du processus d'analyse LTP en boucle ouverte, le domaine [rmin, rmax] des retards possibles est subdivisé en nz sous-intervalles ayant par exemple la même longueur (nz=3 typiquement), et les premières optimisations effectuées à 1 ' étape 90 relativement aux différentes sous-trames sont remplacées par nz optimisations dans les différents sous-intervalles portant chacune sur l'ensemble de la trame. On obtient ainsi nz retards de base R1 ,.... Knz I en résolution entière. La décision voisé/non voisé (étape 92) est prise sur la base de celui des retards de base Ki qui procure la plus grande valeur pour la première estimation en boucle ouverte du gain de prédiction à long terme. Ensuite, si la trame est voisée, on détermine les retards de base en résolution fractionnaire par le même processus qu'à l'étape 100, mais en autorisant seulement les valeurs de retard quantifiées. L'examen 101 des sous-multiples et des multiples n'est pas effectué. Pour la phase 132 de calcul de la seconde estimation du gain de prédiction, on prend comme retards candidats les nz retards de base précédemment déterminés. Cette seconde variante permet de se dispenser de l'examen systématique des sous-multiples et des multiples qui sont en général pris en considération grâce à la subdivision du domaine des retards possibles. Suivant une troisième variante du processus d'analyse LTP en boucle ouverte, la phase 132 est modifiée en ce que, aux étapes d'optimisation 148, on détermine d'une part l'index ist, qui maximise CSt,2(ri)/Gst^(ri) pour ISt (j ) -N1 /2si Ist(i)-N3/2si ZP1, mais une quantité Ymax' définie de la même manière que Ymax mais en référence à l'intervalle de longueur réduite Dans cette troisième variante, la détermination 158 du mode de voisement conduit à sélectionner plus souvent le degré de voisement MV=3. On prend également en compte, en plus du gain Gp précédemment décrit, une troisième estimation en boucle ouverte du gain LTP, correspondant à Ymax': Gp'=20.log10[RO/(RO-Ymax')]. Le degré de voisement est MV=1 si GpSS1, MV=3 si Gp&commat; > S2 et MV=2 si aucune de ces deux conditions n'est vérifiée. En augmentant ainsi la proportion de trames de degré MV=3, on réduit la complexité moyenne de l'analyse en boucle fermée et on améliore la robustesse aux erreurs de transmission. Une quatrième variante du processus d'analyse LTP en boucle ouverte concerne surtout les trames faiblement voisées (MV=1). Ces trames correspondent souvent à un début ou à une fin d'une zone de voisement. Fréquemment, ces trames peuvent comporter de une à trois sous-trames pour lesquelles le coefficient de gain du filtre de synthèse à long terme est nul voire-négatif I1 est proposé de ne pas effectuer l'analyse LTP en boucle fermée pour les sous-trames en question, afin de réduire la complexité moyenne du codage. Ceci peut être réalisé en mémorisant à l'étape 152 de la figure 6 nst pointeurs indiquant pour chaque sous-trame st' si l'autocorrélation Cst. correspondant au retard d'index istl est négative ou encore très petite.Une fois que tous les intervalles référencés dans les listes Zist, les sous-trames pour lesquelles le gain de prédiction est négatif ou négligeable peuvent être identifiées en consultant les nst pointeurs. Le cas échéant le module 38 est désactivé pour les sous-trames correspondantes. Ceci n'affecte pas la qualité de l'analyse LTP puisque le gain de prédiction correspondant à ces sous-trames sera de toutes façons quasiment nul. Un autre aspect de l'inventIon concerne le module 42 de calcul de la réponse impulsionnelle du filtre de synthèse pondéré. Le module 38 d'analyse LTP en boucle fermée a besoin de cette réponse impulsionnelle h sur la durée d'une soustrame pour calculer les convolutions YT(i) selon la formule (1). Le module 40 d'analyse stochastique en a également besoin pour calculer des convolutions comme on le verra plus loin. Le fait d'avoir à calculer des convolutions avec une réponse h s'étendant sur la durée dlune sous-trame (lst=40 typiquement) implique une relative complexité du codage, qu'il serait souhaitable de réduire notamment pour augmenter l'autonomie de la station mobile. Dans certains cas il a été proposé de tronquer la réponse impulsionnelle à une longueur inférieure à la longueur d'une sous-trame (par exemple à 20 échantillons), mais ceci peut dégrader la qualité du codage. On propose selon l'invention de tronquer la réponse impulsionnelle h en tenant compte d'une part de la distribution énergétique de cette réponse et d'autre part du degré de voisement MV de la trame considérée, déterminé par le module 36 d'analyse LTP en boucle ouverte. Les opérations effectuées par le module 42 sont par exemple conformes à l'organigramme de la figure 7. La réponse impulsionnelle est d'abord calculée à l'étape 160 sur une longueur pst supérieure à la longueur dune sous-trame et suffisamment grande pour qu'on soit assuré de prendre en compte toute l'énergie de la réponse impulsionnelle (par exemple pst=60 pour nst=4 et lst=40 si la prédiction linéaire à court terme est d'ordre q=10).A l'étape 160, on calcule également les énergies tronquées de la réponse impulsionnelle Les composantes h(i) de la réponse impulsionnelle et les énergies tronquées Eh(i) peuvent être obtenues en filtrant une impulsion unitaire au moyen d'un filtre de fonction de transfert W(z)/A(z) d'états initiaux nuls, ou encore par récurrence Eh (i) = Eh (h (i)] pour O Ensuite le module 42 détermine la plus petite longueur La telle que l'énergie Eh(La-l) de la réponse impulsionnelle tronquée à La échantillons soit au moins égale à une proportion a de son énergie totale Eh(pst-l) estimée sur pst échantillons. Une valeur typique de a est 98%. Le nombre La est initialisé à pst à l'étape 162 et décrémenté d'une unité en 166 tant que Eh a.Eh Pour tenir compte du degré de voisement MV, un terme correcteur (XV) est ajouté à la valeur de La qui a été obtenue (étape 168). Ce terme correcteur est de préférence une fonction croissante du degré de voisement. On peut par exemple prendre A(0)=-5, A(1)=0,(2)=+5 et A(3)=+7. De cette façon, la réponse impulsionnelle h sera déterminée de façon d'autant plus précise que le voisement de la parole est important. La longueur de troncature Lh de la réponse impulsionnelle est prise égale à La si Lasnst et à nst sinon. Les échantillons restants de la réponse impulsionnelle (h(i)=O avec i2Lh) peuvent être annulés. Avec la troncature de la réponse impulsionnelle, le calcul (1) des convolutions yT(i) par le module 38 d'analyse LTP en boucle fermée est modifié de la façon suivante : L'obtention de ces convolutions, qui représente une part importante des calculs effectués, nécessite donc sensiblement moins de multiplications, d'additions et d'adressages dans le répertoire adaptatif lorsque la réponse impulsionnelle est tronquée. La troncature dynamique de la réponse impulsionnelle faisant intervenir le degré de voisement MV permet d'obtenir une telle réduction de complexité sans affecter la qualité du codage. Les mêmes considérations s'appliquent pour les calculs de convolutions effectués par le module 40 d'analyse stochastique.Ces avantages sont particulièrement appréciables lorsque le filtre de pondération per ceptuelle a une fonction de transfert de la forme W(z)=A(z/'yl)/A(z/&gamma;2) avec 0 Un troisième aspect de l'invention concerne le module 40 d'analyse stochastique servant à modéliser la partie non prédictible de l'excitation. L'excitation stochastique considérée ici est de type multi-impulsionnelle. L'excitation stochastique relative à une sous-trame est représentée par np impulsions de positions p(n) et d'amplitudes, ou gains, g(n) (lsnsnp). Le gain gp de prédiction à long terme peut également être calculé au cours du même processus. De façon générale, on peut considérer que la séquence d'excitation relative à une sous-trame comporte nc contributions associées respectivement à nc gains. Les contributions sont des vecteurs îst échantillons qui, pondérés par les gains associés et sommés correspondent à la séquence d'excitation du filtre de synthèse à court terme.Une des contributions peut être prédictible, ou plusieurs dans le cas dun filtre de synthèse à long terme à plusieurs prises ("multi-tap pitch synthesis filter"). Les autres contributions sont dans le cas présent np vecteurs ne comportant que des 0 sauf une impulsion d'amplitude 1. On a donc nc=np si MV=O, et nc=np+l si Mv=l, 2 ou 3. L'analyse multi-impulsionnelle incluant le calcul du gain gp=g(0) consiste, de façon connue, à trouver pour chaque sous-trame des positions p(n) (lsnsnp) et des gains g(n) (0 nnp) qui minimisent 1' erreur quadratique pondérée perceptuellement E entre le signal de parole et le signal synthétisé, donnée par : les gains étant solution du système linéaire g.B=b. Dans les notations ci-dessus : - X désigne un vecteur-cible initial composé des îst échantillons du signal de parole pondéré SW sans mémoire X=(x(O),x(l),...,x(lst-1)), les x(i) ayant été calculés comme indiqué précédemment lors de l'analyse LTP en boucle fermée - g désigne le vecteur ligne composé des np+1 gains g=(g(0)=gpS g(l),...,g(np)) - les vecteurs-ligne Fp(n) (OSn - b désigne le vecteur ligne composé des nc produits scalaires entre le vecteur X et les vecteurs ligne Fpln);; - B désigne une matrice symétrique à nc lignes et nc colonnes dont le terme Bi.j = Fp(i)= Fp(j)T(0#i,j - ()T désigne la transposition matricielle. Pour les impulsions de l'excitation stochastique (lKnsnp=nc-l) les vecteurs Fp(n) sont simplement constitués par le vecteur de la réponse impulsionnelle h décalée de p(n) échantillons. Le fait de tronquer la réponse impulsionnelle comme décrit précédemment permet donc de réduire sensiblement le nombre d'opérations utiles au calcul des produits scalaires faisant intervenir ces vecteurs Fp(n). Pour la contribution prédictible de l'excitation, le vecteur Fp(o)=YTp a pour composantes Fp(0) (i) (OSi Minimiser l'erreur quadratique E définie ci-dessus revient à trouver 1ensemble des positions p(n) qui maximisent la corrélation normalisée b.B-l.bT puis à calculer les gains selon g=b.B -I Mais une recherche exhaustive des positions d'impul sion'nécessiterait un volume de calculs excessif. Pour atténuer ce problème, l'approche multi-impulsionnelle applique en général une procédure sous-optimale consistant à calculer successivement les gains et/ou les positions d'impulsion pour chaque contribution.Pour chaque contribution n (OSn A 1' issue de la dernière itération nc-l, les gains gnc-1(i) (I) sont les gains sélectionnés et l'erreur quadratique minimisée E est égal à l'énergie du vecteur-cible enC-l La méthode ci-dessus donne des résultats satisfaisants, mais elle nécessite l'inversion d'une matrice B n à chaque itération. Dans leur article "Amplitude Optimization and Pitch Prediction in Multipulse Coders" (IEEE Trans. on Acoustics, Speech, and Signal Processing, Vol.37, N 3, Mars 1989, pages 317-327), S. Singhal et B.S. Atal ont proposé de simplifier le problème de l'inversion des matrices Bn en utilisant la décomposition de Cholesky : Bn=Mn.MnT où Mn est une matrice triangulaire inférieure.Cette décomposition est possible parce que B n est une matrice symétrique à valeurs propres positives. L'avantage de cette approche est que l'inversion d'une matrice triangulaire est relativement peu complexe, Bn-1 pouvant être obtenue par B 1=(Mn-)T.Mn- La décomposition de Cholesky et l'inversion de la matrice Mn nécessitent toutefois d'effectuer des divisions et des calculs de racines carrées qui sont des opérations exigeantes en termes de complexité de calcul. L'invention propose de simplifier considérablement la mise en oeuvre de l'optimisation en modifiant la décomposition des matrices Bn de la façon suivante Bn = Ln.RnT = Ln.(Ln.Kn-)T où Kn est une matrice diagonale et Ln est une matrice triangulaire inférieure n'ayant que des 1 sur sa diagonale principale (soit Ln=Mn.Kn1/2 avec les notations précédentes). Compte-tenu de la structure de la matrice Bn, les matrices Ln=Rn.Kn, Rn, Kn et Ln 1 sont construites chacune par simple adjonction dDune ligne aux matrices correspondantes de l'itération précédente B(n,O) o Bn = Bn-1 Ln = Ln 1 B(n,n-1) O B(n,O) . . B(n,n-l) B(n,n) L(n,O) . . . L(n,n-l) 1 o o Rn= Rn-1 o o R(n,0). . R(n,n-1) R(n,n) O . . O K(n o Ln 1 Ln11 o L-1(n.O)... L-(n,n-1) 1 Dans ces conditions, la décomposition de Bn, l'InversIon de LnS obtention de Bn 1=Kn.(Ln-1)T.Ln-1 et le recalcul des gains ne nécessitent qu'une seule division par itération et aucun calcul de racine carrée. L'analyse stochastique relative à une sous-trame d'une trame voisée (MV=1,2 ou 3) peut dès lors se dérouler comme indiqué sur les figures 8 à 11. Pour calculer le gain de prédiction à long terme, l'index de contribution n est initialisé à O à l'étape 180 et le vecteur Fp(0) est pris égal à la contribution à long terme YTp fournie par le module 38. Si n > O, l'itération n commence par la détermination 182 de la position p(n) de l'impulsion n qui maximise la quantité où e=(e(O),...,e(lst-1)) est un vecteur-cible calculé lors de l'itération précédente. Différentes contraintes peuvent être apportées au domaine de maximisation de la quantité cidessus inclus dans l'intervalle [O,lst[.L'invention utilise de préférence une recherche segmentaire dans laquelle la sous-trame d'excitation est subdivisée en ns segments de même longueur (par exemple ns=lO pour lst=40). Pour la première impulsion (n=l), la maximisation de (Fp.eT)2/(Fp.FpT) est effectuée sur l'ensemble des positions possibles p dans la sous-trame. A l'itération n > 1, la maximisation est effectuée à l'étape 182 sur l'ensemble des positions possibles à l'exclusion des segments dans lesquels ont été respectivement trouvées les positions p(l),...,p(n-l) des impulsions lors des itérations précédentes. Dans le cas où la trame courante a été détectée comme non voisée, la contribution n=O est également constituée par une impulsion de position p(O). L'étape 180 comprend alors seulement l'initialisation n=O, et elle est suivie par une étape de maximisation identique à l'étape 182 pour trouver p(O), avec e=e1=x comme valeur initiale du vecteurcible. On remarque que lorsque la contribution n=O est prédictible (MV=1, 2 ou 3), le module 38 d'analyse LTP en boucle fermée a effectué une opération de nature semblable à la maximisation 182, puisqu'il a déterminé la contribution à long terme, caractérisée par le retard TP, en maximisant la quantité (YT.eT)2/(YT.YTT) dans l'intervalle de recherche des retards T, avec e=e~l=X comme valeur initiale du vecteurcible. On peut également, lorsque l'énergie de la contribution LTP est très faible, ignorer cette contribution dans le processus de recalcul des gains. Après l'étape 180 ou 182, le module 40 procède au calcul 184 de la ligne n des matrices L, R et K intervenant dans la décomposition de la matrice B, ce qui permet de compléter les matrices Ln, Rn et Kn définies ci-dessus. La décomposition de la matrice B permet d'écrire pour la composante située à la ligne n et à la colonne j. On peut donc écrire, pour j croissant de O à n-l : L(n,i) =R(n,j),A(j) et, pour j=n : L(n,n) = 1 Ces relations sont exploitées dans le calcul 184 détaillé sur la figure 9. L'index de colonne j est d'abord -initialisé à 0, à l'étape 186.Pour l'index de colonne j, la variable tmp est d'abord initialisée à la valeur de la composante B(n,j), soit : tmp = Fp(n) .Fp T A l'étape 188, l'entier k est en outre initialisé à O. On effectue alors une comparaison 190 entre les entiers k et j. Si k K(n) est pris égal à 1/tmp si tmp*O (étape 198) et à 0 sinon. On constate que le calcul 184 ne requiert qu'au plus une division 198, pour obtenir K(n). En outre, une éventuelle singularité de la matrice B n n'entraîne pas d'instabilités puisqu'on évite les divisions par 0. En référence à la figure 8, le calcul 184 des lignes n de L, R et K est suivi par l'inversion 200 de la matrice Ln constituée des lignes et des colonnes O à n de la matrice L. Le fait que L soit triangulaire avec des 1 sur sa diagonale principale en simplifie grandement l'inversion comme le montre la figure 10. On peut en effet écrire pour O Linv(j')=L-(n,j').L'inversion 200 commence alors par une initialisation 202 de l'index de colonne j' à n-l. A l'étape 204, le terme Linv L(k',j').Linv(k') à Linv(j'), puis on incrémente d'une unité l'entier k (étape 208) avant de réexécuter la comparaison 206. Quand la comparaison 206 montre que k'=n, on compare j' à O (test 210). Si j' > O, on décrémente l'entier j' d'une unité (étape 212) et on revient à l'étape 204 pour calculer la composante suivante. L'inversion 200 est terminée lorsque le test 210 montre que j'=O. En référence à la figure 8 l'inversion 200 est suivie par le calcul 214 des gains réoptimisés et du vecteur-cible E pour itération suivante. Le calcul des gains réoptimisés est également très simplifié par la décomposition retenue pour la matrice B. On peut en effet calculer le vecteur gn=(gn(O),...,gn(n)) g et g,(i')=g(il ) +L )+L~l(n,it).gn(n) pour o#i b(n) sert de valeur d'initialisation pour la variable tmq. A l'étape 216, on initialise également l'index i à o. On effectue ensuite la comparaison 218 entre les entiers i et n. Si i Le calcul 214 est suivi par une incrémentation 228 de l'index n de la contribution, puis par une comparaison 230 entre l'index n et le nombre de contributions nc. Si n L'optimisation des positions et des gains est terminée lorsque n=nc au test 230. La recherche segmentaire des impulsions diminue sensiblement le nombre de positions d'impulsion à évaluer au cours des étapes 182 de la recherche de l'excitation stochastique. Elle permet en outre une quantification efficace des positions trouvées. Dans le cas typique où la sous-trame de lst=40 échantillons est divisée en ns=lO segments de ls=4 échantillons, l'ensemble des positions d'impulsion possibles peut prendre ns!,ls"P/Inp!(ns-np)!1=258 048 valeurs si np=5 (MV=1, 2 ou 3) ou 860 160 si np=6(MV=O), au lieu de lst!/[np!(lst-np)!]=658 008 valeurs si np=5 ou 3 838 380 si np=6 dans le cas où on impose seulement que deux impulsions ne puissent pas avoir la même position.En d'autres termes, on peut quantifier les positions sur 18 bits au lieu de 20 bits si np=5, et sur 20 bits au lieu de 22 si np=6. Le cas particulier où le nombre de segments par soustrame est égal au nombre d'impulsions par excitation stochastique (ns=np) conduit à la plus grande simplicité de la recherche de 1'excitation stochastique, ainsi qu'au plus faible débit binaire (si lst=40 et np=5, il y a 85=32768 ensembles de positions possibles, quantifiables sur 15 bits seulement au lieu de 18 si nus=10). Mais en réduisant à ce point le nombre de séquences d'innovation possibles, on peut appauvrir la qualité du codage. Pour un nombre d'impulsions donné, le nombre des segments peut être optimisé selon un compromis visé entre la qualité du codage et sa simplicité de mise en oeuvre (ainsi que le débit requis). Le cas où ns > np présente en outre l'avantage qu'on peut obtenir une bonne robustesse aux erreurs de transmission en ce qui concerne les positions des impulsions, grâce à une quantification séparée des numéros d'ordre des segments occupés et des positions relatives des impulsions dans chaque segment occupé. Pour une impulsion n, le numéro d'ordre sn du segment et la position relative prn sont respectivement le quotient et le reste de la division euclidienne de p(n) par la longueur ls d'un segment : p(n)=sn.ls+prn (OSsn Au niveau du décodeur, les mots binaires possibles sont stockés dans une table de quantification dans laquelle les adresses de lecture sont les index de quantification reçus. L'ordre dans cette table, déterminé une fois pour toutes, peut être optimisé de façon qu'une erreur de transmission affectant un bit de l'index (le cas d'erreur le plus fréquent, surtout lorsqu'un entrelacement est mis en oeuvre dans le codeur canal 22) ait, en moyenne, des conséquences minimales suivant un critère de voisinage. Le critère de voisinage est par exemple qu'un mot de ns bits ne puisse être remplacé que par des mots voisins", éloignés d'une distance de Hamming au plus égale à un seuil np-28, de façon à conserver toutes les impulsions sauf 6 d'entre elles à des positions valides en cas d'erreur de transmission de l'index portant sur un seul bit.D'autres critères seraient utilisables en substitution ou en complément, par exemple que deux mots soient considérés comme voisins si le remplacement de l'un par l'autre ne modifie pas l'ordre d'affectation des gains associés aux impulsions. A des fins d'illustration. on peut considérer le cas simplifié où ns=4 et np=2, soit 6 mots binaires possibles quantifiables sur nb=3 bits. Dans ce cas, on peut vérifier que la table de quantification présentée au tableau II permet de conserver np-l=l impulsion bien positionnée pour toute erreur affectant un bit de l'index transmis. I1 y a 4 cas d'erreur (sur un total de 18), pour lesquels on reçoit un index de quantification qu'on sait être erroné (6 au lieu de 2 ou 4 ; 7 au lieu de 3 ou 5), mais le décodeur peut alors prendre des mesures limitant la distorsion, par exemple répéter la séquence d'innovation relative à la sous-trame précédente ou encore affecter des mots binaires acceptables aux index 'impossibles' (par exemple 1001 ou 1010 pour l'index 6 et 1100 ou 0110 pour l'index 7 conduisent encore à np-l=l impulsion bien positionnée en cas de réception de 6 ou 7 avec une erreur binaire). Dans le cas général, l'ordre dans la table de quantification des mots peut être déterminé à partir de considérations arithmétiques ou, si cela est insuffisant, en simulant sur ordinateur les scénarios d'erreurs (de façon exhaustive ou par un échantillonnage statistique de type Monte-Carlo suivant le nombre de cas d'erreurs possibles). Pour sécuriser la transmission de l'index de quantification des segments occupés, on peut en outre tirer parti des différentes catégories de protection offertes par le codeur canal 22, notamment si le critère de voisinage ne peut être vérifié de façon satisfaisante pour tous les cas index de quantification mot d'occupation des segments décimal binaire binaire décimal naturel naturel O 000 0011 3 1 001 0101 5 2 010 1001 9 3 011 1100 12 4 100 1010 10 5 101 0110 6 (6) (110) (1001 ou 1010) (9 ou 10) (7) (111) (1100 ou 0110) (12 ou 6) TABLEAU II d'erreurs possibles affectant un bit de l'index. Le module d'ordonnancement 46 peut ainsi mettre dans la catégorie de protection minimale, ou dans la catégorie non protégée, un certain nombre nx des bits de l'index qui, s'ils sont affectés par une erreur de transmission, donnent lieu à un mot erroné mais vérifiant le critère de voisinage avec une probabilité jugée satisfaisante, et mettre dans une catégorie plus protégée les autres bits de l'index. Cette façon de procéder fait appel à un autre ordonnancement des mots dans la table de quantification. Cet ordonnancement peut également être optimisé au moyen de simulations si on souhaite maximiser le nombre nx des bits de l'index affectés à la catégorie la moins protégée. Une possibilité est de commencer par constituer une liste de mots de ns bits par comptage en code de Gray de O à 2ns1, et d'obtenir la table de quantification ordonnée en supprimant de cette liste les mots n'ayant pas un poids de Hamming de np. La table ainsi obtenue est telle que deux mots consécutifs ont une distance de Hamming de np-2.Si les index dans cette table ont une représentation binaire en code de Gray, toute erreur sur le bits de poids le plus faible fait varier 1 index de +1 et entraîne donc le remplacement du mots d'occupation effectif par un mot voisin au sens du seuil np-2 sur la distance de Hamming, et une erreur sur le i-ième bit de poids le plus faible fait aussi varier l'index de +1 avec une probabilité d'environ 21-i En plaçant les nx bits de poids faible de index en code de Gray dans une catégorie non protégée, une éventuelle erreur de transmission affectant un de ces bits conduit au remplacement du mot d'occupation par un mot voisin avec une probabilité au moins égale à (1+1/2+...+1/2nX 1)/nx. Cette probabilité minimale décroît de 1 à (2/nb)(1-1/2nb) pour nx croissant de 1 à nb.Les erreurs affectant les nb-nx bits de poids fort de l'index seront le plus souvent corrigées grâce à la protection que leur applique le codeur canal. La valeur de nx est dans ce cas choisie selon un compromis entre la robustesse aux erreurs (petites valeurs) et un encombrement réduit des catégories protégées (grandes valeurs) Au niveau du codeur, les mots binaires possibles pour représenter l'occupation des segments sont rangés en ordre croissant dans une table de recherche. Une table d'indexage associe à chaque adresse le numéro d'ordre, dans la table de quantification stockée au décodeur, du mot binaire ayant cette adresse dans la table de recherche. Dans l'exemple simplifié évoqué ci-dessus, le contenu de la table de recherche et de la table d'indexage est donné dans le tableau III (en valeurs décimales). La quantification du mot d'occupation des segments déduit des np positions fournies par le module d'analyse stochastique 40 est effectuée en deux étapes par le module de quantification 44. Une recherche dichotomique est d'abord effectuée dans la table de recherche pour déterminer l'adresse dans cette table du mot à quantifier. L'index de quantification est ensuite obtenu à 1 'adresse déterminée dans la table d'indexage puis fourni au module 46 d'ordonnancement des bits. Adresse Table de recherche Table d'indexage 0 3 O 1 5 1 2 6 5 3 9 2 4 10 4 5 12 3 TABLEAU III Le module 44 effectue en outre la quantification des gains calculés par le module 40. Le gain gTp est par exemple quantifié dans l'intervalle [0;1,6], sur 5 bits si MV=1 ou 2 et sur 6 bits si MV=3 pour tenir compte de la plus grande importance perceptuelle de ce paramètre pour les trames très voisées. Pour le codage des gains associés aux impulsions de l'excitation stochastique, on quantifie sur 5 bits la plus grande valeur absolue Gs des gains g(l), . .. ,g(np), en prenant par exemple 32 valeurs de quantification en progression géométrique dans l'intervalle [0;32767], et on quantifie chacun des gains relatifs g(1)/Gs,...,g(np)/Gs dans l'intervalle [-1;+1], sur 4 bits si MV=1, 2 ou 3, ou sur 5 bits si MV=O. Les bits de quantification de Gs sont placés dans une catégorie protégée par le codeur canal 22, de même que les bits de poids fort des index de quantification des gains relatifs. Les bits de quantification des gains relatifs sont ordonnés de façon à permettre leur affectation aux impulsions associées appartenant aux segments localisés par le mot d'occupation. La recherche segmentaire selon l'invention permet en outre de protéger de manière efficace les positions relatives des impulsions associées aux plus grandes valeurs de gain. Dans le cas où np=5 et ls=4, dix bits par sous-trame sont nécessaires pour quantifier les positions relatives des impulsions dans les segments. On considère le cas où 5 de ces 10 bits sont placés dans une catégorie peu ou pas protégée (II) et où les 5 autres sont placés dans une catégorie plus protégée (IB). La distribution la plus naturelle est de placer le bit de poids fort de chaque position relative dans la catégorie protégée IB, de sorte que les éventuelles erreurs de transmission affectent plutôt les bits de poids fort et ne provoquent donc qu'un décalage d'un échantillon pour l'impulsion correspondante. I1 est toutefois judicieux, pour la quantification des positions relatives, de considérer les impulsions dans ordre décroissant des valeurs absolues des gains associés et de placer dans la catégorie IB les deux bits de quantification de chacune des deux premières positions relatives ainsi que le bit de poids fort de la troisième. De cette façon, les positions des impulsions sont protégées préférentiellement lorsqu'elles sont associées à des gains importants, ce qui améliore la qualité moyenne particulièrement pour les sous-trames les plus voisées. Pour reconstituer les contributions impulsionnelles de l'excitation, le décodeur 54 localise d'abord les segments au moyen du mot d'occupation reçu ; il attribue ensuite les gains associés ; puis il attribue les positions relatives aux impulsions sur la base de l'ordre d'importance des gains. On comprendra que les différents aspects de l'invention décrits ci-dessus procurent chacun des améliorations propres, et qu'il est donc envisageable de les mettre en oeuvre indépendamment les uns des autres. Leur combinaison permet de réaliser un codeur de performances particulièrement intéressantes. Dans l'exemple de réalisation décrit dans ce qui précède, le codeur de parole à 13 kbits/s requiert de l'ordre de 15 millions d'instructions par seconde (Mips) en virgule fixe. On le réalisera donc typiquement en programmant un processeur de signal numérique (DSP) du commerce, de même que le décodeur qui ne requiert que de l'ordre de 5 Mips. REVENDICATIONS 1. Procédé de codage à analyse par synthèse d'un signal de parole (S) numérisé en trames successives divisées en nst sous-trames, comprenant les étapes suivantes : - analyse par prédiction linéaire du signal de parole pour déterminer des paramètres d'un filtre (60) de synthèse à court terme - analyse en boucle ouverte du signal de parole pour détecter les trames voisées du signal et pour déterminer, pour chaque trame voisée, un degré de voisement du signal (MV) et un intervalle de recherche d'un retard de prédiction à long terme - analyse prédictive en boucle fermée du signal de parole pour sélectionner, pour certaines au moins des soustrames des trames voisées, un retard de prédiction à long terme contenu dans l'intervalle de recherche et constituant un paramètre d'un filtre de synthèse à long terme (66); et - détermination d'une excitation stochastique pour chaque sous-trame, de façon à minimiser un écart pondéré perceptuellement entre le signal de parole et l'excitation stochastique filtrée par les filtres de synthèse à long terme et à court terme, caractérisé en ce qu'à l'étape d'analyse en boucle ouverte, on détermine l'intervalle de recherche relatif à chaque trame voisée de façon qu'il contienne un nombre de retards (N1, N3) dépendant du degré de voisement de ladite trame. 2. Procédé selon la revendication 1, caractérisé en ce que l'intervalle de recherche du retard de prédiction à long terme contient moins de retards pour les trames ayant le plus grand degré de voisement que pour les autres trames voisées. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'analyse en boucle ouverte relative à une trame comprend la détermination de nst retards de base (kit) qui maximisent chacun une estimation en boucle ouverte du gain de prédiction à long terme sur une sous-trame respective de ladite trame, puis la comparaison entre un premier seuil prédéterminé (SO) et une première estimation en boucle ouverte du gain de prédiction à long terme sur la trame obtenue sur la base des nst retards de base relativement aux sous-trames correspondantes pour détecter si la trame est voisée, en ce que, si la trame est détectée comme voisée, l'analyse en boucle ouverte comprend en outre pour chaque sous-trame la détermination d'une liste (Ist) de retards candidats pour lesquels 1'estimation en boucle ouverte du gain de prédiction sur la sous-trame est supérieure à une fraction déterminée () de l'estimation relative au retard de base pour la sous-trame, en ce qu'on sélectionne dans lesdites listes le retard candidat pour lequel une seconde estimation en boucle ouverte du gain de prédiction à long terme sur la trame est maximale, la seconde estimation en boucle ouverte sur la trame associée à un retard candidat étant obtenue sur la base de nst retards optimaux, compris dans un intervalle de N1 retards centré sur ledit retard candidat, qui, respectivement, maximisent sur ledit intervalle l'estimation en boucle ouverte du gain de prédiction sur les nst sous-trames, en ce que la détermination du degré de voisement de la trame comprend une comparaison entre la seconde estimation maximisée du gain de prédiction sur la trame et au moins un autre seuil prédéterminé (S1,S2), et en ce que l'intervalle de recherche déterminé à l'issue de l'analyse en boucle ouverte est centré sur ledit retard sélectionné. 4. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'analyse en boucle ouverte relative à une trame comprend la détermination d'un retard de base (K) qui maximise une première estimation en boucle ouverte du gain de prédiction à long terme sur ladite trame, puis la comparaison entre un premier seuil prédéterminé (SO) et la première estimation maximisée du gain de prédiction à long terme sur la trame pour détecter si la trame est voisée, en ce que, si la trame est détectée comme voisée, l'analyse en boucle ouverte comprend en outre la détermination d'une liste (I) de retards candidats pour lesquels l'estimation en boucle ouverte du gain de prédiction sur la trame est supérieure à une fraction déterminée (B > de de l'estimation relative au retard de base, en ce qu'on sélectionne dans ladite liste le retard candidat pour lequel une seconde estimation en boucle ouverte du gain de prédiction à long terme sur la trame est maximale, la seconde estimation en boucle ouverte sur la trame associée à un retard candidat étant obtenue sur la base de nst retards optimaux, compris dans un intervalle de N1 retards centré sur ledit retard candidat, qui, respectivement, maximisent sur ledit intervalle l'estimation en boucle ouverte du gain de prédiction sur les nst soustrames, en ce que la détermination du degré de voisement de la trame comprend une comparaison entre la seconde estimation maximisée du gain de prédiction sur la trame et au moins un autre seuil prédéterminé (S1,S2), et en ce que l'intervalle de recherche déterminé à l'issue de l'analyse en boucle ouverte est centré sur ledit retard sélectionné. 5. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'analyse en boucle ouverte relative à une trame comprend la détermination d'un nombre nz de retards de base (K1,..., Knz ) qui maximisent chacun, sur un sous-intervalle respectif de valeurs de retard possibles, une première estimation en boucle ouverte du gain de prédiction à long terme sur ladite trame, puis la comparaison entre un premier seuil prédéterminé (SO) et la plus grande des nz premières estimations maximisées du gain de prédiction à long terme sur la trame pour détecter si la trame est voisée, en ce que, si la trame est détectée comme voisée, on sélectionne, parmi nz retards candidats obtenus à partir des nz retards de base, le retard candidat pour lequel une seconde estimation en boucle ouverte du gain de prédiction à long terme sur la trame est maximale, la seconde estimation en boucle ouverte sur la trame associée à un retard candidat étant obtenue sur la base de nst retards optimaux, compris dans un intervalle de -N1 retards- centré sur ledit retard candidat, qui, respectivement, maximisent sur ledit intervalle l'estimation en boucle ouverte du gain de prédiction sur les nst soustrames, en ce que la détermination du degré de voisement de la trame comprend une comparaison entre la seconde estimation maximisée du gain de prédiction sur la trame et au moins un autre seuil prédéterminé (S1,S2), et en ce que l'intervalle de recherche déterminé à l'issue de l'analyse en boucle ouverte est centré sur ledit retard sélectionné. 6. Procédé selon l'une quelconque des revendications 3 à 5, caractérisé en ce que si la seconde estimation maximisée du gain de prédiction sur une trame voisée est supérieure à un des seuils (S2), on détermine si les nst retards optimaux sont compris dans un intervalle centré sur le retard sélectionné et contenant un nombre de retards N3 inférieur à N1 et, dans laffirmative, on attribue à la trame un degré de voisement pour lequel l'intervalle de recherche du retard de prédiction à long terme contient N3 retards, l'intervalle de recherche contenant N1 retards pour au moins un autre degré de voisement. 7. Procédé selon l'une quelconque des revendications 3 à 5, caractérisé en ce que lors de la maximisation de la seconde estimation en boucle ouverte du gain de prédiction à long terme sur une trame voisée, on calcule en outre une troisième estimation en boucle ouverte du gain sur la trame sur la base de nst retards, compris dans un intervalle centré sur le retard sélectionné et contenant un nombre N3 de retards inférieur à N1, qui, respectivement, maximisent sur ledit intervalle de N3 retards leestimation en boucle ouverte du gain de prédiction sur les nst sous-trames, et en ce qu'on attribue à la trame un degré de voisement pour lequel l'intervalle de recherche contient N3 retards si ladite troisième estimation dépasse un seuil prédéterminé (S2), l'intervalle de recherche contenant N1 retards pour au moins un autre degré de voisement. 8. Procédé selon la revendication 3 ou 4, caractérisé en ce que les retards candidats d'une liste sont choisis parmi les sous-multiples du retard de base associé à ladite liste et parmi les multiples du plus petit desdits sousmultiples pour lequel l'estimation en boucle ouverte du gain de prédiction est supérieure à ladite fraction déterminée de l'estimation relative au retard de base. 9. Procédé selon la revendication 8, caractérisé en ce que les retards de prédiction à long terme peuvent correspondre à des nombres entiers ou fractionnaires d'échantillons du signal de parole, en ce qu'on détermine les retards de base (rbf) en résolution fractionnaire pour chercher les sous-multiples et les multiples à inclure dans une liste des retards candidats, et en ce que les retards de base sont déterminés en résolution entière pour évaluer les premières estimations en boucle ouverte du gain de prédiction sur une trame. 10. Procédé selon lune quelconque des revendications 3 à 9, caractérisé en ce qu'on n'effectue pas l'analyse prédictive en boucle fermée relativement à chaque sous-trame pour laquelle lautocorrélation