En mathématiques, et plus précisément en algèbre linéaire, la notion de pseudo-inverse (ou inverse généralisé) généralise celle d’inverse d’une application linéaire ou d’une matrice[1] aux cas non inversibles en lui supprimant certaines des propriétés demandées aux inverses, ou en l’étendant aux espaces non algébriques plus larges.

En général, il n’y a pas unicité du pseudo-inverse. Son existence, pour une application linéaire entre espaces de dimension éventuellement infinie, est équivalente à l'existence de supplémentaires du noyau et de l'image. Selon les propriétés demandées, le pseudo-inverse défini permet toutefois de généraliser la notion d'inverse en se restreignant au semi-groupe associatif multiplicatif seul, même s'il ne respecte pas les autres contraintes du corps ou de l'algèbre (en particulier les propriétés de distributivité ou de commutativité ne sont plus vraies dans le cas général, là où le véritable inverse peut les respecter).

Ont été étudiés en particulier les types de pseudo-inverses suivants :

  • le pseudo-inverse de Moore-Penrose dans le cas des matrices carrées non inversibles, mais généralisable à toute algèbre de matrices à valeurs dans un corps.
  • le pseudo-inverse de Drazin qui détermine la matrice qui constitue un point fixe dans la multiplication par l'exponentiation de matrices carrées au-delà d'un degré fini.
  • le pseudo-inverse à gauche et le pseudo-inverse à droite, utiles dans le cas des matrices non carrées qui ne sont jamais inversibles pour déterminer la factorisation en valeurs singulières, et qui ne sont pas nécessairement égaux non plus dans le cas de transformées non commutatives comme les opérateurs fonctionnels et distributions non discrètes.

Le pseudo-inverse se calcule à l’aide d’une généralisation du théorème spectral aux matrices non carrées.

Il est notamment utile dans le calcul de régressions (méthode des moindres carrés) pour un système d'équations linéaires.

Historique modifier

Pour une matrice à coefficients réels ou complexes (pas nécessairement carrée), ou pour une application linéaire entre espaces euclidiens ou hermitiens, en ajoutant certaines conditions supplémentaires, on a unicité du pseudo-inverse, appelé pseudo-inverse de Moore-Penrose (ou simplement « pseudo-inverse »), décrit par Eliakim Hastings Moore[2] dès 1920 et redécouvert indépendamment par Roger Penrose[3] en 1955. Erik Ivar Fredholm avait déjà introduit le concept de pseudo-inverse pour un opérateur intégral en 1903.

Cas général pour une application linéaire modifier

Définition et premières propriétés modifier

Soient   une application linéaire entre deux espaces vectoriels   et   et   une application linéaire de   dans  . Ces deux applications sont pseudo-inverses l'une de l'autre si les deux conditions suivantes sont satisfaites :

  et  .

Dans ce cas, les propriétés suivantes sont vérifiées :

  • l'espace   est la somme directe du noyau de   et de l'image de   ;
  • l'espace   est la somme directe du noyau de   et de l'image de   ;
  • les applications   et   induisent des isomorphismes réciproques entre leurs images ;
  • si l'application   est inversible, alors son inverse est l'application  .

Cette définition se traduit naturellement sous forme matricielle dans le cas d'espaces vectoriels de dimension finie.

Existence et construction modifier

Réciproquement, soit   une application linéaire entre deux espaces vectoriels   et  , dont le noyau admette un supplémentaire   dans   et dont l'image admette un supplémentaire   dans  . Alors la restriction de   à   induit un isomorphisme entre   et son image. L'application réciproque de l'image de   vers   s'étend de façon unique par l'application nulle sur  , en une application linéaire   de   dans   qui est par construction pseudo-inverse de  .

Il y a donc correspondance biunivoque entre les pseudo-inverses d'une application linéaire et les couples de supplémentaires pour son noyau et son image.

Remarque : ceci s'applique évidemment aux cas où l'un des supplémentaires   et   est réduit à l'origine ou à l'espace vectoriel tout entier, ce qui a lieu en particulier lorsque   est inversible :   est alors égal à   et   est réduit à l'origine.

Choix des supplémentaires modifier

Il n'y a pas de choix canonique d'un supplémentaire en général, mais une structure d'espace euclidien ou hermitien sur les espaces vectoriels source et but permet d'en déterminer un par la définition de l'orthogonal. Cette définition du pseudo-inverse correspond au « pseudo-inverse de Moore-Penrose » pour les matrices.

Cas matriciel modifier

Définition modifier

Étant donné une matrice   à coefficients réels ou complexes avec   lignes et   colonnes, son pseudo-inverse   est l'unique matrice à   lignes et   colonnes vérifiant les conditions suivantes :

  1.   ;
  2.         (  est un inverse pour le semi-groupe multiplicatif) ;
  3.         (  est une matrice hermitienne) ;
  4.         (  est également hermitienne).

Ici, la notation   désigne la matrice adjointe à  , donc la transposée pour le cas réel.

Cette matrice peut s'obtenir comme une limite :

 

qui existe même si les matrices produits ( ) et ( ) ne sont pas inversibles.

Propriétés modifier

 

Identités valables pour toute matrice  
(à coefficients réels ou complexes)

L’opération de pseudo-inversion :

  • est involutive ;
  • commute avec la transposition et la conjugaison ;
  • est un antimorphisme sur le produit, sous certaines hypothèses : soient   et   deux matrices dont le produit   existe. Si l’une au moins est unitaire, ou si les deux matrices sont de rang égal à leur dimension commune, alors   ;
  • n’est pas continue en 0 (la matrice nulle). En effet, elle est homogène de degré –1 : pour tout scalaire   non nul,  . Plus généralement, elle n’est continue qu’au voisinage des matrices inversibles (ou de rang maximum pour des matrices non carrées[4]).

Projections modifier

Dans le cas matriciel,   et   sont des projecteurs orthogonaux, c'est-à-dire des matrices hermitiennes ( ,  ) et idempotentes (  et  ), et l'on a les résultats suivants[5] :

  •   et   ;
  •   est le projecteur orthogonal sur l'image de   (égale à l'orthogonal du noyau de  ) ;
  •   est le projecteur orthogonal sur l'image de   (égal à l'orthogonal du noyau de  ) ;
  •   est le projecteur orthogonal sur le noyau de   ;
  •   est le projecteur orthogonal sur le noyau de  .

Calcul effectif modifier

Si la matrice  , avec   lignes et   colonnes, est de rang  , alors elle peut s'écrire comme un produit de matrices de même rang  , où   possède   lignes et   colonnes et   possède   lignes et   colonnes. Dans ce cas les produits ( ) et ( ) sont inversibles et la relation suivante est vérifiée :

 .

Des approches optimisées existent pour le calcul de pseudo-inverses de matrices par blocs.

Algorithmiquement, le pseudo-inverse s'obtient à partir de la décomposition en valeurs singulières : muni de cette décomposition  , on calcule

 ,

 est le pseudo-inverse de  . La matrice   est constituée de deux blocs : une matrice diagonale et une matrice nulle. Son pseudo-inverse est une matrice dont les éléments non nuls sont obtenus en inversant les éléments non nuls (de la diagonale) de  , et en prenant le transposé de la matrice ainsi obtenue.

À partir d'une matrice dont le pseudo-inverse est connu, il existe des algorithmes spécialisés[Lesquels ?] qui effectuent le calcul plus rapidement pour des matrices en rapport avec la première. En particulier, si la différence n'est que d'une ligne ou colonne changée, supprimée ou ajoutée, des algorithmes itératifs peuvent exploiter cette relation.

Cas particuliers modifier

 

Pseudo-inverse d'un vecteur colonne

  • Le pseudo-inverse d’une matrice nulle est sa transposée (également nulle).
  • Le pseudo-inverse d’un vecteur colonne non nul est son vecteur adjoint divisé par sa norme au carré.
  • Plus généralement, si le rang de   est égal à son nombre de lignes[1], la matrice   peut être choisie égale à l’identité et dans ce cas :
     .
  • De même, si le rang de   est égal à son nombre de colonnes,
     .
  • A fortiori si la matrice   est inversible, son pseudo-inverse est son inverse.
  • Si le pseudo-inverse de   est connu, on peut en déduire   par l’égalité :
      ;
  • De même, si   est connu, le pseudo-inverse de   est donné par :
     .

Exemple d’utilisation modifier

Le pseudo-inverse donne une solution approchée   à un système   d’équations linéaires : celle que donnerait la méthode des moindres carrés[6].

Le système   a des solutions si et seulement si  , et les solutions sont alors[7] tous les vecteurs de la forme

 , où le vecteur   est arbitraire (si ce n'est sa dimension).

En effet (voir supra), l'image de   est égale au noyau de  , et le noyau de   est égal à l'image de  .

Notes et références modifier

(en)/(de) Cet article est partiellement ou en totalité issu des articles intitulés en anglais « Moore–Penrose inverse » (voir la liste des auteurs) et en allemand « Pseudoinverse » (voir la liste des auteurs).
  1. a et b (en) Adi Ben-Israel et Thomas N. E. Greville, Generalized Inverses : Theory and Applications, Springer-Verlag, , 2e éd. (1re éd. 1974) (ISBN 0-387-00293-6)
  2. (en) Eliakim Hastings Moore, « On the reciprocal of the general algebraic matrix », Bull. AMS, vol. 26,‎ , p. 394-395 (lire en ligne, consulté le )
  3. (en) Roger Penrose, « A generalized inverse for matrices », Proceedings of the Cambridge Philosophical Society, vol. 51,‎ , p. 406-413
  4. (en) Vladimir Rakočević, « On continuity of the Moore–Penrose and Drazin inverses », Matematički Vesnik, vol. 49,‎ , p. 163-72 (lire en ligne).
  5. (en) Gene H. Golub et Charles F. Van Loan, Matrix Computations, Baltimore/London, Johns Hopkins, , 3e éd. (ISBN 0-8018-5414-8), p. 257-258.
  6. (en) Roger Penrose, « On best approximate solution of linear matrix equations », Proceedings of the Cambridge Philosophical Society, vol. 52,‎ , p. 17-19.
  7. (en) M. James, « The generalised inverse », Mathematical Gazette, vol. 62, no 420,‎ , p. 109-114 (DOI 10.1017/S0025557200086460).

Voir aussi modifier

Article connexe modifier

Régularisation Tychonoff

Liens externes modifier

Bibliographie modifier

  • (de) W. Mackens et H. Voß , « Mathematik I für Studierende der Ingenieurwissenschaften »
  • (de) A. Kielbasinski et H.Schwetlick, « Numerische lineare Algebra », Deutscher Verlag der Wissenschaften, 1988