圆拟合程序matlab

一个高效的matlab代码，用来拟合图像中的圆。

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