UNIVERSIDAD DE CONCEPCION

FACULTAD DE CIENCIAS
FISICAS Y MATEMATICAS
DEPARTAMENTO DE INGENIERIA MATEMATICA

ECUACIONES DIFERENCIALES ORDINARIAS (521218)

Listado N ◦ 5 (Método de Variación de Parámetros y Problemas de Vibraciones).

Problemas a resolver en práctica

1. Resolver aplicando método de variación de parámetros:

ex
a) y 00 (x) − 3y 0 (x) + 2y(x) = , y(0) = 0, y 0 (0) = 0
1 + ex
b) 2y 00 (x) − 4y 0 (x) + 2y(x) = x−1 ex , y(1) = 0, y 0 (1) = 0
d
En lo que sigue, D =
dx
Solución a): La EDO se puede re – escribir en forma de operador como
ex
(D2 − 3 D + 2)[y](x) = ,
1 + ex
la cual ya se encuentra normalizada, y se identifica aquı́
ex
f (x) := .
1 + ex
Aplicando el Principio de Superposición, la solución de la EDO se descompone en
y = yh + yp , donde yp es una solución particular de la EDO, e yh es la solución
general de la EDO homogénea asociada,

(D2 − 3 D + 2)[yh ] = 0.

La ecuación caracterı́stica de la EDO precedente está dada por

λ2 − 3 λ + 2 = 0 ⇐⇒ (λ − 1) (λ − 2) = 0 ⇐⇒ λ1 = 1, λ2 = 2,

de donde, la solución general de la EDO homogénea asociada está dada por

yh (x) = C1 ex + C2 e2 x , ∀ x ∈ R, con C1 , C2 ∈ R.

Para encontrar una solución particular por el método de variación de parámetros,

se propone como solución particular a

yp (x) = A1 (x) ex + A2 (x) e2 x ,

1
donde las funciones A01 y A02 son soluciones del siguiente sistema de ecuaciones,
escrito en forma matricial,

0
   
0
A01 (x) A01 (x)
 x       
e e2 x 1 ex
= ⇐⇒ =
   
A02 (x) A02 (x)
 ex 
ex 2 e2 x 1 2 ex 1 
1 + ex 1 + ex
 
0
A01 (x)
   
1 ex
⇐⇒ =
 
x 0 1
0 e A2 (x)

1 + ex
donde primero se dividió cada ecuación por ex y luego se realizó la operación elemen-
tal por filas f2 ← f2 − f1 . De esta manera, la única ecuación del sistema precedente
es
1 e−x e−x e−2 x
A01 (x) = − = − , A02 (x) = = , x ∈ R.
1 + ex 1 + e−x 1 + ex 1 + e−x

Ahora, integrando con respecto a x las expresiones precedentes, y considerando el

cambio de variable e−x = t, se deduce que
e−x
Z Z
dt
A1 (x) = − −x
dx = = ln(1 + t) + K1 = ln(1 + e−x ) + K1 ,
1+e 1+t

e−2 x
Z Z Z  
t 1
A2 (x) = dx = − dt = − 1 dt = ln(1 + t) − t + K2
1 + e−x 1+t 1+t
= ln(1 + e−x ) − e−x + K2 .

Podemos eligir las constantes de integración tales que K1 = K2 = 0 pues buscamos

una solución particular. Por lo tanto,
 
−x x −x −x
yp (x) = ln(1 + e ) e + ln(1 + e ) − e e2 x
= ln(1 + e−x ) ex + ln(1 + e−x ) e2 x − ex .

Ası́, la solución general de la EDO planteada está dada por

y(x) = C1 ex + C2 e2 x + ln(1 + e−x ) ex + ln(1 + e−x ) e2 x − ex

=Cf1 ex + C2 e2 x + ln(1 + e−x ) ex + ln(1 + e−x ) e2 x , ∀ x ∈ R, con C
f1 , C2 ∈ R.

f1 := C1 − 1. Ahora, notemos que la derivada de la solución

donde se ha definido C
general de la EDO original está dada por
x
f1 ex + 2 C2 e2 x − 1 + e + ln(1 + e−x ) ex + 2 ln(1 + e−x ) e2 x ,
y 0 (x) = C ∀ x ∈ R.
1 + e−x

2
Imponiendo las condiciones iniciales, se deduce que
 ( 
y(0) = 0 Cf1 + C2 + 2 ln(2) = 0 f1 = −1 − ln(2)
C
0 ⇐⇒ ⇐⇒
y (0) = 0 f1 + 2 C2 − 1 + 3 ln(2) = 0
C C2 = 1 − ln(2)

Por lo tanto, la única solución de PVI está dada por

y(x) = −[1 + ln(2)] ex + [1 − ln(2)] e2 x + ln(1 + e−x ) ex + ln(1 + e−x ) e2 x , ∀ x ∈ R.

Solución b): La EDO se puede re – escribir en forma de operador y a la vez

normalizada como
ex
(D2 − 2 D + 1)[y](x) = , x > 0,
2x
y se identifica aquı́
ex
f (x) := .
2x
Aplicando el Principio de Superposición, la solución de la EDO se descompone en
y = yh + yp , donde yp es una solución particular de la EDO, e yh es la solución
general de la EDO homogénea asociada,

(D2 − 2 D + 1)[yh ] = 0.

La ecuación caracterı́stica de la EDO precedente está dada por

λ2 − 2 λ + 1 = 0 ⇐⇒ (λ − 1)2 = 0 ⇐⇒ λ1 = λ2 = 1,

de donde, la solución general de la EDO homogénea asociada está dada por

yh (x) = C1 ex + C2 x ex , ∀ x ∈ R, ∀ C1 , C2 ∈ R.

Para encontrar una solución particular por el método de variación de parámetros,

se propone como solución particular a

yp (x) = A1 (x) ex + A2 (x) x ex ,

donde las funciones A01 y A02 son soluciones del siguiente sistema de ecuaciones,
escrito en forma matricial,
   
0 0
A01 (x) A01 (x)
 x       
e x ex 1 x
= ⇐⇒ =
   
x x x 0 x 0
e e + xe A2 (x) e 1 1+x A2 (x) 1
  
2x  2x
  0  0
1 x A1 (x)
⇐⇒ =
 
0 1 A02 (x) 1

2x
3
 donde primero se dividió cada ecuación por ex y luego se realizó la operación elemen-
tal por filas f2 ← f2 − f1 . De esta manera, la única ecuación del sistema precedente
es
1 1
A01 (x) = − , A02 (x) = , x>0
2 2x
Ahora, integrando con respecto a x las expresiones precedentes, se deduce que
Z
1 x
A1 (x) = − dx = − + K1 , x > 0
2 2

Z
1 1
A2 (x) = dx = ln(x) + K2 , x > 0.
2x 2
Sin perdida de generalidad podemos eligir las constantes K1 = K2 = 0. Por lo tanto,
x x
yp (x) = − ex + ln(x) ex .
2 2

Ası́, la solución general de la EDO planteada está dada por

x x x
y(x) = C1 ex + C2 x ex −e + ln(x) ex
2 2
x x x
f2 x e + ln(x) ex ,
= C1 e + C ∀ x > 0, f2 ∈ R,
con C1 , C
2
f2 := C2 − 1/2. Ahora, notemos que la derivada de la solución
donde se ha definido C
general de la EDO original está dada por
x
y 0 (x) = C1 ex + C f2 x ex + ln(x) e + 1 ex + x ln(x) ex ,
f2 ex + C ∀ x > 0.
2 2 2

Imponiendo las condiciones iniciales, se deduce que


1
 C1 = 2
 
e C1 + e C
f2 = 0
 

y(1) = 0

⇐⇒ ⇐⇒
y 0 (1) = 0
f2 + 1 e = 0 f2 = − 1
 
e C1 + 2 e C 2

 C

2

Por lo tanto, la única solución de PVI está dada por

1 x 1 x x
y(x) = e − x e + ln(x) ex , ∀ x > 0.
2 2 2

2. Movimiento libre sin amortiguamiento: Un cuerpo se acopla a un resorte que

cuelga del techo, produciendo un estiramiento de 12 [cm]. De repente el cuerpo es

4
soltado desde un punto a 32 [cm] por debajo de la posición de equilibrio con una ve-
locidad inicial hacia arriba de 34 [cm/s]. Despreciando el amortiguamiento del medio,
determine la ecuación de movimiento del cuerpo, indicando su amplitud, periodo y
frecuencia. ¿En qué instante pasará el cuerpo por su posición de equilibrio, por se-
gunda vez? ¿En qué instante el resorte se elonga (comprime) lo máximo posible, por
primera vez? Considerar g = 980 [cm/s2 ].

Solución: Consideremos un sistema de referencia orientado positivamente en la

dirección de la gravedad, con origen en el punto de equilibrio del centro de masa del
cuerpo. Además, se definen las variables de trabajo

t : instante de tiempo, medido en segundos,

x(t) : posición del centro de masa en centı́metros, en el instante t.

Del enunciado, se obtiene que el cuerpo produce un estiramiento l = 21 [cm] al

equilibrio. Usando la condición de equilibrio k l = m g, la frecuencia angular natural
de vibración del sistema verifica que

k g 980 [cm/s2 ] √
ω02 = = = 1 = 1960 [(rad/s)2
] ⇐⇒ ω0 = 14 10 [rad/s].
m l 2
[cm]

2 4
A partir de aquı́, como x0 = [cm] y v0 = − [cm/s] (el signo menos refleja el hecho
3 3
que la velocidad apunta hacia arriba), se deduce que x(t) satisface el siguiente PVI
(se recuerda que se está ante un movimiento libre sin amortiguamiento),



 x00 (t) + 1960 x(t) = 0, t ≥ 0,

 2
x(0) = x0 = [cm], (1)
3

 0 4

 x (0) = v0 = − [cm/s]
3

La solución de este PVI está dada, para todo t ≥ 0, por

" √ #
v0 2 √ 10 √
x(t) = x0 cos(ω0 t) + sen(ω0 t) = cos(14 10 t) − sen(14 10 t) [cm].
ω0 3 105

A partir de aquı́, se deduce que la amplitud (A), periodo (T ) y frecuencia (f ) del

movimiento están dadas por
s
2
r √
v0 1 982 ω 0 7 10 1 π
A = x20 + 2 = [cm], f = = [Hz], T = = √ [s].
ω0 21 5 2π π f 7 10

5
Para saber los instante donde el cuerpo pasará por su posición de equilibrio, se debe
resolver la ecuación
√
2 √ 10 √
x(t) = 0, t ≥ 0 ⇐⇒ cos(14 10 t) − sen(14 10 t) = 0, t ≥ 0
3 105
√ 2 × 105 √
⇐⇒ tg(14 10 t) = √ = 7 10, t ≥ 0
3 10
√ √
⇐⇒ 14 10 tk = Arctg(7 10) + k π, k ∈ N
" √ #
Arctg(7 10) π
⇐⇒ tk = √ +k √ [s], k ∈ N.
14 10 14 10

Ası́, el instante en que pasará el cuerpo por su posición de equilibrio

√ por primera vez
Arctg(7 10)
(sin considerar la posición inicial) está dado por t0 = √ [s] ≈ 0,034 [s]
14 10
y el instante en que pasará por segunda vez está dado por
" √ #
Arctg(7 10) π
t1 = √ + √ [s] ≈ 0,105 [s].
14 10 14 10

Por otro lado, dado que al cuerpo se le ha inyectado una velocidad inicial hacia arri-
ba, la máxima elongación/compresión ocurre cuando la velocidad es nula y además,
primero se da la compresión máxima. Recordemos que la velocidad en todo instante
de tiempo t ≥ 0 está dada por
" √ #
28 10 √ 4 √
x0 (t) = − sen(14 10 t) − cos(14 10 t) [cm/s], ∀ t ≥ 0,
3 3

se debe resolver la ecuación,

√
28 √10 4 √
x0 (e
t) = 0, t ≥ 0 ⇐⇒ −
e t ) − cos(14 10 e
sen(14 10 e t≥0
t ) = 0, e
3 3 √
√ 4 10 e
⇐⇒ tg(14 10 e t)=− √ =− , t≥0
28 10 70
√ !
√ 10
⇐⇒ 14 10 ten = −Arctg + n π, n ∈ N∗
70
√ √ ! √
10 10 n 10 π
⇐⇒ ten = − Arctg + [s], k ∈ N∗ .
140 70 140

Notamos que puesto que para n = 0 se obtiene un t̃0 negativo, la solución corres-
pondiente no se toma en cuanto. De esta manera, los primeros instantes de máxima

6
 compresión y elongación están dados, respectivamente, por
" √ √ ! √ #
10 10 10 π
tcom,1 = te1 = − Arctg + [s] ≈ 0,070 [s],
140 70 140
" √ √ ! √ #
10 10 10 π
telon,1 = te2 = − Arctg + [s] ≈ 0,141 [s]
140 70 70

Observación: En este apunte, se denota

N = {0, 1, · · · }, N∗ = N − {0}.

3. Movimiento libre con amortiguamiento: Un cuerpo que pesa 19, 6 [N ] se une

a un resorte de 5 [m] de longitud, el cual cuelga de cierto techo. En la posición
de equilibrio, el resorte mide 9, 9 [m]. Si el cuerpo se eleva y se suelta desde el
reposo 2 [m] por arriba de la posición de equilibrio, determine los máximos relativos
del desplazamiento, x(t). Considere que el medio que rodea al sistema ofrece una
resistencia µ = 1 [N · s/m] ¿Cuál es su máximo desplazamiento con respecto del
equilibrio? Considerar g = 9,8 [m/s2 ].

Solución: Usaremos las mismas notaciones y sistema de referencia que en el pro-

blema anterior. Del enunciado se obtiene que el cuerpo produce al equilibrio un
estiramiento l = 9,9 [m] − 5 [m] = 4,9 [m] y la masa está dada por

mg 19,6 [N ]
m= = = 2 [kg].
g 9,8 [m]

Usando la condición de equilibrio k l = m g, tenemos

r s
9,8 [m/s2 ] √
r
k g
ω0 = = = = 2 [rad/s].
m l 4,9 [m]

Se recuerda que se está ante un movimiento libre con amortiguamiento, la EDO que
modela el movimiento del cuerpo es x00 (t) + 2 λ x0 (t) + ω02 x(t) = 0 con
µ 1
λ= = [s−1 ].
2m 4
A partir de aquı́, como x0 = −2 [m] y v0 = 0 [m/s] (se suelta el cuerpo desde el
reposo), se deduce que x satisface el siguiente PVI
 00
 x (t) + 21 x0 (t) + 2 x(t) = 0, t ≥ 0,
x(0) = −2 [m], (2)
0
x (0) = 0 [m/s]


7
La ecuación caracterı́stica asociada a la EDO para x está dada por
√ √
1 1 31 1 31
r2 + r + 2 = 0 ⇐⇒ r1 = − + i, r2 = − − i,
2 4 4 4 4
de donde la solución general de la EDO es
" √ ! √ !#
31 31
x(t) = e−t/4 C1 cos t + C2 sen t [m], ∀ t ≥ 0, con C1 , C2 ∈ R.
4 4
Notemos que la primera derivada de la solución general está dada por
(√ √ ! √ √ !)
31 C 2 − C 1 31 31 C 1 + C 2 31
x0 (t) = e−t/4 cos t − sen t [m/s]
4 4 4 4
Imponiendo las condiciones iniciales, se obtiene el sistema de ecuaciones
 

 C 1 = −2 [m] 
 C1 = −2 [m]
 
√ ⇐⇒ √
31 C 2 − C 1 C 1 2 31
 C2 = √ = − [m]
 

 = 0 [m/s] 
4 31 31
Por lo tanto, la ecuación de movimiento del cuerpo está dada por
" √ ! √ √ !#
31 2 31 31
x(t) = e−t/4 −2 cos t − sen t [m], ∀ t ≥ 0.
4 31 4
Ahora bien, notemos que los máximos relativos del desplazamiento (que está dada
por la función |x|) ocurren cuando la velocidad es nula, pues tanto los puntos de
máximo relativo y de mı́nimo relativo de x son puntos de máximo relativo de |x|; o
bien, en t = 0 (esto es porque al optimizar una función de clase C 1 en un intervalo
cerrado del tipo [a, +∞[, t = a es candidato a un punto de extremo relativo de la
función). Recordemos que la velocidad en todo instante de tiempo t > 0 está dada
por (usando los valores de C1 y C2 ya conocidos)
√ √ !
16 31 31
x0 (t) = e−t/4 sen t [m/s], ∀ t ≥ 0,
31 4
de donde, se debe resolver la ecuación,
√ √ !
16 31 31
x0 (t) = 0, t ≥ 0 ⇐⇒ e−t/4 sen t = 0, t≥0
31 4
√ !
31
⇐⇒ sen t = 0, t ≥ 0
4
√
31
⇐⇒ tk = k π, k ∈ N
4 √
4 k 31 π
⇐⇒ tk = [s], k ∈ N.
31
8
 siendo estos los puntos de extremos relativos de x. Finalmente, se concluye que
el máximo desplazamiento ocurre en el instante t0 = 0, pues la amplitud va de-
creciendo a medida que avanza el tiempo. Esto no nos sorprende, pues al poseer
el cuerpo velocidad inicial nula y ser el sistema amortiguado, necesariamente el
máximo desplazamiento ocurre en t = 0.

4. Movimiento libre con amortiguamiento: Un cuerpo de masa m se une a un

resorte de rigidez k que cuelga del techo. Se supone que la resistencia al movimiento
del medio que rodea al cuerpo es proporcional a la velocidad de este, con una
constante de amortiguamiento µ.
(a) Muestre que en el régimen de vibraciones subamortiguadas, la posición
√ vertical
−λt 2 2
x(t) del centro de masa del cuerpop esta dada por x(t) = Ae sen( ω − λ t +
θ), donde λ = µ/(2m) y ω0 = k/m. Determine la amplitud A y fase θ en
términos de λ, ω0 y de la posición y velocidad iniciales x0 y v0 del cuerpo.
(b) Determine los instantes de tiempo t0 , t1 , . . ., para los cuales la curva gráfica de
x(t) toca las curvas exponenciales ±Ae−λt .
(c) Muestre que para t > 0, x(t) alcanza sus máximos y mı́nimos relativos a
instantes de tiempos t∗k distintos de tk , con k = 0, 1, ....
(d) Si la velocidad inicial √
del cuerpo v0 es nula, pero no ası́ su posición inicial x0 ,
∗
muestre que tk = kπ/ ω 2 − λ2 , donde k ≥ 0 es un entero.

Solución a): Bajo las condiciones del enunciado, x(t) es solución del siguiente PVI,
 00
 x (t) + 2 λ x0 (t) + ω02 x(t) = 0, t ≥ 0,
x(0) = x0 , ,
0
x (0) = v0


siendo r
µ k
λ= , ω0 = .
2m m
La única solución del PVI está dada por
" #
v + λ x
q q
0 0
x(t) = e−λ t x0 cos( ω02 − λ2 t) + p 2 sen( ω02 − λ2 t) , t ≥ 0.
ω0 − λ 2

Esta solución se puede re-escribir como

√
x(t) = Ae−λt sen( ω 2 − λ2 t + θ)
donde la amplitud de la función sinusoidal que forma parte de x está dada por
s
(v0 + λ x0 )2
A = x20 + .
ω02 − λ2

9
Para determinar el ángulo de fase θ, dado que
q q q
A sen( ω02 − λ2 t + θ) = A sen(θ) cos( ω02 − λ2 t) + A cos(θ) sen( ω02 − λ2 t),

comparando con x, se deduce que


 A sen(θ) = x0
v0 + λ x0 (3)
 A cos(θ) = p 2
ω0 − λ2

luego, si v0 + λ x0 6= 0, entonces θ es tal que

p
x0 ω02 − λ2
tg(θ) = ,
v0 + λ x0
siendo elegido de tal forma que sea consistente con (3).

Si v0 + λ x0 = 0, entonces cos(θ) = 0, y es posible elegir,

 π
 2 , x0 > 0,


θ=
 − π , x0 < 0


2

Solución b): Para que el gráfico de x corte al gráfico

p de las exponenciales t 7→
−λ t 2
±A e , usando el ı́tem a), debe ocurrir que sen( ω0 − λ2 t + θ) = ±1. Ahora,
q q
2 2
sen( ω0 − λ t + θ) = ±1 ⇐⇒ ω02 − λ2 t + θ = k π + π/2, k ∈ Z
kπ π/2 − θ
⇐⇒ t = p 2 +p 2 , k ∈ Z.
ω0 − λ2 ω0 − λ2

Puesto que t > 0, el entero k debe satisfacer

π
kπ > − + θ.
2

Solución c): Para t > 0, como x es diferenciable, t∗ es un extremo relativo de x sı́,

y sólo sı́, x0 (t∗ ) = 0. Ahora, dado que
h q q q i
0 −λ t 2 2 2 2 2 2
x (t) = A e − λ sen( ω0 − λ t + θ) + ω0 − λ cos( ω0 − λ t + θ) ,

10
se deduce que
p
ω02 − λ2
q
x0 (t∗ ) = 0 ⇐⇒ tg( ω02 − λ2 t∗ + θ) =
λ !
p
2 2
ω − λ
q
⇐⇒ ω02 − λ2 t∗ + θ = Arctg 0
+ k π, k ∈ Z
λ
p !
ω02 − λ2
Arctg −θ
λ kπ
⇐⇒ t∗ = p
2
+p 2 , k ∈ Z,
ω0 − λ2 ω0 − λ2
Puesto que t∗ > 0, el entero k debe satisfacer
p !
2 2
ω0 − λ
k π > −Arctg + θ.
λ

Finalmente, dado que para cada k ∈ Z,

p
tg2 ( ω02 − λ2 t∗k + θ) ω02 − λ2
q
2 2 2 ∗
sen ( ω0 − λ tk + θ) = =
ω02
p
1 + tg2 ( ω02 − λ2 t∗k + θ)
q
2
6= 1 = sen ( ω02 − λ2 tk + θ)

se concluye que si el amortiguamiento no es nulo, esto es, λ > 0, los puntos tk > 0
y t∗k > 0 no coinciden.

Solución d): Si la velocidad inicial del cuerpo es nula, y además su posición inicial
x0 es no nula, del Ãtem a) se deduce que
" #
λ x
q q
0
x(t) = e−λ t x0 cos( ω02 − λ2 t) + p 2 sen( ω02 − λ2 t) ,
ω0 − λ 2

siendo su primera derivada,

!
2
λ x
q q
0 2 0 −λ t
2
x (t) = − x0 ω0 − λ + p 2 e sen( ω02 − λ2 t),
ω0 − λ 2

y ası́, los extremos relativos de x son las soluciones de la ecuación

q
x (t ) = 0, t ≥ 0 ⇐⇒ sen( ω02 − λ2 t∗ ) = 0, t∗ ≥ 0
0 ∗ ∗

kπ
⇐⇒ t∗k = p 2 , k ∈ N.
ω0 − λ2

11
 5. Movimiento forzado sin amortiguamiento: ¿Para qué valores de m el sistema
masa – resorte, modelado por la ecuación diferencial m y 00 (t) + 64y(t) = 15 sen(wt)
manifiesta resonancia, si sen(wt) tiene una frecuencia de 12 [Hz] ?

Solución: Dado que la frecuencia de la fuerza excitante es de f = 12 [Hz], su

frecuencia angular está dada por

ω = 2 π f = 24 π [rad/s].

Además, dado la EDO propuesta, deducimos que la constante de rı́gidez del resorte
es k = 64 [N/m]. Luego la frecuencia angular natural del sistema masa – resorte
está dada por r
64
ω0 = [rad/s].
m
Por lo tanto, el sistema forzado presentará resonancia sı́, y sólo sı́,
r
64 64 1
ω = ω0 ⇐⇒ 24 π [rad/s] = ⇐⇒ m = 2
= [kg] ≈ 0,011 [kg].
m (24π) 9 π2

Problemas propuestos para el estudiante:

1. Usando el método de variación de parámetros, resuelva:

a) y 00 − 2y 0 + 5y = 1/x, y(0) = 0, y 0 (0) = 0;

b) y 000 + 3y 00 + 3y 0 + y = x−3 e−x ;
c) y 00 (x) + 2y 0 (x) + y(x) = e−x ln(x), y(1) = 0, y 0 (1) = 0.
e−x
d ) y 00 (x) + 2y 0 (x) + y(x) = .
1 + x2
2. a) Encuentre las soluciones de la EDO

x3 y (3) (x) − 2x2 y (2) (x) − 2xy (1) (x) + 8y(x) = 0 ∀x > 0

con la forma y(x) = xr con r ∈ R.

b) Encuentre la solución general de la EDO:

x3 y (3) (x) − 2x2 y (2) (x) − 2xy (1) (x) + 8y(x) = 0 ∀x > 0.

c) Encuentre la solución general de la EDO:

x3 y (3) (x) − 2x2 y (2) (x) − 2xy (1) (x) + 8y(x) = x ∀x > 0.

12
3. Determine la solución general de
(9 − 2x)y 00 (x) − 4(x − 5)y 0 (x) + 4y(x) = 4x2 − 36x + 81
sabiendo que el kernel (o núcleo) del operador asociado a la EDO tiene una base
dada por {y1 (x) = x − 5, y2 (x) = e−2x }.
4. Combinar los métodos de anuladores y de variación de parámetros para resolver el
PVI:

3y 00 − 6y 0 + 30y = 15 sin(x) + ex tan(3x)

y(0) = 0
y 0 (0) = 1
5. Movimiento horizontal libre sin amortiguamiento: Una masa de 2[kg] se des-
plaza en el plano horizontal y está unida a un resorte con rigidez k = 50[N/m] sujeto
a una pared. La masa se desplaza (1/4)[m] a la izquierda del punto de equilibrio y
recibe una velocidad de 1[m/s] hacia la izquierda. Desprecie el amortiguamineto y
determine la ecuación de movimiento de la masa junto con su amplitud, periodo y
frecuencia ¿Cuánto tiempo después de su liberación pasa la masa por su posición
de equilibrio por primera vez?
6. Movimiento forzado sin/con amortiguamiento: Un cuerpo que pesa 490 [dina]
queda suspendido de un resorte, alargándolo 49 [cm]. De pronto, el cuerpo es des-
plazado 3 [cm] por debajo del punto de equilibrio y es liberado. En este instan-
te se aplica al sistema una fuerza externa F (t) := 8 cos(6t) [dina]. Considerando
g = 980 [cm/s2 ], determine y esboce una gráfica de la función que describe el des-
plazamiento del sistema, suponiendo:
(a) que no hay amortiguamiento
(b) que existe amortiguamiento igual a 4 veces la velocidad del cuerpo.
7. Movimiento forzado con amortiguamiento: Un cuerpo, que tiene una ma-
sa de 2 [kg], se une a un resorte que cuelga del techo, haciendo que el resorte se
estire 20 [cm] hasta llegar en reposo al equilibrio. En el instante t = 0, el cuerpo
es desplazado 5 [cm] por debajo de la posición de equilibrio y es liberado. En ese
mismo instante se aplica una fuerza externa F (t) = 0, 3 cos(t) [N ] al sistema. Si la
constante de amortiguamiento para el sistema es 5 [N · s/m], encuentre la ecuación
de movimiento para el sistema. Determine el desplazamiento de la masa.
8. Resonancia sin amortiguamiento: Considere el sistema de masa-resorte mode-
lado por
10 X 00 (t) + k X(t) = f (t),
donde la la función de forzamiento es f (t) = cos(ωt) + 2 sen(ωt/2). Conociendo que
cos(ωt) oscila con una frecuencia entre 20 y 80 [Hz]. ¿Qué valores de k pueden llevar
al sistema a la resonancia?

13
9. Cicuito eléctrico RLC. Un circuito eléctrico RLC (Resistencia - Capacitor - In-
ductor) conectado en serie a un generador es gobernado por la ecuación diferencial
di
L + Ri + q/C = E (4)
dt
donde L representa a la constante de inductancia medida en Faradios [F ], R la
constante de resistencia medida en Ohms [Ω], C es la capacidad eléctrica del ca-
pacitor medida en Henrios [H], i la intensidad de corriente eléctrica en el circuito
medida en Amperes [A], q es la carga eléctrica medida en Coulombs [C], E es la
fuerza electromotriz del sistema medida en Volts [V ] y t es el tiempo medido en
[s]. Además, como la corriente eléctrica es el flujo de carga se tiene que i = dq
dt
. De
donde la ecuación (4) se puede escribir como
q
Lq̈ + Rq̇ + =E (5)
C
o bien, derivando con respecto al tiempo
d2 i di i dE
L 2
+R + = (6)
dt dt C dt
(a) Utilice la ecuación (5) o (6) para plantear una equivalencia con el sistema
mecánico masa – amortiguador – resorte.
(b) Se tiene un sistema eléctrico que consta de una fuerza electromotriz E(t) =
100 sin(ωt)[V ], una resistencia de R = 2[Ω], un inductor de L = 0,1[H] y un
condencador de C = 1/260[F ]. Si la corriente y la carga inicial del condensador
son cero, determinar la carga del condensador en el instante t > 0. ¿Para que
valor de ω el sistema manifiesta resonancia?

14
20/10/2020 MDT/JMS/CMG//cmg

15