fitcurve2dpoints.m

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

function fitCurve2Dpoints(X,err,rng)
 i = size(X,2);
 x = 0:0.01:rng;
 while(true)
 Mat = makeMat(X,i);
 A = Mat\X(:,2);
 y = polyval(A,x);
 diffy = (X(:,2) - polyval(A,X(:,1))).^2;
 meanErry = sum(diffy)/size(X,1);
 if meanErry < err
    break;
 end
 i = i + 1;
 end
plot(X(:,1),X(:,2),'o',x,y,'g')
%axis ([ 0 rng 0 0.5*rng^2])
end