vấn đề trọng số SVD?


12

ABxy

minij(AijxiyjBij)2.
Adiag(x)Bdiag(y)=AB(xy)

Nói chung, tôi muốn tìm nhiều vectơ đơn vị và ở dạng trong đó là các hệ số thực dương.xy

minij(Aijk=1nsixi(k)yj(k)Bij)2.
si

Điều này tương đương với phân rã giá trị số ít (SVD) khi .(B)ij=1

Có ai biết vấn đề này được gọi là gì không? Có một thuật toán nổi tiếng như SVD cho giải pháp cho vấn đề như vậy không?

(di chuyển từ math.SE)


Tôi tin rằng đây là SVD tổng quát . Mục nhập Wikipedia không chi tiết lắm, vì vậy bạn có thể nên kiểm tra các nguồn được liên kết. Đặc biệt, trang 466 của liên kết sách Google này có thể hữu ích.
ely

1
Đối với tôi, nó không giống bất cứ thứ gì như SVD tổng quát. Đặc biệt vì B không nhất thiết phải là đường chéo hoặc đối xứng, vì vậy mỗi hoặc có thể xuất hiện nhiều lần trong tổng. xy
Victor Liu

B không cần phải là đường chéo cũng như đối xứng trong SVD tổng quát. Cả hai liên kết tôi cung cấp chỉ ra rằng A và B có thể là các ma trận có giá trị phức tạp chung có thứ nguyên M-by-N và P-by-N tương ứng.
ely

Cảm ơn lời đề nghị @EMS. Tôi sẽ đánh giá cao nếu bạn có thể xây dựng kết nối.
Ghi nhớ

Câu trả lời:


8

Điều này là xa SVD tổng quát.

Nếu B là một ma trận dương, bạn có thể sử dụng gói BIRSVD của tôi http://www.mat.univie.ac.at/~neum/software/birsvd/

Bài viết http://www.mat.univie.ac.at/~neum/software/birsvd/svd_incomplete_data.pdf mô tả phương pháp đó cũng cung cấp các tài liệu tham khảo mà bạn có thể xem xét để thực hiện tìm kiếm tài liệu.


Ah, chuyển vấn đề thành xấp xỉ thứ hạng thấp có trọng số! Cảm ơn rất nhiều!
Ghi nhớ

Để thêm chi tiết vào câu trả lời của @ Arnold, vấn đề này có thể được chuyển thành vấn đề xấp xỉ thứ hạng thấp có trọng số trong đó mục tiêu là giảm thiểu định mức có trọng số thay vì chỉ tiêu Frobenius. nơi và biểu thị sản phẩm Schur (aka Sản phẩm Hadamard). ||Csixiyi||W2||C||W=||CW||F
Ghi nhớ

Đúng. Điều này cho một cái tên đẹp cho vấn đề của bạn. Làm thế nào để giải quyết nó là một vấn đề khác nhau. Đây không phải là một vấn đề tiêu chuẩn và khá khó để tìm ra một thuật toán vừa nhanh vừa đáng tin cậy.
Arnold Neumaier

@ArnoldNeumaier điều này thật tuyệt, cảm ơn bạn. có thể nhận được giấy phép và thông báo bản quyền với mã của bạn không? Như bây giờ nó là phần mềm độc quyền. Nếu bạn phát hành nó theo GPLv3 hoặc tương thích, nó có thể tìm đường đến gói đại số tuyến tính của GNU Octave.
JuanPi
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.