Fast Kronecker matrix multiplication: 任何大小矩阵的快速 Kronecker matrix multiplication-matlab开发

快速 Kronecker 矩阵乘法，适用于全矩阵和稀疏矩阵任何大小。 从不计算实际的 Kronecker 矩阵并省略乘以单位矩阵。 y = kronm(Q,x) 计算y = (Q{k} kron ... Q{2} kron Q{1})*x 如果 Q 仅包含两个矩阵且 x 是向量，则代码使用身份( Q{2} 克朗 Q{1} )*vec(X) = vec(Q{1}*X*Q{2}'), 其中 vec(X)=x。 如果 Q 包含两个以上的矩阵和/或如果 x 有更多除了一列之外，该算法使用此身份的广义形式。 该算法的思想是将 x 视为一个多维数组并分别为每个维度 i 应用线性映射 Q{i}。 致谢： 此代码遵循与 Paul G. Constantine & 的“kronmult”相同的想法David F. Gleich（斯坦福大学，2009 年）。 但是，我避免循环并允许非平方输入 Q{i}。 我还