PROCESAMIENTO DIGITAL DE SEÑALES
S0x:
Ing. Luis Gutiérrez Magán
[Link]@[Link]
PROCESAMIENTO DE SEÑALES
CONTENIDO:
Transformada Z Inversa 9.1
TZI por expansión en fracciones parciales 9.2
Propiedades de la Transformada Z 9.3
LOGROS DE LA SESIÓN
Al término de la sesión, el estudiante conoce los procedimientos para
obtener la Transformada Z Inversa y las propiedades y aplicaciones de
la Transformada Z
Transformada z inversa
Expresión general: Recuerde que, para 𝑧 =∞𝑟𝑒 −𝑗𝑤 , la Transformada
∞
Z G(z) dada por
𝐺 𝑧 = 𝑔[𝑛]𝑧 −𝑛 = 𝑔[𝑛]𝑟 −𝑛 𝑒 −𝑗𝑤𝑛
𝑛=−∞ 𝑛=−∞
−𝑛
es simplemente la DTFT de la secuencia modificada 𝑔 𝑛 𝑟
En consecuencia, la DTFT inversa viene dada por
1 𝜋
𝑔]𝑛]𝑟 −𝑛 = න 𝐺(𝑟𝑒 𝑗𝑤 )𝑒 𝑗𝑤𝑛 𝑑𝑤
2𝜋 −𝜋
Al hacer un cambio de la variable 𝑧 = 𝑟𝑒 −𝑗𝑤 , la ecuación anterior se puede convertir en una integral de
contorno dada por
1
𝑔𝑛 = ර 𝐺(𝑧)𝑧 𝑛−1 𝑑𝑧
2𝜋𝑗 𝐶 ′
Donde 𝐶 ′ es un contorno de integración en sentido antihorario definido por |z|=r
Transformada z inversa
Pero la integral permanece sin cambios cuando se reemplaza con cualquier contorno C que rodea el
punto z = 0 en la ROC de G (z)
La integral de contorno se puede evaluar usando el teorema del residuo de Cauchy que resulta en:
𝑟𝑒𝑠𝑖𝑑𝑢𝑜 𝑑𝑒 𝐺(𝑧)𝑧 𝑛−1
𝑔 𝑛 =
𝑒𝑛 𝑙𝑜𝑠 𝑝𝑜𝑙𝑜𝑠 𝑑𝑒𝑛𝑡𝑟𝑜 𝑑𝑒 𝐶
La ecuación anterior debe evaluarse en todos los valores de n.
Una transformada z racional G(z) con una transformada inversa causal g[n] tiene una ROC que es
exterior a un círculo
Es más conveniente expresar G(z) en una forma de expansión de fracción parcial y luego determinar
g[n] sumando la transformada inversa de los términos individuales en la expansión.
Transformada z inversa por expansión en fracción parcial
Una G(z) racional puede ser expresada como:
𝑃(𝑧) σ𝑀
𝑖=0 𝑝𝑖 𝑧
−𝑖
𝐺 𝑧 =
𝐷(𝑧) σ𝑁𝑗=0 𝑑𝑗 𝑧
𝑗
Si 𝑀 ≥ 𝑁 entonces G(z) puede ser expresado como:
𝑀−𝑁
𝑃1 (𝑧)
𝐺 𝑧 = 𝜂𝑙 𝑧 −𝑙 +
𝐷(𝑧)
𝑙=0
donde el grado de 𝑃1 (𝑧) es menor que N. La función racional 𝑃1 (𝑧)/D(z) es llamada fracción propia.
Ejemplo: Considerar
2 + 0.8𝑧 −1 + 0.5𝑧 −2 + 0.3𝑧 −3
𝐺 𝑧 =
1 + 0.8𝑧 −1 + 0.2𝑧 −2
Luego de una larga división obtenemos:
−1
5.5 + 2.1𝑧 −1
𝐺 𝑧 = −3.5 + 1.5𝑧 +
1 + 0.8𝑧 −1 + 0.2𝑧 −2
Transformada z inversa por expansión en fracción parcial
Polos Simples: En la mayoría de los casos prácticos, la transformada z racional de interés
G(z) es una fracción propia con polos simples
Sean los polos de G(z) en 𝑧 = 𝜆𝑘 1 ≤ 𝑘 ≤ 𝑁
Entonces, una expansión de fracción parcial de G(z) tiene la forma:
𝑁
𝜌𝑙
𝐺 𝑧 =
1 − 𝜆𝑙 𝑧 −1
𝐿=1
Las constantes en la expansión de fracciones parciales se denominan residuos y están
dadas por:
𝜌𝑙 = 1 − 𝜆𝑙 𝑧 −1 𝐺(𝑧)|𝑧=𝜆𝑙
Cada término de la suma en la expansión de fracciones parciales tiene una
𝑛
ROC dada por
𝑧 > 𝜆𝑙 y, por tanto, tiene una transformada inversa de la forma 𝜌(𝜆) 𝜇[𝑛]
Transformada z inversa por expansión en fracción parcial
Por lo tanto, la transformada inversa g[n] de G(z) viene dada por:
𝑔 𝑛 = 𝜌𝑙 (𝜆𝑙 )𝑛 𝜇[𝑛]
𝑙=1
Nota: El enfoque anterior con una ligera modificación también se puede utilizar para determinar la inversa de una transformada z racional de una
secuencia no causal.
Ejemplo: sea la Transformada Z H(z) de una secuencia causal h[n] será dado por
𝑧(𝑧 + 2) 1 + 2𝑧 −1
𝐻 𝑧 = =
(𝑧 − 0.2)(𝑧 + 0.6) (1 − 0.2𝑧 −1 )(1 + 0.6𝑧 −1 )
Entonces, una expansión de fracción parcial de H(z) tiene la forma:
𝜌1 𝜌2
𝐻 𝑧 = +
1 − 0.2𝑧 −1 1 + 0.6𝑧 −1
Transformada z inversa por expansión en fracción parcial
Ahora
1 + 2𝑧 −1
𝜌1 = 1 − 0.2𝑧 −1 𝐻 𝑧 |𝑧=0.2 = −1
|𝑧=0.2 = 2.75
1 + 0.6𝑧
1 + 2𝑧 −1
𝜌2 = 1 + 0.6𝑧 −1 𝐻 𝑧 |𝑧=0.6 = | = −1.75
1 − 0.2𝑧 −1 𝑧=−0.6
Por lo tanto
2.75 1.75
𝐻 𝑧 = −1 −
1 − 0.2𝑧 1 + 0.6𝑧 −1
Por tanto, la transformada inversa de lo anterior viene dada por:
ℎ 𝑛 = 2.75(0.2)𝑛 𝜇 𝑛 − 1.75 −0.6 𝑛 𝜇[𝑛]
Transformada z inversa por expansión en fracción parcial
Polos múltiples: Si G (z) tiene varios polos, la expansión de la fracción parcial es de forma ligeramente
diferente
Sea el polo en z = ν de multiplicidad L y el resto N
o Los polos L sean simples y en z = λ, 1 ≤ l ≤ N – L
Entonces la expansión de fracción parcial de G(z) es de la forma
𝑀−𝑁 𝑁−𝐿 𝐿
−𝑙
𝜌𝑙 𝛾𝑖
𝐺 𝑧 = 𝜂𝑙 𝑧 + +
1 − 𝜆𝑙 𝑧 −1 (1 − 𝑣𝑧 −1 )𝑖
𝑙=0 𝑙=1 𝑖=1
donde las constantes se calculan usando
1 𝑑 𝐿−𝑖
𝛾𝑖 = 𝐿−𝑖 −1 𝐿−𝑖
1 − 𝑣𝑧 −1 𝐿 𝐺(𝑧) 𝑧=𝑣 1≤𝑖≤𝐿
𝐿 − 𝑖 ! (−𝑣) 𝑑(𝑧 )
Los residuos de 𝜌𝑙 se calculan como antes
Transformada z inversa por expansión en fracción parcial
La Transformada Z G(z) de una secuencia causal {g [n]} puede ser expandido en una serie de potencias en 𝑧 −1
En la expansión de la serie, el coeficiente que multiplica el término 𝑧 −𝑛 es entonces la enésima muestra g[n]
Para una Transformada Z racional expresada como una relación de polinomios en 𝑧 −1, la expansión de la serie
de potencias se puede obtener por división larga
Ejemplo: considere
1 + 2𝑧 −1
𝐻 𝑧 =
1 + 0.4𝑧 −1 − 0.12𝑧 −2
La división larga del numerador por el denominador produce:
𝐻 𝑧 = 1 + 1.6𝑧 −1 − 0.52𝑧 −2 + 0.4𝑧 −3 − 0.2224𝑧 −4 + ⋯
Como resultado ℎ 𝑛 = 1 1.6 − 0.52 0.4 − 0.2224 … 𝑛 ≥ 0
Propiedades de la transformada z
Propiedad Secuencia Transformada – Z ROC
𝑔[𝑛] 𝐺 𝑧 ℛ𝑔
ℎ[𝑛] 𝐻(𝑧) ℛℎ
Conjugación 𝑔∗ [𝑛] 𝐺 ∗ (𝑧 ∗ ) ℛ𝑔
Inversión del tiempo 𝑔[−𝑛] 𝐺 1ൗ𝑧 1ൗ
ℛ𝑔
Linealidad 𝛼𝑔 𝑛 + 𝛽ℎ[𝑛] 𝛼𝐺 𝑧 + 𝛽𝐻(𝑧) Incluye ℛ𝑔 ∩ ℛℎ
Desplazamiento en el 𝑔 𝑛 − 𝑛0 𝑧 −𝑛0 𝐺(𝑧) ℛ𝑔 , excepto posiblemente el
tiempo punto z = 0 o ∞
Multiplicación por una 𝛼 𝑛 𝑔[𝑛] 𝐺(𝑧Τ𝛼) 𝛼 ℛ𝑔
secuencia exponencial
Diferenciación de 𝐺(𝑧) 𝑛𝑔[𝑛] 𝑑𝐺(𝑧) ℛ𝑔 , excepto posiblemente el
−𝑧
𝑑𝑧 punto z = 0 o ∞
Convolución 𝑔[𝑛]⨂ℎ[𝑛] 𝐺 𝑧 𝐻(𝑧) Incluye ℛ𝑔 ∩ ℛℎ
Modulación 𝑔 𝑛 ℎ[𝑛[ 1 𝑧 Incluye ℛ𝑔 ℛℎ
ර 𝐺 𝑣 𝐻( )𝑣 −1 𝑑𝑣
2𝜋 𝑣
∞
Relación de Parseval 1
𝑔 𝑛 ℎ∗ 𝑛 = ර 𝐺 𝑣 𝐻 ∗ (1/𝑣 ∗ )𝑣 −1 𝑑𝑣
2𝜋𝑗
𝑛=−∞
Nota: si ℛ𝑔 denota la región ℛ𝑔− < 𝑧 < ℛ𝑔+ y ℛℎ denota la región ℛ ℎ− < 𝑧 < ℛ ℎ+ , luego 1ൗℛ𝑔 denota la
región 1/ℛ𝑔− < 𝑧 < 1/ℛ𝑔+ y ℛ𝑔 ℛℎ denota la regiónℛ𝑔− ℛℎ− < ℛ𝑔+ ℛℎ+
Propiedades de la transformada z
Ejemplo: Considerar la secuencia de dos lados
𝑣 𝑛 = 𝛼 𝑛 𝜇 𝑛 − 𝛽𝑛 𝜇[−𝑛 − 1]
Dejar 𝑥 𝑛 = 𝛼 𝑛 𝜇 𝑛 y y 𝑛 = 𝛽𝑛 𝜇[−𝑛 − 1] con X(z) y Y(z) que tienen sus
transformadas
1 1
Ahora 𝑋 𝑧 = −1 , 𝑧 > 𝛼 y 𝑌 𝑧 = −1 , 𝑧 < 𝛽
1−𝛼𝑧 1−𝛽𝑧
Usando la propiedad de linealidad, podemos decir
1 1
𝑉 𝑧 =𝑋 𝑧 +𝑌 𝑧 = −1
+
1 − 𝛼𝑧 1 − 𝛽𝑧 −1
La ROC de V(z) esta dada por la superposición de las regiones de 𝑧 > 𝛼 y 𝑧 <
𝛽
Propiedades de la transformada z
Ejemplo: Determinar la Transformada Z y la ROC de la secuencia causal
𝒗 𝒏 = 𝒓𝒏 𝒄𝒐𝒔𝒘𝟎 𝒏 𝝁[𝒏]
Podemos expresar 𝑥 𝑛 = 𝑣 𝑛 + 𝑣 ∗ [𝑛] donde
𝟏 𝒏 𝒋𝒘 𝒏 𝟏
𝒗𝒏 = 𝒓 𝒆 𝟎 𝝁 𝒏 = 𝜶𝒏 𝝁[𝒏]
𝟐 𝟐
La transformada Z de 𝒗 𝒏 esta dada por
1 1 1 1
𝑉 𝑧 = . = . , 𝑧 > 𝛼 =𝑟
2 1 − 𝛼𝑧 −1 2 1 − 𝑟𝑒 𝑗𝑤0 𝑧 −1
Usando la propiedad de la conjugada obtenemos la Transformada Z de 𝑣 ∗ [𝑛] como:
1 1 1 1
𝑉∗ 𝑧∗ = . = . , 𝑧 > 𝛼
2 1 − 𝛼 ∗ 𝑧 −1 2 1 − 𝑟𝑒 −𝑗𝑤0 𝑧 −1
Propiedades de la transformada z
Usando la propiedad de linealidad, obtenemos:
1 1 1
𝑋 𝑧 = 𝑉 𝑧 + 𝑉∗ 𝑧∗ = +
2 1 − 𝑟𝑒 𝑗𝑤0 𝑧 −1 1 − 𝑟𝑒 −𝑗𝑤0 𝑧 −1
o
1 − (𝑟𝑐𝑜𝑠𝑤0 )𝑧 −1
𝑋 𝑧 = , 𝑧 >𝑟
1 − 2𝑟𝑐𝑜𝑠𝑤0 𝑧 −1 + 𝑟 2 𝑧 −2
Ejemplo: Determinar la Transformada Z Y(z) y la ROC de la secuencia
𝑦 𝑛 = 𝑛 + 1 𝛼 𝑛 𝜇[𝑛]
Podemos escribir 𝑦 𝑛 = 𝑛𝑥 𝑛 + 𝑥[𝑛] donde 𝑥 𝑛 = 𝛼 𝑛 𝜇[𝑛]
Propiedades de la transformada z
Ahora, la Transformada Z de X(z) de 𝑥 𝑛 = 𝛼 𝑛 𝜇[𝑛] esta dada por
1
𝑋 𝑧 = , 𝑧 > 𝛼
1 − 𝛼𝑧 −1
Usando la propiedad de Diferenciación, decimos que la Transformada Z de 𝑛𝑥 𝑛 es:
𝑑𝑋 𝑧 𝛼𝑧 −1
−𝑧 = , 𝑧 > 𝛼
𝑑𝑧 1 − 𝛼𝑧 −1
Usando la propiedad de linealidad obtenemos:
1 𝛼𝑧 −1
𝑌 𝑧 = + , 𝑧 > 𝛼
1 − 𝛼𝑧 −1 1 − 𝛼𝑧 −1
Propiedades de la transformada z
Ahora la Transformada Z X(z) de 𝑥 𝑛 = 𝛼 𝑛 𝜇[𝑛] esta dada por:
1
𝑋 𝑧 = −1
, 𝑧 > 𝛼
1 − 𝛼𝑧
Usando la propiedad de la diferenciación, arribamos a la transformada Z de nx[n]
como:
𝑑𝑋 𝑧 𝛼𝑧 −1
−𝑧 = −1
, 𝑧 > 𝛼
𝑑𝑧 1 − 𝛼𝑧
Usando la propiedad de linealidad podemos obtener finalmente
1 𝛼𝑧 −1 1
𝑌 𝑧 = −1
+ −1 2
= −1 2
, 𝑧 > 𝛼
1 − 𝛼𝑧 (1 − 𝛼𝑧 ) (1 − 𝛼𝑧 )
MUCHAS
GRACIAS!!!
por su atención…
Ing. Luis Gutiérrez Magán
[Link]@[Link]