若要将得到的数组进行圆的拟合,改用什么程序编写呢?请大家帮帮我吧 ,最好有具体的程序,谢谢..

一段文章，供参考：
Matlab中如何作圆回归？
:#Peter Boettcher (boettcher@ll.mit.edu),2002/5/16, comp.soft-sys.matlab#

Q5.5: How can I fit a circle to a set of XY data?
=================================================

An elegant chunk of code to perform least-squares circle fitting was
written by Bucher Izhak and has been floating around the newgroup for
some time. The first reference to it that I can find is in:

function [xc,yc,R,a] = circfit(x,y)
%CIRCFIT Fits a circle in x,y plane
%
% [XC, YC, R, A] = CIRCFIT(X,Y)
% Result is center point (yc,xc) 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));

Tom Davis provided a more sophisticated approach that works for more
cases in and Code included.

一个实例：
x=[38.849, 37.706, 31.188, 20.238, 6.8660, -6.5940, -18.056, -26.135, -30.328, -30.814, -28.035, -22.542, -14.782, -5.2590, 5.4750, 16.587, 26.893, 34.852];
y=[ 0., 13.724, 26.107, 35.053, 38.934, 37.402, 31.276, 21.931, 11.039, .100e-2, -10.204, -18.915, -25.603, -29.825, -31.049, -28.729, -22.565, -12.684];
b=[35.06010729, 4.016838807, .2587265862e-2]
x1=linspace(min(x)-1,max(x)+.5,120);
y1=sqrt(b(1)*b(1)-(x1-b(2)).^2)+b(3)+4;
plot(x,y,'o','markerfacecolor','k','markeredgecolor','k','markersize',8)
hold on
plot(x1,y1,'b-','linewidth',3)
y2=-sqrt(b(1)*b(1)-(x1-b(2)).^2)+b(3)+4;
plot(x1,y2,'b-','linewidth',3)
plot(b(2),4+b(3),'or','markerfacecolor','r','linewidth',2)
axis tight
在本例中，b(2)=4.0168,4+b(3)=4.0026即为圆心位置，b(1)=35.06为半径。
若能提供数据，也可进行类似拟合。