Discussion :

[Newenda]

Bonjour,

J'aimerais faire un mix de plotyy et plotxx, j'ai en effet 2 valeurs parfaitement linéaires entre elles en abscisses (le temps et un valeur d'angle) et pour chaque un valeur de vitesse (1er ordonnée) et de hauteur (2eme ordonnée) que je mets en ordonnée.

Mais je ne vois pas comment faire et pas trouver d'aide non plus sur le net.

De plus, comme une des abscisses est le temps, j'aimerais bien les afficher en utilisant Datetick.

Si vous pouviez me filer un coup de main, ça serait très agréable.

Merci d'avance.

[duf42]

Bonjour,

La solution est de superposer 2 axes, tu peux t'inspirer du code des fonctions PLOTYY et/ou PLOTXX.
Une fois que tu as créé 2 axes superposés, tu peux les régler à ta guise pour obtenir ce que tu veux (et utiliser DATETICK sans problème).

Duf

[Newenda]

Bonjour,

La solution est de superposer 2 axes, tu peux t'inspirer du code des fonctions PLOTYY et/ou PLOTXX.
Une fois que tu as créé 2 axes superposés, tu peux les régler à ta guise pour obtenir ce que tu veux (et utiliser DATETICK sans problème).

Duf

Merci Duf pour ta réponse. Mais j'aimerais simplement ajouter un axe en 'x' un peu au dessus de l'axe des abscisses. J'aimerais également qu'il réagisse avec le zoom.

En sommes j'ai une série de mesure où j'ai, à chaque fois :
le temps, un angle, 2 valeurs de hauteur, 1 valeur de vitesse

J'essaye de représenter la hauteur et la vitesse en fonction du temps (donc utilisation de plotty) mais j'aimerais aussi rajouter l'angle en abscisse. Le petit problème étant que l'angle a des valeurs de 0 à 25° puis de 25° à 0° (en gros il monte et descend ensuite)

Bref voilà à quoi je pensais (avec des données au pif):

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
clear all
close all
 
%% data temps
deb = datenum('March 28, 2012  3:37:00 PM')
fin = datenum('March 28, 2012  5:37:00 PM')
time = linspace(deb,fin,41) ;
 
%% data angle
anglee = [0:1:20  19:-1:0] ; % monté et descente
 
%% data hauteur (2 vecteurs)
hauteur1 = rand(1,41);
hauteur2 = rand(1,41);
hauteur = [hauteur1 ; hauteur2 ]
 
%% data vitesse
vitesse = rand(1,41);
 
figure
[AX,H1,H2] = plotyy(time,vitesse, time,hauteur) ;
    set(get(AX(1),'Ylabel'),'String','vitesse')
    set(get(AX(2),'Ylabel'),'String','hauteur')
    set(H1,'Marker','s','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',3)
    set(H2,'LineStyle','-','LineWidth',1.6)
 
pos = get(gca, 'position')+[0 0.12 0 0];
axes('position', pos);
 
set(gca, 'color', 'none', 'ytick', [],'xtick',  [0 10 20 30 40],'XTickLabel',{'0','10','20','10','0'});

Mais voilà, comme vous pouvez l'essayer, le résultat est foireux : - impossible de zoomer, l'axe de l'angle est vide (?), et il y a un espèce de trait qui part du sommet de l'ordonnée de gauche..

J'ai aussi tenté un test avec hline mais bon..c'est pire

Code :

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

1
2
3
4
5
6
7
8
9
10
  coul = hsv(2)
figure
for ii = 1:2
    hl1 = line(time,hauteur(ii,:),'Color',coul(ii,:));
end
ax1 = gca;
 
 
hl2 = line(anglee,vitesse,'Color','g');
ax2 = axes('Position',get(ax1,'Position')+[0 0 0 0],'XAxisLocation','top','YAxisLocation','right','Color','none','XColor','k','YColor','k');

Bref je vois pas trop...

Merci !

[Newenda]

Alors comme ça :

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
35
36
37
clear all
close all
 
%% data temps
deb = datenum('March 28, 2012  3:37:00 PM')
fin = datenum('March 28, 2012  5:37:00 PM')
time = linspace(deb,fin,41) ;
 
%% data angle
anglee = [0:1:20  19:-1:0] ; % monté et descente
 
%% data hauteur (2 vecteurs)
hauteur1 = rand(1,41);
hauteur2 = rand(1,41);
hauteur = [hauteur1 ; hauteur2 ]
 
%% data vitesse
vitesse = rand(1,41);
 
 
 coul = hsv(2)
figure
for ii = 1:2
hl1 = line(time,hauteur(ii,:),'Color',coul(ii,:));
end
ax(1)=gca;
 
 
set(ax(1),'Position',[0.12 0.12 0.75 0.70])
set(ax(1),'XColor','k','YColor','k');
 
ax(2)=axes('Position',get(ax(1),'Position'),'XAxisLocation','top','YAxisLocation','right','Color','none','XColor','k','YColor','k');
 
set(ax,'box','off')
 
hl2 = line(anglee,vitesse,'Color','g','Parent',ax(2));
set(ax(2),'Xtick',[0 5 10 15 20 25 30 35 40],'XTickLabel',{'0','5','10','15','20','15','10','5','0'});

ça marche pas mal mais :

- Le zoom ne marche que sur un axe
- Je vois pas toujours comment faire l'axe des angles de 0° > 25° > 0°
- Pas moyen de placer ce deuxième axe des abscisses non pas en haut mais en bas, au dessus du premier ? (en pouvant ajouter des labels)

[Invité]

Bonjour,

Si j'ai bien compris, tu peux reprendre le même principe que je t'avais conseillé ici.
Par contre, je n'ai pas saisi le problème des angles

[Newenda]

Exactement.. mais à l'époque j'avais eu des problèmes d'axe (je jouais avec pcolor aussi)

Mais sinon oui c'est ça :

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
  yWidth = 0.04;                      %# y axes spacing, 
    xLimit = [min(time) max(time)];         %# Range of time values
    xOffset = -yWidth*diff(xLimit)/axesPosition(3);
 
    %# Create the figure and axes:
 
    % VITESSE
    h1 = axes('Position',axesPosition,...
        'Color','w','XColor','k','YColor','r',...
        'XLim',xLimit,'YLim',[0 20],'NextPlot','add','FontSize',12,'FontWeight','bold');
          set(h1, 'Xtick', [xLimit(1):1/24/5:xLimit(2)]);
    datetick('x','HH:MM','keeplimits','keepticks')    
 
    % HAUTEUR
    h2 = axes('Position',axesPosition+yWidth.*[-1 0 1 0],...
        'Color','none','XColor','k','YColor','k',...
         'XLim',xLimit+[xOffset 0],'YLim',[-30 30],...
   'XTick',[],'XTickLabel',[],'NextPlot','add','FontSize',12,'FontWeight','bold');
 
    % ANGLE
    h3 = axes('Position',axesPosition+yWidth.*[-2 0 2 0],...
        'Color','none','XColor','k','YColor','b',...
         'XLim',xLimit+[2*xOffset 0],'YLim',[-2 30],...
       'XTick',[],'XTickLabel',[],'NextPlot','add','FontSize',12,'FontWeight','bold');
 
       set(h1, 'Xtick', [xLimit(1):1/24/5:xLimit(2)]);
    datetick('x','HH:MM','keeplimits','keepticks')    %
 
    xlabel(h1,'temps','FontSize',12,'FontWeight','bold');
    ylabel(h1,'vitesse (m/s)','FontSize',12,'FontWeight','bold');
    ylabel(h2,'hauteur (mm)','FontSize',12,'FontWeight','bold');
    ylabel(h3,'angle (°)','FontSize',12,'FontWeight','bold');

Pour l'angle, comme ce n'est pas une échelle classique linéaire je l'ai ploté comme une donnée et non un axe.

Merci !