拟合以下数据
[0.00016 0.00031 0.00049 0.00067 0.00082];
[0.03025 0.0605 0.09075 0.125 0.15625];
代码：
X=[0.00016 0.00031 0.00049 0.00067 0.00082];
Y=[0.03025 0.0605 0.09075 0.125 0.15625];
n=5;
x2=sum(X.^2); % 求Σ(xi^2)
x1=sum(X); % 求Σ(xi)
x1y1=sum(X.Y); % 求Σ(xiyi)
y1=sum(Y); % 求Σ(yi)
a=(nx1y1-x1y1)/(nx2-x1x1); %解出直线斜率b=(y1-ax1)/n
b=(y1-ax1)/n; %解出直线截距
%作图
% 先把原始数据点用蓝色十字描出来
figure
plot(X,Y,'+');
hold on
grid on
% 用红色绘制拟合出的直线
px=linspace(0,0.001,45);
py=a*px+b;
plot(px,py,'r');
运行结果
