如何比较两个数组的所有元素？

Question

我有两个大约 1000 行和 1000 列的大数组。我需要比较这些数组的每个元素，如果相应的元素相等，则将 1 存储在另一个数组中。

我可以用 for 循环来做到这一点，但这需要很长时间。我怎样才能更快地做到这一点？

Answer 1

给出的答案都是正确的。我只是想详细说明gnovice的评论关于浮点测试。

比较浮点数是否相等时，需要使用容差值。常用的公差比较有两种：绝对公差和相对公差。 （来源）

绝对容差a 和 b 的比较如下所示：

|a-b| < tol

相对容差比较如下所示：

|a-b| < tol*max(|a|,|b|) + tol_floor

您可以将上述两个函数实现为匿名函数：

%# absolute tolerance equality
isequalAbs = @(x,y,tol) ( abs(x-y) <= tol );

%# relative tolerance equality
isequalRel = @(x,y,tol) ( abs(x-y) <= ( tol*max(abs(x),abs(y)) + eps) );

然后您可以将它们用作：

%# let x and y be scalars/vectors/matrices of same size
x == y
isequalAbs(x, y, 1e-6)
isequalRel(x, y, 1e-6)

Answer 2

如果您的两个矩阵 A 和 B 大小相同，那么您可以这样做：

index = A == B;

并且 index 将是 逻辑数组，其中 1 处处都是 A 的元素和 B 相等且为零。

警告...

如果A和B包含整数，上面的应该没问题。但是，如果它们包含浮点值，您可能会得到不需要的结果。对于完全相等的元素，上面的代码只会有值1。即使是最小的差异也会导致元素被视为不相等。

您可以查看此问题的答案 有关处理“浮点运算的危险”的更多信息。一种解决方案是检查数组元素是否在彼此给定的容差范围内，如下所示：

tolerance = 0.0001;
index = abs(A-B) <= tolerance;

上面将为您提供一个逻辑数组 index ，其中 A 和 B 彼此相差在 0.0001 以内，否则为零。

Answer 3

只需使用普通的 == 运算符即可：

>> [1 2; 3 4] == [1 5; 6 4]      

ans =

     1     0
     0     1