M Etodos Num Ericos: Gu Ia de Trabajos PR Acticos Nro. 7

Métodos Numéricos
Cursada 2do Cuatrimestre 2012

Guı́a de Trabajos Prácticos Nro. 7
Temas: Problemas de valores iniciales. Método de Euler, Método de Heun, Método de la Serie
de Taylor, Método de Runge-Kutta. Concepto de métodos predictores-correctores.
Problemas de valores de frontera. Sistemas de ecuaciones diferenciales. Ecuaciones diferenciales
de orden superior. Introducción a las ecuaciones diferenciales parciales.
dy
1. Sea la EDO dt = 2ty 2 con condicion inicial y(0) = 1,
a) Calcule manualmente y(0.4) usando el método de Euler, con pasos h=0.04, h=0.02
y h=0.01
b) Grafique las soluciones.
1
c) Sabiendo que la solución de la EDO es y(t) = 1−t2
, calcule el error absoluto de las
soluciones obtenidas en el inciso anterior.
2. Sea la EDO y ′ = (t − y)/2, con condicion unicial y(0) = 1,
a) Encuentre manualmente la solución y(t) en el intervalo [0 3] con tamaños de paso

h = 1, h = 1/2 y h = 1/4.
b) Genere una tabla con las soluciones para cada método y cada tamaño de paso, mas
la solución exacta y(t) = 3e−t/2 − 2 + t.
c) Grafique las soluciones.
d) Calcule el error absoluto y relativo para cada caso cuando t = 3.
3. Realice un programa MATLAB que implemente el método de Euler simple.
4. Realice un programa MATLAB que implemente el método de Euler mejorado.
5. Realice un programa MATLAB que implemente el método de Runge-Kutta de cuarto

orden.
6. Encuentre la solución y(t) para el ejercicio 1 usando los métodos previamente programados,
para los mismos tamaños de paso.
7. Encuentre la solución y(t) para el ejercicio 2 usando los métodos previamente programados,
para los mismos tamaños de paso.
8. Se le aplica una fuerza de 10 Newtons a un cuerpo de masa m = 5g originalmente en

reposo. Se quiere poder graficar la distancia recorrida x como una función del tiempo t.
Plantear el problema como una ecuación diferencial de segundo orden, y resolver por Euler
entre t = 0 y t = 5, con intervalos de 1 segundo.
9. Un proyectil de masa m=0.11Kg se lanza verticalmente hacia arriba con una velocidad
inicial v0 = 80m/seg. Dicha velocidad irá disminuyendo por efecto de la fuerza de gravedad
y de la resistencia ofrecida por el aire. La ecuación diferencial que rige la velocidad del
proyectil es:
dv
m∗ = −m ∗ g − k ∗ v 2
dt
Calcule por medio del método de Euler mejorado los diferentes valores de velocidad que to-
ma el proyectil hasta alcanzar su altura máxima. Considere g = 9.8m/seg2 y k=0.002Kg/m
Cátedra de Métodos Numéricos - Departamento de Matemática - Facultad de Ingenierı́a - U.N.M.D.P. 1

a) Utilice un paso h=0.1seg.

b) Compare con los resultados obtenidos con h=0.01seg.
10. Aplicar el método de Adams-Bashforth-Moulton al problema del ejercicio 2 de esta guia,

con h = 0.25. Usar Runge-Kutta para calcular los primeros valores y1 , y2 , y3 (y0 = y(0)).
Calcular la primera iteracion manualmente. Calcular el error para t = 3. Comparar los
resultados con el comando ode113 de Matlab.
11. Dado el siguiente sistema de ecuaciones:
dx
= x + 2y
dt con x(0) = 6
dy y(0) = 4
= 3x + 2y
dt
a) Usar los métodos de Euler y Runge-Kutta para resolver el sistema de ecuaciones: para
t = 0.20, con intervalo h = 0.02.
b) Graficar y comparar el resultado con la solucion exacta data por
x(t) = 4e4t + 2e−t
y(t) = 6e4t − 2e−t
c) Comparar los resultados con el comando ode45 de Matlab.
12. Un problema de simulación denominado cazador-presa es modelado por el siguiente sistema

de ecuaciones diferenciales.
dp
= k1 ∗ p − k2 ∗ c ∗ p
dt
dc
= k3 ∗ c ∗ p − k4 ∗ c
dt
Calcule como varı́an las poblaciones de cazadores y presas por medio de un método de
Runge-Kutta sabiendo que k1 = 0.4, k2 = 0.02, k3 = 0.001, k4 = 0.03. Adopte como
condiciones iniciales p(0)=30 y c(0)=3.
13. Dado el problema

x

u′′ − 1 − u=x
5
con condiciones de frontera u(1) = 2, u(4) = −1, y tamaño de intervalo h = 1.
a) Resolver usando el método del disparo.

b) Resolver usando diferencias finitas.
14. Dada la ecuación diferencial u′′ = u, con u(1) = 1.175201 y u(3) = 10.017875

a) Resolver por el método de las diferencias finitas, usando h = 1 y h = 0.5.

b) Comparar la solucion obtenida en todos los puntos con la solución analı́tica u(x) =
sinh(x)
c) En que punto x se genera el mayor error absoluto respecto de la solución analı́tica?
d) Como se comporta el error con respecto a la disminución por la mitad del intervalo
h (en el punto x = 2)?
15. Se desea conocer el flujo de calor en un alambre aislado, de 2cm de longitud, cuyo modelo
responde a
∂T ∂2T
=K 2
∂t ∂x
donde K = 0.1515 es el coeficiente de conductividad. La temperatura inicial está dada por
T (x, t0 ) = 100x si x < 1 (1)

T (x, t0 ) = 200 − 100x si x ≥ 1 (2)
y las condiciones de contorno estan dadas por
T (0, t) = 0 (3)
T (2, t) = 0 (4)
Usar k = ∆t = 0.2 y h = ∆x = 0.25
a) Calcular las temperaturas manualmente para Ti,2 (t=0.2 seg) usando el método
explı́cito
b) Calcular las temperaturas para todo t hasta 10 seg, usando el método explı́cito
c) Calcular las temperaturas para todo t hasta 10 seg, usando el método implı́cito
d) Repetir los dos cálculos para k = ∆t = 0.25
e) Graficar y comparar los resultados (usar la función mesh or surf )
16. Analizar la vibración de una cuerda atada en los extremos usando la ecuación de onda
∂2y 2
2∂ y
= c
∂t2 ∂x2
c=2
0 < x < 1, 0 < t < 0.5
y(0, t) = 0, y(1, t) = 0 para 0 ≤ t ≤ 0.5
Posición inicial: y(0, t) = f (x) = sin(πx) + sin(2πx), para 0 < x < 1
Velocidad inicial: y ′ (0, t) = g(x) = 0, para 0 < x < 1

h = ∆x = 0.1
a) Analizar la vibración hasta t = 0.5 usando k = ∆t = 0.05

b) Analizar la vibración hasta t = 0.5 usando k = ∆t = 0.04
c) Graficar y comparar los resultados
17. Se desea conocer la distribución de temperaturas en estado de equilibrio en el interior de

una placa rectangular de acero, cuyo modelo responde a la siguiente ecuación diferencial.
∂2T ∂2T
2
+ =0
∂x ∂y 2
Halle los valores de temperatura en los nodos, por medio de diferencias finitas.

M Etodos Num Ericos: Gu Ia de Trabajos PR Acticos Nro. 7

Cargado por

Copyright:

Formatos disponibles

M Etodos Num Ericos: Gu Ia de Trabajos PR Acticos Nro. 7

Cargado por

Información del documento

Descripción original:

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

M Etodos Num Ericos: Gu Ia de Trabajos PR Acticos Nro. 7

Cargado por

Copyright:

Formatos disponibles

Métodos Numéricos

Cursada 2do Cuatrimestre 2012

2. Sea la EDO y ′ = (t − y)/2, con condicion unicial y(0) = 1,

a) Encuentre manualmente la solución y(t) en el intervalo [0 3] con tamaños de paso

3. Realice un programa MATLAB que implemente el método de Euler simple.

4. Realice un programa MATLAB que implemente el método de Euler mejorado.

5. Realice un programa MATLAB que implemente el método de Runge-Kutta de cuarto

8. Se le aplica una fuerza de 10 Newtons a un cuerpo de masa m = 5g originalmente en

Cátedra de Métodos Numéricos - Departamento de Matemática - Facultad de Ingenierı́a - U.N.M.D.P. 1

a) Utilice un paso h=0.1seg.

10. Aplicar el método de Adams-Bashforth-Moulton al problema del ejercicio 2 de esta guia,

11. Dado el siguiente sistema de ecuaciones:

x(t) = 4e4t + 2e−t

y(t) = 6e4t − 2e−t

c) Comparar los resultados con el comando ode45 de Matlab.

12. Un problema de simulación denominado cazador-presa es modelado por el siguiente sistema

13. Dado el problema

con condiciones de frontera u(1) = 2, u(4) = −1, y tamaño de intervalo h = 1.

a) Resolver usando el método del disparo.

Cátedra de Métodos Numéricos - Departamento de Matemática - Facultad de Ingenierı́a - U.N.M.D.P. 2

a) Resolver por el método de las diferencias finitas, usando h = 1 y h = 0.5.

donde K = 0.1515 es el coeficiente de conductividad. La temperatura inicial está dada por

T (x, t0 ) = 100x si x < 1 (1)

y las condiciones de contorno estan dadas por

Usar k = ∆t = 0.2 y h = ∆x = 0.25

Cátedra de Métodos Numéricos - Departamento de Matemática - Facultad de Ingenierı́a - U.N.M.D.P. 3

a) Analizar la vibración hasta t = 0.5 usando k = ∆t = 0.05

17. Se desea conocer la distribución de temperaturas en estado de equilibrio en el interior de

Cátedra de Métodos Numéricos - Departamento de Matemática - Facultad de Ingenierı́a - U.N.M.D.P. 4

También podría gustarte