Ecuaciones Diferenciales – Apuntes de clase

Lic. Adriana Valverde Calderón

SISTEMAS NO LINEALES

¿CÓMO PUEDE APROXIMARSE UN SISTEMA NO LINEAL CERCA DE UN PUNTO

DE EQUILIBRIO POR MEDIO DE UN SISTEMA LINEAL?
Análisis del punto de equilibrio

1. Ecuación de Van Der Pol:

𝑑2𝑥 𝑑𝑥
2
− (1 − 𝑥 2 ) +𝑥 =0
𝑑𝑡 𝑑𝑡
Expresar la ecuación por medio de un sistema de primer orden:
𝑑𝑥
=𝑦
{ 𝑑𝑡
𝑑𝑦
= (1 − 𝑥 2 )𝑦 − 𝑥
𝑑𝑡
El sistema no lineal tiene como único punto de equilibrio el punto (0 , 0)
¿Cómo es el término no lineal cerca del punto de equilibrio?
𝑥 𝑦 𝑥2𝑦
0,1 0,1 10−3
0,01 0,01 10−6

Se observa que el término no lineal 𝑥 2 𝑦 tiende más rápidamente a cero. Lo que

significa que estudiar el sistema no lineal cerca de cero será lo mismo que estudiar el
sistema lineal:
𝑑𝑥
=𝑦 1 ± √3𝑖
{ 𝑑𝑡 con valores propios
𝑑𝑦 2
=𝑦−𝑥
𝑑𝑡

2. Ejemplo no polinomial
𝑑𝑥
=𝑦
{ 𝑑𝑡 ; tiene infinitos puntos de equilibrio( (0,0); (π, 0); (−π, 0) … )
𝑑𝑦
= −𝑦 − 𝑠𝑒𝑛𝑥
𝑑𝑡

1
Ecuaciones Diferenciales – Apuntes de clase
Lic. Adriana Valverde Calderón

Usando el desarrollo en serie de potencias de la función 𝑠𝑒𝑛𝑥

𝑥 3 𝑥 5 𝑥 7 𝑥 9 𝑥11
𝑠𝑒𝑛𝑥 = 𝑥 − + − + − +⋯
3! 5! 7! 9! 11!
podemos escribir:
𝑑𝑥
=𝑦
𝑑𝑡
𝑑𝑦 𝑥3 𝑥5 𝑥7
= −𝑦 − (𝑥 − + − + ⋯ )
{ 𝑑𝑡 3! 5! 7!

cerca del origen cancelamos los términos no lineales

𝑑𝑥
=𝑦
{ 𝑑𝑡 ; tiene un solo punto de equilibrio
𝑑𝑦
= −𝑦 − 𝑥
𝑑𝑡
Con valores propios de la matriz de coeficientes:
−1 ± √3𝑖
𝜆=
2
Plano fase de:

𝑦 𝑦

𝑥 𝑥
0 0
−𝜋 𝜋

𝑑𝑥 𝑑𝑥
=𝑦 =𝑦
{ 𝑑𝑡 { 𝑑𝑡
𝑑𝑦 𝑑𝑦
= −𝑦 − 𝑠𝑒𝑛𝑥 = −𝑦 − 𝑥
𝑑𝑡 𝑑𝑡

3. Un modelo de especies en competencia

Sean 𝑥 e 𝑦 las poblaciones de dos especies que compiten por recurso. Un incremento
en cualquier especie tienen un efecto adverso sobre la razón de crecimiento de la otra.

2
Ecuaciones Diferenciales – Apuntes de clase
Lic. Adriana Valverde Calderón

𝑑𝑥 𝑥
= 2𝑥 (1 − ) − 𝑥𝑦
{ 𝑑𝑡 2 ; ¿ Cuántos puntos de equilibrio tiene?
𝑑𝑦 𝑦
= 3𝑦 (1 − ) − 2𝑥𝑦
𝑑𝑡 3

de las ecuaciones: 𝑥(2 − 𝑥 − 𝑦) = 0 ; 𝑦(3 − 𝑦 − 2𝑥) = 0 se obtienen cuatro puntos

de equilibrio: (0
⏟ , 0) ; (0,
⏟ 3) ; (2 (1 ,1)
⏟ , 0) 𝑦 ⏟
(𝑥,𝑦) (𝑥,𝑦) (𝑥,𝑦) (𝑥,𝑦)

¿Cuál es el significado de cada punto?

Es de particular interés el punto (1 ,1) que indica que es posible para estas dos especies
coexistir en equilibrio.
¿Será posible que el análisis cualitativo de los sistemas lineales sirva para analizar el
comportamiento de un sistema no lineal cerca de los puntos de equilibrio (P.E)?
¿Cómo se comporta el sistema cerca de (1 ,1)?
Aproximar: Sistema no lineal → Sistema lineal
P.E (1,1) → (0 , 0) siempre de P.E
¿Cómo realizar esta aproximación del sistema NO LINEAL al sistema LINEAL?
1º- Mover el punto (1 ,1) al (0 , 0) por medio de la introducción de dos variables
𝑢 =𝑥−1
{
𝑣 =𝑦−1
Observe que (𝑢 , 𝑣) → (0 , 0) cuando (𝑥 , 𝑦) → (1 , 1)
2º- Obtener el sistema de las nuevas variables:
Calculamos:
𝑑𝑢 𝑑(𝑥 − 1) 𝑑𝑥 𝑑𝑣 𝑑(𝑦 − 1) 𝑑𝑦
= = ; = =
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
Sustituyendo en el sistema dado:
𝑑𝑢 𝑢+1
= 2(𝑢 + 1) (1 − ) − (𝑢 + 1)(𝑣 + 1)
𝑑𝑡 2
𝑑𝑣 (𝑣 + 1)
= 3(𝑣 + 1) (1 − ) − 2(𝑢 + 1)(𝑣 + 1)
{ 𝑑𝑡 3

En términos de las nuevas variables tenemos el sistema:

3
Ecuaciones Diferenciales – Apuntes de clase
Lic. Adriana Valverde Calderón

𝑑𝑢
= −𝑢 − 𝑣 − 𝑢2 − 𝑢𝑣
{ 𝑑𝑡
𝑑𝑣
= −2𝑢 − 𝑣 − 2𝑢𝑣 − 𝑣 2
𝑑𝑡
Observe que el punto (0 , 0) es Punto de Equilibrio para este último sistema no lineal,
pero cerca de (0 , 0) los términos no lineales tienden más rápidamente a cero.
3º- Aproximar el sistema no lineal cerca de (0 , 0) a un sistema lineal:
𝑑𝑢
= −𝑢 − 𝑣
{ 𝑑𝑡
𝑑𝑣
= −2𝑢 − 𝑣
𝑑𝑡
Para este sistema lineal, los valores propios son: 𝜆 = −1 ± √2 , por lo que el punto de
equilibrio (0 , 0) se comporta como un punto de silla.

𝑦
Planos Fase
𝑣
3

2
𝑢
0
1

𝑥
0 2
1

𝑑𝑥 𝑥 𝑑𝑢
= 2𝑥 (1 − ) − 𝑥𝑦 = −𝑢 − 𝑣
{ 𝑑𝑡 2 { 𝑑𝑡
𝑑𝑦 𝑦 𝑑𝑣
= 3𝑦 (1 − ) − 2𝑥𝑦 = −2𝑢 − 𝑣
𝑑𝑡 3 𝑑𝑡

¿Cuál es el proceso analítico para aproximar un sistema no lineal cualesquiera a un sistema

lineal?

4
Ecuaciones Diferenciales – Apuntes de clase
Lic. Adriana Valverde Calderón

PROCESO DE LINEALIZACIÓN

Se desea usar un sistema lineal para predecir el comportamiento de un sistema no lineal

cerca de un punto de equilibrio.
Considere el sistema no lineal:
𝑑𝑥
= 𝑓(𝑥, 𝑦)
{ 𝑑𝑡 que tiene a (𝑥0 , 𝑦0 ) como Punto de equilibrio
𝑑𝑦
= 𝑔(𝑥, 𝑦)
𝑑𝑡
Para linealizar el sistema cerca del punto de equilibrio, introducimos nuevas variables:
𝑢 = 𝑥 − 𝑥0
{𝑣 = 𝑦 − 𝑦
0

Que mueven al punto de equilibrio (𝑥0 , 𝑦0 ) en XY al origen (0 , 0) en UV, de manera que:

Cuando (𝑥 , 𝑦) → (𝑥0 , 𝑦0 ) ; (𝑢 , 𝑣) → (0 , 0).
El sistema en términos de 𝑢 y 𝑣 será:

𝑑𝑢 𝑑𝑥 𝑑𝑢
= = 𝑓(𝑥, 𝑦) = 𝑓(𝑥0 + 𝑢 , 𝑦0 + 𝑣) = 𝑓(𝑥0 + 𝑢 , 𝑦0 + 𝑣)
{ 𝑑𝑡 𝑑𝑡 → { 𝑑𝑡
𝑑𝑣 𝑑𝑦 𝑑𝑣
= = 𝑔(𝑥, 𝑦) = 𝑔(𝑥0 + 𝑢 , 𝑦0 + 𝑣) = 𝑔(𝑥0 + 𝑢 , 𝑦0 + 𝑣)
𝑑𝑡 𝑑𝑡 𝑑𝑡
Con único punto de equilibrio: (𝑢 , 𝑣) = (0 , 0) , verificándose:
0 = 𝑓(𝑥0 , 𝑦0 )
0 = 𝑔(𝑥0 , 𝑦0 )
Lo que significa que el punto de equilibrio (𝑥0 , 𝑦0 ) en XY, se ubica en (0 , 0) en UV.
Analizando las funciones 𝑓 y 𝑔 : la mejor aproximación lineal de una función de dos
variables está dada por el plano tangente.
La ecuación del plano tangente en el punto (𝑥, 𝑦, 𝑧) es:
𝑑𝑧 = ∇𝑓(𝑥, 𝑦)(𝑑𝑥, 𝑑𝑦)
En el punto (𝑥0 , 𝑦0 , 𝑧0 ) se tendrá:
𝜕𝑓(𝑥0 , 𝑦0 ) 𝜕𝑓(𝑥0 , 𝑦0 )
𝑧 − 𝑧0 = ( , ) (𝑥 − 𝑥0 , 𝑦 − 𝑦0 )
𝜕𝑥 𝜕𝑦
∴ 𝐿a mejor aproximación lineal de 𝑓 en el punto (𝑥0 , 𝑦0 , 𝑧0 ) será:
𝜕𝑓(𝑥0 , 𝑦0 ) 𝜕𝑓(𝑥0 , 𝑦0 )
⏟ 0 , 𝑦0 ) +
𝑓(𝑥, 𝑦) = 𝑓(𝑥 (𝑥 − 𝑥0 ) + (𝑦 − 𝑦0 )
𝜕𝑥 𝜕𝑦
0

5
Ecuaciones Diferenciales – Apuntes de clase
Lic. Adriana Valverde Calderón

Considerando:
𝜕𝑓(𝑥0 , 𝑦0 ) 𝜕𝑓(𝑥0 , 𝑦0 )
𝑓(𝑥, 𝑦)
⏟ 𝑓(𝑥0 , 𝑦0 ) +
=⏟ (𝑥 − 𝑥0 ) +
⏟ (𝑦 − 𝑦0 )
⏟
𝜕𝑥 𝜕𝑦
𝑓(𝑥0 +𝑢 , 𝑦0 +𝑣) 0 𝑢 𝑣

Reescribiendo el sistema en términos de UV:

𝑑𝑢 𝜕𝑓(𝑥0 , 𝑦0 ) 𝜕𝑓(𝑥0 , 𝑦0 )
= 𝑢+ 𝑣
𝑑𝑡 𝜕𝑥 𝜕𝑦
En forma análoga obtenemos la aproximación para la función 𝑔(𝑥, 𝑦):
𝑑𝑣 𝜕𝑔(𝑥0 , 𝑦0 ) 𝜕𝑔(𝑥0 , 𝑦0 )
= 𝑢+ 𝑣
𝑑𝑡 𝜕𝑥 𝜕𝑦
El sistema linealizado en las variables 𝑢, 𝑣:

𝑑𝑢 𝜕𝑓(𝑥0 , 𝑦0 ) 𝜕𝑓(𝑥0 , 𝑦0 )
𝜕𝑥 𝜕𝑦 𝑢
( 𝑑𝑡 ) = ( )
𝑑𝑣 𝜕𝑔(𝑥0 , 𝑦0 ) 𝜕𝑔(𝑥0 , 𝑦0 ) 𝑣
𝑑𝑡 ⏟
( 𝜕𝑥 𝜕𝑦 )
𝑀𝑎𝑡𝑟𝑖𝑧 𝐽𝑎𝑐𝑜𝑏𝑖𝑎𝑛𝑎

Ejemplo:
𝑑𝑥
= −2𝑥 + 2𝑥 2
{ 𝑑𝑡 con Pts. de E: (0,0) 𝑦 (1,0)
𝑑𝑦
= −3𝑥 + 𝑦 + 3𝑥 2
𝑑𝑡
Calculamos la matriz Jacobiana
𝜕𝑓(𝑥 , 𝑦) 𝜕𝑓(𝑥 , 𝑦)
𝜕𝑥 𝜕𝑦 −2 + 4𝑥 0
=( )
𝜕𝑔(𝑥 , 𝑦) 𝜕𝑔(𝑥 , 𝑦) −3 + 6𝑥 1
( 𝜕𝑥 𝜕𝑦 )
Cerca de (0,0) el plano fase para el sistema no lineal debe parecerse al plano fase del
sistema lineal:
𝑑𝑌 −2 0
=( )𝑌
𝑑𝑡 −3 1
1 0
Valores propios: 𝜆1 = −2 ; 𝜆2 = 1 y vectores propios 𝑉𝜆1 = ( ) ; 𝑉𝜆2 = ( )
1 1
Cerca de (1,0) el plano fase para el sistema no lineal debe parecerse al plano fase del
sistema lineal:
𝑑𝑌 2 0
=( )𝑌
𝑑𝑡 3 1

6
Ecuaciones Diferenciales – Apuntes de clase
Lic. Adriana Valverde Calderón

1 0
Valores propios: 𝜆1 = 2 ; 𝜆2 = 1 y vectores propios 𝑉𝜆1 = ( ) ; 𝑉𝜆2 = ( )
3 1

Cerca de (0,0) Cerca de (1,0)

𝑣 𝑣

𝑢 𝑢
𝑦

Punto de silla Fuente

1
𝑥

Separatriz: Separa comportamientos