Computing determinant of matrix A using LU factorization

27 ビュー (過去 30 日間)
Emily Gallagher
Emily Gallagher 2019 年 10 月 8 日
コメント済み: the cyclist 2019 年 10 月 8 日
I must compute the determinant of an nxn matrix A using LU factorization. I have already calculated L and U in a different 'lufact' function. For some reason though, my funcion for computing the determinant of A, which is just equal to the product of the diagonal entries of matrix U is not working.
seed = 7;
rng(seed)
A = randn(n,n);
[L,U] = lufact(A);
function f = determ(A)
%Computing determinant of A
f = prod(diag(U));
end
  1 件のコメント
the cyclist
the cyclist 2019 年 10 月 8 日
One weirdness here. If your function definition is really
function f = determ(A)
%Computing determinant of A
f = prod(diag(U));
end
then that is going to throw an error, because you don't pass in U. So, maybe you are not really calculating with the inputs you expect?

サインインしてコメントする。

回答 (1 件)

the cyclist
the cyclist 2019 年 10 月 8 日
編集済み: the cyclist 2019 年 10 月 8 日
If I use the lu function from MATLAB (instead of your lufact function), I get the expected result. So, I'd double-check your LU factorization code.
seed = 7;
n = 5;
rng(seed)
A = randn(n,n);
[L,U] = lu(A);
f = prod(diag(U))
d = det(A)
f =
-8.0745
d =
-8.0745

カテゴリ

Help Center および File ExchangeFunction Creation についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by