Interpolación y Ajuste de curvas

Aproximación Polinomial Simple:

La interpolación es de gran importancia en el campo de la ingeniería ya que, al consultar fuentes
de información presentadas en forma tabular, es frecuente no encontrar el valor buscado como
un punto en la tabla.

En general, sí se desea aproximar una función con un polinomio de grado ¨n¨, se necesitan n+1
puntos que sustituidos en el polinomio general de grado ¨n¨ de la forma:

𝑃𝑛 (𝑥 ) = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 … + 𝑎𝑛 𝑥 𝑛

Se genera un sistema de n+1 ecuaciones lineales con incógnitas a0, a1,a2,…an. Una vez resuelto
el sistema de n+1 en el polinomio propuesto, con la cual se obtiene un polinomio de
aproximación que se utilizará para realizar la interpolación en un punto especifico.

Polinomios de Interpolación de Lagrange

Existen métodos de aproximación polinomial en los cuales no siempre se requiere resolver un
sistema de ecuaciones lineales, sino que los cálculos se realizan a través de fórmulas como este
método.
Para obtener el polinomio de interpolación de primer grado, se propone el polinomio de
aproximación de primer grado que de acuerdo a Lagrange tiene la siguiente forma:

𝑃(𝑥 ) = 𝑎0 (𝑥 − 𝑥1 ) + 𝑎1 (𝑥 − 𝑥0 )

Para obtener el valor de a0 se hace que 𝑥 = 𝑥0

𝑓(𝑥0 ) = 𝑎0 (𝑥0 − 𝑥1 )

𝑓(𝑥0 )
𝑎0 =
(𝑥0 − 𝑥1 )
Para obtener el valor de a1 se hace que 𝑥 = 𝑥1

𝑓(𝑥1 ) = 𝑎1 (𝑥1 − 𝑥0 )

𝑓(𝑥1 )
𝑎1 =
(𝑥1 − 𝑥0 )

Sustituyendo a a0 y a1 en el polinomio propuesto, se tiene:

𝑓 (𝑥0 ) 𝑓 (𝑥1 )
𝑓 (𝑥 ) = (𝑥 − 𝑥1 ) + (𝑥 − 𝑥0 )
(𝑥0 − 𝑥1 ) (𝑥1 − 𝑥0 )

Reacomodando términos y dejando el polinomio de Lagrange de primer grado:

(𝑥 − 𝑥1 ) (𝑥 − 𝑥0 )
𝑓 (𝑥 ) = 𝑓 (𝑥0 ) + 𝑓(𝑥1 )
(𝑥0 − 𝑥1 ) (𝑥1 − 𝑥0 )

Polinomio de Lagrange de segundo grado:

𝑓(𝑥0 )(𝑥 − 𝑥1 ) (𝑥 − 𝑥2 ) (𝑥 − 𝑥0 ) (𝑥 − 𝑥2 ) (𝑥 − 𝑥0 ) (𝑥 − 𝑥1 )
𝑓(𝑥) = + 𝑓(𝑥1 ) + 𝑓(𝑥2 )
(𝑥0 − 𝑥1 ) (𝑥0 − 𝑥2 ) (𝑥1 − 𝑥0 ) (𝑥1 − 𝑥2 ) (𝑥2 − 𝑥0 ) (𝑥2 − 𝑥1 )

Polinomio de Langrange de tercer grado:

𝑓(𝑥0 )(𝑥 − 𝑥1 ) (𝑥 − 𝑥2 ) (𝑥 − 𝑥3 ) (𝑥 − 𝑥0 ) (𝑥 − 𝑥2 ) (𝑥 − 𝑥3 )
𝑓(𝑥) = + 𝑓(𝑥1 )
(𝑥0 − 𝑥1 ) (𝑥0 − 𝑥2 ) (𝑥0 − 𝑥3 ) (𝑥1 − 𝑥0 ) (𝑥1 − 𝑥2 ) (𝑥1 − 𝑥3 )
(𝑥 − 𝑥0 ) (𝑥 − 𝑥1 ) (𝑥 − 𝑥3 ) (𝑥 − 𝑥0 ) (𝑥 − 𝑥1 ) (𝑥 − 𝑥2 )
+ 𝑓(𝑥2 ) + 𝑓(𝑥3 )
(𝑥2 − 𝑥0 ) (𝑥2 − 𝑥1 ) (𝑥2 − 𝑥3 ) (𝑥3 − 𝑥0 ) (𝑥3 − 𝑥1 ) (𝑥3 − 𝑥2 )
 MINIMOS CUADRADOS
(Ajuste de curva)

A veces la información dada en tablas tiene errores de significativos sobretodo cuando provienen
de medidas experimentales. En este caso no tiene sentido pasar un polinomio de aproximación
por los puntos dados sino solo cerca de ellos.

En este método el grado del polinomio no tiene relación con el numero de puntos utilizados.

Ajuste lineal.

Proponiendo el polinomio de primer grado de la forma:

𝑦 = 𝑎0 + 𝑎1 𝑥

∑ 𝑦 = 𝑛𝑎0 + 𝑎1 ∑ 𝑥

∑ 𝑥𝑦 = 𝑎0 ∑ 𝑥 + 𝑎1 ∑ 𝑥 2

Ajuste polinomial.

De la misma manera como se genera el sistema de ecuaciones para el polinomio de primer grado
puede obtenerse para cualquier polinomio de grado ¨n¨ de la forma:

𝑦 = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + ⋯ 𝑎𝑛 𝑥 𝑛

Sistema general de ecuaciones para un polinomio ¨n¨.

∑ 𝑦 = 𝑛𝑎0 + 𝑎1 ∑ 𝑥 + 𝑎2 ∑ 𝑥 2 + ⋯ 𝑎𝑛 ∑ 𝑥 𝑛

∑ 𝑥𝑦 = 𝑎0 ∑ 𝑥 + 𝑎1 ∑ 𝑥 2 + 𝑎2 ∑ 𝑥 3 + ⋯ 𝑎𝑛 ∑ 𝑥 𝑛+1

∑ 𝑥 2 𝑦 = 𝑎0 ∑ 𝑥 2 + 𝑎1 ∑ 𝑥 3 + 𝑎2 ∑ 𝑥 4 + ⋯ 𝑎𝑛 ∑ 𝑥 𝑛+2
…
…
∑ 𝑥 𝑛 𝑦 = 𝑎0 ∑ 𝑥 𝑛 + 𝑎1 ∑ 𝑥 𝑛+1 + 𝑎2 ∑ 𝑥 𝑛+2 + ⋯ 𝑎𝑛 ∑ 𝑥 2𝑛
Ajuste no polinamial.

Se puede tener cualquier función y el sistema de ecuaciones se puede generar de la misma

manera que para una función polinomial.

Las funciones típicas son:

- Logaritmica (Log): 𝑦 = 𝑎0 + 𝑎1 ln 𝑥

- Exponencial (Exp): 𝑦 = 𝑎0 𝑒 𝑎1 𝑥 ó ln 𝑦 = ln 𝑎0 + 𝑎1 𝑥

- Potencia (pwr): 𝑦 = 𝑎0 𝑥 𝑎1 ó ln 𝑦 = ln 𝑎0 + 𝑎1 ln 𝑥

1
- Inversa (Inv): 𝑦 = 𝑎0 + 𝑎1 𝑥

Otra forma de ajuste:

𝑦 = 𝑎0 𝑥 + 𝑎1 sin 𝑥

∑ 𝑥𝑦 = 𝑎0 ∑ 𝑥 2 + 𝑎1 ∑ 𝑥 sin 𝑥

∑ 𝑦 sin 𝑥 = 𝑎0 ∑ 𝑥 sin 𝑥 + 𝑎1 ∑ sin^2( 𝑥)

 Regresión lineal múltiple por mínimos cuadrados

𝑦 = 𝑎0 + 𝑎1 𝑥1 + 𝑎2 𝑥2

𝑛 ∑ 𝑥1𝑖 ∑ 𝑥2𝑖 𝑎0 ∑ 𝑦𝑖
∑ 𝑥1𝑖 ∑ 𝑥 2 1𝑖 ∑ 𝑥1𝑖 𝑥2𝑖 𝑎1 = ∑ 𝑥1𝑖 𝑦𝑖
∑ 𝑥2𝑖 ∑ 𝑥1𝑖 𝑥2𝑖 ∑ 𝑥 2 2𝑖 𝑎2 ∑ 𝑥2𝑖 𝑦𝑖

Diferencias Divididas
𝑓(𝑥)−𝑓(𝑥0 )
Por definición la derivación en el punto X0 de una función analítica f(x) es 𝑓 ′ (𝑥) = , sin
(𝑥−𝑥0 )
embargo, cuando la función esta en forma tabular las derivadas solo pueden obtenerse
aproximadamente como sigue:

Para la tabla siguiente

Puntos X F(x)
0 𝑥0 𝑓(𝑥0 )
1 𝑥1 𝑓(𝑥1 )

2 𝑥2 𝑓(𝑥2 )
3 𝑥3 𝑓(𝑥3 )
4 𝑥4 𝑓(𝑥4 )
5 𝑥5 𝑓(𝑥5 )
Las diferencias divididas quedan como:

Primera

𝑓 (𝑥1 ) − 𝑓(𝑥0 )
𝑓[𝑥0 , 𝑥1 ] =
𝑥1 − 𝑥0

𝑓 (𝑥2 ) − 𝑓(𝑥1 )
𝑓[𝑥1 , 𝑥2 ] =
𝑥2 − 𝑥1

𝑓 (𝑥3 ) − 𝑓(𝑥2 )
𝑓[𝑥2 , 𝑥3 ] =
𝑥3 − 𝑥2

𝑓 (𝑥4 ) − 𝑓(𝑥3 )
𝑓[𝑥3 , 𝑥4 ] =
𝑥4 − 𝑥3

𝑓 (𝑥5 ) − 𝑓(𝑥4 )
𝑓[𝑥4 , 𝑥5 ] =
𝑥5 − 𝑥4

Segunda

𝑓[ 𝑥1 , 𝑥2 ] − 𝑓[𝑥0 , 𝑥1 ]
𝑓[ 𝑥0 , 𝑥1 , 𝑥2 ] =
𝑥2 − 𝑥0

𝑓[ 𝑥2 , 𝑥3 ] − 𝑓[𝑥1 , 𝑥2 ]
𝑓[ 𝑥1 , 𝑥2 , 𝑥3 ] =
𝑥3 − 𝑥1

𝑓[ 𝑥3 , 𝑥4 ] − 𝑓[𝑥2 , 𝑥3 ]
𝑓[ 𝑥2 , 𝑥3 , 𝑥4 ] =
𝑥4 − 𝑥2

𝑓[ 𝑥4 , 𝑥5 ] − 𝑓[𝑥3 , 𝑥4 ]
𝑓[ 𝑥3 , 𝑥4 , 𝑥5 ] =
𝑥5 − 𝑥3
Tercera

𝑓[ 𝑥1 , 𝑥2 , 𝑥3 ] − 𝑓[𝑥0 , 𝑥1 , 𝑥2 ]
𝑓[ 𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 ] =
𝑥3 − 𝑥0

𝑓[ 𝑥2 , 𝑥3 , 𝑥4 ] − 𝑓[𝑥1 , 𝑥2 , 𝑥3 ]
𝑓[ 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 ] =
𝑥4 − 𝑥1

𝑓[ 𝑥3 , 𝑥4 , 𝑥5 ] − 𝑓[𝑥2 , 𝑥3 , 𝑥4 ]
𝑓[ 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ] =
𝑥5 − 𝑥2

Cuarta

𝑓[ 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 ] − 𝑓[𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 ]
𝑓[ 𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 ] =
𝑥4 − 𝑥0

𝑓[ 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ] − 𝑓[𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 ]
𝑓[ 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ] =
𝑥5 − 𝑥1

Quinta

𝑓[ 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ] − 𝑓[𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 ]
𝑓[ 𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 , 𝑥4 , 𝑥5 ] =
𝑥5 − 𝑥0
 Polinomio de Newton en diferencias divididas
(Aproximación polinomial de Newton)
Sí se desea aproximar una función dada en forma tabular con un polinomio de grado ¨n¨, se
puede establecer en general que este polinomio puede escribirse de la forma:

𝑃𝑛 (𝑥) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) + ⋯ . 𝑎𝑛 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥𝑛 )

Donde los coeficientes a0,a1, a2, an, están dados por:

𝑎0 = 𝑓(𝑥0 )

𝑎1 = 𝑓[𝑥0 , 𝑥1 ]

𝑎2 = 𝑓[𝑥0 , 𝑥1 , 𝑥2 ]
.
.
.

𝑎2 = 𝑓[𝑥0 , 𝑥1 , 𝑥2 , 𝑥𝑛 ]

Polinomio de Newton en Diferencias Finitas

Cuando la distancia ¨h¨ entre dos argumentos consecutivos cualesquiera es la misma a lo largo
de la tabla, el polinomio de Newton en diferencias divididas puede expresarse con mas sencillez.

Definiendo las diferencias finitas hacia adelante:

Primera

∆f(𝑥𝑖 ) = 𝑓(𝑥𝑖+1 ) − 𝑓(𝑥𝑖 )

Segunda

∆2 𝑓(𝑥𝑖 ) = ∆𝑓 (𝑥𝑖+1 ) − ∆𝑓(𝑥𝑖 )Type equation here.

Tercera

∆3 𝑓(𝑥𝑖 ) = ∆2 𝑓(𝑥𝑖+1 ) − ∆2 𝑓(𝑥𝑖 )

.
.
.
∆𝑛 𝑓(𝑥𝑖 ) = ∆𝑛−1 𝑓(𝑥𝑖+1 ) − ∆𝑛−1 𝑓(𝑥𝑖 )
𝑥−𝑥
0
Definiendo el nuevo parámetro 𝑠 = , entonces el polinomio de Newton en diferencias finitas
ℎ
hacia adelante puede expresarse como:

𝑠(𝑠 − 1) 2 𝑠(𝑠 − 1)(𝑠 − 2) 3

𝑝𝑛 (𝑥) = 𝑓(𝑥0 ) + 𝑠∆𝑓(𝑥0 ) + ∆ 𝑓(𝑥0 ) + ∆ 𝑓(𝑥0 ) + ⋯
2! 3!
𝑠(𝑠 − 1)(𝑠 − 2) … (𝑠 − (𝑠 − 1)) 𝑛
+ ∆ 𝑓(𝑥0 )
𝑛!