¡Descarga Metodos numericos y más Apuntes en PDF de Ingeniería Matemática solo en Docsity!
INTEGRACIÓN NUMÉRICA
En los cursos de Cálculo Integral, nos enseñan como calcular una integral definida de una función contínua mediante una aplicación del Teorema Fundamental del Cálculo:
Teorema Fundamental del Cálculo
Sea una función contínua en el intervalo y sea una antiderivada
de. Entonces:
El problema en la práctica, se presenta cuando nos vemos imposibilitados de encontrar la antiderivada requerida, aún para integrales aparentemente sencillas como:
la cual simplemente es imposible de resolver con el Teorema Fundamental del Cálculo.
En este capítulo estudiaremos diversos métodos numéricos que nos permitirán obtener aproximaciones bastante exactas a integrales como la mencionada anteriormente. Esencialmente, veremos dos tipos de integración numérica: las fórmulas de Newton-Cotes y el algoritmo de Romberg.
Las fórmulas de Newton-Cotes están conformadas por las bien conocidas reglas del trapecio y de Simpson (regla de un tercio y de tres octavos). El algoritmo de Romberg forma parte de un método conocido como método de extrapolación de Richardson.
Haciendo uso de algunos programas computacionales (por ejemplo, en Mathematica ) es posible discernir sobre las cualidades y defectos de cada uno de los métodos mencionados arriba.
FORMULAS DE INTEGRACION DE NEWTON-COTES
Estas fórmulas se basan en la idea de integrar una función polinomial en vez
de :
donde es un polinomio de interpolación de grado
para ciertos datos de que se escogen apropiadamente.
Es importante observar que estas fórmulas se pueden aplicar inclusive a una tabla de datos, ya que lo que se usa es un polinomio de interpolación, el cual puede ser calculado con la tabla.
Dentro de las fórmulas de Newton-Cotes, existen las formas cerradas y abiertas.
En las formas cerradas se conocen los valores de y ; en caso contrario, se llaman formas abiertas.
Nosotros nos remitiremos a estudiar únicamente las formas cerradas, y por lo tanto,
siempre suponemos que conocemos los valores y.
REGLA DEL TRAPECIO
Corresponde al caso donde , es decir :
donde es un polinomio de interpolación (obviamente de grado 1) para los datos:
Del capítulo anterior, sabemos que este polinomio de interpolación es:
Integrando este polinomio, tenemos que:
Por lo tanto, tenemos que:
Que es la conocida Regla del Trapecio. Este nombre se debe a la interpretación geométrica que le podemos dar a la fórmula. El polinomio de interpolación para una tabla que contiene dos datos, es una línea recta. La integral, corresponde al área
bajo la línea recta en el intervalo , que es precisamente el área del trapecio que se forma.
Ejemplo 1: Utilizar la regla del trapecio para aproximar la integral:
Sustituyendo el valor de h y usando la notación sigma, tenemos finalmente:
Esta es la regla del trapecio para n subintervalos. Obviamente, esperamos que entre más subintervalos usemos, mejor sea la aproximación a la integral.
Ejemplo 1: Aplicar la regla del trapecio para aproximar la integral
si subdividimos en 5 intervalos.
Solución****.
En este caso, identificamos , y la partición generada es:
Así, aplicando la fórmula tenemos que:
Cabe mencionar que el valor verdadero de esta integral es de 1.4626…
Así, vemos que con 5 intervalos, la aproximación no es tan mala. Para hacer cálculos con más subintervalos, es conveniente elaborar un programa que aplique la fórmula con el número de subintervalos que uno desee. El lector debería hacer su propio programa y checar con 50, 500, 1000, 10000 y 20000 subintervalos, para observar el comportamiento de la aproximación.
REGLA DE SIMPSON DE UN TERCIO
Suponemos que tenemos los datos:
donde es el punto medio entre y.
En este caso se tiene que:
donde es el polinomio de interpolación para los datos en la tabla anterior. Usaremos el polinomio de Lagrange.
Así, tenemos que:
Si denotamos , entonces:
Simplificando términos:
Vemos que cada uno de los términos anteriores, es esencialmente de la misma
forma, es decir, una constante por
Así, calculamos la siguiente integral por partes:
Sea:
por lo tanto,
Usamos esta fórmula para calcular la integral de cada uno de los tres términos de
.
Ejemplo 2. Usar la regla de Simpson de 1/3, para aproximar la siguiente integral:
Solución****. Igual que en el ejercicio anterior, sustituímos datos adecuadamente:
Al igual que con la regla del trapecio, podemos extender la regla de Simpson de 1/3,
si subdividimos el intervalo en subintervalos de la misma longitud
Sea la partición que se forma al hacer la subdivisión, y
denotemos por el punto medio en cada subintervalo.
Aplicamos primero propiedades básicas de la integral definida:
Ahora, aplicamos la regla de Simpson de 1/3, en cada una de las integrales de arriba:
Sustituímos y usamos la notación sigma:
Ejemplo 1.
Aproximar la siguiente integral, aplicando la regla de Simpson de y subdividiendo en 5 intervalos.
Solución****.
En este caso, tenemos que , y la partición que se genera es:
Además, los puntos medios de cada subintervalo son:
Por lo tanto, sustituímos los datos en la fórmula para obtener:
Nótese que esta aproximación ya es exacta hasta el cuarto decimal!
Ejemplo 2.
Aproximar la siguiente integral, utilizando la regla de Simpson de y subdividiendo en 4 intervalos.
Solución****.
En este caso, tenemos que , y la partición que se genera es:
Además, los puntos medios de cada subintervalo son:
Sustituyendo todos estos datos en la fórmula obtenemos la siguiente aproximación:
REGLA DE SIMPSON DE TRES OCTAVOS
Este caso corresponde a , es decir,
donde es un polinomio de interpolación para los siguientes datos:
Y donde , y , son los puntos que dividen en tres partes iguales
al intervalo.
Igual que en el caso anterior, se usa el polinomio de interpolación de Lagrange, y usando el método de integración por partes se llega a la siguiente fórmula:
Esta última, es la regla de Simpson de 3/8 para n subintervalos todos de la misma longitud.
Ejemplo 2. Aproximar la siguiente integral:
aplicando la regla de Simpson de 3/8, y subdiviendo en 3 intervalos.
Solución****.
Identificamos y la partición correspondiente:
Al considerar los puntos que dividen en tres partes iguales a cada subintervalo, tenemos los siguientes datos:
Sustituyendo todos los datos en la fórmula, obtenemos:
De acuerdo a los ejemplos vistos, resulta evidente que la regla de Simpson de 3/8, es más exacta que la de 1/3 y a su vez, ésta es más exacta que la regla del trapecio. En realidad, pueden establecerse cotas para los errores que se cometen en cada uno de estos métodos.
Puesto que no es nuestra intención justificar formalmente cada uno de los teoremas, los siguientes resultados se mencionan para completar la información, pero omitimos las demostraciones correspondientes.
REGLA F O R M U L A E R R O R D O N D E****... Trapecio
Simpson
Simpson
INTEGRACIÓN EN INTERVALOS DESIGUALES
Cuando la longitud de los subintervalos no es igual, se usa una combinación de la regla Trapezoidal y las reglas de Simpson, procurando seguir el siguiente orden jerárquico:
1 .- Simpson
Esta se aplica, si contamos con 4 puntos igualmente espaciados.
2 .- Simpson
Esta se aplica si falla (1) y contamos con 3 puntos igualmente espaciados.
3 .- Regla Trapezoidal
Solo se aplica si no se cumple y
Ejemplo 1.
Evaluar , usando la siguiente tabla :
Solución****.
Vemos que en el intervalo podemos aplicar la regal del trapecio, en el
intervalo la regal de Simpson de 3/8 y en el intervalo la regal de Simpson de 1/3. Así, tenemos las siguientes integrales:
Finalmente, la integral buscada es la suma de las tres integrales anteriores:
Ejemplo 2.
Calcula la integral , usando la siguiente tabla de datos:
donde es un promedio de la doble derivada entre ciertos valores que pertenecen a cada uno de los subintervalos.
Ahora bien, si suponemos que el valor de es constante, entonces :
Sustituyendo esto último en nuestra primera igualdad, tenemos que:
De aquí podemos despejar :
En el caso especial cuando (que es el algoritmo de Romberg), tenemos :
Esta fórmula es solo una parte del algoritmo de Romberg. Para entender el método, es conveniente pensar que se trabaja en niveles de aproximación. En un primer nivel, es cuando aplicamos la regla del Trapecio, y para poder usar la fórmula anterior, debemos de duplicar cada vez el número de subintervalos: así, podemos comenzar con un subintervalo, luego con dos, cuatro, ocho, etc, hasta donde se desee.
Posteriormente, pasamos al segundo nivel de aproximación, que es donde se usa la fórmula anterior, tomando las parejas contiguas de aproximación del nivel anterior,
y que corresponden cuando.
Después pasamos al nivel tres de aproximación, pero aquí cambia la fórmula de Romberg, y así sucesivamente hasta el último nivel, que se alcanza cuando solo contamos con una pareja del nivel anterior.
Desde luego, el número de niveles de aproximación que se alcanzan, depende de las aproximaciones que se hicieron en el nivel 1. En general, si en el primer nivel, iniciamos con n aproximaciones, entonces alcanzaremos a llegar hasta el nivel de aproximación n.
Hacemos un diagrama para explicar un poco más lo anterior.
Ejemplo 1. Usar el algoritmo de Romberg, para aproximar la integral
usando segmentos de longitud.
Solución****. Primero calculamos las integrales del nivel 1, usando la regla del trapecio para las longitudes de segmentos indicadas:
Con estos datos, tenemos:
Ejemplo 2. Usar el algoritmo de Romberg para aproximar la integral:
Agregando a la tabla anterior donde.
Solución****.
Calculamos con la regla del trapecio:
Tenemos entonces la siguiente tabla:
De donde concluímos que la aproximación buscada es:
Ejemplo 3. Aproximar la siguiente integral:
usando el método de Romberg con segmentos de longitud
Solución****. Igual que arriba, primero usamos la regla del trapecio (con los valores de h indicados) para llenar el nivel 1. Tenemos entonces que:
A continuación, usamos las fórmulas de Romberg para cada nivel y obtenemos la siguiente tabla:
De donde concluímos que la aproximación buscada es:
Podemos escribir una fórmula general para calcular las aproximaciones en cada uno de los niveles como sigue:
ALGORITMO DE INTEGRACIÓN DE ROMBERG
Los coeficientes en cada una de las fórmulas en el método de Romberg, deben sumar 1. Así se tiene la siguiente fórmula recursiva:
donde:
es la integral más exacta
es la integral menos exacta
y el indice k indica el nivel de integración o de aproximación. Por ejemplo, digamos
que , entonces tenemos:
que es nuestra fórmula del nivel 2 de aproximación.
Como todo proceso iterativo, éste se detiene cuando se obtiene una aproximación suficientemente buena. En este caso se pide que:
donde es la cota suficiente.
Ejemplo 1. Aplicar el algoritmo de integración de Romberg a la integral:
i) Dividiendo en un solo intervalo. ii) Dividiendo en 4 intervalos. Soluciones : i) 82.60511 ii) 76.
3. Usar la regla de Simpson 3/8 para aproximar,
i) Dividiendo en un solo intervalo. ii) Dividiendo en 4 intervalos. Soluciones : i) 2.76591 ii) 2.
4. Integrar las siguientes tablas de datos:
i)
ii) Soluciones : i) -17.11458 ii) 9.
5. Usar el algoritmo de integración de Romberg para aproximar,
i) Usando 1, 2 y 4 intervalos. ii) Agregando al inciso anterior, 8 intervalos. Soluciones : i) 9.156626413 ii) 9.
6. Aproxime la integral del ejercicio anterior, tomando como cota suficiente.
Solución****. 9.
