Taller 6 - Metodos Numericos

UNIVERSIDAD AUSTRAL DE CHILE Bachillerato en Ciencias de la Ingeniería
FACULTAD DE CIENCIAS Semestre 4

Instituto de Matemáticas
Taller 6
Bain 053 — Método Numéricos
1 Derivadas y ecuaciones diferenciales

1. Obtenga las siguientes fórmulas de aproximación de derivadas
−f(x + 2h) + 4f(x + h) − 3f(x)

(a) f (x) = + O(h2 )
2h
f(x + 2h) − 2f(x + h) + f(x)
(b) f (x) = + O(h2 )
h2
2. La tabla siguiente corresponde a la función de Bessel de primer orden, J1 (x):
x 0 0, 1 0, 2 0, 3 0, 4 0, 5
J1 (x) 0, 0000 0, 0499 0, 0995 0, 1483 0, 1960 0, 2423
x 0, 6 0, 7 0, 8 0, 9 1, 0
J1 (x) 0, 2867 0, 3290 0, 3688 0, 4059 0, 44018
Complétela con los valores aproximados de J1 (x), J1 (x)

3. Se busca una aproximación para f (x) que utilice los valores de f en tres puntos, digamos: x−1 , x0 , x1 ,
con x− 1− < x0 < x1 y x0 − x−1 = x1 − x0 = h. La fórmula debe ser del tipo
f (x0 ) ≈ Af(x−1 ) + Bf(x0 ) + Cf(x1 )
Se asume, además, que ésta debe ser exacta sobre los polinomios de grado menor o igual a dos.
(a) Calcule los coeficientes A, B, C.

h2
(b) Si |f | < M, se puede probar que el error de truncación verifica |e| < M. ¿Es esta fórmula
6
más exacta que otras que sólo usan dos puntos?
(c) De manera similar, deduzca la fórmula siguiente, la cual utiliza cinco puntos(de hecho, sólo cuatro).
1
f (x0 ) ≈ [f(x−2 − 8f(x−1 ) + 8f(x1 ) − f(x2 )]
12h
cuyo error de truncación es del orden de h4 .
(d) Finalmente, obtenga la fórmula de cinco puntos
1
f (x0 ) ≈ [−f(x−2 + 16f(x1 ) − 30f(x0 ) + 16f(x1 ) − f(x2 )]
12h2
4. Considere la función
x2 sen 792π si x = 0
f(x) =

100x
0 si x = 0
(a) Muestre que es derivable en R y calcule su derivada.

(b) Calcule f (0, 1) usando las fórmulas más usuales.
(c) Use un programa gráfico para trazar la curva de escuación y = f(x) y, a partir de esto, explique
los resultados obtenidos.
1
5. Un sistema resorte-masa tiene resistencia al movimiento proporcional al cuadrado de la velocidad:
d2 x dx
2
+ 0, 1 + 0, 6x = 0 x(0) = x (0) = 0
dt2 dt
Investigue sobre el uso de series de Taylor para resolver esta clase de ecuaciones y encuentre una
solución aproximada.
6. Dada la ecuación y = 0, 04y, y(0) = 100
(a) Calcule y(1) mediante el método de Euler, con h = 1, h = 0, 5, h = 0, 25.

(b) Calcule y(1) a partir de la solución y = 1000e 0,04t
(c) Interprete los resultados anteriores.
7. Resuelva la ecuación
y = −10y
y(0) = 1
usando RK2 y paso h = 0, 2. ¿Qué fenómeno se observa?
8. Use RK2 para calcular una solución aproximada en [0, 1] de la ecuación diferencial:
dy 1
= , y(0) = 0
dx x+y
9. Un paracaidista salta de un aeroplano. Durante el tiempo que tarda en abrir su paracaídas, la resisten-
cia del aire es proporcional a la raíz cuadrada de su velocidad. Se sabe que, bajo estas condiciones,
la velocidad máxima de caída es 80 mph. Determine por RK2 su velocidad durante los primeros dos
segundos (use ∆t = 0, 2).
10. La tasa de emisión radiactiva de una substancia es proporcional a la cantidad de materia restante. La
ecuación diferencial es, entonces y = −λy. Suponga que λ = 0, 01 y que se tienen 100 g de material
en el instante t = 0. ¿Cuánto material restará cuando t = 100?
La solución de la ecuación es y = 100e−λt y la respuesta exacta es 36,788 g. Resuelva la ecuación
numéricamente, usando
(a) Euler, con h = 25, h = 10, h = 5, h = 1

(b) Euler mejorado, con h = 20, h = 10
(c) Euler modificado, con h = 20, h = 10
(d) RK4, con h = 100, h = 50
(e) Predictor-corrector de orden 2, con h = 20.
Compare los resultados anteriores.
11. Un cuerpo con una masa inicial de 200 slugs es acelerado con una fuerza constante de 2000 lb; su masa
decrese a una tasa de 1 slug/s. Supongamos que éste está sujeto a una resistencia igual al doble de su
velocidad y que, en el instante t = 0, se encuentra en reposo. La ecuación diferencial es
dv 2000 − 2v
=
dt 200 − t
y la solución es v = 10t − t2 /40. Calcule v(50), usando
(a) Euler, Euler mejorado y Euler modificado, con h = 10

(b) RK4, con h = 10
(c) Predictor corrector, con h = 10.
Compare los resultados obtenidos.
2
12. Escriba programas Matlab para resolver y = f(t, y), y(x0 ) = y0 , por los métodos de Euler, Euler
modificado, RK4 y predictor-corrector. Cada programa debe
(a) Suponer la definición de f(t, y) en un programa externo.

(b) Leer valores x0 , y0 (condición inicial), h (el paso), xF (último valor de x).
(c) El programa para el predictor-corrector debe usar el método RK4 para encontrar los valores
requeridos.
13. Para la ecuación diferencial:

dy
= t2 − t, y(0) = 1
dt
se sabe que la solución es y = t3 /3 − t2 /2 +1/6. Si usamos un método de Adams_Bashford de orden
2, ¿cuál es el error esperado? Confirme esta predición realizando los cálculos correspondientes.
x
14. Para la ecuación diferencial y = se tienen los siguientes valores
y
x 0 1
√ 2 √ 3√
y 1 2 5 10
¿Cuántos decimales exactos dará el método de Milne en x = 4? ¿Con cuántos decimales exactos deben
escribirse los valores iniciales para asegurar esta precisión?
15. Resuelva el problema de valores en la frontera, tomando h = 0, 2:
y + xy − xy = 2x, y(0) = 1, y(1) = 0
16. Dado el problema de valores en la frontera:
y + y = 0, y(0) = 0, y(π/2) = 1,
normalice al intervalo [0, 1] y luego resuelva el problema por un método de diferencias finitas al orden
2. Use h = 0, 25.
17. Resuelva el problema de condiciones en la frontera
y + y = 0, y(0) + y (0) = 1, y(π/2) + y (π/2) = 0,
subdividiendo el intervalo en cuatro partes iguales.

18. Escriba el algoritmo RK2 para el sistema diferencial siguiente (asuma valores iniciales x0 , y0 ) :
x = x + 2y
y = 5x + 4y, x(0) = x0 , y(0) = y0
19. Reduzca la ecuación y + y = F (x) a un par de ecuaciones de primer orden y muestre que el lado
derecho de cada una es independiente de la variable derivada a la izquierda. ¿Esto significa que la
fórmula correctora no requiere iteraciones?
20. Una viga de longitud L en cantilever, que se curva bajo la acción de una fuerza uniforme de w lb/pies,
sufre también la acción de una fuerza axial P en su lado libre. La ecuación de la curva elástica es
solución de
d2 y wx2
EI = Py −
dx 2 2
dy
y(0) = (0) = 0
dx
Resuelva este problema por diferencias finitas, en el caso de una viga de madera de 2 × 2 × 10 pies

para la cual E = 12 × 105 lb/pul2 . Suponga que la deflexión vertical es de 4 , que w = 25 lb/pies y

que P = 500 lb.
3
21. Considere la ecuación diferencial elíptica
∂ 2 u ∂ 2 u ∂u
2 + 2 − =2
∂x2 ∂y ∂x
Utilice el método de diferencias finitas sobre el cuadrado unitario, usando como condición de borde
u = 0 y h = 0, 2 para escribir un sistema de ecuaciones lineales que genere una solución aproximada.
Aproxime la derivada de primer orden por diferencias centrales.

Taller 6 - Metodos Numericos

Cargado por

Copyright:

Formatos disponibles

Taller 6 - Metodos Numericos

Cargado por

Información del documento

Descripción 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

Taller 6 - Metodos Numericos

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD AUSTRAL DE CHILE Bachillerato en Ciencias de la Ingeniería

FACULTAD DE CIENCIAS Semestre 4

1 Derivadas y ecuaciones diferenciales

−f(x + 2h) + 4f(x + h) − 3f(x)

Complétela con los valores aproximados de J1 (x), J1 (x)

f (x0 ) ≈ Af(x−1 ) + Bf(x0 ) + Cf(x1 )

(a) Calcule los coeficientes A, B, C.

(a) Muestre que es derivable en R y calcule su derivada.

(a) Calcule y(1) mediante el método de Euler, con h = 1, h = 0, 5, h = 0, 25.

(c) Interprete los resultados anteriores.

(a) Euler, con h = 25, h = 10, h = 5, h = 1

(a) Euler, Euler mejorado y Euler modificado, con h = 10

(a) Suponer la definición de f(t, y) en un programa externo.

13. Para la ecuación diferencial:

y + xy − xy = 2x, y(0) = 1, y(1) = 0

16. Dado el problema de valores en la frontera:

y + y = 0, y(0) = 0, y(π/2) = 1,

y + y = 0, y(0) + y (0) = 1, y(π/2) + y (π/2) = 0,

subdividiendo el intervalo en cuatro partes iguales.

que P = 500 lb.

También podría gustarte