3個特定的嵌套for循環的向量化

Question

我有3循環，我想矢量化它。

我嘗試了很多，但我做不到。

這是簡單的版本。

我無法解釋整個情況，但：

A是100 * 82。

B是203 * 2000。

C是1 * 10。

D是100 * 5 * 2000。

for nset=1:100 
     for nindy=1:5 
      for nimp=1:2000 
    if A(nset,2+(nindy-1).*3)==1 
     if B(33.*6+nindy,nimp)<= C(1+(nindy-1).*2)+... 
      ((C(2+(nindy-1).*2)-C(1+(nindy-1).*2))./11).*A(nset,1+(nindy-1).*3) 

     D(nset,nindy,nimp)=1; 
     end 
    end 

     end 
     end 
    end 

我想更換3與量化的循環......

感謝....

Answer 1

是的，這是可能的！這些線路應該這樣做：

nset=1:100; 
nindy=1:5; 
nimp=1:2000; 
D(bsxfun(@and, (A(nset,2+(nindy-1).*3)==1),bsxfun(@le,shiftdim(B(33.*6+nindy,nimp),-1),bsxfun(@plus,C(1+(nindy-1).*2),bsxfun(@times,((C(2+(nindy-1).*2)-C(1+(nindy-1).*2))./11),A(nset,1+(nindy-1).*3))))))=1; 

通過我的電腦改變了運行時間約1秒到1/20秒。