matlab拟合圆的方程代码

以下是使用MATLAB拟合圆的方程代码示例：

% 生成示例数据
theta = linspace(0, 2*pi, 100); % 角度范围
r = 5; % 圆的半径
x = r*cos(theta); % x坐标
y = r*sin(theta); % y坐标

% 添加噪声
noise = randn(size(x))*0.1;
x = x + noise;
y = y + noise;

% 使用最小二乘法拟合圆
A = [x' y' ones(size(x))']; % 构建系数矩阵
b = -x'.^2 - y'.^2; % 构建常数向量
coeff = A\b; % 线性回归求解

% 提取圆心和半径
center_x = -coeff(1)/2;
center_y = -coeff(2)/2;
radius = sqrt(center_x^2 + center_y^2 - coeff(3));

% 绘制拟合结果
figure;
plot(x, y, 'ro'); % 绘制原始数据点
hold on;
theta_fit = linspace(0, 2*pi, 100); % 拟合圆的角度范围
x_fit = center_x + radius*cos(theta_fit); % 拟合圆的x坐标
y_fit = center_y + radius*sin(theta_fit); % 拟合圆的y坐标
plot(x_fit, y_fit, 'b-'); % 绘制拟合圆

% 添加标题和标签
title('Circle Fitting');
xlabel('x');
ylabel('y');
axis equal; % 设置坐标轴比例相同
legend('Original Data', 'Fitted Circle');

这段代码生成了一个半径为5的圆，并添加了一些噪声。然后使用最小二乘法拟合圆，并绘制拟合结果。你可以根据需要修改半径、噪声等参数。