8

我有一个 Nx1 值向量。我想做的是创建一个 NxN 矩阵，其中每个值代表第 i 个和第 j 个值之间的差异 - 有点像一个大的相关矩阵。我已经用循环完成了这个，但我正在寻找一种更优雅的方法来使用 MATLAB 的矢量化功能，因为这个矢量可能会变得非常大。

4

3 回答 3

12

关于什么

    diff__ = bsxfun(@minus,repmat(A,N,1),A');

这绝对可以通过做来改善

    diff__ = bsxfun(@minus,A,A');

?

一点性能检查：

   N = 1000;
   v = rand(N,1);

   tic
   diff__ = bsxfun(@minus,repmat(v,N,1),v');
   toc

   tic
   diff__ = bsxfun(@minus,v,v');
   toc

结果

  Elapsed time is 105.343344 seconds.
  Elapsed time is 1.124946 seconds.

（蒂姆的数据检查：

差异__ =

 0     2     6     4
-2     0     4     2
-6    -4     0    -2
-4    -2     2     0

）。

于 2012-11-16T18:55:17.597 回答
7

meshgrid可以生成适合此目的的矩阵。获得差分矩阵

meshgrid(v) - meshgrid(v)'

例子：

>> v = [1 3 7 5]

v =

     1     3     7     5

>> meshgrid(v)

ans =

     1     3     7     5
     1     3     7     5
     1     3     7     5
     1     3     7     5

>> meshgrid(v) - meshgrid(v)'

ans =

     0     2     6     4
    -2     0     4     2
    -6    -4     0    -2
    -4    -2     2     0

>>
于 2012-11-16T18:53:04.070 回答
0

已经给出了很好的答案。但是为了凑热闹，这里有另一种方式（使用 Tim 数据）

v=[1 3 7 5];
 cell2mat(arrayfun(@(i) (v(i)-v)',1:size(v,2), 'UniformOutput',false))

答案=

 0     2     6     4
-2     0     4     2
-6    -4     0    -2
-4    -2     2     0
于 2012-11-16T19:08:32.323 回答