Interpolación polinomial y ajustes de curvas

Actividad evaluativa Eje 3

Análisis numérico

Presentado Por:

Luisa Fernanda Jaimes Hernández

Tutor:

Richard. Fabian reyes ramos

Fundación Universitaria Del Área Andina.

Ingeniería industrial virtual

2024
Objetivo de aprendizaje

Aplicar los diferentes métodos de interpolación y ajuste de curvas desarrollados en

el eje 3 para resolver ejercicios. prácticos.

Instrucciones:

1. Realice el estudio del referente de pensamiento del eje 3 profundizando

sobre cada una de las lecturas complementarias, video cápsulas y el video

resumen.

2. Si lo requiere, revise otros documentos que crea le pueden ayudar a tener

una visión diferente del tema tratado, puede incluir videos.

3. Tenga a la mano las diferentes expresiones de la interpolación lineal,

cuadráticas, de Lagrange, de sumas divididas de Newton y las de mínimos

cuadrados a manera de una ficha nemotécnica.

4. En el desarrollo de cada ejercicio, trate de ser lo más ordenado posible en la

argumentación de cada paso y en la aplicación de cada expresión de los métodos

de interpolación.
 1. Determinar la función lineal de interpolación que pasa por los puntos

(-3, 0), (3,1). Luego interpole el valor a=1 y extrapole el valor de b=4

La interpolación lineal nos permite determinar el conjunto de valores que

toma una función en un determinado intervalo dados dos puntos o valores iniciales.

Para hallar dicha función, basta con aplicar de forma correcta la expresión:

𝑦1− 𝑦0
𝑓(𝑥 ) = 𝑦0+ (𝑥 − 𝑥 0 )
𝑥1 − 𝑥0

Definimos los puntos P (x 0, y 0) = (−3,0) y Q (x1, y 1) = (3,1). Reemplazando

estos puntos en la expresión de arriba se obtiene:

1−0
f(x) = 0 + (x − (−3)
3 − (−3)

1
f(x) = (x + 3)
6

1 1
f(x) = x
6 2
En la siguiente figura se puede apreciar la gráfica de la función encontrada:

Una vez que hallemos la función de interpolación procedemos a

encontrar los valores de interpolación y extrapolación.

1 1 4 2
Interpolando a = 1 obtenemos f (1) =6 (1+3) ⤍ 6 (4) ⤍6 =3⤍ 0.6667

1 1 7
Extrapolando b = 4 obtenemos f (4) = 6 (4+3) ⤍ 6 (7) ⤍ 6⤍ 1.16667

Interpolación cuadrática

2. Mediante la interpolación cuadrática, determine la función cuadrática

que contiene las coordenadas (0,-3), (1,0), (3, 0). Luego, interpole el
valor a= 2 y extrapole el valor b= -1

A continuación, tenemos los puntos:

(x0−y0) = (0,−3)

(x1−y1) = (1,−0)

(x2−y2) = (3,0)
 Para el caso de la interpolación cuadrática, se desea encontrar un polinomio

de la forma:

f (x)=a+b(x−x1) +c(x−x1)(x−x2)

Reemplazamos el valor de x1, x2 para armar un sistema de ecuaciones y

hallar sus tres incógnitas

f(x) = a + b(x − 0) + c(x − 0)(x − 1)

f(x) = a + b(0 − 0) + c(0 − 0)(0 − 1)

Ahora reemplazamos la imagen de x1 y nos queda:

−3 = 𝑎 + 𝑏(0 − 0) + 𝑐 (0 − 0)(0 − 1)

−3 = 𝑎

Seguimos el mismo procedimiento anterior para las siguientes incógnitas:

0 = −3+b (1−0) +c (1−0(1−1)

0 = − 3+b (1) +c (1)

0 = − 3+b b =3
Hallamos c teniendo

0=−3+3(3−0) +c (3−0) (3−1)

0=−3+9+6c

6
0= 6+ 6 c ⤍ −6 = ¿6 c ⤍ c = = −1 ¿
6

Teniendo ya los valores para completar el polinomio procedemos a

reemplazar

f (x) = − 3 +3 (x−0) =1 (x−0) (x−1)

f (x) = −3+3 x − x2 + x

Finalmente tenemos f (x)=−x2+4 x−3⤍ función de interpolación

Interpolando a = 2

f (x)=− ¿¿

¿−4+8−3

f (x)=−7+8=1

Extrapolando b = -1f (x) = −x2+4 x−3

f (x) = − ¿¿

¿−1−4−3

f (x) =−8
En la siguiente figura se puede apreciar la gráfica de la función encontrada:

3.Interpolación de Lagrange
 Dados los puntos:

(x0 − y0 ) = (0,0)

(x1 − y1 ) = (−3, −84)

(𝑥2 − 𝑦2 ) = (2,26)

(𝑥3 − 𝑦3 ) = (−4, −196)

 Para hallar los polinomios de LaGrange tenemos

(𝑥 − 𝑥1)(𝑥 − 𝑥2)(𝑥 − 𝑥3)

𝐿0(𝑥) =
(𝑥0 − 𝑥1)(𝑥0 − 𝑥2)(𝑥0 − 𝑥3)

(𝑥 + 3)(𝑥 − 2)(𝑥 + 4)
𝐿0(𝑥) =
(0 + 3)(0 − 2)(0 + 4)

(𝑥 + 3)(𝑥 − 2)(𝑥 + 4)
𝐿0(𝑥) =
−24
 = −1/24𝑥3 − 5/242 + 1/12𝑥 + 1

 Hallamos L1(x)

(𝑥 − 𝑥0)(𝑥 − 𝑥2)(𝑥 − 𝑥3
𝐿1(𝑥) =
(𝑥1 − 𝑥0)(𝑥1 − 𝑥2)(𝑥1 − 𝑥3)
(𝑥 − 0)(𝑥 − 2)(𝑥 + 4)
𝐿1(𝑥) =
(−3 − 0)(−3 − 2)(−3 + 4)
(𝑥 − 0)(𝑥 − 2)(𝑥 + 4)
𝐿1(𝑥) = = 1/15𝑥3 + 2/15𝑥2 − 8/15
15

 Hallamos L2(x)
(𝑥 − 𝑥0)(𝑥 − 𝑥1)(𝑥 − 𝑥3)
𝐿2(𝑥) =
(𝑥2 − 𝑥0)(𝑥2 − 𝑥1)(𝑥2 − 𝑥3)

(𝑥 − 0)(𝑥 + 3)(𝑥 + 4)
𝐿2(𝑥) =
(2 + 0)(2 + 3)(2 + 4)

(𝑥 − 0)(𝑥 + 3)(𝑥 + 4)
𝐿2(𝑥) = = 1/60𝑥3 + 7/60𝑥2/+1/5𝑥
60

 Hallamos L3(x)
 (𝑥 − 𝑥0)(𝑥 − 𝑥1)(𝑥 − 𝑥2)
𝐿3(𝑥) =
(𝑥3 − 𝑥0)(𝑥3 − 𝑥1)(𝑥3 − 𝑥2)

(𝑥 − 0)(𝑥 + 3)(𝑥 − 2)
𝐿3(𝑥) =
(−4 − 0)(−4 + 3)(−4 − 2)

(𝑥 − 0)(𝑥 + 3)(𝑥 − 2)
𝐿3(𝑥) = = −1/24𝑥3 − 1/24𝑥2 + 1/4
−24

 Tomando las expresiones y reemplazándolas en el

polinomio de LaGrange, tenemos:

𝑃(X) = ∑3k=0 yklk (x) = 0L(x) − 84L1(X) + 26L2(X) − 196L3(X)

P(X) = X + 3X 3

 Cuando se destruyen paréntesis y se realizan las operaciones

indicadas, utilizando herramientas computacionales obteniendo como

resultado.

2 3 2 8 1 7 1 1
− 84 ( 𝑋 + 𝑋 2 − 𝑋) + 26 ( 𝑋 3 + 𝑋 2 + + 𝑋)
15 15 15 60 60 5 5

1 3 1 2 1
− 196 (− 𝑋 − 𝑋 + 𝑋)
24 24 4

 Desarrollar

1 3 2 2 8 28𝑋 3 56𝑋 2 224𝑋

84 ( 𝑋 + 𝑋 − 𝑋) : − − +
15 15 15 5 5 5
 28𝑥 3 56𝑥 2 224𝑥 1 7 1
− − + + 26 ( 𝑋 3 + 𝑋 2 − 𝑋)
5 5 5 60 60 5
1 3 1 2 1
− 196 (− 𝑋 + 𝑋 + )
24 24 4

1 3 7 2 1 13𝑋 3 91𝑋 2 26𝑋

26 ( 𝑋 + 𝑋 − 𝑋) : − − +
60 60 5 5 5 5

28𝑥 3 56𝑥 2 224𝑥 13𝑥 3 91𝑥 2 26𝑥

=− − + + + +
5 5 5 30 30 5
1 3 1 2 1
− 196 (− 𝑋 − 𝑋 + 𝑥)
24 24 4

1 3 1 1 49𝑋 3 49𝑋 2 49𝑋

−196 (− 𝑋 − 𝑋 2 − 𝑥) : − − + − 49𝑥
24 24 24 6 6 6

28𝑥 3 56𝑥 2 224𝑥 13𝑥 3 91𝑥 2 26𝑥 49𝑥 3 49𝑥 2

=− − + + + + + + + − 49𝑥
5 5 5 30 30 5 6 6

 Factorizacion

28𝑥 3 56𝑥 2 224𝑥 13𝑥 3 91𝑥 2 26𝑥 49𝑥 3 49𝑥 2

=− − + + + + + + + − 49𝑥
5 5 5 30 30 5 6 6
∶ 3𝑥 3 + 0 + 𝑥

= 3𝑥 3 + 0 + 𝑥

 Aplicar la propiedad: a + 0 = a

= 3𝑥 3 + 𝑥

f (x)=x+3 x3
 Para obtener la interpolación en el punto x = -2 debemos reemplazar

en todas las ecuaciones de L0 (x), L1 (x), L2(x), L3 (x)para obtener el

polinomio resultante.

𝐿0(𝑥)
(−2 + 3)(−2 − 2)(−2 + 4) 1
= =
−24 3

𝐿1(𝑥)
(−2 − 0)(−2 − 2)(−2 + 4)
=
15
16
=
15

𝐿2(𝑥)
(−2 − 0)(−2 + 3)(−2 + 4)
=
60
−1
=
15

𝐿3(𝑥)
(−2 − 0)(−2 + 3)(−2 − 2) 1
= =
−24 3

 Por tanto:
 ∞

𝑝(−2) = ∑ 𝑦𝑘𝑙𝑘(−2)
𝑛=1
= 0L0(−2) − 84L1(−2) + 26L2(−2) − 196L3(−2)
= −26

 Solución computacional:

16 1 1
−84 ( ) + 26 (− ) − 196 (− )
15 15 3

16 448
−84 ( )=−
15 5

1 26
26 (− )=−
15 15

1 196
196 (− ) = −
3 3

448 26 196
=− − − (− )
5 15 3

448 26 196
=− − − (− ) = 26
5 15 3
= − 26
  En la siguiente figura se puede apreciar la gráfica de la función

encontrada:

4.

pn(x)=f [x 0]+f [x0, x1](x−x0) +f [x 0, x 1, x 2] (x−x 0) (x−x1) +…+f [x 0, x1, …, xn]

(x−x 0)(x−x1)
 f[xk ] = f(xk )

f[xk+1 ] − f[xk ]
f[xk xk+1 ] =
xk+1 − xk

f[xk+1 xk+2 ] − f[xk xk+1 ]

f[xk xk+1 xk+2 ] =
xk+2 − xk

f[xk+1 xk+2 xk+3 ] − f[xk xk+1 xk+2 ]

f[xk xk+1 xk+2 xk+3 ] =
xk+3 − xk

En las que aparecen las diferencias divididas f [x0,..., xi], obtenidas a partir

de los valores proporcionados por la tabla inicial. Calculamos entonces la tabla de

diferencias divididas:

xk Yk F [xk || F [xk || F [xk || F [xk || F [xk ||

xk+1] xk+2] xk+3] xk+4] xk+5]
3 -37

1 -7 -15

-7 -127 15 -3

-1 -1 21 -3 0

6 -127 -18 -3 0 0

5 -91 -36 -3 0 0 0

En donde se ha expresado la diferencia dividida f [xk, xk+1, …, xk+p] como f

[xk || xk + p]. La diagonal de la tabla de diferencias divididas, es entonces:

[−37,−15,−3,0,0,0], que se corresponde exactamente con el conjunto de valores que

aparece en la fórmula y, por tanto, los polinomios de Newton son los siguientes:
 P 0 (x) = −37(interpola en el primer punto)

P 1(x) = −15 (x−3) +p 0(x) =−15 x+8(interpola en los2 primeros puntos)

P 2(x)= −3(x−3) (x−1) +p1(x) = −1−3 x2−3 x (interpola en

todos los puntos)

O también:

P (x)=−37−15(x−3) −3(x−3) (x−1) =−1−3 x2−3 x

La gráfica del polinomio de interpolación:

p(x)=−1−3 x2−3 x
Luego para interpolar reemplazamos x = - 3 en el polinomio
encontrado.

p (−3) =−1−3 (−3)2−3(−3) = − 19

En la siguiente figura se puede apreciar la gráfica de la función encontrada:

 Primeramente, calculamos las medias de los valores de x y los valores de y,

la suma de los cuadrados de los valores de x, y la suma de cada valor de x

multiplicado por su valor correspondiente y.

x y Xy x2
-3 14 -42 9
-1 4 -4 1
1 2 2 1
3 8 24 9
5 22 110 25
 7 44 308 49
∑x=12 y=94 ∑x y=398 ∑ x2=94
Calculamos la pendiente aplicando la formula siguiente:

(∑ 𝑥 )(∑ 𝑦)
𝑀 = ∑ 𝑥𝑦 − 𝑛
( ∑ 𝑥)2
∑ 𝑧2
𝑛
(12) (94)
398 −
6
m = 94− ¿

Procedemos a calcular la intercepción en y

Primero, se calcula la media de los valores de x y la media de los valores de

∑ 𝑥 12
𝑥 𝑚𝑒𝑑𝑖𝑎 = = =2
𝑛 6

∑ 𝑥 94 47
𝑦 𝑚𝑒𝑑𝑖𝑎 = = =
𝑛 6 3

Usando la fórmula para calcular la intercepción en y.

 47
𝑏= − (3 ∗ 2)
3
47
𝑏= − (6)
3

𝑏 = 9.666667

Usamos la pendiente y la intercepción en y para formar la ecuación de la recta que

mejor se ajusta.

La pendiente de la recta es 3 y la intercepción en y es 9.666667

Luego la función lineal ajustada y=3 x+9.666667

En la siguiente figura se puede apreciar la gráfica de la función encontrada línea

ajustada
Conclusiones
 • El polinomio de interpolación suele usarse para estimar valores de una

función. tabulada, en las abscisas que no aparecen en la tabla.

• El aumento de grado no siempre mejora la aproximación.

• El polinomio es muy sensible a los errores de los datos.

Referencias bibliográficas

• https://www.academia.edu/36380579/

Interpolacion_y_ajuste_de_curvas

• Burden, R., Faires, J. y Solorio, P. (2011).

Análisis numérico. México, México: Cengage Learning.

• Chapra, S. y Canale, R. (2006). Métodos

numéricos para ingenieros. México, México:

McGraw-Hil

 https://www.geogebra.org/?lang=es