mínimos cuadrados orden

superior

ANDRES FELIPE DIAZ SILVA

UNIVERSIDAD SURCOLOMBIANA - INGENIERIA

 Asignatura:
Métodos Numéricos

Trabajo rubrica sobre mínimos cuadrados de orden superior

Presenta:
Andres Felipe diaz silva
Código: 20201186235

Docente
Ing. YAMIL ARMANDO CERQUERA ROJAS MSc.

Neiva, marzo 01 de 2022

 Tabla de contenido

1. Planteamiento del Problema......................................................................................................4

2. Análisis de la situación planteada...............................................................................................5
3. Revisión bibliográfica..................................................................................................................6
4. Planteamiento de la solución......................................................................................................7
5. Solución del problema................................................................................................................8
6. Análisis de Resultados...............................................................................................................14
7. Conclusiones.............................................................................................................................15
8. Bibliografía................................................................................................................................16
9. Anexos: Códigos........................................................................................................................17
1. Planteamiento del Problema

Se pone a cargar un teléfono móvil a partir desde que se descarga completamente y se

toman los valores de la batería cada 10 min para así efectuar el ajuste de curvas o de
función y saber cual es la que mejor se ajusta
2. Análisis de la situación planteada
Al obtener los datos de carga según el tiempo se proponen hallar ciertos valores para
hallar la funcion que mejor se ajuste a eso valores según que ecuación se dese utilizar
3. Revisión bibliográfica
No se tomó ninguna.
4. Planteamiento de la solución
Se utiliza la ecuación de mínimos cuadrados con ajuste lineal, parabólico, cubico y
de mayor orden para obtener el mínimo error y saber cuál es la función que mejor
se ajusta a los datos; todos estos pasos tienen una solución fácil utilizando scilab.
5. Solución del problema

Batería% 1.6 4. 7.7 13 1 24 2 36 4 52

7 9 8 4
tiempo 1 2 3 4 5 6 7 8 9 10

Grafica 1. Porcentaje de batería vs Tiempo (min)

Mínimos cuadrados de forma manual.

[línea recta]

Ecuaciones para hallar a y b.

● 𝑎Σ𝑥2+𝑏Σ𝑥=Σ𝑌𝑥
● 𝑎Σ𝑥+𝑛𝑏=Σ𝑌

Procedimiento

A continuación, se presenta la información necesaria para formar el sistema de ecuaciones.

Tabla 1. Resumen de los datos.

n x y0 X2 y0x
1 1 1.6 1 1.6
2 2 4.7 4 9.4
3 3 7.7 9 23.1
 4 4 13 16 52
5 5 19 25 95
6 6 24 36 144
7 7 28 47 196
8 8 36 64 288
9 9 44 81 396
10 10 52 100 520
∑ 55 230 385 1725.1

Formando así las siguientes ecuaciones:

385a + 55b =1725.1
55a + 10b =230
Se realiza el sistema 2x2, mediante matrices.
𝑚𝑝= [
385 55
55 10 ]
a
𝑣𝑝=[ ]
b

𝑡𝑖𝑝= [ 1725.1
230 ]

Entonces 𝑚𝑝∗𝑣𝑝=𝑡𝑖𝑝, obteniendo mediante algebra, que:

𝑣𝑝=𝑖𝑛𝑣(𝑚𝑝)∗𝑡𝑖𝑝
Es decir
𝑣𝑝=𝑖𝑛𝑣
([
385 55
55 10 ])[ 1725.1
230 ]

Al resolver el sistema a través de Scilab se obtiene que:

𝑎= 5.5769697
𝑏= -7.6733333
Por lo tanto, se obtiene que la función calculada es:

y=5.5769697 x−7.6733333
Grafica 2. Porcentaje de batería vs Tiempo (min). Datos observados, datos calculados y
función.

Finalmente se procede a realizar el cálculo de errores para los valores calculados respecto de los
valores observados.

Tabla 2. Valores Observados vs Valores Calculados + Valor de A y B.

n y0 ye
A B
1 1,6 -2,0963636 629,8274659 13,66310386
2 4,7 3,4806061 381,0067382 1,486921483
3 7,7 9,0575758 194,3911926 1,843012053
4 13 14,634545 69,98083736 2,671737357
5 19 20,211515 7,775648595 1,467768595
6 24 25,788485 7,775648595 3,198678595
7 28 31,365455 69,98083736 11,32628736
8 36 36,942424 194,391187 0,888162996
9 44 42,519394 381,0067421 2,192194127
10 52 48,096364 629,827486 15,23837402
∑ 230 230,000 2566 53,97624045
ymed=
∑ y 0 = 23 =¿
n 10

A=∑ ( y e − ymed ) =2566

2

B=∑ ( y − y ) =53,97624045
e 0 2

2 A 2566
r= = =0.9793979
A+ B 2566+53,97624045

𝑟2= 0.9793979 [Ajuste lineal]

[Parábola]

Ecuaciones para hallar a y b y c.

● a ∑ x 4 + b ∑ x3 + c ∑ x 2=∑ y 0 x 2
● a ∑ x 3 +b ∑ x 2 +c ∑ x =∑ y 0 x
● a ∑ x 2 +b ∑ x +c∗n=∑ y 0

Procedimiento

A continuación, se presenta la información necesaria para formar el sistema de ecuaciones.

Tabla 3. Resumen de los datos.

n x y0 X2 X3 X4 y0x y0x2
1 1 1.6 1 1 1 1.6 1.6
2 2 4.7 4 8 16 9.4 18.8
3 3 7.7 9 27 81 23.1 69.3
4 4 13 16 64 256 52 208
5 5 19 25 125 625 95 475
6 6 24 36 216 1296 144 864
7 7 28 47 343 2401 196 1372
8 8 36 64 512 4096 288 2304
9 9 44 81 729 6561 396 3564
 10 10 52 100 1000 10000 520 5200
∑ 55 230 385 3025 25333 1725.1 14076.7

Formando así las siguientes ecuaciones:

25333 a+3025 b+ 385 c=14076.7

3025 a+385 b+55 c=1725.1
385 a+55 b+ 15 c=230

Se realiza el sistema 2x2, mediante matrices.

[ ]
25333 3025 385
𝑚𝑝= 3025 385 55
385 55 15

[]
a
𝑣𝑝= b
c

[ ]
14076.7
𝑡𝑖𝑝= 1725.1
230

Entonces 𝑚𝑝∗𝑣𝑝=𝑡𝑖𝑝, obteniendo mediante algebra, que:

𝑣𝑝=𝑖𝑛𝑣(𝑚𝑝)∗𝑡𝑖𝑝
Es decir

([ ] [ ])
25333 3025 385 14076.7
𝑣𝑝=𝑖𝑛𝑣 3025 385 55 1725.1
385 55 15 230
Al resolver el sistema a través de Scilab se obtiene que:
𝑎= 0.3041667
b=2.2311364
c= -0.9816667
Por lo tanto, se obtiene que la función calculada es:
2
y=0.3041667 x +2.2311364 x−−0.9816667

Grafica 3. Porcentaje de batería vs Tiempo (min). Datos observados, datos calculados y función.
Finalmente se procede a realizar el cálculo de errores para los valores calculados respecto
de los valores observados.

Valores Observados vs Valores Calculados

n y0 yc
A B
1 1,6 2 459,9465117 0,002149583
2 4,7 5 334,9898266 7,43817E-06
3 7,7 8 211,7245467 0,561364174
4 13 13 103,8453731 0,036273107
5 19 18 27,26738323 1,492839225
6 24 23 0,126132943 0,415828943
7 28 30 42,77755161 2,373001607
8 36 36 177,7979828 0,111616796
9 44 44 429,9842258 0,069663796
10 52 52 826,353443 0,06433122
∑ 230 230 2614,812978 5,127075890

ymed=
∑ y 0 = 230 =23
n 10

A=∑ ( y c − ymed ) =2614,812978

2
 B=∑ ( y − y ) =¿ 5,127075890 ¿
c 0 2

2 A 2614,812978
r= = =¿0.9980431
A+ B 2614,812978+5,127075890

𝑟2=0.9980431 [Ajuste parabólico]

[Cúbico]

Ecuaciones para hallar a, b, c y d.

● a ∑ x 6 +b ∑ x 5 +c ∑ x 4 +d ∑ x3 =∑ y 0 x 3
● a ∑ x 5 +b ∑ x 4 + c ∑ x 3+ ¿ d ∑ x 2=∑ y 0 x 2 ¿
● a ∑ x 4 + b ∑ x3 + c ∑ x 2+ d ∑ x=∑ y 0 x
● a ∑ x 3 +b ∑ x 2 +c ∑ x +d∗n=∑ y 0

Procedimiento

A continuación, se presenta la información necesaria para formar el sistema de

ecuaciones.

Tabla 5. Resumen de los datos

n x y0 X2 X3 X4 X5 X6 y0x y0x2 y0x3
1 1 1.6 1 1 1 1 1 1.6 1.6 1.6
2 2 4.7 4 8 16 32 64 9.4 18.8 37.6
3 3 7.7 9 27 81 243 729 23.1 69.3 207.9
4 4 13 16 64 256 1024 4096 52 208 832
5 5 19 25 125 625 3125 15625 95 475 2375
6 6 24 36 216 1296 7776 46656 144 864 5184
7 7 28 47 343 2401 16807 117649 196 1372 9604
8 8 36 64 512 4096 32768 262144 288 2304 18432
9 9 44 81 729 6561 59049 531441 396 3564 32076
10 10 52 100 1000 1000 10000 100000 520 5200 52000
0 0 0
∑ 55 230 385 3025 2533 22082 197840 1725. 14076.7 12075
3 5 5 1 0.1

Formando así las siguientes ecuaciones:

1978405 a+220825 b+25333 c +3025 d=120750.1

220825 a+25333 b+ 3025 c+385 d=14076.7
25333 a+3025 b+ 385 c+55 d =1725.1
3025 a+385 b+55 c+10 d =230

Se realiza el sistema 4x4, mediante el mediante matrices.

[ ]
1978405 220825 25333 3025
mp= 220825 25333 3025 385
25333 3025 385 55
3025 385 55 10

[]
a
b
vp=
c
d

[ ]
120750.1
14076.7
tip=
1725.1
230

Entonces mp∗vp=tip , obteniendo mediante algebra, que:

vp=inv ( mp )∗tip

Es decir

([ ]) [ ]
1978405 220825 25333 3025 120750.1
vp=inv 220825 25333 3025 385 ∗ 14076.7
25333 3025 385 55 1725.1
3025 385 55 10 230

Al resolver el sistema a través de Scilab se obtiene que:

a=0.0099262

b=0.1403846
c=2.9865190

d=−1.8333333

Por lo tanto, se obtiene que la función calculada es:

3 2
y=0.0099262 x + 0.1403846 x +2.9865190 x−1.8333333
Grafica 4. Porcentaje de batería vs Tiempo (min). Datos observados, datos calculados y
función.

Finalmente se procede a realizar el cálculo de errores para los valores calculados respecto
de los valores observados.
Tabla 6. Valores Observados vs Valores Calculados
n y0 yc
A B
1 1,6 1,30350 470,7382641 0,087914326
2 4,7 4,78065 331,944616 0,006504858
3 7,7 8,65769 205,7017902 0,917174541
4 13 12,99417 100,116594 3,39656E-05
5 19 17,84965 26,52610512 1,323305123
6 24 23,28368 0,080476044 0,513110044
7 28 29,35583 40,39654957 1,838269566
8 36 36,12564 172,2824517 0,015785661
9 44 43,65268 426,5332325 0,120630488
10 52 51,99650 840,797186 1,2229E-05
∑ 230 230 2615,117265 4,822740801
 ymed=
∑ y 0 = 230 =23
n 10

A=∑ ( y − ymed) =2615,117265

c 2

B=∑ ( y − y ) =4,822740801
c 0 2

A 2615,117265
r 2= = =0.9981592
A+ B 2615,117265+4,822740801

2
r =0.9981592[ Ajuste Cúbico]

[Polinomio de orden 4]

Ecuaciones para hallar a, b, c,d y e.

● a ∑ x 8 +b ∑ x 7 +c ∑ x 6+ d ∑ x 5 +e ∑ x 4 =∑ y 0 x 4
● a ∑ x 7 +b ∑ x 6 +c ∑ x 5+ d ∑ x 4 + e ∑ x 3=∑ y 0 x3
● a ∑ x 6 +b ∑ x 5 +c ∑ x 4 +d ∑ x3 +e ∑ x2 =∑ y 0 x 2
● a ∑ x 5 +b ∑ x 4 + c ∑ x 3+ d ∑ x 2 +e ∑ x =∑ y 0 x
● a ∑ x 4 + b ∑ x3 + c ∑ x 2+ d ∑ x +e∗n=∑ y 0

Procedimiento

A continuación, se presenta la información necesaria para formar el sistema de

ecuaciones.

Tabla 7. Resumen de los datos

n x y0 X2 X3 X4 X5 X6 X7 X8 y0x y0x2 y0x3
1 1 1.6 1 1 1 1 1 1 1 1.6 1.6 1.6
2 2 4.7 4 8 16 32 64 128 256 9.4 18.8 37.6
3 3 7.7 9 27 81 243 729 2187 6561 23.1 69.3 207.9
4 4 13 16 64 256 1024 4096 1638 6553 52 208 832
4 6
5 5 19 25 125 625 3125 15625 7812 3906 95 475 2375
5 25
6 6 24 36 216 1296 7776 46656 2799 1679 144 864 5184
 36 616
7 7 28 47 343 2401 1680 11764 8235 5764 196 1372 9604
7 9 43 801
8 8 36 64 512 4096 3276 26214 2097 1677 288 2304 1843
8 4 152 7216 2
9 9 44 81 729 6561 5904 53144 4782 4304 396 3564 3207
9 1 969 6721 6
1 1 52 10 1000 1000 1000 10000 1000 1000 520 5200 5200
0 0 0 0 00 00 0000 0000 0
0
∑ 5 230 38 3025 2533 2208 19784 1808 1.67 1725 14076. 1207
5 5 3 25 05 0425 7D+ .1 7 50.1
08

Formando así las siguientes ecuaciones:

1.677 D+08 a+18080425 b+1978405 c +220825 d +25333=1070375.5

18080425 a+1978405 b+220825 c +25333 d +3025 e=120750.1
1978405 a+220825 b+25333 c +3025 d +385 e=14076.7
220825 a+25333 b+ 3025 c+385 d +55 e=¿ 1725.1
25333 a+3025 b+ 385 c+55 d +10 e=230

Se realiza el sistema 4x4, mediante el mediante matrices.

[ ]
1.677e+8 18080425 1978405 220825 25333
18080425 1978405 220825 25333 3025
mp= 197840 220825 25333 3025 385
220825 25333 3025 385 55
25333 3025 385 55 10

[]
a
b
vp= c
d
e
 [ ]
1070375.5
120750.1
tip= 14076.7
1725.1
230

Entonces mp∗vp=tip , obteniendo mediante algebra, que:

vp=inv ( mp )∗tip

Es decir

([ ]) [ ]
1.677e+8 18080425 1978405 220825 25333 1070375.5
18080425 1978405 22082525333 3025 120750.1
vp=inv 197840 220825 25333 3025 385 ∗ 14076.7
220825 25333 3025 385 55 1725.1
25333 3025 385 55 10 230

Al resolver el sistema a través de Scilab se obtiene que:

a=0.0069202

b=−0.1423174

c=1.2545309

d=−0.0583528

e=0.5416667

Por lo tanto, se obtiene que la función calculada es:

4 3 2
y=0.0069202 x −0.1423174 x +1.2545309 x −0.0583528 x+ 0.5416667

Grafica 5. Porcentaje de humedad vs Tiempo (h). Datos observados, datos calculados y

función.
Finalmente se procede a realizar el cálculo de errores para los valores calculados respecto
de los valores observados.

Tabla 8. Valores Observados vs Valores Calculados

n y0 yc
A B
1 1,6 1,60245 457,8552487 5,99075E-06
2 4,7 4,41527 345,3922598 0,081072255
3 7,7 8,37535 213,8803964 0,456097217
4 13 13,04400 99,12197582 0,001935824
5 19 18,14860 23,53607226 0,724880257
6 24 23,58263 0,339462378 0,174194378
7 28 29,40565 41,03239036 1,975860356
8 36 35,84330 164,9503035 0,024555517
9 44 43,28730 411,574379 0,507946992
10 52 52,29546 858,223684 0,087293657
∑ 230 230 2615,906172 4,033842444

ymed=
∑ y 0 230
= =23
n 1

A=∑ ( y − ymed) =¿2615,906172

c 2
B=∑ ( y − y ) =4,033842444
c 0 2

2 A 2615,906172
r= = =0.9984603
A+ B 2615,906172+ 4,033842444

r 2=0.9984603[ Ajuste Polinomio de orden 4 ]

6. Análisis de Resultados
Se puede observar que al utilizar los métodos entre más sean el número de pasos menor
será el error efectuado en el problema a solucionar y que para dichos problemas no tienen
que dar unos valores específicos iniciales también que la grafica que mejor se efectúa y
toma la mayor cantidad de números es el de polinomio de orden 4 aunque la parabólica
también coge muchos puntos de la gráfica, pero con mayor exactitud es la de orden
superior
7. Conclusiones
De manera que, a través del conjunto de datos pares y ordenados inicialmente
presentados, se halló una función lineal, parabólica, cubica y de orden superior
que se ajustaron lo mejor posible, con mínimas diferencias y errores entre los
valores observados y calculados, los cuales nos proporcionaron además una
demostración visual entre los puntos como se puede también observar al unir
todas las graficas algunas toman punto pero no los suficientes
8. Bibliografía
No se hizo consulta exteriores
9. Anexos: Códigos
Códigos estudiantiles

20201186235