Есть ли более быстрый способ для этого сценария MATLAB?

Question

Оригинальный сценарий

Y = [1 2 3 3 2 1 1 2 3]'; 
    n = length(Y); 
    Ym = zeros(n, n); 
    for i=1:n 
     index = find(Y==Y(i)); 
     Ym(i, index') = 1; 
    end 

пор Ym является

Ym = 

1  0  0  0  0  1  1  0  0 
0  1  0  0  1  0  0  1  0 
0  0  1  1  0  0  0  0  1 
0  0  1  1  0  0  0  0  1 
0  1  0  0  1  0  0  1  0 
1  0  0  0  0  1  1  0  0 
1  0  0  0  0  1  1  0  0 
0  1  0  0  1  0  0  1  0 
0  0  1  1  0  0  0  0  1 

Answer 1

Да! С bsxfun -

Ym = bsxfun(@eq,Y,Y.') 

Answer 2

Я считаю, что логическая индексация работает быстрее, чем bsxfun на моем компьютере. Здесь времена выборки для различных методов:

tic; 
for j=1:10000 
    Y = [1 2 3 3 2 1 1 2 3]'; 
    n = length(Y); 
    Ym = zeros(n, n); 
    for i=1:n 
     index = find(Y==Y(i)); 
     Ym(i, index') = 1; 
    end 
end 
disp('Method 1:'); 
toc; 

tic; 
for j=1:10000 
    Y = [1 2 3 3 2 1 1 2 3]'; 
    n = length(Y); 
    Ym = zeros(n, n); 
    for i=1:n 
     Ym(i, Y==Y(i)') = 1; 
    end 
end 
disp('Method 2:'); 
toc; 
tic; 
for j=1:10000 
    Y = [1 2 3 3 2 1 1 2 3]'; 
    n = length(Y); 
    Ym = zeros(n, n); 
    a=repmat(Y,1,n); 
    b=repmat(Y',n,1); 
    Ym(a==b)=1; 
end 
disp('Method 3:'); 
toc; 

tic; 
for j=1:10000 
    Y = [1 2 3 3 2 1 1 2 3]'; 
    Ym = bsxfun(@eq,Y,Y.'); 
end 
disp('Method 4'); 
toc 

ВЫВОД:

Method 1: 
Elapsed time is 0.111412 seconds. 
Method 2: 
Elapsed time is 0.069617 seconds. 
Method 3: 
Elapsed time is 0.246780 seconds. 
Method 4 
Elapsed time is 0.103120 seconds.