UNIVERSIDAD TECNOLÓGICA METROPOLITANA

FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.

DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

UNIDAD 5: ECUACIONES DIFERENCIALES LINEALES DE ORDEN N .

Forma general

dny d n 1 y d n2 y dy
a 0 ( x) n  a1 ( x) n 1  a 2 ( x) n 2      a n 1 ( x)  a n ( x) y  F ( x)
dx dx dx dx

si ai son constantes i , la ecuación se llama ecuación diferencial lineal con coeficientes

constante. En caso contrario ecuación diferencial lineal de coeficientes variables.

Ejemplos
3y’’’+ 2y’’–4y’+8y = 3e-x-5x2 ec. diferencial lineal con coeficientes constantes, de orden 3
xy’’ – 2y’ + x2 y = x3 – 2 ec. diferencial lineal con coeficiente variables, de orden 2.

dy d2y
Usando los operadores : D, D2 , D3 , ….., donde Dy = , D2y = , …...
dx dx 2
Anotaremos

ao Dny + a1 Dn-1y + ---+ an-1 Dy + any = F o bien (ao Dn + a1 Dn-1 + --- + an-1 D + an) y = F

Si  (D) = ao Dn + a1 Dn-1 + --- + an , anotaremos  (D) y = F

Propiedades de los operadores

1) Dn (u+v) = Dn (u) + Dn (v) , u,v funciones diferenciables
2) Dn ( u) =  Dn (u) ,  = constante
3)  (D) (u + v) =  (D)(u) +  (D)(v)
4)  (D) ( u) =   (D) u,  = constante
 Dn y  (D) son operadores lineales
Definición :
Sea  (D) y = F(x) , F(x)  0 ecuación diferencial general de orden n.

Se puede asociar con la ecuación  (D) y = 0 llamada ecuación complementaria o ecuación

homogénea de la ecuación diferencial lineal general.

Teorema fundamental I
Si y = u (x) es cualquier solución de la ecuación  (D) y = F(x) e y = v(x) es cualquier
solución de la ecuación  (D) y = 0, entonces y = u(x) + v(x) es solución de la ecuación  (D)
y = F(x).
Si v(x) tiene n constantes arbitrarias, donde n es el orden de la ecuación, entonces es la solución
general de la ecuación complementaria.
Si u(x) no tiene constantes arbitrarias, es una solución particular de  (D) y = F(x).

63
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Definición : Solución complementaria de la ecuación  (D) y = F(x)

Corresponde a la solución general de la ecuación complementaria y la denotaremos por yc.

Definición : Solución particular de la ecuación  (D) y = F(x)

Es una solución seleccionada de la ec. dada (sin constantes arbitrarias) y la denotaremos por yp

Teorema fundamental II
La solución general de la ecuación diferencial  (D) y = F(x) viene dada por y = yc + y p

Teorema de existencia – unicidad

Si ai (x) y F(x) son funciones continuas en a, b y ao (x)  0 y po, …., pn, son constantes
dadas.
Entonces ! solución y(x), tal que ao(x) Dn + a1(x) Dn-1 + …. + an(x) y = F(x), cumple las
condiciones y(c) = po , y’(c) = p1, ….., y(n-1) (c) = pn, c  a, b

Nota: El teorema proporciona sólo condiciones suficientes, es decir, aún si las condiciones
enunciadas no se satisfacen todas, pueden existir soluciones únicas.

Solución general de la ecuación complementaria de una ecuación diferencial de coeficientes

constantes
Considerar una ecuación auxiliar, que se obtiene al reemplazar el operador D por m en la
ecuación diferencial, se obtiene un polinomio en m y se calculan sus raíces. Se pueden presentar
tres casos.

Caso 1: Raíces reales distintas

Si m1, m2,...., mk son raíces reales distintas de la ecuación auxiliar en m, entonces
y1  e m1x , y 2  e m2 x ......... y n  e mn x son soluciones linealmente independientes de la ecuación

diferencial complementaria  (D) y = 0, y también y = c1 e 1  c 2 e 2    c k e k es solución

mx m x m x

de la ecuación complementaria.

Si la ecuación diferencial es de orden n y se obtiene n raíces reales distintas m 1, m2, ..., mn,
entonces y = c1 e m1 x + ......... + cn e mn x es la solución general de la ecuación complementaria, es
decir, es la solución complementaria.

Ejemplos: Resolver las siguientes ecuaciones diferenciales:

1) y’’ – 3y’ + 2y = 0
Solución:
m2 – 3m + 2 = 0
(m – 1) (m – 2) = 0  m = 1  m = 2  y = ex , y = e2x sea solución
y = c1 ex + c2 e2x es la solución general.

64
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

2) y’’’ – 6y’’ + 11 y’ – 6y = 0
Solución:
m3 – 6m2 + 11 m – 6 = 0
(m – 1) (m – 2) (m – 3) = 0  m = 1, 2, 3.
y = c1 ex + c2 e2x + c3 e3x es la solución general.

Caso 2: Raíces reales iguales

Si la ecuación auxiliar tiene m k como raíz con multiplicidad k, entonces
y1  e mk x , y2  xe mk x ......... yk  x k 1e mk x son soluciones linealmente independientes de la ecuación
diferencial complementaria  (D) y = 0, y también y = (c1 + c2 x + .... + ck xk-1) e mk x es solución
de la ecuación diferencial.

Ejemplos: Resolver las siguientes ecuaciones diferenciales:

1) (D2 – 6D + 9)y = 0
Solución:
m2 – 6 m + 9 = 0, luego (m – 3)2 = 0, es decir m = 3 con multiplicidad 2.
 y = (c1 + c2 x ) e3x solución general.

2) (D3 – 6D2 + 12 D - 8)y = 0

Solución:
m3 – 6m2 + 12 m - 8 = 0, luego (m – 2)3 = 0, es decir m = 2, con multiplicidad 3.
y = (c1 + c2 x + c3 x2) e2x es la solución general.

3) (D6 – 6D5 + 12 D4 – 6D3 – 9D2 + 12D – 4) y = 0

Solución:
m6 – 6m5 + 12m4 – 6m3 – 9m2 + 12m – 4 = 0
(m – 1)3 (m – 2)2 (m + 1) = 0
 m = 1 con multiplicidad 3
m = 2 con multiplicidad 2
m = - 1 con multiplicidad 1

 y = (c1 + c2 x + c3 x2 ) ex + ( c4 + c5 x ) e2x + (c6 e-x) es la solución general.

65
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Caso 3: Raíces complejas

Si la ecuación auxiliar presenta raíces complejas conjugadas, de la forma: m1 = a +ib, m2 = a – ib,
entonces una solución de la ecuación diferencial es

y = eax c1 sen (b x) + c2 cos (b x)

pues e(a + ib)x = eax  eibx = eax (cos b x + i sen bx)

x x2 x3
ex  1    ....
1! 2! 3!
ibx (ibx) 2 (ibx) 3 b 2 x 2 ib 3 x 3
Nota: e
ibx
 1    ....  1  ibx    ....
1! 2! 3! 2! 3!
 b2 x2   b3 x3 
e ibx  1   ....  i bx   ....  cos(bx)  isen(bx)
 2!   3! 
y1 = eax cos(bx) y2 = eax sen(bx) son soluciones independientes

Ejemplos

1) (D2 + 2D + 5) y = 0
m2 + 2m + 5 = 0
m=-1 2i
y = e-x c1 sen 2x + c2 cos 2x

2) (D4 – 5D2 + 12 D + 28) y = 0

m4 – 5m2 + 12 m + 28 = 0
raíces : -2 , -2, 2  3 i
y = (c1 + c2x) e-2x + e2x (c3 sen 3 x + c4 cos 3 x)

3) (D2 – 2D + 5)2 y = 0
(m2 – 2 m + 5)2 = 0
m2 – 2m + 5 = 0  m = 1  2i con multiplicidad 2.
 y = ex (c1 cos 2x + c2 sen 2x) + x ex (c3 cos 2x + c4 sen 2x).

Nota:
2  4  20 2  4i
m2 – 2m + 5 = 0  m    1  2i
2 2

66
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Ejercicios propuestos:

1. Encuentre las soluciones generales de cada una de las siguientes ecuaciones:

a) y’’ + 4y’ – 5y = 0 b) (4D2 – 25)y = 0 c) y’’ = 4y
d) 2y’’’ – 5y’’ + 2y’ = 0 e) I’’(t) – 4I’(t) + 2I(t) = 0 f) (D3 +2D2 – 5D – 6)y = 0

2. Encuentre las soluciones que satisfagan las condiciones dadas.

a) y’’ – y = 0, y(0) = 2, y’(0) = -3 b) (D2 –3D +2)y = 0, y(0) = -1, y’(0) = 0
c) (D3 –16D)y = 0, y(0) = 0, y’(0) = 0, y’’(0) = 16

3. Encuentre las soluciones generales de cada una de las siguientes ecuaciones:

a) (D2 – 4D + 4)y = 0 b) 16y’’ – 8y’ + y = 0 c) 4I’’(t) – 12I’(t) + 9I(t) = 0
d) (D6 – 4D4)y = 0 e) (D4 – 2D3 + D2)y = 0 f) 4y(IV) – 20y’’ + 25y = 0

4. Encuentre las soluciones que satisfagan las condiciones dadas:

a) (D2 – 2D + 1)y = 0; y(0) = 1, y’(0) = -2
b) (D3– D2)y = 0; y(0) = 1, y’(0) = y’’(0) = 0

2
d s ds ds
c)
2 16  64s; s  0,  4 donde t  0
dt dt dt
5. Encuentre la solución general de cada una de las siguientes ecuaciones:
d 2s
a) y’’ + 4y = 0 b) (D2 + 4D + 5)y = 0 c) 4 2  9 s
dt
d) 4y’’ – 8y’ + 7y = 0 e) y(IV) = -16y’’ f) (D3 + D2 - 2)y = 0

6. Encuentre las soluciones que satisfagan las condiciones dadas.

a) (D2 + 1)y = 0; y(0)=4, y’(0)=0 b)U’’(t) = -16U(t); U(0)=0, U’(0) = 4
c) I’’(t) + 2I’(t) + 5I(t) = 0; I(0) = 2, I’(0) = 0

67
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Independencia lineal y Wronskianos

Definición: Llamaremos wronskiano de las funciones y1 e y2 al determinante:

y1 y2
W (y1, y2) = = y1 y’2 – y2 y’1
y1 ' y2 '

Teorema 1
Sean y1, y2 dos soluciones de la ecuación diferencial
(ao D2 + a1 D + a2) y = 0
ao  0, a1, a2 funciones continuas de x en algún intervalo. Entonces
a ( x)
y1 y2  1 dx
W (y1, y2) =  ce ao ( x ) (identidad de Abel)
y '1 y'2

Teorema 2
Sean y1 , y2 soluciones de la ecuación diferencial
(ao D2 + a1 D + a2) y = 0 en algún intervalo I. Entonces
a) y1, y2 son linealmente dependientes (l.d.) en I  W (y1, y2) = 0 en I.
b) y1, y2 son linealmente independiente (l.i.) en I  W (y1, y2)  0 en I.

Teorema 3
Sea y1 una solución de la ecuación diferencial (ao D2 + a1 D + a2) y = 0.
Entonces una solución l. i con y1 está dada por
a1
  ao
dx
e
y2 = y1
 dx
y12

Nota
Este teorema permite determinar una solución y2 conocida otra solución y1.

Ejemplo
Encontrar la solución general de x2 y’’ + xy’ – y = 0, si se sabe que y = x es una solución de la
ecuación diferencial.

Respuesta
ao = x2 , a1 = x, y1 = x
dx
 
= x e x
e  ln x x 1 x 2  1
 dx  x   x2 
3
y2  x dx  x x dx  x 
x2 x2  2 2x
1
x
1 x   2  0, x
 y = c1x + c2  w(y1 , y2) = 1
x 1 x
x2

68
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Teorema 4
Si y1 , y2 son soluciones l.i. de la ecuación diferencial (aoD2 + a1 D + a2)y = 0 e yp es una
solución particular de (a0 D2 + a1 D + a2) y = F(x), entonces y = c1 y1 + c2 y2 es una solución de
la ecuación complementaria e y = c1 y1 + c2 y2 + yp es una solución de la ecuación general.

Nota
Los teoremas anteriores se pueden generalizar al caso de n funciones y ecuación diferencial de
orden n.
(ao Dn + a1 Dn-1 + ….. + an ) y = 0 (*)

y1  yn
y '1  y'n
1) w (y1 ,....., yn) =
 
y1( n 1)  y n( n 1)
a1
2) w (y1,....., yn) = c   ao dx identidad de Abel.
e

3) y1 ,...., yn son l.d. en I  w = 0 en I.

y1 ,...., yn son l.i. en I  w  0 en I.
4) Si y1, ...., yn son soluciones l.i. de la ecuación diferencial () e yp es una solución particular de
la ecuación general (ao Dn + .... + an) y = F(x) entonces y = c1 y1 + .... + cn yn es una solución de
la ecuación diferencial () e y = c1 y1 + ..... + cn yn + yp es una solución de la ecuación diferencial
general.

Ejercicio:
Determine si las siguientes funciones son linealmente independientes en IR.
a) x + 2 , 2x – 3 w = 7  0 l.i.
b) sen 2x , cos2x , 2 w = 0 l.d.
2 3
c) 1, x, x , x . w = 12  0 l.i.

69
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Solución particular de la ecuación diferencial  (D)y = F(x).

Método de los coeficientes indeterminados

Se aplica a la ecuación diferencial (D)y = F(x), donde F(x) contiene un polinomio, términos de
la forma sen(rx), cos(rx), erx , donde r es constante o combinaciones de sumas y productos de
estos.

Caso a : si F(x) contiene la forma erx

Se trata de buscar una función y cuya derivada n-ésima produzca F(x). Luego yp = a erx

Ejemplo:
Resolver y’’ + 4y = 4 e2x
Respuesta:
Sea y = a e2x una posible solución
y’ = 2a e2x
y’’ = 4a e2x

reemplazando en la ecuación, se tiene

4a e2x + 4a e2x = 4 e2x,
luego 8a e2x = 4 e2x,
1
:. a=
2
1 2x
 yp = e es la solución particular.
2

La ecuación complementaria y’’ + 4y = 0 tiene solución general yc = c1 cos 2x + c2 sen 2x

En efecto: m2 + 4 = 0  m =  2i raíces complejas conjugadas.

e2x
la solución general de la ecuación dada es y = yc + yp = c1 cos 2x + c2 sen 2x +
2

Caso b : si F(x) contiene la forma sen(rx) o cos(rx)

Se trata de buscar una función y cuya derivada n-ésima produzca sen(rx) o cos(rx).
Luego
yp = a sen(rx) + b cos(rx)
______________________________________________________________________________

Ejemplo:
Resolver (D2 + 4D + 4) y = 6 sen(3x)
Respuesta:
Sea yp = a sen(3x) + b cos(3x)
yp’ = 3a cos(3x) - 3b sen(3x)
yp’’ = - 9a sen(3x) - 9b cos(3x).

70
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Reemplazando en la ecuación dada, se tiene

- 9a sen(3x) – 9b cos(3x) + 4 (3a cos(3x) – 3b sen (3x)) + 4 (a sen(3x) + b cos(3x)) = 6 sen(3x)

(-5a – 12b) sen(3x) + (12a – 5b) cos(3x) = 6 sen(3x).
 (- 5a – 12b) = 6  (12a – 5b) = 0
 30  72
Resolviendo el sistema se tiene: a = , b=
169 169
 30  72
 yp = sen 3x + cos 3x solución particular
169 169

La ecuación complementaria (D2 + 4D + 4)y = 0 tiene solución general yc = (c1 + c2 x) e-2x

En efecto: m2 + 4m + 4 = 0  (m + 2)2 = 0  m = -2 con multiplicidad 2.

 la solución general de la ecuación dada es :

 30  72
y = yc + yp = (c1 + c2 x) e-2x + sen 3x + cos 3x
169 169

Caso c: si F(x) contiene la forma de un polinomio de grado n.

Se trata de buscar una función y, cuya derivada n-ésima produzca F(x).
Luego yp = an xn + ….. + ao

Ejemplo:
Resolver (D2 + 4D + 9)y = x2 + 3x.
Respuesta:
Sea y = ax2 + bx + c
y’ = 2ax + b
y’’ = 2a.

Reemplazando en la ecuación dada, se tiene: 2a + 8ax + 4b + 9ax2 + 9bx + 9c = x2 + 3x.

Resolviendo el sistema se tiene:
9a = 1
8a + 9b = 3
2a + 4b + 9c = 0

1 19  94
Luego: ,b 
a= ,c 
9 81 729
1 2 19 94
 yp = x  x  solución particular
9 81 729

La ecuación complementaria (D2 + 4D + 9) y = 0 tiene solución general

yc = e-2x (c1 cos 5 x + c2 sen 5 x)

71
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

En efecto:
m2 + 4m + 9 = 0  m = - 2  5 i raíces complejas conjugadas

 la solución general de la ecuación dada es :

1 2 19 94
y = yc + yp = e-2x (c1 cos ( 5 x) + c2 sen ( 5 x) + x  x
9 81 729

Ejercicio
Resolver (D2 + 2D + 1) y = 2 cos(2x) + 3x + 2 + 3 ex
Respuesta:
Sea yp = a sen (2x) + b cos( 2x) + cx + d + f ex
y’ = 2a cos (2x) – 2b sen (2x )+ c + f ex
y’’ = 4a sen( 2x) – 4 b cos( 2x) + f ex

Reemplazando en la ecuación dada.

(-3a - 4b) sen 2x + (4a - 3b) cos 2x + cx + d + 2c + 4 f ex = 2 cos 2x + 3x + 2 + 3 ex

 - 3a - 4b = 0
4a - 3b = 2
8 6
c=3  a= , b , c = 3, d = - 4, f =
3
4
25 25
 8 6 3
yp = sen 2x - cos 2x + 3x – 4 + ex
25 25 4

La ecuación complementaria (D2 + 2D + 1)y = 0 tiene solución general yc = (c1 + c2 x ) e-x

En efecto:
m2 + 2m + 1 = 0  (m + 1)2 = 0  m = -1 con multiplicidad 2.

 la solución general de la ecuación dada es :

8 6 3 x
y = yc + yp = (c1 + c2 x) e-x + sen 2x - cos 2x + 3x – 4 + e
25 25 4

El método aniquilador (Justificación al método de los coeficientes indeterminados)

Dada la ecuación diferencial  (D) y = F(x) se trata de buscar un operador tal que aplicado a la
ecuación dada, se obtenga una ecuación con el lado derecho igual a cero. El operador requerido
se llama operador de aniquilación o aniquilador.

Ejemplos
1) Resolver (D2 + 4) y = 4 e2x

Respuesta: (D2 + 4) y = 4 e2x /D

72
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

D (D2 + 4) y = 8 e2x

73
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

 2 (D2 + 4) y = 8 e2x
D (D2 + 4) y = 8 e2x

Aniquilador: (D – 2)

Restando: (D – 2) (D2 + 4)y = 0  (m – 2) (m2 + 4) = 0 ecuación auxiliar

m = 2, m =  2i

y = c1 cos 2x + c2 sen 2x + c3 e2x solución general de 

yc = c1 cos 2x + c2 sen 2x , yp = c3 e2x

1
Reemplazando yp en la ecuación dada, se tiene c3 =
2

2) Resolver (D2 + 4D + 4) y = 6 sen 3x

Respuesta:
Aniquilador: (D2 + 9)
(D2 + 9) (D2 + 4D + 4) y = (D2 + 9) (6 sen 3x)
(D2 + 9) (D2 + 4D + 4) y = 0
(m2 + 9) (m2 + 4m) = 0  m =  3i , m = -2 con multiplicidad 2

y = (c1 + c2 x) e-2x + c3 sen 3x + c4 cos 3x

yp = c3 sen 3x + c4 cos 3x
 30  72
Reemplazando yp en la ecuación dada.: c3 = , c4 =
169 169

3) Resolver (D2 + 2D + 1) y = 2 cos 2x + 3 x + 2 + 3 ex

Respuesta
para 2 cos 2x usamos D2 + 4
para 3x + 2 usamos D2
x
para 3e usamos D - 1

 D2 (D – 1) (D2 + 4) (D2 + 2D + 1) y = 0
 y = (c1 + c2 x) e-x + c3 sen 2x + c4 cos 2x + c5 ex + c6 x + c7
yp = c3 sen 2x + c4 cos 2x + c5 ex + c6 x + c7

Reemplazando yp en la ecuación dada, se tiene

8 6 3
c3 = c4 = c5 = c6 = 3 c7 = - 4.
25 25 4

Ejercicios propuestos:

74
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Resolver las siguientes ecuaciones diferenciales

1) (D3 + 4D) y = ex + sen x

2) s’’(t) – 3s’ (t) + 2s (t) = 8t2 + 12 e-t , s(o) = 0, s’(o) = 2
3) y’’ + y = 6 cos2x , y(0) = 0 y (/2) = 0
4) y’’ – 3y’ + 2y = 4 sen3 3x
5) Demuestre el uso del método de aniquilación para
a) y’’ + y = 2 e3x
b) (D2 + 2D + 1) y = 4 sen 2x
c) (D2 – 4) y = 8x2
Excepciones en el método de los coeficientes indeterminados

Ejemplo:
Resolver (D2 + 3D + 2)y = 4 e-2x
a) Usando el método de los coeficientes indeterminados
b) Usando aniquiladores
Respuesta
a) Usando el método de los coeficientes indeterminados
ecuación auxiliar m2 + 3 m + 2 = 0  m = - 1 , m = - 2
 yc = c1 e-ex + c2 e-x
sea y = a e-2x
y’ = - 2a e-2x
y’’ = + 4 a e-2x
reemplazando
4a e-2x + 3 (-2a e-2x ) + 2ae-2x = 4e-2x
4a e-2x - 6a e-2x + 2ae-2x = 4e-2x
0 = 4e-2x contradicción
Mejor considerar y = ax e , si y = a e-2x está en yc
-2x

(D2 + 3D + 2) (ax e-2x ) = - a e-2x = 4 e-2x  a = - 4

 y = c1 e-2x + c2 e-x - 4x e-2x

b) Por el método de aniquilación

(D2 + 3D + 2) y = 4 e-2x / (D + 2)
(D + 2) (D2 + 3D + 2) y = - 8 e-2x + 8 e-2x = 0
(m + 2) (m + 2) (m + 1) = 0
 y = c1 e-2x + c2 e-x + c3 x e-2x , c3 = -4

Observación
Si un término de yp está en yc, se debe multiplicar el término de y p por una potencia mínima de x
para que esto no ocurra.

Ejemplos
1) Resolver (D2 + 4) y = 6 sen 2x + 3x2
Respuesta :
ecuación auxiliar m2 + 4 = 0  m =  2i
 yc = c1 cos (2x) + c2 sen (2x), f(x) = 6 sen(2x) + 3x2

75
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

como sen (2x) aparece en yc.

sea yp = x (a sen (2x) + b cos (2x) ) + cx2 + dx + f.
y’p = (a sen (2x) + b cos (2x) ) + x (2a cos (2x) - 2b sen (2x)) + 2cx + d
y’’p = 2a cos(2x) - 2b sen(2x) + 2a cos(2x) - 2b sen(2x) + x(-4a sen(2x) - 4b cos(2x)) + 2c
= 4a cos (2x) - 4b sen (2x) + x (-4a sen (2x) - 4b cos (2x)) + 2c
al sustituir en la ecuación diferencial queda:
4a cos(2x) - 4b sen(2x) + x(-4a sen(2x) - 4b cos(2x)) + 2c + 4x(a sen(2x) + b cos(2x)) + cx 2 + dx
+ f) = 6 sen 2x + 3 x2
4a cos 2x – 4b sen 2x + 4cx2 + 4dx + 4f + 2c = 6 sen 2x + 3 x2
 4a = 0 , - 4b = 6, 4c = 3 , 4d = 0 , 4f + 2c = 0
3 3 3
 a=0, b=- , c= , d=0 f= -
2 4 8
3 3 2 3
 yp = - x cos 2x + x -
2 4 8
 3 3 2 3
y = c1 sen 2x + c2 cos 2x - x cos 2x + x -
2 4 8

2) Resolver (D3 – 3D2 + 3D – 1) y = 2 ex

Respuesta:
ecuación auxiliar (m3 – 3 m2 + 3m – 1) = 0
(m – 1)3 = 0
2 x
yc = (c1 + c2 x + c3 x ) e
como ex está en yc y x ex está en yc y x2 ex está en yc , sea yp = a x3 ex
1
 (D3 – 3D2 + 3D – 1) (a x3 ex) = 6a ex = 2 ex  a =
3
 1 3 x
y = (c1 + c2 x + c3 x2 ) ex + x e.
3

Método de variación de parámetro

Dada la ecuación diferencial

y’’ + P(x) y’ + Q(x) y = F(x)

76
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

sabemos que tiene como solución general

y (x) = yc + yp

donde yc = c1 y1(x) + c2 y2(x) e yp es una solución particular

El método consiste en reemplazar las constantes c1 y c2 por las funciones v1 (x) y v2 (x), tal
que y = v1 y1 + v2 y2 sea una solución particular de la ecuación dada.

derivando : y’ = v1’ y1 + v1 y1’ + v2’ y2 + v2 y2’

ordenando : y’ = ( v1 y1’ + v2 y2’) + (v1’ y1 + v2’ y2)

para simplificar, hacemos v1’ y1 + v2’ y2 = 0

 y’ = v1 y1’ + v2 y2’

derivando y’
queda y’’ = v1 y1’’ + v1’ y1’ + v2 y2’’ + v2’ y2’

Reemplazando en la ecuación dada, se tiene

(v1 y1’’ + v1’ y1’ + v2 y2’’ + v2’ y2’) + P(x) (v1 y1’ + v2 y2’) + Q(x) (v1 y1 + v2 y2) = F(x).

Ordenado queda

v1 (y1’’ + P(x) y1’ + Q(x)y1) + v2 (y2’’ + P(x) y2’ + Q(x)y2) + v1’ y1’ + v2’ y2’ = F(x)
 v1  0 + v2  0 + v1’ y1’ + v2’ y2’ = F(x)

pues y1, y2 son soluciones de la ecuación diferencial complementaria.

Luego v1’y1’ + v2’y2’ = F(x)

y como agregamos la condición v1’y1 + v2’y2 = 0 se tiene el sistema

v1’ y1’ + v2’ y2’ = F(x) v1’ y1 + v2’ y2 = 0


v1’ y1 + v2’ y2 = 0 v1’ y1’ + v2’ y2’ = F(x)

0 y2 y1 0
por Regla de Cramer v1’ = F y2 ' , v2’ = y1 ' F
w w
 y 2 F ( x) y1 F ( x )
de donde v1’ = , v2’ =
w( y1 , y 2 ) w( y1 , y 2 )
 y 2 F ( x) y F ( x)
:. v1 =  dx v2 =  1 dx
w w

77
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

 y 2 F ( x) y1 F ( x )
 yp = y1  w
dx + y2  w
dx es la solución particular buscada.

Ejemplo
Resolver y’’ + 4y’ + 4y = x e2x
Respuesta
Solución complementaria
y’’ + 4y + 4y = 0
m2 + 4m + 4 = 0  (m + 2)2 = 0  m = - 2 con multiplicidad 2.

yc = (c1 + c2 x) e-2x
Solución particular
sea y1 = e-2x , , y2 = x e-2x
-2x
y1’ = -2e , y2’ = e-2x (1 – 2x)

yp = v1 (x) y1 (x) + v2 (x) y2(x)

e 2 x xe 2 x
w (y1 , y2) = = e-4x
 2e  2 x e  2 x (1  2 x)
 xe 2 x xe 2 x e 2 x xe 2 x
 yp = e-2x  e 4 x dx + x e-2x  e 4 x dx
x 1
 x e 4 x dx  xe 2 x  xe 4 x dx = e2x   
2
yp = e-2x
16 32 
x 1
 y = yc + yp = e-2x (c1 + c2 x ) + e2x   
16 32 

Ejercicios propuestos:

1) Encuentre la solución general de las siguientes ecuaciones diferenciales con coeficientes

constantes.
a) y’’ + 4y = tg 2x
b) y’’ – y’- 6y = e-x
c) y’’ + y = cosec x
d) y’’ + 2y’ + y = e-x ln x
e) y’’ – 2y’ – 3y = 64x e-x
f) y’’ + 2y’ + 5y = e-x sec 2x

2) Encuentre la solución general de las siguientes ecuaciones diferenciales con coeficientes

variables.
a) (x2 – 1) y’’ – 2xy’ + 2y = (x2 - 1)2 , si y = x es solución de la ec. complementaria.
b) (x2 + x)y’’ + (2 – x2)y’ – (2 + x)y = x(x + 1)2, si y = ex es sol. de la ec. complementaria

Observación
El método de variación de parámetros puede extenderse a ecuaciones diferenciales de orden
arbitrario.

Dada la ecuación diferencial de orden n.

78
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

y(n) + a 1 (x) y(n-1) + ............+ an(x) = F(x)

y supongamos que se conoce la solución complementaria.

yc = c1 y1(x) + ...... + cn yn (x)

se busca una solución particular de la forma.

yp = c1(x)y1(x) + ...... + cn(x) yn (x)

tal que
c1’ y1 + .................. + cn’ yn = 0

c1’ y1’ + .................. + cn’ yn’ = 0



c1‘ y 1( n  2 ) + ........... + cn’yn(n-2) = 0
c1’ y 1( n 1) + ...... + c n ' y n
( n 1)
= F(x)

sea W y1, ....., yn = Wronskiano de y1, ....., yn

y sea Vk ( x) = determinante obtenido de W al reemplazar su columna k por

0
 
0
 
 
0
1
 
Vk ( x ) F ( x ) Vk ( x)  F ( x)
Entonces c k' ( x ) 
w
 ck ( x)   w
dx

V1 ( x ) F ( x ) V ( x ) F ( x)
Luego: yp = y1(x)  dx  .......  y n ( x )  n dx es la solución particular
w w
buscada.

Ejemplos
1) Hallar una solución particular yp para la ecuación 3 y’’’ + 5y’’ – 2y’ = ex
Respuesta
5 2 ex ex
y ' ' ' y' ' y'  , luego F(x) =
3 3 3 3

Ecuación complementaria: 3y’’’ + 5y’’ – 2y’ = 0

79
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

Ecuación auxiliar: 3 m3 + 5 m2 – 2m = 0
m (3 m2 + 5 m- 2) = 0
3 m (m + 2) (m – 1/3) = 0

 yc = c1 + c2 e-2x + c3 ex/3

y1 = 1 , y2 = e-2x , y3 = ex/3

V1 ( x) F ( x) V ( x) F ( x)dx V ( x) F ( x)dx
yp = y1( x )  dx  y 2  2  y3  3
w w w

x x

1 e 2 x e3 0 e 2 x e3
x x
5 x 5 x
2 x e3
14 2 x e3
7
w = 0  2e  e 3
V1 = 0  2e  e 3
3 9 3 3
x x

e3 e3
0 4e  2 x 1 4e  2 x
9 9

1 0 e 3
x x
1 e 2 x 0
e3 e3
V2 = 0 0  V3 = 0  2e  2 x 0  2e  2 x
3 3 0 4e  2 x 1
x

e 3
0 1
9
x
5
7 3x ex e3 ex ex
e   x  2e  2 x 
 yp = 1 3 3 dx  e  2 x 3 3 dx  e 3 3 dx
  14 3
5 x   14 3
5
x
  14 3
5 x
e e e
9 9 9

1 x 1 x 9 x 1 x
yp = - e  e  e  e
2 42 14 6

2) Resolver y’’’ – 2y’’ – y’ + 2y = ex

Respuesta:
Ec. complementaria (D3 – 2D2 – D + 2) y = 0
Ec. auxiliar m3 – 2m2 – m + 2 = 0  m = 2, 1, - 1.

 yc = c1 e2x + c2 ex + c3 e-x

sean y1 = e2x , y2 = ex , y3 = e-x

80
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

e 2x ex ex 1 1 1 1 1 1
1 3
W= 2e 2 x ex  ex  e2x  e xe x 2 1  1  e2x 0 1  3  e2x  e 2 x (3  9) 
3 3
4e 2 x ex ex 4 1 1 0 3 3
 6e 2 x
0 ex ex
ex e x
V1 = 0 ex  e x  e  2
e  ex
1 ex ex
e 2x 0 ex
e2x ex
V2 = 2e 2 x 0  e x     (  e x  2e x )  3e x
2e 2 x  ex
4e 2 x 1 e x

e2x ex 0
e2x ex
V3 = 2e 2 x ex 0   e 3 x  2e 3 x   e 3 x
2e 2 x ex
4e 2 x ex 1

 2e x 3e x e x  e 3x e x
  6e 2 x   6e 2 x   6e 2 x dx
x
yp = e2x dx  e x
dx  e

2e 2 x  e x x ex 1 x 1
yp =  e x     ex  ex  ex
6 2 6 3 2 12

1 x x x x x
yp = 
4
e  e
2
como ex  yc , yp = -
2
e

x x
y = c1 e2x + c2 ex + c3 e-x - e
2

Ejercicios propuestos de ecuaciones diferenciales lineales no homogéneas

Resuelva cada una de las siguientes ecuaciones diferenciales sujetas a condiciones dadas si las
hay.
d2y dy
2
1) y’’ + 3y = x + 1; y(0) = 0, y’(0) = 2 2) 2  3  2 y  senx
dx dx
3) (D2 + 2D + 1)y = e x  e  x 4) y' ' ' - 4y  4x  2  3e -2x
d 2I dI d 4x
5) 2  2  5I  34 cos 2t 6) 4
 x  8e t
dt dt dt
7) y ' '4 y  xe 2 x ; y (0)  y ' (0)  0 8) x 2 y ' '6 y  0; y (1)  2, y ' (1)  0

81
 UNIVERSIDAD TECNOLÓGICA METROPOLITANA
FACULTAD DE C. NAT, MATEMÁTICAS Y DEL M. AMB.
DEPARTAMENTO DE MATEMÁTICAS
LIDIA ORTEGA SILVA

d3y d2y
9)  2 1 10) y ( IV )  16 y ' '  64 cos 4 x
dx 3 dx 2
d 2r
11) y’’  4 y  x(1  cos x ) 12)  2r  e  2
d 2
13) y’’’ - 4y’’ + 4y’ = 12 e 2 x  24x 2 14) y’’ + y = secx, y(0) = 1, y’(0) = 2
d 4s d 2s
15) x2y’’ – 4xy’ + 4y = 24(x+1) 16) 4  2 2  s  100 cos 3t
dt dt
17) 4y’’ – 4y’ + y = ln x 18) D(D – 1)(D2 – 4)y = X2 – X + ex
2

4 2
d I d I ln x
19)  9 2  20e .t ; I (0)  I ' (0)  0 20) x2y’’’ – xy’’ + y’ =
dt 4
dt x

82