Transformada Z
# Transformada Z

En matemáticas y en el procesamiento de señales, la transformada Z convierte una señal real o compleja definida en el dominio del tiempo discreto en una representación en el dominio de la frecuencia compleja.

El nombre de transformada Z procede de la variable del dominio, al igual que se podría llamar "Transformada S" a la Transformada de Laplace. Un nombre más adecuado para la TZ podría haber sido "Transformada de Laurent", ya que está basada en la serie de Laurent. La TZ es a las señales de tiempo discreto lo mismo que Laplace a las señales de tiempo continuo.

## Definición

La transformada Z, igual que otras transformaciones integrales, puede ser definida como una transformada unilateral o bilateral.

### Transformada Z bilateral

La TZ bilateral de una señal definida en el dominio del tiempo discreto $x[n]$ es una función $X(z)$ que se define

    $X(z)=Z\{x[n]\}=\sum _{n=-\infty }^{\infty }x[n]z^{-n}\$

donde $n\in \mathbb {Z}$ y en general $z\in \mathbb {C}$, es decir, es un número complejo de la forma

    $z=Ae^{j\omega }$

donde $A$ es el módulo de $z$, y $\omega$ es el argumento de ese complejo que bien podría representar la frecuencia angular (pulsación) en radianes por segundo (rad/s).

### Transformada Z unilateral

De forma alternativa, en los casos en que $x[n]$ está definida únicamente para $n\geq 0$, la transformada Z unilateral se define como

    $X^{+}(z)=Z^{+}\{x[n]\}=\sum _{n=0}^{\infty }x[n]z^{-n}\$

En el procesamiento de señales, se usa esta definición cuando la señal es causal. En este caso, la Transformada Z resulta una serie de Laurent, con Región de Convergencia (ROC) del tipo $|z|>R$ ; es decir que converge "hacia afuera".

Un ejemplo interesante de la TZ unilateral es la función de generadora de probabilidades, donde $x[n]$ es la probabilidad que toma una variable discreta aleatoria en el instante $n$, y la función $X(z)$ suele escribirse como $X(s)$, ya que $s=z^{-1}$. Las propiedades de las transformadas Z son útiles en la teoría de la probabilidad.

## Transformada Z inversa

La Transformada Z inversa se define

    $x[n]=Z^{-1}\{X(z)\}={\frac {1}{2\pi j}}\oint _{C}X(z)z^{n-1}dz\$

donde $C\$ es un círculo cerrado que envuelve el origen y la región de convergencia (ROC). El contorno, $C\$, debe contener todos los polos de $X(z)\$.

Un caso especial y simple de esta integral circular es que cuando $C\$ es el círculo unidad (que también puede usarse cuando la ROC incluye el círculo unidad), obtenemos la transformada inversa de tiempo discreto de Fourier:

    $x[n]={\frac {1}{2\pi }}\int _{-\pi }^{+\pi }X(e^{j\omega })e^{j\omega n}d\omega \$

La TZ con un rango finito de $n$ y un número finito de $z$ separadas de forma uniforme puede ser procesada de forma eficiente con el algoritmo de Bluestein. La transformada discreta de Fourier (DFT) es un caso especial de la TZ, y se obtiene limitando z para que coincida con el círculo unidad.

## Región de convergencia (ROC)

La región de convergencia, también conocida como ROC, define la región donde la transformada-z existe. La ROC es una región del plano complejo donde la TZ de una señal tiene una suma finita. La ROC para una x[n] es definida como el rango de z para la cual la transformada-z converge. Ya que la transformada–z es una serie de potencia, converge cuando $x[n]z^{-n}$ es absolutamente sumable.

    $ROC=\{z:\sum _{n=-\infty }^{\infty }x[n]z^{-n}<\infty \}\$

Propiedades de la Región de Convergencia:

La región de convergencia tiene propiedades que dependen de las características de la señal, $x[n]$.

* La ROC no tiene que contener algún polo. Por definición un polo es donde x[z] es infinito. Ya que x[z] tiene que ser finita para todas las z para tener convergencia, no puede existir ningún polo para ROC.
* Si $x[n]$ es una secuencia de duración finita, entonces la ROC es todo el plano-z, excepto en |z|=0 o |z|=∞.
* Si $x[n]$ es una secuencia del lado derecho entonces la ROC se extiende hacia fuera en el último polo desde x[z].
* Si $x[n]$ es una secuencia del lado izquierdo, entonces la ROC se extiende hacia dentro desde el polo más cercano en x[z].
* Si $x[n]$ es una secuencia con dos lados, la ROC va a ser un anillo en el plano-z que está restringida en su interior y exterior por un polo.

### Ejemplo 1 (Sin ROC)

Sea $x[n]=0.5^{n}\$. Expandiendo $x[n]\$ en $(-\infty ,\infty )\$ obtenemos

    $x[n]=\{...,0.5^{-3},0.5^{-2},0.5^{-1},1,0.5,0.5^{2},0.5^{3},...\}=\{...,2^{3},2^{2},2,1,0.5,0.5^{2},0.5^{3},...\}\$

Siendo la suma

    $\sum _{n=-\infty }^{\infty }x[n]z^{-n}<\infty \$

No hay ningún valor de $z\$ que satisfaga esta condición.

### Ejemplo 2 (ROC causal)

Sea $x[n]=0.5^{n}u[n]\$ (donde $u$ es la función escalón). Expandiendo $x[n]\$ en $(-\infty ,\infty )\$ obtenemos

    $x[n]=\{...,0,0,0,1,0.5,0.5^{2},0.5^{3},...\}\$

Siendo la suma

    $\sum _{n=-\infty }^{\infty }x[n]z^{-n}=\sum _{n=0}^{\infty }0.5^{n}z^{-n}=\sum _{n=0}^{\infty }\left({\frac {0.5}{z}}\right)^{n}={\frac {1}{1-0.5z^{-1}}}\$

La última igualdad se obtiene con la fórmula del sumatorio para series geométricas, y la igualdad solo se conserva si $\left|0.5z^{-1}\right|<1\$, lo cual puede ser reescrito para definir $z\$ de modo $\left|z\right|>0.5\$. Por lo tanto, la ROC es $\left|z\right|>0.5\$. En este caso la ROC es el plano complejo exterior al círculo de radio 0,5 con origen en el centro.

### Ejemplo 3 (ROC anticausal)

Sea $x[n]=-(0.5)^{n}u[-n-1]\$ (donde $u$ es la función escalón). Expandiendo $x[n]\$ entre $(-\infty ,\infty )\$ obtenemos

    $x[n]=\{...,-(0.5)^{-3},-(0.5)^{-2},-(0.5)^{-1},0,0,0,...\}\$

Siendo la suma

    $\sum _{n=-\infty }^{\infty }x[n]z^{-n}=-\sum _{n=-\infty }^{-1}0.5^{n}z^{-n}=-\sum _{n=-\infty }^{-1}\left({\frac {z}{0.5}}\right)^{-n}\$
    $=-\sum _{m=1}^{\infty }\left({\frac {z}{0.5}}\right)^{m}=-{\frac {0.5^{-1}z}{1-0.5^{-1}z}}={\frac {z}{z-0.5}}={\frac {1}{1-0.5z^{-1}}}\$

De nuevo, usando la fórmula de sumatorio para series geométricas, la igualdad solo se mantiene si $\left|0.5^{-1}z\right|<1\$, de modo que podemos definir $z\$ como $\left|z\right|<0.5\$.

Aquí, la ROC es $\left|z\right|<0.5\$, es decir, el interior de un círculo centrado en el origen de radio 0,5.

### Conclusión de los ejemplos

Los ejemplos 2 y 3 muestran claramente que la transformada $X(z)\$ de $x[n]\$ es única si y solo si se especifica cuál es la ROC. Dibujando los gráficos de polos y ceros para los casos causal y anticausal, comprobaríamos como la ROC de ambos casos no incluye el polo que está en 0,5. Esto se extiende a los casos con múltiples polos: la ROC nunca contiene polos.

En el ejemplo 2, el sistema causal tiene una ROC que incluye $\left|z\right|=\infty \$, mientras que al sistema anticausal del ejemplo 3 le pertenece una ROC que incluye $\left|z\right|=0\$.

En los sistemas con múltiples polos, es posible tener una ROC que no incluya ni $\left|z\right|=\infty \$ ni $\left|z\right|=0\$. La ROC crea una región circular. Por ejemplo, $x[n]=0.5^{n}u[n]-0.75^{n}u[-n-1]\$ tiene dos polos en 0,5 y 0,75. La ROC será $0.5<\left|z\right|<0.75\$, la cual no incluye ni el origen ni el infinito. Este tipo de sistemas se conoce como sistemas de causalidades mezcladas, ya que contiene un término causal $0.5^{n}u[n]\$ y otro anticausal $-(0.75)^{n}u[-n-1]\$.

La estabilidad de un sistema se puede determinar simplemente conociendo su ROC. Si esta ROC contiene el círculo unidad (p. ej. $\left|z\right|=1\$) entonces el sistema es estable. En los sistemas anteriores, el sistema causal es estable porque $\left|z\right|>0.5\$ contiene el círculo unidad.

Si tenemos la TZ de un sistema sin su ROC (p.ej., un $x[n]\$ ambiguo) podemos determinar una única señal $x[n]\$ en función de que queramos o no las siguientes propiedades:

* Estabilidad
* Causalidad

Si queremos un sistema estable, la ROC debe contener el círculo unidad. Si queremos un sistema causal, la ROC debe contener al infinito. Si queremos un sistema anticausal, la ROC debe contener al origen.

De este modo, podemos encontrar una señal en el tiempo $x[n]\$ que sea única.

## Propiedades

Propiedades de la Transformada Z :

| | Dominio del tiempo | Espacio Z | Demostración | ROC |
| - | - | - | - | - |
| Notación | $x[n]={\mathcal {Z}}^{-1}\{X(z)\}$ | $X(z)={\mathcal {Z}}\{x[n]\}$ | | ROC: $r_{2}<|z|<r_{1}$ |
| Linealidad | $a_{1}x_{1}[n]+a_{2}x_{2}[n]$ | $a_{1}X_{1}(z)+a_{2}X_{2}(z)$ | ${\begin{aligned}X(z)&=\sum _{n=-\infty }^{\infty }(a_{1}x_{1}(n)+a_{2}x_{2}(n))z^{-n}\\&=a_{1}\sum _{n=-\infty }^{\infty }x_{1}(n)z^{-n}+a_{2}\sum _{n=-\infty }^{\infty }x_{2}(n)z^{-n}\\&=a_{1}X_{1}(z)+a_{2}X_{2}(z)\end{aligned}}$ | Al menos la intersección de $ROC_1$ y $ROC_2$ |
| Expansión en el tiempo | $x_{(k)}[n]={\begin{cases}x[r],&n=rk\\0,&n\not =rk\end{cases}}$ r: integral | $X(z^{k})$ | ${\begin{aligned}X_{k}(z)&=\sum _{n=-\infty }^{\infty }x_{k}(n)z^{-n}\\&=\sum _{r=-\infty }^{\infty }x(r)z^{-rk}\\&=\sum _{r=-\infty }^{\infty }x(r)(z^{k})^{-r}\\&=X(z^{k})\end{aligned}}$ | $R^{\frac {1}{k}}$ |
| Desplazamiento en el tiempo | $x[n-k]$ | $z^{-k}X(z)$ | ${\begin{aligned}Z\{x[n-k]\}&=\sum _{n=0}^{\infty }x[n-k]z^{-n}\\&=\sum _{j=-k}^{\infty }x[j]z^{-(j+k)}&&j=n-k\\&=\sum _{j=-k}^{\infty }x[j]z^{-j}z^{-k}\\&=z^{-k}\sum _{j=-k}^{\infty }x[j]z^{-j}\\&=z^{-k}\sum _{j=0}^{\infty }x[j]z^{-j}&&x[\beta ]=0,\beta <0\\&=z^{-k}X(z)\end{aligned}}$ | ROC, excepto z =0 si k > 0 y z = ∞ si k < 0 |
| Escalamiento en el espacio Z | $a^{n}x[n]\$ | $X(a^{-1}z)\$ | ${\begin{aligned}{\mathcal {Z}}\{a^{n}x[n]\}&=\sum _{n=-\infty }^{\infty }a^{n}x(n)z^{-n}&\\&=\sum _{n=-\infty }^{\infty }x(n)(a^{-1}z)^{-n}\\&=X(a^{-1}z)\end{aligned}}$ | $|a|r_{2}<|z|<|a|r_{1}\$ |
| Inversión en el tiempo | $x[-n]$ | $X(z^{-1})$ | ${\begin{aligned}{\mathcal {Z}}\{x(-n)\}&=\sum _{n=-\infty }^{\infty }x(-n)z^{-n}\\&=\sum _{m=-\infty }^{\infty }x(m)z^{m}\\&=\sum _{m=-\infty }^{\infty }x(m){(z^{-1})}^{-m}\\&=X(z^{-1})\\\end{aligned}}$ | ${\frac {1}{r_{1}}}<|z|<{\frac {1}{r_{2}}}\$ |
| Conjugación compleja | $x^{*}[n]\$ | $X^{*}(z^{*})\$ | ${\begin{aligned}{\mathcal {Z}}\{x^{*}(n)\}&=\sum _{n=-\infty }^{\infty }x^{*}(n)z^{-n}\\&=\sum _{n=-\infty }^{\infty }\left[x(n)(z^{*})^{-n}\right]^{*}\\&=\left[\sum _{n=-\infty }^{\infty }x(n)(z^{*})^{-n}\right]^{*}\\&=X^{*}(z^{*})\end{aligned}}$ | ROC |
| Parte Real | $\operatorname {Re} \{x[n]\}$ | ${\frac {1}{2}}\left[X(z)+X^{*}(z^{*})\right]$ | | ROC |
| Parte Imaginaria | $\operatorname {Im} \{x[n]\}$ | ${\frac {1}{2j}}\left[X(z)-X^{*}(z^{*})\right]$ | | ROC |
| Diferenciación | $nx[n]$ | $-z{\frac {dX(z)}{dz}}$ | ${\begin{aligned}{\mathcal {Z}}\{nx(n)\}&=\sum _{n=-\infty }^{\infty }nx(n)z^{-n}\\&=z\sum _{n=-\infty }^{\infty }nx(n)z^{-n-1}\\&=-z\sum _{n=-\infty }^{\infty }x(n)(-nz^{-n-1})\\&=-z\sum _{n=-\infty }^{\infty }x(n){\frac {d}{dz}}(z^{-n})\\&=-z{\frac {dX(z)}{dz}}\end{aligned}}$ | ROC |
| Convolución | $x_{1}[n]*x_{2}[n]$ | $X_{1}(z)X_{2}(z)$ | ${\begin{aligned}{\mathcal {Z}}\{x_{1}(n)*x_{2}(n)\}&={\mathcal {Z}}\left\{\sum _{l=-\infty }^{\infty }x_{1}(l)x_{2}(n-l)\right\}\\&=\sum _{n=-\infty }^{\infty }\left[\sum _{l=-\infty }^{\infty }x_{1}(l)x_{2}(n-l)\right]z^{-n}\\&=\sum _{l=-\infty }^{\infty }x_{1}(l)\left[\sum _{n=-\infty }^{\infty }x_{2}(n-l)z^{-n}\right]\\&=\left[\sum _{l=-\infty }^{\infty }x_{1}(l)z^{-l}\right]\!\!\left[\sum _{n=-\infty }^{\infty }x_{2}(n)z^{-n}\right]\\&=X_{1}(z)X_{2}(z)\end{aligned}}$ | Al menos la intersección de $ROC_1$ y $ROC_2$ |
| Correlación cruzada | $r_{x_{1},x_{2}}=x_{1}^{*}[-n]*x_{2}[n]$ | $R_{x_{1},x_{2}}(z)=X_{1}^{*}({\tfrac {1}{z^{*}}})X_{2}(z)$ | | Al menos la intersección de ROC de $X_{1}({\tfrac {1}{z^{*}}})$ y $X_{2}(z)$ |
| Primera diferencia | $x[n]-x[n-1]$ | $(1-z^{-1})X(z)$ | | Al menos la intersección de ROC de $X_1(z)$ y $|z|>0$ |
| Acumulación | $\sum _{k=-\infty }^{n}x[k]$ | ${\frac {1}{1-z^{-1}}}X(z)$ | ${\begin{aligned}\sum _{n=-\infty }^{\infty }\sum _{k=-\infty }^{n}x[k]z^{-n}&=\sum _{n=-\infty }^{\infty }(x[n]+\cdots +x[-\infty ])z^{-n}\\&=X[z]\left(1+z^{-1}+z^{-2}+\cdots \right)\\&=X[z]\sum _{j=0}^{\infty }z^{-j}\\&=X[z]{\frac {1}{1-z^{-1}}}\end{aligned}}$ | |
| Multiplicación | $x_{1}[n]x_{2}[n]$ | ${\frac {1}{j2\pi }}\oint _{C}X_{1}(v)X_{2}({\tfrac {z}{v}})v^{-1}\mathrm {d} v$ | | - |
| Relación de Parseval | $\sum _{n=-\infty }^{\infty }x_{1}[n]x_{2}^{*}[n]$ | ${\frac {1}{j2\pi }}\oint _{C}X_{1}(v)X_{2}^{*}({\tfrac {1}{v^{*}}})v^{-1}\mathrm {d} v$ | | |

   Teorema de valor inicial Si x[n] causal
       $x[0]=\lim _{z\to \infty }X(z).$
   Teorema de valor final Si los polos de $(z-1)X(z)$ están dentro del círculo unitario, entonces
       $x[\infty ]=\lim _{z\to 1}(z-1)X(z).$

## Tabla con los pares más habituales de la transformada Z

| | Señal, $x(n)$ | Transformada Z, $X(z)$ | ROC |
| - | - | - | - |
| 1 | $\delta (n)\,$ | $1\,$ | ${\mbox{todo }}z\,$ |
| 2 | $u(n)\,$ | ${\frac {1}{1-z^{-1}}}$ | $|z|>1\,$ |
| 3 | $a^{n}u(n)\,$ | ${\frac {1}{1-az^{-1}}}$ | $|z|>|a|\,$ |
| 4 | $na^{n}u(n)\,$ | ${\frac {az^{-1}}{(1-az^{-1})^{2}}}$ | $|z|>|a|\,$ |
| 5 | $-a^{n}u(-n-1)\,$ | ${\frac {1}{1-az^{-1}}}$ | $|z|<|a|\,$ |
| 6 | $-na^{n}u(-n-1)\,$ | ${\frac {az^{-1}}{(1-az^{-1})^{2}}}$ | $|z|<|a|\,$ |
| 7 | $\cos(\omega _{0}n)u(n)\,$ | ${\frac {1-z^{-1}\cos(\omega _{0})}{1-2z^{-1}\cos(\omega _{0})+z^{-2}}}$ | $|z|>1\,$ |
| 8 | $\sin(\omega _{0}n)u(n)\,$ | ${\frac {z^{-1}\sin(\omega _{0})}{1-2z^{-1}\cos(\omega _{0})+z^{-2}}}$ | $|z|>1\,$ |
| 9 | $a^{n}\cos(\omega _{0}n)u(n)\,$ | ${\frac {1-az^{-1}\cos(\omega _{0})}{1-2az^{-1}\cos(\omega _{0})+a^{2}z^{-2}}}$ | $|z|>|a|\,$ |
| 10 | $a^{n}\sin(\omega _{0}n)u(n)\,$ | ${\frac {az^{-1}\sin(\omega _{0})}{1-2az^{-1}\cos(\omega _{0})+a^{2}z^{-2}}}$ | $|z|>|a|\,$ |

## Relación con Laplace

La TZ bilateral es simplemente la transformada de Laplace bilateral de la señal muestreada

       $x(t)\sum _{n=-\infty }^{\infty }\delta (t-nT)=\sum _{n=-\infty }^{\infty }x[n]\delta (t-nT)\$

donde $x(t)\$ es la señal continua muestreada, $x[n]=x(nT)\$ la n-ésima muestra, $T\$ el período de muestreo, y con la sustitución $z=e^{sT}\$.

Del mismo modo, la TZ unilateral es simplemente la transformada de Laplace unilateral de la señal ideal muestreada. En ambas se asume que la señal muestreada vale cero para todos los índices negativos en el tiempo.

## Relación con Fourier

La TZ es una generalización de la transformada de Fourier de tiempo discreto (DTFT). La DTFT puede hallarse evaluando la TZ $X(z)\$ en $z=e^{j\omega }\$ o, lo que es lo mismo, evaluada en el círculo unidad. Para determinar la respuesta en frecuencia del sistema, la TZ debe ser evaluada en el círculo unidad.

## Ecuación en diferencias de coeficientes lineales constantes

La ecuación diferencial de coeficientes lineales constantes (en inglés LCCD) es una representación de un sistema lineal, basada en la ecuación de la media autorregresiva.

    $\sum _{p=0}^{N}y[n-p]\alpha _{p}=\sum _{q=0}^{M}x[n-q]\beta _{q}\$

Ambos términos de esta ecuación pueden dividirse por $\alpha _{0}\$, si no es cero, normalizando $\alpha _{0}=1\$ la ecuación LCCD puede ser escrita

    $y[n]=\sum _{q=0}^{M}x[n-q]\beta _{q}-\sum _{p=1}^{N}y[n-p]\alpha _{p}\$

Esta forma de la ecuación LCCD es más explícita para comprobar que la salida actual $y[{n}]\$ se define en función de las salidas anteriores $y[{n-p}]\$, la entrada actual $x[{n}]\$, y las entradas anteriores $x[{n-q}]\$.

### Función de transferencia

Se calcula haciendo la TZ de la ecuación

    $Y(z)\sum _{p=0}^{N}z^{-p}\alpha _{p}=X(z)\sum _{q=0}^{M}z^{-q}\beta _{q}\$

y dividiendo

    $H(z)={\frac {Y(z)}{X(z)}}={\frac {\sum _{q=0}^{M}z^{-q}\beta _{q}}{\sum _{p=0}^{N}z^{-p}\alpha _{p}}}={\frac {\beta _{0}+z^{-1}\beta _{1}+z^{-2}\beta _{2}+...+z^{-M}\beta _{M}}{\alpha _{0}+z^{-1}\alpha _{1}+z^{-2}\alpha _{2}+...+z^{-N}\alpha _{N}}}\$

### Ceros y polos

Gracias al teorema fundamental del álgebra sabemos que el numerador tiene M raíces (llamadas ceros) y el denominador tiene N raíces (llamadas polos). Factorizando la función de transferencia

    $H(z)={\frac {(1-q_{1}z^{-1})(1-q_{2}z^{-1})...(1-q_{M}z^{-1})}{(1-p_{1}z^{-1})(1-p_{2}z^{-1})...(1-p_{N}z^{-1})}}\$

donde $q_{k}\$ es el k-ésimo cero y $p_{k}\$ es el k-ésimo polo. Los ceros y polos son por lo general complejos, y por tanto se pueden dibujar en el plano complejo.

En definitiva, los ceros son las soluciones de la ecuación obtenida de igualar el numerador a cero, mientras que los polos son las de la ecuación que se obtiene al igualar a cero el denominador.

Se puede factorizar el denominador mediante la descomposición en fracciones simples, las cuales pueden ser transformadas de nuevo al dominio del tiempo. Haciendo esto obtenemos la respuesta al impulso y la ecuación diferencial de coeficientes lineales constantes del sistema.

### Salida del sistema

Si por un sistema $H(z)\$ pasa una señal $X(z)\$ entonces la salida será $Y(z)=H(z)X(z)\$. Haciendo una descomposición en fracciones simples de $Y(z)\$ y la TZ inversa de cada una de ellas puede encontrarse entonces la salida $y[n]\$.

