MATEMÁTICAS III

TEMA 8: Resolución numérica de E.D.O. de primer orden. Métodos de un

paso: Euler, Taylor, Runge-Kutta. Estabilidad, consistencia y convergencia.

Araceli Queiruga Dios

Departamento de Matemática Aplicada
E-mail: [email protected]

GRADO EN INGENIERÍA ELECTRÓNICA INDUSTRIAL Y AUTOMÁTICA

GRADO EN INGENIERÍA ELÉCTRICA
GRADO EN INGENIERÍA MECÁNICA
CONTENIDOS
(MATEMÁTICAS III. TEMA 8)

1. Resolución numérica de E.D.O. de primer orden. Métodos de un paso:

o Euler.
o Taylor.
o Runge-Kutta.
2. Estabilidad, consistencia y convergencia.
RESOLUCIÓN NUMÉRICA DE EDOs
(MATEMÁTICAS III. TEMA 8)

Definición:
Se llama problema de Cauchy o problema del valor inicial (P.V.I.) al
conjunto formado por una EDO y una condición inicial:
 y ' = f ( x, y ), x  [a , b]

 y0 = y ( x0 )
➢ Métodos de paso simple o unipaso: se calcula y(xi+1) a partir de la

información proporcionada por y(xi).

➢ Métodos de paso múltiple o multipaso: se calcula y(xi+1) a partir de la
información proporcionada por los valores y(xi), …, y(xi-p), y del teorema
fundamental del cálculo, según el cual:
xj
y ( x j ) = y ( xi ) +  f ( x, y )dx
xi
RESOLUCIÓN NUMÉRICA DE EDOs: NOTACIÓN
(MATEMÁTICAS III. TEMA 8)

x0 = a, x1 ,..., xn −1 , xn = b es el conjunto discreto de puntos equidistantes

del intervalo [a, b]
xi = x0 + hi; i  0  n
b−a
h=
n
y(x) = curva solución
y(xi) = valor exacto; y' (xi) = f (xi, y(xi))
yi = valor aproximado; f = f (xi, yi) = yi
Error global: en el paso i-ésimo viene dado por Ei = |y(xi) - yi|.
Error local: (error cometido en cada paso) para el paso i-ésimo es
ei =|ŷ (xi) – yi|, donde ŷ(xi) es la solución exacta del PVI:
 y ' = f ( x, y )

 yi −1 = y ( xi −1 )
Matemáticas III. Ecuaciones Diferenciales Ordinarias
RESOLUCIÓN NUMÉRICA DE EDOs
(MATEMÁTICAS III. TEMA 8)

Definición:
Un método numérico es una ecuación en diferencias que contiene un cierto número de
aproximaciones sucesivas y consecutivas yn+j que permiten calcular secuencialmente la
sucesión yn.
Es decir, discretizamos la variable independiente t. De esta forma aproximaremos las
solución y(t) en un conjunto finito de puntos: t0, t1, ..., tn ∈ [a, b].
Llamamos y(ti) al valor exacto e yi al aproximado.
f (xi, yi) = fi

Matemáticas III. Ecuaciones Diferenciales Ordinarias

RESOLUCIÓN NUMÉRICA DE EDOs
(MATEMÁTICAS III. TEMA 8)

Clasificación:
Métodos de paso simple: se basan en el cálculo de yi+1 a partir de la información
proporcionada por yi.
Métodos multipaso: el cálculo de yi+1 se realiza a partir de la información
proporcionada por los valores de los puntos previamente calculados yi, … yi+1-k y por
el uso del Teorema Fundamental del Cálculo.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

RESOLUCIÓN NUMÉRICA DE EDOs
(MATEMÁTICAS III. TEMA 8)

Ejemplo:

yn + 2 + yn +1 − 2 yn =
h
4
( f ( xn + 2 , yn + 2 ) + 8 f ( xn +1 , yn +1 ) + 3 f ( xn , yn ) )

1 1 3 h
a) yn +3 + yn + 2 − yn +1 − yn = (19 f n + 2 + 5 f n )
2 2 4 8

b) yn + 2 − yn = h ( f ( xn + 2 , yn*+ 2 ) + f n )
h
yn*+ 2 = 3 yn +1 − 2 yn + ( f n +1 − 3 f n )
2

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODOS DE PASO SIMPLE
(MATEMÁTICAS III. TEMA 8)

Los métodos unipaso son aquellos en los que para obtener una aproximación de
y(xi+1), se utiliza la información proporcionada por y(xi).
1. Métodos de Taylor (según el grado del polinomio de Taylor que se elija)
a) Método de Euler.
b) Método de Taylor de orden k (con k >1).
2. Métodos de Runge-Kutta.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

 MÉTODOS DE TAYLOR
(MATEMÁTICAS III. TEMA 8)

Se basan en el hecho de que se puede obtener y(xi+1) mediante el desarrollo de Taylor

de y(x) en el punto xi

h2 hn (n
y ( xi +1 ) = y ( xi ) + hy '( xi ) + y ''( xi ) + + y ( xi ) + Rn , h = xi +1 − xi .
2! n!

Se sustituyen y, y', y'',... por los datos que nos dan en el enunciado

h2 h ( n −1
y ( xi +1 ) ≈ yi + hf ( xi , yi ) + f '( xi , yi ) + + f ( n ( xi , yi ).
2! ( n − 1)!

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODO DE EULER
(MATEMÁTICAS III. TEMA 8)

Se utiliza el polinomio de Taylor de primer orden.

Dado el PVI:
 y ' = f ( x, y ), x  [a , b]

 y0 = y ( x0 )

con solución única y(x). El método de Euler consiste en aproximar el valor de y en los
puntos xi = xi + ih por los valores de yi obtenidos a partir de y0 por el algoritmo:

yi +1 = yi + hf ( xi , yi ).
La cota del error local se calcula:
h2 M
ei  ; M = sup y '' ( )
2   ( xi −1 , xi )

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODO DE EULER: EJEMPLO
(MATEMÁTICAS III. TEMA 8)

Usar el método de Euler para aproximar la solución del P.V.I.:

4𝑥
𝑦 ′ = −𝑥𝑦 + , 𝑦 0 = 1, en 0 ≤ 𝑥 ≤ 1, con ℎ = 0,25
𝑦

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODO DE TAYLOR DE ORDEN k > 1
(MATEMÁTICAS III. TEMA 8)

Consiste en aproximar y(xi+1) por yi+1, donde

ℎ2 ′ ℎ𝑛−1
𝑦 𝑥𝑖+1 ≈ 𝑦𝑖 + ℎ𝑓 𝑥𝑖 , 𝑦𝑖 + 𝑓 𝑥𝑖 , 𝑦𝑖 + ⋯ + 𝑓 (𝑛 𝑥𝑖 , 𝑦𝑖 ,
2! 𝑛−1 !

para todo 0 ≤ 𝑖 ≤ 𝑛. Una cota para el error cometido con este método es

ℎ𝑘+1 𝑀
𝑒𝑖 ≤ , 𝑀= sup 𝑦 (𝑘+1 (𝜉) .
𝑘+1 ! 𝜉∈(𝑥𝑖−1 ,𝑥𝑖 )

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 1
(MATEMÁTICAS III. TEMA 8)

Resuelve por un método numérico la siguiente ecuación:

y' = x

 y (0) = 0
con h = 0.5, en el intervalo [0,5].
1. Con el siguiente método numérico: yn+1 – 2 yn = 0

2. Con el siguiente método numérico: yn+1 – 2 yn = h xn

3. Con el siguiente método numérico: yn+1 – yn = h xn

Calcula la solución exacta y compárala con los métodos anteriores.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 2
(MATEMÁTICAS III. TEMA 8)

Dada la siguiente ecuación, cuya solución exacta se conoce, obtener el valor

aproximado de y(1) mediante el método de Euler con h = 0,1. Buscad la razón del mal
comportamiento de la aproximación:

 y ' = 3( y + x )
1 4

 y (0) = 1
solución exacta: y ( x ) = − x − + e3 x
3 3

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 3
(MATEMÁTICAS III. TEMA 8)

Usa el método de Taylor de orden 2, con h = 0,1, para aproximar la solución del
siguiente PVI:

 y ' = 1 + x sen( xy ) 0  x  0.5


 y (0) = 0

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODOS DE PASO SIMPLE: RUNGE-KUTTA
(MATEMÁTICAS III. TEMA 8)

La idea general de Runge-Kutta es mejorar la aproximación

yi +1 = yi + hf ( xi , yi ).
del método de Euler, sin necesidad de calcular derivadas de orden superior. Esto puede
lograrse con aproximaciones de la forma

yi +1 = yi + h ( xi , yi , h )

Siendo  ( xi , yi , h ) una conveniente combinación lineal de derivadas de y (o sea, de

valores de f ) evaluadas en puntos “apropiados”.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODOS DE PASO SIMPLE: RUNGE-KUTTA
(MATEMÁTICAS III. TEMA 8)

(El orden indica el número de puntos en los que se evalúa f en cada paso).
➢ Métodos de Runge Kutta de 1 etapa: (s = 1) yn +1 = yn + hf
➢ Métodos de Runge Kutta de 2 etapas: (s = 2)
o Método modificado de Euler:

k1 = f ( xn , yn )

yn +1 = yn + hk2  h h
k
 2 = f ( x n + , y n + k1 )
2 2
o Método mejorado de Euler:

h k1 = f ( xn , yn )
yn +1 = yn + ( k1 + k2 ) 
2 k2 = f ( xn + h, yn + hk1 )

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODOS DE PASO SIMPLE: RUNGE-KUTTA
(MATEMÁTICAS III. TEMA 8)

(El orden indica el número de puntos en los que se evalúa f en cada paso).
➢ Método de Runge Kutta de orden 4:

h
yn +1 = yn + ( k1 + 2k2 + 2k3 + k4 ) , 0  i  n − 1,
6
Donde

k1 = f ( xi , yi ),
h h
k2 = f ( xi + , yi + k1 ),
2 2
h h
k3 = f ( xi + , yi + k2 ),
2 2
k4 = f ( xi + h, yi + hk3 ).

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJERCICIO 1
(MATEMÁTICAS III. TEMA 8)

Dado el PVI
1 y
y' = 2
− − y 2
, 1 x  2
x x
y (1) = −1
cuya solución exacta es y(x) = -1/x, usa un método de Runge-Kutta de orden 2
(Euler mejorado) con h = 0,1 para aproximar la solución y compárala con los valores
reales de y.

Resuelve la misma ecuación mediante otro método de Runge Kutta de orden 2 (Euler
modificado) y compara el resultado con el obtenido antes.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJERCICIO 3
(MATEMÁTICAS III. TEMA 8)

Resuelve mediante un método de Runge Kutta de orden 4 con h = 0.1, el siguiente

P.V.I.:
y ' = x 2 − 3 y, 0  x  0, 4
y (0) = 1

h
yn +1 = yn + ( k1 + 2k2 + 2k3 + k4 ) , 0  i  n − 1,
6
k1 = f ( xi , yi ),
 h h
k2 = f ( xi + , yi + k1 ),
 2 2

k3 = f ( xi + h , yi + h k2 ),
 2 2
k = f ( x + h, y + hk ).
 4 i i 3

Matemáticas III. Ecuaciones Diferenciales Ordinarias

MÉTODOS DE ADAMS-BASHFORD
(MATEMÁTICAS III. TEMA 8)

➢ Método Adams-Bashford de orden 2

h
yi + 2 = yi +1 + ( 3 fi +1 − fi )
2
➢ Método Adams-Bashford de orden 3
h
yi + 3 = yi + 2 + ( fi +2 + 4 fi +1 + fi )
3

➢ Método Adams-Bashford de orden 4

h
yi + 4 = y i + 3 + ( 55 fi +3 − 59 fi +2 + 37 f i +1 − 9 f i )
24

(El orden indica el número de puntos en los que se evalúa f en cada paso).

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 3
(MATEMÁTICAS III. TEMA 8)

Usa el método de Adams-Bashforth de orden 2, para resolver aproximadamente el

problema de PVI con h = 0,1:

 y ' = 1 + t sen(ty ) 0  x  0.5


 y (0) = 0

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLOS
(MATEMÁTICAS III. TEMA 8)

Sea y(t) la solución del problema de Cauchy y' = -2 t y2, y(0) = 1.

a) Calcula y(1), usando un método de Runge-Kutta de orden 4 con

h = 0.2.
b) Calcula y(1), usando el método de Adams-Bashford de orden 4 con h =
0.2, iniciándolo con los primeros valores del apartado anterior
c) Obtener la solución exacta de la EDO y comparar los resultados.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLOS
(MATEMÁTICAS III. TEMA 8)

dT
La ley de enfriamiento de Newton establece que = −k (T − T0 ), donde T(t) es la
dt
temperatura de un objeto en el instante t y T0 es la temperatura ambiente. Una taza de
café tiene una temperatura de 80ºC y se encuentra en una habitación con una
temperatura ambiente de 20ºC. En este caso k = 0,080. Obtened una tabla de
temperaturas del café en los siguientes 10 minutos, a intervalos de 30 segundos.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLOS
(MATEMÁTICAS III. TEMA 8)

Error local: (error cometido en cada paso) para el paso i-ésimo es

ei =|y (xi) – yi|, donde y(xi) es la solución exacta del PVI.
Error relativo local: (en el paso i-ésimo )
y ( xi ) − yi
ei =
y ( xi )

Ejemplo: Utilizar el método de Euler para aproximar el valor de la EDO

y' = x + y, y(0) = 1,
utilizando incrementos de longitud h = 0,2 y h = 0,1. Obtener el error relativo local en
cada paso (como %).
Sabiendo que la solución exacta es: y(x) = -1 - x + 2ex

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 3
(MATEMÁTICAS III. TEMA 8)

Usa el método de Taylor de orden 2, con h = 0.1 (n = 10), para aproximar la solución
del siguiente PVI:
 y ' = y, 0  x  0.5

 y (0) = 1,

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 3
(MATEMÁTICAS III. TEMA 8)

Un proyectil de masa m = 0,11 kg se lanza verticalmente hacia arriba con una velocidad
inicial v(0) = 8 m/s y se va frenando debido a la fuerza de la gravedad Fg = – mg y a la
resistencia del aire Fr = – kv2, donde g = 9,8 m/s2 y k = 0,002 kg/m.
La ecuación diferencial para la velocidad está dada por
mv' = – mg – kv2
Usar el método de Runge-Kutta de orden 4 para obtener la velocidad y determinar con
un error menor que una décima de segundo el instante en el que el proyectil alcanza su
altura máxima y empieza a caer.

Matemáticas III. Ecuaciones Diferenciales Ordinarias

EJEMPLO 3
(MATEMÁTICAS III. TEMA 8)

En un circuito eléctrico se dispone de un condensador con capacidad constante de

C = 1,1 faradios y una resistencia constante de R = 2,1 ohmios. Se le aplica un voltaje
E (t ) = e −0,06 tsen(2t −  )
Sabiendo que la intensidad en el instante inicial es I(0) = 1 amperio, obtener el valor
aproximado de la intensidad cada décima de segundo durante los 2 primeros segundos,
comparando las aproximaciones ofrecidas por los métodos de Euler y Runge-Kutta de
orden 2.

La ecuación que rige el comportamiento de un circuito eléctrico sin inducción es

dI I dE
R + =
dt C dt

Matemáticas III. Ecuaciones Diferenciales Ordinarias

BIBLIOGRAFÍA
(MATEMÁTICAS III. TEMA 8)

1. http://ocw.upm.es/matematica-aplicada/programacion-y-metodos-numericos/contenidos/
TEMA_2/Presentaciones/Errores_1.pdf/view
2. A. García, et al. Ecuaciones diferenciales ordinarias. Teoría y Problemas. Ed. Clagsa.
3. Á. Martín del Rey. Apuntes Cálculo Numérico, 2010.
4. J. H. Mathews, K. D. Fink, Métodos Numéricos utilizando Matlab, Prentice Hall, 3ª Edición,
2000.

Matemáticas III. Ecuaciones Diferenciales Ordinarias