function [P,newxy] = myline(x,y)

% 調(diào)用形式：
%       [P]=myline(x,y)
% 功能：
%       考慮到壞點(diǎn)的剔除的最小二乘法的直線擬合 f(x) = A + Bx
%       其中：調(diào)用子程序submyline1, submyline2
% 輸入：
%        -- x    行向量，實(shí)驗(yàn)點(diǎn)數(shù)據(jù)橫坐標(biāo)
%        -- y    行向量，實(shí)驗(yàn)點(diǎn)數(shù)據(jù)縱坐標(biāo)
% 輸出：
%        -- P    行向量，依次代表：A,B,C,S,R,dA,dB,dC
% 袁怡圃 2003/04/24

a=x;
b=y;
[A,B,S]=submyline1(a,b);

N=length(x);
for j=1:N
    k=1; 
    for i=1:N
        if abs(y(i)-A-B*x(i))<2*S
            u(k)=x(i);
            v(k)=y(i);
            k=k+1;
        end
    end
    u=u(1:k-1);
    v=v(1:k-1);
    [A,B,S]=submyline1(u,v);
    x=u;
    y=v;
    N=length(x);
end 
newxy(1,:) = u;
newxy(2,:) = v;
% xmin=min(a);
% xmax=max(a);
% x=xmin:0.01:xmax;
% y=B*x+A;
% plot(a,b,'*r',u,v,'.g',x,y,'b');
% disp('A B C S R dA dB dC分別為：')
[C,R,dA,dB,dC]=submyline2(u,v);
P=[A B C S R dA dB dC];