Asignatura: Ecuaciones diferenciales II.

Carrera: Licenciatura en matemáticas

Alumno: Raúl Ibáñez Couoh

Matrícula: ES172001745

Docente: Rafael Pacheco Espinosa

Unidad 2. Sistemas de ecuaciones no lineales

Evidencia de aprendizaje. Solución de sistemas autónomos y

análisis de estabilidad

20/02/2020, Zihuatanejo, Guerrero, México.

 dx ⎡2 5⎤
―＝ 2 x + 5 y A＝⎢
dt ⎣ 1 6 ⎥⎦
dy
―＝ x + 6 y
dt
|⎡ 2 5 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ 1 6 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ - 8 ⋅ λ + 7
2

solve , λ ⎡ 7 ⎤
λ 2 - 8 ⋅ λ + 7 ＝ 0 ――― →⎢ ⎥
⎣1⎦
⎛⎡ 2 5 ⎤⎞ ⎡ 1 ⎤
eigenvals ⎜⎢ →
⎝⎣ 1 6 ⎥⎦⎟⎠ ⎢⎣ 7 ⎥⎦

⎛⎡ 2 5 ⎤⎞ ⎡ 1 -5 ⎤ ⎡ 1 -5 ⎤
eigenvecs ⎜⎢ →
⎝⎣ 1 6 ⎥⎦⎟⎠ ⎢⎣ 1 1 ⎥⎦ ⎢⎣ 1 1 ⎥⎦

⎡1⎤ ⎡ -5 ⎤ 7 t
XC ＝ c1 ⎢ ⎥ e t + c2 ⋅ ⎢ ⋅e
⎣1⎦ ⎣ 1 ⎥⎦

x ((t)) ＝ c1 ⋅ e t - 5 c2 ⋅ e 7 t

y ((t)) ＝ c1 ⋅ e t + c2 ⋅ e 7 t
⎛⎡ 2 5 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ 1 6 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦

Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦
 Mediante pplane8:

Tenemos una fuente nodal.

dx
―＝ 2 x - y ⎡ 2 -4 ⎤
dt A＝⎢
⎣ -6 1 ⎥⎦
dy
―＝ -6 x + y
dt
|⎡ 2 -4 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ -6 1 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ - 3 ⋅ λ - 22
2

⎡ ‾‾ ⎤
⎢ ―――97 + 3 ⎥
solve , λ ⎢ 2 ⎥
λ 2 - 3 ⋅ λ - 22 ＝ 0 ――― →
⎢ ⎥
- ‾‾
⎢ ――――97 + 3 ⎥
⎢⎣ 2 ⎥⎦
⎡ ⎤
- ‾‾
⎢ ―――― 97 + 3 ⎥ ⎡ ⎤
- ‾‾
97 - 1 ‾‾
97 - 1
⎛⎡ 2 -4 ⎤⎞ ⎢ 2 ⎥ ⎛ ⎡ 2 -4 ⎤ ⎞ ⎢ ―――⎥
eigenvals ⎜⎢ ⎥ ⎟ → eigenvecs ⎜⎢ ⎥ ⎟ → ⎢ ――――
12 12 ⎥
⎝⎣ -6 1 ⎦⎠ ⎢ ‾‾ ⎥ ⎝⎣ -6 1 ⎦⎠
⎢ ―――⎥97 + 3 ⎢⎣ ⎥⎦
1 1
⎢⎣ 2 ⎥⎦

⎛⎡ 2 -4 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ -6 1 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto Silla.

dx
―＝ 5 x + 6 y ⎡5 6⎤
dt A＝⎢
⎣ 1 4 ⎥⎦
dy
―＝ x + 4 ⋅ y
dt
|⎡ 5 6 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ 1 4 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ - 9 ⋅ λ + 14
2

solve , λ ⎡ 7 ⎤
λ 2 - 9 ⋅ λ + 14 ＝ 0 ――― →⎢ ⎥
⎣2⎦
⎛⎡ 5 6 ⎤⎞ ⎡ 2 ⎤
eigenvals ⎜⎢ →
⎝⎣ 1 4 ⎥⎦⎟⎠ ⎢⎣ 7 ⎥⎦

⎛⎡ 5 6 ⎤⎞ ⎡ 3 -2 ⎤
eigenvecs ⎜⎢ →
⎝⎣ 1 4 ⎥⎦⎟⎠ ⎢⎣ 1 1 ⎥⎦

⎛⎡ 5 6 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ ⎥⎟ → ⎢ ⎥
 ⎛⎡ 5 6 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ 1 4 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦

Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto espiral inestable, o repulsor fuente o fuente nodal.

dx
―＝ 5 x - 3 y ⎡ 5 -3 ⎤
dt A＝⎢
⎣ -1 -1 ⎥⎦
dy
―＝ -x - y
dt
|⎡ 5 -3 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ -1 -1 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ - 4 ⋅ λ - 8
2

solve , λ ⎡ 2 ⋅ ‾‾ 3 +2 ⎥
⎤
λ 2 - 4 ⋅ λ - 8 ＝ 0 ――― →⎢
⎢⎣ -⎛⎝2 ⋅ ‾‾ ⎞
3 ⎠ + 2 ⎥⎦

⎛⎡ 5 -3 ⎤⎞ ⎡ -⎛⎝2 ⋅ ‾‾ ⎞ ⎤
3⎠ + 2 ⎥
eigenvals ⎜⎢ → ⎢
⎝⎣ -1 -1 ⎥⎦⎟⎠ ⎢⎣ 2 ⋅ ‾‾
3 + 2 ⎥⎦
 ⎛⎡ 5 -3 ⎤⎞ ⎡ -⎛⎝2 ⋅ ‾‾ ⎞
3 ⎠ - 3 2 ⋅ ‾‾
⎤
3 -3⎥
eigenvecs ⎜⎢ ⎥ ⎟ →⎢
⎝⎣ -1 -1 ⎦⎠ ⎣ 1 1 ⎦
⎛⎡ 5 -3 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ -1 -1 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦

Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto silla

dx
―＝ 2 x + 9 y ⎡ 2 9 ⎤
dt A＝⎢
⎣ -1 -4 ⎥⎦
dy
―＝ -x - 4 y
dt
|⎡ 2 9 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ -1 -4 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ + 2 ⋅ λ + 1
2

solve , λ ⎡ -1 ⎤
λ 2 + 2 ⋅ λ + 1 ＝ 0 ――― →⎢
⎣ -1 ⎥⎦

⎛⎡ 2 9 ⎤⎞ ⎡ -1 ⎤
eigenvals ⎜⎢ ⎥⎟ → ⎢ ⎥
 solve , λ ⎡ -1 ⎤
λ 2 + 2 ⋅ λ + 1 ＝ 0 ――― →⎢
⎣ -1 ⎥⎦

⎛⎡ 2 9 ⎤⎞ ⎡ -1 ⎤
eigenvals ⎜⎢ →
⎝⎣ -1 -4 ⎥⎦⎟⎠ ⎢⎣ -1 ⎥⎦

⎛⎡ 2 9 ⎤⎞ ⎡ -3 ⎤
eigenvecs ⎜⎢ →
⎝⎣ -1 -4 ⎥⎦⎟⎠ ⎢⎣ 1 ⎥⎦
⎛⎡ 2 9 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ -1 -4 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦

Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto vórtice.
 dx
―＝ x + y ⎡ 1 1⎤
dt A＝⎢
⎣ -1 3 ⎥⎦
dy
―＝ -x + 3 y
dt
|⎡ 1 1 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ -1 3 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ - 4 ⋅ λ + 4
2

solve , λ ⎡ 2 ⎤
λ 2 - 4 ⋅ λ + 4 ＝ 0 ――― →⎢ ⎥
⎣2⎦
⎛⎡ 1 1 ⎤⎞ ⎡ 2 ⎤
eigenvals ⎜⎢ →
⎝⎣ -1 3 ⎥⎦⎟⎠ ⎢⎣ 2 ⎥⎦

⎛⎡ 1 1 ⎤⎞ ⎡ 1 ⎤
eigenvecs ⎜⎢ →
⎝⎣ -1 3 ⎥⎦⎟⎠ ⎢⎣ 1 ⎥⎦

⎛⎡ 1 1 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ -1 3 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto espiral inestable o repulsor fuente

 dx
―＝ 7 x - 4 y ⎡ 7 -4 ⎤
dt A＝⎢
⎣ 3 0 ⎥⎦
dy
―＝ 3 x
dt
|⎡ 7 -4 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ 3 0 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ - 7 ⋅ λ + 12
2

solve , λ ⎡ 4 ⎤
λ 2 - 7 ⋅ λ + 12 ＝ 0 ――― →⎢ ⎥
⎣3⎦
⎛⎡ 7 -4 ⎤⎞ ⎡ 3 ⎤
eigenvals ⎜⎢ →
⎝⎣ 3 0 ⎥⎦⎟⎠ ⎢⎣ 4 ⎥⎦
⎡4 ⎤
⎛⎡ 7 -4 ⎤⎞ ⎢ ―1 ⎥
eigenvecs ⎜⎢ → 3
⎝⎣ 3 0 ⎥⎦⎟⎠ ⎢ ⎥
⎣ 1 1⎦
⎛⎡ 7 -4 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ 3 0 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto espiral inestable o repulsor fuente

 dx
―＝ -x + y ⎡ -1 1 ⎤
dt A＝⎢
⎣ -6 -6 ⎥⎦
dy
―＝ -6 x - 6 y
dt
|⎡ -1 1 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ -6 -6 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ + 7 ⋅ λ + 12
2

solve , λ ⎡ -3 ⎤
λ 2 + 7 ⋅ λ + 12 ＝ 0 ――― →⎢
⎣ -4 ⎥⎦
⎛⎡ -1 1 ⎤⎞ ⎡ -4 ⎤
eigenvals ⎜⎢ →
⎝⎣ -6 -6 ⎥⎦⎟⎠ ⎢⎣ -3 ⎥⎦
⎡ 1 1⎤
⎛⎡ -1 1 ⎤⎞ ⎢ -―-― ⎥
eigenvecs ⎜⎢ ⎥ ⎟ → 2 3⎥
⎝⎣ -6 -6 ⎦⎠ ⎢
⎣ 1 1 ⎦
⎛⎡ -1 1 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ -6 -6 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Nodo estable.
 dx
―＝ -x + y ⎡ -1 1 ⎤
dt A＝⎢
⎣ 0 -1 ⎥⎦
dy
―＝ -y
dt
|⎡ -1 1 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ 0 -1 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ + 2 ⋅ λ + 1
2

solve , λ ⎡ -1 ⎤
λ 2 + 2 ⋅ λ + 1 ＝ 0 ――― →⎢
⎣ -1 ⎥⎦
⎛⎡ -1 1 ⎤⎞ ⎡ -1 ⎤
eigenvals ⎜⎢ →
⎝⎣ 0 -1 ⎥⎦⎟⎠ ⎢⎣ -1 ⎥⎦

⎛⎡ -1 1 ⎤⎞ ⎡ 1 ⎤
eigenvecs ⎜⎢ →
⎝⎣ 0 -1 ⎥⎦⎟⎠ ⎢⎣ 0 ⎥⎦

⎛⎡ -1 1 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ 0 -1 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto vórtice.
 dx
―＝ -x ⎡ -1 0 ⎤
dt A＝⎢
⎣ 0 -1 ⎥⎦
dy
―＝ -y
dt
|⎡ -1 0 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ 0 -1 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ + 2 ⋅ λ + 1
2

solve , λ ⎡ -1 ⎤
λ 2 + 2 ⋅ λ + 1 ＝ 0 ――― →⎢
⎣ -1 ⎥⎦
⎛⎡ -1 0 ⎤⎞ ⎡ -1 ⎤
eigenvals ⎜⎢ →
⎝⎣ 0 -1 ⎥⎦⎟⎠ ⎢⎣ -1 ⎥⎦

⎛⎡ -1 0 ⎤⎞ ⎡ 0 1 ⎤
eigenvecs ⎜⎢ →
⎝⎣ 0 -1 ⎥⎦⎟⎠ ⎢⎣ 1 0 ⎥⎦

⎛⎡ -1 0 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ 0 -1 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto Vórtice.
 dx
―＝ -3 x ⎡ -3 0 ⎤
dt A＝⎢
⎣ 0 2 ⎥⎦
dy
―＝ 2 y
dt
|⎡ -3 0 ⎤ ⎡ 1 0 ⎤|
||⎢⎣ 0 2 ⎥⎦ - λ ⋅ ⎢⎣ 0 1 ⎥⎦|| → λ + λ - 6
2

solve , λ ⎡ 2 ⎤
λ 2 + λ - 6 ＝ 0 ――― →⎢
⎣ -3 ⎥⎦
⎛⎡ -3 0 ⎤⎞ ⎡ -3 ⎤
eigenvals ⎜⎢ →
⎝⎣ 0 2 ⎥⎦⎟⎠ ⎢⎣ 2 ⎥⎦

⎛⎡ -3 0 ⎤⎞ ⎡ 0 1 ⎤
eigenvecs ⎜⎢ →
⎝⎣ 0 2 ⎥⎦⎟⎠ ⎢⎣ 1 0 ⎥⎦

⎛⎡ -3 0 0 ⎤⎞ ⎡ 1 0 0 ⎤
rref ⎜⎢ →
⎝⎣ 0 2 0 ⎥⎦⎟⎠ ⎢⎣ 0 1 0 ⎥⎦
Como x ＝ 0 y y ＝ 0 , tenemos que los puntos críticos del sistema autónomo son: [[ 0 0 ]]
⎡ x0 ⎤ ⎡ 0 ⎤
⎢ ⎥＝⎢ ⎥
⎣ y0 ⎦ ⎣ 0 ⎦

Punto silla.
Tenemos un Sumidero nodal o un nodo estable.
Sin punto de equilibrio específico.

𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽:
−2 0
� �
0 1
Eigenvalores e eigenvectores:
−2 1 0
� �� �
1 01
 𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽𝐽:
2 0
� �
0 −1
Eigenvalores e eigenvectores:
−1 0 1
� �� �
2 10
Puntos de equilibrio en 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠.
(−1,0)
(1,0)

a) Siendo 𝑎𝑎 = 𝑐𝑐 y𝑏𝑏 = 𝑑𝑑
Existe una línea de equilibrio donde cada línea se aproxima a un punto de equilibrio.
Se cumple 𝑎𝑎𝑑𝑑 − 𝑐𝑐𝑐𝑐 = 0.
Donde todo punto es un punto de equilibrio
 Conclusiones
El trabajo mostrado anteriormente, nos otorga la capacidad de analizar la estabilidad
de sistemas ordinarios de ecuaciones diferenciales, dónde es importante conocer
las matemáticas que rodean el análisis mismo.
Por medio de 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝8 elaboré los diagramas de fase de los distintos sistemas a
estudiar en esta Evidencia de aprendizaje, lo cual me permitió tener un mejor
entendimiento de la unidad.
Sin lugar a duda Matlab es una de las mejores herramientas para el estudio de las
matemáticas y sus métodos numéricos que emplean para la obtención de resultados
de sistemas de ecuaciones diferenciales lineales, es de gran ayuda.

Bibliografía
Ayres, F. (1952 ). Ecuaciones Diferenciales . México: Mc GRaw-Hill.
Blanchard, P. (1999). Ecuaciones Diferenciales . México : Thomson .
Braun, M. (1983). Ecuaciones Diferenciales y sus aplicaciones . México : Grupo
Editorial Iberoamérica.
Bronson, R. (2006). Ecuaciones Diferenciales . México : McGraw-Hill .
Carmona, I. (2011). Ecuaciones diferenciales. México: Pearson.
Carmona, I. (2011). Ecuaciones Diferenciales . México : Pearson .
David, A. (2005). Ecuaciones Diferenciales y problemas con valores en la frontera .
México : Pearson Educación .
Edwards, H. (2009). Ecuaciones Diferenciales y problemas con valores en la
frontera Cómputo y modelado . México : Pearson .
Espinoza, E. (1996). Ecuaciones Diferenciales Aplicaciones. Lima Peru : Lima.
Gabiria, J. A. (1997 ). Ecuaciones Diferenciales Ordinarias para Estudiantes de
Física. Madrid: Argitalpen Zerbitzua .
[Link]. (1969 ). Ecuaciones Diferenciales y cálculo variacional . Moscu : MIR .