Blogia

ANÁLISIS NUMÉRICO

UNIVERSIDAD NACIONAL DE LOJA

 

CARRERA DE INGENIERÍA EN SISTEMAS

 

MÓDULO 8

 

ANÁLISIS NUMÉRICO

 

 

INTEGRACIÓN NUMÉRICA

 

 

El presenta tema tiene objetivo mostrar las utilidades de una herramienta de programación como es Maple11 para resolver problemas de cálculo numérico, esta fabulosa herramienta nos permite realizar desde cálculos de operaciones básicas como operaciones con matrices, vectores, funciones trigonométricas, funciones de integración y más. Además nos permite representar los resultados de estas operaciones en gráficas unidimensionales y 3D.

 

El presente artículo muestra como podemos aplicar esta herramienta para resolver problemas de integración numérica mediante los métodos de Simpson y la de los trapecios.

 

Primeramente ponemos a disposición la fundamentación teórica del tema.

Esperamos que esta información les sirva y esperamos sus comentarios.

 

 

Autoras: Jhomara Luzuriaga, Doris Viñamagua, Johanna Torres

 

 

Integración Numérica

Integración Numérica

La integración numérica constituye una amplia gama de algoritmos para calcular el valor numérico de una integral definida y, por extensión, el término se usa a veces para describir algoritmos numéricos para resolver ecuaciones diferenciales.

El problema básico considerado por la integración numérica es calcular una solución aproximada a la integral definida:


Este problema también puede ser enunciado como un problema de valor inicial para una ecuación diferencial ordinaria, como sigue:


Encontrar y(b) es equivalente a calcular la integral. Los métodos desarrollados para ecuaciones diferenciales ordinarias, como el método de Runge-Kutta, pueden ser aplicados al problema reformulado. En este artículo se discuten métodos desarrollados específicamente para el problema formulado como una integral definida.

Hay varias razones para llevar a cabo la integración numérica. La principal puede ser la imposibilidad de realizar la integración de forma analítica. Es decir, integrales que requerirían un gran conocimiento y manejo de matemática avanzada pueden ser resueltas de una manera más sencilla mediante métodos numéricos. Incluso existen funciones integrables pero cuya primitiva no puede ser calculada, siendo la integración numérica de vital importancia. La solución analítica de una integral nos arrojaría una solución exacta mientras que la solución numérica nos daría una solución aproximada. El error de la aproximación, que depende del método que se utilice y de qué tan fino sea, puede llegar a ser tan pequeño que es posible obtener un resultado idéntico a la solución analítica en las primeras cifras decimales.

Método del Trapecio

Las integrales son el area debajo de una curva delimitada por 2 puntos sobre el eje de las x, un punto inicial o punto a y un punto final o punto b.
Al derivar una funcion y evaluarla entre esos 2 puntos se obtiene el area que esta debajo de la curva, como lo muestra la imagen, se puede conocer el area sombreada que se muestra integrando la funcion y evaluandola entre los 2 puntos.




(Se supone que la base inferior de todos los trapecios mide lo mismo)
y lo que se hace es sacar el area de cada uno de los trapecios y luego sumarlas todas para obtener mas o menos el area sombreada.
Mientras mas trapecios se coloquen mas precisa sera la aproximacion, y tambien mientras mas trapecios se coloquen su base inferior tiende a ser 0.


 

 

 

El area de un trapecio se obtiene como muestra la imagen de arriba, asi que solo hay que hacer un pequeño analisis matematico para obtener una ecuacion que nos permita obtener el area de los trapecios en la curva:

Se separa al trapecio como si estuviera compuesto de un triangulo y un rectangulo.
La serie de ecuaciones que se muestran abajo son un desarrollo en el cual se trata al trapecio priero sacando el area del rectangulo y luego la del triangulo y sumandolas. Luego se simplifican las ecuaciones y se obtiene la que estamos buscando.

El area del rectangulo es bxh y la base es la resta entre los puntos x1-x0 la altura es el punto y1.
El area del triangulo es (bxh)/2
la base es la misma que la del rectangulo y la altura es la resta del punto y0-y1.
Se factorizan y simplifican las ecuaciones.

Por ultimo queda la ecuacion del area de cada trapecio.


 

 


Pero no todas las funciones son faciles de integrar, de hecho algunas ni siquiera se pueden integrar, asi que se usan otros metodos que se aproximan a una integral.

 

REGLA DEL TRAPECIO

Es la primera de las fórmulas de integración cerrada de Newton–Cotes. Corresponde al caso donde el polinomio en la ecuación de integración es de primer orden:

 

Geométricamente, es equivalente a aproximar el área del trapezoide bajo la línea recta que conecta f(a) y f(b).

La integral se representa como:

I ≈ ancho x altura promedio

Error de la regla trapezoidal:

Una estimación para el error de truncamiento local de una sola aplicación de la regla trapezoidal es:

 

 

 


 

 

 

El metodo del trapecio es uno de los mas antiguos y sirvio de base para desarrollar las integrales(junto con otros metodos parecidos como el de los rectangulos), pero nunca seran exactos, solo se aproximan.
Este metodo consiste en "llenar" el area que se quiere conocer con trapecios como lo muestra la imagen:

 

 

http://www.scribd.com/doc/4755474/mat31

Método de Simpson

En este procedimiento, se toma el intervalo de anchura 2h, comprendido entre xi y xi+2, y se sustituye la función f(x) por la parábola que pasa por tres puntos (xi, yi), (xi+1, yi+1), y (xi+2, yi+2). El valor del área aproximada, sombreada en la figura, se calcula con un poco más de trabajo y el resultado es

El método de Simpson deberá ser mucho más exacto que el procedimiento del trapecio. El área aproximada en el intervalo [a, b] es

bien, agrupando términos

El primer paréntesis, contiene la suma de los extremos, el segundo, la suma de los términos de índice impar, y el tercero la suma de los términos de índice par. En el método de Simpson, el número de divisiones n debe de ser par. En el caso de que el usuario introduzca un número impar el programa lo convierte en el número par siguiente.

 REGLAS DE SIMPSON

Regla de Simpson 1/3: Resulta cuando una interpolación polinomial de segundo orden es sustituida en una ecuación de aproximación con integral:

 

Después de la integración y manejo algebraico, resulta:

Donde h = (b – a)/2

Regla de Simpson 3/8: Resulta cuando una interpolación polinomial de tercer orden es sustituida en la ecuación de aproximación:

Para obtener:

 

http://www.scribd.com/doc/5784437/Metodo-Simpson

VIDEO INTEGRACION NUMERICA

VIDEO REGLA DE LOS TRAPECIOS

VIDEO REGLA DE SIMPSON