Discussion :

[assma.m]

Bonjour à tous,

J'ai besoin de découper des grandes matrices creuses carrées en sous matrices.
J'ai cru que en écrivant de cette manière : M1=M(debL:finL,debC:finC) (avec M une matrice creuse, debL : indice de début au niveau des lignes dans la matrice M , finL : indice de fin au niveau des lignes dans la matrice M, debC: indice de début au niveau des colonnes dans la matrice M , finC: indice de début au niveau des colonnes dans la matrice M)
je peux extraire une partie de la matrice.
Mais le résultat était M1=la matrice extraite de M avec le nombre de lignes = finL-debL+1 et le nombre de colonnes = = finC-debC+1 et l'indice de début des lignes et des colonnes est toujours 1.

Est ce qu'il ya une fonction de Matlab qui peut extraire une partie de la matrice?

Cordialement
Asma

[Invité]

Bonjour,

Mais le résultat était M1=la matrice extraite de M avec le nombre de lignes = finL-debL+1 et le nombre de colonnes = = finC-debC+1 et l'indice de début des lignes et des colonnes est toujours 1.

Ce résultat est tout à fait normal, je ne vois pas le problème...
Contrairement à d'autres langages comme le C, MATLAB fait débuter ses indices à 1.
De plus si tu fais M1 = M(debL:finL , debC:finC); il est tout à fait normal que tu aies fin-deb+1 éléments.
1:10 %=> 1 2 3 4 5 6 7 8 9 10 compte bien 10 valeurs égal à 10-1+1.

[assma.m]

Le problème est qu'il me donne pas les éléments entre les indices (debL et finL ; debC et finC ) Il me donne toujours le résultat des premières lignes et colonnes.

[Clotho]

Tu pourrais donner un code succint avec une matrice 5,5 par exemple? Je vois mal quel est le soucis, sans résultats concrets.

J'ai utilisé ce type d'écriture récemment pour extraire une sous matrice d'une image et la superposer dans une autre sans soucis de valeurs. Et ce que tu as employé ou ce qui a été donné par Winjerome marche bien.

[Jerome Briot]

Je ne suis pas sûr de bien comprendre le problème mais ne chercherais-tu pas à faire ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
M = rand(4);
M = (M>.3).*M
 
S = sparse(M);
 
debL = 2;
finL = 3;
debC = 2;
finC = 4;
 
[i,j,v] = find(S(debL:finL,debC:finC));
 
newS = sparse(i+debL-1,j+debC-1,v)