Was kann ich über die Kovarianzmatrix sagen, wenn eine inverse Kovarianzmatrix spärlich ist?


Antworten:


4

Wie bereits von Yair kommentiert, gibt es keine spezifische Sparsity-Bedingung der inversen Kovarianzmatrix, die die tatsächliche Kovarianzmatrix beeinflusst oder umgekehrt. Alles andere als triviale Sparsity-Matrixmuster (dh Diagonale) haben keine Garantie dafür, dass sie sich sowohl auf eine bestimmte Matrix als auch auf deren Inverse auswirken. Sogar tridiagonale Matrizen können leicht nicht spärliche Inversen aufweisen.

In bestimmten Fällen, in denen die Sparsity der Matrix in Blöcken auftritt, können Sie möglicherweise einige Ergebnisse ableiten, die sich aus dem pseudoinversen Blockmatrix- Algorithmus ergeben, der Folgendes besagt:

[EINB.C.D.]]- -1=[(EIN- -B.D.- -1C.)- -1- -EIN- -1B.(D.- -C.EIN- -1B.)- -1- -D.- -1C.(EIN- -B.D.- -1C.)- -1(D.- -C.EIN- -1B.)- -1]]

aber das ist wahrscheinlich auch so (rein anekdotisch habe ich versucht, Sparsity-Muster durch die Cholesky-Zerlegung einer PSD-Matrix aufzuerlegen, aber ich habe meinen Versuch-und-Irrtum-Streifzug nicht bestanden). Möglicherweise möchten Sie auch den Cuthill-McKee-Algorithmus (CM) untersuchen, wenn Sie erwarten, dass sich ein Adjazenzmerkmal in der Kovarianzmatrix widerspiegelt. Der CM-Algorithmus permutiert eine dünn besetzte Matrix mit einem symmetrischen Sparsity-Muster in eine Bandmatrixform mit einer kleinen Bandbreite. Dies könnte dazu beitragen, eine gewisse Sparsity gegenüber den nicht diagonalen Einträgen der inversen Matrix beizubehalten, dies ist jedoch nicht garantiert. (Das Anwenden von CM - wenn dies sinnvoll ist - kann für bestimmte Anwendungen (z. B. in 2D-Glättungsroutinen) sehr hilfreich sein und Ihre Berechnungen erheblich beschleunigen.)


2
(+1) Da die blockdiagonalen Matrizen einer bestimmten Form ein Ring sind, haben ihre Inversen (wann immer sie existieren) dieselbe blockdiagonale Struktur und bewahren so viel Sparsity-Muster. Als extremes Beispiel sind diagonale Matrizen blockdiagonal, wodurch der von @Aksakal aufgezeigte Fall veranschaulicht wird. Am weitesten kann man in diese Richtung gehen, indem man blockdiagonale Matrizen durch Permutationsmatrizen konjugiert (wodurch offensichtlich alle Null- und Nicht-Null-Einträge erhalten bleiben, sondern nur verschoben werden).
whuber

1
(+1) Danke für diesen Kommentar (kurze Antwort wirklich). Es ist wirklich aufschlussreich. Ich werde es auf jeden Fall in Zukunft in Betracht ziehen.
usεr11852
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.