Voir la source de Règle de Cramer

{{Voir homonymes|Cramer}}

La '''règle de Cramer''' (ou '''méthode de Cramer''') est un [[théorème]] en [[algèbre linéaire]] qui donne la solution d'un '''système de Cramer''', c'est-à-dire un [[système d'équations linéaires]] avec autant d'équations que d'inconnues et dont le [[déterminant (mathématiques)|déterminant]] de la matrice de coefficients est non nul, sous forme de quotients de déterminants.

En calcul, la méthode est moins efficace que la [[élimination de Gauss-Jordan|méthode de résolution de Gauss]] pour des grands systèmes (à partir de quatre équations) dont les coefficients dans le premier membre sont explicitement donnés. Cependant, elle est d'importance théorique, car elle donne une expression explicite pour la solution du système, et elle s'applique dans des systèmes où par exemple les coefficients du premier membre dépendent de paramètres, ce qui peut rendre la méthode de Gauss inapplicable.

Elle est nommée d'après le [[mathématicien]] [[République de Genève|genevois]] [[Gabriel Cramer]] (1704-1752).

== Description ==

Le [[système d'équations linéaires|système de ''n'' équations à ''n'' inconnues]], de forme générale :
:<math>\left\{\begin{matrix} 
a_{1,1}x_1+a_{1,2}x_2+...+a_{1,n}x_n = \lambda_1 \\ 
a_{2,1}x_1+a_{2,2}x_2+...+a_{2,n}x_n = \lambda_2 \\ 
\vdots \\
a_{n,1}x_{1}+a_{n,2}x_{2}+...+a_{n,n}x_n = \lambda_n
\end{matrix}\right.</math>
est représenté sous la forme d'un [[produit matriciel]] :
:<math>\begin{pmatrix} 
a_{1,1} & a_{1,2} & \cdots & a_{1,n}\\ 
a_{2,1} & a_{2,2} & \cdots & a_{2,n}\\ 
\vdots & \vdots & \ddots & \vdots\\ 
a_{n,1} & a_{n,2} & \cdots & a_{n,n}\\ 
\end{pmatrix} \times 
\begin{pmatrix} 
x_1\\ 
x_2\\ 
\vdots\\ 
x_n\\ 
\end{pmatrix} = \begin{pmatrix} 
\lambda_1\\ 
\lambda_2\\ 
\vdots\\ 
\lambda_n\\ 
\end{pmatrix} \Leftrightarrow A \cdot X=\Lambda</math>

où la [[Matrice (mathématiques)|matrice]] <math>A</math>, carrée, contient les coefficients des inconnues, le [[vecteur colonne]] <math>X</math> contient ces inconnues et le vecteur colonne <math>\Lambda</math> contient les membres de droite des équations du système ; les coefficients et les inconnues font partie d'un même [[corps commutatif]].

Le théorème affirme alors que le système admet une unique solution [[Équivalence logique|si et seulement si]] sa matrice <math>A</math> est [[matrice inversible|inversible]] (déterminant non nul), et cette solution est alors donnée par :

:<math>x_k = { \det(A_k) \over \det(A) }</math>

où <math>A_k</math> est la matrice carrée formée en remplaçant la ''k''-ième colonne de <math>A</math> par le vecteur colonne <math>\Lambda</math>.
:<math>A_k = ( a_{k|i,j} ) \mbox{ avec } a_{k|i,j} = \left\{\begin{matrix} a_{i,j} & \mbox{si } j \ne k \\ \lambda_{i} & \mbox{si }j = k\end{matrix}\right.</math>

Un système carré (i.e. avec autant d'équations que d'inconnues) est dit de Cramer si le déterminant de sa matrice est non nul.

Lorsque le système (toujours carré) n'est pas de Cramer (i.e. lorsque le déterminant de ''A'' est nul) :
* si le déterminant d'une des matrices <math>A_k</math> est non nul, alors le système n'a pas de solution ;
* la réciproque est fausse : il peut arriver que le système n'ait pas de solution bien que les déterminants <math>\det(A_k)</math> soient tous nuls. Un exemple en est donné par :
:<math>\left\{\begin{matrix}x+y+z=1\\
x+y+z=2\\
x+y+z=3\end{matrix}\right.</math>
Pour plus de précisions, voir [[Théorème de Rouché-Fontené]].

Le nombre d'opérations à effectuer pour résoudre un système linéaire à l'aide de la règle de Cramer dépend de la méthode utilisée pour calculer le déterminant. Une méthode efficace pour les calculs de déterminant est l'élimination de Gauss-Jordan ([[Théorie de la complexité (informatique théorique)|complexité polynomiale]]). Cependant, la règle de Cramer demandera d'avoir recours à un nombre de calculs de déterminants égal à la taille du système, une élimination de Gauss-Jordan appliquée directement au système résout donc le problème plus efficacement.

== Démonstrations==

Si ''A'' est inversible, calculons la solution ''X'' (dont on sait qu'elle existe et est unique).

;Méthode directe<ref>{{Ouvrage|auteur1=J.-P. Marco|auteur2=L. Lazzarini|responsabilité2=dir.|titre=Mathématiques L1|sous-titre=Cours complet avec fiches de révision, {{Unité|1000|tests}} et exercices corrigés|éditeur=[[Pearson (maison d'édition)|Pearson]]|année=2013|numéro d'édition=2|passage=479|isbn=|lire en ligne=https://books.google.fr/books?id=r1bzeQoVx_MC&pg=PA479}}.</ref>
:Notons ''C''{{ind|1}}, … , ''C{{ind|n}}'' les colonnes de ''A''. L'égalité ''AX'' = Λ se réécrit{{Retrait|<math>\Lambda=x_1C_1+\ldots+x_nC_n.</math>}}Par linéarité du déterminant suivant la ''k''-ième colonne, on en déduit{{Retrait|<math>\det(A_k)=x_1\det(B_{k,1})+\cdots+x_n\det(B_{k,n})</math>}}où ''B{{ind|k,j}}'' désigne la matrice ''A'' dans laquelle la ''k''-ième colonne est remplacée par ''C{{ind|j}}''. Or pour tout ''j'' ≠ ''k'', la matrice ''B{{ind|k,j}}'' a deux colonnes égales donc son déterminant est nul. Il reste alors{{Retrait|<math>\det(A_k)=x_k\det(B_{k,k})=x_k\det(A),</math>}}d'où le résultat, en divisant par det(''A'') qui par hypothèse est non nul.
;Méthode utilisant la [[Comatrice#Formules de Laplace|formule de Laplace]]<ref>{{Ouvrage|auteur1=[[Jean-Pierre Ramis]]|auteur2=[[André Warusfel]]|responsabilité2=dir.|titre=Mathématiques Tout-en-un pour la Licence|sous-titre=Niveau 1|éditeur=[[Éditions Dunod|Dunod]]|année=2013|numéro d'édition=2|passage=382|isbn=|lire en ligne=https://books.google.fr/books?id=K3gcey1WUXMC&pg=PA382}}.</ref>{{,}}<ref>L. Thomas, [http://imbsrv1.epfl.ch/csag/cours/alglin0910/PolycopAlgLin0910.pdf Algèbre Linéaire Bachelor {{1re|année}} 2009-2010] (polycopié initial élaboré par E. B. Fluckiger et P. Chabloz), [[EPFL]], septembre 2009, {{p.|47-48}}.</ref>
:La [[matrice inverse]] ''A''{{-1}} est donnée par la formule de Laplace{{Retrait|<math>A^{-1}=\frac1{\det A} \, {}^t{{\rm com} A}</math>}}où <math>{}^t{{\rm com} A}</math> est la [[Matrice transposée|transposée]] de la [[comatrice]] de ''A''. Exprimons chacune des coordonnées de l'unique solution ''X = A''{{-1}}Λ, pour ''k'' variant de 1 à ''n'' :{{Retrait|<math>x_k=\frac{\sum_{j=1}^{n} \lambda_j{\rm Cof}_{j,k}}{\det(A)}.</math>}}À nouveau d'après la formule de Laplace, la somme au numérateur est le développement de <math>\det(A_k) </math> par rapport à sa ''k''-ième colonne donc{{Retrait|<math>x_k={ \det(A_k)\over \det(A) }.</math>}}
;Remarque
{{Refsou|date=juin 2015|La formule Cramer permet, inversement, de démontrer celle de Laplace.}}

== Exemples ==
=== Système d'ordre 2 ===
Si <math>ad-bc\ne0</math>, le système
:<math>\left\{\begin{matrix}
ax+by = e\\
cx+dy = f\end{matrix}\right.</math>
a pour unique solution :
:<math>x = { \begin{vmatrix}e&b\\f&d\end{vmatrix} \over \begin{vmatrix}a&b\\c&d\end{vmatrix} } = { ed - bf \over ad - bc},\quad y = { \begin{vmatrix}a&e\\c&f\end{vmatrix} \over \begin{vmatrix}a&b\\c&d\end{vmatrix} } =  { af - ec \over ad - bc}.</math>

Exemple numérique :
:<math>\left.\begin{matrix}
4x + 2y = 24\\
2x + 3y = 16\end{matrix}\right\}\Leftrightarrow\left\{\begin{matrix}x = {24 \cdot 3-16 \cdot 2 \over 8 } = {40 \over 8} = 5\\y = {4 \cdot 16-2 \cdot 24 \over 8} = {16 \over 8} = 2\end{matrix}\right.</math>

=== Système d'ordre 3 ===

:<math>\left\{\begin{matrix}a_1x_1 + b_1x_2 + c_1x_3 = d_1\\
a_2x_1 + b_2x_2 + c_2x_3 = d_2\\
a_3x_1 + b_3x_2 + c_3x_3 = d_3\end{matrix}\right.</math>

Posons :
:<math>A = \begin{pmatrix}a_1&b_1&c_1\\a_2&b_2&c_2\\a_3&b_3&c_3\end{pmatrix},\quad X= \begin{pmatrix}x_1\\x_2\\x_3\end{pmatrix}\quad\text{et}\quad
\Lambda = \begin{pmatrix} 
d_1\\ 
d_2\\ 
d_3
\end{pmatrix}.</math>

Le système admet une solution unique si et seulement si <math>\det(A) \ne 0</math> :
:<math>x_1 = \frac{\det(A_1)}{\det(A)} = \frac{\begin{vmatrix}d_1&b_1&c_1\\d_2&b_2&c_2\\d_3&b_3&c_3\end{vmatrix}}{\det(A)}</math>
:<math>x_2 = \frac{\det(A_2)}{\det(A)} = \frac{\begin{vmatrix}a_1&d_1&c_1\\a_2&d_2&c_2\\a_3&d_3&c_3\end{vmatrix}}{\det(A)}</math>
:<math>x_3 = \frac{\det(A_3)}{\det(A)} = \frac{\begin{vmatrix}a_1&b_1&d_1\\a_2&b_2&d_2\\a_3&b_3&d_3\end{vmatrix}}{\det(A)}</math>

Ou plus simplement :
:<math>X=\begin{pmatrix}x_1\\x_2\\x_3\end{pmatrix} = \frac1{\det(A)} \cdot \begin{pmatrix}
\det(A_1)\\
\det(A_2)\\
\det(A_3)\end{pmatrix}.</math>

====Cas du déterminant nul====
Pour que le système n'admette aucune solution, il ''[[condition suffisante|suffit]]'' que :
:<math>\det(A) = 0\quad\text{et}\quad\Big( \det(A_1) \ne 0\quad\text{ou}\quad\det(A_2) \ne 0\quad\text{ou}\quad\det(A_3) \ne 0 \Big)\,</math>
En effet, cela signifie que les 3 lignes du système sont linéairement dépendants lorsqu'on ne considère que le membre de gauche, mais elle ne le sont plus lorsqu'on inclut le membre de droite. Il ne peut donc y avoir de solution.

Dans le cas 
:<math>\det(A) = \det(A_1) = \det(A_2) = \det(A_3) = 0\,</math> 
on peut avoir soit une infinité de solutions, soit aucune:
Comme remarqué précédemment, det(''A'') = 0 signifie que les lignes de coefficients du système sont linéairement dépendants, l'une est [[combinaison linéaire]] des autres. Il faut dans ce cas extraire un système de Cramer, en choisissant ''r'' = rang ''A'' inconnues principales et r equations qui donnent un système de Cramer lorsqu'on met les inconnues non principales dans le membre de droite, et donc une solution "unique" exprimant les inconnues principales en termes des autres inconnues qui restent paramètres libres. (Cela correspond à un [[Espace affine|sous-espace affine]] de dimension ''n'' - ''r''.) Enfin, on vérifie si les équations supplémentaires sont satisfaites pour la solution obtenue ainsi. Dans le cas contraire, il n'y a aucune solution.

<!--hors sujet en l'état
=== Système d'ordre 4 ===
Le calcul du déterminant du système d'ordre 4 est basé sur la méthode du système d'ordre 3, utilisant les coordonnées du 1er vecteur vertical comme coefficients de signe alterné.
Soit :
:<math>\left\{\begin{matrix}a_1x_1 + b_1x_2 + c_1x_3 + d_1x_4 = e_1\\
a_2x_1 + b_2x_2 + c_2x_3 + d_2x_4 = e_2\\
a_3x_1 + b_3x_2 + c_3x_3 + d_3x_4 = e_3\\
a_4x_1 + b_4x_2 + c_4x_3 + d_4x_4= e_4\end{matrix}\right.</math>

Posons :
:<math>A = \begin{pmatrix}a_1&b_1&c_1&d_1\\a_2&b_2&c_2&d_2\\a_3&b_3&c_3&d_3\\a_4&b_4&c_4&d_4\end{pmatrix} \mbox{ et }X=\begin{pmatrix}x_1\\x_2\\x_3\\x_4\end{pmatrix}</math>.

:<math>\det(A) = a_1\det(A_1)-a_2\det(A_2)+a_3\det(A_3)-a_4\det(A_4)\,</math>
-->

== Références==
{{Références|colonnes=1}}

{{Portail|algèbre}}

[[Catégorie:Déterminant]]
[[Catégorie:Théorème d'algèbre|Cramer]]