TP4 - Aides : 3. 4 - boucle en temps

TP4 - Aides : Etape 3 : schéma en temps : 3. 4 - boucle en temps

3. 4 - boucle en temps

Faire une boucle sur $\color{mycolor}{n}$ (temps) dans laquelle on résout le système linéaire et on stocke la solution $\color{mycolor}{(n+1)}$ dans la $\color{mycolor}{(n+2)}$ ième colonne du tableau $\color{mycolor}{U}$ :

dans cette boucle (de n=1 à nt-1), il faut : (commandes sous scilab en brun)

construire le second membre : $\color{brown}F=dt^2* Fc+M*(2*U(:,n)-U(:,n-1))-(1-\theta )*c^2*dt^2 *K*U(:,n-1) \; ;$
modifier le second membre pour tenir compte des CL de Dirichlet (FCLdi du TP1)
résoudre le système linéaire $A U^{n+1}=F$ : U(:,n+1)=lusolve(A,F)

Sous scilab : écrire en complétant les zones grises du fichier ondes1DP1.sci ,

fin de l’aide sur l’étape 3

$\includegraphics[width=8cm]{./prog-images/ondes1DP1trou.png}$