使用重复的 Kronecker 乘积矩阵实现快速高效的矩阵向量乘法计算 - Matlab开发

计算矩阵向量乘积 y = (Q1 kron Q2 kron ... kron Qm) * x 可以在不形成 Kronecker 产品的大矩阵的情况下完成。 这个 m 文件为 Fernandes 等人的第 394 页中的这个任务实现了一个算法。 1998, JACM 45(3): 381--414 (doi:10.1145/278298.278303)。 该实现也适用于 X 是矩阵的情况。 不要被 for 循环吓到，这段代码与 Matlab JIT 编译器配合得很好，适用于超过 5000 万个条目的向量。