4. 3 -  Appels, tracés et test

Appeler les sous-programmes de résolution statique () et instationnaire () :

par les commandes (sous scilab) :

  Te=chaleur1DP1temps( nx,x,fk,fs,cldi,clne,rc,theta ,dt,nt,To)
  Ts=chaleur1DP1( nx,x,fk,cldi,clne,fs)

Tracer les courbes des au cours du temps et tracer les résultats en 3D () : (commandes sous scilab en brun) :

• Tracer au cours du temps :  
  on initialise la figure :  figure(1)  ;  clf()  ;
  on itère sur , toutes les 10 itérations :   for n = 1 : 10 : nt
  on trace en bleu (avec * et tirets) :  plot(x,Te(:,n), ’b*–’ )
  on impose un temps de pause :  sleep(5000)  ;
  on ferme la boucle :  end ;
  On peut de plus tracer la solution initiale en rouge et la solution statique en vert :  plot(x,To, ’r’,x,Ts, ’g’ ) ;
  et mettre titre et labels :
    title(’ theta=’+string(theta)...
       +’, dt=’+string(dt) +’, nb de points=’+string(nx) ) ;
    ylabel(’Temperature’) ;  xlabel(’abscisse le long du barreau’) ;

• tracer les résultats en 3D () :
  on initialise la figure :  figure(2)  ;  clf()  ;
  on crée le tableau des temps :   t = 0 : dt : (nt-1)*dt ;
  on précise la table des couleurs choisie :
    xset("colormap",jetcolormap(64)) ;
  on utilise la commande plot3D1 (voir le help de Scilab) :
   plot3d1(x,t,Te,110,45,"X@Y@Z",flag=[-2,2,4]) ; et on met labels et titre :
    xlabel(’abscisse sur le barreau’) ;
    ylabel(’temps’)  ;   zlabel(’Temperature’) ;
    title(’nb de points=’+string(nx)...
       +’, theta=’+string(theta)+’, dt=’+string(dt))
  

Test :  

Tester ce programme en prenant , et .
Après exécution du programme, taper dans la console Scilab; Vous obtenez les valeurs :

vous pouvez continuer, sinon vérifiez vos programmes.

Prenez maintenant et ;
Choisir différents et faire varier le pas de temps ; retrouver les propriétés de ce schéma (cf. cours 10)
Comparer la solution atteinte au bout d’une longue période avec la solution statique.

fin de l’aide sur l’étape 4