Computing the determinant of a block matrix

5 views (last 30 days)
I am having trouble using a well-known formula for computing the determinant of a block matrix. That is, if
in which A and D are square matrices and exists then
The code snippet below should achieve this but it returns two different results. I'm sure it's a very trivial mistake I've made but I've probably been staring at it too long now to find it.
As a follow up, does anyone know if there is a suitably adjusted version of this formula for the case when it is B and C which are the square matrices?
function [detM,detBlock] = BlockDeterminant()
M = magic(4);
detM = det(M);
A = M(1:3,1:3);
B = M(1:3,4);
C = M(4,1:3);
D = M(4,4);
detBlock = det(A)*det(D-C*inv(A)*B);

Accepted Answer

KSSV on 1 Aug 2022
A = rand(2) ;
B = rand(2) ;
C = rand(2) ;
D = rand(2) ;
M = [A B; C D] ;
ans = -0.0221
ans = -0.0221
In your case, D is a scalar.
Bruno Luong
Bruno Luong on 1 Aug 2022
Good. I don't get why contribution is so helpful, but if you say so...

Sign in to comment.

More Answers (1)

Walter Roberson
Walter Roberson on 1 Aug 2022
Just round-off error. You are calculating two different order of operations and that affects the results .
If you use M = sym(magic(4)) you will get exact zero for each




Community Treasure Hunt

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

Start Hunting!

Translated by