Discussion :

[Snurgz]

Bonjour à tous!
(je m'excuse déjà pour le pavé que ca va être ) je suis totalement débutant sur MatLab et j'ai un exercice que j'essaye de finir depuis pas mal d'heures déjà^^.

Je dois donc tracer deux courbes(Y1,Y2) représentant le "gain" de chaleur dans la terre en fonction de la profondeur.

Suite à des calculs on en déduit que la Temperature est liée à la profondeur selon cette équation : T(z) = -(A/2k)*z^2 +((Q+Ad)/k)*z)

ou A vaut A1 = 1.25 x 10-6 W.m-3
A2 = 2.50 x 10-6 W.m-3 => Les 2 courbes se font en fonction de ces deux paramètres.

k = 2.5 W.m-1
d = 50 km
Q = 0.021 W.m-2
z = 0 à 50 km.

J'ai donc écrit mon code sur MatLab :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
A1=1.25*10^-6;
A2=2.5*10^-6;
k=2.5;
Qd=0.021;
d=50000;
z=0:1000:50000;
y1=(-(A1/2*k))*z.^2 + ((Qd+A1*d)/k)*z;
y2=(-(A2/2*k))*z.^2 + ((Qd+A2*d)/k)*z;
plot(z,y1,'b',z,y2,'g');
legend('A1=1.25*10^-6 w.m^-1.c^1','A2=02.5*10^-6 w.m^-1.c^1')
title('Temp?rature en fonction de la profondeur')
xlabel('Profondeur en m');
ylabel('Temp?rature en Celsius');
legend('A1=1.25*10^-6 w.m^-1.c^1','A2=02.5*10^-6 w.m^-1.c^1').

Ce qui me donne ceci:



Et le petit problème que j'ai c'est que de
1) J'ai une perte de chaleur alors que je descend ce qui devrait être le contraire^^
et de 2 ) Si j'enleve le signe moins de l'axe Y il me semble bizarre qu'a seulement 50 km il fasse deja 5000 degré celsius( Donc plus chaud que dans le noyau de la terre^^)

Voilà si qqn a une petite idée de la faute que j'ai pu commettre c'est volontiers
Merci bye

[Jerome Briot]

Un indice :

Code :

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

1
2
3
4
5
>> 1/2*3
 
ans =
 
    1.5000

[Snurgz]

salut Dut, merci de ton indice

Je crois que j'ai compris ton indice , lorsqu'on écrit 1/2*3 on s'attend à un 1/6 alors que la machine le comprend comme étant 1/2 +1/2+1/2= 1.5?