如何在Matlab中对这个for循环进行矢量化?

时间:2015-04-18 14:49:25

标签: matlab for-loop matrix vector vectorization

这可能是一个简单的问题。我有一个向量和一个矩阵,并希望根据一些操作创建一个新的矩阵。我使用for循环构建了新矩阵,我想知道如何用Vector运算符编写它可能更快。

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

我想要新的Delta矩阵[n x n]如下:

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)]

1 个答案:

答案 0 :(得分:3)

您可以使用高效的bsxfun -

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