function?[xcycRa]=circlefit（xy）
%?最小二乘法圓的擬合函數
%?CIRCLEFIT?fits?a?circle?in?xy?plane
%?x^2+y^2+a（1）*x+a（2）*y+a（3）=0
%?（xcxy）:圓心坐標；R：半徑；a：方程系數

n=length（x）;
xx=x.*x;
yy=y.*y;
xy=x.*y;

A=[sum（x）?sum（y）?n;sum（xy）?sum（yy）?sum（y）;sum（xx）?sum（xy）?sum（x）];
B=[-sum（xx+yy）;-sum（xx.*y+yy.*y）;-sum（xx.*x+xy.*y）];
a=A\B;
xc?=?-0.5*a（1）;
yc?=?-0.5*a