Control Predictivo
Control Predictivo
Control Predictivo
El propósito de las unidades de procesos es, cómo no, producir una serie de productos
de forma que se maximicen los beneficios. En muchos casos, esta maximización viene
dictada bien por la máxima capacidad de producción o bien por las condiciones del
mercado en el que han de introducir los productos obtenidos.
A continuación se indican las ideas comunes a todos los controladores que se engloban
en el control predictivo:
- Utilización de manera explícita de un modelo para predecir la salida del proceso
en instantes de tiempo futuros, llamado horizonte temporal.
- Cálculo de las acciones de control minimizando una cierta función objetivo.
- Uso de una estrategia deslizante de tal modo que en cada instante el horizonte es
desplazado hacia delante en el tiempo.
Las diferencias más reseñables entre los controles predictivos radican en el modelo
utilizado para el proceso, las perturbaciones y la función objetivo a minimizar. Estos
contrastes, aunque pequeños, pueden causar grandes diferencias en el funcionamiento
en bucle cerrado.
Entre las ventajas del MPC frente a otros tipos de control caben destacar:
- Conceptos intuitivos, por lo tanto no es necesario personal muy cualificado en
control.
- Es aplicable a multitud de procesos, incluyendo aquellos con grandes retardos
(posee compensación implícita del retardo), fase no mínima o inestables.
- Fácilmente aplicable a procesos multivariables.
- Muy útil en procesos en los que son conocidas las referencias futuras.
Asimismo el MPC presenta algunos inconvenientes como son el aumento de carga de
cálculo, resuelta con los potentes ordenadores de hoy en día, y la necesidad de disponer
del modelo apropiado del proceso, el cual es crítico para las prestaciones del
controlador.
Como último concepto básico se expone a continuación la estrategia que siguen todos
los controladores de la familia MPC:
u(t+k|t)
u(t)
ŷ(t+k|t)
y(t)
N
1.- En cada instante t utilizamos el modelo del proceso para predecir las futuras salidas
para un determinado horizonte temporal u horizonte de predicción N. Se obtienen así las
salidas predichas ŷ(t+k|t) durante todo el horizonte, las cuales dependen de los valores
en el instante t y de las futuras acciones de control que se vayan a aplicar u(t+k|t) y que
hay que calcular.
2.- Las señales de control futuras se calculan de manera que se optimice un criterio
determinado para mantener el proceso lo más próximo posible a la trayectoria de
referencia w(t+k), que será o bien el set point o una aproximación suave.
3.- Se envía al proceso la señal de control u(t|t), desechándose el resto de las calculadas,
puesto que en el siguiente instante de muestreo ya se conoce y(t+1).
4.- Se repite el proceso desde el punto 1, con los valores de las señales actualizados.
Trayectoria
de referencia
Entradas y salidas
pasadas
Salidas predichas +
MODELO
-
Controles
futuros
Errores futuros
OPTIMIZADOR
Función
de coste Restricciones
obtienen tras aplicar un impulso como entrada al proceso. Sin embargo estos modelos
suelen truncarse de forma que nos quedamos con un modelo de la forma
N
y (t ) = ∑ hi ⋅ ∆u (t − i )
i =1
• Subspace Identification:
Se trata de una nueva tecnología de identificación paramétrica introducida por
DMCplus, que ofrece las siguientes ventajas:
- Utiliza un modelo en el espacio de estados para representar relaciones internas
entre variables.
- Especialmente indicada para procesos MIMO (Multiple In/Multiple Out),
pudiendo obtener la mínima parametrización y aumentando su eficiencia.
- Es un proceso no iterativo que usa álgebra lineal.
Su mayor inconveniente es su poca implantación hasta el momento en la industria, por
lo que no está demasiado depurado su funcionamiento.
g1 0 L 0
g g1 L 0
2
en donde M M O M se denomina matriz dinámica del ensayo en escalón y
G=
g m g m −1 L g1
M M O M
g p g p −1 L g p −m +1
lleva los coeficientes del modelo para todas las variables formada por m (horizonte de
∧
control) columnas, f es la respuesta libre del sistema, u son las señales de control e y
son las predicciones de salida (vector de dimensión p, número de variables del sistema).
j =1 j =1
(
u = G T G + λI ) −1
G T (w − f )
siendo ésta la acción de control a lo largo de todo el horizonte de control, pero hay que
recordar que sólo ha de aplicarse la primera de las acciones.
Además podemos añadir restricciones al problema del tipo matricial
R ⋅u ≤ c
aunque ésto nos lleva a soluciones no analíticas y sí numéricas.
Por último resaltar que las expresiones anteriores son fácilmente aplicables a sistemas
multivariables sin más que convertir las variables en vectores de forma que recojan
todas las señales de entrada y salida que intervengan en el proceso. La matriz G
quedaría entonces como:
G11 G12 L G1nu
G G 22 L G2 nu
G=
21
M M O M
Gny1 Gny 2 L Gnynu
donde cada submatriz Gij contiene los coeficientes de la respuesta ante escalón i-ésima