circfit.m

本人编写的基于最小二乘法的圆拟合程序

function   [xc,yc,R,a] = circfit(x,y)

%   [xc yx R] = circfit(x,y)

%  圆心为 (yc,xc) 半径为 R

%   x^2+y^2+a*x+b*y+c=0

%  用最小二乘法，xc=-0.5a(1),yc=-0.5a(2)

    x=x(:); y=y(:);
   abc=[x y ones(size(x))]\[-(x.^2+y.^2)];
   a=abc(1);
   b=abc(2);
   c=abc(3);
   xc = -.5*a;
   yc = -.5*b;
   R  =  sqrt((a^2+b^2)/4-c)