Discussion :

[nant44]

Bonjour,

J'ai deux questions et j'aimerai que qqn pourra me répondre.
1) comment afficher la courbe en bleu (dans le code ci-dessous) uniquement qd th varie entre pi/2 et pi ?

2) Le 2eme problème c'est je n'arrive pas à créer un fichier .avi de l'animation !! [j'ai mis la commande mov(n-1) = getframe; juste avant drawnow et j'ai mis movie2avi(mov,'exemple_anim.avi') après end. mais cela ne marche pas (voila l'erreur Warning: The frame height has been padded to be a multiple of four as required by the specified codec.)

merci d'avance

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
th=0:0.02:2*pi;
y=sin(th);
fig=figure;
set(fig,'doublebuffer','on')
p=plot(th(1),y(1),'r-');
pp=line([0 th(1)],[-1 y(1)]);
xlim([min(th) max(th)])
ylim([min(y) max(y)])
for n=2:numel(th)
  set(p,'xdata',th(1:n),'ydata',y(1:n));
  set(pp,'xdata',[0 th(n)],'ydata',[-1 y(n)]);
  drawnow 
end

[Jerome Briot]

Pour la première question.... j'ai pas compris

Pour la seconde, il semble que ce soit un problème de compatibilité du codec utilisé. As-tu essayé avec un autre codec ?

[nant44]

Merci Dut.
c'est vrai que la première question n'est pas bien formulée !
Donc je reformule : j'ai modifié le code précédent pour éclaircir les choses.

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
clear all
close all
clc
 
th=0:0.005:2*pi;
y=sin(th);
th_begin = pi/2;
th_end = 3*pi/2;
begin_index = find(abs(th-th_begin)==min(abs(th-th_begin)));
end_index = find(abs(th-th_end)==min(abs(th-th_end)));
fig=figure;
hold on
plot(th,y)
set(fig,'doublebuffer','on')
p=plot(th(1),y(1),'rp');
p1=line([th(begin_index) th(begin_index)],[-1 y(begin_index)],'color','green')
xlim([min(th) max(th)])
ylim([min(y) max(y)])
for n=1:numel(th)
  set(p,'xdata',th(n),'ydata',y(n));
  if (n>begin_index)&(n<end_index)
      set(p1,'xdata',[th(n) th(n)],'ydata',[-1 y(n)]);
  end
      drawnow 
end

Dans ce nouveau code, j'ai une courbe sinusoïdale sur laquelle se déplace un point. Ce que je souhaite faire c'est de marquer le passage de mon point entre deux points de la courbe. C'est pour cela, j'ai ajouté un segment de droite dont l'une des extrémités est le point en question et la seconde extrémité sera la projection de ce point sur l'axe des abscisses.

Mon problème c'est je ne vois pas de moyen pour synchroniser l'affichage du segment avec le passage de mon point dans la zone d'intérêt (comme vous pouvez le voir l'affichage de segment commence dès le début de l'animation ce qui normal car c'est que je lui demande de faire dans la ligne "p1=line ....")


Pour la deuxième question je ne vois pas comment choisir un autre codec. Peux-tu me donner une piste ?

Merci

[Jerome Briot]

Tout d'abord :

Code :

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

1
2
begin_index = find(abs(th-th_begin)==min(abs(th-th_begin)))
end_index = find(abs(th-th_end)==min(abs(th-th_end)))

se simplifie comme ceci :

Code :

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

1
2
[pasbesoin,begin_index] = min(abs(th-th_begin));
[pasbesoin,end_index] = min(abs(th-th_end));

Et ensuite, tu peux faire 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
26
27
28
29
30
31
clear all
close all
clc
 
th=0:0.005:2*pi;
y=sin(th);
th_begin = pi/2;
th_end = 3*pi/2;
 
[pasbesoin,begin_index] = min(abs(th-th_begin));
[pasbesoin,end_index] = min(abs(th-th_end));
 
fig=figure;
hold on
plot(th,y)
set(fig,'doublebuffer','on')
p=plot(th(1),y(1),'rp');
p1=line([nan nan],[nan nan],'color','green');
xlim([min(th) max(th)])
ylim([min(y) max(y)])
 
vis = {'off' 'on'};
for n=1:numel(th)
 
	set(p,'xdata',th(n),'ydata',y(n));
 
	v = ((n>begin_index)&(n<end_index)) + 1;
    set(p1,'xdata',[th(n) th(n)],'ydata',[-1 y(n)],'visible',vis{v});
 
    drawnow 
end

Pour le problème de codec, regarde le paramètre Compression de la fonction MOVIE2AVI

[Jerome Briot]

Tu peux aussi faire directement 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
26
27
28
clear all
close all
clc
 
th=0:0.005:2*pi;
y=sin(th);
th_begin = pi/2;
th_end = 3*pi/2;
 
fig=figure;
hold on
plot(th,y)
set(fig,'doublebuffer','on')
p=plot(th(1),y(1),'rp');
p1=line([nan nan],[nan nan],'color','green');
xlim([min(th) max(th)])
ylim([min(y) max(y)])
 
vis = {'off' 'on'};
v = ((th>th_begin) & (th<th_end))+1;
 
for n=1:numel(th)
 
	set(p,'xdata',th(n),'ydata',y(n));
    set(p1,'xdata',[th(n) th(n)],'ydata',[-1 y(n)],'visible',vis{v(n)});
 
    drawnow 
end

[nant44]

merci Dut.