拟合圆matlab代码

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