stringtranslate.com

Transformada Z

En matemáticas y procesamiento de señales , la transformada Z convierte una señal de tiempo discreto , que es una secuencia de números reales o complejos, en una representación de dominio de frecuencia de valor complejo (el dominio z o plano z ). [1] [2]

Puede considerarse un equivalente en tiempo discreto de la transformada de Laplace (el dominio s o plano s ). [3] Esta similitud se explora en la teoría del cálculo de escala de tiempo .

Mientras que la transformada de Fourier de tiempo continuo se evalúa en el eje vertical del dominio s (el eje imaginario), la transformada de Fourier de tiempo discreto se evalúa a lo largo del círculo unitario del dominio z. El semiplano izquierdo del dominio s se asigna al área dentro del círculo unitario del dominio z, mientras que el semiplano derecho del dominio s se asigna al área fuera del círculo unitario del dominio z.

En el procesamiento de señales, uno de los medios para diseñar filtros digitales es tomar diseños analógicos, someterlos a una transformación bilineal que los mapea del dominio s al dominio z y luego producir el filtro digital mediante inspección, manipulación o aproximación numérica. Estos métodos tienden a no ser precisos excepto en la proximidad de la unidad compleja, es decir, a frecuencias bajas.

Historia

El concepto fundacional que ahora se reconoce como la transformada Z, que es una piedra angular en el análisis y diseño de sistemas de control digital, no era completamente novedoso cuando surgió a mediados del siglo XX. Sus principios embrionarios se remontan al trabajo del matemático francés Pierre-Simon Laplace , quien es más conocido por la transformada de Laplace , una técnica matemática estrechamente relacionada. Sin embargo, la formulación explícita y la aplicación de lo que ahora entendemos como la transformada Z fueron significativamente avanzadas en 1947 por Witold Hurewicz y colegas. Su trabajo fue motivado por los desafíos presentados por los sistemas de control de datos muestreados, que se estaban volviendo cada vez más relevantes en el contexto de la tecnología de radar durante ese período. La transformada Z proporcionó un método sistemático y efectivo para resolver ecuaciones diferenciales lineales con coeficientes constantes, que son omnipresentes en el análisis de señales y sistemas de tiempo discreto. [4] [5]

El método se perfeccionó aún más y obtuvo su nomenclatura oficial, "la transformada Z", en 1952, gracias a los esfuerzos de John R. Ragazzini y Lotfi A. Zadeh , quienes formaban parte del grupo de control de datos muestreados en la Universidad de Columbia. Su trabajo no solo solidificó el marco matemático de la transformada Z, sino que también amplió su alcance de aplicación, particularmente en el campo de la ingeniería eléctrica y los sistemas de control. [6] [7]

El desarrollo de la transformada Z no se detuvo con Ragazzini y Zadeh. Una extensión notable, conocida como la transformada Z modificada o avanzada, fue introducida más tarde por Eliahu I. Jury . El trabajo de Jury amplió la aplicabilidad y robustez de la transformada Z, especialmente en el manejo de condiciones iniciales y proporcionando un marco más completo para el análisis de sistemas de control digital. Esta formulación avanzada ha desempeñado un papel fundamental en el diseño y análisis de estabilidad de sistemas de control de tiempo discreto, contribuyendo significativamente al campo del procesamiento de señales digitales. [8] [9]

Curiosamente, los fundamentos conceptuales de la transformada Z se cruzan con un concepto matemático más amplio conocido como el método de funciones generadoras , una poderosa herramienta en combinatoria y teoría de la probabilidad. Esta conexión fue insinuada ya en 1730 por Abraham de Moivre , una figura pionera en el desarrollo de la teoría de la probabilidad. De Moivre utilizó funciones generadoras para resolver problemas de probabilidad, sentando las bases para lo que eventualmente evolucionaría en la transformada Z. Desde una perspectiva matemática, la transformada Z puede verse como una instancia específica de una serie de Laurent , donde la secuencia de números bajo investigación se interpreta como los coeficientes en la expansión (de Laurent) de una función analítica . Esta perspectiva no solo resalta las profundas raíces matemáticas de la transformada Z, sino que también ilustra su versatilidad y amplia aplicabilidad en diferentes ramas de las matemáticas y la ingeniería. [10]

Definición

La transformada Z se puede definir como una transformada unilateral o bilateral . (Al igual que tenemos la transformada de Laplace unilateral y la transformada de Laplace bilateral ). [11]

Transformada Z bilateral

La transformada Z bilateral o de dos lados de una señal de tiempo discreto es la serie de potencia formal definida como:

donde es un entero y es, en general, un número complejo . En forma polar , puede escribirse como:

donde es la magnitud de , es la unidad imaginaria , y es el argumento complejo (también conocido como ángulo o fase ) en radianes .

Transformada Z unilateral

Alternativamente, en los casos donde se define solo para , la transformada Z unilateral o de un solo lado se define como:

En el procesamiento de señales , esta definición se puede utilizar para evaluar la transformada Z de la respuesta al impulso unitario de un sistema causal de tiempo discreto .

Un ejemplo importante de la transformada Z unilateral es la función generadora de probabilidad , donde el componente es la probabilidad de que una variable aleatoria discreta tome el valor . Las propiedades de las transformadas Z (enumeradas en el § Propiedades) tienen interpretaciones útiles en el contexto de la teoría de la probabilidad.

Transformada Z inversa

La transformada Z inversa es:

donde es una trayectoria cerrada en sentido antihorario que rodea el origen y se encuentra completamente en la región de convergencia (ROC). En el caso en que la ROC sea causal (ver Ejemplo 2), esto significa que la trayectoria debe rodear todos los polos de .

Un caso especial de esta integral de contorno se produce cuando es el círculo unitario. Este contorno se puede utilizar cuando la ROC incluye el círculo unitario, lo que siempre se garantiza cuando es estable, es decir, cuando todos los polos están dentro del círculo unitario. Con este contorno, la transformada Z inversa se simplifica a la transformada de Fourier inversa de tiempo discreto , o serie de Fourier , de los valores periódicos de la transformada Z alrededor del círculo unitario:

La transformada Z con un rango finito de y un número finito de valores espaciados uniformemente se puede calcular de manera eficiente a través del algoritmo FFT de Bluestein . La transformada de Fourier de tiempo discreto (DTFT), que no debe confundirse con la transformada de Fourier discreta (DFT), es un caso especial de una transformada Z de este tipo que se obtiene al restringir que se encuentre en el círculo unitario.

Región de convergencia

La región de convergencia (ROC) es el conjunto de puntos en el plano complejo para el cual la suma de la transformada Z converge (es decir, no aumenta en magnitud hasta el infinito):

Ejemplo 1 (sin ROC)

Ampliando el intervalo se obtiene

Mirando la suma

Por lo tanto, no hay valores de que satisfagan esta condición.

Ejemplo 2 (ROC causal)

ROC (azul), | z | = .5 (círculo negro punteado) y el círculo unitario (círculo gris punteado).

Sea (donde es la función escalonada de Heaviside ). Desarrollando el intervalo se obtiene

Mirando la suma

La última igualdad surge de la serie geométrica infinita y la igualdad solo se cumple si que puede reescribirse en términos de como Por lo tanto, la ROC es En este caso, la ROC es el plano complejo con un disco de radio 0,5 en el origen "perforado".

Ejemplo 3 (ROC anticausal)

ROC (azul), | z | = .5 (círculo negro punteado) y el círculo unitario (círculo gris punteado).

Sea (donde es la función escalonada de Heaviside ). Desarrollando el intervalo se obtiene

Mirando la suma

y usando nuevamente la serie geométrica infinita , la igualdad solo se cumple si lo cual puede reescribirse en términos de como Por lo tanto, la ROC es En este caso, la ROC es un disco centrado en el origen y de radio 0,5.

Lo que diferencia este ejemplo del anterior es únicamente el ROC. Esto tiene la intención de demostrar que el resultado de la transformación por sí solo no es suficiente.

Ejemplos de conclusión

Los ejemplos 2 y 3 muestran claramente que la transformada Z de es única cuando y solo cuando se especifica la ROC. La creación del gráfico de polos y ceros para el caso causal y anticausal muestra que la ROC para ambos casos no incluye el polo que está en 0,5. Esto se extiende a los casos con múltiples polos: la ROC nunca contendrá polos.

En el ejemplo 2, el sistema causal produce un ROC que incluye, mientras que el sistema anticausal del ejemplo 3 produce un ROC que incluye

ROC se muestra como un anillo azul 0,5 < | z | < 0,75

En sistemas con múltiples polos es posible tener un ROC que no incluya ni El ROC crea una banda circular. Por ejemplo,

Tiene polos en 0,5 y 0,75. La ROC será 0,5 < | z | < 0,75, que no incluye ni el origen ni el infinito. Un sistema de este tipo se denomina sistema de causalidad mixta, ya que contiene un término causal y un término anticausal.

La estabilidad de un sistema también se puede determinar conociendo únicamente la ROC. Si la ROC contiene el círculo unitario (es decir, | z | = 1), entonces el sistema es estable. En los sistemas anteriores, el sistema causal (ejemplo 2) es estable porque | z | > 0,5 contiene el círculo unitario.

Supongamos que se nos proporciona una transformada Z de un sistema sin una ROC (es decir, una ambigua ). Podemos determinar una única siempre que deseemos lo siguiente:

Para que haya estabilidad, la ROC debe contener el círculo unitario. Si necesitamos un sistema causal, entonces la ROC debe contener el infinito y la función del sistema será una secuencia del lado derecho. Si necesitamos un sistema anticausal, entonces la ROC debe contener el origen y la función del sistema será una secuencia del lado izquierdo. Si necesitamos estabilidad y causalidad, todos los polos de la función del sistema deben estar dentro del círculo unitario.

Lo único entonces se puede encontrar.

Propiedades

Teorema de Parseval

Teorema del valor inicial : Sies causal, entonces

Teorema del valor final : Si los polos deestán dentro del círculo unitario, entonces

Tabla de pares de transformadas Z comunes

Aquí:

es la función escalonada unitaria (o de Heaviside) y

es la función de impulso unitario de tiempo discreto (cf. función delta de Dirac , que es una versión de tiempo continuo). Las dos funciones se eligen juntas de modo que la función de escalón unitario sea la acumulación (total acumulado) de la función de impulso unitario.

Relación con la serie de Fourier y la transformada de Fourier

Para los valores de en la región , conocida como círculo unitario , podemos expresar la transformada como una función de una única variable real definiendo Y la transformada bilateral se reduce a una serie de Fourier :

que también se conoce como la transformada de Fourier de tiempo discreto (DTFT) de la secuencia. Esta función periódica es la suma periódica de una transformada de Fourier , lo que la convierte en una herramienta de análisis ampliamente utilizada. Para entender esto, sea la transformada de Fourier de cualquier función, , cuyas muestras en algún intervalo sean iguales a la secuencia. Entonces la DTFT de la secuencia se puede escribir de la siguiente manera.

donde tiene unidades de segundos, tiene unidades de hercios . La comparación de las dos series revela que es una frecuencia normalizada con unidad de radianes por muestra . El valor corresponde a . Y ahora, con la sustitución, la ecuación 1 se puede expresar en términos de (una transformada de Fourier):

A medida que cambia el parámetro T , los términos individuales de la ecuación 2 se alejan o se acercan a lo largo del eje f . Sin embargo, en la ecuación 3 , los centros permanecen separados 2 π , mientras que sus anchos se expanden o contraen. Cuando la secuencia representa la respuesta al impulso de un sistema LTI , estas funciones también se conocen como su respuesta en frecuencia . Cuando la secuencia es periódica, su DTFT es divergente en una o más frecuencias armónicas y cero en todas las demás frecuencias. Esto a menudo se representa mediante el uso de funciones delta de Dirac variantes en amplitud en las frecuencias armónicas. Debido a la periodicidad, solo hay un número finito de amplitudes únicas, que se calculan fácilmente mediante la transformada de Fourier discreta (DFT), mucho más simple. (Véase Transformada de Fourier de tiempo discreto § Datos periódicos ).

Relación con la transformada de Laplace

Transformación bilineal

La transformada bilineal se puede utilizar para convertir filtros de tiempo continuo (representados en el dominio de Laplace) en filtros de tiempo discreto (representados en el dominio Z), y viceversa. Se utiliza la siguiente sustitución:

convertir alguna función en el dominio de Laplace en una función en el dominio Z ( transformación de Tustin ), o

del dominio Z al dominio de Laplace. A través de la transformación bilineal, el plano s complejo (de la transformada de Laplace) se mapea al plano z complejo (de la transformada z). Si bien esta asignación es (necesariamente) no lineal, es útil porque mapea todo el eje del plano s sobre el círculo unitario en el plano z. Como tal, la transformada de Fourier (que es la transformada de Laplace evaluada en el eje) se convierte en la transformada de Fourier de tiempo discreto. Esto supone que existe la transformada de Fourier; es decir, que el eje está en la región de convergencia de la transformada de Laplace.

Transformación con estrella

Dada una transformada Z unilateral de una función muestreada en el tiempo, la transformada con asterisco correspondiente produce una transformada de Laplace y restaura la dependencia de (el parámetro de muestreo):

La transformada de Laplace inversa es una abstracción matemática conocida como función de muestreo de impulso .

Ecuación diferencial lineal con coeficientes constantes

La ecuación de diferencia de coeficientes constantes lineales (LCCD) es una representación de un sistema lineal basada en la ecuación de promedio móvil autorregresivo :

Ambos lados de la ecuación anterior se pueden dividir por si no es cero. Al normalizar con la ecuación LCCD se puede escribir

Esta forma de la ecuación LCCD es favorable para hacer más explícito que la salida "actual" es una función de las salidas pasadas, la entrada actual y las entradas previas.

Función de transferencia

Tomando la transformada Z de la ecuación anterior (utilizando leyes de linealidad y desplazamiento temporal) se obtiene:

donde y son la transformada z de y respectivamente. (Las convenciones de notación normalmente utilizan letras mayúsculas para referirse a la transformada z de una señal indicada por una letra minúscula correspondiente, similar a la convención utilizada para notar las transformadas de Laplace).

La reorganización da como resultado la función de transferencia del sistema :

Ceros y polos

Del teorema fundamental del álgebra, el numerador tiene raíces (que corresponden a los ceros de ) y el denominador tiene raíces (que corresponden a los polos). Reescribiendo la función de transferencia en términos de ceros y polos

donde es el cero y es el polo. Los ceros y los polos son comúnmente complejos y cuando se grafican en el plano complejo (plano z) se denomina diagrama de polos y ceros .

Además, también pueden existir ceros y polos en y Si tomamos en consideración estos polos y ceros así como los ceros y polos de orden múltiple, el número de ceros y polos es siempre igual.

Al factorizar el denominador, se puede utilizar la descomposición en fracciones parciales , que luego se puede transformar nuevamente al dominio del tiempo. De esta manera, se obtendría la respuesta al impulso y la ecuación diferencial de coeficientes constantes lineales del sistema.

Respuesta de salida

Si un sistema de este tipo es impulsado por una señal , entonces la salida es Al realizar una descomposición en fracciones parciales en y luego tomar la transformada Z inversa, se puede encontrar la salida. En la práctica, a menudo es útil descomponer fraccionariamente antes de multiplicar esa cantidad por para generar una forma de que tiene términos con transformadas Z inversas fácilmente computables.

Véase también

Referencias

  1. ^ Mandal, Jyotsna Kumar (2020). " Codificación reversible basada en transformada Z ". Esteganografía reversible y autenticación mediante codificación por transformada . Estudios en inteligencia computacional. Vol. 901. Singapur: Springer Singapur. págs. 157–195. doi :10.1007/978-981-15-4397-5_7. ISBN 978-981-15-4396-8. ISSN  1860-949X. S2CID  226413693. Z es una variable compleja. La transformada Z convierte la señal del dominio espacial discreto en una representación compleja del dominio de frecuencia. La transformada Z se deriva de la transformada de Laplace.
  2. ^ Lynn, Paul A. (1986). "La transformada de Laplace y la transformada z ". Electronic Signals and Systems . Londres: Macmillan Education UK. págs. 225–272. doi :10.1007/978-1-349-18461-3_6. ISBN 978-0-333-39164-8La transformada de Laplace y la transformada z están estrechamente relacionadas con la transformada de Fourier. La transformada z es especialmente adecuada para trabajar con señales y sistemas discretos. Ofrece una notación más compacta y conveniente que la transformada de Fourier de tiempo discreto.
  3. ^ Palani, S. (26 de agosto de 2021). "El análisis de la transformada z de señales y sistemas de tiempo discreto". Señales y sistemas . Cham: Springer International Publishing. págs. 921–1055. doi :10.1007/978-3-030-75742-7_9. ISBN 978-3-030-75741-0. S2CID  238692483. La transformada z es la contraparte discreta de la transformada de Laplace. La transformada z convierte ecuaciones diferenciales de sistemas de tiempo discreto en ecuaciones algebraicas, lo que simplifica el análisis de sistemas de tiempo discreto. La transformada de Laplace y la transformada z son comunes, excepto que la transformada de Laplace se ocupa de señales y sistemas de tiempo continuo.
  4. ^ ER Kanasewich (1981). Análisis de secuencia temporal en geofísica. Universidad de Alberta. págs.186, 249. ISBN 978-0-88864-074-1.
  5. ^ ER Kanasewich (1981). Análisis de secuencia temporal en geofísica (3ª ed.). Universidad de Alberta. págs. 185-186. ISBN 978-0-88864-074-1.
  6. ^ Ragazzini, JR; Zadeh, LA (1952). "El análisis de sistemas de datos muestreados". Transacciones del Instituto Americano de Ingenieros Eléctricos, Parte II: Aplicaciones e Industria . 71 (5): 225–234. doi :10.1109/TAI.1952.6371274. S2CID  51674188.
  7. ^ Cornelius T. Leondes (1996). Implementación de sistemas de control digital y técnicas computacionales. Academic Press. p. 123. ISBN 978-0-12-012779-5.
  8. ^ Eliahu Ibrahim Jury (1958). Sistemas de control de datos muestreados . John Wiley & Sons.
  9. ^ Eliahu Ibrahim Jury (1973). Teoría y aplicación del método de la transformada Z. Krieger Pub Co. ISBN 0-88275-122-0.
  10. ^ Eliahu Ibrahim Jury (1964). Teoría y aplicación del método de la transformada Z. John Wiley & Sons. pág. 1.
  11. ^ Jackson, Leland B. (1996). "La transformada z". Filtros digitales y procesamiento de señales . Boston, MA: Springer US. págs. 29-54. doi :10.1007/978-1-4757-2458-5_3. ISBN 978-1-4419-5153-3La transformada z es para los sistemas de tiempo discreto lo que la transformada de Laplace es para los sistemas de tiempo continuo. z es una variable compleja. A veces se la denomina transformada z bilateral , siendo la transformada z unilateral la misma excepto por una suma desde n = 0 hasta el infinito. El uso principal de la transformada unilateral ... es para secuencias causales, en cuyo caso las dos transformadas son las mismas de todos modos. Por lo tanto, no haremos esta distinción y nos referiremos a ... simplemente como la transformada z de x ( n ).
  12. ^ Bolzern, Paolo; Scattolini, Ricardo; Schiavoni, Nicola (2015). Fondamenti di Controlli Automatici (en italiano). Educación de MC Graw Hill. ISBN 978-88-386-6882-1.
  13. ^ abc AR Forouzan (2016). "Región de convergencia de la derivada de la transformada Z". Electronics Letters . 52 (8): 617–619. Bibcode :2016ElL....52..617F. doi :10.1049/el.2016.0189. S2CID  124802942.

Lectura adicional

Enlaces externos