Polinomio de Newton

format_list_bulleted Contenido keyboard_arrow_down
ImprimirCitar
Expresión matemática

En el campo matemático del análisis numérico, un polinomio de Newton, llamado así por su inventor Isaac Newton, es un polinomio de interpolación para un conjunto dado de puntos de datos. El polinomio de Newton a veces se denomina polinomio de interpolación de diferencias divididas de Newton porque los coeficientes del polinomio se calculan utilizando el método de diferencias divididas de Newton.

Definición

Dado un conjunto de k + 1 puntos de datos

()x0,Sí.0),...... ,()xj,Sí.j),...... ,()xk,Sí.k){displaystyle (x_{0},ldots(x_{j},y_{j}),ldots(x_{k},y_{k})}

donde no hay dos xj iguales, el polinomio de interpolación de Newton es una combinación lineal de polinomios de base de Newton

N()x):=.. j=0kajnj()x){displaystyle N(x):=sum _{j=0}{k}a_{j}n_{j}(x)}

con los polinomios de base de Newton definidos como

nj()x):=∏ ∏ i=0j− − 1()x− − xi){displaystyle n_{j}(x):=prod _{i=0}{j-1}(x-x_{i}}

para j " 0 " n0()x)↑ ↑ 1{displaystyle n_{0}(x)equiv 1}.

Los coeficientes se definen como

aj:=[Sí.0,...... ,Sí.j]{displaystyle a_{j}:=[y_{0},ldotsy_{j}}

dónde

[Sí.0,...... ,Sí.j]{displaystyle [y_{0},ldotsy_{j}}

es la notación para diferencias divididas.

Por lo tanto, el polinomio de Newton se puede escribir como

N()x)=[Sí.0]+[Sí.0,Sí.1]()x− − x0)+⋯ ⋯ +[Sí.0,...... ,Sí.k]()x− − x0)()x− − x1)⋯ ⋯ ()x− − xk− − 1).{displaystyle N(x)=[y_{0}]+[y_{0},y_{1}](x-x_{0})+cdots +[y_{0},ldotsy_{k}](x-x_{0})(x-x_{1})cdots (x-x_{k-1}). }

Fórmula de diferencia dividida hacia delante de Newton

El polinomio de Newton se puede expresar en forma simplificada cuando x0,x1,...... ,xk{displaystyle x_{0},x_{1},dotsx_{k} se organizan consecutivamente con igual espaciamiento. Introducción de la notación h=xi+1− − xi{displaystyle h=x_{i+1}-x_{i} para cada uno i=0,1,...... ,k− − 1{displaystyle i=0,1,dotsk-1}y x=x0+sh{displaystyle x=x_{0}+sh}, la diferencia x− − xi{displaystyle x-x_{i}} puede ser escrito como ()s− − i)h{displaystyle (s-i)h}. Así que el polinomio Newton se convierte

N()x)=[Sí.0]+[Sí.0,Sí.1]sh+⋯ ⋯ +[Sí.0,...... ,Sí.k]s()s− − 1)⋯ ⋯ ()s− − k+1)hk=.. i=0ks()s− − 1)⋯ ⋯ ()s− − i+1)hi[Sí.0,...... ,Sí.i]=.. i=0k()si)i!hi[Sí.0,...... ,Sí.i].{displaystyle {begin{aligned}N(x) ventaja=[y_{0}+[y_{0},y_{1}]sh+cdots [y_{0},ldotsy_{k}]s(s-1)cdots (s-k+1){h}^{k}\\=sum ¿Por qué? {fnMicrosoft Sans Serif}

Esto se llama la fórmula de diferencia dividida hacia delante de Newton.

Fórmula de diferencia dividida hacia atrás de Newton

Si los nodos son reordenados como xk,xk− − 1,...... ,x0{displaystyle {x},{x}_{k-1},dots{x}, el polinomio de Newton se convierte

N()x)=[Sí.k]+[Sí.k,Sí.k− − 1]()x− − xk)+⋯ ⋯ +[Sí.k,...... ,Sí.0]()x− − xk)()x− − xk− − 1)⋯ ⋯ ()x− − x1).{displaystyle N(x)=[y_{k}]+[{y},{y}_{k-1}](x-{x}_{k})+cdots (x-{x})(x-{x}_{k-1})cdots (x-{x}_{1}). }

Si xk,xk− − 1,...... ,x0{displaystyle {x}_{k},;{x}_{k-1},;dots;{x}_{0}} are equally spaced with x0=xk+sh{displaystyle {x}}={x}_{k}+sh} y xi=xk− − ()k− − i)h{displaystyle {fnMicrosoft Sans Serif} para i = 0, 1,... k, entonces,

N()x)=[Sí.k]+[Sí.k,Sí.k− − 1]sh+⋯ ⋯ +[Sí.k,...... ,Sí.0]s()s+1)⋯ ⋯ ()s+k− − 1)hk=.. i=0k()− − 1)i()− − si)i!hi[Sí.k,...... ,Sí.k− − i].{displaystyle {begin{aligned}N(x) ventaja=[{y}_{y}_{k},{y}_{y}_{k-1}]sh+cdots [{y}_{k},ldots{y}_{0}s(s+1)cdots (s+k-1) {h} {k}\\cH00} - ¿Qué? {fnMicrosoft Sans Serif}

se llama la fórmula de diferencia dividida hacia atrás de Newton.

Importancia

La fórmula de Newton es interesante porque es la versión en diferencias directa y natural del polinomio de Taylor. El polinomio de Taylor indica adónde irá una función, según su valor y y sus derivadas (su tasa de cambio y la tasa de cambio de su tasa de cambio, etc.) en un valor particular de x. La fórmula de Newton es el polinomio de Taylor basado en diferencias finitas en lugar de tasas de cambio instantáneas.

Adición de nuevos puntos

Al igual que con otras fórmulas de diferencias, el grado de un polinomio de interpolación de Newton se puede aumentar agregando más términos y puntos sin descartar los existentes. La forma de Newton tiene la simplicidad de que los nuevos puntos siempre se agregan en un extremo: la fórmula directa de Newton puede agregar nuevos puntos a la derecha y la fórmula inversa de Newton puede agregar nuevos puntos a la izquierda.

La precisión de la interpolación polinomial depende de cuán cerca esté el punto interpolado de la mitad de los valores x del conjunto de puntos utilizados. Obviamente, a medida que se agregan nuevos puntos en un extremo, ese medio se aleja cada vez más del primer punto de datos. Por lo tanto, si no se sabe cuántos puntos se necesitarán para la precisión deseada, la mitad de los valores de x podría estar lejos de donde se realiza la interpolación.

Gauss, Stirling y Bessel desarrollaron fórmulas para remediar ese problema.

La fórmula de Gauss agrega alternativamente nuevos puntos en los extremos izquierdo y derecho, manteniendo así el conjunto de puntos centrados cerca del mismo lugar (cerca del punto evaluado). Al hacerlo, utiliza términos de la fórmula de Newton, con puntos de datos y valores de x renombrados de acuerdo con la elección de qué punto de datos se designa como x 0 punto de datos.

La fórmula de Stirling permanece centrada en un punto de datos en particular, para usarse cuando el punto evaluado está más cerca de un punto de datos que del medio de dos puntos de datos.

La fórmula de Bessel permanece centrada en un medio particular entre dos puntos de datos, para usarse cuando el punto evaluado está más cerca de un medio que de un punto de datos.

Bessel y Stirling logran eso usando a veces el promedio de dos diferencias, y a veces usando el promedio de dos productos de binomios en x, donde la de Newton o la de Gauss use solo una diferencia o producto. Stirling's usa una diferencia promedio en términos de grado impar (cuya diferencia usa un número par de puntos de datos); Bessel's usa una diferencia promedio en términos de grado par (cuya diferencia usa un número impar de puntos de datos).

Fortalezas y debilidades de varias fórmulas

Para cualquier conjunto finito de puntos de datos, solo hay un polinomio de menor grado posible que pasa por todos ellos. Así, es apropiado hablar de la "forma de Newton", o forma de Lagrange, etc., del polinomio de interpolación. Sin embargo, diferentes métodos para calcular este polinomio pueden tener diferentes eficiencias computacionales. Existen varios métodos similares, como los de Gauss, Bessel y Stirling. Se pueden derivar de los de Newton cambiando el nombre de los valores x de los puntos de datos, pero en la práctica son importantes.

Bessel contra Stirling

La elección entre Bessel y Stirling depende de si el punto interpolado está más cerca de un punto de datos o más cerca de un punto medio entre dos puntos de datos.

El error de interpolación de un polinomio se aproxima a cero, a medida que el punto de interpolación se aproxima a un punto de datos. Por lo tanto, la fórmula de Stirling trae su mejora de precisión donde menos se necesita y Bessel trae su mejora de precisión donde más se necesita.

Por lo tanto, se podría decir que la fórmula de Bessel es la fórmula de diferencia más consistentemente precisa y, en general, la más consistentemente precisa de las fórmulas familiares de interpolación de polinomios.

Métodos de diferencias divididas frente a Lagrange

A veces se dice que Lagrange requiere menos trabajo y, a veces, se recomienda para problemas en los que se sabe, de antemano, por experiencia previa, cuántos términos se necesitan para una precisión suficiente.

Los métodos de diferencias divididas tienen la ventaja de que se pueden agregar más puntos de datos para mejorar la precisión. Se pueden seguir utilizando los términos basados en los puntos de datos anteriores. Con la fórmula ordinaria de Lagrange, resolver el problema con más puntos de datos requeriría volver a hacer todo el problema.

Hay un "baricéntrico" versión de Lagrange que evita la necesidad de rehacer todo el cálculo al agregar un nuevo punto de datos. Pero requiere que se registren los valores de cada término.

Pero la capacidad de Gauss, Bessel y Stirling de mantener los puntos de datos centrados cerca del punto interpolado les da una ventaja sobre Lagrange, cuando no se sabe de antemano cuántos puntos de datos serán necesario.

Además, suponga que desea averiguar si, para algún tipo particular de problema, la interpolación lineal es lo suficientemente precisa. Eso se puede determinar evaluando el término cuadrático de una fórmula de diferencia dividida. Si el término cuadrático es insignificante, lo que significa que el término lineal es lo suficientemente preciso sin agregar el término cuadrático, entonces la interpolación lineal es lo suficientemente precisa. Si el problema es lo suficientemente importante, o si el término cuadrático es lo suficientemente grande como para tener importancia, entonces uno podría querer determinar si la suma de los términos cuadráticos y cúbicos es lo suficientemente grande como para tener importancia en el problema.

Por supuesto, solo se puede usar un método de diferencias divididas para tal determinación.

Para ello, se debe elegir la fórmula de la diferencia dividida y/o su punto x0 de modo que la fórmula utilice, para su término lineal, los dos puntos de datos entre los que se realizaría la interpolación lineal de interés.

Las fórmulas de diferencias divididas son más versátiles, útiles en más tipos de problemas.

La fórmula de Lagrange es óptima cuando toda la interpolación se realiza en un valor de x, con solo los puntos de datos' valores de y que varían de un problema a otro, y cuando se sabe, por experiencia pasada, cuántos términos se necesitan para una precisión suficiente.

Con la forma de Newton del polinomio de interpolación existe un algoritmo compacto y efectivo para combinar los términos para encontrar los coeficientes del polinomio.

Precisión

Cuando, con Stirling's o Bessel's, el último término usado incluye el promedio de dos diferencias, entonces se está usando un punto más que el que usaría Newton's u otras interpolaciones polinómicas para el mismo grado polinomial. Entonces, en ese caso, Stirling's o Bessel's no están poniendo un polinomio de N−1 grado a través de N puntos, sino que están intercambiando equivalencia con la de Newton para un mejor centrado y precisión, dando a esos métodos a veces una precisión potencialmente mayor, para un grado polinomial dado, que otras interpolaciones polinómicas.

Caso general

Para el caso especial xi = i, hay un conjunto estrechamente relacionado de polinomios, también llamados los polinomios de Newton, que son simplemente los coeficientes binomiales para el argumento general. Es decir, uno también tiene los polinomios de Newton pn()z){displaystyle p_{n}(z)} dado por

pn()z)=()zn)=z()z− − 1)⋯ ⋯ ()z− − n+1)n!{displaystyle p_{n}(z)={z choose n}={frac {z(z-1)cdots (z-n+1)}{n!}}

De esta forma, los polinomios de Newton generan la serie de Newton. Estos son a su vez un caso especial de los polinomios en diferencias generales que permiten la representación de funciones analíticas a través de ecuaciones en diferencias generalizadas.

Idea principal

Resolver un problema de interpolación lleva a un problema de álgebra lineal donde tenemos que resolver un sistema de ecuaciones lineales. Usando una base monomio estándar para nuestro polinomio de interpolación, obtenemos la muy complicada matriz de Vandermonde. Al elegir otra base, la base de Newton, obtenemos un sistema de ecuaciones lineales con una matriz triangular inferior mucho más simple que se puede resolver más rápido.

Para k + 1 puntos de datos, construimos la base de Newton como

n0()x):=1,nj()x):=∏ ∏ i=0j− − 1()x− − xi)j=1,...... ,k.{displaystyle n_{0}(x):=1,qquad n_{j}(x):=prod _{i=0}^{j-1}(x-x_{i})qquad j=1,ldotsk.}

Utilizar estos polinomios como base para ▪ ▪ k{displaystyle Pi _{k} tenemos que resolver

[1...... 01x1− − x01x2− − x0()x2− − x0)()x2− − x1)⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱ ⋱ 1xk− − x0...... ...... ∏ ∏ j=0k− − 1()xk− − xj)][a0⋮ ⋮ ak]=[Sí.0⋮ ⋮ Sí.k]{displaystyle {begin{bmatrix}1 ventajaldots {01} {0} {0}cH00} {begin{bmatrix}y_{0}\\\vdots {\\\\\\\\\\\\\\\\\\\\\\\\cH}\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

para resolver el problema de interpolación de polinomios.

Este sistema de ecuaciones se puede resolver iterativamente resolviendo

.. i=0jaini()xj)=Sí.jj=0,...... ,k.{displaystyle sum ¿Qué?

Derivación

Si bien la fórmula de interpolación se puede encontrar resolviendo un sistema lineal de ecuaciones, hay una pérdida de intuición en lo que muestra la fórmula y no es evidente por qué funciona la fórmula de interpolación de Newton. Para comenzar, primero necesitaremos establecer dos hechos:

Hechos 1. Revertir los términos de una diferencia dividida deja sin cambios: [Sí.0,...... ,Sí.n]=[Sí.n,...... ,Sí.0].{displaystyle [y_{0},ldotsy_{n}=[y_{n},ldotsy_{0}.}

La prueba de esto es una inducción fácil: para n=1{displaystyle n=1} nosotros computamos

[Sí.0,Sí.1]=[Sí.1]− − [Sí.0]x1− − x0=[Sí.0]− − [Sí.1]x0− − x1=[Sí.1,Sí.0].{displaystyle [y_{0},y_{1}={frac [y_{1}]- [y_{0} {x_{1}-x_{0}={frac} {f}}= {f} [y_{0}= [y_{1}]} {x_{0}= [y_{1},y_{0}

Paso de inducción: Supongamos que el resultado es válido para cualquier diferencia dividida en la mayoría n+1{displaystyle n+1} términos. Luego utilizando la hipótesis de inducción en la siguiente 2a igualdad vemos que para una diferencia dividida involucrando n+2{displaystyle n+2} términos que tenemos

[Sí.0,...... ,Sí.n+1]=[Sí.1,...... ,Sí.n+1]− − [Sí.0,...... ,Sí.n]xn+1− − x0=[Sí.n,...... ,Sí.0]− − [Sí.n+1,...... ,Sí.1]x0− − xn+1=[Sí.n+1,...... ,Sí.0].{displaystyle [y_{0},ldotsy_{n+1}={frac {[y_{1},ldotsy_{n+1}]-[y_{0},ldotsy_{n}}{x_{n+1}-x_{0}}}}={frac [y_{n},ldotsy_{0}]-[y_{n+1},ldotsy_{1}}{x_{0}-x_{n+1}=[y_{n+1},ldotsy_{0}}}

Formulamos el siguiente dato 2 que para fines de inducción y claridad también llamamos Declaración n{displaystyle n}()Stmn{displaystyle {text{Stm}_{n}})

Hechos 2. ()Stmn{displaystyle {text{Stm}_{n}}) Si ()x0,Sí.0),...... ,()xn− − 1,Sí.n− − 1){displaystyle (x_{0},ldots(x_{n-1},y_{n-1})} cualquiera n{displaystyle n} puntos con distintos puntos x{displaystyle x}-coordinados y P=P()x){displaystyle P=P(x)} es el único polinomio de grado (en la mayoría) n− − 1{displaystyle n-1} cuyo gráfico pasa por estos n{displaystyle n} puntos entonces allí tiene la relación

[Sí.0,...... ,Sí.n]()xn− − x0)⋅ ⋅ ...... ⋅ ⋅ ()xn− − xn− − 1)=Sí.n− − P()xn){displaystyle [y_{0},ldotsy_{n}](x_{n}-x_{0})cdot ldots cdot (x_{n}-x_{n-1})=y_{n}-P(x_{n})}

Prueba. (Será útil para una lectura fluida de la prueba tener la declaración precisa y su sutileza en mente: P{displaystyle P} se define por pasar ()x0,Sí.0),...,()xn− − 1,Sí.n− − 1){displaystyle (x_{0},y_{0}),...(x_{n-1},y_{n-1}} pero la fórmula también habla a ambos lados de un punto arbitrario adicional ()xn,Sí.n){displaystyle (x_{n},y_{n}} con x{displaystyle x}-coordinado distinto del otro xi{displaystyle x_{i}}.)

De nuevo demostramos estas declaraciones por inducción. Para mostrar Stm1,{displaystyle {text{Stm}_{1}} Deja ()x0,Sí.0){displaystyle (x_{0},y_{0}} ser cualquier punto y dejar P()x){displaystyle P(x)} ser el polinomio único del grado 0 pasando por ()x0,Sí.0){displaystyle (x_{0},y_{0}}. Entonces evidentemente P()x)=Sí.0{displaystyle P(x)=y_{0} y podemos escribir

[Sí.0,Sí.1]()x1− − x0)=Sí.1− − Sí.0x1− − x0()x1− − x0)=Sí.1− − Sí.0=Sí.1− − P()x1){displaystyle [y_{0},y_{1}](x_{1}-x_{0}={frac {y_{1}-y_{0} {x_{1}-x_{0}(x_{1}-x_{0})=y_{1}-y_{0}=y_{1}-P(x_{1})}

Prueba Stmn+1,{displaystyle {text{Stm}_{n+1} suposición Stmn{displaystyle {text{Stm}_{n}} ya establecido: Vamos P()x){displaystyle P(x)} ser el polinomio de grado (en la mayoría) n{displaystyle n} pasando ()x0,Sí.0),...... ,()xn,Sí.n).{displaystyle (x_{0},ldots(x_{n},y_{n}).

Con Q()x){displaystyle Q(x)} siendo el único polinomio de grado (en la mayoría) n− − 1{displaystyle n-1} pasando por los puntos ()x1,Sí.1),...... ,()xn,Sí.n){displaystyle (x_{1},y_{1}),ldots(x_{n},y_{n}}, podemos escribir la siguiente cadena de iguales, donde usamos en penúltima igualdad que Stmn{displaystyle ¿Qué? aplicable Q{displaystyle Q}:


[Sí.0,...... ,Sí.n+1]()xn+1− − x0)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)=[Sí.1,...... ,Sí.n+1]− − [Sí.0,...... ,Sí.n]xn+1− − x0()xn+1− − x0)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)=()[Sí.1,...... ,Sí.n+1]− − [Sí.0,...... ,Sí.n])()xn+1− − x1)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)=[Sí.1,...... ,Sí.n+1]()xn+1− − x1)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)− − [Sí.0,...... ,Sí.n]()xn+1− − x1)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)=()Sí.n+1− − Q()xn+1))− − [Sí.0,...... ,Sí.n]()xn+1− − x1)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)=Sí.n+1− − ()Q()xn+1)+[Sí.0,...... ,Sí.n]()xn+1− − x1)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)).##### {0}{_____# -[y_{1},n]

La hipótesis de inducción para Q{displaystyle Q} también se aplica a la segunda igualdad en el siguiente cálculo, donde ()x0,Sí.0){displaystyle (x_{0},y_{0}} se añade a los puntos que definen Q{displaystyle Q}:

Q()x0)+[Sí.0,...... ,Sí.n]()x0− − x1)⋅ ⋅ ...... ⋅ ⋅ ()x0− − xn)=Q()x0)+[Sí.n,...... ,Sí.0]()x0− − xn)⋅ ⋅ ...... ⋅ ⋅ ()x0− − x1)=Q()x0)+Sí.0− − Q()x0)=Sí.0=P()x0).{} {} {cH00}cH00}cH00}}cdots}cdotsy_{0} {} {0}}cdotcdotscdot (x_{0}} {0}{0}{0}{0} {}{0} {c}}}c}}c}c}c}c}c}}c}cdot}c}c}c}cdot}c}c]

Ahora mira. Q()x)+[Sí.0,...... ,Sí.n]()x− − x1)⋅ ⋅ ...... ⋅ ⋅ ()x− − xn).{displaystyle Q(x)+[y_{0},ldotsy_{n}(x-x_{1})cdot ldots cdots (x-x_{n}). } Por definición Q{displaystyle Q} este polinomio pasa por ()x1,Sí.1),...,()xn,Sí.n){displaystyle (x_{1},y_{1}),...,(x_{n},y_{n}} y, como acabamos de demostrar, también pasa a través de ()x0,Sí.0).{displaystyle (x_{0},y_{0}). } Así es el polinomio único de grado ≤ ≤ n{displaystyle leq n} que pasa por estos puntos. Por lo tanto este polinomio es P()x);{displaystyle P(x);} i.e.: P()x)=Q()x)+[Sí.0,...... ,Sí.n]()x− − x1)⋅ ⋅ ...... ⋅ ⋅ ()x− − xn).{displaystyle P(x)=Q(x)+[y_{0},ldotsy_{n}](x-x_{1})cdot ldots cdots (x-x_{n}). }

Así podemos escribir la última línea en la primera cadena de igualdades como:Sí.n+1− − P()xn+1){displaystyle Y...' y por lo tanto han establecido que

[Sí.0,...... ,Sí.n+1]()xn+1− − x0)⋅ ⋅ ...... ⋅ ⋅ ()xn+1− − xn)=Sí.n+1− − P()xn+1).{displaystyle [y_{0},ldotsy_{n+1}](x_{n+1}-x_{0})cdot ldots cdots (x_{n+1}-x_{n})=y_{n+1}-P(x_{n+1}).}
Stmn+1{displaystyle {text{Stm}_{n+1}

Ahora mira el Hechos 2: Se puede formular de esta manera: Si P{displaystyle P} es el único polinomio de grado en la mayoría n− − 1{displaystyle n-1} cuyo gráfico pasa por los puntos ()x0,Sí.0),...,()xn− − 1,Sí.n− − 1),{displaystyle (x_{0},y_{0}),...(x_{n-1},y_{n-1}),} entonces P()x)+[Sí.0,...... ,Sí.n]()x− − x0)⋅ ⋅ ...... ⋅ ⋅ ()x− − xn− − 1){displaystyle P(x)+[y_{0},ldotsy_{n}(x-x_{0})cdot ldots cdots (x-x_{n-1})} es el único polinomio de grado en la mayoría n{displaystyle n} pasando a través de puntos ()x0,Sí.0),...,()xn− − 1,Sí.n− − 1),()xn,Sí.n).{displaystyle (x_{0},y_{0}),...,(x_{n-1},y_{n-1}),(x_{n},y_{n}).} Así que vemos que la interpolación de Newton permite de hecho añadir nuevos puntos de interpolación sin destruir lo que ya se ha calculado.

Polinomio de Taylor

El límite del polinomio de Newton si todos los nodos coinciden es un polinomio de Taylor, porque las diferencias divididas se convierten en derivadas.

lim()x0,...... ,xn)→ → ()z,...... ,z)f[x0]+f[x0,x1]⋅ ⋅ ().. − − x0)+⋯ ⋯ +f[x0,...... ,xn]⋅ ⋅ ().. − − x0)⋅ ⋅ ⋯ ⋯ ⋅ ⋅ ().. − − xn− − 1)=f()z)+f.()z)⋅ ⋅ ().. − − z)+⋯ ⋯ +f()n)()z)n!⋅ ⋅ ().. − − z)n################################################################################################################################################################################################################################################################

Solicitud

Como se puede ver en la definición de las diferencias divididas, se pueden agregar nuevos puntos de datos al conjunto de datos para crear un nuevo polinomio de interpolación sin volver a calcular los coeficientes anteriores. Y cuando cambia un punto de datos, normalmente no tenemos que volver a calcular todos los coeficientes. Además, si las xi se distribuyen equidistantemente, el cálculo de las diferencias divididas se vuelve significativamente más fácil. Por lo tanto, las fórmulas de diferencias divididas suelen preferirse a la forma de Lagrange para fines prácticos.

Ejemplos

Las diferencias divididas se pueden escribir en forma de tabla. Por ejemplo, para una función f debe interpolarse en puntos x0,...... ,xn{displaystyle x_{0},ldotsx_{n}. Escriba

x0f()x0)f()x1)− − f()x0)x1− − x0x1f()x1)f()x2)− − f()x1)x2− − x1− − f()x1)− − f()x0)x1− − x0x2− − x0f()x2)− − f()x1)x2− − x1x2f()x2)⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ xnf()xn)#### {x}#########################################################################################################################################################################################################################################################

Luego, el polinomio de interpolación se forma como se indicó anteriormente utilizando las entradas superiores de cada columna como coeficientes.

Por ejemplo, supongamos que vamos a construir el polinomio de interpolación a f(x) = tan(x) usando diferencias divididas, en los puntos

n{displaystyle n}xn{displaystyle x_{n}f()xn){displaystyle f(x_{n}}
0{displaystyle 0}− − 32{displaystyle -{tfrac {3}{2}}− − 14.1014{displaystyle -14.1014}
1{displaystyle 1}− − 34{displaystyle -{tfrac {3}{4}}− − 0,931596{displaystyle -0.931596}
2{displaystyle 2}0{displaystyle 0}0{displaystyle 0}
3{displaystyle 3}34{fnMicroc} {3}{4}}0,931596{displaystyle 0.931596}
4{displaystyle 4}32{fnMicroc} {3}{2}}14.1014{displaystyle 14.1014}

Usando seis dígitos de precisión, construimos la tabla

− − 32− − 14.101417.5597− − 34− − 0,931596− − 10.87841.242134.8348400001.242134.83484340,93159610.878417.55973214.1014{displaystyle {begin{matrix}-{tfrac {3}{2} unos-14.1014 unos cuantos tarden\\\fnMicrosoft Sans Serif} {3}{4} {0.931596 tendríamos que-10.8784 correspondía\fnMicrosoft Sansículo1.24213 tendrían 4.834 millones de personas tendrían un problema con un cambio de sentido. {3}{4} {0.931596 tendrían que ser 10.8784 correspondían\fnMicrosoft Sans Serif} {3}{2}} {3}}} {3}} {3}}}} {3}} {3}}}}} {3} {3}} {3}}}}} {3}}}}}}}}}} {3}}}}}}}}}}}}}}}}}}}}}}} {3}}}}} {3} {3}} {3}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}} {3}}}}}} {3}}}}} {3}}}}} {3}}}} {3}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}

Por lo tanto, el polinomio de interpolación es

− − 14.1014+17.5597()x+32)− − 10.8784()x+32)()x+34)+4.83484()x+32)()x+34)()x)+0()x+32)()x+34)()x)()x− − 34)=− − 0,00005− − 1.4775x− − 0,00001x2+4.83484x3{displaystyle {begin{aligned} [3}{2})-10.8784(x+{tfrac {3}{2})(x+{tfrac [3}{4})+4.83484(x+{tfrac {3}{2})(x+{tfrac {3}{4}})(x)+0(x+{tfrac {3}})(x+{tfrac {3} {3}{4}}})(x)(x-{tfrac} {3}{4}}}}})(x)(x)(x {3}{4})={3}={3}d}} {3}} {3} {3}}}} {3}}}}}

Con más dígitos de precisión en la tabla, se encontrará que el primer y el tercer coeficiente son cero.

Otro ejemplo:

La secuencia f0{displaystyle f_{0} tales que f0()1)=6,f0()2)=9,f0()3)=2{displaystyle f_{0}(1)=6,f_{0}(2)=9,f_{0}=2} y f0()4)=5{displaystyle f_{0}(4)=5}, es decir, son 6,9,2,5{displaystyle 6,9,2,5} desde x0=1{displaystyle x_{0}=1} a x3=4{displaystyle x_{3}=4}.

Usted obtiene la pendiente de orden 1{displaystyle 1} de la siguiente manera:

  • f1()x0,x1)=f0()x1)− − f0()x0)x1− − x0=9− − 62− − 1=3{displaystyle f_{1}(x_{0},x_{1}={frac {f_{0}(x_{1})-f_{0}(x_{0}{x_{1}-x_{0}={0}={frac}={0}={0} {9-6}{2-1}=3}
  • f1()x1,x2)=f0()x2)− − f0()x1)x2− − x1=2− − 93− − 2=− − 7{displaystyle f_{1}(x_{1},x_{2}={frac {f_{0}(x_{2})-f_{0}(x_{1}}{x_{2}-x_{1}={1}}={frac}={frac} {2-9}{3-2}=-7}
  • f1()x2,x3)=f0()x3)− − f0()x2)x3− − x2=5− − 24− − 3=3{displaystyle f_{1}(x_{2},x_{3}={frac {f_{0}(x_{3})-f_{0}(x_{2}}{x_{3}-x_{2}}={frac}={frac} {5-2}{4-3}=3}}

Como tenemos las pendientes de orden 1{displaystyle 1}, es posible obtener el siguiente pedido:

  • f2()x0,x1,x2)=f1()x1,x2)− − f1()x0,x1)x2− − x0=− − 7− − 33− − 1=− − 5{displaystyle f_{2}(x_{0},x_{1},x_{2}={frac {f_{1}(x_{1},x_{2})-f_{1}(x_{0},x_{1})}{x_{2}-x_{0}}}}}={frac}}={frac {-7-3}{3-1}=-5}
  • f2()x1,x2,x3)=f1()x2,x3)− − f1()x1,x2)x3− − x1=3− − ()− − 7)4− − 2=5{fnMicrosoft Sans Serif} {2}(x_{2},x_{2})={frac {f_{1}(x_{2},x_{3})-f_{1}(x_{1},x_{2}} {x_{3}-x_{1}}}={4} {}=}{4} {}}}} {}}} {}}}}}}}{2} {}}}}} {}}}}} {}}}} {} {}}}}} {} {} {}}}}}}}}}}}} {}}} {}}}}} {}}}}}}}} {}} {} {} {}}}} {}}}}} {}}}}}}} {}}}}}}}}}}}}}}}}}}}}}} {}}}}}}}}}}}}}}}}}}}}}}}}}}}}}

Finalmente, definimos la pendiente del orden 3{displaystyle 3}:

  • f3()x0,x1,x2,x3)=f2()x1,x2,x3)− − f2()x0,x1,x2)x3− − x0=5− − ()− − 5)4− − 1=103[displaystyle f_{3}(x_{0},x_{1},x_{2},x_{3})={frac {f_{2}(x_{1},x_{2},x_{3})-f_{2}(x_{0},x_{2} {0} {0} {0} {0} {0} {0} {0}}}}} {5-(-5)}{4-1}={frac {}{3}}

Una vez que tenemos la pendiente, podemos definir los polinomios consecuentes:

  • p0()x)=6{displaystyle p_{0}(x)=6}.
  • p1()x)=6+3()x− − 1){displaystyle p_{1}(x)=6+3(x-1)}
  • p2()x)=6+3()x− − 1)− − 5()x− − 1)()x− − 2){displaystyle p_{2}(x)=6+3(x-1)-5(x-1)(x-2)}.
  • p3()x)=6+3()x− − 1)− − 5()x− − 1)()x− − 2)+103()x− − 1)()x− − 2)()x− − 3){displaystyle p_{3}(x)=6+3(x-1)-5(x-1)(x-2)+{frac {10}{3}}(x-1)(x-2)(x-3)}

Contenido relacionado

Yarda cúbica

Una yarda cúbica es una unidad de volumen imperial/usual estadounidense utilizado en Canadá y los Estados Unidos. Se define como el volumen de un cubo con...

William Jones (matemático)

William Jones, FRS fue un matemático galés, más conocido por su uso del símbolo π para representar la relación entre la circunferencia de un círculo y...

Vladimir vapnik

Vladimir Naumovich Vapnik es un informático, investigador y académico. Es uno de los principales desarrolladores de la teoría de aprendizaje estadístico...
Más resultados...
Tamaño del texto:
undoredo
format_boldformat_italicformat_underlinedstrikethrough_ssuperscriptsubscriptlink
save