Discussion :

[pablot]

Bonjour a tous,

Le but de mon bout de programme est d'afficher 3 angles évoluant dans le temps, en direct. J'ai donc 2 boucles for imbriquées les unes dans les autres :

La première (grande) parcourt les différentes valeurs d'angles
La deuxième parcourt un certain nombre de rayon pour un angle donné par la première boucle for (afin de tracer le rayon)

Mes problèmes :

- Est il possible de tracer un trait continu plutôt que des cercles comme j'ai fais ?
- Mes deux premiers cercles ne prennent pas en compte le hold off, et je conserve donc la "trace" des rayons, censée s'effacer a chaque changement d'angle, pour suivre l'évolution des angles dans le temps clairement.

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
 
 
for k=1:20:(length(roll)) %boucle de parcours des valeurs d'angle
 
    for i=0:0.1:1; %boucle de tracage de rayon
 
        subplot(3,2,2)
        polar(yaw(1,k),i,'--or') %plot de la valeur numero k de l'angle yaw, pour un rayon i
 
        hold on
 
        subplot(3,2,4)
        polar(pitch(1,k),i,'--or') %plot de la valeur numero k de l'angle pitch, pour un rayon i
 
        hold on
 
        subplot(3,2,6)
        polar(roll(1,k),i,'--or') %plot de la valeur numero k de l'angle roll, pour un rayon i
 
        hold on
 
    end;
 
    pause(0.0001)
    hold off
 
end

Voici ce que j'obtiens apres avoir parcouru tous les angles :



Merci a tous pour votre aide !

[Jerome Briot]

Voici deux variantes selon que les valeurs sont calculées avant ou pendant les itérations :

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
r = 0:0.1:1;
th = 2*pi*r;
[x,y] = pol2cart(th, r);
 
figure
 
p = polar(NaN,max(r),'ro-');
 
for n = 1:numel(x)
 
    set(p, 'XData', x(1:n), 'YData', y(1:n));    
 
    pause(0.2)
 
end

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
figure
 
p = polar(NaN,1,'ro-');
 
for r = 0:0.1:1
 
    th = 2*pi*r;
 
    [x,y] = pol2cart(th, r);
 
    xp = get(p, 'XData');
    yp = get(p, 'YData');
 
    set(p, 'XData', [xp x], 'YData', [yp y]);    
 
    pause(0.2)
 
end

[pablot]

Merci beaucoup je vais regarder ça de plus près !

[pablot]

J'ai du mal m'exprimer,

En fait je cherche a représenter l'angle calculé avant comme une aiguille de compteur, c'etait le but de ma seconde fonction for, tracer le rayon.
Donc a chaque changement d'angle, l'angle précédent doit disparaitre, ce qui n'est pas le cas.

[Jerome Briot]

Comme 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
% Données à tracer
r = 0:0.1:1;
th = 2*pi*r;
 
% Conversion en coordonnées cartésiennes
[x,y] = pol2cart(th, r);
 
figure
 
% Tracé "vide" et stockage de l'identifiant de la ligne dans la variable p
p = polar(NaN,max(r),'r-');
 
for n = 1:numel(x)
 
    % Mise à jour des coordonnées des deux extrémités de la ligne
    set(p, 'XData', [0 x(n)], 'YData', [0 y(n)]);    
 
    pause(0.2)
 
end

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
figure
 
% Tracé "vide" et stockage de l'identifiant de la ligne dans la variable p
p = polar(NaN,1,'r-');
 
for r = 0:0.1:1
 
    % Valeur de theta variable
    th = 2*pi*r;
 
    % Conversion en coordonnées cartésiennes
    [x,y] = pol2cart(th, r);
 
    % Mise à jour des coordonnées des deux extrémités de la ligne
    set(p, 'XData', [0 x], 'YData', [0 y]);    
 
    pause(0.2)
 
end

[pablot]

Oui c'est l'idée merci !
Tu pourrais commenter ton code stp car je ne comprends pas tout loin de la...

Merci !

[Jerome Briot]

Tu pourrais commenter ton code stp car je ne comprends pas tout loin de la...

C'est fait.

N'hésite pas à poser des questions sur les lignes que tu ne comprends pas.

[pablot]

Merci !
J'ai finis par m'en sortir et voila que ca fonctionne au top comme ceci (pour ceux que ca peut aider) :

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
 
 
subplot(3,2,2)
p1 = polar(NaN,1,'r-')
 
subplot(3,2,4)
p2 = polar(NaN,1,'r-')
 
subplot(3,2,6)
p3 = polar(NaN,1,'r-') 
 
[x1,y1] = pol2cart(yaw, 1);
[x2,y2] = pol2cart(pitch, 1);
[x3,y3] = pol2cart(roll, 1);
 
for k=1:(length(roll)-1) %boucle de parcours des valeurs d'angle
 
 
set(p1, 'XData', [0 x1(k)], 'YData', [0 y1(k)]);    
set(p2, 'XData', [0 x2(k)], 'YData', [0 y2(k)]);   
set(p3, 'XData', [0 x3(k)], 'YData', [0 y3(k)]);   
 
    pause(0.1)
 
end

Je ne connaissais pas la fonction set mais c'est exactement ce qu'il me fallait !

Est il possible de changer l'épaisseur du trait ? Car j'ai eu l'impression que les parametres LineSpec ne passaient pas comme avec la fonction plot ?

Merci encore

[pablot]

Au passage voici le resultat :

[Jerome Briot]

Est il possible de changer l'épaisseur du trait ?

Comme 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
figure
subplot(3,2,2)
p(1) = polar(NaN, 1, 'r-');
 
subplot(3,2,4)
p(2) = polar(NaN, 1,'r-');
 
subplot(3,2,6)
p(3) = polar(NaN, 1, 'r-');
 
set(p, 'LineWidth', 2)
 
[x1,y1] = pol2cart(yaw, 1);
[x2,y2] = pol2cart(pitch, 1);
[x3,y3] = pol2cart(roll, 1);
 
for k = 1:(length(roll)-1) %boucle de parcours des valeurs d'angle
 
    set(p(1), 'XData', [0 x1(k)], 'YData', [0 y1(k)]);
    set(p(2), 'XData', [0 x2(k)], 'YData', [0 y2(k)]);
    set(p(3), 'XData', [0 x3(k)], 'YData', [0 y3(k)]);
 
    pause(0.1)
 
end

Tu remarqueras qu'il est plus avantageux d'utiliser un vecteur d'identifiants ici pour appliquer directement l'épaisseur à toutes les lignes

[pablot]

C'est vraiment parfait merci beaucoup !
Je vais prendre des infos sur la fonction set car elle a l'air tres utile !

[tachmou]

Salut,

Je te conseille de tester :

Code :

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

set(p)

Code :

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

get(p)

Cela te permettra d'accéder à la liste de toutes les propriétés disponibles (via get) et modifiables (via set) de ta structure "p".

De plus, des petits :

Code :

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

doc set

Code :

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

doc get

... t'aideront sûrement :-)

Cordialement,

[pablot]

Merci bien

[pablot]

Je reviens a vous car j'ai adapté mon troisieme graphe (une représentation 3D d'un vecteur) avec la fonction set, ca fonctionne nikel et c'est plus léger.
Mais je ne vois pas comment afficher 3 vecteurs sur le meme graphe en utilisant cette fonction...

Voici mon code :

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
 
 
subplot(3,2,[1 3 5])
g = plot3(NaN,NaN,NaN)
 
grid on
axis ([-100 100 -100 100 -100 100])
 
set(p, 'LineWidth', 2)
set(g, 'LineWidth', 2)
 
pause(2)
 
for k=1:(length(roll)-1) %boucle de parcours des valeurs d'angle
 
    set(p(1), 'XData', [0 x1(k)], 'YData', [0 y1(k)]);
    set(p(2), 'XData', [0 x2(k)], 'YData', [0 y2(k)]);
    set(p(3), 'XData', [0 x3(k)], 'YData', [0 y3(k)]);
    set(g, 'XData', [0 vector_new_x(k,1)] ,'YData', [0 vector_new_x(k,2)], 'ZData', [0 vector_new_x(k,3)]) %affichage du vecteur 3D
 
 
    pause(0.1)
 
end

Merci !!!

[tachmou]

Il te suffit de créer d'autres vecteurs de la même manière, mais lorsque tu les crées, il faut que tu utilises :

Code :

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

hold on

... ou bien, dit autrement, que tu "set" le "hold" du handle de ta figure à "on".

Cordialement,

[Jerome Briot]

... ou bien, dit autrement, que tu "set" le "hold" du handle de ta figure à "on".

Un survol du code de hold permet de s'en convaincre :

Code :

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

edit hold

[pablot]

Merci pour la reponse mais je n'ai pas compris.

La pour un vecteur, j'ai set les valeurs de Xdata Ydata et Zdata de mon vecteur 1, sur g, qui est ma figure.
D'apres ce que vous me dites, je laisse tel quel, sauf que je rajoute hold on a la fin comme pour un plot classique, et ensuite je re set g de la meme maniere avec mes nouveaux vecteurs ?

[tachmou]

En gros, c'est ça, oui.

"Hold on" te permet d'afficher des nouveaux objets dans la même zone graphique (via p(..) = polar...), sans effacer les précédents.

Ensuite, le fait d'utiliser "set" te permet de mettre à jours ces mêmes objets, donc tu n'en afficheras pas d'autres avec "set".

[pablot]

J'ai tenté ceci mais je n'ai qu'un vecteur qui apparait :

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
 
 
for k=1:(length(roll)-1) %boucle de parcours des valeurs d'angle
 
    set(p(1), 'XData', [0 x1(k)], 'YData', [0 y1(k)]); %plot cercles trigos
    set(p(2), 'XData', [0 x2(k)], 'YData', [0 y2(k)]);
    set(p(3), 'XData', [0 x3(k)], 'YData', [0 y3(k)]);
 
    set(g, 'XData', [0 vector_new_x(k,1)] ,'YData', [0 vector_new_x(k,2)], 'ZData', [0 vector_new_x(k,3)]), hold on %plot vecteur x
    set(g, 'XData', [0 vector_new_y(k,1)] ,'YData', [0 vector_new_y(k,2)], 'ZData', [0 vector_new_y(k,3)]), hold on %plot vecteur y
    set(g, 'XData', [0 vector_new_z(k,1)] ,'YData', [0 vector_new_z(k,2)], 'ZData', [0 vector_new_z(k,3)]), hold on %plot vecteur z
 
    pause(0.1)
 
end

[pablot]

J'ai egalement essayé de le mettre apres mon g = plot3[...] mais toujours rien.

La chose importante aussi est que je dois afficher 3 vecteurs différents a chaque itération de ma boucle for, le hold ne doit se faire donc que pendant une itération...