FR 2696297 A1 19940401 FR 9211467 A 19920925 "Procédé pour déterminer la transmittance d'un circuit de filtrage prévu pour transformer la réponse impulsionnelle d'un filtre en une réponse à phase minimale et filtre mettant en oeuvre ce procédé" La présente invention concerne un procédé pour déterminer la transmittance d'un circuit de filtrage prévu pour transformer la réponse impulsionnelle d'un filtre en une réponse à phase minimale. Elle vise également un circuit de filtrage obtenu par mise en oeuvre de ce procédé. Une réponse est dite à phase minimale lorsque sa transformée en z n' a de zéro qu'à l'intérieur du cercle unité dans le plan complexe. Jusqu'à présent, on effectuait soit une recherche systématique de tous les zéros de la fonction de transfert associée à la réponse impulsionnelle à transformer et on réalisait l'extraction de ceux dont le module est supérieur ou égal à 1, soit une recherche des seuls zéros de module supérieur ou égal à 1, comme cela est divulgué dans l'article "adaptive adjustment of receiver for distorted digital signals" (réglage adaptatif d'un récepteur de signaux numériques distordus) IEE Proc Part F (Août 1984) Vol 131, pp 526-536. Cependant, ces méthodes ont une précision limitée et leur complexité augmente avec la longueur du filtre En outre, le nombre d'opérations à effectuer est une variable aléatoire et le temps de convergence est variable. Le but de l'invention est de remédier à ces inconvénients en proposant un procédé pour déterminer la transmittance d'un circuit de filtrage associé à un filtre présentant une réponse impulsionnelle H(z), ce circuit de filtrage étant prévu pour que la réponse impulsionnelle globale G(z) du filtre et du circuit de filtrage soit à phase minimale, ce procédé présentant un temps de convergence fixe et étant de complexité réduite. Suivant l'invention, le procédé comprend: une détermination de la réponse en fréquence H(f) du filtre par transformée de Fourier de la réponse impulsionnelle H(z) de ce filtre, une détermination de la partie réelle a(f) du cepstre H(f) = -ln(H(f)) de la réponse impulsionnelle H(f) du filtre pour extraire de cette partie réelle a(f) une partie paire p(f) et une partie impaire q(f), une détermination du cepstre G(f) = -ln(G(f)) de la réponse globale en fréquence à partir de ces parties paire p(f) et impaire q(f) du cepstre H(f) de la réponse impulsionnelle du filtre, une détermination de la réponse globale en fréquence G(f) à partir du cepstre G(f) de cette réponse, une détermination de la réponse théorique en fréquence C(f) du circuit de filtrage au moyen de la réponse globale en fréquence G(f) et de la réponse en fréquence H(f) du filtre, un calcul de la transmittance théorique C(z) du circuit de filtrage par transformée de Fourier inverse de la réponse en fréquence C(f) de ce circuit et, une détermination de la transmittance P(z) du circuit de filtrage par troncature de la transmittance C(z) dont ne sont conservés qu'un nombre prédéterminé de coefficients. Ainsi, le procédé selon l'invention permet de concentrer l'énergie dans les premiers échantillons de la réponse impulsionnelle fournie en sortie et la transforme en réponse à phase minimale. Cé procédé met en oeuvre le principe du cepstre qui est fondé sur la causalité de la réponse impulsionnelle G(z)= - ln(G(z)) qui est obtenue lorsque la réponse impulsionnelle G(z) est à phase minimale. Selon une forme préférée de mise en oeuvre du procédé selon l'invention, on procède à une estimation de la réponse impulsionnelle globale G(z) à partir d'un développement limité du cepstre G(z). Le procédé selon l'invention peut comprendre en outre l'obtention d'une première estimation Go(z) de la réponse globale à phase minimale G(z) comportant: un calcul des coefficients de la réponse impulsionnelle à phase minimale G(z) par transformée de Fourier discrète inverse de la réponse en fréquence à phase minimale G(f), et la conservation d'un nombre prédéterminé de coefficients significatifs pour obtenir ladite estimation Go(z). Suivant un autre aspect de l'invention, il est proposé un circuit de filtrage numérique transversal correcteur de phase, prévu pour transformer une réponse impulsionnelle d'un fitre en une réponse à phase minimale Des moyens de calcul algorithmique sont associés au cicuit de filtrage pour en adapter les coefficients par mise en oeuvre du procédé selon l'invention. D'autres particularités et avantages de l'invention apparaîtront encore dans la description ci-après. Aux dessins annexés donnés à titre d'exemples non limitatifs: La figure 1 est un schéma-bloc d'un dispositif de réception incluant un circuit de filtrage selon l'invention, la figure 2 illustre la variation de l'énergie de la réponse impulsionnelle globale d'un filtre et du circuit de filtrage selon l'invention en fonction du temps, les figures 3 à 5 illustrent l'amélioration apportée par un circuit de filtrage selon l'invention en représentant l'évolution de taux d'erreur binaire à la sortie d'un récepteur en fonction du rapport de l'énergie moyenne par bit en sortie du filtre sur la densité spectrale unilatérale de puissance d'un bruit blanc additif gaussien pour des configurations de canal suivantes: * figure 3: canal à deux trajets d'égales puissances, * figure 4: canal à deux trajets de puissances inégales, * figure 5: canal à trois trajets d'égales puissances. on va maintenant décrire un exemple de mise en oeuvre du procédé selon l'invention au sein d'un dispositif de réception. Le récepteur 10, représenté sous la forme d'un schéma- bloc en figure 1, est agencé pour recevoir des symboles Sn émis en entrée d'un filtre 2, par exemple un canal de transmission 2, présentant une réponse impulsionnelle H(z) échantillonnée à une fréquence prédéterminée 1/T, o T est une période d'échantillonnage A ces symboles est ajouté au cours de la transmission, du bruit B Le récepteur 10 selon l'invention comprend un circuit de filtrage 3 prévu pour transformer la réponse impulsionnelle globale en une réponse à phase minimale et situé en amont d'un détecteur 4 de type DFSE (Decision feedback Sequence Estimator: Estimateur de séquence à Réaction de Décision)qui génère des symboles de réception Sn, des moyens 5 pour estimer le filtre et calculer une estimation H(z) de la réponse impulsionnelle du filtre, et des moyens 1 pour calculer un algorithme d'adaptation des coefficients du circuit de filtrage 3 qui matérialise le procédé selon l'invention Les moyens de calcul 1 génèrent une estimation P(z) de la transmittance P(z) du circuit de filtrage 3 qui est un filtre numérique transversal de réponse impulsionnelle finie approchant un filtre correcteur de phase idéal de réponse impulsionnelle infinie, et une estimation G(z) de la réponse impulsionnelle à phase minimale en sortie du circuit de filtrage 3. Soient: H(z) la transformée en z de la réponse impulsionnelle du filtre échantillonnée à la fréquence 1/T, G(z) la transformée en z de la réponse globale à phase minimale, C(z) la transmittance en z du circuit de filtrage. K K H h z G(z)= Zg hn ZG(z)= gn C(z)=G(z) n=O n= O o K est la longueur du filtre. Si on considère dans un premier temps le cas o H(z) n'a pas de zéro sur le cercle unitaire, p zéros à l'intérieur de ce cercle notés ai, q zéros à l'extérieur notés Tj: ( 1) H(z)=ho 1 ( 1 aiz-') f I ( 1-Pj 6 ') jk I Mj>l K=p+q i=l j=l La propriété de phase minimale sur G(z) impose pour G(z) et C(z) les expressions suivantes: ( 2) G(z)= hg H( 1-z-') ( 1 z-') j=li=l jl O avec le sigle * représentant l'opération de conjugaison. La transmittance théorique C(z) du circuit de filtrage passe-tout et de gain unitaire peut s'exprimer de la façon suivante: q r I ( 1 -Cj Z) fla à ( 1z) j=l j=l pi Les pôles de la transmittance C(z) étant à l'extérieur du cercle unitaire, on peut développer C(z) sous la forme suivante: q lI (FEZ) ( 4) C(z)= " f (l+z -q+ Z j=l qui peut être représenté sous la forme d'une série: ( 5) C(Z)=Zcz n=O Les coefficients Cn constituent une suite qui tend vers O lorsque N tend vers l'infini et C(z) est la transmittance d'un filtre idéal anticausal à réponse impulsionnelle infinie stable Une troncature du développement de la transmittance C(z) à un ordre prédéterminé L permet une réalisation pratique avec un filtre numérique transversal de longueur finie La transmittance de ce filtre peut s'exprimer sous la forme polynômiale suivante: ( 6) P(z)=co+c,z++c Lz L La réponse obtenue après passage dans le filtre de transmittance P(z) s'écrit: L ( 7) H(z) P(z)=G(z)+ +en zn pour LK n=L+I-K En effet, H(z)P(z)=H(z)(C(z) cz 9 n=L+ 1 =G(z)-H(z) 1 cn zn n=L+I Les termes en z-n de H(z) sont nuls pour m>K tandis que les termes en zk de H(z)P(z) sont nuls pour k>L On peut donc transformer le deuxième terme du deuxième membre de l'équation précédente comme suit: + co L H(z) cnzn X E Snzn n= L+ 1 n=L+ 1-K Si, par exemple, H(z)=ho( 1-pz-1), alors l'expression ( 7) peut s'écrire de la façon suivante: 1/z L ( 8) H(z) P(z)=G(z)+h ( 1) ZL ( 8) el A La troncature du filtre peut être adoucie en la pondérant par une fenêtre. Lorsque la réponse impulsionnelle H(z) possède des zéros sur le cercle unitaire, ceux-ci ne sont pas affectés par le filtre L'expression de la transmittance théorique C(z) exprimée par l'équation ( 3) précitée reste inchangée mais les équations ( 1) et ( 2) deviennent: P q r H(z)=ho ( 1-aiz-l) H ( 1-jz-l,) H ( 1-Yiz-l) ( 9) l i||a;y= l j= 1 i=l Ioi 11 3 > jj=i K=p+q+ r et la réponse à phase minimale s'exprime de la façon suivante: ( 10) G(z)=hof-I ( 1 z-') ( 1 z-') ( 1-ri Z') j=l i=l j=l i= Le procédé selon l'invention met en oeuvre la méthode du cepstre fondée sur le fait que la causalité de la réponse impulsionnelle du filtre de transmittance G(z) = -ln(G(z)) est obtenue lorsque le filtre de transmittance G(z) est à phase minimale. Pour lzoil O n-0 On considère maintenant la réponse en fréquence H(f) du filtre En supposant que la réponse impulsionnelle n'a pas de zéro sur le cercle unitaire, le module lH(f)l ne peut pas s'annuler et on peut écrire, pour les expressions respectives de la réponse en fréquence H(f) du filtre, de la transmittance en fréquence à phase minimale G(f) et de la transmittance en fréquence C(f) du circuit de filtrage idéal: ( 13) H(f) = ea(f)J(f) ( 14) G( O =e-&( -"- ( 15) C(f)=ejf(f - A A avec gn = g 1,n + jg 2,n et il est possible de séparer a(f) en une partie paire p(f) et une partie impaire q(f) qui peuvent être exprimées de la façon suivante: ( 16) P = a(f) +a(f) =g ,,cos( 2,fn T) a(f)=a(Oa(-f) ( 17)q(f= a(af) =_ g 2 sin( 27 ún T) La partie imaginaire O (f) du cepstre peut alors s'exprimer ainsi: ( 1) 9 (f) = -É (g 1 sin( 2 a T) gi, cos( 2 z T)) ( 18),-o = Tsin( 2 fn'r) l_p(v) cos( 2 =zi 1) du -ô Tcos( 2 zn T) | q(v) sin( 2 a Ui T) dv A partir de la phase O (f) qui ne dépend que des parties paire et impaire de a(f), on peut calculer la transmittance en fréquence C(f) et la réponse à phase minimale G(f). Les développements mathématiques qui viennent d'être présentés ont pour seul but de fournir les éléments nécessaires à la compréhension du fonctionnement de l'algorithme de détermination des coefficients du circuit de filtrage concrétisant le procédé selon l'invention, et l'exploitation pratique de cet algorithme va maintenant être expliquée. Dans le cas o la réponse impulsionnelle H(z) du filtre n'a pas de zéro sur le cercle unitaire, la transformée H(f) est obtenue par transformée de Fourier, éventuellement transformée de Fourier rapide, de la réponse impulsionnelle H(z) Ce premier traitement qui comporte le calcul de N= 2 m échantillons nécessite 2 m N multiplications et 3 m N additions, l'obtention des quantités p(f) et q(f) met en oeuvre 2 N multiplications et 3 N additions, {gn} est une suite décroisante qui tend vers zéro lorsque N tend vers l'infini En limitant le développement de G(z) à l'ordre N/2, les coefficients g 1,n et gz 2,n sont respectivement obtenus à l'aide de transformées en cosinus et en sinus inverses: Pour g N 2 p()cos(r p(O)k Pourn= 0, n:N ^ I 2 2 p() cos ( 2 r-) + p( 0)+ N ^ 2 2 k Pour n El 1,N/2 -1 l gj 2 = 2-( ri 2 p( 1 T)cos( 271- N-)+p(O)+p(-)(-l)n) k=l N 1 22 k nk Pourn El 1,N/2-1 l g 2,:=- 1 ( E 2 q( -)sn( 2 r-7) et 2 N=o k= 1 2 La fixation arbitraire de g 20 à zéro modifie les expressions ( 2) et ( 3) de la réponse à phase minimale G(z) et de la transmittance théorique C(z) qui deviennent: ( 19) G(z)=Ihl I(l-aiz-') f J( 1-,z-') j-1 j-I q h O FI ( 1-P z) ( 20) j= 1 Ihol li j( 1 z) j= 1 pi Ces calculs mettent en oeuvre 4 m N multiplications et 6 m N additions. La phase O (f) est la partie imaginaire de la transformée de Fourier rapide de la séquence complexe gn = gl, n + Jg 2,n; son calcul nécessite 2 m N multiplications et 3 m N additions. L'expression de la transformée de Fourier C(f) de la transmittance théorique du filtre peut se mettre sous la forme suivante: il ( 21) C(f) = H*(f)ea(f)(cos O(f)-jsin O(f)) Son calcul met en oeuvre 6 N multiplications et 2 N additions La transmittance C(z) est obtenue par transformée de Fourier rapide inverse et nécessite 2 m N multiplications et 3 m N additions L'estimation P(z) de la transmittance du circuit de filtrage selon l'invention est obtenue à partir de la transmittance théorique C(z) dont on ne garde que les L+ 1 coefficients significatifs. Au total, l'algorithme d'adaptation met en oeuvre: l Om N+ 8 N multiplications m N+ 4 N additions entre nombres réels, et aucune division. On est bien en présence d'un algorithme non itératif et le nombre d'opérations n'est pas variable et ne dépend pas de la durée K de la réponse impulsionnelle à transformer, ni de l'ordre L de troncature de la transmittance C(z). Cependant, le nombre d'échantillons N doit satisfaire la relation suivante: ( 22) N> L> K Dans le cas o la réponse impulsionnelle H(z) du filtre a des zéros sur le cercle unitaire, il n'est plus possible de mettre en oeuvre la méthode du cepstre, mais il est encore possible d'utiliser l'algorithme qui vient d'être présenté par une modification heuristique consistant à remplacer l'expression de a(f) par l'expression suivante: ( 23) a(f) = -1/2 ln(H(f)H*(f) + E) o 8 doit être choisi avec soin pour éviter soit la divergence de l'algorithme ( 8 trop petit), soit une diminution des performances de l'algorithme (ú trop grand). Les modifications mineures apportée à l'algorithme n'augmentent sa complexité que de N additions réelles et de N comparaisons réelles. Les moyens de calcul algorithmique 1 délivrent une réponse estimée G(z) de la réponse à phase minimale G(z). Cette réponse estimée peut revêtir deux formes: Une première forme d'estimation Go(z) est obtenue en ne gardant que les termes significatifs des coefficients du polynôme G(z) obtenu par transformée de Fourier discrète inverse de la transformée G(f) calculée à partir de a(f) et de 0 (f) Le calcul de cette estimation nécessite 2 m N+ 2 N multiplications et 3 m N additions. Une seconde forme d'estimation G 1 (z) est constituée par la partie causale du produit H(z)P(z) et nécessite 2 (K+ 1)(K+ 2)multiplications et K(K+ 1) additions entre nombres réels. Dans un exemple d'application du procédé selon l'invention à un canal téléphonique, on a mis en évidence le transfert de l'énergie E de la réponse impulsionnelle globale du canal et du circuit de filtrage vers l'origine des temps, en référence à la figure 2 o les courbes 20, 21 représentent respectivement la variation de l'énergie E en fonction du temps échantillonné n T en présence et en l'absence du circuit de filtrage selon l'invention. La mise en oeuvre du procédé selon l'invention dans un dispositif complet de réception tel que représenté sous forme d'un schéma-bloc en figure 1 a été simulée Un tel dispositif de réception peut être utilisé dans des transmissions numériques avec des mobiles à accès multiple à répartition dans le temps (AMRT) Les messages sont émis sous la forme d'une succession d'intervalles de temps, au centre desquels est placée une séquence d'apprentissage permettant d'estimer la réponse du filtre et d'adapter les coefficients du circuit de filtrage en mettant en oeuvre le procédé selon l'invention. A titre d'exemple non limitatif, on peut considérer sur le canal de transmission 2 une modulation d'amplitude en quadrature à quatre états, une mémoire du filtre prise en compte par le récepteur 10 de 6 symboles, un nombre d'états pris en compte dans l'algorithme de Viterbi modifié du détecteur DFSE 4 égal à 16, correspondant à une mémoire de 2 symboles seulement L'intervalle de temps est constitué de 22 symboles d'apprentissage et de 144 symboles d'information. Les performances du récepteur 10 sont évaluées par un taux d'erreur par bit t en fonction du rapport de l'énergie moyenne Eb par bit en sortie du canal 2 sur la densité spectrale unilatérale de puissance No d'un bruit blanc additif gaussien Différents profils de canaux ont été simulés: un premier canal ci est un canal à deux trajets d'égales puissances distants de 6 échantillons de temps, un second canal c 2 est un canal à deux trajets de puissances inégales, distants de 6 échantillons de temps, un troisième canal c 3 est un canal à trois trajets d'égales puissances, situés respectivement aux instants 0, 3 Ts et 6 Ts, o Ts est la période d'échantillonnage. En référence à la figure 3 qui correspond au cas du premier canal ci, les courbes 31 et 30 représentent respectivement l'évolution du taux d'erreur binaire t en fonction du rapport Eb/No dans un dispositif récepteur non muni d'un circuit de filtrage selon l'invention (*) et dans un dispositif récepteur équipé d'un circuit de filtrage selon l'invention (o) On observe que le taux d'erreur est systématiquement diminué en présence d'un circuit de filtrage et que l'amélioration croît significativement avec l'énergie moyenne. Dans le cas du canal c 2 à deux trajets de puissances inégales illustré en figure 4, l'amélioration du taux d'erreur binaire est encore très significative, l'écart entre les deux courbes 40, 41 étant sensiblement inférieur à celui observé dans le canal ci. Dans le cas du canal c 3 à trois trajets d'égales puissances illustré en figure 5, l'écart entre les courbes , 51 de taux d'erreur binaire correspondant respectivement à la présence et à l'absence d'un circuit de filtrage selon l'invention est particulièrement important, en particulier pour des niveaux élevés d'énergie moyenne. Bien sûr, l'invention n'est pas limitée aux exemples qui viennent d'être décrits et de nombreux aménagements peuvent être apportés à ces exemples sans sortir du cadre de l'invention. Ainsi, le procédé selon l'invention peut s'appliquer à un nombre quelconque de réponses impulsionnelles Le nombre de zéros de la réponse impulsionnelle d'un canal ne constitue pas non plus une entrave à l'application de ce procédé, qui ne peut être éventuellement limitée que par les capacités de calcul rapide mises en oeuvre pour exécuter l'algorithme d'adaptation. REVENDICATIONS 1 Procédé pour déterminer la transmittance P(z) d'un circuit de filtrage ( 3) associé à un filtre ( 2) présentant une réponse impulsionnelle H(z) caractérisé en ce que, ledit circuit de filtrage ( 3) étant prévu pour que la réponse impulsionnelle globale G(z) dudit filtre ( 2) et de ce circuit ( 3) soit à phase minimale, il comprend: une détermination de la réponse en fréquence H(f) dudit filtre ( 2) par transformée de Fourier de la réponse impulsionnelle H(z) de ce filtre ( 2), une détermination de la partie réelle a(f) du cepstre H(f) = -ln(H(f)) de la réponse impulsionnelle H(f) dudit filtre ( 2) pour extraire de cette partie réelle a(f) une partie paire p(f) et une partie impaire q(f), une détermination du cepstre G(f) = -ln(G(f)) de la réponse globale en fréquence à partir desdites parties paire p(f) et impaire q(f) du cepstre H(f) de la réponse impulsionnelle dudit filtre ( 2), une détermination de la réponse globale en fréquence G(f) à partir du cepstre G(f) de cette réponse, une détermination de la réponse théorique en fréquence C(f) dudit circuit de filtrage ( 3) au moyen de la réponse globale en fréquence G(f) et de la réponse en fréquence H(f) dudit filtre ( 2), un calcul de la transmittance théorique C(z) dudit circuit de filtrage ( 3) par transformée de Fourier inverse de la réponse en fréquence C(f) de ce circuit ( 3) et, une détermination de la transmittance P(z) dudit circuit de filtrage ( 3) par troncature de la transmittance C(z) dont ne sont conservés qu'un nombre prédéterminé de coefficients. 2 Procédé selon la revendication 1, caractérisé en ce qu'il permet une estimation de la réponse impulsionnelle globale G(z) à partir d'un développement limité du cepstre G(z). 3 Procédé selon l'une des revendications 1 ou 2, caractérisé en ce qu'il comprend en outre l'obtention d'une première estimation Go(z) de la réponse à phase minimale G(z) comportant: un calcul des coefficients de la réponse impulsionnelle à phase minimale G(z) par transformée de Fourier discrète inverse de la réponse en fréquence à phase minimale G(f), et la conservation d'un nombre prédéterminé de coefficients significatifs pour obtenir ladite estimation GO(z). 4 Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre l'obtention d'une deuxième estimation G 1 (z) sous la forme de la partie causale du produit de la réponse impulsionnelle du filtre H(z) et de l'estimation P(z) de la transmittance P(z) du filtre ( 3). Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend un nombre d'opérations à effectuer et un temps de calcul fixes et indépendants de la longueur du filtre H(z), et en ce que lesdits développements limités et lesdites transformées ont des tailles prédéterminées. 6 Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend en outre un traitement heuristique du cas o la réponse impulsionnelle H(z) du filtre ( 2) a des zéros sur le cercle unité, par la modification suivante: _a(f) = - 1/2 ln(H(f)H*(f) + E), e représentant un paramètre de réglage. 7 Circuit de filtrage numérique transversal correcteur de phase ( 3) associé à un filtre ( 2) et prévu pour que la réponse impulsionnelle globale dudit filtre ( 2) et dudit circuit de filtrage ( 3) soit à phase minimale, comprenant des moyens de calcul algorithmique ( 1) associés audit circuit de filtrage correcteur ( 3) pour en adapter les coefficients par mise en oeuvre du procédé selon l'une quelconque des revendications précédentes.