Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

7.1 Metodo Euler, Euler-Meorado, Runge-Kuttia

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 30

Ecuaciones Diferenciales

Tema 1. Parte 2: Métodos Numéricos para


Ecuaciones Diferenciales
Ester Simó Mezquita
Matemática Aplicada IV

1 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


Tema 1. Parte 2: Métodos numéricos para
Ecuaciones Diferenciales

1. Introducción
2. El método de Euler
3. El término de error
4. Método de Euler mejorado
5. Método de Runge-Kutta

2 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


1. Introducción
A estas alturas del curso un estudiante podría pensar que la
mayoría de las ecuaciones diferenciales pueden resolverse
explícitamente, con la solución de una fórmula dada

Aunque es posible demostrar de forma abstracta que casi cualquier


EDO posee una solución, por lo menos localmente, en general
resulta muy difícil expresar explícitamente de que solución se
trata

Pero lo primordial es que muchas de las ecuaciones que debemos


resolver en ingeniería no poseen soluciones de forma cerrada

Por ejemplo, las ecuaciones que rigen la forma del ala de un avión
no pueden resolverse. Y sin embargo, se vuela a diario

3 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


1. Introducción
La llegada de los ordenadores de alta velocidad ha hecho viable
y fácil llevar a cabo aproximaciones numéricas de las
soluciones

Las soluciones se obtienen con cualquier grado de exactitud, se


trazan gráficas y se lleva a cabo cualquier análisis que se desee

Pero los métodos numéricos jamás deben emplearse de forma


aislada. Siempre que sea posible, el usuario de estos métodos
debería utilizar técnicas cualitativas e intentar determinar si la
solución está acotada, si es estable, ¿cómo son sus asíntotas en el
infinito? ¿cómo se relacionan las diferentes soluciones entre sí?

De esta manera, los ingenieros no utilizan los métodos numéricos a


ciegas, sino, más bien, lo hace para brindar argumentos a su
entendimiento
4 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales
1. Introducción

Los métodos numéricos para resolver EDO tienen dos características


que se han de tener en cuenta:

1. Sólo permiten hallar soluciones particulares. Por lo tanto, para


poderlos aplicar, hará falta dar un conjunto completo de condiciones
iniciales

2. Necesitamos que las EDO o el conjunto de EDO que se les pasa


sean todas de primer orden

5 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


2. El método de Euler
Consideremos un problema de valores iniciales

Podemos integrar de a para obtener

Ya que la función desconocida se presenta en el integrando a la derecha,


no podemos proceder , a menos que contemos con un método de aproximación
de la integral.

El método de Euler se obtiene a partir de la técnica más simple para aproximar


la integral

6 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


2. El método de Euler
Problema de valores iniciales

Supongamos que el integrando no varía mucho en el intervalo 


resultará un error muy pequeño si reemplazamos por su valor en el
punto extremo izquierdo.

Colocando en su lugar una partición


del intervalo que se estudia. Supongamos que cada intervalo
tiene longitud

7 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


2. El método de Euler
Problema de valores iniciales

Sobre la base de estos cálculos definimos

Continuando de esta manera, y estableciendo que definimos

En general, estableciendo que , definimos

8 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


2. El método de Euler

Método de Euler

Fijado , es posible obtener aproximaciones de la solución


del problema de valores iniciales

en los puntos donde

mediante el método recurrente

son los puntos de nuestra solución aproximada


de la ecuación diferencial

9 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
2. El método de Euler
son los puntos de
nuestra solución aproximada de la ecuación diferencial

Solución exacta

Solución aproximada

10 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
2. El método de Euler

Para obtener una forma cómoda de medir el comportamiento de la técnica


numérica que se emplea definiremos el error relativo local en el n-esimo
paso

Normalmente esta cantidad vendrá dada como un porcentaje

11 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
2. El método de Euler
Ejemplo
Apliquemos la técnica de Euler a la EDO

utilizando incrementos de longitud y

Nuestro punto de referencia será calcular numéricamente y compararlo


con el valor exacto de que obtendremos aplicando el correspondiente
método de EDO

12 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
2. El método de Euler

La tabla muestra los cálculos para

El porcentaje de error

13 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
2. El método de Euler

La tabla muestra los cálculos para

Los datos muestran que cuando se reduce el incremento, la exactitud


mejora, pero el inconveniente es que se requieren más cálculos
14 Tema 1 Métodos Numéricos para
Ecuaciones Diferenciales
2. El método de Euler
El método de Euler se puede utilizar para resolver un número arbitrario
de EDO de primer orden.

Por ejemplo, si tenemos el sistema

los valores que se han de calcular de forma recurrente son

15 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
2. El método de Euler
Todos los métodos para resolver EDO trabajan con EDO de primer
orden.

Si tenemos una de orden superior  la tenemos que reescribir como un


sistema de EDO de orden 1.

Por ejemplo, sea la segunda ley de Newton en una dimensión con fuerza
arbitraria

Considerando la velocidad como una nueva variable


dependiente, además de , tendremos

M. de
Euler

16 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
3. El término de error

La noción de error es fundamental en cualquier técnica numérica.

Asociado al hecho de hacer muchas operaciones si es grande, está el


problema del error de redondeo

Los números reales no pueden representarse exactamente en un ordenador y


se han de redondear. Eso quiere decir que, cada vez que se hace una
operación, es posible que se pierdan dígitos del resultado, y en principio,
cuantas más operaciones más información se va perdiendo

Además, el método de Euler, introduce por sí mismo un error, que se llama


error de truncamiento. Los dos tipos de errores se mezclan, y de hecho, el
error total se puede amplificar

Analicemos a continuación el error que se comete al aproximar aplicando el


método de Euler

17 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
3. El término de error

El error de truncamiento local en el n-ésimo paso se define como

donde es el valor exacto en la de la ecuación diferencial e


es la aproximación de Euler

Podemos emplear la fórmula de Taylor para obtener una aproximación útil


de este término de error

para algún valor entre y .

Sabemos por la EDO que 

18 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
3. El término de error
El error de truncamiento local en el n-ésimo paso se define como

donde es el valor exacto en la de la ecuación diferencial y


es la aproximación de Euler

Podemos concluir

el error es proporcional a

19 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
3. El término de error

El error total de truncamiento para ir de a en pasos del


método de Euler será

Donde es la media de los y

Por lo tanto, el error total de truncamiento al aplicar el método de


Euler es proporcional al paso

20 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
4. Método de Euler mejorado
Mejoremos el método de Euler

Recordemos que nuestra antigua ecuación era

La idea del método de Euler consistía en reemplazar el integrando por


(aproximar la integral por medio de un área de un rectángulo)

Ahora proponemos reemplazar el integrando por

Así tenemos

El problema con la ecuación propuesta consiste en que se desconoce, debido


a que desconocemos la solución exacta

Lo que podemos hacer es reemplazar por su valor aproximado


determinado por el método de Euler
21 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales
4. Método de Euler mejorado
Representemos este nuevo valor por medio de
Así la ecuación se convierte

En general nuestro esquema de recurrencia es

en el cual

Este método que recibe el nombre de método mejorado de Euler o método de


Heun, primero predice y luego corrige una aproximación de

22 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


4. Método de Euler mejorado
Método de Euler mejorado

Fijado , es posible obtener aproximaciones de la solución


del problema de valores iniciales

en los puntos donde

mediante el método recurrente

Donde

El error global de truncamiento del método es del orden de


23 Tema 1 Métodos Numéricos para
Ecuaciones Diferenciales
4. Método de Euler mejorado
Ejemplo
Apliquemos el método de Euler mejorado a la EDO

utilizando un incrementos de longitud y midamos la mejora en


exactitud con respecto al método ordinario de Euler

24 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
4. Método de Euler mejorado

El valor aproximado obtenido para es 3.40542 . El error relativo es del 1%,


mientras que con el método de Euler es del 13%
25 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales
5. Método de Runge-Kutta
Mejoremos el método de Euler mejorado

Recordemos que nuestra antigua ecuación era

La idea del método de Runge-Kutta consiste en aproximar la integral


sustituyendo el integrando por una parábola

26 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


5. Método de Runge-Kutta
Método de Runge-Kutta de cuarto orden

Fijado , es posible obtener aproximaciones de la solución


del problema de valores iniciales

en los puntos donde

mediante el método recurrente

Donde

El error global de truncamiento del método es del orden de


27 Tema 1 Métodos Numéricos para
Ecuaciones Diferenciales
5. Método de Runge-Kutta
Ejemplo
Apliquemos el método de Runge-Kutta a la EDO

utilizando un incrementos de longitud

28 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales
5. Método de Runge-Kutta

29 Tema 1 Métodos Numéricos para Ecuaciones Diferenciales


6. Bibliografía

1. Simmons, G.F., Krantz, S.G., Ecuaciones diferenciales. Teoría, técnica


y práctica. McGraw-Hill Interamericana, 2007. ISBN 978-0-07-286315-4

2. Batlle, C., Massana, I., Zaragozá, M., Àlgebra i Equacions diferencials,


Edicions UPC, 2000. ISBN 84-8301-405-X

3. Batlle, C, Apunts tema 5 – Mètodes numèrics per a equacions


diferencials, Atenea-Campus Digital, 2010

30 Tema 1 Métodos Numéricos para


Ecuaciones Diferenciales

También podría gustarte