Khatri-Rao积（KR积）
输入：A,B两个矩阵，要求两矩阵列数一致
输出：维度为（A行数*B行数）*列数的矩阵。
计算过程：A,B矩阵相同列做kron积运算，并逐列摆放组成结果矩阵。
举例：
A=[1 2;3 4],B=[5 6;7 8]
其中AB均为2*2的矩阵，得到（2*2）*2=4*2的矩阵。
第一列：用A的第一列和B的第一列做kron积运算
[1*5=5
1*7=7
3*5=15
3*7=21]
第二列：用A的第二列和B的第二列做kron积运算
[2*6=12
2*8=16
4*6=24
4*8=32]
最终结果:
[5 12
7 16
15 24
21 32]
Matlab代码：
for k=1:K %K为输入矩阵的列数
C(:,k)=kron(A(:,k),B(:,k));
end