Discussion :

[francois_S]

Bonjour,
J'ai un simple problème qui m'énerve.
J'ai une matrice assez conséquente A , et je dois mettre certains indices de sa diagonale à une certaine valeur, disons à 1, par exemple.

J'ai un vecteur d'indice ind, qui contient les indices de la diagonale.

J'ai donc par exemple,

Code :

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

A= rand(1000,1000); ind=[1 4 100 330];

Et quand j'essaie quelque chose de simple comme A(ind,ind) =1.. bah ca me met tout le carré correspondant à 1, ce qui n'est pas ce que je désire, vu que je veux JUSTE les éléments de la diagonale correspondant aux indices donnés à 1...
je veux donc A(ind(1),ind(1)) = 1 puis A(ind(2),ind(2)) =1, etc...
J'aimerais me passer de boucler for pour le faire.

Quelqu'un serait m'aider?
Un grand merci d'avance

[Jean Dumoncel]

Bonsoir,

une solution en bidouillant les matrices :

Code :

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

1
2
3
d = diag(A);
d(ind) = 1;
result = A-diag(diag(A))+diag(d);

ou un peu plus propre avec les indices :

Code :

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

A(sub2ind(size(A),ind,ind)) = 1;

[francois_S]

Un grand merci!

la solution avec le sub2ind est exactement le genre de solution que je cherchais!

[Jerome Briot]

Code :

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

A(sub2ind(size(A),ind,ind)) = 1;

La même solution sans SUB2IND :

Code :

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

1
2
3
4
A = rand(7);
ind = [1 2 5];
 
A((ind-1)*size(A,1)+ind) = 1

Voir la technique d'indexage linéaire dans la