Octave 的默認 LU 分解函數錯誤

Question

這個論壇的新手。

我正在嘗試以完整的旋轉方式運行八度音程的 LU 分解函數：

[L, U, p, q] = lu(A)

對於矩陣 AI，我不斷收到此錯誤：

“返回列表中未定義元素編號 4”

元素 4 是列排列矩陣 問：這是怎么回事？ 為什么不顯示？ 提前致謝

Answer 1

如果矩陣A已滿，則lu函數不會在 Octave 中執行列交換（強調我的）：

當使用兩個或三個輸出參數和一個備用 [原文如此] 輸入矩陣調用時， lu不會嘗試執行保持稀疏性的列排列。 使用第四個輸出參數調用，返回保留稀疏性的列轉換Q ，使得P * A * Q = L * U 。

因此，如果為稀疏矩陣提供第四個輸出參數，則僅對稀疏矩陣執行完全旋轉以最大化稀疏性。 上面的引用使用“ A ”，但根據鏈接的 Octave 文檔部分頂部提供的函數簽名，我相信他們的意思是寫“ S ”：“ [L, U, P, Q] = lu (S) ” .

默認情況下，似乎沒有完整矩陣的完整旋轉選項。

我會注意到 MATLAB 對其lu的 第四個輸出具有相同的行為：

列置換... . 使用此輸出來減少稀疏矩陣因子中的填充（非零值的數量）。