krp.m

来自「基于卷积信号的MIMO系统盲信号估计」· M 代码 · 共 15 行

function [AkrpB] = krp(A,B);
%%%% calculate the khatri-rao product of A and B

[I F] = size(A);
[J F1] = size(B);

if (F1 ~= F)
 disp('krp.m: column dimensions do not match!!! - exiting matlab');
 exit;
end

AkrpB = [];
for f=1:F,
 AkrpB = [AkrpB kron(A(:,f),B(:,f))];
end