如何在Matlab中选择子矩阵(不是任何特定的模式)

时间:2012-10-26 16:55:42

标签: matlab matrix matrix-indexing

如何在Matlab中选择子矩阵(不是任何模式)?例如,对于大小为10乘10的矩阵,如何选择由第1和第9行与第4和第6列的交点组成的子矩阵?

感谢您提供任何有用的答案!

3 个答案:

答案 0 :(得分:40)

TLDR:简答

至于你的问题,假设你有一个10×10的任意矩阵A。提取所需子矩阵的最简单方法是使用索引向量:

B = A([1 2 9], [4 6]);


在MATLAB中建立索引

官方文档中有一个interesting article,可以全面解释MATLAB中的索引。 基本上,有几种方法可以提取值的子集,我将为您总结它们:

1。索引向量

索引向量指示要提取的元素的索引。它们可以包含一个或多个索引,如下所示:

A = [10 20 30 40 50 60 70 80 90]

%# Extracts the third and the ninth element
B = A([3 9])  %# B = [30 90]

可以分别为每个维度指定索引向量,例如:

A = [10 20 30; 40 50 60; 70 80 90];

%# Extract the first and third rows, and the first and second columns
B = A([1 3], [1 2])  %# B = [10 30; 40 60]

还有两个特殊的下标:end和冒号(:):

  • end只是表示该维度中的最后一个索引。
  • 冒号只是“1:end”的简写符号。

例如,您可以编写A([1 2 3], [2 3]),而不是撰写A(:, 2:end)。这对于大型矩阵特别有用。

2。线性索引

线性索引通过将列连接到一个列向量并分别为元素分配索引,将任何矩阵视为列向量。例如,我们有:

A = [10 20 30; 40 50 60; 70 80 90];

我们要计算b = A(2)。等效列向量是:

A = [10;
     40;
     70;
     20;
     50;
     80;
     30;
     60;
     90]

因此b等于40。

当然也允许使用特殊冒号和end下标。因此,A(:)会将任何矩阵A转换为列向量。

使用矩阵下标进行线性索引: 也可以使用另一个矩阵进行线性索引。下标矩阵简单地转换为列向量,并用于线性索引。然而,得到的矩阵总是与下标矩阵具有相同的尺寸 例如,如果I = [1 3; 1 2],则A(I)与撰写reshape(A(I(:)), size(I))相同。

从矩阵下标转换为线性索引,反之亦然: 为此,您分别有sub2indind2sub。例如,如果要将矩阵[1, 3]中的下标A(对应于元素30)转换为线性索引,则可以编写sub2ind(size(A), 1, 3)(在这种情况下,结果应为7 ,当然)。

3。逻辑索引

在逻辑索引中,下标是二进制的,其中逻辑1表示选择了相应的元素,而0表示不是。下标向量必须与原始矩阵具有相同的维度,或者具有相同数量的元素。例如,如果我们有:

A = [10 20 30; 40 50 60; 70 80 90];

我们希望使用逻辑索引提取A([1 3], [1 2]),我们可以这样做:

Ir = logical([1 1 0]);
Ic = logical([1 0 1]);
B = A(Ir, Ic)

或者这个:

I = logical([1 0 1; 1 0 1; 0 0 0]);
B = A(I)

或者这个:

I = logical([1 1 0 0 0 0 1 1 0]);
B = A(I)

请注意,后两种情况是一维向量,如有必要,应重新形成矩阵(例如,使用reshape)。

答案 1 :(得分:13)

让我用一个例子来解释:

让我们定义一个6x6矩阵

A = magic(6)

A = 
35     1     6    26    19    24
 3    32     7    21    23    25
31     9     2    22    27    20
 8    28    33    17    10    15
30     5    34    12    14    16
 4    36    29    13    18    11

从这个矩阵中,您需要第1,2和5行以及第4和第6列中的元素

B = A([1 2 5],[4 6])

B = 

26    24
21    25
12    16

希望这有帮助。

答案 2 :(得分:0)

function f = sub(A,i,j)
[m,n] = size(A);

row = 1:m;
col = 1:n;

x = row;
x(i) = [];

y=col;
y(j) = [];

f= A(x,y);

返回矩阵A,删除了第i th 行和第j th 列。