我有以下问题:
我有一个给定的矩阵,让我们说4x4
。
如何获得以下组合的索引:
例如,如果我有矩阵:
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
此外,是否可以组合这些操作,例如我得到了“奇数行”和“奇数列”的索引。和'行甚至列甚至'?
谢谢!
答案 0 :(得分:11)
使用索引很容易:
B = A(1:2:end, 1:2:end);
B = A(1:2:end, 2:2:end);
B = A(2:2:end, 1:2:end);
B = A(2:2:end, 2:2:end);
以上假设您自己想要实际的矩阵值。这有点令人困惑,因为您的矩阵元素与线性索引值本身相同。如果要确定要访问矩阵的实际列主索引,可以从1 to N
生成向量,其中N
是矩阵中元素的总数,然后将此矩阵重新整形为所需大小你要的那个。之后,使用上面相同的逻辑来获得实际的线性指数:
N = numel(A);
B = reshape(1:N, size(A,1), size(A,2));
ind = B(1:2:end, 1:2:end); %// For odd rows, odd columns
%// Repeat for the other ones...
现在,根据您的评论,您想要创建一个新矩阵,该矩阵将仅存储这些提取的矩阵值,同时使所有其他元素为零。如果要执行此操作,只需预先分配一个零矩阵,然后复制这些值以使用计算的索引提取到新矩阵中。换句话说:
N = numel(A);
B = reshape(1:N, size(A,1), size(A,2));
ind = B(1:2:end, 1:2:end); %// For odd rows, odd columns - Change to suit your tastes
out = zeros(size(A));
out(ind(:)) = A(ind(:));
如果你想组合像奇数行 - 奇数列甚至行 - 偶数列这样的索引,只需计算两组索引,将它们连接成一个向量,并像以前一样执行相同的语法。因此:
N = numel(A);
B = reshape(1:N, size(A,1), size(A,2));
ind = B(1:2:end, 1:2:end); %// For odd rows, odd columns
ind2 = B(2:2:end, 2:2:end); %// For even rows, even columns
ind = [ind(:); ind2(:)];
out = zeros(size(A));
out(ind) = A(ind);
答案 1 :(得分:2)
<强>代码强>
N = size(A,1); %// Get size of input matrix A
case1_ind = bsxfun(@plus,[1:2:N]',(0:N/2-1)*2*N)
case2_ind = case1_ind + N
case3_ind = case1_ind + 1
case4_ind = case3_ind + N
注意:这些输出是索引。因此,要获得实际输出,请将它们用作索引。
要合并case 1
和case 4
的索引,只需连接 -
case14comb_ind = [case1_ind ; case4_ind]
修改:
%// To copy onto some other matrix of the same size as A, do this for case 1
new_matrix = zeros(size(A))
new_matrix(case1_ind(:)) = A(case1_ind(:))
对其他情况也重复此操作。