Matlab没有优化以下内容吗？

Question

我有一个很长的向量 1xr v，和一个很长的向量 w 1xs，还有一个矩阵 A rxs，它是稀疏的（但尺寸很大）。

我希望 Matlab 会优化以下内容，这样我就不会遇到内存问题：

 A./(v'*w)

但似乎 Matlab 实际上正在尝试生成完整的 v'*w 矩阵，因为我遇到了内存不足的问题。有没有办法克服这个问题？请注意，不需要计算所有v'*w，因为A 的许多值都是0。

编辑：如果可能的话，一种方法是A(find(A))./(v'*w)(find(A));

但是您不能选择矩阵的子集（在这种情况下为v'*w）而不先计算它并将其放入变量中。

Answer 1

• 您可以使用bsxfun。这给出了与A./(v'*w) 相同的结果，但不生成矩阵v.'*w：

  bsxfun(@rdivide, bsxfun(@rdivide, A, v'), w)
  

• 另一种可能性：如果您只想要非零值，请使用：

  [ii jj Anz] = find(A);
  Anz./v(ii)'./w(jj).'
  

  这给出了一个对应于您的A(find(A))./(v'*w)(find(A)) 的列向量，同样不会生成v.'*w。如果您需要稀疏矩阵A./(v'*w)（而不是其非零值的列向量），请使用sparse(ii,jj,Anz./v(ii)'./w(jj).')。