curvefit03.m

数值分析与计算的PPT幻灯片与MAT源码

%curvefit03.m

x0=0:0.1:1;
y0=[-.447,1.978,3.11,5.25,5.02,...
     4.66,4.01,4.58,3.45,5.35,9.22];
plot(x0,y0,'or')              %画出离散点
pause
n=7;
a=polyfit(x0,y0,n);           %拟合多项式
p=poly2sym(a)                 %求出多项式的表达式
xi=0:0.01:1;                  %细分区间
yi=polyval(a,xi);             %求多项式的值
hold on
plot(xi,yi,'-b','MarkerSize',20)%绘制拟合曲线图
xlabel('x') 
legend('离散点','拟合曲线',2)