Scribd
Cargar
122 vistas19 páginas

Interpolación Spline Cúbico: Teoría y Ejemplos

Este documento explica la interpolación spline cúbica. Define un spline cúbico como una función polinómica cúbica por tramos que es dos veces continuamente diferenciable. Explica cómo determinar los coeficientes de los polinomios spline para que satisfagan las condiciones de interpolación en los puntos de datos e igualdad de derivadas primera y segunda en los puntos de unión. Distingue entre spline cúbico natural, periódico y forzado dependiendo de las condiciones impuestas en los extremos. Incluye un ejemplo numérico para ilustrar

Cargado por

Tu Persona Favorita
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PPTX, PDF, TXT o lee en línea desde Scribd
122 vistas19 páginas

Interpolación Spline Cúbico: Teoría y Ejemplos

Este documento explica la interpolación spline cúbica. Define un spline cúbico como una función polinómica cúbica por tramos que es dos veces continuamente diferenciable. Explica cómo determinar los coeficientes de los polinomios spline para que satisfagan las condiciones de interpolación en los puntos de datos e igualdad de derivadas primera y segunda en los puntos de unión. Distingue entre spline cúbico natural, periódico y forzado dependiendo de las condiciones impuestas en los extremos. Incluye un ejemplo numérico para ilustrar

Cargado por

Tu Persona Favorita
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PPTX, PDF, TXT o lee en línea desde Scribd

Interpolación Spline

Por Rosa Mercedes Garrido Juárez

Prof. Ing. R.M.G.J. 1


Interpolación Spline
Sea x una subdivisión del intervalo [a b]

x : a  x0  x1  ...  xn 1  xn  b
Con intervalos [x j x j 1 ] j=0,..,n-1
Definición: Una función S:[a b]  Es decir spline cúbico, si
2
• S C [ a b] es decir, S es dos veces continuamente diferenciable en [a, b];
• S Es un polinomio cúbico en cada parte del intervalo [x j x j 1 ] 0  j  n  1
S x j   f x x j x j1   s j ( x)  d j  c j ( x  x j )  b j ( x  x j ) 2  a j ( x  x j )3
Objetivo: Interpolación de datos x , f , 0  j  n, con un spline cúbico
j j S,

Tal que:
S x j   f x j , para 0  j  n
Prof. Ing. R.M.G.J. 2
Interpolación con Spline Cúbico
Observaciones:
• Un spline cúbico posee 4n parámetros, los cuales son determinados como
siguen:
• Propiedades de la Interpolación
S x j 1   f j 1 y S x j   f j para todo 1  j  n;
• Continuidad en la derivada
S x j   S j 1 x j  para todo 1  j  n  1;
• Continuidad en la segunda derivada
S x j   S j1 x j  para todo 1  j  n  1;

• Se trata de un total 4N - 2 ecuaciones para los parámetros 4N

• OBS! Faltan dos condiciones por determinar

Prof. Ing. R.M.G.J. 3


Las otras dos condiciones

Definición: Clases de Spline Cúbico


• Spline cúbico natural, si S a   S b   0
• Spline cúbico periódico, si S (i )
a   S (i )
b , para i  0,1,2;
• Spline cúbico forzado o completo, si S a  f (a) y S b  f (b)
Nota: cada una de las tres condiciones antes descritas proporciona las otras
dos condiciones (ecuaciones).
2
Proposición: Entre todas las funciones de interpolación C el spline cúbico
natural reduce al mínimo el funcional
2
I ( y )    y( x)  dx
b

Nota El funcional I mide la curvatura de y aproximado


Prof. Ing. R.M.G.J. 4
Calculo de los splines cúbicos naturales

Si S en la parte del intervalo [ x j x j 1 ] es dado por:


S x j   f x x j x j1   s j ( x)  d j  c j ( x  x j )  b j ( x  x j ) 2  a j ( x  x j )3

Tal que: dj  fj
f j 1  f j 2M j  M j 1
cj   hj
hj 6
Mj
bj 
2 0  j  n 1
M j 1  M j
aj 
6h j
Donde h j  x j 1  x j para todo
Los momentos M j  S ( x j ) se obtienen a partir de un sistema tridiagonal.

Prof. Ing. R.M.G.J. 5


Deducción de los coeficientes

Prof. Ing. R.M.G.J. 6


Deducción de los coeficientes

Prof. Ing. R.M.G.J. 7


Deducción de la ecuación de Momentos

Con la condición de S j 1 x j   S j x j  para xj, 1 j  n-1,

M 1 ,..M n 1

Prof. Ing. R.M.G.J. 8


Deducción de la ecuación de Momentos
Caso: Spline Natural

Para 1 j  n-1, y para las n-1 incógnitas M 1 ,..M n 1


Nota: Los momentos M0= 0 y Mn = 0, son conocidos.
Prof. Ing. R.M.G.J. 9
Sistema tridiagonal de los momentos

Prof. Ing. R.M.G.J. 10


Ejemplo de spline natural

Prof. Ing. R.M.G.J. 11


Ejemplo de spline natural
Solución

Coeficientes de los splines

Polinomio interpolante o spline cúbico

Para el primer spline

Prof. Ing. R.M.G.J. 12


Spline Forzado (frontera Sujeta)

S0 x0   A y Sn xn   B , con lo cual se agregan dos ecuaciones:

2h0 M 0  h0 M 1  6 yx0 , x1   A
hn 1 M n 1  2hn 1 M n  6B  yxn 1 , x n 
Y llegamos a tener n+1 ecuaciones con n+1 incógnitas:

2h0 h0 0  0  M 0   yx0 , x1   A 
h 2h0  h1  h1     M 1   yx , x   f x , x  
 0  1 2 0 1 
 0    0      6  
    
   hn 2 2hn 2  hn 1  hn 1   M n 1   y x n 1 , x n   f x n  2 , x n 1 
 0  0 hn1 2hn 1   M n   B  yx n 1 , x n  

Prof. Ing. R.M.G.J. 13


Ejemplo Spline Forzado

 Obtener una interpolación por Spline Cúbico forzado


para f(x)=(x-1/2)^4 en x=-0.5 , 0, 1. Se pide:
 Mostrar las funciones Spline S(x) para cada intervalo.
 Demuestre que las funciones Spline cumple las
condiciones mínimas.
 Interpole para x=0.25 y x=0.75 y determine el error
cometido en cada caso.

Prof. Ing. R.M.G.J. 14


Ejemplo Spline Forzado
x -1/2 0 1
y 1 1/16 1/16

h0=0.5 h1=1
y[x0,x1]=-15/8 y[x1,x2]=0

A  f  1 / 2  4 B  f 1  1 / 2
2h0 h0 0  M 0   yx0 x1   A 
h 2 h  h  h   M   6 yx x   yx x 
 0 0 1 1  1   1 2 0 1 

 0 h1 2h1   M 2   B  yx1 x2  
 1 1 / 2 0  M 0  17 / 8  51 / 4 
1 / 2 3 1  M   615 / 8  45 / 4
  1     
 0 1 2  M 2   1 / 2   3 
M 0  12 M 1  3 / 2 M 2  3 / 4
Prof. Ing. R.M.G.J. 15
Ejemplo Spline Forzado

M i 1  M i
a0 = -7/2 a1=-1/8 ai 
6hi
b0 = 6 b1=3/4 Mi
bi 
c0 = -4 c1=-5/8 2
d0 = 1 d1=1/16 M i 1  2M i
ci  yxi xi 1   hi
6
d i  yi
f’(-1/2)=-4 f ‘(1)=1/2

Prof. Ing. R.M.G.J. 16


Ejemplo Spline Forzado

a) S0(x)= -7/2(x+1/2)3+6(x+1/2)2-4(x+1/2)+1 -1/2<= x<=0

S1(x)=-1/8 x3+3/4 x2-5/8x+1/16 0<= x<=1

b) S0(0) = S1(0)= 1/16


S0(xj)=yj S1(xj)=yj j=0,1,2

S’0(x) = -21/2(x+1/2)2+12(x+1/2)-4 S’1(x)= 3/8x2+3/2x-5/8


S’0(0) = S’1(0)=-5/8

S’’0(x) = -21(x+1/2)+12 S’’1(x)= 3/4x+3/2


S’’0(0) = S’’1(0)=3/2

c)
S1(0.25)= -0.0923 f(0.25)= 0.0039 Error1= |f(0.25)- S1(0.25)| =0.088
S1(0.75)=-0.0449 f(0.75)= 0.0039 Error2= |f(0.75)- S1(0.75)| =0.041

Prof. Ing. R.M.G.J. 17


Problema de Spline Natural

Prof. Ing. R.M.G.J. 18


Solución

Prof. Ing. R.M.G.J. 19

También podría gustarte