基于最小二乘法的工程测量圆拟合的Matlab代码，并画出散点和拟合的圆，并求出各个点到该圆圆心的距离

下面是基于最小二乘法的工程测量圆拟合的Matlab代码：

% 输入数据
x = [1.5 2.3 3.1 3.8 4.5 5.2 6.1 6.8 7.6 8.3];
y = [0.7 1.6 2.7 3.5 4.5 5.5 6.7 7.4 8.2 9.1];

% 定义目标函数
fun = @(a,x) (x(1,:)-a(1)).^2 + (x(2,:)-a(2)).^2 - a(3)^2;

% 初始化参数
a0 = [0,0,1];

% 使用最小二乘法进行拟合
a = lsqcurvefit(fun,a0,[x;y],zeros(size(x)),[],[],[],[],optimset('Display','off'));

% 画出散点和拟合的圆
plot(x,y,'o');
hold on;
th = linspace(0,2*pi,100);
xc = a(1) + a(3)*cos(th);
yc = a(2) + a(3)*sin(th);
plot(xc,yc);
title('圆拟合');
xlabel('x');
ylabel('y');

% 计算各个点到该圆圆心的距离
d = sqrt((x-a(1)).^2 + (y-a(2)).^2);
disp(['各个点到该圆圆心的距离为：', num2str(d)]);

运行以上代码，可以得到散点和拟合的圆，以及各个点到该圆圆心的距离。