2009-07-20 2 views
2

불충분 한 방정식 시스템 (예 : 시스템을 설명하는 행렬이 키가 크다)에 단일 솔루션을 생성하는 함수를 작성하려고합니다. 이렇게하기 위해서 LAPACK 문서에서 Mathematica와 TI 계산기의 rref() 함수와 비슷한 행렬을 줄이는 방법을 찾고있었습니다. 내가 만난 가장 가까운 것은 http://software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/이 작은 실이었다. 그러나이 스레드는 단순히 "U"위 삼각 행렬 (각 행을 대각선으로 나눔)을 사용하는 것이 행렬의 축소 된 형태와 동일하다는 것을 의미하는 것으로 보입니다. 내가 직접 rref() 코드를 작성할 수는 있지만 LAPACK이 유명한 성능을 얻을 수 있다고는 생각하지 않는다.Lapack의 행 축소

1) 제대로 표현되지 않은 시스템에 대한 특정 해결책을 얻는 더 좋은 방법이 있습니까?
2) 그렇지 않다면 LAPACK이 행렬을 줄이는 방법이 있습니까?

감사합니다. 이것

+0

LU 분해의 U 행렬은 'ref'이다 sgelsx LAPACK 참조. 각 행을 선행 항목으로 나누면'ref'를'rref'로 변환합니다. – Sharpie

답변

3

한 자주 사용되는 방법이 최소 제곱 해법이다 이제

+0

나는 이것을 시도하고 그것이 작동하는 것 같습니다. 감사. – Litherum