3,1 EDO - (Orden 2, Superior & Soluciones in SP) - (Teoría) - 1er&2do Parcial
3,1 EDO - (Orden 2, Superior & Soluciones in SP) - (Teoría) - 1er&2do Parcial
3,1 EDO - (Orden 2, Superior & Soluciones in SP) - (Teoría) - 1er&2do Parcial
Profesor: Antonio Chong Escobar, Ph.D. Referencia Bibliográfica: Ecuaciones diferenciales y problemas con valores en la frontera, Boyce – Diprima, 4ta edición
__________________________________________________________________EDO (ORDEN 2 Y SUPERIOR & SOLUCIONES EN SERIES DE POTENCIAS)
৺DEDUCCIÓN DE UNA SOLUCIÓN LINEALMENTE INDEPENDIENTE (𝑳. 𝑰.) A UNA SOLUCIÓN CONOCIDA PARA UNA EDO LINEAL DE 2DO
ORDEN HOMOGÉNEA
Sea 𝑦1 (𝑥) una solución de la ecuación: 𝑦 ′′ + 𝑝(𝑥)𝑦 ′ + 𝑞(𝑥)𝑦 = 0, es decir, 𝑦1 ′′ + 𝑝(𝑥)𝑦1 ′ + 𝑞(𝑥)𝑦1 = 0.
Una solución 𝐿. 𝐼. a 𝑦1 (𝑥) está dada por 𝑦2 (𝑥) = 𝑣(𝑥)𝑦1 (𝑥), tal que 𝑣(𝑥) es una función no constante por hallar.
Por lo tanto, 𝑦2′ (𝑥) = 𝑣 𝑦1 ′ + 𝑣′𝑦1
𝑦2′′ (𝑥) = 𝑣 𝑦1′′ + 𝑣 ′ 𝑦1′ + 𝑣 ′′ 𝑦1 + 𝑣 ′ 𝑦1′ .
2
Sustituyendo 𝑦2 (𝑥) ; 𝑦2′ (𝑥) y 𝑦2′′ (𝑥) en la ecuación diferencial se tiene
𝑦2′′ + 𝑝(𝑥)𝑦2′ + 𝑞(𝑥)𝑦2 = 0
`
𝑣 𝑦1′′ + 𝑣 ′ 𝑦1′ + 𝑣 ′′ 𝑦1 + 𝑣 ′ 𝑦1′ + 𝑝(𝑥)(𝑣 𝑦1′ + 𝑣 ′ 𝑦1 ) + 𝑞(𝑥)𝑣 𝑦1 = 0
𝑣⏟ [𝑦1′′ + 𝑝(𝑥)𝑦1′ + 𝑞(𝑥)𝑦1 ] + 𝑣 ′ [2𝑦1′ + 𝑝(𝑥)𝑦1 ] + 𝑣 ′′ 𝑦1 = 0
𝑖𝑔𝑢𝑎𝑙 𝑎 0 𝑑𝑎𝑑𝑜 𝑞𝑢𝑒 𝑦1 𝑒𝑠 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛
𝑣 ′ [2𝑦1′ + 𝑝(𝑥)𝑦1 ] + 𝑣 ′′ 𝑦1 = 0 (EDO de 2do orden en la que está ausenta la variable dependiente 𝑣).
Al realizar las sustituciones 𝑣 ′ = 𝑤 ; 𝑣 ′′ = 𝑤 ′ para reducir del orden de la EDO se tiene
𝑑𝑤
𝑤[2𝑦1′ + 𝑝(𝑥)𝑦1 ] + 𝑦1 =0
𝑑𝑥
𝑑𝑤 𝑑𝑤 2𝑦′1 +𝑝(𝑥)𝑦1 𝑦′1
→ 𝑦1 = −𝑤[2𝑦1′ + 𝑝(𝑥)𝑦1 ] → ∫ = −∫ 𝑑𝑥 → 𝑙𝑛(𝑤) = − ∫ (2 + 𝑝(𝑥)) 𝑑𝑥
𝑑𝑥 𝑤 𝑦1 𝑦1
−2 −2
→ 𝑙𝑛(𝑤) = −2𝑙𝑛(𝑦1 ) − ∫ 𝑝(𝑥)𝑑𝑥 → 𝑤 = 𝑒 𝑙𝑛(𝑦1 ) −∫ 𝑝(𝑥)𝑑𝑥
→ 𝑤(𝑥) = (𝑦1 (𝑥)) 𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 .
Regresando a la variable original se tiene:
𝑑𝑣 𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑒− ∫ 𝑝(𝑥)𝑑𝑥
= 2 → ∫ 𝑑𝑣 = ∫ 2 𝑑𝑥 → 𝑣(𝑥) = ∫ 2 𝑑𝑥.
𝑑𝑥
(𝑦1 (𝑥)) (𝑦1 (𝑥)) (𝑦1 (𝑥))
Entonces una solución 𝐿. 𝐼. a la solución conocida 𝑦1 (𝑥) es
𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑒− ∫ 𝑝(𝑥)𝑑𝑥
𝑦2 (𝑥) = 𝑦1 (𝑥) ∫ 2 𝑑𝑥 . Así, la solución general de la EDO es: 𝑦(𝑥) = 𝑐1 𝑦1 (𝑥) + 𝑐2 𝑦1 (𝑥) ∫ 2 𝑑𝑥 ; 𝑐1 , 𝑐2 ∈ ℝ.
(𝑦1 (𝑥)) (𝑦1 (𝑥))
৺TEOREMA DE ABEL
Si 𝑦1 y 𝑦2 son soluciones de la EDO de 2do orden lineal homogénea 𝑦 ′′ + 𝑝(𝑥)𝑦 ′ + 𝑞(𝑥)𝑦 = 0, tal que 𝑝 y 𝑞 son continuas en un intervalo
abierto 𝐼, entonces el Wronskiano de 𝑦1 y 𝑦2 está dado por 𝑊(𝑦1 , 𝑦2 )(𝑥) = 𝑐𝑒 − ∫ 𝑝(𝑥)𝑑𝑥, tal que 𝑐 ∈ ℝ y 𝑐 depende de 𝑦1 y 𝑦2 pero no
depende de 𝑥.
Demostración:
𝐻1 : 𝑦1 es solución de la EDO 𝑦1 ′′ + 𝑝(𝑥)𝑦1 ′ + 𝑞(𝑥)𝑦1 = 0.
𝐻2 : 𝑦2 es solución de la EDO 𝑦2 ′′ + 𝑝(𝑥)𝑦2 ′ + 𝑞(𝑥)𝑦2 = 0.
Al multiplicar la ecuación de 𝐻1 por (−𝑦2 ) se tiene −𝑦2 𝑦1′′ − 𝑦2 𝑝(𝑥)𝑦1′ − 𝑦2 𝑞(𝑥)𝑦1 = 0.
Al multiplicar la ecuación de 𝐻2 por (𝑦1 ) se tiene 𝑦1 𝑦2′′ + 𝑦1 𝑝(𝑥)𝑦2′ + 𝑦1 𝑞(𝑥)𝑦2 = 0.
Sumando estas 2 ecuaciones: 𝑦1 𝑦2′′ − 𝑦1′′ 𝑦2 + 𝑝(𝑥) [𝑦 ′ ′
⏟1 𝑦2 − 𝑦1 𝑦2 ] = 0.
𝑊(𝑦1 ,𝑦2 )
Se puede verificar que los 2 primeros términos, esto es, 𝑦1 𝑦2′′ − 𝑦1′′ 𝑦2 , corresponden a la derivada de 𝑊(𝑦1 , 𝑦2 ) como sigue:
𝑦1 𝑦2
𝑊(𝑦1 , 𝑦2 )(𝑥) = |𝑦 ′ 𝑦 ′ | = 𝑦1 𝑦2′ − 𝑦1′ 𝑦2 → 𝑊 ′ (𝑦1 , 𝑦2 )(𝑥) = (𝑦1 𝑦2′ − 𝑦1′ 𝑦2 )′ = (𝑦1′ 𝑦2′ + 𝑦1 𝑦2′′ ) − (𝑦1′′ 𝑦2 + 𝑦1′ 𝑦2′ ) = 𝑦1 𝑦2′′ − 𝑦1′′ 𝑦2 .
1 2
Así, representando a 𝑊(𝑦1 , 𝑦2 )(𝑥) con 𝑊 y representando a 𝑊 ′ (𝑦1 , 𝑦2 )(𝑥) con 𝑊 ′ se tiene la ecuación 𝑊 ′ + 𝑝(𝑥)𝑊 = 0.
Resolviendo esta última EDO:
𝑑𝑊 𝑑𝑊
= −𝑝(𝑥)𝑊 → ∫ = − ∫ 𝑝(𝑥)𝑑𝑥 → 𝑙𝑛|𝑊| = − ∫ 𝑝(𝑥)𝑑𝑥 + 𝑘 → |𝑊| = 𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 𝑒 𝑘 → 𝑊 = ±𝑐𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 ; 𝑐 > 0.
𝑑𝑥 𝑊
Como 𝑊 = 0 también satisface la ecuación 𝑊 ′ + 𝑝(𝑥)𝑊 = 0, entonces 𝑊(𝑦1 , 𝑦2 )(𝑥) = 𝑐𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 ; 𝑐 ∈ ℝ.
Por lo tanto, el Wronskiano de conjuntos fundamentales de soluciones de una EDO lineal de 2do orden homogénea son múltiplos
escalares. Además, sin resolver este tipo de EDO, el Wronskiano de cualquier conjunto fundamental de soluciones puede determinarse
en términos de una constante.
৺DEDUCCIÓN DE UNA SOLUCIÓN LINEALMENTE INDEPENDIENTE (𝑳. 𝑰.) A UNA SOLUCIÓN CONOCIDA PARA UNA EDO LINEAL DE 2DO
ORDEN HOMOGÉNEA USANDO EL TEOREMA DE ABEL
Considerando el teorema de Abel con 𝑐 ≠ 0, una solución 𝑦1 conocida y una solución 𝑦2 por determinar se tiene
𝑦1 𝑦2 𝑦′ 𝑐𝑒− ∫ 𝑝(𝑥)𝑑𝑥
𝑊(𝑦1 , 𝑦2 ) = 𝑐𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 → |𝑦 ′ 𝑦 ′ | = 𝑐𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 → 𝑦1 𝑦2′ − 𝑦1′ 𝑦2 = 𝑐𝑒 − ∫ 𝑝(𝑥)𝑑𝑥 → 𝑦2′ − 1 𝑦2 = (EDO lineal para 𝑦2 ).
1 2 𝑦1 𝑦1
𝑦′
∫ −𝑦1 𝑑𝑥
Para esta última ecuación, el factor integrante está dado por 𝑢(𝑥) = 𝑒 1 = 𝑒 −𝑙𝑛(𝑦1) = 𝑦1 −1 .
Multiplicando la EDO lineal por el factor integrante 𝑢(𝑥):
𝑦′1 𝑐𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑐𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑐𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑒− ∫ 𝑝(𝑥)𝑑𝑥
𝑦1 −1 𝑦2′ − 2 𝑦2 = 2 → 𝑑(𝑦2 𝑦1 −1 ) = 2 𝑑𝑥 → 𝑦2 𝑦1 −1 = ∫ 2 𝑑𝑥 → 𝑦2 = 𝑐𝑦1 ∫ 2 𝑑𝑥.
(𝑦1 ) (𝑦1 ) (𝑦1 ) (𝑦1 ) (𝑦1 )
Por consiguiente, una solución 𝐿. 𝐼. a una solución 𝑦1 conocida para la EDO 𝑦 ′′ + 𝑝(𝑥)𝑦 ′ + 𝑞(𝑥)𝑦 = 0, tal que 𝑝 y 𝑞 son continuas en un
𝑒− ∫ 𝑝(𝑥)𝑑𝑥
intervalo abierto 𝐼, está dada por 𝑦2 (𝑥) = 𝑐𝑦1 (𝑥) ∫ 2 𝑑𝑥 ; 𝑐 ≠ 0 (esta expresión se conoce como la identidad de Abel). Por
(𝑦1 (𝑥))
𝑒− ∫ 𝑝(𝑥)𝑑𝑥
ejemplo, con 𝑐 = 1: 𝑦2 (𝑥) = 𝑦1 (𝑥) ∫ 2 𝑑𝑥.
(𝑦1 (𝑥))
Profesor: Antonio Chong Escobar, Ph.D. Referencia Bibliográfica: Ecuaciones diferenciales y problemas con valores en la frontera, Boyce – Diprima, 4ta edición
__________________________________________________________________EDO (ORDEN 2 Y SUPERIOR & SOLUCIONES EN SERIES DE POTENCIAS)
৺EDO LINEAL DE 2DO ORDEN HOMOGÉNEA DE COEFICIENTES CONSTANTES
Es una EDO de la forma 𝑎𝑦 ′′ (𝑥) + 𝑏𝑦 ′ (𝑥) + 𝑐𝑦(𝑥) = 0, tal que 𝑎, 𝑏, 𝑐 ∈ ℝ ∧ 𝑎 ≠ 0. Para resolver esta EDO se plantea la solución de la
forma 𝑦(𝑥) = 𝑒 𝑟𝑥 tal que 𝑟 es una constante por determinar.
Para determinar el valor de 𝑟 se sustituye 𝑦(𝑥) = 𝑒 𝑟𝑥 y sus derivadas (𝑦 ′ (𝑥) = 𝑟𝑒 𝑟𝑥 y 𝑦 ′′ (𝑥) = 𝑟 2 𝑒 𝑟𝑥 ) en la EDO como sigue:
𝑎𝑟 2 𝑒 𝑟𝑥 + 𝑏𝑟𝑒 𝑟𝑥 + 𝑐𝑒 𝑟𝑥 = 0 → 𝑒 𝑟𝑥 (𝑎𝑟 2 + 𝑏𝑟 + 𝑐) = 0 → (𝑒 𝑟𝑥 = 0) ∨ (𝑎𝑟 2 + 𝑏𝑟 + 𝑐 = 0).
La ecuación exponencial obtenida no tiene solución, mientras que la ecuación cuadrática obtenida tiene dos soluciones y presenta tres
3
casos posibles. Así, el planteamiento es válido siempre que 𝑟 sea solución de esta ecuación cuadrática, a la cual se la denomina ecuación
`
característica o auxiliar.
Soluciones de acuerdo con la ecuación auxiliar
Caso 1: Si el discriminante 𝑏 2 − 4𝑎𝑐 > 0, entonces la ecuación auxiliar tiene dos soluciones reales distintas 𝑟1 y 𝑟2 , con lo cual las funciones
𝑦1 (𝑥) = 𝑒 𝑟1𝑥 y 𝑦2 (𝑥) = 𝑒 𝑟2𝑥 son soluciones de la EDO. Estas soluciones forman un conjunto fundamental de soluciones para la EDO, lo
cual se muestra a continuación calculando 𝑊(𝑦1 , 𝑦2 ) y verificando que es diferente de cero como se muestra a continuación.
𝑒 𝑟1𝑥 𝑒 𝑟2𝑥
𝑊(𝑦1 , 𝑦2 )(𝑥) = | 𝑟1𝑥 | = 𝑟2 𝑒 (𝑟1+𝑟2)𝑥 − 𝑟1 𝑒 (𝑟1+𝑟2)𝑥 = (𝑟1 − 𝑟2 )𝑒 (𝑟1+𝑟2)𝑥 ≠ 0, puesto que 𝑟1 ≠ 𝑟2 .
𝑟1 𝑒 𝑟2 𝑒 𝑟2𝑥
Entonces la solución general de la EDO está dada por: 𝑦(𝑥) = 𝑐1 𝑒 𝑟1𝑥 + 𝑐2 𝑒 𝑟2𝑥 ; 𝑐1 , 𝑐2 ∈ ℝ.
Caso 2: Si el discriminante 𝑏 2 − 4𝑎𝑐 = 0, la ecuación auxiliar tiene solución real repetida 𝑟 y la función 𝑦1 (𝑥) = 𝑒 𝑟𝑥 es solución de la
EDO. Para hallar una segunda solución linealmente independiente a 𝑦1 se puede utilizar la identidad de Abel, como se muestra a
continuación.
𝑒− ∫ 𝑝(𝑥)𝑑𝑥 𝑏 𝑐 𝑏
𝑦2 (𝑥) = 𝑦1 (𝑥) ∫ 2 𝑑𝑥 , considerando la EDO de forma canónica: 𝑦 ′′ (𝑥) + 𝑦 ′ (𝑥) + 𝑦(𝑥) = 0, tal que 𝑝(𝑥) =
(𝑦1 (𝑥)) 𝑎 𝑎 𝑎
𝑏
𝑒− ∫ 𝑎 𝑑𝑥 −𝑏±√𝑏2 −4𝑎𝑐 𝑏 𝑏
→ 𝑦2 (𝑥) = 𝑒 𝑟𝑥 ∫ [ 𝑟𝑥 ]2 𝑑𝑥 y además se conoce que para la ecuación auxiliar 𝑟 = =− , con lo cual 2𝑟 = − .
𝑒 2𝑎 2𝑎 𝑎
𝑒∫ 2𝑟𝑑𝑥 𝑒2𝑟𝑥
→ 𝑦2 (𝑥) = 𝑒 𝑟𝑥 ∫ [ 𝑑𝑥 = 𝑒 𝑟𝑥 ∫ 𝑑𝑥 = 𝑒 𝑟𝑥 ∫ 𝑑𝑥 = 𝑥𝑒 𝑟𝑥 .
𝑒𝑟𝑥 ]2 𝑒2𝑟𝑥
Entonces la solución general de la EDO está dada por: 𝑦(𝑥) = 𝑐1 𝑒 𝑟𝑥 + 𝑐2 𝑥𝑒 𝑟𝑥 ; 𝑐1 , 𝑐2 ∈ ℝ.
Caso 3: Si el discriminante 𝑏 2 − 4𝑎𝑐 < 0, la ecuación auxiliar tiene soluciones complejas conjugadas: 𝑟1,2 = 𝛼 ± 𝜆𝑖 y las funciones
𝑦1 (𝑥) = 𝑒 𝑟1𝑥 y 𝑦2 (𝑥) = 𝑒 𝑟2𝑥 son soluciones de la EDO, es decir, las funciones
𝑦1 (𝑥) = 𝑒 (𝛼+𝜆𝑖)𝑥 = 𝑒 𝛼𝑥 𝑒 𝜆𝑥𝑖 = 𝑒 𝛼𝑥 (𝑐𝑜𝑠(𝜆𝑥) + 𝑖𝑠𝑒𝑛(𝜆𝑥)),
𝑦2 (𝑥) = 𝑒 (𝛼−𝜆𝑖)𝑥 = 𝑒 𝛼𝑥 𝑒 −𝜆𝑥𝑖 = 𝑒 𝛼𝑥 (𝑐𝑜𝑠(𝜆𝑥) − 𝑖𝑠𝑒𝑛(𝜆𝑥)).
Así, la combinación lineal 𝑦3 (𝑥) = 𝑐1 𝑦1 (𝑥) + 𝑐2 𝑦2 (𝑥) también es una solución de la EDO, pero no es la solución general debido a la
dependencia lineal que se puede probar entre 𝑦1 y 𝑦2 , así como a los valores imaginarios presentes. Sin embargo, se obtiene la solución
general de la EDO, a partir de los siguientes artificios:
𝑦3 (𝑥) = 𝑐1 𝑒 𝛼𝑥 (𝑐𝑜𝑠(𝜆𝑥) + 𝑖𝑠𝑒𝑛(𝜆𝑥)) + 𝑐2 𝑒 𝛼𝑥 (𝑐𝑜𝑠(𝜆𝑥) − 𝑖𝑠𝑒𝑛(𝜆𝑥))
𝑦3 (𝑥) = 𝑐1 𝑒 𝛼𝑥 𝑐𝑜𝑠(𝜆𝑥) + 𝑖𝑐1 𝑒 𝛼𝑥 𝑠𝑒𝑛(𝜆𝑥) + 𝑐2 𝑒 𝛼𝑥 𝑐𝑜𝑠(𝜆𝑥) − 𝑖𝑐2 𝑒 𝛼𝑥 𝑠𝑒𝑛(𝜆𝑥)
⏟1 + 𝑐2 ) 𝑒 𝛼𝑥 𝑐𝑜𝑠(𝜆𝑥) + 𝑖(𝑐
𝑦3 (𝑥) = (𝑐 ⏟ 1 − 𝑐2 ) 𝑒 𝛼𝑥 𝑠𝑒𝑛(𝜆𝑥).
𝑘1 𝑘2
Se puede mostrar que las funciones 𝑒 𝛼𝑥 𝑐𝑜𝑠(𝜆𝑥) y 𝑒 𝛼𝑥 𝑠𝑒𝑛(𝜆𝑥) son soluciones de la EDO y que además son linealmente independientes.
Entonces la solución general de la EDO está dada por: 𝑦(𝑥) = 𝑘1 𝑒 𝛼𝑥 𝑐𝑜𝑠(𝜆𝑥) + 𝑘2 𝑒 𝛼𝑥 𝑠𝑒𝑛(𝜆𝑥) ; 𝑘1 , 𝑘2 ∈ ℝ.
ii) exponencial: 𝑎𝑒 𝜆𝑥 ⏞ 𝑒 𝜆𝑥 ] 𝑥 𝑆
[𝐴
? ?
tal que 𝑝𝑛 (𝑥) y 𝑞𝑚 (𝑥) son polinomios de grados 𝑛 y 𝑚 tal que el grado de los polinomios 𝑃 y 𝑄 es: 𝑁 = 𝑚𝑎𝑥{𝑛, 𝑚}
? ?
tal que 𝑝𝑛 (𝑥) y 𝑞𝑚 (𝑥) son polinomios de grados 𝑛 y 𝑚 tal que el grado de los polinomios 𝑃 y 𝑄 es: 𝑁 = 𝑚𝑎𝑥{𝑛, 𝑚}
𝑆: es el menor entero no negativo (positivo o cero), tal que cada término de la solución particular planteada 𝑦𝑝 (𝑥) sea linealmente
independiente a cada término de la solución de la EDO homogénea correspondiente (la solución denominada complementaria).
Profesor: Antonio Chong Escobar, Ph.D. Referencia Bibliográfica: Ecuaciones diferenciales y problemas con valores en la frontera, Boyce – Diprima, 4ta edición
__________________________________________________________________EDO (ORDEN 2 Y SUPERIOR & SOLUCIONES EN SERIES DE POTENCIAS)
৺MÉTODO DE VARIACIÓN DE PARÁMETROS (PARA LAS EDO LINEALES DE 2DO ORDEN NO HOMOGÉNEAS)
Este método se puede aplicar tanto a las EDO no homogéneas de coeficientes constantes como a las de coeficientes variables. Además,
no es necesario que el término no homogéneo sea de algún tipo específico.
Considere la EDO no homogénea 𝑦 ′′ (𝑥) + 𝑝(𝑥)𝑦 ′ (𝑥) + 𝑞(𝑥)𝑦(𝑥) = 𝑔(𝑥), tal que 𝑔(𝑥) ≠ 0. El método de variación de parámetros
consiste en plantear la solución particular de la forma 𝑦𝑝 (𝑥) = 𝑣1 (𝑥)𝑌1 (𝑥) + 𝑣2 (𝑥)𝑌2 (𝑥), tal que 𝑣1 y 𝑣2 son parámetros variables
(funciones no constantes) que deben hallarse, mientras que 𝑌1 (𝑥) y 𝑌2 (𝑥) son soluciones linealmente independientes de la EDO
5
homogénea correspondiente.
`
De esta forma, se tiene 𝑦𝑝′ (𝑥) = 𝑣1′ 𝑌1 + 𝑣1 𝑌1′ + 𝑣2′ 𝑌2 + 𝑣2 𝑌2′ → 𝑦𝑝′ (𝑥) = 𝑣1 𝑌1′ + 𝑣2 𝑌2′ + [𝑣1′ 𝑌1 + 𝑣2′ 𝑌2 ].
Se hace la suposición de que 𝑣1′ 𝑌1 + 𝑣2′ 𝑌2 = 0 (1era condición que deben satisfacer 𝑣1 y 𝑣2 )
Por lo tanto, se tiene 𝑦𝑝′ (𝑥) = 𝑣1 𝑌1′ + 𝑣2 𝑌2′ ,
𝑦𝑝′′ (𝑥) = 𝑣1′ 𝑌1′ + 𝑣1 𝑌1′′ + 𝑣2′ 𝑌2′ + 𝑣2 𝑌2′′ .
Sustituyendo 𝑦𝑝 (𝑥), 𝑦𝑝 ′(𝑥) y 𝑦𝑝 ′′(𝑥) en la EDO no homogénea se obtiene
𝑦𝑝 ′′(𝑥) + 𝑝(𝑥)𝑦𝑝 ′(𝑥) + 𝑞(𝑥)𝑦𝑝 (𝑥) = 𝑔(𝑥)
[𝑣1′ 𝑌1′ + 𝑣1 𝑌1′′ + 𝑣2′ 𝑌2′ + 𝑣2 𝑌2′′ ] + 𝑝(𝑥)[𝑣1 𝑌1′ + 𝑣2 𝑌2′ ] + 𝑞(𝑥)[𝑣1 𝑌1 + 𝑣2 𝑌2 ] = 𝑔(𝑥)
𝑣1 (𝑥) [𝑌⏟1′′ + 𝑝(𝑥)𝑌1′ + 𝑞(𝑥)𝑌1 ] + 𝑣2 (𝑥) [𝑌 ′′ ′ ′ ′ ′ ′
⏟2 + 𝑝(𝑥)𝑌2 + 𝑞(𝑥)𝑌2 ] + 𝑣1 𝑌1 + 𝑣2 𝑌2 = 𝑔(𝑥).
0 0
(𝑌1′′ + 𝑝(𝑥)𝑌1′ + 𝑞(𝑥)𝑌1 = 0 y 𝑌2′′ + 𝑝(𝑥)𝑌2′ + 𝑞(𝑥)𝑌2 = 0 puesto que 𝑌1 y 𝑌2 son soluciones de la EDO homogénea correspondiente).
Entonces se tiene: 𝑣1′ 𝑌1′ + 𝑣2′ 𝑌2′ = 𝑔(𝑥) (2da condición que deben satisfacer 𝑣1 y 𝑣2 ).
Por lo tanto, para que todo lo planteado y supuesto sea válido, 𝑣1 y 𝑣2 deben satisfacer el sistema
𝑣 ′ 𝑌 + 𝑣2′ 𝑌2 = 0 𝑌1 𝑌2 𝑣1′ 0
{ 1′ 1′ , cuya representación matricial es [ ′ ][ ] = [ ].
𝑣1 𝑌1 + 𝑣2′ 𝑌2′ = 𝑔(𝑥) 𝑌1 𝑌2′ 𝑣2′ 𝑔(𝑥)
Así, las soluciones del sistema, utilizando la regla de Cramer, son las siguientes
0 𝑌2 𝑌1 0
| | | ′ |
𝑔(𝑥) 𝑌′2 −𝑌2 𝑔(𝑥) −𝑌 𝑔(𝑥) 𝑌
1 𝑔(𝑥) 𝑌 𝑔(𝑥) 𝑌 𝑔(𝑥)
𝑣1′ = 𝑌1 𝑌2 = → 𝑣1 = ∫ 𝑊(𝑌2 ,𝑌 ) 𝑑𝑥 ; 𝑣2′ = 𝑌1 𝑌2
1
= 𝑊(𝑌 → 𝑣2 = ∫ 1 𝑑𝑥.
𝑊(𝑌1 ,𝑌2 ) 1 2 1 ,𝑌2 ) 𝑊(𝑌
1 ,𝑌2 )
| ′ | | ′ |
𝑌′2
𝑌1 𝑌 1 𝑌′2
−𝑌2 𝑔(𝑥) 𝑌 𝑔(𝑥)
Finalmente, 𝑦𝑝 (𝑥) = 𝑌1 ∫ ( 𝑑𝑥 + 𝑌2 ∫ 1 𝑑𝑥.
𝑊 𝑌1 ,𝑌2 ) 𝑊(𝑌1 ,𝑌2 )
৺MÉTODO DE VARIACIÓN DE PARÁMETROS (PARA LAS EDO LINEALES DE ORDEN SUPERIOR NO HOMOGÉNEAS)
Sea {𝑌1 , 𝑌2 , … , 𝑌𝑛 } un conjunto fundamental de soluciones de la EDO homogénea
𝑦 (𝑛) + 𝑎𝑛−1 (𝑥)𝑦 (𝑛−1) + 𝑎𝑛−2 (𝑥)𝑦 (𝑛−2) + ⋯ + 𝑎1 (𝑥)𝑦 ′ + 𝑎0 (𝑥)𝑦 = 0,
entonces la solución particular de la EDO no homogénea: 𝑦 (𝑛) + 𝑎𝑛−1 (𝑥)𝑦 (𝑛−1) + 𝑎𝑛−2 (𝑥)𝑦 (𝑛−2) +. . . +𝑎1 (𝑥)𝑦 ′ + 𝑎0 (𝑥)𝑦 = 𝑔(𝑥),
usando el método de variación de parámetros, se plantea de la forma 𝑦𝑝 (𝑥) = 𝑣1 (𝑥)𝑌1 (𝑥) + 𝑣2 (𝑥)𝑌2 (𝑥) + ⋯ + 𝑣𝑛 (𝑥)𝑌𝑛 (𝑥), tal
que 𝑣1 (𝑥), 𝑣2 (𝑥), … 𝑣𝑛 (𝑥) son parámetros variables (funciones no constantes) que deben hallarse y que de manera análoga a lo deducido
para el caso de las EDO de 2do orden deben satisfacer el sistema de ecuaciones
𝑌1 𝑣1′ + 𝑌2 𝑣2′ + … + 𝑌𝑛 𝑣𝑛′ = 0 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 𝑣1′ 0
𝑌1′ 𝑣1′ + 𝑌2′ 𝑣2′ + … + 𝑌𝑛′ 𝑣𝑛′ = 0 𝑌1′ 𝑌2′ … 𝑌𝑛−1 ′
𝑌𝑛′ 𝑣2′ 0
⋮ ⋮ , cuya forma matricial es ⋮ ⋮ ⋮ ⋮ ⋮ = ⋮ .
(𝑛−2) ′ (𝑛−2) ′ (𝑛−2) ′ (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) ′
𝑌1 𝑣1 + 𝑌2 𝑣2 + ⋯ + 𝑌𝑛 𝑣𝑛 = 0 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 𝑣𝑛−1 0
(𝑛−1) ′ (𝑛−1) ′ (𝑛−1) ′ (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) [ 𝑣 ′ ] [𝑔(𝑥)]
{𝑌1 𝑣1 + 𝑌2 𝑣2 + ⋯ + 𝑌𝑛 𝑣𝑛 = 𝑔(𝑥) [ 𝑌1 𝑌2 … 𝑌 𝑌𝑛 𝑛−1 ] 𝑛
Así, las soluciones del sistema, utilizando la regla de Cramer, son las siguientes
0 𝑌2 … 𝑌𝑛−1 𝑌𝑛 𝑌1 0 … 𝑌𝑛−1 𝑌𝑛 𝑌1 𝑌2 … 𝑌𝑛−1 0
0| 𝑌′2
… 𝑌′𝑛−1 𝑌′𝑛 | 0 … | 𝑌′1 𝑌′𝑛−1 𝑌′𝑛 | | 𝑌′1 𝑌′2 … 𝑌′𝑛−1 0 |
⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮
(𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2)
| 0 𝑌2 … 𝑌𝑛−1 𝑌𝑛 | | 𝑌1 0 … 𝑌𝑛−1 𝑌𝑛 | | 𝑌1 𝑌2 … 𝑌𝑛−1 0 |
(𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1)
𝑔(𝑥) 𝑌2 … 𝑌𝑛−1 𝑌𝑛 𝑌 𝑔(𝑥) … 𝑌𝑛−1 𝑌𝑛 𝑌 𝑌2 … 𝑌𝑛−1 𝑔(𝑥)
𝑣1′ (𝑥) = 𝑌 𝑌 … 𝑌 𝑌 , 𝑣2′ (𝑥) = 𝑌1 𝑌 … 𝑌 𝑌 , … , 𝑣𝑛′ (𝑥) = 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 .
1 2 𝑛−1 𝑛 1 2 𝑛−1 𝑛 1
| 𝑌′1 𝑌′2 … 𝑌′𝑛−1 𝑌′𝑛 | | 𝑌′1 𝑌′2 … 𝑌′𝑛−1 𝑌′𝑛 | | 𝑌′1 𝑌′2 … 𝑌′𝑛−1 𝑌′𝑛 |
⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮
(𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2) (𝑛−2)
| 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 | | 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 | | 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 |
(𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1)
𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛 𝑌1 𝑌2 … 𝑌𝑛−1 𝑌𝑛
Profesor: Antonio Chong Escobar, Ph.D. Referencia Bibliográfica: Ecuaciones diferenciales y problemas con valores en la frontera, Boyce – Diprima, 4ta edición
__________________________________________________________________EDO (ORDEN 2 Y SUPERIOR & SOLUCIONES EN SERIES DE POTENCIAS)
SOLUCIONES DE EDO UTILIZANDO SERIES DE POTENCIAS
Se considera las EDO lineales de 2do orden de la forma 𝑃(𝑥)𝑦 ′′ (𝑥) + 𝑄(𝑥)𝑦 ′ (𝑥) + 𝑅(𝑥)𝑦(𝑥) = 𝐺(𝑥).
Una gran cantidad de problemas de la física matemática llevan a ecuaciones de la forma 𝑃(𝑥)𝑦 ′′ (𝑥) + 𝑄(𝑥)𝑦 ′ (𝑥) + 𝑅(𝑥)𝑦(𝑥) = 0.
Como ejemplos se puede citar
1) la ecuación de Bessel (de orden 𝑣): 𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ + (𝑥 2 − 𝑣 2 )𝑦 = 0, tal que 𝑣 ∈ ℝ.
7
2) la ecuación de Legendre: (1 − 𝑥 2 )𝑦 ′′ − 2𝑥𝑦 ′ + 𝛼(𝛼 + 1)𝑦 = 0, tal que 𝛼 ∈ ℝ.
3) la ecuación de Hermite: 𝑦 ′′ − 2𝑥𝑦 ′ + 𝜆𝑦 = 0, tal que 𝜆 ∈ ℝ. `
4) la ecuación de Chebyshev: (1 − 𝑥 )𝑦 ′′ − 𝑥𝑦 ′ + 𝛼 2 𝑦 = 0,
2
tal que 𝛼 ∈ ℝ.
A continuación, se considera sólo a las ecuaciones homogéneas puesto que las no homogéneas se analizan de manera similar.
Tipos de puntos alrededor de los cuales se buscan soluciones en series de potencias
Considerando que 𝑃(𝑥), 𝑄(𝑥) y 𝑅(𝑥) son polinomios y que no hay factor que sea común para ellos se tiene que:
i) Un punto 𝑥0 tal que 𝑃(𝑥0 ) ≠ 0 se denomina punto ordinario. Dado que 𝑃(𝑥) es continua, existe un intervalo alrededor de 𝑥0 en el que
𝑃(𝑥) nunca es cero. En dicho intervalo es factible dividir a la EDO 𝑃(𝑥)𝑦 ′′ (𝑥) + 𝑄(𝑥)𝑦 ′ (𝑥) + 𝑅(𝑥)𝑦(𝑥) = 0 para 𝑃(𝑥) como sigue:
𝑄(𝑥) 𝑅(𝑥)
𝑦 ′′ (𝑥) + 𝑦 ′ (𝑥) + 𝑦(𝑥) = 0,
⏟
𝑃(𝑥) ⏟
𝑃(𝑥)
𝑝(𝑥) 𝑞(𝑥)
tal que 𝑝(𝑥) y 𝑞(𝑥) son continuas, por lo que, según el teorema de existencia y unicidad, en ese intervalo existe una solución única para
la EDO sujeta a las condiciones iniciales 𝑦(𝑥0 ) = 𝑦0 , 𝑦 ′ (𝑥0 ) = 𝑦0′ , para valores arbitrarios de 𝑦0 , 𝑦0′ .
?
Alrededor de un punto ordinario se buscan soluciones de la forma 𝑦(𝑥) = ∑+∞ ⏞𝑛 (𝑥 − 𝑥0 )𝑛 = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )2 + ⋯ ,
𝑛=0 𝑎
tal que se supone que la serie converge en un intervalo de la forma |𝑥 − 𝑥0 | < 𝛿, para algún 𝛿 > 0. Por consiguiente, la solución general
tiene la forma 𝑦(𝑥) = 𝑎0 𝑦1 (𝑥) + 𝑎1 𝑦2 (𝑥) ; 𝑎0 , 𝑎1 ∈ ℝ.
ii) Un punto 𝑥0 tal que 𝑃(𝑥0 ) = 0 se denomina punto singular. En este caso, 𝑄(𝑥0 ) ≠ 0 ó 𝑅(𝑥0 ) ≠ 0 y al menos uno de los coeficientes
𝑝(𝑥) o 𝑞(𝑥) de la ecuación de forma canónica se vuelve no acotado (discontinuo) cuando 𝑥 se acerca a 𝑥0 . Por lo tanto, el teorema de
existencia y unicidad no es aplicable en este caso y se utiliza la teoría de Frobenius.
TEORÍA DE FROBENIUS
Esta teoría hace uso de desarrollos en series de potencias alrededor de puntos singulares denominados regulares.
Dada la ecuación diferencial ordinaria 𝑃(𝑥)𝑦 ′′ + 𝑄(𝑥)𝑦′ + 𝑅(𝑥)𝑦 = 0,
(𝑥−𝑥0 )𝑄(𝑥) (𝑥−𝑥0 )2 𝑅(𝑥)
1) se dice que un punto singular 𝑥0 es regular si: lím existe y lím existe (más adelante, el resultado del primer
𝑥→𝑥0 𝑃(𝑥) 𝑥→𝑥0 𝑃(𝑥)
límite es denominado 𝑓0 y el resultado del segundo límite es denominado 𝑔0 ).
2) cualquier punto singular que no sea regular se denomina punto singular irregular.
Alrededor de un punto singular regular se buscan soluciones linealmente independientes de la forma
? ?
𝑦1 (𝑥) = (𝑥 − 𝑥0 )𝑟1 ∑+∞ +∞ ⏞ (𝑥
⏞𝑛 (𝑥 − 𝑥0 )𝑛 = ∑𝑛=0
𝑛=0 𝑎 𝑎𝑛 − 𝑥0 )𝑛+𝑟1 ,
? ?
𝑦2 (𝑥) = (𝑥 − 𝑥0 )𝑟2 ∑+∞ ⏞ 𝑛 +∞ ⏞
𝑛=0 𝑏𝑛 (𝑥 − 𝑥0 ) = ∑𝑛=0 𝑏𝑛 (𝑥 − 𝑥0 )
𝑛+𝑟2 .
?
Para esto, primero se plantea una solución de la forma 𝑦(𝑥) = ∑+∞ ⏞𝑛 (𝑥 − 𝑥0 )𝑛+𝑟 .
𝑛=0 𝑎
Luego, se deduce una ecuación cuadrática para 𝑟, llamada ecuación indicial. Se puede demostrar que la ecuación indicial tiene la forma
𝑟(𝑟 − 1) + 𝑓0 𝑟 + 𝑔0 = 0. Una vez hallados los valores de 𝑟 se plantean las soluciones 𝑦1 (𝑥) y 𝑦2 (𝑥), para las cuales se debe encontrar
los coeficientes 𝑎𝑛 y 𝑏𝑛 respectivamente.
Así, la solución general de la EDO analizada está dada por 𝑦(𝑥) = 𝑎0 𝑦1 (𝑥) + 𝑏0 𝑦2 (𝑥); 𝑎0 , 𝑏0 ∈ ℝ.
Casos de acuerdo con la ecuación indicial
Caso 1: Si se tienen dos raíces reales 𝑟1 y 𝑟2 tal que 𝑟1 > 𝑟2 y (𝑟1 − 𝑟2 ) ∉ ℤ, entonces las soluciones linealmente independientes de la
EDO se plantean de la forma
𝑦1 (𝑥) = ∑+∞𝑛=0 𝑎𝑛 (𝑥 − 𝑥0 )
𝑛+𝑟1 ; 𝑦 (𝑥) = ∑+∞ 𝑏 (𝑥 − 𝑥 )𝑛+𝑟2 .
2 𝑛=0 𝑛 0
Caso 2: Si se tiene raíz real repetida 𝑟, entonces las soluciones linealmente independientes de la EDO se plantean de la forma
𝑦1 (𝑥) = ∑+∞𝑛=0 𝑎𝑛 (𝑥 − 𝑥0 )
𝑛+𝑟 ; 𝑦 (𝑥) = 𝑦 (𝑥) 𝑙𝑛(𝑥 − 𝑥 ) + ∑+∞ 𝑏 (𝑥 − 𝑥 )𝑛+𝑟 .
2 1 0 𝑛=0 𝑛 0
Caso 3: Si se tienen dos raíces reales 𝑟1 y 𝑟2 tal que 𝑟1 > 𝑟2 y (𝑟1 − 𝑟2 ) ∈ ℤ, entonces las soluciones linealmente independientes de la
EDO se plantean de la forma
𝑦1 (𝑥) = ∑+∞𝑛=0 𝑎𝑛 (𝑥 − 𝑥0 )
𝑛+𝑟1 ; 𝑦 (𝑥) = 𝑎𝑦 (𝑥) 𝑙𝑛(𝑥 − 𝑥 ) + ∑+∞ 𝑏 (𝑥 − 𝑥 )𝑛+𝑟2 , tal que 𝑎 es una constante por hallar.
2 1 0 𝑛=0 𝑛 0
Observación: Si la primera solución obtenida 𝑦1 converge a una función conocida, entonces para hallar una segunda solución 𝑦2
linealmente independiente a 𝑦1 se puede utilizar la identidad de Abel.
Profesor: Antonio Chong Escobar, Ph.D. Referencia Bibliográfica: Ecuaciones diferenciales y problemas con valores en la frontera, Boyce – Diprima, 4ta edición