¿Hay alternativas a la transformación bilineal?

26

Cuando diseñamos un filtro digital basado en un filtro analógico, usualmente usamos la transformación bilineal . Para aproximar una función de transferencia discreta de la función de transferencia analógica (continua) , sustituimosDa(z)A(s)

z=1+sT/21sT/2

donde es el período de muestreo. Alternativamente, para aproximar una función de transferencia continua de la función de transferencia discreta , sustituimosTAa(s)D(z)

s=2Tz1z+1

¿Existen métodos alternativos para realizar tales conversiones? ¿Hay mejores aproximaciones?

Phonon
fuente

Respuestas:

16

Los filtros analógicos son estables si los polos están en la mitad izquierda del plano s (figura a la izquierda) y los filtros digitales son estables si los polos están dentro del círculo unitario (figura a la derecha). Entonces, matemáticamente, todo lo que se necesita para convertir de analógico a digital es un mapeo (¿conforme?) Desde el medio espacio al disco de la unidad y el al círculo de la unidad . Cualquier transformación que haga esto es un posible candidato para ser una alternativa a la transformación bilateral.ȷΩ|z|=1

ingrese la descripción de la imagen aquí

Dos de los métodos bien conocidos son el método de invariancia de impulso y el método de transformación Z emparejado . Conceptualmente, ambos son similares al muestreo de una forma de onda continua con la que estamos familiarizados. Denotando la transformada inversa de Laplace por y la transformada Z como , ambos métodos implican calcular la respuesta al impulso del filtro analógico comoL1Z

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

y muestrear en un intervalo de muestreo que es lo suficientemente alto como para evitar el aliasing. La función de transferencia del filtro digital se obtiene de la secuencia muestreada comoa(t)Ta[n]

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

Sin embargo, hay diferencias clave entre los dos.

Método de invariancia de impulso:

En este método, expande la función de transferencia analógica como fracciones parciales (no en la transformación Z coincidente como lo menciona Peter ) como

A(s)=mCmsαm

donde es una constante y son los polos. Matemáticamente, cualquier función de transferencia con un numerador de menor grado que el denominador se puede expresar como una suma de fracciones parciales . Solo los filtros de paso bajo satisfacen este criterio (el paso alto y el paso de banda / parada de banda tienen al menos el mismo grado) y, por lo tanto, el método de impulso invariante no se puede utilizar para diseñar otros filtros.Cmαm

La razón por la que falla también es bastante clara. Si tuviera un polinomio en el numerador del mismo grado que en el denominador, tendrá un término constante independiente, que al transformarse inversamente, dará una función delta que no se puede muestrear.

Si lleva a cabo la transformación inversa de Laplace y Z hacia adelante, verá que los polos se transforman como que significa que si su filtro analógico es estable, el filtro digital también será estable . Por lo tanto, preserva la estabilidad del filtro.αmeαmT

Transformada Z combinada

En este método, en lugar de dividir la respuesta al impulso como fracciones parciales, realiza una transformación simple de los polos y los ceros de manera similar (coincidente) como y (también preservando la estabilidad), dandoβmeβmTαmeαmT

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

Puede ver fácilmente la limitación de ambos métodos. Impulse invariante es aplicable solo si su filtro es de paso bajo y el método de transformación z combinado es aplicable a los filtros de paso de banda y paso de banda (y paso alto hasta la frecuencia Nyquist). También están limitados en la práctica por la frecuencia de muestreo (después de todo, solo puede llegar a un cierto punto) y sufren los efectos del aliasing.

La transformación bilineal es, con mucho, el método más utilizado en la práctica y los dos anteriores son más bien para intereses académicos. En cuanto a la conversión a analógico, lo siento, pero no lo sé y no puedo ser de mucha ayuda, ya que casi nunca uso filtros analógicos.

Lorem Ipsum
fuente
Wow Wow ..... esta es la mejor explicación que he visto sobre este tema. Muchas gracias por compartir. Hermoso trabajo.
la transformación z coincidente es mejor para los filtros Bessel porque la característica importante de los filtros Bessel es su retraso de grupo plano, no su respuesta de frecuencia
endolito el
9

Hay muchas formas de hacer el mapeo de a . La comunidad de control tiene algunas cosas que decir al respecto.sz

Algunos ejemplos son:

La transformación Z combinada

Aquí, la función de transferencia de dominio se escribe como una expansión de fracción parcial:s

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

Y la conversión de cada parte de la expansión de fracción parcial se realiza directamente usando:

s+sn=1z1exp(snT)

La regla de Simpson

Una interpretación de la transformación bilineal es que es una forma de transformar de tiempo continuo a tiempo discreto mediante integración aproximada usando la regla trapezoidal .

Una técnica más precisa para la integración aproximada utiliza la regla de Simpson. Si se utiliza esta aproximación, la asignación resultante es:

s=3Tz21z2+4z+1
Peter K.
fuente
1
¿La regla de Simpson, esencialmente interpolación cuadrática (donde la regla trapezoidal es lineal)?
Peter Mortensen
1
@ Peter Mortensen: ¡Sí, más o menos!
Peter K.
¿Es su transformación Z combinada diferente de la de Lorem Ipsum? No veo descomposición de fracción parcial en ningún otro lugar.
endolito el
@endolith ver el enlace de wikipedia en mi respuesta. De ahí es de donde lo saqué. 😂 Respondí antes de Lorem, y no lo he editado.
Peter K.