Observador de Estado Control III
Observador de Estado Control III
Observador de Estado Control III
Narvis Camacho
Reimer Bolivar
Andres Marquez
Variable de estado:
Es una cantidad de variables que describe completamente el estado actual de un sistema
dinámico
Figura 1
Figura 2
Observabilidad
Observador de Estado
(1.2)
(1.3)
(1.4)
(1.5)
El problema de diseñar un observador de orden completo está en determinar la matriz de ganancias del observador
Ke de forma que la dinámica de error definida mediante la Ecuación (1.5) sea asintóticamente estable con una
velocidad de respuesta suficiente.
(La estabilidad asintótica y la velocidad de respuesta de la dinámica de error se determinan mediante los valores
característicos de la matriz A.KeC.) Por tanto, el diseño del observador de orden completo se convierte en
determinar un Ke apropiado tal que A .KeC tenga los valores propios deseados.
sea de estado completamente controlable. La condición de controlabilidad completa del estado para este sistema
dual es que el rango de
sea n.
la ecuación para el observador de estado de orden completo es
Observador de Luenberger
El observador de Luenberger permite estimar el estado interno de un sistema dinámico lineal
utilizando una aproximación basada en un sistema clon y correcciones de entrada.
Esquema de realimentación lineal de salida para un sistema no lineal, utilizando un observador dinámico de estado
Convergencia
Tomaremos como nuestro observador de estados, o reconstructor, al sistema dinámico lineal
dado por
(4.1)
donde la matriz Φ y los vectores Γ y L serán determinados en lo sucesivo, de tal manera que sea posible
obtener un estimador asintótico de los estados del sistema lineal.
(4.2)
con valor inicial Puesto que el valor inicial del estado verdadero
es desconocido, no podemos hacer suposición alguna acerca del valor inicial del error de estimación. En
general tendremos que suponer que tal error inicial es diferente de cero pero no podremos dar
calificaciones adicionales sobre la magnitud de esta cantidad vectorial.
(4.3)
La dinámica del error de observación (4.3) no debe depender ni de los valores de las acciones de control ni de los
valores correspondientes que éstas producen en los estados del sistema.
Si hemos de manipular uδ para lograr, digamos, la estabilización del sistema linealizado hacia el origen de coordenadas,
no parece lógico que el error de observación dependiera de cuan grandes o pequeñas deban ser tales acciones de control.
El mismo razonamiento nos conduce a concluir que es igualmente inconveniente que el error de estimación dependa de
los valores que va tomando el vector de estado en su camino hacia la estabilización, pues la existencia de sus transientes
alejaría, aún más, el valor del estado estimado del real, actuando como perturbaciones en la dinámica
del error. En consecuencia, debemos anular, con una escogencia apropiada de las matrices y vectores Φ, Γ y L, la
influencia de los estados y controles del sistema original linealizado.
De esta forma, el error de estimación eδ solo dependerá de si mismo y, en particular, de su valor inicial.
La solución de la ecuación diferencial lineal vectorial de primer orden, la cual representa al vector de error de
estimación, estaría dada aún en términos del vector columna desconocido L mediante la siguiente expresión:
En virtud del resultado anterior, estamos en capacidad de predecir el comportamiento cualitativo de eδ, el cual depende
principalmente de la naturaleza de los autovalores de la matriz Φ = A − LC.
Nuestro objetivo será entonces el de imponer, sobre la dinámica autonoma del error de estimación, un comportamiento
asintóticamente estable a cero a través de la escogencia juiciosa del vector L. Es decir, debemos escoger L de tal
manera que los autovalores de A − LC se encuentren en el semiplano izquierdo del plano complejo.
Al hacer esto se garantiza, finalmente, que el valor del estado estimado converge, efectivamente, de manera asintótica
al valor del estado verdadero del sistema linealizado.
el mismo. A partir de (4.2) y (4.4), tenemos
(4.6)
Si hacemos
Es decir, si damos el nombre de salida estimada al valor que produce el estado estimado cuando es afectado
por el mapa de salida C, entonces podemos reescribir las ecuaciones del observador (4.6)
(4.7)
(4.12)
El sistema de ecuaciones (4.12) lo podemos reescribir en una sola ecuación vectorial de estado aumentado o estado
compuesto de la planta y del observador, de la manera siguiente
(4.13)
Es difícil ver el efecto de los vectores (o matrices) de diseño K y L en la ecuación lineal (4.13) pues los autovalores de
la matriz ampliada en lazo cerrado no son fáciles de discernir.
Para resolver este problema, aprovecharemos el hecho de que dos matrices similares, es decir, aquellas que se
relacionan mediante una transformación similar, tienen los mismos autovalores.
(4.14)
El sistema transformado a partir de (4.14) resulta entonces:
(4.15)
Observador de Orden Reducido
Supóngase que el vector de estado x es un vector de dimensión n y que el vector de salida y es un vector de dimensión
m medible. Como las m variables de salida son combinaciones lineales de las variables de estado, no necesitan
estimarse m variables de estado, sino sólo n .m variables de estado. Así, el observador de orden reducido se vuelve un
observador de (n .m)-ésimo orden. Tal observador de (n . m)-ésimo orden es el observador de orden mínimo. La Figura
siguiente muestra el diagrama de bloques de un sistema con un observador de orden mínimo.
Sin embargo, es importante considerar que, si la medida de las variables de salida contiene ruido significativo y es
relativamente imprecisa, el uso del observador de orden completo puede dar un mejor comportamiento del sistema.
Para presentar la idea básica del observador de orden mínimo, sin complicaciones matemáticas innecesarias, se
expondrá el caso en el que la salida es un escalar (es decir, m =1) y se derivará la ecuación de estado para el observador
de orden mínimo.
Sea el sistema:
(4.2)
(4.21)
donde el vector de estado x se particiona en dos partes (un escalar) y [un vector de dimensión (n . 1)]. Aquí la
variable de estado es igual a la salida y de modo que se mide directamente, y es la parte no medible del vector de
estado. De esta forma, el estado particionado y las ecuaciones de salida son
(4.23) (4.24)
Donde:
A partir de la Ecuación (4.23), la ecuación para la parte medida del estado es
O bien
(4.25)
Los términos del miembro izquierdo de la Ecuación (4.25) se pueden medir. La Ecuación (4.25) funciona como la
ecuación de salida. A1 diseñar el observador de orden mínimo, se considera que el lado izquierdo de la Ecuación
(4.25) contiene cantidades conocidas. Por tanto, la ecuación (4.25) relaciona las cantidades medibles y no medibles
del estado.
(4.26)
Si se considera que los términos son cantidades conocidas, la Ecuación (4.26) describe la dinámica de la
parte no medida del estado.
A continuación se presentará un método para diseñar un observador de orden mínimo. El procedimiento de diseño se
simplifica si se utiliza la técnica de diseño desarrollada para el observador de estado de orden completo.
Compárese la ecuación de estado para el observador de orden completo con la del observador de orden
mínimo. La ecuación de estado para el observador de orden completo es
El diseño del observador de orden mínimo se realiza del modo siguiente. Primero, obsérvese que la
ecuación del observador para el observador de orden completo se obtuvo a partir de la Ecuación:
(4.27)
(4.29)
Se define
y (4.30)
(4.31)
Se define
(4.32)
La Ecuación (4.32), junto con la Ecuación (4.30), definen el observador de orden mínimo como
A continuación se obtendrá la ecuación del error del observador. Utilizando la Ecuación (4.25), la
Ecuación (4.28) se modifica a
(4.34)
Sistema con realimentación de estado observado donde el observador es el observador de orden mínimo.
Restando la Ecuación (4.34) de la Ecuación (4.26), se obtiene
(4.35)
Se define
(4.36)
Esta es la ecuación del error para el observador de orden mínimo. Obsérvese que e es un vector
de dimensión (n - 1).
La dinámica del error se puede escoger como se desee siguiendo la técnica desarrollada para
el observador de orden completo, siempre y cuando el rango de la matriz
(4.37)
donde son valores propios deseados para el observador de orden mínimo. La matriz de
ganancias del observador Ke se determina seleccionando primero los valores propios deseados para el
observador de orden mínimo [es decir, colocando las raíces de la ecuación característica, Ecuación (4.37),
en las posiciones deseadas] y después empleando el procedimiento desarrollado para el observador de
orden completo con las modificaciones adecuadas.
debe modificarse a
(4.38)
donde Ke es una matriz de (n - 1) x 1 y
debe modificarse a
Donde
Observador de Orden Completo Vs Observador de Orden Mínimo
Donde
x2 = X
(4.4)
(4.41)
Linealizando:
(4.42)
(4.43)
(4.44)
(4.45)
(4.46)
(4.47)
No existe, entonces, problema alguno en proponer un valor de la ganancia L1 de tal forma que el error de
estimación sea asintóticamente estable a cero. De hecho, cualquier valor positivo
de tal ganancia lo lograría. Sin embargo, este estimador aún requiere del cálculo de la derivada temporal
de la altura incremental en el segundo tanque (recordemos que
(4.49)
(4.50)
Lo cual reescribimos:
(4.51)
(4.52)
Observador de orden reducido caso general
Esquema de control realimentado lineal de la salida para un sistema de tanques mediante el uso de un observador de orden
reducido
Podemos, entonces, particionar las matrices que intervienen en el sistema linealizado y reescribir
estas ecuaciones de estado incremental como:
(4.60)
(4.61)
Reescribamos las ecuaciones de estado linealizadas (4.38) en función de las variables conocidas (asumiremos como
antes que y ̇δ es medible):
Puesto que yδ y uδ son conocidas, es concebible que también podamos conocer como señal del tiempo al vector n1
dimensional dado por la segunda ecuación anterior. Designemos a este vector como y2δ:
(4.62)
El problema de estimación de x2δ puede verse como el problema de reconstruir el estado del subsistema dado por:
donde el vector v1δ = A21x1δ +B2uδ es completamente conocido, al igual que la señal y2δ.
La reconstrucción del estado incremental desconocido x2δ, con error de estimación asintóticamente estable a cero, es
factible si el par (A22, A12) es observable o, en su defecto, simplemente detectable.
Teorema 4.1: Un sistema lineal z ̇ = Az, y = Cz es observable (es decir, el par (A, C) es observable) si, y solamente si,
el único vector η para el cual se cumple que
es para η = 0.
Prueba: Supongamos, por un momento, que el sistema lineal dado es inobservable. Entonces, existe una
transformación del vector de estado z que exhibe tal inobservabilidad. Tal transformación lleva al sistema a su forma
canónica inobservable de Kalman:
Puesto que la exponencial matricial de este sistema está dada por:
entonces es claro que existen vectores η distintos de cero que hacen cero el siguiente producto:
(4.63*)
es decir:
(4.64)
Si hacemos
(4.65)
Entonces obtenemos
(4.66)
Estructura de control por realimentación lineal de la salida, a base de un observador dinámico de orden reducido
El vector de estado x(t) es de dimensión n y la salida y(t) es un escalar medible que se expresa como una
combinación lineal de las variables de estado (y(t)=Cx(t)). Por lo que no necesitan estimarse las n
variables de estado, sino sólo n-1. Así, el observador de orden reducido se vuelve un observador de (n-1)-
ésimo orden