Métodos Numéricos Splines

5.8 INTERPOLACIÓN POR SPLINES

Particionamos el intervalo [a,b] en n sub-intervalos [x j , x j1 ] , j  1, n  1 . La función
interpolante es ahora en cada sub-intervalo un polinomio de grado m tal que los
polinomios se encuentran en los xj , j  1, n  1 , (m-1)-veces diferenciables.

5.8.1 SPLINES LINEALES

Dados los puntos interpolantes a  x 0  x1  x 2   xn  b

Con los valores correspondientes y0 , , yn . Utilizando la interpolación por splines

lineales se conecta los puntos (x 0 , y0 ), , (x n , y n ) solamente por segmentos
formando una función continua sobre [a,b].

Sn= 1

S1

S0

X0 Xn

Los splines lineales se calculan mediante la fórmula

y j1  y j
S j (x)  y j  (x  x j ) x j  x  x j1 j  0, n  1
x j1  x j

El error es acotado por

h2
max  f (x)  S(x)   8 max  f (x)  h  max x j1  x j
j 0,n 1
x[a,b] x[a,b]

5.8.2 SPLINES CÚBICOS:

Dados los puntos a  x 0  x1  x 2   x n  b con sus valores y 0 , y1 , , yn

S n=1

S1
.
S0 Sj
S j+1

hj hj+1

X1 Xj Xj+1 Xj+2 Xn Xn
X0

Dra. Pilar Nuñez Blas Página 1

Métodos Numéricos Splines

Definamos un polinomio cúbico S, denotado por Sj en el sub-intervalo [x j , x j1 ]

j  1, n  1 tal que
a) S j (x j )  y j
b) S j (x j1 )  S j1 (x j1 )
c) S j (x j1 )  S j1 (x j1 )
d) S''j (x j1 )  S''j1 (x j1 )
e) Que satisfaga una de las condiciones de frontera

1) S(x 0 )  S(x n )  0 (FRONTERA LIBRE)

2) S(x 0 )  m 0 , S(x n )  m1 (FRONTERA SUJETA)
3) S(x 0 )  S(x n ) , S(x 0 )  S(x n ) (FRONTERA PERIÓDICA)

DEDUCCIÓN DE LOS SPLINES CÚBICOS

Consideremos

Sj (x)  a j  b j (x  x j )  c j (x  x j )2  d j (x  x j )3 , j  0, n  1 (1)

Según a) S j (x j )  a j  y j  a j  y j

Aplicando b)

a j1  S j1 (x j1 )  S j (x j1 ) 

a j1  Sj (x j1 )  a j  b j (x j1  x j )  c j (x j1  x j )2  d j (x j1  x j )3 , j  0, n  2

Definamos h j  x j1  x j 

a j1  a j  b jh j  c jh 2j  d jh 3j (2)

De (1)

S'j (x)  b j  2c j (x  x j )  3d j (x  x j )2

S'j (x j )  b j

Según c)

b j1  S'j1 (x j1 )  S'j (x j1 )  b j  2c j (x j1  x j )  3d j (x j1  x j )2

 b j1  b j  2c jh j  3d jh 2j (3)

De (1)

Dra. Pilar Nuñez Blas Página 2

Métodos Numéricos Splines

S''j (x)  2c j (x  x j )  6d j (x  x j )

 S''j (x j )  2c j

Aplicando d)

2c j1  S"j1 (x j1 )  S"j (x j1 )  2c j  6d jh j

c j1  c j
 dj  (4)
3h j

Al reemplazar (4) en (2) y (3) se tiene

hj
a j1  a j  b jh j  (2c j  c j1 ) (5)
3

b j1  b j  h j (c j  c j1 ) (6)

De (5)

1 hj
bj  (a j1  a j )  (2c j  c j1 ) (7)
hj 3

De (7)

1 h j1
b j1  (a j  a j1 )  (2c j1  c j ) (8)
h j1 3

De (6)

b j  b j1  h j1 (c j1  c j ) (9)

Al sustituir (7) y (8) en (9) se tiene

3 3
h j1c j1  2(h j1  h j )c j  h jc j1  (a j1  a j )  (a j  a j1 ) j  1, n  1 (10)
hj h j 1

FRONTERA LIBRE

2c n  S(x n )  0  cn=0

2c0  S(x 0 )  0  c0=0

Luego, de esto y de (10) se produce el sistema Ax  B , donde

Dra. Pilar Nuñez Blas Página 3

Métodos Numéricos Splines

1 0 0 0 
h 2(h 0  h1 ) h1 0 0 
 0
0 h1 2(h1  h 2 ) h1 0 
 
A0 0 
 0 
 
0 hn2 2(h n  2  h n 1 ) h n 1 
 0 0 0 0 1 


es una matriz estrictamente diagonal dominante.

 0 
 3 3 
 (a 2  a1 )  (a 1  a 0 ) 
 h1 h0 
 
B  
 3 3 
 (a n  a n 1 )  (a n 1  a n  2 ) 
 h n 1 hn2 
 0 
 

 c0  0 
 c 
 1 
X 
 
 c n 1 
 c n  0 

De este sistema se obtiene los cj , de (7) los bj y de (4) los dj.

FRONTERA SUJETA

m 0  S(x 0 )  b 0

Según la ec. (7)

a1  a 0 h 0
m0  b0   (2c0  c1 )
h0 3

3
 2h 0 c0  h 0c1  (a1  a 0 )  3m0 (11)
h0

Similarmente,

m1  S(x n )  b n (++)

Pero, de (6) se obtiene

b n  b n 1  h n 1 (c n 1  c n ) (+++)

Dra. Pilar Nuñez Blas Página 4

Métodos Numéricos Splines

y de (7) se tiene

1 h
bn 1  (a n  a n 1 )  n 1 (2cn 1  cn ) (++++)
h n 1 3

Por lo tanto, de (++)-(++++) se tiene

3
h n 1cn 1  2h n 1cn  3m1  (a n  a n 1 ) (12)
h n 1

Luego de (10) , (11) y (12) se tiene el sistema Ax  B donde

 2h 0 0 0 0 
h 2(h 0  h1 ) h1 0 0 
 0
 0 h1 2(h1  h 2 ) h1 0 
 
A 0 0 
 0 
 
 0 hn2 2(h n  2  h n 1 ) h n 1 
 0 0 0 0 2h n 1 


 3 
 (a1  a 0 )  3m 0 
h0
 
 3 3 
 (a 2  a1 )  (a 1  a 0 ) 
 h1 h0 
B  
 
 3 (a  a )  3 (a  a ) 
 h n 1 n n 1
hn2
n 1 n2 
 
 3 
 3 m 1  (a n  a n 1 ) 
 h n 1 

 c0 
 c 
 1 
X 
 
 c n 1 
 c n 

ERROR

Si S(a)  m0 , S(b)  m1 , f  C 4 ([a, b]) a  x 0  x1   xn  b 

5h 4
max  f (x)  S(x)  
x[a,b]

max f ( 4) (x)
384 [a,b]
x

donde h  max  x j1  x j 
0 j n 1

Dra. Pilar Nuñez Blas Página 5

Métodos Numéricos Splines

Dra. Pilar Nuñez Blas Página 6