我有一个选择矩阵W
。我需要这个矩阵转置的kronecker:
% Illustrative setup
n = 100;
ind = rand(1, n) > 0.5;
% Current solution
W = sparse(eye(n));
W = W((ind==1),:);
kronWW = kron(W', W');
然而,这似乎是非常低效的,因为它做了很多1 * 1
操作。由于我知道这基本上是一个花哨的repmat
,我想看看运行时是否因repmat
而改善了 - 但我到目前为止还没有想到办法如果没有显式循环W
的元素(我已经解除了可能的速度,可能会更慢),就可以做到这一点。
有没有办法根据另一个矩阵的逻辑输入有条件地复制矩阵,这个矩阵会在这里复制kronecker产品?