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.
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]
Una extensión notable, conocida como la transformada Z modificada o avanzada , fue introducida posteriormente 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 integral 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]
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]
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 .
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.
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 da 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.
Los tres métodos siguientes se utilizan a menudo para la evaluación de la transformada inversa:
Este método implica aplicar el teorema de residuos de Cauchy para evaluar la transformada inversa. Al integrar alrededor de un contorno cerrado en el plano complejo, se suman los residuos en los polos de la función de transformada dentro de la ROC. Esta técnica es particularmente útil cuando se trabaja con funciones expresadas en términos de variables complejas.
En este método, la transformada se expande en una serie de potencias. Este enfoque es útil cuando la función de transformada es racional, lo que permite la aproximación de la inversa mediante la expansión en una serie y la determinación de los coeficientes de la señal término por término.
Esta técnica descompone la transformada en una suma de fracciones más simples, cada una correspondiente a pares de transformadas conocidos. Luego, la transformada inversa se determina buscando cada término en una tabla estándar de pares de transformadas. Este método se usa ampliamente por su eficiencia y simplicidad, especialmente cuando la función original se puede descomponer fácilmente en componentes reconocibles.
A) Determine la transformada inversa de la siguiente ecuación mediante el método de expansión en serie,
Solución:
Caso 1:
República Checa:
Dado que la ROC es el exterior de un círculo, es causal (señal existente para n≥0).
de este modo,
(la flecha indica el término en x(0)=1)
Tenga en cuenta que en cada paso del proceso de división larga eliminamos el término de menor potencia de .
Caso 2:
República Checa:
Dado que la ROC es el interior de un círculo, es anticausal (señal existente para n<0).
Realizando una división larga obtenemos,
(la flecha indica el término en x(0)=0)
Tenga en cuenta que en cada paso del proceso de división larga eliminamos el término de menor potencia de .
Nota:
B) Determine la transformada inversa de la siguiente ecuación mediante el método de expansión en serie,
Eliminando potencias negativas si y dividiendo por ,
Por expansión de fracciones parciales,
=
Caso 1:
República Checa:
Ambos términos son causales, por lo tanto es causal.
Caso 2:
República Checa:
Ambos términos son anticausales, por lo tanto es anticausal.
Caso 3:
República Checa:
Uno de los términos es causal (p=0,5 proporciona la parte causal) y el otro es anticausal (p=1 proporciona la parte anticausal), por lo tanto, es bidireccional.
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):
Ampliando el intervalo se obtiene
Mirando la suma
Por lo tanto, no hay valores de que satisfagan esta condición.
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".
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.
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
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.
Teorema del valor inicial : Sies causal, entonces
Teorema del valor final : Si los polos deestán dentro del círculo unitario, entonces
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.
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 ).
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.
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 .
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.
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 :
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. Al hacerlo, se obtendría la respuesta al impulso y la ecuación diferencial de coeficientes constantes lineales del sistema.
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.
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.
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.
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.
transformada z es a los sistemas de tiempo discreto lo que la transformada de Laplace es a 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 ).