Matlab中正方形平面挖正方形孔的代码异常问题排查

问题根源 #

你遇到的两个问题本质是掩码创建逻辑的精度与边界判断方式错误：

• 角落斜角：大概率用了欧氏距离（如sqrt((x-x0)^2 + (y-y0)^2)）来筛选孔区域，这种方式会生成圆形过渡边缘，而非严格直角；
• 壁厚偏差：物理坐标与网格索引的转换未对齐，浮点计算的累积误差导致实际偏移量偏离预期值。

修复方案与更优实现 #

以下两种方法都能彻底解决问题，优先推荐基于网格索引的精确实现：

方法1：基于网格索引的精确矩形掩码（无浮点误差）

直接通过行列索引定义区域，完全对齐网格点，确保直角和壁厚绝对精确：

a = 3;          % 外正方形物理尺寸
b = 2;          % 内孔物理尺寸
N = 101;        % 网格点数量（建议取奇数，中心对齐；N越大边缘越平滑）

% 创建物理坐标网格
[x, y] = meshgrid(linspace(0, a, N), linspace(0, a, N));

% 计算孔的边界索引（将物理偏移转换为网格索引）
offset_phys = (a - b)/2;
grid_step = a/(N-1);
idx_start = round(offset_phys / grid_step) + 1;
idx_end = round((a - offset_phys) / grid_step) + 1;

% 生成掩码：1为保留区域，0为孔区域
mask = ones(N, N);
mask(idx_start:idx_end, idx_start:idx_end) = 0;

% 可视化验证
figure;
imagesc(x, y, mask);
axis equal tight;
xlabel('X'); ylabel('Y');
title('带精确正方形孔的网格掩码');

方法2：基于物理坐标的严格范围判断

如果必须用物理坐标，直接用矩形的x/y范围逻辑与，避免距离判断：

a = 3;
b = 2;
N = 101;

[x, y] = meshgrid(linspace(0, a, N), linspace(0, a, N));

offset = (a - b)/2;
% 严格判断：仅保留x/y在孔外的区域
mask = (x <= offset | x >= offset + b) & (y <= offset | y >= offset + b);

% 可视化
figure;
imagesc(x, y, mask);
axis equal tight;

关键注意事项 #

• 网格点数量N：取奇数可让中心网格点对齐正方形中心，避免偏移；N越大，边缘过渡越平滑，壁厚精度越高；
• 避免欧氏距离判断：这类逻辑会生成圆形边缘，必然导致角落斜角；
• 壁厚验证：当a=3、b=2时，两种方法的偏移量均为精确的0.5，完全符合预期。

内容的提问来源于stack exchange，提问作者Fish1996