func4p4.m

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

%func4p4.m二次曲线拟合
x=[2,3,4,5,7,8,10,11,14,15,16,18,19]
y=[106.42,108.26,109.58,109.50,110.00,109.93,110.49,110.59,110.60,110.90,110.76,111.00,111.20]
p=polyfit(x,y,2)
xx=linspace(2,19,100)
yy=polyval(p,xx)
plot(xx,yy,'-',x,y,'or')
f=polyval(p,x)
disp(abs(norm(y,inf)-norm(f,inf)))
%此时残量(residue)为:0.1338
%多项式各项系数：c2=-0.0201    c1=0.6069  c0=106.4826