가중 SVD 문제?


12

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

min∑ij(Aij−xiyjBij)2.
.A−diag(x)⋅B⋅diag(y)=A−B∘(xy⊤)

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

min∑ij(Aij−∑k=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

||C−∑sixiyi⊤||W2||C||W=||C∘W||F∘

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

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