ce6_1b.m

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

clear all
n = 100;
z = zeros(n,2);
x = [2.4,4.32]';
v = zeros(2,n);
for i = 1:n;
    v(:,i) = [0.05,0;0.6*0.015,0.015*0.8]*randn(2,1);
    z(i,:) = ([0,i;0.5,i/10]*x + v(:,i))';
end

[a,b] = line_fit(z,n)

plot(z(:,1),z(:,2),'+R')
hold
tempX = [0,450];
tempY = a.*tempX +b;
line(tempX,tempY)