인증된 실수 고유값 위치 추정: Gershgorin 디스크와 Hermite 행렬의 융합

초록

본 논문은 실수 고유값을 정확히 구간으로 인증하는 새로운 혼합 기법을 제시한다. Gershgorin 디스크를 이용해 후보 영역을 설정하고, 일반화된 Hermite 행렬을 통해 해당 구간에 실근이 존재함을 증명한다. 구간은 이진 탐색 방식으로 정밀히 축소 가능하며, 구현은 Julia로 제공된다.

상세 분석

이 연구는 고유값 문제를 다항식 근 근사와 구간 인증이라는 두 축으로 재구성한다. 먼저, 행렬 A의 특성 다항식 p(x)를 La Budde 방법으로 안정적으로 계산한다. La Budde는 행렬을 상삼각형 형태(Hessenberg)로 변환한 뒤, 실수 연산만으로 계수를 재귀적으로 구해 O(n³) 복잡도를 유지한다는 장점이 있다. 특성 다항식이 확보되면, 고유값은 p(x)의 실근과 동일함을 이용한다.

다음 단계는 Gershgorin 디스크 정리를 활용해 고유값이 존재할 가능성이 높은 복소 평면 영역을 빠르게 추정하는 것이다. 각 행 i에 대해 중심 a_{ii}와 반지름 R_i=∑{j≠i}|a{ij}| 로 정의된 디스크 D_i는 모든 고유값이 ∪_i D_i 안에 포함된다는 보장을 제공한다. 여기서 중요한 점은 실수 고유값은 실축 위에 놓이므로, 각 디스크와 실축의 교집합이 실제 후보 구간이 된다.

핵심 기여는 이 후보 구간에 대해 Hermite 행렬을 이용해 “실근 존재 여부”를 인증하는 절차다. 일반화된 Hermite 행렬 H_q(p)는 보조 다항식 q(x)에 따라 정의되며, q(x)=1인 경우 H_1(p)의 서명 σ(H_1) 은 p(x)의 전체 실근 개수를 바로 제공한다(Hermite 정리). q(x)=(x−a)(x−b) 로 설정하면 구간