curvefit02.m

来自「数值分析与计算的PPT幻灯片与MAT源码」· M 代码 · 共 24 行

% curvefit02.m

t=[1:16];
y=[4 4.6 8 8.4 9.28 9.5 9.7 9.86 10 10.2 ...
   10.32 10.30 10.24 10.18 10.00 9.40];
plot(t,y,'.k','markersize',20)   %描点
pause

a=polyfit(t,y,2)              %求出二次多项式的系数
p=poly2sym(a)                   %给出多项式的符号表达式
xi=[0:0.1:16];                %将区间[1,16]160等分
yi=polyval(a,xi);             %求多项式在分点的值
hold on
plot(xi,yi,'b-','markersize',20)%绘制拟合曲线图
pause

yii=interp1(t,y,xi,'liner');    %求线性插值函数值
hold on
plot(xi,yii,'r-','markersize',20)%绘制插值曲线图
x0=15.5;
y0=polyval(a,x0)
hold on
plot(x0,y0,'*r')
legend('已知数据点','拟合曲线','样条插值曲线','计算值',0)