Integración Numérica
Integración Numérica
Integración Numérica
Los mtodos de integracin numrica nos permiten integrar funciones que estn de finidas analticamente o de los que solo conocemos su tabla en un numero finito de puntos.
Teorema Fundamental del Clculo Sea de una funcin contnua en el intervalo . Entonces: y sea una antiderivada
El problema en la prctica, se presenta cuando nos vemos imposibilitados de encontrar la antiderivada requerida, an para integrales aparentemente sencillas como:
La cual simplemente es imposible de resolver con el Teorema Fundamental del Clculo. FORMULAS DE INTEGRACION DE NEWTON-COTES Las frmulas de Newton-Cotes estn 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 mtodo conocido como mtodo de extrapolacin 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 mtodos mencionados arriba. Estas frmulas se basan en la idea de integrar una funcin polinomial en vez de :
donde
para ciertos datos de que se escogen apropiadamente. Es importante observar que estas frmulas se pueden aplicar inclusive a una tabla de datos, ya que lo que se usa es un polinomio de interpolacin, el cual puede ser calculado con la tabla. Dentro de las frmulas de Newton-Cotes, existen las formas cerradas y abiertas. En las formas cerradas se conocen los valores de llaman formas abiertas. REGLA DEL TRAPECIO Corresponde al caso donde , es decir : y ; en caso contrario, se
donde datos:
Que es la conocida Regla del Trapecio. Este nombre se debe a la interpretacin geomtrica que le podemos dar a la frmula. El polinomio de interpolacin para una tabla que contiene dos datos, es una lnea recta. La integral, corresponde al rea bajo la lnea recta en el intervalo forma. , que es precisamente el rea del trapecio que se
donde es el polinomio de interpolacin para los datos en la tabla anterior. Usaremos el polinomio de Lagrange. As, tenemos que:
Si denotamos
, entonces:
Simplificando trminos:
Vemos que cada uno de los trminos 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 frmula para calcular la integral de cada uno de los tres trminos de .
Debido al factor
se le conoce como la regla de Simpson de un tercio. para obtener nuestra frmula final:
donde
Y donde
al intervalo . Igual que en el caso anterior, se usa el polinomio de interpolacin de Lagrange, y usando el mtodo de integracin por partes se llega a la siguiente frmula:
donde . Debido al factor es que se le di el nombre de Regla de Simpson de 3/8. En la prctica, se sustituye el valor de h para obtener:
Sea
donde es el error de truncamiento que se comete al aplicar la regla. El mtodo de extrapolacin de Richardson combina dos aproximaciones de integracin numrica, para obtener un tercer valor ms exacto.
El algoritmo ms eficiente dentro de ste mtodo, se llama Integracin de Romberg, la cual es una frmula recursiva. Supongamos que tenemos dos aproximaciomnes : e
Se puede demostrar que el error que se comete con la regla del trapecio para n subintervalos est dado por las siguientes frmulas:
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 :
Esta frmula es solo una parte del algoritmo de Romberg. Para entender el mtodo, es conveniente pensar que se trabaja en niveles de aproximacin. En un primer nivel, es cuando aplicamos la regla del Trapecio, y para poder usar la frmula anterior, debemos de duplicar cada vez el nmero de subintervalos: as, podemos comenzar con un subintervalo, luego con dos, cuatro, ocho, etc, hasta donde se desee. Posteriormente, pasamos al segundo nivel de aproximacin, que es donde se usa la frmula anterior, tomando las parejas contiguas de aproximacin del nivel anterior, y que corresponden cuando . Despus pasamos al nivel tres de aproximacin, pero aqu cambia la frmula de Romberg, y as sucesivamente hasta el ltimo nivel, que se alcanza cuando solo contamos con una pareja del nivel anterior. Desde luego, el nmero de niveles de aproximacin 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 aproximacin n. Hacemos un diagrama para explicar un poco ms lo anterior.
Ejemplos: TRAPECIO Ejemplo 1: Utilizar la regla del trapecio para aproximar la integral:
Solucin. Igual que en el ejemplo anterior, sustitumos los datos de manera directa en la frmula del trapecio. En este caso, tenemos los datos:
La regla del trapecio se puede ampliar si subdividimos el intervalo subintervalos, todos de la misma longitud .
en
Sea la particin que se forma al hacer dicha subdivisin. Usando propiedades de la integral tenemos que:
Ahora bien, ya que todos los subintervalos tienen la misma longitud h, tenemos que:
Esta es la regla del trapecio para n subintervalos. Obviamente, esperamos que entre ms subintervalos usemos, mejor sea la aproximacin a la integral. Ejemplo 1: Aplicar la regla del trapecio para aproximar la integral
si subdividimos en 5 intervalos.
= 1.48065 Cabe mencionar que el valor verdadero de esta integral es de 1.4626 As, vemos que con 5 intervalos, la aproximacin no es tan mala. Para hacer clculos con ms subintervalos, es conveniente elaborar un programa que aplique la frmula con el nmero de subintervalos que uno desee. El lector debera hacer su propio programa y checar con 50, 500, 1000, 10000 y 20000 subintervalos, para observar el comportamiento de la aproximacin. SIMPSON DE UN TERCIO Ejemplo 1. Usar la regla de Simpson de 1/3 para aproximar la siguiente integral:
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
por el punto medio en cada subintervalo. Aplicamos primero propiedades bsicas de la integral definida:
Ahora, aplicamos la regla de Simpson de 1/3, en cada una de las integrales de arriba:
Sustitumos
Adems, los puntos medios de cada subintervalo son: Por lo tanto, sustitumos los datos en la frmula para obtener:
Ntese que esta aproximacin ya es exacta hasta el cuarto decimal! Ejemplo 2. Aproximar la siguiente integral, utilizando la regla de Simpson de en 4 intervalos. y subdividiendo
SIMPSON DE TRES OCTAVOS Ejemplo 1. Aproximar la siguiente integral, usando la regla de Simpson de :
Al igual que en los dos casos anteriores, la regla de Simpson de 3/8, se puede extender si subdividimos el intervalo . Sea la particin determinada de esta forma. Cada subintervalo y los puntos determinados as: en intervalos de la misma longitud
Aplicando la regla de
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. Solucin. Identificamos y la particin correspondiente:
Al considerar los puntos que dividen en tres partes iguales a cada subintervalo, tenemos los siguientes datos:
Ejemplo 1.
Evaluar
Solucin. Vemos que en el intervalo podemos aplicar la regal del trapecio, en el intervalo la regal de Simpson de
la regal de Simpson de 3/8 y en el intervalo 1/3. As, tenemos las siguientes integrales:
Ejemplo 2.
Calcula la integral
Solucin. En este caso, vemos que podemos aplicar la regla de Simpson de 1/3 en el intervalo , la regla del trapecio en el intervalo intervalo y la regla de Simpson de 3/8 en el
Vale la pena comentar que no siempre tiene que suceder que se apliquen exactamente las tres reglas. En realidad, esto depende de cmo se encuentran espaciados los intervalos de la tabla de datos. INTEGRACIN DE ROMBERG Ejemplo 1. Usar el algoritmo de Romberg, para aproximar la integral
usando segmentos de longitud . Solucin. Primero calculamos las integrales del nivel 1, usando la regla del trapecio para las longitudes de segmentos indicadas:
Ahora pasamos al segundo nivel de aproximacin donde usaremos la frmula que se dedujo anteriormente:
donde es la integral menos exacta (la que usa menos subintervalos) e es la ms exacta (la que usa el doble de subintervalos). En un diagrama vemos lo siguiente:
Para avanzar al siguiente nivel, debemos conocer la fmula correspondiente. De forma similar a la deduccin de la frmula,
se puede ver que la frmula para el siguiente nivel de aproximacin (nivel 3) queda como sigue:
donde: es la integral ms exacta es la integral menos exacta En el siguiente nivel (nivel 4) se tiene la frmula
As, podemos concluir que el valor de la aproximacin, obtenido con el mtodo de Romberg en el ejemplo 1, es:
donde
usando el mtodo de Romberg con segmentos de longitud , , , Solucin. Igual que arriba, primero usamos la regla del trapecio (con los valores de h indicados) para llenar el nivel 1. Tenemos entonces que:
A continuacin, usamos las frmulas de Romberg para cada nivel y obtenemos la siguiente tabla: