line_fit.m

利用最小二乘的方法进行直线拟合,生成100组带有高斯噪声的数据进行模拟

function [a,b]=line_fit(z,n)
meanX = mean(z(:,1));
meanY = mean(z(:,2));
upTemp = 0;
lowTemp = 0;
for i = 1:n
    upTemp = upTemp + (z(i,1)-meanX)*z(i,2);
    lowTemp = lowTemp + (z(i,1)-meanX)^2;
end
a = upTemp/lowTemp;
b = meanY - a*meanX;