w

w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 5

Captulo 2: INTERPOLACIN.

Captulo 2.1. INTRODUCCIN.

Supongamos que partimos de una tabla de valores de la funcin seno. Por ejemplo:

x ( ) ( ) x sen

60 866025 , 0
65 906308 , 0
70 939693 , 0
75 965926 , 0
80 984808 , 0


Cmo podramos calcular, por ejemplo, el valor de ( ) 66 sen ?

Para determinar dicho valor, de una manera aproximada, podramos construir un
polinomio ( ) x P tal que el polinomio coincida con la funcin ( ) ( ) x sen x f = en los valores de la
tabla suministrada, tal que:
( )
( )
( ) ( )
i i
x sen x P
P
P
=
=
=
...
906308 , 0 65
866025 , 0 60

As, ( ) ( ) 66 66 P sen

Pero, qu tipo de polinomio verificara est condicin?

Veamos unos ejemplos ms sencillos.

En primer lugar, suministramos una tabla con dos valores:


0
x ( )
0 0
y x f =
1
x ( )
1 1
y x f =

Representando grficamente los puntos ( )
0 0
, y x y ( )
1 1
, y x vemos como definen una
recta. En este caso el polinomio que buscamos sera de la forma ( ) x c c x P
1 0
+ = , esto es, un
polinomio de grado menor o igual a uno.
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 6



Ahora, suministraremos una tabla con tres valores:


0
x ( )
0 0
y x f =
1
x ( )
1 1
y x f =
2
x ( )
2 2
y x f =

Representando grficamente los puntos ( )
0 0
, y x , ( )
1 1
, y x y ( )
2 2
, y x vemos como definen
una parbola. En este caso el polinomio que buscamos sera de la
forma ( )
2
2 1 0
x c x c c x P + + = , esto es, un polinomio de grado menor o igual a dos.


x
y
y
0

x
2
x
0

y
2

(x
0
,y
0
)

(x
2
,y
2
)
x
1

y
1
(x
1
,y
1
)

x
y
y
0

x
1

x
0

y
1

(x
0
,y
0
)

(x
1
,y
1
)
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 7

Observamos, pues, que para una muestra de 1 + n trminos es necesario un polinomio
de grado menor o igual a n .

Captulo 2.1.1. EL PROBLEMA DE LA INTERPOLACIN:

Sean
n
x x x x ,..., , ,
2 1 0
, 1 + n distintos pertenecientes a R y sea f una funcin real
definida en el intervalo [ ] b a I , = , con I x x x x
n
,..., , ,
2 1 0
. Queremos construir un
polinomio ( ) x P de grado menor o igual a n que interpola a f en los puntos
n
x x x x ,..., , ,
2 1 0
, es
decir, que ( ) ( )
0 0
x f x P = , ( ) ( )
1 1
x f x P = , ( ) ( )
2 2
x f x P = , , ( ) ( )
n n
x f x P =
( 1 + n condiciones).

Demostraremos que tal polinomio existe y es nico, y lo denominaremos polinomio
interpolante.

En primer lugar demostraremos la unicidad y a continuacin la existencia,
construyndolo.

Captulo 2.1.2. TEOREMA DE EXISTENCIA Y UNICIDAD:

Vamos a verificar en primer lugar que existe, a lo sumo, un polinomio de grado menor o
igual a n que interpola a f en los 1 + n puntos distintos
n
x x x x ,..., , ,
2 1 0
:

Suponemos la existencia de dos polinomios ( ) x P y ( ) x Q de grado menor o igual a n que
interpolan a f en los 1 + n puntos distintos
n
x x x x ,..., , ,
2 1 0
, es decir,
( ) ( ) ( ) ( ) [ ] n i x f x Q x f x P
i i i i
, 0 , y = = .

Definimos ahora un polinomio ( ) ( ) ( ) x Q x P x R = de grado menor o igual a n . De este
modo tenemos:
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) 0
...
0
0
0
2 2 2 2 2
1 1 1 1 1
0 0 0 0 0
= = =
= = =
= = =
= = =
n n n n n
x f x f x Q x P x R
x f x f x Q x P x R
x f x f x Q x P x R
x f x f x Q x P x R


Es decir, el polinomio ( ) x R es nulo en los 1 + n puntos
n
x x x x ,..., , ,
2 1 0
, o lo que es igual,
tiene 1 + n races. Sin embargo, de acuerdo con el teorema fundamental del lgebra, sabemos
que todo polinomio de grado menor o igual a n tiene n races (reales o imaginarias, iguales o
repetidas). Cmo podemos justificar esta aparente contradiccin? La nica posible solucin
es que ( ) x R sea el polinomio nulo, ( ) 0 = x R . En tal caso, tenemos ( ) ( ) ( ) 0 = = x Q x P x R y
por tanto, ( ) ( ) x Q x P = y queda demostrada la unicidad del polinomio interpolante.

Probaremos ahora la existencia de dicho polinomio mediante su construccin. Para
llevarla a cabo vamos a aplicar dos mtodos: la Frmula de Lagrange y la Frmula de Newton.


w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 8

Captulo 2.2. POLINOMIO INTERPOLANTE. FRMULAS DE LAGRANGE Y
NEWTON

Captulo 2.2.1. FRMULA DE LAGRANGE:

Si
n
x x x x ,..., , ,
2 1 0
son 1 + n puntos distintos, el polinomio ( ) x P
n
, de grado menor o igual
a n , que verifica ( ) ( ) [ ] n i x f x P
i i n
, 0 , = (polinomio de interpolacin) viene dado en la
frmula de Lagrange por la siguiente expresin:
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )

=
= + + + + =
n
i
i i n n n
x l x f x l x f x l x f x l x f x l x f x P
0
. . . 2 2 . 1 1 . 0 0
... ,
siendo ( ) ( ) ( ) ( ) x l x l x l x l
n
,..., , ,
2 1 0
un polinomio de grado n que se anula en todos los
puntos ( ) [ ] x f x, excepto en uno en que su valor es la unidad.
( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( )
( )
( )
( )
( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( )
( )
( )
( )
( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( )
( )
( )
( )
( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( )
( )
( )
( )

=
=
=
=


=

=
=
=
=


=

=
=
=
=


=

=
=
=
=


=

0
...
1
0
0
...
...
...
0
...
1
0
0
...
...
0
...
0
1
0
...
...
0
...
0
0
1
...
...
2
2 2
1 2
0 2
1 2 1 0
1 2 1 0
2
2 2
1 2
0 2
2 3 2 1 2 0 2
3 1 0
2
1
2 1
1 1
0 1
1 3 1 2 1 0 1
3 2 0
1
1
2 0
1 0
0 0
0 3 0 2 0 1 0
3 2 1
0
n
n n n n n
n
n
n
n
n
n
n
n
n
n
n
x l
x l
x l
x l
x x x x x x x x
x x x x x x x x
x l
x l
x l
x l
x l
x x x x x x x x
x x x x x x x x
x l
x l
x l
x l
x l
x x x x x x x x
x x x x x x x x
x l
x l
x l
x l
x l
x x x x x x x x
x x x x x x x x
x l

( )
( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( )

=
+
+

=


=
n
i k
k
k i
k
n i i i i i i i i
n i i
i
x x
x x
x x x x x x x x x x x x
x x x x x x x x x x x x
x l
0
1 1 2 1 0
1 1 2 1 0
... ...
... ...

w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 9

Vamos a verificar que ( ) x P
n
es el polinomio de interpolacin:
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( )
0
2 1 0
. . 2 2 . 1 1 . 0 0 0
0 ... 0 0 1
...
x f
x f x f x f x f
x l x f x l x f x l x f x l x f x P
n
n n n
=
= + + + + =
= + + + + =


La identidad se verifica sucesivamente para los distintos valores
n
x x x x ,..., , ,
2 1 0
.

Ejercicio:
Se dispone de la siguiente tabla de valores:



2 23
1 7
0 1
1 1
3 17


Calcular el polinomio de interpolacin aplicando la frmula de Lagrange.
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( )
120
1 2 1
17
12
3 1 2
1
6
3 1 1 2
1
10
3 1 2
7
30
3 1 1
23
3 1 0 1 2
. 4 . 3 . 2 . 1 . 0 4
+ +
+

+ +
+
+
+ +

+
=
= + + + + =
x x x x x x x x
x x x x x x x x x x x x
x l f x l f x l f x l f x l f x P

Un problema que tiene la frmula de Lagrange es que si queremos aadir un nuevo
punto,
1 + n
x , el polinomio se tiene que volver a escribir desde cero ya que los
polinomios ( ) ( ) ( ) ( ) x l x l x l x l
n
,..., , ,
2 1 0
y ahora, ( ) x l
n 1 +
se tienen que volver a calcular. Este
problema se solventa con la frmula de Newton.

Captulo 2.2.2. FRMULA DE NEWTON:

La idea bsica es construir el polinomio en pasos sucesivos: primero construiremos el
polinomio ( ) x P
0
, de grado menor o igual a cero, que coincide con la funcin en el punto
0
x ;
luego construiremos el polinomio ( ) x P
1
, de grado menor o igual a uno, que coincide con la
funcin en los puntos
1 0
, x x , as hasta llegar al polinomio ( ) x P
n
, de grado menor o igual
a n que coincide con la funcin en los puntos
n
x x x x ,..., , ,
2 1 0
. De este modo, cada polinomio
viene definido a partir del anterior, es decir, viene dado por una relacin de recurrencia.
( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
n n n n n n n
x f x P x f x P x f x P x f x P x P
x f x P x f x P x f x P x P
x f x P x f x P x P
x f x P x P
= = = =
= = =
= =
=
,..., , ,
...
, ,
,
2 2 1 1 0 0
2 2 2 1 1 2 0 0 2 2
1 1 1 0 0 1 1
0 0 0 0

w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 10

Vamos a calcular estos polinomios:
( ) ( )
0 0 0
c x f x P = =

( ) ( ) ( )
( ) ( ) ( )
( ) ( )
( ) ( ) ( )
( ) ( ) ( ) ( ) ( )
( ) ( )
( ) ( )
( ) ( )
0 1 0 1
0 1
0 1
1
1 1 1
0 1 0 1 1 0 1 1
0 1 1 0 1
0 1
0 1 0 0 1
1 0 1

0
x x c c x P
x x
x f x f
c
x f x P
x x c x f x Q x f x P
x x c x Q x Q
x f x P
x Q x f x P
x Q x P x P
+ =

=
)
`

=
+ = + =
= =
)
`

=
+ =
+ =


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

0
0
x x x x c x x c c x P
x x x x
x x c c x f
c
x f x P
x x x x c x x c c x Q x P x P
x x x x c x Q
x Q
x f x P
x Q x f x Q x P x P
x Q
x f x P
x Q x f x Q x P x P
x Q x P x P
+ + =


=

)
`

=
+ + = + =
=

=
)
`

=
+ = + =
=
)
`

=
+ = + =
+ =

As tenemos:
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
1 2 1 0 1 0 2 0 1 0
... ...

+ + + + =
n n n
x x x x x x x x c x x x x c x x c c x P

De este modo solventamos el inconveniente que tenamos al aadir ms puntos. Por
ejemplo, si aadimos un punto
1 + n
x :
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
n n n
x x x x x x x x c x x x x c x x c c x P + + + + =
+ +
... ...
2 1 0 1 1 0 2 0 1 0 1

Pero, cmo calculamos los coeficientes
n
c c c c ,..., , ,
2 1 0
?

Captulo 2.2.3. DIFERENCIAS DIVIDIDAS DE UNA FUNCIN:

Dados 1 + n puntos distintos
n
x x x x ,..., , ,
2 1 0
y una funcin f definida en tales puntos, se
llama diferencia dividida de la funcin f en los puntos
n
x x x x ,..., , ,
2 1 0
y se
representa [ ]
n
x x x x f ,..., , ,
2 1 0
, al coeficiente de
n
x en el desarrollo del correspondiente
polinomio interpolador.

As, [ ] [ ] [ ] [ ]
n n
c x f c x f c x f c x f = = = = ,..., , ,
2 2 1 1 0 0
, siendo n el orden de la diferencia
dividida.

w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 11

El valor de una diferencia dividida es independiente del orden en que se escriban sus
argumentos.

El polinomio de Newton, usando diferencias divididas, vendr dado por:
( ) [ ] [ ]( ) [ ]( ) ( )
[ ]( ) ( ) ( ) ( )
1 2 1 0
1 0 2 0 1 0
...
...

+
+ + + + =
n n
n
x x x x x x x x x f
x x x x x f x x x f x f x P


Verificndose que [ ]
[ ] [ ]
0
1 2 1 0 3 2 1
2 1 0
,..., , , ,..., , ,
,..., , ,
x x
x x x x f x x x x f
x x x x f
n
n n
n

=



De esta manera, las diferencias divididas de cualquier orden se definen a partir de las de
rdenes precedentes y su clculo se sintetiza en forma de tabla triangular:


0
x [ ]
0
x f


[ ]
[ ] [ ]
0 1
0 1
1 0
,
x x
x f x f
x x f

=

1
x [ ]
1
x f


[ ]
[ ] [ ]
1 2
1 2
2 1
,
x x
x f x f
x x f

=

2
x [ ]
2
x f


[ ]
[ ] [ ]
2 3
2 3
3 2
,
x x
x f x f
x x f

=

3
x [ ]
3
x f








[ ]
[ ] [ ]
0 2
1 0 2 1
2 1 0
, ,
, ,
x x
x x f x x f
x x x f

=


[ ]
[ ] [ ]
0 3
2 1 0 3 2 1
3 2 1 0
, , , ,
, , ,
x x
x x x f x x x f
x x x x f

=

[ ]
[ ] [ ]
1 3
2 1 3 2
3 2 1
, ,
, ,
x x
x x f x x f
x x x f

=





w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 12

Ejercicio:
Se dispone de la siguiente tabla de valores:



2 23
1 7
0 1
1 1
3 17


Calcular el polinomio de interpolacin aplicando la frmula de Newton, empleando
diferencias divididas.
( ) ( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
3 2 1 0 4
2 1 0 3 1 0 2 0 1 0 4


x x x x x x x x c
x x x x x x c x x x x c x x c c x P
+
+ + + + =



2

23



[ ] 16 ,
1 0
= x x f


1 7

[ ] 5 , ,
2 1 0
= x x x f



[ ] 6 ,
2 1
= x x f

[ ] 1 , , ,
3 2 1 0
= x x x x f


0 1

[ ] 2 , ,
3 2 1
= x x x f


[ ] 0 , , , ,
4 3 2 1 0
= x x x x x f


[ ] 2 ,
3 2
= x x f


[ ] 1 , , ,
4 3 2 1
= x x x x f


1 1

[ ] 2 , ,
4 3 2
= x x x f



[ ] 8 ,
4 3
= x x f


3 17


( ) ( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )
1 3 2
1 1 2 0 1 2 1 2 5 2 16 23


2 3
3 2 1 0 4
2 1 0 3 1 0 2 0 1 0 4
+ =
= + + + + + + + + + + =
= +
+ + + + =
x x x
x x x x x x x x x x
x x x x x x x x c
x x x x x x c x x x x c x x c c x P

w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 13

Captulo 2.3. ERROR EN LA INTERPOLACIN.

Sea la funcin [ ] R b a f , : , 1 + n veces derivable en( ) b a, y sea ( ) x P
n
el
polinomio de grado menor o igual a n que interpola a la funcin en los 1 + n puntos distintos,
[ ] b a x x x x
n
, ,..., , ,
2 1 0
tales que b x x x x a
n
= < < < < = ...
2 1 0
. Si [ ] b a x , , se verifica que
existe un punto
x
c tal que el error en x viene dado por la expresin:
( ) ( ) ( )
( )
( )
( ) ( )( )( ) ( )
( )
( )
( ) ( )

=
+
+

+
=
=
+
= =
n
i
i
n
n
n
n n
x x c f
n
x x x x x x x x c f
n
x P x f x e
0
1
2 1 0
1

! 1
1
...
! 1
1




Si
( )
( ) ( ) b a x M x f M
n
, ,
1

+
, entonces ( )
( )
( )

n
i
i n
x x
n
M
x e
0

! 1
.

En este caso tomaremos
( )
( )
( ) x f M
n
b a x
1
,
sup
+

= .
f(x) P(x)
x
0

x
y
x
1

x
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 14

Ejercicio:
Se dispone de los siguientes valores de la funcin logaritmo neperiano:

x
x ln

9 1972 . 2
5 . 9 2513 . 2
11 3979 . 2


a) Aproximar el valor de 2 . 9 ln utilizando interpolacin lineal. Dar una cota del error
cometido.
Vamos a construir el polinomio de interpolacin lineal aplicando la frmula de
Newton empleando diferencias divididas.
( ) ( )
0 1 0 1
x x c c x P + =
La tabla de diferencias divididas ser:

9 1972 . 2


[ ] 1082 . 0 ,
1 0
= x x f

5 . 9 2513 . 2

[ ] 0052 . 0 , ,
2 1 0
= x x x f

[ ] 0977 . 0 ,
2 1
= x x f

11 3979 . 2

As, el polinomio de interpolacin vendr dado por:
( ) ( ) ( ) 9 1082 . 0 1972 . 2
0 1 0 1
+ = + = x x x c c x P
Evaluando el valor de dicho polinomio en 2 . 9 = x obtendremos el valor aproximado
de 2 . 9 ln
( ) ( ) 2 . 9 ln 2188 . 2 9 2 . 9 1082 . 0 1972 . 2 2 . 9
1
= + = P
A continuacin vamos a evaluar el error cometido en la aproximacin:
( ) ( ) ( ) ( ) ( ) ( )
( ) ( )
( )
( )
81
03 . 0
9
03 . 0
5 . 9
1 1
9
1
5 . 9
1 1
9
1
9 5 . 9 5 . 9 , 9
03 . 0 1
03 . 0 03 . 0
2
06 . 0
5 . 9 2 . 9 9 2 . 9
2
1
2 . 9
! 2
1
2
2 2 2
2 2
' ' ' '
' '
0
' '
=
< < < < < <
= = = =
= = =

=
x e
c c
c c
c c
c f c f
c f x c f x e
n
n
i
i n

b) Aproximar el valor de 2 . 9 ln utilizando interpolacin cuadrtica. Dar una cota del error
cometido.

w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 15

Captulo 2.3.1. ERROR EN LA INTERPOLACIN LINEAL:

Sea una funcin f , dos veces derivable en ( )
1 0
, x x , y sea ( ) x P
1
el polinomio que
interpola a f en esos dos puntos. Sabemos que si [ ]
1 0
, x x x ,
( ) ( ) ( ) ( ) ( ) ( )( )
1 0
' '
1 1 0

! 2
1
, x x x x c f x P x f x e x x c
n
= = .

Vamos a hallar, a partir de esta expresin, una cota para el error (en valor absoluto) en la
interpolacin lineal, vlida para cualquier punto [ ]
1 0
, x x x .

( )( ) ( )( )
( ) ( )
( ) ( )
2 0
2 0 0
0 0 0 1 1 0
1 1 , 0

h t t
h
x x
t
h
h
h x x
h
x x
h x x x x x x h x x x x
=

= =
=

= = = =




( ) ( ) ( ) c f h t t x e
' ' 2
1
1
2
1
=
( ) ( ) ( )
( )
( )
4
1
2
1
(mximo)
2
1
0 2 1
0
2 1
1 1
'
'
2
=
|
.
|

\
|
= =

=
=
= = =
g
t t
t g
t t g
t t t t t t t g

( ) ( ) c f h x e
' ' 2
1

4
1

2
1

( ) ( ) [ ]
1 0
' '
2
1
, ,
8
x x x c f
h
x e
x
0
x
1

h
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 16

Captulo 2.3.2. ERROR EN LA INTERPOLACIN CUADRTICA:

Sea una funcin f , tres veces derivable en( )
2 0
, x x , y sea ( ) x P
2
el polinomio que
interpola a f en los puntos
2 1 0
, , x x x , tal que
2 1 0
x x x < < . Sabemos que si [ ]
2 0
, x x x ,
( ) ( ) ( ) ( ) ( ) ( )( )( )
2 1 0
' ' '
1 2 0

! 3
1
, x x x x x x c f x P x f x e x x c
n
= = .

Vamos a hallar, a partir de esta expresin, una cota para el error (en valor absoluto) en la
interpolacin cuadrtica, vlida para cualquier punto [ ]
1 0
, x x x y supuesto, adems, que las
abscisas de los nodos son equidistantes.

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

2
h t t t
h
x x
t h
h
h x x
h
h x x
h
x x
h x x h x x x x x x h x x x x x x
=

= =

=
= = = =




( ) ( ) ( ) ( ) c f h t t t x e
' ' ' 3
2
2 1
! 3
1
=
( ) ( ) ( )
( ) ( ) [ ]
( ) ( ) ( ]
[ ]
( ]
( )
( )
( )
( )
( )
( ) ( )
9
3 2
2 1 max
9
3 2
3
3

3
2
3
3

3
1
1
3
3
1
3
3

3
3
1 2
3
3
1 1
3
3
1
3
3
1
3
3
1
(mximo)
3
3
1 0 2 6 3
0
2 , 1 2 3 3
1 , 0 2 6 3
2 , 1 2 3
1 , 0 2 3
2 , 1 2 1
1 , 0 2 1
2 1
2 , 0
2
'
2
2
'
2 3
2 3
=
= =
|
.
|

\
|
=
=
|
|
.
|

\
|

|
|
.
|

\
|
+ =
|
|
.
|

\
|
+
|
|
.
|

\
|
+
|
|
.
|

\
|
+ =
|
|
.
|

\
|
+
= = +

+
+
=

+
+
=

= =

t t t
g
t t t
t g
t si t t
t si t t
t g
t si t t t
t si t t t
t si t t t
t si t t t
t t t t g
t

x
0
x
1

h
x
2

2h
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 17

As, ( ) ( ) ( ) ( ) c f h t t t x e
' ' ' 3
2
2 1
! 3
1
=
( ) ( ) [ ]
2 0
' ' '
3
2
, ,
3 9
x x x c f
h
x e

Ejercicio:
Acotar el error de interpolacin para la funcin ( )
x
e x f = en el intervalo[ ] 1 , 0 , cuando se
interpola por un polinomio de grado 2 en los nodos 1 ,
2
1
, 0 .
( ) ( )
( )
( )
3 72
1 1 0 1 , 0

2
1

3 9
1

3 9
2
1 0
3
' ' '
3
2
e
x e
e e e e e c c
e c f
h
x e
c c
c

< < < < < <

+

+
= =
n
i
i
n
i
i
n
n n
x x
n
x x c f
n
x P x f x e
0 0
1

! 1
1

! 1
1


Por otra parte:



( ) ( ) ( ) ( ) a b a b a b a b x x x x x x x x
n
... ...
2 1 0


As, se verifica que:
( )
( )
( )! 1
1
+

+
n
a b
x e
n
n

x
0
=a x
1
x
2
x x
n
=b
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 18

Qu ocurre cuando aumentamos los nodos?
( )
( )
0
! 1
lim
1
=
+

+

n
a b
n
n


Vemos, pues, que a medida que incrementamos el nmero de nodos el error cometido
en la interpolacin tiende a cero.

Pero en general, no se consigue disminuir el error incrementando el nmero de nodos y
por lo tanto, el grado del polinomio.

Para algunas funciones, la cota del error aumenta a medida que aumenta n . Un elevado
nmero de puntos puede dar lugar a polinomios con muchas oscilaciones que aproximan mal a
la funcin que se interpola, especialmente en los extremos del intervalo.

Por ejemplo:
La funcin ( )
2
1
1
x
x f
+
= en[ ] 5 , 5
( ) x P
4
en 5 , 3 , 0 , 3 , 5



w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 19

( ) x P
10
en 5 , 4 ,..., 0 ,..., 4 , 5



( ) x P
20
en 5 , 5 . 4 ,..., 0 ,..., 5 . 4 , 5



Para solucionar estos problemas introduciremos el concepto de la interpolacin a trozos.


w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 20

Captulo 2.4. INTERPOLACIN POLINMICA A TROZOS.

Una buena forma de obtener, por interpolacin, una funcin es la interpolacin a trozos o
interpolacin segmentaria que consiste en dividir el intervalo[ ] b a, en subintervalos lo
suficientemente pequeos y aproximar la funcin en esos subintervalos mediante polinomios
de grado pequeo.

Captulo 2.4.1. INTERPOLACIN LINEAL A TROZOS:

Sea f dos veces derivable en( ) b a, y dados 1 + n puntos b x x x x a
n
= < < < < = ...
2 1 0
,
se aproxima la funcin f en cada intervalo de la forma[ ] 1 ,..., 0 , ,
1
=
+
n i x x
i i
mediante la recta
que une los puntos ( ) ( )
i i
x f x , y ( ) ( )
1 1
,
+ + i i
x f x .



x
0

x
y
x
1
x
2
x
3

f(x
0
)
f(x
2
)
f(x
1
)
f(x
3
)
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 21

Captulo [Link]. ERROR EN LA INTERPOLACIN LINEAL A TROZOS:

Queremos calcular una cota para el error en la interpolacin lineal a trozos independiente
del punto escogido.
Si [ ]
1 0
, x x x , ( ) ( ) ( ) ( )
1 0 1 1
' '
2
0 1 1
, ,
8
1
x x c c f x x x e
Si [ ]
2 1
, x x x , ( ) ( ) ( ) ( )
2 1 1 2
' '
2
1 2 1
, ,
8
1
x x c c f x x x e

Si [ ]
n n
x x x ,
1
, ( ) ( ) ( ) ( )
n n n n n
x x c c f x x x e , ,
8
1
1 1
' '
2
1 1

As, si [ ]
n
x x x ,
0
, ( ) ( ) ( ) { } [ ] M x x x P x f x e
i i
n i
, max
8
1
2
1
1 0
1 1 +

= ,
siendo
( )
( ) ( ) b a x M x f
n
, ,
1

+
.

Si las abscisas de los nodos son equidistantes, entonces, el error cometido en la
interpolacin lineal a trozos vendr dado por:
( ) M h x e
8
1
2
1


Ejercicio:
Sea la funcin ( )
3
x x f = en el intervalo[ ] 1 , 1 .
a) Se desea generar una tabla de valores de f con abscisas equiespaciadas de forma
que la interpolacin lineal a trozos entre dos puntos consecutivos tenga un error
menor que
4
10 3

. Determinar el nmero mnimo de puntos.
x
( ) x f
1

h + 1

h 2 1+

...

h n 1+


El error cometido viene dado por:
( ) M h x e
8
1
2
1

Calculamos el valor de M :
( )
( ) ( ) ( )
( ) [ ] 1 , 1 , 6
6 3
' '
' ' 2 ' 3
' '

= = =

c c f
x x f x x f x x f
M c f

w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 22

As:
( )
( )
( )( )
2
max
2 2 4 2
4 2 4 2 4 2
4
1
2 2
1
10 2
0 10 2 10 2 0 10 4
10 4 10
4
1
10 3
4
3
10 3

4
3
6
8
1

=
+

=
h
h h h
h h h
x e
h h x e

Por otra parte sabemos que:
( )
n n n
a b
h
2 1 1
=

=

=
Y por tanto:
100
1
10
2
10 2
2
10 2
2 2
2
= = =

=
=

n
n n
n
h
h

El mnimo nmero de puntos es 101 1= + n puntos.

b) Se desea utilizar un nico polinomio interpolante (con nodos equidistantes) en todo
el intervalo, manteniendo la cota del error
4
10 3

. Determinar el nmero mnimo de
puntos y el polinomio de interpolacin.
Anteriormente hemos probado que con interpolacin lineal es necesario
emplear101puntos para que el error cometido sea menor que
4
10 3

.
Vamos a probar ahora con la interpolacin cuadrtica.
x
( ) x f
1

0

1


El error cometido viene dado por:
( ) M h x e
3 9
1
3
2

Calculamos el valor de M :
( )
( ) ( ) ( ) ( )
( ) [ ] 1 , 1 , 6
6 6 3
' ' '
' ' ' ' ' 2 ' 3
' ' '
=
= = = =

c c f
x f x x f x x f x x f
M c f

As:
( )
( )
? 10 3
3 3
2

10 3
6
3 9
1

3 9
1
4
4
1
1

=
x e
M x e

Con tres nodos equidistantes no se puede asegurar que el error sea menor
que
4
10 3

.
w
w
w
.
.
c
o
m
M
a
t
e
m
a
t
i
c
a
1
Universidad de Oviedo
MTODOS NUMRICOS
Escuela Universitaria de Ingeniera Tcnica en Informtica de Oviedo (E.U.I.T.I.O)

Alberto Surez Lpez
Pgina 23

Vamos a probar ahora la interpolacin cbica.
x
( ) x f
1

3 1

3 1

1


Tenemos que calcular un polinomio de grado menor o igual a 3 y en el que se
verifique el valor de la funcin en los nodos anteriores. Deducimos, pues, que el
polinomio de interpolacin coincide con la propia funcin, luego el error es cero.
As, el mnimo nmero de puntos necesarios es cuatro y el polinomio de
interpolacin es ( ) ( ) x f x x P = =
3
3
.