Suite d'opérations élémentaires sur les lignes (et éventuellement les colonnes) d'une matrice pour la transformer en forme échelonnée ou forme échelonnée réduite, utilisée pour résoudre des systèmes linéaires, calculer le rang, trouver des inverses et obtenir des bases des espaces de solutions par élimination systématique des variables.

Principe

Les opérations élémentaires sur les lignes (échanger deux lignes, multiplier une ligne par un scalaire non nul, ajouter un multiple scalaire d'une ligne à une autre) préservent l'ensemble des solutions du système linéaire correspondant ; les appliquer systématiquement supprime les inconnues étape par étape pour triangulariser la matrice des coefficients.

Démonstration

Résoudre un système 3×3 : constituer la matrice augmentée, utiliser des opérations sur les lignes pour obtenir une forme échelonnée supérieure, effectuer la substitution arrière pour obtenir la solution unique si les pivots sont non nuls ; en cas singulier la procédure révèle des dépendances et paramétrise l'espace des solutions.

Mauvaise application

Diviser par un pivot nul sans permuter d'abord les lignes, ou traiter négligemment les pivots en arithmétique flottante (omettant pivot partiel/total) peut provoquer des divisions par zéro ou une instabilité numérique sévère et des solutions numériques incorrectes.

Conséquence

Une application correcte donne une forme échelonnée (ou réduite) à partir de laquelle on lit le rang, on détermine la solvabilité, on calcule les inverses en rang plein et on obtient des descriptions paramétriques explicites des espaces de solutions ; elle sous-tend la factorisation LU et de nombreux solveurs directs linéaires.

Inversion

La conversion inverse consiste à effectuer des opérations sur les lignes non inversibles (mise à zéro par multiplication, transformations non réversibles) qui détruisent l'équivalence des solutions ; annuler l'élimination exige d'enregistrer les opérations sur les lignes ou de stocker des facteurs (LU) car la transformation n'est pas récupérable de façon unique à partir de la forme triangulaire seule.

Limite

S'applique aux systèmes linéaires sur des corps (ou anneaux divisionnaires avec précaution) ; sur des anneaux sans division (comme les entiers) l'application directe échoue sans adaptations (forme normale de Smith). En contexte numérique, la stabilité et la stratégie de pivotage sont des limites importantes à l'usage direct.

Tension sémantique

Tension avec les méthodes itératives : l'élimination est une procédure directe et finie adaptée aux systèmes denses de petite à moyenne taille, tandis que les solveurs itératifs (Gradient Conjugué, GMRES) conviennent mieux aux grands systèmes creux ou mal conditionnés pour des raisons de mémoire et complexité.

Synthèse

L'élimination de Gauss est la procédure par opérations élémentaires sur les lignes qui réduit un système linéaire en forme triangulaire (ou réduite), révélant le rang et la solvabilité et permettant une solution directe, l'inversion matricielle et la factorisation LU lorsque c'est applicable.