LABORATORIO

PRIMER LABORATORIO DE MÉTODOS NUMÉRICOS II
- ALUMNOS:
Medina Sanchez, Brian (20181040)
Mendieta Pacotaype, Patricia (20181041)
Mendoza Yapias, Martha (20181042)
Elegir el nú mero de matrícula de uno de los integrantes del grupo e identificar a,b,c y d:
2 0 1 8 1 0 4 2
2 0 ? ? a b c d
P= (b+c+d)/3
Donde: a=1, b=0, c=4, d=2 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
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

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

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) = 60 + 3t
 V2(t) = V1(0) + t (3-2-4)

V2(t) = 60 - 3t
Hallamos las concentraciones:
 C1 = X1(t) / V1(t) = X1(t) / 60+3t

 C2 = X2(t) / 60 - 3t
Para el tanque A:
 dX1/dt = (4 gal/min) (0 lb) +(2 gal/min) (X2(t) / 60 - 3t) - (3 gal/min) (X1(t) / 60 +

3t)
dX1/dt = (2*X2(t)/60 - 3t) - (3*X1(t) / 60+3t)
Para el tanque B
 dX2/dt = (3 gal/min) (X1(t) / 60 + 3t) -(2 gal/min) (X2(t) / 60 - 3t) - (4 gal/min)

(X1(t) / 60 + 3t)
dX2/dt = (3*X1(t)/60 + 3t) - (6*X2(t) / 60-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)/60 - 3t) - (3*X1(t) / 60+3t)
dX2/dt = (3*X1(t)/60 + 3t) - (6*X2(t) / 60-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/(60-3*t))*x(2) - ((3/(60+3*t))*x(1)); (3/(60+3*t))*x(1)-(6/(60-

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
{
2 2
=1+ y + sen(t + y ) , 4 ≤t ≤ 9
dt
y ( 4 )=2
a) Runge Kutta clásico
b) Aplicando el método predictor corrector que se obtiene a partir de los siguientes
métodos multipasos de Adams-Moulton y Adams-Bashfort respectivamente
h
y p +1= y p + 9 f ( t p+1 , y p +1 ) +19 f ( t p , y p ) −5 f ( t p−1 , y p−1 ) + f ( t p−2 , y p−2) )
24 (
h
y p +1= y p + (55 f ( t p , y p ) −59 f ( t p−1 , y p−1 ) +37 f ( t p−2 , y p−2 ) −9 f ( t p−3 , y p−3 ) )
24
c) Compare los resultados obtenidos en a) y b) Graficando i) t vs y
Solución
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

5.
6. Aproxime la solución del problema
∂2 u ∂2 u
=3 , 0< x< 2 ,0< t<0.4
∂ x2 ∂ t2
u ( 0 , t ) =0 ,u ( 2 , t )=0 , 0 ≤t ≤0.4
❑
∂u
u ( x , 0 ) =x ( 2−x ) , ∫ ¿ x ,0 ≤ x ≤ 2
∂ t t=0
Considerando el paso espacial de 0.5 y el paso temporal de 0.2.

Solución
∂2 U 2
2∂ U ∂2 U 1 ∂2 U
=c → = , 0< x <2 , 0<t <0.4
∂t 2 ∂ x2 ∂ t 2 √32 ∂ x 2
U ( 0 , t ) =0 , para 0 ≤ t ≤ 0.4
U ( 2 ,t )=0 , para 0 ≤t ≤0.4
U ( x , 0 ) =f ( x )=x ( 2−x ) , para 0 ≤ x ≤2
U ( 0.5,0 )=U ( 1.5,0 ) =0.5 × ( 2−0.5 )=0.75 ; U ( 1,0 )=1 × ( 2−1 )=1
∂u
( x , 0 )=g ( x ) =x , para 0 ≤ x ≤2
∂t
∆t 0.2
h=∆ x=0.5 , k =∆ t=0.2 , r=c = =0.23
∆ x 0.5 × √3
r =0.23 ≤1 , implicaque hay estabilidad numérica
Reemplazando
U i , j+1 =( 2−2 r 2 ) ×U i , j +r 2 × ( U i +1, j +U i−1 , j )−U i , j −1
U i , j+1 =1.8942×U i , j +0.0529 × ( U i+1 , j +U i−1 , j )−U i , j−1
Para calcular U i , j+1 se requiere los valores de U i , j ,U i+1 , j , U i−1 , j ,U i , j−1
A partir de los valores iniciales se obtiene

2
2 r
U i , 1=( 1−r ) × f i +k × g i+ × ( f i +1+ f i−1 ) para i=2,3 ,… . , n−1
2
Donde f i=f ( x )=x i ( 2−xi ) , gi=g( x ) =xi , r=0.23

i i
U i , 1=0.9471× f ( x ) +0.2 × g( x ) +0.02645 × ( f ( x ) + f ( x

i i i+ 1 i−1 ) )
 Para i=1
U 1,1 =0.9471× f ( x ) +0.2 × g( x )+ 0.02645× ( f ( x )+ f ( x
1 1 2 0 ) )
U 1,1 =0.9471×(0.5 × (2−0.5 ) )+ 0.2× 0.5+0.02645 × ( 1 × ( 2−1 )+ 0 ×(2−0) )
U 1,1 =0.836775
 Para i=2
U 2,1 =0.9471× f ( x ) +0.2 × g ( x )+ 0.02645× ( f ( x ) +f ( x
2 2 3 1 ) )
U 2,1 =0.9471×(1 × ( 2−1 ))+0.2 ×1+0.02645 × ( 1.5 × (2−1.5 ) +0.5 ×(2−0.5) )
U 2,1 =1.186775
 Para i=3
U 3,1 =0.9471× f ( x )+ 0.2× g ( x )+ 0.02645× ( f ( x ) + f ( x ) )

3 3 4 3
U 3,1 =0.9471×(1.5 × ( 2−1.5 ))+0.2 ×1.5+ 0.02645× ( 2× ( 2−2 ) +1 ×(2−1) )
U 3,1 =1.036775
Con esto se completa la segunda fila:

Aplicando:
U i , j+1 =1.8942×U i , j +0.0529 × ( U i+1 , j +U i−1 , j )−U i , j−1
 Para j=1
U i , 2=1.8942× U i ,1 +0.0529 × ( U i +1,1 +U i−1,1 )−U i ,0
- Para i=1
U 1,2 =1.8942×U 1,1+ 0.0529× ( U 2,1 +U 0,1 ) −U 1,0
U 1,2 =1.8942× 0.836775+ 0.0529× ( 1.186775+0 ) −0.75
U 1,2 =0.8977996025
- Para i=2
U 2,2 =1.8942×U 2,1+ 0.0529 × ( U 3,1 +U 1,1 ) −U 2,0
U 2,2 =1.8942×1.186775+ 0.0529× ( 1.036775+ 0.836775 )−1
U 2,2 =1.3471
- Para i=3
U 3,2 =1.8942×U 3,1+ 0.0529× ( U 4,1 +U 2,1 )−U 3,0
U 3,2 =1.8942×1.036775+ 0.0529× ( 0+1.186775 )−0.75
U 3,2 =1.276639603
Con los datos podemos completar la tercera fila:
U ( 0.5,0 .2 )=U ( x 1 ,t 1)
≈U 1,1=0.836775
U ( 1,0.2 )=U ( x ,t ) ≈U 2,1=1.186775

2 1
U ( 1.5,0.2 )=U ( x ,t ) ≈ U 3,1=1.036775

3 1
U ( 0.5,0 .4 ) =U ( x , t ) ≈ U 1,2 =0.8977996025

1 2
U ( 1,0.4 ) =U ( x ,t ) ≈ U 2,2 =1.3471

2 2
U ( 1.5,0.4 )=U ( x , t ) ≈ U 3,2=1.276639603

3 2

LABORATORIO

Cargado por

Copyright:

Formatos disponibles

LABORATORIO

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

LABORATORIO

Cargado por

Copyright:

Formatos disponibles

PRIMER LABORATORIO DE MÉTODOS NUMÉRICOS II

Donde f ( t , y )=cos ⁡(7 y−2−3 t 2+ 3)

Supongamos que f ( t , y )=cos ⁡(7 y−2−3 t 2+ 3) es continua sobre

b) Resuelva el siguiente p.v.i. para δ =0 , ε =± 10−q , q=2,3,4, aplicando el método de

c) Resuelva el siguiente p.v.i. para ε =0 , δ=± 10−q , q=2,3,4, aplicando el método de

Se observa en los problemas perturbados cuando y0=0 da como resultado NaN

2. Resuelva el siguiente p.v.i usando longitud de paso h=0.05

b) Runge Kutta clásico

Compare los resultados obtenidos en a) y b).

d) ¿Cuál de estos métodos recomendaría?

 V1(t) = V1(0) + t (4+2-3)

 V2(t) = V1(0) + t (3-2-4)

Hallamos las concentraciones:

 C1 = X1(t) / V1(t) = X1(t) / 60+3t

 dX1/dt = (4 gal/min) (0 lb) +(2 gal/min) (X2(t) / 60 - 3t) - (3 gal/min) (X1(t) / 60 +

 dX2/dt = (3 gal/min) (X1(t) / 60 + 3t) -(2 gal/min) (X2(t) / 60 - 3t) - (4 gal/min)

dX1/dt = (2*X2(t)/60 - 3t) - (3*X1(t) / 60+3t)

dX2/dt = (3*X1(t)/60 + 3t) - (6*X2(t) / 60-3t)

 b. Para encontrar la cantidad de sal en el tanque B a los 5 minutos, aplicamos Runge-

 >> f= @(t,x)[(2/(60-3*t))*x(2) - ((3/(60+3*t))*x(1)); (3/(60+3*t))*x(1)-(6/(60-

 Los resultados se muestran a continuación:

k tk (x1)k ~ x1(tk) (x2)k ~ x2(tk)

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)

a) Runge Kutta clásico

c) Compare los resultados obtenidos en a) y b) Graficando i) t vs y

Considerando el paso espacial de 0.5 y el paso temporal de 0.2.

U i , j+1 =( 2−2 r 2 ) ×U i , j +r 2 × ( U i +1, j +U i−1 , j )−U i , j −1

U i , j+1 =1.8942×U i , j +0.0529 × ( U i+1 , j +U i−1 , j )−U i , j−1

Para calcular U i , j+1 se requiere los valores de U i , j ,U i+1 , j , U i−1 , j ,U i , j−1

A partir de los valores iniciales se obtiene

Donde f i=f ( x )=x i ( 2−xi ) , gi=g( x ) =xi , r=0.23

U i , 1=0.9471× f ( x ) +0.2 × g( x ) +0.02645 × ( f ( x ) + f ( x

U 3,1 =0.9471× f ( x )+ 0.2× g ( x )+ 0.02645× ( f ( x ) + f ( x ) )

U 3,1 =0.9471×(1.5 × ( 2−1.5 ))+0.2 ×1.5+ 0.02645× ( 2× ( 2−2 ) +1 ×(2−1) )

Con esto se completa la segunda fila:

U i , j+1 =1.8942×U i , j +0.0529 × ( U i+1 , j +U i−1 , j )−U i , j−1

U 2,2 =1.8942×U 2,1+ 0.0529 × ( U 3,1 +U 1,1 ) −U 2,0

U 2,2 =1.8942×1.186775+ 0.0529× ( 1.036775+ 0.836775 )−1

U 3,2 =1.8942×U 3,1+ 0.0529× ( U 4,1 +U 2,1 )−U 3,0

U 3,2 =1.8942×1.036775+ 0.0529× ( 0+1.186775 )−0.75

Con los datos podemos completar la tercera fila:

U ( 1,0.2 )=U ( x ,t ) ≈U 2,1=1.186775

U ( 1.5,0.2 )=U ( x ,t ) ≈ U 3,1=1.036775

U ( 0.5,0 .4 ) =U ( x , t ) ≈ U 1,2 =0.8977996025

U ( 1,0.4 ) =U ( x ,t ) ≈ U 2,2 =1.3471

U ( 1.5,0.4 )=U ( x , t ) ≈ U 3,2=1.276639603

También podría gustarte

dX1/dt = (2X2(t)/60 - 3t) - (3X1(t) / 60+3t)

dX2/dt = (3X1(t)/60 + 3t) - (6X2(t) / 60-3t)

 >> f= @(t,x)[(2/(60-3t))x(2) - ((3/(60+3t))x(1)); (3/(60+3t))x(1)-(6/(60-