Integracion Numerica

1.
INTEGRACIÓN NUMÉRICA
La integración numérica es una herramienta muy útil para obtener valores aproxima-
dos de las integrales definidas que no pueden calcularse analíticamente, ya sea por que la
integral no tenga primitivas expresados analíticamente ó que dicha integral no se conoce
en forma analítica ,si no en forma discreta tabulada tales como datos procedentes de un
experimento.
Sabiendo que una integral es el limite de una suma infinita ,es natural que su aproxi-
mación consista en una suma finita de muestras ponderadas con pesos wi del integrando
f (xi ), dicha suma denominada formula de cuadratura.
Definición: Sea f : [a, b] ⊂ R −→ R una función y una partición P = a = x0 < x1 < x2 < · · · < xn = b
del intervalo [a, b] con los xi los nodos en el eje X. Una expresión de tipo Q[ f (x)] =
n
∑ wi f (xi) = w1 f (x1) + w2 f (x2) + · · · + wn f (xn) de modo que .
i=1
Z b
f (x) dx = Q[ f (x)] + E[ f (x)] (1)
a
donde E[ f (x)] representa error de truncadura de la ecuación (1) Los valores {xi }ni=1 son
los nodos en el eje x de integración o cuadratura
Los nodos {xi }ni=1 = {x0 < x1 < x2 · · · < xn } se eligen de diferentes maneras ,dependiendo
que apliquemos en los criterios.
Para la regla del trapecio o la de simpson ,los nodos se toman igualmente espaciados o sea
x1 − x0 = x2 − x1 = x3 − x2 = · · · = xn − Xn−1 = ∆x = h para el método de Gauss-Legendre
los nodos que se forman son las raíces del polinomio de Legengre .
Un detalle importante en todas las aplicaciones es la estimación del error cometido y el
conocimiento del orden de precision del criterio de la regla utilizada .
En las reglas de cuadratura es frecuente que estén basados en la interpolación polinómica
,esto es :se sustituyen los valores de la función f (x) por el polinomio que interpola a f (x)
en los nodos y su integral es la aproximación buscada.

MÉTODO DE NEWTON COTES
Si un numero determinado n de nodos en el intervalo [a, b], esto es particionar en n subin-
tervalos .Siempre es posible determinar un polinomio de grado n y generar una formula
de cuadratura de n puntos.
si x1 − x0 = x2 − x1 = · · · = xn − xn−1 = ∆x = h osea cada subintervalo de [a, b] tiene la
misma longitud.
Las reglas resultantes se conocen como la fórmula cuadratura de Newton.
si se incluye los extremos de [a, b] de integración de dice fórmula cerrada, caso contrario
la fórmula es abierta.
a continuación nos remitimos al gráfico siguiente
Del gráfico se tiene que los xi −xi−1 = h para todo i = 0, 1, 2, · · · , n y que f (x0 ) , f (x1 ), · · · , f (xn )
son las bases del trapecio en sus respectivos nodos xi . Así la suma de las área de todos
los trapecios se aproximan al área acotado por el eje X, la recta x = a,la recta x = b y por
último la curva del la función f .
Si extraemos el i-ésimo trapecio de la figura, se tiene como área:
h
Ai = [ f (xi−1 ) + f (xi )]
2
luego en general la suma total de las áreas del trapecio es:
h h h h
[ f (x0 ) + f (x1 )] + [ f (x1 ) + f (x2 )] + · · · + [ f (xi−1 ) + f (x1 )] + · · · + [ f (xn−1 ) + f (xn )]
2 2 2 2
seguidamente factorizando h/2 se tiene:
h h
[ f (x0 ) + f (x1 ) + f (x1 ) + f (x2 ) + f (x2 ) + · · · + f (xn−1 ) + f (xn−1 ) + f (xn )] = [ f (x0 )+
2 2
2 f (x1 ) + 2 f (x2 ) + · · · + 2 f (xn−1 ) + f (xn )]
Luego agrupando se tiene:

" # Z
n−1 b
h
f (x0 ) + 2 ∑ f (xi ) + f (xn ) ≈ f (x) dx
2 i=1 a
Ejemplo Aplicamos la regla del trapecio para aproximar la integral:

Z 1
2
ex dx
0
solución En este caso, identificamos n = 5, y la partición general del intervalo [0], 1 es

1−0 1
h= = = 0,2.
5 5
Por tanto tenemos los nodos en los siguientes puntos p = {0 , 0.2 , 0.4 , 0.6 , 0.8 , 1} luego
Z 1
2 0.2
ex dx ≈ [ f (0) + 2( f (0.2) + f (0.4) + f (0.6) + f (0.8)) + f (1)]
Z0 2
1 2 2 2 2 2
ex dx ≈ (0.1[1 + 2(e0.2 + e0.4 + e0.6 + e0.8 + e)] = 1.48065
0
Cabe mencionar que el valor verdadero de esta integral es 1.4626 · · · .
Así, vemos con cinco intervalos, la aproximación solo permite una cifra decimal. Para
hacer cálculos con mas subintervalos , es conbeniente elaborar un programa en python
que aplique la fórmula con el numero de subintervalos que uno desee .

Integracion Numerica

Cargado por

Copyright:

Formatos disponibles

Integracion Numerica

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Integracion Numerica

Cargado por

Copyright:

Formatos disponibles

1.

integral no tenga primitivas expresados analíticamente ó que dicha integral no se conoce

f (xi ), dicha suma denominada formula de cuadratura.

los nodos en el eje x de integración o cuadratura

que apliquemos en los criterios.

x1 − x0 = x2 − x1 = x3 − x2 = · · · = xn − Xn−1 = ∆x = h para el método de Gauss-Legendre

Un detalle importante en todas las aplicaciones es la estimación del error cometido y el

conocimiento del orden de precision del criterio de la regla utilizada .

En las reglas de cuadratura es frecuente que estén basados en la interpolación polinómica

en los nodos y su integral es la aproximación buscada.

Si un numero determinado n de nodos en el intervalo [a, b], esto es particionar en n subin-

tervalos .Siempre es posible determinar un polinomio de grado n y generar una formula

si x1 − x0 = x2 − x1 = · · · = xn − xn−1 = ∆x = h osea cada subintervalo de [a, b] tiene la

Las reglas resultantes se conocen como la fórmula cuadratura de Newton.

a continuación nos remitimos al gráfico siguiente

último la curva del la función f .

Si extraemos el i-ésimo trapecio de la figura, se tiene como área:

Luego agrupando se tiene:

Ejemplo Aplicamos la regla del trapecio para aproximar la integral:

solución En este caso, identificamos n = 5, y la partición general del intervalo [0], 1 es

hacer cálculos con mas subintervalos , es conbeniente elaborar un programa en python

que aplique la fórmula con el numero de subintervalos que uno desee .

También podría gustarte