Discussion :

[SEA123]

Bonjour,

j'ai un petit problème en ce qui concerne cette équation :
je ne sais pas quelle valeur prendre pour N

Sinon est ce que vous pensez que j'ai commis une erreur pour le programme:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
% Initialisation des paramètres
phi=1000;L=14.8;lambdaz=98.6;
t=0:1:300;
z=0;
Cp=200;
ro=14965;
alphaz=(lambdaz/(ro*Cp));
N=100;som=0;
for p=1:N
    som=som+(2*phi/(L*lambdaz))*cos(p*pi*z/L)/((p*pi/L)^2)*(1-exp(-alphaz*((p*pi/L)^2)*t))
end
T=som+((phi/(ro*Cp*L))*t);
figure;plot(t,T)

[Mat32]

Je pense qu'il y a une erreur dans la somme :

Code :

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

1
2
3
4
for p=1:N
    som=som+cos(p*pi*z/L)/((p*pi/L)^2)*(1-exp(-alphaz*((p*pi/L)^2)*t))
end
T=(2*phi/(L*lambdaz))+som+((phi/(ro*Cp*L))*t)

Par curiosité à quoi correspond cette formule ?

[SEA123]

Je pense qu'il y a une erreur dans la somme :

Code :

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

1
2
3
4
for p=1:N
    som=som+cos(p*pi*z/L)/((p*pi/L)^2)*(1-exp(-alphaz*((p*pi/L)^2)*t))
end
T=(2*phi/(L*lambdaz))+som+((phi/(ro*Cp*L))*t)

Par curiosité à quoi correspond cette formule ?

c'est une réponse thermique d'un échantillon,

il n y a pas d'addition entre les 2 termes de la formules donc pourquoi tu as écris:

Code :

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

T=(2*phi/(L*lambdaz))+som+((phi/(ro*Cp*L))*t)

[Mat32]

Désolé je voulais écrire :

Code :

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

1
2
3
4
for p=1:N
    som=som+cos(p*pi*z/L)/((p*pi/L)^2)*(1-exp(-alphaz*((p*pi/L)^2)*t))
end
T=(2*phi/(L*lambdaz))*som+((phi/(ro*Cp*L))*t)

A mon avis tu ne devrais pas mettre le code ci-dessous dans la boucle

Code :

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

(2*phi/(L*lambdaz))