Oct 062016
 
CALIFICA EL POST

Para sistemas lineales e invariantes en el tiempo el análisis de estabilidad es una tarea relativamente simple, sin embargo los métodos de estabilidad vistos anteriormente, no pueden ser aplicados en el análisis de sistemas lineales variantes en el tiempo ni tampoco en sistemas no lineales.

 

Una forma genérica para analizar la estabilidad de sistemas es a través de la teoría de estabilidad de Lyapunov, donde dicho análisis es determinado de forma indirecta. Porque no se usa directamente la dinámica del sistema si no que por el contrario se analiza como el sistema se está comportando.

 

Esta idea surgió en el siglo 19 donde todo se basa en energía, en aquella época solo existían sistemas mecánicos, entonces si la energía de ese sistema mecánico va para cero al momento del tiempo tender a infinito, se podria caracterizar un sistema como estable.

 

Se sabe que la energia siempre es una función escalar construida a partir de la dinámica del sistema

 

Las expresiones para definir la energía son infinitas, dependiendo del tipo de aplicación con la cual se desee trabajar (energia potencial cinetica,…)

 

Básicamente se asocia al sistema una función escalar V(||x||) de tal forma que capture la energía asociada a la trayectoria del sistema y voy a analizar si dicha función va para cero en el transcurrir del tiempo.

 

Esta función escalar asocia la norma de los estados del sistema.

 

Si para nuestro analisis comenzamos recordando cual es la energia de una señal en el tiempo:

 

E_x=\sqrt{\int_{0}^{\infty}x(t)^2dt}

E_x^2=\int_{0}^{\infty}x(t)^2dt

Ahí se analiza si esa función va para cero a lo largo del tiempo. Una característica de esa definición de energía es que ella siempre es positiva

Energia Lyapunov

Para la integral poder ir para cero, se analiza que la característica de la gráfica es siempre positiva y que la derivada es siempre negativa porque ella siempre decrece. Y si se cumplen esas dos características, ese sistema siempre va ir para cero.

 

Entonces se verifica que un sistema es internamente estable si la energía es

E_x>0

\dfrac{dE_x}{dt}<0

Ahora si se tiene más de un estado, osea que no es solo un escalar, se puede usar la definición de producto interno

E_x^2=\int_{0}^{\infty}x(t)'x(t)dt

El único problema es que si por ejemplo mi función de energía no tiene una derivada que sea siempre negativa

Energia de Lyapunov

Dicha función es estable porque va para cero, pero no cumple con las dos propiedades, entonces yo no consigo capturar la dinámica de esa energía. El problema de esa abordaje, es que la función que se vaya a utilizar debe ser monotonicamente decreciente. Para conseguir eso se coloca una matriz P en los estados y se intenta hallar la matriz tal que la función sea monotonicamente decreciente

E_x^2=\int_{0}^{\infty}x(t)'Px(t)dt

En conclusión, para determinar la estabilidad por Lyapunov, consiste en poder determinar una función que presente un comportamiento de energía monotono para decir que el sistema es estable, de lo contrario si no se consigue determinar esta función, en términos generales no quiere decir que el sistema sea inestable.

Propiedades de Matrices Simetricas

Para entender como analizar la estabilidad por Lyapunov, es necesario rever algunas propiedades de las matrices simétricas:

 

Si M es una matriz cuadrada (m\times m)

 

Si M=M^T se dice que es una matriz Simetrica.

Si M=-M^T se dice que es una matriz Antisimetrica.

1. Cualquier matriz cuadrada equivale a la suma de una matriz simetrica y una matriz antisimetrica

M=\dfrac{M+M^T}{2}+\dfrac{M-M^T}{2}

2. La función cuadrática x^TMx asociada a una matriz antisimetrica es siempre CERO.

x^TMx=-x^TM^Tx=(-x^TM^Tx)^T=-x^TMx=0

Luego, la función x^TMx con M simétrica o no, es igual a una función cuadrática con matriz simétrica.

Matrices definidas positivas:

Una matriz cuadrada (m\times m) es definida positiva si

 

\forall x \neq 0 \rightarrow x^TMx>0

 

Lo que quiere decir que la Matriz M es definida positiva si la funcion x^TMx es definida positiva.

 

Aplicación a sistemas Lineales

Si tenemos un sistema lineal \dot{x}=Ax y probamos con la función de Lyapunov V(x)=x^TMx>0

 

Primero derivamos la función de Lyapunov

 

\dot{V}(x)=\dot{x}^TMx+x^TM\dot{x}=x^TA^TMx+x^TMAx

\dot{V}(x)=x^T (A^TM+MA)x<0

 

La derivada tiene que ser menor que cero

Como (A^TM+MA) es una forma cuadrática si yo tomo ese resultado que va a ser una matriz que es Simetrica basta con analizar los autovalores de esa matriz, si todos los autovalores son negativos, la matriz es definida negativa lo que indica que la derivada de Lyapunov siempre dará negativa.

 

Por otro lado si en la ecuación cuadratica (x^TMx>0) si M tiene todos los autovalores positivos, esa ecuación siempre va a dar positivo garantizando la estabilidad

 

Ahora esto solo vale si yo conozco la matriz M,si yo no conozco dicha matriz tendria que encontrar una Matriz M simetrica con todos los autovalores positivos de forma que (A^TM+MA) tenga todos los autovalores negativos

 

\exists M=M^T>0 : A^TM+MA<0

 

Para poder resolver ese problema analíticamente, en lugar de trabajar con esa desigualdad, se puede transformar esa expresión en una ecuación lineal.

Si tengo esta ecuación:

 

\exists M=M^T>0, N=N^T>0:A^TM+MA=-N

 

Si esa matriz N tiene todos los autovalores positivos y colocamos un signo (-) alfrente de N, quiere decir que esa matriz es Negativa. De esa forma consigo demostrar que el termino (A^TM+MA) es menor que cero resolviendo una simple ecuación. Esto es conocido como la ECUACIÓN DE LYAPUNOV.

 

De esa forma consigo por medio de un sistema de ecuaciónes encontrar una matriz M. Solo que cuando hallamos esa M, ademas de ser una matriz Simetrica ella debe tener todos los autovalores Positivos.

 

Teorema 1

Si todos los autovalores de A tienen parte real negativa si y solo si para una matriz dada que sea simetrica positiva N a la ecuación de Lyapunov A^TM+MA=-N se consigue una unica solución simetrica y definida positiva M. Indicando que el sistema es Estable caso contrario es Inestable (esto es solo valido para sistemas LTI), para No lineales o variantes en el tiempo, simplemente no se puede concluir.

Teorema 2

Si todos los autovalores de A tienen parte real negativa entonces la ecuación de Lyapunov A^TM+MA=-N tiene una unica solución para una dada matriz N simetrica y definida positiva dada por:

M=\int_{0}^{\infty}e^{A^Tt}Ne^{At}dt

Se define a P=A^TQ+QA como la ecuación de Lyapunov. Y la estabilidad asintotica del sistema está asegurada si siendo P una matriz simetrica definida negativa, se puede encontrar Q como una matriz simetrica definida positiva (Solución de la ecuación de Lyapunov)

Si los autovalores de A cumplen lo siguiente:

– Q es definida positiva si Re(\lambda_i)>0 \forall \lambda_i de A

– Q es definida negativa si Re(\lambda_i)<0 \forall \lambda_i de A

– Q es indefinida si \exists Re(\lambda_i)>0 y \exists Re(\lambda_i)<0 de A

Si te está sirviendo de algo la información de la pagina, podrías invitarme a un café y ayudarme a seguir manteniendo en pie el sitio WEB. Solo cuesta $2USD y me ayudarías enormemente a seguir con mi trabajo. Muchas Gracias por tu visita.




EJEMPLO

 

\dot{x}(t)=\begin{bmatrix}  0 & 1\\  -4 & -5  \end{bmatrix}x(t)

 

Sabemos que el sistema es asintoticamente estable si

 

V(||x||)>0 y [/latex]\dfrac{dV(||x||)}{dt}<0[/latex]

Tomamos como función de Lyapunov la expresión cuadratica:

V(x)=x^TMx>0 y [/latex]\dot{V}(x)=x^T (A^TM+MA)x<0[/latex]

Determinamos los Autovalores de A.

det(\lambda I-A)=0 \rightarrow (\lambda_1=-1, \lambda_2=-4)

Sea N=I y M=\begin{bmatrix}  a & b\\  b & c  \end{bmatrix}

A^TM+MA=-N

\begin{bmatrix}0 & 1\\-4 & -5\end{bmatrix}^T\begin{bmatrix}a & b\\b & c\end{bmatrix}+\begin{bmatrix}a & b\\b & c\end{bmatrix}\begin{bmatrix}0 & 1\\-4 & -5\end{bmatrix}=\begin{bmatrix}-1 & 0\\0 & -1\end{bmatrix}

\begin{bmatrix}0 & -4\\1 & -5\end{bmatrix}\begin{bmatrix}a & b\\b & c\end{bmatrix}+\begin{bmatrix}a & b\\b & c\end{bmatrix}\begin{bmatrix}0 & 1\\-4 & -5\end{bmatrix}=\begin{bmatrix}-1 & 0\\0 & -1\end{bmatrix}

\begin{bmatrix}-4b & -4c\\a-5b & b-5c\end{bmatrix}+\begin{bmatrix}-4b & a-5b\\-4c & b-5c\end{bmatrix}=\begin{bmatrix}-1 & 0\\0 & -1\end{bmatrix}

\begin{bmatrix}-8b & a-5b-4c\\a-5b-4c & 2b-10c\end{bmatrix}=\begin{bmatrix}-1 & 0\\0 & -1\end{bmatrix}

Ecuaciones:

-8b=-1 \rightarrow b=1/8=0.125

2b-10c=-1 \rightarrow c=\dfrac{2(1/8)+1}{10}=0.125

a-5b-4c=0 \rightarrow a=5(1/8)+4(1/8)=9/8=1.125

M=\begin{bmatrix}  1.125 & 0.125\\  0.125 & 0.125  \end{bmatrix}

V(x)=x^TMx>0\rightarrow \begin{bmatrix}x_1^T\\x_2^T\end{bmatrix}\begin{bmatrix}1.125 & 0.125\\0.125 & 0.125\end{bmatrix}\begin{bmatrix}x_1\\x_2\end{bmatrix}=1.125x_1^2+0.250x_1x_2+0.125x_2^2>0

Calculando los autovalores de la matriz M tenemos que:

det(\lambda I-M)=0 \rightarrow (\lambda_1=0.1096, \lambda_2=1.1404)

Como los autovalores son positivos, quiere decir que la Matriz M es simetrica y definida positiva, por lo tanto el sistema es estable.

En MATLAB:

> M=lyap(A,N)

 

Para Sistemas Discretos

Para el caso discreto la función de la energia de una señal viene dada por:

E_x^2=\sum_{k=0}^{\infty}x(k)^TPx(k)

Energia de Lyapunov

Para verificar si la función está decreciendo, siempre la diferencia entre un instante futuro con relación a un instante actual debe dar negativo.

V(x[k])>0 y [/latex]\Delta V(x[k])=V(x[k+1])-V(x[k])<0[/latex]

La ecuación de Lyapunov para el caso discreto viene dado por:

Si se tiene un sistema discreto

x[k+1]=Ax[k]

Y se utiliza la función de Lyapunov mas simple, conocida como función de Lyapunov Cuadrada:

V(x[k])=x[k]^TMx[k]>0

Y ya sabemos que como es una ecuación cuadrática, la matriz M tiene que tener todos los autovalores positivos y también ser simétrica.

Se calcula entonces la variación de la función:

\Delta V(x[k])=V(x[k+1])-V(x[k])<0

\Delta V(x[k])=x[k+1]^TMx[k+1]-x[k]^TMx[k]<0

\Delta V(x[k])=(Ax[k])^TM(Ax[k])-x[k]^TMx[k]<0

\Delta V(x[k])=x[k]^TA^TMAx[k]-x[k]^TMx[k]<0

\Delta V(x[k])=x[k]^T(A^TMA-M)x[k]<0

Y esta también es una forma cuadratica.

Para determinar si es estable basta con analizar si todos los autovalores de la matriz (A^TMA-M<0) son menosres que cero.

Para resolver analíticamente, se iguala la matriz a otra matriz negativa N, igual que en el caso continuo

A^TMA-M=-N

Si consigo determinar una matriz M definida positiva puedo decir que mi sistema dinámico Lineal e invariante en el tiempo es estable, recordando que si no es definida positiva seria inestable.

Para el caso de Lineal y Variante en el tiempo y No lineal en el caso de la matriz M ser definida negativa, no se puede determinar nada, y se tendria que seguir buscando otra función de Lyapunov.

 

Si te sirvió de algo la información de la pagina, podrías invitarme a un café y ayudarme a seguir manteniendo en pie el sitio WEB. Solo cuesta $2USD y me ayudarías enormemente a seguir con mi trabajo. Muchas Gracias por tu visita.




 Leave a Reply

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code class="" title="" data-url=""> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong> <pre class="" title="" data-url=""> <span class="" title="" data-url="">

(required)

(required)