fitline2dpoints.m

来自「Line & Curve Fitting To 2D Points」· M 代码 · 共 20 行

function fitLine2Dpoints(X,n)
[x,y,beta,alpha] = fitLine(X,n);
plot(X(:,1),X(:,2),'o',x,y);
i1 = 0;
i2 = 0;
for i = 1:size(X,1)
    if X(i,2) > (beta * X(i,1) + alpha)
        i1 = i1 + 1;
        X1(i1,:) = X(i,:);
    else
        if X(i,2) < (beta * X(i,1) + alpha)
            i2 = i2 + 1;
            X2(i2,:) = X(i,:);
        end
    end
end
[x1,y1,beta1,alpha1] = fitLine(X1,size(X1,1));
[x2,y2,beta2,alpha2] = fitLine(X2,size(X2,1));
plot(X(:,1),X(:,2),'o',x,y,'g',x1,y1,x2,y2)
end