Khatri-Rao积(KR积)

Khatri-Rao积(KR积)_第1张图片

输入: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

你可能感兴趣的:(MATLAB)