Medi2 U3 Ea Jeru
Medi2 U3 Ea Jeru
Medi2 U3 Ea Jeru
(MEDI2)
Licenciatura en Matemáticas
Evidencia de aprendizaje
1. Utilice un análisis de linealización para trazar las órbitas fase de los siguientes sistemas.
Utilizando algún programa matemático grafique las soluciones numéricas y compare las
soluciones obtenidas.
𝑥′ = 𝑥 + 𝑦
a)
𝑦′ = 𝑥2 + 𝑦
Puntos críticos
𝑥 ′ = 𝑥 + 𝑦; 𝑥 = 𝑦 = 0, 𝑥 = −𝑦
𝑦 ′ = 𝑥 2 + 𝑦; 𝑦 ′ = (−𝑦)2 + 𝑦 = 𝑦 2 + 𝑦 = 𝑦(𝑦 + 1) = 0; 𝑦1 = 0, 𝑦2 = −1
Por lo tanto, (0,0) es punto de equilibrio
𝑥′ = 𝑥 + 𝑦 = 0 + 0 = 0
𝑦 ′ = 𝑥 2 + 𝑦 = 02 + 0 = 0
Por lo tanto, (1, −1) es punto de equilibrio
𝑥 ′ = 𝑥 + 𝑦 = 1 + (−1) = 0
𝑦 ′ = 𝑥 2 + 𝑦 = (1)2 + (−1) = 0
Cerca de los puntos de equilibrio
𝜕 𝜕
𝑥 ′ = 𝑥 + 𝑦 = 𝑓 (𝑥, 𝑦), 𝑥 ′ ≈ 𝑓 (𝑥, 𝑦)(𝑥 − 𝑥0 ) + 𝑓(𝑥, 𝑦)(𝑦 − 𝑦0 )
𝜕𝑥 𝜕𝑦
𝜕 𝜕
𝑦 ′ = 𝑥 2 + 𝑦 = 𝑔(𝑥, 𝑦), 𝑦 ′ ≈ 𝑔(𝑥, 𝑦)(𝑥 − 𝑥0 ) + 𝑔(𝑥, 𝑦)(𝑦 − 𝑦0 )
𝜕𝑥 𝜕𝑦
Por lo tanto:
𝜕𝑓 𝜕𝑓
′ 𝜕𝑦 𝑥 − 𝑥0
𝑥 𝜕𝑥 1 1 𝑥 − 𝑥0
( ′) ≈ (𝑦 − 𝑦 ) = ( )( )
𝑦 𝜕𝑔 𝜕𝑔 0 2𝑥 1 𝑦 − 𝑦0
(𝜕𝑥 𝜕𝑦 )
𝜆 = 1 ∓ √2(0) = 1
Las raíces son iguales y positivas por lo tanto es un nodo inestable
Jesús Abraham Rojas Úrzulo Licenciatura en Matemáticas
Matricula: ES1821013126 Ecuaciones Diferenciales II (MEDI2)
Por lo tanto, toda solución que comience fuera de los puntos críticos es inestable.
𝑥 ′ = 𝑦 2 − (𝑥 + 1) 2
b)
𝑦 ′ = −𝑥
Puntos críticos
𝑦 = 1+𝑥
𝑥 ′ = 𝑦 2 − (𝑥 + 1 )2 ; 𝑦 2 = (𝑥 + 1 )2 , |𝑦 | = |𝑥 + 1 |,
𝑦 = −1 − 𝑥
𝑦1 = 1
𝑦 ′ = −𝑥; 𝑥 = 0;
𝑦2 = −1
Puntos críticos (0, −1) y (0,1)
𝜕𝑓 𝜕𝑓
′ 𝜕𝑦 𝑥 − 𝑥0
𝑥 𝜕𝑥 −2(𝑥 + 1) 2𝑦 𝑥 − 𝑥0
( ′) ≈ (𝑦 − 𝑦 ) = ( )( )
𝑦 𝜕𝑔 𝜕𝑔 0 −1 0 𝑦 − 𝑦0
(𝜕𝑥 𝜕𝑦 )
𝜆2 + 𝜆(2(0) + 2) + 2(−1) = 𝜆2 + 2𝜆 − 2 = 0; 𝜆 = −1 ± √3
Es un punto silla, 𝜆1 < 0, 𝜆2 > 0
𝑥′ = 1 − 𝑦2
c)
𝑦′ = 𝑦 − 𝑥
Puntos críticos
𝑥 ′ = 1 − 𝑦 2 ; 𝑦1 = 1, 𝑦2 = −1
𝑦 ′ = 𝑦 − 𝑥; 𝑦 = 𝑥; 𝑥1 = 1, 𝑥2 = −1
(−1, −1) y (1,1)
1 √7
𝜆2 − 𝜆 − 2𝑦0 = 𝜆2 − 𝜆 − 2(−1) = 𝜆2 − 𝜆 + 2, 𝜆 = + 𝑖
2 2
Es un foco inestable
En el punto de equilibrio (1,1)
𝜆2 − 𝜆 − 2𝑦0 = 𝜆2 − 𝜆 − 2(1) = 𝜆2 − 𝜆 − 2, 𝜆1 = −1, 𝜆2 = 2
Es un punto silla
𝑥 ′ = 𝑥 (2 − 𝑦 )
d)
𝑦 ′ = 𝑦(2 − 𝑥)
Jesús Abraham Rojas Úrzulo Licenciatura en Matemáticas
Matricula: ES1821013126 Ecuaciones Diferenciales II (MEDI2)
Puntos críticos
𝑥 ′ = 𝑥(2 − 𝑦); 𝑥 = 0, 𝑦 = 2
𝑦 ′ = 𝑦(2 − 𝑥 ); 𝑦 = 0, 𝑥 = 2
(0,0) y (2,2)
La matriz Jacobina cerca de (𝑥0 , 𝑦0 )
(2 − 𝑦0 ) −𝑥0 𝑦0 (2 − 𝑦0 ) −𝑥0 𝑦0 1 0 (2 − 𝑦0 ) − 𝜆 −𝑥0 𝑦0
( ) ; det [( )−𝜆( )] = det [( )]
−𝑥0 𝑦0 (2 − 𝑥0 ) −𝑥0 𝑦0 (2 − 𝑥0 ) 0 1 −𝑥0 𝑦0 (2 − 𝑥0 ) − 𝜆
2 2
= [(2 − 𝑦0 ) − 𝜆][(2 − 𝑥0 ) − 𝜆] − 𝑥0 𝑦0
= 𝜆2 − 2𝜆((2 − 𝑦0 ) + (2 − 𝑥0 )) + (2 − 𝑥0 )(2 − 𝑦0 ) − 𝑥02 𝑦02
2. Analizar la estabilidad del punto crítico del siguiente sistema, proponiendo como función de
Lyapunov 𝑉 (𝑥, 𝑦) = 𝑎𝑥 2𝑚 + 𝑏𝑦 2𝑛
𝑥 ′ = −2𝑥𝑦
a)
𝑦′ = 𝑥2 − 𝑦3
𝜕𝑉 𝑑𝑥 𝜕𝑉 𝑑𝑦
𝑉 ′ (𝑥, 𝑦) = + = [2𝑚𝑎𝑥 2𝑚−1 ][−2𝑥𝑦] + [2𝑛𝑏𝑦 2𝑛−1 ][𝑥 2 − 𝑦 3 ]
𝜕𝑥 𝑑𝑡 𝜕𝑦 𝑑𝑡
−4 ∗ 𝑚 ∗ 𝑎 ∗ 𝑥 2𝑚 𝑦 + 2 ∗ 𝑛 ∗ 𝑏 ∗ 𝑥 2 𝑦 2𝑛−1 − 2 ∗ 𝑛 ∗ 𝑏 ∗ 𝑦 2𝑛+2
Los términos posibles a eliminar, son:
2 ∗ 𝑛 ∗ 𝑏 ∗ 𝑥 2 𝑦 2𝑛−1 = 4 ∗ 𝑚 ∗ 𝑎 ∗ 𝑥 2𝑚 𝑦; 𝑚 = 1, 𝑛 = 1
2𝑏 = 4𝑎 → 𝑏 = 2𝑎
𝑉 (𝑥, 𝑦) = 𝑎𝑥 2 + 2𝑎𝑦 2
𝜕𝑉 𝑑𝑥 𝜕𝑉 𝑑𝑦
𝑉 ′ (𝑥, 𝑦) = + = [2𝑎𝑥 ][−2𝑥𝑦] + [4𝑎𝑦][𝑥 2 − 𝑦 3 ] = −4𝑎𝑥 2 𝑦 + 4𝑎𝑥 2 𝑦 − 4𝑎𝑦 4 = −4𝑎𝑦 4
𝜕𝑥 𝑑𝑡 𝜕𝑦 𝑑𝑡
𝑥3
𝑥 ′ = 𝑥𝑦 2 − 2
b) 𝑥 2𝑦 𝑦3
𝑦′ = −
5 2
′(
𝜕𝑉 𝑑𝑥 𝜕𝑉 𝑑𝑦 2𝑚−1 ] [𝑥𝑦 2
𝑥3 2𝑛−1
𝑥 2𝑦 𝑦 3
𝑉 𝑥, 𝑦) = + [
= 2𝑚𝑎𝑥 [
− ] + 2𝑛𝑏𝑦 ] [ − ]
𝜕𝑥 𝑑𝑡 𝜕𝑦 𝑑𝑡 2 5 2
2
2𝑚𝑎𝑥 2𝑚 𝑦 2 − 𝑚𝑎𝑥 2𝑚+2 + 𝑛𝑏𝑥 2 𝑦 2𝑛 − 𝑛𝑏𝑦 2𝑛+2
5
𝑚 = 2, 𝑛 = 2, 𝑎 = 1, 𝑏 = 5
Jesús Abraham Rojas Úrzulo Licenciatura en Matemáticas
Matricula: ES1821013126 Ecuaciones Diferenciales II (MEDI2)
−𝑥 4 + 4𝑥 2 𝑦 2 − 5𝑦 4 → 𝑢 = 𝑥 2 , 𝑣 = 𝑦 2 → −𝑢2 + 4𝑢𝑣 − 5𝑣 2
𝜕 𝜕
(−𝑢2 + 4𝑢𝑣 − 5𝑣 2 ) = −2𝑢 + 4𝑣 = 0; (−𝑢2 + 4𝑢𝑣 − 5𝑣 2 ) = 4𝑢 − 10𝑣 = 0
𝜕𝑢 𝜕𝑣
El máximo se la función se da en 𝑢 = 𝑣 = 0, 𝑥 = 𝑦 = 0, por lo tanto
Es decir,
𝐹 (𝑥, 𝑦) = −𝑥 4 + 4𝑥 2 𝑦 2 − 5𝑦 4 < 0, (𝑥, 𝑦) ≠ (0,0), 𝐹 (0,0) = 0
Por lo tanto, 𝑉 (𝑥, 𝑦) > 0, 𝑉(0,0) = 0 y 𝑉 ′ (𝑥, 𝑦) < 0, 𝑉′(0,0) = 0
Por lo tanto, es asintóticamente estable, para, 𝑉 (𝑥, 𝑦) = 𝑥 2 + 5𝑦 2
𝑥 ′ = −2𝑥 + 𝑥𝑦 3
a)
𝑦 ′ = −𝑥 2 𝑦 2 − 𝑦 3
𝜕𝑉 𝑑𝑥 𝜕𝑉 𝑑𝑦
𝑉 ′ (𝑥, 𝑦) = + = [2𝑚𝑎𝑥 2𝑚−1 ][−2𝑥 + 𝑥𝑦 3 ] + [2𝑛𝑏𝑦 2𝑛−1 ][−𝑥 2 𝑦 2 − 𝑦 3 ]
𝜕𝑥 𝑑𝑡 𝜕𝑦 𝑑𝑡
−4𝑚𝑎𝑥 2𝑚 + 2𝑚𝑎𝑥 2𝑚 𝑦 3 − 2𝑛𝑏𝑥 2 𝑦 2𝑛+1 − 2𝑛𝑏𝑦 2𝑛+2
Los términos a eliminar:
2𝑚𝑎𝑥 2𝑚 𝑦 3 = 2𝑛𝑏𝑥 2 𝑦 2𝑛+1
𝑛 = 1, 𝑚 = 1
2𝑎𝑥 2 𝑦 3 = 2𝑏𝑥 2 𝑦 3 ; 𝑎 = 𝑏
𝑉 (𝑥, 𝑦) = 𝑎𝑥 2 + 𝑎𝑦 2
𝜕𝑉 𝑑𝑥 𝜕𝑉 𝑑𝑦
𝑉 ′ (𝑥, 𝑦) = + = [2𝑎𝑥][−2𝑥 + 𝑥𝑦 3 ] + [2𝑎𝑦][−𝑥 2 𝑦 2 − 𝑦 3 ]
𝜕𝑥 𝑑𝑡 𝜕𝑦 𝑑𝑡
Jesús Abraham Rojas Úrzulo Licenciatura en Matemáticas
Matricula: ES1821013126 Ecuaciones Diferenciales II (MEDI2)
3. Mediante un programa científico obtenga las órbitas para el oscilador de Var Der Pol
𝑑2𝑥 𝑑𝑥
2
− μ (1 − 𝑥 2 ) +𝑥 = 0
𝑑𝑡 𝑑𝑡
𝑑𝑥 𝑑 2 𝑥 𝑑𝑦
= 𝑦, =
𝑑𝑡 𝑑𝑡 2 𝑑𝑡
𝑑𝑦
= μ ( 1 − 𝑥 2 )𝑦 − 𝑥
𝑑𝑡
Conclusiones:
El análisis cualitativo de los sistemas de ecuaciones diferenciales es de suma importancia, en la
primera unidad de la materia el propósito fue resolver sistemas lineales de ecuaciones diferenciales,
es decir, encontrar la solución explicita para cada variable, sin embargo, en la segunda unidad el
objetivo era analizar cualitativamente dichos sistemas, es decir, encontrar los puntos de equilibrio y
clasificarlos, así como definir si el sistema era estable o inestable.
En esta tercera unidad, el objetivo fue realizar el análisis cualitativo a sistemas de ecuaciones no
lineales, ya que dichos sistemas, en la mayoría de las ocasiones es muy difícil encontrar una solución
explicita, sin embargo, analizando con diferentes técnicas es posible entender la evolución del sistema,
Jesús Abraham Rojas Úrzulo Licenciatura en Matemáticas
Matricula: ES1821013126 Ecuaciones Diferenciales II (MEDI2)
esto mediante la linealización en los puntos de equilibrio, la ecuación de Lyapunov para determinar la
estabilidad, los ciclos límites e inclusive los atractores extraños; en esta evidencia de aprendizaje no
se resolvió explícitamente ninguna ecuación diferencial, sin embargo, en todas aplicando los métodos
anteriormente mencionados se tuvo un entendimiento de la evolución del sistema.
Bibliografía
Escutia, J. I. (2013). Matemáticas 5. México D.F.: Mc Graw Hill.
Unadm. (s.f.). Ecuaciones Diferenciales II (Universidad Abierta y a Distancia de México). Universidad
Abierta y a Distancia de México.
Zill, D. G. (2006). Ecuaciones diferenciales con aplicaciones de modelado. México D.F.: THOMSON.