Programmation de l'équation de la chaleur
Résolu
Discrétisation et matlab
-
fuyehuc -
fuyehuc -
Bonjour,
je veux programmer l'équation de la chaleur discrétisée en coordonnées cylindriques mais j'ai des problémes aprés avoir déclaré mes données. je ne peut plus continuer.
Je mets ci joint cette équation et la programmation débutée.
T_i^(n-1)=-p.(1-1/(2.i)).T_(i-1)^n+(1+2.p).T_i^n-p.(1+1/(2.i)).T_(i+1)^n
B.T_(n+1)=T_n+A.S_(n+1)
clear;clf;clc;hold off;
k=80.2;ro=7870;Cp=447e3;TL=0;TR=200;
alpha=k/ro/Cp;dx=(50-1)/49;
%*******************************
% Construction de la matrice A *
%*******************************
A(1,1:2)=[-2 1];A(1,3:50)=0;
A(50,1:48)=0;A(50,49:50)=[1 -2];
for i=2:49
for j=1:50
if i<j-1 & j>i+1
A(i,j)=0;
end
if i==j
A(i,j)=-2;
A(i,j-1)=1;
A(i,j+1)=1;
end
end
end
je veux programmer l'équation de la chaleur discrétisée en coordonnées cylindriques mais j'ai des problémes aprés avoir déclaré mes données. je ne peut plus continuer.
Je mets ci joint cette équation et la programmation débutée.
T_i^(n-1)=-p.(1-1/(2.i)).T_(i-1)^n+(1+2.p).T_i^n-p.(1+1/(2.i)).T_(i+1)^n
B.T_(n+1)=T_n+A.S_(n+1)
clear;clf;clc;hold off;
k=80.2;ro=7870;Cp=447e3;TL=0;TR=200;
alpha=k/ro/Cp;dx=(50-1)/49;
%*******************************
% Construction de la matrice A *
%*******************************
A(1,1:2)=[-2 1];A(1,3:50)=0;
A(50,1:48)=0;A(50,49:50)=[1 -2];
for i=2:49
for j=1:50
if i<j-1 & j>i+1
A(i,j)=0;
end
if i==j
A(i,j)=-2;
A(i,j-1)=1;
A(i,j+1)=1;
end
end
end
Configuration: Windows XP Internet Explorer 6.0
A voir également:
- Programmation de l'équation de la chaleur
- Test chaleur pc - Guide
- Editeur d'équation - Télécharger - Vie quotidienne
- Application de programmation - Guide
- Télécharger le programme de pmu - Télécharger - Médias et Actualité
- Fermer un programme de force - Guide
j'ai un problème de programmation car je n'avait jamais itulisé ce logiciel
c'est exactement ce que je suis en train d'essayer de faire! T'en es où, est-ce que tu a résolu le problème. j'ai presque la même équation que toi.
Par ailleur dans ta première équation, ça ne serait pas Ti^(n+1) plutôt que (n-1)?