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).')。