Discussion :

[cozo22]

Bonjour,


J'aimerais savoir si vous seriez regler mon petit probleme...

quand je lance ce code, il me donne l evolution de la courbe point par point mais j aimerais qu il me trace la courbe finale d evolution

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;
K=1.63;
Tp=122;
w=0.0982543280127309;
Q=490.91;
do=0.1;
Rho=0.0233;
Cp=712;
ro=2
n=1;
clear F
 
for r=0:1/40:4;
    t=50;
    l=quad(@(x)(exp(-((r-ro).^2)./(4*w*x'))./((pi*x').^(1/2)))',1,t);
    y=Q'*sqrt(w)./K'*l'
plot(r,y);
axis([0 4 0 1600]);
title('Surface temperature profile without conduction')
xlabel('x in mm')
ylabel('Temperature in Kelvin')
set(gcf, 'DoubleBuffer', 'on')
drawnow
pause(0.05);
F(n)=getframe;
n=n+1;
end

pour mieux comprendre lance le code tel qu il est ....
et moi je voudrais le trace de la courbe total plutot que de l avoir point par point...

[Caro-Line]

Au lieu de faire la commande plot dans la boucle il faut la sortir (avec ce qui va avec).
Et pour créer ton vecteur il faut indexer y avec r.
Par contre du coup tu fais quoi de ton getframe ?
sans le getframe :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
for r=0:1/40:4;
    t=50;
    l=quad(@(x)(exp(-((r-ro).^2)./(4*w*x'))./((pi*x').^(1/2)))',1,t);
    y(r)=Q'*sqrt(w)./K'*l'
end
r=0:1/40:4;
plot(r,y);
axis([0 4 0 1600]);
title('Surface temperature profile without conduction')
xlabel('x in mm')
ylabel('Temperature in Kelvin')

[cozo22]

maintenant quand je lance mon code ca fait une erreur

Code :

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

1
2
 
??? Subscript indices must either be real positive integers or logicals.

et pour repondre a ta question je n ai plus besoin du getfram

j ai 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
 
clear all;
close all;
K=1.63;
Tp=122;
w=0.0982543280127309;
Q=490.91;
do=0.1;
Rho=0.0233;
Cp=712;
ro=2
n=1;
clear F
for r=0:1/40:4;
    t=50;
    l=quad(@(x)(exp(-((r-ro).^2)./(4*w*x'))./((pi*x').^(1/2)))',1,t);
    y(r)=Q'*sqrt(w)./K'*l'
end
r=0:1/40:4;
plot(r,y);
axis([0 4 0 1600]);
title('Surface temperature profile without conduction')
xlabel('x in mm')
ylabel('Temperature in Kelvin')

[tug83]

c'est du à la solution de caro
y(0) ce n'est pas possible
en R2008a, vous auriez ce message clair:
??? Attempted to access y(0); index must be a positive integer or logical.

[Jerome Briot]

maintenant quand je lance mon code ca fait une erreur

Code :

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

1
2
 
??? Subscript indices must either be real positive integers or logicals.

Voir la : Subscript indices must either be real positive integers or logicals.

[Caro-Line]

Oups, OK j'ai été trop vite (pas la 1ère fois que ça arrive).
Mais c'est facile à résoudre.
Donc :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
r=0:1/40:4;
for k=1:length(r)
    t=50;
    l=quad(@(x)(exp(-((r(k)-ro).^2)./(4*w*x'))./((pi*x').^(1/2)))',1,t);
    y(k)=Q'*sqrt(w)./K'*l'
end
 
plot(r,y);
axis([0 4 0 1600]);
title('Surface temperature profile without conduction')
xlabel('x in mm')
ylabel('Temperature in Kelvin')

Et désolée je ne peux pas tester...

[mr_samurai]

Salut,

Matlab a raison, tu utilises y(r) alors que for r=0:1/40:4 .

Essayes 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
 
clear all;
close all;
K=1.63;
Tp=122;
w=0.0982543280127309;
Q=490.91;
do=0.1;
Rho=0.0233;
Cp=712;
ro=2;
n=1;
r=0:1/40:4;
 
for i=1:length(r)
    t=50;
    l=quad(@(x)(exp(-((r(i)-ro).^2)./(4*w*x'))./((pi*x').^(1/2)))',1,t);
    y(i)=Q*sqrt(w)./K*l';
end
 
plot(r,y);
axis([0 4 0 1600]);
title('Surface temperature profile without conduction')
xlabel('x in mm')
ylabel('Temperature in Kelvin')

+++ bonne chance,

PS :
transposer ne sers a rien quand ta variable est scalaire

[cozo22]

Ca c est une vrai reponse caro....merci beaucoup et faut pas etre gene par la petitte remarque de TUG...

a bientot.....

[cozo22]

merci aussi a samourai....