线性拟合：对于y=a*x+b的形式

    a=(N*ΣXi*Yi-ΣXi*ΣYi)/(N*ΣXi*Xi-(ΣXi)2)

    b=((ΣXi*Xi)*(ΣYi)-(ΣXi)*(ΣXi*Yi))/(N*ΣXi*Xi-(ΣXi)2)

    利用Matlab自带函数可实现：

    拟合函数：pn=polyfit(x,y,n)  返回pn系数向量，降阶排列，n为阶数

    函数：yy=polyval(pn,x) pn为降阶排列的多项式系数，x为向量或者矩阵，返回yy：将x带入pn后计算得到的向量或者矩阵。

%初始数据
x=[1 2 3 4];
y=[1.1 2.2 2.7 3.8];
%获得线性拟合系数
pn=polyfit(x,y,1);

yy=polyval(pn,x);
subplot(1,2,1);
plot(x,y);
subplot(1,2,2);
plot(x,yy);

结果：pn=0.86,0.3