function?[abr]=circfit（xy）????????????????????????????%??????????ab為圓心坐標，r為半徑
n=length（x）;
c1=sum（x）;c2=sum（y）;c3=sum（x.*x）;c4=sum（y.*y）;?????????????%最小二乘估計
c5=sum（x.*y）;c6=sum（x.*x.*x）;c7=sum（y.*y.*y）;
c8=sum（x.*y.*y）;c9=sum（x.*x.*y）;
d1=2*（c3-c1.*c1/n）;d2=2*（c5-c1*c2/n）;d3=c6+c8-c1*（c3+c4）/n;
d4=2*（c4-c2*c2/n）;d5=c9+c7-c2*（c3+c4）/n;
a=（d3*d4-d2*d5）/（d1*d4-d2*d2）;
b=（d1*d5-d2*d3）/（d1*d4-d2*d2）;
r=sqrt（（（c3+c4-2*a*c1-2*b*c2）/n）+a*a+b*b）;

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

?????文件????????470??2009-12-01?20:18??light_spot_center\circfit.m

?????文件????????581??2009-12-03?09:37??light_spot_center\gray1.m

?????文件????????729??2009-12-03?09:37??light_spot_center\gray11.m

?????文件??????28603??2009-06-17?17:50??light_spot_center\guangban23.jpg

?????文件???????1703??2009-11-25?17:13??light_spot_center\optimumthreshold.m

?????目錄??????????0??2012-07-06?09:27??light_spot_center

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

????????????????32086????????????????????6