MATLAB: creare una matrice diagonale a blocchi con lo stesso blocco ripetuto

Question

Ho una matrice K di dimensioni n x n. Voglio creare un nuovo blocco matrice diagonale M di dimensioni N x N, tale da contenere d blocchi di matrice K come sua diagonale.

avrei usato direttamente M = blkdiag (K, K, K) ecc avuto d stato inferiore. Sfortunatamente, d è molto grande e non voglio scrivere manualmente la formula con d esattamente gli stessi argomenti per la funzione blkdiag().

Esiste un modo più breve e intelligente per farlo?

Answer 1

per questo è possibile utilizzare kron.

M = kron(X,Y) 

restituisce il prodotto tensoriale Kronecker di X e Y. Il risultato è una vasta gamma formata adottando tutte le possibili prodotti tra gli elementi di X e quelli di Y. Se x è m-by-n e Y è p-by-q, quindi kron (X, Y) è m * p-by-n * q. Quindi nel tuo caso qualcosa come questo farà:

M = kron(eye(L),K) 

con L il # di blocchi.

Answer 2

Quanto segue dovrebbe funzionare:

d = 5; K = occhio (3); T = cell (1, d);

per j = 1: d T {j} = K; fine

M = blkdiag (T {:})

Answer 3

s = 'A,'; 
s = repmat(s,[1,n2]); 
s = ['B=blkdiag(', s(1:end-1),');']; 
eval(s); 

Può essere più veloce rispetto all'utilizzo di Kron-eye.

Answer 4

tmp = repmat({K},d,1); 
M = blkdiag(tmp{:}); 

Non si dovrebbe mai usare eval, o andare dentro per loop inutilmente. Kron è un modo molto elegante. Volevo solo condividere questo come funziona anche.

Answer 5

Un ciclo "for" potrebbe essere d'aiuto. Ad esempio:

M = k; 
for i=1:N/n - 1 
    M=blkdiag(M,k); 
end