Hora 6 8 10 12 14 16 18 20

Grados 7 9 12 18 21 19 15 10
Un problema de Aproximacin
Evolucin de la temperatura diurna
4
8
20
6 8 10 12 14 16 18 20 22
6
10
12
14
16
18
22
Hora
G
r
a
d
o
s

Interpolacion
Interpolacin Polinomial
Polinomios Osculadores: Interpolacin de
Hermite
Interpolacin Racional: Aproximaciones de
Pade
Interpolacin segmentaria: Splines
Otros
Ajuste
Polinomios de Taylor
Mnimos Cuadrados
Minimizacin de normas
Aproximacin Racional
Series de Fourier
Curvas de Bezier
B-Splines
Interpolacin Polinmica
Segmentaria
Limitaciones de la interpolacin polinmica
Grado del polinomio
Carcter de la funcin a interpolar
Alternativa propuesta: Splines.
Numricamente estable
Matrices dispersas
Agradable a la vista
Interpolacin Polinomica
Segmentaria: Splines
Interpolacin Segmentaria
Interpolacin Segmentaria Lineal
Interpolacin Segmentaria Cbica
Condiciones Naturales
Condiciones sobre la derivada
Interpolacin Segmentaria
Lineal: Funcin de Runge
-1 0 1
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Spline lineal
-1 0 1
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
Polinomio grado 4
y
x
=
+
1
1 25
2
Perfil para un diseo

Polinomio interpolador
Aplicaciones

Ingeniera y Diseo (CAD/CAM, CNCs)
Geologa
Aeronutica y automocin
Economa
Procesamiento de seales e imgenes (Reconocimiento
de patrones, recuperacin de imgenes)
Robtica
Medicina (Aparatos auditivos, mapas cerebrales)
Meteorologa (Mapas climticos, deteccin de
inundaciones,...)
Mundo Virtual Distribuido Multiusuario



Interpolacin Polinmica
Segmentaria
s
Dados n+1 puntos (x
0
,y
0
), (x
1
,y
1
), ..., (x
n
,y
n
) con
x
0
<x
1
<x
n
, una funcin spline de orden k (k-Spline)
sobre dichos puntos es una funcin S verificando:
(i) S(x) = q
k
(x) polinomio de grado k, x e[x
k
,x
k+1
],
k=0,1,...,n-1
(ii) S(x
k
) = y
k
, k=0,1,...,n
(iii)
| |
1
0 1
,
k
S C x x

e
Splines Lineales
Polinomio de Lagrange
Polinomio de Newton
q x
x x
x x
y
x x
x x
y
k
k
k k
k
k
k k
k
( ) =

+
+ +
+
1
1 1
1
q x f x f x x x x
y
y y
x x
x x
k k k k k
k
k k
k k
k
( ) [ ] [ , ]( )
( )
= + =
= +

+
+
+
1
1
1
Splines Lineales
Interpolacin Segmentaria
Lineal: Funcin de Runge
-1 0 1
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Spline lineal
-1 0 1
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
Polinomio grado 4
y
x
=
+
1
1 25
2
Splines Cbicos

Spline cbico

4n incgnitas
Condiciones de interpolacin
n+1 ecuaciones
Condiciones de conexin



3(n-1) ecuaciones
q x a b x x c x x d x x
k k k k k k k k
( ) ( ) ( ) ( ) = + + +
2 3
( )
k k
S x y =
1 1 1
' '
1 1 1
'' ''
1 1 1
( ) ( )
( ) ( )
( ) ( )
k
k
k k k k
k k k
k k k
q x q x
q x q x
q x q x
+ + +
+ + +
+ + +
=
=
=
h
a a
h
a a
k
k k
k
k k +


=
1
1
1
3 3
( ) ( )
h c h h c h c
k k k k k k k +
+ + +
1 1 1 1
2 ( ) =
a f x k n
k k
= = ( ), , ,..., 01
b
h
a a
h
c c k n
k
k
k k
k
k k
= + =
+ +
1
3
2 0 1 1
1 1
( ) ( ), , ,...,
d c c h k n
k k k k
= =
+
( ) / ( ), , ,
1
3 01 1





h x x
k k k
=
+1
q x a b x x c x x d x x
k k k k k k k k
( ) ( ) ( ) ( ) = + + +
2 3
n-1 ecuaciones y n+1 incgnitas
Condiciones Naturales
Teorema 1
Sea f(x) una funcin definida en [x
0
,x
n
]. Entonces
existe un nico s(x) spline interpolante cbico
para f(x) en [x
0
,x
n
] tal que
s(x
0
) = 0 y s(x
n
) = 0.
c
n
= s(x
n
)/2 = 0
s(x
0
) = 2c
0
= 0 c
0
= 0.

Matriz del sistema
M
h h h
h h h h
h h h
h h h
h h h h
h h h
n n n
n n n n
n n n
=
+
+
+
+
+
+

(
(



2 0 0 0 0
2 0 0 0
0 2 0 0 0
0 0 0 2 0
0 0 0 2
0 0 0 0 2
0 1 1
1 1 2 2
2 2 3
4 3 3
3 3 2 2
2 2 1
( )
( )
( )
( )
( )
( )








(
(
(
(
(
(
(
(
(
p
h
a a
h
a a
h
a a
h
a a
n
n n
n
n n
=

(
(
(
(
(
(
(
(
(
(
(

3 3
3 3
1
2 1
0
1 0
1
1
2
1 2
( ) ( )
( ) ( )

Trminos independientes
Ejemplo de la temperatura
5 10 15 20
6
8
10
12
14
16
18
20
22
Hora
G
r
a
d
o
s

Spline cbico
5 10 15 20
6
8
10
12
14
16
18
20
22
Hora
G
r
a
d
o
s

Polinomio interpolador
Condiciones sobre la derivada
Teorema 2
Sea f(x) una funcin definida en [x
0
,x
n
]. Entonces existe un nico
s(x) spline cbico interpolante para f(x) en [x
0
,x
n
].tal que
s(x
0
) = f(x
0
) y s(x
n
) = f(x
n
).
2
3
3
0 0 0 1
0
1 0 0
h c h c
h
a a f x + = ( ) ' ( )
h c h c f x
h
a a
n n n n n
n
n n

+ =
1 1 1
1
1
2 3
3
' ( ) ( ).



Matriz del sistema
M
h h
h h h h
h h h h
h h h
h h h
h h h h
h h
n n n
n n n n
n n
=
+
+
+
+
+

(
(
(
(
(



2 0 0 0 0 0
2 0 0 0 0
0 2 0 0 0
0 0 2 0 0 0
0 0 0 0 2 0
0 0 0 0 2
0 0 0 0 0 2
0 0
0 0 1 1
1 1 2 2
2 2 3
3 2 2
2 2 1 1
1 1

( )
( )
( )
( )
( )
(
(
(
(
(
(
Trminos independientes
p
h
a a f x
h
a a
h
a a
h
a a
h
a a
f x
h
a a
n
n n
n
n n
n
n
n n
=

(
(
(
(
(
(
(
(
(
(

3
3
3 3
3 3
3
3
0
1 0 0
1
2 1
0
1 0
1
1
2
1 2
1
1
( ) ' ( )
( ) ( )
( ) ( )
' ( ) ( )

Splines Cbicos
Interpolacin segmentaria con
MATLAB
Interpolacin segmentaria cbica
ps = spline(x,y) % Devuelve el Spline, no los
coeficientes
[x,s] = unmkpp(ps) % Devuelve los coeficientes
ps = mkpp(x,s)
syy = spline(x,y,xx) = ppval(ps,xx)
Interpolacin segmentaria lineal
lyy = interp1(x,y,xx)
-1 0 1
0
0.5
1
Spline Natural
-1 0 1
0
0.5
1
Spline Derivada
-1 0 1
0
0.5
1
Interpolacin Lineal
-1 0 1
0
0.5
1
Spline de MATLAB