가중 SVD 문제?


12

두 행렬 주어 및 B , I는 벡터 찾으려 XY 와 같은 것을 분간 Σ I J ( I J - X Y j 개의 B 용 I의 J ) 2 . 행렬 형태로, 나는 A - diag ( x ) B diag ( y ) = A B ( x y 의 Frobenius 규범을 최소화하려고합니다.ABxy

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

일반적으로, I는 복수 개의 단위 벡터 찾으려 Y를 "형태이야 Σ I J ( I J - N Σ의 K = 1I X ( K ) Y ( K ) J의 B의 I , J ) 2 . 이야 의는 양의 실수 계수이다.xy

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

이는 때 특이 값 분해 (SVD)와 같습니다 .(B)ij=1

이 문제가 무엇인지 아는 사람이 있습니까? 그러한 문제의 해결을 위해 SVD와 같은 잘 알려진 알고리즘이 있습니까?

(math.SE에서 마이그레이션)


나는 이것이 일반화 된 SVD 라고 믿는다 . Wikipedia 항목은 매우 상세하지 않으므로 링크 된 소스를 확인해야합니다. 특히이 Google 도서 링크의 466 페이지 가 도움이 될 수 있습니다.
ely

1
xy

일반화 된 SVD에서 B는 대각선이나 대칭 일 필요는 없습니다. 내가 제공 한 두 링크는 ​​A와 B가 각각 MxN 및 PxN 차원의 일반 복소수 행렬 일 수 있음을 나타냅니다.
ely

@EMS 제안에 감사드립니다. 연결을 자세히 설명해 주시면 감사하겠습니다.
Memming

답변:


8

이것은 일반화 된 SVD와는 거리가 멀다.

B가 양의 행렬이면 내 패키지 BIRSVD를 사용할 수 있습니다 http://www.mat.univie.ac.at/~neum/software/birsvd/

이 방법을 설명하는 http://www.mat.univie.ac.at/~neum/software/birsvd/svd_incomplete_data.pdf 문서 는 문헌 검색을 고려할 수있는 참고 자료도 제공합니다.


아, 문제를 가중 낮은 순위 근사치로 변환! 고마워요!
Memming

||Csixiyi||W2||C||W=||CW||F

예. 이것은 당신의 문제에 좋은 이름을줍니다. 그것을 해결하는 방법은 다른 문제입니다. 표준 문제는 아니며 빠르고 안정적인 알고리즘을 찾는 것이 매우 까다 롭습니다.
Arnold Neumaier

@ArnoldNeumaier 감사합니다. 코드로 라이센스 및 저작권 고지를받을 수 있습니까? 지금은 독점 소프트웨어입니다. GPLv3 또는 호환 버전으로 배포하면 GNU Octave의 선형 대수 패키지로가는 길을 찾을 수 있습니다.
JuanPi
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.