UNIVERSIDAD DE CONCEPCION

FACULTAD DE CIENCIAS
FISICAS Y MATEMATICAS
DEPARTAMENTO DE INGENIERIA MATEMATICA

Evaluación 1 Ecuaciones Diferenciales Ordinarias (521218) 2018-1

PROBLEMA 1. Considere la forma diferencial

e2x−y (x + 1)y − y dx + x e2x−y + y + 1 dy = 0 .

(a) Muestre que la forma dada no es exacta.

(b) Compruebe que la forma diferencial dada admite el factor integrante µ(x, y) = ey−x .

(c) Resuelva la ecuación diferencial asociada a la forma diferencial dada.

Desarrollo:

(a) La forma diferencial no es exacta pues

∂  2x−y 
e (x + 1) y − y = e2x−y (x − y − xy) − 1 , y sin embargo, (1)
∂y
∂  2x−y 
xe + y + 1 = e2x−y (1 + 2x). (2)
∂x
(2 puntos)

(b) Multiplicando la forma dada por el factor ey−x , se obtiene la forma

   
ex (x + 1) y − yey−x dx + x ex + (y + 1)ey−x dy = 0, (3)

que resulta exacta puesto que

∂ h x i ∂ h x i
e (x + 1) y − yey−x = x e + (y + 1)ey−x .
∂y ∂x

(3 puntos)

(c) Como la forma auxiliar es exacta, existe F = F (x, y) de modo que

∂  
F (x, y) = ex (x + 1) y − yey−x (4)
∂x
∂  
F (x, y) = x ex + (y + 1)ey−x (5)
∂y

de modo que integrando en (4) respecto de x , sigue

F (x, y) = xy ex + y ey−x + C(y)

 donde C(y) es una constante que depende de y. (5 puntos)
Ahora, derivando parcialmente la F respecto de y , sigue que

∂
F (x, y) = xex + ey−x (1 + y) + C 0 (y).
∂y
Igualando con (5), se obtiene que C(y) es constante.
Finalmente, la solución de la forma diferencial dada es F (x, y) = constante, esto es,

xy ex + y ey−x = c con c constante en R.

(5 puntos)

PROBLEMA 2.

(2.1) Considere el conjunto de funciones S := {e−x cos(2x) , x sen(3x)}. Determine la EDO

lineal homogénea a coeficientes constantes reales, de menor orden posible, que tenga a
las funciones dadas en S como soluciones particulares.
Desarrollo: Un aniquilador de e−x cos(2x), de menor orden posible, es L1 = (D +
1)2 + 4. A su vez, L2 = (D2 + 9)2 aniquila a x sen(3x), y es del menor orden posible.
Finalmente, la EDO solicitada es Ly(x) = 0, siendo

L := L1 L2 = [(D+1)2 +4](D2 +9)2 = D6 +2D5 +23D4 +36D3 +171D2 +162D+405 .

(5 puntos)

(2.2) Resuelva el siguiente PVI:

 2 00
x y (x) − x y 0 (x) + 5 y(x) = 0 , x>0
y(1) = 2 , y 0 (1) = −1 .

Desarrollo: Reconociendo que la EDO dada es tipo Euler-Cauchy, realizamos el

cambio de variable: x = et ⇔ t = ln(x). Esto permite transformar el operador
diferencial lineal a coeficientes variables

L := x2 D2 − x D + 5

en
L̃ := D(D − 1) − D + 5 = D2 − 2 D + 5 ,
d d
siendo D := dx
, y D := dt
.
De esta manera, queda por resolver ahora la EDO lineal homogénea a coeficientes
constantes:
L̃ỹ(t) = 0 ⇔ (D2 − 2 D + 5)ỹ(t) = 0 .
Ecuación caracterı́stica:

λ2 − 2 λ + 5 = 0 ⇔ (λ − 1)2 + 22 = 0 ,
 de donde las raı́ces son λ1 = 1 + 2 i, λ2 = 1 − 2 i. Las funciones asociadas son

ỹ1 (t) = et cos(2t) , ỹ2 (t) = et sen(2t) ,

que en la variable original corresponde con

y1 (x) = x cos(2 ln(x)) , y2 (x) = x sen(2 ln(x)) .

La solución general de la Ecuación de Euler-Cauchy es

y(x) = C1 x cos(2 ln(x)) + C2 x sen(2 ln(x)) ,

siendo C1 y C2 constantes reales arbitrarias. (7 puntos)

Ahora determinamos los valores de C1 y C2 para descubrir aquella solucón que satisface
las condiciones iniciales dadas. Para esto, se requiere calcular la derivada de y(x), la
cual resulta ser
   
0
y (x) = cos(2 ln(x)) − 2 sen(2 ln(x)) C1 + sen(2 ln(x) + 2 cos(2 ln(x)) C2 .

Imponiendo las condiciones iniciales, se deduce que C1 = 2 y C2 = −3/2. Finalmente,

la solución del PVI planteado es
3
y(x) = 2 x cos(2 ln(x)) − x sen(2 ln(x)) .
2
(3 puntos)
PROBLEMA 3. Determine la solución general de la ecuación diferencial ordinaria
1
y 00 (t) + 4 y (t) = ∀t ∈ ]π/2, π[ .
sen (2 t)
Desarrollo: Tenemos que
y (t) = yh (t) + yp (t) ,
donde t ∈ ]0, π/2[, yh es la solución general de

yh00 (t) + 4 yh (t) = 0 (6)

y yp es una solución cualquiera de

1
yp00 (t) + 4 yp (t) = ∀t ∈ ]π/2, π[ . (7)
sen (2 t)

(1 punto)
Primeramente resolveremos (6). El polinomio caracterı́stico de (6) es

p (λ) = λ2 + 4,

que tiene dos raı́ces simples ±2i. Por lo tanto

yh (t) = A cos (2 t) + B sen (2 t) ,

con A, B ∈ R constantes arbitrarias.
(4 puntos)
A continuación buscaremos la solución particular de (7) con la forma

yp (t) = A (t) cos (2 t) + B (t) sen (2 t) ,

donde las derivadas de A (t) y B (t) satisfacen el sistema de ecuaciones

A0 (t) cos (2 t) + B 0 (t) sen (2 t) = 0



1 . (8)
 A0 (t) (−2 sen (2 t)) + B 0 (t) (2 cos (2 t)) =
sen (2 t)

Resolviendo (8) obtenemos

 
0
 0 sen (2 t) 
A (t) =  1 /∆ = −1/∆
sen(2 t)
2 cos (2 t)

y  
 cos (2 t) 0  1 cos (2 t)
B 0 (t) =  1  /∆ = ∆ sen (2 t)
−2 sen (2 t) sen(2 t)
con  
 cos (2 t) sen (2 t)  = 2 cos (2 t)2 + 2 sen (2 t)2 = 2.

∆ = 
−2 sen (2 t) 2 cos (2 t)
R
Integrando se llega a A (t) = (−1/2)dt = −t/2 y
Z
1 cos (2 t) 1 1
B (t) = dt = ln (|sen (2 t)|) = ln (− sen (2 t))
2 sen (2 t) 4 4

pues t ∈ ]π/2, π[. Note que las constantes de integración se han tomado igual a 0 sin pérdida
de generalidad ya que solo buscamos una solución de (7). Por lo tanto
t 1
yp (t) = − cos (2 t) + ln (− sen (2 t)) sen (2 t)
2 4
(9 puntos)
Finalmente
t 1
y (t) = A cos (2 t) + B sen (2 t) − cos (2 t) + ln (− sen (2 t)) sen (2 t)
2 4
donde A, B ∈ R son constantes arbitrarias. (1 punto)
PROBLEMA 4. Un cuerpo de masa 1 [kg] se acopla a un resorte de rigidez 25 [N/m]
que cuelga del techo. El cuerpo se mueve solamente en la dirección vertical y está sujeto
a las acciones de la fuerza del resorte en el régimen de validez de la ley de Hooke y a
una fuerza de amortiguamiento proporcional a su velocidad. No se le aplica ninguna fuerza
externa (vibraciones amortiguadas libres). La constante de amortiguamiento es 8 [N seg/m].
Consideramos el eje z en la dirección vertical apuntando hacia abajo.

(4.1) Escribir la ecuación diferencial para la posición z(t) del centro de masa del cuerpo según
el eje z, medida respecto de la posición de equilibrio (z = 0), en cualquier instante de
tiempo t. Luego, determine la solución general de dicha ecuación.
Nota: No se pide deducir la EDO a partir de la ley de Newton).

(4.2) Ahora, si el cuerpo está inicialmente a z0 = 2 [cm] debajo de su posición de equilibrio,

y se suelta al instante t0 = 0 con una velocidad inicial de 2 [cm/seg] apuntando hacia
arriba, determine la posición del centro de masa del cuerpo en todo instante.

(4.3) ¿En qué instante pasará el cuerpo por su posición de equilibrio por primera vez?

(4.4) ¿Pasará el cuerpo a 2 [cm] arriba del punto de equilibrio en algún instante? Justificar
su respuesta (considere ln(2) ' 0,69).

(4.5) Bosquejar la gráfica de la función z(t).

Desarrollo:
Comentario: se deberı́a especificar las unidades después de todos los resultados.
Observación: Cada item tiene un puntaje máximo de (3 puntos)

(a) La ecuación del movimiento es

z 00 (t) + 2λz 0 (t) + ω 2 z(t) = 0 , (9)

donde
r s
8 N seg m−1 k 25 N m−1
λ= = 4 seg−1 , ω= = = 5 seg−1
2 × 1 kg M 1 kg

(recuerda que 1 N = 1 kg m seg−2 ). Dado que ω > λ, estamos en el régimen de vibraciones

subamortiguadas con el cuasiperiodo
2π 2π 2π
T =√ =p = seg .
ω 2 − λ2 (25 − 16)seg−2 3

El polinomio caracterı́stico P (r) = r2 + 8r + 25 de la EDO lineal homogénea de segundo

orden (9) tiene raı́ces complejas r± = −4 ± 3 i. Luego la solución general de (9) es

z(t) = e−4t (C1 cos (3t) + C2 sen (3t)) = Ae−4t sen (3t + φ) , (10)

donde las constantes reales C1 , C2 , la amplitud A y la fase φ son arbitrarias, y serán

fijadas con las condiciones iniciales.
(b) Las condiciones iniciales son
z(0) = z0 = 2 cm , z0 (0) = v0 = −2 cm seg−1 . (11)
Note que es necesario introducir un signo menos antes de 2 cm seg−1 , pues esta especifi-
cado que la velocidad inicial apunta hacia arriba y el eje z apunta hacia abajo.
Primera forma: Determinamos las constantes C1 y C2 .
Usando (10), obtenemos z(0) = C1 = 2 cm. Luego, calculando la derivada de z(t)
z 0 (t) = e−4t (−4C1 cos (3t) − 3C1 sen (3t) − 4C2 sen (3t) + 3C2 cos (3t)) ,
deducimos z 0 (0) = −4C1 + 3C2 = −2 cm seg−1 de la segunda condición inicial. Por lo
tanto, C2 = (−2 + 4 × 2)/3 = 2 cm. La posición z(t) del centro de masa del cuerpo está
dada por
z(t) = 2e−4t (cos (3t) + sen (3t)) cm (t en seg).
Aplicando fórmulas trigonométricas, se puede transformar dicha solución en
√  π
z(t) = 2 2e−4t sen 3t + cm (t en seg). (12)
4
Segunda forma: Determinamos directamente las amplitud A y fase φ en función de la
condiciones iniciales.
Usando (10) y (11), tenemos
z(0) = A sen φ = 2 cm , z0 (0) = −4A sen φ + 3A cos φ = −2 cm seg−1 . (13)
Substituyendo la primera ecuación en la segunda, se obtiene
A cos φ = (−2 + 4 × 2)/3 = 2 cm .
Para hallar A, tomamos el cuadrado de esa expresión y adicionamos el cuadrado de la
primera ecuación en (13):
p √ √
A = A2 cos2 φ + A2 sen2 φ = 22 + 22 = 2 2 cm
Para hallar φ, dividimos estas dos mismas ecuaciones, llegando a
A sen φ 2 cm π
tan(φ) = = =1 ⇒ φ= + kπ .
A cos φ 2 cm 4
Dado que sen(φ) > 0 y cos(φ) > 0, los números enteros k son pares y podemos eligir
k = 0 (el resultado no cambia si k = 2, 4, ...). Luego, la solución del PVI formado por la
EDO (9) y las condiciones iniciales (11) está dada por (12).
(c) El cuerpo pasa por su posición de equilibrio z = 0 en los instantes τm tales que
sen(3τm + π/4) = 0 ⇔ 3τm + π/4 = mπ ⇔ τm = (4m − 1)π/12 ,
donde m es un número entero arbitrario. Si m ≤ 0 los tiempos τm son negativos, si m = 1
obtenemos un tiempo τ1 positivo. Ası́, el cuerpo pasa por su posición de equilibrio por
primera vez en el instante
τ1 = π/4 seg
(d) Se puede mostrar con el siguiente argumento que el cuerpo nunca no pasará a 2 [cm]
arriba del punto de equilibrio, esto es, por la posición z(t) = −2 cm. Comenzamos
determinando el primer instante t∗ tal que Ae−4t < z0 para cada t > t∗ . Dado que
√ ln 2
2 2e−4t < 2 ⇔ 2 < e8t ⇔ <t,
8
es claro que t∗ = ln 2/8. Notamos que t∗ < τ1 (de hecho ln 2 ' 0,693 < 2π). Luego, se
deduce de (12) que |z(t)| < 2 cm para cada t > t∗ y en particular para cada t > τ1 .
Además, claramente 0 ≤ z(t) ≤ z0 = 2 cm si 0 ≤ t ≤ τ1 . Entonces, el cuerpo núnca
pasará por la posición z = −2 cm, es decir, a 2 cm mas arriba que su punto de equilibrio.
Note que si la velocidad inicial v0 hubiera sido mucho mayor en valor absoluto, es-
to no hubiera sido cierto. De hecho, si |v0 | es muy grande la cuasi-amplitud A =
p
z02 + (v0 + λz0 )2 de las oscilaciones puede ser mucho más grande que la posición inicial
z0 , de modo que el tiempo t∗ serı́a mayor que τ2 = 7π/12. En este caso, por continuidad
de la función solución, z(t∗ ) = −2 cm en algún instante t∗ ∈ [τ1 , τ2 ]. Hay que señalar que
este instante t∗ no es fácil de determinar de manera explı́cita.

(e) El bosquejo de la gráfica de z(t) deberı́a mostrar claramente:

a) el decaimiento√exponencial de la amplitud de las oscilaciones (dado por las funciones

±Ae−λt = ±2 2e−4t cm);
b) el punto (0, 2) que corresponde a la posición inicial (aquı́ en centı́metros);
c) la tangente z 0 (0) = −2 cm/seg que apunta hacia el semi-eje z negativo;
d ) el instante τ1 = π/4 en el cual z(t) se anula por primera vez.

3
exp(-4*x)*2*sqrt(2)*sin(3*x+pi/4)
2*sqrt(2)*exp(-4*x)
-2*sqrt(2)*exp(-4*x)
2.5 0

2
z(t) , A exp(-lambda t) [cm]

1.5

0.5

-0.5
0 0.5 1 1.5 2
tiempo [seg]