Manual IngCtrlAuto
Manual IngCtrlAuto
Manual IngCtrlAuto
Tecnólogo en Mecatrónica
MEC-05 Ingenierı́a de Control Automático
Septiembre, 2023
Resumen
Debido a que los avances en la teorı́a y práctica del control automático apor-
tan los medios para obtener un desempeño óptimo de los sistemas dinámicos,
mejorar la productividad, aligerar la carga de muchas operaciones manua-
les repetitivas y rutinarias, ası́ como de otras actividades, casi todos los
ingenieros y cientı́ficos deben tener un buen conocimiento de este campo.
No se pretende hacer un texto enciclopédico que cubra todos los temas que
hoy en dı́a constituyen la disciplina de Ingenierı́a de Control Automático.
Ese objetivo estarı́a por encima de nuestras posibilidades y capacidades. Se
trata de una selección de temas en las que forman un papel preponderante
las experiencias acumuladas y formación del participante. En otras palabras,
se han seleccionado los temas más significativos para los objetivos del curso
y la exclusión de algún tema no significa necesariamente un juicio de valor
sobre el mismo.
Índice general
Índice de figuras V
Índice de tablas IX
i
ÍNDICE GENERAL
ii
ÍNDICE GENERAL
iii
ÍNDICE GENERAL
iv
Índice de figuras
v
ÍNDICE DE FIGURAS
vi
ÍNDICE DE FIGURAS
vii
ÍNDICE DE FIGURAS
viii
Índice de tablas
ix
GLOSARIO
x
1
Introducción a la Ingenierı́a de
Control
1.1. Introducción
1
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
del buque real en maniobras, y también se usan modelos estructurales para predecir el
funcionamiento en carga estática o dinámica de grandes construcciones como son los
puentes o los rascacielos.
Estos modelos a escala sirven para estudiar algunas caracterı́sticas muy concretas
de los sistemas reales. Por ejemplo, en el diseño de un avión, una de las etapas más
importantes es la realización de un modelo aerodinámico que permita contemplar sus
caracterı́sticas de vuelo en un túnel de viento, y no interesan en absoluto otros aspec-
tos particulares como por ejemplo la distribución interna del cableado eléctrico, o el
funcionamiento de los sistemas hidráulicos.
Aquı́ trataremos otro tipo de modelos llamados modelos matemáticos, que tam-
bién son versiones simplificadas de la realidad y que apoyándose en las leyes que gobier-
nan los sistemas fı́sicos, intentan describir y predecir el comportamiento de aspectos
particulares del sistema real. De la misma manera que en el modelo a escala del avión
fabricado para estudiar su comportamiento aerodinámico no es necesario incorporar los
sistemas hidráulicos, un modelo matemático elaborado para describir las deformacio-
nes en régimen elástico de una estructura no tiene por que incorporar informaciones
relativas al comportamiento térmico de la misma; el comportamiento térmico, si es
necesario estudiarlo, deberá hacerse utilizando otro modelo de aplicación de las leyes
fı́sicas correspondientes.
El primer paso para obtener un modelo matemático de un sistema fı́sico real es la
obtención de un modelo fı́sico y es un prerrequisito básico en el desarrollo de casi
todas las estrategias de control. En la práctica es difı́cil representar un sistema fı́sico
con toda su complejidad y por ello en el análisis y diseño de sistemas de control es
común realizar algunas idealizaciones que faciliten su estudio.
Establecido el modelo, se resuelven las ecuaciones obtenidas para diversas condicio-
nes de entrada, y el resultado representa la respuesta dinámica del sistema, esto
es, la predicción de su comportamiento que permitirá establecer posteriormente cual es
el sistema de control más adecuado.
La obtención de un modelo matemático presupone establecer un equilibrio adecua-
do entre simplicidad y precisión en los resultados del análisis. La precisión se puede
aumentar si se aumenta su complejidad, pero si no se necesita una precisión extre-
ma, es preferible obtener un modelo simplificado que sea adecuado al problema que se
pretende resolver, aunque sea a costa de ignorar ciertas propiedades fı́sicas inherentes
2
1.3 Sistemas Lineales Invariantes en el Tiempo y Causales
al sistema (por ejemplo el rozamiento, las no linealidades, los efectos de segundo or-
den...), dado que ello redundará en una simplificación de los cálculos (téngase en cuenta
que en algunos sistemas reales es necesario plantear y resolver cientos de ecuaciones
simultáneamente para describir su comportamiento).
Como regla general, cuando se soluciona un problema nuevo es conveniente desa-
rrollar primero un modelo simplificado para obtener una idea general de la solución, y
después uno más completo para realizar un análisis detallado.
Desde el punto de vista matemático la primera simplificación que vamos a realizar
para abordar el problema es suponer que los sistemas (en general) y la planta (en
particular) son lineales, invariantes con el tiempo y causales.
3
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
Supondremos también en lo sucesivo que además de lineales, los sistemas que nos in-
teresa estudiar pertenecen a una subclase de éstos denominada invariantes en el tiempo,
haciendo este término referencia a que los parámetros (a0 , a1 , a2 , · · · , b0 , b1 , · · ·) presen-
tes en la ecuación diferencial, se mantienen fijos a lo largo del tiempo y por tanto no
son funciones del tiempo. Conceptualmente, un sistema es invariante en el tiempo si
se cumple que manteniendo fijas las condiciones iniciales y la entrada, la respuesta
obtenida es la misma, independientemente del instante de tiempo considerado.
Para el caso más simple de que la planta esté gobernada por una única ecuación
diferencial, exista una única excitación r(t) y una única variable independiente (que
habitualmente suele ser el tiempo (t)) se tiene que el sistema responderá a una ecuación
diferencial de este tipo:
d d2 d d2
a0 y(t) + a1 y(t) + a2 2 y(t) + · · · = b0 r(t) + b1 r(t) + b2 2 r(t) + · · · (1.2)
dt dt dt dt
4
1.3 Sistemas Lineales Invariantes en el Tiempo y Causales
Superposición: Si una entrada r1 (t) produce una salida y1 (t), y una entrada r2 (t)
produce una salida y2 (t), entonces la entrada r1 (t) + r2 (t) dará como resultado
una salida y1 (t) + y2 (t). Expresado matemáticamente:
Las funciones lineales de una variable tienen la expresión y = kx, donde k es una
constante. De modo similar, las funciones lineales de varias variables tienen la forma:
X
y = f (x1 , x2 , ...) = a1 x1 + a2 x2 + ... = ai xi (1.5)
f (x1 + x2 ) = x1 + x2 + 4
f (x1 ) + f (x2 ) = x1 + x2 + 8
f (x1 + x2 ) ̸= f (x1 ) + f (x2 )
5
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
6
1.3 Sistemas Lineales Invariantes en el Tiempo y Causales
7
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
8
1.4 La Transformada de Laplace
Linealidad: Si L[f1 (t)] = F1 (s) y L[f2 (t)] = F2 (s), entonces L[af1 (t) + bf2 (t)] =
aF1 (s) + bF2 (s). Esta propiedad indica que la transformada de una suma es la
suma de transformadas y se mantiene la multiplicación por una constante en la
transformación.
Traslación en el dominio del tiempo: Si L[f (t)] = F (s), entonces L[f (t−a)] =
e−as F (s). Ası́, si se retarda una función en a segundos, el retardo aparece reflejado
en la transformada como una multiplicación por la función e−as .
Teorema del valor final: Si L[f (t)] = F (s) y todos los polos de sF(s) están
en el semiplano izquierdo del plano s, entonces se verifica que: lı́mt→∞ f (t) =
lı́ms→0 sF (s). La aplicabilidad del teorema debe considerarse con cuidado, por-
que la evaluación de sF (s) en s = 0 en ciertas circunstancias proporciona un valor
9
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
Teorema del valor inicial: Si L[f (t)] = F (s) y f (0) existe, entonces se verifica
que lı́mt→0 f (t) = lı́ms→∞ sF (s). Este teorema nos permite encontrar el valor de
f(t) en 0+ directamente a partir de la transformada de Laplace de f (t) sin nece-
sidad de tener que calcular la respuesta en el tiempo y tomar lı́mites. Obsérvese
que el teorema del valor inicial no proporciona el valor de f(t) en 0 sino en 0+ y
que para aplicar este teorema no es necesario comprobar previamente la posición
de los polos de sF (s).
10
1.4 La Transformada de Laplace
11
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
Si F (s) tiene únicamente polos reales de primer orden, F (s) puede ponerse como:
N (s) k1 k2 kn
F (s) = = + + ··· + (1.12)
(s + p1 )(s + p2 ) · · · (s + pn ) s + p1 s + p2 s + pn
kj = (s + pj ) F (s)s=−pj (1.13)
Una vez conocidos los residuos de las fracciones, la L−1 es inmediata, ya que cada
sumando representará un término exponencial:
16s + 16
F (s) = (1.14)
s3 + 6s2 + 8s
16s + 16 16(s + 1) A B C
F (s) = = = + + (1.15)
s3 2
+ 6s + 8s s(s + 2)(s + 4) s s+2 s+4
A = [F (s).s]s=0 = 2
B = [F (s).(s + 2)]s=−2 = 4
C = [F (s).(s + 4)]s=−4 = −6
2 4 −6
F (s) = + + (1.16)
s s+2 s+4
y por tanto la L−1 :
f (t) = 2u(t) + 4e−2t u(t) − 6e−4t u(t) = (2 + 4e−2t − 6e−4t )u(t) (1.17)
12
1.4 La Transformada de Laplace
Supóngase ahora que F (s) tuviera en pj un polo de orden r (es decir, r raı́ces
repetidas). El desarrollo de este tipo de término serı́a de la forma:
k1 k2 kr
F (s) = r
+ r−1
+ ··· + (1.18)
(s + pj ) (s + pj ) s + pj
Al hacer s = −pj se anulan todos los términos del segundo miembro excepto k1 ,
que se puede evaluar, resultando:
k1 = (s + pj )r F (s)s=−pj (1.20)
1 d2
k3 = [(s + pj )r F (s)]s=−pj (1.23)
2! ds2
Para el término general kj se llega a la siguiente expresión general:
1 dj−1
kj = [(s + pj )r F (s)]s=−pj (1.24)
(j − 1)! dsj−1
Una vez conocidos los residuos de este tipo de polo, la transformada inversa se puede
obtener ası́:
−1 k1 k2 kr
f (t) = L r
+ r−1
+ ··· +
(s + p1 ) (s + p2 ) (s + pj )
k1 tr−1 k2 tr−2 kj tr−j
= + + ··· + + · · · + kr e−pj t
(r − 1)! (r − 2)! (r − j)!
Ejemplo 2: Sea la función:
8s + 8
F (s) = (1.25)
s3 + 4s2 + 4s
13
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
d
B= F (s).(s + 2)2 s=−2 = −2 (1.29)
ds
y entonces:
2 2 4
F (s) = − + (1.30)
s s + 2 (s + 2)2
y la L−1 :
f (t) = 2u(t) − 2e−2t u(t) + 4te−2t u(t) = 2 − 2e−2t + 4te−2t e(t) (1.31)
Los dos métodos anteriores son válidos para cualquier tipo de polos, bien sean
reales o complejos, sin embargo, si existen polos complejos, los residuos en general
serán complejos y puede hacerse una simplificación debido a que siempre que aparecen
raı́ces complejas en una ecuación lo hacen en forma de pares conjugados. En este caso
resulta más cómodo el desarrollo en fracciones parciales de la forma:
1 k1 s + k2
F (s) = = 2 (1.32)
s2 + as + b s + as + b
e identificar los coeficientes k1 y k2 .
Ejemplo 3: Sea la función:
16s + 26
F (s) = (1.33)
s (s2+ 4s + 13)
que tiene tres raı́ces: una real y dos complejas conjugadas, por tanto descomponemos
ası́:
A Bs + C
F (s) = + 2 (1.34)
s s + 4s + 13
14
1.4 La Transformada de Laplace
2 −2s + 8
F (s) = + 2 (1.37)
s s + 4s + 13
para calcular la L−1 hay que hacer alguna transformación en el segundo término para
asimilarlo a una suma de funciones seno y coseno: primero buscamos un cuadrado
perfecto en el denominador.
2 −2s + 8 2 −2s + 8
F (s) = + = + (1.38)
s s2 + 4s + 13 s (s + 2)2 + 9
ahora sacamos factor común s + 2 en el numerador:
2 (s + 2)(−2) + 12
F (s) = + (1.39)
s (s + 2)2 + 9
2 s+2 12 2 (s + 2) 3
F (s) = −2 + = −2 +4 (1.40)
s (s + 2)2 + 9 (s + 2)2 + 9 s (s + 2)2 + 9 (s + 2)2 + 9
y la L−1 queda: f (t) = 2 − 2e−2t cos 3t + 4e−2t sin 3t u(t).
Una de las formas más sencillas de resolver una ecuación diferencial de coeficien-
tes constantes es la utilización de la transformada de Laplace. Las razones son dos:
la primera es que mediante la transformación, los operadores derivada e integral se
convierten en productos y cocientes. La segunda es que no tenemos que preocuparnos
por la forma de la función empleada como entrada o excitación (habrı́a que tenerlo en
cuenta para resolver la ecuación por los métodos clásicos).
15
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
Para obtener la evolución de la salida del sistema y(t) al aplicar una entrada deter-
minada r(t) con condiciones iniciales dadas, necesitamos resolver la ecuación diferencial
lineal:
dy n drm
an + · · · + a0 y = bm + · · · + b0 r (1.41)
dtn dtm
Los pasos a seguir para su resolución aplicando la transformada de Laplace son los
siguientes:
Como los sistemas que nos ocupan son reales (causales, realizables), la entrada r(t)
es causal, ası́ que:
r(t) = 0; ∀t < 0 (1.44)
Esto quiere decir que, sea r(t) continua o discontinua en t = 0 siempre se cumple
que:
r(0) = 0 (1.45)
dr
=0 (1.46)
dt t=0
y ası́ sucesivamente
drm−1
=0 (1.47)
dtm−1 t=0
con lo que
an sn Y (s) − sn−1 y(0) + · · · + y n−1 (0) + · · · + a0 Y (s) = bm sm R(s) + · · · + b0 R(s)
(1.48)
16
1.4 La Transformada de Laplace
Si las condiciones iniciales son nulas y la entrada es no nula nos queda únicamente
la respuesta forzada, es decir, la salida cuando el sistema está en equilibrio y se le aplica
una determinada entrada (por ejemplo si a un péndulo se le desplaza bruscamente de su
posición de equilibrio). Si por el contrario la entrada es nula y las condiciones iniciales
son no nulas, la respuesta es la respuesta libre del sistema, es decir, como se comporta
el sistema de forma autónoma debido a unas determinadas condiciones iniciales (por
ejemplo si se suelta un péndulo desde un punto distinto a la vertical).
En resumen, la respuesta de un sistema es la suma de dos términos: su respuesta
libre ante condiciones iniciales no nulas y su respuesta forzada ante condiciones iniciales
nulas.
Ejemplo:
Encontrar la solución para la siguiente Ecuación Diferencial.
d2 x(t) dx(t)
+5 + 6x(t) = 0
dt2 dt
x (0) = 0
dx
(0) = 7
dt
dx
s2 X(s) − sx(0) − (0) + 5[sX(s) − x(0)] + 6X(s) = 0
dt
s+5 1 dx
X(s) = x(0) + 2 (0)
s2 + 5s + 6 s + 5s + 6 dt
7 7 7
X(s) = = −
(s + 2)(s + 3) s+2 s+3
17
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
El segundo sumando (respuesta libre) es nulo si las condiciones iniciales son nulas.
En muchas situaciones prácticas, en t = 0, antes de aplicar entrada, el sistema está en
equilibrio, esto se traduce en condiciones iniciales nulas, y en lo sucesivo, salvo que se
diga lo contrario, se supondrán condiciones iniciales nulas.
Ante condiciones iniciales nulas, la solución es únicamente el primer sumando, esto
es, la solución forzada. Al cociente de polinomios que multiplica a R(s) se le denomina
función de transferencia del sistema, y al segundo término (si existiera), término
de condiciones iniciales.
Por tanto, se denomina función de transferencia G(s) al cociente entre la transfor-
mada de Laplace de la salida, y la transformada de Laplace de la entrada de un sistema
ante condiciones iniciales nulas
Y (s)
G(s) = (1.52)
R(s) c.i.=0
L [y(t)] Y (s)
G(s) = = (1.53)
L [r(t)] R(s)
La función de transferencia tiene las siguientes caracterı́sticas:
18
1.5 Concepto de Función de Transferencia.
Permite caracterizar el comportamiento del sistema ante señales tı́picas sin necesi-
dad de obtener analı́ticamente la salida del mismo. Para ello, la única información
relevante son los polos de la función de transferencia, es decir, las raı́ces de su
denominador: en efecto, ellas determinarán el tipo de fracciones simples en que
podrá descomponerse y, por tanto, el tipo de funciones (exponenciales o/y trigo-
nométricas) de que constará la referida respuesta.
Ejemplo:
dy
Encontrar y(t) cuando x(t) = u(t), y(0) = 0 y dt = 0.
d2 y(t) dy(t)
2
+4 + 3y(t) = 3x(t)
dt dt
19
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
3A = 3; A = 1
−3
B(−1)(2) = 3; B =
2
1
C(−3)(−2) = 3; C =
2
1 3/2 1/2
Y (s) = − +
s s+1 s+3
3 1
y(t) = u(t)(1 − e−t + e−3t )
2 2
u(t) = ejwt
Z ∞
y(t) = h(τ )ejw(t−τ ) dτ
Z−∞
∞
= h(τ )e−jwτ dτ ejwt
−∞
u (t) = e−jwt
Z ∞
y(t) = h(τ )e−jw(t−τ ) dτ
−∞
Z ∞
= h(τ )ejwτ dτ e−jwt
−∞
20
1.6 Respuesta en Frecuencia de Sistemas Lineales Continuos
Z ∞
λ(−w) = h(τ )e−(−jwτ ) dτ
−∞
Z ∞
= h(τ )e−jwτ dτ
−∞
= λ(w)
λ(w) = Aejϕ
λ(w) = λ(−w) = Ae−jϕ
Resumen
21
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
Ejemplo:
Encontar la respuesta del sistema F(s) a la señal de entrada
s
F (s) =
(s + 0,7)(s + 0,8)
x(t) = 3 sen 50t + 15 cos(100t + 30)
Para w1 = 50rad/seg.
j50 ∼ 1
|F (jw1 )| = =
(j50 + 0,7)(j50 + 0,8) 50
arg [F (jw1 )] = arg(j50) − arg(j50 + 0,7) − arg(j50 + 0,8)
∼
= 90 − 90 − 90
∼
= −90
Para w2 = 100rad/seg
j100 ∼ 1
|F (jw2 )| = =
(j100 + 0,7)(j100 + 0,8) 100
arg [F (jw2 )] = arg(j100) − arg(j100 + 0,7) − arg(j100 + 0,8)
∼
= 90 − 90 − 90
∼
= −90
ÿ + 2ẏ + y = 3u̇ + 2u
22
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
A continuación estudiaremos los modelos, las leyes fı́sicas y las funciones de transfe-
rencia de los elementos más sencillos y más ampliamente utilizados en ingenierı́a. Estos
sistemas sencillos pueden reducirse a componentes básicos que son versiones matemáti-
cas idealizadas de elementos reales tales como resistencias, válvulas, masas, muelles,
etc. Como si se tratara de un kit de construcción, a partir de estos elementos simples
puede llegarse a modelar otros de mayor complejidad. Por ejemplo, a partir del cono-
cimiento de las relaciones entrada-salida de resistencias, capacidades e inductancias, y
combinando sus efectos puede modelarse el comportamiento de un motor.
Las plantas (y los sistemas de control completos) a modelar son muy variados y
abarcan campos como la mecánica, electricidad, electrónica, hidráulica, etc. Para rea-
lizar un estudio global y coherente es conveniente elegir para cada elemento básico un
par de variables que permita manejar con facilidad sistemas integrados por tecnologı́as
diferentes. Por ejemplo, un sistema de control puede incluir elementos mecánicos y
eléctricos al mismo tiempo (un motor eléctrico con una carga acoplada). Aquı́ emplea-
remos las variables de potencia, llamadas ası́ porque el producto de ambas resulta,
desde un punto de vista dimensional, una potencia.
23
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
En mecánica el par de variables que usamos son la fuerza y la velocidad para los
sistemas de traslación y el par y la velocidad angular para los de rotación; para los
sistemas eléctricos el voltaje y la corriente y para los hidráulicos la presión y el caudal.
En todos los sistemas que vamos a estudiar, los elementos idealizados que los forman
son únicamente de tres tipos: de carácter resistivo, inductivo o capacitivo, haciendo esta
clasificación referencia al tipo de relación matemática existente entre las dos variables
elegidas.
En los elementos de carácter resistivo, o lo que es lo mismo, elementos disipati-
vos (rozamiento mecánico y viscoso, resistencia eléctrica, y resistencia hidráulica), la
relación entre las variables fuerza y velocidad es una constante. En los inductivos se
almacena energı́a (masa mecánica, inductancia eléctrica, inertancia hidráulica), la re-
lación entre las variables de potencia es una derivada en el tiempo. En los capacitivos
(muelle, condensador eléctrico, capacidad hidráulica), la relación entre las variables es
una integral.
Por otra parte, las distintas fuentes de energı́a empleadas son: la gravedad en los sis-
temas mecánicos, los generadores de intensidad de corriente y de voltaje en los sistemas
eléctricos, y las bombas en los hidráulicos.
24
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
Los sistemas mecánicos están gobernados por la segunda ley de Newton, la cual
establece para sistemas mecánicos de traslación que “la suma vectorial de las fuerzas
exteriores aplicadas sobre un cuerpo de masa m, es igual al producto de su masa por
su aceleración ”. Cuando se trata de sistemas mecánicos de rotación la segunda ley de
Newton establece de forma dual que “la suma vectorial de momentos exteriores aplica-
dos sobre un cuerpo es igual a su momento de inercia multiplicado por su aceleración
angular ”.
X X
F ext = ma y Mext = Jα (1.54)
25
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
una cierta constante elástica k lo que les confiere su comportamiento caracterı́stico. Los
muelles de torsión también se construyen en acero con su correspondiente constante k.
El tercer elemento es el amortiguador que modela el rozamiento viscoso. Está com-
puesto por un cilindro hueco dentro del cual se sitúa un pistón móvil. La fricción viscosa
es el rozamiento entre superficies móviles separadas por un fluido viscoso o también se
define como la fuerza entre un cuerpo sólido y un medio fluido. El cilindro está lleno
con un fluido (habitualmente aire comprimido o aceite) y dispone de un camino que le
permite retornar al lado opuesto del pistón a través del pequeño espacio existente entre
el pistón y el cilindro o bien mediante orificios hechos en el pistón.
La diferencia entre el amortiguador de traslación y el de rotación es la dirección del
movimiento (traslación y giro) y la resistencia al mismo en un medio viscoso (fuerza
y par). Los amortiguadores producen una fuerza (o un par) que es proporcional a la
velocidad lineal (o angular), y la fuerza resistente (o el par resistente), idealmente,
es proporcional a la velocidad lineal (o angular). Estos dispositivos no son capaces
de almacenar energı́a y sólo pueden disiparla en forma de calor (son elementos de
carácter resistivo). Ambos tipos de amortiguadores se caracterizan por un coeficiente b
(coeficiente de rozamiento viscoso).
El amortiguador en movimientos de rotación representa el rozamiento que impide
el giro. Por ejemplo el movimiento de un disco en el aire, que ofrece una resistencia
aerodinámica al movimiento de rotación del disco; o el giro del eje de un motor dentro
de un cojinete de apoyo.
Ejemplo 1: Sistema masa-muelle-amortiguador (traslación).
Combinemos la acción de los tres elementos básicos en un sistema de traslación como
el de la figura 1.3 y calculemos, por ejemplo, la función de transferencia X(s)/F (s).
Antes de nada es necesario especificar el sentido y dirección tanto de las fuerzas que
actúan como de los desplazamientos (lo mejor es dibujar el sentido real de las fuerzas
y los desplazamientos y especificar su módulo de manera que quede positivo). Después
se aı́sla el cuerpo que tiene masa (el resorte y el amortiguador son elementos ideales sin
masa) y se plantea la segunda ley de Newton sobre él. Obsérvese que algunas fuerzas
son del tipo acción-reacción por lo que tienen el mismo módulo y sentido contrario.
Como el modelo fı́sico incluye un muelle en posición vertical, es conveniente elegir
una referencia de desplazamiento que coincida con la posición relajada del resorte. El
efecto de fuerza de gravedad (mg) hará que el resorte, en su posición relajada esté algo
26
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
k
b
kx(t) bdx(t)/dt
0 0
m m
x(t) x(t)
f (t) f (t)
extendido. Entonces, para simplificar los cálculos, conviene tomar como referencia de
desplazamiento esa posición de reposo con masa incluida.
Dibujamos el diagrama de cuerpo libre (aislamos la masa y dibujamos todas las
fuerzas que actúan sobre ella con su sentido) y planteamos la ecuación dinámica del
movimiento, resultando:
dx(t) dx(t)2
f (t) − kx(t) − b =m (1.55)
dt dt2
Tomando L:
X(s) 1
F (s) − kX(s) − bsX(s) = ms2 X(s) =⇒ = 2
(1.56)
F (s) ms + bs + k
b
f
m1 m2
k1 k2
27
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
b
f
m1 m2
k1 k2
x1 x2
m1 m2
k1 (x1 ) k2 (x1 − x2 ) k2 (x1 − x2 )
x1 x2
d(x1 − x2 ) d2 x1
f − k l xl − b − k2 (x1 − x2 ) = m1 (1.57)
dt dt
d(x1 − x2 ) d2 x2
b + k2 (x1 − x2 ) = m2 (1.58)
dt dt
Tomando transformada de Laplace en ambas ecuaciones:
F (s) − k1 X1 (s) + bs [X1 (s) − X2 (s)] + k2 [X1 (s) − X2 (s)] = m1 s2 X1 (s) (1.59)
28
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
X1 (s) m2 s2 + bs + k2
= (1.61)
F (s) (m1 s2 + bs + k + k2 ) (m2 s2 + bs + k2 ) − (bs + k2 )2
Ejemplo 3: Sistema masa-amortiguador (rotación).
Supongamos el sistema de la figura 1.7 donde un disco de inercia J gira en un medio
viscoso de coeficiente b. Al aplicar un par T (t) al sistema, se obtiene un desplazamiento
θ(t) del eje. Calculemos la función de transferencia θ(s)/T (s).
T(t) θ(t)
b
J
dθ(t) d2 θ(t)
T (t) − b =J (1.62)
dt dt2
θ(s) 1
G(s) = = (1.63)
T (s) s(sJ + b)
En los circuitos eléctricos hay tres elementos básicos: resistencia eléctrica, induc-
tancia y condensador, y dos variables de potencia: voltaje e intensidad de corriente.
Si al suministrar energı́a al elemento, éste la disipa en forma de campo magnético
es una bobina, y si lo hace en forma de campo eléctrico se trata de un condensador.
En la práctica, los componentes reales de un circuito se comportan de más de una de
dichas formas, y muchas veces de las tres simultáneamente, pero lo normal es que pre-
domine uno de los tres efectos citados sobre los otros dos. Por ejemplo, para construir
una bobina (elemento real) es necesario arrollar un hilo conductor en espiral, lo cual
dará a la bobina la inductancia L que le es caracterı́stica, pero a la vez; el hilo conduc-
tor presentará una cierta resistencia R al paso de la corriente, lo cual será un efecto
secundario que sólo se tendrá en cuenta si se quiere construir un modelo muy exacto
29
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
de la bobina real. Esto también ocurre con los elementos mecánicos reales: un muelle
tiene una masa, pero no la consideramos porque el efecto que nos interesa estudiar es
su elasticidad y no su inercia.
La ley fı́sica que gobierna el funcionamiento de los elementos resistivos es la ley de
Ohm: “El voltaje v(t) en los terminales de un elemento resistivo puro es directamente
proporcional a la corriente i(t) que circula por él”. La constante de proporcionalidad
R se llama resistencia eléctrica.
La ley fı́sica que gobierna el funcionamiento de una inductancia es: “El voltaje v(t)
que se origina en los terminales de una bobina cuando circula por ella una corriente
variable con el tiempo i(t), es proporcional a dicha variación”. La constante de propor-
cionalidad se denomina inductancia L.
La ley fı́sica que gobierna el funcionamiento de un condensador eléctrico es: “El
voltaje v(t) que se origina en los terminales de un condensador cuando circula por el
una corriente i(t) es proporcional a la integral de dicha corriente”. La constante de
proporcionalidad es el inverso de la capacidad C.
Si queremos calcular funciones de transferencia de circuitos eléctricos en los cuales
existan combinaciones más o menos complejas de los tres elementos estudiados necesi-
taremos el auxilio de las leyes de Kirchhoff, las cuales gobiernan los circuitos eléctricos.
La primera enuncia que “la suma de las corrientes que confluyen en un nudo es igual
a cero”. La segunda dice que “la suma de los voltajes a lo largo de un trayectoria
cerrada o malla es cero”. Ambas leyes pueden ser usadas de manera combinada para
determinar el conjunto de ecuaciones integro-diferenciales necesarias para predecir el
comportamiento de un sistema eléctrico. La elección de una ley u otra para resolver un
circuito es únicamente una cuestión de comodidad, pues cualquiera de ellas conduce a
la misma solución.
30
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
R L
+ bc b bc +
C
i(t)
e(t) vC (t)
− bc b bc −
y teniendo en cuenta que la corriente i(t) es la misma: vR (t) = Ri(t); vL (t) = L di(t)
dt y
R
1 ∞ di(t) R
1 ∞
vC (t) = C 0 i(t)dt sustituyendo: e(t) = Ri(t) + L dt + C 0 i(t)dt
Tomando L:
I(s)
E(s) = RI(s) + LsI(s) + (1.65)
Cs
31
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
I(s)
por otra parte: V c(s) = Cs sustituyendo I(s) en la ecuación anterior:
1
E(s) = CsV c(s) R + Ls + (1.66)
Cs
de donde
VC (s) 1 1
= = 2
(1.67)
E(s) Cs [R + Ls + 1/Cs] LCs + RCs + 1
Ejemplo 2: Sistema RC-RC.
Calculemos la función de transferencia V (s)/E(s) en el sistema eléctrico de la figura
1.9. Como la estructura es de tipo mallas, es más cómodo utilizar la segunda ley de
Kirchoff.
R1 R2
+ bc b b bc +
C1 C2
e(t) v(t)
i1 (t) i2 (t)
- bc b b bc -
32
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
Los sistemas que involucran mecánica de fluidos están entre los más difı́ciles de
modelar y analizar en cuanto a su comportamiento dinámico. Aquı́ nos quedaremos
con una subclase denominada sistemas hidráulicos, caracterizados por trabajar con altas
presiones, bajos caudales y fluidos incompresibles (aunque se incremente la presión, el
volumen se mantiene constante).
La dinámica de los sistemas hidráulicos puede representarse mediante ecuaciones
diferenciales lineales, sólo si el fluido es incompresible y el régimen de flujo es laminar
(el fluido se desplaza de forma ordenada). La mayorı́a de los sistemas hidráulicos reales
no se comportan ası́, sino que se comprimen si se les somete a presión y el régimen de
flujo es turbulento (caracterizado porque el desplazamiento del fluido consiste en un
movimiento aleatorio e irregular superpuesto al movimiento principal).
Regimen Laminar
Regimen turbulento
33
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
P2 − P1
R= (1.75)
q
siendo p la presión (a uno y otro lado del elemento resistivo) y q el caudal que lo
atraviesa. La capacitancia hidráulica es el término que se emplea para describir el
almacenamiento de energı́a en el lı́quido en forma de energı́a potencial. Se define como:
Z
1
R= (q1 − q2 )dt (1.76)
p
siendo p la presión a la que está sometida el depósito y q el caudal (de entrada y salida).
La inertancia está ligada con la inercia del fluido, esto, es su capacidad para permanecer
inmóvil o para mantener su movimiento. Consideremos un flujo no estacionario, es decir,
que no ha alcanzado el régimen permanente, de un fluido incompresible. La fuerza f
f q(t)
necesaria para acelerar el fluido en la tuberı́a serı́a: f = ma, como p = A y v(t) = A ,
entonces p1 A − p2 A = m dv(t) dq(t)
dt = ρl dt por tanto:
ρl dv(t) dq(t)
p1 − p2 = =I (1.77)
A dt dt
es decir, que la pérdida o caı́da de presión es proporcional a la derivada del caudal
respecto al tiempo, donde I representa la inertancia del fluido. Si las tuberı́as son
cortas y anchas este efecto (inertancia) no suele tenerse en cuenta.
Ejemplo 1: Sistema tanque en régimen laminar.
Calculemos la función de transferencia Qs (s)/Qe (s) del sistema de la figura 1.11. En
los sistemas de nivel de lı́quido la condición de equilibrio se enuncia ası́: “La diferencia
entre el caudal de entrada y el caudal de salida por unidad de tiempo es igual a la
cantidad de lı́quido acumulado ”.
Si en régimen permanente entra en el depósito una cierta cantidad de lı́quido por
unidad de tiempo (caudal en m3 /s) y qe es el incremento respecto a ese caudal de
entrada en equilibrio, a la salida se dispondrá de otro caudal en equilibrio más un
incremento qs correspondiente a la variación. Lo que interesa es la dinámica, es decir,
las variaciones sobre la situación de equilibrio.
34
1.7 Sistemas Fı́sicos: Modelos, Ecuaciones y Funciones de Transferencia.
Por otra parte, el caudal de salida qs dependerá del grado de apertura de la válvula
de salida (efecto que contempla el parámetro R resistencia hidráulica, ya que las válvu-
las se modelan como elementos resistivos: la resistencia al paso del flujo está controlada
por la posición del mando de la válvula), y la diferencia entre el caudal de entrada y el
de salida provocará un incremento en la cantidad de lı́quido acumulado en el tanque.
qe
h0 + ∆h
qs
35
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
1 dh(t)
qe (t) − h(t) = A (1.83)
R dt
1 H(s) R
Qe (s) − H(s) = AsH(s) ⇒ = (1.84)
R Qe (s) sAR + 1
Ejemplo 2: Sistema tanque en régimen turbulento.
Para el mismo sistema tanque del ejemplo anterior, calculemos Qs(s)/Qe(s). Plan-
teamos la ecuación de la dinámica del sistema:
dV (t) dh(t)
qe (t) − qs (t) = =A (1.85)
dt dt
dqs (t)
qe (t) − qs (t) = 2AB 2 qs (t) (1.87)
dt
q0
h
h0
Figura 1.12: Relación no lineal entre la altura y el caudal. - Sistemas hidráulicos
36
1.8 Diagramas de Bloques
En el entorno del punto de operación (h0 , qs0 ), la relación entre ambas variables
puede suponerse lineal si no nos alejamos mucho, por tanto, la ecuación que liga ambas
variables se puede aproximar por una recta tangente a la curva en el punto. La ecuación
de la recta tangente viene dada por:
dqs B
donde la pendiente m es la derivada de la función en el punto: m = dh q0 ,h0 = √
2 h0
.
Por tanto, la relación ahora es:
B
qs − qs0 = √ (h − h0 ) (1.89)
2 h0
lo cual nos lleva ahora a una ecuación diferencial lineal:
√
h0 dqs (t)
qe (t) − qs (t) = 2A (1.90)
B dt
Tomando transformadas de Laplace:
√
h0
Qe (s) − Qs (s) = 2A sQs (s (1.91)
B
quedando la función de transferencia:
Qs (s) 1
= √ (1.92)
Qe (s) 1 + 2AB h0 s
37
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
ellas de los cuales se ha despejado alguna variable), mostrando los diferentes modelos de
subsistemas y los caminos que recorren las señales, y proporcionando un conocimiento
rápido de las caracterı́sticas de los mismos al permitir observar como se interrelacionan
las diferentes partes.
Las flechas se usan para representar las direcciones en las que fluyen las señales.
Cuando las señales son funciones del tiempo se representan por letras minúsculas se-
guidas por (t), por ejemplo i(t), aunque con frecuencia (t) se omite cuando es obvio
que las señales son funciones del tiempo. Cuando las señales están en el dominio s se
representan con letras mayúsculas seguidas por (s), por ejemplo I(s).
Un sumador es un elemento en el que las señales se suman algebraicamente. Si una
de las señales que entran a dicho punto se indica como positiva y la otra como negativa,
entonces dicha suma es la diferencia entre las dos señales.
Si ambas señales se indican como positivas entonces la suma es la adición de las dos
señales. Cuando de algún punto de la trayectoria de la señal se toma una derivación de
una señal, se habla de punto de bifurcación y se representa igual que en los circuitos
eléctricos, mediante un punto.
Los bloques se dibujan con la función de transferencia en su interior.
V(s) k θ(s) b
s(sT +1)
+
-
El término trayectoria directa se usa para los elementos a través de los cuales
pasa la señal en la dirección entrada → salida a lo largo del sistema. Las funciones
de transferencia para los elementos en esta trayectoria directa en general se designan
mediante G o G(s). El término trayectoria de realimentación se usa para los elementos
para los cuales pasa la señal cuando se inyecta de regreso desde la salida hacia la
entrada. Las funciones de transferencia para los elementos de esa trayectoria por lo
común se designan por H o H(s). El término trayectoria de prealimentación se usa para
los elementos que están en paralelo con la trayectoria directa y a través de los cuales
la señal se mueve en la misma dirección, es decir, entrada → salida.
38
1.8 Diagramas de Bloques
G1 (s) G2 (s) b b
G1 (s)
H(s) G2 (s)
Vf Km θs
s(Rf +sLf )(sJ+ B)
O se puede utilizar otro en el que se detallen los subsistemas de que consta y las
relaciones entradas salida entre cada uno de ellos: campo, constante del motor, carga
del motor, conversión de velocidad angular a ángulo girado por el eje.
Vf 1
If Tm ωs θs
1 1
Rf +sLf Km sJ+B s
39
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
i(t)
ei (t) C eo (t)
− −
bc b bc
40
1.8 Diagramas de Bloques
Ei (s) I(s)
1/R
+ I(s) Eo (s)
- 1/Cs
Eo (s)
+
-
Eo (s)
41
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
G2 (s)
G1 (s)
+
Y(s)
X(s)
-
G2 (s)
42
1.9 Sistemas Lineales en Tiempo Discreto
tabla siguiente se resumen las reglas más frecuentes, la mayorı́a de ellas se demuestran
por simple inspección.
X(s) Y(s)b
G(s)
Eliminación de un lazo + X(s) G(s) Y(s)
∓ 1 ± G(s)H(s)
de realimentación H(s)
to de bifurcación des-
X(s) 1
X(s) G(s)
pués de un bloque
43
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
unitarios
∞
X
δT (t) = δ(t − kT ) (1.99)
k=0
Este muestreador ideal está definido como un muestreador que abre y cierra de
manera instantánea, en tiempo cero, cada T segundos.
La salida del muestreador ideal se expresa como
∞
X
f ∗ (t) = f (kT )δ(t − kT ) = f (t)δT (t) (1.100)
k=0
∞
X
yn = xk fn−k (1.102)
k=−∞
ó
∞
X
yn = xn−k fk (1.103)
k=−∞
44
1.10 Transformada Z
f (t)
1
(a)
7T t
0 T 2T 3T 4T 5T 6T
δT (t)
1
(b) b b b
t
0 T 2T 3T 4T 5T 6T 7T
∗
f (t)
f (2T )
1 f (T )
f (3T )
(c) f (0)
f (4T ) t
f (7T )
f (5T )
f (6T )
Figura 1.21: Muestreador ideal: (a) Señal de entrada; (b) Tren de impulsos
unitarios; (c) Señal de salida. - f ∗ (t) = f (t)δT (t).
xk = {0; 1↑ ; 1; 1; 1; · · ·} (1.105)
1.10. Transformada Z
Euler
∞
X
z zk
e = (1.106)
k!
k=0
45
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
n x(n) y(n)
−1 0 0
0 1 1
1 1 1,6
2 1 1,76
3 1 1,736
4 1 1,6896
5 1 1,6666
6 1 1,662
7 1 1,6639
8 1 1,6659
9 1 1,6668
10 1 1,6669
1.8
1.6
1.4
1.2
Amplitud
0.8
0.6 x(n)
y(n)
0.4
0.2
0
−2 0 2 4 6 8 10 12
Número de muestra [n]
46
1.10 Transformada Z
∞
X z (2k+1)
sin z = (−1)k (1.107)
(2k + 1)!
k=0
∞
X z 2k
cos z = (−1)k (1.108)
(2k)!
k=0
iz
e = cos z + i sin z (1.109)
∞
X
S= ak = lı́m Sn (1.113)
n→∞
k=0
n
X
Sn = ak = 1 + a + a2 + · · · + an (1.114)
k=0
n+1
X
aSn = ak = a + a2 + a3 + · · · + an + an+1 (1.115)
k=1
Sn (1 − a) = (1 − an+1 ) (1.118)
1− an+1
Sn = (1.119)
1−a
1
1−a si |a| < 1
1 − an+1
S = lı́m Sn = lı́m = (1.120)
n→∞ n→∞ 1 − a
indef inido si |a| ≥ 1
Ejemplo:
9 9 9
0,99 · · · = + + + ··· (1.121)
10 " 100 1000 #
2
9 1 1
= 1+ + + ··· (1.122)
10 10 10
∞
9 X 1 k
= (1.123)
10 10
k=0
47
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
!
9 1
= 1 (1.124)
10 1 − 10
9 1
= 9 (1.125)
10 10
9 10
= (1.126)
10 9
= 1 (1.127)
1.10.2. Definición
∞
X
F (z) = f (kT )z −k (1.128)
k=−∞
donde si
f (kT ) = 0 ∀ t < 0
∞
X
F (z) = f (kT )(z −1 )k (1.129)
k=0
Región de convergencia
( ∞
)
X
ROC = z ∈ C| xk ∥z∥−k < ∞ (1.130)
k=0
Ejemplos:
f (kT ) = 1
X∞
1
F (z) = (z −1 )k = 1
k=0
1− z
ROC
1
< 1
z
|z| > 1
48
1.10 Transformada Z
∞
X ∞ αT k
X
αkT −k a
F (z) = a z =
z
k=0 k=0
z
F (z) =
z − aαT
ROC
aαT
< 1
z
|z| > aαT
Operador lineal
49
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
Prueba de Convolución
∞
!
X
Z ({gd∗ xk }) = Z gk xn−k
k=−∞
∞ ∞
!
X X
= gk xn−k z −n
n=−∞ k=−∞
X∞ ∞
X
−k
= gk z xn−k z −(n−k)
k=−∞ n=−∞
= G(z)X(z)
∞
X
Z (xk−n ) = (xk−n ) z −k
k=−∞
∞
X
z −n x(k−n) z −(k−n) = z −n X(z)
k=−∞
1. Por Definición. I
1
xk = X(z)z k−1 dz
2πi C
50
1.10 Transformada Z
2. División Larga.
b0 z m + b1 z m−1 + · · · + bm−1 z + bm
F (z) =
z n + a1 z n−1 + · · · + an−1 z + an
b0 z m + b1 z m−1 + · · · + bm−1 z + bm
F (z) =
z n + a1 z n−1 + · · · + an−1 z + an
el sistema es causal si n ≥ m. Para efectos prácticos se considera m = n.
Ejemplo de Implementación:
Considerar el sistema discreto
0,003115z + 0,003033
F (z) =
z 2 − 1,917z + 0,9231
Simular la respuesta del sistema a una señal senoidal de amplitud 3 y frecuencia de
2Hz y frecuencia de muestreo de 100Hz.
La ecuación de diferencias es:
Programa en MATLAB
clc
clear all
fs=1000;
Ts=1/fs;
f=20;
N=4*fs/f;
x(1)=3*sin(2*pi*f*Ts*1);
51
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
x(2)=3*sin(2*pi*f*Ts*2);
y(1)=0;
y(2)=0.003115*x(1);
kTs(1)=1*Ts;
kTs(2)=2*Ts;
for k=3:N;
kTs(k)=k*Ts;
x(k)=3*sin(2*pi*f*Ts*k);
y(k) = 1.917*y(k-1) - 0.9231*y(k-2) + 0.003115*x(k-1) + 0.003033*x(k-2);
end
n1=N/2;
n2=N;
plot(kTs(n1:n2),x(n1:n2),’r’, kTs(n1:n2),y(n1:n2),’k’)
∞
X
y(k) = f (n)u(k − n) (1.131)
n=−∞
u(kT ) y(kT )
F (z)
52
1.11 Respuesta en Frecuencia de Sistemas Lineales Discretos
.
Ahora consideramos que la entrada al sistema como
∞
!
X
y(k) = F (n)ejw0 nT e−jw0 kT
n=−∞
u (kT ) = sen(w0 kT )
F (ejw0 T ) jw0 kT F (e−jw0 T ) −jw0 kT
y (kT ) = e − e
2j 2j
En sistemas Reales
r = F (ejw0 T ) = F (e−jw0 T )
ϕ = arg F (ejw0 T ) = −arg F (e−jw0 T )
53
1. INTRODUCCIÓN A LA INGENIERÍA DE CONTROL
x(kT ) y(kT )
F (z)
= r sen (w0 kT + ϕ)
Resumen
x(kT ) = A sen(w0 kT + ϕ)
y(kT ) = F ejw0 T A sen w0 kT + ϕ + arg F ejw0 T
⇒ Realizar Práctica 1
54
2
Métodos Clásicos de
Identificación de Procesos
Objetivo:
Al término del tema los participantes compararán los métodos clásicos de identifi-
cación de sistemas para obtener modelos de manera rápida siguiendo los lineamientos
indicados durante el estudio del tema, incluyendo por lo menos el 80
2.1. Introducción
2. Ajuste sobre la marcha de los parámetros del regulador sobre la base de una
identificación recursiva de los parámetros del modelo (Identificación on-line).
55
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
56
2.2 Métodos de Identificación
En el presente tema vamos a hacer una rápida revisión de algunos métodos clásicos
de identificación de procesos que han sido exitosamente utilizados en la práctica y
que aún hoy siguen siendo útiles en muchos casos, sobre todo por su simplicidad. El
contenido del capı́tulo está organizado en la forma siguiente: En 2.2 se discuten los
principales métodos o enfoques utilizados en la identificación de procesos. En 2.3 se
revisan las técnicas experimentales de identificación basadas en la respuesta a la señal
escalón que hasta nuestros dı́as sigue siendo la opción más utilizada en la práctica
industrial. También en esta sección se encuentran equivalencias útiles entre los modelos
continuos y discretos de primer y segundo orden. En 1.4 se incursiona por el uso de las
señales aleatorias para la identificación, no sin antes hacer un breve repaso de algunos
conceptos básicos de la dinámica estadı́stica como las funciones de auto correlación y
correlación cruzada.
Teniendo en cuenta los objetivos que se persiguen en la identificación, ası́ como los
recursos e información disponibles, se han utilizado distintos enfoques o métodos para
este propósito, siendo los principales: Identificación Analı́tica, Identificación Experimen-
tal utilizando señales especiales de prueba e Identificación Recursiva. A continuación
haremos una breve discusión de cada uno de ellos.
57
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
La señal más simple que puede utilizarse para la identificación es, sin duda, la
función escalón. Esta señal se aproxima, por ejemplo, mediante la apertura o cierre
súbitos de una válvula, un cambio rápido de voltaje o corriente que alimenta a algún
tipo de actuador eléctrico, etc. Un escalón ideal es una señal cuyo tiempo de crecimiento
inicial es cero. Resulta obvio que tal señal no puede ser realizada fı́sicamente, pues serı́a
necesaria una energı́a infinita, sin embargo, si el tiempo de crecimiento inicial es varias
veces más pequeño que el perı́odo correspondiente a la frecuencia máxima de interés,
el error que se introduce en la identificación puede considerarse despreciable.
58
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
Ke−θs
Gm (s) = (2.1)
T1 s + 1
Ke−θs
Gm (s) = (2.2)
(T1 s + 1)(T2 s + 1)
Ke−θs
Gm (s) = (2.3)
s2 /ωn2 + 2ξs/ωn + 1
59
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
KU
0.632KU
T +θ
−(t−θ)
y(t) = KU (1 − e T ) (2.4)
t−θ =T (2.5)
tenemos que:
y(t = θ + T ) = KU (1 − e−1 ) = 0,632KU (2.6)
Y est
K= (2.7)
U
donde: Yest - valor estacionario de la respuesta. U - amplitud del escalón.
60
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
KU
0.632KU
0.284KU
θ
T
θ+ 3
θ+T
t0,284 = θ + T /3 (2.8)
t0,632 = θ + T (2.9)
61
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
e−5s
G(s) = (2.10)
(5s + 1)(10s + 1)
0.9
0.8
0.7
0.632
0.6
0.5
0.4
0.3
0.284
0.2
0.1
0
0 10 20 30 40 50 60 70 80
12.6 20.8 t[seg]
Vamos a tratar de aproximar este sistema por uno de primero orden usando el
método descrito en la figura 2.2. Su respuesta ante una señal escalón de amplitud 5 se
muestra en la figura 2.3, en la que aparecen también los tiempos en que la respuesta
alcanza el 28.4 y el 63.2 % del nuevo valor estacionario. Es posible entonces escribir las
ecuaciones:
T1 + θ = 20,8 (2.11)
T1
+ θ = 12,6 (2.12)
3
de donde se obtienen los valores de los parámetros correspondientes al sistema de primer
orden equivalente:
T1 = 12,3s, θ = 8,5s (2.13)
62
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 10 20 30 40 50 60 70 80
t[seg]
Figura 2.4: Respuesta del sistema original vs sistema de primer orden aproxi-
mado - Ejemplo 1.
⇒ Realizar Práctica 2
63
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
Ti
θ Tb Tc
Ta
Figura 2.5: Sistemas de segundo orden sobre amortiguado - Respuesta al escalón
unitario.
64
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
y entonces:
θ = Ti + Tc − Ta − Tb (2.18)
KU
0.632KU
0.284KU
θ
θ + t2
θ + t6
65
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
K ∼ K
G(s) = = (2.19)
(T1 s + 1)(T2 s + 1) . . . (Tn s + 1) (τ s + 1)n
66
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
Td Tc
K
ηK
Q
Ta Tb
Tc
Figura 2.9: Sistema aperiódico de alto orden - Respuesta al escalón unitario
Tabla 2.1: Relaciones de n con los cocientes Ta /Tb y Te /Tb - Sistemas aperiódicos
n 1 2 3 4 5 6 7 8 9 10
Ta /Tb 0 0.104 0.218 0.319 0.410 0.493 0.570 0.642 0.709 0.773
Te /Tb 1 0.736 0.677 0.647 0.629 0.616 0.606 0.599 0.593 0.587
η 0 0.264 0.323 0.353 0.371 0.384 0.394 0.401 0.407 0.413
n 1 2 3 4 5 6 7 8 9 10
Ta /τ 0 0.282 0.805 1.425 2.1 1.811 3.549 4.307 5.081 5.869
Tb /τ 1 2.718 3.695 4.463 5.119 5.699 6.226 6.711 7.164 7.59
Td /τ 0 1 2 3 4 5 6 7 8 9
Te /τ 1 2 2.5 2.888 3.219 3.51 3.775 4.018 4.245 4.458
67
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
Ejemplo 2. En este ejemplo vamos a ilustrar el uso del método de Strejc descrito en
2.3.5 que permite identificar un sistema con varias constantes de tiempo mediante un
modelo del mismo orden, pero con una sola constante de tiempo repetida.
Consideremos, en efecto, un sistema descrito mediante la función transferencial:
1
G(s) = (2.20)
(3s + 1)(5s + 1)(8s + 1)(12s + 1)
68
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
ese punto se determinó evaluando la primera derivada del polinomio en el punto Q, con
lo que se obtiene su pendiente y también su intercepto.
La ecuación de la tangente obtenida para este ejemplo fue la siguiente:
que evaluada en y(t) = 0 y y(t) = 5 permite completar los parámetros requeridos para
utilizar las tablas 2.1 y 2.2 del método de Strejc y que se muestran en la figura 2.10.
Entrando en la tabla 2.1 con los valores calculados de los cocientes Ta /Tb y Te /Tb
ası́ como con el valor del parámetro η, puede determinarse que la mejor estimación
para el orden del modelo aproximado es 4. Al mismo tiempo, utilizando la tabla 2.2
con el valor de n = 4 y los cocientes Ta /τ , Tb /τ , Td /τ y Te /τ , se obtienen cuatro valores
posibles de τ , que se muestran en la figura 2.11. En esta figura se presentan también la
curva de respuesta del sistema original (lı́nea sólida) y las de los sistemas aproximados,
correspondientes a los cuatro valores distintos de τ . Evidentemente, el mejor ajuste se
obtiene para τ = 6,94 cuya respuesta (curva 4) prácticamente no se distingue de la del
sistema original.
69
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
u(t)
1
(a)
5T 6T 7T
t
0 T 2T 3T 4T
u∗ (t)
u(2T )
1 u(T )
u(3T )
(b) u(0)
u(4T ) t
u(7T )
u(5T )
m(t) u(6T )
1
(c)
5T 6T 7T
t
0 T 2T 3T 4T
Figura 2.13: Proceso de muestreo - (a) variable continua, (b) variable discreta, (c)
salida retenedor de orden cero.
70
2.3 Identificación Experimental Utilizando la Respuesta a la Señal Escalón
1 − e−Ts s
H(s) = (2.23)
s
La función transferencial discreta equivalente de un proceso continuo con función
transferencial continua F (s) viene dada por:
1 − e−Ts s
F (z −1 ) = Z( F (s)) (2.24)
s
La función transferencial (2.24) puede simplificarse si se considera un perı́odo de
discretización muy pequeño comparado con las constantes de tiempo del proceso. En
efecto, si consideramos el desarrollo en serie de la exponencial e−Ts s :
Ts2 2
e−Ts s = 1 − Ts s + s + ... (2.25)
2
y si se desprecian los términos de grado mayor que 1, tenemos que:
e−Ts s ≈ 1 − Ts s (2.26)
y
1 − (1 − Ts s)
H(s) ≈ = Ts (2.27)
s
Entonces la función transferencial discreta (2.24) se simplifica a:
z −P
F1 (z −1 ) = K1 (2.29)
1 − αz −1
donde:
71
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
θ
p = Round (2.30)
Ts
KTs
K1 = (2.32)
T1
y la ecuación de diferencias equivalente:
z −(p+1)
F2 (z −1 ) = K2 (2.34)
1 + a1 z −1 + a2 z −2
donde:
θ
p = Round (2.37)
Ts
a1 = −(α1 + α2 ) (2.38)
a2 = α1 α2 (2.39)
α1 − α2
K2 = KTs (2.40)
T1 − T2
Ecuación de diferencias equivalente:
72
2.4 Conclusiones.
z −(p+1)
F3 (z −1 ) = K3 (2.42)
1 + a1 z −1 + a2 z −2
donde:
θ
p = Round (2.43)
Ts
α = eξωn Ts (2.44)
p
q = ωn 1 − ξ 2 Ts (2.45)
Kωn Ts αsen(q)
K3 = p (2.46)
1 − ξ2
a1 = −2αcos(q) (2.47)
a2 = α2 (2.48)
Ke−θs
F (s) = (2.50)
(T1 s + 1)(T2 s + 1)
considerando el retenedor de orden cero sin simplificación.
2.4. Conclusiones.
73
2. MÉTODOS CLÁSICOS DE IDENTIFICACIÓN DE PROCESOS
Hasta aquı́ hemos obtenido varios modelos genéricos que son capaces de aproximar
casi cualquier sistema lineal con una exactitud aceptable utilizando la respuesta a la
señal escalón se usan modelos de primer y segundo orden con retardo y sus equivalentes
discretizados los cuales nos son de gran utilidad cuando no disponemos de la información
suficiente y el tiempo es crı́tico.
74
3
Identificación de Sistemas
Mediante el Método de Mı́nimos
Cuadrados
Objetivo:
Al finalizar el estudio del tema de identificación mediante mı́nimos cuadrados los
participantes representaran un sistema de control simulado usando el método estudiado
e incluyendo por lo menos el 80
3.1. Introducción
Se hará una exposición detallada del método de mı́nimos cuadrados (MMC) para
la identificación de sistemas discretos lineales. El MMC resulta sin dudas, uno de los
enfoques básicos y más utilizados en la teorı́a y la práctica recientes de la identificación
de procesos. Numerosos autores, tales como Eykhoff (1974), Peterka(1975) y muchos
otros han dedicado un esfuerzo considerable al estudio de este método y sus diversas
extensiones y modificaciones. Aparte de su indudable atracción intuitiva, el MMC posee
una serie de propiedades estadı́sticas muy convenientes y sobre todo, puede encontrarse
una forma recursiva de este método suficientemente simple. Esta última posibilidad ha
venido a ser la base de los llamados métodos de identificación en lı́nea (on-line) de
procesos, que han alcanzado gran popularidad como resultado del desarrollo de los
sistemas de control digital.
75
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
Supongan ahora que disponen de un horno que conocen sus caracterı́sticas dinámicas
y que construimos un controlador para que mantenga una temperatura que nosotros
deseamos, y que funciona satisfactoriamente en el verano, pero cuando llega el invierto
tenemos un error en el control de la temperatura que nosotros no deseamos. Lo que
aprenderemos en este capı́tulo será un método de identificación que nos permita ajustar
los cambios que sufren los sistemas por agentes externos y evitar problemas de control.
Este método es de gran utilidad ya que todos los sistemas o procesos cambian con el
tiempo y esos cambios deben ser registrados en el momento para realizar los ajustes
pertinentes.
En el tema se discutirá en detalle el método de mı́nimos cuadrados, tanto en su
versión convencional, como en varias modificaciones orientadas a mejorar la precisión
del MMC cuando los datos están contaminados con ruido. El capı́tulo esta organizado
de la siguiente manera: En 3.2 se introduce el modelo ARX en sus varias formas posibles,
que es la estructura básica utilizada para aplicar los métodos de identificación en lı́nea.
En 3.3 se aplica el método general de mı́nimos cuadrados al modelo ARX y se obtienen
expresiones para la estimación óptima de los parámetros ası́ como de la matriz de
covarianza de la estimación. En 3.4 se discute una técnica heurı́stica de extrema utilidad,
la llamada técnica del olvido exponencial. En 3.5 se deriva en detalle el algoritmo
recursivo mı́nimo cuadrado convencional y se discuten las bondades y limitaciones de
este método.
v(t)
+
u(t) + y(t)
P
76
3.2 Estructura del Modelo ARX
E{e(k)} = 0 (3.2)
77
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
es decir, e(k) se considera un proceso de ruido blanco, con media cero y varianza σ 2 ,
no correlacionado con la entrada o salida del proceso. El modelo (3.1) puede escribirse
también en la forma:
y(k) = ŷ(k) + e(k) (3.6)
donde:
n
X n
X
ŷ(k) = aj y(k − j) + bj u(k − j) (3.7)
j=1 j=0
A(z −1 ) = a1 z −1 + a2 z −2 + · · · + an z −n (3.10)
B(z −1 ) = b0 + b1 z −1 + b2 z −2 + · · · + bn z −n (3.11)
o también
B(z −1 ) 1
Y (z) = U (z) + E(z) (3.13)
1 − A(z −1 ) 1 − A(z −1 )
78
3.2 Estructura del Modelo ARX
B(z −1 )
Y (z) = U (z) + V (z) (3.14)
1 − A(z −1 )
donde:
1
V (z) = E(z) (3.15)
1 − A(z −1 )
También:
n
X
v(k) = ai v(k − i) + e(k) (3.16)
i=1
donde P T y z(k) son respectivamente los vectores de los parámetros y de las medi-
ciones. Ambos vectores pueden definirse de varias formas; una muy conveniente desde
el punto de vista de los algoritmos de identificación en lı́nea, es la siguiente:
P T = [b0 , a1 , b1 , · · · , an , bn ] (3.18)
z T (k) = [u(k), y(k − 1), u(k − 1), · · · , y(k − n), u(k − n)] (3.19)
79
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
con
Y T (t) = [y(1), y(2, ), · · · , y(t)] (3.21)
y
u(1) ···
u(2) u(t)
y(0) · · · y(t − 1)
y(1)
u(0) · · · u(t − 1)
u(1)
Z T (t) = .. .... .. (3.22)
. . . .
y(1 − n) y(2 − n) · · · y(t − n)
u(1 − n) u(2 − n) · · · u(t − n)
es la matriz de las mediciones, de dimensión (2n + 1) × t.
También:
E T (t) = [e(1), e(2), · · · , e(t)] (3.23)
donde:
y(k) - ν vector de las salidas, u(k) - µ vector de las entradas, Aj , Bj - matrices ν × ν y
ν × µ de los coeficientes de regresión, e(k) - ν vector del ruido.
Un modelo compacto equivalente a (3.24) puede escribirse en la misma forma que
el correspondiente a los sistemas monovariables. En efecto:
80
3.2 Estructura del Modelo ARX
es un vector de dimensión:
ρ = [µ + (µ + ν)n] (3.27)
P T = [B0 , A1 , B1 , · · · , An , Bn ] (3.28)
Por último, el modelo generalizado para t mediciones dado en (3.20), puede escri-
birse para el caso multivariable en la misma forma:
Nótese que Y (t) y E(t) son ahora matrices de dimensión t × ν y Z(t) es una matriz
t × ρ.
Antes de finalizar esta sección, quisiéramos hacer dos observaciones que considera-
mos de suma importancia. La primera se refiere al orden del modelo ARX considerado
en (3.1). Por orden del modelo entendemos al número entero n que aparece como lı́mite
superior de las sumatorias en la expresión (3.1). Por comodidad y en aras de la claridad
de la exposición hemos supuesto que el lı́mite de ambas sumatorias es el mismo, dejan-
do abierta la posibilidad de que algunos de los coeficientes aj o bj sean ceros. El caso
frecuente de sistemas con retraso por transporte queda incluido si se hacen iguales a
cero algunos de los primeros coeficientes bj , por ejemplo b0 , b1 , etc. No obstante lo an-
terior, en capı́tulos posteriores será conveniente utilizar modelos en los que los órdenes
de las dos sumatorias pueden ser distintos y también donde aparezca eventualmente el
retraso por transporte en forma explı́cita.
La segunda observación se refiere a la inclusión de b0 en el modelo que proponemos,
es decir, la admisión implı́cita de que el valor de la variable de entrada u(k) influye
sobre la salida en el mismo tiempo y(k).
Muchos autores utilizan un modelo en el que no se considera esa influencia y se
supone la existencia siempre, como mı́nimo, de un perı́odo de retardo. Desde nuestro
81
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
Tm Tm Tm
Consideremos el modelo ARX para sistemas de una entrada y una salida, corres-
pondiente a t mediciones:
donde P, Y (t), Z(t) y E(t) se definen en (3.18), (3.21), (3.22) y (3.23). El vector de las
salidas según (3.31) puede considerarse como un vector aleatorio que tiene las siguientes
propiedades estadı́sticas:
E Y T (t) = P T Z T (t) (3.32)
82
3.3 Método de Mı́nimos Cuadrados para la Estimación de Parámetros del
Modelo ARX
donde W es una matriz definida positiva arbitraria llamada matriz de peso, que puede
asumir distintas formas especiales, dando lugar a diversas modificaciones del MMC.
Nótese que para darle mayor generalidad al problema estamos considerando que los
parámetros pueden variar con el tiempo. El caso más simple y que será objeto de
nuestro interés inmediato, consiste en considerar W = I. En este caso, el criterio J
puede también escribirse:
t
X
J = E T (t)E(t) = e(k)2 (3.35)
k=1
J = Y T (t)Y (t) − P T (t)Z T (t)Y (t) − Y T (t)Z(t)P (t) + P T (t)Z T (t)Z(t)P (t) (3.36)
pero P T (t)Z T (t)Y (t) = Y T (t)Z(t)P (t) pues se trata de un escalar y entonces:
83
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
La necesidad de que la matriz Z(t) sea de rango completo, nos lleva a profundizar
un tanto en los requerimientos que deben cumplir las señales de entrada u(k) a fin
de que el sistema sea identificable. Astrom y Eykhof (1975), introducen el concepto
de excitación persistente y lo consideran un requerimiento indispensable para que sea
posible obtener un estimador mı́nimo cuadrático único de los parámetros de la planta.
Este concepto puede expresarse en la forma siguiente:
Un sistema se encuentra sometido a una excitación persistente de orden n si existen
los lı́mites:
t
1X
ū = lı́m u(k) (3.40)
t→∞ t
k=1
t
1X
Rui = lı́m [u(k) − ū][u(k + i) − ū] (3.41)
t→∞ t
k=1
es definida positiva.
La condición de excitación persistente asegura no solamente que la matriz Z(t)
sea de rango completo sino también que la estimación de los parámetros del modelo
mediante el MMC sea estadı́sticamente consistente. Esta condición, por otra parte,
no impone serias restricciones sobre las señales de entrada u admisibles, cumpliéndose
en la gran mayorı́a de los casos prácticos, bajo condiciones normales de operación.
Debe señalarse, sin embargo, que en adición a la condición de excitación persistente
generalmente se requiere que las señales de entrada utilizadas sean perturbantes, es
decir, que su espectro de frecuencias sea suficientemente amplio, de manera que las
frecuencias fundamentales del sistema sean perturbadas. Esta última condición, sin
embargo, no resulta necesaria desde el punto de vista de la existencia de un estimador
mı́nimo cuadrático de los parámetros del proceso.
Regresando al problema fundamental que nos ocupa, resulta pertinente hacer un
análisis de algunas propiedades del estimador mı́nimo cuadrático de los parámetros,
cuya expresión se obtuvo en (3.39). La esperanza matemática del estimador, viene
dada por:
n o
E P̂ (t) = E (Z T (t)Z(t))−1 Z T (t)Y (t) (3.43)
84
3.3 Método de Mı́nimos Cuadrados para la Estimación de Parámetros del
Modelo ARX
La matriz de mediciones Z(t), contiene valores de u(k − i), y(k − i), para k =
1, 2, · · · , t; i = 1, 2, · · · , n tal como puede apreciarse en (3.22). Las variables de salida
y(k-i) a su vez contienen términos de ruido e(k − i) y por tanto, el producto de matrices
(Z T (t)Z(t))−1 Z(t), que es una pseudo-inversa de Z(t), contiene funciones de la secuen-
cia e(k − i), para i = 1, 2, · · · , n. Si se parte de la suposición de que el ruido e(k) es
blanco (ecuaciones (3.2) y (3.3) ), puede demostrarse que la esperanza matemática de
los términos aleatorios presentes en la pseudo-inversa de Z(t) es cero. Estos términos
tienen en general la forma f (e(k − i)e(k − j)), donde f es una función no lineal de la
secuencia e(k − i) y además i ̸= j, i ̸= 0. Por otra parte, admitiendo la no correlación
de los términos aleatorios presentes en (Z T (t)Z(t))−1 Z(t) e Y (t), consecuencia también
de la suposición de ruido blanco para e(k), la expresión (3.43) puede escribirse:
n o
E P̂ (t) = (Z T (t)Z(t))−1 Z(t)E {Y (t)} (3.44)
y entonces:
n o
E P̂ (t) = (Z T (t)Z(t))−1 Z T (t)Z(t)P (t) = P (t) (3.46)
R(t) = cov (Z T (t)Z(t))−1 Z T (t)Y (t) = (Z T (t)Z(t))−1 Z T (t)cov {Y (t)} Z(t)(Z T (t)Z(t))−1
(3.47)
Teniendo en cuenta a (3.33), se llega al resultado:
85
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
Pt h i
1 1
σ̂ 2 (t) = t
2
k=1 e(k) = t (Y (t) − Z(t)P̂ (t))T (Y (t) − Z(t)P̂ (t))
h i
1
= t (Y T (t)Y (t) − 2P̂ T (t)Z T (t)Y (t) + P̂ T (t)Z T (t)Z(t)P̂ (t)
h i
1
= t (Y T (t)Y (t) − 2P̂ T (t)Z T (t)(Z(t)P̂ (t) + E(t)) + P̂ T (t)Z T (t)Z(t)P̂ (t)
h i
1
= t (Y T (t)Y (t) − P̂ T (t)Z T (t)Z(t)P̂ (t) − 2t P̂ T (t)Z T (t)E(t)
1h T i
σ̂ 2 (t) = (Y (t)Y (t) − P̂ T (t)Z T (t)Z(t)P̂ (t) (3.49)
t
Hasta aquı́ hemos obtenido expresiones para los estimadores mı́nimos cuadráticos
de los parámetros del sistema, de la matriz de covarianza de la estimación y de la
varianza del ruido que pueden calcularse a partir de t mediciones de entrada-salida
del sistema que se identifica mediante el modelo ARX definido en la forma (3.31). Los
resultados obtenidos son válidos también para el caso multivariable, si P T (t), Y (t) y
Z(t) se definen como en (3.28) y (3.30). La única diferencia consiste en que en este
caso σ 2 (t) es una matriz ν × ν o sea, la estimación de la matriz de covarianza del ruido
y además, es necesario definir un arreglo tetra-dimensional para la covarianza de la
estimación de los parámetros, en la forma:
h i h i T
E P̂ (t) − P (t) P̂ (t) − P (t) 2
= σij (Z (t)Z(t))−1 (3.50)
ir js
86
3.4 Técnica del Olvido Exponencial.
87
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
E e2 (τ ) = φ−2(t−τ ) σ 2 (3.52)
donde v(τ ) es un proceso aleatorio con propiedades de ruido blanco y varianza cons-
tante:
E v 2 (τ ) = σ 2 (3.54)
Nótese que la diferencia entre (3.55) y (3.51) consiste en que los datos y(τ ) y z(τ )
son ahora pasados por el factor de olvido φ(t−τ ) . La introducción del factor de olvido
exponencial no altera la forma de la ecuación del estimador mı́nimo cuadrático de los
parámetros dados en (3.39), salvo que ahora las matrices Z(t) y Y (t) se definen en la
forma ligeramente modificada:
φ(t−1) z T (1) φ(t−1) y T (1)
φ(t−2) z T (2) φ(t−2) y T (2)
.. ..
Z(t) = . Y (t) = . (3.56)
φz T (t − 1) φy T (t − 1)
z T (t) y T (t)
88
3.4 Técnica del Olvido Exponencial.
φ n(φ2n = 0,1)
0.999 1155
0.998 580
0.997 385
0.996 290
0.995 230
0.994 195
0.993 165
0.992 145
0.991 134
0.990 125
89
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
En este epı́grafe serán derivadas las formas recursivas necesarias para calcular los
valores de los estimadores P̂ (t+1), σ̂ 2 (t+1) y la matriz C(t+1) que se define mediante:
−1
C(t + 1) = Z T (t + 1)Z(t + 1) (3.58)
a partir de P̂ (t), σ̂ 2 (t) y C(t) que se consideran conocidos, ası́ como de los datos obteni-
dos en el tiempo t + 1, es decir y(t + 1), u(t + 1). Según las definiciones de la matriz Z(t)
y el vector Y (t), dadas en (3.56), pueden obtenerse fácilmente las siguientes ecuaciones:
Resulta ahora necesario acudir a una bien conocida identidad matricial, cuya expli-
cación puede encontrarse, por ejemplo en (Householder, 1964):
−1 −1 T −1
A + BDB T = A−1 − A−1 B D−1 + B T A−1 B B A (3.64)
A = φ2 C −1 (t) (3.65)
B = z(t + 1) (3.66)
D=I (3.67)
90
3.5 Algoritmo Recursivo para la Estimación Mı́nimo Cuadrático.
se llega fácilmente a:
1
C(t + 1) = 2
C(t) − C(t)z(t + 1)(φ2 + z T (t + 1)C(t)z(t + 1))−1 z T (t + 1)C(t)
φ
(3.68)
Nótese que la expresión que hay que invertir en (3.68) es un escalar. Podemos
escribir entonces la siguiente expresión recursiva para actualizar la matriz C(t):
1 1 T
C(t + 1) = 2 C(t) − 2 C(t)z(t + 1)z (t + 1)C(t) (3.69)
φ α (t + 1)
donde:
α2 (t + 1) = φ2 + z T (t + 1)C(t)z(t + 1) (3.70)
91
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
h i
P̂ (t + 1) = P̂ (t) + C(t + 1)z(t + 1) y(t + 1) − z T (t + 1)P̂ (t) (3.78)
donde:
t
X
v(t + 1) = φ2i = 1 + φ2 v(t) (3.82)
i=0
α2 (t + 1) = φ2 + z T (t + 1)g(t + 1) (3.84)
92
3.5 Algoritmo Recursivo para la Estimación Mı́nimo Cuadrático.
1 1 T
C(t + 1) = 2 C(t) − 2 g(t + 1)g (t + 1) (3.89)
φ α (t + 1)
93
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
E1 (s)
+
U1 (s) + Y1 (s)
b
F11 (s)
+
F12 (s)
F21 (s)
+
U2 (s) + Y2 (s)
b
F22 (s)
+
E2 (s)
0,33 + 0,17z −1
F11 (z −1 ) = (3.91)
1 − 0,2z −1 − 0,3z −2
0,40 − 0,20z −1
F21 (z −1 ) = (3.92)
1 − 0,2z −1 − 0,3z −2
0,35 − 0,15z −1
F12 (z −1 ) = (3.93)
1 − 0,4z −1 − 0,2z −2
0,28 + 0,12z −1
F22 (z −1 ) = (3.94)
1 − 0,4z −1 − 0,2z −2
El modelo de regresión general expresado en la forma
94
3.5 Algoritmo Recursivo para la Estimación Mı́nimo Cuadrático.
B0 A1 B1 A2
0.331 0.389 0.205 0.001 0.169 -0.201 0.303 0.002
0.998 0.283 0.001 0.401 -0.153 -0.121 0.002 0.201
Puede apreciarse que las estimaciones obtenidas son muy exactas, a pesar del núme-
ro relativamente pequeño de etapas de tiempo simuladas, considerando que el número
de parámetros a estimar es en este caso de 16. Este resultado ilustra la fuerte relación
existente entre la exactitud de la identificación y el nivel de ruido en la salida.
95
3. IDENTIFICACIÓN DE SISTEMAS MEDIANTE EL MÉTODO DE
MÍNIMOS CUADRADOS
3.6. Conclusiones.
96
4
Objetivo: Después del tema cuatro los participantes explicarán las ventajas y des-
ventajas de la técnica de raı́z cuadrada en identificación de sistemas discretos con res-
pecto a métodos clásicos de identificación mediante el método de mı́nimos cuadrados
convencional, e incluyendo por lo memos el 80 % de los conceptos de algoritmos de raı́z
cuadrada.
4.1. Introducción
97
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
Ax = b (4.1)
AT Ax − AT b = AT e (4.3)
o sea, la solución mı́nimo-cuadrática del sistema (4.1) (Smuk, 1970) viene dada por la
expresión:
−1 T
x = AT A A b (4.5)
98
4.2 Solución Mı́nimo-Cuadrática de Sistemas sobre Determinados de
Ecuaciones Lineales
S = AT A = F T F (4.6)
F T F x = AT b (4.7)
donde F es la matriz n×n triangular superior o inferior según las conveniencias, definida
en (4.7). A partir de (4.8) se obtiene en forma simple la solución mı́nimo-cuadrática
deseada. Como se verá en lo que sigue, el método no requiere la factoración explı́cita
de la matriz AT A, ni tampoco la inversión de F T .
Consideremos, en efecto, la matriz de transformación:
j k
1
1
..
.
T =
c s
(4.9)
−s c
..
.
1
1
99
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
En la que se señalan solamente los elementos distintos de cero. Es fácil comprobar que
si s y c se escogen de manera que se cumpla:
c2 + s2 = 1 (4.10)
Entonces:
T T = T −1 (4.11)
Supongamos que queremos triangularizar a D haciendo cero los elementos d21 , d31
y d32 señalados en (4.14). Para anular el 1, escogemos la matriz de transformación:
c s 0
T ′ = −s c 0 (4.15)
0 0 1
100
4.2 Solución Mı́nimo-Cuadrática de Sistemas sobre Determinados de
Ecuaciones Lineales
−2s + c = 0 (4.17)
s2 + c2 = 1 (4.18)
Si se escoge ahora:
c′ 0 s′
T ′′ = 0 1 0 (4.22)
−s′ 0 c′
y procedemos en forma similar se logra anular el 2.
Finalmente, aplicando la transformación:
1 0 0
T ′′′ = 0 c′′ s′′ (4.23)
0 −s′′ c′′
∆ = T ′′′ T ′′ T ′ D = T D (4.24)
101
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
A partir de (4.25) y (4.26) y teniendo en cuenta además (4.18), se obtiene para d′11 , s
y c, las siguientes expresiones:
q
d21 d11
r = d211 + d221 , d′11 = r, s= , c= (4.27)
r r
Puede comprobarse fácilmente que en general, para cualquier etapa, se obtienen los
resultados:
q dji dii
r= d2ii + d2ji , d′ii = r, s= , c= (4.28)
r r
para i = 1, 2, · · · , n − 1, j = i + 1, i + 2, · · · , n. Un elemento de la fila k en (4.13) se
anula el cual corresponde a dji .
El algoritmo que ofrecemos a continuación da la posibilidad de triangularizar una
matriz cuadrada de dimensión n, utilizando el método de las transformaciones ortogo-
nales de rotación, aunque sin hacer explı́citamente el producto de matrices definido en
(4.24):
for i = 1, 2, · · · , n − 1 do
for j = i + 1, i + 2, · · · , n do
q
r = d2ii + d2ji
dji
s= r
dii
c= r
dii = r
dji = 0
for k = i + 1, · · · , n do
d = dik
dik = cd + sdjk
djk = cdjk − sd
end for
end for
end for
El método de las transformaciones ortogonales de rotación, puede emplearse tam-
bién ventajosamente para hallar la solución mı́nimo-cuadrática de un sistema sobrede-
terminado de ecuaciones lineales, lo cual es nuestro objetivo inmediato.
Consideremos, en efecto, el sistema sobredeterminado:
Ax − b = e (4.29)
102
4.2 Solución Mı́nimo-Cuadrática de Sistemas sobre Determinados de
Ecuaciones Lineales
El producto T D será una matriz triangular de dimensión n+1, por lo que el sistema
transformado (4.31) tendrá la forma:
n+1
x1
x2
..
. = Te
m
xn
−1
0, 0, 0, 0, · · · , 0
0, 0, 0, 0, · · · , 0
0, 0, 0, 0, · · · , 0
Figura 4.1: Matriz triangular de dimensión n + 1 - sistema transformado
103
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
for i = 1, 2, · · · , n do
q
r = d2ii + fi2
fi
s= r
dii
c= r
dii = r
for k = i + 1, · · · , n do
d = dik
dik = cd + sfk
fk = cfk − sd
end for
end for
La solución mı́nimo-cuadrática xT = (x1 , x2 , · · · , xn ) de (4.32) resulta evidentemen-
te de resolver las primeras n ecuaciones de manera que se anulen e′1 , e′2 , · · · , e′n . De esta
forma todo el error aparece concentrado en el elemento b′n+1 y
X 2 2
min e′ i = b′ n+1 (4.33)
Es fácil demostrar que el valor mı́nimo del criterio cuadrático (4.33) no se ha mo-
dificado por la transformación ortogonal. En efecto, se tiene que:
e′ = T e (4.34)
y también:
T
e′ e′ = eT T T T e (4.35)
T
eT e = e′ e′ (4.36)
D′ x̂ = b′ (4.37)
104
4.3 Algoritmo REDCO para la Identificación en Linea de Sistemas
Lineales Discretos
o sea, que la matriz triangular obtenida, no es otra cosa que la raı́z cuadrada de AT A.
La solución del sistema de ecuaciones:
F x̂ = b′ (4.39)
105
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
donde Z(t) es la matriz de las mediciones, de dimensión t × (2n + 1), siendo t el número
de observaciones del par entrada-salida y n el orden del modelo. La traspuesta de esta
matriz se da en (4.21). Y (t) es el vector de las salidas de dimensión t. P es el vector de
parámetros y E(t) es el vector de ruido.
El problema de identificación de los parámetros del modelo (4.40), contenidos en el
vector P , puede interpretarse como el de la solución de un sistema sobredeterminado
de ecuaciones, ya que en general, el número de observaciones t puede hacerse arbitra-
riamente grande y siempre mayor que el número de parámetros 2n + 1. La solución
mı́nimo cuadrática de este sistema, es como se sabe:
Z T (t)Z(t) P = Z T (t)Y (t) (4.41)
F (t)
q(t)
2n + 2
106
4.3 Algoritmo REDCO para la Identificación en Linea de Sistemas
Lineales Discretos
donde
z T (t + 1) = [u(t + 1), y(t), u(t), · · · , y(t − n + 1), u(t − n + 1)] (4.44)
2n + 2
F (t)
q(t)
2n + 3
minE T (t)E(t)
z T (t + 1) y(t + 1)
Este sistema puede ser de nuevo triangularizado anulando la última fila, correspon-
diente a las mediciones en t + 1, utilizando para ello el método de las transformaciones
ortogonales de rotación. Se obtiene entonces la matriz:
La nueva estimación de los parámetros, puede calcularse fácilmente, resolviendo el
sistema triangular de ecuaciones:
Nótese que si los parámetros no varı́an rápidamente con el tiempo, caso muy fre-
cuente en la práctica, la actualización correspondiente al tiempo t + 1, puede consistir
solamente en la anulación de la fila [z T (t + 1), y(t + 1)], calculándose explı́citamente
los parámetros después de un número arbitrario de etapas de tiempo, con lo que se
107
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
2n + 2
F (t + 1)
q(t + 1)
minE T (t + 1)E(t + 1)
F (0) = 0 (4.46)
q(0) = 0 (4.47)
108
4.3 Algoritmo REDCO para la Identificación en Linea de Sistemas
Lineales Discretos
El método descrito conocido como algoritmo REDCO, puede ser utilizado también
en el caso de sistemas multivariables con µ entradas y ν salidas. En este caso el modelo:
donde F (t) es la raı́z cuadrada de la matriz Z T (t)Z(t) con dimensión ρ y q(t) es una
matriz ρ × ν.
La matriz compuesta que se obtiene como resultado de la aplicación del método
REDCO tendrá la estructura:
ρ+ν
F (t)
q(t) ρ
S(t) ν
El pequeño triángulo S(t), contiene información acerca del mı́nimo del criterio
cuadrático, que en el caso de sistemas multivariables puede definirse:
min traza E T (t)E(t) = traza S T (t)S(t) (4.53)
109
4. ALGORITMOS DE RAÍZ CUADRADA PARA LA
IDENTIFICACIÓN DE MODELOS DISCRETOS
110
5
Objetivo:
Al término del tema tres los participantes aplicarán las diferentes técnicas de control
clásico según los lineamientos del tema, incluyendo al menos el 80
5.1. Introducción
111
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
112
5.2 Sistemas de Control Automático
tránsito mediante señales operadas en una base de tiempo es otro ejemplo de control
en lazo abierto.
Los sistemas de control en lazo abierto constituyen el tipo más sencillo y económico
de los sistemas de control. En la figura 5.1 se muestra un sistema de este tipo repre-
sentado en forma de diagrama de bloques.
Entrada Salida
CON T ROLADOR P LAN T A
r(t) y(t)
113
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
114
5.2 Sistemas de Control Automático
Detector de error
Salida
Referencia
Controlador Elemento F inal de Control P lanta b
+
Error
- Control Actuador
Bucle de realimentacion
Sensor
Está claro que de este modo, se obtiene un sistema de control con un comporta-
miento totalmente automático, ya que no es precisa la intervención de un operador
para adaptar la salida a la entrada, de ahı́ el nombre con que también se conocen los
sistemas re-alimentados o de lazo cerrado: sistemas de control automático.
115
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
Una ventaja del sistema de control en lazo cerrado es que el uso de la re-alimentación
vuelve la respuesta del sistema relativamente insensible a las perturbaciones externas
y a las variaciones internas en los parámetros del sistema. Por tanto, es posible usar
componentes relativamente precisos y baratos para obtener el valor adecuado de una
planta determinada, en tanto que hacer eso es imposible en el caso de un sistema en
lazo abierto.
Desde el punto de vista de la estabilidad, el sistema de control en lazo abierto es más
fácil de desarrollar, porque la estabilidad del sistema no es un problema importante.
Por otra parte, la estabilidad es una función principal en el sistema de control en lazo
cerrado, lo cual puede conducir a corregir en exceso errores que producen oscilaciones
de amplitud constante o cambiante.
Debe señalarse que, para los sistemas en los que se conocen con anticipación las en-
tradas y en los cuales no hay perturbaciones, es aconsejable emplear un control en lazo
abierto. Los sistemas de control en lazo cerrado sólo tienen ventajas cuando se presen-
tan perturbaciones impredecibles y/o variaciones impredecibles en los componentes del
sistema. Observe que la valoración de la energı́a de salida determina en forma parcial
el costo, el peso y el tamaño de un sistema de control. La cantidad de componentes
usados en un sistema de control en lazo cerrado es mayor que la que se emplea para un
sistema de control equivalente en lazo abierto. Por tanto, el sistema de control en lazo
cerrado suele tener costos y potencias más grandes. Para disminuir la energı́a requerida
de un sistema, se emplea un control en lazo abierto cuando puede aplicarse. Por lo
general, una combinación adecuada de controles en lazo abierto y en lazo cerrado es
menos costosa y ofrecerá un desempeño satisfactorio del sistema general.
En los inicios de la era industrial, el control de los procesos se llevó a cabo mediante
tanteos basados en la intuición y en la experiencia acumulada. Un caso tı́pico es el
control de acabado de un producto en el horno donde el operario era realmente el
“instrumento de control”que juzgaba la marcha del proceso por el color de la llama,
por el tipo de humo, el tiempo transcurrido y el aspecto del producto y decidı́a ası́ el
116
5.3 Acciones de Control.
momento de retirar la pieza; en esta decisión influı́a muchas veces la suerte, de tal modo
que no siempre la pieza se retiraba en las mejores condiciones de fabricación. Más tarde
el mercado empezó a exigir mayor calidad en las piezas fabricadas lo que condujo al
desarrollo de teorı́as para explicar el funcionamiento del proceso, de las que derivaron
estudios analı́ticos que a su vez permitieron realizar el control de la mayor parte de las
variables de interés en los procesos.
En la figura 5.3 se muestra un esquema general del bucle de control re-alimentado, la
función de transferencia del proceso normalmente incluye, además de su propia dinámi-
ca, los bloques correspondientes a actuadores y/o bloques de amplificación necesarios
para accionarlos. La mayorı́a de los controladores automáticos industriales son de tipo
electrónico, hidráulico, neumático o alguna combinación de éstos.
U (s) Planta
R(s) E(s) Y (s)
Gc (s) Amplif. Actuador Proceso b
+
-
G(s)
Sensor
117
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
cierto tiempo hasta que el operario empiece a detectar que el agua está más frı́a. Cuan-
do éste note la disminución de temperatura deberá compararla con la temperatura que
desea y calcular mentalmente cuantas vueltas debe darle a la válvula de vapor y en que
sentido, y a continuación realizar la corrección, necesitando para toda esta operación,
un cierto tiempo. También pasa un cierto tiempo hasta que los efectos de la corrección
de la posición de la válvula se notan en la temperatura de salida y son sentidos por
el operador. Sólo entonces es capaz de saber si su primera corrección ha sido escasa o
excesiva. En este punto efectuará una segunda corrección, que al cabo de un tiempo
dará lugar a otro cambio de temperatura. El operador observará los resultados de esta
segunda corrección y realizará otra tercera, y ası́ sucesivamente, constituyendo el lazo
cerrado de control manual.
Es evidente que si el tanque contiene una gran cantidad de agua de tal modo que
esa masa considerable de lı́quido estabiliza y resiste los cambios que puedan causarse
a la temperatura por variaciones en el caudal de entrada o de salida; en la presión
del vapor, o en la temperatura ambiente, el operario deberı́a estar muy atento y le
serı́a casi imposible mantener la temperatura de salida del agua en un valor constante.
Se llega ası́ a la necesidad de introducir un dispositivo dentro del lazo con objeto de
automatizar el control. El controlador es el elemento encargado de procesar la
señal de error y generar otra capaz de disminuir su valor para conseguir la
máxima precisión.
Pero el controlador aún tiene más tareas. Cuando se plantea un problema de control,
por lo general se especifican un conjunto de restricciones que debe cumplir el sistema
para considerar que está operando adecuadamente. Estas restricciones están dadas en
términos de estabilidad, velocidad de respuesta y precisión, las tres premisas
básicas del control. Se dice que un sistema está controlado cuando todas ellas se
cumplen satisfactoriamente.
Desafortunadamente suele suceder que cuando se mejora alguna de ellas, una de las
dos restantes o ambas se degradan, por lo que se hace necesario establecer un compromi-
so que satisfaga, en la medida de lo posible, las especificaciones funcionales del sistema;
ası́ el controlador debe también garantizar una determinada estabilidad re-
lativa y una velocidad de respuesta dentro de los márgenes especificados.
El procedimiento mediante el cual el controlador genera la señal de control se de-
nomina acción de control.
118
5.3 Acciones de Control.
En general, el controlador debe mejorar las prestaciones del sistema, lo que supone:
Dado que la planta Gp (s) ası́ como los elementos de medición de un sistema H(s)
son fijos, es decir, no los podemos alterar para ejecutar nuestra acción de control, sólo
nos queda modificar la manera en que manipulamos las señales dentro del lazo. El
estudio de los controladores, se aplica suponiendo una configuración de sistema de un
único lazo (o un sistema que se pueda reducir a esta configuración).
La acción de control se refiere a como se emplea la señal de actuación u(t) en los
elementos de control (controlador o regulador) para lograr la corrección de la salida.
En el control manual descrito en el ejemplo, el operador puede hacer las correcciones
en la válvula de vapor de diversas formas:
Puede abrir la válvula un número constante de vueltas, por cada unidad de des-
viación
119
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
Este modo de control se caracteriza por tener dos posiciones fijas. Si u(t) es la salida
de control y e(t) la señal de error, su funcionamiento queda descrito tanto por la figura
5.4. como por su expresión matemática:
u(t)
u1
−u1
120
5.4 Control Todo-Nada (On-Off )
u(t)
u1
+u1 para e(t) > bm bm
u(t) = 0 para |e(t)| ≤ bm e(t)
−u1 para e(t) < bm
−u1
121
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
u(t)
u1
u para e(t) > bm
u(t) = 1 e(t)
0 para e(t) < −bm
En este tipo de control la relación entre la salida del controlador u(t) y la señal de
error e(t) es:
u(t) = Kp e(t) (5.1)
U (s)
= Kp (5.2)
E(s)
donde Kp se denomina ganancia proporcional.
Cualquiera que sea el mecanismo real y la forma de aplicar la potencia a los ac-
cionamientos, el controlador proporcional es en esencia un amplificador de ganancia
ajustable, ya que existe una relación lineal continua entre el valor de la variable con-
trolada y la posición del elemento final de control.
Veamos cómo reacciona un sistema al ser excitado por un escalón de amplitud A.
Analizaremos la respuesta en régimen estacionario a entrada escalón, tanto para el
sistema en bucle abierto, como para el sistema en bucle cerrado.
En bucle abierto, a una señal de entrada en escalón de amplitud A, el proceso
responde con una salida:
122
5.5 Control Proporcional (P).
A
Y (s) = U (s)G(s) = G(s) (5.3)
s
G(s)
yss = lims→0 sY (s) = lims→0 sA = A · G(0) (5.4)
s
A Kp G(s) Kp G(0)
yss = lims→0 sY (s) = lims→0 s = .A (5.5)
s 1 + Kp G(s) 1 + Kp G(0)
Se observa que para valores altos de Kp la salida del sistema prácticamente seguirá
a la señal de referencia, en efecto:
123
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
A Ki G(s)
yss = lims→0 sY (s) = lims→0 s =A (5.7)
s s + Ki G(s)
y entonces ess → 0.
Sin embargo, la acción de control integral puede empeorar de un modo substancial la
estabilidad relativa del sistema, aumentando el sobre impulso de la respuesta transitoria,
y pudiendo llegar a hacer que se vuelva inestable, debido al desplazamiento de los polos
de lazo cerrado hacia la derecha. Por ello puede resultar conveniente que la acción
integral se acompañe de otras acciones de control.
124
5.6 Control Proporcional-Integral (PI).
Dar un valor elevado a Ti (por ejemplo 100 minutos) de forma que a efectos
prácticos, el controlador PI se comporte como P .
Disminuir poco a poco Ti hasta el mı́nimo posible, de forma que se anule el ess pero
manteniendo la forma de respuesta deseada. Es decir, debe buscarse el mı́nimo
valor de Ti que no afecte de forma significativa al coeficiente de amortiguamiento.
125
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
126
5.8 Control Proporcional Integral Derivativo (PID).
La acción de control PID genera una señal resultado de la combinación de las tres
acciones ya estudiadas.
Z t
de(t)
u(t) = Kp e(t) + Kd + Ki e(t)dt
d(t) 0
Z
de(t) 1 t
= Kp e(t) + Td + e(t)dt (5.13)
dt Ti 0
127
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
128
5.8 Control Proporcional Integral Derivativo (PID).
Ke−sθ
Geq (s) = (5.15)
Ts + 1
siendo θ el retardo y T la constante de tiempo del sistema de primer orden equivalente.
b) Calcular los parámetros del sistema equivalente de primer orden más retardo
aplicando las expresiones empı́ricas: K = valor final para entrada escalón unitario,
T = 1,5(t6 − t2 ), (5.16)
θ = t6 − T. (5.17)
Aplicar las fórmulas del primer método de Ziegler-Nichols para calcular los coefi-
cientes del controlador PID:
T
Kp = 1,2 (5.18)
Kθ
Ti = 2θ (5.19)
Td = 0,5θ (5.20)
129
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
estos son los valores de los parámetros correspondientes al sistema de primer orden
equivalente.
Finalmente, los coeficientes del controlador PID son:
12,3
Kp = 1,2 = 4,2, (5.24)
3,5
Ti = 2 × 3,5 = 7, (5.25)
El análisis del sistema es un método para predecir la respuesta del sistema con base
al modelo matemático. El principal objetivo del análisis del sistema es la respuesta de
las siguientes dos cuestiones:
Es el sistema estable y si lo es, ¿cómo lo es?
Que tan rápida es la respuesta del sistema.
La estabilidad del sistema indica la respuesta natural. Como se mueve el motor de
una posición a otra, esto puede ocasionar un sobre impulso en la tarjeta ocasionando
tiempos crı́ticos antes de estabilizarse. La importancia del sobre impulso está directa-
mente relacionada con el grado de estabilidad, el cual es llamado amortiguamiento. Del
más bajo amortiguamiento, al más largo sobre impulso. Un sistema inestable causa que
130
5.9 Análisis del Sistema de Control
1. Iniciar con el modelo matemático del sistema y derivar L(s), la función de trans-
ferencia en lazo abierto. En un sistema simple con un simple lazo de retroalimen-
tación, L(s) es el producto de todos los elementos en el lazo. En sistemas más
complejos, derive la función de transferencia entre la entrada y la salida. Reorde-
ne los términos en el denominador para tenerla de la forma 1 + L(s). La función
resultante, L(s), es la función de transferencia en lazo abierto correspondiente.
2. Reemplace los parámetros s en L(s) por los parámetros jω. La función resultante
L(jω) es llamada respuesta en frecuencia en lazo abierto
|L(jω)| en ω = ωc (5.27)
4. Encontrar la fase de L(jωc ). Dado que L(jω) es una función compleja, se encon-
trara su argumento, y se llamará el argumento como Φ.
131
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
Φ = arg[L(jωc )] (5.29)
Θ = 180 + Φ (5.30)
+
-
Encoder
Kf
El motor que se controlará tiene los parámetros que se muestran en la tabla 5.1.
Se asume por adelantado que los parámetros del filtro son un controlador PD:
P = 100 (5.31)
132
5.9 Análisis del Sistema de Control
Tabla 5.1: Parámetros del motor - Servomotor sin escobillas de imán permanente
D = 0,4 (5.32)
133
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
considerando los parámetros dados en la tabla 5.1 y la ecuación 5.39, L(s) es igual a
541435,6(s + 250)
L(s) = (5.41)
s(s + 187,6)(s + 2246,3)
|L(jωc )| = 1 (5.42)
El resultado del análisis indica que el sistema es estable. El margen de fase de 75.2
grados representa un buen grado de amortiguamiento y la frecuencia de cruce de 268.1
rad/seg representa el tiempo de respuesta que es de aproximadamente de 3.7 ms.
⇒ Tarea 5.1.
134
5.10 Diseño del Control PID
Se iniciará asumiendo que todos los elementos del sistema son fijos y que su modelo
matemático es conocido. Como el amplificador, el sensor y la ganancia del DAC son
considerados fijos en este punto, se definirá una sola ganancia combinada por K. La
función de transferencia del motor es considerando la posición como salida H(s). Los
elementos del controlador que deben ser considerados ahora son el filtro F (s), y el
retenedor de orden cero representado por G(s). El diagrama a bloque resultante se
presente en la siguiente figura 5.8
Filtro Ganacia Motor
F (s)G(s) K H(s) b
+
-
El objetivo del diseño es seleccionar los parámetros del filtro de tal manera que
produzca una respuesta deseada. La respuesta del servosistema está caracterizada por
la velocidad de respuesta y el nivel de estabilidad. La velocidad de respuesta está de-
terminada por la frecuencia de cruce ωc . Ésta es la frecuencia en la cual la magnitud
de la función de transferencia en lazo abierto es uno. El nivel de estabilidad está de-
terminado por el margen de fase Θ. Donde si Θ está entre 30 y 45 grados, la respuesta
del sistema es amortiguada y estable. El margen de fase más pequeño indica respuesta
de bajo amortiguamiento, y por lo tanto, grandes sobre impulsos. Pero si es mayor de
45 grados, la respuesta puede ser sobre amortiguada y lenta.
Se asumirá que el filtro es de tipo PID con coeficientes P , I y D.
I
F (s) = P + sD + (5.47)
s
Para un diseño preciso, se considerará la naturaleza digital del controlador. Cuando
el controlador opera con un periodo de muestreo T , el efecto del muestreo en el sistema
dinámico es aproximado a la de un retardo de medio periodo de muestreo. Esto está
expresado por G(s)
−sT
G(s) = e 2 (5.48)
135
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
|L(jwc )| = 1 (5.50)
y que
Θ = 180 + Arg [L(jwc )] (5.51)
136
5.10 Diseño del Control PID
ecuación (5.50). Como quiera que sea el desvió de fase asociada con el retardo se
encuentra como
−jwc T
G(jwc ) = e 2 (5.53)
h −jwc T i −w T −180wc T o
c
Arg e 2 = rad = (5.54)
2 2π
La simplificación nos permite iniciar el proceso de diseño con P y D, y se escribe
la ecuación (5.50) y (5.51) como
Θ = 180 + Arg[L(jwc )]
−1 wc D 180wc T
= 180 + Arg[H(jwc )] + tan − −5 (5.57)
P 2π
Porque H(s) es conocida, la magnitud del argumento de H(jwc ) puede ser encon-
cos(B)
trada, llevando estas dos simples ecuaciones algebraicas (5.58) y (5.59) en P = KA
sen(B)
yD= ωc KA puede ser resueltas relativamente fácil.
Para sumar la constante I, se asumirá que el desvı́o causado por la constante
integral es de -5 grados, por lo que el arg(−5o ) = I/jωc P , despejando I tenemos
I = P ωc tan(5o ).
Considerando el sistema de control para el motor, sensor y amplificador utilizado
en la sección anterior, se toman los parámetros dados en la tabla 5.1.
rad
El objetivo del diseño es lograr una frecuencia de cruce de wc = 200 seg con un
margen de fase de Θ = 45o y un Ts = 0,001 seg.
La función de transferencia del motor es
4,1322 1741654,5
H(s) = = (5.60)
s(0,0053s + 1)(0,000445s + 1) s(s + 187,6)(s + 2246,3)
137
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
La ganancia combinada es
K = Kd Kv Kf = 0,7764 (5.61)
rad
La magnitud y el argumento de H(s) a la frecuencia de cruce de 200 seg son
encontrados de sustituir s = j200 en (5.60)
1741654,5
H(j200) = (5.62)
j200(j200 + 187,6)(j200 + 2246,3)
resultando en
A = |H(j200)| = 0,0141 (5.63)
El correspondiente ángulo B
B = 17,65o (5.65)
⇒ Realizar Práctica 5
5.11. Conclusiones
138
5.11 Conclusiones
139
5. TÉCNICAS CLÁSICAS DE CONTROL DE PROCESOS
140
6
Objetivo:
Al finalizar el estudio del cuarto tema los participantes manejarán técnicas de con-
trol discreto basándose por lo menos en tres técnicas revisadas en el curso, incluyendo
al menos el 80
6.1. Introducción
En este capı́tulo vamos a hacer una incursión por algoritmos y métodos de control
digital bien establecidos en la práctica industrial tales como los reguladores PID, las
conexiones anticipativas (feed-forward), predictores de Smith, PID autoajustables, etc.
La razón de incluir estos temas en un libro básicamente dedicado a la identifica-
ción y control adaptable; es nuestro convencimiento de que cualquier nuevo método o
enfoque de control, si pretende ganar un espacio en la práctica y no ser simplemente
un tema académico, debe compararse ventajosamente en algún sentido con las solu-
ciones clásicas. Esta comparación debe ser, además, lo más justa posible; es decir, se
trata de tomar como patrón la mejor variante clásica posible, en la que se incorporan
reguladores PID profesionales y un conjunto de adiciones o extensiones que en su con-
junto permiten dar respuestas satisfactorias a una muy amplia variedad de problemas
de control industrial.
También en este capı́tulo se va a tratar un tema más cercano al espı́ritu de los
restantes y que presenta un interés especial: algoritmos PID auto ajustables. Se trata
de una clase de control adaptable que ha recibido mucha atención en la literatura de
141
6. TÉCNICAS DE CONTROL DISCRETO
142
6.2 Algoritmos PID Discretos.
donde:
∆u(k) = u(k) − u(k − 1) (6.3)
ó
u(k) = u(k − 1) + ∆u(k) (6.4)
Una serie de estudios comparativos han sido publicados sobre las dos formas prin-
cipales del algoritmo PID discreto. Puede consultarse, por ejemplo, el libro de Astrom
y Witenmark (1984).
A continuación haremos algunos comentarios sobre el tema:
143
6. TÉCNICAS DE CONTROL DISCRETO
144
6.3 Versión Profesional del Algoritmo PID.
Ru Td ∆y(k)
u(k) = en (k) − + int(k) (6.5)
Bp Tc Ry
donde Ru = umax − umin es el intervalo del actuador a menudo 0 − 100 %. Ry =
ymax − ymin es el intervalo del transmisor con que se mide la variable en unidades
fı́sicas, por ejemplo: 0 − 100 C ◦ , 4,2 − 20 Kg/cm2 , etc.
yr (k) − yf (k)
en (k) = (6.6)
Ry
es el error normalizado, donde yr (k) es el valor de referencia, eventualmente filtrado
mediante un filtro exponencial de primer orden. yf (k) es el valor de la salida filtrada.
Bp es la banda proporcional. Se relaciona con la ganancia mediante Kp = Ru /Bp . En
muchas ocasiones Kp = 100/Bp .
Tc Ru
int(k) = int(k − 1) + en (k) (6.7)
Ti Bp
int(0) = u(0) (6.8)
y
e(k − 1) = yr (k − 1) − y(k − 1) (6.10)
145
6. TÉCNICAS DE CONTROL DISCRETO
sea en este caso el uso del negativo del incremento de la salida es perfectamente equi-
valente al del incremento del error para conformar el modo derivativo del PID. No
obstante, si se produce un cambio de la referencia, tendremos que:
A diferencia del cambio en la salida, que por naturaleza es lento dada la dinámica
del proceso, el cambio de la referencia es inmediato y brusco, lo cual significa que
el incremento del error y por ende la componente derivativa del regulador también
cambiará bruscamente. Este es precisamente el efecto desestabilizador que se pretende
evitar cuando se sustituye al incremento del error por el incremento de la salida en el
modo derivativo. La variable u(k) generalmente se encuentra en un rango normalizado
en la computadora (0 − 100) ó (0 − 1). Después, en el convertidor D/A se transforma en
una señal normalizada, p. ej. de 4 a 20 ma. Si el convertidor de salida es de 12 bits y el
rango interno de u(k) es 0 − 100, el valor de entrada al convertidor, se calcula mediante
la expresión:
u(k)
usal = × 4095 (6.13)
100
Es necesario tener en cuenta dos detalles que a menudo se confunden y que sin
embargo tienen un tratamiento distinto. El primero consiste en el signo de la ganancia
del proceso. Si este signo es positivo, es decir, a un aumento de la variable de control
corresponde un aumento de la salida, entonces el signo de Bp (o de Kp ) en el PID
debe escogerse positivo. En caso contrario, es decir, cuando la ganancia del proceso es
negativa, puede entonces escogerse Bp negativo y se obtiene el comportamiento estable
deseado.
El segundo aspecto que hay que tomar en consideración es la naturaleza de la válvula
de control. La expresión (6.13) es válida en el caso de que la válvula sea normalmente
cerrada (Acción Directa), que es la más frecuente. Existen casos en la práctica, no
obstante, de válvulas normalmente abiertas (Acción inversa), que cierran a medida que
aumenta la señal de control. En este caso, el valor que debe entregarse al convertidor
es:
u(k)
usal = 4095 × 1 − (6.14)
100
Cualquier regulador digital y en particular los PID, deben tener al menos 3 modos
de operación: 1- Respaldo (Back-up), 2- Manual y 3- Automático. En el modo Respaldo,
146
6.4 Filtraje de las Variables
el regulador sólo recibe las mediciones y el control se ejerce mediante otro dispositivo
externo. La salida del regulador no está fı́sicamente conectada a la entrada del proceso.
En el modo Manual, desde el regulador se puede actuar sobre la válvula cerrándola
o abriéndola. Es frecuente que para este fin se utilicen 2 teclas, una de incremento y
otra de decremento.
En el modo Automático, el regulador digital asume la función de control, calculando
el valor de la variable de control mediante la expresión (6.5) u otra similar. En este
modo el operador puede generalmente modificar el valor de referencia (set-point) de la
variable controlada.
Deben tomarse algunas previsiones en la inicialización y cambio del modo de un
regulador digital PID si se requiere un comportamiento suave y libre de saltos bruscos.
A continuación detallamos las más importantes:
Cuando se pasa del modo Respaldo al modo Manual, debe hacerse el valor
int(0) = u(0). Esto puede lograrse automáticamente si se mide la variable de
control u(t) (posición de una válvula, p. ej.) o manualmente, leyendo el valor de
u(t) que tiene el equipo de respaldo e insertándolo en la computadora.
147
6. TÉCNICAS DE CONTROL DISCRETO
Tm 12 3 12 3 12 3 12 3 12 3
El tipo de filtro para ruidos de frecuencias baja y media más utilizados en la práctica
es el llamado filtro exponencial de primer orden cuya expresión es:
donde α es un valor real en el intervalo 0-1. yf (k) y yf (k − 1) son los valores filtrados
en los perı́odos de muestreo k y k − 1 respectivamente. ym (k) es el valor medido en el
perı́odo k.
La ecuación (6.15) puede también ser expresada mediante:
n
X
yf (k) = α (1 − α)i ym (k − i) + (1 − α)n+1 yf (k − n − 1) (6.16)
i=0
148
6.4 Filtraje de las Variables
Yf (z) α
= (6.21)
Ym (z) 1 − (1 − α)z −1
cuando z → 1 se obtiene:
Yf (z)
=1 (6.22)
Ym (z)
O sea, que el valor filtrado converge al valor de la medición si este se mantiene
constante durante un número suficientemente alto de perı́odos de muestreo.
Es importante tener en cuenta que el filtraje de las mediciones introduce un retardo
adicional en el lazo de control, cuyo efecto crece a medida que el parámetro α es menor.
Por otra parte se conoce que el efecto pernicioso del ruido es más notable sobre el modo
derivativo en un regulador PID. Debe existir entonces un compromiso entre el filtraje
de la señal y el retardo introducido que tiende a hacer inestable al lazo de control. Una
buena solución parece ser trabajar con 2 valores diferentes del coeficiente de filtraje,
uno mayor para los modos P e I y otro más pequeño para el modo D. Es decir, se
calculan 2 valores distintos de la variable de salida filtrada:
149
6. TÉCNICAS DE CONTROL DISCRETO
La elección de los valores concretos del coeficiente α, debe hacerse en cada caso
particular según la siguiente heurı́stica: Para perı́odos de control largos (100 o más
perı́odos de muestreo) y señales ruidosas, se seleccionan valores pequeños de α (0.1
o menores). En situaciones intermedias (T c ≤ 5Tm ) y señales filtradas por hardware,
usualmente pueden escogerse valores de α > 0,4. En última instancia, en cada lazo de
control debe ajustarse experimentalmente el valor más adecuado de α.
ke−θs
F (s) = (6.27)
τs + 1
150
6.6 Análisis de un Sistema en Tiempo Discreto.
orden cero (Sampler and Hold). Entonces es posible utilizar los parámetros que se dan
en la tabla 6.1, donde Kp es la ganancia del proceso, K la ganancia del regulador y Ti ,
Td son respectivamente los tiempos integral y derivativo.
151
6. TÉCNICAS DE CONTROL DISCRETO
e(h − 1) + e(h)
f (h) = , f (0) = 0 (6.31)
2
Entonces se escribe
( k ) ( k )
X e(h − 1) + e(h) X 1
Z =Z f (h) = F (z) (6.32)
2 1 − z −1
h=1 h=1
Pk
⇒ Tarea 6.1 Calcular la transformada Z de h=1 f (h).
Además se tiene que la transformada Z de f (h) es
1 + z −1
F (z) = Z {f (h)} = E(z) (6.33)
2
Por lo tanto ( k )
X e(h − 1) + e(h) 1 + z −1
Z = E(z) (6.34)
2 2(1 − z −1 )
h=1
152
6.6 Análisis de un Sistema en Tiempo Discreto.
o
Ki
U (z) = [Kp + + Kd (1 − z −1 )]E(z) (6.37)
1 − z −1
donde
KT Ki
Kp = K − =K− (6.38)
2Ti 2
es la ganancia proporcional,
KT
Ki = (6.39)
Ti
es la ganancia integral y
KTd
Kd = (6.40)
T
es la Ganancia derivativa.
A esta función de transferencia del PID discreto se le conoce comúnmente como
forma posicional del esquema PID. Considerando los resultados de la sección anterior
se tiene que
153
6. TÉCNICAS DE CONTROL DISCRETO
154
6.6 Análisis de un Sistema en Tiempo Discreto.
similar a la del ruido aleatorio. Y en efecto, el proceso de cuantificación actúa como una
fuente de ruido aleatorio. Según (Ogata, K.1996), si el nivel cuantificación Q es pequeño
comparado con la amplitud promedio de la señal de entrada, entonces la varianza del
ruido de cuantificación es un doceavo del cuadrado de cuantificación. En nuestro caso
se tiene que esta varianza es igual a un doceavo de cuenta, es decir 0.00013 radianes
de varianza del error de la posición, el cual no es significativo para el intervalo que se
maneja.
como sigue:
1. Para que el sistema sea estable, los polos en lazo cerrado o las raı́ces de la ecuación
caracterı́stica deben presentase en el plano z dentro del cı́rculo unitario. Cualquier
polo en lazo cerrado exterior al cı́rculo unitario hace inestable al sistema.
3. Los ceros en lazo cerrado no afectan la estabilidad absoluta y por lo tanto pueden
quedar en cualquier punto del plano z.
155
6. TÉCNICAS DE CONTROL DISCRETO
1 − e−T s 1741654,5
Gp (s) = (6.51)
s s(s + 187,6)(s + 2246,3)
1−e−T s
F (z) s H(s) b
+
-
156
6.6 Análisis de un Sistema en Tiempo Discreto.
z1 = 0,9812
z2 = 0,8612 + j0,2120
z3 = 0,8612 − j0,2120
z4 = 0,2324
z5 = −0,03991 (6.59)
En vista que las magnitudes de los polos del sistema en lazo cerrado están dentro
del cı́rculo unitario del plano z el sistema es estable. Se debe de tener especial cuidado
con los errores de cuantificación, errores por redondeo, dado que estos pueden mover
la posición de las raı́ces y ubicarlas fuera del cı́rculo unitario e inestabilizar el sistema.
En este caso no se tiene este problema debido a la resolución que se maneja.
Las caracterı́sticas de la respuesta transitoria del sistema de control en tiempo
discreto dependen del periodo de muestreo T . Un periodo de muestreo grande tiene
efectos dañinos sobre la estabilidad relativa del sistema. Como se ve en el análisis
precedente, para los valores dados de las constantes del PID digital, aumentar el periodo
de muestreo T hará que el sistema de control en tiempo discreto sea menos estable y que
eventualmente se convierta en inestable. De manera alternativa, al reducir el periodo de
muestreo T , se permite que los valores crı́ticos de las ganancias respecto a la estabilidad
sean mayores. De hecho, reducir el periodo de muestreo más y más tiende a hacer que
el sistema se comporte muy parecido a un sistema en tiempo continuo.
En el sistema mostrado en la figura 6.2 se puede observar que las raı́ces del sistema
dependen del periodo T . Si se aumenta el periodo T a 0.01 y se aplica el procedimiento
anterior se obtienen las raı́ces.
157
6. TÉCNICAS DE CONTROL DISCRETO
z1 = −0,5593 + j1,0950
z2 = −0,5593 − j1,0950
z3 = 0,8270
z4 = 0,1529
z5 = −0,00186 (6.60)
En donde, se puede ver que z1 y z2 están fuera del cı́rculo unitario del plano z y
por lo tanto el sistema es inestable en lazo cerrado.
Para comprender mejor el efecto del control en cascada, es conveniente recordar una
propiedad general del feed-back negativo. Supongamos un bloque de primer orden, con
constantes de tiempo T1 y ganancia K:
U (s) K Y (s)
T1 s+1
Yr (s) K Yb (s)
K1 T1 s+1
+
-
1 KK KK1
y(s) KK1
= T1 s+1
KK
= = 1+KK1
T1
(6.61)
yr (s) 1 + T1 s+1
1 T1 s + (1 + KK1 ) 1+KK1 s + 1
Los efectos sobre la función transferencial original son dos:
K1
Una reducción de la ganancia en el factor 1+KK1
158
6.7 Control en Cascada
Una reducción aún mayor de la constante de tiempo, que queda dividida por
(1 + KK1 ).
Pv T0
D1 D2
+ +
Td u Fv Ta
R G1 G2 b
+
-
Pv T0
D1 D2
+ +
Td Fd u + Fv + Ta
R2 R1 G1 b G2 b
+ +
- -
159
6. TÉCNICAS DE CONTROL DISCRETO
El perı́odo de control del regulador esclavo debe ser igual o menor que el del
regulador maestro. Esto es consecuencia de la suposición inicial de que el sub-
proceso G1 tiene una dinámica más rápida que G2 .
El ajuste de los parámetros se hace siempre comenzando por el lazo esclavo, con
la cascada desconectada. Una vez lograda la estabilización del lazo esclavo, se
pasa al ajuste del maestro.
160
6.8 Control de Relación
Lazo esclavo manual: Se puede actuar sobre la válvula (abrir o cerrar) desde
el regulador. La cascada está desconectada. Para pasar a este modo se debe
inicializar el modo integral del lazo interno en caso de que éste lo posea.
Lazo esclavo automático: El lazo interno se cierra, la referencia de este lazo
se puede variar normalmente. Para pasar a este modo es conveniente hacer
previamente la referencia igual a la salida del sub-proceso 4.
Lazo maesto en manual, esclavo en automático: En este modo, la salida
del regulador maestro se hace igual a la referencia del regulador esclavo. Se
puede actuar sobre la variable controlada principal, moviendo el valor de la
referencia, pero en este modo no necesariamente coinciden ambos valores.
Lazo en cascada: Se cierra el regulador maestro. Previamente debe hacerse
el valor de la referencia 2 igual a la variable controlada principal, para evitar
un cambio brusco de la posición de la válvula.
SPv1 v1
R v2 + B R G1
161
6. TÉCNICAS DE CONTROL DISCRETO
v1 = Rv2 + B (6.63)
si B = 0, entonces v1 = Rv2 .
El operador puede modificar el valor de los parámetros R y B. El control de relación
tiene en común con el control en cascada que la referencia de un regulador es calculado,
en este caso a través de un bloque lineal. En el caso de las cascadas, la referencia de
uno de los reguladores se calcula también, pero mediante un segundo regulador.
El regulador de relación tiene generalmente 4 modos de operación:
4. Ratio: Se calcula el SPv1 a través del bloque de relación. Debe inicializarse esta
conexión haciendo:
v1m
B=0 R= (6.64)
v2m
162
6.9 Conexiones Feed-Forward
v b
FF G2
- +
yr + + + y
R G1 b
∆y G2 (z) F F G1 (z)
= + (6.65)
∆v 1 + R(z)G1 (z) 1 + R(z)G1 (z)
∆y G2 (z) + F F G1 (z)
= (6.66)
∆v 1 + RG1 (z)
de donde se deduce:
G2 (z)
FF = (6.67)
G1 (z)
Si se aproximan G1 (z) y G2 (z) como dos funciones transferenciales de 1er orden
con retardo:
K1 z −p1
G1 (z) = (6.68)
1 − α1 z −1
K2 z −p2
G2 (z) = (6.69)
1 − α2 z −1
entonces:
K2 z −p2 (1 − α1 z −1 )
F F (z) = (6.70)
K1 z −p1 )(1 − α2 z −1 )
donde:
α1 = e−Ts /T1 , α2 = e−Ts /T2 (6.71)
163
6. TÉCNICAS DE CONTROL DISCRETO
finalmente:
∆u 1 − α1 z −1
= F F (z) = −KF F z (p1 −p2 ) (6.72)
∆v 1 − α2 z −1
K2
KF F = (6.73)
K1
En muchas ocasiones es suficiente sólo una compensación del valor estacionario y
entonces se hace F F = KF F .
Si se utiliza una compensación dinámica tal como lo expresamos mediante (6.72),
entonces hay que tomar en consideración la siguiente condición: Para que la compen-
sación tenga un sentido real es necesario que p1 ≤ p2 , es decir, el retardo en el canal
principal debe ser menor o igual que el retardo en el canal de la perturbación, pues
de otra manera habrı́a que utilizar predicciones de valores futuros de la perturbación.
Esto se comprende mejor observando la ecuación de diferencia asociada:
Como es conocido, uno de los factores que más afecta la calidad y estabilidad de
un sistema de control, ya sea este basado en un algoritmo PID u otro cualquiera,
es la presencia de “retardo puro” o “retardo por transporte” en el canal de control.
Muchas soluciones han sido intentadas para compensar el efecto negativo del retardo,
destacándose entre ellas por su sencillez y efectividad el llamado “Predictor de Smith”
conocido desde finales de los años 50. El uso generalizado de este concepto de control
tropezó, no obstante, con la dificultad de su realización mediante la tecnologı́a analógica
disponible en los años 50 y 60.
164
6.10 Predictor de Smith
-
Sp1 u1 + y1
R1 b
G11 b
+ +
- +
b
G12
F F1 F F2
b
G21
- +
Sp2 + + y2
R2 b
G22 b
+ u2
-
Ke−θs
G(s) = (6.76)
Ts + 1
165
6. TÉCNICAS DE CONTROL DISCRETO
+
yr Kp ypi + yp
R(s) b
Tp s+1 e−θp s +
b
+ +
- -
Km
ymi ym -
Tm s+1
b
e−θm s
en el modelo. La conexión interna es la que tiende a compensar el efecto del retardo, pues
la acción de control se anticipa como resultado de la predicción de la salida mediante
el modelo utilizado.
La función transferencial en lazo cerrado correspondiente a la estructura de la figura
6.10, resulta:
R(s)Kp −θp s
Yp (s) Tp s+1 e
= R(s)Kp −θp s
(6.77)
Yr (s) + R(s)K m R(s)Km −θm s
1+ Tp s+1 e Tm s+1 − Tm s+1 e
θp = θm (6.78)
Kp = Km (6.79)
Tp = Tm (6.80)
RKp
Yp (s) Tp s+1
= RKp
e−θp s (6.81)
Yr (s) 1 + Tp s+1
166
6.11 PID Autoajustable
En esta sección, sin pretender agotar el tema, lo que no serı́a posible dada la diver-
sidad de enfoques propuestos, presentaremos dos variantes totalmente distintas en su
concepción y que tienen en común su simplicidad y facilidad de realización.
Este método utiliza el criterio clásico de margen de fase como criterio de diseño
(Banyasz y Keviczky, 1982) y parte de la descripción del proceso mediante la función
de transferencia discreta:
Y (z) b1 z −(k+1)
Gp (z) = = (6.82)
U (z) 1 − a1 z −1 − a2 z −2
p0 − p1 z −1 − p2 z −2 P (z)
Gr (z) = −1
= (6.83)
1−z 1 − z −1
donde:
p1
p′1 = = a1 (6.88)
p0
p2
p′2 = = a2 (6.89)
p0
167
6. TÉCNICAS DE CONTROL DISCRETO
b1 p0 z −(k+1) b1 p0 −k
L(z) = Gp (z)Gr (z) = −1
= z −1 z (6.90)
1−z 1 − z −1
Nótese que se trata de la función de transferencia de un integrador en serie con un
elemento de retardo puro y con una ganancia:
KI = b1 p0 (6.91)
donde:
1 1 − e−T s
Lf (s) = (6.93)
T s
corresponde a un fijador de orden cero y T es el perı́odo de muestreo.
KI 1
LI (s) = (6.94)
T s
y
Lk (s) = e−sT k (6.95)
son las funciones de transferencia del integrador y del elemento de retardo puro.
⇒ Tarea 6.2
A partir de (6.93), (6.94) y (6.95), haciendo s = jw y con algunas transformaciones
simples, pueden determinarse las caracterı́sticas de amplitud y fase vs. frecuencia del
sistema en lazo abierto. Ellas son:
KI sen wT2
A(w) = wT
(6.96)
wT 2
y
wT π
φ(w) = − (2k + 1) − (6.97)
2 2
De acuerdo con el criterio de Nyquist, el sistema en lazo cerrado correspondiente
tendrá un determinado margen de fase φa , si se cumple:
φ(w) = −π + φa (6.98)
168
6.11 PID Autoajustable
para
A(w) = 1 (6.99)
wT
Utilizando ahora los resultados (6.96) y (6.97), con x = 2 , se obtiene de (6.99):
2x2 x
KI = = 2x (6.100)
|sen(x)| |sen(x)|
x
El cociente |sen(x)| puede aproximarse muy bien por 1 para valores de x pequeños
y por tanto:
KI ≈ 2x (6.101)
π − 2φa
x= (6.102)
2(2k + 1)
Sustituyendo (6.102) en (6.101) obtenemos entonces:
p2 = â2 p0 (6.107)
donde por b̂1 , â1 , â2 deben interpretarse como las estimaciones de estos parámetros
obtenidas en tiempo real mediante cualquier variante del método de mı́nimos cuadrados
recursivos.
De acuerdo con (6.83), la ley de control se calcula mediante:
169
6. TÉCNICAS DE CONTROL DISCRETO
y(t) b1 z −1
Gp (z −1 ) = = (6.109)
u(t) 1 − a1 z −1 − a2 z −2
r1 + r2 z −1 + r3 z −2
Gr (z −1 ) = (6.110)
1 − z −1
a1 − a2 + b1 r2 −2 b1 r1 − a1 − 1 −1 1
z −3 + z + z + =0 (6.111)
b1 r3 + a2 b1 r3 + a2 b1 r3 + a2
Imponemos ahora la condición de que el sistema en lazo cerrado tenga un polo triple,
de valor a, que como es lógico, debe estar en la zona de estabilidad. Esto equivale a
suponer que el polinomio caracterı́stico del sistema tenga la forma:
1 3
r1 = (1 + a1 − ) (6.113)
b1 a
1 3
r2 = (a2 − a1 + 2 ) (6.114)
b1 a
1 1
r3 = − (a2 + 3 ) (6.115)
b1 a
170
6.11 PID Autoajustable
Ki = r1 + r2 + r3 (6.120)
Kd = r3 (6.121)
171
6. TÉCNICAS DE CONTROL DISCRETO
o también:
0≤κ≤1 (6.124)
Ti = βTc (6.128)
γ
Td = Tc (6.129)
β
donde Tc es el perı́odo de control.
En el trabajo de Marsik se da el siguiente procedimiento recursivo para el cálculo
del coeficiente α:
0,5
∆α(t) = 0,5fv α(t) −1 (6.130)
κ(t)
que debe converger al valor óptimo de α cuando κ ≈ 0,4.
Un problema clave del algoritmo es el cálculo en tiempo real del ı́ndice de oscilación
κ(t). Para ello fue necesario acudir a una vieja expresión dada por Rice (1944) que
calcula el número de cruces por un determinado nivel de una señal gaussiana continua:
s
1 σv2
fe = (6.131)
π σe2
donde σv2 y σe2 son respectivamente las varianzas del error de regulación y de su primera
derivada. En una forma análoga:
s
1 σa2
fv = (6.132)
π σv2
172
6.11 PID Autoajustable
t
1 X
σv2 = ∆e(i)2 (6.134)
t−1
i=1
t
1 X 2
σa2 = ∆e e(i)2 (6.135)
t−1
i=1
173
6. TÉCNICAS DE CONTROL DISCRETO
s
1 a2 (t)
fv (t) = (6.142)
π v 2 (t)
v 2 (t)
κ(t) = q (6.143)
a2 (t)e2 (t)
Para determinar las reglas de adaptación de los parámetros γ(t) y β(t), vamos a
escribir la ecuación (6.125) en la siguiente forma ligeramente modificada:
t
X
u(t) = [α(i)(e(i) + β(i)∆e(i) + γ(i)∆2 e(i))] (6.145)
i=0
De acuerdo con Marsik (1983), la adaptación de los parámetros γ(i) y β(i), puede
hacerse en una forma eficiente imponiendo la condición de que la energı́a efectiva de
control aportada por las tres componentes del regulador sea igual, lo cual se traduce a
la condición: q q q
β(i) ∆e(i)2 = γ(i) ∆2 e(i)2 = e(i)2 (6.146)
El algoritmo se usa como una herramienta opcional para el ajuste del algorit-
mo PID posicional descrito mediante la ecuación (6.1), usando las relaciones de
transformación dadas en (6.127), (6.128) y (6.129) y no como un algoritmo autóno-
mo de control como fue concebido por el autor. Entonces, dado un conjunto de
parámetros del PID, es necesaria una inicialización adecuada de los parámetros
α, β y γ que garantice una transición suave cuando la opción de auto-ajuste se
conecte.
174
6.12 Conclusiones
El comportamiento del lazo cerrado cuando se usan los parámetros del PID cal-
culados mediante el algoritmo de Marsik, resulta en algunos casos excesivamente
oscilatorio. Esto se debe, en nuestra opinión, a la relación demasiado rı́gida que
se establece entre los tiempos integral y derivativo calculados. Hemos obtenido
un comportamiento más flexible mediante la introducción de un nuevo parámetro
en el cálculo del tiempo integral:
T i = ξβTc (6.148)
6.12. Conclusiones
En este capı́tulo hemos visto con detalle los PID discretos haciendo énfasis en las
consideraciones prácticas que se deben de tomar para garantizar el éxito de la aplicación,
también hemos tratado algunas variaciones en el lazo de control como lo son el control
en cascada, las conexiones feed-forward y el predictor de Smith los cuales estudian
algunas alteraciones de los sistemas de control en la práctica y como tratarlos para
disminuir su efecto en el controlador PID. Finalmente se trata con PID auto ajustable
los cuales con apoyo con el método de mı́nimos cuadrados se logra una forma de control
adaptable.
En los últimos años numerosos enfoques han sido utilizados con el fin de encon-
trar un procedimiento de auto-ajuste suficientemente confiable para los parámetros de
los reguladores PID. Distintos criterios de ajuste se han empleado con mayor o menor
éxito: garantizar un margen de fase prefijado, uso de la regla de Ziegler-Nichols par-
tiendo de oscilaciones provocadas en el proceso, minimización de un criterio cuadrático,
asignación de polos, auto ajuste basado en razonamientos heurı́sticos que han sido los
175
6. TÉCNICAS DE CONTROL DISCRETO
que han arrojado mejores resultados. En numerosas simulaciones el algoritmo PID auto
ajustable mediante la asignación de un polo triple en la función de transferencia de lazo
cerrado, ha mostrado su eficacia para la regulación de sistemas de hasta tercer orden.
El algoritmo de auto-ajuste, con algunas modificaciones, ha sido probado exitosa-
mente en una gran variedad de ejemplos simulados, incluyendo modelos de hasta cuarto
orden con retardo, ası́ como en casos prácticos reales. En el transcurso de este capı́tulo
tenemos herramienta de gran utilidad para controlar sistemas dinámicos que cambian
con el tiempo y las condiciones del medio ambiente.
176
7
Conclusiones Generales
177
7. CONCLUSIONES GENERALES
la identificación de sistemas, sobre este tema hay mucho trabajo por realizar. Con lo
aprendido hasta aquı́ podemos aplicar los métodos clásicos de la ingenierı́a de control
automático.
178