ECUACIONES DIFERENCIALES LINEALES

DE SEGUNDO ORDEN
YOEL GUTIRREZ
UNEXPO-PUERTO ORDAZ

Introduccin

En la ingeniera se desarrollan modelos que producen ED lineales de segundo orden.

Algunos ejemplos de estos modelos son las vibraciones mecnicas.
Consideremos un resorte que resiste compresin y estiramiento, sujeto a un soporte
rgido. Un cuerpo de masa m se sujeta en el otro extremo del resorte y lo estira una
lingitud s, llegando a una posicin de equilibrio (ver Figura 1)
Soporte rgido

Resorte
sin
estirar

m
Posicin de
equilibrio

Figura 1: Sistema resorte-masa

Denotemos con x la distancia del cuerpo a su posicin de equilibrio. Tomamos
x > 0 cuando el resorte est por debajo de la posicin de equilibrio y x < 0 cuando
est por encima.
La segunda ley de Newton, fuerza es igual a masa por aceleracin
m

d2 x
= F;
dt2

es una EDO de segundo que parece con mucha frecuencia en la prctica. Aplicaremos
esta ley al sistema masa-resorte para obtener una EDO de segundo orden ms general,
1

que gobierne el movimiento del cuerpo, tomando en cuenta las fuerzas que actan sobre
l debido a la elasticidad del resorte, el peso del cuerpo, la friccin (amortiguamiento)
y las posibles fuerzas externas.
Al desplazar la masa m con respecto de la posicin de equilibrio, el resorte se
estira o se comprime y ejerce una fuerza que reiste al desplazamiento. Para la mayor
parte de los resortes, esta fuerza es directamente proporcional a la distancia que el
resorte ha sido alargado o comprimido, por lo que est dada por
Fr =

(1)

k (s + x) ;

donde la constante positiva k es la rigidez del resorte y el signo negativo indica la

naturaleza de oposicin de la fuerza. La ecuacin (1), conocida como la ley de
Hooke, slo es vlida para desplazamientos sucientemente pequeos.
El peso del cuerpo
W = mg
est equilibrado por Fr : En la posicin de equilibrio
mg

(2)

ks = 0

Prcticamente todos los sistemas mecnicos experimentan la fuerza de friccin;

por lo general, para el movimiento de vibracin, esta fuerza se modela mediante un
trmino proporcional a la velocidad:
Ff =

dx
;
dt

(3)

donde b 0 es el coeciente de amortiguamiento, y el signo negativo tiene el mismo

signicado que en la ecuacin (1).
Las otras fuerzas que actan sobre el oscilador se consideran por lo general como
externas al sistema. Aunque estas pueden ser gravitacionales, elctricas o magnticas,
lo comn es que las fuerzas externas ms importantes sean transmitidas a la masa
sacudiendo la base de la que cuelga el resorte. Representaremos todas las fuerzas
externas por una funcin conocida Fe (t):
Si no hay otras fuerzas de retardo que acten sobre el sistema, entonces la fuerza
neta o resultante es
Fn = Fr + W + Ff + Fe (t)
Aplicando la segunda Ley de Newton obtenemos
m

d2 x
=
dt2
=

Como mg

k(s + x) + mg
ks

kx + mg

dx
+ Fe (t)
dt
dx
b + Fe (t)
dt
b

ks = 0;
m

d2 x
=
dt2

kx

b
2

dx
+ Fe (t)
dt

As, la EDO para el sistema mecnico masa-resorte es

m

d2 x
dx
+ b + kx = Fe (t):
2
dt
dt

(4)

Nuestro objetivo es obtener una descripcin cualitativa de las soluciones de la EDO

de la forma (4), llamada ED lineal de segundo orden con coecientes constantes.
Para comenzar el estudio de estas ecuaciones diferenciales necesitamos ciertas
deniciones y terminologas comnes.
Denicin 1.1 Una ecuacin diferencial ordinaria es lineal de segundo orden
si se puede escribir de la forma
a2 (x)
o ms brevemente

d2 y
dy
+ a1 (x) + a0 (x)y = f (x);
2
dx
dx
00

(5)

a2 (x)y + a1 (x)y + a0 (x)y = f (x);

donde a2 ; a1 ; a0 y f dependen slo de x y no de y: Cuando a2 ; a1 ; a0 son constantes,
diremos que la ecuacin (5) tiene coecientes constantes, en caso contrario, tiene
coecientes variables.
Modelaremos problemas mediantes un PVI de segundo orden que involucra un
ED lineal, estos PVI son de la forma
00

a2 (x)y + a1 (x)y + a0 (x)y = f (x);

y(x0 ) = y0 ;

y 0 (x0 ) = y1 :

(6)

Recurdese que, para un problema como ste, se busca una funcin denida en
algn intervalo I que contenga a x0 , y satisfaga la ecuacin diferencial y las 2 condiciones iniciales especicadas en x0 : y(x0 ) = y0 ; y 0 (x0 ) = y1 .
Existe un teorema de existencia y unicidad para el PVI anteior, similar al correspondiente teorema para el PVI que involucra una ED lineal de primer orden.
Teorema 1.1 (Teorema de existencia y unicidad) Sean a2 (x); a1 (x); a0 (x) y f (x)
funciones continuas en un intervalo I, y sea a2 (x) 6= 0 para cada x del intervalo. Si
x0 es cualquier punto en I y si y0 ; y1 ; son nmeros arbitrarios, el PVI (6) tiene una
y slo una solucin y(x) en el intervalo I.
La cuestin de existencia de una solucin no constituye en realidad un problema.
Sin embargo, es muy importante saber que la solucin es nica.

Ecuaciones lineales homogneas: Caso general

Denicin 2.1 Una ecuacin lineal de segundo orden de la forma

d2 y
dy
+ a1 (x) + a0 (x)y = 0
2
dx
dx
se llama homognea, mientras que una ecuacin
a2 (x)

dy
d2 y
+ a1 (x) + a0 (x)y = f (x)
2
dx
dx
donde f (x) nos es idnticamente cero, se llama no homognea.
a2 (x)

(7)

(8)

Al estudiar la ecuacin no homognea (8), es necesario considerar a la par la

ecuacin homognea (7). Bajo estas condiciones se habla de (8) como una ecuacin
completa y de (7) como la ecuacin reducida asociada a ella.
En el siguiente teorema veremos que la suma o superposicin de dos soluciones de
una ecuacin diferencial lineal homognea de segundo orden tambin es una solucin.
Teorema 2.1 (Principio se superposicin, ecuaciones homogneas)
Si
y1 ; y2 son 2 soluciones de la ED lineal homognea (7) de orden 2 en un intervalo
I, la combinacin lineal
y(x) = c1 y1 (x) + c2 y2 (x)
en donde c1 ; c2 son constantes arbitrarias, tambin es una solucin cuando x est en
I.
Prueba. Como y1 ; y2 son 2 soluciones de la ED lineal homognea (7) en I; entonces
a2 (x)

dy1
d2 y1
+ a1 (x)
+ a0 (x)y1 = 0
2
dx
dx

y
dy2
d2 y2
+
a
(x)
+ a0 (x)y2 = 0
1
dx2
dx
en I: Por lo tanto, si c1 ; c2 son constantes arbitrarias y
a2 (x)

y = c1 y1 + c2 y2 ;
entonces
d2 y
dy
a2 (x) 2 + a1 (x) + a0 (x)y
dx
dx
d2
d
= a2 (x) 2 (c1 y1 + c2 y2 ) + a1 (x) (c1 y1 + c2 y2 ) + a0 (x) (c1 y1 + c2 y2 )
dx
dx
d2 y1
dy1
d2 y2
dy2
= c1 a2 (x) 2 + a1 (x)
+ a0 (x)y1 + c2 a2 (x) 2 + a1 (x)
+ a0 (x)y2
dx
dx
dx
dx
= c1 :0 + c2 :0
= 0
en I; y el teorema queda demostrado.
4

Ntese que:
1. Un mltiplo constante y(x) = c1 y1 (x) de una solucin y1 (x) de una ED lineal
homognea tambin es una solucin.
2. Una ED lineal homognea siempre tiene la solucin trivial y(x) = 0.
Sabemos que la solucin general de una ED lineal homognea de segundo orden
es biparamtrica (tiene dos parmetros). Ahora, por el principio de superposicin, si
y1 ; y2 son 2 soluciones de la ED lineal homognea (7) en un intervalo I, entonces
y(x) = c1 y1 (x) + c2 y2 (x)

(9)

en donde c1 ; c2 son constantes arbitrarias, tambin es una solucin en I. Entonces

nos preguntamos, es la ecuacin (9) la solucin general de la ED lineal homognea
(7) en I?, esto es, toda solucin de la ecuacin homognea (7) se puede escribir de
la forma (9), para una eleccin adecuada de las constantes c1 ; c2 ?
Para dar respuesta a esta pregunta ntese que:
1. Si y2 (x) es la solucin idnticamente nula en I, entoces
c1 y1 (x) + c2 y2 (x) = c1 y1 (x)
tiene en realidad una sola constante y no puede esperarse que sea la solucin
general de la EDO (7) en I:
2. Si y2 (x) es un mltiplo constante de y1 (x) en I; digamos y2 (x) = ky1 (x) en I;
entonces de nuevo
c1 y1 (x) + c2 y2 (x) = c1 y1 (x) + ky1 (x) = (c1 + k)y1 (x) = Cy1 (x)
tiene en realidad una sola constante.
Luego, una de las condiciones que necesitamos es la independencia lineal.
Denicin 2.2 Dos funciones f1 (x) y f2 (x) se dice que son linealmente dependientes (l.d.) en un intervalo I si existen dos constantes c1 ; c2 no todas cero, tales que
para toda x en I
c1 f1 (x) + c2 f2 (x) = 0:
(10)
En caso contrario se dice que el conjunto es linealmente independiente (l.i.).
En otras palabras, dos funciones son l.i. en un intervalo si las nicas constantes
para las que se cumple (10) para toda x en el intervalo son c1 = c2 = 0:

Ntese que si las dos funciones f1 (x) y f2 (x) son l.d. en un intervalo, existen
constantes, c1 y c1 , que no son cero a la vez, tales que, para toda x en el intervalo se
cumple (10); por consiguiente, si suponemos que c1 6= 0, entonces
c2
f2 (x);
c1

f1 (x) =

esto es, si dos funciones son l.d., entonces una es mltiplo constante de la
otra. Al revs, si f1 (x) = c2 f2 (x) para alguna constante c2 , entonces
( 1)f1 (x) + c2 f2 (x) = 0
para toda x en algn intervalo. As, las funciones son l.d. porque al menos unas de
las constantes no es cero. Llegamos a la conclusin de que dos funciones son l.i.
en un intervalo I cuando ninguna es mltiplo constante de la otra en I.
(x)
Esto es, el cociente ff21 (x)
no es constante en un intervalo en que f1 (x) y f2 (x) son l.i..
Para estudiar la dependencia lineal de 2 soluciones particulares de una ED lineal
homognea podemos recurrir a un determinante, como se muestra en el siguiente
teorema
Teorema 2.2 Sean y1 ; y2 ; 2 soluciones de la ecuacin diferencial homognea de orden
2 (7) , donde a2 ; a1 y a0 son continuas en algn intervalo dado I y an (x) 6= 0 para
toda x en I.
1. y1 ; y2 ; son linealmente dependientes en I si y slo si el Wronskiano de y1 e y2
dado por
y1 y2
W (y1 ; y2 ) =
y10 y20
es idnticamente nulo en el intervalo.
2. y1 ; y2 ; son linealmente independientes en I si y slo si W (y1 ; y2 ) 6= 0 para cada
x en I.
3. W (y1 ; y2 ) = 0 para cada x 2 I, o bien, W (y1 ; y2 ) 6= 0 para cada x 2 I.
Cualquier conjunto y1 ; y2 de dos soluciones l.i. de la ED lineal homognea (7)
de segundo orden en un intervalo I se dice que es un conjunto fundamental de
soluciones en el intervalo.
Semejante al hecho de que cualquier vector en dos dimensiones puede expresarse como una combinacin lineal de los vectores linealmente independientes i y
j; Cualquier solucin de una ED lineal homegnea de segundo orden en un intervalo
I puede expresarse como una combinacin lineal de dos soluciones l.i en I:

Teorema 2.3 (Solucin general, ecuaciones homogneas) Si y1 ; y2 ; son 2 soluciones linealmente independientes de la ecuacin homognea de orden 2 (7), en un
intervalo I, entonces
y = c1 y1 + c2 y2
(11)
es una solucin de (7) para cualesquiera constantes c1 ; c2 . Recprocamente, toda solucin de (7) tiene la forma (11) para selecciones apropiadas de las constantes c1 ; c2 .
Prueba. Sean y1 ; y2 ; 2 soluciones linealmente independientes de la EDO lineal homognea
dy
d2 y
a2 (x) 2 + a1 (x) + a0 (x)y = 0
dx
dx
en un intervalo I; por el principio de superposicin es inmediato que
y = c1 y1 + c2 y2
tambin es una solucin en I para cualesquiera constantes c1 ; c2 .
Sea (x) una solucin de la EDO lineal homognea en I; sea x = t un punto en I
y consideremos el sistema de ecuaciones
c1 y1 (t) + c2 y2 (t) = (t)
0
0
0
c1 y1 (t) + c2 y2 (t) = (t):

(12)

Como y1 ; y2 ; son soluciones linealmente independientes de la ED lineal homognea

de segundo orden en I, entonces W (y1 ; y2 ) 6= 0 para cada x 2 I, por lo tanto, como
t2I
y1 (t) y2 (t)
6= 0:
0
0
y1 (t) y2 (t)
Esto garantiza que el sistema (12) tiene solucin nica, esto es, podemos determinar
c1 y c2 de manera nica. Luego, toda solucin de (7) se puede escribir de la forma
(11) para selecciones apropiadas de las constantes c1 ; c2 .

Ecuaciones
lineales
coecientes constantes

homogneas

con

Estudiaremos con detalle las soluciones de la ED lineal homognea

ay 00 + by 0 + cy = 0

(13)

en el caso especial en que a, b y c son constantes reales.

Nuestro punto de partida es la propiedad de la funcin exponencial emx de que sus
derivadas son todas mltiplos de la propia funcin, lo que nos conduce a considerar
y = emx
7

como una posible solucin de (13) si la constante m se escoge adecuadamente. Como

y 0 = memx e y 00 = m2 emx , sustituyendo en (13) vemos que
(am2 + bm + c)emx = 0;

(14)

y puesto que emx nunca se anula, (14) se cumple si y slo m satisface la ecuacin
am2 + bm + c = 0;

(15)

llamada ecuacin auxiliar o ecuacin caracterstica de la EDO (13). Examinaremos tres casos: Las soluciones de la ecuacin caracterstica que corresponden a races
reales distintas, races reales e iguales y races complejas conjugadas.
3.0.1

Caso 1: Races reales distintas.

Si la ecuacin (15) tiene dos races reales distintas, m1 y m2 , llegamos a las soluciones
y1 = em1 x y y2 = em2 x . Estas funciones son linealmente independientes en ( 1; 1),
en efecto, como m2 6= m1
W (y1 ; y2 ) =

em1 x
em2 x
m1 x
m1 e
m2 em2 x

= (m2

m1 ) e(m1 +m2 )x

es no nulo para cada x 2 ( 1; 1); en consecuencia, la solucin general de la ecuacin

(13) en ese intervalo es
y = c1 em1 x + c2 em2 x :
3.0.2

Caso 2: Races reales e iguales.

Si la ecuacin (15) tiene un nica raz real, entonces b2 4ac = 0 y la raz es m =

en este caso
b
y1 = emx = e 2a x

b
:
2a

es una solucin de la ED lineal homognea (13) en ( 1; 1). Se demuestra que

y2 = xemx = xe

b
x
2a

es una segunda solucin de la ecuacin (13) en ( 1; 1), en efecto, sustituyendo y2 ;

0
00
y2 y y2 en (13) obtenemos que

ay200 + by20 + cy2 = a

=
=
=
=

b
b
b2
b
+ 2 x e 2a x + b 1
x e
a 4a
2a
b
b2
b2
x + cx e 2a x
b+ x+b
4a
2a
2
2
b
b
2b + 4ac
xe 2a x
4a
2
b
b
4ac
xe 2a x
4a
b
0
xe 2a x
4a

b
x
2a

+ c xe

b
x
2a

= 0:
Por otra parte, como
e

w(y1 ; y2 )(x) =
=
= e

b
x
2a

b
e
2a

b
x
2a

b
x e
2a

xe

b
x
2a

b
x
2a

1
b
x
a

bx
e
2a

b
x
2a

b
x
a

b
x
a

no se anula en ( 1; 1), las dos soluciones, y1 y y2 son l.i en en ( 1; 1): Por lo

tanto (15) tiene a
y = c1 emx + c2 xemx = (c1 + c2 x) emx :
como su solucin general en ( 1; 1).
3.0.3

Caso 3: Races complejas conjugadas.

Si la ecuacin (15) tiene dos races complejas conjungadas, digamos m1 =

m2 =
i , donde y son nmeros reales; entonces,
y = c1 e(

+i )x

+ c2 e(

i )x

+i y
(16)

es solucin de la ecuacin (13) para cada eleccin de las constantes c1 y c2 : Sin

embargo, en la prctica se preere trabajar con funciones reales y no con exponenciales
complejas. Con este objeto se usa la frmula de Euler:
ei = cos + isen ;
en que

es un nmero real.

Como (16) es una solucin de (13) para cualquier eleccin de las constantes c1 y
c2 , si c1 = c2 = 1, obtenemos la solucin
y1 = e(
y si c1 = 1 y c2 =

+i )x

i )x

+ e(

= e x (ei

+e

i x

i x

) = 2e

cos( x);

1, obtenemos la solucin

y2 = e(

+i )x

i )x

e(

= e x (ei

) = 2ie x sen( x):

En consecuencia, por el principio de superposicin para las ecuaciones homogneas

1
y1 = e
2

1
y2 = e x sen( x)
2i

cos( x) y

tambin son soluciones de la ecuacin (13). Adems, se demuestra que esas soluciones
son linealmente independientes en ( 1; 1); por tanto, la solucin general de la
ecuacin (13) en ( 1; 1) es
y = (c1 cos( x) + c2 sen( x))e x :
Teorema 3.1 (Existencia y unicidad: Caso homogneo) Para
cualesquiera
nmeros reales a2 ;,a1 ;,a0 ; x0 ; y0 ; y1 ; con a2 6= 0; existe una nica solucin del PVI
00

a2 y + a1 y + a0 y = 0; y(x0 ) = y0 ;

y 0 (x0 ) = y1 ;

vlida para toda x en ( 1; 1):

Ejercicios 1
1. Verique que la ecuacin
m

d2 x
+ kx = 0
dt2

tiene una solucin de la forma x(t) = sen!t; donde ! =

2. Para la ecuacin
m

k
:
m

d2 x
dx
+
b
+ kx = 0;
dt2
dt

verique lo siguiente:
(a) Si x(t) es una solucin, tambin lo es cx(t), para cualquier constante c,
(b) si x1 (t) y x2 (t) son soluciones, tambin lo es su suma.
3. Verique que x(t) = 3sen3t + cos 3t es una solucin del PVI
00

2x + 18x = 0;

x(0) = 1;

x0 (0) = 6:

Determine el mximo de jx(t)j para -1 < t < 1:

10

4. Determinar si las funciones f (x) y g(x) son linealmente independientes en el

intervalo (0; 1):
(a) f (t) = e t sen2t y g(t) = e

cos 2t

(b) f (t) = tan2 t

sec2 t y g(t) = 3

(c) f (t) = 0

g(t) = et

(d) f (t) = x y g(t) = x2

5. Sean f (t) = t3 y g(t) = jt3 j :
(a) Son f (t) y g(t) linealmente independientes en los siguientes intervalos?
i) (0; 1); ii) ( 1; 0); iii) ( 1; 1)

(b) Calcules el wronskiano de esta dos funciones en ( 1; 1)

6. Hallar una solucin general de cada ecuacin
00

(a) x + x

x=0

00

(b) 2w + 7w
(c) x

00

4w = 0
11x = 0

00

(d) 4y + 20y + 125y = 0

00

(e) y + y = 0
(f) y

00

10y + 26y = 0

(g) y

00

4y + 7y = 0

00

(h) y + 4y + 8y = 0
7. Resuelva cada PVI
0

(a) y + y = 0; y(0) = 2; y (0) = 1

(b) y

00

4y + 3y = 0; y(0) = 1; y (0) =

(c) y

00

6y + 9y = 0; y(0) = 2; y (0) =

00

1
3
25
3

(d) y + 2y + 17y = 0; y(0) = 1; y (0) =

00

(e) y + 9y = 0; y(0) = 1; y (0) = 1

11

Ecuaciones lineales no homogneas

El siguiente teorema ser de mucha utilidad cuando consideremos los mtodos para
encontrar soluciones particulares de EDO lineales no homogneas con coecientes
constantes.
Teorema 4.1 (Principio de superposicin, ecuaciones no homogneas) Si ypi
representa una solucin particular de ED lineal no homognea
a2 (x)

dy
d2 y
+ a1 (x) + a0 (x)y = fi (x);
2
dx
dx

en donde i = 1; 2; : : : ; k. Entonces
yp = c1 yp1 + c2 yp2 + : : : + ck ypk
es una solucin particular de
a2 (x)

d2 y
dy
+ a1 (x) + a0 (x)y = c1 f1 (x) + c2 f2 (x) + : : : + cfk (x);
2
dx
dx

para cualquier eleccin de las constantes c1 ; c2 ; :::; ck :

Qu forma tiene las soluciones generales de las ED lineales no homogneas? La
respuesta a esta pregunta la mostramos a continuacin.
Teorema 4.2 (Solucin general, ecuaciones no homogneas) Si y1 ; y2 son 2 soluciones linealmente independientes de la ecuacin
a2 (x)

d2 y
dy
+
a
(x)
+ a0 (x)y = 0
1
dx2
dx

(17)

y yp es una solucin particular de

a2 (x)

dy
d2 y
+ a1 (x) + a0 (x)y = f (x)
2
dx
dx

(18)

donde f es continua en I, entonces

y = c1 y1 + c2 y2 + yp

(19)

es una solucin de (18) para cualesquieras constantes c1 ; c2 . Recprocamente, toda

solucin de (18) tiene la forma (19) para selecciones apropiadas de las constantes
c 1 ; c2 .

12

Usulamente, la solucin general de la homognea asociada c1 y1 + c2 y2 se denota

por yc y se llama solucin complementaria. As, la solucin general de la ecuacin
no homognea (18), en un intervalo I es de la forma
y = yc + yp
Para encontrar la solucin general de la ecuacin diferencial lineal no homognea
(18), debemos pasar por dos momentos
1. Determinar la solucin general de la ecuacin homognea asociada, yc .
2. Establecer cualquier solucin particular, yp , de la ecuacin no homognea.
Entonces, la solucin general de (18) en un intervalo es y = yc + yp .
Estudiaremos las ecuaciones lineales en donde los coecientes de la homognea
asociada son constantes, en este caso no existir problema alguno para hallar yc :
Pero, en tales circunstancias, cmo hallar yp ? Existen muchos mtodos por medio
de los cuales se pueden obtener soluciones particulares, a continuacin nos referiremos
a dos de ellos.

4.1

El mtodo de los coecientes indeterminados

Un mtodo a menudo usado en Fsica e Ingeniera es el mtodo de los coecientes

indeterminados. La idea bsica es una conjetura o propuesta coherente acerca de la
forma de yp originada por los tipos de funciones que forman el dato g(x). El mtodo
es bsicamente directo, pero est limitado a ecuaciones lineales no homogneas, como
la ecuacin (18), en que:
1. Los coecientes a0 ; a1 ; a2 son constantes.
2. f (x) es una funcin polinomial, una funcin exponencial erx donde r es constante, una funcin seno o coseno como senrx, cos rx, o sumas y productos
nitos de estas funciones.
A continuacin mostramos la forma de yp originada por algunos tipos de funciones
que forman el dato f (x):
1. Para un polinomio de grado n, asuma un polinomio de grado n.
2. Para trminos como senrx o cos rx, asuma Asenrx + B cos rx.
3. Para trminos como erx , asuma Aerx .
4. Para un producto nito de trminos como los anteriores, asuma una combinacin lineal de todas las funciones linealmente independientes generadas por
diferenciaciones repetidas de f (x)
13

5. Para una suma o diferencia nita de trmininos como los anteriores, aplique el
principio de superposicin.
Si algunos de los trminos asumidos ocurre en la solucin complementaria, debemos multiplicar estos trminos asumidos por una potencia de la forma xn ; con n el
entero positivo ms pequeo, de modo que ninguno de los trminos asumidos aparezca
en la solucin complementaria.

4.2

El mtodo de variacin de parmetros

Consideremos la ED lineal no homonea

a2 (x)y 00 + a1 (x)y 0 + a0 (x)y = f (x):

(20)

El mtodo de los coecientes indeterminados, para hallar una solucin particular de la

ecuacin (20) tiene dos serias limitaciones: puede usarse slo cuando los coecientes
a2 , a1 y a0 son constantes, e incluso entonces slo funciona si el trmino de la derecha
f (x) tiene una forma particularmente sencilla. Dentro de estas limitaciones, ese
mtodo suele ser el ms simple de aplicar.
Ahora desarrollaremos otro mtodo ms potente que funciona siempre, sean cuales
sean a2 , a1 y a0 y f (x), y supuesto slo que la solucin general de la correspondiente
ecuacin homognea
a2 (x)y 00 + a1 (x)y 0 + a0 (x)y = 0;
(21)
se conoce de antemano.
La ED (20) se puede escribir de la forma
y 00 + P (x)y 0 + Q(x)y = g(x);

(22)

(x)
; Q(x) = aa02 (x)
y g(x) = af2(x)
Supongamos que P (x); Q(x) g(x)
donde P (x) = aa21 (x)
(x)
(x)
son continuas en algn intervalo I; y que se ha encontrado de algn modo la solucin
general
y = c1 y1 (x) + c2 y2 (x)

de la homognea asociada.
El mtodo de variacin de parmetros permite encontrar dos funciones desconocidas, v1 (x) y v2 (x); de manera tal que
yp = v1 y1 + v2 y2

(23)

sea una solucin particular de (22). Para encontrar estas dos funciones desconocidas,
empecemos calculando la derivada de (23), agrupada como sigue:
yp0 = (v1 y10 + v2 y20 ) + (v10 y1 + v20 y2 ):

14

(24)

Otra derivacin introducira segundas derivadas de las incgnitas v1 y v2 . Evitamos

esa complicacin exigiendo que
v10 y1 + v20 y2 = 0:

(25)

yp0 = v1 y10 + v2 y20 :

(26)

yp00 = v1 y100 + v10 y10 + v2 y200 + v20 y20 :

(27)

Esto da
As que
Sustituyendo (23), (26) y (27) en (22) y reordenando se llega a
v1 (y100 + P (x)y10 + Q(x)y1 ) + v2 (y200 + P (x)y20 + Q(x)y2 ) + v10 y10 + v20 y20 = g(x):
Como y1 e y2 son soluciones de la homognea asociada a la ecuacin (22) , las dos
expresiones entre parntesis son cero, y se obtiene
v10 y10 + v20 y20 = g(x):

(28)

Teniendo en cuenta (25) y (28) conjuntamente tenemos dos ecuaciones con dos incgnitas v10 y v20
v10 y1 + v20 y2 =
0
0 0
0 0
v1 y1 + v2 y2 = g(x);
cuya solucin es
v10 =

y2 g(x)
W (y1 ; y2 )

v20 =

y1 g(x)
:
W (y1 ; y2 )

(29)

Hay que hacer notar que estas frmulas son legtimas, ya que el Wronskiano de los
denominadores es no nulos por la independencia lineal de y1 y y2 . Integrando cada
ecuacin en (29) obtenemos.
Z
Z
y2 g(x)
y1 g(x)
v1 =
dx y v2 =
dx:
W (y1 ; y2 )
W (y1 ; y2 )
Resumiendo toda la informacin podemos armar que
Z
Z
y2 g(x)
y1 g(x)
yp = y1
dx + y2
dx
W (y1 ; y2 )
W (y1 ; y2 )
es la solucin particular de (20) que buscbamos.

15

Ejercicios 2
1. Hallar una solucin general para cada ED lineal, aplicando el mtodo de los
coecientes indeterminados
00

6y = e4x

00

6y = e2x

00

(a) y + y
(b) y + y

(c) y + y
6y = 5te2t ;
Sol : y = c1 e2t + c2 e
00

(d) 2y + y =

1 + 2e

3t
1
x
2

1 2
t
2

1
t
5

e2t

;
1

1
x
2

+ 2x2 9x 2xe 2 x
Sol : y = c1 + c2 e
p
00
(e) y + 3y = 12 sen( 3x);
p
p
p
Sol : y = c1 cos( 3x) + c2 sen( 3x) 4p1 3 t cos( 3x)
(f) y

00

4y + 13y = 3 cos(3x)
00

(g) 2y + y

6y = senx + 3e

Sol : y = c1 e

2x

+ c2 e

3
x
2

2x

8
senx
65

1
65

cos x

3
xe 2x
7

00

(h) y + 9y = 2sen(3x) + senx;

Sol : y = c1 cos(3x) + c2 sen(3x)
(i) y

00

(j) y

00

10

2y = 12x

1
x cos(3x)
3

+ 12 senx

2y + 5y = e senx

00

(k) y + y + y = xsenx
00

(l) y + 4y = 4sen(2x);
Sol : y = c1 cos(2x) + c2 sen(2x)
00

x cos(2x)

(m) y + y
12y = 2x + e3x ;
1
Sol : c1 e3x c2 e 4x 61 x 72
+ 17 xe3x
2. Hallar la solucin general de la homognea asociada y determine la forma, aplicando el mtodo de los coecientes indeterminados, de una solucin particular
de cada ecuacin. (No evale los coecientes)
00

(a) 2y + 3y = 3 + e

3
t
2

00

(b) y + 9y = 4t3 sen(3t)

00

(c) y + y
00

6y = 5te2t
0

(d) y + 2y + 2y = 8t2 e t sent

16

3. Resuelva por el mtodo de variacin de parmetros el PVI

y 00
Sol: (x4

x2

4y 0 + 4y = 12x2

6x e2x ; y(0) = 1; y (0) = 0:

2x + 1) e2x

4. Hallar una solucin general de cada ED lineal aplicando el mtodo de variacin

de parmetros
(a) y 00 + 4y 0 + 4y = e 2x ln x
Sol : y = c1 + c2 x + x2

1
2

ln x

3
4

2x

(b) 4y 00 + 4y 0 + y = x 2 e 2
Sol : y = c1 + c2 x 14 (1 + ln x) e
(c)

3y 0 + 2y = (1 + e x )

(d) y 00 + 2y 0 + y = (ex
Sol : y = (c1 + c2 x

1) 2
ln (1

e x )) e

(e) y 00 + 2y 0 + y = (ex + 1) 2
Sol : y = (c1 + c2 x + ln (1 + e x )) e

1
x
2

Vibraciones mecnicas

Regresaremos al sistema masa-resorte descrito en la introduccin de este tema y

analizaremos su movimiento con ms detalle. La ecuacin que lo describe es
d2 x
dx
m 2 + b + kx = Fe (t);
dt
dt

(30)

donde m > 0 es la masa del cuerpo sujeta en uno de los extremos del resorte, b 0
es la constante de amortiguaiento, k > 0 es la rigidez del resorte y Fe (t) es la fuerza
externa que acta sobre el sistema.

5.1

Movimiento libre no amortiguado

Primero nos centraremos en un caso sencillo, cuando b = 0 y Fe (t) = 0; llamado el

movimiento libre sin amortiguamiento. Entonces la ecuacin (30) se reduce a
d2 x
+ kx = 0:
dt2
q
k
Dividiendo por m y haciendo ! = m
, la ecuacin anterior se puede escribir como
m

d2 x
+ !2x = 0
dt2
17

(31)

Esta ecuacin describe el movimiento libre no amortiguado de un cuerpo. La ecuacin

auxiliar de (31) es r2 + ! 2 = 0; que tiene races complejas conjugadas !i: Por lo
tanto, su solucin general es
(32)

x(t) = c1 cos(!t) + c2 sen(!t)

Dos condiciones iniciales asociadas con la ecuacin (31) son x(0) = a, la cantidad
0
de desplazamiento inicial, y x (0) = b; la velocidad inicial del cuerpo. Por ejemplo,
si a > 0 y b < 0; el cuerpo parte de un punto debajo de la posicin de equilibrio con
una velocidad hacia arriba. Si a < 0 y b = 0; el cuerpo se suelta partiendo del reposo
desde un punto jaj unidades arriba de la posicin de equilibrio.
Para analizar el movimiento descrito por la ecuacin (31) conviene escribir su
solucin general (32) de la forma
(33)

x(t) = Csen(!t + );
con C

0: Aplicando identidades trigonomtricas obtenemos que

Csen(!t + ) = C (sen(!t) cos( ) + sen( ) cos(!t))
= (Csen( )) cos(!t) + (C cos( )) sen(!t)

Luego, como (32) y (33) son iguales, entonces

C1 cos(!t) + C2 sen(!t) = (Csen( )) cos(!t) + (C cos( )) sen(!t)
Por lo tanto
C1 = Csen( )

C2 = C cos( )

C1
C

cos( ) =

Esto es
sen( ) =

C2
c

donde
C12 + C22
C2
C12 + C22
1 =
2
qC

sen2 ( ) + cos2 ( ) =

C =
y

8
>
arctag cc12 ;
>
>
>
>
>
>
+ arctag cc21 ;
>
>
>
>
<
2 + arctag cc12 ;
=
>
0;
>
>
>
>
;
>
>
>
>
;
>
>
: 32
;
2

18

c21 + c22

si C1 > 0 y C2 > 0
si C2 < 0
si C1 < 0 y
si C1 = 0
si C1 = 0
si C1 > 0
si C1 < 0

C2 > 0
y C2 > 0
y C2 > 0
y C2 = 0
y C2 = 0

Es evidente de (33) que el movimiento de una masa en un sistema libre de amortiguamiento es un onda sinusoidal, o lo que se llama un movimiento armnico
simple (ver Figura 2). La constante C es la amplitud del movimiento y ; que
vara de 0 a 2 ; es el ngulo fase. El movimiento es peridico con periodo T = 2! ;
y frecuencia natural 2! : El periodo se mide en unidades de tiempo y la frecuencia
natural en periodos (o ciclos) por unidades de tiempo. La constante ! es la frecuencia circular o angular para la funcin seno en (33) y se mide en radianes por
unidad de tiempo.

C
Periodo
2/

-/

t
2/-/

Figura 2: Movimiento libre y no amortiguado

Ntese que la amplitud y el ngulo fase dependen de las constantes C1 y C2 ; que
a su vez quedan determinadas por la posicin y velocidad inicial de la masa. Sin
embargo, el periodo y la frecuencia dependen slo de k y m y no de las condiciones
iniciales.
Las unidades con las cuales trabajaremos en los problemas que involucran vibraciones mecnicas se resumen en la siguiente tabla
Unidades
Longitud
Masa
Tiempo
Velocidad
Aceleracin
Fuerza
Gtavedad

cgs
cm
g
s
cm=s
cm=s2
dina
980 cm=s2

mks
m
kg
s
m=s
m=s2
Newton (N )
9; 8 m=s2

fps
Pie (f t)
slug
s
f t=s
f t=s2
Libra (lb)
32 f t=s2

Ejercicios 3
1. Se ja una masa de 81 de Kg: a un resorte con rigides k = 16 N=m: La
masa se desplaza 21 m: hacia abajo de la posicin de equilibrio y luego se le
19

imparte una velocidad hacia abajo de 12 m=s: Despreciando cualquier fuerza

de amortiguamiento externa que pueda esar presente, determinar la ecuacin
del movimiento de la masa, junto con su amplitud, periodo y frecuencia natural.
Despus de cunto tiempo pasa la masa por primera vez, por la posicin de
equilibrio?
2. Una masa de 1 Kg: suspendida de un resorte lo estira 3; 5 cm. Si la masa se
desplaza 7 cm: por debajo de la posicin de equilibrio y se la aplica una velocidad hacia abajo de 7 cm=s: Establezca una ecuacin diferencial y condiciones
iniciales que describan el movimiento. Encuentre la posicin y velocidad de la
masa en cada tiempo t:Encuentre la amplitud, perodo, ngulo fase y frecuencia del movimiento. Determina la posicin y velocidad 1 s: despus de soltar la
masa.
3. Una masa de 30 gr: se une a un resorte. En equilibrio el resorte se alargo20
cm: El resorte se desplaza hacia abajo otros 10 cm: y se suelta. Establezca
la ecuacin diferencial para el movimiento y resulvala para determinar el
movimiento resultante ignorando las fuerzas externas y de amortiguamiento.
4. Una masa de 6 gr: se une a un sistema de masa-resorte con una constante de
resorte de 30gr=s2 : Cules deben ser las condiciones iniciales para obtener una
respuesta con amplitud 3 y ngulo de fase de 4 ? Suponga que el movimiento
es libre y no amortiguado.
5. Una masa de 3 Kg: est unida a un resorte con rigidez K = 48 N=m: La masa
se desplaza 1=2 m: hacia arriba del punto de equilibrio y recibe una velocidad
de 2m=s: hacia abajo. La fuerza de amortiguamiento es despreciable. Determine la ecuacin del movimiento de la masa, junto con su amplitud, perodo y
frecuencia. cunto tiempo despus de la liberacin pasa la masa por primera
vez por la posicin de equilibrio?
6. Una masa de 2 Kg: est unida a un resorte con rigidez k = 50 N=m:La masa se
desplaza 1=4 m hacia arribas de la posicin de equilibrio y recibe una velocidad
de 1 m=s: hacia arriba. Desprecie la fuerza de amortiguamiento y determine la
ecuacin del movimiento de la masa, junto con su amplitud, perodo y frecuencia. cunto tiempo despus de la liberacin pasa la masa por primera vez por
la posicin de equilibrio?
7. Una masa unida a un resorte oscila con un periodo de 3 s: Despus de agregar
2 Kg:; el perodo se convierte en4 s: Si se desprecian las fuerzas de amortiguamiento o externas, determine cunta masa se adjunto originalmente al
resorte.

20

5.2

Movimiento libre amortiguado

En la mayora de las aplicaciones de las vibraciones mecnicas existe cierto tipo de

fuerza de friccin o de amortiguamiento que afecta las vibraciones. Esta fuerza puede
deberse a un medio que rodea el sistema o bien un amortiguador. Ahora estudiaremos
los efectos del amortiguamiento sobre las vibraciones libres, de modo que la ecuacin
(30) se generaliza a
d2 x
dx
(34)
m 2 + b + kx = 0
dt
dt
La ecuacin (34) se puede escribir como
d2 x
dx
+ !2x = 0
+2
2
dt
dt

(35)

donde

b
k
y !2 =
m
m
La ecuacin (35) describe el movimiento libre amortiguado de un cuerpo. Su ecuacin
auxiliar es
r2 + 2 r + ! 2 = 0
(36)
2 =

y las races correspondientes son

p
2
r1 =
+

!2

r2 =

!2

Ahora distinguiremos tres casos posible que dependen del discriminante

!2:

Caso 1. Movimiento sobreamortiguado ( 2 ! 2 > 0) Existe dos races

reales distintas de la ecuacin caracterstica, por lo tanto, la solucin general de la
ecuacin (35) es
x(t) = c1 er1 t + c2 er2 t
(37)
p
2
! 2 < ; por tanto tambin r1 < 0: As
Es claro que r2 < 0 y que
l{m x(t) = 0

t!+1

Adems, como
0

x (t) = c1 r1 er1 t + c2 r2 er2 t =1 er1 t c1 r1 + c2 r2 e(r2

r1 )t

vemos que la derivada es idnticamente nula cuando c1 = c2 = 0; o se anula a lo ms

en un valor de t que es solucin de la ecuacin
c1 r1 + c2 r2 e(r2

21

r1 )t

Si se ignora la solucin trivial x(t) = 0, implica que x(t) tiene a lo ms un

mximo o un mnimo local para t > o. Por tanto x(t) no oscila. Esto deja solo tres
posibilidades para el movimiento de x(t); como se muestra en la Figura 3, de acuerdo
con las condiciones iniciales.

Hay un
mximo

No hay mximo
ni mnimo

Hay un
Mnimo

Figura 3: Movimiento sobreamortiguado

Caso 2. Movimiento crticamente amortiguado ( 2 ! 2 = 0) Existe una
nica raz de la ecuacin caracterstica, por tanto la solucin general de la ecuacin
(35) es
x(t) = (c1 + c2 t) e t
Para comprender el movimiento descrito por x(t); primero consideraremos el comportamiento de x(t) cuano t ! +1: Por la ragla de LHopital
l{m x(t) = l{m

t!+1

(Recuerde que

t!+1

c2
c1 + c2 t
= l{m
=0
t
t!+1
e
e t

> 0). A continuacin, como

0

x (t) = (c2

c1

c2 t) e

vemos de nuevo que una solucin no trivial puede tener a lo ms un mximo o un

mnimo local para t > 0; de modo que el movimiento no es oscilatorio. As se tiene un
movimiento crticamente amortiguado. Grcamente, los movimientos crticamente
amortiguados son similares a los movimientos sobreamortiguados.
Caso 3. Movimiento subamortiguado u oscilatorio ( 2 ! 2 = 0) En este
caso existe dos races complejas conjugadas de la ecuacin caracterstica. Estas races
son
i ; donde
p
2
=
y
= !2
22

Por tanto, a solucin general de la ecuacin (35) es

x(t) = e t (C1 cos( t) + C2 sen( t))
Como en el caso del movimiento libre no amortiguado, podemos expresar esta
solucin en la forma alternativa
p
2
x(t) = Ce t sen
!2
t+

Ahora es evidente que x(t) es el producto de un factor de amortiguamiento exponencial

ce t
y un factor sinosoidal
sen

p
!2

t+

que produce el movimiento oscilatorio. Debido a que el factor sinusoidal vara entre
1 y 1 con periodo p 22 2 ; la solucin x(t) vara entre ce t y ce t con cuasiperi!

odo
P =p

2
!2

y cuasifrecuencia p1 : Adems como a =

< 0; el factor exponencial tiende a cero
cuando t ! +1: En la Figura 4 aparece la grca de una solucin tpica. El sistema se
llama subamortiguado, porque no hay un amortiguamiento suciente ( es demaciado
pequeo) para evitar que el sistema oscile.
c

-c

Figura 4: Movimiento subamortiguado u oscilatorio

Es fcil ver que cuando ! 0; el factor de amortiguamiento tiende a la constante
c y la cuasifrecuencia tiende a la frecuencia natural del movimiento libre subamortiguado correspondiente. Por otra parte, los valores de t donde la grca de x(t) toca
a las curvas exponenciales ce t no siempre son iguales, pero son cercanos, a los
mismos valores de t donde x(t) alcanza, los valores mximos y mnimos.
23

Ejercicios 4
1. Una masa de 1=4 Kg: est unida a un resorte con una rigidez de4 N=m: La
constante de amortiguamiento para el sistema es de 1 N s=m: Si la masa se
desplaza 1=2 m: hacia arriba y recibe una velocidad inicial de 1 m=s: hacia arriba, determine la ecuacin del movimiento.Cul es el mximo desplazamiento
que alcanzar la masa?
2. Una masa de 1=8 Kg: se une a un resorte con rigidez 16 N=m: La constante
de amortiguamiento para el sistema es 2 N s=m: Si la masa se mueve 3=4
m hacia arriba de la posicin de equilibrio y recibe una velocidad inicial hacia
arriba de 2 m=s:, determine la ecuacin del movimiento del cuerpo y de su
factor de amortiguamiento, cuasiperiodo y cuasifrecuencia.
3. Una masa de 20 Kg: se une a un resorte con rigidez 200 N=m. La constante
de amortiguamiento para el sistema es 140 N s=m: Si la masa se mueve 25
cm:hacia abajo de la posicin de equilibrio y recibe una velocidad inicial hacia
arriba de 1 m=s:; cundo regresar por primera vez a su posicin de equilibrio?
4. Una masa de 2 Kg. se une a un resorte
p con rigidez 40N=m: La constante de
amortiguamiento para el sistema es 8 5 N s=m: Si la masa se mueve 10 cm
hacia abajo de la posicin de equilibrio y recibe una velocidad inicial hacia abajo
de 2 m=s:, Cul ser el desplazamiento mximo con respecto de la posicin de
equilibrio?
5. Una masa de1=4 Kg: se une a un resorte con rigidez 8 N=m: La constante de
amortiguamiento para el sistema es 1=4 N s=m: Si la masa se mueve 1 m:
hacia arriba de la posicin de equilibrio y se libera, Cul ser el desplazamiento
mximo hacia abajo?
6. Una masa de 1 Kg: se une a un resorte con rigidez 100 N=m: La constante de
amortiguamiento para el sistema es 0; 2 N s=m: Si la masa se empuja hacia
abajo de la posicin de equilibrio con una velocidad de 1 m=seg:, Cundo
alcanzar su desplazamiento mximo hacia abajo?
7. Una masa de 1=4 Kg: se une a un resorte con rigidez 8 N=m: La constante de
amortiguamiento para el sistema es 2 N s=m: Si la masa se mueve 50 cm:
hacia arriba de la posicin de equilibrio y recibe una velocidad inicial hacia
arriba de 2 m=s:, Cundo alcanzar su desplazamiento mximo hacia arriba?
8. Un resorte con una constante de 8g=seg 2 tiene un objeto unido que lo alarga 245
cm. El coeciente de amortiguamiento es 8 g=seg: En el tiempo t = 0; la masa
se encuentra en posicin de equilibrio y tiene una velocidad de 3 cm=seg: hacia
abajo. Suponga que el movimiento es libre y: (a) Establezca la ecuacin diferencial que describe el movimiento, (b) hallar la solucin de la ecuacin que satisfaga las condiciones iniales consideradas, (c) establezca el tipo de movimiento,
24

(d) bosqueje la representacin grca del movimiento del objeto, (d) determine
que tanto se aleja el objeto de su posicin de equilibrio. (Sol : x(t) = 3te 2t ; el
objeto se aleja o,55 cm.)
9. Un resorte con una constante de 2 kg=seg 2 tiene un objeto unido con una
masa de 10 kg. El coeciente de amortiguamiento es 9 kg=seg: En el tiempo
t = 0 el cuerpo se encuentra en la posicin de equilibrio y tiene una velocidad de 5 m=seg: hacia abajo. Suponga que el movimiento e libre y amortiguado y: (a) Establezca el PVI que describe el movimiento, (b) hallar la
posicin del cuerpo en cada tiempo t, (c) bosqueje la representacin grca del
movimiento, (d) encuentre que tan lejos se mueve la masa de la posicin de
2
1
equilibrio. (Sol : x(t) = 50 e 5 t e 2 t ; la masa se mueve 4,095 m: de la
posicin de equilibrio.)
10. Un resorte con una constante de 226 gr=seg 2 tiene un objeto unido con una
masa de 25 gr. El coeciente de amortiguamiento es 10 gr=seg: En el tiempo
t = 0 el cuerpo se encuentra 20 cm: por debajo de la posicin de equilibrio y
tiene una velocidad de 41 cm=seg: hacia abajo. Suponga que el movimiento e
libre y amortiguado y: (a) Establezca el PVI que describe el movimiento, (b)
hallar la posicin del cuerpo en cada tiempo t, (c) bosqueje la representacin
grca de la funcin posicin del cuerpo, (d) determine el tiempo que tarda el
cuerpo en pasar por primera vez por la posicin de equilibrio. (Sol : x(t) =
1
e 5 t (20 cos(3t) + 15sen(3t)) ; el cuerpo tarda 0,74 seg en pasar por primera vez
por la posicin de equilibrio.)

5.3

Movimiento forzado

Ahora tomaremos en cuenta una fuerza externa Fe (t), que acta sobre una masa
oscilatoria en un resorte. Esto es, investigaremos la EDO
m

dx
d2 x
+
b
+ kx = Fe (t);
dt2
dt

(38)

que puede escribirse de la forma

d2 x
dx
+2
+ ! 2 x = F (t);
2
dt
dt
donde

(39)

b
k
Fe (t)
; !2 =
y F (t) =
m
m
m
Es muy comn que la fuerza externa de un sistema masa resorte sea armnica
simple.
Fe (t) = F0 cos( t) o Fe (t) = F0 sen( t);
(40)
2 =

25

donde F0 es la amplitud de la fuerza peridica y es su frecuencia angular.

A continucin mostraremos un ejemplo de un movimiento forzado con amortiguamiento.
Ejemplo 5.1 Consideremos un sistema msa-resorte formado por una masa de 1=5
Kg:;unida a un resorte con rigidez de 2 N=m: La masa parte del reposo a 1=2 m:
por debajo de su posicin de equilibrio. El movimiento tiene una constante de amortiguamiento de 1; 2 N s=m: y est impulsada por una fuerza externa f (t) = 5 cos(4t)
que se inicia cuando t = 0:
Solucin. El PVI que modela el problema es
1 d2 x 6 dx
+ 2x = 5 cos(4t);
+
5 dt2
5 dt

1
x(0) = ;
2

x (0) = 0

La solucin del PVI es

x(t) = xc (t) + xp (t)
86
25
50
38
cos(t)
sen(t)
cos(4t) + sen(4t)
= e 3t
51
51
102
51
2p
25 p
=
2210e 3t sen (t + 2; 73) +
17sen (4t + 6; 04)
51
102
Ntese que la solucin complementaria xc (t) del ejemplo anterior tiene la propiedad
de que
2p
2210e 3t sen (t + 2; 73) = 0
t!+1
t!+1 51
Como xc (t) se vuelve insignicante cuando t ! +1, se dice que es un trmino
transitorio o solucin transitoria. As, cuando el tiempo es grande, los desplazamientos de la masa del problema anterior son muy bien aproximados por la solucin
particcular xp (t): Esta lmina funcin se llama tambin solucin de estado estable,
de estado estacionario o de estado permamente.
Cuando F (t) es una funcin peridica como (40), la solucin general de la ecuacin
(39) es de la forma
l{m xc (t) = l{m

x(t) = parte transitoria + parte estable

Cuado se ejerce una fuerza peridica y no existe fuerza de amortiguamiento, no
hay parte transitoria en la solucin del problema. Si la fuerza peridica que se ejerce
tiene una frecuencia que es igual o casi igual a la de las vibraciones no amortiguadas
libres, se puede originar un grave problema en un sistema mecnico oscilatorio

26

Ejemplo 5.2 Resuelva el PVI

x00 + 9x = 5sen( t);

x(0) = 0;

x0 (0) = 0;

j j=
6 3

Solucin. La solucin de la ecuacin homognea asociada a la EDO es

xc (t) = c1 cos(3t) + c2 sen(3t)
Aplicando el mtodo den los coecientes indeterminados obtenemos la solucin particular
5
xp (t) =
2 sen( t)
9
Luego, la solucin general de la EDO es
x(t) = c1 cos(3t) + c2 sen(3t) +

5
9

2 sen(

Aplicando las condiciones iniciales obtenemos c1 = 0 y c2 =

t)
5
3(9

: Por tanto, la

solucin del PVI es

x(t) =

Para

5
2

3 9

sen(3t) + 3sen( t)) ;

j j=
6 3

= 3; la solucin del PVI se dene como

5
!3 3

x(t) = l{m

sen(3t) + 3sen( t)
5
= sen(3t)
2
18
9

Ntese que cuando tn =

n
3

5
t cos(3t)
6

; n = 1; 2; 3; :::;

jx(tn )j =

5n
! +1 si n ! +1
8

Por tanto, los desplazamientos crecen cuando t ! +1: Este fenmeno se llama
resonancia pura. La Figua 5 muestra un movimiento caracterstico de este caso

27

Figura 5: Resonancia pura

La solucin anterior es consecuencia de resolver el PVI
x00 + 9x = 5sen(3t);

x(0) = 0;

x0 (0) = 0

directamente por los mtodos convencionales.

Si la gura anterior describe en realidad los desplazamientos de un sistema de
resorte y masa, este sistema se destruira. Las oscilaciones grandes de la masa forzarn
al resorte a rebasar su lmite elstico. Se podra decir que el modelo descrito es
irreal por completo, porque no tiene en cuenta los efectos retardantes de las siempre
presentes fuerzas de amortiguamiento. No se puede tener resonancia pura cuando
se considera un amortiguamiento mnimo, pero si se pueden desarrollar amplitudes
grandes e igualmente destructivas de vibracin, pero acotadas cuando t ! +1:

Ejercicios 5
1. Un objeto con na masa de 2 kg. est suspendido del extremo de un resorte
que tiene una constante de 128 kg=seg 2 , suponiendo que el movimiento es no
amortiguado y que se aplica una fuerza externa dada por f (t) = 40 cos(8t) para
t > 0 : (a) encuentre la posicin del objeto en cada tiempo t si ste se mueve 16 m:
hacia abajo de la posicin de equilibrio y se suelta, (b) dar una interpretacin
fsica de lo que sucede con el movimiento del objeto cuando t crece.
Sol : x(t) = 61 cos(8t) + 54 tsen(8t)
2. Una masa de 8 kg: se une a un resorte que cuelga desde el techo, haciendo que
el resorte se estire 1.96 m: hasta llegar al reposo en equilibrio. En el instante
t = 0 se aplica una fuerza externa f (t) = cos(2t) N: al istema. La constante de
amortiguamiento del sistema es 3 N seg=m: Determine la solucin de estado
estable para el sistema.
Sol : xp (t) =

2
25

cos(2t) +

3
sen(2t):
50

3. La respuesta de un sistema sobreamortiguado a una fuerza constante queda

descrita mediante la ecuacin
00

2x + 8x + 6x = 18:
0

Si el sistema parte del reposo (x(0) = x (0) = 0), calcule x(t) y bosqueje su
representacin grca. Cul es el valor de l{m x(t)?
t!+1

1
2

Sol : x(t) = (3e

3t

9e ) + 3:

28

Circuitos en serie LRC

Si i(t) representa la corriente en un circuito elctrico en serie LRC, de acuerdo con

la segunda ley de Kirchho, la suma de las cadas de voltaje a travs del inductor,
resistor y capacitor es igual al voltaje aplicado al circuito; esto es
L

di
1
+ Ri + q = E(t)
dt
C

Pero i(t) = dq
, de manera que la ecuacin anterior se transforma en la ecuacin
dt
diferencial lineal de segundo orden
1
d2 q
dq
(41)
+ R + q = E(t)
2
dt
dt C
La nomenclatura que se emplea en el anlisis de circuitos, es similar a la que se
usa en los sistemas de resorte y masa.
L

Sistema mecnico
Masa m
constante de amortiguamiento
Constante del resorte k
Posicin x
Fuerza externa f

Sistema elctrico
Inductancia L
Resistencia R
Recproco de la capacitancia
Carga q
Fuerza electromotriz E

1
C

En consecuencia, la mayora de los resultados deducidos para los sistemas mecnicos pueden aplicarse de inmediato a circuitos elctricos. El hecho de que la misma
ecuacin diferencial sirva como modelo matemtico para sistemas fsicos tan diferentes
es una ilustracin poderosa del papel unicador de las matemticas en la investigacin
de fenmenos naturales.
Si E(t) = 0, las vibraciones elctricas del circuito se llaman libres. Como la
ecuacin auxiliar de la ecuacin homognea asociado a la ecuacin anterior es
Lm2 + Rm + 1=C = 0;
habrn tres formas de la solucin cuando R 6= 0, dependiendo del valor del discriminante R2 4L
. Se dice que el circuito es
C
4L
C

1. Sobreamortiguado si R2

>0

2. Crticamente amortiguado si R2
3. Subamortiguado si R2

4L
C

<0

29

4L
C

=0

En cada uno de los tres casos, la solucin de la ecuacin lineal contiene el factor
R
e 2L , as que q(t) ! 0 cuando t ! +1. Cuando E(t) = 0 y R = 0, se dice que el
circuito es no amortiguado, y las vibraciones elctricas no tienden a cero cuando t
aumenta sin lmite; la respuesta al circuito es armnica simple.
Cuando hay un voltaje E(t) aplicado en el circuito, se dice que las vibraciones
elctricas son forzadas. Cuando R 6= 0, la solucin de la homognea asociada se
llama solucin transitoria. Si E(t) es peridico o constante. La solucin particular
de la ecuacin lineal, es una solucin de estado estable, solucin permanente
o solucin remanente.

Ejercicios 6
1. Encuentre la carga en el capacitor de un circuito LRC en serie en cuando t =
0; 01s: cuando L = 0:05 h:; R = 2 :; C = 0:01 f:; E(t) = 0 v:; q(0) = 5 c: e
i(0) = 0 A: Determine la primera vez que la carga en capacitor es igual a cero.
2. Calcule la carga en el capacitor de un circuito en serie LC cuando L = 0:01 h:;
C = 0:1 f:; E(t) = 100sent v:; q(0) = 0 c: e i(0) = 0 A: Cul es la carga en el
capacitor cuando ha transcurrido mucho tiempo?
Sol : q(t) = 23 e

10t

(cos(10t) + sen(10t)) +

3
2

3. Determine la carga y la corriente en un circuito en serie LRC cuando L = 21 h:;

R = 10 :; C = 0:01 f:; E(t) = 150 v:; q(0) = 1 c: e i(0) = 0 A: Cul es la
carga en el capacitor cuando ha transcurrido mucho tiempo?
Sol : q(t) =

100
(10sen(t)
99

sen(10t))

4. Encuentre la carga permanente y la corriente remanete en un circuito LRC en

serie cuando L = 1 h:; R = 2 :; C = 0:25 f: y E(t) = 50 cos t v:

30