如何有效地从（非稀疏）矩阵中删除零？

Question

我有一个矩阵：

x = [0 0 0 1 1 0 5 0 7 0];

我需要删除所有零，如下所示：

x = [1 1 5 7];

我使用的矩阵很大（1x15000），我需要多次执行此操作（5000+），所以效率是关键！

Answer 1

一种方法：

x(x == 0) = [];

关于时间的说明：

正如木片，与 狐YMG。 Loren 在她的一篇 MathWorks 博客文章。 由于您提到必须执行此操作数千次，因此您可能会注意到差异，在这种情况下，我会首先尝试 x = x(x~=0); 。

警告：如果您使用非整数，请小心。 例如，如果您有一个非常小的数字，您希望将其视为足够接近零以便将其删除，则上述代码将不会删除它。 仅删除精确的零。 以下内容还将帮助您删除“足够接近”零的数字：

tolerance = 0.0001;  % Choose a threshold for "close enough to zero"
x(abs(x) <= tolerance) = [];

Answer 2

只是为了与众不同：

x=x(x~=0);

或者

x=x(abs(x)>threshold);

这也有处理复数的好处

Answer 3

这是三种常见的解决方案。 它有助于看到差异。

x = round(rand(1,15000));

y = x;

tic,y(y==0) = [];toc

Elapsed time is 0.004398 seconds.

y = x;

tic,y = y(y~=0);toc

Elapsed time is 0.001759 seconds.

y = x;

tic,y = y(find(y));toc

Elapsed time is 0.003579 seconds.

正如您所看到的，最便宜的方法是直接逻辑索引，选择要保留的元素。 查找的成本更高，因为 matlab 查找这些元素，返回它们的列表，然后索引到向量中。

Answer 4

这是另一种方法

y = x(find(x))

我将让您来计算您尝试的各种方法的相对效率 - 请写下来并让我们都知道。

Answer 5

虽然我的计时结果并不能确定它是否明显更快，但这似乎是最快和最简单的方法：

y = nonzeros(y) 

Answer 6

x = [0 0 0 1 1 0 5 0 7 0]
y = [0 2 0 1 1 2 5 2 7 0]

然后 x2 和 y2 可以通过以下方式获得：

x2=x(~(x==0 & y==0))
y2=y(~(x==0 & y==0))

x2 = [0     1     1     0     5     0     7]
y2 = [2     1     1     2     5     2     7]

希望这有帮助！