y=191,321,442,565,686,819,930,1032,1153,1252;%因变量xmean=mean(x);ymean=mean(y);sumx2=(x-xmean)*(x-xmean);sumxy=(y-ymean)*(x-xmean)a=sumxy/sumx2;%解出直线斜率a(即传感器灵敏度)b=ymean-a*xmean;%解出直线截距bz=(a*(x(1,10)+b-(y(1,10)/(y(1,10);%“10”是自变量的个数,z为非线性误差(即线性度)abz%作图,先把原始数据点用蓝色十字描出来figureplot(x,y,+); hold on% 用红色绘制拟合出的直线px=linspace(0,6,50);%(linspace语法(从横坐标负轴起点0画到横坐标正轴终点6,50等分精度)py=a*px+b;plot(px,py,r运行结果: a =236.9818 b =87.4000另一种简单一点的方法:p=polyfit(x,y,1);p p = 236.9818 87.4000