Creazione di una matrice diagonale a blocchi sparsi in Matlab

Question

Supponiamo che B sia una matrice di celle di matrici sparse in Matlab e voglio formare una matrice diagonale a blocchi sparsi M i cui blocchi diagonali sono le matrici archiviate in B. Qual è il più semplice/più efficiente modo di farlo?

Answer 1

Uso blkdiag su comma-separated list generato dalla matrice di celle:

result = blkdiag(B{:}); 

Ad esempio, con

B = {sparse([1 0 0; 2 2 0; 3 3 3]), 4*speye(2)}; 

questo produce

>> result 
result = 
    (1,1)  1 
    (2,1)  2 
    (3,1)  3 
    (2,2)  2 
    (3,2)  3 
    (3,3)  3 
    (4,4)  4 
    (5,5)  4 
>> full(result) 
ans = 
    1  0  0  0  0 
    2  2  0  0  0 
    3  3  3  0  0 
    0  0  0  4  0 
    0  0  0  0  4