Existe-t-il des alternatives à la transformation bilinéaire?

26

Lors de la conception d'un filtre numérique basé sur un filtre analogique, nous utilisons généralement la transformation bilinéaire . Pour approximer une fonction de transfert discrète de la fonction de transfert analogique (continue) nous substituonsA ( s )une(z)UNE(s)

z=1+sT/21-sT/2

où est la période d'échantillonnage. Alternativement, pour approximer une fonction de transfert continue de la fonction de transfert discrète nous substituonsA a ( s ) D ( z )TUNEune(s)(z)

s=2Tz-1z+1

Existe-t-il d'autres méthodes pour effectuer de telles conversions? Y a-t-il de meilleures approximations?

Phonon
la source

Réponses:

16

Les filtres analogiques sont stables si les pôles sont dans la moitié gauche du plan s (figure à gauche) et les filtres numériques sont stables si les pôles sont à l'intérieur du cercle unitaire (figure à droite). Donc, mathématiquement, tout ce qui est nécessaire pour convertir l'analogique en numérique est un mappage (conforme?) Du demi-espace au disque unitaire et l' axe au cercle unitaire . Toute transformation qui le fait est un candidat possible pour être une alternative à la transformation bilatérale.| z | = 1ȷΩ|z|=1

entrez la description de l'image ici

Deux des méthodes bien connues sont la méthode d'invariance d'impulsion et la méthode de transformée en Z appariée . Conceptuellement, les deux sont similaires à l'échantillonnage d'une forme d'onde continue que nous connaissons. Indiquant la transformée de Laplace inverse par et la transformée Z en , ces deux méthodes impliquent de calculer la réponse impulsionnelle du filtre analogique comme ZL1Z

a(t)=L1{A(s)}

et échantillonnage à un intervalle d'échantillonnage suffisamment élevé pour éviter le repliement. La fonction de transfert du filtre numérique est alors obtenue à partir de la séquence échantillonnée commea(t)Ta[n]

Da(z)=Z{a[n]}

Cependant, il existe des différences clés entre les deux.

Méthode d'invariance d'impulsion:

Dans cette méthode, vous développez la fonction de transfert analogique sous forme de fractions partielles (pas dans la transformation Z correspondante comme mentionné par Peter ) comme

A(s)=mCms-αm

où est une constante et sont les pôles. Mathématiquement, toute fonction de transfert avec un numérateur de degré inférieur au dénominateur peut être exprimée comme une somme de fractions partielles . Seuls les filtres passe-bas satisfont à ce critère (passe-haut et passe-bande / coupe-bande ont au moins le même degré), et donc la méthode invariante d'impulsion ne peut pas être utilisée pour concevoir d'autres filtres.Cmαm

La raison de l'échec est également assez claire. Si vous aviez un polynôme au numérateur du même degré qu'au dénominateur, vous aurez un terme constant libre qui, lors de la transformation inverse, donnera une fonction delta qui ne peut pas être échantillonnée.

Si vous effectuez les transformations inverses de Laplace et Z en avant, vous verrez que les pôles sont transformés en ce qui signifie que si votre filtre analogique est stable, le filtre numérique sera également stable . Il préserve donc la stabilité du filtre.αmeαmT

Transformation Z adaptée

Dans cette méthode, au lieu de diviser la réponse impulsionnelle en fractions partielles, vous effectuez une simple transformation des pôles et des zéros de la même manière (appariée) que et (également préservant la stabilité), donnantβmeβmTαmeαmT

UNE(s)=m(s-βm)n(s-αn)m(1-z-1eβmT)n(1-z-1eαnT)

Vous pouvez facilement voir la limitation de ces deux méthodes. L'invariant d'impulsion n'est applicable que si votre filtre est passe-bas et que la méthode de transformée z adaptée est applicable aux filtres coupe-bande et passe-bande (et passe-haut jusqu'à la fréquence de Nyquist). Ils sont également limités dans la pratique par le taux d'échantillonnage (après tout, vous ne pouvez aller jusqu'à un certain point) et souffrent des effets de repliement.

La transformation bilinéaire est de loin la méthode la plus couramment utilisée dans la pratique et les deux ci-dessus sont plutôt plus pour les intérêts académiques. Quant à la conversion en analogique, je suis désolé mais je ne sais pas et je ne peux pas être d'une grande aide car je n'utilise presque jamais de filtres analogiques.

Lorem Ipsum
la source
Wow Wow ..... ce sont les meilleures explications que j'ai vues sur ce sujet. Merci beaucoup d'avoir partagé. Beau travail.
la transformée z adaptée est meilleure pour les filtres de Bessel parce que la caractéristique importante des filtres de Bessel est leur retard de groupe plat, pas leur réponse en fréquence
endolith
9

Il existe de nombreuses façons de faire le mappage de à . La communauté du contrôle a quelque chose à dire à ce sujet.sz

Certains exemples sont:

La transformée en Z correspondante

s

Oui(s)=une0s+s0+une1s+s1+...

Et la conversion de chaque partie de l'expansion de fraction partielle se fait directement en utilisant:

s+sn=1-z-1exp(-snT)

La règle de Simpson

Une interprétation de la transformée bilinéaire est qu'elle est un moyen de passer du temps continu au temps discret par intégration approximative en utilisant la règle trapézoïdale .

Une technique plus précise pour une intégration approximative utilise la règle de Simpson. Si cette approximation est utilisée, le mappage résultant est:

s=3Tz2-1z2+4z+1
Peter K.
la source
1
Règle de Simpson, interpolation essentiellement quadratique (où la règle trapézoïdale est linéaire)?
Peter Mortensen
1
@Peter Mortensen: Oui, à peu près!
Peter K.
Votre transformation Z assortie est-elle différente de celle de Lorem Ipsum? Je ne vois aucune décomposition partielle de fraction ailleurs.
endolith
@endolith voir le lien wikipedia dans ma réponse. C'est de là que je l'ai obtenu. 😂 J'ai répondu avant Lorem et je ne l'ai pas édité.
Peter K.