Discussion :

[bernard6]

Bonjour,

J'ai un petit problème pour combiner deux matrices de manière spéciale.

Je m'explique...

Soient :

Code :

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

1
2
3
 
A=[92 107; 107 122; 122 13; 13 92]
B=[3; 1;0;47]

J'aimerai bien avoir comme résultat, en les combinant :

Code :

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

1
2
 
C=[92 3; 3 107; 107 1; 1 122; 122 13; 13 47; 47 92]

c'est à dire que j'intercale toutes les deux lignes une ligne de B dans A SAUF si la valeur de B est égale à 0.

Je cherche une technique pour que ce soit fonctionnel pour des matrices de toute taille (A et B possédant toujours le même nombre de lignes)

Voilà, je ne sais pas si c'était très clair...

Merci d'avance pour votre aide !

[Jerome Briot]

Une première solution :

Code :

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

1
2
3
4
5
A=[92 107; 107 122; 122 13; 13 92]
B=[3; 1;0;47]
 
C=reshape([A(:,1) B B A(:,2)].',2,[]).';
C=[C(C(:,1)~=0,1) C(C(:,2)~=0,2)]

[bernard6]

Bonjour,

Votre première solution fonctionne parfaitement !

Merci !

[Jerome Briot]

une deuxième solution:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
A=[92 107; 107 122; 122 13; 13 92]
B=[3; 1;0;47]
 
C=kron(A,[1;1]);
idx=2*find(B);
 
[C(idx,1),C(idx-1,2)]=deal(B(B~=0));
 
idx=2*find(~B);
C(idx,:)=[];
 
C

Comparer le temps d'exécution entre les deux codes

[bernard6]

Bonjour, en reprenant des données plus conséquentes (A et B ont 401 lignes) et en exécutant le code suivant plusieurs fois pour chacune des techniques, je me rends compte que, de façon générale, le deuxième code est légèrement plus rapide...

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
 
 
PointsVitesse = csvread('donnees vitesse.txt')
t0=cputime; 
 
A=[PointsVitesse(:,1),PointsVitesse(:,2)]
B=PointsVitesse(:,3)
 
C=reshape([A(:,1) B B A(:,2)].',2,[]).';
C=[C(C(:,1)~=0,1) C(C(:,2)~=0,2)]
 
disp(['le temps d''execution est ; ',num2str(cputime -t0)])
 
 
% 
% C=kron(A,[1;1]);
% idx=2*find(B);
% 
% [C(idx,1),C(idx-1,2)]=deal(B(B~=0));
% 
% idx=2*find(~B);
% C(idx,:)=[];
% 
% C
% disp(['le temps d''execution est ; ',num2str(cputime -t0)])

Mode pour code 1 = 0.046875
Mode pour code 2 = 0.03125

Il faudra donc préférer la deuxième technique

Merci !