Integracion Numerica

Integracin Numrica
JULIO FUENTEALBA V.
MS. CS. U. DE CHILE
DR. U.DE CHILE
Integracin Numrica
La integracin numrica es una herramienta en la

ingeniera para obtener valores aproximados de
integrales definidas que no pueden calcularse
analticamente.
Integracin Numrica
Matemticamente la integracin se representa

por:
Ec 1
que se tiene para la integral de la funcin f(x) con

respecto a la variable independiente x, evaluada en
los lmites x=a y x=b
El significado de la Ec 1 es el valor total o
sumatoria de f(x) sobre el rango x=a a x=b.
Integracin Numrica
el proceso representado en la Ec 1 y en la Fig 1 es

llamado integracin definida
a b
Fig 1
Formula de Newton Cotes
Las frmulas de integracin de Newton-Cotes son los

esquemas de integracin numrica ms comunes.
Se basan en la estrategia de reemplazar una funcin
complicada o datos tabulados con una funcin
aproximada que sea fcil de integrar:
b b
I f ( x)dx f n ( x)dx Ec 2
a a
donde fn(x) es igual a un polinomio de la forma:

Ec 3
donde n es el orden del polinomio.

Por ejemplo, en la Fig. 2 se usa el polinomio de

primer orden (una lnea recta) como una
aproximacin. Mientras que en la Fig. 3 se emplea
una parbola para el mismo propsito.
Fig 2 Fig 3
Por ejemplo, en la Fig. 4

se usan tres segmentos
de lnea recta para
aproximar la integral.
Pueden utilizarse
polinomios de orden
superior para los mismos
propsitos.
Fig 4
La regla del Trapecio o Trapezoidal
La Regla Trapezoidal es la primera de las frmulas

de integracin cerrada de Newton-Cotes.
Corresponde al caso donde el polinomio de la Ec 2 es
de primer orden.
b f (a ) f (b)
I a
f ( x)dx (b a )
2

Ec 2
en MatLab
syms x
f = input('ingrese formula =' );
a = input('integrar desde =' );
b = input('hasta =' );
imagenDeA = subs(f,x,a);
imagenDeB = subs(f,x,b);
trapecio = (b-a)*((imagenDeA + imagenDeB)/2);
disp(trapecio);
Geomtricamente, la regla del trapecio es

equivalente a aproximar el rea del trapecio bajo la
lnea recta que conecta a f(a) y f(b) como se muestra
en Fig. 7.
Recuerde, que la frmula para calcular el rea de un

trapezoide es la altura por el promedio de las bases,
tal y como se muestra en la figura.
Al recostar el trapecio la formula es la misma.
Ejercicios
Use la Regla del Trapecio para aproximar los valores

de las siguientes integrales:
Regla del Trapecio se puede ampliar si subdividimos

el intervalo [a,b] en n subintervalos, todos de la
misma longitud.
Usando las propiedades de la integral, tenemos que
Aplicando la Regla del Trapecio a cada una de las

integrales, obtenemos :
En MAtLab
syms x
%trapecio con n intervalos
n = input('numero d intervalos para trapacio =' );
xi = (b - a)/n;
resultado = 0 ;
primero = a ;
for i=1 : n
segundo = primero+xi;
disp(segundo);
ImagPrim = subs(f,x,primero);
ImagSegua = subs(f,x,segundo);
resultado = resultado + (segundo-primero)*((ImagPrim +ImagSegua)/2);
primero = primero + xi;
end
disp(resultado);
Aplicando la Ahora bien, ya que los subintervalos

tienen la misma longitud h, tenemos que :
Sustituyendo el valor de h y haciendo uso de la

notacin de sumatoria, tenemos finalmente:
Ec. 6
Esta es la regla del trapecio para n subintervalos. Entre ms subintervalos

usemos, mejor ser la aproximacin a la integral.
Aplicacin mltiple de la R. del T.
Ilustracin de la Regla
Trapezoidal de aplicacin
mltiple:
a) dos segmentos
b) tres segmentos
Aplicacin mltiple de la R. del T.
Ilustracin de la Regla
Trapezoidal de aplicacin
mltiple:
a) cuatro segmentos
b) cinco segmentos
Ejercicios
Use la Regla del Trapecio para aproximar los valores

de las siguientes integrales:
Si subdividimos en 5 partes
Regla de Simpson
Adems de aplicar la regla trapezoidal con

segmentacin ms fina, otra forma de obtener una
estimacin ms exacta de la integral es con el uso de
polinomios de orden superior para conectar los
puntos.
Regla de Simpson
Por ejemplo, si hay un punto extra a la mitad del

camino f(a) y f(b), los tres puntos se pueden conectar
en una parbola, tal y como se muestra en la Fig. 12.
Fig. 12
Regla de Simpson
Si hay dos puntos igualmente espaciados entre f(a) y

f(b), los cuatro puntos se pueden conectar con un
polinomio de tercer orden, tal y como se muestra en
la Fig. 13.
Fig. 13
Regla de Simpson
Las frmulas que resultan al tomar las integrales

bajo estos polinomios son conocidas como Regla de
Simpson.
La Regla de Simpson 1/3 resuelta cuando una
interpolacin polinomial de segundo orden es
sustituida en la ecuacin:
Ec. 7
Regla de Simpson
Si a y b se designan como xo y x2 y f2(x) es

representada por un polinomio de Lagrange de
segundo orden y la integral se transforma en:
( x x1 )( x x2 ) ( x x0 )( x x2 ) ( x x0 )( x x1 )
x2
I f ( x0 ) f ( x1 ) f ( x2 )dx
x0
( x0 x1 )( x0 x2 ) ( x1 x0 )( x1 x2 ) ( x2 x0 )( x2 x1 )
Despus de la integracin y manejo algebraico,

resulta la siguiente frmula:
(b a )
I f ( x0 ) 4 f ( x1 ) f ( x2 ) f ( x0 ) 4 f ( x1 ) f ( x2 )
h 2
3 3
Regla de Simpson
Recuerde que x1 es el punto medio entre a y b.
f ( xo ) 4 f ( x1 ) f ( x2 )
I (b a) Ec. 8
6
Esta ecuacin es conocida como Regla de Simpson

de 1/3. Es la segunda frmula de integracin cerrada
de Newton-Cotes.
La especificacin 1/3 surge del hecho de que h est
dividida entre 3 en la ecuacin anterior.
syms x
disp('USO DE LA REGLA DE SIMPSON 3/8')
% esta es la regla de simspson 3/8
valordelaintegral = (b-a)*((subs(f,x,a)+
3*subs(f,x,a+(b-a)/3)+3*subs(f,x,a+2*(b-a)/3)+
subs(f,x,b))/8);
% muestra el resultado
disp(valordelaintegral);
syms x
b = input('hasta=);
%ahora se hace un grafico de la funcion uando 100 puntos
incremento = (b-a)/100 ;
A=(1:100);
B=(1:100);
for i=1 : 100
A(i) = a ;
B(i)= subs(f,x,a);
a = a + incremento ;
end
plot(A,B);
Ejercicios
Use la Regla de Simpson de 1/3 para aproximar el

valor de las siguientes integrales:
Regla de Simpson
La Regla de Simpson se puede ampliar si

subdividimos el intervalo [a,b] en n intervalos de
igual longitud. h= (b-a)/n
. Sea P={x0 , x1, x2 ., xn} la particin que se forma al
hacer dicha subdivisin y sea Pm={ p1, p2 ., pn} el
conjunto de los untos medios de los subintervalos.
Usando las propiedades de la integral, tenemos que:
Regla de Simpson
Al sustituir la Regla de Simpson de 1/3 a cada una de

las integrales, obtenemos:
f ( xo ) 4 f ( x1 ) f ( x2 ) f ( x1 ) 4 f ( x2 ) f ( x3 ) f ( xn 2 ) 4 f ( xn 1 ) f ( xn )
I 2h 2h ... 2h
6 6 6
Regla de Simpson
Al Combinando trminos y
sustituyendo nos
queda f ( x ) 4 n f ( x ) 2 n1 f ( x ) f ( x )
o m i 2
I (b a) i 1 i 1 Ec. 9
6n
Representacin grfica de la
Regla de Simpson 1/3 de
aplicacin mltiple. Observe
que el mtodo se puede
emplear slo si el nmero de
segmentos es par
Ejercicios
En los ejercicios 1 a 3, use (a) la Regla del Trapecio y (b) la

Regla de Simpson, con el valor de n indicado para estimar
las integrales definidas.
Aplique valores aproximados de f (xk) que tengan una
precisin de cuatro decimales y redondee las respuestas a
dos decimales.
Ninguna de las integrales definidas de los ejercicios 4 a 6
puede ser evaluada exactamente en trminos de funciones
elementales. Utilice la Regla de Simpson, con el valor de n
que se indica, para determinar un valor aproximado de la
integral definida dada. Exprese el resultado con tres cifras
decimales.
Ejercicios
Ejercicios
Solucin 1
Solucin 2
Solucin 3
Regla de Simpson 3/8
En una manera similar a la derivacin de la Regla

Trapezoidal y Regla de Simpson 1/3, un polinomio
de Lagrange de tercer orden se puede ajustar a
cuatro puntos e integrarse:
Ec. 10
Para obtener:
Donde h = (b-a)/3 . Esta ecuacin se llama Regla

e Simpson de 3/8 debido a que h se multiplica por
3/8.
Notese que x1 y x2 son los puntos que dividen en tres
partes iguales el intervalo [a,b]
sta es la tercera frmula de integracin cerrada de

Newton-Cotes. La Regla de Simpson 3/8 se puede
expresar tambin de la forma:
Ec. 11
Ilustracin de cmo se puede

usar en conjunto las Reglas
de Simpson de 1/3 y 3/8 para
menejar aplicaciones
mltiples con nmeros nones
de intervalos.
Al igual que en los casos anteriores, la Regla de

Simpson de 3/8 se puede extender si subdividimos
el intervalo [a.b] en n intervalos de la misma
longitud h.
Sea la particin determinada de esta forma. Cada subintervalo
lo dividimos en tres partes iguales, y sean y
los puntos determinados as:

Aplicado la Regla de Simpson de 3/8 en cada uno

de los intervalos, tenemos
ba n n 1

b
f ( x)dx f ( xo ) 3 f ( yi ) f ( z i ) 2 f ( xi ) f ( xn )
a
8n i 1 i 1
Taller
CALCULAR EL VALOR DE LA INTEGRAL:
HACIENDO USO DE:

1. REGLA DEL TRAPECIO SIMPLE
2. REGLA DEL TRAPACIO COMPUESTO EN n=3
3. REGLA DE SIMPSON DE 1/3 SIMPLE
4. REGLA DE SIMPSON DE 1/3 COMPUESTO CON n=3
5. REGLA DE SIMPSON DE 3/8 SIMPLE
. no me pilla huichipirichi

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

Integracin Numrica

La integracin numrica es una herramienta en la

Matemticamente la integracin se representa

que se tiene para la integral de la funcin f(x) con

el proceso representado en la Ec 1 y en la Fig 1 es

Las frmulas de integracin de Newton-Cotes son los

donde fn(x) es igual a un polinomio de la forma:

donde n es el orden del polinomio.

Por ejemplo, en la Fig. 2 se usa el polinomio de

Por ejemplo, en la Fig. 4

La Regla Trapezoidal es la primera de las frmulas

Geomtricamente, la regla del trapecio es

Recuerde, que la frmula para calcular el rea de un

Use la Regla del Trapecio para aproximar los valores

Regla del Trapecio se puede ampliar si subdividimos

Aplicando la Regla del Trapecio a cada una de las

Aplicando la Ahora bien, ya que los subintervalos

Sustituyendo el valor de h y haciendo uso de la

Esta es la regla del trapecio para n subintervalos. Entre ms subintervalos

Use la Regla del Trapecio para aproximar los valores

Adems de aplicar la regla trapezoidal con

Por ejemplo, si hay un punto extra a la mitad del

Si hay dos puntos igualmente espaciados entre f(a) y

Las frmulas que resultan al tomar las integrales

Si a y b se designan como xo y x2 y f2(x) es

Despus de la integracin y manejo algebraico,

Recuerde que x1 es el punto medio entre a y b.

Esta ecuacin es conocida como Regla de Simpson

Use la Regla de Simpson de 1/3 para aproximar el

La Regla de Simpson se puede ampliar si

Al sustituir la Regla de Simpson de 1/3 a cada una de

En los ejercicios 1 a 3, use (a) la Regla del Trapecio y (b) la

En una manera similar a la derivacin de la Regla

Donde h = (b-a)/3 . Esta ecuacin se llama Regla

sta es la tercera frmula de integracin cerrada de

Ilustracin de cmo se puede

Al igual que en los casos anteriores, la Regla de

los puntos determinados as:

Aplicado la Regla de Simpson de 3/8 en cada uno

CALCULAR EL VALOR DE LA INTEGRAL:

HACIENDO USO DE:

También podría gustarte