clc;
close all;
x = [2064.12,1932.77,1885.74,1718.55,1657.15,1569.56,1471.99,1418.19,1262.71];
y = [19.12685537,21.45685537,21.92685537,25.67685537,27.85685537,31.47685537,37.12685537,40.49685537,48.75685537];
plot(x,y,'r*')
p = polyfit(x,y,2)
y1 = polyval(p,x)
hold on;
plot(x,y1,'b.')
poly2sym(p)
结果:
y1 =
19.3798 20.8708 21.7292 26.1662 28.3386 31.9424 36.6557 39.5696 49.2495
ans =
(356785675701145*x^2)/9223372036854775808 - (5979427441899477*x)/36028797018963968 + 6936059580474597/35184372088832
拟合效果: