Ecuaciones lineales no homogéneas

7.2 Ecuaciones lineales no homogéneas

Objetivos de aprendizaje

  • 7.2.1 Escribir la solución general de una ecuación diferencial no homogénea.
  • 7.2.2 Resolver una ecuación diferencial no homogénea por el método de los coeficientes indeterminados.
  • 7.2.3 Resolver una ecuación diferencial no homogénea por el método de variación de parámetros.

En esta sección, examinamos cómo resolver ecuaciones diferenciales no homogéneas. La terminología y los métodos son diferentes de los que utilizamos para las ecuaciones homogéneas, así que vamos a empezar por definir algunos términos nuevos.

Solución general para una ecuación lineal no homogénea

Consideremos la ecuación diferencial lineal no homogénea

a2 (x)y″+a1(x)y′+a0(x)y=r(x).a2 (x)y″+a1(x)y′+a0(x)y=r(x).

La ecuación homogénea asociada

a2 (x)y″+a1(x)y′+a0(x)y=0a2 (x)y″+a1(x)y′+a0(x)y=0
(7.3)

se denomina ecuación complementaria. Veremos que la resolución de la ecuación complementaria es un paso importante para resolver una ecuación diferencial no homogénea.

Definición

Una solución yp(x)yp(x) de una ecuación diferencial que no contiene constantes arbitrarias se llama solución particular de la ecuación.

Teorema 7.4

Solución general para una ecuación no homogénea

Supongamos que yp(x)yp(x) es cualquier solución particular de la ecuación diferencial no homogénea lineal

a2 (x)y″+a1(x)y′+a0(x)y=r(x).a2 (x)y″+a1(x)y′+a0(x)y=r(x).

Además, supongamos que c1y1(x)+c2 y2 (x)c1y1(x)+c2 y2 (x) denota la solución general de la ecuación complementaria. Entonces, la solución general de la ecuación no homogénea viene dada por

y(x)=c1y1(x)+c2 y2 (x)+yp(x).y(x)=c1y1(x)+c2 y2 (x)+yp(x).
(7.4)

Prueba

Para demostrar que y(x)y(x) es la solución general, primero debemos demostrar que resuelve la ecuación diferencial y, segundo, que cualquier solución de la ecuación diferencial puede escribirse en esa forma. Al sustituir y(x)y(x) en la ecuación diferencial, tenemos

a2 (x)y″+a1(x)y′+a0(x)y=a2 (x)(c1y1+c2 y2 +yp)″+a1(x)(c1y1+c2 y2 +yp)′+a0(x)(c1y1+c2 y2 +yp)=[a2 (x)(c1y1+c2 y2 )″+a1(x)(c1y1+c2 y2 )′+a0(x)(c1y1+c2 y2 )]+a2 (x)yp″+a1(x)yp′+a0(x)yp=0+r(x)=r(x).a2 (x)y″+a1(x)y′+a0(x)y=a2 (x)(c1y1+c2 y2 +yp)″+a1(x)(c1y1+c2 y2 +yp)′+a0(x)(c1y1+c2 y2 +yp)=[a2 (x)(c1y1+c2 y2 )″+a1(x)(c1y1+c2 y2 )′+a0(x)(c1y1+c2 y2 )]+a2 (x)yp″+a1(x)yp′+a0(x)yp=0+r(x)=r(x).

Así que y(x)y(x) es una solución.

Ahora, supongamos que z(x)z(x) es cualquier solución a a2 (x)y″+a1(x)y′+a0(x)y=r(x).a2 (x)y″+a1(x)y′+a0(x)y=r(x). Entonces

a2 (x)(z−yp)″+a1(x)(z−yp)′+a0(x)(z−yp)=(a2 (x)z″+a1(x)z′+a0(x)z)−(a2 (x)yp″+a1(x)yp′+a0(x)yp)=r(x)−r(x)=0,a2 (x)(z−yp)″+a1(x)(z−yp)′+a0(x)(z−yp)=(a2 (x)z″+a1(x)z′+a0(x)z)−(a2 (x)yp″+a1(x)yp′+a0(x)yp)=r(x)−r(x)=0,

por lo que z(x)−yp(x)z(x)−yp(x) es una solución a la ecuación complementaria. Pero, c1y1(x)+c2 y2 (x)c1y1(x)+c2 y2 (x) es la solución general de la ecuación complementaria, por lo que hay constantes c1c1 y c2 c2 tal que

z(x)−yp(x)=c1y1(x)+c2 y2 (x).z(x)−yp(x)=c1y1(x)+c2 y2 (x).

Por lo tanto, vemos que z(x)=c1y1(x)+c2 y2 (x)+yp(x).z(x)=c1y1(x)+c2 y2 (x)+yp(x).

□

Ejemplo 7.11

Verificar la solución general

Dado que yp(x)=xyp(x)=x es una solución particular de la ecuación diferencial y″+y=x,y″+y=x, escriba la solución general y compruebe que la solución satisface la ecuación.

Solución

La ecuación complementaria es y″+y=0,y″+y=0, que tiene la solución general c1cosx+c2 senx.c1cosx+c2 senx. Entonces, la solución general de la ecuación no homogénea es

y ( x ) = c 1 cos x + c 2 sen x + x . y ( x ) = c 1 cos x + c 2 sen x + x .

Para comprobar que esta es una solución, sustitúyala en la ecuación diferencial. Tenemos

y ′ ( x ) = − c 1 sen x + c 2 cos x + 1 y y ″ ( x ) = − c 1 cos x − c 2 sen x . y ′ ( x ) = − c 1 sen x + c 2 cos x + 1 y y ″ ( x ) = − c 1 cos x − c 2 sen x .

Entonces

y ″ ( x ) + y ( x ) = − c 1 cos x − c 2 sen x + c 1 cos x + c 2 sen x + x = x . y ″ ( x ) + y ( x ) = − c 1 cos x − c 2 sen x + c 1 cos x + c 2 sen x + x = x .

Así que, y(x)y(x) es una solución para y″+y=x.y″+y=x.

Punto de control 7.10

Dado que yp(x)=−2yp(x)=−2 es una solución particular para y″−3y′−4y=8,y″−3y′−4y=8, escriba la solución general y verifique que la solución general satisface la ecuación.

En el apartado anterior aprendimos a resolver ecuaciones homogéneas con coeficientes constantes. Por lo tanto, para las ecuaciones no homogéneas de la forma ay″+by′+cy=r(x),ay″+by′+cy=r(x), ya sabemos cómo resolver la ecuación complementaria, y el problema se reduce a hallar una solución particular para la ecuación no homogénea. A continuación examinamos dos técnicas para ello: el método de los coeficientes indeterminados y el método de la variación de los parámetros.

Coeficientes indeterminados

El método de los coeficientes indeterminados consiste en hacer conjeturas sobre la forma de la solución particular basándose en la forma de r(x).r(x). Cuando tomamos derivadas de polinomios, funciones exponenciales, senos y cosenos, obtenemos polinomios, funciones exponenciales, senos y cosenos. Así que cuando r(x)r(x) tiene una de estas formas, es posible que la solución de la ecuación diferencial no homogénea tenga esa misma forma. Veamos algunos ejemplos para ver cómo funciona.

Ejemplo 7.12

Coeficientes indeterminados cuando r(x)r(x) es un polinomio

Halle la solución general de y″+4y′+3y=3x.y″+4y′+3y=3x.

Solución

La ecuación complementaria es y″+4y′+3y=0,y″+4y′+3y=0, con la solución general c1e–x+c2 e−3x.c1e–x+c2 e−3x. Dado que r(x)=3x,r(x)=3x, la solución particular podría tener la forma yp(x)=Ax+B.yp(x)=Ax+B. Si este es el caso, entonces tenemos yp′(x)=Ayp′(x)=A y yp″(x)=0.yp″(x)=0. Para ypyp sea una solución de la ecuación diferencial, debemos encontrar valores para AA y BB de manera que

y ″ + 4 y ′ + 3 y = 3 x 0 + 4 ( A ) + 3 ( A x + B ) = 3 x 3 A x + ( 4 A + 3 B ) = 3 x . y ″ + 4 y ′ + 3 y = 3 x 0 + 4 ( A ) + 3 ( A x + B ) = 3 x 3 A x + ( 4 A + 3 B ) = 3 x .

Al igualar los coeficientes de los términos similares, tenemos

3 A = 3 4 A + 3 B = 0 . 3 A = 3 4 A + 3 B = 0 .

Entonces, A=1A=1 y B=–43,B=–43, así que yp(x)=x−43yp(x)=x−43 y la solución general es

y ( x ) = c 1 e – x + c 2 e −3 x + x − 4 3 . y ( x ) = c 1 e – x + c 2 e −3 x + x − 4 3 .

En el Ejemplo 7.12, observe que aunque r(x)r(x) no incluía un término constante, era necesario que incluyéramos el término constante en nuestra conjetura. Si hubiéramos asumido una solución de la forma yp=Axyp=Ax (sin término constante), no habríamos podido hallar una solución. (¡Verifique esto!). Si se grafica la función r(x)r(x) es un polinomio, nuestra conjetura para la solución particular debe ser un polinomio del mismo grado, y debe incluir todos los términos de orden inferior, independientemente de que estén presentes en r(x).r(x).

Ejemplo 7.13

Coeficientes indeterminados cuando r(x)r(x) es un exponencial

Halle la solución general de y″−y′−2 y=2 e3x.y″−y′−2 y=2 e3x.

Solución

La ecuación complementaria es y″−y′−2 y=0,y″−y′−2 y=0, con la solución general c1e–x+c2 e2 x.c1e–x+c2 e2 x. Dado que r(x)=2 e3x,r(x)=2 e3x, la solución particular podría tener la forma yp(x)=Ae3x.yp(x)=Ae3x. Entonces, tenemos yp′(x)=3Ae3xyp′(x)=3Ae3x como yp″(x)=9Ae3x.yp″(x)=9Ae3x. Para ypyp sea una solución de la ecuación diferencial, debemos encontrar un valor para AA tal que

y ″ − y ′ − 2 y = 2 e 3 x 9 A e 3 x − 3 A e 3 x − 2 A e 3 x = 2 e 3 x 4 A e 3 x = 2 e 3 x . y ″ − y ′ − 2 y = 2 e 3 x 9 A e 3 x − 3 A e 3 x − 2 A e 3 x = 2 e 3 x 4 A e 3 x = 2 e 3 x .

Así que, 4A=2 4A=2 y A=1/2 .A=1/2 . Entonces, yp(x)=(12 )e3x,yp(x)=(12 )e3x, y la solución general es

y ( x ) = c 1 e – x + c 2 e 2 x + 1 2 e 3 x . y ( x ) = c 1 e – x + c 2 e 2 x + 1 2 e 3 x .

Punto de control 7.11

Halle la solución general de y″−4y′+4y=7sent−cost.y″−4y′+4y=7sent−cost.

En el punto de control anterior, r(x)r(x) incluyó los términos del seno y del coseno. Sin embargo, aunque r(x)r(x) incluyó un término de seno solamente o un término de coseno solamente, ambos términos deben estar presentes en la conjetura. El método de los coeficientes indeterminados también funciona con productos de polinomios, exponenciales, senos y cosenos. Algunas de las principales formas de r(x)r(x) y las conjeturas asociadas para yp(x)yp(x) se resumen en la Tabla 7.2.

r(x)r(x) Conjetura inicial para yp(x)yp(x) grandes.
kk (una constante) AA (una constante)
ax+bax+b Ax+BAx+B (Nota: La conjetura debe incluir ambos términos aunque b=0.b=0.)
ax2 +bx+cax2 +bx+c Ax2 +Bx+CAx2 +Bx+C (Nota: La conjetura debe incluir los tres términos aunque bb o cc son cero).
Polinomios de orden superior Polinomio del mismo orden que r(x)r(x) grandes.
aeλxaeλx AeλxAeλx
acosβx+bsenβxacosβx+bsenβx Acosβx+BsenβxAcosβx+Bsenβx (Nota: La conjetura debe incluir ambos términos, incluso si cualquiera de ellos sean a=0a=0 o b=0.b=0.)
aeαxcosβx+beαxsenβxaeαxcosβx+beαxsenβx Aeαxcosβx+BeαxsenβxAeαxcosβx+Beαxsenβx
(ax2 +bx+c)eλx(ax2 +bx+c)eλx (Ax2 +Bx+C)eλx(Ax2 +Bx+C)eλx
(a2 x2 +a1x+a0)cosβx+(b2 x2 +b1x+b0)senβx(a2 x2 +a1x+a0)cosβx+(b2 x2 +b1x+b0)senβx (A2 x2 +A1x+A0)cosβx+(B2 x2 +B1x+B0)senβx(A2 x2 +A1x+A0)cosβx+(B2 x2 +B1x+B0)senβx
(a2 x2 +a1x+a0)eαxcosβx+(b2 x2 +b1x+b0)eαxsenβx(a2 x2 +a1x+a0)eαxcosβx+(b2 x2 +b1x+b0)eαxsenβx (A2 x2 +A1x+A0)eαxcosβx+(B2 x2 +B1x+B0)eαxsenβx(A2 x2 +A1x+A0)eαxcosβx+(B2 x2 +B1x+B0)eαxsenβx
Tabla 7.2 Formas clave del método de los coeficientes indeterminados

Hay que tener en cuenta que este método tiene un inconveniente importante. Consideremos la ecuación diferencial y″+5y′+6y=3e−2x.y″+5y′+6y=3e−2x. Con base en la forma de r(x),r(x), adivinamos una solución particular de la forma yp(x)=Ae−2x.yp(x)=Ae−2x. Pero cuando sustituimos esta expresión en la ecuación diferencial para encontrar un valor para A,A, nos encontramos con un problema. Tenemos

yp′(x)=−2Ae−2xyp′(x)=−2Ae−2x

y

yp″=4Ae−2x,yp″=4Ae−2x,

por lo que queremos

y″+5y′+6y=3e−2x4Ae−2x+5(−2Ae−2x)+6Ae−2x=3e−2x4Ae−2x−10Ae−2x+6Ae−2x=3e−2x0=3e−2x,y″+5y′+6y=3e−2x4Ae−2x+5(−2Ae−2x)+6Ae−2x=3e−2x4Ae−2x−10Ae−2x+6Ae−2x=3e−2x0=3e−2x,

que no es posible.

Si observamos detenidamente, vemos que, en este caso, la solución general de la ecuación complementaria es c1e−2x+c2 e−3x.c1e−2x+c2 e−3x. La función exponencial en r(x)r(x) es en realidad una solución de la ecuación complementaria, por lo que, como acabamos de ver, todos los términos del lado izquierdo de la ecuación se cancelan. En este caso podemos seguir utilizando el método de los coeficientes indeterminados, pero tenemos que modificar nuestra conjetura multiplicándola porx.x. Usando la nueva conjetura, yp(x)=Axe−2x,yp(x)=Axe−2x, tenemos

yp′(x)=A(e−2x−2 xe−2x)yp′(x)=A(e−2x−2 xe−2x)

y

yp″(x)=−4Ae−2x+4Axe−2x.yp″(x)=−4Ae−2x+4Axe−2x.

La sustitución da como resultado

y″+5y′+6y=3e−2x(−4Ae−2x+4Axe−2x)+5(Ae−2x−2 Axe−2x)+6Axe−2x=3e−2x−4Ae−2x+4Axe−2x+5Ae−2x−10Axe−2x+6Axe−2x=3e−2xAe−2x=3e−2x.y″+5y′+6y=3e−2x(−4Ae−2x+4Axe−2x)+5(Ae−2x−2 Axe−2x)+6Axe−2x=3e−2x−4Ae−2x+4Axe−2x+5Ae−2x−10Axe−2x+6Axe−2x=3e−2xAe−2x=3e−2x.

Así que, A=3A=3 y yp(x)=3xe−2x.yp(x)=3xe−2x. Esto nos da la siguiente solución general

y(x)=c1e−2x+c2 e−3x+3xe−2x.y(x)=c1e−2x+c2 e−3x+3xe−2x.

Observe que si xe−2xxe−2x fuera también una solución de la ecuación complementaria, tendríamos que multiplicar por xx de nuevo, y trataríamos yp(x)=Ax2 e−2x.yp(x)=Ax2 e−2x.

Estrategia de resolución de problemas

Estrategia para la resolución de problemas: Método de los coeficientes indeterminados

  1. Resuelva la ecuación complementaria y escriba la solución general.
  2. Con base en la forma de r(x),r(x), haga una conjetura inicial para yp(x).yp(x).
  3. Compruebe si algún término de la conjetura de yp(x)yp(x) es una solución a la ecuación complementaria. Si es así, multiplique la estimación porx.x. Repita este paso hasta que no haya términos en yp(x)yp(x) que resuelven la ecuación complementaria.
  4. Sustituya yp(x)yp(x) en la ecuación diferencial e iguale los términos similares para encontrar los valores de los coeficientes desconocidos en yp(x).yp(x).
  5. Sume la solución general de la ecuación complementaria y la solución particular que acaba de encontrar para obtener la solución general de la ecuación no homogénea.

Ejemplo 7.14

Resolver ecuaciones no homogéneas

Halle las soluciones generales de las siguientes ecuaciones diferenciales.

  1. y″−9y=−6cos3xy″−9y=−6cos3x
  2. x″+2 x′+x=4e−tx″+2 x′+x=4e−t
  3. y″−2 y′+5y=10x2 −3x−3y″−2 y′+5y=10x2 −3x−3
  4. y″−3y′=−12ty″−3y′=−12t

Solución

  1. La ecuación complementaria es y″−9y=0,y″−9y=0, que tiene la solución general c1e3x+c2 e−3xc1e3x+c2 e−3x (paso 1). Con base en la forma de r(x)=−6cos3x,r(x)=−6cos3x, nuestra suposición inicial para la solución particular es yp(x)=Acos3x+Bsen3xyp(x)=Acos3x+Bsen3x (paso 2). Ninguno de los términos de yp(x)yp(x) resuelven la ecuación complementaria, por lo que es una suposición válida (paso 3).
    Ahora queremos encontrar valores para AA y B,B, así que sustituimos ypyp en la ecuación diferencial. Tenemos
    yp′(x)=−3Asen3x+3Bcos3xyyp″(x)=−9Acos3x−9Bsen3x,yp′(x)=−3Asen3x+3Bcos3xyyp″(x)=−9Acos3x−9Bsen3x,

    por lo que queremos encontrar valores de AA y BB de manera que
    y″−9y=−6cos3x−9Acos3x−9Bsen3x−9(Acos3x+Bsen3x)=−6cos3x−18Acos3x−18Bsen3x=−6cos3x.y″−9y=−6cos3x−9Acos3x−9Bsen3x−9(Acos3x+Bsen3x)=−6cos3x−18Acos3x−18Bsen3x=−6cos3x.

    Por lo tanto,
    −18A=−6−18B=0.−18A=−6−18B=0.

    Esto da A=13A=13 y B=0,B=0, así que yp(x)=(13)cos3xyp(x)=(13)cos3x (paso 4).
    Poniendo todo junto, tenemos la solución general
    y(x)=c1e3x+c2 e−3x+13cos3x.y(x)=c1e3x+c2 e−3x+13cos3x.
  2. La ecuación complementaria es x″+2 x′+x=0,x″+2 x′+x=0, que tiene la solución general c1e−t+c2 te−tc1e−t+c2 te−t (paso 1). Con base en la forma r(t)=4e−t,r(t)=4e−t, nuestra suposición inicial para la solución particular es xp(t)=Ae−txp(t)=Ae−t (paso 2). Sin embargo, vemos que esta conjetura resuelve la ecuación complementaria, por lo que debemos multiplicar port,t, lo que da una nueva suposición xp(t)=Ate−txp(t)=Ate−t (paso 3). Comprobando esta nueva conjetura, vemos que también resuelve la ecuación complementaria, por lo que debemos volver a multiplicar por t, lo que da xp(t)=At2 e−txp(t)=At2 e−t (paso 3, de nuevo). Ahora, comprobando esta conjetura, vemos que xp(t)xp(t) no resuelve la ecuación complementaria, por lo que es una conjetura válida (paso 3, de nuevo).
    Ahora queremos hallar un valor para A,A, por lo que sustituimos xpxp en la ecuación diferencial. Tenemos
    xp(t)=At2 e−t,por lo quexp′(t)=2 Ate−t−At2 e−txp(t)=At2 e−t,por lo quexp′(t)=2 Ate−t−At2 e−t

    y xp″(t)=2 Ae−t−2 Ate−t−(2 Ate−t−At2 e−t)=2 Ae−t−4Ate−t+At2 e−t.xp″(t)=2 Ae−t−2 Ate−t−(2 Ate−t−At2 e−t)=2 Ae−t−4Ate−t+At2 e−t.
    Sustituyendo en la ecuación diferencial, queremos hallar un valor de AA para que
    x″+2 x′+x=4e−t2 Ae−t−4Ate−t+At2 e−t+2 (2 Ate−t−At2 e−t)+At2 e−t=4e−t2 Ae−t=4e−t.x″+2 x′+x=4e−t2 Ae−t−4Ate−t+At2 e−t+2 (2 Ate−t−At2 e−t)+At2 e−t=4e−t2 Ae−t=4e−t.

    Esto da A=2 ,A=2 , por lo que xp(t)=2 t2 e−txp(t)=2 t2 e−t (paso 4). Juntando todo, tenemos la solución general
    x(t)=c1e−t+c2 te−t+2 t2 e−t.x(t)=c1e−t+c2 te−t+2 t2 e−t.
  3. La ecuación complementaria es y″−2 y′+5y=0,y″−2 y′+5y=0, que tiene la solución general c1excos2 x+c2 exsen2 xc1excos2 x+c2 exsen2 x (paso 1). Con base en la forma r(x)=10x2 −3x−3,r(x)=10x2 −3x−3, nuestra suposición inicial para la solución particular es yp(x)=Ax2 +Bx+Cyp(x)=Ax2 +Bx+C (paso 2). Ninguno de los términos de yp(x)yp(x) resuelven la ecuación complementaria, por lo que es una suposición válida (paso 3). Ahora queremos hallar valores para A,A, B,B, y C,C, por lo que sustituimos ypyp en la ecuación diferencial. Tenemos yp′(x)=2 Ax+Byp′(x)=2 Ax+B y yp″(x)=2 A,yp″(x)=2 A, por lo que queremos hallar valores de A,A, B,B, y CC de manera que
    y″−2 y′+5y=10x2 −3x−32 A−2 (2 Ax+B)+5(Ax2 +Bx+C)=10x2 −3x−35Ax2 +(5B−4A)x+(5C−2 B+2 A)=10x2 −3x−3.y″−2 y′+5y=10x2 −3x−32 A−2 (2 Ax+B)+5(Ax2 +Bx+C)=10x2 −3x−35Ax2 +(5B−4A)x+(5C−2 B+2 A)=10x2 −3x−3.

    Por lo tanto,
    5A=105B−4A=−35C−2 B+2 A=−3.5A=105B−4A=−35C−2 B+2 A=−3.

    Esto da A=2 ,A=2 , B=1,B=1, y C=−1,C=−1, así que yp(x)=2 x2 +x–1yp(x)=2 x2 +x–1 (paso 4). Juntando todo, tenemos la solución general
    y(x)=c1excos2 x+c2 exsen2 x+2 x2 +x–1.y(x)=c1excos2 x+c2 exsen2 x+2 x2 +x–1.
  4. La ecuación complementaria es y″−3y′=0,y″−3y′=0, que tiene la solución general c1e3t+c2 c1e3t+c2 (paso 1). Con base en la forma r(t)=−12t,r(t)=−12t, nuestra suposición inicial para la solución particular es yp(t)=At+Byp(t)=At+B (paso 2). Sin embargo, vemos que el término constante de esta conjetura resuelve la ecuación complementaria, por lo que debemos multiplicar por t,t, lo que da una nueva suposición yp(t)=At2 +Btyp(t)=At2 +Bt (paso 3). Comprobando esta nueva conjetura, vemos que ninguno de los términos de yp(t)yp(t) resuelve la ecuación complementaria, por lo que es una conjetura válida (paso 3, de nuevo). Ahora queremos hallar valores para AA y B,B, por lo que sustituimos ypyp en la ecuación diferencial. Tenemos yp′(t)=2 At+Byp′(t)=2 At+B y yp″(t)=2 A,yp″(t)=2 A, por lo que queremos hallar valores de AA y BB de manera que
    y″−3y′=−12t2 A−3(2 At+B)=−12t−6At+(2 A−3B)=−12t.y″−3y′=−12t2 A−3(2 At+B)=−12t−6At+(2 A−3B)=−12t.

    Por lo tanto,
    −6A=−122 A−3B=0.−6A=−122 A−3B=0.

    Esto da A=2 A=2 y B=4/3,B=4/3, así que yp(t)=2 t2 +(4/3)typ(t)=2 t2 +(4/3)t (paso 4). Juntando todo, tenemos la solución general
    y(t)=c1e3t+c2 +2 t2 +43t.y(t)=c1e3t+c2 +2 t2 +43t.

Punto de control 7.12

Halle la solución general de las siguientes ecuaciones diferenciales.

  1. y″−5y′+4y=3exy″−5y′+4y=3ex
  2. y″+y′−6y=52cos2 ty″+y′−6y=52cos2 t

Variación de los parámetros

A veces, r(x)r(x) no es una combinación de polinomios, exponenciales o senos y cosenos. Cuando este es el caso, el método de los coeficientes indeterminados no funciona y tenemos que utilizar otro enfoque para hallar una solución particular a la ecuación diferencial. Utilizamos un enfoque denominado método de variación de los parámetros.

Para simplificar un poco nuestros cálculos, vamos a dividir la ecuación diferencial entre a,a, por lo que tenemos un coeficiente principal de 1. Entonces la ecuación diferencial tiene la forma

y″+py′+qy=r(x),y″+py′+qy=r(x),

donde pp y qq son constantes.

Si la solución general de la ecuación complementaria viene dada por c1y1(x)+c2 y2 (x),c1y1(x)+c2 y2 (x), vamos a buscar una solución particular de la forma yp(x)=u(x)y1(x)+v(x)y2 (x).yp(x)=u(x)y1(x)+v(x)y2 (x). En este caso, utilizamos las dos soluciones linealmente independientes de la ecuación complementaria para formar nuestra solución particular. Sin embargo, estamos asumiendo que los coeficientes son funciones de x, en vez de constantes. Queremos encontrar funciones u(x)u(x) y v(x)v(x) de manera que yp(x)yp(x) satisfagan la ecuación diferencial. Tenemos

yp=uy1+vy2 yp′=u′y1+uy1′+v′y2 +vy2 ′yp″=(u′y1+v′y2 )′+u′y1′+uy1″+v′y2 ′+vy2 ″.yp=uy1+vy2 yp′=u′y1+uy1′+v′y2 +vy2 ′yp″=(u′y1+v′y2 )′+u′y1′+uy1″+v′y2 ′+vy2 ″.

Sustituyendo en la ecuación diferencial, obtenemos

yp″+pyp′+qyp=[(u′y1+v′y2 )′+u′y1′+uy1″+v′y2 ′+vy2 ″]+p[u′y1+uy1′+v′y2 +vy2 ′]+q[uy1+vy2 ]=u[y1″+py1′+qy1]+v[y2 ″+py2 ′+qy2 ]+(u′y1+v′y2 )′+p(u′y1+v′y2 )+(u′y1′+v′y2 ′).yp″+pyp′+qyp=[(u′y1+v′y2 )′+u′y1′+uy1″+v′y2 ′+vy2 ″]+p[u′y1+uy1′+v′y2 +vy2 ′]+q[uy1+vy2 ]=u[y1″+py1′+qy1]+v[y2 ″+py2 ′+qy2 ]+(u′y1+v′y2 )′+p(u′y1+v′y2 )+(u′y1′+v′y2 ′).

Observe que y1y1 y y2 y2 son soluciones de la ecuación complementaria, por lo que los dos primeros términos son cero. Por lo tanto, tenemos

(u′y1+v′y2 )′+p(u′y1+v′y2 )+(u′y1′+v′y2 ′)=r(x).(u′y1+v′y2 )′+p(u′y1+v′y2 )+(u′y1′+v′y2 ′)=r(x).

Si simplificamos esta ecuación imponiendo la condición adicional u′y1+v′y2 =0,u′y1+v′y2 =0, los dos primeros términos son cero, y esto se reduce a u′y1′+v′y2 ′=r(x).u′y1′+v′y2 ′=r(x). Así, con esta condición adicional, tenemos un sistema de dos ecuaciones en dos incógnitas:

u′y1+v′y2 =0u′y1′+v′y2 ′=r(x).u′y1+v′y2 =0u′y1′+v′y2 ′=r(x).

La resolución de este sistema nos da u′u′ y v′,v′, que podemos integrar para hallar u y v.

Entonces, yp(x)=u(x)y1(x)+v(x)y2 (x)yp(x)=u(x)y1(x)+v(x)y2 (x) es una solución particular de la ecuación diferencial. Resolver este sistema de ecuaciones, a veces, es un reto, así que aprovechemos para repasar la regla de Cramer, que nos permite resolver el sistema de ecuaciones utilizando determinantes.

Regla: regla de Cramer

El sistema de ecuaciones

a1z1+b1z2 =r1a2 z1+b2 z2 =r2 a1z1+b1z2 =r1a2 z1+b2 z2 =r2

tiene una solución única si y solo si el determinante de los coeficientes es distinto de cero. En este caso, la solución viene dada por

z1=|r1b1r2 b2 ||a1b1a2 b2 |yz2 =|a1r1a2 r2 ||a1b1a2 b2 |.z1=|r1b1r2 b2 ||a1b1a2 b2 |yz2 =|a1r1a2 r2 ||a1b1a2 b2 |.

Ejemplo 7.15

Usar la regla de Cramer

Utilice la regla de Cramer para resolver el siguiente sistema de ecuaciones.

x2 z1+2 xz2 =0z1−3x2 z2 =2 xx2 z1+2 xz2 =0z1−3x2 z2 =2 x

Solución

Tenemos

a 1 ( x ) = x 2 a 2 ( x ) = 1 b 1 ( x ) = 2 x b 2 ( x ) = −3 x 2 r 1 ( x ) = 0 r 2 ( x ) = 2 x . a 1 ( x ) = x 2 a 2 ( x ) = 1 b 1 ( x ) = 2 x b 2 ( x ) = −3 x 2 r 1 ( x ) = 0 r 2 ( x ) = 2 x .

Entonces,

| a 1 b 1 a 2 b 2 | = | x 2 2 x 1 −3 x 2 | = −3 x 4 – 2 x | a 1 b 1 a 2 b 2 | = | x 2 2 x 1 −3 x 2 | = −3 x 4 – 2 x

y

| r 1 b 1 r 2 b 2 | = | 0 2 x 2 x −3 x 2 | = 0 − 4 x 2 = −4 x 2 . | r 1 b 1 r 2 b 2 | = | 0 2 x 2 x −3 x 2 | = 0 − 4 x 2 = −4 x 2 .

Por lo tanto,

z 1 = | r 1 b 1 r 2 b 2 | | a 1 b 1 a 2 b 2 | = −4 x 2 −3 x 4 – 2 x = 4 x 3 x 3 + 2 . z 1 = | r 1 b 1 r 2 b 2 | | a 1 b 1 a 2 b 2 | = −4 x 2 −3 x 4 – 2 x = 4 x 3 x 3 + 2 .

Además,

| a 1 r 1 a 2 r 2 | = | x 2 0 1 2 x | = 2 x 3 − 0 = 2 x 3 . | a 1 r 1 a 2 r 2 | = | x 2 0 1 2 x | = 2 x 3 − 0 = 2 x 3 .

Por lo tanto,

z 2 = | a 1 r 1 a 2 r 2 | | a 1 b 1 a 2 b 2 | = 2 x 3 −3 x 4 – 2 x = –2 x 2 3 x 3 + 2 . z 2 = | a 1 r 1 a 2 r 2 | | a 1 b 1 a 2 b 2 | = 2 x 3 −3 x 4 – 2 x = –2 x 2 3 x 3 + 2 .

Punto de control 7.13

Utilice la regla de Cramer para resolver el siguiente sistema de ecuaciones.

2 xz1−3z2 =0x2 z1+4xz2 =x+12 xz1−3z2 =0x2 z1+4xz2 =x+1

Estrategia de resolución de problemas

Estrategia para la resolución de problemas: Método de variación de los parámetros

  1. Resuelva la ecuación complementaria y escriba la solución general
    c1y1(x)+c2 y2 (x).c1y1(x)+c2 y2 (x).
  2. Utilice la regla de Cramer u otra técnica adecuada para hallar funciones u′(x)u′(x) y v′(x)v′(x) que satisfacen
    u′y1+v′y2 =0u′y1′+v′y2 ′=r(x).u′y1+v′y2 =0u′y1′+v′y2 ′=r(x).
  3. Integre u′u′ y v′v′ para calcular u(x)u(x) y v(x).v(x). Entonces, yp(x)=u(x)y1(x)+v(x)y2 (x)yp(x)=u(x)y1(x)+v(x)y2 (x) es una solución particular de la ecuación.
  4. Sume la solución general de la ecuación complementaria y la solución particular encontrada en el paso 3 para obtener la solución general de la ecuación no homogénea.

Ejemplo 7.16

Usar el método de variación de los parámetros

Halle la solución general de las siguientes ecuaciones diferenciales.

  1. y″−2 y′+y=ett2 y″−2 y′+y=ett2
  2. y″+y=3sen2 xy″+y=3sen2 x

Solución

  1. La ecuación complementaria es y″−2 y′+y=0y″−2 y′+y=0 con la solución general asociada c1et+c2 tet.c1et+c2 tet. Por lo tanto, y1(t)=ety1(t)=et y y2 (t)=tet.y2 (t)=tet. Calculando las derivadas, obtenemos y1′(t)=ety1′(t)=et y y2 ′(t)=et+tety2 ′(t)=et+tet (paso 1). Entonces, queremos hallar funciones u′(t)u′(t) y v′(t)v′(t) para que
    u′et+v′tet=0u′et+v′(et+tet)=ett2 .u′et+v′tet=0u′et+v′(et+tet)=ett2 .

    Si aplicamos la regla de Cramer, tenemos
    u′=|0tetett2 et+tet||ettetetet+tet|=0−tet(ett2 )et(et+tet)−ettet=−e2 tte2 t=−1tu′=|0tetett2 et+tet||ettetetet+tet|=0−tet(ett2 )et(et+tet)−ettet=−e2 tte2 t=−1t

    y
    v′=|et0etett2 ||ettetetet+tet|=et(ett2 )e2 t=1t2 (paso 2).v′=|et0etett2 ||ettetetet+tet|=et(ett2 )e2 t=1t2 (paso 2).

    Integrando, obtenemos
    u=−∫1tdt=−ln|t|v=∫1t2 dt=−1t(paso 3).u=−∫1tdt=−ln|t|v=∫1t2 dt=−1t(paso 3).

    Entonces tenemos
    yp=−etln|t|−1ttet=−etln|t|−et(paso 4).yp=−etln|t|−1ttet=−etln|t|−et(paso 4).

    La intersección etet es una solución de la ecuación complementaria, por lo que no necesitamos llevar ese término a nuestra solución general explícitamente. La solución general es
    y(t)=c1et+c2 tet−etln|t|(paso 5).y(t)=c1et+c2 tet−etln|t|(paso 5).
  2. La ecuación complementaria es y″+y=0y″+y=0 con la solución general asociada c1cosx+c2 senx.c1cosx+c2 senx. Así que, y1(x)=cosxy1(x)=cosx como y2 (x)=senxy2 (x)=senx (paso 1). Entonces, queremos hallar funciones u′(x)u′(x) y v′(x)v′(x) de manera que
    u′cosx+v′senx=0−u′senx+v′cosx=3sen2 x.u′cosx+v′senx=0−u′senx+v′cosx=3sen2 x.

    Si aplicamos la regla de Cramer, tenemos
    u′=|0senx3sen2 xcosx||cosxsenx−senxcosx|=0−3sen3xcos2 x+sen2 x=−3sen3xu′=|0senx3sen2 xcosx||cosxsenx−senxcosx|=0−3sen3xcos2 x+sen2 x=−3sen3x

    y
    v′=|cosx0−senx3sen2 x||cosxsenx−senxcosx|=3sen2 xcosx1=3sen2 xcosx(paso 2).v′=|cosx0−senx3sen2 x||cosxsenx−senxcosx|=3sen2 xcosx1=3sen2 xcosx(paso 2).

    Integrando primero para hallar u, obtenemos
    u=∫−3sen3xdx=−3[−13sen2 xcosx+2 3∫senxdx]=sen2 xcosx+2 cosx.u=∫−3sen3xdx=−3[−13sen2 xcosx+2 3∫senxdx]=sen2 xcosx+2 cosx.

    Ahora, integramos para encontrar v. Usando la sustitución (con w=senxw=senx), obtenemos
    v=∫3sen2 xcosxdx=∫3w2 dw=w3=sen3x.v=∫3sen2 xcosxdx=∫3w2 dw=w3=sen3x.

    Entonces,
    yp=(sen2 xcosx+2 cosx)cosx+(sen3x)senx=sen2 xcos2 x+2 cos2 x+sen4x=2 cos2 x+sen2 x(cos2 x+sen2 x)(paso 4).=2 cos2 x+sen2 x=cos2 x+1yp=(sen2 xcosx+2 cosx)cosx+(sen3x)senx=sen2 xcos2 x+2 cos2 x+sen4x=2 cos2 x+sen2 x(cos2 x+sen2 x)(paso 4).=2 cos2 x+sen2 x=cos2 x+1

    La solución general es
    y(x)=c1cosx+c2 senx+1+cos2 x(paso 5).y(x)=c1cosx+c2 senx+1+cos2 x(paso 5).

Punto de control 7.14

Halle la solución general de las siguientes ecuaciones diferenciales.

  1. y″+y=secxy″+y=secx
  2. x″−2 x′+x=ettx″−2 x′+x=ett

Sección 7.2 ejercicios

Resuelva las siguientes ecuaciones utilizando el método de los coeficientes indeterminados.

54.

2 y ″ − 5 y ′ − 12 y = 6 2 y ″ − 5 y ′ − 12 y = 6

55.

3 y ″ + y ′ − 4 y = 8 3 y ″ + y ′ − 4 y = 8

56.

y ″ − 6 y ′ + 5 y = e – x y ″ − 6 y ′ + 5 y = e – x

57.

y ″ + 16 y = e −2 x y ″ + 16 y = e −2 x

58.

y ″ − 4 y = x 2 + 1 y ″ − 4 y = x 2 + 1

59.

y ″ − 4 y ′ + 4 y = 8 x 2 + 4 x y ″ − 4 y ′ + 4 y = 8 x 2 + 4 x

60.

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

61.

y ″ + 2 y ′ + y = sen x + cos x y ″ + 2 y ′ + y = sen x + cos x

62.

y ″ + 9 y = e x cos x y ″ + 9 y = e x cos x

63.

y ″ + y = 3 sen 2 x + x cos 2 x y ″ + y = 3 sen 2 x + x cos 2 x

64.

y ″ + 3 y ′ − 28 y = 10 e 4 x y ″ + 3 y ′ − 28 y = 10 e 4 x

65.

y ″ + 10 y ′ + 25 y = x e −5 x + 4 y ″ + 10 y ′ + 25 y = x e −5 x + 4

En cada uno de los siguientes problemas,

  1. Escriba la forma de la solución particular yp(x)yp(x) para el método de los coeficientes indeterminados.
  2. [T] Utilice un sistema de álgebra computacional para hallar una solución particular a la ecuación dada.
66.

y ″ − y ′ − y = x + e – x y ″ − y ′ − y = x + e – x

67.

y ″ − 3 y = x 2 − 4 x + 11 y ″ − 3 y = x 2 − 4 x + 11

68.

y ″ − y ′ − 4 y = e x cos 3 x y ″ − y ′ − 4 y = e x cos 3 x

69.

2 y ″ − y ′ + y = ( x 2 − 5 x ) e – x 2 y ″ − y ′ + y = ( x 2 − 5 x ) e – x

70.

4 y ″ + 5 y ′ − 2 y = e 2 x + x sen x 4 y ″ + 5 y ′ − 2 y = e 2 x + x sen x

71.

y ″ − y ′ − 2 y = x 2 e x sen x y ″ − y ′ − 2 y = x 2 e x sen x

Resuelva la ecuación diferencial utilizando el método de los coeficientes indeterminados o el de la variación de los parámetros.

72.

y ″ + 3 y ′ − 4 y = 2 e x y ″ + 3 y ′ − 4 y = 2 e x

73.

y ″ + 2 y ′ = e 3 x y ″ + 2 y ′ = e 3 x

74.

y ″ + 6 y ′ + 9 y = e – x y ″ + 6 y ′ + 9 y = e – x

75.

y ″ + 2 y ′ − 8 y = 6 e 2 x y ″ + 2 y ′ − 8 y = 6 e 2 x

Resuelva la ecuación diferencial utilizando el método de variación de los parámetros.

76.

4 y ″ + y = 2 sen x 4 y ″ + y = 2 sen x

77.

y ″ − 9 y = 8 x y ″ − 9 y = 8 x

78.

y ″ + y = sec x , 0 < x < π / 2 y ″ + y = sec x , 0 < x < π / 2

79.

y ″ + 4 y = 3 csc 2 x , 0 < x < π / 2 y ″ + 4 y = 3 csc 2 x , 0 < x < π / 2

Halle la única solución que satisface la ecuación diferencial y las condiciones iniciales dadas, donde yp(x)yp(x) es la solución particular.

80.

y″−2 y′+y=12ex,y″−2 y′+y=12ex, yp(x)=6x2 ex,yp(x)=6x2 ex, y(0)=6,y′(0)=0y(0)=6,y′(0)=0

81.

y″−7y′=4xe7x,y″−7y′=4xe7x, yp(x)=2 7x2 e7x−449xe7x,yp(x)=2 7x2 e7x−449xe7x, y(0)=−1,y′(0)=0y(0)=−1,y′(0)=0

82.

y″+y=cosx−4senx,y″+y=cosx−4senx, yp(x)=2 xcosx+12 xsenx,yp(x)=2 xcosx+12 xsenx, y(0)=8,y′(0)=–4y(0)=8,y′(0)=–4

83.

y″−5y′=e5x+8e−5x,y″−5y′=e5x+8e−5x, yp(x)=15xe5x+425e−5x,yp(x)=15xe5x+425e−5x, y(0)=–2,y′(0)=0y(0)=–2,y′(0)=0

En cada uno de los siguientes problemas, dos soluciones linealmente independientes —y1y1 y y2 y2 —-que satisfacen la ecuación homogénea correspondiente. Utilice el método de variación de los parámetros para hallar una solución particular a la ecuación no homogénea dada. Supongamos que x > 0 en cada ejercicio.

84.

x2 y″+2 xy′−2 y=3x,x2 y″+2 xy′−2 y=3x, y1(x)=x,y2 (x)=x−2y1(x)=x,y2 (x)=x−2

85.

x2 y″−2 y=10x2 –1,x2 y″−2 y=10x2 –1, y1(x)=x2 ,y2 (x)=x−1y1(x)=x2 ,y2 (x)=x−1