function - Octave 的默认LU分解函数错误

Question

该论坛的新手。

我正在尝试像这样完全旋转来运行 Octave 的LU分解功能:

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

对于我拥有的矩阵A，我不断收到此错误:

"element number 4 undefined in return list"

元素4是列置换矩阵Q。这是怎么回事？为什么不显示？提前致谢

Answer 1

如果矩阵A已满，则lu函数does not perform column exchanges in Octave(强调我的):

When called with two or three output arguments and a spare[sic] input matrix, lu does not attempt to perform sparsity preserving column permutations. Called with a fourth output argument, the sparsity preserving column transformation Q is returned, such that P * A * Q = L * U.

因此，如果为稀疏矩阵提供了第四个输出参数，则仅对稀疏矩阵执行 full pivoting以使稀疏性最大化。上面的引号使用“A”，但是根据链接的Octave文档部分顶部提供的功能签名，我认为它们的意思是写“S”:“ [L, U, P, Q] = lu (S)”。

默认情况下，似乎没有用于完整矩阵的完整旋转选项。

我会注意到MATLAB的 fourth output of its lu 具有相同的行为:

Column permutation ... . Use this output to reduce the fill-in (number of nonzeros) in the factors of a sparse matrix.