function [A,B,S]=submyline1(x,y)

% [A,B,S]=submyline1(x,y)  最小二乘法的直線擬合 f(x) = A + Bx
% 輸入：   -- x  1xN 的向量
%          -- y  1xN 的向量 (x,y)是實驗點數(shù)據(jù)
% 輸出：   -- A  擬合直線的截距
%          -- B  擬合直線的斜率
%          -- S  擬合結(jié)果的剩余標準誤差
%袁怡圃 2003/04/24

xmean=mean(x);
ymean=mean(y);
sumx2=(x-xmean)*(x-xmean)';
sumxy=(x-xmean)*(y-ymean)';
B=sumxy/sumx2;
A=ymean-B*xmean;
N=length(x);
Q=0;
for i=1:N
    Q=Q+(y(i)-A-B*x(i))^2;
end
S=sqrt(Q/(N-2));