Boucle For

[invite032bc05b]

Bonjour,

Je souhaite réaliser une itération d'un calcul, mais je n'y parviens pas.

Voici l'exemple de calcul que je voudrai réalisé en boucle for

clear all
clc
x=0.01;
e=0.01;
Esoleil=400;
alpha=0.5;
epsi=0.26;
Tint=0;
Text=0;
Tcarro=(Tint+Text)/2;
Te=273;
hint=4;
hext=4;
k=46;
sigma=5.67E-8;

phiint=alpha*epsi*sigma*Te^4;
phiemis=epsi*sigma*Tcarro^4;
phis=alpha*Esoleil;
B=(hext*Text+phiint+phis-2*phiemis-hint*((e/k)*(-phis-hext*Text+phiemis)-Tint))*(1/(hext+(hint*hext*e/k)+hint));
A=(-phis-hext*(Text-B)+phiemis)/k;

T1=A*x+B

Tcarro=T1+273.15
phiint=alpha*epsi*sigma*Te^4;
phiemis=epsi*sigma*Tcarro^4;
phis=alpha*Esoleil;
B=(hext*Text+phiint+phis-2*phiemis-hint*((e/k)*(-phis-hext*Text+phiemis)-Tint))*(1/(hext+(hint*hext*e/k)+hint));
A=(-phis-hext*(Text-B)+phiemis)/k;

T2=A*x+B

Tcarro=T2+273.15
phiint=alpha*epsi*sigma*Te^4;
phiemis=epsi*sigma*Tcarro^4;
phis=alpha*Esoleil;
B=(hext*Text+phiint+phis-2*phiemis-hint*((e/k)*(-phis-hext*Text+phiemis)-Tint))*(1/(hext+(hint*hext*e/k)+hint));
A=(-phis-hext*(Text-B)+phiemis)/k;

T3=A*x+B

Tcarro=T3+273.15
phiint=alpha*epsi*sigma*Te^4;
phiemis=epsi*sigma*Tcarro^4;
phis=alpha*Esoleil;
B=(hext*Text+phiint+phis-2*phiemis-hint*((e/k)*(-phis-hext*Text+phiemis)-Tint))*(1/(hext+(hint*hext*e/k)+hint));
A=(-phis-hext*(Text-B)+phiemis)/k;

T4=A*x+B

Tcarro=T4+273.15
phiint=alpha*epsi*sigma*Te^4;
phiemis=epsi*sigma*Tcarro^4;
phis=alpha*Esoleil;
B=(hext*Text+phiint+phis-2*phiemis-hint*((e/k)*(-phis-hext*Text+phiemis)-Tint))*(1/(hext+(hint*hext*e/k)+hint));
A=(-phis-hext*(Text-B)+phiemis)/k;

T5=A*x+B



J'ai essayé de faire une boucle for comme suit, mais scilab me note que la fonction T n est pas définit. Or je veux seulement faire varier le nom de la variable avec n.
for n=1:10
T(0)=Tcarro;
T(n-1)=T(n-1);
Tcarro=T(n-1);
phiint=alpha*epsi*sigma*Te^4;
phiemis=epsi*sigma*Tcarro^4;
phis=alpha*Esoleil;
B=(hext*Text+phiint+phis-2*phiemis-hint*((e/k)*(-phis-hext*Text+phiemis)-Tint))*(1/(hext+(hint*hext*e/k)+hint));
A=(-phis-hext*(Text-B)+phiemis)/k;"

T(n)=A*x+B;
end

disp ("Température parois soumise au rayonnement solaire: "+string(T5)+"°C")


Merci de votre aide
-----

[Jack]

Je constate qu'une bonne lecture de nos règle participatives est nécessaire!
http://forums.futura-sciences.com/pr...eau-forum.html

[umfred]

ton code fonctionne bien (en tout cas, ne génère pas d'erreurs) si:
- tu retires le " sur la ligne de calcul de A
- tu remplaces T5 dans le disp par T(5)
- tu pars de T(1) et donc n va de 2 à 11 (pour conserver ton nombre d'itérations)

[invite032bc05b]

Merci bien, cela fonctionne

[A voir en vidéo sur Futura]

[danyvio]

Je ne connais pas spécifiquement scilab, mais je me demande ce que l'instruction suivante puisée dans la boucle for de la demande de ce post donne :

T(n-1)=T(n-1);