masalah SVD tertimbang?


12

Diberikan dua matriks A dan B , saya ingin mencari vektor x dan y , sehingga,

minij(AijxiyjBij)2.
Dalam bentuk matriks, saya mencoba untuk meminimalkan norma Frobenius dari Adiag(x)Bdiag(y)=AB(xy) .

Secara umum, saya ingin mencari beberapa vektor satuan x dan y dalam bentuk

minij(Aijk=1nsixi(k)yj(k)Bij)2.
di mana si adalah koefisien nyata positif.

Ini sama dengan dekomposisi nilai singular (SVD) ketika (B)ij=1 .

Adakah yang tahu apa sebutan masalah ini? Apakah ada algoritma terkenal seperti SVD untuk solusi masalah seperti itu?

(dimigrasikan dari math.SE)


Saya percaya ini adalah SVD Generalized . Entri Wikipedia tidak terlalu detail, jadi Anda mungkin harus memeriksa sumber yang tertaut. Secara khusus, halaman 466 tautan Google buku ini mungkin bermanfaat.
ely

1
Bagi saya, ini tidak terlihat seperti SVD umum. Terutama karena B tidak harus diagonal atau simetris, sehingga setiap atau dapat muncul berkali-kali dalam jumlah. xy
Victor Liu

B tidak perlu diagonal atau simetris dalam SVD umum. Kedua tautan yang saya berikan menunjukkan bahwa A dan B masing-masing dapat berupa matriks umum yang bernilai kompleks dengan dimensi M-by-N dan P-by-N.
ely

Terima kasih atas sarannya @EMS. Saya akan sangat menghargai jika Anda dapat menguraikan koneksi.
Memming

Jawaban:


8

Ini jauh dari SVD umum.

Jika B adalah matriks positif, Anda dapat menggunakan paket saya BIRSVD http://www.mat.univie.ac.at/~neum/software/birsvd/

Makalah ini http://www.mat.univie.ac.at/~neum/software/birsvd/svd_incomplete_data.pdf yang menjelaskan metode di sana juga memberikan referensi yang dapat Anda pertimbangkan untuk melakukan pencarian literatur.


Ah, mengubah masalah menjadi pendekatan peringkat rendah tertimbang! Terima kasih banyak!
Memming

Untuk menambahkan rincian ke jawaban @ Arnold, masalah ini dapat diubah menjadi masalah pendekatan peringkat rendah tertimbang di mana tujuannya adalah untuk meminimalkan norma tertimbang, bukan norma Frobenius. mana dan menunjukkan produk Schur (alias Produk Hadamard). ||Csixiyi||W2||C||W=||CW||F
Memming

Iya. Ini memberi nama yang bagus untuk masalah Anda. Bagaimana mengatasinya adalah masalah yang berbeda. Ini bukan masalah standar dan cukup sulit untuk menemukan algoritma yang cepat dan andal.
Arnold Neumaier

@ArnoldNeumaier ini bagus, terima kasih. apakah mungkin untuk mendapatkan lisensi dan pemberitahuan hak cipta dengan kode Anda? Seperti sekarang ini adalah perangkat lunak berpemilik. Jika Anda melepaskannya di bawah GPLv3 atau yang kompatibel, ia mungkin menemukan jalannya ke paket aljabar linier GNU Octave.
JuanPi
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.