Bonjour,

j ai un soucis sur ce 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
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.