Discussion :

[cozo22]

Bonjour,

j ai ce code qui fonctionne mais qui ne me donne pas ce que je veux.

au lieu de faire une animation qui fasse varier la fonction y (temperature ) en fonction de x (distance) j aimerais faire une animation afin que pour chaque t (temps) donne, il y ai un trace de la courbe y=f(x)...

Merci de votre aide

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;
Q=490.916900919889;
K=1.63;
w=0.0982543280127309;
tp=600;
R=0;
Y0=[20];
t=120;
x=0:1/300:0.4;
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))));
y=A*[(Y0+(2*Q*sqrt(w*t)/K))*cos(x/B)-(2*Q)/K*B*1/A*sin(x/B)];
 
plot(x,y)
axis([0 0.4 20 1500]);
title('temperature evolution on the surface of the susbtrate')
xlabel('x in cm')
ylabel('Temperature in degres celsius')
set(gcf, 'DoubleBuffer', 'on')
for i = 1 : length(x)
  p1 = plot(x(i), y(i), 'bo')
hold on
drawnow
end

[cozo22]

j ai trouve tout seul....

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;
Q=490.916900919889;
K=1.63;
w=0.0982543280127309;
tp=600;
R=0;
Y0=[20];
for t=1:1:122;
x=0:1/300:0.4;
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))));
y=A*[(Y0+(2*Q*sqrt(w*t)/K))*cos(x/B)-(2*Q)/K*B*1/A*sin(x/B)];
 
plot(x,y)
axis([0 0.4 20 1500]);
title('temperature evolution on the surface of the susbtrate')
xlabel('x in cm')
ylabel('Temperature in degres celsius')
set(gcf, 'DoubleBuffer', 'on')
%for i = 1 : length(x)
 % p1 = plot(x(i), y(i), 'bo')
%hold on
drawnow
%end
 
end