1.

ECUACIONES EN DIFERENCIA Y TRANSFORMADA Z

1.1 ECUACIONES EN DIFERENCIA
1.2 DEFINICIN DE LA TRANSF_Z
1.3 TRANSF_ Z DE FUNCIONES ELEMENTALES
1.4 PROPIEDADES Y TEOREMAS
1.5 TRANSFORMADA Z INVERSA
1.6 SOLUCIN DE ECUACIONES EN DIFERENCIA

1.1 ECUACIONES EN DIFERENCIA

Una ecuacin en diferencias es una ecuacin que relaciona varias secuencias con ellas mismas
desplazadas (o sea, una ecuacin de recurrencia).

a 0 y k a1 y k 1 a 2 y k 2 a n y k n b0 u k b1 u k 1 b2 u k 2 bm u k m
n m

a
i 0
i y k i b j u k j
j 0

1.1.1 DIFERENCIA FINITA

Sea y f [k ] una funcin para k Z, se llama primera diferencia, o diferencia finita de primer
grado de y f [k ] a la expresin dada por
f [k ] f [ k 1] f [k ] , grficamente es

De esta manera vemos que f [ k ] corresponde al incremento que sufre y f [k ] cuando la

variable k se incrementa en una unidad.
CONTROL DIGITAL 2

Ejemplo: Si y [k ] 1 2k 3k 2 ; hallar la diferencia de primer orden.

Solucin:
Yk f [k ] 1 2k 3k 2
Yk 1 f [k 1] 1 2(k 1) 3(k 1) 2 1 2k 2 3(k 2 2k 1) 6 8k 3k 2
f [k ] f [k 1] f [k ] (6 8k 3k 2 ) (1 2k 3k 2 ) 5 6k

Como se observa, a partir del resultado del ejemplo anterior, f [ k ] es otra funcin que
depende de la misma variable k, por tanto se puede hablar de la segunda diferencia de f [k ]
(o primera diferencia de f [ k ] ). Esta se escribe:
2 f [k ] f [ k 1] f [k ]
y as sucesivamente para diferencias de orden mayor.

EJERCICIOS
Hallar la primera y segunda diferencia de cada una de las siguientes funciones.

1. f [k ] k 2 2k 1
2. g [k ] 1 / k
3. h [k ] ln (k )

1.1.2. ECUACIONES DE DIFERENCIA FINITA DE PRIMER ORDEN

Una ecuacin que relacione los valores de una funcin y f [k ] con una o varias de sus
diferencias finitas, se llama ecuacin de diferencia finita en f [k ] .
En adelante se llamar simplemente ecuaciones de diferencia. El orden de estas ecuaciones lo
determina la diferencia de mayor grado que se encuentre en la ecuacin.

Ejemplos:
1. 3Yk 1 Yk 0

2. Yk 1 4Yk k
3. Yk 1 Yk 3 (2)
k

1.1.3. MODELAJE DE SISTEMAS DE TIEMPO DISCRETO MEDIANTE ECUACIONES EN

DIFERENCIAS

Como los sistemas de tiempo discreto que se usan aqu no son tomados del mundo fsico real
(la naturaleza) sino que son una invencin del hombre, normalmente no parten de un modelo
fsico, como ocurra en los sistemas de tiempo continuo, sino que directamente se representan
mediante algn tipo de modelo matemtico.

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 3

Lo ms normal es que los sistemas de tiempo discreto se representen grficamente mediante

diagramas de bloques, donde los elementos fundamentales son: multiplicacin por una
constante, sumador y retardo. As, un sistema de tiempo discreto puede representarse
mediante un diagrama de bloques con estos elementos bsicos y de l se puede obtener la
ecuacin en diferencias correspondiente. Se asume que u k es la entrada, yk es la salida y n es el
orden del sistema de tiempo discreto.

Ejemplo: Obtener un modelo matemtico para los voltajes de nodo V k (k= 1, 2, 3, , N) en la

red de resistencias R-2R.

Solucin: Planteando la ecuacin de suma de corrientes en l nodo k-1 (voltajes de nodo) se

obtiene:
1 1 1 1 1 5 1 1
R 2 R R V k 1 R V k 2 R V k 0 V k 1 V k 2 V k 0
2R R R
O mejor
2V k 5V k 1 2V k 2 0 para k=2, 3, 4, , N
Que es la ecuacin en diferencias que modela a los voltajes de nodo para la red R-2R. En el
nodo N se tiene la restriccin impuesta por la ecuacin
1
V N V N 1
2
Ejemplo: Obtener la ecuacin en diferencias que modela al sistema de tiempo discreto de la
siguiente figura:

Solucin:
Del diagrama de bloques se obtiene: y k u k ay k 1
Y reorganizando queda la ecuacin de diferencias as: y k ay k 1 u k

Ejemplo: Obtener la ecuacin en diferencias que permite modelar al sistema de tiempo

discreto representado por el diagrama de bloques de la siguiente figura:

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 4

Solucin:
1 2
Del diagrama de bloques se tiene que: yk y k 2 y k 3 u k 2u k 2
4 8
1 2
quedando la ecuacin de diferencias as: y k y k 2 y k 3 u k 2u k 2
4 8

1.1.4. SOLUCIN DE ECUACIONES DE DIFERENCIA FINITA

Una funcin Y f [ k ] es una solucin de la ecuacin de diferencia finita si est definida para
todo k= 0, 1, 2, 3,... y satisface la ecuacin dada. , .

Ejemplo: Probar que Yk f [ k ] k 1 es una solucin de la ecuacin de diferencia:

2

Yk 1 Yk 2k 1
Solucin:
Yk 1 Yk [(k 1) 2 1] [ k 2 1] 2k 1
o sea que la funcin Y f [ k ] cumple con la ecuacin de diferencia dada en el ejemplo.

Existen dos clases de soluciones para una ecuacin de diferencia finita: la solucin particular y
la solucin general.
Se llama solucin particular a aquella que no tiene constantes arbitrarias y cumple con la
definicin anterior de solucin, como en el caso del ejemplo.
Se llama solucin general a aquella que teniendo una o varias constantes arbitrarias, cumple
con la definicin anterior de solucin.

Ejemplo: Yk f [ k ] k 1 C , donde C es una constante arbitraria, es tambin una solucin

2

de Yk 1 Yk 2k 1 .

EJERCICIOS.

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 5

Dada la funcin Yk f [ k ] averiguar si es la solucin de la ecuacin correspondiente

1. Yk k ; Yk 1 Yk 1
2. Yk 4 (2) ; Yk 1 2Yk 0
k

3. Yk 1 (k 1) * k ; Yk 1 Yk 2k 3
4. Yk 3 ; Yk 1 Yk 0
(1 k ) * k
5. Yk C ; Yk 1 Yk k
2
1
6. Yk ; Yk 1 3Yk 1
1 k
1 3k 1 Yk
7. k
Y ; Yk 1
2 1 Yk
k
1
8. Yk 1 ; 2Yk 1 Yk 3
2
k 1
9. Yk 3(2 1) ; Yk 1 2Yk 1
10. Yk e ; Yk 1 Yk (1 e)Yk
k

1.1.5. ECUACIN DE DIFERENCIA LINEAL DE PRIMER ORDEN

Se llama ecuacin de diferencia lineal de primer orden con coeficientes constantes en Yk a toda
expresin de la forma:
a1 Yk 1 a0 Yk f [k ]
con a0, a1 constantes y f [ k ] una funcin tal que k = 0, 1, 2, 3,

Ejemplos:

1. 5Yk 1 3Yk k 1
2. 2Yk 1 Yk 2
k

La ecuacin en diferencia de orden n es

an Yk n an 1 Yk n 1 a1 Yk 1 a0 Yk U k

Se recalca que estas ecuaciones son aplicables si y slo si el valor de la funcin en un perodo
cualquiera, est relacionado o depende del valor de dicha funcin en el perodo inmediatamente
anterior, adems la variable independiente debe tomar los valores enteros 0, 1, 2, 3,

Hasta aqu siempre se ha considerado el intervalo [k, k+1] y sus valores correspondientes de Y k
y Yk+1 en los extremos de este intervalo. Sin embargo, se puede presentar otra notacin,
considerando el intervalo [k-1, k] , cuyo valor de Yk f [ k ] sern Yk-1 y Yk respectivamente, por
tanto
f [ k ] Yk 1 Yk
y la ecuacin de diferencia lineal de primer orden con coeficientes constantes se escribir:

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 6

a1 Yk a0 Yk 1 f1 [k ]

Donde a0, a1 son constantes y f1 [ k ] es una funcin que depende de k relacionada con f [ k ] .
El manejo de estas ecuaciones es idntico al de las anteriores, simplemente que Y k de la
primera ecuacin corresponde a Y k-1 de la segunda y Yk+1 de la primera ecuacin corresponde a
Yk de la segunda.

Ejemplo: En un instante cualquiera, una seal es igual a las tres cuartas partes de la recibida en
el periodo anterior aumentada en 0.2. Plantear una ecuacin de diferencia que relacione dos
seales consecutivas.

Solucin:
Yk 1 0.75Yk 0.2 tambin: Yk 0.75Yk 1 0.2

A manera de ejemplo, consideremos una interesante poblacin de canbales a la cual se le toma

muestras de la poblacin cada 15 aos, esta crece la mitad en este periodo. Adems el
muestreo ha demostrado que dichos canbales siempre se comen (literalmente) el 80% de la
poblacin que haba hace dos periodos de muestreo, es decir, 30 aos y que adems, canbales
inmigrantes de otros pueblos llegan a esta poblacin.

Este es un tpico problema de sistemas discretos y queremos modelarlo mediante un ecuacin

de diferencia, para ello suponemos un periodo de muestreo de 15 aos y llamaremos a y[k ] la
salida, es decir, el periodo en que queremos conocer la poblacin de canbales, donde k puede
tomar valores de 0, 1, 2, 3, etc., y u[k ] la asumimos como la entrada, es decir los canbales
inmigrantes de otros pueblos, as la poblacin de canbales en cualquier periodo k (recordar
que un periodo son 15 aos) y[k ] , est dada por:

y[k ] 1.5 y[ k 1] 0.8 y[ k 2] u[k ]

Donde y[ k 1] y y[ k 2] representan la poblacin de canbales de hace uno y dos periodos
respectivamente.

Se puede organizar la ecuacin de tal forma que:

y[k ] 1.5 y[k 1] 0.8 y[k 2] u[k ]
O tambin
y[ k 2] 1.5 y[k 1] 0.8 y[ k ] u[ k 2]

Las anteriores ecuaciones de diferencia modelan el sistema discreto propuesto.

1.1.6. SOLUCIN DE LAS ECUACIONES DE DIFERENCIA

Uso de la tabla

Un mtodo para resolver una ecuacin de diferencia, no muy prctico pero interesante para
conocer la dinmica del sistema es la construccin de una tabla en la cual se van registrando las

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 7

salidas y entradas para los respectivos periodos en el tiempo conforme ste avanza. Para
explicar este mtodo, retomemos el ejemplo de la poblacin de canbales y calculemos cual
sera la oblacin en 60 aos, es decir 4 periodos de muestreo, suponiendo que en un instante
inicial la poblacin existente son los 2 canbales que cada periodo llegan, inmigrantes de otros
pueblos que se toman como entrada al sistema. Conociendo la ecuacin de diferencia:

y[k ] 1.5 y[ k 1] 0.8 y[ k 2] u[k ]

Construimos nuestra tabla as:

k y[k ] y[ k 1] y[ k 2] u[k ]
0 2 0 0 2
1 5 2 0 2
2 7.9 5 2 2
3 9.85 7.9 5 2
4 10.45 9.85 7.9 2

De la tabla anterior, se llega a la conclusin de que la poblacin de canbales en 60 aos es de

10.45.

Como se puede observar, este mtodo es muy conveniente para mirar la dinmica del sistema,
pero que tal si la pregunta fuera cual es la poblacin en 1000 aos?. La construccin de la tabla
para resolver esta inquietud es muy tediosa, de ah que mejor se empleen mtodos
matemticos para resolver este tipo de ecuaciones de diferencia

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 8

Caso 1: Si todas las races ri son diferentes, la solucin general es:

n
y[ k ] c i ri , k 0
k

i 1
De las condiciones iniciales se determinan las constantes c i:
n
y[ k ] c i ri , k= 0, 1, 2, ,n -1
k

i 1

y[0] 1 1 1 c1
y[1] r r2 rn c 2
1
y[2]
n 1 n 1 n 1
y[3] r1 r2 rn c n

Si todas las ri son diferentes, el conjunto de ecuaciones resultante de esta matriz siempre se
puede resolver. Es posible que algunas races as como algunas constantes sean complejas
conjugadas, pero para condiciones iniciales reales, la solucin siempre ser real.

Caso 2: en caso de que se presenten races reales repetidas con multiplicidad m, se tiene que
adicionar a la solucin y a la matriz, trminos como
k k k
ri , kri , , k m 1 ri

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 9

Ejemplo:

Resolver la siguiente ecuacin homognea:

y[ k 3] 7 y[ k 2] 16 y[ k 1] 12 y[ k ] 0

Con condiciones iniciales y[0]=1, y[1]=3, y[2]=0.

La ecuacin caracterstica r 3 7 r 2 16r 12 0 tiene las siguientes races:

r1=3, r2 = r3 =2

La solucin general est dada por:

y[k ] c1 3 k c 2 2 k kc 3 2 k y[k ] c1 3 k (c 2 kc 3 )2 k

Las constantes se calculan de:

1 1 1 0 c1
3 3 2 2 c
2
0 9 4 8 c 3

De donde c1=-8, c2=9 y c3=4.5. Obtenindose la solucin:

y[ k ] (8)3 k (9 4.5k ) 2 k

Caso 3: En el caso de tener races conjugadas complejas, se puede probar que para
condiciones iniciales reales, las constantes c i son tambin conjugadas complejas. Esto resultar
en una solucin real.

Por ejemplo si:

ri Re j , c i R0 e j0
ri 1 Re ( j ) , c i 1 R0 e ( j0 )

De la solucin general:
n
y[k ] c ri 1
i i
k
, k 0

Se tiene que:

ci ri ci 1ri k1 R0 R k e j ( k 0 ) R0 R k e j ( k 0 )
k

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 10

k
ci ri ci 1ri k1 2 R0 R k cos(k 0 )

Ejemplo: Supngase la ecuacin de diferencia de segundo orden:

y[k 2] 2ay[k 1] y[k ] 0

Esta ecuacin tiene como ecuacin caracterstica asociada a r 2 2ar 1 0 cuyas races son:
r1 , r2 a a 2 1
Lo que implica una solucin homognea de la forma:
k k
y[ k ] c1 r1 c 2 r2
As:

y[k ] c1 a a 2 1 k

c2 a a 2 1 k

Dependiendo de a, la solucin puede tener diferentes formas:

Si a <1, las races son complejas y significa que:

a 2 1 0 1 a 2 0
Por tanto:
r1, 2 a j 1 a 2

De donde:
R r1 a 2 (1 a 2 ) 1

Las races se pueden escribir como:

r1 Re ( j ) r2 Re ( j )
(1 a 2 )
R=1 tan 1
a

La solucin es y[k ] 2 R0 cos(k 0 ) donde R 0 , 0 depende de las condiciones iniciales.
Si a=1 se tiene una raz que se repite: r1 r2 1 , de donde la solucin es:
y[k ] (c1 c 2 k )1k
Si a= -1 las races son: r1 r2 1 , de donde la solucin es:
y[k ] (c1 c 2 k )(1) k

Si a>1 entonces tenemos races reales y se puede aplicar el caso 1.

La ecuacin no homognea

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 11

La solucin de la ecuacin general se obtiene mediante la adicin de una solucin particular a la

solucin de la ecuacin homognea, dicha solucin particular no necesariamente satisface las
condiciones iniciales. Una solucin particular puede encontrarse con el mtodo de variacin de
parmetros o el mtodo de parmetros indeterminados. Este ltimo mtodo propone una
solucin particular semejante al lado derecho de la ecuacin pero con un parmetro
indeterminado .
Algunas de las soluciones propuestas dependiendo del lado derecho de la ecuacin se muestran
en el siguiente cuadro:

Lado derecho de la ecuacin Solucin propuesta

en un instante k
k kn
ak k nak
k n a k cos k 1k 0

Polinomio p(k) de
grado m
n

k
i 0
i
i
ak

n n
i k i a k i k i a k cos(k 0 )
i 0 i 0
1[ k ] 1[ k ]

Si a e j son races de multiplicidad m de la ecuacin caracterstica, entonces la solucin

propuesta debe ser multiplicada por k m .

Ejemplo: Resolver la siguiente ecuacin:

y[k 3] 7 y[k 2] 16 y[ k 1] 12 y[ k ] (1) k

Con condiciones iniciales y[0]=1, y[1]=3, y[2]=0.

El lado derecho de la ecuacin es una serie de potencias y por tanto seleccionamos una solucin
particular de la misma serie de potencias. As, se puede seleccionar:
y[ k ] ( 1) k
Sustituyendo:
(1) k 3 7 ( 1) k 2 16 (1) k 1 12 (1) k ( 1) k

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 12

Lo que resulta en:

(1) k (1) 3 7 (1) k (1) 2 16 (1) k (1) 1 12 (1) k (1) k
De donde
36 (1) k ( 1) k
1
Por tanto .
36
Teniendo en cuenta que la solucin a la ecuacin homognea ya fue calculada en el ejemplo
anterior, la solucin general es:
1 1
y[k ] c1 3 k (c 2 kc 3 )2 k (1) k y[k ] (1) k c1 3 k (c 2 kc 3 )2 k
36 36

Ahora calculamos las constantes ci.

37 / 36 1 1 0 c1
107 / 36 3 2 2 c
2
1 / 36 9 4 8 c3

279 316 156

Resultando c1 , c2 , c3 ,
36 36 36

Retomando ahora el problema de canbales y conociendo su ecuacin de diferencia dada por:

y[ k 2] 1.5 y[k 1] 0.8 y[ k ] u[ k 2]

Podemos encontrar la solucin de la anterior ecuacin y de esa forma conocer la poblacin de

canbales en cualquier periodo k sin necesidad de construir una tabla. Para ello, empezamos
por la solucin a la ecuacin homognea y por tanto necesitamos conocer las races de la
ecuacin caracterstica:
r 2 1.5r 0.8 0

Cuyas races son: r1 , r2 0.75 0.4873 j o expresadas de otra forma:

0.487 rad
R (0.75) 2 (0.487) 2 0.894 y tan 1 33.14 0.576
0.75 s
De donde:
r1 (0.894)e j 0.576 r2 (0.894)e j 0.576

Lo que implica una solucin homognea de la forma:

y[k ] c1 (0.894) k e ( j 0.576) k c 2 (0.894) k e ( j 0.576) k

Ahora solucionamos la ecuacin particular y[ k ] :

El lado derecho de la ecuacin es una entrada paso de la forma 2(1[k]) y por tanto la solucin
propuesta tendr la forma (1[k]), as en la ecuacin:

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 13

y[ k 2] 1.5 y[k 1] 0.8 y[ k ] u[ k 2]

Se sustituye la solucin propuesta y por tanto: 1.5 0.8 2
De donde 6.666
Por consiguiente la solucin total ser:
y[k ] c1 (0.894) k e ( j 0.576 ) k c 2 (0.894) k e ( j 0.576 ) k 6.666

Teniendo en cuenta las races complejas la solucin se puede escribir como:

y[ k ] 2 R 0 (0.894) k cos(k 0.576 0 ) 6.666

Ahora con condiciones iniciales y[0]=2, y[1]=5 se tiene que:

2 2 R 0 cos 0 6.666 5 2 R 0 (0.894) cos(0.576 0 ) 6.666

As tenemos que:
2 R 0 cos 4.666 2 R0 cos( 0.576) 1.864
De las anteriores ecuaciones se obtiene que:
2 R 0 cos 4.666 2 R 0 sen 3.763

Y as:
4.666
2 R0 5.994 cos 1 3.82
5.994
De donde la solucin total estar dada por:

y[ k ] 5.994(0.894) k cos(k 0.576 3.82) 6.666

1.2 DEFINICIN DE LA TRANSF_Z

Sea x(t) : Seal en tiempo continuo

x(kT) : Seal en tiempo discreto
La transformada z se define de la siguiente manera:

Z[x(kT)] = X(z) =
k 0
x ( kT ) z k

= x(0) + x(T) z -1 + x(2T) z -2

+ + x(kT) z k
+

1.3 TZ DE FUNCIONES ELEMENTALES

1.3.1 ESCALN UNITARIO

AGUSTIN SOTO JORGE A. POLANA
CONTROL DIGITAL 14

x(kT) =1, kT 0 , y x(kT) = 0, kT 0

Reemplazando en la definicin

X ( z ) x (kT ) * z k ; donde x (kT ) 1
k 0

X ( z ) 1* z k z k s k 1 z 1 z 2 z 3 z k Ec. (1)
k 0 k 0

Si multiplicamos la ecuacin (1) en ambos miembros por z

zs k z 1 z 1 z 2 z 3 z k 1 Ec. (2)

Si restamos la ecuacin (2) - (1)

zsk sk z 1 z 1 z 2 z 3 z k 1 1 z 1 z 2 z 3 z k
s k ( z 1) z z k
z z k
sk
z 1
z 1 z
lim sk lim
k k z 1 z 1

1 z
X ( z) 1

1 z z 1

La anterior relacin se obtiene por Matlab de dos formas:

Utilizando la definicin de la transformada Z

>> syms z k %se crean las variables simblicas z,k

>> f=symsum(z^(-k),k,0,inf);
>> pretty(f) %se arregla la presentacin
z
------
-1 + z

Utilizando el comando ztrans

>> syms z k %se crean las variables simblicas z,k

>> g=1^(k); %se crea la funcin a transformar
>> f=ztrans(g); %se calcula la transformada z
>> pretty(f) %se arregla la presentacin
z
------
-1 + z

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 15

Para graficar la seal escaln unitaria discreta por Matlab, se hace:

% GENERACIN DE ESCALN UNITARIO DISCRETO

x = ones (1,11); % define once valores de 1's
v = [0 10 0 2]; % define valores de ejes
axis (v);
plot (x,'ro') % grafica crculos de color rojo
xlabel ('k') % asigna rotulo al eje x
ylabel ('x(k)') % asigna rotulo al eje y
title (ESCALON UNITARIO DISCRETO)

1.3.2 RAMPA UNITARIA

% GENERACIN DE LA RAMPA UNITARIA DISCRETA

k = 0:10; % define valores de k
x = k; % funcin rampa para x
axis([0 10 0 10]); % define ejes
grid % rejilla para grafica
plot(k, x,'ro') % grafica x en funcin de k
xlabel('k'); % rotulo para eje x
ylabel('x(k)'); % rotulo para eje y
title('RAMPA UNITARIA DISCRETA')

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 16

kT, para k 0 , 1, 2 , 3,
x (kT) seal discreta
0, para k 0

X(z) Z[x(kT)] Z[kT]

k 0
kT * z
k 0
-k

s k 0 z 1 2 z 2 3z 3 kz k Ec. (1)

Si multiplicamos la ecuacin (1) en ambos miembros por z

zs k 1 2 z 1 3 z 2 4 z 3 (k 1) z k Ec. (2)

Si restamos la ecuacin (2) - (1)

zs k s k 1 2 z 1 3z 2 4 z 3 kz k 1 z 1 2 z 2 3z 3 kz k
s k ( z 1) 1 z 1 z 2 z 3 z k
z
Del ejercicio anterior se tiene que 1 z 1 z 2 z 3 z k
z 1
z 1 Tz
sk *
z 1 z 1 ( z 1) 2
Tz 1 Tz
X (z)
1 z
1 2
z 1 2

La anterior relacin se obtiene por Matlab de dos formas:

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 17

Utilizando la definicin de la transformada Z

>> syms z k %se crean las variables simblicas z,k

>> f=symsum(z^(-k)*k,k,0,inf);
>> pretty(f) %se arregla la presentacin
z
---------
2
(-1 + z)

Utilizando el comando ztrans

>> syms z k %se crean las variables simblicas z,k

>> g=k; %se crea la funcin a transformar
>> f=ztrans(g); %se calcula la transformada z
>> pretty(f) %se arregla la presentacin
z
---------
2
(-1 + z)

1.3.3 POTENCIAL: a k (a = constante)

%GENERACION DE LA FUNCION POTENCIAL x(k) = 2 k

k=linspace(0,5,20); % define valores de k
x=2.^ k; % funcin potencial
grid % rejilla para grfica
plot(k, x,'ro') % grfica x en funcin de k
xlabel('k'); % rotulo para eje x
ylabel('x(k)'); % rotulo para eje y
title('POTENCIAL DISCRETA')

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 18

a k , para k 0 , 1, 2 , 3,
x (k) seal discreta
0, para k 0

X ( z ) a k z - k sk 1 az 1 a 2 z 2 a 3 z 3 a k z k Ec. (1)
k 0

Si multiplicamos la ecuacin (1) en ambos miembros por az 1

az 1 s k az 1 a 2 z 2 a 3 z 3 a k z k a k 1 z ( k 1) Ec. (2)

Si restamos la ecuacin (2) - (1)

az 1 sk sk az 1 a 2 z 2 a 3 z 3 a k z k a k 1 z ( k 1) 1 az 1 a 2 z 2 a 3 z 3 a k z k

1 a k 1 z ( k 1)
sk (az 1 1) 1 a k 1 z ( k 1) sk
az 1 1
1 a k a z 1 z k 1 a a z 1 1 z 1 a a z 1 (0)
lim sk lim
k k az 1 1 az 1 1 az 1 1
1 (0) 1 1 z z
lim sk 1
1 a
k az 1 az 1 z 1 a z z a
1 z
X ( z)
1 az 1
za

La anterior relacin se obtiene por Matlab de dos formas:

Utilizando la definicin de la transformada Z

>> syms z k a %se crean las variables simblicas z, k, a

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 19

>> f=symsum(z^(-k)*a^(k),k,0,inf);
>> pretty(f) %se arregla la presentacin
z
- -----
a - z

Utilizando el comando ztrans

>> syms z k a %se crean las variables simblicas z, k, a

>> g=a^(k); %se crea la funcin a transformar
>> f=ztrans(g); %se calcula la transformada z
>> pretty(f) %se arregla la presentacin
z
------------
a (-1 + z/a)

1.3.4 EXPONENCIAL: e -akT (a = constante)

%GENERACION DE LA FUNCION EXPONENCIAL x(k) = e -2k

k = linspace (1,5,20); % define valores de k con
% espaciamiento lineal
x = exp(-2* k); % funcin exponencial
grid % rejilla para grfica
plot(k, x,'bo') % grfica x en funcin de k
xlabel('k'); % rotulo para eje x
ylabel('x(k)'); % rotulo para eje y
title('EXPONENCIAL DISCRETA')

AGUSTIN SOTO JORGE A. POLANA

 CONTROL DIGITAL 20

e akT , para k 0 , 1, 2 , 3,
x (k) seal discreta
0, para k 0

X ( z ) e akT z - k
k 0

sk 1 e aT z 1 e 2 aT z 2 e 3aT z 3 e akT z k Ec. (1)

Si multiplicamos la ecuacin (1) en ambos miembros por e aT z 1

e aT z 1 sk e aT z 1 e 2 aT z 2 e 3aT z 3 e 2 akT z ( k 1) Ec. (2)

Si restamos la ecuacin (2) - (1)

aT
e z s k s k e aT z 1 e 2 aT z 2 e 3aT z 3 e 2 akT z ( k 1) 1 e aT z 1 e 2 aT z 2 e 3aT z 3 e akT z k
1

e akT z 1 1
sk (e aT z 1 1) e akT z 1 1 sk aT 1
e z 1

e akT z 1 1 e aT 1 e z 1 1 e aT (0) z 1 1
lim sk lim aT 1
k k e aT z 1 1 e aT z 1 1 e z 1
1 1 z
lim s k aT 1
k e z 1 1 e aT z 1 z e aT

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 21

1 z
X ( z) aT 1

1 e z z e aT

La anterior relacin se obtiene por Matlab de dos formas:

Utilizando la definicin de la transformada Z

>> syms z k a T %se crean las variables simblicas z k a T

>> f=symsum(z^(-k)*exp(-a*k*T),k,0,inf);
>> pretty(f) %se arregla la presentacin
z exp(a T)
---------------
-1 + z exp(a T)

Utilizando el comando ztrans

>> syms z k a T %se crean las variables simblicas z k a T

>> g=exp(-a*k*T); %se crea la funcin a transformar
>> f=ztrans(g); %se calcula la transformada z
>> pretty(f) %se arregla la presentacin
z
--------------------------
/ z \
exp(-a T) |-1 + ---------|
\ exp(-a T)/

1.3.5 SENOIDAL : sen(wkT)

%GENERACION DE LA FUNCION SENO:

x(k) = sen(wkT)
k = linspace(1,20); % define valores de k con espaciamiento
lineal

x = sin(k); % funcin exponencial

grid % rejilla para grafica
plot(k, x,'bo') % grafica x en funcin de k
xlabel('k'); % rotulo para eje x
ylabel('x(k) =seno(k)'); % rotulo para eje y
title('SENOIDAL DISCRETA')

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 22

sen(kT ) , para k 0 , 1, 2 , 3,
x (k) seal discreta
0, para k 0

X(z) = Z[x(k)] = Z[sen(wkT)]

Por al ecuacin de Euler:

Sen(wkT) = (1/2j) ( e jwkT - e jwkT), reemplazando,

Z[sen(wkT)] = Z[(1/2j)( e jwkT - e jwkT)],

aplicando la transf_z de la exponencial,

1 1 1 1
X (z)
2 j 1 e jT z 1 2 j 1 e jT z 1

reemplazando las exponenciales :

z 1 sen(T )
X ( z)
1 2 z 1 cos(T ) z 2

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 23

1.4 PROPIEDADES Y TEOREMAS.

1.4.1 MULTIPLICACIN POR UNA CONSTANTE

Z[a x(k)] = a Z[x(k)] = a X(z)

1.4.2 LINEALIDAD
Si x(kT) = a f(kT) + b g(kT), entonces,
X(z) = a F(z) + b G(z)

1.4.3 MULTIPLICACIN POR a k

Si y(kT) = a k x(kT), entonces,

Z[a k y(kT)] = a k x(kT) z k
k 0

k 0
x(kT) (a -1
z) k = X(a -1
z)

1.4.4 TEOREMA DE TRASLACIN

Si y(kT) = e - akT
x(kT), entonces,

Z[e - akT
x(kT)] =
k 0
e -akT
x(kT) z k

k 0
x(kT) (e aT
z) k = X(e aT
z)

1.4.5 TEOREMA DEL CORRIMIENTO

Corrimiento hacia atrs:

Z[x(k-n)T] = z n Z[x(k)] = z n
X(z)

Corrimiento hacia adelante:

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 24

n 1
Z[x(k+n)T] = z n
[ X(z)
k 0
x(kT)*z k ]

= z n
X(z) - z n
x(0) - z n-1
x(1) - z n-2
x(2) - - z x(n-1)

Ejemplo:
Z[x(k+3)T] = z 3
X(z) - z 3 x(0) - z 2 x(1) - z x(2)

1.4.6 SUMA DE FUNCIONES

k
Sea y(k) =
h0
x(h) , para k = 0,1, 2,

y(k) = x(0) + x(1) + x(2) + + x(k-1) + x(k)

y(k-1) = x(0) + x(1) + x(2) + + x(k-1), restando estas dos expresiones,
y(k) - y(k-1) = x(k), sacando Transf._Z,
Y(z) z 1
Y(z) = X(z), entonces despejando Y(z) se tiene que:
1
Y ( z) X ( z)
1 z 1

.4.7 TEOREMA DEL VALOR INICIAL

Si el lmite lim X(z) existe, entonces el valor inicial de x(k) = x(0) es igual a:
x(0) lim X ( z )
z

1.4.8 TEOREMA DEL VALOR FINAL

El valor final de x(k), o sea, cuando k (Si X(z) es estable) , es:
x() lim x( k ) lim (1 z 1 ) X ( z )
k z 1

EJEMPLO 1-1

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 25

Encontrar la transformada Z de una funcin escaln de amplitud 4 y desfase en atraso de 3

periodos de muestreo.
Solucin:
x(kT) = 4*u(kT- 3T), asumiendo T = 1 por simplicidad,
x(k) = 4u(k-3)
Z[4u(k-3)] = 4Z[u(k-3)] = 4 z -3
Z[u(k)]

aplicando teorema de corrimiento en atraso

X(z) = 4 z -3
(1/ (1-z 1
)) = 4 / (z 3 z 2)

EJEMPLO 1-2
Obtener la transformada Z de y(k) = 5 k2
para k = 1, 2, 3, .... e igual a cero para k 0.
Solucin:
Sea x(k) = 5 k, entonces y(k) = x(k 2) = 5 k2

Z[y(k)] = Z[5 k2
] = Z[x(k -2)] = z 2
Z[x(k)] = z -2
Z[5 k ] = z 2 * 1/(1 5 z 1
)
Z[5 k2
] = 1 / (z 2 5 z)

EJEMPLO 1-3
Obtener la transformada Z de y(k) =k e 5k
para k = 1, 2, 3, .... e igual a cero para k 0.

Solucin:
Sea x(k) = k, entonces, X(z) = z 1
/ (1 z ) 2 , y adems, y(k) = e
1 5k
x(k ) ,
Aplicando teorema de traslacin,
Z[y(k)] = Z[k e 5k
] = X(e 5k
z )] , reemplazando , en X(z) se tiene:
Y(z) = (e 5k
z) 1
/ (1 (e 5k
z) 1
) 2
=e - 5k
z 1
/ (1 - e - 5k
z 1) 2

EJEMPLO 1-4
Determinar el valor inicial x(0) de una seal si su transformada Z es igual a :
(1 e 5 k ) z 1
X ( z)
(1 z 1 )(1 e 5 k z 1 )
Aplicando el Teorema de valor inicial,

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 26

(1 e 5 k ) 1
x (0) lim X ( z ) 0
z (1 1 )(1 e 5 k 1 )

EJEMPLO 1-5
Determinar el valor final x() de una seal si su transformada Z es igual a:
1 1
X ( z)
(1 z ) (1 e 5 z 1 )
1

1 z 1

x() lim x(k ) lim (1 z 1 ) X ( z ) lim 1 5 1
1
k z 1 z 1
1 e z

EJEMPLO 1-6
Obtener la transformada Z de la figura dada. Tiempo de muestreo = 1.0

Si x(k) = (1/3)k (rampa de pendiente 1/3)

y(k) = x(k) x(k- 3), entonces,

Y(z) = z[y(k)] = z[(1/3)k] z -3
z[(1/3)k]
1 z 1
1 z *z3 1
1 z 1 (1 z 3 )
Y ( z)
3 (1 z 1 ) 2 3 (1 z 1 ) 2 3 (1 z 1 ) 2

1.5 TRANSFORMADA Z INVERSA

Con la transformada Z inversa se obtiene la seal discreta en los instantes de muestreo x(kT).
Los siguientes son los mtodos ms utilizados para obtener la transformada Z inversa.

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 27

1.5.1 MTODO DE DIVISIN DIRECTA

El mtodo consiste en arreglar la funcin X(z) en trminos de z 1 tanto el numerador como el

denominador, dividir algebraicamente el numerador entre el denominador y su cociente
mediante comparacin con la definicin de X(z) encontrar la seal x(kT).

x(k ) z
k 0
k
x(0) + x(1) z 1
+ x(2) z 1 + x(3) z 3 + + x(k) z k

EJEMPLO 1-7
Obtener la transformada Z inversa x(k) de:
5 z 10 5 z 10
X ( z) 2
( z 0.8)( z 0.2) z z 0.16
multiplicando numerador y denominador por z 2
,
5 z 1 10 z 2
X ( z)
1 z 1 0.16 z 2
dividiendo numerador entre denominador, se tiene,
X(z) = 5 z 1
+ 15 z 2
+ 14.2 z 3
+ 11.8 z 4
+
Comparando con la definicin de X(z), se obtiene que:
X(0) = 0, x(1) = 5, x(2) =15, x(3) = 14.2, x(4) = 11.8,
Si se quiere obtener ms muestras de la seal, es mejor aplicar Matlab de esta forma:
% EJEMPLO 1-6: DE TRANSFORMADA Z INVERSA
x = [1 zeros(1,40)]; % para k = 40 muestras
num = [0 5 10]; % coeficientes del numerador
den = [1 -1 0.16]; % coeficientes del denominador
y = filter(num,den,x) % obtencin de las 40 muestras
k = 0:40;
plot(k,y,'ro',k,y,'-')
xlabel('k')
ylabel('y(k)')

Los 40 resultados obtenidos de x(0) hasta x(39) son:

0 5.0000 15.0000 14.2000 11.8000 9.5280 7.6400
6.1155 4.8931 3.9146 3.1317 2.5054 2.0043 1.6035
1.2828 1.0262 0.8210 0.6568 0.5254 0.4203 0.3363
0.2690 0.2152 0.1722 0.1377 0.1102 0.0882 0.0705

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 28

0.0564 0.0451 0.0361 0.0289 0.0231 0.0185 0.0148

0.0118 0.0095 0.0076 0.0061 0.0048 0.0039
Cuya representacin es la siguiente:

1.5.2 MTODO DE FRACCIONES PARCIALES

Consiste el mtodo en expandir la funcin X(z) en fracciones parciales con el fin de que
queden trminos ms simples y luego encontrar a cada fraccin la transformada Z
inversa.
EJEMPLO 1-8
Obtener la transformada Z Inversa de:
5 z 3 26 z 2 44 z 29
X ( z)
z 3 6 z 2 11 z 6
Para representar esta funcin en fracciones parciales usamos el comando Matlab residue que
encuentra los valores del vector r, del vector p y del trmino independiente k segn la siguiente
expresin:
r1 r r
X ( z) 2 n k
z p1 z p 2 z pn
Usando Matlab:
num = [5 26 44 29];
den = [1 6 11 6];
[r, p, k] = residue(num,den)
% Los resultados son:
r = [-2 -5 3], p = [-3 -2 -1], k = 5,

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 29

por tanto, las fracciones parciales quedan :

2 5 3
X ( z) 5
z 3 z 2 z 1
Multiplicando por z-1 :
2 z 1 5 z 1 3 z 1
X ( z) 5
1 3z 1 1 2 z 1 1 z 1
Si se supone que:
v(k) = ak
Z[a k1
] = Z[v(k -1)] = z 1
Z[v(k)] = z -1
Z[a k ] = z 1 * 1/(1 a z 1
)
entonces ,
x(k) = -2(-3) k1
-5(-2)k 1 +3(-1)k 1 + 5k
donde k es el delta kronecker (delta de Dirac en control continuo) cuya transformada Z es igual
a 1.

1.5.3 MTODO DE LOS RESIDUOS

El mtodo plantea que:
m
x(kT) = i 1
(residuos de X(z) z k - 1 en el polo z = zi )

x(kT) = k1 + k2 + + km

(a) Si X(z) z k 1 tiene un polo simple en z = zi , entonces, el residuo es :

ki lim ( z zi ) * X ( z ) *z k 1
z zi

(b) Si X(z) z k 1 tiene un polo mltiple en z = zi de orden q, entonces,

el residuo es:

ki
1
lim
q 1
( q 1)! z zi z q 1

( z zi ) q * X ( z ) *z k 1
EJEMPLO 1-9

CASO POLO SIMPLE.

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 30

Se halla

tiene tres polos simples en z=2, z=3, z=1

Calculo de

POR MATLAB

num=[0 0 1 0];% coeficientes del numerador hay que colocar los ceros para
hallar Xz1
den=[1 -6 11 -6];% coeficientes del denominador
Xz=tf(num,den,-1); % funcin de transferencia en trminos de z
pole(Xz)% obtencin de polos
% la funcin tiene polos simple en z =2, z=3 z=1
% se debe reconstruir la funcin para aplicar el
% comando limit de Matlab
syms z
Xz1 = sum(num.*[z^3 z^2 z 1])/sum(den.*[z^3 z^2 z 1])
% a) clculo del residuo para k1,k2,k3 plos simples

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 31

syms z k
k1 =limit ((z-2)*Xz1*z ^(k-1),z,2)
k2=limit((z-3)*Xz1*z^(k-1),z,3)
k3=limit((z-1)*Xz1*z^(k-1),z,1)
%X(k)=k1+k2+k3

CASO POLO MULTIPLE.

3z 2 9 z
X ( z)
z 3 1.8 z 2 1.05 z 0.2

ki
1
lim
q 1
( q 1)! z zi z q 1

( z zi ) q * X ( z ) *z k 1

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 32

Ahora resolvamos el mismo ejercicio por MAT-LAB

% EJEMPLO 1-9 : PROGRAMA EN MATLAB
num=[0 3 -9 0]; % coeficientes del numerador
den=[1 -1.8 1.05 -0.2];% coeficientes del denominador
Xz=tf(num,den,-1); % funcin de transferencia en
% trminos de z
pole(Xz) % obtencin de polos

% la funcin tiene un polo simple en z = 0.8 y un polo

% doble en z = 0.5
% se debe reconstruir la funcin para aplicar el
% comando limit de Matlab
syms z
Xz1 = sum(num.*[z^3 z^2 z 1])/sum(den.*[z^3 z^2 z 1]);

% a) clculo del residuo para el polo simple

syms z k
k1 =limit ((z-0.8)*Xz1*z ^(k-1),z,0.8)
% k1 = -220/3*4^ k/(5^ k)

% b) clculo del residuo para el polo doble

h=diff((z-0.5)^2*Xz1*z^(k-1),1) % primera derivada
k2=limit(h,z,0.5)
%K2=1/3*(220+150*k)/(2^k)

x(k) = k1 + k2 = -220/3*4^ k/(5^ k) + 1/3*(220+150*k)/(2^k)

1.5.4 UTILIZANDO MATLAB PARA HALLAR Z INVERSA

Para calcular una transformada Z inversa se utiliza el comando iztrans

>> syms z k %se crean las variables simblicas z,k

>> g=z/(z-1); %se crea la funcin a transformar
>> f=iztrans(g);%se calcula la transformada z inversa
>> pretty(f) %se arregla la presentacin

1
>> syms z k %se crean las variables simblicas z,k

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 33

>> g=z/(z-1)^2; %se crea la funcin a transformar

>> f=iztrans(g); %se calcula la transformada z inversa
>> pretty(f) %se arregla la presentacin
n
>> syms z k a %se crean las variables simblicas z,k
>> g=z/(z-a); %se crea la funcin a transformar
>> f=iztrans(g); %se calcula la transformada z
>> pretty(f) %se arregla la presentacin

1.6 ECUACIONES EN DIFERENCIA

Considrese un sistema discreto LTI (Lineal e invariante en el tiempo) dado por la ecuacin en
diferencias:
x(k) + a1 x(k-1) + a2 x(k-2) + +an x(k-n) =
b0 u(k) + b1 u(k-1) + b2 u(k-2) + +bn u(k-n)
donde u(k) es la entrada al sistema y x(k) es la salida.
La forma de solucionar esta ecuacin en diferencia consiste en calcular la transformada Z,
luego aplicar las condiciones iniciales dadas y por ltimo obtener la transformada Z inversa. Se
debe recordar que:
Z[x(k-n)T] = z n Z[x(k)] = z n
X(z), y,
n 1
Z[x(k+n)T] = z n
[ X(z)
k 0
x(kT)*z k ] =

z n
X(z) - z n
x(0) - z n-1
x(1) - z n-2
x(2) - - z x(n-1)

EJEMPLO 1-10
Resolver la siguiente ecuacin en diferencias.
x( k +2) + 5x(k +1) + 6x(k) = 0
Condiciones iniciales: x(0) = 0, x(1) = 1
Solucin :
a) Aplicando Transf._Z, se tiene:
[z 2
X(z) - z 2 x(0) - z x(1)] + 5[z X(z) - z x(0)] + 6 X(z) = 0
b) Sustituyendo condiciones iniciales:

[z 2
X(z) - z 2 (0) - z (1)] + 5[z X(z) - z (0)] + 6 X(z) = 0
=z 2
X(z) - z + 5z X(z) + 6 X(z) = 0
AGUSTIN SOTO JORGE A. POLANA
CONTROL DIGITAL 34

despejando X(z):
z z z
X ( z) (fracciones parciales)
z 5z 6 z 2 z 3
2

c) Obtener la transf_z inversa :

Sabiendo que, z[ ak ] = 1/ (1 a z 1
), entonces ,
1 1
X ( z) 1

1 2z 1 3 z 1
Por tanto su Transf._Z inversa es:
x(k) = (-2)k (-3)k

EJEMPLO 1-11
Resolver la siguiente ecuacin en diferencias.
x(k) + 5x(k-1) + 6x(k-2) = u(k), donde u(k) es el escaln unitario
Solucin :
a) Aplicando Transf._Z, se tiene:
X(z) + 5 z 1
X(z) + 6 z -2
X(z) = U(z), pero U(z) = 1 / (1 z- 1)
Despejando X(z) :
1 z3
X ( z)
(1 z 1 )(1 5 z 1 6 z 2 ) ( z 1)( z 2 5 z 6)

b) Obtener la transf_z inversa:

% Aplicando Matlab para encontrar las fracciones parciales
Xz = zpk([0 0 0 ], [1 -2 -3], 1, -1);
Xz = tf(Xz)
pole(Xz)
% tiene polos simples en : -3.0000 -2.0000 1.0000
[num,den]=tfdata(Xz,'v')
[r,p,k]=residue(num,den)
% r = -6.7500 2.6667 0.0833
% p = -3.0000 -2.0000 1.0000
% k = 1

z 6.75 z 2.6667 z 0.0833

Con base en lo anterior, X ( z) 1
z3 z2 z 1
Multiplicando por z-1 :

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 35

1 6.75 z 1 1 2.6667 z 1 1 0.0833z 1

X ( z) 1
1 3z 1 1 2 z 1 1 z 1
Sabiendo que, z[ak] = 1/ (1 a z 1
), entonces ,
x(k) = (-3)k + 6.75(-3)k - 1 + (-2)k - 2.6667(-2)k 1 + (1)k 0.0833(1)k 1 + k
= (-3)k + 6.75(-1/3)(-3)k + (-2)k 2.6667(-1/2) (-2)k + 1 - 0.0833 + k
x(k) = -5/4*(-3)k + 7/3*(-2)k + 0.9167 + k

EJEMPLO 1-12
Resolver la siguiente ecuacin en diferencias.
x (k + 2) + 0.5x (k + 1) + 0.2x(k) = u(k + 1) + 0.3u(k), (1)
Condiciones iniciales:
x(k) = 0, para k 0 y u(k) = 0, para k 0 y adems,
u(0) = 1.5, u(1) = 0.5, u(2) = -0.5, u(k) =0 para k = 3, 4, 5,
Solucin :
a) Aplicando Transf._Z, se tiene:
[ z 2 X(z) z2 x(0) z x(1) ] + 0.5[ z X(z) z x(0)] + 0.2 X(z)
= [z U(z) z u(0)] + 0.3 U(z) (2)
b) Sustituyendo condiciones iniciales:
Como no se conoce x(1), se debe encontrar reemplazando k por -1 en la ecuacin (1),
x(-1+2) + 0.5x(-1+1) + 0.2x(-1) = u(-1+1) + 0.3u(-1), entonces,
x(1) + 0.5 x(0) + 0.2 x(-1) = u(0) + 0.3 u(-1),
de las condiciones iniciales, x(0)=0, x(-1)= 0, u(0) = 1.5, u(-1) = 0,
reemplazando se tiene que,
x(1) = 1.5
encontrado x(1) se sustituyen condiciones iniciales en (2),
[ z 2 X(z) z (1.5) ] + 0.5[ z X(z) ] + 0.2 X(z) = [z U(z) - z (1.5)] + 0.3 U(z) (2)
Despejando X(z):
z 0 .3
X ( z) U ( z)
z 0 .5 z 0 .2
2

U (z) u (k ) z
k 0
k
u (0) u (1) z 1 u ( 2) z 2 u (3) z 3

k=0
= 1.5 + 0.5 z 1 - 0.5 z 2 , reemplazando,

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 36

( z 0.3)(1.5 z 2 0.5 z 0.5) 1.5 z 3 0.95 z 2 0.35 z 0.15

X ( z)
( z 2 0.5 z 0.2) z 2 z 4 0.5 z 3 0.2 z 2
Para obtener la seal discreta por Matlab:
% OBTENER DIEZ VALORES DE LA SEAL
num = [0 1.5 0.95 -0.35 -0.15];
den = [1 0.5 0.2 0 0];
u = [1 zeros(1,10)];
xk = filter(num,den,u)
k = 0:10;
stem(xk, k) % grafica la seal muestreada

Valores x(k) dados por Matlab:

x(0)= 0, x(1)= 1.5000, x(2)= 0.2000, x(3)=-0.7500, x(4)= 0.1850, x(5)= 0.0575
x(6)= -0.0658, x(7)= 0.0214, x(8)= 0.0025, x(9)= -0.0055, x(10)= 0.0023

Ejemplo 1-13 canibales

Considerando la transformada Z de una secuencia desplazada en el tiempo, la TZ de

(1) ser,

De esta manera, se tiene que,

De donde,

Ahora bien, definiendo la transformada Z de la funcin de excitacin (entrada),

Considerando (3) en (2),

Ahora bien, para que la funcin racional (4) sea una fraccin propia, y se puede realizar
su descomposicin en fracciones parciales, se transformara (4) como sigue:

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 37

Donde y sern numero complejos conjugados, considerando el

orden del polinomio de z.
Realizando la descomposicin a travs de MatLab se tiene,

>> num = [2 0 0 ];
>> den = [1 -2.5 2.3 -0.8];

>> [r ,p, k] = residue(num,den);

r =
6.6667
-2.3333 - 1.8810i
-2.3333 + 1.8810i

p =

1.0000
0.7500 + 0.4873i
0.7500 - 0.4873i

k =
[]

Donde la estructura [r, p, k] corresponde a los coeficientes r-, los polos de

la fraccin racional p-, y los trminos enteros k-.De esta manera se tiene
que la expresin (5) se transforma en:

Determinando la transformada inversa Z de la expresin (7) se tiene:

M-FILE asociado

function y = tz(x)
num = [2 0 0 ];
den = [1 -2.5 2.3 -0.8];
AGUSTIN SOTO JORGE A. POLANA
CONTROL DIGITAL 38

[r ,p, k] = residue(num,den);

y = (r(1))*(p(1).^x) + (r(2))* (p(2).^x) + (r(3))*(p(3).^x);

end

EJECUCION EN EL WORKSPACE

>> x = 0:1:15;
>> tz(x)

TRANSFORMADA Z MODIFICADA.

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 39

La transformada Z modificada se utiliza cuando la funcin de transferencia de un

sistema, presenta un tiempo muerto o retardo .
Asumiendo que la funcin de transferencia del sistema est dada por:

En donde no contiene tiempo muerto y es el tiempo muerto, el procedimiento

para evaluar la transformada z de esta funcin es el siguiente:
Sea

En donde T es el periodo de muestreo y N es la parte entera del cociente:

Reemplazando en la funcin de transferencia de la planta se tiene:

Aplicando la transformada z

El termino se define como la transformada z modificada de y se

denota por . Entonces:

En donde

EJEMPLO

Hallar la transformada z de la funcin

Asumir que el periodo de muestreo es T=1s.

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 40

SOLUCION

Por tabla transformada de Z modificada, se tiene:

Ahora con N=1 y y utilizando la ecuacin

Se tiene:

Por MATLAB se tiene:

sys = tf(5,[1 6 9],InputDelay,1.3)

gpz=c2d(sys,1,zoh)

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 41

TABLA DE LA TRANSFORMADA Z MODIFICADA

AGUSTIN SOTO JORGE A. POLANA

CONTROL DIGITAL 42

EJERCICIOS CAPTULO I

1) Clasifique los sistemas definidos por cada una de las relaciones entrada salida

especificaciones.

a) y k x k 1 2 x k d) y k x k sen( x k k )

b) y k kx k 1 e) y k y k 1 y k 2 x k 2 x k 2

c) y k x k

2) Resuelva la siguiente ecuacin es diferencias: 2 x(k ) 2 x(k 1) x(k 2) u (k )

1, k 0,1,2
Donde x(k)=0 para k<0 y u ( k )
0, k 0

3) Considere la ecuacin en diferencias

x( k 2) 1.3679x( k 1) 0.3679 x( k ) 0.3679u ( k 1) 0.2642u ( k )

Donde x(k) es la salida y x(k)=0 para k 0 y donde u(k) es la entrada y est dada por

u(k)=0, k 0 u(k)=0, k=3, 4, 5,

u(0)=1 u(1)=0.2142 u(2)=-0.2142

determine la salida x(k).

4) Obtenga la transformada z de la curva x(t) que se muestra en la figura

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 43

z 1
5) Dada la transformada z X ( z )
(1 z 1 )(1 1.3 z 1 0.4 z 2 )

Determine los valores inicial y final de x(k). Tambin encuentre x(k), la transformada z inversa
de X(z), en una forma cerrada

z 3
6) Obtenga la transformada z inversa de X ( z ) en una forma cerrada.
(1 z 1 )(1 0.2 z 1 )

7) Resuelva la siguiente ecuacin en diferencias: x (k 2) x(k 1) 0.25 x(k ) u (k 2)

Donde x(0)=1 y x(1)=2. La funcin de entrada u(k) est dada por

u(k)=1, k= 0, 1, 2,
Resuelva este problema tanto de manera analtica como por MATLAB.

0.5 z
8) Desarrollar por fracciones parciales y por series de potencias F ( z )
( z 0.5)( z 0.7)

9) Deducir la transformada z para f (t ) sent u (t )

z ( z 1)( z 2)
10) Encuentre f(kT) si F ( z )
( z 0.5)( z 0.7)( z 0.9)

11) Para cada F(z), encuentre f(kT) usando fracciones parciales

z ( z 3)( z 5)
a) F ( z )
( z 0.4)( z 0.6)( z 0.8)

( z 0.2)( z 0.4)
b) F ( z )
( z 0.1)( z 0.5)( z 0.9)

( z 1)( z 0.2)
c) F ( z )
( z 0.5)( z 0.6)

12) Repita el punto anterior usando Matlab.

BIBLIOGRAFIA BASICA
SISTEMAS DE CONTROL EN TIEMPO DISCRETO.
AUTOR: Katsuhito Ogata.
EDITORIAL: Pearson

SISTEMAS DE CONTROL DIGITAL

AUTOR: BENJAMIN C. KUO
EDITORIAL: CECSA
AGUSTIN SOTO JORGE A. POLANA
CONTROL DIGITAL 44

BIBLIOGRAFIA COMPLEMENTARIA

DIGITAL CONTROL using digital signal processing.

AUTOR: FARZAD NEKOOGAR
EDITORIAL: Prentice Hall

SISTEMAS DE CONTROL CONTINUO Y DISCRETO.

AUTOR: JHON DORSEY
EDITORIAL: Mc GRAW Hill.

CONTROL SYSTEM DESIGN.

AUTOR: GRAHAM C. GOODWIN
EDITORIAL: Prentice Hall.

CONTROL EN EL ESPACIO DE ESTADO.

AUTOR: SERGIO DOMINGUEZ Y OTROS
EDITORIAL: Prentice Hall.

CONTROL DE PROCESOS.
AUTOR: ALFRED ROCA CUSIDO.
EDITORIAL: ALFAOMEGA.

SENALES Y SISTEMAS.
AUTOR: ALAN V. OPPENHEIM
EDITORIAL: Pearson

MATLAB.
AUTOR: MAYH WORKS
EDITORIAL: Prentice Hall
MATEMATICAS AVANZADAS PARA CONTROL(APUNTES).
AUTOR: P.H.D JOSE ALDEMAR MUNOZ H.

UNIVERSIDAD SURCOLOMBIANA
CAPTULO I 45

AGUSTIN SOTO JORGE A. POLANA