Soluciones de ecuaciones diferenciales mediante series

7.4 Soluciones de ecuaciones diferenciales mediante series

Objetivos de aprendizaje

  • 7.4.1 Utilizar las series de potencias para resolver ecuaciones diferenciales de primer y segundo orden.

En Introducción a series de potencias estudiamos cómo se pueden representar las funciones como series de potencias, y(x)=∑n=0∞anxn.y(x)=∑n=0∞anxn. También vimos que podemos hallar representaciones en serie de las derivadas de dichas funciones diferenciando la serie de potencias término a término. Esto da y′(x)=∑n=1∞nanxn–1y′(x)=∑n=1∞nanxn–1 y y″(x)=∑n=2 ∞n(n–1)anxn–2 .y″(x)=∑n=2 ∞n(n–1)anxn–2 . En algunos casos, estas representaciones en serie de potencias pueden utilizarse para hallar soluciones a ecuaciones diferenciales.

Observe que este tema se trata muy brevemente en este texto. La mayoría de los libros de texto de introducción a las ecuaciones diferenciales incluyen un capítulo entero sobre soluciones de series de potencias. Este texto solo tiene una sección sobre el tema, por lo que no se abordan aquí varias cuestiones importantes, sobre todo las relacionadas con la existencia de soluciones. Se han elegido los ejemplos y ejercicios de esta sección para los que existen soluciones de potencia. Sin embargo, no siempre existen soluciones de potencias. Aquellos que estén interesados en un tratamiento más riguroso de este tema deberán consultar un texto de ecuaciones diferenciales.

Estrategia de resolución de problemas

Estrategia para la resolución de problemas: Buscar soluciones en serie de potencia para ecuaciones diferenciales

  1. Suponga que la ecuación diferencial tiene una solución de la forma y(x)=∑n=0∞anxn.y(x)=∑n=0∞anxn.
  2. Diferencie la serie de potencias término a término para obtener y′(x)=∑n=1∞nanxn–1y′(x)=∑n=1∞nanxn–1 y y″(x)=∑n=2 ∞n(n–1)anxn–2 .y″(x)=∑n=2 ∞n(n–1)anxn–2 .
  3. Sustituya las expresiones de la serie de potencias en la ecuación diferencial.
  4. Reajuste las sumas según sea necesario para combinar los términos y simplificar la expresión.
  5. Iguale los coeficientes de las potencias similares de xx para determinar los valores de los coeficientes anan en la serie de potencias.
  6. Sustituya los coeficientes en la serie de potencias y escriba la solución.

Ejemplo 7.25

Soluciones en serie para las ecuaciones diferenciales

Halle una solución mediante series de potencias para las siguientes ecuaciones diferenciales.

  1. y″−y=0y″−y=0
  2. (x2 –1)y″+6xy′+4y=–4(x2 –1)y″+6xy′+4y=–4

Solución

  1. Supongamos que y(x)=∑n=0∞anxny(x)=∑n=0∞anxn (paso 1). Entonces, y′(x)=∑n=1∞nanxn–1y′(x)=∑n=1∞nanxn–1 y y″(x)=∑n=2 ∞n(n–1)anxn–2 y″(x)=∑n=2 ∞n(n–1)anxn–2 (paso 2). Queremos hallar valores para los coeficientes anan de manera que
    y″−y=0∑n=2 ∞n(n–1)anxn–2 –∑n=0∞anxn=0(paso 3).y″−y=0∑n=2 ∞n(n–1)anxn–2 –∑n=0∞anxn=0(paso 3).

    Queremos que los índices de nuestras sumas coincidan para poder expresarlos mediante una única suma. Es decir, queremos reescribir la primera suma para que empiece por n=0.n=0.
    Para volver a indexar el primer término, sustituya n por n+2 n+2 dentro de la suma y cambie el límite inferior de la suma por n=0.n=0. Obtenemos
    ∑n=2 ∞n(n–1)anxn–2 =∑n=0∞(n+2 )(n+1)an+2 xn.∑n=2 ∞n(n–1)anxn–2 =∑n=0∞(n+2 )(n+1)an+2 xn.

    Esto da
    ∑n=0∞(n+2 )(n+1)an+2 xn−∑n=0∞anxn=0∑n=0∞[(n+2 )(n+1)an+2 −an]xn=0(paso 4).∑n=0∞(n+2 )(n+1)an+2 xn−∑n=0∞anxn=0∑n=0∞[(n+2 )(n+1)an+2 −an]xn=0(paso 4).

    Dado que las expansiones en serie de potencias de las funciones son únicas, esta ecuación solo puede ser cierta si los coeficientes de cada potencia de x son cero. Así que tenemos
    (n+2 )(n+1)an+2 −an=0paran=0,1,2 ,….(n+2 )(n+1)an+2 −an=0paran=0,1,2 ,….

    Esta relación de recurrencia nos permite expresar cada coeficiente anan en términos del coeficiente de dos términos anteriores. Se obtiene una expresión para valores pares de n y otra para valores impares de n. Observando primero las ecuaciones que implican valores pares de n, vemos que
    a2 =a02 a4=a2 4⋅3=a04!a6=a46⋅5=a06!⋮.a2 =a02 a4=a2 4⋅3=a04!a6=a46⋅5=a06!⋮.

    Así, en general, cuando n es par, an=a0n!an=a0n! (paso 5).
    Para las ecuaciones que implican valores impares de n, vemos que
    a3=a13⋅2 =a13!a5=a35⋅4=a15!a7=a57⋅6=a17!⋮.a3=a13⋅2 =a13!a5=a35⋅4=a15!a7=a57⋅6=a17!⋮.

    Por lo tanto, en general, cuando n es impar, an=a1n!an=a1n! (continuación del paso 5).
    Uniendo todo esto, tenemos
    y(x)=∑n=0∞anxn=a0+a1x+a02 x2 +a13!x3+a04!x4+a15!x5+⋯.y(x)=∑n=0∞anxn=a0+a1x+a02 x2 +a13!x3+a04!x4+a15!x5+⋯.

    Volviendo a indexar las sumas para tener en cuenta los valores pares e impares de n por separado, obtenemos
    y(x)=a0∑k=0∞1(2 k)!x2 k+a1∑k=0∞1(2 k+1)!x2 k+1(paso 6).y(x)=a0∑k=0∞1(2 k)!x2 k+a1∑k=0∞1(2 k+1)!x2 k+1(paso 6).

    Análisis para la parte a.
    Como es de esperar para una ecuación diferencial de segundo orden, esta solución depende de dos constantes arbitrarias. Sin embargo, observe que nuestra ecuación diferencial es una ecuación diferencial de coeficiente constante, pero la solución de la serie de potencias no parece tener la forma común (que contiene funciones exponenciales) que estamos acostumbrados a ver. Además, como y(x)=c1ex+c2 e–xy(x)=c1ex+c2 e–x es la solución general de esta ecuación, debemos ser capaces de escribir cualquier solución en esta forma, y no está claro si la solución de la serie de potencias que acabamos de hallar puede, de hecho, escribirse en esa forma.
    Afortunadamente, después de escribir las representaciones en serie de potencias de exex y e–x,e–x, y haciendo un poco de álgebra, encontramos que si elegimos
    c0=(a0+a1)2 , c1=(a0−a1)2 ,c0=(a0+a1)2 , c1=(a0−a1)2 ,

    tenemos entonces a0=c0+c1a0=c0+c1 ya1=c0−c1,a1=c0−c1, y
    y(x)=a0+a1x+a02 x2 +a13!x3+a04!x4+a15!x5+⋯=(c0+c1)+(c0−c1)x+(c0+c1)2 x2 +(c0−c1)3!x3+(c0+c1)4!x4+(c0−c1)5!x5+⋯=c0∑n=0∞xnn!+c1∑n=0∞(−x)nn!=c0ex+c1e–x.y(x)=a0+a1x+a02 x2 +a13!x3+a04!x4+a15!x5+⋯=(c0+c1)+(c0−c1)x+(c0+c1)2 x2 +(c0−c1)3!x3+(c0+c1)4!x4+(c0−c1)5!x5+⋯=c0∑n=0∞xnn!+c1∑n=0∞(−x)nn!=c0ex+c1e–x.

    Así que, de hecho, encontramos la misma solución general. Observe que esta elección de c1c1 y c2 c2 no es evidente. Este es un caso en el que sabemos cuál debería ser la respuesta, e hicimos una "ingeniería inversa" de nuestra elección de coeficientes.
  2. Supongamos que y(x)=∑n=0∞anxny(x)=∑n=0∞anxn (paso 1). Entonces, y′(x)=∑n=1∞nanxn–1y′(x)=∑n=1∞nanxn–1 y y″(x)=∑n=2 ∞n(n–1)anxn–2 y″(x)=∑n=2 ∞n(n–1)anxn–2 (paso 2). Queremos hallar valores para los coeficientes anan de manera que
    (x2 –1)y″+6xy′+4y=−4(x2 –1)∑n=2 ∞n(n–1)anxn–2 +6x∑n=1∞nanxn–1+4∑n=0∞anxn=−4x2 ∑n=2 ∞n(n–1)anxn–2 –∑n=2 ∞n(n–1)anxn–2 +6x∑n=1∞nanxn–1+4∑n=0∞anxn=−4.(x2 –1)y″+6xy′+4y=−4(x2 –1)∑n=2 ∞n(n–1)anxn–2 +6x∑n=1∞nanxn–1+4∑n=0∞anxn=−4x2 ∑n=2 ∞n(n–1)anxn–2 –∑n=2 ∞n(n–1)anxn–2 +6x∑n=1∞nanxn–1+4∑n=0∞anxn=−4.

    Tomando los factores externos dentro de los sumandos, obtenemos
    ∑n=2 ∞n(n–1)anxn−∑n=2 ∞n(n–1)anxn–2 +∑n=1∞6nanxn+∑n=0∞4anxn=−4(paso 3).∑n=2 ∞n(n–1)anxn−∑n=2 ∞n(n–1)anxn–2 +∑n=1∞6nanxn+∑n=0∞4anxn=−4(paso 3).

    Ahora, en la primera suma, vemos que cuando n=0n=0 o n=1,n=1, el término se evalúa a cero, por lo que podemos añadir estos términos a nuestra suma para obtener
    ∑n=2 ∞n(n–1)anxn=∑n=0∞n(n–1)anxn.∑n=2 ∞n(n–1)anxn=∑n=0∞n(n–1)anxn.

    Del mismo modo, en el tercer término, vemos que cuando n=0,n=0, la expresión se evalúa a cero, por lo que también podemos añadir ese término. Tenemos
    ∑n=1∞6nanxn=∑n=0∞6nanxn.∑n=1∞6nanxn=∑n=0∞6nanxn.

    Entonces, solo tenemos que desplazar los índices en nuestro segundo término. Obtenemos
    ∑n=2 ∞n(n–1)anxn–2 =∑n=0∞(n+2 )(n+1)an+2 xn.∑n=2 ∞n(n–1)anxn–2 =∑n=0∞(n+2 )(n+1)an+2 xn.

    Así, tenemos
    ∑n=0∞n(n–1)anxn−∑n=0∞(n+2 )(n+1)an+2 xn+∑n=0∞6nanxn+∑n=0∞4anxn=−4(paso 4).∑n=0∞[n(n–1)an−(n+2 )(n+1)an+2 +6nan+4an]xn=−4∑n=0∞[(n2 −n)an+6nan+4an−(n+2 )(n+1)an+2 ]xn=−4∑n=0∞[n2 an+5nan+4an−(n+2 )(n+1)an+2 ]xn=−4∑n=0∞[(n2 +5n+4)an−(n+2 )(n+1)an+2 ]xn=−4∑n=0∞[(n+4)(n+1)an−(n+2 )(n+1)an+2 ]xn=–4∑n=0∞n(n–1)anxn−∑n=0∞(n+2 )(n+1)an+2 xn+∑n=0∞6nanxn+∑n=0∞4anxn=−4(paso 4).∑n=0∞[n(n–1)an−(n+2 )(n+1)an+2 +6nan+4an]xn=−4∑n=0∞[(n2 −n)an+6nan+4an−(n+2 )(n+1)an+2 ]xn=−4∑n=0∞[n2 an+5nan+4an−(n+2 )(n+1)an+2 ]xn=−4∑n=0∞[(n2 +5n+4)an−(n+2 )(n+1)an+2 ]xn=−4∑n=0∞[(n+4)(n+1)an−(n+2 )(n+1)an+2 ]xn=–4

    Observando los coeficientes de cada potencia de x, vemos que el término constante debe ser igual a −4,−4, y los coeficientes de todas las demás potencias de x deben ser cero. Entonces, mirando primero el término constante,
    4a0−2 a2 =−4a2 =2 a0+2 (paso 3).4a0−2 a2 =−4a2 =2 a0+2 (paso 3).

    Para n≥1,n≥1, tenemos
    (n+4)(n+1)an−(n+2 )(n+1)an+2 =0(n+1)[(n+4)an−(n+2 )an+2 ]=0.(n+4)(n+1)an−(n+2 )(n+1)an+2 =0(n+1)[(n+4)an−(n+2 )an+2 ]=0.

    Dado que n≥1,n≥1, n+1≠0,n+1≠0, vemos que
    (n+4)an−(n+2 )an+2 =0(n+4)an−(n+2 )an+2 =0

    y, por lo tanto,
    an+2 =n+4n+2 an.an+2 =n+4n+2 an.

    Para valores pares de n, tenemos
    a4=64(2 a0+2 )=3a0+3a6=86(3a0+3)=4a0+4⋮.a4=64(2 a0+2 )=3a0+3a6=86(3a0+3)=4a0+4⋮.

    En general, a2 k=(k+1)(a0+1)a2 k=(k+1)(a0+1) (paso 5).
    Para valores impares de n, tenemos
    a3=53a1a5=75a3=73a1a7=97a5=93a1=3a1⋮.a3=53a1a5=75a3=73a1a7=97a5=93a1=3a1⋮.

    En general, a2 k+1=2 k+33a1a2 k+1=2 k+33a1 (continuación del paso 5).
    Uniendo todo esto, tenemos
    y(x)=∑k=0∞(k+1)(a0+1)x2 k+∑k=0∞(2 k+33)a1x2 k+1(paso 6).y(x)=∑k=0∞(k+1)(a0+1)x2 k+∑k=0∞(2 k+33)a1x2 k+1(paso 6).

Punto de control 7.22

Halle una solución mediante series de potencias para las siguientes ecuaciones diferenciales.

  1. y′+2 xy=0y′+2 xy=0
  2. (x+1)y′=3y(x+1)y′=3y

Cerramos esta sección con una breve introducción a las funciones de Bessel. El tratamiento completo de las funciones de Bessel va mucho más allá del alcance de este curso, pero aquí tenemos una pequeña muestra del tema para que podamos ver cómo se utilizan las soluciones en serie de las ecuaciones diferenciales en aplicaciones del mundo real. La ecuación de Bessel de orden n viene dada por

x2 y″+xy′+(x2 −n2 )y=0.x2 y″+xy′+(x2 −n2 )y=0.

Esta ecuación surge en muchas aplicaciones físicas, en particular las que implican coordenadas cilíndricas, como la vibración de una cabeza de tambor circular y el calentamiento o enfriamiento transitorio de un cilindro. En el siguiente ejemplo, encontramos una solución en serie de potencias para la ecuación de Bessel de orden 0.

Ejemplo 7.26

Solución de la serie de potencias de la ecuación de Bessel

Halle una solución en serie de potencias de la ecuación de Bessel de orden 0 y grafique la solución.

Solución

La ecuación de Bessel de orden 0 viene dada por

x 2 y ″ + x y ′ + x 2 y = 0 . x 2 y ″ + x y ′ + x 2 y = 0 .

Suponemos una solución de la forma y=∑n=0∞anxn.y=∑n=0∞anxn. Entonces y′(x)=∑n=1∞nanxn–1y′(x)=∑n=1∞nanxn–1 y y″(x)=∑n=2 ∞n(n–1)anxn–2 .y″(x)=∑n=2 ∞n(n–1)anxn–2 . Sustituyendo esto en la ecuación diferencial, obtenemos

x 2 ∑ n = 2 ∞ n ( n – 1 ) a n x n – 2 + x ∑ n = 1 ∞ n a n x n – 1 + x 2 ∑ n = 0 ∞ a n x n = 0 Sustitución. ∑ n = 2 ∞ n ( n – 1 ) a n x n + ∑ n = 1 ∞ n a n x n + ∑ n = 0 ∞ a n x n + 2 = 0 Lleve los factores externos dentro de las sumas. ∑ n = 2 ∞ n ( n – 1 ) a n x n + ∑ n = 1 ∞ n a n x n + ∑ n = 2 ∞ a n – 2 x n = 0 Vuelva a indexar la tercera suma. ∑ n = 2 ∞ n ( n – 1 ) a n x n + a 1 x + ∑ n = 2 ∞ n a n x n + ∑ n = 2 ∞ a n – 2 x n = 0 Separe el término n = 1 de la segunda suma. a 1 x + ∑ n = 2 ∞ [ n ( n – 1 ) a n + n a n + a n – 2 ] x n = 0 Recoja los términos de la suma. a 1 x + ∑ n = 2 ∞ [ ( n 2 − n ) a n + n a n + a n – 2 ] x n = 0 Multiplique por en el primer término. a 1 x + ∑ n = 2 ∞ [ n 2 a n + a n – 2 ] x n = 0, Simplifique. x 2 ∑ n = 2 ∞ n ( n – 1 ) a n x n – 2 + x ∑ n = 1 ∞ n a n x n – 1 + x 2 ∑ n = 0 ∞ a n x n = 0 Sustitución. ∑ n = 2 ∞ n ( n – 1 ) a n x n + ∑ n = 1 ∞ n a n x n + ∑ n = 0 ∞ a n x n + 2 = 0 Lleve los factores externos dentro de las sumas. ∑ n = 2 ∞ n ( n – 1 ) a n x n + ∑ n = 1 ∞ n a n x n + ∑ n = 2 ∞ a n – 2 x n = 0 Vuelva a indexar la tercera suma. ∑ n = 2 ∞ n ( n – 1 ) a n x n + a 1 x + ∑ n = 2 ∞ n a n x n + ∑ n = 2 ∞ a n – 2 x n = 0 Separe el término n = 1 de la segunda suma. a 1 x + ∑ n = 2 ∞ [ n ( n – 1 ) a n + n a n + a n – 2 ] x n = 0 Recoja los términos de la suma. a 1 x + ∑ n = 2 ∞ [ ( n 2 − n ) a n + n a n + a n – 2 ] x n = 0 Multiplique por en el primer término. a 1 x + ∑ n = 2 ∞ [ n 2 a n + a n – 2 ] x n = 0, Simplifique.

Entonces, a1=0,a1=0, y para n≥2 ,n≥2 ,

n 2 a n + a n – 2 = 0 a n = − 1 n 2 a n – 2 . n 2 a n + a n – 2 = 0 a n = − 1 n 2 a n – 2 .

Dado que a1=0,a1=0, todos los términos impares son cero. Entonces, para valores pares de n, tenemos

a 2 = − 1 2 2 a 0 a 4 = − 1 4 2 a 2 = 1 4 2 ⋅ 2 2 a 0 . a 6 = − 1 6 2 a 4 = − 1 6 2 ⋅ 4 2 ⋅ 2 2 a 0 a 2 = − 1 2 2 a 0 a 4 = − 1 4 2 a 2 = 1 4 2 ⋅ 2 2 a 0 . a 6 = − 1 6 2 a 4 = − 1 6 2 ⋅ 4 2 ⋅ 2 2 a 0

En general,

a 2 k = ( –1 ) k ( 2 ) 2 k ( k ! ) 2 a 0 . a 2 k = ( –1 ) k ( 2 ) 2 k ( k ! ) 2 a 0 .

Por lo tanto, tenemos

y ( x ) = a 0 ∑ k = 0 ∞ ( –1 ) k ( 2 ) 2 k ( k ! ) 2 x 2 k . y ( x ) = a 0 ∑ k = 0 ∞ ( –1 ) k ( 2 ) 2 k ( k ! ) 2 x 2 k .
(7.10)

El gráfico aparece a continuación.

Esta figura es el gráfico de una función. El gráfico oscila con la mayor amplitud por encima del origen. El eje horizontal está marcado en incrementos de 2,5. El eje vertical está marcado en incrementos de 0,2.

Punto de control 7.23

Compruebe que la expresión encontrada en el Ejemplo 7.26 es una solución de la ecuación de Bessel de orden 0.

Sección 7.4 ejercicios

Halle una solución mediante series de potencias para las siguientes ecuaciones diferenciales.

104.

y ″ + 6 y ′ = 0 y ″ + 6 y ′ = 0

105.

5 y ″ + y ′ = 0 5 y ″ + y ′ = 0

106.

y ″ + 25 y = 0 y ″ + 25 y = 0

107.

y ″ − y = 0 y ″ − y = 0

108.

2 y ′ + y = 0 2 y ′ + y = 0

109.

y ′ − 2 x y = 0 y ′ − 2 x y = 0

110.

( x − 7 ) y ′ + 2 y = 0 ( x − 7 ) y ′ + 2 y = 0

111.

y ″ − x y ′ − y = 0 y ″ − x y ′ − y = 0

112.

( 1 + x 2 ) y ″ − 4 x y ′ + 6 y = 0 ( 1 + x 2 ) y ″ − 4 x y ′ + 6 y = 0

113.

x 2 y ″ − x y ′ − 3 y = 0 x 2 y ″ − x y ′ − 3 y = 0

114.

y ″ − 8 y ′ = 0 , y ( 0 ) = –2 , y ′ ( 0 ) = 10 y ″ − 8 y ′ = 0 , y ( 0 ) = –2 , y ′ ( 0 ) = 10

115.

y ″ − 2 x y = 0 , y ( 0 ) = 1 , y ′ ( 0 ) = −3 y ″ − 2 x y = 0 , y ( 0 ) = 1 , y ′ ( 0 ) = −3

116.

La ecuación diferencial x2 y″+xy′+(x2 –1)y=0x2 y″+xy′+(x2 –1)y=0 es una ecuación de Bessel de orden 1. Utilice una serie de potencias de la forma y=∑n=0∞anxny=∑n=0∞anxn para hallar la solución.