Transformada Z

Contents
1 Introduccin 5
2 Transformada Z 7
2.1 Propiedades de la transformada Z . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2 La transformada Z inversa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2.1 Mtodo de la divisin directa . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.2.2 Mtodo de expansin en fracciones parciales . . . . . . . . . . . . . . . . . 15
2.3 Ecuaciones en diferencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3 Anlisis de sistemas de control en tiempo discreto 21
3.1 Muestreo mediante impulsos y retencin de datos . . . . . . . . . . . . . . . . . . 21
3.1.1 La funcin de transferencia pulso. . . . . . . . . . . . . . . . . . . . . . . 22
3.1.2 Funcin de transferencia pulso de sistemas en lazo cerrado . . . . . . . . . 25
4 Lugar de las races 28
4.0.3 Reglas generales para la construccin de los lugares geomtrico de la raz. 28
5 Diseo de controladores 32
5.1 Anlisis en el dominio continuo . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
5.1.1 Controlador PI (Proporcional-Integral) . . . . . . . . . . . . . . . . . . . . 33
5.2 Anlisis en el dominio discreto . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
5.3 Equivalente discreto de un controlador continuo . . . . . . . . . . . . . . . . . . . 38
5.4 Controladores PID . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2
Dr. Virgilio Vsquez Lpez
Depto. de Mecatrnica y Automatizacin
ITESM-CEM
5.4.1 Sintonizacin de controladores PID utilizando los mtodos de Ziegler-
Nichols . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
6 Modelo en variables de estado 44
7 Sistemas electroneumticos 45
7.1 Mtodo paso a paso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
A 1er regla de Ziegler-Nichols 51
3
ITESM-CEM
Chapter 1
Introduccin
Este curso trata del uso de computadoras digitales para el control de procesos. Primero tratare-
mos el anlisis de sistemas en lazo abierto para establecer los fundamentos tericos necesarios
para el diseo de sistemas de control por computadora. El uso de las computadoras digitales
como elementos de control han crecido rpidamente en los pasados 30 aos. Las principales
razones de este uso extensivo radican en su alto rendimiento, versatilidad y conabilidad. Ms
aun, gracias a las computadoras, se han podido implementar nuevas aplicaciones que requieren
de clculos muy complicados, imposibles de realizar por otros medios. Ejemplo de estas apli-
caciones se pueden encontrar en robtica, anlisis de seales, optimizacin de procesos y en el
rea del control adaptable.
El tipo de seales sobre las cuales esta basado el curso de control continuo es seales con-
tinuas, esto signica que ellas estn denidas sobre un rango continuo de tiempo y pueden
cambiar de valor en cualquier momento. Las seales discretas, por otro lado, estn denidas
nicamente en ciertos instantes de tiempo, y pueden cambiar de valor nicamente en esos in-
stantes de tiempo. Muy frecuentemente las seales discretas son el resultado de un proceso
de muestreo de seales continuas. Los instantes de muestreo estn, generalmente, igualmente
espaciados por un tiempo denominado periodo de muestreo T. Cuando una seal continua es
muestreada, los valores muestreados resultantes forman una seal en tiempo discreto, o simple-
mente una seal discreta. Grcamente este fenmeno se representa en la gura 1.
5
ITESM-CEM
Para que una seal pueda ser administrada en una computadora, esta tiene que ser digitalizada
previamente. Para efectuar esta operacin se hace uso de los dispositivos electrnicos llama-
dos convertidores AD (AnlogoDigital). Tambin existen dispositivos que convierten seales
digitales a continuas conocidos como convertidores DA ( DigitalAnalgica). Dado que estos
dispositivos requieren de un tiempo de conversin, la accin de mandar una seal a una com-
putadora se efecta con un dispositivo de muestreo y retencin junto con un convertidor AD,
(vase la gura 2). Como lo indica la Figura 2, la conversin de la seal consiste en convertir
los valores muestreados x(i) a nmeros binarios, los cuales ya pueden ser suministrados a la
computadora.
6
ITESM-CEM
Chapter 2
Transformada Z
La transformada Z es una herramienta clsica para el anlisis y sntesis de sistemas discretos.
El papel de la transformada Z en sistemas en tiempo discreto, es similar al de la transformada
de Laplace en sistemas en tiempo continuo y se obtiene aplicando al transformada de Laplace
en seales discretas. Las seales en tiempo discreto surgen si el sistema involucra la operacin
de muestreo de seales en tiempo continuo, por lo que la trasnformada Z est relacionada
inherentemente a un proceso de muestreo.
La transformada Z (transformada Z unilateral)de una seal arbitraria es:
X(Z) = z {x(kT)} = z {x(t)} =

X
k=0
x(kT)z
k
(2.1)
= x(0) + x(T)z
1
+ x(2T)z
2
+ x(3T)z
3
+ ...
Ejemplo 1 Obtenga la transformada Z de la funcin escaln unitario
u(t) =
1 t 0
0 t < 0
Note que la seal x(kT) = 1 k 0, por lo tanto utilizando la denicin de tarnsformada Z
(2.1), se tiene que
z {u(t)} =

X
k=0
z
k
= 1 + z
1
+ z
2
+ z
3
+ .... (2.2)
7
ITESM-CEM
Por otro lado, utilizando un resultado de series
1 + a
1
+ a
2
+ a
3
+ ... =
1
1 a
1
(2.3)
La expresin (2.2) se puede reescribir como
z {u(t)} =

X
k=0
z
k
= 1 + z
1
+ z
2
+ z
3
+ .... =
1
1 z
1
Ejemplo 2 Obtenga la transformada Z de la funcin exponencial
x(t) =
e
at
t 0
0 t < 0
De la denicin de transformada Z
z {x(t)} =

X
k=0
e
akT
z
k
= 1 +
e
aT
z
1
+
e
aT
z
2
+
e
aT
z
3
+ ....
Utilizando nuevamente el resultado (2.3), se tiene que
z {x(t)} =
1
1 (e
aT
z)
1
=
z
z e
aT
Ejemplo 3 Obtenga la transformada Z de la siguiente expresin
X(s) =
1
s (s + a)
aplicando fracciones parciales
X(s) =
c
1
s
+
c
2
s + a
donde
c
1
= s
1
s (s + a)
s=0
=
1
a
, c
1
= (s + a)
1
s (s + a)
s=a
=
1
a
8
ITESM-CEM
por lo tanto
x(t) =
1
a

1
a
e
at
x(kT) =
1
a

1
a
e
akT
X(Z) =
1
a
z
z1

1
a
z
ze
aT
2.1 Propiedades de la transformada Z
El uso de la transformada Z se puede facilitar sustancialmente al utilizar algunas de las
propiedades de la transformada Z, las cuales se pueden obtener directamente de la denicin.
En las siguientes propiedades listadas a continuacin, se asume que
F
i
(z) = z {f
i
(t)} = z {f
i
(kT)}
Propiedad de la linealidad
Una funcin f(x) es lineal si f (x
1
+ x
2
) = f (x
1
) + f (x
2
). Aplicando este resultado a la
denicin de la transformada Z, se obtiene inmediatamente que
Z {f
1
(kT) + f
2
(kT)} =

X
k=0
[f
1
(kT) + f
2
(kT)] z
k
= Z {f
1
(kT)} + Z {f
2
(kT)}
= F
1
(z) + F
2
(z)
Entonces, la transformada Z es una funcin lineal. La propiedad de la linealidad hace posible
que se pueda aplicar la tcnica de fracciones parciales.
Teorema de la traslacin real
Z {f(t nT)} = z
n
F(z) (2.4)
y
Z {f(t + nT)} = z
n
"
F(z)
n1
X
k=0
f(kT)z
k
#
(2.5)
9
ITESM-CEM
A partir de esta ltima expresin se tiene que
Z {f(t + T)} = z [F(z) f(0)] = zF(z) zf(0)
Z {f(t + 2T)} = z
2
F(z) f(0) f(T)z

1
= z
2
F(z) z
2
f(0) zf(T)
Z {f(t + 3T)} = z
3
F(z) f(0) f(T)z

1
f(2T)z
2
= z
3
F(z) z
3
f(0) z
2
f(T) zf(2T)
Esta propiedad es una herramienta indispensable en la solucin de ecuaciones en diferencias.
De las propiedades anteriores note que la multiplicacin de F(z) por z tiene el efecto de avanzar
la seal f(kT) un perodo de muestreo y que la multiplicacin de F(z) por z
1
tiene el efecto
de retrasar la seal f(kT) un perodo de muestreo.
Ejemplo 4 Determine la transformada Z de la siguiente funcin
Figure 2-1: Escalon desplazado en el tiempo
Solucin 5 La funcin que describe la grca anterior es
f(t) = u(t T)
Aplicando la propiedad de la traslacin real (2.4)
Z {u(t T)} = z
1
z
z 1
=
1
z 1
10
ITESM-CEM
Figure 2-2:
Figure 2-3:
Teorema de la tralacin compleja
Si F(z) es la transformada Z de f(t), entonces,
Z

e
at
f(t)
= F(ze
aT
)
la transformada. Esto se conoce como el teorema de la traslacin compleja
Ejemplo 8 Determine la transforma Z de g(t) = te
t
Solucin 9 Note que f(t) = t, a partir de tablas, se tiene que
F(z) =
Tz
1
(1 z
1
)
2
por lo que
G(z) =
T

ze
T
1 (ze
T
)
1
2
=
T

z
1
e
T
(1 z
1
e
T
)
2
=
Tze
T
(z
2
e
T
)
2
11
ITESM-CEM
ITESM-CEM
Teorema del valor inicial
Si X(z) = z {x(t)} = z {x(k)} y si el lim
z
X(z) existe, entonce el valor inicial x(0) de x(t) o de
x(k) est dado por
x(0) = lim
z
X(z) (2.6)
Para probar este teorema, note que
X(z) =

X
k=0
x(k)z
k
= x(0) + x(1)z
1
+ x(2)z
2
+ ....
al tomar el lmite cuando z tiende a innito
lim
z
X(z) = x(0) +
x(1)
z
+
x(2)
z
2
+ .... = x(0)
Ejemplo 10 Determine el valor inicial x(0) si X(z) est dada por
X(z) =
1 z
1
+ Tz
1
(1 z
1
)
2
Solucin 11
lim
z
X(z) = lim
z
1
1
z
+ T
1
z
1
1
z
2
= 1
X(z) es la transformada z de x(t) = 1 + t, por lo que x(0) = 1
Teorema del valor nal
Si X(z) = z {x(k)}, donde x(k) = 0 para k < 0 y que X(z) es estable, es decir, x(k) permanesca
nita (k = 0, 1, 2, ...). Entonces el valor nal de x(k) puede darse mediante
lim
k
x(k) = lim
z1
1 z
1
X(z)
(2.7)
Ejemplo 12 Determine el valor nal x() de
X(z) =
1
1 z
1

1
1 e
aT
z
1
, a > 0
12
ITESM-CEM
mediante el uso del teorema del valor nal
Solucin 13
x() = lim
z1
1 z
1
e
aT
z
1
+ z
1
(1 z
1
) (1 e
aT
z
1
)
= 1
Note que X(z) es la transformada z de x(t) = 1 e
at
, por lo tanto x() = 1
2.2 La transformada Z inversa
Como ya se menciono al inicio de este captulo, la transformada Z en sistemas de control en
tiempo discreto juega el mismo papel que la transformada de Laplace en sistemas de control
continuo, por lo que es necesario, al igual que en Laplace, obtener la transformada Z inversa
para que esta transformada sea til.
La notacin para la transformada Z inversa es Z
1
. La transformada Z
1
de X(z) da como
resultado la correspondiente secuencia en el tiempo x(k). Note que a partir de la transfor-
mada Z
1
slo se obtiene la secuencia de tiempo en los instantes de muestreo, por lo que la
transformada Z
1
de X(z) da una nica x(k) pero no da una nica x(t), es decir, se obtiene
una secuencia de tiempo que especica los valores de x(t) solamente en los valores discretos de
tiempo, t = 0, T, 2T, ..., y no dice nada acerca de los valores de x(t) en todos los otros tiempos.
Esto es, muchas funciones del tiempo x(t) diferentes pueden tener la misma x(kT).
Existen diversos mtodos para obtener la transformada Z
1
que no implican el uso de tablas:
1. Mtodo de la divisin directa
2. Mtodo de expansin en fracciones parciales
3. Mtodo de la integral de inversin
4. Mtodo computacional
Polos y ceros en el plano Z
La ubicacin de los polos y ceros de X(z) determina las caractersticas de x(k), la secuencia de
valores o nmeros. Para encontrar los polos y los ceros de X(z), es conveniente expresar X(z)
13
ITESM-CEM
como un cociente de polinomios en z. Por ejemplo, cuntos polos y ceros tiene la siguiente
funcin X(z)
X(z) =
1 + z
1
z
2
(1 + 0.2z
1
) (1 + 0.5z
1
)
al reescribirla en potencias positivas
X(z) =
1 +
1
z
1
z
2
1 +
0.2
z

1 +
0.5
z
=
z+1
z
1
z
2
z+0.2
z
z+0.5
z
=
z + 1
z(z + 0.2)(z + 0.5)
por lo que X(z) tiene polos ubicados en z = 0, z = 0.2 y z = 0.5 y tiene un cero, ubicado
en z = 1.
2.2.1 Mtodo de la divisin directa
En el mtodo de la divisin directa, la transformada z inversa se obtiene mediante la expansin
de X(z) en una serie innita de potencias de z
1
. Este mtodo es til cuando es dcil obtener
una expresin en forma cerrada o se desea encontrar slo algunos de los primeros trminos de
la secuencia x(k).
El mtodo de la divisin directa proviene del hecho de que si X(z) est expandida en una
serie de potencias de z
1
, esto es, si
X(z) =

X
k=0
x(kT)z
k
= x(0) + x(T)z
1
+ x(2T)z
2
+ x(3T)z
3
+ ....
entonces, x(kT) es el coeciente del trmino z
k
. Por lo tanto, los valores de x(kT) para
k = 0, 1, 2, ... ,se pueden determinar por inspeccin.
Ejemplo 14 Encuentre x(k) para k = 0, 1, 2, 3, 4 cuando X(z) est dada por
X(z) =
z + 1
(z + 0.2)(z + 0.5)
Solucin 15 Primero, X(z) se reescribe como un cociente de polinomios en z
1
X(z) =
1 + z
1
(1 + 0.2z
1
) (1 + 0.5z
1
)
14
ITESM-CEM
y realizando una divisin algebraica
z
1
+ z
2
1 + 0.7z
1
+ 0.1z
2
= z
1
+ 0.3z
2
0.31z
3
+ 0.187z
4
+ ...
por inspeccin se obtiene la solucin para x(k), con k = 0, 1, 2, 3, 4, es decir,
x(0) = 0
x(1) = 1
x(2) = 0.3
x(3) = 031
x(4) = 0.187
.
.
.
Ejemplo 16 Encuentre x(k) para k = 0, 1, 2, 3, 4 cuando X(z) est dada por
X(z) =
1
z 1
2.2.2 Mtodo de expansin en fracciones parciales
El mtodo de expansin en fraciones parciales que se presenta aqu es idntico al mtodo de
expansin en fracciones parciales que se utiliza en la Transformada de Laplace. Para encontrar
la transformada Z inversa de X(z) por fracciones parciales, primero se factoriza el polinomio
denominador de X(z) y se encuentran los polos
X(z) =
N(z)
D(z)
=
b
0
z
m
+ b
1
z
m1
+ ... + b
m1
z + b
m
(z + p
1
)(z + p
2
)...(z + p
n
)
luego se expande X(z) en fracciones parciales de manera que cada uno de lo trminos sea
identicado faclmente utilizando las tablas de transformada Z. Si X(z) tiene uno o ms ceros
en el origen (z = 0), entonces
X(z)
z
se debe expandir en lugar de X(z).
15
ITESM-CEM
Polos diferentes
Considrese que todos los polos son diferentes y que hay por lo menos un cero en el origen,
entonces se aplicar fracciones parciales a
X(z)
z
y se expandir de la forma
X(z)
z
=
a
1
z + p
1
+
a
2
z + p
2
+ ... +
a
n
z + p
n
donde las a
i
son constantes y se denominan residuos de la raz z = p
i
. La frmula para obtener
el residuo es
a
i
=

(z + p
i
)
X(z)
z
z=p
i
Ejemplo 17 Determine la transformada Z inversa de
X(z) =
z
(z 1) (z 0.1)
Solucin 18 (a) Note que existe un cero en el origen, por lo que se expande en fracciones
parciales
X(z)
z
, es decir
X(z)
z
=
1
(z 1) (z 0.1)
=
a
z 1
+
b
z 0.1
donde
a =
1
(z 1) (z 0.1)
(z 1)
z=1
= 1.1
b =
1
(z 1) (z 0.1)
(z 0.1)
z=0.1
= 1.1
por lo tanto
X(z) = 1.1
z
z 1
1.1
z
z 0.1
= 1.1
1
1 z
1
1.1
1
1 0.1z
1
x(k) = 1.1 (1)
k
1.1 (0.1)
k
Qu sucedera si no se "guarda" el cero en el origen?
Solucin 19 (b)
X(z) =
z
(z 1) (z 0.1)
=
a
z 1
+
b
z 0.1
16
ITESM-CEM
donde
a =
z
(z 1) (z 0.1)
(z 1)
z=1
= 1.1
b =
z
(z 1) (z 0.1)
(z 0.1)
z=0.1
= 0.11
sustituyendo los valores de los residuos
X(z) = 1.1
1
z 1
0.11
1
z 0.1
= 1.1
z
1
1 z
1
0.11
z
1
1 0.1z
1
por lo tanto
x(k) = 1.1 (1)
k1
0.11 (0.1)
k1
Note que al parecer las dos soluciones son distintas, ser cierto?, Por qu la diferencia
entre las dos soluciones? Verique si en realidad son diferentes.
Ejemplo 20 Determine la transforma z inversa de
X(z) =
1
z
2
(z 1)
Reeescribiendo en potencias negativas
X(z) =
z
3
1 z
1
Note que s F(z) =
1
1z
1
, entonces f(k) = u(k) = (1)
k
. Utilizando el teorema de la traslacin
real
Z {x(t nT)} = z
n
X(z)
se tiene que
Z
1
z
3
F(z)
= Z
1
z
3
1
1 z
1
= u(k 3) = (1)
k3
17
ITESM-CEM
Polos mltiples
Sea X(z) escrita en forma factorizada
X(z) =
N(z)
D(z)
=
N(z)
(z + p
1
)
r
(z + p
r+1
) + ... + (z + p
n
)
es decir, existen r races mltiples y n r races diferentes. la expansin de X(z) en fracciones
parciales es
X(z) =
b
r
(z + p
1
)
r
+
b
r1
(z + p
1
)
r1
+ ... +
b
1
z + p
1
+
a
r+1
z + p
r+1
+
a
n
z + p
n
donde b
r
, b
r1
, b
1
estn dados por
b
r
=
h
(z + p
1
)
r
N(z)
D(z)
i
z=p
1
b
r1
=
n
d
dz
h
(z + p
1
)
r
N(z)
D(z)
io
z=p
1
.
.
.
b
rj
=
n
1
j!
d
j
dz
j
h
(z + p
1
)
r
N(z)
D(z)
io
z=p
1
Ejemplo 21 Determine la transformada Z inversa de
X(z) =
z
(z 1)
2
(z 2)
X(z)
z
=
1
(z 1)
2
(z 2)
=
a
z 2
+
b
2
(z 1)
2
+
b
1
z 1
donde
a =
1
(z1)
2
(z2)
(z 2)
z=2
= 1
b
2
=
1
(z1)
2
(z2)
(z 1)
2
z=1
= 1
b
1
=
d
dz
1
(z1)
2
(z2)
(z 1)
2
z=1
=
d
dz
1
(z2)
z=1
=
1
(z2)
2
z=1
= 1
Por lo tanto
X(z) =
z
z 2

z
(z 1)
2

z
z 1
=
1
1 2z
1

z
1
(1 z
1
)
2

1
1 z
1
18
ITESM-CEM
x(k) = (2)
k
k (1)
k1
(1)
k
Polos complejos y conjugados
Esta metodologa es ilustrada por medio de una serie de ejemplos.
Ejemplo 22 Determine la transformada z inversa de
X(z) =
z
2
z
2
z + 1
Note que las races son complejas y conjugadas. Reescribiendo la expresin anterior en potencias
negativas
X(z) =
1
1 z
1
+ z
2
y utilizando la tabla de transformada z, se identica que
2 cos(T) = 1 cos(T) =
1
2
T = cos
1
(T)
Por lo tanto
X(z) =
1
1
2
z
1
+
1
2
z
1
1 z
1
+ z
2
=
1
1
2
z
1
1 z
1
+ z
2
+
1
2 sin(T)
sin(T) z
1
1 z
1
+ z
2
x(k) = cos(kT) +
1
2 sin(T)
sin(kT)
Ejemplo 23 Determine la transformada z inversa de
X(z) =
z
2
2z
2
z + 1
19
ITESM-CEM
X(s) x(t) x(kT) o x(k) X(z)
s
2
+
2
sin(t) sin(kT)
z
1
sinT
12z
1
cos T+z
2
s
s
2
+
2
cos (t) cos (kT)
1z
1
cos T
12z
1
cos T+z
2
(s+)
2
+
2
e
t
sin(t) e
kT
sin(kT)
e
T
z
1
sinT
12e
T
z
1
cos T+e
2T
z
2
s+
(s+)
2
+
2
e
t
cos (t) e
kT
cos (kT)
1e
T
z
1
cos T
12e
T
z
1
cos T+e
2T
z
2
Tabla de Transformada Z
Ejercicio 24 Determine la transformada z inversa de las siguientes funciones
X(z) =
1 e
aT
z
(z 1) (z e
aT
)
X(z) =
1 e
aT
z
(z 1)
2
(z e
aT
)
X(z) =
z
2
+z + 2
(z 1) (z
2
z + 1)
2.3 Ecuaciones en diferencias
Un sistema en tiempo discreto, lineal e invariante en el tiempo es caracterizado por la ecuacin
x(k) +a
1
x(k 1) +... +a
n
x(k n) = b
0
u(k) +b
1
u(k 1) +... +b
m
u(k m)
donde u(k) y x(k) son la entrada y la salida del sistema respectivamente, en la k-sima
iteracin. El problema es obtener una expresin para x(k). Denase
X(z) = z {x(kT)} = z {x(k)}
entonces, aplicando el teorema de la translacin real se obtienen la Transformada Z de cada
uno de los trminos que conforman la ecuacin en diferencias. Posteriomente, se despeja la
variable dependiente X(z) y se aplica la transformada Z
1
para obtener la secuencia x(k)
Ejemplo 25 [1] Obtenga la solucin de la siguiente ecuacin en diferencias en trminos de
20
ITESM-CEM
x(0) y x(1)
x(k + 2) + (a +b) x(k + 1) +abx(k) = 0
donde a y b son constantes y k = 0, 1, 2, ...
Ejemplo 26 [1]Resuelva la siguiente ecuacin en diferencias
2x(k) 2x(k 1) +x(k 2) = u(k)
donde x(k) = 0 para k < 0 y
u(k) =
1, k = 0, 1, 2, ...
0, k < 0
21
ITESM-CEM
Chapter 3
Anlisis de sistemas de control en
tiempo discreto
3.1 Muestreo mediante impulsos y retencin de datos
Considrese un muestreador cticio cuya salida se considera como un tren de impulsos que
comienza en t = 0, con el perodo de muestreo igual a T y la magnitud de cada impulso igual al
valor muestreado de la seal en tiempo continuo en el instante de muestreo correspondiente (ver
gura 3-1). La seal muestreada x
(t), se puede representar mediante una sumatoria innita

1
x
(t) =

X
k=0
x(kT)(t kT)
= x(0)(t) +x(T)(t T) +... +x(kT)(t kT) +... (3.1)
1
A lo largo de estas notas, se supone que la operacin de muestreo es uniforme; esto es, slo existe un perodo
de muestreo en el sistema el cual es constante. Si un sistema de control en tiempo discreto incluye dos o ms
muestreadores en el sistema, se supone que los muestreadores estn sincronizados y tienen la misma frecuencia
de muestreo.
22
ITESM-CEM
La transformada de Laplace de la ecuacin (3.1)
X
(s) = L[x
(t)] = x(0)L[(t)] +x(T)L[(t T)] +x(2T)L[(t 2T)] +...

= x(0) +x(T)e
Ts
+x(2T)e
2Ts
+...
=

X
k=0
x(kT)e
kTs
Figure 3-1: Muestreador mediante impulsos
Si se dene
z = e
Ts
o
s =
1
T
lnz
entonces
X
(s)|
s=
1
T
lnz
= X(z) (3.2)
La transformada de Laplace de la seal muestreada mediante impulsos x
(t) es la transformada
Z de la seal x(t) si e
Ts
se dene como z, es decir z = e
Ts
.
3.1.1 Retenedor de orden cero
En un muestreador ideal, un interruptor se cierra cada perodo de muestreo T para admitir una
seal de entrada. Un muestreador convierte una seal de tiempo continuo en un tren de pulsos
que se presenta en los instantes de muestreo t = 0, T, 2T, ...
23
ITESM-CEM
La retencin de datos es un proceso de generacin de una seal de tiempo continuo h(t) a
partir de una secuencia en tiempo discreto x(kT). Un circuito de retencin convierte la seal
muestreada en una seal de tiempo continuo, que reproduce aproximadamente la seal aplicada
al muestreador. El circuito de retencin ms simple es el Retenedor de Orden Cero (ROC), este
Figure 3-2:
circuito retiene la amplitud de la muestra en un instante de muestreo al siguiente. La funcin
de transferencia del ROC, G
ROC
(s) es
G
ROC
(s) =
1 e
Ts
s
(3.3)
3.1.2 La funcin de transferencia pulso.
La funcin de transferencia relaciona las transformada de Laplace de la seal de salida con la
correspondiente entrada del sistema, mientras que la funcin de transferencia pulso relaciona
las transformadas Z de salida en los instantes de muestreo con la correspondiente entrada
muestreada.
Considere la respuesta de un sistema continuo excitado por una seal muestreada como se
muestra en la gura 3-3. La seal muestreada mediante impulsos x
(t) es la entrada al sistema

cuya funcin de transferencia es G(s). Se supone que la salida del sistema es una seal en
24
ITESM-CEM
Figure 3-3: Sistema en tiempo continuo G(s) excitado con una seal muestreada mediante
impulsos
tiempo continuo y(t)
1
. Entonces la salida Y (s) es
Y (s) = G(s)X
(s)
Considere que en la salida hay otro muestreador, sincronizado en fase con el muestreador de la
entrada, y ambos operan con el mismo perodo de muestreo, entonces la salida es
Y
(s) = [G(s)X
(s)]
= G
(s)X
(s)
De este modo, utilizando la relacin (3.2) se obtiene
Y (z) = G(z)X(z)
Ejemplo 27 Considere los sistemas que se muestran en la gura 3-4. Donde G
1
(s) =
1
s
,
G
2
(s) =
1
s+1
. Obtenga la funcin de transferencia pulso
Y (z)
X(z)
para cada uno de estos sistemas
Clculo de la transformada Zque involucran un retenedor de orden cero (ROC)
Suponga que la funcin de transferencia G(s) sigue de un ROC. Entonces el producto de la
funcin de transferencia del ROC y de G(s) se convierte en
X(s) =
1 e
Ts
s
G(s) (3.4)
1
Se supone que x(t) < 0 e y(t) < 0 para t < 0.
25
ITESM-CEM
Figure 3-4: Sistemas muestreados
Note que (3.4) se puede escribir como
X(s) =
1 e
Ts
G(s)
s
=
1 e
Ts
G
1
(s) = G
1
(s) e
Ts
G
1
(s) (3.5)
la transformada Z de (3.5) es dada por
X(z) = G
1
(z) z
1
G
1
(z) =
1 z
1
G
1
(z)
Por lo tanto, para obtener la transformada Z de X(s), el trmino 1 e
Ts
= 1 z
1
y
nicamente hay que obtener la transformada Z
n
G(s)
s
o
Ejemplo 28 Obtenga la transformada Z de
G(s) =
1 e
Ts
s
1
s
Solucin 29
G(z) =
1 z
1
1
s
2
1 z
1
Tz
1
(1 z
1
)
2
=
Tz
1
1 z
1
=
T
z 1
26
ITESM-CEM
3.1.3 Funcin de transferencia pulso de sistemas en lazo cerrado
Para determinar la funcin de transferencia pulso del sistema de control en lazo cerrado que se
muestra en la gura (3-5), considere el siguiente algoritmo:
Figure 3-5: Sistema de control en lazo cerrado
Figure 3-6: Grafo del sistema en lazo cerrado
1. Denomine la entrada de los muestreadores como E(s) y la salida de los muestreadores
como E
(s)
2. Trace el grafo del sistema (vase la gura 3-6)
3. Escriba las salidas del grafo en funcin de las entradas del grafo (se considera que las
entradas del grafo son las entradas del sistema y las salidas de los muestreadores y que
las salidas del grafo son las salidas del sistema y las entradas al muestreador). En este
caso, las ecuaciones quedan de la siguiente forma
Y
(s) = G
1
(s)G
2
(s)E
(s)
E(s) = R(s) G
1
(s)G
2
(s)E
(s)
(3.6)
27
ITESM-CEM
4. Muestree las ecuaciones obtenidas en el paso anterior
Y
(s) = G
1
G
2
(s)E
(s)
E
(s) = R
(s) G
1
G
2
(s)E
(s)
(3.7)
5. Resuelva las ecuaciones por el mtodo ms conveniente:
(a) Masson
Figure 3-7: Grafo de lazo cerrado muestreado
Y
(s)
R
(s)
=
G
1
G
2
(s)
1 +G
1
G
2
(s)
Y (z)
R(z)
=
G
1
G
2
(z)
1 +G
1
G
2
(z)
(b) Sustitucin
E
(s) = R
(s) G
1
G
2
(s)E
(s)
E
(s) =
R
(s)
1+G
1
G
2
(s)
Y
(s) = G
1
G
2
(s)E
(s) =
G
1
G
2
(s)
1+G
1
G
2
(s)
R
(s)
por lo tanto
Y (z)
R(z)
=
G
1
G
2
(z)
1 +G
1
G
2
(z)
Ejercicio 30 Obtenga la funcin de transferencia pulso de lazo cerrado del siguiente diagrama
a bloques
28
ITESM-CEM
Figure 3-8: Diagrama a bloques discreto
29
ITESM-CEM
Chapter 4
Diseo de controladores digitales
4.1 Correspondencia entre el plano s y el plano z
En el plano complejo s, la ubicacin de los polos y los ceros nos permitian predecir el compor-
tamiento dinmico el sistema, de aqu la importancia de estudiar la relacin entre los requisitos
de diseo (por ejemplo tiempo de establecimiento t
s
, mximo sobreimpulso M
p
) con la ubi-
cacin de los polos en el plano s. De igual manera, en los sistemas discretos es muy importante
la ubicacin de los polos y los ceros en el plano z. A continuacin determinaremos la relacin
existente entre el plano s y el plano z.
En un proceso donde se encuentre involucrado un muestreo por impulsos, las variables
complejas s y z se encuentran relacionadas por
z = e
Ts
dado que la variable compleja s est formada por una parte real y una parte imaginaria ,
es decir
s = +j
s
se tiene que
z = e
T(+j)
= e
T
e
jT
Si se considera un punto representativo en el eje j en el plano s, y conforme este punto se
30
ITESM-CEM
Figure 4-1: Plano SPlano Z
mueve desde j
s
2
hasta j
s
2
siendo
s
la frecuencia de muestreo, tenemos que |z| = 1 y ]z
vara de a en direccin contraria a las manecillas del reloj en el plano z. Conforme
el punto representativo se mueve desde j
s
2
hasta j
3s
2
, el punto correspondiente en el plano
z traza un criculo unitario en direccin contraria a las menecillas del reloj. Por lo tanto,
conforme el punto se mueve en el eje j del plano s dibujaremos un crculo unitario en el plano
z un nmero innito de veces.
Tiempo de establecimiento t
s
El tiempo de establecimiento queda determinado por el valor de atenuacin de los polos
dominantes en lazo cerrado. Si se especica el tiempo de establecimiento, se puede dibujar una
lnea =
1
corresponde en el plano z a la parte interior de un crculo de radio e
1
T
31
ITESM-CEM
Lugar geomtrico de frecuencias constantes
Un lugar geomtrico de frecuencia constante =
1
en el plano s corresponde en el plano z a
una lnea radial de ngulo constante T
1
(en radianes)
Figure 4-2:
Lugar geomtrico de amortiguamiento constante
Una lnea de factor de amortiguamiento constante (una lnea radial) en el plano s cofrresponde
a una espiral en el plano z
s =
n
+j
n
q
1
2
=
n
+j
d
z = e
Ts
= e
T(n+j
d
)
Por lo tanto
|z| = exp(
n
T) = exp
2
n
p
1
2
s
p
1
2
!
= exp
2
d
s
p
1
2
!
y
]z =
d
T =
2
d
s
Entonces, la magnitud de z se reduce y el ngulo de z se aumenta linealmente conforme
d
se
incremente y el lugar geomtrico en el plano z se convierte en una espiral logartmica
32
ITESM-CEM
Figure 4-3:
33
ITESM-CEM
Chapter 4
Diseo de controladores digitales
4.1 Correspondencia entre el plano s y el plano z
En el plano complejo s, la ubicacin de los polos y los ceros nos permitian predecir el compor-
tamiento dinmico el sistema, de aqu la importancia de estudiar la relacin entre los requisitos
de diseo (por ejemplo tiempo de establecimiento t
s
, mximo sobreimpulso M
p
) con la ubi-
cacin de los polos en el plano s. De igual manera, en los sistemas discretos es muy importante
la ubicacin de los polos y los ceros en el plano z. A continuacin determinaremos la relacin
existente entre el plano s y el plano z.
En un proceso donde se encuentre involucrado un muestreo por impulsos, las variables
complejas s y z se encuentran relacionadas por
z = e
Ts
dado que la variable compleja s est formada por una parte real y una parte imaginaria ,
es decir
s = +j
s
se tiene que
z = e
T(+j)
= e
T
e
jT
Si se considera un punto representativo en el eje j en el plano s, y conforme este punto se
30
ITESM-CEM
Figure 4-1: Plano SPlano Z
mueve desde j
s
2
hasta j
s
2
siendo
s
la frecuencia de muestreo, tenemos que |z| = 1 y ]z
vara de a en direccin contraria a las manecillas del reloj en el plano z. Conforme
el punto representativo se mueve desde j
s
2
hasta j
3s
2
, el punto correspondiente en el plano
z traza un criculo unitario en direccin contraria a las menecillas del reloj. Por lo tanto,
conforme el punto se mueve en el eje j del plano s dibujaremos un crculo unitario en el plano
z un nmero innito de veces.
Tiempo de establecimiento t
s
El tiempo de establecimiento queda determinado por el valor de atenuacin de los polos
dominantes en lazo cerrado. Si se especica el tiempo de establecimiento, se puede dibujar una
lnea =
1
corresponde en el plano z a la parte interior de un crculo de radio e
1
T
31
ITESM-CEM
Lugar geomtrico de frecuencias constantes
Un lugar geomtrico de frecuencia constante =
1
en el plano s corresponde en el plano z a
una lnea radial de ngulo constante T
1
(en radianes)
Figure 4-2:
Lugar geomtrico de amortiguamiento constante
Una lnea de factor de amortiguamiento constante (una lnea radial) en el plano s cofrresponde
a una espiral en el plano z
s =
n
+j
n
q
1
2
=
n
+j
d
z = e
Ts
= e
T(n+j
d
)
Por lo tanto
|z| = exp(
n
T) = exp
2
n
p
1
2
s
p
1
2
!
= exp
2
d
s
p
1
2
!
y
]z =
d
T =
2
d
s
Entonces, la magnitud de z se reduce y el ngulo de z se aumenta linealmente conforme
d
se
incremente y el lugar geomtrico en el plano z se convierte en una espiral logartmica
32
ITESM-CEM
Figure 4-3:
Ejemplo 31 Determine la regin en el plano complejo z de tal forma que se cumplan con las
siguientes condiciones: Tiempo de establecimiento t
s
< 4s, Tiempo pico t
p
< 3s y un mximo
sobreimpulso M
p
< 10%. Considere que el perodo de muestreo es T = 0.5seg
4.2 Anlisis de estabilidad en el plano z
Considere la siguiente funcin de transferencia pulso en lazo cerrado
Y (z)
R(z)
=
G(z)
1 +GH(z)
(4.1)
La estabilidad del sistema en lazo cerrado se puede determinar por las ubicaciones de los polos
de lazo cerrado en el plano z, es decir, por las races de la ecuacin caracterstica
1 +GH(z) = 0 (4.2)
como sigue:
1. Para que el sistema sea estable los polos de lazo cerrado deben estar dentro del crculo
unitario. Cualquier polo en lazo cerrado, fuera del crculo unitario hace inestable al
sistema
2. Si un polo simple se presenta en z = 1, el sistema es criticamente estable. Si existe un par
de polos complejos y conjugados sobre el crculo unitario, tambin hace que el sistema sea
33
ITESM-CEM
crticamente estable
3. Los ceros de lazo cerrado no afectan la estabilidad y por lo tanto pueden quedar localizados
en cualquier parte del plano complejo z.
Ejemplo 32 La funcin de transferencia de lazo cerrado de un sistema de control digital, es
expresada por
G(z) =
z
z
2
z + 2
El sistema es estable o inestable?
Solucin 33 Las races del sistema son
z = 0.5 j1.232
Por lo tanto la magnitud de los polos son 1.32, por lo que los polos estan fuera del crculo
unitario, lo cual implica que el sistema es inestable
Ejemplo 34 Considere el siguiente sistema de control
donde
Figure 4-4:
G
ROC
(s) =
1 e
Ts
s
y G(s) =
1
s + 1
Si T = 0.5s, determine el rango de k que garantice la estabilidad del sistema en lazo cerrado
Solucin 35 Sea
G
T
= z
1 e
Ts
s
1
s + 1
1 z
1
1
s (s + 1)
= k
1 e
T
z e
T
34
ITESM-CEM
Por lo tanto, la ecuacin caracterstisca es
z e
T
+k
1 e
T
= 0
La estbilidad del sistema es vericada al considerar que la magnitud del polo es menor que uno,
es decir
|z| =
e
T
k
1 e
T
< 1
de lo anterior se tiene que
e
T
k
1 e
T
< 1 k > 1
y
e
T
k
1 e
T
< 1 k <
1 +e
T
1 e
T
4.2.1 La prueba de estabilidad de Jury
Esta prueba de estabilidad nos indica la existencia de cualquier raz inestable (races en el plano
z fuera del crculo unitario). Sin embargo, no indica nada sobre las ubicaciones de las races
inestable. Para aplicar la prueba de estabilidad de Jury a una ecuacin caracterstica dada
P(z) = 0, se construye una tabla cuyos elementos se basan en los coecientes de P(z) = 0.
Suponga que la ecuacin caracterstica es de la forma
P(z) = a
0
z
n
+a
1
z
n1
+a
2
z
n2
+... +a
n1
z +a
n
(4.3)
35
ITESM-CEM
donde a
0
> 0. Entonces la tabla de Jury se construye como sigue
Rengln z
0
z
1
z
2
z
3
z
n2
z
n1
z
n
1 a
n
a
n1
a
n2
a
n3
a
2
a
1
a
0
2 a
0
a
1
a
2
a
3
a
n2
a
n1
a
n
3 b
n1
b
n2
b
n3
b
n4
b
1
b
0
4 b
0
b
1
b
2
b
3
b
n2
b
n1
.
.
.
2n 5 p
3
p
2
p
1
p
0
2n 4 p
0
p
1
p
2
p
3
2n 3 q
2
q
1
q
0
donde
b
k
=
a
n
a
n1k
a
0
a
k+1
, k = 0, 1, 2, ..., n 1
.
.
.
q
k
=
p
3
p
2k
p
0
p
k+1
, k = 0, 1, 2
Note que el ltimo rengln de la tabla est formado por tres elementos (para sistemas de segundo
orden, 2n 3 = 2(2) 3 = 1 la tabla estar formada por un rengln)
Un sistema con la ecuacin caracterstica dad por (4.3), el cual por comodidad se reescribe
P(z) = a
0
z
n
+a
1
z
n1
+a
2
z
n2
+... +a
n1
z +a
n
donde a
0
> 0, es estable, si todas las condiciones siguientes se satisfacen:
1. |a
n
| < a
0
2. P(z)|
z=1
> 0
3. P(z)|
z=1
> 0 para n par

< 0 para n impar
36
ITESM-CEM
4.
|b
n1
| > |b
0
|
|c
n2
| > |c
0
|
.
.
.
|q
2
| > |q
0
|
Ejemplo 36 Examine la estabilidad de la ecuacin caracteristica siguiente
P(z) = z
3
1.5z
2
0.1z + 0.5 = 0
Solucin 37 Antes de construir la tabla se verifcan las condiciones 1, 2 y 3. Note que a
0
= 0.5
y a
n
= a
3
= 1, por lo tanto |0.5| < 1 por lo que la primera condicin se cumple.
P(z)|
z=1
= (1)
3
1.5 (1)
2
0.1 (1) + 0.5 = 0.1 < 0
Esta condicin no se satisface por lo que el sistema es inestable
Ejemplo 38 Examine la estabilidad de la ecuacin caracteristica siguiente
P(z) = z
4
z
3
+ 0.5z
2
+ 0.1z 0.5 = 0
Solucin 39 Nuevamente, antes de construir la tabla se verifcan las condiciones 1, 2 y 3.
Note que a
0
= 0.5 y a
n
= a
4
= 1, por lo tanto |0.5| < 1 por lo que la primera condicin se
cumple.
P(z)|
z=1
= (1)
4
(1)
3
+ 0.5 (1)
2
0.1 (1) + 0.5 = 0.9 > 0
P(z)|
z=1
= (1)
4
(1)
3
+ 0.5 (1)
2
0.1 (1) + 0.5 = 3.1 > 0
37
ITESM-CEM
Ahora construimos la tabla de estabilidad de Jury
Rengln z
0
z
1
z
2
z
3
z
4
1 0.5 0.1 0.5 1 1
2 1 1 0.5 0.1 0.5
3 0.75 0.95 0.75 0.4
4 0.4 0.75 0.95 0.75
5 0.4 0.41 0.18
Vericando las condiciones
|b
n1
| > |b
0
|
que en este caso es
|0.75| > |0.4|
la cual se cumple. La siguiente condicin es
|q
2
| > |q
0
|
al sustituir los valores
|0.4| > |0.18|
Por lo tanto, la ecuacin caracterstica dad es estable, o lo que es lo mismo, todas las races
estn dentro del crculo unitario.
Ejemplo 40 Considere un sistema de control en tiempo discreto con retroalimentacin unitaria
cuya funcin de transferencia pulso es
GH(z) =
k (0.5z + 1)
(z 1) (z 0.36)
considere que T = 1seg. Determine el rango de valores de la ganancia k para garantizar la
estabilidad
38
ITESM-CEM
4.2.2 Anlisis de estabilidad mediante la transformacin bilineal y el criterio
de Routh
La transformacin bilineal denida por
z =
w + 1
w 1
misma que, al ser resuelta en funcin de w, resulta en
w =
z + 1
z 1
hace corresponder el interior del crculo unitario del plano z con el semiplano izquierdo del
plano w.
Considere la ecuacion caracterstica (4.3), es decir
P(z) = a
0
z
n
+a
1
z
n1
+a
2
z
n2
+... +a
n1
z +a
n
y sustituyendo
w+1
w1
en lugar de z
a
0
w + 1
w 1
n
+a
1
w + 1
w 1
z
n1
+a
n2
2
w + 1
w 1
+... +a
n1
w + 1
w 1
+a
n
= 0
Al multiplicar ambos miembros de la ecuacin por (w 1)
n
se obtiene
Q(w) = b
0
w
n
+b
1
w
n1
+b
2
w
n2
+... +b
n1
w +b
n
= 0
es posible aplicar el criterio de estabilidad de Routh de la misma forma que en los sistemas de
tiempo continuo.
4.3 Lugar de las races
El mtodo del lugar de las races desarrollados para sistemas en tiempo continuo puede ser
extendido sin modicaciones a sistemas discretos en el tiempo, excepto por la estabilidad, la
cual se modica, del eje j en el plano S al crculo unitario en el plano Z.
39
ITESM-CEM

Transformada Z

Cargado por

Copyright:

Formatos disponibles

Transformada Z

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Transformada Z

Cargado por

Copyright:

Formatos disponibles

Contents

F(z) f(0) f(T)z

F(z) f(0) f(T)z

(t), se puede representar mediante una sumatoria innita

(t)] = x(0)L[(t)] +x(T)L[(t T)] +x(2T)L[(t 2T)] +...

(t) es la entrada al sistema

> 0 para n par

También podría gustarte