Departamento de
produção e Sistemas – Métodos numéricos
Método dos
mínimos quadrados para a função do segundo grau em Matlab
>> x=[46.000000;
20.000000;
52.000000;
30.000000;
57.000000;
25.000000;
28.000000;
36.000000;
57.000000;
44.000000;
31.000000;
31.000000;
52.000000;
23.000000;
60.000000;
48.000000;
34.000000;
51.000000;
50.000000;
34.000000;
46.000000;
23.000000;
37.000000;
40.000000;
30.000000]
x =
46
20
52
30
57
25
28
36
57
44
31
31
52
23
60
48
34
51
50
34
46
23
37
40
30
>> y=[354.000000;
190.000000;
405.000000;
263.000000;
451.000000;
302.000000;
288.000000;
385.000000;
402.000000;
365.000000;
209.000000;
290.000000;
356.000000;
254.000000;
395.000000;
434.000000;
220.000000;
374.000000;
308.000000;
220.000000;
311.000000;
181.000000;
274.000000;
303.000000;
244.000000]
y =
354
190
405
263
451
302
288
385
402
365
209
290
356
254
395
434
220
374
308
220
311
181
274
303
244
>> p= polyfit (x,y,2)
p =
0.0106
4.5554 113.6752
>> f= polyval (p,x)
f =
345.7159
209.0354
379.2984
259.9042
407.8683
234.2040
249.5603
291.4458
407.8683
334.6918
265.1080
265.1080
379.2984
224.0728
425.2654
356.8250
280.8469
373.6481
368.0192
280.8469
345.7159
224.0728
296.7772
312.8988
259.9042
>> r=sum((y-f).^2)
r =
4.5756e+004