clear all
% 生成离散点
t = [0 1 2 3 4 5 6 7 8 9];
y = [0 2 4 7 8 9 12 14 15 18];
% 构造参数矩阵
A = ones(10,2);
A(:,2) = t';
% 参数初值
a = 0; %平均速度
b = 0; %初始位置
% 误差值
l = y - (b + a*t);
% 解算参数
for i = 1:10
dX = inv(A'*A)*A'*l';
b = b + dX(1);
a = a + dX(2);
l = y - (b + a*t);
end
% 成图，红色‘+’为生成的离散点，黑色直线为最小二乘拟合直线
plot(t,y,'r+',t,b+a*t,'k')
disp(a)
disp(b)
