Certamen 2 - Ecuaciones Diferenciales Ordinarias (2008-2)

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 6

UNIVERSIDAD DE CONCEPCION

FACULTAD DE CIENCIAS
FISICAS Y MATEMATICAS
DEPARTAMENTO DE INGENIERIA MATEMATICA

EVALUACION 2 (521218)

Problema 1. Resolver la siguiente EDO

x2 y 00 (x) + 2xy 0 (x) − 2y(x) = x−2 + sen(ln(x))

Solución Problema 1.

Notamos que la EDO Ly(x) = f (x) , con L = x2 D2 + 2xD − 2 y f (x) = x−2 + sen(ln(x)),
es del tipo Euler-Cauchy. Hacemos entonces el cambio de variable x = ez , con lo que nuestra
EDO se escribe, respecto a la variable independiente z, como una de coeficientes constantes
L̃ỹ(z) = f˜(z) , siendo

L̃ = D(D − 1) + 2D − 2 = D2 + D − 2 , f˜(z) = f (ez ) = e−2z + sen(t)


d d
Se recuerda que D = dx
yD= dz
, y que la solución general de la EDO transformada, es :

y(z) = yh (z) + yp (z)


Resolviendo L̃ỹ = 0.
La ecuación caracterı́stica asociada a la EDO es

λ2 + λ − 2 = 0 (λ + 2)(λ − 1)

de donde se tiene λ1 = −2, λ2 = 1.


En consecuencia, la solución de la EDO homogénea (en z! ) , es:

ỹh (z) = C1 e−2z + C2 ez ,

siendo C1 , C2 ∈ R constantes arbitrarias.


Hallando ỹp tal que L̃ỹp , = f˜.
Aplicando el método de los aniquiladores, se busca un aniquilador para f˜ de menor orden
posible, esto es,
L̃1 f˜ = (D + 2)(D2 + 1)f˜ = 0
.
Ası́ , la solución particular buscada es de la forma

ỹp = Aze−2z + Bcos(z) + Csen(z) + De−2z + Eez

1
que al reemplazar en la EDO se reduce a

ỹp (z) = Aze−2z + Bcos(z) + Csen(z)

con A, B y C reales a determinar. Derivando ỹp (z) y reemplazando en L̃ỹ = f˜, resulta

−3Ae−2z + (C − B)cos(z) − (3C + B)sen(z) = e−2z + sen(z)

de donde A = −1/3 , B = 1/2 , C = −1/2 ,


con lo cual se obtiene
1 1 1
ỹp (z) = − ze−2z + cos(z) − sen(z).
3 2 2
De esta manera, la solución general de L̃ỹ = f˜ es
1 1 1
ỹ(z) = C1 e−2z + C2 ez − ze−2z + cos(z) − sen(z).
3 2 2
Finalmente, la solución general de la EDO original, es
1 1n o
y(x) = C1 x−2 + C2 x − x−2 ln(x) + cos(ln(x)) − sen(ln(x)) .
3 2
siendo C1 , C2 , C3 ∈ R constantes arbitrarias.

2
Pregunta 2.

(i) Dos tanques A y B cada uno de 55 litros de capacidad, contienen respectivamente 12 y


15 lts de agua. Agua con sal a una concetración de 0,5 Kg/lt se bombea al tanque A a
razón de 2lt/min. La mezcla fluye del tanque A al tanque B a razón de 4 lt/min, y del
tanque B al tanque A a razón de 2 lt/min. Si del tanque B fluye al externo 1 lt/min de
mezcla, escriba (SIN RESOLVER) el sistema de ecuaciones diferenciales que determina
la cantidad de masa de sal hasta el instante del derrame. (No se pide resolver).

(ii) Considere el sistema masa-resorte dado por

m x00 (t) + k x(t) = h(t)

Si h(t) = 3cos(wt), m = 25 , ¿Para qué valores de k el sistema manifiesta resonancia si


cos(wt) tiene una frecuencia de 15 Hz ?

Solución Problema 2.

(i) Sean x(t) la masa de sal en el tanque A e y(t) la masa de sal en el tanque B.
Observe que el volumen en el tanque A es siempre constante en 12 lt; en cambio en el
tanque B la variación de volumen es ∆VB = 1, de donde VB (t) = t + 15. Por lo tanto el
derrame en el tanque B inicia a los 40 minutos. Ası́, el sistema se reduce a:

 0 y(t) x(t)
 x (t) = 1 + 2 t+15 − 4 12 ; t ∈ [0, 40]
y 0 (t) = 4 x(t)
12
y(t)
x(t) − 3 t+15 ; t ∈ [0, 40]

x(0) = 0, y(0) = 0

(ii) Tenemos que las soluciones del sistema homogéneo es


√ √
k k
xh (t) = Acos( ) + Bsen(( ).
5 5

Ası́ , el sistema manifiesta resonancia , si w = 5k donde 2π
w
= 15, es decir, si w = 30π.
Por lo tanto, existe resonancia si

k √
= 30π, esto es, si k = 150π
5

3
Problema 3. Resuelva el siguiente PVI

y 00 (t) + 4y 0 (t) + 6y(t) = H(t − 2)δ(t − 1) + δ(t − 2)


y 0 (0) = y(0) = 0
Solución:
Primero aplicamos T. de L. a H(t − 2)δ(t − 1) ; tenemos:

L[H(t − 2)δ(t − 1)](s) = e−2s L[f (t)]


donde f (t − 2) = δ(t − 1), esto es, f (t) = δ(t + 1). Ası́,

L[H(t − 2)δ(t − 1)](s) = e−2s L[δ(t + 1)](s) = e−2s es = e−s

Ahora, aplicamos Transformada de Laplace a la EDO dada, tenemos

s2 L[y(t)] + 4sL[y(t)] + 6L[y(t)] = e−s + e−2s

(s2 + 4s + 6)L[y(t)] = e−s + e−2s

−s
 1 
−2s
 1 
L[y(t)] = e +e .
(s + 2)2 + 2 (s + 2)2 + 2
Aplicando Transformada inversa, obtenemos:
    
−1 −s 1 −1 −2s
 1
y(t) = L e (t) + L e (t),
(s + 2)2 + 2 (s + 2)2 + 2

esto es, y(t) = H(t − 1)f1 (t − 1) + H(t − 2)f1 (t − 2)


donde √
2 −2t √ 
n o
1
f1 (t) = L−1 (t) = e sen 2 t .
(s + 2)2 + 2 2
Por lo tanto,
√ √
2 −2(t−1) √  2 −2(t−2) √ 
y(t) = H(t − 1) e sen 2 (t − 1) + H(t − 2) e sen 2 (t − 2) .
2 2

4
Problema 4. Usando valores propios determine la solución del siguiente PVI:
 0
 x (t) = x(t) ; x(0) = 0
0
y (t) = 2x(t) + y(t) ; y(0) = 0
 0
z (t) = 2x(t) + y(t) + z(t) ; z(0) = 0

SOLUCION.
De p(λ) := |A − λI| = 0, se obtiene que: λ1 = λ2 = λ3 = 1
Las soluciones del sistema se obtienen de
n o
λt 2 2 3 3
φ(t) = e I + t(A − λI) + 1/2 t (A − λI) + 1/6 t (A − λI) u

para u autovector genearalizado asociado a λ = 1.


De (A − λI)u1 = (0, 0, 0)t obtenemos el primer autovector generalizado , que resulta ser
u1 = (0, 0, 1) Ası́, la primera solución para el sistema ,es

X1 (t) = et (0, 0, 1).

Ahora de (A − λI)2 u2 = (0, 0, 0)t se obtiene u2 = (0, 1, 0) ( que es l.i. con u1 !); la correspon-
diente solución al sistema se obtiene de
n o
λt
φ(t) = e I + t(A − λI) u2 ,

realizando las respectivas operaciones se obtiene

X2 (t) = et (0, 1, t).

Finalmente de (A − λI)3 u3 = (0, 0, 0)t , se obtiene

u3 = (1, 0, 0) ( que es l.i. con u1 y con u2 ! )


y la tercera solución linealmente independientes con las otras dos, es
n o
λt 2 2
X3 (t) = e I + t(A − λI) + 1/2 t (A − λI) u3

que al realizar los cálculos, se obtiene

X3 (t) = (et , 2tet , (2t + t2 )et ).

Ası́, toda solución del sistema es dada por

Xh (t) = c1 X1 (t) + c2 (t)X2 (t) + c3 X3 (t)

5
Teniendo presente que Xh (0) = (x(0), y(0), z(0)) = (1, 1, 1), se obtiene
     
0 0 et
Xh (t) = c1  0  + c2  et  , +c3  2tet ,
et tet (2t + t2 )et

c1 = c2 = c3 = 1
de donde la única solución al PVI dado, es.

(x(t), y(t), z(t)) = (0, 0, et ) + (0, et , tet ) + (et , 2tet , (2t + t2 )et )

es decir,

   
x(t) et
 y(t)  =  et + 2tet 
t 2
z(t) e (1 + 3t + t )

JMS/LNB//jms.
26/11/2008.

También podría gustarte