zhixian.m

这是根据点拟合圆和直线的matlab程序

function [a,b,res,time]=zhixian(x,y)
c1=clock;
n=length(x);
x1=sum(x);
y1=sum(y);
x2=sum(x.*x);
y2=sum(y.*y);
xy=sum(x.*y);
a=(n*xy-x1*y1)/(n*x2-x1*x1);
b=(y1-a*x1)/n;
res=sum((a*x+b-y).*(a*x+b-y))/n;
plot(x,y,'*',x,a*x+b)
c2=clock;
time=c2-c1;
time=[0 0 86400 3600 60 1]*time';