matlab中已知最小二乘拟合圆圆心和半径求铆钉孔垂直度代码？

以下是基于最小二乘拟合圆的铆钉孔垂直度计算的MATLAB代码：

% 已知圆心和半径 xc = 0; % 圆心x坐标 yc = 0; % 圆心y坐标 r = 10; % 圆半径

% 生成随机的铆钉孔坐标 n = 50; % 铆钉孔数量 theta = linspace(0, 2pi, n+1); % 角度范围 theta = theta(1:end-1); % 去掉重复的最后一个 x = xc + rcos(theta) + 0.1randn(1,n); % x坐标加入随机噪声 y = yc + rsin(theta) + 0.1*randn(1,n); % y坐标加入随机噪声

% 最小二乘拟合圆 A = [x', y', ones(n,1)]; % 设计矩阵 b = -(x.^2 + y.^2); % 右侧向量 c = A\b; % 系数向量

% 圆心和半径 xc_fit = -c(1)/2; yc_fit = -c(2)/2; r_fit = sqrt(xc_fit^2 + yc_fit^2 - c(3));

% 绘制原始数据和拟合圆 figure; plot(x, y, 'bo'); hold on; theta_fit = linspace(0, 2pi, 100); x_fit = xc_fit + r_fitcos(theta_fit); y_fit = yc_fit + r_fit*sin(theta_fit); plot(x_fit, y_fit, 'r-'); axis equal;

% 计算铆钉孔垂直度 d = abs(sqrt((x-xc_fit).^2 + (y-yc_fit).^2) - r_fit); % 铆钉孔到拟合圆的距离 verticality = max(d) - min(d); % 铆钉孔垂直度

disp(['铆钉孔垂直度为：', num2str(verticality)])