八度：比較兩個向量

Question

我有這樣的載體：

t = 1: 10 % t = 1 2 3 ..10 

A= [3 4 5] % a column vector 

如果鍵入：

(3 == t) 

我得到的結果是：

0 0 1 0 0 0 0 0 0 0 % it means: 1 at location equals, and 0 at others 

我想爲矢量做到這一點a，這意味着它將採用向量A中的每個元素並比較並返回另一個向量。所以在這種情況下，結果將是一個3×10的矩陣。

但是這行會導致錯誤：A==t。

當然，我可以通過使用for循環來做到這一點，但我想要矢量化這個操作。

Answer 1

你正在尋找的是功能ismember

octave> t = 1:10 
t = 
    1 2 3 4 5 6 7 8 9 10 

octave> A = ismember (t, [2 3 4]) 
A = 
    0 1 1 1 0 0 0 0 0 0 

Answer 2

使用倍頻的廣播能力，它存在於3.6.3（不知道它被引入時），你可以簡單地這樣說：

A'==t 

如果你想讓它產生相同的結果carandraug的ismember命令，你只需添加一個「任何」，像這樣：

any(A'==t) 

對於較小的向量，此方法比ismember方法快得多。

octave:209> tic; for i=1:10000 B=ismember(t,A); end; toc; 
Elapsed time is 1.5 seconds. 
octave:211> tic; for i=1:10000 B=any(A'==t); end; toc; 
Elapsed time is 0.2 seconds. 

注意：如果您的倍頻的版本不支持廣播，或者你想保持它與舊版本兼容，A「== t可通過bsxfun更換（@當量，A」，T） 。