Moore–Penroses pseudoinvers

Moore–Penroses pseudoinvers är inom linjär algebra en generalisering av vissa egenskaper hos matrisinversen för icke-kvadratiska matriser, uppkallad efter Eliakim Hastings Moore och Roger Penrose, som beskrev den oberoende av varandra 1920 respektive 1955.

Definition

Moore–Penroses pseudoinvers till en matris A {\displaystyle A} är en matris A + {\displaystyle A^{+}} som uppfyller:

  1. A A + A = A {\displaystyle AA^{+}A=A}       ( A A + {\displaystyle AA^{+}} behöver inte vara en enhetsmatris, men ska avbilda alla kolonnvektorer i A {\displaystyle A} på sig själva);
  2. A + A A + = A + {\displaystyle A^{+}AA^{+}=A^{+}}       ( A + {\displaystyle A^{+}} is är en svag invers för den multiplikativa semigruppen);
  3. ( A A + ) = A A + {\displaystyle (AA^{+})^{*}=AA^{+}}       ( A A + {\displaystyle AA^{+}} är en hermitesk matris)
  4. ( A + A ) = A + A {\displaystyle (A^{+}A)^{*}=A^{+}A}       ( A + A {\displaystyle A^{+}A} är också hermitesk).

A {\displaystyle A^{*}} är det hermiteska konjugatet till A {\displaystyle A} . För reella matriser är detta samma sak som transponatet.

Egenskaper

Givet en matris A {\displaystyle A} med Moore–Penroses pseudoinvers A + {\displaystyle A^{+}} , gäller följande:

  • A + {\displaystyle A^{+}} är unik.
  • Om A {\displaystyle A} är en inverterbar matris, är A 1 = A + {\displaystyle A^{-1}=A^{+}} .
  • Pseudoinversen av pseudoinversen är den ursprungliga matrisen, ( A + ) + = A {\displaystyle (A^{+})^{+}=A} .
  • A A + {\displaystyle AA^{+}} är en ortogonal projektion på A {\displaystyle A} s värderum.
  • A + A {\displaystyle A^{+}A} är en ortogonal projektion på A {\displaystyle A^{*}} s värderum.
  • Pseudoinversen till en nollmatris är dess transponat.

Specialfall

Ortonormala rader och kolonner

Om A {\displaystyle A} har ortonormala kolonnvektorer ( A A = I {\displaystyle AA^{*}=I} ) eller ortonormala radvektorer ( A A = I {\displaystyle A^{*}A=I} så är A + = A {\displaystyle A^{+}=A^{*}} ).

Linjärt oberoende kolonner och rader

Om kolonnerna i A {\displaystyle A} är linjärt oberoende är A A {\displaystyle A^{*}A} inverterbar och Moore–Penroses pseudoinvers kan beräknas med:

A + = ( A A ) 1 A {\displaystyle A^{+}=(A^{*}A)^{-1}A^{*}} .

Det följer då att A + {\displaystyle A^{+}} är vänsterinvers till A {\displaystyle A} .

Om raderna i A {\displaystyle A} är linjärt oberoende är A A {\displaystyle AA^{*}} inverterbar och Moore–Penroses pseudoinvers kan beräknas med:

A + = A ( A A ) 1 {\displaystyle A^{+}=A^{*}(AA^{*})^{-1}} .

Det följer då att A + {\displaystyle A^{+}} är högerinvers till A {\displaystyle A} .

Beräkning

Singulärvärdesfaktorisering

Om matrisen A {\displaystyle A} har singulärvärdesfaktoriseringen A = U Σ V {\displaystyle A=U\Sigma V^{*}} så fås A + = V Σ + U {\displaystyle A^{+}=V\Sigma ^{+}U^{*}} . Pseudoinversen av Σ {\displaystyle \Sigma } , som är en "nästan diagonal" matris med matrisens singulärvärden i diagonalen, genom att ersätta varje element σ i {\displaystyle \sigma _{i}} i diagonalen med 1 σ i {\displaystyle {\frac {1}{\sigma _{i}}}} . Exempel:

Σ = ( σ 1 0 0 0 0 σ 2 0 0 0 0 σ 3 0 )     Σ + = ( 1 σ 1 0 0 0 0 1 σ 2 0 0 0 0 1 σ 3 0 ) {\displaystyle \Sigma ={\begin{pmatrix}\sigma _{1}&0&0&0\\0&\sigma _{2}&0&0\\0&0&\sigma _{3}&0\end{pmatrix}}~~\Sigma ^{+}={\begin{pmatrix}{\frac {1}{\sigma _{1}}}&0&0&0\\0&{\frac {1}{\sigma _{2}}}&0&0\\0&0&{\frac {1}{\sigma _{3}}}&0\end{pmatrix}}}

Tillämpningar

Moore–Penroses pseudoinvers ger en minsta kvadrat-lösning till system av linjära ekvationer. Om systemet ges av A x = b {\displaystyle Ax=b} ges minsta kvadrat-lösningen av x = A + b {\displaystyle x=A^{+}b} .