从多维矩阵中提取块

Question

考虑多维矩阵A，其中size(A) 具有相同的偶数元素N。应该如何找到矩阵B 和size(B)=size(A)/2，这样：

B(1,1,...,1)=A(1,1,...,1),
B(1,1,...,2)=A(1,1,...,2),
...
B(N/2,N/2,...,N/2)=A(N/2,N/2,...,N/2).

Answer 1

我一般不喜欢arrayfun（或循环函数），但如果维度数不是以千计，那么这应该没问题：

Nv = size(A)/2;
S = arrayfun(@(x){1:x},Nv);
B = A(S{:});

也应该适用于不同尺寸的尺寸。只需决定如何处理mod(size(A),2)~=0 所在的维度。