Esistono alternative alla trasformazione bilineare?


26

Quando si progetta un filtro digitale basato su un filtro analogico di solito si utilizza la trasformazione bilineare . Per approssimare una funzione di trasferimento discreta dalla funzione di trasferimento (continua) analogica A ( s ) che sostituiamoDa(z)A(s)

z=1+sT/21sT/2

dove è il periodo di campionamento. In alternativa, per approssimare una funzione di trasferimento continuo A a ( s ) dalla funzione di trasferimento discreta D ( z ) sostituiamoTAa(s)D(z)

s=2Tz1z+1

Esistono metodi alternativi per eseguire tali conversioni? Ci sono approssimazioni migliori?

Risposte:


16

I filtri analogici sono stabili se i poli si trovano nella metà sinistra del piano s (figura a sinistra) e i filtri digitali sono stabili se i poli si trovano all'interno del cerchio dell'unità (figura a destra). Quindi matematicamente tutto ciò che serve per convertire da analogico a digitale è una mappatura (conforme?) Dal al disco dell'unità e dall'asse ȷ Ω al cerchio dell'unità | z | = 1 . Qualsiasi trasformazione che faccia questo è un possibile candidato per essere un'alternativa alla trasformazione bilaterale.ȷΩ|z|=1

inserisci qui la descrizione dell'immagine

L1Z

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

a(t)Ta[n]

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

Tuttavia, ci sono differenze chiave tra i due.

Metodo di invarianza dell'impulso:

In questo metodo, si espande la funzione di trasferimento analogico come frazioni parziali (non nella trasformata Z corrispondente come menzionato da Peter ) come

A(s)=mCmsαm

Cmαm

Il motivo per cui fallisce è anche abbastanza chiaro. Se avessi un polinomio nel numeratore dello stesso grado del denominatore, avrai un termine permanente indipendente, che alla trasformazione inversa darà una funzione delta che non può essere campionata.

αmeαmT

Trasformata Z abbinata

βmeβmTαmeαmT

A(s)=m(sβm)n(sαn)m(1z1eβmT)n(1z1eαnT)

Puoi facilmente vedere la limitazione di entrambi questi metodi. L'impulso invariante è applicabile solo se il filtro è passa basso e il metodo di trasformazione z adattato è applicabile ai filtri passa banda e passa banda (e passa alto fino alla frequenza di Nyquist). In pratica sono anche limitati dalla frequenza di campionamento (dopotutto, puoi salire fino a un certo punto) e soffrono degli effetti dell'aliasing.

La trasformazione bilineare è di gran lunga il metodo più comunemente usato nella pratica e le due sopra menzionate sono piuttosto più per gli interessi accademici. Per quanto riguarda la conversione in analogico, mi dispiace ma non lo so e non posso esservi di grande aiuto in quanto non uso quasi mai i filtri analogici.


Wow Wow ..... questa è la migliore spiegazione che ho visto su questo argomento. Grazie mille per la condivisione. Bel lavoro.

la trasformazione z abbinata è migliore per i filtri Bessel perché la caratteristica importante dei filtri Bessel è il loro ritardo di gruppo piatto, non la loro risposta in frequenza
endolith

9

sz

Alcuni esempi sono:

Trasformazione Z abbinata

s

Y(s)=a0s+s0+a1s+s1+...

E la conversione di ogni parte dell'espansione della frazione parziale viene eseguita direttamente utilizzando:

s+sn=1z1exp(snT)

La regola di Simpson

Un'interpretazione della trasformata bilineare è che è un modo di trasformare da tempo continuo a tempo discreto mediante integrazione approssimativa usando la Regola trapezoidale .

Una tecnica più accurata per l'integrazione approssimativa utilizza la regola di Simpson. Se viene utilizzata questa approssimazione, la mappatura risultante è:

s=3Tz21z2+4z+1

1
Regola di Simpson, interpolazione essenzialmente quadratica (dove la regola trapezoidale è lineare)?
Peter Mortensen,

1
@Peter Mortensen: Sì, praticamente!
Peter K.

La tua trasformazione Z abbinata è diversa da quella di Lorem Ipsum? Non vedo decomposizione parziale della frazione da nessun'altra parte.
endolito

@endolith vedi il link di Wikipedia nella mia risposta. Ecco da dove l'ho preso. 😂 Ho risposto prima a Lorem e non l'ho modificato.
Peter K.
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.