5 votos

Neocercanía positivamente invariante mediante la función de Lyapunov

Dado el siguiente sistema de EDO no lineal,

$$x_1'=-x_1-x_2$$ $$x_2'=2x_1-x_2^3$$

Necesito usar la función cuadrática de Lyapunov

$$V(x) = x^TQx$$

où $Q$ es una matriz definida positiva tal que

$$A^TQ+QA=-I$$

y donde $A=Df(0,0)$ para encontrar un barrio $U$ alrededor del origen, lo más grande posible, tal que $U$ es positivamente invariante y todas las soluciones que empiezan en $U$ tienden al origen como $t\to\infty$ .

He examinado detenidamente los apuntes del curso y el libro de Perko ( Sistemas dinámicos y ecuaciones diferenciales ), pero no pude encontrar una pista de cómo hacer este problema.

Agradecería mucho que alguien me lo explicara:

  1. ¿Cómo se encuentra este barrio?

  2. ¿Cómo se $V$ relacionado con todo esto?

1 votos

Seguro que no puede encontrar $(a,b,c)$ tal que $V(x)=ax_1^2+bx_2^2+cx_1x_2$ es positiva definida y $\dot V(x)<0$ en un barrio de $(0,0)$ menos el punto $(0,0)$ y $\dot V(0,0)=0$ ? Normalmente se empieza por comprobar si algunos $V(x)=ax_1^2+bx_2^2$ no haría el trabajo...

1voto

CTNT Puntos 1718

Escriba a $$\dot{x}=Ax+g(x)=\left[\matrix{-1 & -1\\ 2 & 0}\right]x+\left[\matrix{0\\-x_2^3}\right]$$ La solución simétrica definida positiva de $$QA+A^TQ=-I$$ es (por cálculo directo) $$Q=\frac{1}{2}\left[\matrix{3 & 1\\ 1 & 2}\right]$$ Ahora para la función de Lyapunov candidata $V=x^TQx$ tenemos $$\dot{V}=\dot{x}^TQx+x^TQ\dot{x}=(Ax+g(x))^TQx+x^TQ(Ax+g(x))\\ =x^T(QA+A^TQ)x+2x^TQg(x)\\ =-\|x\|^2+2x^TQg(x)\\ \leq -\|x\|^2+2\lambda_{\max}(Q)\|x\|\|g(x)\|$$ Tenga en cuenta que $\|g(x)\|=|x_2|^3\leq \|x\|^3$ y por lo tanto $$\dot{V}\leq -\|x\|^2+2\lambda_{\max}(Q)\|x\|^4\\ \leq -\|x\|^2+2\frac{\lambda_{\max}(Q)}{\lambda_{\min}(Q)}V\|x\|^2\\ \leq -\|x\|^2\left(1-2\frac{\lambda_{\max}(Q)}{\lambda_{\min}(Q)}V\right)$$ Así $\dot{V}< 0$ para todos $$V< \frac{\lambda_{\min}(Q)}{2\lambda_{\max}(Q)}$$ con $V\neq 0$ y por lo tanto todas las trayectorias que comienzan dentro de la vecindad $$U:=\left\{x: V(x)< \frac{\lambda_{\min}(Q)}{2\lambda_{\max}(Q)}\right\}$$ permanecer en $U$ para todo el tiempo posterior y tienden al origen como $t\rightarrow\infty$ .

0 votos

¿Dónde está el $g(x) = \begin{bmatrix} 0\\ -x_2^3 \end{bmatrix}$ ¿De dónde viene? Además, ¿por qué $\dot{V}=x^T(QA+A^TQ)x+2x^TQg(x)$ ? Gracias

0 votos

Aún no he entendido del todo tu planteamiento, pero ¿podrías echarle un vistazo a lo que he hecho, tal vez esta forma también sea aceptable? $V(x) = 3/2 x_1^2 + x_1x_2 + x_2^2 > (1/2 x_1+x_2)^2 \ge 0$ . $\dot{V} = 3x_1 x_1' + x_1'x_2 + x_1x_2' +2x_2x_2' =-x_1^2-x_2^2-x_1x_2^3-2x_2^4 \le -\lvert x \lvert^2 +3\lvert x\lvert^4= \lvert x\lvert^2 (-1 + 3\lvert x\lvert^2) \implies \lvert x \lvert \le 1/\sqrt{3}.$ [Dónde $\lvert \cdot \lvert$ denota la norma vectorial]. @CTNT

1 votos

@secuencia El razonamiento anterior no es del todo correcto. Tenga en cuenta que $V\leq c_1$ es el área de una elipse, mientras que $|x| \leq c_2$ es un disco. Lo que has probado es que $V$ es decreciente para todos los puntos de algún disco. Así, la región positivamente invariante $U$ es la mayor superficie definida por $V \leq c_1$ e incluido dentro del disco no el disco en sí .

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X