Unidad 6.1-Ecuaciones Diferenciales Ordinarias
Unidad 6.1-Ecuaciones Diferenciales Ordinarias
Unidad 6.1-Ecuaciones Diferenciales Ordinarias
dy
f (t , y ) (1)
dt
yi 1 yi * h (2)
f (ti , yi ) (3)
donde f(ti,yi) es la ecuación diferencial evaluada en ti y yi. Esta estimación puede ser sustituida en 2:
yi 1 yi f (ti , yi ) * h
(4)
- Ejemplo. Use el método de Euler para integrar y ' 4e0.8t 0.5 y desde t=0 a t=4 con un paso de 1. La
condición inicial en t=0 es y=2. Note que la solución exacta puede ser determinada analíticamente
como:
4 0.8t
y (e e 0.5t ) 2e 0.5t
1.3
Por eso,
y(1) = 2 + 3*(1) = 5
t yanal yEuler Er
0 2.00000 2.00000
1 6.19463 5.00000 19.28
2 14.84392 11.40216 23.19
3 33.67717 25.51321 24.24
4 75.33896 56.84931 24.54
MÉTODOS NUMÉRICOS - INGENIERÍA CIVIL - PROF: HARVETH GIL
Note que aunque el cálculo captura la tendencia general de la solución verdadera, el error es
considerable. Este error puede ser reducido utilizando un paso más pequeño.
- El Método de Runge-Kutta. Existen diferentes variaciones de los métodos RK, pero el más común es el
siguiente:
yi 1 yi * h (2)
donde φ es llamada la función de incremento, la cual puede ser interpretada como una representación
de la pendiente sobre el intervalo. La función incremento puede ser escrita de forma general como:
k1 f (ti , yi )
...
kn f (ti pn1h, yi qn1,1k1h qn1,2k2h ... qn1,n1kn1h)
donde las p y las q son constantes. Note que las k tienen relaciones recurrentes, es decir, k1 aparece
en la ecuación de k2, la cual aparece en la ecuación de k3, y así sucesivamente. Debido a que cada k es
una evaluación funcional, esta recurrencia hace el método RK eficiente para cálculos
MÉTODOS NUMÉRICOS - INGENIERÍA CIVIL - PROF: HARVETH GIL
computacionales. Note que el método RK de primer orden con n=1, es de hecho el método de Euler.
- Método clásico de RK de cuarto orden. El método más popular de RK es de cuarto orden aunque
existen numerosas versiones. La siguiente, es la forma más común utilizada y es por ello llamada el
método clásico:
1
yi 1 yi (k1 2k 2 2k3 k 4 ) * h
6 (6)
donde
k1 f (ti , yi )
1 1
k 2 f (ti h, yi k1h)
2 2
1 1
k3 f (ti h, yi k 2 h)
2 2
k4 f (ti h, yi k3h)
Solución. Para este caso, la pendiente al principio del intervalo es calculada como:
Este valor es usado para calcular un valor de y además de una pendiente en el punto medio:
Esta pendiente es usada para calcular otro valor de y además de otra pendiente en el punto medio:
Ahora, esta pendiente es usada para calcular un valor de y además de una pendiente al final del
intervalo:
Finalmente, las cuatro pendientes estimadas son combinadas para originar una pendiente promedio.
Esta pendiente promedio es posteriormente usada para hacer una última predicción al final del
intervalo:
Ejercicio. Resuelva el siguiente problema para un intervalo desde t=0 hasta 2, usando un paso de 0.01
donde y(0)=1.
MÉTODOS NUMÉRICOS - INGENIERÍA CIVIL - PROF: HARVETH GIL
dy
2 y t 2
dt
Obtenga la solución mediante el método de Runge-Kutta de cuarto orden, con cinco cifras significativas.
Muestre la solución en una gráfica.
- Ejemplo: Resuelva para la velocidad y la posición de un saltador de Bunge en caída libre mediante el
método de Euler. Asuma que a t=0, x==0, e integre hasta t=10s con un paso de 2s. Tome la
aceleración de la gravedad como g=9.81 m/s2, y la masa del saltador como 68.1kg con un coeficiente
de arrastre de 0.25 kg/m.
dx
dt
d C
g d 2
dt m
La solución analítica es:
gm gm
(t ) tanh t
Cd Cd
MÉTODOS NUMÉRICOS - INGENIERÍA CIVIL - PROF: HARVETH GIL
m gCd
x(t ) ln cosh t
Cd m
Calcule también el error.
Solución. Las ODE pueden ser usadas para calcular la pendiente en t=0, así:
dx
0
dt
d 0.25 2
9.81 0 9.81
dt 68.1
El método de Euler es entonces usado en t=2s:
x 0 0(2) 0
0 9.81(2) 19.62
el proceso puede ser repetido en t=4s:
x 0 19.62(2) 39.24
0.25
19.62 9.81 19.622 2 36.41368
68.1
Siguiendo el procedimiento obtenemos:
Referencias:
- J. Kiusalaas, Numerical methods in engineering with Matlab, Cambridge University Press (2005).
- S.C. Chapra, Applied numerical methods with Matlab, Mc Graw Hill (2008).
- C. Woodford y C. Phillips, Numerical methods with worked examples, Springer (1997).