使這個matlab片段運行沒有循環

Question

我想要一個代碼更高效的時間下面的代碼。最好沒有迴路。

參數：

t % time values vector 
t_index = c % one of the possible indices ranging from 1:length(t). 
A % a MXN array where M = length(t) 
B % a 1XN array 

代碼：

m = 1; 
for k = t_index:length(t) 
     A(k,1:(end-m+1)) = A(k,1:(end-m+1)) + B(m:end); 
     m = m + 1; 
end 

非常感謝。

Answer 1

我從B構建了一個大小爲NxM的矩陣（稱之爲B2），根據條件在正確的位置填充零點，然後根據條件填充三角形，然後您需要做的就是A+B2。

是這樣的：

N=size(A,2); 
B2=zeros(size(A)); 
k=c:length(t); 
B2(k(1):k(N),:)=hankel(B) 
ans=A+B2; 

注意，事實上，它是「矢量」並不意味着它是更快的這些日子。 Matlab的JIT使循環具有可比性，有時甚至比內置矢量化選項更快。