Guest User

Untitled

a guest
Jan 21st, 2019
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.41 KB | None | 0 0
  1. T = 5;
  2. dim = 3;
  3.  
  4. dKmat = Array[dKm, {T, T}];
  5. M1mat = Array[M1m, {dim, dim}];
  6. Kron3 = KroneckerProduct[M1mat, dKmat];
  7.  
  8. ArrayReshape[Kron3, {Times @@ Dimensions[Kron3], 1}] ===
  9. KroneckerProduct[
  10. Flatten[Table[
  11. KroneckerProduct[IdentityMatrix[T, SparseArray],
  12. Transpose[{M1mat[[;; , i]]}]], {i, dim}], 1],
  13. IdentityMatrix[T, SparseArray]].ArrayReshape[
  14. dKmat, {Times @@ Dimensions[dKmat], 1}]
Add Comment
Please, Sign In to add comment