line_fit.m

来自「给定带有高斯噪声的三组样本,对其进行最小二乘拟合,并且绘图,求出误差协方差」· M 代码 · 共 11 行

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;