MATLAB: Sparse knonecker product of large and several matrices

kroneckersparse

Hi,
I'm trying to create a sparse kronecker product without actually generating the kron product and then generating the sparse matrix. Is there a way to do this? For example
E=eye(100);
I=ones(20);
K=sparse(kron(E,I));
Right now the code first generates the kron product and then keeps on the information of non-zero elements. As I need to do this lots of time because I need to generate and keep in memory different K matrices (K1, K2, … K200, for E1,…E200, and I1,…I200), the generation of the kron product uses all my memory. Is there anyway I could go from the E's and I's to the K's without generating the full matrices?
Thanks,

Best Answer

E=speye(100);
I=ones(20);
K=kron(E,I);