MATLAB网格网格错误

Question

我在我的一个循环调用的matlab函数中使用meshgrid，特定的代码行如下（对不起，我不能在这里粘贴完整的代码）

rect
[m,n] = meshgrid(rect(1):rect(3),rect(2):rect(4));
size(m)

现在我正在打印 rect 和 size(m)，在输出中显示为 ans 以及 meshgrid 代码。我在最后两次迭代下面发布

1) 第一个给出 size(m) 为 46x48 的尺寸是正确的尺寸，并且不会在我稍后进行的矩阵加法中引发错误。

2) 第二个显示 size(m) 为 46x47 的尺寸是错误的，尽管 rect 值表明不是这样。

rect =

  119.9290   69.3807  166.9290  114.3807


ans =

    46    48


rect =

  120.0508   69.5380  167.0508  114.5380


ans =

    46    47

Error using  + 
Matrix dimensions must agree.

现在我取同样的数字，在命令窗口做一个meshgrid，结果正确如下

>> [lm,kl] = meshgrid(120.0508:167.0508,69.5380:114.5380);
>> size(lm)

ans =

    46    48

这个异常的原因可能是什么？

Answer 1

您可能遇到浮点精度错误。为了有一个对小精度误差具有鲁棒性的解决方案，我建议切换到 linspace。

[m,n] = meshgrid(linspace(rect(1),rect(3),round(rect(3)-rect(1)+1),linspace(rect(2),rect(4),round(rect(4)-rect(2)+1))

作为替代方案，您可以在第二个值中添加一些小数字。这将确保差异始终比整数距离略大，而不是略小。

[m,n] = meshgrid(rect(1):(rect(3)+0.1),rect(2):(rect(4)+0.1));