对于空间中的三维离散坐标点，如何拟合空间圆得到空间圆心坐标和半径，请给出具体方法并用matlab代码实现一下。

拟合空间圆的问题可以转化为拟合平面和求解圆心和半径的问题。下面是一个具体的方法：

1. 选取三个不共线的点，将它们看做是圆上的三个点。

2. 计算这三个点所在平面的法向量。

3. 将所有点投影到该平面上，得到它们在平面上的二维坐标。

4. 在平面上拟合一个圆，得到圆心坐标和半径。

5. 将圆心坐标和法向量反投影到空间中，得到圆心坐标。

下面是Matlab代码实现：

% 生成测试数据
x = rand(100, 1);
y = rand(100, 1);
z = rand(100, 1);

% 选取三个不共线的点
p1 = [x(1), y(1), z(1)];
p2 = [x(2), y(2), z(2)];
p3 = [x(3), y(3), z(3)];

% 计算平面法向量
v1 = p2 - p1;
v2 = p3 - p1;
n = cross(v1, v2);
n = n / norm(n);

% 投影到平面上
d = dot(n, p1);
P = [x, y, z];
P = P - (d + P * n')' * repmat(n, size(P, 1), 1);
P2D = P(:, 1:2);

% 拟合圆
[x0, y0, r] = fit_circle(P2D(:,1), P2D(:,2));

% 反投影到空间中
p0 = [x0, y0, 0];
p0 = p0 * n' / n;
p0 = p0 + (d / norm(n)) * n;

其中，fit_circle 函数的实现可以参考以下链接：https://ww2.mathworks.cn/matlabcentral/fileexchange/22643-circle-fit/content/fit_circle.m