ACTIVIDAD EVALUATIVA EJE 3

ANALISIS NUMERICO

FACULTAD DE INGENIERIAS

INGENIERIA INDUSTRIAL

PRESENTADO POR:
CARLOS ANDRES ROLDAN
ISABEL ANDREA MORENO
WENDY YOHANA LOPEZ

FUNDACION UNIVERSITARIA DEL AREA ANDINA

SEPTIEMBRE 2023
COLOMBIA
 INTRODUCCION

En el vasto campo de la matemática y la ciencia de datos, la capacidad de analizar y modelar

datos experimentales desempeña un papel crucial en la toma de decisiones informadas y la
comprensión de fenómenos del mundo real. Para lograrlo, se han desarrollado herramientas
y técnicas poderosas, entre las cuales destacan la interpolación polinomial y el ajuste de
curvas. En este trabajo, exploraremos a fondo estos dos conceptos fundamentales, que
constituyen la columna vertebral de la aproximación y modelado de datos en diversos campos
de la ciencia y la ingeniería.

La interpolación polinomial nos permite conectar puntos de datos discretos con suaves curvas
polinómicas que pasan exactamente por cada punto, proporcionando una representación
continua de los datos en el espacio de las variables independientes. Por otro lado, el ajuste de
curvas nos permite encontrar modelos matemáticos que se adapten de manera óptima a un
conjunto de datos, incluso cuando la relación subyacente entre las variables no es puramente
polinómica.
PUNTO 1
a.
𝑥0 = 0 𝑓(𝑥0 ) = 1
𝑥1 = 2 𝑓(𝑥1 ) = 3
𝑥2 = 3 𝑓(𝑥2 ) = 0

Los polinomios de Lagrange están definidos por:

𝑛−1 𝑥−𝑥𝑗 𝑖 = 0,1, . . , 𝑛 − 1

𝑃𝑛 (𝑥) = ∑𝑖=0 𝑓(𝑥𝑖 )𝐿𝑖 (𝑥) Donde 𝐿𝑖 (𝑥) = ∏𝑛−1
𝑖=0
𝑗≠𝑖
𝑥𝑖 −𝑥𝑗 𝑗 = 0,1, … , 𝑛 − 1

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

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

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

𝑃(𝑥) = 𝑓(𝑥0 ) 𝐿0 (𝑥) + 𝑓(𝑥1 ) 𝐿1 (𝑥) + 𝑓(𝑥2 ) 𝐿2 (𝑥)

(𝑥 − 2)(𝑥 − 3) (𝑥 − 0)(𝑥 − 3) (𝑥 − 0)(𝑥 − 2)

𝑃(𝑥) = 1 ∗ +3∗ +0∗
(0 − 2)(0 − 3) (2 − 0)(2 − 3) (3 − 0)(3 − 2)
 (𝑥 − 2)(𝑥 − 3) (𝑥)(𝑥 − 3)
𝑃(𝑥) = +3∗
(−2)(−3) (2)(−1)
(𝑥 2 − 5𝑥 + 6) (𝑥 2 − 3𝑥) 1 2 3
= +3∗ = (𝑥 − 5𝑥 + 6) − (𝑥 2 − 3𝑥)
6 (−2) 6 2
1 5 3 9
𝑃(𝑥) = ( 𝑥 2 − 𝑥 + 1) + (− 𝑥 2 + 𝑥)
6 6 2 2

4 11
𝑃(𝑥) = − 𝑥 2 + 𝑥+1
3 3

b.
𝑥0 = −1 𝑓(𝑥0 ) = 0
𝑥1 = 2 𝑓(𝑥1 ) = 1
𝑥2 = 3 𝑓(𝑥2 ) = 1
𝑥3 = 5 𝑓(𝑥3 ) = 2
Se aplica el mismo procedimiento del ejercicio anterior.
(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )(𝑥 − 𝑥3 ) (𝑥 − 2)(𝑥 − 3)(𝑥 − 5)
𝐿0 (𝑥) = =
(𝑥0 − 𝑥1 )(𝑥0 − 𝑥2 )(𝑥0 − 𝑥3 ) (−1 − 2)(−1 − 3)(−1 − 5)

(𝑥 − 𝑥0 )(𝑥 − 𝑥2 )(𝑥 − 𝑥3 ) (𝑥 + 1)(𝑥 − 3)(𝑥 − 5)

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

(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥3 ) (𝑥 + 1)(𝑥 − 2)(𝑥 − 5)

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

(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 ) (𝑥 + 1)(𝑥 − 2)(𝑥 − 3)

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

𝑃(𝑥) = 𝑓(𝑥0 ) 𝐿0 (𝑥) + 𝑓(𝑥1 ) 𝐿1 (𝑥) + 𝑓(𝑥2 ) 𝐿2 (𝑥)

(𝑥 − 2)(𝑥 − 3)(𝑥 − 5) (𝑥 + 1)(𝑥 − 3)(𝑥 − 5)

𝑃(𝑥) = 0 ∗ +1∗ +1
(−1 − 2)(−1 − 3)(−1 − 5) (9)
(𝑥 + 1)(𝑥 − 2)(𝑥 − 5) (𝑥 + 1)(𝑥 − 2)(𝑥 − 3)
∗ +2∗
(−8) 36
 (𝑥 3 − 7𝑥 2 + 7𝑥 + 15) (𝑥 3 − 6𝑥 2 + 3𝑥 + 10) (𝑥 3 − 4𝑥 2 + 𝑥 + 6)
𝑃(𝑥) = + +2∗
(9) (−8) 36

1 7 7 5 1 3 3 5 1 2 1 1
𝑃(𝑥) = ( 𝑥 3 − 𝑥 2 + 𝑥 + ) + (− 𝑥 3 + 𝑥 2 − 𝑥 − ) + ( 𝑥 3 − 𝑥 2 + 𝑥+ )
9 9 9 3 8 4 8 4 18 9 18 3
1 3 1 2 11 3
𝑃(𝑥) = 𝑥 − 𝑥 + 𝑥+
24 4 24 4

c.
𝑥0 = 0 𝑓(𝑥0 ) = −2
𝑥1 = 2 𝑓(𝑥1 ) = 1

𝑥2 = 4 𝑓(𝑥2 ) = 4

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

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

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

𝑃(𝑥) = 𝑓(𝑥0 ) 𝐿0 (𝑥) + 𝑓(𝑥1 ) 𝐿1 (𝑥) + 𝑓(𝑥2 ) 𝐿2 (𝑥)

(𝑥 − 2)(𝑥 − 4) (𝑥 − 0)(𝑥 − 4) (𝑥 − 0)(𝑥 − 2)

𝑃(𝑥) = −2 ∗ +1∗ +4∗
(0 − 2)(0 − 4) (2 − 0)(2 − 4) (4 − 0)(4 − 2)
(𝑥 − 2)(𝑥 − 4) (𝑥)(𝑥 − 4) (𝑥)(𝑥 − 2)
𝑃(𝑥) = −2 ∗ + +4∗
8 (−4) 8
1 1 1
𝑃(𝑥) = − (𝑥 2 − 6𝑥 + 8) + (− 𝑥 2 + 𝑥) + ( 𝑥 2 − 𝑥)
4 4 2
1 3 1 1
𝑃(𝑥) = (− 𝑥 2 + 𝑥 − 2) + (− 𝑥 2 + 𝑥) + ( 𝑥 2 − 𝑥)
4 2 4 2
 3
𝑃(𝑥) = 𝑥 − 2
2
PUNTO 2.
a.
𝑥0 = 0 𝑓(𝑥0 ) = 1
𝑥1 = 2 𝑓(𝑥1 ) = 3
𝑥2 = 3 𝑓(𝑥2 ) = 0

El polinomio de interpolación de Newton por diferencias divididas está definido por:

𝑃𝑛 (𝑥) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) + ⋯ + 𝑎𝑛 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 ) … (𝑥 − 𝑥𝑛−2 )

𝑎0 , 𝑎1 , 𝑎2 , … . , 𝑎𝑛

𝑗 = 0,1,2, … , 𝑛 − 1

𝑗 𝑥𝑗 𝑓(𝑥𝑗 )
0 0 1
1 2 3 1
2 3 0 -3 _- 4/3

𝑓(𝑥1 ) − 𝑓(𝑥0 ) 3 − 1 2
𝑓(𝑥0 , 𝑥1 ) = = = =1
𝑥1 − 𝑥0 2−0 2
𝑓(𝑥2 ) − 𝑓(𝑥1 ) 0 − 3 −3
𝑓(𝑥1 , 𝑥2 ) = = = = −3
𝑥2 − 𝑥1 3−2 1
𝑓(𝑥1 , 𝑥2 ) − 𝑓(𝑥0 , 𝑥1 ) −3 − 1 4
𝑓(𝑥0 , 𝑥1 , 𝑥2 ) = = =−
𝑥2 − 𝑥0 3−0 3

𝑃𝑛 (𝑥) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
4
𝑃(𝑥) = 1 + 1(𝑥 − 0) + (− ) (𝑥 − 0)(𝑥 − 2)
3
4
𝑃(𝑥) = 1 + 1(𝑥) + (− ) (𝑥 2 − 2𝑥)
3
4
𝑃(𝑥) = 1 + 𝑥 + (− ) (𝑥 2 − 2𝑥)
3

4 11
𝑃(𝑥) = − 𝑥 2 + 𝑥+1
3 3
b.
𝑥0 = −1 𝑓(𝑥0 ) = 0
𝑥1 = 2 𝑓(𝑥1 ) = 1
𝑥2 = 3 𝑓(𝑥2 ) = 1
𝑥3 = 5 𝑓(𝑥3 ) = 2

𝑗 𝑥𝑗 𝑓(𝑥𝑗 )
0 -1 0
1 2 1 1/3
2 3 1 0 -1/12
3 5 2 1/2 1/6 1/24

𝑓(𝑥1 ) − 𝑓(𝑥0 ) 1 − 0 1
𝑓(𝑥0 , 𝑥1 ) = = =
𝑥1 − 𝑥0 2+1 3
𝑓(𝑥2 ) − 𝑓(𝑥1 ) 1 − 1 0
𝑓(𝑥1 , 𝑥2 ) = = = =0
𝑥2 − 𝑥1 3−2 1
𝑓(𝑥3 ) − 𝑓(𝑥2 ) 2 − 1 1
𝑓(𝑥2 , 𝑥3 ) = = =
𝑥3 − 𝑥2 5−3 2
1
𝑓(𝑥1 , 𝑥2 ) − 𝑓(𝑥0 , 𝑥1 ) 0 − 3 1
𝑓(𝑥0 , 𝑥1 , 𝑥2 ) = = =−
𝑥2 − 𝑥0 3+1 12
1
𝑓(𝑥2 , 𝑥3 ) − 𝑓(𝑥1 , 𝑥2 ) 2 − 0 1
𝑓(𝑥1 , 𝑥2 , 𝑥3 ) = = =
𝑥3 − 𝑥1 5−2 6
1 1 1
𝑓(𝑥1 , 𝑥2 , 𝑥3 ) − 𝑓(𝑥0 , 𝑥1 , 𝑥2 ) 6 + 12 4 1
𝑓(𝑥0 , 𝑥1 , 𝑥2 , 𝑥3 ) = = = =
𝑥3 − 𝑥0 5+1 6 24

𝑃𝑛 (𝑥) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) + 𝑎3 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )

1 1 1
𝑃(𝑥) = 0 + (𝑥 + 1) + (− ) (𝑥 + 1)(𝑥 − 2) + ( ) (𝑥 + 1)(𝑥 − 2)(𝑥 − 3)
3 12 24
1 1 1 1
𝑃(𝑥) = 𝑥 + + (− ) (𝑥 2 − 𝑥 − 2) + ( ) (𝑥 3 − 4𝑥 2 + 𝑥 + 6)
3 3 12 24
1 3 1 2 11 3
𝑃(𝑥) = 𝑥 − 𝑥 + 𝑥+
24 4 24 4
c.

𝑥0 = 0 𝑓(𝑥0 ) = −2
𝑥1 = 2 𝑓(𝑥1 ) = 1

𝑥2 = 4 𝑓(𝑥2 ) = 4

𝑗 𝑥𝑗 𝑓(𝑥𝑗 )
0 0 -2
1 2 1 3/2
2 4 4 3/2 _0

𝑓(𝑥1 ) − 𝑓(𝑥0 ) 1 + 2 3
𝑓(𝑥0 , 𝑥1 ) = = =
𝑥1 − 𝑥0 2−0 2
𝑓(𝑥2 ) − 𝑓(𝑥1 ) 4 − 1 3
𝑓(𝑥1 , 𝑥2 ) = = =
𝑥2 − 𝑥1 4−2 2
3 3
𝑓(𝑥1 , 𝑥2 ) − 𝑓(𝑥0 , 𝑥1 ) 2 − 2
𝑓(𝑥0 , 𝑥1 , 𝑥2 ) = = =0
𝑥2 − 𝑥0 4−0

𝑃𝑛 (𝑥) = 𝑎0 + 𝑎1 (𝑥 − 𝑥0 ) + 𝑎2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
3
𝑃(𝑥) = −2 + (𝑥) + 0(𝑥)(𝑥 − 2)
2
3
𝑃(𝑥) = 𝑥 − 2
2
PUNTO 3
Interpolación lineal:
𝑦2 − 𝑦1
𝑦= ∗ (𝑥 − 𝑥1 ) + 𝑦1
𝑥2 − 𝑥1
𝑦2 − 𝑦1
𝑦= ∗ (𝑥 − 𝑥1 ) + 𝑦1
𝑥2 − 𝑥1
Para x= 0.2:
−3 − 1
𝑦= ∗ (0.2 − 0) + 1
1−0
𝑦 = 0.2
Para x= 3:
5+3
𝑦= ∗ (3 − 1) − 3
4−1
𝑦 = 2.33

Interpolación cuadrática:

𝑦 = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐
Se calculan los coeficientes de la ecuación reemplazando los puntos dados

1 = 𝑎 ∗ 02 + 𝑏 ∗ 0 + 𝑐 1=𝑐
−3 = 𝑎 ∗ (1)2 + 𝑏(1) + 𝑐 } → −3 = 𝑎 + 𝑏 + 𝑐 }
5 = 𝑎 ∗ (4)2 + 𝑏 ∗ 4 + 𝑐 5 = 16𝑎 + 4𝑏 + 𝑐

Resolviendo el sistema de ecuaciones:

𝑎 + 𝑏 = −4
16𝑎 + 4𝑏 = 4
De la ecuación 1:

𝑏 = −4 − 𝑎
Reemplazando en la ecuación 2:

16𝑎 + 4(−4 − 𝑎) = 4

16𝑎 − 16 − 4𝑎 = 4
12𝑎 = 20
5
𝑎=
3
 5 17
𝑏 = −4 − 𝑎 = −4 − =−
3 3
La ecuación cuadrática sería entonces:
5 2 17
𝑦= 𝑥 − 𝑥+1
3 3

Para x= 0.2:
5 17
𝑦 = (0.2)2 − (0.2) + 1
3 3
1
𝑦=−
15
Para x=3:
5 17
𝑦= (3)2 − (3) + 1
3 3
𝑦 = −1
PUNTO 4.

𝑦 = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐
Se calculan los coeficientes de la ecuación reemplazando los puntos dados, para simplificar los datos
se hace un reajuste así:

Años 0 5 10 15
Empleados 24 30 38 43

Se toman solo 3 puntos:

24 = 𝑎 ∗ 02 + 𝑏 ∗ 0 + 𝑐 24 = 𝑐
30 = 𝑎 ∗ (5)2 + 𝑏(5) + 𝑐 } → 30 = 25𝑎 + 5𝑏 + 𝑐 }
38 = 𝑎 ∗ (10)2 + 𝑏 ∗ 10 + 𝑐 38 = 100𝑎 + 10𝑏 + 𝑐

Resolviendo el sistema de ecuaciones:

25𝑎 + 5𝑏 = 6
100𝑎 + 10𝑏 = 14
De la ecuación 1:
6 − 5𝑏
𝑎=
25
Reemplazando en la ecuación 2:
6 − 5𝑏
100 ( ) + 10𝑏 = 14
25
24 − 20𝑏 + 10𝑏 = 14

−10𝑏 = −10
𝑏=1
6−5 1
𝑎= =
25 25
La ecuación cuadrática sería entonces:
1 2
𝑦= 𝑥 + 𝑥 + 24
24

Para encontrar el numero de empleados en el año 2000 se reemplaza x=20:

1
𝑦= (20)2 + 20 + 24 = 60
24

Para el año 2000 había aproximadamente 60 empleados.

 CONCLUSION

En este trabajo, hemos explorado en profundidad dos herramientas esenciales en el análisis

de datos y el modelado matemático: la interpolación polinomial y el ajuste de curvas. Estas
técnicas, aunque distintas en enfoque y aplicabilidad, comparten un objetivo común:
transformar datos dispersos en modelos comprensibles y predictivos.

A lo largo del trabajo, hemos aprendido cómo la interpolación polinomial nos permite
conectar puntos de datos de manera precisa mediante polinomios, brindando una
representación continua de información discreta. También hemos descubierto cómo el ajuste
de curvas nos capacita para encontrar modelos matemáticos que se ajusten óptimamente a
datos reales, incluso cuando las relaciones subyacentes son complejas y no se limitan a
polinomios.