matlab中提取不规则圆周上几个数据点并进行最小二乘拟合圆代码？

以下是一个示例代码，可以提取不规则圆周上的数据点并进行最小二乘拟合圆：

% 生成一个不规则圆周的数据点
theta = linspace(0, 2*pi, 100);
x = 2*cos(theta) + randn(size(theta))*0.1;
y = 2*sin(theta) + randn(size(theta))*0.1;

% 选择一部分数据点进行拟合
idx = 1:10:100;
x_fit = x(idx);
y_fit = y(idx);

% 构造矩阵A和向量b
A = [x_fit', y_fit', ones(size(x_fit'))];
b = -x_fit'.^2 - y_fit'.^2;

% 最小二乘拟合圆
c = pinv(A)*b;
xc = -0.5*c(1);
yc = -0.5*c(2);
R = sqrt((xc^2 + yc^2) - c(3));

% 绘制结果
figure;
plot(x, y, '.', 'MarkerSize', 10);
hold on;
theta_fit = linspace(0, 2*pi, 100);
x_circle = xc + R*cos(theta_fit);
y_circle = yc + R*sin(theta_fit);
plot(x_circle, y_circle);
axis equal;

上述代码首先生成一个不规则圆周的数据点，然后选择其中一部分数据点进行拟合。拟合过程使用最小二乘法，构造一个线性方程组，其中矩阵A和向量b的形式为：

$$ A = \begin{bmatrix} x_1 & y_1 & 1 \ x_2 & y_2 & 1 \ \vdots & \vdots & \vdots \ x_n & y_n & 1 \end{bmatrix}, \quad b = \begin{bmatrix} -(x_1^2 + y_1^2) \ -(x_2^2 + y_2^2) \ \vdots \ -(x_n^2 + y_n^2) \end{bmatrix} $$

其中$n$为选择的数据点数量。解出线性方程组的解向量$c$后，圆心坐标为$(-0.5c_1, -0.5c_2)$，半径为$\sqrt{(c_1^2 + c_2^2) - c_3}$。最后绘制出拟合的圆和原始数据点的散点图