Interpolaciòn
Podemos definir la interpolciòn como el “arte de leer entre los
valores de una tabla”. Supongamos que se conocen los valores y0, y1,…yn
correspondientes a los (n +1) valores distintos de x0, x1,… xn de la variable
independiente x. La pregunta es ¿qué valor le corresponde a un valor x distinto
de los xn+1 correspondidos entre x0 y xn?
Por los (n+1) puntos dados podemos hacer pasar una curva
representativa de un polinomio de grado n. Así como, por dos puntos pasa una
recta; por tres una parábola, etc., por n+1 puntos pasa una curva de orden n.
Esa curva representativa permitirá asignar a cada x por un valor “interpolante”
de y.
1º Caso: Valores equidistantes
Se cumple la condición:
x0 - x1 = x1 - x2 = x2 - x3 = … = xn - xn+1= h
De esta manera
x1 = x0 + h
x2 = x1 + h = x0 + 2h
…
xn = xn-1+ h = …= x0 + nh
Formaremos una tabla de diferencias de la siguiente manera:
x y ∆y ∆2 y ∆3 y
x0 y0
x1 = x0 + h y1 ∆ y0 = y1 - y0
y2 ∆2 y0 = ∆y1 - ∆y0
x2 = x1 + h y3 ∆ y1 = y2 – y1 ∆3 y0 = ∆2y1 - ∆2y0
. ∆2 y1= ∆y2 - ∆y1
x3 = x2 + h . ∆ y2 = y3 – y2
. .
. yn
.
xn = xn-1 + h
Observamos que
∆ y 0 = y 1 - y0
∆2 y0 = ∆y1 - ∆y0 = y2 - 2y1 + y0
∆3 y0 = ∆2y1 - ∆2y0 = y3 - 3y2 + 3y1 - y0
-
-
-
Las diferencias n-simas de y0 se forman a partir de la ordenada n-sima y
sus n ordenadas anteriores afectadas por los coeficientes del desarrollo del
binomio de Newton: (a - b)n .
Fórmula de Gregory-Newton:
El polinomio propuesto por los autores del método tiene la forma:
Pn(x) = a0 + a1(x - x0) + a2 (x - x0)(x - x1) + …+ an(x - x0)(x - x1) … (x - xn-1)
Esta fórmula es un polinomio de grado n que pasa por los (n+1) valores o
puntos dados como datos.
Para determinar los coeficientes a0…an imponemos el valor y0 al
polinomio cuando x sea igual a x0.
Pn(x0)= a0 + a1.0 +a2.0 + … + an.0 = y0
De modo que a0 = y0
Tomamos el par (x1; y1) y reemplazamos en el polinomio
Pn(x1) = a0 + a1(x1 - x0) + a2.0 + … + an.0 = y1
Del cual despejamos:
a1 = = =
Hacemos lo mismo para a2 reemplazando x e y por el par (x2; y2)
Pn(x2) = y2 = a0 + a1(x2 - x0) + a2 (x2 - x0) (x2 - x1)
a2 = = = =
Deducimos por analogía que el coeficiente a3 será
a3 =
y por lo tanto el coeficiente n-ésimo tendrá el valor:
an =
Ejemplo: Formemos la tabla de diferencias para los siguientes valores:
x y ∆y ∆2 y ∆3 y
0 0 - - -
1 1 1 - -
2 4 3 2 -
3 9 5 2 0
4 16 7 2 0
Vemos que a partir de ∆3 las diferencias son nulas.
Esto es así porque la función que hemos tomado para el ejemplo es
y = x2
Continuemos:
Pn(x0) = y0 Pn (0) = 0 = a0 ; a0 = 0
Pn(x1) = y1 a1= =
Pn(x2)= y2 a2=
Pn(x) = a0 + a1 (x - x0) + a2 (x - x0) (x - x1)
= 0 +1 ( x - 0) +1(x - 0) (x -1)
= x + x2 - x = x2
Pn(x) = x2 , que es la función buscada.
Cuando la función que queremos aproximar es una función trascendente
como ex; ln x; sen x; etc. Las diferencias terceras, cuartas, quintas, etc., no se
anulan pero sus valores podrán hacerse despreciables de acuerdo al grado de
exactitud que se desee obtener.
Formula de Lagrange:
Aplicable a intervalos no equidistantes.
Tiene la forma:
Pn(x) = a0( x -x1)(x - x2)…(x - xn) + a1 (x - x0)(x - x2)…(x – xn)+…
+ an( x –x0)(x – x1)…(x – xn-1) = y
Imponiendo a esta fórmula los (n +1) condiciones o pares de valores
(x ; y) obtendremos (n+ 1) coeficientes. Así se obtiene:
y0 = a0 (x0 - x1)…(x0 - xn)
y1 = a1 (x1 - x0)…(x1 - xn)
...
yn = an (xn - x0)…(xn - xn-1)
La fórmula de Lagrange quedará finalmente:
Pn (x) =
Ejemplo:
En un ensayo se obtiene la siguiente tabla de valores:
x = 0 ; 490 ; 840 ; 1003
y = 0 ; 3150 ; 5036 ; 5773
¿Qué valor corresponde a x = 700?
Calculemos primero los coeficientes a0 : a1 ; a2 ; y a3
a0=
a1=
a2 =
a3 =
Finalmente:
P3(x) =3,58 x 10-5 x (x-840) (x-1003) -10,51 x 10-5 x(x-490) (x-1003) + 6,88 x
10-5 x(x-490) (x-840) =
Luego para x = 700
P3(700) = 4328
Ejemplo: Calcular sen 0,42 teniendo como dato
Sen 0,3 = 0,295520
Sen 0,4 = 0,389418
Sen 0,5 = 0,479426
Sen 0,6 = 0,564642
Usando los coeficiente lagrangianos dados como ejemplo y por ser x0 =0,4
x0 + hµ = 0,5 hµ = 0,1 µ= 0,1
El punto buscado es x = 0,42 x0 + hµ = 0,42
reemp: 0,4 + h . 0,1 = 0,42
despejo h: h= = 0,2
Como h = 0,2 es el caso del ejemplo podemos usar esos coeficientes:
P(0,42) = - 0,048 . (sen 0,3) + 0,864 (sen 0,4) + 0,216 (sen 0,5)
- 0,032 (sen 0,6)
P(0,42) = 0,407760
La interpolaciòn tiene 6 cifras exactas. Si aplicáramos interpolación lineal
tendría solo 3 cifras exactas.
Interpolación Parabólico Progresiva:
Esta fórmula es aplicable a intervalos no equidistantes. Su
deducción se debe a Newton y el polinomio que se obtiene tiene la forma:
Pn (x) = a0 + a1(x-x0) + a2(x-x0)(x-x1) + ….
+ an(x-x0)(x-x1)…(x-xn-1)
Imponiendo luego la condición de pasar por los (n+1) puntos
conocidos x0 ··· Pn. Tomemos primero para x = x0 ; y = y0
P0(x0)= a0 = y0
y la ecuación obtenida es una recta paralela al eje x.
La aproximación es muy mala.
Si imponemos la siguiente condición x = x1 ; y = y1
P1(x1) = y0 +a1(x1 – x0)
a1 =
Y el polinomio será una recta de pendiente a1 que pasa por (x0 : y0) y por
(x1 ; y1).
Si ahora imponemos un punto más obtenemos una parábola de
segundo orden
P2(x2)= a0 + a1 (x2-x0) + a2 (x2-x0)(x2-x1) = y2
Donde el coeficiente será
a2 =
Vemos que el numerador de a2 es el error que se cometería si se
tomara como valor de y2 el valor calculado por la aproximación de primer orden
para el punto y2.
a2 =
Así para n + 1 puntos
Pn(x) =a0 +a1 (x-x0)+…..+ an (x-x0)…(x-xn-1)
an =
Ejemplo:
Se ha determinado la presión del vapor de agua para diversos
valores de la temperatura. Se desea calcular la presión para t = 10 ºC
a0 = P0 = 63
t P a1 =
0 ºC 63
P1 = 63 + 5,2 t P1(t2) = 141
5 ºC 89
15 ºC 173
30 ºC 429 a2 =
P2 = 63 + 5,2 t + 0,21 t (t – 5) P2(t3) = 376,5
a3 =
El polinomio de interpolación de 3º grado será:
P3(t) = 63 +5,2 t + 0,21 t(t -5) + 0,0047 t(t -5)(t -15)
Reemplazamos t = 10 ºC
P3 (10 ºC) = 124
Aproximación de Funciones
Método de los Mínimos Cuadrados
El método consiste en buscar un conjunto de coeficientes a0 …. am
del polinomio de manera que se haga mínimo el error de aproximación del
mismo a la función dato.
La cantidad de puntos datos del problema será n.
Se debe cumplir que m < n.
Veamos el siguiente grafico:
x1 f1 f3 P1, P2, P3 son
x2 f2 puntos del polinomio
x3 f3 f1 f1 P2 3 que aproxima
… … P3 P(x)
fi
xn fn 1 2
P1 f2
+
hf h h2
i
x1 x2 x3
ξi = ƒi – P(xi) = Error de orden i [1]
El polinomio tendrá la forma
P(x) = a0 + a1x + a2 x2 + …. + am xm
Que también se puede escribir
P(x) =
Entonces: reemplazado en [1]
ξi =ƒi -
La suma de los cuadrados de los n puntos será la función Error cuadrático
total:
ƒ(a0 ; a1 ;….am) = (ƒi - )2
Derivando esta función respecto de cada uno de los ak coeficientes
obtendremos m+1 ecuaciones que igualadas a cero y resueltas en forma
simultanea nos darán las m+1 soluciones a0 …am para formar el polinomio
buscado:
[ (ƒi - )2] = 0
[ (ƒi - )] = 0
ƒi - =0
ƒi = ecuación Nº K del sistema
Hacemos k = 0
=
Hacemos k = 1
=
Finalmente llegamos a k = m
=
Obtenemos así las m+ 1 ecuaciones que resueltas en forma simultánea
nos permitirán calcular las incógnitas a0, a1 . . . am.
Ejemplo:
Dados los Valores de la tabla construir un polinomio interpolante
de grado m = 3 y n = 5 .
x f
0 -2
1 0,04
1,5 2,91
3 -1,6
4,5 1,1
Calculamos:
Σ xi =10 Σ fi = 0,45
Σ xi2= 32,5 Σ fixi = 4,555
Σ xi3= 122,5 Σ fixi2 = 14,4625
Σ xi4= 497,125 Σ fixi3 = 66,8987
Σ xi5= 2096,87
Σ xi6= 9045,156
Formamos el sistema de ecuaciones:
5 a0 + 10 a1 + 32,5 a2 + 122,5 a3 = 0,45
10 a0 + 32,5 a1 + 122,5 a2 + 497,125 a3 = 4,555
32,5 a0 +122,5 a1 + 497,125 a2 + 2096,87 a3 = 14,4625
122,5 a0 +497,125 a1 + 2096,87 a2 + 9045,156 a3 = 66,8987
Resulto el sistema tenemos:
a0= -2,246447
a1= 7,067418
a2= -3,9111967
a3= 0,556097
Finalmente el polinomio que aproxima a la función será:
P3(x) = -2,246447 + 7,067418 x + (-3,9111967) x2 + 0,556097 x3
