animation 3D sur matlab

Version imprimable

Bonjour,

j ai un soucis sur ce code

Code:

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
38
39
40
41
42
43
44
 
 %-----------------------------------------
% Animation
% d'aprés 'Apprendre et maitriser Matlab'
%-----------------------------------------
echo on;
%------------------------
%  Animation           
%------------------------
echo off;
clear all;
close all;
Q=490.916900919889;
K=1.63;
w=0.0982543280127309;
tp=600;
R=0;
Y0=[20];
 
 
t=1:10:600;
x=1:10:600;
[T,X]=meshgrid(t,x);
 
  A=(K*sqrt(pi)*Y0+2*Q*sqrt(w*t))/(K*sqrt(pi)*Y0+2*Q*sqrt(w*t)*sqrt(pi));
  B=log((K*sqrt(pi)*Y0+2*Q*sqrt(w*t)/((K*sqrt(pi)*Y0+2*Q*sqrt(w*t)*sqrt(pi))/(-w*t))));   
  z= A'*[(Y0'+(2*Q'*sqrt(w'*t')./K))'*cos(x./B)'-(2*Q)'./K'*B'*1/A'*sin(x./B)'];
mesh(T,X,Z); % filaire
 
lim=axis;% vecteur ligne
 
n=20;
M=moviein(n);%réserve la mémoire
 
for i=1:n
   mesh(X,Y,Z); %amortit la surface
   axis(lim);             % les axes ne changent pas
   M(:,i)=getframe;    %capture 1 image
end
if 0                  % pb en fin d'exécution si ce code est activé!
 loops=3;         % boucles
 ips=20;           % nbr images/s
 movie(M,loops,ips);
end

je souhait faire un film en 3D de la fonction z (temperature) en fonction de t (temps) et x ( profondeur)....
est ce que quelqu un pourrait m aider

merci.

26/03/2008, 12h42
mr_samurai

Je ne comprend pas ce que tu veux faire :aie: .

Si tu veux représenter Z = f(t,x) , c'est quoi la 3eme variable pour faire l'animation ?

Tu ne veux pas dire plutôt : Z=f(x) , et faire une animation pour t=1:tmax ?

Sinon , si une variable p est scalaire, p' perturbe uniquement le lecteur ;).

Tu peux essayer ceci :

Code:

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
38
39
 
%-----------------------------------------
% Animation
% d'aprés 'Apprendre et maitriser Matlab'
%-----------------------------------------
echo on;
%------------------------
%  Animation           
%------------------------
echo off;
clear all;
close all;
Q=490.916900919889;
K=1.63;
w=0.0982543280127309;
tp=600;
R=0;
Y0=[20];
 
t=1:10:600;
x=1:10:600;
[T,X]=meshgrid(t,x);
 
A=(K*sqrt(pi)*Y0+2*Q*sqrt(w*t))/(K*sqrt(pi)*Y0+2*Q*sqrt(w*t)*sqrt(pi));
B=log((K*sqrt(pi)*Y0+2*Q*sqrt(w*t)/((K*sqrt(pi)*Y0+2*Q*sqrt(w*t)*sqrt(pi))/(-w*t))));   
z= A'*[(Y0+(2*Q'*sqrt(w*t')./K))*cos(x./B) - (2*Q)./K*B'*1/A'*sin(x./B)];
 
 plot(x,z(1,:),'tag','xzPlot') 
 xlim([min(min(x))-1 max(max(x))+1])
 ylim([min(min(z))-1 max(max(z))+1])
 gui_data=guihandles(gcf);
 guidata(gcf,gui_data);
 
 for i=1:length(t)
      %set(gui_data.xzPlot,'xdata',x);
      set(gui_data.xzPlot,'ydata',z(i,:));
      guidata(gcf,gui_data);
      pause(0.05)
 end

26/03/2008, 12h53
cozo22

z est la temperature qui varie en fonction du temps t et de la profondeur x.

je veux donc faire une animation en 3D (axe t , axe x, et axe z) pour visualiser l evolution de la temperature en fonction de la profondeur x a chaque t (temps) donne.
26/03/2008, 13h02
cozo22

effectivement je viens de comprendre il me manque un axe ...

ce qu il se passe c est qu on a une evolution de la temperature en fonction du temps suivant deux directions x et y (ici x=y). il faut donc rajouter la direction y..
26/03/2008, 13h03
mr_samurai

Citation:

Envoyé par cozo22

z est la temperature qui varie en fonction du temps t et de la profondeur x.

je veux donc faire une animation en 3D (axe t , axe x, et axe z) pour visualiser l evolution de la temperature en fonction de la profondeur x a chaque t (temps) donne.

C'est bien ce que je pensais :nono: .

Une animation 3D prend en compte 4 paramétres : 3 pour le graphique 3D et le 4eme pour l'evolution .
Dans ton cas tu dispose de 3 paramètres et pas 4 .

Citation:

Envoyé par cozo22

visualiser l'evolution de la temperature en fonction de la profondeur x a chaque t

Autrement dit, tu trace z=ft(x) et tu fai varier t.

Si, tu ne vois pas ce que je veux dire, je peux rien pour toi :( .
26/03/2008, 13h05
cozo22

si si c est bon...

il fat juste que je retouche mon equation et la j en ai pour un petit moment....

merci bien