如何在Matlab中向量化此for循環？

Question

這也許是一個簡單的問題。 我有一個向量和一個矩陣，並希望基於某種操作來制作一個新的矩陣。 我使用for循環構造了新矩陣，我想知道如何用可能更快的Vector運算符編寫它。

d=[n x 1];
t= [n x n];

我想要的新增量矩陣是[nxn]，如下所示：

for i=1:39
    for j=1:39
  Delta(i,j)=d(i)-d(j)-t(i,j);
    end
end

結果

[d (1) - d (1) - t( 1 ,1),d (1) - d (2) - t( 1 ,2), ...   d(1) - d (39) -  t( 1 ,39)
 d (2) - d (1) - t( 2 ,1),d (2) - d (2) - t( 2 ,2), .... ,d (2) - d (39) - t( 2 ,39)  
.                                                                     
.
.
d (38) - d (1) - t( 38 ,1),d (38) - d (2) - t( 38 ,2), ... , d(38) -d (39)-t(38,39)  
d (39) - d (1) - t( 39 ,1),d (39) - d (2) - t( 39 ,2), ...,  d(39)- d (39)- t(39 ,39)]

Answer 1

您可以使用高效的bsxfun

Delta = bsxfun(@minus,d,d.') - t