Comment sélectionner une sous-matrice (pas dans un modèle particulier) dans Matlab

TLDR: Réponse Courte

Quant à votre question, supposons que vous avez un arbitraire de 10 par 10 de la matrice A. La façon la plus simple d'extraire les sous-matrice serait avec un indice de vecteur:

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

Indexation dans MATLAB

Il y a un article intéressant dans la documentation officielle que de manière exhaustive l'indexation dans MATLAB.
Fondamentalement, il existe plusieurs façons d'extraire un sous-ensemble de valeurs, je vais résumer pour vous:

1. L'Indexation Des Vecteurs

Indexation des vecteurs indiquent les indices de l'élément à extraire. Ils peuvent contenir un index unique ou à plusieurs, comme suit:

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

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

Indexation des vecteurs peut être spécifié pour chaque dimension séparément, par exemple:

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]

Il y a aussi deux indices: end et le côlon (:):

• end indique simplement le dernier indice de cette dimension.
• Le côlon est juste un court de main de notation pour les "1:fin".

Par exemple, au lieu d'écrire A([1 2 3], [2 3])vous pouvez écrire A(:, 2:end). Ceci est particulièrement utile pour les grandes matrices.

2. Linéaire Indexation

Linéaire indexation traite toute matrice comme s'il s'agissait d'un vecteur colonne par la concaténation des colonnes dans un vecteur colonne et l'attribution des indices pour les éléments respectivement. Par exemple, nous avons:

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

et nous voulons calculer b = A(2). L'équivalent vecteur colonne est:

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

et donc b est égal à 40.

La spéciale du côlon et end indices sont aussi autorisés, bien sûr. Pour cette raison, A(:) convertit une matrice A dans un vecteur colonne.

Linéaire de l'indexation avec la matrice des indices:
Il est également possible d'utiliser une autre matrice linéaire de l'indexation. L'indice de la matrice est tout simplement converti en un vecteur colonne, et utilisé pour l'indexation linéaire. La matrice obtenue est cependant toujours de la même dimension que l'indice de la matrice.

Par exemple, si I = [1 3; 1 2], puis A(I) est la même que l'écriture de reshape(A(I(:)), size(I)).

De la conversion à partir de la matrice des indices linéaire des indices et vice versa:
Pour cela, vous avez sub2ind et ind2sub, respectivement. Par exemple, si vous voulez convertir les indices [1, 3] dans la matrice A (correspondant à l'élément 30) dans un linéaire d'index, vous pouvez écrire sub2ind(size(A), 1, 3) (le résultat est dans ce cas devrait être de 7, bien sûr).

3. Logique D'Indexation

Dans la logique d'indexation les indices sont binaires, où une logique 1 indique que l'élément correspondant est sélectionné, et 0 signifie qu'il n'est pas. L'indice du vecteur doit être de mêmes dimensions que l'original de la matrice ou d'un vecteur avec le même nombre d'éléments. Par exemple, si nous avons:

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

et nous voulons extraire A([1 3], [1 2]) à l'aide de logiques d'indexation, nous pouvons faire ceci:

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

ou ceci:

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

ou ceci:

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

Noter que dans les deux derniers cas est une dimension du vecteur, et devraient être réorganisées de retour dans une matrice si nécessaire (par exemple, à l'aide de carré).

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);

Renvoie la matrice A, jeth ligne et jth colonne supprimée.