func4p3.asv

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

x=linspace(1960,1968,9)
y=[]
p=polyfit(x,y,2)
xx=linspace(0.5,3,100)
yy=polyval(p,xx)
plot(xx,yy,'-')
f=polyval(p,x)
disp(abs(norm(y,inf)-norm(f,inf)))
%此时残量(residue)为: 0.0950
%多项式各项系数：c2=0.5614，c1=0.8287，c0=1.1560
%如提高为用三次多项式拟合，residue=0.0517,
%更为精确；一次：residue=0.3729；四次：residue=0.0238,但波动较大
%can also see the figure