clear
clc
t = 1 : 16;
y = [1 4 9 16 25 36 49 64 81 100 121 144 169 800 900 1000];
plot(t, y, '*')
p = polyfit(t, y, 10); % p是拟合的多项式系数向量
t1 = 1 : 0.1 : 16;
y1 = polyval(p, t1);
figure
plot(t1, y1);
figure
plot(t, y, '*', t1, y1, 'red');
结果为: