가능한 한 간단한 패키지로 합리적인 수준의 선형 대수 기능을 제공하기 위해 헤더 전용 행렬 라이브러리에서 작업하고 있으며 현재의 최신 상태가 무엇인지 조사하려고합니다. 복소수 행렬. 나는 2 단계 분해, 이중 대각선 화와 단일 값 계산을하고 있습니다. 지금 나는 bidiagonalization (나는 LAPACK이이를 사용하여 생각), 그리고 생각의 호주 방법을 사용하고 그 약의가 현재수록 …
다음 최소 제곱 문제에 대한 방법을 추천 할 수 있습니까? 다음 을 최소화하는 을 . 여기서 R 은 단일 (회전)입니다. 매트릭스.R∈R3×3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i=0N(Rxi−bi)2→min∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i - b_i)^2 \rightarrow \minRRR ∑i=0N(Axi−bi)2→min∑i=0N(Axi−bi)2→min\sum\limits_{i=0}^N (Ax_i - b_i)^2 \rightarrow \min (임의의 A∈R3×3A∈R3×3A \in \mathbb{R}^{3 \times 3} ) 을 최소화하여 대략적인 해결책을 얻을 수 있습니다. …
한 번에 하나씩 특이 값을 계산하는 잘린 SVD 알고리즘이 있습니까? 내 문제 : 큰 밀도 행렬 M 의 첫 kkk 특이 값 (및 특이 벡터) 을 계산하고 싶지만 적절한 k 값이 무엇인지 모르겠습니다 . M 은 크므로 효율성 때문에 전체 SVD를 평가하지 않고 나중에 가장 작은 SV를 잘라 버릴 것입니다.MMMkkkMMM …
모형 축소를 위해 행렬 의 20 개 가장 큰 특이 값과 연관된 왼쪽 특이 벡터를 계산하려고합니다 . 여기서 N ≈ 10 6 및 k ≈ 10 3 입니다. 불행히도 내 행렬 A 는 아무런 구조없이 밀도가 높습니다.∈ R엔, kㅏ∈아르 자형엔,케이A \in \mathbb R^{N,k}엔≈ 106엔≈106N\approx 10^6k ≈ 10삼케이≈10삼k\approx 10^3ㅏㅏA 이 크기의 …
나는 (인텔 MKL의 SVD를 사용하고 dgesvdSciPy를 통해) 나는 사이의 정밀도를 변경하면 결과가 크게 다른 것으로 나타났습니다 float32과 float64내 행렬이 심하게 조절한다 /하지 전체 순위. 결과를 float32-> float64변경에 둔감하게 만들기 위해 추가해야하는 최소한의 정규화에 대한 가이드가 있습니까? 특히 하면 의 규범이 와 사이의 정밀도를 변경할 때 약 1만큼 이동합니다 . 규범 …