Unidad II. Interpolacion e Integracion Numerica

ESCUELA DE INGENIERIA
METODOS NUMERICOS
UNIDAD II. INTERPOLACION E
INTEGRACION NUMERICA
Dr. Deny González

Contenido
Introducción.
INTERPOLACION POLINOMIAL DE NEWTON EN DIFERENCIAS
DIVIDIDAS.
 INTERPOLACION LINEAL
 INTERPOLACION CUADRATICA
 FORMA GENERAL DE LOS POLINOMIOS DE INTERPOLACION
DE NEWTON.
ESQUEMAS DE INTEGRACION NUMERICA
METODO DE TRAPECIO SIMPLE Y MULTIPLE
METODO DE SIMPSON 1/3 SIMPLE Y MULTIPLE
METODO DE SIMPSON 3/8 SIMPLE Y MULTIPLE
UNIDAD II. INTERPOLACION E INTEGRACION NUMERICA
La Interpolación como función es aquella que pasa a través de puntos dados como datos mostrados en
tablas o tomados de una función dada. Con el propósito de estimar valores intermedios entre datos
definidos por puntos.
Interpolación Lineal:
Es la mas simple interpolación, se basa en la unión de dos puntos por una línea recta.
Formula de la Interpolación Lineal:
f ( x1 )  f ( x0 )
f1 ( x)  f ( x0 )   ( x  x0 )
x1  x0
f1 (X) = Notación de Interpolación de primer grado.
Nota: Mientras menor sea la distancia entre los puntos mejor será la aproximación.
Ejemplo 1 :
Estime el valor del Logaritmo natural de 2.
Seleccionamos los intervalos, ln (2) = 0.6931472
Pto. 1, X = 1 , ln (1) = 0
Pto. 2, X = 3 , ln (3) = 1.098612 f ( x1 )  f ( x0 )
f1 ( x)  f ( x0 )   ( x  x0 )
Pto. 3, X = 5 , ln (5) = 1.609437 x1  x0
Seleccionamos entre el punto 1 y 3.

1.609437 0
f1 ( x)  0   (2  1)  0.402359 Error = 41.95 %
5 1
1.098612 0
f1 ( x)  0   (2  1)  0.549306 Error = 26.18 %
3 1
Y = Ln X
Valor verdadero
Gráficamente =
1 2 3 4 5
Ejemplo 2:
Utilice la tabla de vapor de agua sobrecalentada a 200 Mpa dada a continuación para
encontrar la entropía, s, mediante interpolación lineal, correspondiente al volumen
especifico v = 0.118m3/kg.
Puntos 1 2 3
v (m3/kg) 0.10377 0.11144 0.12547
s (Kj / kg K) 6.4147 6.5453 6.7664
f ( x1 )  f ( x0 )
f1 ( x)  f ( x0 )   ( x  x0 )
x1  x0
6.7664  6.4147
f1 ( x)  6.4147   (0.118  0.10377)  6.6453
0.12547 0.10377
Valor de entropía, s, = 6.6453 kj / kg k
Estime usted, ahora entre los puntos 2 y 3. Analice el resultado.

Interpolación Cuadrática:
Una estrategia para mejorar la estimación en la interpolación consiste en introducir alguna

curvatura a la línea que une los puntos. Si se tienen conocidos tres puntos como datos, estos
pueden ajustarse en un polinomio de segundo grado ( también conocido como polinomio
cuadrático o parábola)
f 2 ( x)  b0  b1( x  x0 )  b2 ( x  x0 )( x  x1)
Donde,
b0  f ( x0 )
f ( x1 )  f ( x0 ) f ( x1 )  b0
b1  
x1  x0 x1  x0
f ( x2 )  f ( x1 ) f ( x1 )  f ( x0 ) f ( x2 )  f ( x1 )
  b1
x2  x1 x1  x0 x2  x1
b2  
x2  x0 x2  x0
Ejemplo 1:
Pto. 1, X = 1 , ln (1) = 0
Pto. 2, X = 3 , ln (3) = 1.098612
Pto. 3, X = 5 , ln (5) = 1.609437
f 2 ( x)  b0  b1( x  x0 )  b2 ( x  x0 )( x  x1)
Donde,
b0  f ( x0 )
f ( x1 )  f ( x0 ) f ( x1 )  b0
b1  
x1  x0 x1  x0
f ( x2 )  f ( x1 ) f ( x1 )  f ( x0 ) f ( x2 )  f ( x1 )
  b1
x2  x1 x1  x0 x2  x1
b2  
x2  x0 x2  x0
x= 2,0000 b0 = 0 b1 = 0,54930614 b2 = -0,07347333
Solucion = 0,622779 Valor real = 0,693142 Error = 10,1512

Ejemplo 2:
Pto. 1, X = 1.5 , ln (1.5) = 0.405465

Pto. 2, X = 1.7 , ln (1.7) = 0.530628
Pto. 3, X = 2.1 , ln (2.1) = 0.741937
f 2 ( x)  b0  b1( x  x0 )  b2 ( x  x0 )( x  x1)
Donde,
b0  f ( x0 )
f ( x1 )  f ( x0 ) f ( x1 )  b0
b1  
x1  x0 x1  x0
f ( x2 )  f ( x1 ) f ( x1 )  f ( x0 ) f ( x2 )  f ( x1 )
  b1
x2  x1 x1  x0 x2  x1
b2  
x2  x0 x2  x0
x= 2,0000 b0 = 0,40546511 b1 = 0,62581571 b2 = -0,16257163
Solucion = 0,693987 Valor real = 0,693142 Error = 0,1219

Ejemplo 2: Mediante programa en python

Ejemplo 3:
Utilice la tabla de vapor de agua sobrecalentada a 200 Mpa dada a continuación para
encontrar la entropía, s, mediante interpolación lineal, correspondiente al volumen
especifico v = 0.118m3/kg.
Puntos 0 1 2
v (m3/kg) 0.10377 0.11144 0.12547
s (Kj / kg K) 6.4147 6.5453 6.7664
Solución,
b0  f ( x0 )  f (0.10377)  6.4147
f ( x1 )  f ( x0 ) 6.5453 6.4147
b1    17.0273
x1  x0 0.11144 0.10377
f ( x2 )  f ( x1 ) f ( x1 )  f ( x0 ) 6.7664  6.5453
  17.0273
x2  x1 x1  x0 0.12547  0.11144
b2    58.4466
x2  x0 0.12547  0.10377
f 2 ( x)  b0  b1( x  x0 )  b2 ( x  x0 )( x  x1)
f 2 ( x)  6.4147  17.0273(0.118  0.10377)  58.4466(0.118  0.10377)(0.118  0.11144)
f 2 ( x)  6.651542
FORMA GENERAL DE LOS POLINOMIOS DE INTERPOLACION DE NEWTON
El polinomio de n-ésimo grado es
Los puntos asociados a con datos se utilizan para evaluar los coeficientes bo, b1,…,bn. Para un
polinomio de n-ésimo grado se requieren n+1 puntos [xo, f(xo)], [x1, f(x1)], …, [xn, f(xn)].
Usamos esos datos y las ecuaciones para evaluar los coeficientes;
Donde las evaluaciones de la función colocadas entre

paréntesis son diferencias divididas finitas
Por ejemplo, la primera diferencia dividida finita en forma general se presenta como
La segunda diferencia dividida finita, que representa la diferencia de las dos primeras
diferencias divididas, se expresa en forma general como
En forma similar la n-ésima diferencia dividida finita es,
El polinomio de interpolación, se conoce como polinomio de interpolación de Newton en

diferencias divididas.
Tabla. Representación Gráfica de la naturaleza de las diferencias divididas finitas

EJEMPLO POLINOMIOS DE INTERPOLACION DE NEWTON
Estime el valor del Logaritmo natural de 2., Seleccionamos los intervalos, ln (2) = 0.6931472
Pto. 0, X0= 1.5 , ln (1.5) = 0.405465 / Pto. 2, X2 = 2.1 , ln (2.1) = 0.741937

Pto. 1, X1= 1.7 , ln (1.7) = 0.530628 / Pto. 3, X3 = 3, ln(3) = 1.098612
Con n = 3, el polinomio de tercer grado es;

EJEMPLO POLINOMIOS DE INTERPOLACION DE NEWTON

Interpolación de Lagrange:
El polinomio de interpolación de lagrange es simplemente una reformulación del polinomio de

Newton que evita el calculo de las diferencias divididas, y se representa de manera concisa como
n
f n ( x)   Li ( x) f ( xi )
Donde,
i 0
n x xj
Li ( x)  
x  xj
j 0 i
j 1
Donde, ∏ designa el “producto de”. Por ejemplo, la versión final (n=1) es
x  x1 x  x0
f1 ( x)  f ( x0 )  f ( x1 )
x0  x1 x1  x0
Y la de segundo grado es,
( x  x1 )( x  x2 ) ( x  x0 )( x  x2 ) ( x  x0 )( x  x1 )
f 2 ( x)  f ( x0 )  f ( x1 )  f ( x2 )
( x0  x1 )( x0  x2 ) ( x1  x0 )( x1  x2 ) ( x2  x0 )( x2  x1 )
Ejemplo 1 (Interpolación de Lagrange)

Pto. 0, X0= 1.5 , ln (1.5) = 0.405465

Pto. 1, X1= 1.7 , ln (1.7) = 0.530628
Pto. 2, X2 = 2.1 , ln (2.1) = 0.741937
Para Primer grado, se toman los ptos 1 y 2

x  x1 x  x0 2  2.1 2  1.7
f1 ( x)  f ( x0 )  f ( x1 ) f1 ( x)  0.530628 0.741937  0.689109
x0  x1 x1  x0 1.7  2.1 2.1  1.7
Para Segundo Grado

( x  x1 )( x  x2 ) ( x  x0 )( x  x2 ) ( x  x0 )( x  x1 )
f 2 ( x)  f ( x0 )  f ( x1 )  f ( x2 )
( x0  x1 )( x0  x2 ) ( x1  x0 )( x1  x2 ) ( x2  x0 )( x2  x1 )
(2  1.7)(2  2.1) (2  1.5)(2  2.1) (2  1.5)(2  1.7)
f 2 ( x)  0.405465 0.530628 0.741937
(1.5  1.7)(1.5  2.1) (1.7  1.5)(1.7  2.1) (2.1  1.5)(2.1  1.7)
f 2 ( x)  0.1013662 0.3316425 0.463710 0.693986

Estime el valor del Logaritmo natural de 2., Seleccionamos los intervalos, ln (2) = 0.6931472
Pto. 0, X0= 1.5 , ln (1.5) = 0.405465

Pto. 1, X1= 1.7 , ln (1.7) = 0.530628
Pto. 2, X2 = 2.1 , ln (2.1) = 0.741937
REGLA DEL TRAPECIO
El método de los trapecios es muy simple y se puede explicar fácilmente a partir

de la siguiente figura.
REGLA DEL TRAPECIO SIMPLE
f (a)  f (b)
I  (b  a) 
2
Eligiendo un espaciado
se divide el intervalo [a, b] por medio de puntos igualmente espaciados
tenemos que, las ordenadas de dichos puntos son

REGLA DEL TRAPECIO
Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el

intervalo de integración de “a” a “b” en varios segmentos, y aplicar el metodo a cada
uno de ellos, según se muestra;
f(x)
f(x)
X0 X1 X2 X3
X0 X1 X2 X0 X1 X2 X3 X4
X0 X1 X2 X3 X4 X5
REGLA DEL TRAPECIO
El área total aproximada es la suma de

las áreas de los n pequeños trapecios de
anchura h
o bien, agrupando términos
Cuanto mayor sea el número de divisiones del intervalo [a, b] que hagamos,
menor será h, y más nos aproximaremos al valor exacto de la integral. Sin
embargo, no podremos disminuir h tanto como queramos, ya que el ordenador
maneja números de precisión limitada.
REGLA DEL TRAPECIO. EJEMPLO

Se tiene la siguiente integral, determine su valor numérico por la regla de trapecio simple
1 𝑓 𝑎 +𝑓 𝑏 0+2
𝐼= 𝑥 2 + 𝑥 𝑑𝑥 𝐼 = 1−0 ∗ = 1 =1
0 2 2
Valor real = 5/6 = 0.833333
Error = 0.16667
REGLA DEL TRAPECIO MULTIPLE

Las áreas de los segmentos se suman después para obtener la integral en todo el
intervalo. Las ecuaciones resultantes se llaman formulas de integración, de
aplicación múltiple o compuestos.
La formula general de la integral se presenta como:
n 1
f ( x0 )  2 f ( xi )  f ( xn )
I  (b  a ) i 1
2n
Se puede calcular el error con la regla del trapecio de aplicación múltiple al sumar
los errores individuales de cada segmento, la ecuación general es:
b  a 3 n
f '' i 
Et 
12n 3 
i 1
Ejemplo de Trapecio Múltiple

1
𝐼= 𝑥 2 + 𝑥 𝑑𝑥 Donde , n = 4
0
n 1
f ( x0 )  2 f ( xi )  f ( xn ) h= (1-0) / 4 = 0.25
I  (b  a ) i 1
2n
f(a) f( a+h) f ( a + 2h ) f ( a + 3h ) f(b)
f(0) f(0,25) f(0,5) f(0,75) f(1)
0 0,3125 0,75 1,3125 2
Valor real = 5/6 = 0.833333

Integral = 0.84375
Error = 0.01042
Ejemplo de Trapecio Múltiple

1
𝐼= 𝑥 2 + 𝑥 𝑑𝑥
0
Ejemplo de Trapecio Múltiple n 1

f ( x0 )  2 f ( xi )  f ( xn )
5 I  (b  a ) i 1
ln xdx
 x 2n
3
Para n = 4
f(a) f ( xo+ h ) f ( xo+ 2h ) f ( xo+ 3h ) f(b)
0,634284 0,66963 0,693147 0,709029 0,719763
f(x)
Integral = 1.3758
1
0,9
0,8
0,7
0,6
0,5 f(x)
0,4
0,3
0,2
0,1
0
3 3,5 4 4,5 5
REGLA DE SIMPSON 1/3
Este método consiste en la aproximación del cálculo del área plana bajo una
curva utilizando trapecios curvilíneos a partir una interpolación con una función
cuadrática:
f ( x0 )  4 f ( x1 )  f ( x2 )
I  (b  a )
6
Designando a y b como X0 y X2, f2(x)
se representa por un polinomio de
Lagrange de segundo grado,
obtendremos la integral
REGLA DE SIMPSON 1/3, MULTIPLE
A modo de interpretación, las sumatorias corresponden a los valores de la función

en los puntos impares y pares respectivamente sin contar los extremos. Esto es
debido a que, en los puntos impares, se encuentra el factor de 4 que introduce la
fórmula de Simpson 1/3 simple para los puntos intermedios y, en los puntos pares,
el factor de 2 se debe a que ese punto es compartido por los trapecios curvilíneos
adyacentes.
𝑛−1 𝑛−2
𝑏−𝑎 𝑓 𝑥0 + 4 𝑖=1,3,5 𝑓 𝑥𝑖 + 2 𝑗=2,4,6 𝑓 𝑥𝑗 + 𝑓 𝑥𝑛
𝐼≈ ∗
𝑛 3
𝑏 − 𝑎 5 (4)
Error → 𝐸𝑎 = − 𝑓
180𝑛4
Ejemplo, Regla de Simpson 1/3, para n =4
4
 x 3 ln xdx
2
f ( li ) f(h) f ( 2h ) f ( 3h ) f ( ls )
5,545177444 14,31704 29,66253 53,71221 88,72284
100
Sum. Par = 59,32506 80
60
Sum. Impar = 272,117 Serie1
40
20
Integral = 70,95168
0
2 2,5 3 3,5 4
REGLA DE SIMPSON 3/8

La derivación de la regla de los tres octavos de simpson es similar a la
regla de un tercio, excepto que se determina el área bajo una parábola
de tercer grado que conecta 4 puntos sobre una curva dada. La forma
general de la parábola de tercer grado es:
I   f ( xo )  3 f ( x1 )  3 f ( x2 )  f ( x3 )
3h
8
ba
donde h 
3
Se considera que a = Xo y b = X3
X1 = a+h
X2 = a+2h
REGLA DE SIMPSON 3/8. MULTIPLE

La formula de simpson 3/8 se puede transformar para utilizar n segmentos, múltiplos de 3
(ya que con tres segmentos queda definida una ecuación cubica). Así, la fórmula 3/8 con
segmentos múltiples es,
𝑏 𝑛−2 𝑛−1
3ℎ
𝐼≈ 𝑓 𝑥 𝑑𝑥 ≈ 𝑓 𝑥0 + 2 𝑓 𝑥𝑖 + 3 𝑓 𝑥𝑗 + 𝑓 𝑥𝑛
8
𝑎 𝑖=3,6,9 𝑗=1,2,4,5
Donde,
𝑥𝑛 − 𝑥0 𝑏 − 𝑎 “N debe ser múltiplo de 3”

ℎ= =
𝑛 𝑛
REGLA DE SIMPSON 3/8. MULTIPLE

METODO DE ROMBERG.
La integración de Romberg es una técnica diseñada para obtener integrales numéricas

de funciones manera eficiente. Se basa en aplicaciones sucesivas de la regla del
trapecio, que mediante manipulaciones matemáticas, se alcanzan mejores resultados.
4 1
𝐼≅ 𝐼𝑚 − 𝐼𝑙
3 3
La ecuación usada para una exactitud O(h6) es
16 1
𝐼≅ 𝐼 − 𝐼
15 𝑚 15 𝑙
Donde Im e Il son las estimaciones mayor y menor, respectivamente. De manera similar,

dos resultados O(h6) se combinan para calcular una integral que es O(h8) utilizando;
64 1
𝐼≅ 𝐼 − 𝐼
63 𝑚 63 𝑙
METODO DE ROMBERG.
Estimar la siguiente integral aplicando el método de Romberg;
1
𝑥 2
𝑒 𝑑𝑥
0 4 1
𝐼 ≅ 𝐼𝑚 − 𝐼𝑙
3 3
16 1
𝐼≅ 𝐼𝑚 − 𝐼𝑙
15 15

Unidad II. Interpolacion e Integracion Numerica

Cargado por

Copyright:

Formatos disponibles

Unidad II. Interpolacion e Integracion Numerica

Cargado por

Información del documento

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

Unidad II. Interpolacion e Integracion Numerica

Cargado por

Copyright:

Formatos disponibles

ESCUELA DE INGENIERIA

Dr. Deny González

Formula de la Interpolación Lineal:

f1 (X) = Notación de Interpolación de primer grado.

Seleccionamos entre el punto 1 y 3.

Seleccionamos entre el punto 1 y 3.

Estime usted, ahora entre los puntos 2 y 3. Analice el resultado.

Una estrategia para mejorar la estimación en la interpolación consiste en introducir alguna

Solucion = 0,622779 Valor real = 0,693142 Error = 10,1512

Pto. 1, X = 1.5 , ln (1.5) = 0.405465

Solucion = 0,693987 Valor real = 0,693142 Error = 0,1219

Ejemplo 2: Mediante programa en python

FORMA GENERAL DE LOS POLINOMIOS DE INTERPOLACION DE NEWTON

El polinomio de n-ésimo grado es

Donde las evaluaciones de la función colocadas entre

FORMA GENERAL DE LOS POLINOMIOS DE INTERPOLACION DE NEWTON

En forma similar la n-ésima diferencia dividida finita es,

El polinomio de interpolación, se conoce como polinomio de interpolación de Newton en

FORMA GENERAL DE LOS POLINOMIOS DE INTERPOLACION DE NEWTON

Tabla. Representación Gráfica de la naturaleza de las diferencias divididas finitas

EJEMPLO POLINOMIOS DE INTERPOLACION DE NEWTON

Pto. 0, X0= 1.5 , ln (1.5) = 0.405465 / Pto. 2, X2 = 2.1 , ln (2.1) = 0.741937

Con n = 3, el polinomio de tercer grado es;

EJEMPLO POLINOMIOS DE INTERPOLACION DE NEWTON

El polinomio de interpolación de lagrange es simplemente una reformulación del polinomio de

Y la de segundo grado es,

Ejemplo 1 (Interpolación de Lagrange)

Pto. 0, X0= 1.5 , ln (1.5) = 0.405465

Para Primer grado, se toman los ptos 1 y 2

Para Segundo Grado

f 2 ( x)  0.1013662 0.3316425 0.463710 0.693986

Pto. 0, X0= 1.5 , ln (1.5) = 0.405465

REGLA DEL TRAPECIO

El método de los trapecios es muy simple y se puede explicar fácilmente a partir

REGLA DEL TRAPECIO SIMPLE

se divide el intervalo [a, b] por medio de puntos igualmente espaciados

tenemos que, las ordenadas de dichos puntos son

REGLA DEL TRAPECIO

Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el

REGLA DEL TRAPECIO

El área total aproximada es la suma de

o bien, agrupando términos

REGLA DEL TRAPECIO. EJEMPLO

Valor real = 5/6 = 0.833333

REGLA DEL TRAPECIO MULTIPLE

Ejemplo de Trapecio Múltiple

Valor real = 5/6 = 0.833333

Ejemplo de Trapecio Múltiple

Ejemplo de Trapecio Múltiple n 1

f(a) f ( xo+ h ) f ( xo+ 2h ) f ( xo+ 3h ) f(b)

0,634284 0,66963 0,693147 0,709029 0,719763

REGLA DE SIMPSON 1/3

REGLA DE SIMPSON 1/3, MULTIPLE

A modo de interpretación, las sumatorias corresponden a los valores de la función

Ejemplo, Regla de Simpson 1/3, para n =4

REGLA DE SIMPSON 3/8

REGLA DE SIMPSON 3/8. MULTIPLE

𝑥𝑛 − 𝑥0 𝑏 − 𝑎 “N debe ser múltiplo de 3”

REGLA DE SIMPSON 3/8. MULTIPLE

La integración de Romberg es una técnica diseñada para obtener integrales numéricas