2

我想从 MATLAB 中的 LUP 分解中计算矩阵的行列式。行列式可以从公式中找到:

在此处输入图像描述

P 是一个置换矩阵,S 是将 P 转换为单位矩阵所需的行交换数。如何在 MATLAB 的上述公式中找到 S?它是否有任何预定义的功能等?

4

1 回答 1

3

如果您解释P为邻接矩阵,并且向量cycles包含由 描述的图中所有循环的长度P,则S=sum(cycles) - length(cycles).

现在剩下的就是找到所有循环的长度,文件交换上有几个函数,比如这个

顺便说一句:[L, U, P] = lu(A)det(A) = det(inv(P))*det(L)*det(U)

于 2014-01-10T14:20:14.110 回答