Compte Rendu1 RT4 GR 3 Calcul Scientifique
Compte Rendu1 RT4 GR 3 Calcul Scientifique
Compte Rendu1 RT4 GR 3 Calcul Scientifique
2013/2014
Calcul Scientifique
Compte rendu TP N1
RT 4 Gr 3 Ralis par :
Omar Khayem Cherif
Mehdi Fajraoui
Calcul Scientifique
Compte Rendu N1
RT 4 Groupe 3
2
I ntroduction
La premire sance de Travaux Pratiques, a t consacre pour se familiariser la
console scilab et aux diffrentes commandes qui sont similaires celles utilises sur
MATLAB que nous nous sommes familiariss avec.
Les commandes excutes sont des commandes basiques qui permettent
deffectuer des oprations sur les matrices et de dterminer certaines caractristiques
comme le dterminant dune matrice [det(X)] ou linverse dune matrice [inv(X ] ou
encore les valeurs propres dune matrice [[V,D]=spec(X)] ainsi que certaines
messages derreurs et davertissements.
Introduction de Scilab:
Scilab: est un logiciel open source gratuit de calcul numrique qui
fournit un puissant environnement de dveloppement pour les
applications scientifiques et lingnierie. Scilab Enterprises est l'diteur et le
fournisseur de services professionnels de Scilab.
Cet outil prsente plusieurs fonctions mathmatiques et graphiques en [2D et 3D]
ainsi qun environnement programmation ce qui lui permet de fournir plusieurs
services dans diffrentes domaines tels que le domaine scientifique et de
lingnierie, comptable, de dveloppement etc
Toutes les informations propos cet outil sont disponible sur son site web :
http://www.scilab.org/
Travaux Pratiques
Travail raliser: Annulation des lments A
ij
de la matrice A en dessous du pivot
A
i-1 i-1.
Pour aboutir notre objectif, on procde comme suivant:
I . Premire tape
1- Dterminer les lignes de codes seulement pour annuler le premier
lment de la deuxime ligne de la matrice A
12
//Prenons une matrices A2
A=[5 4 ;1 1];
//Determinant le pivot
pivo=A(1,1);
//Determinant le coefficient d'annulation
coef=A(2,1)/pivo;
//Annulons le premiers lments sous le pivot
A=[A(1,:);A(2,:)-coef*A(1,:)]
Calcul Scientifique
Compte Rendu N1
RT 4 Groupe 3
3
2- Gnraliser le code pour tous les lments en dessous du premier Pivot
Pour ceci on prend une matrice A
3
.
Attention : llment traiter, doit tre diffrent de nul.
I I . Deuxime Etape :
1. Dtermination de tous les pivots et introductions des boucles pour obtenir une
matrice triangulaire :
//Prenons une matrices A3
A=[5 4 2; 1 1 1; 5 5 5];
//Determinant le pivot
pivo=A(1,1);
//Determinant les coefficients d'annulation
coef=[A(:,1)/pivo];
//Annulons les lments sous le pivot
for i=2:3,if A(i,1) <> 0 then
A(i,:)=A(i,:)-coef(i)*A(1,:),end, end
//Prenons une matrices A3
A=[5 4 2; 1 1 1; 5 5 5];
for j=1:3,
//Determinons les pivots
pivo=A(j,j);
//Determinant les coefficients d'annulation
coef=[A(:,j)/pivo];
//Annulons les lments sous le pivot
for i=j+1:3,if A(i,j) <> 0 then
A(i,:)=A(i,:)-coef(i)*A(j,:),end, end,
end
Calcul Scientifique
Compte Rendu N1
RT 4 Groupe 3
4
2. Test : Prenons la matrice A
5
Prenons la matrice:
A =
5. 4. 2. 2. 2.
1. 1. 1. 1. 1.
5. 5. 5. 10. 8.
- 1. - 2. 5. 8. 0.
0.5 6. 4. 8. - 1.
Aprs modification de code prcdent on obtient: