利用MATLAB进行直线拟合的程序

源程序：

% 初始数据

x=[0.00310,0.00313,0.00316,0.00319,0.00322,0.00326,0.00329,0.00332];

y=[9.985,9.893,9.810,9.729,9.646,9.581,9.507,9.451];

% 拟合直线y=p(1)*x+p(2)

p=polyfit(x,y,1);

% 显示拟合前后直线，其中方框为拟合前数据，直线为拟合后的直线

yy=polyval(p,x);

plot(x,y,'s',x,yy)

%拟合的直线方程

poly2sym(p,'x')

结果：

ans =

4648504458470549*x-396243168544487872

>> vpa(ans,8)

ans =

.46485045e16*x-.39624317e18