UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA

Control Digital
________________________________________________________________________________________

Captulo 1
Introduccin a Sistemas de Control Digital
El desarrollo y uso de computadores digitales cada vez ms potentes, veloces y
econmicos, han sido decisivos para el empleo de tcnicas digitales para el control de
sistemas dinmicos, como por ejemplo, el control inteligente en robots industriales, la
optimizacin econmica en el uso de combustible en los automviles, refinamiento en la
operacin de equipos de uso domstico, sistema de pilotaje de aviones, etc.

1.1 Sistemas de Control Digital

El esquema bsico de un sistema de control digital, se muestra en la figura 1.1. El sistema
incluye un control prealimentado y realimentado. La estrategia de control a usar depender
del diseador, en funcin a una determinada aplicacin.
Perturbacin
O
ruido

Controlador
Prealimentado

r
Refer.

Muestreador
y
A/D

Computadora
Digital

D/A
y
retenedor

y
Salida
Actuador

Planta

Controlador digital

Transductor
o
sensor

Filtraje

n
Ruido
Figura 1.1: Diagrama de bloques de un sistema de control digital.
______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

Implementacin del
algoritmo de control

PC
Computadora de Procesos
Microcontrolador
DSP

1.2 Muestreo y Retencin

El proceso de muestreo es sumamente til en sistemas de control donde se usa un
controlador digital. El muestreo de seales consiste en tomar muestras en puntos discretos
de tiempo, que luego son cuantificadas por un proceso de retencin, para luego pasar a un
convertidor analgico/digital, producindose su conversin en cdigos binarios. Tales
cdigos, denominados tambin datos binarios pueden ser procesados por el computador,
para producir una determinada decisicin.
En la gran mayora de las aplicaciones de control realimentado, los procesos o plantas
(por naturaleza) utilizan actuadores que son controlados por seales analgicas, situacin
que hace necesario la reconversin de la seal digital hacia analgico. Tal propsito se
logra por un proceso de conversin digital/analgico (D/A) y luego por un retenedor, que
haga posible la reconstruccin de la seal analgica, la cual excita al actuador controlando
directamente el proceso. En la figura 1.2 (a) se muestra un diagrama de bloques de un
sistema de adquisicin de datos, y en la figura 1.2 (b) se muestra un diagrama de bloques de
un sistema de reconstruccin de la seal analgica.

Transductor

Amplificador

Filtro
pasa-bajo

Multiplexor
Analgico

Muestreador
y retenedor

Al controlador digital
o PC

(a)

Variable fsica

Del
controlador
o PC
Registro

Convertidor
D/A

Demultiplexor

Convertidor
A/D

Al
actuador
Retenedor

(b)
Figura 1.2: a) Diagrama de bloques de un sistema de adquisicin de datos;
b) diagram de bloques de un sistema de reconstruccin de seal analgica.

______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

Es importante anotar que la seleccin del periodo de muestreo T es determinante en los

resultados del sistema de control digital. Debemos recordar, por Nyquist, que la frecuencia
de muestreo debe ser

fs 2 fm

(1.1)

donde Fs es la frecuencia de muestreo y fm es la frecuencia mxima de la seal a ser

muestreada. Informacin amplia sobre sistemas de adquisicin de datos y reconversin
analgica puede encontrarse en cualquier texto de control digital o control en tiempo
discreto (vase [5]).
La figura 1.3 muestra el proceso de muestreo y retencin de la seal analgica u ,
permitiendo la reconstruccin de la seal por aproximaciones rectangulares.
u(t)

u(kT)

t
u(t)
u(s)

0 T 2T

Muestreador

(t)

6T

0 T 2T
Retenedor de
orden cero

u(kT)
u*(s)

6T
(t)
(s)

Figura 1.3: Muestreador y retenedor de orden cero.

La funcin de transferencia del retenedor de orden cero viene dada por

Gr 0

u ( s) 1 e Ts

s
u * ( s)

(1.2)

donde T es el periodo de muestreo.

A continuacin se presenta una lista de mtodos de retencin usados con Matlab.

______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

1.3 Discretizacin Directa

Es bastante til discretizar directamente expresiones que contengan integrales y
derivadas.En el desarrollo del curso nos encontraremos frecuentemente con funciones que
contengan derivadas e integrales, las cuales pueden ser discretizadas usando ecuaciones en
diferencias. Veamos:

d
f (kT ) f (kT ) f (kT T )
f (t )

dt
T
T
2
2
d
f (kT ) f (kT ) f (kT T )
f (t )

2
dt
T2
T2

f (t )dt Te(iT T )

(1.3)
(1.4)

i 0

Es necesario anotar que en este caso los perodos de muestreo deben ser pequeos para
que el resultado del proceso discreto sea bastante aproximado al obtenido mediante el
mtodo de discretizacin general (que es el que usa Matlab).

Ejemplo 1.1: Considerando las ecuaciones de estado y de salida del circuito RLC
siguiente:

______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

1
x2
c
1
R
1
x 2 x1 x 2 Vin
L
L
L
y x1
x1

(1.5)
(1.6)
(1.7)

Determine:
a) La discretizacin directa del sistema de ecuaciones
b) La grfica de su respuesta a un escaln unitario
Asuma un periodo de muestreo de 0.01 segundos.

Solucin
a) Del sistema de ecuaciones de estado y de salida del circuito RLC podemos
encontrar el equivalente discreto en adelanto de la derivada (un corrimiento a la
derecha), entonces, el sistema discreto para cada ecuacin diferencial ser:
T
x1 (kT 1) x1 (kT ) x 2 (kT )
(1.8)
c
1
R
1
x 2 (kT 1) x 2 (kT ) T [ x1 (kT ) x 2 (kT ) Vin (kT )]
(1.9)
L
L
L
y (kT ) x1 (kT )
(1.10)
Las ecuaciones (1.8), (1.9) y (1.10) pueden reescribirse implcitamente as:

T
x 2 (k )
c
1
R
1
x 2 (k 1) x 2 (k ) T [ x1 (k ) x 2 (k ) Vin (k )]
L
L
L
y (k ) x1 (k )
x1 (k 1) x1 (k )

(1.11)
(1.12)
(1.13)

b) La grfica de la respuesta al escaln unitario se obtiene mediante la ejecucin del

siguiente programa en MATLAB (ejem3_1):

______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

1.4

Te n s i n e n e l c a p a c it o r (V )

1.2

0.8

0.6

0.4

0.2

0
0

50

100

150

200

250

N m e ro d e m u e s t ra s (k )

1.4 La Transformada Z
La transformada Z es un mtodo operacional sumamente potente en el trabajo con sistemas
en tiempo discreto.
La transformada Z de una funcin f(t), toma en cuenta slo los valores muestreados
de f(t), es decir la secuencia de valores f(kT): f(0), f(T), f(2T), ... , donde k 0 y T es el
periodo de muestreo. Se define mediante la siguiente ecuacin:
______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

F ( z ) Z [ f (t )] Z [ f (kT )] f (kT ) z k f (0) f (T ) z 1 f (2T ) z 2

(1.14)

k 0

La relacin entre la transformada z y la transformada s viene dada por:

z eTs eT ( j ) eT [cos T jsenT ]

(1.15)

Para una secuencia de nmeros f(k), la transformada z se define como

F ( z ) Z [ f (k )] f (k ) z k

(1.16)

k 0

La transformada z dadas por las ecuaciones (3.14) y (3.6) se denomina transformada z

unilateral.
Si - t , o si k adopta valores enteros (k = 0, 1, 2, ...), entonces la
transformada z de f(t) o de f(k) viene dada por:
F ( z ) Z [ f (t )] Z [ f (kT )]

f (kT ) z

(1.17)

F ( z ) Z [ f (k )]

f (k ) z

(1.18)

Las ecuaciones (1.17) y (1.18) representan la transformada z bilateral.

1.4.1 Transformada z de algunas funciones elementales

La transformada z de algunas funciones elementales, considerando el caso unilateral, se
presenta a continuacin:
Funcin escaln unitario: Dada la funcin escaln unitario definida por

1(t ), t 0
f (t )
0, t 0

(1.19)

la transformada z de f(t) usando la ecuacin (1.14) viene dada por

______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

k 0

k 0

F ( z ) Z [1(t )] 1z k z k
1

1 z z
1
z

1
z 1
1 z

(1.20)

Funcin rampa unitaria: Dada la funcin rampa unitaria definida por

t, t 0
f (t )
0, t 0

(1.21)

la transformada z de f(t) viene dada por

k 0

k 0

k 0

F ( z ) Z [ f (t )] Z [t ] f (kT ) z k kTz k T kz k
1

T (z 2z
T

z
Tz

1 2
(1 z )
( z 1) 2

(1.22)

Funcin exponencial: Dada la funcin exponencial definida por

e at , t 0
f (t )
0, t 0

(1.23)

la transformada z de f(t) viene dada por

F ( z ) Z [ f (t )] Z [e at ] Z [e akT ] e akT z k
k 0

1 e

aT

z e

1
1 e

aT

z 1

2 aT

z )
z

z e aT

(1.24)

En la tabla 3.1 se presenta la transformada z de funciones muy tiles, y su correspondiente

transformada de Laplace.

______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

Tabla 1.1: Transformada de Laplace versus transformada z.

f(t)
(t)

F(s)
1

1(t) o (t)

e at

t2

1
s
1
sa
1
s2
2
s3

1 e at

f(kT) o f(k)
(k)
1,
k=0
0,
k0
1(k) o (k)

e akT
kT

F(z)
1

z
z 1
z
z e aT
Tz
( z 1) 2

(kT ) 2

T 2 z ( z 1)
( z 1)3

a
s ( s a)

1 e akT

z (1 e aT )
( z 1)( z e aT )

e at ebt

ba
( s a)( s b)

e akT ebkT

z (e aT ebT )
( z e aT )( z e bT )

t e at

1
( s a) 2

kTe akT

Tze aT
( z e aT ) 2

t 2 e at

2
( s a )3

(kT )2 e akT

T 2 ze aT ( z e aT )
( z e aT )3

10

(1 at ) e at

s
( s a)2

(1 akT ) e akT

11

sen t

senkT

12

cos t

s
s
2
s 2

z 2 z (1 aT )e aT
( z e aT ) 2
z senT
2
z 2 z cosT 1

cos kT

z 2 z cosT
z 2 2 z cosT 1

13

e at sen t

e akT senkT

e aT z senT
z 2 2e aT z cosT e 2 aT

e akT cos kT

z ( z e aT cosT )
z 2 2e aT z cosT e 2 aT

( s a)
2

14

e at cos t

sa
( s a) 2 2

1.4.2 Propiedades y teoremas de la transformada z

En la tabla 1.2 se presenta un grupo de propiedades y teoremas de la transformada z
unilateral.
______________________________________________________________________
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

Tabla 3.2: Propiedades y teoremas de la transformada z.

1
2
3
4

f(t) o f(k)
a f (t )
a f1 (t ) b f 2 (t )
f (t T ) o f (k 1)
f (t 2T )

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

f (k 2)

z 2 F ( z ) z 2 f (0) z f (1)

f (t kT )

z k F ( z ) z k f (0) z k 1 f (T ) zf (kT T )

f (t kT )

z k F (z )

f (n k )

z k F ( z) z k f (0) z k 1 f (1) zf (k 1)

f (n k )

z k F (z )

10

e at f (t )

11

a k f (k )

F ( ze aT )
z
F
a
z
F ( z)
z 1
F (1)

12

f (k )
k 0

13

f (k )

Z[f(t)] o Z[f(k)]
a F (z )
a F1 ( z ) b F2 ( z )
z F ( z ) z f (0)

k 0

14

k m f (k )

z F (z )
dz

1.4.3 Forma general para obtener la transformada Z de una funcin

laplaciana.
La transformada z de una funcin laplaciana se puede obtener mediante el mtodo de los
residuos, que a continuacin se detalla.
Dada una funcin:
B( s )
G( s)
(1.25)
A( s)
donde el grado del polinomio A(s) es mayor que el de B(s), y suponiendo que todas las
races de B(s) poseen parte real negativa, entonces:

______________________________________________________________________ 10
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________
P
G(s) z

G ( z ) Z [G ( s)] ( s ai )
z e Ts s ai
i 1
m 1
Q
1
d j
m j G(s) z

(
s

b
)

j
m j 1
z e Ts s b
j 1
(m j 1)! ds
j

(1.26)

donde:
P es el nmero de polos ai no repetidos de G(s) .
Q es el nmero de polos bj que se repiten con multiplicidad mj.

Ejemplo 1.2: Determinar la transformada z de la funcin:

G( s)

1
s ( s 1)
2

Solucin
Para este caso: P = 1, a1 = -1, Q = 1, b1 = 0, m1 = 2. Luego:

1
z
G ( z ) ( s 1) 2

s ( s 1) z e Ts s 1

1 d 2
1
z
z
z
zT

s 2
Ts
T
(2 1)! ds s ( s 1) z e s 0 z e
z 1 ( z 1) 2

z
z( z 1 T )

T
ze
( z 1) 2

(1.27)

1.5 La Transformada Z Inversa

As como en sistemas de control en tiempo continuo la transformada s cumple un papel
muy importante, idnticamente, la transformada z juega un papel muy importante en
sistemas de control en tiempo discreto.
La notacin para la transformada z inversa es Z-1. La transformada z inversa de G(z)
da como resultado una nica g(kT), pero no da una nica g(t), debido a que la transformada
z inversa slo obtiene la secuencia de tiempo que especifica los valores de g(t) en los
valores discretos de tiempo t = 0, T, 2T, ..., y no est definido en los otros tiempos. Esto
implica entonces, que est definido en instantes de tiempo de muestreo T, el cual puede ser
por ejemplo 1 seg., 0.1 seg., 2 seg. ,etc.

Z 1[G( z)] g (kT )

(1.28)

______________________________________________________________________ 11
M.Sc., Ing. Ral Benites Saravia

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERA ELCTRICA Y ELECTRNICA
Control Digital
________________________________________________________________________________________

Para obtener g(kT) se recomienda seguir cualquier de los siguientes mtodos:

a) Tabla de transformadas z.
b) Obtener la transformada z inversa sin usar la tabla de transformadas z. En este caso
se puede aplicar los mtodos de divisin directa, computacional, expansin en
fracciones parciales y la integral de inversin (vase [1],[3]).

______________________________________________________________________ 12
M.Sc., Ing. Ral Benites Saravia