title('拉格朗日插值点的函数值图像'); subplot(325)

plot(xp,yf,'b+',xp,yp,'g-',x,y,'r--',x1,y1,'y-','LineWidth',2); title('样条、拉格朗日、最小二乘法与原始数据对比图

');

方法二：

因为不知道其函数形式，可先plot而后确定使用哪种模型比较合适。 函数图形程序：

x=[2 3 4 5 6 7 8 9 10 11 12 13 14 15 16];

y=[6.42 8.2 9.58 9.5 9.7 10 9.93 9.99 10.49 10.59 10.6 10.8 10.6 10.9 10.76];

plot(x,y)

与指数函数趋势类似（但是趋势不同，一个趋于递减，一个趋于递增，使用倒数），故拟合模型为：

y ae

b

x

x=[2 3 4 5 6 7 8 9 10 11 12 13 14 15 16];

y=[6.42 8.2 9.58 9.5 9.7 10 9.93 9.99 10.49 10.59 10.6 10.8 10.6 10.9 10.76]; plot(x,y,'r')

myfunc=inline('beta(1)*exp(beta(2)./x) ','beta','x'); beta=nlinfit(x,y,myfunc,[0.5 0.5]); xt=2:0.01:16;

yt= beta(1)*exp(beta(2)./xt); hold on

plot(xt,yt,'b--');