DISEO DE OBSERVADORES DE ESTADO

CONCEPTO

TIPOS

1. OBSERVADOR DE ORDEN COMPLETO

DISEO DE OBSERVADORES DE ORDEN COMPLETO

1.1. METODO DE DISEO ABREVIADO

1.2. METODO DE DISEO POR LA FORMULA DE ACKERMAN

1.3. METODO DE DISEO COMPLETO

1.4. DISEO MEDIANTE EL SOFTWARE MATLAB

2. OBSERVADOR DE ORDEN REDUCIDO

2.1. DISEO DE OBSERVADORES DE ORDEN REDUCIDO

2.2. METODOLOGA DE DISEO

3. OBSERVADOR PARA SISTEMAS MIMO

REFERENCIAS BIBLIOGRFICAS

REFERENCIAS DE INTERNET

CONCEPTO:

Los observadores de estado, son herramientas virtuales, que permiten estimar las variables o estados de un
sistema en base a mediciones de las seales de salida y seales de control. Estos observadores permiten enviar
informacin estimada acerca del valor que toman dichos estados, permitiendo conocer un aproximado del
valor real, adems cuentan con muy poco margen de diferencia o error.

Se le considera una herramienta virtual, puesto que se desarrolla como software o programa dentro de una
computadora.

TIPOS:

Existen 2 tipos de observadores: observadores de orden completo, y observadores de orden reducido u orden
mnimo.

Los observadores de orden Completo, son aquellos utilizados para observar o estimar todos los
estados de un sistema.

Los observadores de orden Reducido, son aquellos utilizados para observar o estimar solo algunos
estados de un sistema.

1
1. OBSERVADOR DE ORDEN COMPLETO

Dado el sistema:

donde:

x Vector de estado (n x 1)

u Seal de control (escalar)

y Seal de salida (escalar)

A Matriz (n x n)

B Matriz (n x 1)

C Matriz (1 x n)

D Matriz (escalar)

Se puede estimar sus estados mediante la siguiente expresin:

donde:

L Vector de ganancias que permiten la observacin de estados (1 x n)

Vector de estados estimados

Salida estimada

2
Debe notarse que las matrices A, B, C, D son las mismas tanto para un sistema real como para el sistema
estimado. Para los clculos siguientes se asume que el valor de D es cero.

La diferencia existente entre x y

Para determinar el error de observacin restamos

Si se sabe que el error esta definido como la diferencia entre el estado real y el estado estimado, entonces se
tendr:

3
A partir de esta expresin se puede conocer el comportamiento dinmico y la estabilidad del sistema, si la
matriz |ALC| es estable, entonces el observador har bien su trabajo, y dada cualquier condicin inicial, el
sistema tender a un error cero.

La eleccin de correctos valores para el vector de observabilidad L, permitir que el comportamiento

La estabilidad asinttica y la velocidad de respuesta de la dinmica del error se determina mediante los
autovalores de la matriz |ALC|, dados por el polinomio caracterstico |sIA+LC|.

Existe una condicin necesaria, la cual consiste en que el sistema obtenido sea estable y completamente
controlable y observable.

Ejemplo:

Determinar la ecuacin caracterstica del sistema siguiente, si se le agrega un observador de estados L.

Solucin.

Si el sistema es de orden 2, es de suponer que el observador tambin ser de orden 2, y lo podemos definir
como

Luego el polinomio caracterstico estar dado por:

DISEO DE OBSERVADORES DE ORDEN COMPLETO

1.1. METODO DE DISEO ABREVIADO

Analizando la respuesta del ejemplo anterior nos podemos dar cuenta que los valores que toman L1 y L2 estn
condicionadas por las races del polinomio, las cuales a su vez estn condicionadas por las caractersticas con
que queremos que cuente el sistema, por tanto se puede elegir races de tal modo de poder controlar la
respuesta del sistema en lazo cerrado.

Por lo tanto podemos asumir valores para dichas races, a los que llamaremos y de modo tal que el
polinomio tenga una respuesta estable. Luego por simple equivalencia de trminos podemos hallar el valor de
las incgnitas.

Ejemplo:

4
Dado el polinomio caracterstico del ejemplo anterior: s2+(4+L2)s+(2+L1), encontrar el valor de L1 y L2 si
se quiere que los polos deseados del sistema se ubiquen en 4 y 3.

Solucin.

Las races del polinomio son = 4 y = 3

Luego por equivalencia |sI A + LC| = s2 + 7s + 12

Es decir, s2 + (4 + L2)s + (2 + L1) = s2 + 7s + 12

s2 = s2

donde (4 + L2)s = 7s L2 = 3

(2 + L1) = 12 L1 = 10

Podemos generalizar la metodologa seguida anteriormente, de la siguiente manera:

Si tenemos que [ n] son los autovalores deseados para la matriz del observador |ALC|, estos
conforman el polinomio caracterstico:

(s) (s) (s) (sn)

Este polinomio se iguala al polinomio caracterstico original |sIA+LC|, crendose una equivalencia entre
trminos:

|sIA+LC| = (s) (s) (s) (sn)

Resolviendo la equivalencia se podr encontrar el valor del vector L.

NOTA: Este mtodo esta restringido a sistemas de hasta 3er orden, adems el sistema debe estar en la forma
cannica observable.

Es aconsejable que los polos del observador sean de 3 a 5 veces mayores (ms negativos) que los polos del
controlador por realimentacin de estados, pero sin salirse de la regin de estabilidad dada por el lugar
geomtrico de las races. La eleccin de los polos deseados van a determinar las caractersticas de la respuesta
obtenida, por lo que puede existir un conjunto infinito de vectores L como solucin, de las cuales solo un
limitado nmero de soluciones cumplen con las necesidades requeridas para el sistema (como por ejemplo:
sobreimpulso, velocidad de respuesta, etc. del sistema estimado), por lo que se aconseja probar, mediante
simulacin, la respuesta del sistema a diferentes valores de polos escogidos.

1.2. METODO DE DISEO POR LA FORMULA DE ACKERMAN

La formula de Ackerman aplicada al diseo de observadores de estado, esta dada por:

5
En donde
es equivalente a
, que es el polinomio caracterstico deseado, pero en vez de la s se coloca la matriz A.

Ejemplo:

Para el siguiente sistema, determinar el vector de observadores de estados L, si se quiere que los polos
deseados se ubiquen en 3+j y 3j

Solucin.

Si (s) = (s) (s) = (s+3j) (s+3+j) = s2 + 6s + 10

por lo tanto (A) = A2 + 6A + 10I

1.3. METODO DE DISEO COMPLETO

1) Determinar la controlabilidad del sistema y la observabilidad

Controlabilidad:

Observabilidad:

2) Calcular el polinomio caracterstico original |sIA|, el cual ser :

|sIA| = sn + a1sn1 + a2sn2 + + an1s + an = 0

6
3) Es conveniente trabajar con las ecuaciones de estado en su forma cannica observable, si no se encuentra
en esta forma, se debe determinar una matriz de transformacin para llevarla a esta forma, la cual se define
como:

Q = (W x Wo)1

En donde Wo es la matriz de observabilidad, y W se define como:

En donde a1, a2, an2, an3, son los coeficientes del polinomio caracterstico original |sIA|.

4) Se determina el polinomio caracterstico deseado a partir de (s) (s) (s) (sn), donde i es un
polo deseado, obtenindose:

sn + b1sn1 + b2sn2 + + bn1s + bn

5) Finalmente el vector L se encuentra a partir de la siguiente expresin: (*)

Ejemplo:

Para el siguiente sistema, determinar el vector de observadores de estados L, si se quiere que los polos
deseados se ubiquen en 5, 2+j y 2j

Solucin.

1)

Controlabilidad Observabilidad

rank(Wc) = 3 rank(Wo) = 3

det(Wc) = 1 det(Wo) = 1

7
El sistema es controlable El sistema es observable

2)

3) Q = ( W x Wo )1

4)

(s ) (s ) (s )

= (s + 5) (s + 2 + j) (s + 2 j)

= (s + 5) (s2 + 4s + 5)

= s3 + 9s2 + 25s + 25 " s3 + b1s2 + b2s + b3

5)

1.4. DISEO MEDIANTE EL SOFTWARE MATLAB

Se puede hacer uso del software Matlab, para lo cual se emplea el comando acker o el comando place.

Dado el sistema

, y un vector de polos deseados:

Se puede obtener un observador de estados utilizando:

L = place (A',B',P)' o tambin

L = acker (A',B',P)'

8
Ejemplo:

Para el siguiente sistema, determinar el vector de observadores de estados L, si se quiere que los polos
deseados se ubiquen en 2, 1+j y 1j

Solucin.

>> A = [0 1 0; 0 0 1; 3 2 1];

>> C = [2 0 0];

>> P = [2 1+j 1j];

>> L = acker(A',C',P)'

L=

1.5000

0.5000

3.0000

>> L = place(A',C',P)'

L=

1.5000

0.5000

3.0000

Ejemplo:

Para el siguiente sistema, determinar el vector de observadores de estados L aplicando los 4 mtodos antes
descritos, si se quiere que los polos deseados se ubiquen en 2, 3+0.5j y 30.5j

Solucin.

9
1) METODO COMPLETO

Controlabilidad Observabilidad

rank(Wc) = 3 rank(Wo) = 3

det(Wc) = 16 det(Wo) = 1

El sistema es controlable El sistema es observable

Polinomio caracterstico original

Q = ( W x Wo )1

Polinomio caracterstico deseado

(s ) (s ) (s )

10
= (s + 2) (s + 3 + 0.5j) (s + 3 0.5j)

= (s + 2) (s2 + 6s + 9.25)

= s3 + 8s2 + 21.25s + 18.5 " s3 + b1s2 + b2s + b3

Vector Observador

2) METODO ABREVIADO

Polinomio deseado

Luego por equivalencia s3 + (L3)s2 + (1 + L2)s + (4 + L1) = s3 + 8s2 + 21.25s + 18.5

(2+L3)s2 = 8s2 L3 = 6

donde (1 + L2)s = 21.25s L2 = 20.25

(4 + L1) = 18.5 L1 = 14.5

3) METODO POR FORMULA DE ACKERMAN

Si (s) = (s) (s) (s ) = s3 + 8s2 + 21.25s + 18.5

11
por lo tanto (A) = A3 + 8A2 + 21.25A + 18.5I

4) USANDO MATLAB

>> A = [0 0 4; 1 0 1; 0 1 2];

>> C = [0 0 1];

>> P = [2 3+0.5j 30.5j];

>> L = acker(A',C',P)'

L=

14.5000

20.2500

6.0000

>> L = place(A',C',P)'

L=

14.5000

20.2500

6.0000

2. OBSERVADOR DE ORDEN REDUCIDO

12
En la prctica no todas las variables necesitan ser observadas, habr algunas que se podrn medir directamente
y con buena precisin, por tanto no ser necesario un observador que estime todos los estados, sino ms bien
solo algunos de ellos.

Si se cuenta con un vector de estados X de dimensin (n x 1) del cual m estados pueden ser medibles, se
tendr que el orden del observador ser (nm x 1).

El vector X puede ser dividido en 2 vectores:

Xa que corresponde a los estados medidos, de orden (m x 1)

X1 Xm Estados conocidos o medibles

Xm+1 Xn Estados no conocidos que requieren ser observados

Obtenindose:

13
Las dimensiones de las submatrices son:

Aaa m x m

Aab m x nm

Aba nm x m

Abb nm x nm

Ba m x 1

Bb nm x 1

Ca 1 x m

Cb 1 x nm

El sistema queda reducido a la siguiente expresin:

2.1. DISEO DE OBSERVADORES DE ORDEN REDUCIDO

En el diseo de observadores de orden completo, el sistema era descrito por

En cambio para el diseo de observadores de orden reducido el sistema es descrito por

la Ecuacin de Estado:

y por la Ecuacin de Salida:

De estos 2 sistemas se pueden establecer una serie de equivalencias:

Observador Orden Completo Observador Orden Reducido

Cabe mencionar que se busca encontrar un vector de observadores L, de orden (nm x 1)

En el diseo de observadores de orden completo se determino la siguiente ecuacin:

14
La cual puede llevarse a su correspondiente equivalencia para el caso de observadores de orden reducido,
obtenindose la siguiente representacin: (*)

y la Ecuacin del Error es:

La Ecuacin Caracterstica para el observador es la siguiente:

|sI Abb + [Link]| = si + a1si1 + a2si2 + + ai1s + ai = 0

donde i equivale al orden de L, es decir (nm)

2.2. METODOLOGA DE DISEO

Se pueden aplicar los mismos mtodos usados para hallar los observadores de orden completo, pero tenemos
que hacer una variacin la cual consiste en reemplazar los ndices: en vez de considerar el ndice n, se debe
considerar el ndice i. Por ejemplo si un sistema es de orden 4 y tiene un estado medible, entonces m=1, y por
tanto el sistema observado ser de orden 3 (n m = 4 1 = 3 = i).

El otro cambio que hay que hacer es realizar una equivalencia entre las matrices, dada de la siguiente manera:

Aaa D

Aab C

Aba B

Abb A

Lo que resta por hacer, es aplicar la misma metodologa que para los observadores de orden completo,
considerando el nuevo orden del sistema (i) y las nuevas matrices del sistema (Abb, Aba, Aab, Aaa).

Ejemplo:

Para el siguiente sistema, suponga que el estado X1 se puede medir con precisin, disee el observador de
orden reducido L aplicando los 4 mtodos antes descritos, si se quiere que los polos deseados se ubiquen en
3+0.5j y 30.5j

Solucin.

Nuevo orden para el sistema observado: 3 1 i = 2

15
Nuevas matrices del sistema:

D Aaa

C Aab

B Aba

A Abb

1) METODO COMPLETO

Controlabilidad Observabilidad

rank(Wc) = 2 rank(Wo) = 2

det(Wc) = 36 det(Wo) = 1

El sistema es controlable El sistema es observable

Polinomio caracterstico original

Q = ( W x Wo )1

Polinomio caracterstico deseado

16
(s ) (s )

= (s + 2 + 3.4641j) (s + 2 3.4641j)

= s2 + 4s + 16 " s2 + b1s + b2

Vector Observador

2) METODO ABREVIADO

Polinomio deseado

Luego por equivalencia

6 + L1 = 4 L1 = 2

6 L1 + 11 + L2 = 16 L2 = 17

3) METODO POR FORMULA DE ACKERMAN

Si (s) =

17
por lo tanto (A) " (Abb)= Abb2 + 4Abb + 16I

4) USANDO MATLAB

>> Aaa = [0];

>> Aab = [1 0];

>> Aba = [0; 6];

>> Abb = [0 1; 11 6];

>> P = [2 + 3.4641j 23.4641j];

>> L = acker(Abb', Aab',P)'

L=

17

>> L = place(Abb', Aab',P)'

L=

17

3. OBSERVADOR PARA SISTEMAS MIMO

Hasta el momento hemos visto el diseo de observadores para sistemas SISO (single input, single output), es
decir, sistemas que tienen una entrada y una salida. A continuacin vamos a estudiar el diseo de
observadores para el caso de sistemas con varias entradas y varias salidas, llamados sistemas MIMO (multiple
input, multiple output).

Un sistema MIMO puede ser descrito de la siguiente manera:

18
donde:

x Vector de estado (n x 1)

u Seal de control (1 x r)

y Seal de salida (escalar)

A Matriz (n x n)

B Matriz (n x r)

C Matriz (q x n)

D Matriz (q x r)

n Es el orden del sistema (numero de estados)

r Es el nmero de entradas

q Es el nmero de salidas

En el diseo de observadores, se sabe que el nmero de entradas que tenga el sistema, no afectar el diseo
del observador, puesto que estos datos no intervienen en los clculos.

En cambio, el nmero de salidas con que cuente el sistema si afecta el diseo del observador. En este caso se
aplica el Principio de Separabilidad Lineal, el cual consiste en separar las salidas y trabajarlas como si fueran
provenientes de sistemas distintos.

Ejemplo:

Si se tiene la siguiente matriz C (matriz de salidas)

Esta se podr descomponer bajo el Principio de Separabilidad Lineal, en dos submatrices:

Para cada una de estas nuevas matrices se requerir un observador independiente: L1 y L2, cuyo orden ser de
(n x 1), en este caso ser de (3 x 1) c/u. As mismo cada observador requiere de polos deseados, los cuales
pueden ser independientes o tambin comunes, siempre que cumplan con dar estabilidad al sistema y
dependiendo de las condiciones de respuesta que se desee obtener al estimar cada salida.

Ejemplo:

Para el siguiente sistema, determinar la matriz de observadores de estados L, si se quiere que los polos
deseados se ubiquen en 1 y 2 para la primera salida, y 1j, y 1+j para la segunda salida.

19
Solucin:

C1 = [1 0] C2 = [0 1]

METODO COMPLETO

Observabilidad

Controlabilidad rank(Wo) = 2

det(Wo) = 1

rank(Wc) = 2 El sistema es observable

det(Wc) = 1

El sistema es controlable
rank(Wo) = 2

det(Wo) = 5

El sistema es observable

Polinomio caracterstico original

DISEO DEL PRIMER OBSERVADOR (L1)

Q = ( W x Wo1 )1

20
Polinomio caracterstico deseado

(s ) (s )

= (s + 1) (s + 2)

= s2 + 3s + 2 " s2 + b1s + b2

Vector Observador

DISEO DEL SEGUNDO OBSERVADOR (L2)

Q = ( W x Wo2 )1

Polinomio caracterstico deseado

(s ) (s )

= (s + 1+j) (s + 1j)

21
= s2 + 4s + 5 " s2 + b1s + b2

Vector Observador

Se puede obtener los mismos resultados usando Matlab:

>> A = [0 1; 5 5];

>> B = [0; 1];

>> C = [1 0; 0 1];

>> C1 = C(1,:);

>> C2 = C(2,:);

>> P1 = [1 2] ;

>> P2 = [1+j 1j]

>> L1 = acker(A',C1',P1)'

L1 =

2.0000

7.0000

>> L2 = acker(A',C2',P2)'

L2 =

0.6000

3.0000

>> L = [L1 L2];

L=

22
2.0000 0.6000

7.0000 3.0000

REFERENCIAS BIBLIOGRFICAS

Chi Tsong Chen

Linear System. Theory and Design

Oxford University Press. 3ra Edicin. USA. 1998

Katsuhiko Ogata

Ingeniera de Control Moderna

Prentice Hall. 3ra Edicin. Mxico. 1998

Michael O'Flynn

Linear Systems. Time Domain and Transform Analysis

Wiley. Singapore. 1987

Donald M. Wiberg

Espacio de Estado y Sistemas Lineales

[Link]. Mxico. 1971

MathWorks Inc

MATLAB Users Guide, Toolbox deCcontrol.

Versin 5.3. USA. 1998.

Apuntes de clases de Teora de Sistemas Lineales.

Maestra en Ingeniera de Control y Automatizacin

Pontificia Universidad Catlica del Per

REFERENCIAS DE INTERNET

[Link]
[Link]
[Link]

(*) La deduccin de esta expresin puede encontrarse en el libro Ingeniera de Control Moderna, Katsuhiro
Ogata, 3ra Edicin, Cp.12, Pg. 817820.

(*) La deduccin de esta expresin puede encontrarse en el libro Ingeniera de Control Moderna, Katsuhiro

23
Ogata, 3ra Edicin, Cp.12, Pg. 832833.

24