Transformada Z y Matlab - DSP (2018-03)

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 44

Tema 3: La Transformada Z

La Transformada Z

1. Definición de la Transformada Z
2. Propiedades de la Transformada Z
3. La Transformada Z inversa
4. Sistemas LTI y dominio Z
5. Estructuras para la realización de sistemas
discretos
1. Definición de la Transformada Z.
La Transformada Z directa.
La transformada Z de una señal discreta x(n) está definida
como una serie de potencias

X ( z )  Z x(n)   x ( n ) z n

n  

Donde z es una variable compleja z=e-jw.


La transformada es llamada directa por transformar una
señal del dominio del tiempo x(n) al plano complejo X(z).
El proceso inverso es llamado transformada inversa Z.
Al ser la transformada Z una serie infinita de potencias,
existe solo para valores de z donde la serie converge.
La región de convergencia (ROC) de X(z) es el conjunto de
valores de z para el cual X(z) alcanza valores finitos.

Ejemplos:
x1(n) = {1, 2, 5, 7, 0, 1} X1(z) = 1 + 2z-1 + 5z-2 + 7z-3 + z-5
x2(n) = {1, 2, 5, 7, 0, 1} X2(z) = z2 + 2z + 5 +7z-1 + z-3
x3(n) = {0, 0, 1, 2, 5, 7, 0, 1} X3(z) = z-2 + 2z-3 + 5z-4 + 7z-5 + z-7
x4(n) = {2, 4, 5, 7, 0, 1} X4(z) = 2z2 + 4z +5 +7z-1+z-3
x5(n) = δ(n) X5(z) = 1
x6(n) = δ(n - k), k > 0 X6(z) = z-k, k > 0
x7(n) = δ(n + k), k > 0 X7(z) = zk, k > 0

¿Cuál es la ROC en cada caso?


Determina la transformada Z de la señal x(n) = ½n u(n).

En Matlab ans =

syms n z z/(z - 1/2)


f=1/2^n;
ztrans(f, n, z)
Determina la transformada Z de la señal x(n) = αn u(n).

syms a n x
f = a^n;
ztrans(f, x)
Determina la transformada Z de la señal x(n) = - bn u(-n-1).
Determina la transformada Z de la señal x(n) = αn u(n) - bn u(-n-1).
Una señal discreta x(n) queda unívocamente determinada
por su transformada z, X(z), y la región de convergencia de
X(z).

La ROC de una señal anticausal es el interior de una


circunferencia de radio r1 mientras que la ROC de una señal
causal es el exterior de un círculo de radio r2.

La ROC para una señal que se extiende hasta el infinito por


los dos lados es un anillo (región anular) en el plano z.

Transformada Z unilateral: X ( z )   x ( n) z
 n

n 0
Señales de duración finita
Plano z
Causal excepto z = 0

Plano z
Anticausal excepto z = ∞

Bilateral Plano z
excepto z = ∞ y z = 0

Señales de duración infinita

r2
Causal |z| > r2

Anticausal
… r1 |z| < r1

r1 r2
Bilateral
… … r2 < |z| < r1
2. Propiedades de la Transformada Z.
Linealidad.
Si
x1 (n) 
z
X1 ( z) y x2 (n) 
z
X 2 ( z)
Entonces
x(n)  a1 x1 (n)  a2 x2 (n) 
 X ( z)  a1 X1 ( z)  a2 X 2 ( z )
z

Determinar la transformada Z de:

%%
Determina la transformada Z de syms w n z
las señales x(n) = (cos ωn )u(n) y f = sin(w*n);
x(n) = (sen ωn)u(n). ztrans(f, z)
Desplazamiento en el tiempo.
Si
x(n) 
 X ( z )
z

Entonces
k
x( n  k ) 
 z X ( z )
z

La ROC de z-kX(z) es la misma que la de X(z) salvo para


z = 0 si k > 0 y z = ∞ si k < 0.
Escalado en el dominio z.
Si
x(n) 
z
X ( z ) ROC : r1  z  r2
Entonces
1
a x(n) 
n
 X (a z) ROC : a r1  z  a r2
z

Para cualquier constante a real o compleja.

Libro Ed.3; Ej: 3.2.1 y 3.2.2


Inversión temporal.
Si x(n) 
z
X ( z ) ROC : r1  z  r2
Entonces
1
1 1
x (  n) 
 X ( z ) ROC :  z 
z

r2 r1

Determina la TZ de la señal x(n) = u(-n)


Derivar en el dominio z.
Si
x(n) 
 X ( z )
z

Entonces
dX ( z )
nx(n) 
  z z

dz
Determina la señal
TZ de x(n)
la señal
si X(z)
x(n) = nanu(n).
= log(1 + az-1) con |z| > |a|.
Convolución de dos secuencias.
Si
x1 (n) 
 X1 ( z ) x2 (n) 
z
 X 2 ( z ) z

Entonces
x(n)  x1 (n)  x2 (n) 
z
X ( z)  X1 ( z) X 2 ( z)
La ROC de X(z) es, cuando menos, la intersección de las de
X1(z) y X2(z)
Determina la convolución de x1(n) = {1, -2, 1} y x2(n) = {1, 1, 1, 1, 1, 1}.
Ejemplo 1

Determinar
Ejemplo 2

Determinar
Transformadas Z comunes
Algunos pares de transformada Z.
Ejemplo 3
Usar las propiedades y la tabla de transformadas comunes de Z y obtener
la transformada Z

Desplazamiento

Derivación

De la tabla
Ejemplo 3
Usando Matlab
Transformada Z Inversa.
Ejemplo 4
Calcular la TZI de:

X(z) tiene dos polos; ROC sin definir, entonces hay 3


posibilidades
Ejemplo 4
Transformada Z Inversa.
Implementación Matlab
La función residuez calcula la parte del residuo y los términos directos
de una función racional en z-1

Los polinomios tienen potencias ascendentes de z-1 .


[R, p, C]=residuez(a,b) calcula los residuos, polos y términos directos
de X(z).
[a, b]=residuez(R, p, C) calcula los residuos, polos y términos directos
de X(z), retorna la función original.
Transformada Z Inversa.
Ejemplo 5 Matlab
Calcular en fracciones parciales la siguiente función:
Transformada Z Inversa.
Ejemplo 5 Matlab
Transformada Z Inversa.
Ejemplo 5 Matlab
Calcular la transformada Z inversa de la siguiente función:

Despues de Matlab
Transformada Z Inversa.
Ejemplo 5 Matlab
Calcular la transformada Z inversa de la siguiente función:

Verificacion en Matlab
Transformada Z Inversa.
Ejemplo 6 Matlab
Calcular la transformada Z inversa de la siguiente función:
Transformada Z Inversa.
Ejemplo 6 Matlab
Calcular la transformada Z inversa de la siguiente función:

Verificacion en Matlab
5. Estructuras para la realización de sistemas
discretos

Representación con diagramas de bloques.


x1(n)

+ x1(n) + x2(n) x(n) z x(n + 1)

x2(n)

x(n) z-1 x(n - 1)


x1(n)

× x1(n) × x2(n)
x(n) a ax(n)
x2(n)
Dibujar el diagrama de bloques para: y(n) = 0.25 y(n-1) + 0.5 x(n) + 0.5 x(n-1)
Determinar la ecuación correspondiente a:

x(n) + + y(n)

z-1 -½ z-1

z-1 2 z-1

1/3 z-1
Forma directa
La función de transferencia de un sistema discreto es:
M

k
b z k
Y ( z)
H ( z)  k 0
N

1   ak z  k
X ( z)
k 1

La ecuación diferencial correspondiente es:


N M
y (n)   ak y(n  k )   bk x(n  k )
k 1 k 0
Su realización es la siguiente:

x(n) y(n)
x(n) y(n)
b0

x(n-1) y(n-1)
b1 -a1

x(n-2) y(n-2)
b2 -a2

x(n-3)
b3 Σ -a3
y(n-3)

x(n-M) y(n-N)
bM -aN
Si el almacenamiento de las entradas y salidas pasadas y
presentes se representa a través del operador de retardo z-1
tenemos:
X(z)
b0 Σ Σ Y(z)

z-1 z-1

b1 -a1

z-1 z-1
H1(z) H2(z)
b2 -a2
Ceros Polos
z-1 z-1

z-1 z-1

bM -aN
Realización Forma Directa I
Intercambiando H1(z) con H2(z) y rearreglando tenemos:

X(z) Σ b0 Σ Y(z)

z-1

-a1
b1
z-1

H2(z) -a2 H1(z)


b2
Polos z-1 Ceros

z-1

-aM bM
z-1
Realización Forma Directa II
-aN
Cascada y paralelo
La realización en cascada se obtiene reconociendo que los
polos o ceros de la función de transferencia H(z) al ser
números reales o conjugados complejos se pueden escribir
de manera factorizada.

i1 i  j 1 j
N1 1 N2  1 1
(1  a z ) (1  b z )(1  b j z )
H ( z )  Kz  M
k 1 k l 1 l
D1 1 D2  1 1
(1  c z ) (1  d z )(1  bl z )
Donde
N1 son números reales de ceros para z = ai
N2 son pares conjugados complejos de ceros z = bj y z = bj*
D1 son números reales de polos para z = ck
D2 son pares conjugados complejos de ceros z = dl y z = dl*
Los polos y ceros reales son empleados típicamente en la
Forma Directa II.
Por ejemplo, 1  ai z 1
H ( z) 
1  ck z 1
Se efectua con la estructura:

X(z) Σ Σ Y(z)

z-1

ck -ai

Estructura de primer orden


Los polos y ceros conjugados complejos se hallan en pares.
Su forma general es:
(1  b j z 1 )(1  bj z 1 ) 1  (b j  bj ) z 1  b j bj z 2
H ( z)  1  1

(1  d l z )(1  d z )l 1  (d l  d l ) z 1  dl d l z 2
Y se realizan en la siguiente estructura:
X(z) Σ Σ Y(z)

z-1

dl+dl* Estructura de
-(bj-bj*)
segundo orden
z-1
-dldl* b j b j*

Como A+A* y AA* son números reales, siendo A un número


complejo, todos los multiplicadores son números reales.
La forma paralela resulta de expander H(z) en fracciones
parciales. Su forma general para polos simples es:
1
M D1
1 D2
1  e z
H ( z )   Ai z 1   Bk 1
  C 1
l
 1
  
l
i 0 k 1 1 ck z l 1 (1 d l z )(1 d l
z )

Donde la primera sumatoria se encarga de los términos de


la expansión en fracciones parciales que resulta si M > N.
La segunda sumatoria se encarga de los polos reales y la
tercera de los pares conjugados complejos.
Si se ignoran los errores de cuantización no existe problema
en la realización del sistema. Pero si son considerados, los
resultados pueden cambiar al moverse los polos o ceros
produciendo inestabilidad.
Realizar la estructura en cascada y paralelo de

Hallar la ecuación en diferencias de:

x(n) + s + s y(n)

0.25
-0.25
Hallar la ecuación en diferencias de:
2

x(n) + s + y(n)

-2
Bibliografía
 Digital Signal Processing: Principles, algorithms and applications
J. G. Proakis & D. G. Manolakis.
Pearson Education Inc. 3a Ed. 1996.

 Introduction to Signals and Systems,


D. K. Lindner
McGraw Hill, 1999.

 Signals and Systems: Continuous and Discrete.


R. E. Ziemer, W. H. Tranter & D. R. Fannin
Prentice Hall, 4a Ed. 1998

 Principles of Signals and Systems


F. J. Taylor
McGraw Hill, 1a Ed. 1994

 Signals and Systems


A. V. Oppenheim
Prentice Hall, 1a Ed. 1993.

 Analog and Digital Communication Systems


M. S. Roden
Prentice Hall, 4a Ed. 1996.

También podría gustarte