Si tu étais allé voir l'exemple trois en question, tu aurais vu que c'était exactement ton problème. Je te l'écris ici en ne considérant que l'équation 1 de ce que tu nous as donné (je considèrerai que T2 est une constante). Ce sera à toi de l'adapter pour tes 10 équations:
je suppose que T1 ne dépend que de x et que c1 ne dépend aussi que de x. Si c1 dépendait de y, la procédure serait de toute façon sensiblement la même.
c1.dT1/dx = U.W.(T2-T1) (1) T1(X=0)=Tchaud entrée (2)
Définition de la dépendance en x de c1:
1 2
| C1X=linspace(0,L,25) %génération de x pour C1
C1= C1X.^2 - 3 % génération de C1(x), à remplacer par ta propre équation. |
(Si tu ne connais pas l'équation traduisant les variations de ton C1, tu ne pourras pas résoudre grand chose.)
Ecriture d'un fichier .m pour interpoler les données de C1
1 2 3 4
| function dT1dx=myode(x,T1,C1X,C1)
C1=interpl(C1X,C1,x); interpole les données C1X, C1 à l'instant x
dT1dx = (U*W*(T2-T1))/C1 %cette equation est tirée de l'équation (1), il manque des paramètres d'entrée à la fonction. |
Resolution de l'équation
1 2 3 4
| Xspan=[0 L] %résolution pour x allant de 0 à L
CI=[Tchaud_entree] %condition initiale, équation (2)
[X T1]=ode45(@(x,t1) myode(x,t1,C1X,C1),Xspan,CI); |
Affichage
plot(X,T1) %reste à ajouter des titres et des légendes.
Partager