Ejemplo 1

Rendimiento de un proceso productivo en funcin de la temperatura

En una planta qumica se sintetiza un producto que es utilizado posteriormente como conservante de productos enlatados. El rendimiento del proceso depende de la temperatura. Se dispone de los siguientes datos T ( C ) R ( %) 150 160 170 180 190 200 210 35.5 37.8 43.6 45.7 47.3 50.1 51.2

Se considera un rendimiento ptimo el que va de 38.5 a 45, por lo que la planta trabaja a 175 C . Si la temperatura de trabajo cae a 162 C por una avera, ser el proceso satisfactorio hasta que sea reparada?

Ejemplo 2
En una planta se bombea esencia de trementina, 60 C , desde la base de una columna de fraccionamiento hasta un gran tanque de almacenamiento descubierto. La columna opera a 1,29 atmsferas . En la siguiente tabla se representan los datos relativos los litros por hora que puede bombear la bomba en funcin de la potencia en watios a la que es necesario que trabaje: Q (l/h) N (w) 500 700 900 1100 1300 1500 1700 1900 365 361.6 370.64 379.68 384.46 395.5 395.95 397

Se desea saber si la bomba ser capaz de impulsar un caudal de 1000 l/h de trementina hasta el tanque de almacenamiento trabajando a un mximo de 373 w.

Ejemplo 3
El pentxido de dinitrgeno gaseoso puro reacciona en un reactor intermitente segn la reaccin estequiomtrica N2 O5 2N 2 O4 +O2 Calculamos la concentracinde pentxido de dinitrgeno existente en ciertos instantes, obteniendo los siguientes datos: T (s) C 0 200 400 650 1100 1900 2300 5.5 5.04 4.36 3.45 2.37 1.32 0.71

Si lo tenemos en el reactor un tiempo mximo de 35 minutos ( 2100 segundos), cul es la concentracin de pentxido de dinitrgeno que queda sin reaccionar?

Interpolacin
Concepto de interpolacin
Supongamos que hay dos magnitudes x e y de los que se conocen n + 1 valores relacionados {(x0 , y0 ) , (x1 , y1 ) , , (xn , yn )} , por ejemplo, datos obtenidos en una experimentacin. Con la condicin xi 6= xj si i 6= j . Nos planteamos si existe una funcin p tal que p (xk ) = yk k = 0, , n (1)

Este tipo de problemas suele darse cuando tenemos datos obtenidos por experimentacin y sabemos que hay una funcin f que rige el proceso pero que desconocemos y queremos trabajar con una funcin alternativa p que represente bien a esos datos de la muestra. Si f rige el proceso entonces f (xk ) = yk luego exigiremos a la funcin p ese mismo requisito, esto nos proporciona condiciones que imponer a p con las que trataremos de obtenerla y una vez conseguido nos permitira conocer o predecir qu habra pasado en otros x en los que no se ha experimentado. Supongamos que existe la funcin f tal que f (xk ) = yk , k = 0, , m . Caben varias preguntas: i) La funcin p que interpola los datos dados de qu tipo ha de ser? polinmica, trigonomtrica, racional,...? La respuesta vendr dada por los datos yk . ) Si se observa que los datos presentan periodicidad entonces buscaremos a p dentro de la funciones trigonomtricas. ) Si los datos presentan asntotas entonces p debera ser una funcin racional. ) Si los yk presentan un comportamiento polinomial, entonces p se escogera de tipo polinmico. Nos centraremos en cmo resolver este caso.

es decir, queremos una funcin cuya grca "pase"por los puntos del plano dados. Si p verica (1) diremos que p interpola los datos dados p es una funcin de interpolacin para los datos (xk , yk ) , k = 0, , n.

ii) Una vez escogido el tipo de funcin habr que responder dos cuestiones ,existe p del tipo escogido que interpole los datos dados? Y si existe , es nica? iii) Es la funcin polinmica escogida una buena aproximacin de la funcin original f en los puntos x que no son de la muestra? Nota: entendermos como funcin original la que rige el experimento y de la cual slo sabemos qu pasa en los n+1 puntos de la muestra. Vamos a hacer el estudio contestando a estas cuestiones suponiendo que la funcin p es una funcin polinmica.

Interpolacin polinmica
Planteamiento del problema
Dada una tabla de n + 1 puntos (xk , yk ) con k = 0, , n tales que xi 6= xj si i 6= j . Llamaremos interpolacin polinmica a la determinacin de un polinomio p de grado menor o igual que n tal que p (xk ) = yk , k = 0, , n Si p es de grado menor o igual que n entonces se podr expresar p (x) = a0 + a1 x + a2 x2 + + an xn donde los ai se obtendrn a partir de las condiciones de interpolacin, esto es, p (x0 ) = p (x1 ) = p (x2 ) = . . .
n a0 + a1 x0 + a2 x2 0 + + an x0 = n a0 + a1 x1 + a2 x2 1 + + an x1 = 2 a0 + a1 x2 + a2 x2 + + an xn 2 = . . .

y0 y1 y2 . . .

n p (xm ) = a0 + a1 xm + a2 x2 m + + an xm = ym

Aparece un sistema para las variables a0 , a1 , ..., an que podemos escribir matriciamente n a y 1 x0 x2 x 0 0 0 0 1 x1 x2 xn a1 y1 1 1 1 x2 x2 xn a2 y2 0 2 = . . . . . . . . . . . . . . . . . . . . . n an ym 1 xm x2 m xm La matriz del sistema es cuadrada (n + 1) (n + 1) y habr una nica solucin del problema si, y slo si, n 1 x0 x2 0 x0 1 x1 x2 xn 1 1 Y 2 n det 1 x2 x0 x2 = (xj xi ) 6= 0 xi 6= xj . . . . . . . . . . . . . . i<j . 2 1 xn xn xn n

es decir, la solucin es nica si y slo si todos los puntos de la muestra son distintos.

Construccin de la solucin Mtodo directo

Planteamos el problema tal como se ha descrito en el prrafo anterior y pasamos a resolver el sistema. Cuando hayamos obtenido la solucin (a0 , a1 , ..., an ) pasamos a escribir la funcin polinmica solucin de nuestro problema p (x) = a0 + a1 x + a2 x2 + + an xn

Ventaja del mtodo: la resolucin del problema de interpolacin pasa por resolver un sistema que es un procedimiento ya conocido . Inconvenientes: Si aparecen nuevos datos de la experimentacin, la solucin p de grado n que tengamos para los datos previos no es aprovechable. Hay que rehacer todos los clculos para la nueva muestra (los datos anteriores y los nuevos juntos).

Polinomio de interpolacin : mtodo de Lagrange

Los polinomios de Lagrange
Para cada i, i = 0, 1. , n construiremos un polinomio de grado menor o igual que n, al que llamaremos pi de manera que pi (xi ) = 1 pi (xj ) = 0 si j 6= i Puesto que el polinomio pi se anula en x0 , x1 , ..., xi1 , xi+1 , ..., xn dicho polinomio debe ser pi (x) = (x x0 ) (x x1 ) (x xi1 ) (x xi+1 ) (x xn ) =
n Q n Y j =0 j 6=i

(x xj ) 1 (xi xj )

Adems, como pi (xi ) = 1 = 1 = pi (xi ) = luego cada pi (x) es

j =0 j 6=i n Q n Q

(xi xj ) = = Q n j =0
j 6=i

j =0 j 6=i

(x xj ) = (xi xj )

pi (x) =

j =0 j 6=i

n Y (x xj ) (xi xj ) j =0 j 6=i

Cada pi (x) es el i-simo polionomio de Lagrange para los puntos x0 , x1 , ..., xn .

Un ejemplo
Para el caso n = 2 ,tenemos x0 , x1 , x2 , los polinomios de Lagrange son p0 (x) = (x x1 ) (x x2 ) (x0 x1 ) (x0 x2 ) (x x0 ) (x x2 ) (x1 x0 ) (x1 x2 ) (x x0 ) (x x1 ) (x2 x0 ) (x2 x1 )

p1 (x) =

p2 (x) =

Teorema Dados n + 1 puntos (xi , yi ), i = 0, 1, ..., n tales que xi 6= xj i 6= j existe un nico polinomio de grado menor o igual que n, p (x), tal que p (xi ) = yi con i = 0, 1, ..., n. Haremos uso del teorema anterior para mostrar cmo se determina p haciendo uso de los polinomios de Lagrange. n (x x ) n n P P Q j Sea el polinomio p (x) = yi pi (x) = yi . ( x x i j) i=0 i=0 j =0 El polinomio p verica p (xj ) =
n P j 6=i

yi pi (xj ) =

i=0

(y0 0 + + yj 1 0) + yj 1 + (yj +1 0 + + yn 0) = yj con j = 0, 1, ..., n es decir, p, es un polinomio de grado menor o igual que n y que satisface las condiciones impuestas. La formacin del polinomio p slo precisa formar los polinomios de Lagrange y escribir una combinacin lineal de ellos donde los coeciente nos vienen dados, los yi . Casos particulares n=1 Interpolacin lineal Polinomio que interpola los datos (x0 , y0 ), (x1 , y1 )

j 1 P i=0

yi pi (xj ) + yj pj (xj ) +

i=j +1

n P

yi pi (xj ) =

Formamos los polinomios de Lagrange p0 (x) = x x1 x0 x1 x x0 x1 x0

p1 (x) = y el polinomio p (x) = y0

x x1 x x0 + y1 x0 x1 x1 x0

n=2 Interpolacin cuadrtica Polinomio que interpola los datos (x0 , y0 ), (x1 , y1 ),(x2 , y2 ) Formamos los polinomios de Lagrange p0 (x) = (x x1 ) (x x2 ) (x0 x1 ) (x0 x2 ) (x x0 ) (x x2 ) (x1 x0 ) (x1 x2 ) (x x0 ) (x x1 ) (x2 x0 ) (x2 x1 )

p1 (x) =

p2 (x) = y el polinomio p (x) = y0

(x x1 ) (x x2 ) (x x0 ) (x x2 ) (x x0 ) (x x1 ) +y1 +y2 (x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 )

n=3 Interpolacin cubica Polinomio que interpola los datos (x0 , y0 ), (x1 , y1 ),(x2 , y2 ) y (x3 , y3 ) Formamos los polinomios de Lagrange p0 (x) = (x x1 ) (x x2 ) (x x3 ) (x x0 ) (x x1 ) (x x3 ) p2 (x) = (x0 x1 ) (x0 x2 ) (x0 x3 ) (x2 x0 ) (x2 x1 ) (x2 x3 ) (x x0 ) (x x2 ) (x x3 ) (x x0 ) (x x1 ) (x x2 ) p3 (x) = (x1 x0 ) (x1 x2 ) (x1 x3 ) (x3 x0 ) (x3 x1 ) (x3 x2 )

p1 (x) =

(x x1 ) (x x2 ) (x x3 ) (x x0 ) (x x2 ) (x x3 ) +y1 + (x0 x1 ) (x0 x2 ) (x0 x3 ) (x1 x0 ) (x1 x2 ) (x1 x3 ) (x x0 ) (x x1 ) (x x3 ) (x x0 ) (x x1 ) (x x2 ) + y3 y2 (x2 x0 ) (x2 x1 ) (x2 x3 ) (x3 x0 ) (x3 x1 ) (x3 x2 ) y el polinomio p (x) = y0

Ejemplo Calcule el nico polinomio de grado menor o igual que tres que interpola a los datos (1, 1), (0, 2), (1, 1) y (2, 2). Solucin: Formamos los polinomios de Lagrange p0 (x) = (x 0) (x 1) (x 2) ; (1 0) (1 1) (1 2) (x (1)) (x 1) (x 2) (0 (1)) (0 1) (0 2) ; p2 (x) = (x (1)) (x 0) (x 2) (1 (1)) (1 0) (1 2) (x (1)) (x 0) (x 1) (2 (1)) (2 0) (2 1)

p1 (x) =

p3 (x) =

el polinomio es p (x) =1

(x 0) (x 1) (x 2) (x (1)) (x 1) (x 2) (x (1)) (x 0) ( +2 1 (1 0) (1 1) (1 2) (0 (1)) (0 1) (0 2) (1 (1)) (1 0) (

+2

1 1 (x (1)) (x 0) (x 1) = x (x 1) (x 2) + (x + 1) (x 1) (x 2) + x (x (2 (1)) (2 0) (2 1) 6 2

1 1 2 + x (x 1) (x + 1) = x3 x2 2x + 2 3 3 3 Se puede utilizar la interpolacin tambin para trabajar con un polinomio en vez de con una funcin dada f .Slo necesitamos saber qu grado queremos manejar y, una vez decidido el grado, evaluar f en n + 1 puntos (xi , f (xi )) i = 0, 1, , n, de esta manera tenemos que el polinomio que interpola a f en los puntos xi i = 0, 1, , n es p (x) =
n X i=0

fi pi (x) , donde fi = f (xi ) con i = 0, 1, , n

Ejemplo: Polinomio de grado menor o igual que 2 que interpola a la funcin f (x) = ex en los puntos x0 = 1, x1 = 0, x2 = 1 El polinomio es p (x) =
i=0 n P

fi pi (x) = f0

(x x1 ) (x x2 ) (x x0 ) (x x2 ) (x x0 ) (x x1 ) + f1 + f2 (x0 x1 ) (x0 x2 ) (x1 x0 ) (x1 x2 ) (x2 x0 ) (x2 x1 ) ( x 0) ( x 1) ( x ( 1)) ( x 1) ( x (1)) (x 0) + e0 + e1 = e1 (1 0) (1 1) (0 (1)) (0 1) (1 1) (1 0) 1 1 1 = e x (x 1) (x 1) (x + 1) + e x (x + 1) = 2 2 1 1 1 1 1 = e + e 1 x2 + e1 + e x + 1 = 2 2 2 2 = (Ch1 1) x2 + Sh1 x + 1

Si ya hemos formado el polinomio de interpolacin pn para los datos (xk , yk ) con k = 0, , n y surge la necesidad de interpolar, adems en un nuevo punto (xn+1 , yn+1 ) todos los clculos anteriores no seran vlidos. Tendramos que formar los correspondientes polinomios de Lagrange de grado n + 1con lo empezaramos de nuevo desde el principio. Para que podamos .aprovechar"el polinomio de interpolacin pn para formar el polinomio que, adems, interpola un nuevo punto (xn+1 , yn+1 ) pasamos a la interpolacin con el mtodo de Newton.

Polinomio de interpolacin : mtodo de Newton

Supongamos que pn (x) es el polinomio de interpolacin para los datos x0 y0 x1 y1 x2 y2 xn yn

y tenemos un nuevo dato (xn+1 , yn+1 ) queremos contruir un nuevo polinomio pn+1 (x) que interpole a los datos anteriores y a (xn+1 , yn+1 ) .

Idea del mtodo de Newton

Formamos el polinomio pn+1 (x) pn+1 (x) = pn (x) + A (x x0 ) (x x1 ) (x xn ) | {z } trmino de grado n+1 pn+1 (xn+1 ) = yn+1 tendremos que pn+1 (xn+1 ) pn (xn+1 ) (xn+1 x0 ) (xn+1 x1 ) (xn+1 xn ) Esta idea se puede llevar a cabo para formar tambin el polinomio pn . A=

El polinomio as construdo es de grado n + 1, interpola a los datos anteriores e imponiendo que se verique la nueva condicin de interpolacin

Diferencias divididas
Trataremos de construir el polinomio de interpolacin pn de los datos (xk , yk ) con k = 0, , n de una muestra sin tener que recurrir a resolver un sistema. Lo escribiremos siguiendo la idea antes planteada pn (x) = A0 +A1 (x x0 )+A2 (x x0 ) (x x1 )+ +An (x x0 ) (x x1 ) (x xn1 )

Denicin i) Llamaremos diferencia dividida de orden cero de la funcin f , y lo notaremos, f [xi ] a yi . ii) Llamaremos diferencia dividida de orden uno de la funcin f , y lo notaremos f [xi+1 ] f [xi ] yi+1 yi = f [xi , xi+1 ], a el cociente f [xi , xi+1 ] = xi+1 xi xi+1 xi iii) Llamaremos diferencia dividida de orden dos de la funcin f ,y lo notaremos f [xi+1 , xi+2 ] f [xi , xi+1 ] , y as sucesivaf [xi , xi+1 , xi+2 ], a f [xi , xi+1 , xi+2 ] = xi+2 xi mente, el siguiente resultado establece cmo formar las diferencias divididas de cada orden.

Ley de recurrencia de las diferencias divididas

Para k 1 tenemos f [x0 , x1 , , xk ] = y f [xi ] = yi i = 0, 1, , n.. Pasamos a aplicar las anteriores deniciones al clculo de pn Si pn (x0 ) = y0 y0 = A0 diferencia dividida de orden 0 en x0 y1 y0 pn (x1 ) = y1 y1 = A0 + A1 (x1 x0 ) A1 = = f [x0 , x1 ] diferencia x1 x0 dividida de orden 1 en x0 y x1 . En general, al imponer la condicin pn (xi ) = yi sale la condicin yi = A0 +A1 (xi x0 )+A2 (xi x0 ) (xi x1 )+ +Ai (xi x0 ) (xi x1 ) (xi xn1 ) de donde Ai = f [x0 , x1 , , xi ] diferencia dividida de orden i en x0 , x1 , , xi . La ley de recurrencia me permite ir formando las diferencias divididas de orden superior a partir de las de un orden menos. f [x1 , x2 , , xk ] f [x0 , x1 , , xk1 ] xk x0

Tabla de diferencias divididas

xi x0 x1 x2 x3 . . . xn

yi y0 y1 y2 yi . . . yn

f [xi ] f [x0 ] f [x1 ] f [x2 ] f [x2 ] . . .

f [xi , xi+1 ] f [x0 , x1 ] f [x1 , x2 ] f [x2 , x3 ] . . .

f [xi , xi+1 , xi+2 ]

f [x0 , x1 , , xn ]

f [x0 , x1 , x2 ] f [x1 , x2 , x3 ] . . . f [x0 , x1 , , xn ]

f [xn ] f [xn1 , xn ] f [xn2 , xn1 , xn ]

Los elementos de la diagonal son, precisamente, los coecientes Ai del polinomio pn escrito en la forma de Newton. Frmula de Newton para el polinomio de interpolacin pn (x) = f [x0 ] + f [x0 , x1 ] (x x0 ) + f [x0 , x1 , x2 ] (x x0 ) (x x1 ) + + f [x0 , x1 , , xn ] (x x0 ) (x x1 ) (x xn1 )

Ejemplo 1
Calcule el polinomio de interpolacin de Newton para los datos xi yi -2 -1 2 3 4 1 4 9

Solucin: El polinomio que se nos pide se puede escribir p3 (x) = A0 + A1 (x x0 ) + A2 (x x0 ) (x x1 ) + A3 (x x0 ) (x x1 ) (x x2 ) Formamos la tabla de diferencias divididas para obtener los coecientes xi yi 2 4 1 2 3 1 4 9 f [xi , xi+1 ] f [x0 , x1 ] = f [xi , xi+1 , xi+2 ] f [x0 , x1 , x2 , x3 ]

14 = 3 1 (2) 41 1 (3) f [x1 , x2 ] = f [x0 , x1 , x2 ] = =1 =1 2 (1) 2 (2) 94 51 11 =5 = 1 f [x0 , x1 , x2 , x3 ] = =0 f [x2 , x3 ] = f [x1 , x2 , x3 ] = 32 3 (1) 3 (1) Solucin p3 (x) = 4 3 (x + 2) + (x + 2) (x + 1) + 0 (x + 2) (x + 1) (x 2) = x2

Ejemplo 2
Calcule el polinomio de interpolacin para la funcin f (x) = |x| en los nodos xi = 4, 1, 2, 5 y 7 usando la tabla de diferencias divididas. Solucin: formamos la tabla de diferencias divididas xi yi = |xi | 4 4 1 2 5 7 1 2 5 7 f [xi , xi+1 ] f [x0 , x1 ] = f [xi , xi+1 , xi+2 ]

14 = 1 1 (4) 1 21 = f [x1 , x2 ] = f [x0 , x1 , x2 ] = 2 (1) 3 52 =1 f [x2 , x3 ] = f [x1 , x2 , x3 ] = 52 75 =1 f [x3 , x4 ] = f [x1 , x2 , x3 ] = 75

(1) = 2 (4) 1 1 3 = 5 (1) 11 =0 72

1 3

2 9 1 9

f [xi , xi+1 , xi+2 , xi+3 ]

f [x0 , x1 , x2 , x3 , x4 ]

Luego el polinomio de interpolacin en la forma de Newton que nos piden es 2 1 1 p (x) = 4(x + 4)+ (x + 4) (x + 1) (x + 4) (x + 1) (x 2) (x + 4) (x + 1) (x 2) (x 5) 9 81 1728

2 1 9 = 5 (4) 81 1 0 9 1 = f [x1 , x2 , x3 , x4 ] = 7 (1) 72 f [x0 , x1 , x2 , x4 ] =

1 9

1 1 72 + 81 1 = f [x0 , , x4 ] = 7+4 1728

Problemas de la relacin propuesta

1. (c) Calcule el polinomio de interpolacin para la tabla dada usando el sistema y la frmula de Newton xi 0 2 4 6 yi 0.25 0.6 0.9 1 Compruebe que ambas formas dan lugar al mismo polinomio. Solucin: i) Por medio del sistema p (x) = a0 + a1 x + a2 x2 + a3 x3 Imponemos las condiciones de interpolacin p (0) = p (2) = p (4) = p (6) = Matricialmente 1 1 1 1 0 2 4 6 0,25 a0 + 2a1 + 4a2 + 8a3 = 0,6 a0 + 4a1 + 16a2 + 64a3 = 0,9 a0 + 6a1 + 36a2 + 216a3 = 1

p (x) = 0,25 + 0,162 5x + 0,012 5x2 3. 125 103 x3 ii) Por el mtodo de Newton

0,25 a0 0 0 4 8 a1 = 0,6 0,9 16 64 a2 1 36 216 a3 a0 0,25 a1 0,162 5 = Solucin : a2 0,012 5 3 3. 125 10 a3

p (x) = A0 + A1 x + A2 x (x 2) + A3 x (x 2) (x 4) xi 0 2 4 6 yi 0,25 0,6 0,9 1 f [xi , xi+1 ] 0,6 0,25 = 0,175 20 0,9 0,6 = 0,15 42 1 0,9 = 0,05 64 f [xi , xi+1 , xi+2 ] f [x0 , x1 , x2 , x3 ]

0,15 0,175 = 0,006 25 40 0,05 0,15 = 0,025 62

0,025 (0,006 25) = 3. 125 103 60

Solucin: p (x) = 0,25 + 0,175 x 0,006 25x (x 2) 3. 125 103 x (x 2) (x 4) = 0,162 5x + 0,012 5x2 3. 125 103 x3 + 0,25