Discussion :

[Rorchar]

Bonjour,

Voici ce que j'aimerai réalisé :

Considérons une matrice A de 10 colonnes correspondant à 10 variables et de 100 lignes correspondant au temps.
Le plot consiste donc à tracer l'évolution de ces variables selon le temps : plot(t,A) avec t un vecteur contenant les 100 incrémentation de temps
Cependant, j'aimerai pouvoir déterminer par exemple 1 couleur pour chaque groupe de 2 variables, et au sein de chaque groupe de couleur un type de tracé (pointillé, ligne continu...) pour chaque variable.

Je peux obtenir le résultat voulu en faisant un plot(...) à rallonge mais je cherche une solution plus optimisée ?
L'idée étant aussi d'avoir au final une solution plastique, qui peut s'adapter automatiquement à un nombre différent de variables, ou une taille de groupe différente

[phryte]

Bonjour,

Une idée pour 4 colonnes (avec une boucle hélas !) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
clear
Nl=10;
M=zeros(Nl,4);
M(:,1)=1:Nl;
M(:,2)=2*(1:Nl);
M(:,3)=3*(1:Nl);
M(:,4)=4*(1:Nl);
C=['sr';'*r';'.b';'db']
hold on
grid
for i=1:4
    plot(1:Nl,M(:,i),C(i,:))
end

[Rorchar]

Intéressant merci, c'est nettement plus plastique que ce que j'avais pour le moment.

Une boucle sur un plot c'est quelque chose à éviter en principe ? ça peut rapidement devenir lent ?

[Jerome Briot]

Une boucle sur un plot c'est quelque chose à éviter en principe ? ça peut rapidement devenir lent ?

Dans ce cas, tu peux utiliser cette variante :

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
clear
Nl = 10;
M = zeros(Nl,4);
M(:,1) = 1:Nl;
M(:,2) = 2*(1:Nl);
M(:,3) = 3*(1:Nl);
M(:,4) = 4*(1:Nl);
 
lst = {'none' 'none' 'none' 'none'};
col = {'r' 'r' 'b' 'b'};
mrk = {'s' '*' '.' 'd'};
 
figure
 
h = plot(M);
 
grid on
 
for n = 1:numel(h)
    set(h(n), 'linestyle', lst{n}, 'color', col{n}, 'marker', mrk{n})
end

[Jerome Briot]

Une autre solution sans boucle mais qui permet de gérer uniquement la couleur des courbes :

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
clear
Nl = 10;
M = zeros(Nl,4);
M(:,1) = 1:Nl;
M(:,2) = 2*(1:Nl);
M(:,3) = 3*(1:Nl);
M(:,4) = 4*(1:Nl);
 
col = [1 0 0 ; 1 0 0 ; 0 0 1 ; 0 0 1];
 
figure
 
set(groot, 'DefaultAxesColorOrder', col)
 
plot(M);
 
set(groot, 'DefaultAxesColorOrder', 'remove')
 
grid on

L'objet groot est disponible depuis la version R2014b de MATLAB. Pour les versions précédentes, il faut essayer avec :

Code :

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

1
2
3
set(0, 'DefaultAxesColorOrder', col)
 
set(0, 'DefaultAxesColorOrder', 'remove')

[phryte]

Bonjour,

Une boucle sur un plot c'est quelque chose à éviter en principe ? ça peut rapidement devenir lent ?

En général utiliser des boucles est une solution de facilité.

MATALB a été conçu pour les calculs matriciels : MATrix LABoratory.

Utiliser les vecteurs est donc plus élégant et plus rapide. Mais parfois on ne peut pas faire autrement !

[Jerome Briot]

Ne pas confondre le noyau de calcul de MATLAB et le moteur graphique.

Pour optimiser le premier, il faut utiliser la pré-allocation de mémoire.
Pour le deuxième, il faut diminuer le nombre d'appel aux fonctions de traçage.

Exemple d'un code non optimisé :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
figure
 
hold on
 
for n = 1:10000
 
    M(n,:) = rand(1,500);
 
    plot(1:500, M(n,:))
 
end

Si on optimise le calcul :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
figure
 
hold on
 
M = zeros(10000, 500);
 
for n = 1:10000
 
    M(n,:) = rand(1,500);
 
    plot(1:500, M(n,:))
 
end

Si on optimise le tracé :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
M = zeros(10000, 500);
 
for n = 1:10000
 
    M(n,:) = rand(1,500);
 
end
 
figure
 
plot(M)

C'est un exemple théorique. Il faut bien entendu s'adapter à chaque cas spécifique.

[Rorchar]

Merci Jérôme, très pratique !