數組每列中的唯一元素（Matlab）

Question

我可能會用Matlab生鏽一點，也許這個問題的答案比我現在想象的更微不足道。不過，我已經在網上搜索有效的解決方案，我還沒有找到任何，所以我會在這裏嘗試。

我有一個大的矩陣，像ÿ這裏：

N = 1e6; k = 20; n = 100; 
y = ceil(n * rand(k, N)); 

因爲我想算獨特的元素數量的陣列的每列。該循環非常緩慢：

tic 
r = zeros(N, 1); 
for ii = 1:N 
    r(ii) = numel(unique(y(:, ii))); 
end 
toc 

尋找矢量化，更快的版本。

---

大衛的答案似乎是正確的和快。謝謝！

N = 1e6; k = 20; n = 100; 
y = ceil(n * rand(k, N)); 

tic 
r1 = zeros(N, 1); 
for ii = 1:N 
    r1(ii) = numel(unique(y(:, ii))); 
end 
toc 

tic 
r2 = sum(diff(sort(y)) ~= 0) + 1; 
toc 

all(r1' == r2) 

Answer 1

試試這個

sum(diff(sort(y))~=0)+1 

它使用做vectorise功能。對於你的情況，性能似乎比你的for循環更好，但是我會想象更大的問題，由於內存限制，我提出的解決方案會變得更糟。

隨着N=1e5，你的方法在我的電腦上花費約7.5秒，我的建議花費約0.05秒。與N=1e6，時間約爲75s和〜0.75s。