Apuntes de Ecuaciones Diferenciales UNIDAD 3

Unidad 3: ED Lineales de n-ésimo Orden con Coeficientes Constantes
UNIDAD 3
ECUACIONES DIFERENCIALES LINEALES DE n-ESIMO

ORDEN CON COEFICIENTES CONSTANTES
INTRODUCCIÓN
3.1 Solución de la Ecuaciones Diferenciales Homogéneas con

Coeficientes Constantes
3.2 Solución de la Ecuaciones Diferenciales No Homogéneas

con Coeficientes Constantes
3.2.1 Método de los Coeficientes Indeterminados

3.2.2 Método de Variación de Parámetros
3.3 Aplicaciones
3.3.1 Circuito Eléctrico RLC Serie

3.3.2 Sistemas Masa-Resorte
3.4 Sección de Problemas
3.4.1 Auto Evaluación

3.4.2 Solución de la Auto Evaluación
Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

62
Objetivos particulares de la unidad.- El alumno resolverá Ecuaciones Diferenciales de n-

ésimo Orden por diversos métodos.
INTRODUCCIÓN
En el presente capítulo se estudiarán los métodos de solución para las Ecuaciones

Diferenciales de n-ésimo orden, primordialmente enfocados a las de segundo orden. Se
inicia por la solución de ecuaciones diferenciales (Homogéneas y No Homogéneas) de
Coeficientes Constantes. Posteriormente, para las ecuaciones diferenciales no homogéneas,
se emplean los métodos de solución por variación de Parámetros y por Coeficientes
Indeterminados; finalmente se aplican estos métodos de solución a problemas de circuitos
eléctricos de tipo RLC (resistencia-inductancia-capacitancia) conectados en serie a una
fuente de alimentación de corriente alterna o corriente directa. También, al igual que en las
unidades anteriores, se incluyen problemas resueltos en forma detallada y una auto
evaluación que consiste en una serie de problemas con su respectiva solución.
3.1 Solución de la Ecuaciones Diferenciales Homogéneas con

Coeficientes Constantes
Una ecuación diferencial de la forma
dny d n1 y d2y dy

an n
 a n 1 n 1
 ...  a 2 2
 a1  a0 y  f ( x)
dx dx dx dx
Donde los coeficientes a n , a n 1 ,..., a 2 , a1 , a 0 son todos constantes y por lo menos uno de
ellos diferente de cero, se clasifica como una ecuación diferencial de coeficientes
constantes; además, si f ( x)  0 se trata de una ecuación diferencial homogénea, en caso
contrario se clasifica como no homogénea. Si n  2 y f ( x)  0 se tiene una ecuación
diferencial de segundo orden con coeficientes constantes homogénea. Para este caso:
d2y dy
a 2 2  a1  a0 y  0
dx dx
Al multiplicarla por 1/a2 se obtiene la ecuación diferencial equivalente
d2y dy
2
  y  0 (1)
dx dx
a1 a0
Donde   y 
a2 a2

63
Ejemplo. Clasifique las siguientes ecuaciones diferenciales de acuerdo a
a) sus coeficientes (constantes o variables)

b) homogéneas o no homogéneas
c) su orden (primer orden, segundo orden, tercer orden, etc.)
1. y   4 y   3 y  10
2. y  y  0
3. y   3xy   2  0
4. x 2 y   4 xy   4 y  8 x 2
5. 5 y iv  10 y   3  xy
Para el problema 1 se tiene que todos los coeficientes (los términos que multiplican a la
variable dependiente y sus derivadas) son constantes, por lo que se clasifica como ecuación
diferencial de coeficientes constantes; debido a que se encuentra igualada a una expresión
diferente de cero, se trata de una ED no homogénea. Finalmente, se tiene una ED de
segundo orden.
Para el problema 2 los coeficientes, al igual que el ejemplo anterior son todos ellos
constantes; Como f ( x)  0 , la ED es homogénea. De acuerdo a la mayor derivada que
involucra la ED, ésta es de tercer orden.
Para el problema 3 el coeficiente que afecta a la primera derivada no es constante, por lo

tanto es una ED de coeficientes variables; Aparentemente se tiene una ED homogénea, sin
embrago debe notarse que el termino 2 no está multiplicando a la variable dependiente o a
alguna de sus derivadas, por lo tanto deberá escribirse al lado derecho de la igualdad y se
trata de una ED no Homogénea; ED de segundo orden.
Para el problema 4 por lo menos uno de los coeficientes no es constante, en este caso son
dos, el que multiplica a la segunda derivada y el que corresponde a la primera derivada por
lo que se trata de una ED de coeficientes variables; Nótese que el término 8x2 no afecta a la
variable dependiente ni a ninguna de sus derivadas por lo que es correcto que se represente
al lado derecho de la igualdad y con esto, la ED es no homogénea; La ED es de segundo
orden.
Para el problema 5 todos los coeficientes son constantes excepto el que multiplica a y (el
término al lado derecho de la igualdad que en realidad se debe representar del lado
contrario) se tiene una ED de coeficientes variables; La ED es no homogénea debido a que
al lado izquierdo de la igualdad hay un término que no está multiplicando a la variable
dependiente o a alguna de sus derivadas y por lo tanto debe escribirse al lado derecho con
lo que la función f (x) es diferente de cero; la ED es de cuarto orden.

64
La solución de la ecuación diferencial (1) es una función y que cumpla con la igualdad que
se indica. Si se supone que la solución tiene la forma y  ce rx (c es cualquier constante, se
considerará como caso particular c = 1) entonces, para determinar el valor de la constante r
se procede a sustituir en la ecuación diferencial, con lo que se obtiene:
d 2 (e rx ) d (e rx )
2
    (e rx )  0
dx dx
Desarrollando las derivadas
r 2 e rx  re rx  e rx  0
Factorizando el término r2
e rx (r 2  r   )  0
Por lo que:
r 2  r    0 (2)
Otra forma de encontrar la ecuación característica es mediante la siguiente sustitución:
y 1
y  r
y   r 2
 
yn  rn
Es una ecuación de segundo grado, se le conoce como ecuación característica, ecuación

auxiliar o también como polinomio característica de la ED; para su solución se puede
aplicar factorización o la fórmula general, para lo cual a  1, b   , c   . De acuerdo a
esto se tiene que:
    2  4
r1, 2 
2
Al realizar las operaciones aritméticas al interior de la raíz se tienen tres posibles casos:
Caso I: Las raíces son reales y diferentes (  2  4  0 )
Caso II: Las raíces son reales e iguales (  2  4  0 )
Caso III: Las raíces son complejas y conjugadas (  2  4  0 )

65
    2  4     2  4
Caso I:  2  4  0 por lo que las raíces r1  y r2 
2 2
son reales y diferentes. Al existir dos valores para r se tienen las soluciones:
y1  e r1x y y 2  e r2 x
O de manera más general, para cualquier constante c
y1  c1e r1x y y 2  c2 e r2 x
De acuerdo al Principio de Superposición:
y h  c1e r1 x  c 2 e r2 x
 
Caso II:  2  4  0 por lo que las raíces r1  y r2  son reales e iguales. Al ser
2 2
iguales las raíces, se tiene una solamente una solución y se puede, a partir de esta,
determinar una segunda solución para la ED, con lo cual se determina que las soluciones
son:
y1  e rx y y 2  xe rx
y1  c1e rx y y 2  c 2 xe rx
y h  y1  y 2
Y la solución de la ED homogénea:
y h  c1e rx  c 2 xe rx
O también
 
 x  x
y h  c1e 2
 c2 xe 2

66
    2  4     2  4
Caso III:  2  4  0 por lo que las raíces r1  y r2 
2 2
son complejas y conjugadas. Las raíces al ser complejas poseen parte real y parte
imaginaria y se pueden representar como:
r1  a  bi
r2  a  bi
Donde:

a
2
 2  4
b
2
Las soluciones son:
y1  e r1x y y 2  e r2 x
y1  c1e r1x y y 2  c2 e r2 x
y h  c1e r1 x  c 2 e r2 x
De esta última solución al sustituir r1 y r2 se tiene:
y h  c1e ( a bi ) x  c 2 e ( a bi ) x
Eliminando los paréntesis
y h  c1e ax bxi  c 2 e ax bxi
Aplicando propiedades de la función exponencial
y h  c1e ax e bxi  c 2 e ax e bxi

67
Por la fórmula de Euler
e bxi  cos bx  isenbx
Sustituyendo en la solución homogénea
y h  c1e ax (cosbx  isenbx )  c 2 e ax (cosbx  isenbx )
Desarrollando
y h  c1e ax cos bx  c1e ax isenbx  c 2 e ax cos bx  c 2 e ax isenbx
Agrupando parte real y parte imaginaria
y h  (c1  c 2 )e ax cos bx  i (c1  c 2 )e ax senbx
Si C1  c1  c2 y Si C2  i(c1  c2 ) , se llega a la solución equivalente
y h  C1e ax cos bx  C 2 e ax senbx
Finalmente, factorizando la función exponencial
y h  e ax C1 cos bx  C 2 senbx 
Es la solución de la ED para el caso III, no olvidar que

a
2
 2  4
b
2
Ejemplo. En los problemas 1–5 determine la solución de cada una de las ecuaciones
diferenciales que se indican.
1. y   4 y   3 y  10 ; sujeta a las condiciones y(0)  2, y(0)  1
2. y  9 y  0 ; sujeta a la condición inicial y( )  2 , y (0)  1
3. y   3 y   2  0
4. y  4 y  4 y  0 ; sujeta a la condición inicial y(0)  1, y(0)  4
5. y   y   2  0

68
1. y  4 y  3 y  0 ; sujeta a las

condiciones y(0)  2, y(0)  1 2  c1  c2
Se obtiene la ecuación característica c1  c2  2 (1)

sustituyendo
La segunda condición indica que y  1
y 1 cuando x = 0. Para aplicar estas
y  r condiciones se requiere la derivada de la
solución; así, ésta derivada es
y   r 2
y h  3c1e 3 x  c 2 e x
r  4r  3  0
2
Sustituyendo las condiciones

(r  3)(r  1)  0
1  3c1e 3( 0)  c2 e ( 0)
Los valores que cumplen con la
afirmación anterior (se puede factorizar
de manera inmediata o en su caso aplicar 1  3c1  c2
la formula general) son
3c1  c2  1 (2)
r1  3
Por lo tanto se tiene un sistema de 2
r2  1 ecuaciones con dos incógnitas
De acuerdo a las raíces obtenidas la c1  c2  2 (1)

solución corresponde al caso 1: Raíces
reales y distintas. Por lo tanto la solución 3c1  c2  1 (2)
homogénea o complementaria se expresa
de la forma La solución se puede obtener por diversos
métodos como determinantes, suma y
y h  c1e r1 x  c 2 e r2 x resta, sustitución, entre otros, los valores
para c1 y c2 son
y h  c1e 3 x  c 2 e x
c1   3
2
Para comprobar la solución simplemente
se sustituye ésta en la ED y se debe de c2  7
cumplir la igualdad. 2
De acuerdo a la primer condición inicial Por lo tanto la solución de la ED sujeta a

se tiene que y = 2 cuando x = 0; las condiciones indicadas es
sustituyendo en la solución
3 7
yh   e3x  e x
2  c1e 3( 0)  c2 e ( 0) 2 2

69
2. y  9 y  0 ; sujeta a la condición 3. y   3 y   2  0

inicial y( )  2 , y (0)  1
Se obtiene la ecuación característica
Ecuación característica sustituyendo
r2  9  0 y 1
y  r
Resolviendo se obtiene
y   r 2
r1  3i  0  3i y r2  3i  0  3i
Por lo tanto
De acuerdo a las raíces obtenidas la
r 2  3r  2  0
solución corresponde al caso 3: Raíces
complejas y conjugadas. Por lo tanto la
Ya que no es posible factorizar de forma
solución homogénea tiene la forma
directa se resuelve aplicando la formula
general y se llega a
y h  e ax C1 cos bx  C 2 senbx 
 3  17
y h  e ( 0) x C1 cos 3x  C 2 sen 3x  r1 
2
y h  C1 cos 3x  C 2 sen 3x  3  17
r2 
2
Aplicando la p condición y  2 , x  0
Con lo cual las raíces que se obtienen
2  C1 cos 3(0)  C2 sen3(0) corresponden al caso I: Raíces reales y
diferentes. La solución para este caso se
2  C1 expresa de la forma:
Para determinar al valor de la otra y h  c1e r1 x  c 2 e r2 x

constante se sustituye la segunda
condición y  1 , x  0 ; para esto se Sustituyendo los valores de las raíces
obtiene la derivada de la solución
 3 17   3 17 
 x  x
   
y h  c1e  
 c2 e  
2 2
y   3C1 sen3x  3C2 cos 3x
1  3C1 sen3(0)  3C2 cos 3(0) O también
y h  c1e 0.562 x  c 2 e 3.562 x

1  3C 2  C 2  1
3
No se determinan los valores para las
Sustituyendo la solución es constantes, ya que este problema carece
de condiciones iniciales.
1
y h  2 cos 3 x  sen 3 x
3
70
4. y  4 y  4 y  0 ; sujeta a la condición 5. y   y   2  0

inicial y(0)  1, y(0)  4
Ecuación característica sustituyendo
r 2  4r  4  0 y 1
y  r
r1  r2  2 y   r 2
Se tiene el caso II: Raíces reales e iguales. Por lo tanto

La solución es de la forma
r2  r  2  0
y h  c1e rx  c 2 xe rx
Ya que no es posible factorizar de forma
y h  c1e 2 x  c 2 xe 2 x directa se resuelve aplicando la formula
general y se llega a
Para determinar el valor de las constantes,
se aplica la primera condición y  1 ,  1  7i 1 7
r1    i
x0 2 2 2
1  c1e 2( 0)  c2 (0)e 2( 0)
 1  7i 1 7
r2    i
1  c1 2 2 2
Para determinar el valor de la otra Con lo cual las raíces que se obtienen
constante, se aplica la segunda condición corresponden al caso I: Raíces reales y
y  4 , x  2 . Para esto se deriva la diferentes. La solución para este caso se
solución expresa de la forma:
y h  2c1e 2 x  c 2 (2 xe 2 x  e 2 x ) y h  e ax C1 cos bx  C 2 senbx 
Sustituyendo las condiciones iniciales

1
 x  7 7 
yh  e 2
C1 cos x  C2 sen x
 2 2 
4  2c1e 2( 0)  c2 [2(0)e 2( 0)  e 2( 0) ]
No se determinan los valores para las
4  2c1  c2 [0  1] constantes, ya que este problema carece
de condiciones iniciales.
c2  4  2c1  c2  6
Sustituyendo se llega a la solución
y h  e 2 x  6 xe 2 x

71
3.2 Solución de la Ecuaciones Diferenciales No Homogéneas

con Coeficientes Constantes
De acuerdo a la sección anterior, una ecuación diferencial de la forma
dny d n1 y d2y dy

an n
 a n 1 n 1
 ...  a 2 2
 a1  a0 y  f ( x)
dx dx dx dx
Donde los coeficientes a n , a n 1 ,..., a 2 , a1 , a 0 son todos constantes y por lo menos uno de
ellos diferente de cero, se clasifica como una ecuación diferencial de coeficientes
constantes; además, si f ( x)  0 se trata de una ecuación diferencial no homogénea. Si
n  2 y f ( x)  0 se tiene una ecuación diferencial de segundo orden con coeficientes
constantes no homogénea. Para este caso:
d2y dy
a2 2
 a1  a 0 y  f ( x)
dx dx
Al multiplicarla por 1/a2 se obtiene la ecuación diferencial equivalente
d2y dy
2
  y  f ( x) (1)
dx dx
a1 a0
Donde   y 
a2 a2
La solución general de (1) es la suma de dos soluciones: solución homogénea ( y h ) y

solución particular ( y p ) , es decir
y g  yh  y p
La solución homogénea se obtiene haciendo f ( x)  0 (ver sección anterior). Para la

solución particular existen dos métodos de solución:
1) Coeficientes Indeterminados y,
2) Variación de Parámetros
El primero de ellos es muy limitado. Sólo se puede aplicar para ciertas funciones, mientras
que el segundo es más general aunque más laborioso que el primero. El primero de ellos se
resuelve aplicando un poco de álgebra, mientras que en el segundo se involucran
determinantes, derivadas, integrales y un poco más de procedimiento algebraico.

72
3.2.1 Método de los Coeficientes Indeterminados
Este método consiste en proponer una solución particular yp de acuerdo a la forma que
tenga la función f(x) de la ecuación diferencial. Es importante indicar que aplica sólo aplica
para funciones polinomiales de grado n, exponenciales, seno y coseno, así como
combinaciones de las anteriores a través de las operaciones suma, diferencia, producto y
cociente siempre y cuando el divisor sea una función exponencial. En la siguiente tabla se
muestran algunas posibles expresiones para f(x) y la solución particular que se debe
proponer para yp.
f (x) yp
2 A
4x  1 A  Bx
 x2  9 A  Bx  Cx 2
 3e 4 x Ae4 x
e4x Ae 4 x
5 sen 2 x Asen 2 x  B cos 2 x
 3 cos 2 x Asen 2 x  B cos 2 x
2 xe 5 x ( A  Bx)e 5 x
( x 2  1)e 5 x ( A  Bx  cx 2 )e 5 x
( x  4)senx ( A  Bx)senx  (C  Dx) cos x
e  x senx e  x ( Asenx  B cos x)
e  x  senx Ae  x  Bsenx  C cos x)
5x 2  sen 2 x  3e 3 x A  Bx  Cx 2  Dsen 2x  F cos 2x  Ge 3x

5 cos 3 x  sen 2 x Asen3 x  B cos 3 x  Csen 2 x  D cos 2 x
5  ex A  Be x

73
En la primera fila se observa que la función es una constante, por lo que la solución
particular que se propone es un polinomio de grado cero; hay que recordar que un
polinomio de grado n se puede representar de la forma:
A0  A1 x  A2 x 2  ...  An  2 x n  2  An 1 x n 1  An x n
Por lo tanto un polinomio de grado cero es una constante que se puede representar por
cualquier letra. En los renglones 2 y 3 son casos similares; son polinomios de primero y
segundo grado respectivamente.
Para la función del siguiente renglón se tiene una constante por una función exponencial, o
visto de otro modo un polinomio de grado 0 por la exponencial; por lo tanto la solución
particular que se propone tiene la forma de un polinomio de grado cero (una constante) por
la función exponencial.
En las filas 8 y 9 tenemos que f(x) es el producto de dos funciones; para el renglón 8 una
función cuadrática y una exponencial, por lo tanto la solución que se propone tiene la forma
de un polinomio de segundo grado por una función exponencial. Para el renglón 9 debe de
ser obvio porque la solución particular se representa como el producto de un polinomio de
primer grado y una función exponencial.
Cabe mencionar que si la función tiene está representada por la función seno o coseno, la
solución particular tiene la forma de una suma de una constante por la función seno y una
constante diferente por la función coseno. Todos los demás ejemplos que se muestran son
combinaciones de los anteriores.
Una vez que se propone la solución particular, se procede a sustituirla en la ecuación

diferencial no homogénea para calcular los coeficientes de la solución propuesta.
Finalmente, la solución de la ecuación diferencial no homogénea es la suma de las
soluciones homogénea y particular.
En los problemas 6–10 determine la solución de la ecuación diferencial de coeficientes

constantes no homogénea que se indica, por el método de Coeficientes
Indeterminados.
6. y   y   6 y  2e 7 x
7. y   5 y   6 y  5  x 2 ; con y(0)  2, y(0)  1
8. y  4 y  cos 4x
9. y  4 y  10 ; con y(0)  1, y(0)  4
10. y   7 y   6 y  3e 6 x

74
6. y   y   6 y  2e 7 x Para determinar el coeficiente A la

solución particular se sustituye en la ED
Para determinar la solución homogénea, no homogénea
se iguala a cero la ecuación diferencial
y   y   6 y  2e 7 x (4)
y  y  6 y  0
Si
sustituyendo y p  Ae7 x , entonces
y 1
yp  7 Ae7 x
y  r
y   r 2 yp  49 Ae7 x
Por lo tanto Sustituyendo en (4)
r2  r  6  0 49 Ae7 x  7 Ae7 x  6 Ae7 x  2e 7 x

(r  3)(r  2)  0 Sumando término semejantes
r1  3, r2  2 36 Ae7 x  2e 7 x
La solución que se obtiene corresponde al Para que la igualdad se cumpla se

caso I: Raíces reales y diferentes requiere que el coeficiente de la función
exponencial del lado izquierdo sea
y h  c1e 3 x  c 2 e 2 x (1) idéntico al coeficiente del lado derecho de
la ecuación; con esto, se llega a la
ecuación
Para determinar la solución particular ésta
se debe de proponer de acuerdo a la
36 A  2
forma que tenga f (x) ; ya que f (x) tiene
la forma de un polinomio de grado cero Por lo tanto
(una constante) por una función
exponencial, se propone la solución 1
particular A
18
y p  Ae7 x (2) Sustituyendo este valor en (3)
Por lo tanto la solución general de la ED 1 7x

no homogénea es la suma de (1) y (2) y g  c1e 3 x  c 2 e  2 x  e
18
y g  yh  y p
y g  c1e3x  c2 e 2 x  Ae7 x (3)

75
7. y   5 y   6 y  5  x 2 ; con y(0)  2, Sustituyendo en (4)

y (0)  1 2C  5( B  2Cx )  6( A  Bx  Cx 2 )  5  x 2
Primero, se identifica f ( x)  5  x 2 ; si se Eliminando paréntesis y factorizando los

hace f ( x)  0 se tiene una ED términos que multiplican a x2, x y x0 se
homogénea llega a
y  5 y  6 y  0  6Cx 2  (10C  6 B) x  (2C  5B  6 A)  5  x 2
Resolviendo (ver sección anterior) se Para que la igualdad se cumpla se

obtiene la solución homogénea: requiere que los coeficientes de x2, x y x0
del lado izquierdo sean idénticos a los del
y h  c1e 6 x  c 2 e x (1) lado derecho de la ecuación; con esto, se
llega al siguiente sistema de ecuaciones
se debe de proponer de acuerdo a la  6C  1
forma que tenga f (x) ; ya que f (x) tiene
10C  6 B  0
la forma de un polinomio de segundo
grado, se propone la solución particular
2C  5 B  6 A  5
y p  A  Bx  Cx 2 (2) Resolviendo se tiene
Por lo tanto la solución general de la ED  59 5 1

no homogénea es la suma de (1) y (2) A , B , C
108 18 6
y g  yh  y p Sustituyendo los valores en (3)
y g  c1e 6 x  c2 e x  A  Bx  Cx 2 (3) y g  c1e 6 x  c 2 e x 

59 5 1
 x  x2
108 18 6
Para determinar los coeficientes A, B y C
la solución particular se sustituye en la Para conocer los valores de las constantes
ED no homogénea c1 y c2 se sustituye y = 2, x = 0 en la
expresión anterior
y   5 y   6 y  5  x 2 (4)
59 5 1
2  c1e 6 ( 0 )  c 2 e ( 0 )   ( 0 )  ( 0) 2
Si 108 18 6
Se tiene
y p  A  Bx  Cx 2 , entonces
59
y p  B  2Cx 2  c1  c 2 
108
y p  2C

76
Es decir, 8. y  4 y  cos 4x
275 Para determinar la solución homogénea,

c1  c 2 
108 se iguala a cero la ecuación diferencial
Para obtener una segunda ecuación, y así y  4 y  0

un sistema de dos ecuaciones con dos
incógnitas, se procede a aplicar la otra Se obtiene la ecuación característica
condición inicial, para esto se deriva la sustituyendo
solución yg y se sustituye y  1 , x = 0
y 1
y g  6c1e 6 x
5 2
 c2 e x   x y  r
18 6
y   r 2
5 2
1  6c1e 6 ( 0 )  c 2 e ( 0 )   ( 0) Por lo tanto
18 6
r2  4  0
5
1  6c1  c 2 
18 Factorizando
 6c1  c 2 
13 (r  2)(r  2)  0
18
Por lo tanto los valores de las raíces son
La solución para las constantes es
r1  2, r2  2
197
c1 
756
La solución que se obtiene corresponde al
caso I: Raíces reales y diferentes
16
c2 
7 y h  c1e 2 x  c 2 e 2 x (1)
Por lo que la solución queda expresada
como
se debe de proponer de acuerdo a la
forma que tenga f (x) ; ya que f (x) es
197 6 x 16 x 59 5 1
yg  e  e   x  x2 una función cosenoidal de argumento 4x,
756 7 108 18 6
se propone la solución particular
y p  Asen 4 x  B cos 4 x (2)
Por lo tanto la solución general de la ED

no homogénea es la suma de (1) y (2)
y g  yh  y p

77
y g  c1e3x  c2 e 2 x  Asen4x  B cos 4x 9. y  4 y  10 ; con y(0)  1 y(0)  4
Para determinar los coeficientes A y B la Para determinar la solución homogénea,

solución particular se sustituye en la ED se iguala a cero la ecuación diferencial,
no homogénea obteniendo así una ED homogénea
y  4 y  cos 4x (4) y  4 y  0
Si Se obtiene la ecuación característica

sustituyendo
y p  Asen 4 x  B cos 4 x , entonces
y 1
y p  4 A cos 4 x  4 Bsen 4 x y  r
y p  16 Asen 4 x  16 B cos 4 x y   r 2
Por lo tanto
Sustituyendo en (4)
 16 Asen 4 x  16 B cos 4 x r2  4  0
 4( Asen4 x  B cos 4 x)  cos 4 x
Aplicando la fórmula general se obtienen
los valores de las raíces, estos son
Eliminando paréntesis y factorizando los
coeficientes que multiplican a la fucnión
seno y a la función coseno r1  2i, r2  2i
 20 Asen 4 x  20 B cos 4 x  cos 4 x De acuerdo a la solución que se obtiene

las raíces corresponden al caso III: Raíces
Para que la igualdad se cumpla se complejas y conjugadas
requiere que el coeficiente de la función
seno del lado izquierdo sea idéntico al y h  C1 cos 2 x  C 2 sen 2 x (1)
coeficiente del lado derecho de la
ecuación; de manera similar para la Para determinar la solución particular ésta
función coseno; con esto, se llega a se debe de proponer de acuerdo a la
forma que tenga f (x) ; ya que f (x) es
 20 A  0  A  0 una constante, se propone un polinomio
de grado cero
1
 20 B  1  B  
20 yp  A (2)
Sustituyendo este valor en

y g  c1e3x  c2 e 2 x  Asen4x  B cos 4x
y g  yh  y p
1
y g  c1e 3 x  c 2 e  2 x  cos 4 x
20

78
y g  C1 cos 2 x  C 2 sen 2 x  A (3) 5

 1  C1 
2
Para calcular el valor del coeficiente A se
sustituye la solución particular en la 7
C1  
ecuación diferencial no homogénea 2
Si Ahora se deriva la solución general y se

sustituyen los valores de la segunda
y p  A , entonces condición y   4 , x  0
y p  0 y g  2C1 sen 2 x  2C 2 cos 2 x
y p  0 4  2C1 sen 2(0)  2C2 cos 2(0)
Sustituyendo en 4  2C2
y  4 y  10 C2  2
0  4 A  10 Sustituyendo los valores de las constantes

en (4)
Despejando el valor del coeficiente
7 5
10 5 y g   cos 2 x  2 sen 2 x 
A  2 2
4 2
Sustituyendo en (3)
5
y g  C1 cos 2 x  C 2 sen 2 x  (4)
2
Ahora se aplica la primera condición

inicial para calcular los valores de las
constantes involucradas en la solución
general. La primera condición indica que
y  1 y x  0 , de acuerdo a lo anterior
se obtiene
5
y g  C1 cos 2 x  C 2 sen 2 x 
2
5
 1  C1 cos 2(0)  C 2 sen 2(0) 
2

79
10. y   7 y   6 y  3e 6 x Sustituyendo en la ED no homogénea se

llega a una inconsistencia
Para determinar la solución homogénea,
se iguala a cero la ecuación diferencial, 0  3e 6 x
obteniendo así una ED homogénea
La situación que se presenta es que la
y  7 y  6 y  0 solución que se propuso para y p ya
existe como solución de la ED
Se obtiene la ecuación característica homogénea; para evitar la
sustituyendo indeterminación la solución que se
proponga se multiplica por la variable
r 2  7r  6  0 independiente, en este caso x, quedando
y p  Axe 6 x , por lo tanto la solución
r1  6, r2  1 general será
De acuerdo a la solución que se obtiene y g  c1e 6 x  c2 e  x  Axe6 x (4)

las raíces corresponden al caso I: Raíces
reales y diferentes Si
y p  Axe 6 x , entonces
y h  c1e 6 x  c 2 e  x (1)
La solución particular que se propone es yp  6 Axe 6 x  Ae6 x
y p  Ae6 x (2) yp  36 Axe 6 x  12 Ae6 x
Por lo tanto la solución general de la ED Sustituyendo en la ED no homogénea

y   7 y   6 y  3e 6 x
y g  yh  y p
36 Axe 6 x  12 Ae 6 x  7(6 Axe 6 x  Ae 6 x )
6 x x 6 x
y g  c1e  c2 e  Ae (3)  6( Axe 6 x )  3e 6 x
Procediendo de manera similar a los Multiplicando para eliminar paréntesis y

ejemplos anteriores para determinar el sumando términos semejantes
valor del coeficiente A
Si 3
 5 Ae 6 x  3e 6 x  A  
5
y p  Ae6 x , entonces
Sustituyendo en (4)
yp  6 Ae 6 x
3 6 x
y g  c1e 6 x  c 2 e  x  xe
5
yp  36 Ae6 x

80
3.2.2 Método de Variación de Parámetros
Este método, al igual que el de coeficientes indeterminados, se emplea para hallar la

solución particular cuando se tiene una ecuación diferencial homogénea. A partir de la ED
de 2º orden lineal no homogénea y de coeficientes constantes, es decir
y  y  y  f (x) (1)
Se propone la solución particular
y p  v1 y1  v 2 y 2
Donde y1 y y 2 son soluciones de la ecuación diferencial homogénea
y  y  y  0 (2)
Si
y p  v1 y1  v 2 y 2 , entonces,
y p  v1 y1  v1 y1  v 2 y 2  v 2 y 2
Haciendo
v1 y1  v2 y 2  0 (A)
Se tiene una primera ecuación, por lo tanto
y p  v1 y1  v 2 y 2
y p  v1 y1  v1 y1  v 2 y 2  v 2 y 2
Sustituyendo y p , y p y y p en (1)
(v1 y1  v1 y1  v2 y 2  v2 y 2 )   (v1 y1  v1 y1  v2 y 2  v2 y 2 )   (v1 y1  v2 y 2 )  f ( x)
Multiplicando para eliminar los paréntesis
v1 y1  v1 y1  v2 y 2  v2 y 2  v1 y1  v1 y1  v2 y2  v2 y 2  v1 y1  v2 y 2  f ( x)

81
Factorizando v1 y v2 , se obtiene
v1[ y1  y1  y1 ]  v2 [ y 2  y 2  y 2 ]  v1 y1  v2 y 2  f ( x)
Debido a que y1 y y 2 son soluciones de la ecuación diferencial homogénea se tiene que
y1  y1  y1  0
y2  y2  y 2  0
Por lo tanto
v1 y1  v2 y 2  f ( x) (B)
Se llega a un sistema de dos ecuaciones con dos incógnitas
v1 y1  v2 y 2  0 (A)
v1 y1  v2 y 2  f ( x) (B)
Resolviendo por determinantes
0 y2 0 y2
f ( x) y 2 f ( x) y 2
v1   v1   dx
y1 y 2 y1 y 2
y1 y 2 y1 y 2
y1 0 y1 0
y  f ( x) y1 f ( x)
v 2  1  v2   dx
y1 y 2 y1 y 2
y1 y 2 y1 y 2
Y la solución particular
y p  v1 y1  v 2 y 2
El determinante del divisor se le llama wronskiano del sistema y se representa como
w( y1 , y 2 )

82
En los problemas 11-15 resuelva la ecuación diferencial de coeficientes constantes no

homogénea empleando el método de Variación de Parámetros.
11. y  7 y  10 y  2x  1 Como el Wronskiano es diferente de cero,

implica que el sistema es linealmente
Primero se determina la solución de la ED independiente
homogénea, haciendo f ( x)  0
Resolviendo para v1
y  7 y  10 y  0
0 e 5 x
La ecuación característica es (2 x  1)  5e 5 x
v1 ' ( x) 
 3e 7 x
r  7r  10  0
2
 (2 x  1)e 5 x
Esto implica que v1 ' ( x) 
 3e 7 x
r1  2, r2  5
1
v1 ' ( x )  ( 2 x  1)e 2 x
3
Las raíces son reales y diferentes, por lo
tanto la solución homogénea es 1
v1 ( x) 
3  (2 x  1)e 2 x dx
y h  c1e 2 x  c 2 e 5 x (1)
1 2x
Donde v1 ( x)  xe
3
y1  e 2 x , y2  e 5 x Resolviendo para v2
Ahora, se procede a resolver la solución e 2 x 0

particular; primeramente el Wronskiano
 2e  2 x (2 x  1)
del sistema está dado por v 2 ' ( x) 
 3e 7 x
y1 y2
w( y1 , y 2 )  (2 x  1)e 2 x
y1 y 2 v 2 ' ( x) 
 3e 7 x
e 2 x e 5 x 1
w( y1 , y 2 )  v 2 ' ( x )   (2 x  1)e 5 x
 2e 2 x  5e 5 x 3
w( y1 , y 2 )  5e 7 x  2e 7 x 1
v 2 ( x)  
3  (2 x  1)e 5 x dx
w( y1 , y 2 )  3e 7 x
e5x  3
v 2 ( x)   2 x  5 
15

83
Sustituyendo los valores en yp resulta: 12. y  9 y  sen3x
y p  v1 y1  v 2 y 2 Primero se determina la solución de la ED

1 2 x 2 x e 5 x  3  5 x
y p  xe (e )  2 x  5  e y  9 y  0
3 15
La ecuación característica es
1 1  3
y p  x  2 x  
3 15  5 r2  9  0
Desarrollando Esto implica que
yp 
1 2
x x
3 r1  3i, r2  3i
3 15 75
Las raíces son complejas y conjugadas,
Simplificando por lo tanto la solución homogénea es
1 3 y h  C1 sen 3x  C 2 cos 3x (1)
yp  x (2)
5 75
Donde
Finalmente, recordar que la solución de la
ecuación diferencial no homogénea es la
suma de (1) y (2) y1  sen3x, y2  cos3x
y g  yh  y p Ahora, se procede a resolver la solución

del sistema está dado por
Es decir
1 3 y1 y2
y g  c1e  2 x  c 2 e 5 x  x w( y1 , y 2 ) 
5 75 y1 y 2
sen 3x cos 3x
w( y1 , y 2 ) 
3 cos 3x  3sen 3 x
w( y1 , y 2 )  3sen 2 3x  3 cos2 3x
w( y1 , y 2 )  3
Como el Wronskiano es diferente de cero,

independiente

84
Resolviendo para v1 Sustituyendo los resultados en yp
0 cos 3x y p  v1 y1  v 2 y 2
sen 3x  3sen 3x
v1 ' ( x)  1 
3 yp 
1 1
sen 2 3xsen 3x   x  sen 6 x  cos 3x
18 6 36 
 sen 3 x cos 3 x
v1 ' ( x )  1 1 1
3 yp  sen 3 3 x  x cos 3 x  sen 6 x cos 3 x
18 6 36
1
v1 ' ( x)  sen 3 x cos 3 x Descomponiendo
3
sen 6 x  2 sen 3 x cos 3 x
1
v1 ( x)   sen 3 x cos 3 xdx
3 1 1 1
yp  sen 3 3 x  x cos 3 x  (2 sen 3 x cos 3 x) cos 3 x
18 6 36
1
v1 ( x )  sen 2 3 x 1 1 1
18 yp  sen 3 3 x  x cos 3 x  sen 3 x cos 2 3 x
18 6 18
Resolviendo para v2
Aplicando la identidad
sen 3x 0
cos2 3x  1  sen 2 3x
3 cos 3x sen 3x
v 2 ' ( x) 
3 Simplificando
sen 2 3x yp 
1 1
sen 3 3 x  x cos 3 x 
1
sen 3 x(1  sen 2 3 x)
v 2 ' ( x) 
3 18 6 18
1 1 1 1
1 yp  sen 3 3 x  x cos 3 x  sen 3 x  sen 3 3 x
v 2 ' ( x)   sen 2 3 x 18 6 18 18
3
1 1 (2)
yp   x cos 3 x  sen 3 x
1 6 18
v 2 ( x)    sen 2 3 xdx
3
Empleando identidades trigonométricas ecuación diferencial no homogénea es la
suma de (1) y (2)
1 1 1 
v 2 ( x)      cos 6 x dx
3 2 2 
y g  yh  y p
Es decir
1 1 1 
v2 ( x)    x  sen 6 x  y g  C1 sen 3 x  C 2 cos 3 x 
1 1
sen 3 x  x cos 3 x
3 2 12  18 6

85
1 Resolviendo para v1
13. y   3 y   2 y 
1 ex
0 ex
Primero se determina la solución de la ED
homogénea, haciendo f ( x)  0 1
 ex
1 e x
y  3 y  2 y  0 v1 ' ( x) 
e 3 x
ex
v1 ' ( x)  1 3ex
x
r 2  3r  2  0
e
Esto implica que
ex
v1 ' ( x ) 
r1  2, r2  1 e 3 x (1  e x )
Las raíces son reales y diferentes, por lo e2x

tanto la solución homogénea es v1 ( x)   dx
1 ex
y h  c1e 2 x  c 2 e  x (1) Realizando la división
Donde ex
v1 ( x)   e x dx   dx
ex 1
y1  e 2 x , y2  e  x
v1 ( x)  e x  ln e x  1
Ahora, se procede a resolver la solución
Resolviendo para v2
e 2 x 0
y y2
w( y1 , y 2 )  1
y1 y 2
1
 2e  2 x
1 ex
e 2 x ex v 2 ' ( x) 
w( y1 , y 2 )  e 3 x
 2e 2 x  e x
e 2 x
w( y1 , y 2 )  e 3 x  2e 3 x v 2 ' ( x)  1 3ex
x
e
w( y1 , y 2 )  e 3 x
e 2 x
v 2 ' ( x) 
Como el Wronskiano es diferente de cero, e 3 x (1  e x )
independiente

86
ex 14. y   2 y   8 y  e 4 x csc 2 x
v 2 ' ( x) 
1 ex
Primero se determina la solución de la ED
ex
v 2 ( x)   dx y  2 y  8 y  0
1 ex
v2 ( x)  ln 1  e x
r 2  2r  8  0
Sustituyendo los resultados en yp
Esto implica que
y p  v1 y1  v 2 y 2 (2)
r1  4, r2  2
2 x x
y p  (e  ln e  1 )e
x x
 (ln 1  e )e
x
Las raíces son reales y diferentes, por lo

Finalmente, recordar que la solución de la tanto la solución homogénea es
suma de (1) y (2) y h  c1e 4 x  c 2 e 2 x (1)
y g  yh  y p Donde
Donde y1  e 4 x , y2  e 2 x
y h  c1e 2 x  c 2 e  x Ahora, se procede a resolver la solución

y p  (e x  ln e x  1 )e 2 x  (ln 1  e x )e  x del sistema está dado por
y1 y2
w( y1 , y 2 ) 
y1 y 2
e4x e 2 x
w( y1 , y 2 ) 
4e 4 x  2e 2 x
w( y1 , y 2 )  2e 2 x  4e 2 x
w( y1 , y 2 )  6e 2 x

independiente

87
Resolviendo para v1 15. y   8 y   12 y  3e 4 x  e x
0 e 2 x Para determinar la solución de la ED

e 4 x csc 2 x  2e  2 x homogénea, se hace f ( x)  0
v1 ' ( x) 
 6e 2 x
y   8 y  12 y  0
 e 2 x csc 2 x
v1 ' ( x)  Se obtiene la ecuación característica
 6e 2 x sustituyendo
v1 ' ( x) 
1
csc 2 x y 1
6 y  r
1 y   r 2
v1 ( x)   csc 2 xdx
6
Por lo tanto
1
v1 ( x)  ln csc 2 x  cot 2 x r 2  8r  12  0
12
Factorizando
Resolviendo para v2
(r  2)(r  6)  0
e4x 0
4e 4 x e 4 x csc 2 x Esto implica que
v 2 ' ( x) 
 6e 2x
r1  2, r2  6
8x
e csc2 x
v 2 ' ( x) 
 6e 2 x Las raíces son reales y diferentes, por lo
tanto la solución homogénea es
1
v 2 ' ( x)   e 6 x csx 2 x
6 y h  c1e 2 x  c 2 e 6 x (1)
1 6x Donde
6
v 2 ( x)   e csc 2 xdx
y1  e 2 x , y2  e 6 x
Ahora, se procede a resolver la solución
y p  v1 y1  v 2 y 2 (2) particular; primeramente el Wronskiano
1   1 
y g   ln csc 2 x  cot 2 x  e 4 x    e 6 x csc 2 xdx  e  2 x
12   6  y1 y2
w( y1 , y 2 ) 
y1 y 2
suma de (1) y (2)

88
e2x e6x 3 2 x 1
w( y1 , y 2 )  2 x v 2 ( x) 
4  e dx   e 5 x dx
4
2e 6e 6 x
3 1 5 x
w( y1 , y 2 )  6e 8 x  2e 8 x v 2 ( x)   e 2 x  e
8 20
w( y1 , y 2 )  4e 8 x
implica que el sistema es linealmente y p  v1 y1  v 2 y 2
independiente
 3 1   3 1 5 x  6 x
Resolviendo para v1 y p   e 2 x  e  x  e 2 x   e  2 x  e e
 8 4   8 20 
0 e6x Simplificando
3e 4 x  e x 6e 6 x
v1 ' ( x)  3 1 3 1 x
4e 8 x y p   e4x  e x  e4x  e
8 4 8 20
 3e10 x  e 7 x
v1 ' ( x)  6 4 x
4e 8 x y p   e4x  e (2)
8 20
3 1
v1 ' ( x)   e 2 x  e  x Finalmente, recordar que la solución de la
4 4 ecuación diferencial no homogénea es la
suma de (1) y (2)
3 2x 1
v1 ( x )  
4  e dx   e  x dx
4 y g  yh  y p
3 1
v1 ( x)   e 2 x  e  x Donde
8 4
y h  c1e 2 x  c 2 e 6 x
Resolviendo para v2
3 1
e2x 0 y p   e4x  e x
4 5
2e 2 x 3e  e x
4x
v 2 ' ( x) 
4e 8 x
3e 6 x  e 3 x
v 2 ' ( x) 
4e 8 x
3  2 x 1 5 x
v 2 ' ( x)  e  e
4 4

89
3.3 Aplicaciones
En esta sección se modelaran mediante ecuaciones diferenciales de coeficientes constantes

no homogéneas de segundo grado, dos fenómenos físicos: el comportamiento de la
corriente y carga eléctrica en un circuito eléctrico RLC conectado en serie a una fuente de
alimentación de voltaje y el movimiento de un cuerpo en un sistema masa-resorte. Para el
primer fenómeno se requiere conocer la ley de voltajes de Kircchoff y el comportamiento
del voltaje en cada elemento pasivo. Para el segundo es necesario el conocimiento e
interpretación de la ley de Hooke y de la segunda ley de Newton.
3.3.1 Circuito Eléctrico RLC Serie
Suponga que se tiene un circuito eléctrico RLC conectado en serie a una fuente de
alimentación. Según la ley de voltajes de Kircchoff, la suma de las caídas de voltaje en un
circuito eléctrico conectado en serie es igual al voltaje que se aplica a dicho circuito, es
decir
VL  VR  VC  VF (1)
Donde:
V L Es el voltaje o caída de tensión eléctrica en el elemento inductor
V R Es el voltaje o caída de tensión eléctrica en el elemento resistivo
V C Es el voltaje o caída de tensión eléctrica en el elemento capacitivo
V F Es el voltaje o tensión eléctrica suministrada por la fuente de alimentación
Además
di
VL  L
dt
VR  Ri
1
C
VC  idt
dq
i
dt

90
L, es la inductancia del elemento que se define como la capacidad de almacenar energía en

forma de campo magnético; se mide en Henrios (H)
R , es la resistencia eléctrica, que se define como la oposición de un elemento al paso de la

corriente eléctrica, se mide en ohms (  )
C, es la capacitancia, que se define como la capacidad de un elemento para almacenar

energía en forma de campo eléctrico, se mide en farads (F)
i , es la intensidad de corriente eléctrica, que se define como la variación de carga eléctrica

por unidad de tiempo, en amperes
q Es la carga eléctrica, en Coulombs (C)
La ecuación (1) de acuerdo a lo anterior, se puede expresar en función de la intensidad de

corriente eléctrica como:
di 1
L  Ri   idt  V F
dt C
O en término de la carga eléctrica:
d 2q dq 1
L R  q  VF (2)
dt dt C
Esta última expresión es una ecuación diferencial de segundo grado no homogénea. Los
coeficientes, L, R y C, son constantes. El método de solución se estudió en la sección 3.2
del presente capítulo. Se obtiene una solución homogénea haciendo el valor de la fuente de
alimentación igual con cero; posteriormente se determina la solución particular aplicando el
método de coeficientes determinados, o en su defecto, el método de variación de
parámetros. La solución del problema será, en primera instancia, determinar la carga
eléctrica en función del tiempo, para, posteriormente, encontrar la intensidad de corriente
eléctrica que circula a través del circuito. Recuerde que la intensidad de corriente eléctrica
matemáticamente se represente como la derivada de la carga respecto del tiempo; por lo
tanto, una vez que se obtiene la ecuación de carga eléctrica, basta con derivar este resultado
para encontrar la corriente. Las condiciones iniciales son la corriente y la carga en t = 0;
estas se obtienen cuando se hace conmutar un interruptor de una posición inicial a una
segunda posición.

91
Problema 16. Determine la carga y la intensidad de corriente eléctrica que circula a través
de un circuito eléctrico conectado en serie a una fuente de alimentación de voltaje de 10
volts. Los valores para los elementos pasivos son R = 10 Ω, L = 20 mH y C = 50 μF. La
corriente en t = 0 es de 1 Amper y la carga en el mismo instante de tiempo es de 2 μC.
Primero, se plantea la ley de voltajes de Para determinar la solución homogénea,

Kircchoff se iguala a cero la ecuación diferencial,
obteniendo así una ED homogénea
VL  VR  VC  VF
q   500q   1x10 6 q  0
Donde:
di sustituyendo
VL  L
dt
q 1
VR  Ri q  r
q   r 2
1
C
VC  idt
Por lo tanto
Por lo tanto se llega a r 2  500r  1x10 6  0

di 1
L  Ri   idt  V F r1  250  250 15i
dt C
dq r2  250  250 15i

Si, i  , entonces se tiene la ED de 2º
dt
orden De acuerdo a la solución que se obtiene
las raíces corresponden al caso III: Raíces
complejas y conjugadas
d 2q dq 1
L R  q  VF
dt dt C 
qh  e 250t C1 cos 250 15t  C2 sen 250 15t 
Sustituyendo valores
La solución particular que se propone es
d 2q dq 1
0.02  10  q  10
dt dt 50 x10 6 qp  A (2)
O también
d 2q dq
 500  1x10 6 q  10
dt dt q g  qh  q p (3)
 
q g  e 250t C1 cos 250 15t  C 2 sen 250 15t  A

92
Para determinar el valor de la constante Para obtener la otra constante se aplica la

A, en este ejemplo, es posible emplear el segunda condición inicial, la cual indica
método de coeficientes indeterminados, el que en t = 0, la intensidad de corriente es
cual es más sencillo de aplicar que el de 1 Amper. Es decir i(0)  1 , recordando
método de variación de parámetros. que la corriente es la derivada de la carga
por lo que q(0)  1 y se necesita la
q   500q   1x10 6 q  10 (4) derivada de la carga para poder aplicar la
condición establecida; derivando q g
Si
y p  A , entonces

qg  e250t  250 15C1sen 250 15t  250 15C2 cos 250 15t 

 250e250 t C1 cos 250 15t  C2 sen 250 15t 
y p  0
Sustituyendo t = 0 y q  1
y p  0

1  e0  250 15C1sen 0  250 15C2 cos 0 
Sustituyendo en (4)
 250e0 C1 cos 0  C2 sen 0t 
0  500(0)  1x10 A  10
6
Simplificando
Despejando
1  250 15C 2  250C1
10
A  0.0001
1x10 6 Por lo tanto
Sustituyendo en (3) C2  0.001031
 
q g  e 250t C1 cos 250 15t  C 2 sen 250 15t  0.00001 Finalmente la solución que representa la
carga en el circuito en función del tiempo
Para determinar los valores de las
constantes C1 y C2 (de la solución
general) se procede a sustituir las 
q  e 250 t  0.000008 cos 250 15t
condiciones iniciales. De acuerdo al
enunciado del problema se tiene que la
carga para un tiempo t = 0, es
 0.001031sen 250 15t  0.00001
q  2 x10 6 C . Sustituyendo en la última La función que representa la intensidad de
expresión corriente eléctrica en función del tiempo
es la derivada de la expresión anterior.
2 x10 6  e 0 C1 cos 0  C 2 sen 0  0.00001
2 x10 6  C1  0.00001
C1  0.000008

93
3.3.2 Sistemas Masa-Resorte
Suponga que se tiene un resorte conectado al techo de una estructura. El resorte posee una
longitud natural y se encuentra en reposo, es decir, no hay movimiento del mismo. Si en el
extremo inferior del resorte se agrega un peso w, el resorte sufrirá una elongación debido al
peso que se le agregó; este estiramiento (o deformación) es directamente proporcional al
peso que se agrega, es decir
F S
O también
F  kS
A la ecuación anterior se le conoce como ley de Hooke. F es el peso o fuerza que se agrega
al resorte y S es la elongación o deformación que sufre el resorte. Entre mayor peso se le
agregue al resorte, éste sufrirá una elongación mayor. k es la constante de restitución del
resorte; es una característica implícita para cada resorte. Depende del material del resorte,
grosor entre otros aspectos. En otras palabras, si se tienen diferentes tipos de resortes y a
cada uno de ellos se les agrega un objeto del mismo peso, no todos sufrirán la misma
elongación. Al agregarle el peso, las fuerzas actuantes en el sistema, son el peso y la fuerza
que se opone al peso; si el sistema continua en reposo se tiene
wF  0
mg  kS  0 (1)
Ahora suponga que el objeto se jala hacia abajo a una distancia y de la posición de
equilibrio y enseguida se suelta, ocasionando con esto que el sistema rompa su equilibrio y
el resorte empiece a contraerse y a alongarse producto de lo anterior. Por lo tanto al estar en
movimiento se tendría
w  k (S  y)  ma
Desarrollando
w  kS  ky  ma
De acuerdo a la ecuación (1) la última expresión se reduce a
ky  ma
ma  ky  0

94
Escribiendo la aceleración como la segunda derivada de la posición respecto del tiempo
d2y
m  ky  0
dt 2
Dividiendo toda la expresión entre la masa
d2y k
 y0
dt 2 m
k
Si  2  , se tiene
m
d2y
2
2 y  0
dt
La ecuación anterior describe la posición del cuerpo, que se conectó al resorte, en función
del tiempo; se tiene una ecuación diferencial homogénea de segundo orden de coeficientes
constantes. Resolviendo
r2  2  0
Implica que
r1  i, r2  i
Donde
k
 y se le conoce como frecuencia angular. Si se conoce la frecuencia angular es
m
posible determinar el período. Recuerde que el período de una función es el tiempo que
transcurre para que se complete un ciclo.

T , además
2
1
f 
T
La solución corresponde al caso III: Raíces complejas y conjugadas, y la solución a la que

se llega
y(t )  C1 cost  C2 sent

95
Note que al graficar la solución, esta corresponde a un movimiento perpetuo. A este tipo de
movimiento se le conoce como Armónico Simple. Las condiciones iniciales para calcular
los valores de las constantes se determinan de acuerdo a la distancia por debajo de la
posición de equilibrio y a la velocidad del cuerpo cuando inicia el análisis.
Realizando un análisis de manera similar al anterior, considerando que ahora además del
peso del cuerpo y la fuerza de restitución del resorte actúa una fuerza que realiza la acción
de amortiguar el movimiento se tendría la ecuación diferencial
d2y dy
m 2  c  ky  0
dt dy
Donde el segundo término es c veces la velocidad instantánea representando la fuerza

amortiguadora o retardadora que ocasionará que el movimiento se detenga, regresando al
reposo. Dividiendo entre la masa
d 2 y c dy k
  y0
dt 2 m dy m
O también
d2y dy
2
 2 2 y  0
dt dy
A este tipo de movimiento se le conoce como Libre Amortiguado y de acuerdo al valor de

las raíces de la ecuación característica se clasifica como amortiguado, críticamente
amortiguado y subamortiguado.
Un tercer tipo de movimiento que se abordará es el Movimiento Forzado. Para este caso se
tiene un sistema masa resorte conectado a una fuerza externa que actuará sobre el mismo en
todo momento. La ED que lo describe
d2y dy
2
 2   2 y  F (t )
dt dy
Donde
f (t )
F (t ) 
m
f (t ) Es la fuerza externa que actúa sobre el sistema

96
Problema 17. Se tiene un resorte al que se le agrega un peso de 30 N. Debido a éste peso el
resorte sufre una deformación de 20 cm. Una vez en equilibrio, el resorte se estira 10 cm.
por debajo de su posición de equilibrio y su suelta ocasionando que el cuerpo se ponga en
movimiento. Determine la ecuación que describe el movimiento del cuerpo en cualquier
instante así como su frecuencia y su período.
Primero se plantea la ED de acuerdo al La ED diferencial a resolver es

enunciado del problema
d 2 y 150
 y0
d2y k dt 2 3.06
 y0
dt 2 m
k d2y
Si  2  , se tiene  49.02 y  0
m dt 2
d2y Se trata de una ED homogénea, por lo que

2
2 y  0 se resuelve obteniendo la ecuación
dt
característica y de acuerdo a las raíces se
La constante de restitución del resorte (k) plantea su solución
se obtiene con la ley de Hooke
F  kS sustituyendo
F es la fuerza o peso que se aplica al y 1

resorte y S es la elongación o y  r
deformación que sufre el resorte debido a
y   r 2
ese peso, entonces
30  k (0.20) r 2  49.02  0
Despejando
Resolviendo la ecuación característica
k  150
r1  7.001i, r2  7.001i
Y la masa
De acuerdo a las raíces obtenidas la
w solución corresponde al caso III: Raíces
m
g complejas y conjugadas. Por lo tanto la
solución homogénea o complementaria se
30 expresa de la forma
m
9.81
y h  C1 cos 7.001t  C 2 sen 7.001t
m  3.06

97
De acuerdo al enunciado cuando se inicia

el análisis (t = 0) el cuerpo se encuentra a Por lo que
0.10 m de la posición de equilibrio. Esta
es la primera condición y(0)  0.10 ; 0  7.001C 2
sustituyéndola en la solución homogénea
C2  0
0.10  C1 cos 7.001(0)  C2 sen 7.001(0)
Con los valores de las constantes se
0.10  C1 cos 0  C2 sen 0 obtiene la solución en cualquier instante
bajo las condiciones indicadas en el
Por lo tanto problema
C1  0.10 y h  0.10 cos 7.001t
Ahora, para calcular el valor de la Con esta solución se puede determinar la

segunda constante se emplea la segunda posición del cuerpo en cualquier instante.
condición inicial del problema. La comprobación de las condiciones del
Continuando con el enunciado del problema consiste simplemente en
problema, se indica que el cuerpo se asignar el valor de cero a la variable
suelta después de haberlo estirado una independiente (tiempo) y entonces la
distancia por debajo de la posición de posición será de 0.10 m
equilibrio. Al momento de soltar el
cuerpo la velocidad es cero; en otras
palabras cuando se inicia el análisis
(t  0) la velocidad es también cero. Ya
que la velocidad es la derivada de la
posición se determina que
y (0)  0
Aplicando la condición anterior, en la

derivada de la solución homogénea
y h  7.001C1 sen 7.001t  7.001C 2 cos 7.001t
Sustituyendo y   0 y además t  0
0  7.001C1 sen 7.001(0)  7.001C2 cos 7.001(0)
Reduciendo
0  7.001C1 sen 0  7.001C2 cos 0
0  7.001C1 (0)  7.001C 2

98
Problema 18. Un resorte se deforma 3 cm debido a un peso de 8 N. Si el resorte se pone en

movimiento con una velocidad de 5 m/s dirigida hacia arriba desde su posición de
equilibrio y considerando que el medio en el que se desenvuelve actúa como una fuerza
amortiguadora igual a 2 veces la velocidad instantánea, determine la ecuación que describe
su movimiento en cualquier instante de tiempo.
Primero se plantea la ED de acuerdo al

enunciado del problema Y la masa
d2y dy k w
 2  y0 m
dt 2
dt m g
k 8
Si  2  , y 2  c , se tiene m
m 9.81
d2y dy m  0.81
2
 2 2 y  0
dt dt
La ED diferencial a resolver es
Se debe recordar que el segundo término
de la ecuación diferencial anterior d2y dy 266.67
2
2  y0
representa la fuerza o medio que actúa en dt dt 0.81
contra del movimiento (fuerza retardadora
o martiguadota) y esta es c veces la
velocidad instantánea. d2y dy
2
 2  329.22 y  0
dt dt
La constante de restitución del resorte (k)
se obtiene con la ley de Hooke Se trata de una ED homogénea, por lo que
se resuelve obteniendo la ecuación
F  kS
característica y de acuerdo a las raíces se
plantea su solución
F es la fuerza o peso que se aplica al
resorte y S es la elongación o Se obtiene la ecuación característica
deformación que sufre el resorte debido a
ese peso, entonces
r 2  2r  329.22  0
8  k (0.03)
Resolviendo la ecuación característica se
obtienen los valores
Despejando
r1  1  18.12i y r2  1  18.12i
8
k
0.03 De acuerdo a las raíces obtenidas la
solución corresponde al caso III: Raíces
k  266.67 complejas y conjugadas. Por lo tanto la

99
solución homogénea o complementaria se

expresa de la forma y h  e  t [C1 cos18.12t  C 2 sen18.12t ]
y h  e  t [C1 cos18.12t  C 2 sen18.12t ]
Si C1  0
De acuerdo al enunciado cuando se inicia
el análisis (t = 0) el cuerpo se encuentra y h  e t [(0) cos18.12t  C 2 sen18.12t ]
en la posición de equilibrio, por lo que
y(0)  0 ; sustituyéndola en la solución
Se simplifica a
homogénea
0 y h  e t [C 2 sen18.12t ]
0  e [C1 cos18.12(0)  C 2 sen18.12(0)]
Derivando
Por lo tanto, se tiene
y h  e  t 18.12C 2 cos18.12t  e  t C 2 sen18.12t
C1  0
Ahora, para calcular el valor de la Sustituyendo

segunda constante se emplea la segunda
condición inicial del problema.  5  e ( 0) 18.12C 2 cos18.12(0)  e ( 0) C 2 sen18.12(0)
Continuando con el enunciado del
problema, se indica que el cuerpo se pone  5  18.12C 2 cos(0)
en movimiento a partir de la posición de
equilibrio con una velocidad de 5 m / s Despejando
dirigida hacia arriba; en otras palabras
cuando se inicia el análisis (t  0) la C 2  0.28
velocidad (la derivada de la posición) está
dirigida hacia arriba (en sentido contrario Con los valores de las constantes se
a la gravedad) y se determina que obtiene la solución en cualquier instante
bajo las condiciones indicadas en el
y(0)  5 problema
Antes de proceder a la derivada de la y h  0.28e t sen18.12t

función homogénea, conviene sustituir el
valor que ya se calculó para la primera
constante, por lo que se llega a
3.4 Sección de Problemas
A continuación se proponen una serie de problemas con el propósito de complementar los

problemas que se resolvieron de manera detallada en la presente unidad y con esto el

100
alumno aplique los conocimientos adquiridos y le sirvan de guía en su estudio. La solución

de los mismos se encuentra en la última sección de la presente unidad.
3.4.1 Auto Evaluación

y(0)  4 , y (0)  1
En los problemas 1-15 determine la

solución de la ED homogénea que se
indica. Donde sea necesario determine En los problemas 16-20 determine la
el valor de las constantes. solución de la ED no homogénea
empleando para ello el método de
1. y  4 y  3 y  0 coeficientes indeterminados. Donde se
indique, aplique las condiciones
2. y  y   3 y  0 iniciales.
3. y  25 y  0
16. y  3 y  2 y  6
4. y  16 y  0
1
17. y   y   y  x 2  2 x
5. y   4 y   0 4
6. y  6 y  9 y  0 18. y   3 y  48 x 2 e 3 x
1 19. y  y  2 xsenx

7. y   y   y0
4
20. y   2 y   5 y  e x cos 2 x
8. y  2 y  4 y  0
9. y  2 y  y  0 En los problemas 21-25 determine la

solución de la ED no homogénea
10. y  6 y  5 y  0 ; sujeta a empleando para ello el método de
variación de parámetros. Donde se
indique, aplique las condiciones
y(0)  3 , y (0)  1
iniciales.
11. y  2 y  3 y  0
21. y   3 y   2 y  sene x
12. y  10 y  25 y  0
22. 3 y   6 y   30 y  e x tan 3x
13. 2 y  6 y  2 y  0
x
17 23. 4 y   y  xe 2 ; con y(0)  1 , y (0)  0
14. y   3 y   y0
2
24. y  y  cosh x
15. y  12 y  36 y  0 ; sujeta a

101
25. y  2 y  y  e  x ln x
26. Un cuerpo de masa 1.5 kg. estira 3 cm. un resorte. Si se tiene una fuerza retardadora
equivalente al doble de la velocidad instantánea y el cuerpo se suelta 4 cm. por debajo de la
posición de equilibrio con una velocidad dirigida hacia arriba de 6m / s ¿cuál es la ecuación
que define el movimiento en cualquier instante de tiempo?
27. Un cuerpo con masa 0.5kg estira un resorte 50cm , un poco después se suelta en el
instante t  0 de un punto que está a 8cm por debajo de la posición de equilibrio con una
velocidad dirigida de 2m / s . Determine la función que describe el movimiento del cuerpo
en cualquier instante.
28. Una banda elástica está hecha de un material tal que con un peso de 4 N colgado de
ella, se alarga 6cm . Si una fuerza de 0.5 cos16t N está actuando sobre la banda de tal forma
que el cuerpo adherido a ella es sacado de su equilibrio por un movimiento hacia arriba que
tiene una velocidad de 3m / s , ¿cuál es la posición del objeto en para cualquier instante de
tiempo?
29. Una masa de 100 gr se adhiere a un resorte de acero con una longitud original de 50cm .
El resorte se extiende 5cm por efecto de esta masa. Si se hace mover la masa hacia abajo
con una velocidad de 10cm / s determina la expresión que representa el movimiento
provocado en cualquier instante de tiempo.
30. Un circuito RLC tiene una resistencia de 10 ohms, una capacitancia de 10 -2 farads y una
inductancia de 0.5 henrios. El voltaje que se suministra al circuito eléctrico es de 12 volts.
Suponiendo que no hay corriente inicial o carga en el capacitor, encuentre la carga en
función del tiempo.
31. Un costal de 10kg de masa es atado a un resorte que por su efecto lo estira 70cm de su
longitud normal. El costal se pone en movimiento al estirar el resorte y soltarlo, con una
velocidad inicial de 1m / s en dirección hacia arriba. Encuentre la expresión que representa
el movimiento resultante considerando que existe una fuerza retardadora igual a 90 veces la
velocidad instantánea.

102
3.4.2 Respuestas Auto Evaluación 16. y p  c1e  x  c2 e 2 x  3
1. y h  c1e 3 x  c 2 e x 17. y p  c1e 2 x  c 2 xe 2 x  x 2  4 x  7

2
1
 x  11 11  18. y g  y h  y p
2. y h  e 2
C1 cos x  C2 sen x
 2 2 
yh  C1 cos 3x  C2 sen 3x
3. y h  C1 cos 5 x  C 2 sen 5 x
 4
y p    4 x 2  4 x  e 3 x
4. y h  c1e 4x
 c2 e 4 x
 3
5. y h  c1  c 2 e 4 x
19. y g  y h  y p
6. y h  c1e 3 x  c 2 xe 3 x
y h  C1 cos x  C 2 senx
1 1
 x  x
7. y h  c1e 2
 c2 xe 2
1 1
y p   x 2 cos x  senx
2 2
8. yh  c1e 1 5 x  c2 e 1 5 x
20. y g  y h  y p
9. yh  c1e 1 2 x  c2 e 1 2 x
y h  C1 cos 2 x  C 2 sen 2 x
10. y h  4e  x  e 5 x

11. yh  e  x C1 cos 2 x  C2 sen 2 x  yp 
1 x
4
xe sen 2 x
12. y h  2e 5 x  15 xe 5 x
21. y g  c1e 2 x  c2 e  x  e 2 x sene x
13. y h  c1e 2 x  c 2 e x
3
 x  5 5  22. y g  y h  y p
14. y h  e 2
C1 cos 2 x  C 2 sen 2 x 
y h  e x (C1 cos 3x  C 2 sen 3x)
15. y h  4e 6 x  25 xe 6 x
1 x
yp   e cos 3 x ln sec 3 x  tan 3 x
27

103
1 1 
x x x
1  3
23. y g  e 2  e 2   x 2  x e 2
4 4 8 4 
1
24. y g  c1e x  c 2 e  x  xsenhx
2
25. y g  y h  y p
y h  c1e  x  c 2 xe  x
1 2 x 3
yp  x e ln x  x 2 e  x
2 4
26.
y  e 0.66t [0.04 cos18.06t  0.3307 sen18.06t ]
27. y  0.08 cos 4.42t  0.45sen 4.42t
28.
y  0.013 cos12.9t  0.23sen12.9t  0.013 cos16t
29. y  0.71sen14t
30.
3 3  3
q(t )  e 10t  cos10t  sen10t  
 25 250  25
1 1
31. y   e  2t  e 7 t
5 5

104

Apuntes de Ecuaciones Diferenciales UNIDAD 3

Cargado por

Información del documentohacer clic para expandir la información del documento

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Apuntes de Ecuaciones Diferenciales UNIDAD 3

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Apuntes de Ecuaciones Diferenciales UNIDAD 3

Cargado por

Copyright:

Formatos disponibles

Unidad 3: ED Lineales de n-ésimo Orden con Coeficientes Constantes

ECUACIONES DIFERENCIALES LINEALES DE n-ESIMO

3.1 Solución de la Ecuaciones Diferenciales Homogéneas con

3.2 Solución de la Ecuaciones Diferenciales No Homogéneas

3.2.1 Método de los Coeficientes Indeterminados

3.3.1 Circuito Eléctrico RLC Serie

3.4 Sección de Problemas

3.4.1 Auto Evaluación

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

Objetivos particulares de la unidad.- El alumno resolverá Ecuaciones Diferenciales de n-

En el presente capítulo se estudiarán los métodos de solución para las Ecuaciones

3.1 Solución de la Ecuaciones Diferenciales Homogéneas con

dny d n1 y d2y dy

Al multiplicarla por 1/a2 se obtiene la ecuación diferencial equivalente

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

Ejemplo. Clasifique las siguientes ecuaciones diferenciales de acuerdo a

a) sus coeficientes (constantes o variables)

Para el problema 3 el coeficiente que afecta a la primera derivada no es constante, por lo

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

Desarrollando las derivadas

Otra forma de encontrar la ecuación característica es mediante la siguiente sustitución:

Es una ecuación de segundo grado, se le conoce como ecuación característica, ecuación

Caso I: Las raíces son reales y diferentes (  2  4  0 )

Caso II: Las raíces son reales e iguales (  2  4  0 )

Caso III: Las raíces son complejas y conjugadas (  2  4  0 )

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

O de manera más general, para cualquier constante c

De acuerdo al Principio de Superposición:

O de manera más general, para cualquier constante c

De acuerdo al Principio de Superposición:

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

O de manera más general, para cualquier constante c

De acuerdo al Principio de Superposición:

De esta última solución al sustituir r1 y r2 se tiene:

y h  c1e ( a bi ) x  c 2 e ( a bi ) x

Eliminando los paréntesis

y h  c1e ax bxi  c 2 e ax bxi

Aplicando propiedades de la función exponencial

y h  c1e ax e bxi  c 2 e ax e bxi

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

Por la fórmula de Euler

e bxi  cos bx  isenbx

Sustituyendo en la solución homogénea

y h  c1e ax (cosbx  isenbx )  c 2 e ax (cosbx  isenbx )

y h  c1e ax cos bx  c1e ax isenbx  c 2 e ax cos bx  c 2 e ax isenbx

Agrupando parte real y parte imaginaria

y h  (c1  c 2 )e ax cos bx  i (c1  c 2 )e ax senbx

Si C1  c1  c2 y Si C2  i(c1  c2 ) , se llega a la solución equivalente

y h  C1e ax cos bx  C 2 e ax senbx

Finalmente, factorizando la función exponencial

y h  e ax C1 cos bx  C 2 senbx 

Es la solución de la ED para el caso III, no olvidar que

1. y   4 y   3 y  10 ; sujeta a las condiciones y(0)  2, y(0)  1

2. y  9 y  0 ; sujeta a la condición inicial y( )  2 , y (0)  1

4. y  4 y  4 y  0 ; sujeta a la condición inicial y(0)  1, y(0)  4

Elaboró: Profesor Víctor Manuel Salazar ESIME-ICE-Matemáticas

1. y  4 y  3 y  0 ; sujeta a las

Se obtiene la ecuación característica c1  c2  2 (1)