ECUACIONES DIFERENCIALES ORDINARIAS

SISTEMAS DE ECUACIONES DIFERENCIALES

Hidalgo Zurita, X.

Resumen

En este pequeño artículo se presenta un resumen de las técnicas más comunes para en-
contrar soluciones de ecuaciones diferenciales ordinarias, se presenta en cada caso sendos
ejemplos que ilustran el método empleado, intentando que el mismo sea de fácil entendi-
miento para estudiantes de niveles básicos de una carrera universitaria. Este compendio
es una recopilación de notas de aula y resúmenes de varios autores.

Índice

1
 SISTEMAS DE ECUACIONES DIFERENCIALES

Los Sistemas de Ecuaciones Diferenciales se presentan con frecuencia como modelos de sistemas
mecánicos o eléctricos, y en general en diversos problemas de ingeniería.

Método de Solución por Eliminación

El método de Eliminación para resolver sistemas de ecuaciones diferenciales con coeficientes
constantes se basa en el principio de los sistemas de ecuaciones algebraicas lineales, mediante la
eliminación de variables. Consiste en que las funciones desconocidas y sus derivadas se eliminan
sucesivamente hasta que se llega a una ecuación diferencial de orden superior que contiene solo
una función desconocida y sus derivadas. Se resuelve esta ecuación y se encuentra paso a paso
las otras funciones desconocidas. Se aplica a sistemas de la forma:
 dx  2
dx
 = a1 x + b1 y + f1 (t);  2 = a1 x + b1 y + f1 (t);


dt

dt
 dy = a x + b y + f (t)
 

2
 d y = a2 x + b2 y + f2 (t)
2 2 2
dt dt2

Si f1 (t) = f2 (t) = 0 son sistemas de ecuaciones diferenciales homogéneas.

 dx
 = 4x − 2y; (1)
dt

Ejemplo 1. Resolver:
 dy = x + y

(2)
dt
 
1 dx
De la ecuación (1) despejamos y: y= 4x −
2 dt
2
dx dx dy
derivando (1) respecto a t: 2
=4 −2 (3)
dt dt dt
d2 x dx
reemplazando (2) en (3): 2
= 4 − 2(x + y)
dt dt
reemplazando y en la última ecuación se tiene:
d2 x
 
dx dx
= 4 − 2x − 4x −
dt2 dt dt
d2 x dx dx
con lo que: = 4 − 6x +
dt2 dt dt
d2 x dx
es decir: − 5 + 6x = 0
dt2 dt
Resolvemos la ecuación homogénea obtenida:
la ecuación característica es:
r2 − 5r + 6 = 0 ⇒ (r − 3)(r − 2) = 0 ⇒ r1 = 3; r2 = 2
y la solución homogénea: x = c1 e3t + c2 e2t

2
 dx
derivando respecto a t: = 3c1 e3t + 2c2 e2t
dt
dx
reemplazando x y en y:
dt
1 1
y = 2(c1 e3t + c2 e2t ) − (3c1 e3t + 2c2 e2t ) = c1 e3t + c2 e2t
2 2
Entonces, lasolución del sistema de ecuaciones diferenciales es:
 x(t) = c e3t + c e2t
1 2
1
 y(t) = c e3t
2 1
+ c2 e2t

 dx
 + 2x + 4y = 1 + 4t; (1)
dt

Ejemplo 2. Resolver:
 dy + x − y = 3 t2

(2)
dt 2
 
1 dx
Despejando y de (1): y= 1 + 4t − − 2x
4 dt
d2 x dx dy
Derivando (1) respecto a t: 2
+2 +4 = 4; (3)
dt dt dt
dy
reemplazando y y en (3):
dt
2
dx dx 3
2

+ 2 + 4 t − x + y =4
dt2 dt 2
d2 x
 
dx 2 dx
entonces: + 2 + 6t − 4x + 1 + 4t − − 2x = 4
dt2 dt dt
d2 x dx
Es decir: 2
+ − 6x = 3 − 4t − 6t2
dt dt
Resolviendo la ecuación homogénea, la ecuación característica es:
r2 + r − 6 = 0 ⇒ (r + 3)(r − 2) = 0 ⇒ r1 = −3; r2 = 2
y la solución homogénea: xh = c1 e−3t + c2 e2t
para encontrar la solución particular, suponemos:
xp = At2 + Bt + C
dx d2 x
derivando: = 2At + B ⇒ = 2A
dt dt2
Reemplazando en la ecuación diferencial:
2A + (2At + B) − 6(At2 + Bt + C) = −6t2 − 4t + 3
igualando coeficientes de potencias semejantes:
t2 : −6A = −6 ⇒ A=1
t: 2A − 6B = −4 ⇒ B=1
t0 : 2A + B − 6C = 3 ⇒ C=0
por lo tanto: xp = t2 + t

3
 y entonces: x(t) = xh + xp = c1 e−3t + c2 e2t + t2 + t
dx
derivando: = −3c1 e−3t + 2c2 e2t + 2t + 1
dt
dx
reemplazando x y en y, se obtiene:
dt
1 1
− 3c1 e−3t + 2c2 e2t + 2t + 1 − 2 c1 e−3t + c2 e2t + t2 + t



y = +t−
4 4
1 1
y = c1 e−3t − 2c2 e2t − t2
4 2
Entonces, lasolución del sistema de ecuaciones diferenciales es:
 x(t) = c1 e−3t + c2 e2t + t2 + t
 y(t) = 1 c1 e−3t − c2 e2t − 1 t2
4 2

 2
dx
 2 = y + 4e−2t (1)


dt
Ejemplo 3. Resolver: 2
 dy
= x − e−2t (2)


dt2
d2 x
Despejando y de (1): y= − 4e−2t
dt2
d3 x dy
Derivando (1) respecto a t: 3
= − 8e−2t
dt dt
d4 x d2 y
= + 16e−2t (3)
dt4 dt2
reemplazando (2) en (3):
d4 x
4
− x = 15e−2t
dt
Resolviendo la ecuación homogénea: la ecuación característica es:
r4 − 1 = 0 ⇒ (r2 − 1)(r2 + 1) = 0 ⇒ (r − 1)(r + 1)(r2 + 1) = 0
⇒ r1 = −1; r2 = 1; r3,4 = ±i
y la solución homogénea: xh = c1 e−t + c2 et + c3 cos(t) + c4 sen(t)
para encontrar la solución particular, suponemos:
xp = Ae−2t
derivando:
dx −2t d2 x −2t d3 x d4 x
= −2Ae ; 2
= 4Ae ; 3
= −8Ae−2t ; 3
= 16Ae−2t
dt dt dt dt
reemplazando en la ec. diferencial, se obtiene:
16Ae−2t − Ae−2t = 15e−2t ⇒ A=1
por lo tanto: xp = e−2t
y entonces: x(t) = xh + xp = c1 e−t + c2 et + c3 cos(t) + c4 sen(t) + e−2t

4
 dx
derivando: = −c1 e−t + 2c2 et − c3 sen(t) + c4 cos(t) − 2e−2t
dt
d2 x
= c1 e−t + c2 et − c3 cos(t) − c4 sen(t) + 4e−2t
dt2
d2 x
reemplazando en y
dt2
y = c1 e−t + c2 et − c3 cos(t) − c4 sen(t) + 4e−2t − 4e−2t

Entonces, lasolución del sistema de ecuaciones diferenciales es:

 x(t) = c et + c e−t + c cos(t) + c sen(t) + e−2t
1 2 3 4
 y(t) = c et + c e−t − c cos(t) − c sen(t)
1 2 3 4

Uso de Operadores en la eliminación de incógnitas.

d
Usando la notación del operador D = se puede reemplazar las ecuaciones diferenciales en
dt
ecuaciones algebraicas, tratando a D como un simple operador algebraico.
 dx
 + 2x − y = 0 (1)
dt

Ejemplo 4. Resolver:
 dy + 4x − 3y = 10 cos(t)

(2)
dt
Transformando las ecuaciones diferenciales en ecuaciones algebraicas mediante el operador
D se tiene:
(D + 2)x − y = 0 (3)
4x + (D − 3)y = 10 cos(t) (4)
resolvemos con el método de eliminación: ec. (3) ∗ (D − 3) + ec. (4)
(D − 3)(D + 2)x − (D − 3)y = 0
4x + (D − 3)y = 10 cos(t)
entonces:
(D − 3)(D + 2)x + 4x = 10 cos(t)
(D2 − D − 6 + 4)x = 10 cos(t)
(D2 − D − 2)x = 10 cos(t)
por lo tanto la ecuación diferencial obtenida es:
d2 x dx
− − 2x = 10 cos(t) (∗)
dt2 dt
Resolvemos la ecuación homogénea:
ecuación característica: r2 − r − 2 = 0
cuyas raíces son: (r − 2)(r + 1) = 0 ⇒ r1 = 2; r2 = −1
y la solución homogénea: xh = c1 e2t + c2 e−t
para encontrar la solución particular, suponemos:
xp = A cos(t) + B sen(t)

5
 derivando:
dx d2 x
= −A sen(t) + B cos(t); = −A cos(t) − B sen(t)
dt dt2
reemplazando en la ec. diferencial (∗):
−A cos(t) − B sen(t) − (−A sen(t) + B cos(t)) − 2(A cos(t) + B sen(t)) = 10 cos(t)
agrupando en senos y cosenos:
sen(t) : A − 3B = 0 ⇒ A = 3B
cos(t) : −3A − B = 10 ⇒ B = 1; A = −3
por lo tanto: xp = −3 cos(t) − sen(t)
y entonces, x(t) = xh + xp = c1 e2t + c2 e−t − 3 cos(t) − sen(t)
dx
derivando: = 2c1 e2t − c2 e−t + 3 sen(t) − cos(t)
dt
dx
despejando y de la ec. (1) y reemplazando x y :
dt
y = 2c1 e2t − c2 e−t + 3 sen(t) − cos(t) + 2(c1 e2t + c2 e−t − 3 cos(t) − sen(t))
es decir: y(t) = 4c1 e2t + c2 e−t + sen(t) − 7 cos(t)
Entonces,
 la solución del sistema de ecuaciones diferenciales es:
 x(t) = c e2t + c e−t − 3 cos(t) − sen(t)
1 2
 y(t) = 4c e2t + c e−t + sen(t) − 7 cos(t)
1 2


dx dy
+ = cos(t) (1)


dt dt

Ejemplo 5. Resolver:
 d2 y
 x+
 =2 (2)
dt2
Transformo las ecuaciones diferenciales mediante el operador D:
Dx + Dy = cos(t) (3)
x + D2 y = 2 (4)
resolviendo: ec. (3) ∗ D − ec. (4)
D2 x + D2 y = D(cos(t))
−x − D2 y = −2
entonces: D2 x − x = −2 − sen(t)
por lo tanto la ecuación diferencial obtenida es:
d2 x
− x = −2 − sen(t) (∗)
dt2
ecuación característica: r2 − 1 = 0 ⇒ r = ±1
y la solución homogénea: xh = c1 et + c2 e−t
para encontrar la solución particular, suponemos:

6
 xp = A + B cos(t) + C sen(t)
derivando:
dx d2 x
= −B sen(t) + C cos(t); = −B cos(t) − C sen(t)
dt dt2
reemplazando en la ec. diferencial (∗):
−B cos(t) − C sen(t) − (A + B cos(t) + C sen(t)) = −2 − sen(t)
agrupando e igualando coeficientes:
t0 : −A = −2 ⇒ A=2
cos(t) : −B − B = 0 ⇒ B=0
1
sen(t) : −C − C = −1 ⇒ C =
2
1
por lo tanto: xp = 2 + sen(t)
2
1
y entonces, x(t) = c1 et + c2 e−t + 2 + sen(t)
2
dx 1
derivando: = c1 et − c2 e−t + cos(t)
dt 2
dy dx
despejando de la ec. (1) y reemplazando :
dt dt
dy 1 1
= cos(t) − (c1 et − c2 e−t + cos(t)) = −c1 et + c2 e−t + cos(t)
dt 2 2
1
integrando: y(t) = −c1 et − c2 e−t + sen(t)
2
Entonces,
 la solución del sistema de ecuaciones diferenciales es:
 x(t) = c1 et + c2 e−t + 2 + 1 sen(t)

2
1
 y(t) = −c1 et − c2 e−t + sen(t)

2

Métodos abreviados del operador

Para minimizar el trabajo involucrado en la eliminación de variables mediante operadores, se
puede usar determinantes, es decir, el uso de la Regla de Cramer para sistemas de ecuaciones
algebraicas lineales. Es recomendado principalmente para sistemas no homogéneos.
 dx
 − 3x − 6y = t2 (1)
dt

Ejemplo 6. Resolver:
 dy + dx − 3y = et

(2)
dt dt
Transformo las ecuaciones diferenciales mediante el operador D:
(D − 3)x − 6y = t2 (3)
Dx + (D − 3)y = et (4)
Usando la regla de Cramer:
Solución en x:

7
  
2
−6 

t
 
 t
(D − 3)

e (D − 3)(t2 ) + 6et 2t − 3t2 + 6et
x=  = =
(D − 3)

−6 
 (D − 3)2 + 6D D2 + 9
 
(D − 3)
 
 D
entonces: (D2 + 9)x = 2t − 3t2 + 6et
por lo tanto la ecuación diferencial obtenida es:
d2 x
+ 9x = 2t − 3t2 + 6et (∗)
dt2
la ecuación característica: r2 + 9 = 0 ⇒ r = ±3i
solución homogénea: xh = c1 cos(3t) + c2 sen(3t)
para encontrar la solución particular, suponemos:
xp = Aet + Bt2 + Ct + D
derivando:
dx d2 x
= Aet + 2Bt + C, = Aet + 2B
dt dt2
reemplazando en la ec. diferencial (∗):
Aet + 2B + 9(Aet + Bt2 + Ct + D) = 6et + 2t − 3t2
igualando coeficientes:
3
et : A + 9A = 6 ⇒ A=
5
1
t2 : 9B = −3 ⇒ B=−
3
2
t: 9C = 2 ⇒ C=
9
2
t0 : 2B + 9D = 0 ⇒ D = −
27
3 t 1 2 2 2
por lo tanto: xp = e − t + t −
5 3 9 27
3 1 2 2
y entonces, x(t) = c1 cos(3t) + c2 sen(3t) + et − t2 + t −
5 3 9 27
Solución en y:
 
(D − 3)
 2 
t
 
et 
 
 D (D − 3)(et ) − D(t2 ) et − 3et − 2t
y= = =
(D − 3)

−6 
 (D − 3)2 + 6D D2 + 9
 
(D − 3)
 
 D
entonces: (D2 + 9)y = −2t − 2et
por lo tanto la ecuación diferencial obtenida es:
d2 y
+ 9y = −2t − 2et (∗∗)
dt2

8
 la ecuación característica: r2 + 9 = 0 ⇒ r = ±3i
solución homogénea: yh = c1 cos(3t) + c2 sen(3t)
para encontrar la solución particular, suponemos:
yp = Aet + Bt + C
derivando:
dy d2 y
= Aet + B, = Aet
dt dt2
reemplazando en la ec. diferencial (∗∗):
Aet + 9(Aet + Bt + C) = −2t − 2et
igualando coeficientes:
1
et : A + 9A = −2 ⇒ A=−
5
2
t: 9B = −2 ⇒ B=−
9
t0 : 9C = 0 ⇒ C=0
1 2
por lo tanto: yp = − et − t
5 9
1 2
y entonces, y(t) = c1 cos(3t) + c2 sen(3t) − et − t
5 9
Entonces, la solución del sistema de ecuaciones diferenciales es:
 3 1 2 2
 x(t) = c1 cos(3t) + c2 sen(3t) + et − t2 + t −

5 3 9 27
1
 y(t) = c cos(3t) + c sen(3t) − et − t
 2
1 2
5 9

EJERCICIOS PROPUESTOS:
 2
 d x + 3 dy + 3y = 0

dt2 dt

Resolver: 2
 d x + 3y = te−t


dt2