Computing the determinant of a block matrix

4 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


Find more on Sparse Matrices in Help Center and File Exchange




Community Treasure Hunt

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

Start Hunting!