function?[xcycRa]?=?circfit（xy）
%CIRCFIT?Fits?a?circle?in?xy?plane
%
%?[XC?YC?R?A]?=?CIRCFIT（XY）
%?Result?is?center?point?（ycxc）?and?radius?R.A?is?an
%?optional?output?describing?the?circle’s?equation:
%
%?x^2+y^2+a（1）*x+a（2）*y+a（3）=0
%?by?Bucher?izhak?25/oct/1991
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?=?-.5*a（1）;
yc?=?-.5*a（2）;
R?=?sqrt（（a（1）^2+a（2）^2）/4-a（3））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????531??2008-10-30?13:33??圓的擬合\circfit.m

?????文件???????1017??2008-10-30?16:12??圓的擬合\yuannihe.m

?????文件?????????50??2008-12-05?15:49??圓的擬合\說明.txt

?????目錄??????????0??2010-01-30?17:15??圓的擬合

-----------?---------??----------?-----??----

?????????????????1598????????????????????4