在matlab中,矩阵的存储是按列优先,sub2ind函数将矩阵中指定元素的行列下标转换成存储的序号,即线性索引号。下面,我们举例子进行说明。
1 建立一个3*4*2的矩阵
rng(0,'twister'); % Initialize random number generator.
A = rand(3, 4, 2)
A(:,:,1) =
0.8147 0.9134 0.2785 0.9649
0.9058 0.6324 0.5469 0.1576
0.1270 0.0975 0.9575 0.9706
A(:,:,2) =
0.9572 0.1419 0.7922 0.0357
0.4854 0.4218 0.9595 0.8491
0.8003 0.9157 0.6557 0.9340
将元素A(2, 1, 2)的下标转换成线性索引号
linearInd = sub2ind(size(A), 2, 1, 2) linearInd = 14
A(2, 1, 2) A(14) ans = and = 0.4854 0.4854
2 针对这种三维的矩阵,如果缺省第三个下标,默认为1。
linearInd = sub2ind(size(A), 2, 4) ans = 11
linearInd = sub2ind(size(A), 2, 4, 1) ans = 11
3 假如要转换多个元素,我们可以使用一条命令完成。例子是转换5个元素下标为线性索引号。
sub2ind(size(A), 3, 3, 2); sub2ind(size(A), 2, 4, 1); sub2ind(size(A), 3, 1, 2); sub2ind(size(A), 1, 3, 2); sub2ind(size(A), 2, 4, 1);
sub2ind(size(A), [3 2 3 1 2], [3 4 1 3 4], [2 1 2 2 1]) ans = 21 11 15 19 11
%%[3 2 3 1 2], [3 4 1 3 4], [2 1 2 2 1] 分别代表行,列,纵向的下标。%%
[A(3,3,2), A(2,4,1), A(3,1,2), A(1,3,2), A(2,4,1)] ans = 0.6557 0.1576 0.8003 0.7922 0.1576 A([21, 11, 15, 19, 11]) ans = 0.6557 0.1576 0.8003 0.7922 0.1576