Articulo de referencia

Transformación 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...

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 valor complejo en el dominio de la frecuencia (el dominio z o plano z ). [ 1 ] [ 2 ] [ 3 ]

Puede considerarse una contraparte en tiempo discreto de la transformada de Laplace (el dominio s o plano s ). [ 4 ] 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 corresponde con el área dentro del círculo unitario del dominio z, mientras que el semiplano derecho del dominio s se corresponde con el área fuera del círculo unitario del dominio z.

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

Historia

La idea básica ahora conocida como la transformada Z era conocida por Laplace , y fue reintroducida en 1947 por W. Hurewicz [ 5 ] [ 6 ] y otros como una forma de tratar los sistemas de control de datos muestreados utilizados con radar. Proporciona una manera manejable de resolver ecuaciones de diferencias lineales con coeficientes constantes. Posteriormente, Ragazzini y Zadeh, del grupo de control de datos muestreados de la Universidad de Columbia, la denominaron "la transformada Z" en 1952. [ 7 ] [ 8 ]

La transformada Z modificada o avanzada fue desarrollada y popularizada posteriormente por EI Jury . [ 9 ] [ 10 ]

La idea contenida en la transformada Z también se conoce en la literatura matemática como el método de generación de funciones , que se remonta a 1730, cuando fue introducido por de Moivre junto con la teoría de la probabilidad . [ 11 ] Desde un punto de vista matemático, la transformada Z también puede verse como una serie de Laurent , donde se considera la secuencia de números en cuestión como la expansión (de Laurent) de una función analítica.

Definición

La transformada Z puede definirse como una transformada unilateral o bilateral (de forma similar a la transformada de Laplace unilateral y la transformada de Laplace bilateral ). [ 12 ]

Transformación Z bilateral

La transformada Z bilateral o de dos lados de una señal de tiempo discreto.incógnita[norte]{\displaystyle x[n]}es la serie de potencias formalincógnita(z){\displaystyle X(z)}definido como:

incógnita(z)=Z{incógnita[norte]}=norte=incógnita[norte]znorte{\displaystyle X(z)={\mathcal {Z}}\{x[n]\}=\sum _{n=-\infty }^{\infty }x[n]z^{-n}}

dóndenorte{\displaystyle n}es un número entero yz{\displaystyle z}es, en general, un número complejo . En forma polar ,z{\displaystyle z}puede escribirse como:

z=Amiiϕ=A(porqueϕ+ipecadoϕ){\displaystyle z=Ae^{i\phi }=A\cdot (\cos {\phi }+i\sin {\phi })}

dóndeA{\displaystyle A}es la magnitud dez{\displaystyle z} ,i{\displaystyle i}es la unidad imaginaria yϕ{\displaystyle \phi }es el argumento complejo (también denominado ángulo o fase ) en radianes .

Transformación Z unilateral

Alternativamente, en los casos en queincógnita[norte]{\displaystyle x[n]}se define únicamente paranorte0{\displaystyle n\geq 0} , latransformada Z unilateral o de un solo lado se define como:

incógnita(z)=Z{incógnita[norte]}=norte=0incógnita[norte]znorte.{\displaystyle X(z)={\mathcal {Z}}\{x[n]\}=\sum _{n=0}^{\infty }x[n]z^{-n}.}

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

Un ejemplo importante de la transformada Z unilateral es la función generadora de probabilidad , donde el componenteincógnita[norte]{\displaystyle x[n]}es la probabilidad de que una variable aleatoria discreta tome el valor norte{\displaystyle n}Las propiedades de las transformadas Z (enumeradas en la sección  Propiedades ) tienen interpretaciones útiles en el contexto de la teoría de la probabilidad .

Transformada Z inversa

La transformada Z inversa es:

incógnita[norte]=Z1{incógnita(z)}=12πidoincógnita(z)znorte1dz{\displaystyle x[n]={\mathcal {Z}}^{-1}\{X(z)\}={\frac {1}{2\pi i}}\oint _ {C}X(z)z^{n-1}dz}

dóndedo{\displaystyle C}es una trayectoria cerrada en sentido contrario a las agujas del reloj que rodea el origen y se encuentra completamente en la región de convergencia (ROC). En el caso de que la ROC sea causal (véase el Ejemplo 2 ), esto significa que la trayectoriado{\displaystyle C}debe rodear todos los polos de incógnita(z){\displaystyle X(z)}.

Un caso especial de esta integral de contorno ocurre cuandodo{\displaystyle C}es el círculo unitario. Este contorno se puede utilizar cuando la ROC incluye el círculo unitario, lo cual siempre está garantizado cuandoincógnita(z){\displaystyle X(z)}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 discreta inversa , o serie de Fourier , de los valores periódicos de la transformada Z alrededor del círculo unitario:

incógnita[norte]=12πππincógnita(miiω)miiωnortedω.{\displaystyle x[n]={\frac {1}{2\pi }}\int _{-\pi }^{\pi }X(e^{i\omega })e^{i\omega n}d\omega .}

La transformada Z con un rango finito denorte{\displaystyle n}y un número finito de espaciados uniformementez{\displaystyle z}Los valores se pueden calcular de manera eficiente mediante el algoritmo FFT de Bluestein . La transformada discreta de Fourier (DTFT), que no debe confundirse con la transformada discreta de Fourier (DFT), es un caso especial de dicha transformada Z obtenida al restringirz{\displaystyle z}estar situado sobre el círculo unitario.

Los siguientes tres métodos se utilizan a menudo para la evaluación de la transformada inversa:

Evaluación directa mediante integración de contornos

Este método consiste en aplicar el teorema de los residuos de Cauchy para evaluar la transformada Z inversa. Al integrar alrededor de un contorno cerrado en el plano complejo, se suman los residuos en los polos de la función de la transformada Z dentro de la región de convergencia. Esta técnica resulta especialmente útil al trabajar con funciones expresadas en términos de variables complejas.

Expansión en una serie de términos en las variables z y z −1

En este método, la transformada Z se expande en una serie de potencias. Este enfoque resulta útil cuando la función de la transformada Z es racional, lo que permite aproximar la inversa mediante la expansión en una serie y la determinación de los coeficientes de la señal término a término.

Expansión de fracciones parciales y búsqueda en tablas

Esta técnica descompone la transformada Z en una suma de fracciones más simples, cada una correspondiente a pares de transformadas Z conocidos. La transformada Z inversa se determina consultando cada término en una tabla estándar de pares de transformadas Z. Este método se utiliza ampliamente por su eficiencia y simplicidad, especialmente cuando la función original se puede descomponer fácilmente en componentes reconocibles.

Ejemplo

[ 13 ]

A) Determine la transformada Z inversa de la siguiente expresión mediante el método de expansión en serie: incógnita(z)=111.5z1+0,5z2{\displaystyle X(z)={\frac {1}{1-1.5\,z^{-1}+0.5\,z^{-2}}}}

Solución:

Caso 1:

ROC:|z|>1{\displaystyle \left\vert z\right\vert >1}

Dado que la ROC es el exterior de un círculo,incógnita(norte){\displaystyle x(n)}es causal (señal existente paranorte0{\displaystyle n\geq 0} ). incógnita(z)=1132z1+12z2=1+32z1+74z2+158z3+3116z4+.{\displaystyle X(z)={1 \over 1-{3 \over 2}z^{-1}+{1 \over 2}z^{-2}}=1+{{3 \over 2}z^{-1}}+{{7 \over 4}z^{-2}}+{{15 \over 8}z^{-3}}+{{31 \over 16}z^{-4}}+\ldots .} De este modo, incógnita(norte)={1,32,74,158,3116}{\displaystyle {\begin{aligned}x(n)&=\left\{1,{\frac {3}{2}},{\frac {7}{4}},{\frac {15}{8}},{\frac {31}{16}}\ldots \right\}\\&\qquad \!\uparrow \\\end{aligned}}}(la flecha indica el término en incógnita(0)=1{\displaystyle x(0)=1} ).

Tenga en cuenta que en cada paso del proceso de división larga eliminamos el término de menor potencia de z1{\displaystyle z^{-1}}.

Caso 2:

ROC:|z|<0,5{\displaystyle \left\vert z\right\vert <0.5}

Dado que la ROC es el interior de un círculo,incógnita(norte){\displaystyle x(n)}es anticausal (señal existente paranorte<0{\displaystyle n<0} ).

Al realizar la división larga obtenemos incógnita(z)=1132z1+12z2=2z2+6z3+14z4+30z5+{\displaystyle X(z)={\frac {1}{1-{\frac {3}{2}}z^{-1}+{\frac {1}{2}}z^{-2}}}=2z^{2}+6z^{3}+14z^{4}+30z^{5}+\ldots }

incógnita(norte)={30,14,6,2,0,0}  {\displaystyle {\begin{aligned}x(n)&=\{30,14,6,2,0,0\}\\&\qquad \qquad \qquad \quad \ \ \,\uparrow \\\end{aligned}}}(la flecha indica el término en incógnita(0)=0{\displaystyle x(0)=0} ).

Tenga en cuenta que en cada paso del proceso de división larga eliminamos el término de menor potencia de z{\displaystyle z}.

Nota:

  1. Cuando la señal es causal, obtenemos potencias positivas dez{\displaystyle z}y cuando la señal es anticausal, obtenemos potencias negativas de z{\displaystyle z}.
  2. zk{\displaystyle z^{k}}indica término enincógnita(k){\displaystyle x(-k)}yzk{\displaystyle z^{-k}}indica el término en incógnita(k){\displaystyle x(k)}.

B) Determine la transformada Z inversa de lo siguiente mediante el método de expansión en serie:

Eliminar poderes negativos siz{\displaystyle z}y dividiendo por z{\displaystyle z} , incógnita(z)z=z2z(z21.5z+0,5)=zz21.5z+0,5{\displaystyle {\frac {X(z)}{z}}={\frac {z^{2}}{z(z^{2}-1.5\,z+0.5)}}={\frac {z}{z^{2}-1.5\,z+0.5}}}

Mediante expansión en fracciones parciales, incógnita(z)z=z(z1)(z0,5)=A1z0,5+A2z1A1=(z0,5)incógnita(z)z|z=0,5=0,5(0,51)=1A2=(z1)incógnita(z)z|z=1=110,5=2incógnita(z)z=2z11z0,5{\displaystyle {\begin{aligned}{\frac {X(z)}{z}}&={\frac {z}{(z-1)(z-0.5)}}={\frac {A_{1}}{z-0.5}}+{\frac {A_{2}}{z-1}}\\[4pt]&A_{1}=\left.{\frac {(z-0.5)X(z)}{z}}\right\vert _{z=0.5}={\frac {0.5}{(0.5-1)}}=-1\\[4pt]&A_{2}=\left.{\frac {(z-1)X(z)}{z}}\right\vert _{z=1}={\frac {1}{1-0.5}}={2}\\[4pt]{\frac {X(z)}{z}}&={\frac {2}{z-1}}-{\frac {1}{z-0.5}}\end{aligned}}}

Caso 1:

ROC:|z|>1{\displaystyle \left\vert z\right\vert >1}

Ambos términos son causales, por lo tantoincógnita(norte){\displaystyle x(n)}es causal.

incógnita(norte)=2(1)norte(norte)1(0,5)norte(norte)=(2(0,5)norte)(norte){\displaystyle {\begin{aligned}x(n)&=2{(1)^{n}}u(n)-1{(0.5)^{n}}u(n)\\&=(2-(0.5)^{n})u(n)\\\end{aligned}}}

Caso 2:

ROC:|z|<0,5{\displaystyle \left\vert z\right\vert <0.5}

Ambos términos son anticausales, por lo tantoincógnita(norte){\displaystyle x(n)}es anticarial.

incógnita(norte)=2(1)norte(norte1)(1(0,5)norte(norte1))=((0,5)norte2)(norte1){\displaystyle {\begin{aligned}x(n)&=-2{(1)^{n}}u(-n-1)-(-1{(0.5)^{n}}u(-n-1))\\&=((0.5)^{n}-2)u(-n-1)\\\end{aligned}}}

Caso 3:

ROC:0,5<|z|<1{\displaystyle 0.5<\left\vert z\right\vert <1}

Uno de los términos es causal ( pag=0,5{\displaystyle p=0.5} proporciona la parte causal) y el otro es anticausal (pag=1{\displaystyle p=1}proporciona la parte anticausal), por lo tantoincógnita(norte){\displaystyle x(n)}Es una cuestión de dos caras.

incógnita(norte)=2(1)norte(norte1)1(0,5)norte(norte)=2(norte1)(0,5)norte(norte){\displaystyle {\begin{aligned}x(n)&=-2{(1)^{n}}u(-n-1)-1{(0.5)^{n}}u(n)\\&=-2u(-n-1)-(0.5)^{n}u(n)\\\end{aligned}}}

Región de convergencia

La región de convergencia (ROC) es el conjunto de puntos en el plano complejo para los cuales la suma de la transformada Z converge absolutamente :

ROdo={z:norte=|incógnita[norte]znorte|<}{\displaystyle \mathrm {ROC} =\left\{z:\sum _{n=-\infty }^{\infty }\left|x[n]z^{-n}\right|<\infty \right\}}

Ejemplo 1 (sin ROC)

Dejemosincógnita[norte]=(0,5)norte{\displaystyle x[n]=(0.5)^{n}}. Ampliandoincógnita[norte]{\displaystyle x[n]}en el intervalo(,){\displaystyle (-\infty ,\infty )}se convierte

incógnita[norte]={,(0,5)3,(0,5)2,(0,5)1,1,0,5,(0,5)2,(0,5)3,}={,23,22,2,1,0,5,(0,5)2,(0,5)3,}.{\displaystyle x[n]=\left\{\dots ,(0.5)^{-3},(0.5)^{-2},(0.5)^{-1},1,0.5,(0.5)^{2},(0.5)^{3},\dots \right\}=\left\{\dots ,2^{3},2^{2},2,1,0.5,(0.5)^{2},(0.5)^{3},\dots \right\}.}

Al observar la suma

norte=incógnita[norte]znorte.{\displaystyle \sum _{n=-\infty }^{\infty }x[n]z^{-n}\to \infty .}

Por lo tanto, no hay valores dez{\displaystyle z}que cumplan esta condición.

Ejemplo 2 (ROC causal)

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

Dejarincógnita[norte]=(0,5)norte[norte]{\displaystyle x[n]=(0.5)^{n}\,u[n]}(dónde{\displaystyle u}es la función escalón de Heaviside ). Expandiendoincógnita[norte]{\displaystyle x[n]}en el intervalo(,){\displaystyle (-\infty ,\infty )}se convierte

incógnita[norte]={,0,0,0,1,0,5,(0,5)2,(0,5)3,}.{\displaystyle x[n]=\left\{\dots ,0,0,0,1,0.5,(0.5)^{2},(0.5)^{3},\dots \right\}.}

Al observar la suma

norte=incógnita[norte]znorte=norte=0(0,5)norteznorte=norte=0(0,5z)norte=110,5z1.{\displaystyle \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.5\,z^{-1}}}.}

La última igualdad surge de la serie geométrica infinita y la igualdad solo se cumple si |0,5z1|<1{\displaystyle \vert 0.5\,z^{-1}\vert <1} , que puede reescribirse en términos dez{\displaystyle z}como|z|>0,5{\displaystyle \vert z\vert >0.5} . Por lo tanto, el ROC es|z|>0,5{\displaystyle \vert z\vert >0.5} . En este caso, la ROC es el plano complejo con un disco de radio0,5{\displaystyle 0.5} en el origen "perforado".

Ejemplo 3 (ROC anticausa)

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

Dejarincógnita[norte]=(0,5)norte[norte1]{\displaystyle x[n]=-(0.5)^{n}\,u[-n-1]}(dónde{\displaystyle u}es la función escalón de Heaviside ). Expandiendoincógnita[norte]{\displaystyle x[n]}en el intervalo(,){\displaystyle (-\infty ,\infty )}se convierte

incógnita[norte]={,(0,5)3,(0,5)2,(0,5)1,0,0,0,0,}.{\displaystyle x[n]=\left\{\dots ,-(0.5)^{-3},-(0.5)^{-2},-(0.5)^{-1},0,0,0,0,\dots \right\}.}

Al observar la suma

norte=incógnita[norte]znorte=norte=1(0,5)norteznorte=metro=1(z0,5)metro=(0,5)1z1(0,5)1z=10,5z11=110,5z1{\displaystyle {\begin{aligned}\sum _{n=-\infty }^{\infty }x[n]\,z^{-n}&=-\sum _{n=-\infty }^{-1}(0.5)^{n}\,z^{-n}\\&=-\sum _{m=1}^{\infty }\left({\frac {z}{0.5}}\right)^{m}\\&=-{\frac {(0.5)^{-1}z}{1-(0.5)^{-1}z}}\\&=-{\frac {1}{0.5\,z^{-1}-1}}\\&={\frac {1}{1-0.5\,z^{-1}}}\\\end{aligned}}}

y utilizando de nuevo la serie geométrica infinita , la igualdad solo se cumple si|(0,5)1z|<1{\displaystyle |(0.5)^{-1}z|<1}que se puede reescribir en términos dez{\displaystyle z}como|z|<0,5{\displaystyle \vert z\vert <0.5} . Por lo tanto, el ROC es|z|<0,5{\displaystyle \vert z\vert <0.5} . En este caso, la ROC es un disco centrado en el origen y de radio 0,5{\displaystyle 0.5}.

Lo único que diferencia este ejemplo del anterior es la curva ROC. Esto se hizo intencionadamente para demostrar que el resultado de la transformación por sí solo es insuficiente.

Ejemplos de conclusión

Los ejemplos 2 y 3 muestran claramente que la transformada Zincógnita(z){\displaystyle X(z)}deincógnita[norte]{\displaystyle x[n]}es único cuando y solo cuando se especifica la ROC. La creación del diagrama de polos y ceros para el caso causal y anticausal muestra que la ROC para cualquiera de los casos no incluye el polo que está en 0,5{\displaystyle 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 una ROC que incluye|z|={\displaystyle |z|=\infty }mientras que el sistema anticausal del ejemplo 3 produce una ROC que incluye |z|=0{\displaystyle \vert z\vert =0}.

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

En sistemas con múltiples polos es posible tener una ROC que no incluya ninguno|z|={\displaystyle |z|=\infty }ni|z|=0{\displaystyle \vert z\vert =0} . La ROC crea una banda circular. Por ejemplo,

incógnita[norte]=(0,5)norte[norte](0,75)norte[norte1]{\displaystyle x[n]=(0.5)^{n}\,u[n]-(0.75)^{n}\,u[-n-1]}

tiene postes en0,5{\displaystyle 0.5}y0,75{\displaystyle 0.75} . El ROC será0,5<|z|<0,75{\displaystyle 0.5<\vert z\vert <0.75} , que no incluye ni el origen ni el infinito. Dicho sistema se denomina sistema de causalidad mixta, ya que contiene un término causal.(0,5)norte[norte]{\displaystyle (0.5)^{n}\,u[n]}y un término anticausa(0,75)norte[norte1]{\displaystyle -(0.75)^{n}\,u[-n-1]}.

La estabilidad de un sistema también se puede determinar conociendo solo la ROC. Si la ROC contiene el círculo unitario (es decir ,|z|=1{\displaystyle \vert z\vert =1} ) ​​entonces el sistema es estable. En los sistemas anteriores, el sistema causal (Ejemplo 2) es estable porque|z|>0,5{\displaystyle \vert z\vert >0.5}Contiene el círculo unitario.

Supongamos que se nos proporciona una transformada Z de un sistema sin ROC (es decir, una ambigua incógnita[norte]{\displaystyle x[n]}) . Podemos determinar un únicoincógnita[norte]{\displaystyle x[n]}siempre que deseemos lo siguiente:

  • Estabilidad
  • Causalidad

Para que haya estabilidad, la región de convergencia (ROC) debe contener el círculo unitario. Si necesitamos un sistema causal, la ROC debe contener el infinito y la función del sistema será una secuencia unilateral derecha. Si necesitamos un sistema anticausal, la ROC debe contener el origen y la función del sistema será una secuencia unilateral izquierda. Si necesitamos tanto estabilidad como causalidad, todos los polos de la función del sistema deben estar dentro del círculo unitario.

El únicoincógnita[norte]{\displaystyle x[n]}entonces se puede encontrar.

Propiedades

Teorema de Parseval

norte=incógnita1[norte]incógnita2[norte]=12πidoincógnita1(v)incógnita2(1v)v1dv{\displaystyle \sum _{n=-\infty }^{\infty }x_{1}[n]x_{2}^{*}[n]\quad =\quad {\frac {1}{2\pi i}}\oint _{C}X_{1}(v)X_{2}^{*}({\tfrac {1}{v^{*}}})v^{-1}\mathrm {d} v}

Teorema del valor inicial  : Siincógnita[norte]{\displaystyle x[n]}es causal, entoncesincógnita[0]=límitezincógnita(z).{\displaystyle x[0]=\lim _{z\to \infty }X(z).}

Teorema del valor final : Si los polos de(z1)incógnita(z){\displaystyle (z-1)X(z)}están dentro del círculo unitario, entoncesincógnita[]=límitez1(z1)incógnita(z).{\displaystyle x[\infty ]=\lim _{z\to 1}(z-1)X(z).}

Tabla de pares comunes de transformada Z

Aquí:

:norte[norte]={1,norte00,norte<0{\displaystyle u:n\mapsto u[n]={\begin{cases}1,&n\geq 0\\0,&n<0\end{cases}}}

es la función escalón unitario (o de Heaviside) y

δ:norteδ[norte]={1,norte=00,norte0{\displaystyle \delta :n\mapsto \delta [n]={\begin{cases}1,&n=0\\0,&n\neq 0\end{cases}}}

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

Relación con las series de Fourier y la transformada de Fourier

Para valores dez{\displaystyle z}en la región|z|=1{\displaystyle \vert z\vert =1} , conocido como el círculo unitario , podemos expresar la transformación como una función de una sola variable real.ω{\displaystyle \omega }definiendoz=miiω{\displaystyle z=e^{i\omega }} . La transformada bilateral se reduce a una serie de Fourier :

que también se conoce como la transformada discreta de Fourier en tiempo (DTFT) de laincógnita[norte]{\displaystyle x[n]}secuencia. Esta2π{\displaystyle 2\pi }- La 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, veamos...incógnita(F){\displaystyle X(f)}Sea la transformada de Fourier de cualquier función ,incógnita(t){\displaystyle x(t)} , cuyas muestras a cierto intervaloT{\displaystyle T}igual elincógnita[norte]{\displaystyle x[n]}secuencia. Luego la DTFT de laincógnita[norte]{\displaystyle x[n]}La secuencia se puede escribir de la siguiente manera.

dóndeT{\displaystyle T}tiene unidades de segundos,F{\displaystyle f}tiene unidades de hercios . La comparación de las dos series revela queω=2πFT{\displaystyle \omega =2\pi fT}es una frecuencia normalizada con unidad de radianes por muestra . El valorω=2π{\displaystyle \omega =2\pi }corresponde aF=1T{\displaystyle f={\tfrac {1}{T}}} . Y ahora, con la sustituciónF=ω2πT{\displaystyle f={\tfrac {\omega }{2\pi T}}} , La ecuación 1 puede expresarse en términos deincógnita(ω2πk2πT){\displaystyle X({\tfrac {\omega -2\pi k}{2\pi T}})}(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 . En la ecuación 3, sin embargo, los centros permanecen separados por , mientras que sus anchos se expanden o se contraen. Cuando la secuenciaincógnita(norteT){\displaystyle x(nT)}representa la respuesta impulsional de un sistema LTI , estas funciones también se conocen como su respuesta en frecuencia . Cuando elincógnita(norteT){\displaystyle x(nT)}La secuencia es periódica, su transformada discreta de Fourier (DTFT) diverge en una o más frecuencias armónicas y es cero en todas las demás. Esto se suele representar mediante el uso de funciones delta de Dirac con amplitud variable en las frecuencias armónicas. Debido a la periodicidad, solo existe un número finito de amplitudes únicas, que se calculan fácilmente mediante la transformada discreta de Fourier (DFT), mucho más sencilla. (Véase Transformada discreta de Fourier §  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:

s=2T(z1)(z+1){\displaystyle s={\frac {2}{T}}{\frac {(z-1)}{(z+1)}}}

para convertir alguna funciónH(s){\displaystyle H(s)}en el dominio de Laplace a una funciónH(z){\displaystyle H(z)}en el dominio Z ( transformación de Tustin ), o

z=misT1+sT/21sT/2{\displaystyle z=e^{sT}\approx {\frac {1+sT/2}{1-sT/2}}}

del dominio Z al dominio de Laplace. Mediante la transformación bilineal, el plano complejo s (de la transformada de Laplace) se mapea al plano complejo z (de la transformada z). Si bien este mapeo es (necesariamente) no lineal, es útil porque mapea todo el sistema de coordenadas.iω{\displaystyle i\omega }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 eliω{\displaystyle i\omega }eje) se convierte en la transformada de Fourier de tiempo discreto. Esto supone que la transformada de Fourier existe; es decir, que eliω{\displaystyle i\omega }El eje se encuentra en la región de convergencia de la transformada de Laplace.

Transformación estrellada

Dada una transformación Z unilateralincógnita(z){\displaystyle X(z)}de una función muestreada en el tiempo, la transformación marcada con asterisco correspondiente produce una transformada de Laplace y restaura la dependencia deT{\displaystyle T}(el parámetro de muestreo):

incógnita(s)=incógnita(z)|z=misT{\displaystyle {\bigg .}X^{*}(s)=X(z){\bigg |}_{\displaystyle z=e^{sT}}}

La transformada inversa de Laplace es una abstracción matemática conocida como función muestreada por impulsos .

Ecuación lineal en diferencias con coeficientes constantes

La ecuación de diferencias lineales con coeficientes constantes (LCCD) es una representación de un sistema lineal basada en la ecuación de media móvil autorregresiva :

pag=0nortey[nortepag]αpag=q=0METROincógnita[norteq]βq.{\displaystyle \sum _{p=0}^{N}y[n-p]\alpha _{p}=\sum _{q=0}^{M}x[n-q]\beta _{q}.}

Ambos lados de la ecuación anterior se pueden dividir porα0{\displaystyle \alpha _{0}}si no es cero. Normalizando con α0=1{\displaystyle \alpha _{0}=1} , la ecuación LCCD se puede escribir

y[norte]=q=0METROincógnita[norteq]βqpag=1nortey[nortepag]αpag.{\displaystyle 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 favorable para hacer más explícito que la salida "actual"y[norte]{\displaystyle y[n]}es una función de resultados anterioresy[nortepag]{\displaystyle y[n-p]}, entrada actualincógnita[norte]{\displaystyle x[n]}y entradas anterioresincógnita[norteq]{\displaystyle x[n-q]}.

Función de transferencia

Al aplicar la transformada Z a la ecuación anterior (utilizando las leyes de linealidad y desplazamiento temporal) se obtiene:

Y(z)pag=0nortezpagαpag=incógnita(z)q=0METROzqβq{\displaystyle Y(z)\sum _{p=0}^{N}z^{-p}\alpha _{p}=X(z)\sum _{q=0}^{M}z^{-q}\beta _{q}}

dóndeincógnita(z){\displaystyle X(z)}yY(z){\displaystyle Y(z)}son la transformada Z deincógnita[norte]{\displaystyle x[n]}yy[norte]{\displaystyle y[n]} , respectivamente. (Las convenciones de notación suelen usar letras mayúsculas para referirse a la transformada Z de una señal denotada por una letra minúscula correspondiente, de forma similar a la convención utilizada para notar las transformadas de Laplace).

Reorganizar los resultados en la función de transferencia del sistema :

H(z)=Y(z)incógnita(z)=q=0METROzqβqpag=0nortezpagαpag=β0+z1β1+z2β2++zMETROβMETROα0+z1α1+z2α2++znorteαnorte.{\displaystyle 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}+\cdots +z^{-M}\beta _{M}}{\alpha _{0}+z^{-1}\alpha _{1}+z^{-2}\alpha _{2}+\cdots +z^{-N}\alpha _{N}}}.}

Ceros y polos

Del teorema fundamental del álgebra el numerador tieneMETRO{\displaystyle M}raíces (correspondientes a ceros deH{\displaystyle H}) y el denominador tienenorte{\displaystyle N}raíces (correspondientes a polos). Reescribiendo la función de transferencia en términos de ceros y polos.

H(z)=(1q1z1)(1q2z1)(1qMETROz1)(1pag1z1)(1pag2z1)(1pagnortez1),{\displaystyle H(z)={\frac {(1-q_{1}z^{-1})(1-q_{2}z^{-1})\cdots (1-q_{M}z^{-1})}{(1-p_{1}z^{-1})(1-p_{2}z^{-1})\cdots (1-p_{N}z^{-1})}},}

dóndeqk{\displaystyle q_{k}}es elk{\displaystyle k}el cero ypagk{\displaystyle p_{k}}es elk{\displaystyle k}polo th. Los ceros y los polos suelen ser complejos y cuando se representan en el plano complejo (plano z) se denomina diagrama de polos y ceros .

Además, también pueden existir ceros y polos enz=0{\displaystyle z=0}yz={\displaystyle z=\infty } . 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 siempre es igual.

Al factorizar el denominador, se puede utilizar la descomposición en fracciones parciales , que luego se puede transformar de nuevo al dominio del tiempo . De esta forma, se obtendrían la respuesta impulsional y la ecuación de diferencias lineal con coeficientes constantes del sistema.

Respuesta de salida

Si tal sistemaH(z){\displaystyle H(z)}es impulsado por una señalincógnita(z){\displaystyle X(z)}entonces la salida esY(z)=H(z)incógnita(z){\displaystyle Y(z)=H(z)X(z)} . Al realizar la descomposición en fracciones parciales enY(z){\displaystyle Y(z)}y luego aplicando la transformada Z inversa a la saliday[norte]{\displaystyle y[n]}se puede encontrar. En la práctica, a menudo es útil descomponer fraccionariamenteY(z)z{\displaystyle \textstyle {\frac {Y(z)}{z}}}antes de multiplicar esa cantidad porz{\displaystyle z}para generar una forma deY(z){\displaystyle Y(z)}que tiene términos con transformadas Z inversas fácilmente calculables.

Véase también

Referencias

  1. Mandal, Jyotsna Kumar (2020). " Codificación reversible basada en la transformada Z ". Esteganografía reversible y autenticación mediante codificación por transformada . Estudios en inteligencia computacional. Vol.  901. Singapur: Springer Singapore. pp. 157–195 . doi : 10.1007/978-981-15-4397-5_7 . ISBN  978-981-15-4396-8ISSN 1860-949X . S2CID 226413693. Z es una variable compleja. La transformada Z convierte la señal del dominio espacial discreto en una representación del dominio de la frecuencia compleja. La transformada Z se deriva de la transformada de Laplace .  
  2. Lynn, Paul A. (1986). «La transformada de Laplace y la transformada Z ». Señales y sistemas electrónicos . 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 práctica que la transformada de Fourier de tiempo discreto.
  3. Jury, Eliahu Ibrahim (1964). Teoría y aplicación del método de la transformada Z. Nueva York: John Wiley & Sons. págs. XIII, 330 págs. 
  4. Palani, S. (26 de agosto de 2021). « Análisis de la transformada Z de señales y sistemas de tiempo discreto». Señales y sistemas . Cham: Springer International Publishing. pp. 921–1055 . doi : 10.1007/978-3-030-75742-7_9 . ISBN  978-3-030-75741-0S2CID 238692483. La transformada Z es la contraparte discreta de la transformada de Laplace. La transformada Z convierte ecuaciones en diferencias de sistemas de tiempo discreto en ecuaciones algebraicas, lo que simplifica el análisis de dichos sistemas. La transformada de Laplace y la transformada Z son similares, con la salvedad de que la transformada de Laplace trabaja con señales y sistemas de tiempo continuo. 
  5. ^ ER Kanasewich (1981). Análisis de secuencia temporal en geofísica . Universidad de Alberta. págs.186 , 249. ISBN  978-0-88864-074-1.
  6. ER Kanasewich (1981). Análisis de secuencias temporales en geofísica (3.ª ed.). Universidad de Alberta. págs. 185–186 . ISBN   978-0-88864-074-1.
  7. Ragazzini, JR; Zadeh, LA (1952). "El análisis de sistemas de datos muestreados". Transactions of the American Institute of Electrical Engineers, Part II: Applications and Industry . 71 (5): 225– 234. doi : 10.1109/TAI.1952.6371274 . S2CID 51674188 . 
  8. Cornelius T. Leondes (1996). Implementación de sistemas de control digital y técnicas computacionales . Academic Press. pág. 123. ISBN  978-0-12-012779-5.
  9. Eliahu Ibrahim Jury (1958). Sistemas de control de datos muestreados . John Wiley & Sons.
  10. Eliahu Ibrahim Jury (1973). Teoría y aplicación del método de la transformada Z. Krieger Pub Co. ISBN 0-88275-122-0.
  11. Eliahu Ibrahim Jury (1964). Teoría y aplicación del método de la transformada Z. John Wiley & Sons. pág. 1. 
  12. 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 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, salvo por una sumatoria desde n = 0 hasta el infinito. El uso principal de la transformada unilateral es para secuencias causales, en cuyo caso ambas transformadas son iguales. Por lo tanto, no haremos esta distinción y nos referiremos a ella simplemente como la transformada Z de x ( n ).
  13. Proakis, John; Manolakis, Dimitris. Procesamiento digital de señales: principios, algoritmos y aplicaciones (3.ª ed.). PRENTICE-HALL INTERNATIONAL, INC. 
  14. 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.
  15. 1 2 3 A. R. 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 . 

Lecturas adicionales

  • Refaat El Attar, Apuntes de clase sobre la transformada Z , Lulu Press, Morrisville, Carolina del Norte, 2005. ISBN 1-4116-1979-X.
  • Ogata, Katsuhiko, Sistemas de control de tiempo discreto, 2.ª ed. , Prentice-Hall Inc., 1995, 1987. ISBN 0-13-034281-5.
  • Alan V. Oppenheim y Ronald W. Schafer (1999). Procesamiento de señales en tiempo discreto, 2.ª edición, Prentice Hall Signal Processing Series. ISBN 0-13-754920-2.
  • "Transformada Z" . Enciclopedia de Matemáticas . EMS Press . 2001 [1994].
  • Merrikh-Bayat, Farshad (2014). "Dos métodos para la inversión numérica de la transformada Z". arXiv : 1409.1727 [ math.NA ].
  • Tabla de transformadas Z de algunas transformadas de Laplace comunes
  • Entrada de Mathworld sobre la transformada Z. Archivada el 30/01/2013 en la Wayback Machine.
  • Hilos de transformación Z en Comp.DSP Archivados el 15/06/2012 en Wayback Machine
  • Un gráfico de la relación entre el plano s de la transformada de Laplace y el plano Z de la transformada Z.
  • Explicación en vídeo de la transformada Z para ingenieros.
  • ¿Qué es la transformada Z?