LABORATORIO DE MÉTODOS NUMÉRICOS II
Mendieta Pacotaype, Patricia (20181041)
2 0 1 8 1 0 4 1
2 0 ? ? a b c d
P= (b+c+d)/3
Donde: a=1, b=0, c=4, d=1 y p=2
1. Dado el problema de valor inicial (p.v.i)
dy
{ dt
=f ( t , y ) +δ , 1 ≤t ≤ 2
y ( 1 )=0+ ε
Donde f ( t , y )=cos (7 y−2−3 t 2+ 3)
a) Averigüe si el problema de valor inicial está bien planteado
b) Resuelva el siguiente p.v.i. para δ =0 , ε =± 10−q , q=2,3,4, aplicando el método de Runge
Kutta Clásico y displaye su resultado
c) Resuelva el siguiente p.v.i. para ε =0 , δ=± 10−q , q=2,3,4, aplicando el método de Runge
Kutta Clásico y displaye su resultado
Solución
a) Averigüe si el problema de valor inicial está bien planteado
dy
{ dt
=f (t , y ) ,1 ≤ t ≤ 2
y ( 1 )=0
Supongamos que f ( t , y )=cos (7 y−2−3 t 2+ 3) es continua sobre D= { ( t , y ) :1 ≤ t ≤2 ,−∞< y <+∞ }
. Si f verifica la condición de Lipschitz con respecto a su variable y en D , entonces el p.v.i tiene una
única solución para y= y (t) para t ∈[ 1,2] y es un problema bien planteado.
∂ f (t , y )
| ∂y | =¿
≤|14 y−3|≤ L
La condición de Lipschitz se cumple si L>0, entonces cumpliría para todo y ≠ 0, entonces es un
problema bien planteado.
NOTA: Para las preguntas b y c, cambiamos a y0=1, debido a la presencia de NaN en los valores de
y cuando y0=0 en el problema perturbado.
b) Resuelva el siguiente p.v.i. para δ =0 , ε =± 10−q , q=2,3,4, aplicando el método de Runge
Kutta Clásico y displaye su resultado
Presentaremos la resolución primero con y0=0 y luego con y0=1.
1. Con y0=0
- Editor:
- Valores de t e y:
- Gráfico:
2. Con y0=1
- Editor:
- Valores de t e y:
- Gráfica:
c) Resuelva el siguiente p.v.i. para ε =0 , δ=± 10−q , q=2,3,4, aplicando el método de Runge
Kutta Clásico y displaye su resultado
Presentaremos la resolución primero con y0=0 y luego con y0=1.
1. Con y0=0
- Editor:
- Valores de t e y:
- Gráfica:
2. Con y0=1
- Editor:
- Valores de t e y:
- Gráfica:
Se observa en los problemas perturbados cuando y0=0 da como resultado NaN debido a que no
producen un resultado numérico bien definido. Pero cuando utilizamos y0=1 en ambos
problemas perturbados (b y c) nos dan resultados similares a la función no perturbada es por
esto que respalda la conclusión de la pregunta (a), el problema original está bien planteado ya
que las soluciones de los problemas perturbados cuando y ≠ 0 , es decir y0=1 se aproxima con
exactitud a la solución del problema original.
2. Resuelva el siguiente p.v.i usando longitud de paso h=0.05
d y1
{
=2 y 2 −2 y 1−5t
dt
d y2
=3 y 1−4 y 2+ 4 t ,2 ≤ t ≤ 3
dt
y 1 ( 2 )=2
y2 ( 2 )=1
Aplicando los métodos de
a) Taylor de segundo orden
b) Runge Kutta clásico
c) Compare los resultados obtenidos en a) y b). Graficando: i) t vs y 1, ii) t vs y 2, iii) y 1 vs y 2
d) ¿Cuál de estos métodos recomendaría?
Solución
a) Taylor de segundo orden
- Editor:
- Comandos:
- Valores de t e y:
- Gráfico:
b) Runge Kutta clásico
- Editor:
- Comandos:
- Valores de t e y: - Gráfica:
c)
Compare los resultados obtenidos en a) y b).
Graficando:
i) t vs y 1
ii) t vs y 2
iii) y 1 vs y 2
d) ¿Cuál de estos métodos recomendaría?
3.
Hallamos V1 y V2 en el momento t:
V1(t) = V1(0) + t (4+2-3)
V1(t) = 70 + 3t
V2(t) = V1(0) + t (3-2-4)
V2(t) = 65 - 3t
Hallamos las concentraciones:
C1 = X1(t) / V1(t) = X1(t) / 70+3t
C2 = X2(t) / 65 - 3t
Para el tanque A:
dX1/dt = (4 gal/min) (0 lb) +(2 gal/min) (X2(t) / 65 - 3t) - (3 gal/min) (X1(t) / 70 + 3t)
dX1/dt = (2*X2(t)/65 - 3t) - (3*X1(t) / 70+3t)
Para el tanque B
dX2/dt = (3 gal/min) (X1(t) / 70 + 3t) -(2 gal/min) (X2(t) / 65 - 3t) - (4 gal/min) (X1(t) / 70 + 3t)
dX2/dt = (3*X1(t)/70 + 3t) - (6*X2(t) / 65-3t)
Ahora se muestra el modelo matemático que describe los números x1 (t) y x2 (t) de libras de
sal en los tanques A y B, respectivamente, en el momento t :
dX1/dt = (2*X2(t)/65 - 3t) - (3*X1(t) / 70+3t)
dX2/dt = (3*X1(t)/65 + 3t) - (6*X2(t) / 70-3t)
X1 (0) = 30
X2 (0) = 0
b. Para encontrar la cantidad de sal en el tanque B a los 5 minutos, aplicamos Runge-Kutta en
octave:
>> f= @(t,x)[(2/(65-3*t))*x(2) - ((3/(70+3*t))*x(1)); (3/(70+3*t))*x(1)-(6/(65-3*t))*x(2)];
>> [t, x] = runge_kutta_C(f, Tif= [0,5], y0= [30;0], h0=0.1)
Los resultados se muestran a continuación:
k tk (x1)k ~ x1(tk) (x2)k ~ x2(tk)
0 0.0000000 30.000000 0.0000000
1 0.100000000000000 29.8510 0.1485
2 0.200000000000000 29.7040 0.2941
... ... ... ...
48 4.799999999999999 24.6551 4.3556
49 4.899999999999999 24.5753 4.3969
50 5.000000000000000 24.4965 4.4366
Entonces, la cantidad de sal en el tanque B a los 5 minutos es 4.4366 lb
4. Resuelva el problema de valor inicial (p.v.i)
dy
{ dt
=1+ y + sen(t 2 + y 2) , 4 ≤t ≤ 9
y ( 4 )=2
a) Runge Kutta clásico
b) Muestre los resultados en una tabla y grafica
a) Runge Kutta clásico
- Editor:
- Comandos:
- Valores de t e y:
-
- Gráfica:
b) Método predictor corrector
- Editor:
- Comandos:
- Valores de t e y:
-
- Gráfica:
c) Compare los resultados obtenidos en a) y b) Graficando i) t vs y