Complément de Schur - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Introduction

En algèbre linéaire et plus précisément en théorie des matrices, le complément de Schur est défini comme suit. Soit

M=\left[\begin{matrix} A & B \\ C & D \end{matrix}\right]

une matrice de dimension (p+q)×(p+q), où les blocs A, B, C, D sont des matrices de dimensions respectives p×p, p×q, q×p and q×q, avec D inversible. Alors, le complément de Schur du bloc D de la matrice M est constitué par la matrice de dimension p×p suivante :

ABD − 1C.

Lorsque B est la transposée de C, la matrice M est symétrique définie-positive si et seulement si D et son complément de Schur dans M le sont.

Le complément de Schur apparaît en particulier comme le résultat d'une élimination de Gauss «partielle» en multipliant la matrice M à droite avec la matrice «triangulaire inférieure» par blocs suivante

LT=\left[\begin{matrix} I_p & 0 \\ -D^{-1}C & D^{-1} \end{matrix}\right].

Ici, Ip désigne la matrice identité de dimension p×p. Après multiplication par la matrice LT, le complément de Schur apparaît dans le bloc p×p supérieur. La matrice produit est

M\cdot LT=\left[\begin{matrix} A-BD^{-1}C & BD^{-1} \\ 0 & I_q \end{matrix}\right].

L'inverse de M peut ainsi être exprimée en termes de D − 1 et de l'inverse du complément de Schur

 \left[ \begin{matrix} A & B \\ C & D \end{matrix}\right]^{-1} = \left[ \begin{matrix} \left(A-B D^{-1} C \right)^{-1}  &   -\left(A-B D^{-1} C \right)^{-1} B D^{-1} \\ -D^{-1}C\left(A-B D^{-1} C \right)^{-1} & D^{-1}+ D^{-1} C \left(A-B D^{-1} C \right)^{-1} B D^{-1} \end{matrix} \right],

ou encore plus simplement,

 \left[ \begin{matrix} A & B \\ C & D \end{matrix}\right]^{-1} =  \left[ \begin{matrix} I & 0 \\ -D^{-1}C & I \end{matrix}\right] \left[ \begin{matrix} (A-BD^{-1}C)^{-1} & 0 \\ 0 & D^{-1} \end{matrix}\right] \left[ \begin{matrix} I & -BD^{-1} \\ 0 & I \end{matrix}\right].

Application à la résolution d'équations linéaires

Le complément de Schur apparaît naturellement lors de la résolution d'un système d'équations linéaires de la forme

Ax + By = a
Cx + Dy = b

  • x et a sont des vecteurs colonne de dimension p,
  • y et b sont des vecteurs colonne de dimension q,
  • A, B, C, D sont comme précédemment.

En multipliant la seconde équation par BD − 1 puis en la soustrayant de la première, il vient

(A - BD^{-1} C) x = a - BD^{-1} b.\,

Ainsi, la résolution de cette équation en x est possible dès que D et son complément de Schur sont inversibles. Il est ensuite possible d'obtenir y en résolvant l'équation Cx + Dy = b. Cette méthode réduit le problème de l'inversion d'une matrice de dimension (p+q) \times (p+q) à celui de l'inversion de deux matrices de dimensions respectives p×p et q×q. En pratique, la matrice D doit être bien conditionnée pour rendre la méthode précise.

Bibliographie

  • Chapitre 7 du livre de R. A. Horn et C. R. Johnson, intitulé « Matrix Analysis », plusieurs fois édité par Cambridge University Press.
Page générée en 0.160 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise