MA2601-02 Ecuaciones Diferenciales Ordinarias, Examen
5 de Diciembre, 2014
Profesor : Raúl Manasevich
Auxiliar : Sebastian Contreras
P1.
Resuelva las siguientes ecuaciones diferenciales
(a) (3 ptos.)
ex + yexy + (ey + xexy )y 0 = 0.
(b) (3 ptos.)
p
y0 = (x + y).
P2.
(a) (3 ptos.) Considere la ecuación de Bessel
x2 y 00 + xy 0 + (x2 − 1/4)y = 0.
Demuestre que la ecuación tiene una solución de la forma y1 (x) = xα cos(x),
para un cierto α ∈ R y x > 0. De aquı́ encuentre una segunda solución de
esta ecuación para x > 0.
(b) (3 ptos.) Considere la ecuación
Ly = 5et + 10e−t ,
con Ly = y (v) + a4 y (iv) + ... + a0 y, donde a0 , ... , a4 son constantes. Se sabe
que el polinomio caracterı́stico tiene raı́ces λ = 1 y λ = −1, con mult.
algebraicas 3 y 2 respectivamente. Encontrar forma de la sol. particular.
P3.
(a) (3 ptos.) Sea f : [0, ∞) 7→ R continua a trozos y de orden exponencial y tal
f (t)
que lı́m existe y es finito. Si como siempre F (s) = L(f (t)) denota la
t→0+ t
T de Laplace de f se pide demostrar que
f (t) Z ∞
L = F (σ)dσ.
t s
1
2
R∞R∞
Suponga que puede intercambiar integrales de la forma a b
, si es nece-
sario.
(b) (3 ptos.) Usando los resultados de (a), encuentre
1 − cos(2t)
L .
t
P4. Encuentre la solución del problema
x0 = Ax,
en la forma x(t) = etA x(0), si la matriz A es dada por
3 −4 1 0
4 3 0 1
.
0 0 3 −4
0 0 4 3
Tiempo: 3 horas.
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile
MA2601-2. Ecuaciones Diferenciales Ordinarias 2014.
Profesor: Raúl Manasevich
Auxiliar: Sebastián Contreras
Pauta Examen
P1] (a) Tenemos la ecuación
ex + yexy + (ey + xexy ) y 0 = 0.
| {z } | {z }
M (x,y) N (x,y)
Verificamos que es exacta...
∂M
= exy + xyexy (1)
∂y
∂N
= exy + xyexy (2)
∂x
Z
∂Φ
Como (1) = (2) concluı́mos que sı́, de manera que M = ⇒ M dx = Φ(x, y).
∂x
Z
(ex + yexy )dx = ex + exy + C(y)
∂Φ
Notando además que N = tenemos...
∂y
ey + = ∂Φ = ye
xy
xe + C 0 (y)
xy
∂y
de manera que concluı́mos que C(y) = ey y de esta forma la solución está dada por
ex + ey + exy = C̃
(b)
√ √
y0 = x+y CV: z = x+y
0
2
z = x+y /()
2zz 0 = 1+z
z + +1 − 1 1
⇒ =
1+z 2
Z
1 1
⇔ 1− = / ()
z+1 2
x
z − ln |1 + z| = +C
2
√ √ x
x + y − ln |1 + x + y| = +C
2
1
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile
P3] a)
Z ∞ Z ∞ Z ∞
F (σ)dσ = f (t)e−σt dtdσ
s s 0
Z ∞ Z ∞
= f (t)e−σt dtdσ
0 s
s
∞
e−σt
Z
=
t
0
∞
∞
e−st
Z
=
0 t
f (t)
= L (s)
t
b)
1 − cos(2t) Z ∞
L = L (1 − cos 2t)dσ
t
Zs ∞
1 1 2σ
= −
s σ 2 σ2 + 4
n→∞ n→∞
1 2
= ln σ − ln σ + 4
2
s s
n s
= ln √ − ln √
2
2
n + 4 n→∞
s +4
Ordenando términos obtenemos...
r ! r
4 4 1 1 4 1 4
lı́m ln 1 + 2 + ln 1 + 2 = ln 1 + ln 1 + 2 = ln 1 + 2
n→∞ n s 2 2 s 2 s
2
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile
P4] Es buena idea escribir la matriz A por bloques;
3 −4 1 0
= B I B 0 I
4 3 0 1 0
A= = +
0 0 3 −4 0 B 0 B 0 0
0 0 4 3 | {z } | {z }
C D
3 −4
Donde 0 y I denotan la matriz cero y la identidad de 2 × 2, y B = , de modo que A = C + D. Para
4 3
usar la propiedad e(C+D)t = eCt eDt debemos demostrar que C y D conmutan...
B 0 0 I 0 B
BC = ∗ =
0 B 0 0 0 0
0 I B 0 0 B
CB = ∗ =
0 0 0 B 0 0
Verificado, procedemos entonces...
eBt
Dt 0
e =
0 eBt
de modo que primero debemos calcular eBt .
3−λ −4
p(λ) = det = (3 − λ)2 + 42 = (3 − λ)2 − (4i)2 = (3 − 4i − λ)(3 + 4i − λ)
4 3−λ
de modo que λ = 3 ± 4i son los valores propios. Calculamos el subespacio propio asociado a λ = 3 − 4i
4i −4
ker B − (3 − 4i)I = ker
4 4i
Notando que las filas son ponderación una de la otra obtenemos que el vector propio es
v = 1 0 +i 0 1
las soluciones están dadas por
cos 4t 3t
x1 (t) = e
− sin 4t
sin 4t
x2 (t) = e3t
cos 4t
De modo que una matriz fundamental (que en este caso coincide con la exponencial) está dada por
cos 4t sin 4t
e3t
− sin 4t cos 4t
eDt es nilpotente y conviene calcularla por serie de potencias...
∞
X Di t2 2
eDt = = I + tD + D+ · · ·
i=0
i! 2!
. Luego Bt
eBt teBt
e 0 I tI
eAt = eCt eDt = ∗ =
0 eBt 0 I 0 eBt
[email protected]
