我有一个选择矩阵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产品？