Discussion :

[pierrecvvb]

Bonjour à tous !

je suis ravi de venir poser à nouveau mes questions sur ce forum.
Après plusieurs recherches (forum, FEX, etc...) je n'ai pas trouvé mon bonheur.
Je souhaiterai tracer une barre d'erreur sur un graphique 3D en barres. (fonction bar3)
https://fr.mathworks.com/matlabcentr...ith-error-bars
Cette fonction en lien n'est pas des meilleures et contient des erreurs...
Voici mon cas :
J'ai deux matrices (8 lignes x 3 colonnes).

• M = 373,39 460,73 482,15
  339,58 458,79 473,45
  235,40 358,51 333,25
  220,50 273,54 476,94
  189,35 99,230 453,25
  233,99 85,02 296,37
  116,17 101,34 139,77
  79,88 86,76 146,68
  
• M_error = 27,82 57,95 45,13
  33,48 41,80 54,24
  40,33 48,45 44,16
  18,26 14,80 32,60
  18,79 14,09 31,52
  27,9 9,65 24,98
  11,94 14,87 9,23
  4,87 8,32 9,38
  

Je trace la première matrice (M) ce qui donne le résultats ci dessous :

Code :

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

Graph3D = bar3(M)

Je souhaiterai donc ajouter la barre d'erreur exposée dans la seconde matrice à chaque barre de ce graphique. Chaque valeur de la matrice M_error correspond à l'amplitude d'erreur de M.
Connaissez-vous un moyen d'y parvenir ?
Je vous remercie de votre aide.

Cordialement,
Pierre L

[Jerome Briot]

Essaie ceci :

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
26
27
28
29
30
31
32
33
34
M = [373.39 460.73 482.15
    339.58 458.79 473.45
    235.40 358.51 333.25
    220.50 273.54 476.94
    189.35 99.230 453.25
    233.99 85.02 296.37
    116.17 101.34 139.77
    79.88 86.76 146.68];
 
M_error = [27.82 57.95 45.13
    33.48 41.80 54.24
    40.33 48.45 44.16
    18.26 14.80 32.60
    18.79 14.09 31.52
    27.9 9.65 24.98
    11.94 14.87 9.23
    4.87 8.32 9.38];
 
figure
 
bar3(M)
hold on
 
for n = 1:size(M,2)
 
    for m = 1:size(M,1)
        line([n n], [m m], [M(m,n) M(m,n)+M_error(m,n)], ...
            'color', 'r', ...
            'linewidth', 2)
    end
 
end
 
axis tight

[pierrecvvb]

Merci beaucoup pour cette réponse rapide !
C'est en effet bien ce que je cherche à faire.
la fonction :

Code :

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

line

est bien plus évidente que celle que j'ai essayé :

Code :

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

plot3

Le résultat est à peu de choses près identique.
Cependant, je trouve que cela ne ressemble pas à une vraie barre d'erreur comme la fonction 2D :

Code :

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

errorbar

Pour un lecteur extérieur ce n'est pas évident.
J'ai donc pensé à rendre le plot transparent mais la barre d'erreur ne s'affiche alors qu'en ligne sur la borne supérieure.

Code :

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

1
2
3
4
5
6
 
%Transparence des barres
 Graph3D = bar3(M,0.5,'detached') ; 
    set(M(1),'facecolor',[0.298 0.298 0.298],'FaceAlpha',.5); 
    set(M(2),'facecolor',[0 0.1647 0.3275],,'FaceAlpha',.5);
    set(M(3),'facecolor',[0.6267 0.1412 0.4],,'FaceAlpha',.5);

D'où mes deux prochaines questions :
Y a t'il un moyen de délimiter ces lignes d'erreurs de la même façon qu'en 2D avec la fonction errorbar ?
Comment visualiser proprement l'amplitude de la ligne d'erreur dans la barre ?
Encore merci.

Pierre L

[Jerome Briot]

Y a t'il un moyen de délimiter ces lignes d'erreurs de la même façon qu'en 2D avec la fonction errorbar ?

Tu peux ajouter un marqueur à l'extrémité de la barre d'erreur :

Code :

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

1
2
3
4
5
6
7
8
9
10
for n = 1:size(M,2)
 
    for m = 1:size(M,1)
        line([n n], [m m], [M(m,n) M(m,n)+M_error(m,n)/2], ...
            'color', 'r', ...
            'linewidth', 2)
        plot3(n, m, M(m,n)+M_error(m,n)/2, 'or', 'markersize', 6, 'markerfacecolor', 'r')
    end
 
end

Ou bien de lignes comme avec errorbar :

Code :

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

1
2
3
4
5
6
7
8
9
for n = 1:size(M,2)
 
    for m = 1:size(M,1)
        line([n n-c ; n n+c], [m m; m m], [M(m,n) M(m,n)+M_error(m,n)/2 ; M(m,n)+M_error(m,n)/2 M(m,n)+M_error(m,n)/2], ...
            'color', 'r', ...
            'linewidth', 2)
    end
 
end

ou encore :

Code :

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

1
2
3
4
5
6
7
8
9
for n = 1:size(M,2)
 
    for m = 1:size(M,1)
        line([n n-c n-c ; n n+c n+c], [m m+c m-c; m m-c m+c], [M(m,n) M(m,n)+M_error(m,n)/2 M(m,n)+M_error(m,n)/2 ; M(m,n)+M_error(m,n)/2 M(m,n)+M_error(m,n)/2 M(m,n)+M_error(m,n)/2], ...
            'color', 'r', ...
            'linewidth', 2)
    end
 
end

Comment visualiser proprement l'amplitude de la ligne d'erreur dans la barre ?

Il suffit de remplacer ceci :

Code :

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

[M(m,n) M(m,n)+M_error(m,n)]

par ceci :

Code :

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

[M(m,n)-M_error(m,n)/2 M(m,n)+M_error(m,n)/2]

[pierrecvvb]

Un très grand ... MERCI !