Jika matriks kovarians terbalik jarang, apa yang bisa saya katakan tentang matriks kovarians?

8

Bagaimana kondisi sparsity pada matriks kovarians terbalik mempengaruhi matriks kovarians sebenarnya?

Gegat
sumber
Sayangnya, tidak
Yair Daon
Pasti ada sesuatu. Misalnya matriks identitas jarang, dan kebalikannya juga.
Aksakal

Jawaban:

4

Seperti yang sudah dikomentari oleh Yair, tidak ada kondisi sparsity spesifik dari matriks kovarians terbalik yang mempengaruhi matriks kovarians aktual atau sebaliknya. Apa pun selain pola matriks sparsity sepele (mis. Diagonal) tidak memiliki jaminan bahwa mereka akan tercermin pada matriks tertentu dan kebalikannya. Bahkan matriks tridiagonal dapat dengan mudah memiliki inversi yang jarang.

Untuk kasus-kasus tertentu di mana sparsity dari matriks terjadi di blok Anda mungkin dapat memperoleh beberapa hasil yang berasal dari algoritma pseudoinverse matriks blok yang menyatakan bahwa:

[ABCD]1=[(ABD1C)1A1B(DCA1B)1D1C(ABD1C)1(DCA1B)1]

tapi itu mungkin tentang hal itu (murni anekdot, saya telah mencoba untuk memaksakan pola sparsity melalui dekomposisi Cholesky dari sebuah matriks PSD tetapi saya gagal dalam percobaan coba-coba saya yang salah). Anda mungkin juga ingin mempertimbangkan melihat ke dalam algoritma Cuthill-McKee (CM) jika Anda mengharapkan beberapa fitur adjacency tercermin dalam matriks kovarians. Algoritma CM memungkinkan matriks sparse yang memiliki pola sparsity simetris ke dalam bentuk matriks band dengan bandwidth kecil, ini mungkin membantu mempertahankan beberapa sparsitas terhadap entri off-diagonal dari matriks invers tetapi itu tidak dijamin. (Menerapkan CM - jika masuk akal - dapat sangat membantu untuk aplikasi tertentu (misalnya dalam rutinitas perataan 2D) dan secara signifikan dapat mempercepat perhitungan Anda.)

usεr11852
sumber
2
(+1) Karena matriks blok-diagonal dari bentuk apa pun yang diberikan adalah sebuah cincin, inversinya (kapan pun ada) memiliki struktur blok-diagonal yang sama sehingga mempertahankan banyak pola sparsity. Sebagai contoh ekstrem, matriks diagonal adalah blok-diagonal, dengan demikian memberikan contoh kasus yang ditunjukkan oleh @Aksakal. Yang paling jauh yang bisa menuju ke arah ini adalah untuk mengkonjugasikan matriks blok-diagonal dengan matriks permutasi (yang jelas mempertahankan semua entri nol dan bukan nol, tetapi hanya memindahkannya).
Whuber
1
(+1) Terima kasih atas komentar ini (jawaban singkat benar-benar). Ini sangat mendalam. Saya pasti akan mempertimbangkannya di masa depan.
usεr11852