如何获得 unique() 的二进制输出？

Question

[uniq_elements,uniq_indices]=unique(x)

显示独特的元素，例如矢量x。

什么是简单（即单行）和高效（即快速执行时间）的方法来简单且仅指示“此值是否唯一”和“如果不是”输出一系列重复项中的第一个？有没有办法使用逻辑运算符来解决这个问题？简单地让一个函数输出“这个值是否唯一”可能是一个开始。

示例代码：

x=[5 5 2 2 3 4 4 1 1 1]; %example values
[~,uniq_indices]=unique(x); %returns [1,3,5,6,8]
result=false(size(x));
result(uniq_indices)=1
%console ouput:
result =

  1×10 logical array

   1   0   1   0   1   1   0   1   0   0

基本上我的假设是——因为只需要逻辑输出——上述代码可以更快/更容易地实现。

Answer 1

方法 1

这是内存密集型的，因为它创建了一个中间n×n 逻辑矩阵，其中n 是x 的大小：

result = sum(triu(bsxfun(@eq, x, x.')))==1;

从 R2016b 开始，这可以更简洁地表达为

result = sum(triu(x==x.'))==1;

工作原理

bsxfun(@eq, x, x.')（或x==x.'）创建一个相等比较的方阵。

triu 只保留上三角部分，因此每个元素只与之前的元素或自身进行比较。

sum 给出每列的总和。如果总和是1，这意味着该元素不等于任何先前的元素，仅等于它自己。

方法2

这稍微更节省内存。它创建了一个中间m×n 逻辑矩阵，其中m 是唯一元素的数量，n 是元素的总数：

xu = unique(x);
result = any(diff([false(numel(xu),1) bsxfun(@eq, x, xu.')],[],2)==1);

从R2016b开始，第二行可以替换为

result = any(diff([false(numel(xu),1) x==xu.'],[],2)==1);

工作原理

这将创建一个矩阵，例如 A，将 x（A 的行索引）的每个唯一元素与 x 的每个元素（A 的列索引）进行比较。让B 表示将false 列附加到A 的结果。

第一次出现在x 中的元素对应B 对应行中的[0 1] 子序列。为了检测到这一点，diff 沿每一行应用，并将结果与​​1（这是[0 1] 中元素之间的增量）进行比较。