※ 引述《CBET ( )》之銘言：
: - - - -
: 主要是使用 toeplitz function，然後再計算 index
t=[A B C D E F G H];
k=8;
[m n]=size(A);
ind1=toeplitz([1 k:-1:2],1:k);
ind1=(ind1-1)*m*n;
ind2=reshape(1:m*n,[m n]);
% modified kron
[ia,ib]=meshgrid(1:k,1:m);
[ja,jb]=meshgrid(1:k,1:n);
ind=ind1(ia,ja)+ind2(ib,jb);
T=t(ind);
clear ind ind1 ind2 ia ib ja jb
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.54.158