Sous-sections

6.3 Algorithme de GAUSS (directe)

6.3.1 Exemple

englishméthode d'élimination de Gauss pour un système 3*3 avec Maple

6.3.2 Principe

  1. transformation du système $A.x=b$
    en un système triangulaire $U.x=d$
  2. résolution de $U.x=d$

6.3.3 résolution système triangulaire


\begin{displaymath}
U.x=d  \mbox{{  avec }}   U_{ij}=0 \mbox{si}   i>j\end{displaymath}

Pour $i$ de $N$ à $1$

\begin{displaymath}
x_{i}=\frac{d_{i}-\sum_{j=i+1}^{n}U_{ij}*x_{j}}{U_{ii}}\end{displaymath}

Algorithme


\begin{algorithm}
% latex2html id marker 3108\caption{remont\'{e}
}
\par
\beg...
...i]$\leftarrow$(d{[}i]-somme)/U{[}i,i]
\par
finpour\end{list}\par
\end{algorithm}

6.3.4 Algorithme de GAUSS

Détail
de l'algorithme
Algorithme
 

\begin{algorithm}
% latex2html id marker 3204\caption{factorisation
}
\par
\b...
...$B{[}i]-Coef{*}B{[}k]
\par
~~finpour
\par
finpour~\end{list}\par
\end{algorithm}

coût
$\theta(N^{2})$ opérations

6.3.5 Problème du pivotage

2 stratégies si $A_{kk}^{k}=0$

6.3.5.1 Pivotage partiel

permutation de ligne à chaque étape k
remplace la ligne k par la ligne p telle que

\begin{displaymath}
\left\vert A_{pk}^{k}\right\vert=\max_{k\leq i\leq n}\left\vert A_{ik}^{k}\right\vert\end{displaymath}

6.3.5.2 Pivotage total

permutation
ligne p colonne q pour amener en $A_{kk}$ l'élément $A_{pq}$ telle que

\begin{displaymath}
\left\vert A_{pq}^{k}\right\vert=\max_{k\leq i;j\leq n}\left\vert A_{ij}^{k}\right\vert\end{displaymath}

remarque:
si un pivot est nul apres pivotage partiel ou total, alors la matrice $A$ est singulière


Pr. Marc BUFFAT
marc.buffat@univ-lyon1.fr
2008-02-28