Articulo de referencia

Interpolación discreta mediante splines

En el campo matemático del análisis numérico , la interpolación de splines discretos es una forma de interpolación donde el interpolante es un tipo especial de polinomio por par...

En el campo matemático del análisis numérico , la interpolación de splines discretos es una forma de interpolación donde el interpolante es un tipo especial de polinomio por partes llamado spline discreto. Un spline discreto es un polinomio por partes cuyas diferencias centrales son continuas en los nodos, mientras que un spline es un polinomio por partes cuyas derivadas son continuas en los nodos. Los splines cúbicos discretos son splines discretos donde se requiere que las diferencias centrales de órdenes 0, 1 y 2 sean continuas. [ 1 ]

Las splines discretas fueron introducidas por Mangasarin y Schumaker en 1971 como soluciones de ciertos problemas de minimización que involucran diferencias. [ 2 ]

Splines cúbicos discretos

Sea x 1 , x 2 , . . ., x n -1 una sucesión creciente de números reales. Sea g ( x ) un polinomio a trozos definido por

gramo(incógnita)={gramo1(incógnita)incógnita<incógnita1gramoi(incógnita)incógnitai1incógnita<incógnitai para i=2,3,,norte1gramonorte(incógnita)incógnitaincógnitanorte1{\displaystyle g(x)={\begin{cases}g_{1}(x)&x<x_{1}\\g_{i}(x)&x_{i-1}\leq x<x_{i}{\text{ para }}i=2,3,\ldots ,n-1\\g_{n}(x)&x\geq x_{n-1}\end{cases}}}

donde g 1 ( x ), . . ., g n ( x ) son polinomios de grado 3. Sea h > 0. Si

(gramoi+1gramoi)(incógnitai+jh)=0 para j=1,0,1 y i=1,2,,norte1{\displaystyle (g_{i+1}-g_{i})(x_{i}+jh)=0{\text{ para }}j=-1,0,1{\text{ y }}i=1,2,\ldots ,n-1}

entonces g ( x ) se denomina spline cúbico discreto. [ 1 ]

Formulación alternativa 1

Las condiciones que definen una spline cúbica discreta son equivalentes a las siguientes:

gramoi+1(incógnitaih)=gramoi(incógnitaih){\displaystyle g_{i+1}(x_{i}-h)=g_{i}(x_{i}-h)}
gramoi+1(incógnitai)=gramoi(incógnitai){\displaystyle g_{i+1}(x_{i})=g_{i}(x_{i})}
gramoi+1(incógnitai+h)=gramoi(incógnitai+h){\displaystyle g_{i+1}(x_{i}+h)=g_{i}(x_{i}+h)}

Formulación alternativa 2

Las diferencias centrales de órdenes 0, 1 y 2 de una función f ( x ) se definen de la siguiente manera:

D(0)F(incógnita)=F(incógnita){\displaystyle D^{(0)}f(x)=f(x)}
D(1)F(incógnita)=F(incógnita+h)F(incógnitah)2h{\displaystyle D^{(1)}f(x)={\frac {f(x+h)-f(xh)}{2h}}}
D(2)F(incógnita)=F(incógnita+h)2F(incógnita)+F(incógnitah)h2{\displaystyle D^{(2)}f(x)={\frac {f(x+h)-2f(x)+f(xh)}{h^{2}}}}

Las condiciones que definen una spline cúbica discreta también son equivalentes a [ 1 ].

D(j)gramoi+1(incógnitai)=D(j)gramoi(incógnitai) para j=0,1,2 y i=1,2,,norte1.{\displaystyle D^{(j)}g_{i+1}(x_{i})=D^{(j)}g_{i}(x_{i}){\text{ para }}j=0,1,2{\text{ y }}i=1,2,\ldots ,n-1.}

Esto indica que las diferencias centralesD(j)gramo(incógnita){\displaystyle D^{(j)}g(x)}son continuas en x i .

Ejemplo

Sea x 1 = 1 y x 2 = 2 de modo que n = 3. La siguiente función define una spline cúbica discreta: [ 1 ]

gramo(incógnita)={incógnita3incógnita<1incógnita32(incógnita1)((incógnita1)2h2)1incógnita<2incógnita32(incógnita1)((incógnita1)2h2)+(incógnita2)((incógnita2)2h2)incógnita2{\displaystyle g(x)={\begin{cases}x^{3}&x<1\\x^{3}-2(x-1)((x-1)^{2}-h^{2})&1\leq x<2\\x^{3}-2(x-1)((x-1)^{2}-h^{2})+(x-2)((x-2)^{2}-h^{2})&x\geq 2\end{cases}}}

Interpolador de spline cúbico discreto

Sea x 0 < x 1 y x n > x n -1 y f ( x ) una función definida en el intervalo cerrado [ x 0 - h, x n + h]. Entonces existe una única spline cúbica discreta g ( x ) que satisface las siguientes condiciones:

gramo(incógnitai)=F(incógnitai) para i=0,1,,norte.{\displaystyle g(x_{i})=f(x_{i}){\text{ para }}i=0,1,\ldots ,n.}
D(1)gramo1(incógnita0)=D(1)F(incógnita0).{\displaystyle D^{(1)}g_{1}(x_{0})=D^{(1)}f(x_{0}).}
D(1)gramonorte(incógnitanorte)=D(1)F(incógnitanorte).{\displaystyle D^{(1)}g_{n}(x_{n})=D^{(1)}f(x_{n}).}

Esta spline cúbica discreta única es la interpolación spline discreta de f ( x ) en el intervalo [ x 0 - h, x n + h]. Esta interpolación coincide con los valores de f ( x ) en x 0 , x 1 , . . ., x n .

Aplicaciones

  • Las splines cúbicas discretas se introdujeron originalmente como soluciones de ciertos problemas de minimización. [ 1 ] [ 2 ]
  • Tienen aplicaciones en el cálculo de splines no lineales. [ 1 ] [ 3 ]
  • Se utilizan para obtener una solución aproximada de un problema de valores en la frontera de segundo orden. [ 4 ]
  • Se han utilizado splines interpolatorios discretos para construir wavelets biorthogonales. [ 5 ]

Referencias

  1. 1 2 3 4 5 6 Tom Lyche (1979). "Interpolación de spline cúbico discreto". BIT . 16 (3): 281– 290. doi : 10.1007/bf01932270 . S2CID 122300608 . 
  2. 1 2 Mangasarian, OL; Schumaker, LL (1971). "Splines discretos mediante programación matemática". SIAM J. Control . 9 (2): 174– 183. doi : 10.1137/0309015 .
  3. Michael A. Malcolm (abril de 1977). "Sobre el cálculo de funciones spline no lineales". SIAM Journal on Numerical Analysis . 14 (2): 254– 282. Bibcode : 1977SJNA...14..254M . doi : 10.1137/0714017 .
  4. Fengmin Chen, Wong, PJY (dic. 2012). "Resolución de problemas de valores en la frontera de segundo orden mediante splines cúbicos discretos". Control Automation Robotics & Vision (ICARCV), 12.ª Conferencia Internacional de 2012 : 1800–1805 .{{cite journal}}: CS1 maint: varios nombres: lista de autores ( enlace )
  5. Averbuch, AZ, Pevnyi, AB, Zheludev, VA (noviembre de 2001). "Ondículas de Butterworth biorthogonales derivadas de splines interpolatorios discretos". IEEE Transactions on Signal Processing . 49 (11): 2682– 2692. Bibcode : 2001ITSP...49.2682A . CiteSeerX 10.1.1.332.7428 . doi : 10.1109/78.960415 . {{cite journal}}: CS1 maint: varios nombres: lista de autores ( enlace )
Obtenido de " https://en.wikipedia.org/w/index.php?title=Discrete_spline_interpolation&oldid=1279884267 "