Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

PEP2 - Invierno - Pauta

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

Universidad de Santiago de Chile

Facultad de Ciencia, Departamento Matemática y C.C.


Coordinación de Ecuaciones Diferenciales y Métodos Numéricos
Semestre Intensivo 1/2022

PEP No.2
NOMBRE: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . R.U.T. . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1. (2 ptos.) Un sistema mecánico se le aplica una fuerza externa que actúa sólo durante un intervalo
de tiempo t ∈ [0, π]. Si la ecuación que describe el movimiento es:
y ′′ (t) + y(t) = f (t), y(0) = 0, y ′ (0) = 0,
donde 
 cos(2t) 0 ≤ t < π,

f (t) =

 0, t ≥ π.

(a) Encuentre la solución del P.V.I. usando la Transformada de Laplace.


(b) Escriba la solución como una función definida por parte y verifique si la solución es continua.
Indicación: Simplifique su respuesta.
Solución
(a) La función en términos del escalón es: f (t) = cos (2t) − cos (2t)u(t − π). Aplicando transfor-
mada en ambos lados de la E.D.O.:
s e−πs s
s2 Y (s) + Y (s) = 2 − 2
s +4 s +4
s e−πs s
Y (s) = 2 −
(s + 4)(s2 + 1) (s2 + 4)(s2 + 1)
e−πs s
   
−1 s −1
y(t) = L −L
(s2 + 4)(s2 + 1) (s2 + 4)(s2 + 1)
    
−1 s s −1 −πs s s
y(t) = L − −L e −
3(s2 + 1) 3(s2 + 4) 3(s2 + 1) 3(s2 + 4)
 
cos(t) cos(2t) cos(t) cos(2t)
y(t) = − + + u(t − π)
3 3 3 3
(b) La función escrita por tramos es:
 cos(t) cos(2t)

 3 − 3 0 ≤ t < π,

y(t) = 2 cos(t)


 3 , t ≥ π.

Evaluando la continuidad en el extremo t = π:


cos(π) cos(2π) 2
y(π) = − =
3 3 3
2 cos(π) 2
y(π) = =
3 3
sin(π) 2 sin(2π)
y ′ (π) = − + =0
3 3
2 sin(π)
y ′ (π) = − =0
3

1
Las funciones y sus derivadas comparten el mismo valor en t = π, por lo que es continua.

2. (2 ptos.) Considere el sistema de ecuaciones diferenciales dado por


   
′ 3 −5 1
X (t) = 3 X(t) + et/2
4 −1 −1

(a) Encuentre la solución homogenea del sistema.


(b) Calcule una solución particular del sistema no homogéneo, usando variación de parámetros.

Solución:

(a)

3 − λ −5 2
0.75 −1 − λ = 0 =⇒ 4λ − 8λ + 3 = 0 =⇒ λ1 = 1.5, λ2 = 0.5

      10 
1.5 −5 a 0
λ1 = 1.5 =⇒ = =⇒ 0.75a − 2.5b = 0 ⇒ v = 3
0.75 −2.5 b 0 1
      
2.5 −5 a 0 2
λ1 = 0.5 =⇒ = =⇒ 0.75a − 1.5b = 0 ⇒ v =
0.75 −1.5 b 0 1
 10   
2 0.5t
=⇒ Xh = c1 3 e1.5t + c2 e ; c1 , c2 ∈ R
1 1

(b) Solución particular del sistema no homogéneo:


Z
Xp = ϕ (t) ϕ−1 (t) F (t) dt

Como
0.75e−1.5t −1.5e−1.5t
 10
e1.5t 2e0.5t
 0.5t
−2e0.5t
   
3 −1 1 e
ϕ (t) = =⇒ ϕ (t) = 4 2t =
e1.5t e0.5t 3e
−e1.5t 103 e
1.5t −0.75e−0.5t 2.5e−0.5t
t
!
0.75e−1.5t −1.5e−1.5t −2e−t
Z Z    
−1 e2
=⇒ ϕ (t) F (t) dt = t dt =
−0.75e−0.5t 2.5e−0.5t −e 2 −3.25t
 20 0.5t
− 6.5te0.5t − 20
     
−3e 3 −6.5 t
=⇒ Xp = 0.5t 0.5t = + t e2
−2e − 3.25te −2 −3.25

(c) Solución general


 10     20    
3 1.5t 2 0.5t −3 −6.5 t
Xg = c1 e + c2 e + + t e2 ; c1 , c2 ∈ R
1 1 −2 −3.25

3. (2 ptos.) Considere el P.V.F.

y ′′ − y = −x2 + 2; 0 ≤ x ≤ 3; y(0) = 0; y(3) = 9,

(a) Use el método de diferencias finitas centradas para aproximar y(1) e y(2).
(b) Resuelva el PVF de forma analı́tica y entregue el error cometido en (a).

2
Solución

(a) Para que el método de esas aproximaciones se debe utilizar un valor h = 1. La E.D.O.
discretizada es:
yi+1 − 2yi + yi−1
− yi = −x2i + 2
h2
Y las condiciones pueden ser escritas con los siguientes subindices:

x0 = 0; y0 = 0; x1 = 1; y1 ; x2 = 2; y2 ; x3 = 3; y3 = 9

Para los valores i = 1 e i = 2:


y2 − 2y1 + y0
− y1 = −x21 + 2
12
y3 − 2y2 + y1
− y2 = −x22 + 2
12
y2 − 3y1 = 1
y1 − 3y2 = −11

El sistema tiene como soluciones y1 = 1; y2 = 4.


(b) La solución de la E.D.O. puede ser calculada por diversas formas, esta es:

y(x) = C1 ex + C2 e−x + x2

Al aplicar las condiciones iniciales:

y(0) = C1 e0 + C2 e0 + 02 = C1 + C2 = 0
y(3) = C1 e3 + C2 e−3 + 9 = C1 e3 + C2 e−3 = 0

La solución del sistema anterior es: C1 = 0 C2 = 0


Por lo tanto, la solución del P.V.F. es: y(x) = x2 , que al evaluarlo en los puntos aproximados
da:
y(1) = 12 = 1
y(2) = 22 = 4
Dando un error de 0.

También podría gustarte