Matlab 程序：圆形工件质量检验模型

以下是一个基本的 Matlab 程序框架，可以实现对每个工件的 36 个数据进行圆心和半径的计算，并判断是否为正品。

% 读取数据，假设数据存储在 data 文件夹下的 txt 文件中
data = load('data/workpiece1.txt');

% 计算圆心坐标
x = data(:,1);
y = data(:,2);
A = [x y ones(size(x))];
b = -x.^2 - y.^2;
c = A\b;
xc = -0.5 * c(1);
yc = -0.5 * c(2);
r = sqrt((xc - x(1))^2 + (yc - y(1))^2); % 计算半径

% 判断是否为正品
tolerance = 0.1; % 允许的误差
if all(abs(sqrt((x - xc).^2 + (y - yc).^2) - r) <= tolerance)
    disp('正品');
    disp(['圆心坐标为(', num2str(xc), ', ', num2str(yc), ')']);
else
    disp('次品');
end

需要注意的是，程序中假设所有的工件数据都存储在data文件夹下，且每个工件的数据存储在一个以workpiece为前缀的 txt 文件中（如workpiece1.txt、workpiece2.txt等）。实际使用时需要根据具体情况进行修改。另外，程序中的圆心和半径计算公式是基于最小二乘法的，可能存在精度问题，需要根据实际情况进行调整。