Brevet FR2618277A1 — publié le 19890120

Prédicteur adaptatif. La présente invention concerne les prédicteurs adaptatifs de modélisation de signaux, utilisables notamment pour la compression ou l'identification de signaux. L'invention trouve une application particuliè- rement importante dans la compression'des signaux en vue de leur mémorisation ou de leur transmission. De tels prédicteurs adaptatifs sont par exemple prévus dans les installations de transmission utilisant le codage MICDA dit aussi ADPCM. L'invention concerne plus particulièrement les prédicteurs dits autorégressifs à moyenne ajustée, en l'abréviation ARMA, constituant un filtre linéaire qui fournit, en réponse à un signal s(n) de puissance géné- ralement très variable dans le temps, une erreur de prédiction e(n) selon l'équation: p.q e(nY = s(n) - als(n-i) - E bie(n-j) (1) e(n) = ú s(n-) - E benj i=1 j=1 Le signal s(n) est en général fortement non sta- tionnaire et les coefficients a et b sont en conséquence des fonctions du temps. Le premier terme en E de l'équation (1) corres- pond à un prédicteur purement auto-régressif ayant le caractère d'un filtre transverse, souvent suffisant comme dans le cas d'un signal aléatoire à bande étroite qui est assimilable à la réponse d'un filtre auto- régressif linéaire à une excitation aléatoire blanche. Par contre, dans le cas de beaucoup de signaux de nature différente, comme pour des signaux de parole, on se prive d'un gain notable de prédiction en omettant la partie MA. En conséquence, on sera souvent amené à utiliser un prédicteur ARMA répondant à la formule (1). Mais l'utilisation d'un prédicteur ARMA se heurte, dans le cas de signaux à bande étroite ou, plus généralement, quand le spectre présente des flancs de montée raide, à des problèmes de stabilité. Pour mieux faire apparaître ces problèmes, il peut être utile de rappeler tout d'abord quelques données concernant la prédiction adaptative ARMA. Pour simplifier les équations, on utilisera dans la suite une notation vectorielle: AT = (a1, a2,..., ap)T; BT = (b1, b2,..., bq) (2) S(n) = (s(n-1),.., s(n-p))T; E(n) = (e(n-1),..., e(n-q))T (3) o les vecteurs S(n) et E(n) correspondent aux échantil- lons passés des signaux considérés. Le modèle de prédicteur (1) s'écrit alors: s(n) = ATS(n) + BTE(n) + e(n) (4) On obtient une compression du signal en conser- vant comme informations les vecteurs A et B,'éventuel- lement numérisés, et l'erreur e(n), sous une forme qui dépend de l'application, et en oubliant le signal s(n) originel. Cette méthode de compression est notamment appliquée dans les systèmes de codage numérique à débit réduit du signal téléphonique, tel que le codage MICDA. Lorsque le signal s(n) est non stationnaire, les vecteurs A et B et la puissance de l'erreur e(n) sont des fonctions du temps n. La prédiction adaptative doit identifier en per- manence, au cours du temps, les grandeursA(n), B(n), e(n) sur la base du seul signal s(n). La quantité T- T s(n) = A(n)TS(n) + B(n)TE(n) (6) représente une prédiction du signal basée sur des obser- vations du passé et e(n) = s(n) - s(n) (7) est l'erreur de prédiction. Pour identifier les paramètres A et B, des algo- rithmes adaptatifs ont été proposés dans la littérature. Ils peuvent s'écrire sous la forme: A(n+1)=A(n)+a f[e(n); e(n-1);...; s(n) , s(n-1),...](8) B(n+1)=B(n)+R g(e(n); efn-1);...; s(n), s(n-1),...](9) Ces algorithmes visent à minimiser une moyenne quadratique de e(n). La figure 1 donne une représentation très schématique d'une constitution possible d'un prédicteur ARMA. Le signal modélisé s(n) est appliqué à un addi- tionneur 10 dont l'entrée soustractive reçoit l'estima- tion (n). La sortie de l'additionneur, représentant l'erreur de prédiction en, est appliquée à des circuits adaptatifs 12 et 14 d'élaboration des vecteurs B(n) et A(n) qui peuvent avoir l'une quelconque des constitu- tions bien connues. Les signaux s(n) et e(n) sont appliqués, après retard z-1, respectivement aux filtres 14 et 12 dont les vecteurs de coefficients sont A(n) et B(n). Les sorties des deux filtres sont appliquées sur les deux entrées d'un additionneur 16 qui fournit s(n). Lorsque les paramètres du prédicteur sont fixes, le prédicteur (1) est un filtre linéaire récursif dont la sortie e(n) associée à l'entrée s(n) peut s'écrire à l'aide des polynomes: p q A(z) = E az-; B(z) = E bJz- (10) i=1 j=l selon 1-A(z) e(n) = H(z).s(n) avec H(z) = (11) 1+B(z) Le problème de stabilité dans la mise en oeuvre d'un tel prédicteur adaptatif fait intervenir la fonction H(z) et aussi H 1(z): a) la stabilité de la fonction H(z) est liée à la forme de l'algorithme adaptatif (9). Les zéros du polynome (1+B(z)) doivent demeurer à l'intérieur du cercle unité. Faute d'obtenir cette propriété, l'erreur de prédiction e(n) donnée par (11) risque de diverger; b) la stabilité de la fonction de transfert inverse H 1(z) est liée uniquement à la position des zéros du polynome (1-A(z)) qui doivent demeurer à l'intérieur du cercle unité. Faute d'obtenir cette propriété, la restitution du signal's(n) à partir des paramètres A, B et de l'excitation e(n) risque de diverger, dès que cette dernière comporte un peu de bruit, ce qui est le cas général. Ces deux critères sont d'autant plus difficiles à remplir que le signal s(n) est à bande étroite par rapport à la largeur du canal et que son spectre présen- te des flancs de môntée raide. Une solution permettant d'écarter le problème (a) consiste à adopter B=0 mais, comme on l'a vu plus haut, cette solution n'est pas acceptable dans le cas d'un signal de parole car elle se traduit par une perte.de gain de prédiction appréciable. Pour remplir les critères il est possible de prévoir, dans le prédicteur, des moyens de calcul et de test des zéros de 1-An(z) et de l+Bn à chaque itération d'ordre n. Mais le calcul est très lourd dès que les ordres p et q des parties auto-régressive et moyenne ajustée dépassent 2, ce qui est le cas sauf lorsqu'on se contente d'une prédiction rudimentaire. Une autre solution encore consiste à réaliser le prédicteur sous forme d'une cascade de cellules ARMA d'ordre 2,2, dont la stabilité peut facilement être testée par des circuits simples, à l'aide du triangle de stabilité. Mais en contrepartie, ce montage en cascade rend difficile la minimisation quadratique de e(n) de manière algorithmique. D'autres solutions encore existent: toutes se traduisent soit par une complexité généralement inaccep- table des calculs, soit par un effet défavorable sur la valeur résiduelle de e(n) ou le gain de prédiction. 2618Z77 La présente invention vise & fournir un prédic- teur répondant mieux que ceux antérieurement connus aux exigences de la pratique permettant tout & la fois d'arriver à la stabilité de H(z) et de H-1(z) et d'éviter une effet défavorable sur le gain de prédiction et d'éviter un accroissement appréciable de la complexi- té des calculs et des circuits. Ce résultat a été atteint & partir d'une analyse mathématique des algorithmes de prédiction et des effets recherchés et de la constatation qu'il est possible, sans pour autant affecter de façon défavorable les résultats de prédiction, d'adopter, au lieu d'un critère de prédiction pure, un critère résultant d'un compromis entre la recherche d'un gain de prédiction maximale et la nécessité de rendre stable H(z) et H1 (z). L'analyse sera faite à titre d'exemple pour le cas fréquent o les algorithmes de prédiction adaptative (8) et (9) destinés à l'élaboration des paramètres A et B sont destinés à minimiser par itération la quantité: 2(2 Jc (A,B) = E(e 2(n)) (12) pour un signal s(n) donné. On peut citer l'algorithme classique dit LMS récursif qui s'écrit: A(n+1) = A(n) + ae(n)S(n), avec a>0 (13) B(n+1) = B(n) + ge(n)E(n), avec 0>0 (14) Lors de la convergence, les incréments ont une moyenne nulle. L'erreur de sortie e(n) est donc décorrélée des erreurs passées e(n-j) et des entrées passées s(n-i). Mais lorsque s(n) a un spectre "raide" (ce qui est le cas des tonalités ou des sons voisés dans la parole) les paramètres A(n) et B(n) deviennent très voisins de la frontière de stabilité: les fluctations dues à l'algorithme déterminent des instabilités pour la fonction de transfert H(z) et la fonction de transfert inverse H-1(z). Il faut remarquer que, dans le cas des signaux à bande étroite, on démontre aussi que le paramètre B(n) traverse inexorablement et de manière quasi-périodique, la frontière de stabilité: la fonction de transfert directe H(z) est donc elle aussi périodiquement instable. Pour la stabiliser, on introduit conformément à l'invention, une fonction de coût Js similaire à Jc et représentant la puissance du signal prédit: A 2 Js(A,B) = E(s2(n)) (15) En minimisant la fonction d'espérance mathéma- tique Js on stabilise le prédicteur pour une entrée s(n) donnée. Et, puisque le signal restitué par le prédicteur inverse, pour une entrée e(n) donnée, ne diffère de s(n) que par e(n), en minimisant Js on stabilise aussi le prédicteur inverse. Cette constatation conduit à abandonner la recherche d'une puissance Jc minimum en faveur de la recherche d'une valeur minimum du coût JN défini par: JN(A,B) = Jc(A,B) + AJs(A,B), avec A>O (16) On réalise ainsi un compromis entre la maximisation du gain de prédiction (qui exigerait la recherche d'un minimum de Jc) et la stabilité des fonc- tions H-1(z) et H(z) (qui impliquerait un minimum de Js). En conséquence, l'invention propose un prédic- teur auto-régressif à moyenne ajustée de modélisation de signaux, comportant un filtre linéaire fournissant un signal s(n), en réponse à une erreur de prédiction e(n): s(n) = ATS(n) + BTE(n) + e(n) (4) o oT: AT et BT sont les vecteurs des coefficients du filtre, S(n) est le vecteur des échantillons passés de s(n), et E(n) est le vecteur des erreurs passées e(n), caractérisé par des moyens pour élaborer les vecteurs AT et BT par un algorithme d'adaptation faisant intervenir, dans le terme d'incrémentation: A e(n) - ± s(n) 1+X o A est une valeur fixe comprise entre 0,05 et 0,3 et e(n) est l'erreur de prédiction, égale à la différence entre le signal s(n) et sa valeur s(n) prédite. On choisira toujours une valeur de A nettement inférieure à 1, toute valeur approchant de 1 et surtout le dépassant se traduisant par une dégradation des performances du prédicteur, car il y a pratiquement suppression de la prédiction. On cherchera à s'approcher le plus possible de la valeur minimum de A pour laquelle la stabilité est assurée pour le type de signal à tra'i- ter le plus défavorable. A sera toujours pratiquement supérieur à 0,05 et inférieur à 0,25. Pour des prédicteurs purement auto-régressifs, ayant un nombre de coefficients courant, c'est-à-dire compris entre deux et dix, typiquement six, on peut estimer qu'une valeur de AA (c'est-à-dire intervenant dans A) proche de 0,1 constitue un compromis satisfai- sant. Une détermination plus précise de AA est possible à l'aide d'abaques, le coefficient AA diminuant'lorsque le rapport entre la largeur de bande et la fréquence d'échantillonnage augmente. Au delà d'une largeur de bande W = 1/16 de la fréquence d'échantillonnage on peut fixer AA =. Pour des prédicteurs purement MA, ayant un nombre de coefficients assez élevé, supérieur ou égal à 6, la valeur de AB (intervenant dans B) doit être choisie plus grande que dans le casAR (typiquement 0,3). En effet, ces prédicteurs modélisent mal les signaux à bande étroite. Dans la pratique cependant, ce n'est pas un inconvénient car on n'utilise pas de prédicteurs purement MA. Pour des prédicteurs ARMA dont la partie AR a au moins deux coefficients, on conservera avantageuse- ment ia valeur AA précédemment proposée pour stabiliser la partie AR et on réduira AB à la même valeur que AA, typiquement AA = AB = 0,1. L'invention propose donc de substituer, dans le cas o on utilise l'algorithme LMS récursif, les formule (13 bis) ou (14 bis) ci-après aux formules (13) et (14): A(n+1) = A(n) + a[e(n) - A s(n)]S(n) (13 bis) 1A B(n+1) = B(n) + 0[e(n) - A s(n)]E(n) (14 bis) 1+RB Naturellement, on pourrait de même adopter certains autres algorithmes adaptatifs mettant en oeuvre une fonction autre que le gradient (par exemple l'algo- rithme du signe). Dans le cas de l'algorithme du signe, on aurait par exemple: A(n+1) = A(n) + signe[e(n) - s(n)]S(n) (13 ter) 1+ A B(n+1) = B(n) + signe[e(n) - A s(n)]E(n) (14 ter) 1+RB L'invention sera mieux comprise à la lecture de la description qui suit d'un mode particulier de réali- sation donné à titre d'exemple non limitatif. La des- cription se réfère aux dessins qui l'accompagnent, dans lesquels: - la figure 1, déjà mentionnée, est un schéma de principe d'un prédicteur adaptatif de type connu; - la figure 2, similaire & la figure 1, montre une constitution possible d'un prédicteur suivant l'invention. Le prédicteur montré schématiquement en figure 2 ne se différencie de celui de la figure I que par l'adjonction des éléments permettant de passer des algorithmes (13) et (14) aux algorithmes (13 bis) et (14 bis) respectivement. Les éléments identiques à ceux de 3 la figure 1 sont désignés par le même numéro de référen- ce. Le dispositif comporte de plus un soustracteur 20 qui reçoit le signal e(n) sur son entrée additive et qui reçoit, sur son entrée soustractive, le signal s(n) amplifié dans le rapport A/1+A par un circuit 18. On voit que le circuit supplémentaire nécessaire est de constitution extrêmement simple et que le coût du prédicteur dans son ensemble n'est pratiquement pas modifié. Pour déterminer l'effet obtenu par la mise en oeuvre de l'invention, on considérera successivement l'accroissement de stabilité obtenu d'abord par le passage de (13) à (13 bis) puis de (14) à (14 bis). Si on appelle respectivement Ac et AN les points d'équilibre du système classique (13) et du système stabilisé (13 bis), on montre que: N Ac A = 1 + On démontre théoriquement que, pour AA assez grand donc du type de signal s(n)), l'homothétie correspon- dante (1+AA)1 ramène la prédiction à l'intérieur du A domaine de stabilité si Ac était à la limite de ce domaine de stabilité. Avec la même convention et sous les hypothèses que l'ordre q est assez élevé et que AB est très infé- rieur- à un, on démontre encore que: Bc BC BN =, I + ÀB La conclusion théorique est identique: la contraction (1+AB)- a tendance à ramener le prédicteur à l'intérieur du domaine de stabilité si Bc était à la limite de ce domaine. Dans les adaptations (13 bis) et (14 bis), les évolutions des parties AR et MA sont couplées du fait du e(n) qui intervient dans les deux incréments. Ceci permet une meilleure optimisation du prédicteur que celle qui consiste à optimiser séparément les parties AR et MA. Le gain de prédiction s'en trouve amélioré et la valeur \B peut être réduite. 1 1 REVENDICATIONS 1. Prédicteur auto-régressif de modélisation de signaux, comportant un filtre linéaire fournissant, en réponse à une erreur de prédiction e(n), un signal s(n): s(n) = ATS(n) + BTE(n) +.e(n) o AT et BT sont les vecteurs des coefficients du filtre, S(n) est le vecteur des échantillons passés de s(n), et E(n) est le vecteur des erreurs passées e(n), caractérisé par des moyens (12, 14, 18, 20) pour élabo- AT et -T rer les vecteurs AT etBT par un algorithme d'adaptation faisant intervenir, dans les termes d'incrémentation: e(n) - s(n). et e(n) - + s(n) o AA et B sont des valeurs fixes comprises entre 0,05 et 0,3 et e(n) est l'erreur de prédiction, égale à la différence entre le signal s(n) et sa valeur s(n) prédite. 2. Prédicteur selon la revendication 1, purement auto-régressif, ayant de deux à dix coefficients, carac- térisé en ce que AA est proche de 0,1. 3. Prédicteur selon la revendication 1, à moyen- ne ajustée, caractérisé en ce que AA et >B sont supé- rieurs à 0,05 et inférieurs à 0,25. 4. Prédicteur selon la revendication 1 ou 3, caractérisé en ce que AA et ÀB ont la même valeur. 5. Prédicteur selon l'une quelconque des reven- dications précédentes, caractérisé en ce que lesdits moyens mettent en oeuvre l'algorithme du gradient ou du signe.