Discussion :

[linda05]

bonjour a tous,
Voilà mon problème: j'ai une matrice par exemple A d'ordre 4( Sachant que j'ai des matrices de très grandes tailles)

Code :

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

A=[3 4 5 2;1 2 4 6;3 4 8 7;2 4 3 7]

et une matrice B contient des indices de lignes et de colonnes de la matrice A

Code :

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

B=[2 3 4 1;4 3 2 1;3 2 1 4;4 2 3 1;1 2 3 4]

je cherche une matrice C :

Code :

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

C=[C(1,:) B]

tel que:

Code :

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

1
2
3
4
5
C(1,1)= A(2,3)+A(3,4)+A(4,1)
C(2,1)= A(4,3)+A(3,2)+A(2,1)
C(3,1)= A(3,2)+A(2,1)+A(1,4)
C(4,1)= A(4,2)+A(2,3)+A(3,1)
C(5,1)= A(1,2)+A(2,3)+A(3,4)

ie

Code :

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

C=[13 2 3 4 1;8 4 3 2 1;7 3 2 1 4;11 4 2 3 1;15 1 2 3 4]

?

J'aimerai surtout trouver un moyen de faire ça d'une façon propre et optimisée.
merci beaucoup
a bientôt

[Invité]

Bonjour,

Si j'ai bien compris :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
A = [ 3 4 5 2
      1 2 4 6
      3 4 8 7
      2 4 3 7 ];
 
B = [ 2 3 4 1
      4 3 2 1
      3 2 1 4
      4 2 3 1
      1 2 3 4 ];
 
rows = B(:,1:3);
cols = B(:,2:4);
C = [sum( A( sub2ind(size(A),rows,cols) ) ,2) B]

[linda05]

Je tiens à te remercier,
J’ai d’autres questions:

• _ si on a un ensemble des éléments comment faire pour trouver l’élément qui est juste avant un autre par exemple pour

Code :

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

X=1 :5 ;X(4)=[] ;

donc

Code :

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

X=[1 2 3 5]

je veux l'élément qui est juste avant 5 la c’est 3?
mon code:

Code :

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

X(find(X==5)-1)

• pour changer deux éléments ds une matrice J'ai essayé ce code :

pour transformer les valeurs "5" en valeurs "3" et les valeurs "3" en valeurs "5:":

Code :

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

1
2
3
4
A=[1 3 5 ;2 5 6 ;3 7 8]
A(A==3)=NaN ;
A(A==5)=3
A(isnan(A))=5

donc:

Code :

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

A=[1 5 3;2 3 6;5 7 8]

Y'a t il moyen de faire ca plus simple?
merci encore une fois, et bonne fin de semaine.

[Newenda]

Il me semble que ce soit impossible d'intervertir deux nombres.. peut être une fonction dans matlabexange mais j'ai rien vu.

Sinon, pour être plus élégant (plus compréhensible) on peut utiliser "find" mais ça ne reduira pas ton code bien sûr.

[Invité]

1. Oui, je rajouterais X(find(X==5, 1, 'first')-1).
2. Pour éviter de jongler avec les NaN, on peut retenir le tableau logique localisant les 3 comme ceci :
   

   Code :

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

   1
   2
   3
   4
   A=[1 3 5 ;2 5 6 ;3 7 8];
   Aeq3 = A==3;
   A(A==5) = 3;
   A(Aeq3) = 5;

[Jerome Briot]

Deux alternatives pour transformer les 3 en 5 et inversement :

Code :

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

A = A.*(A~=3 & A~=5)+3*(A==5)+5*(A==3);

ou plus simplement :

Code :

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

A = A-2*(A==5)+2*(A==3);

Mais le code perd en lisibilité par rapport aux solutions classiques données précédemment. Il faut donc surement l'accompagner d'un commentaire.

[linda05]

merci bien,
mais pour le code

Code :

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

A = A-2*(A==5)+2*(A==3)

est ce que dans le cas général :

Code :

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

A = A-(a-b)*(A==a)+(a-b)*(A==b)

pour a,b des éléments dans A.
bonne journée

[Jerome Briot]

est ce que dans le cas général :

Code :

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

A = A-(a-b)*(A==a)+(a-b)*(A==b)

pour a,b des éléments dans A.

Oui si a et b sont des entiers

[linda05]

Ok, merciiiiiiiiiiiiiiiiiiiiiii