Entspricht jede halbpositive bestimmte Matrix einer Kovarianzmatrix?

12

Es ist bekannt, dass eine Kovarianzmatrix halbpositiv definitiv sein muss. Ist das Gegenteil der Fall?

Das heißt, entspricht jede halbpositive bestimmte Matrix einer Kovarianzmatrix?

Jingjings
quelle

Antworten:

12

Wenn ich mich hier an die Definitionen von PD und PSD halte, denke ich schon, da wir dies konstruktionsbedingt tun können. Ich gehe für ein etwas einfacheres Argument davon aus, dass Sie es für Matrizen mit realen Elementen meinen, aber mit entsprechenden Änderungen würde es sich auf komplexe Matrizen erstrecken.

Sei eine echte PSD-Matrix; Von der Definition, mit der ich verlinkt habe, wird es symmetrisch sein. Jede reelle symmetrische positive definitive Matrix A kann als A = L L T geschrieben werden . Dies kann durch L = Q √ erfolgenAAA=LLTwennA=QDQTmit orthogonalemQund DiagonaleDundL=QDQTA=QDQTQD als Matrix der Komponente wise QuadratwurzelnD. Es muss also nicht der volle Rang sein.DD

Sei eine Vektor-Zufallsvariable der entsprechenden Dimension mit der Kovarianzmatrix I (die einfach zu erstellen ist).ZI

Dann hat Kovarianzmatrix A .LZA

[Zumindest ist das theoretisch. In der Praxis gibt es verschiedene numerische Probleme, mit denen Sie sich befassen müssen, wenn Sie gute Ergebnisse erzielen möchten, und aufgrund der üblichen Probleme bei der Gleitkommaberechnung erhalten Sie nur ungefähr das, was Sie benötigen. das heißt, die Varianz eines berechneten würde in der Regel nicht genau A . Aber so etwas ist immer ein Problem, wenn wir Dinge tatsächlich berechnen wollen.]LZ A

Glen_b - Monica neu starten
quelle
Es stimmt zwar, dass eine Zerlegung ohne vollen Rang möglich ist, funktioniert der Cholesky - Algorithmus nur mit herkömmlichen A . Ohne vollen Rang kann es also keine Cholesky-Zerlegung sein. Computergestützt könnte man diese Zerlegung im Singularfall durch Diagonalisierung durchführen. (Obwohl dies weitaus teurer ist)A=LLA
Horst Grünbusch
@Horst: Warum sollte L = Q √ sein ?unteres Dreieck sein? L=QDQT
Amöbe sagt Reinstate Monica
1
@amoeba Man könnte es zwar so organisieren, aber es muss nicht dreieckig sein, damit das Argument funktioniert - es ist eine Funktion des Cholesky, aber es ist nicht erforderlich, damit das Ergebnis funktioniert.
Glen_b -Reinstate Monica
1
@Glen Ist Symmetrie eine notwendige Bedingung, um PSD zu sein, oder ist diese Definition eine von vielen?
114
1
@ 114 für die Beziehung zwischen symmetrisch und PSD siehe math.stackexchange.com/questions/516533/…
Frank