func4p3.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 14 行

%func4p3.m,求拟合线性多项式，预测人口
x=linspace(1960,1968,9)
y=[3.3918,3.4213,3.4503,3.34698,3.4763,3.492,3.5133,3.5322,3.5505]
p=polyfit(x,y,1)
xx=linspace(1960,1968,100)
yy=polyval(p,xx)
plot(xx,yy,'-')
f=polyval(p,x)
y2000=exp(polyval(p,2000))
disp(abs(norm(y,inf)-norm(f,inf)))
%此时残量(residue)为:  0.0041,拟合较好
%多项式各项系数：c1=0.0206  c0=-37.0770
%预测当时人口为： 67.1524