[논문 리뷰] The polynomial eigenvalue problem is well conditioned for random inputs

Carlos Beltrán, Diego Armentano|arXiv (Cornell University)|2017. 06. 19.

Matrix Theory and Algorithms참고 문헌 22인용 수 4

한 줄 요약

이 논문은 입력 행렬이 i.i.d. 표준 복소 가우시안 분포를 따를 때 다항식 고유값 문제(PEVP)가 평균적으로 잘 조건화되어 있음을 입증한다. 랜덤 행렬 이론과 대수기하학을 사용하여 정확한 기대 제곱 조건수를 (d+1)n²−1/d로 도출함으로써, 큰 n과 d에 대해서도 수치적 해법이 일반적으로 안정적임을 보여준다.

ABSTRACT

We compute the exact value of the squared condition number for the polynomial eigenvalue problem, when the input matrices have entries coming from the standard complex Gaussian distribution, showing that in general this problem is quite well conditioned.

연구 동기 및 목표

• 랜덤 입력 행렬 하에서 다항식 고유값 문제(PEVP)의 평균 조건수를 결정하기 위해.
• PEVP에 대한 평균 조건수에 관한 이론적 결과가 부족한 바, 관련 문제에 대한 광범위한 연구에도 불구하고 이를 해결하기 위해.
• 기대 제곱 조건수를 계산하여 PEVP가 대규모 문제일지라도 일반적으로 잘 조건화되어 있음을 입증하기 위해.
• 대수기하학과 랜덤 행렬 이론을 활용해 다른 고유값 문제로의 분석을 일반화하기 위해.
• 실제로 수치적 해법이 높은 정밀도로 안정적이고 계산 가능하다는 이론적 근거를 제공하기 위해.

제안 방법

• 논문 [16, Th. 4.2]에서 유도된 PEVP의 조건수 공식을 사용하며, 이는 고유값, 고유벡터 및 행렬 다항식의 도함수에 의존한다.
• 복소 프로젝티브 공간에서 가우시안 입력을 가진 대수적 고유값 문제에 대한 기대 조건수에 관한 일반적 결과인 정리 2.1을 적용한다.
• Coarea 공식과 정규 자코비안을 사용하여 대수다양체 위의 적분을 계산하고 기하학적 구조를 기대값 계산과 연결한다.
• 조건수가 척도에 대해 불변임(μ(tA, (α,β)) = μ(A, (α,β)))이라는 사실을 활용하여 분산을 1로 정규화할 수 있다.
• 해의 다양체와 그 여차원의 구조를 분석하여 이론적 프레임워크의 적용 가능성을 검증한다.
• 동차 다양체와 관련된 핵심 적분을 계산하기 위해 보조정리 A.1을 활용하여 기대 제곱 조건수의 정확한 평가를 가능하게 한다.

실험 결과

연구 질문

• RQ1 입력 행렬이 i.i.d. 표준 복소 가우시안일 때 다항식 고유값 문제의 기대 제곱 조건수는 얼마인가?
• RQ2 특정 경우에서 악조건화될 수 있음에도 불구하고, 큰 n과 d에 대해서도 PEVP가 평균적으로 잘 조건화되어 있는가?
• RQ3 [16]의 이론적 프레임워크를 확장하여 PEVP 및 관련 문제의 정확한 기대 조건수를 계산할 수 있는가?
• RQ4 PEVP의 평균 조건수는 표준 또는 일반화된 고유값 문제와 비교해 어떻게 다른가?
• RQ5 유사한 방법을 고유값에 대수적 제약 조건이 있는 다른 구조적 고유값 문제에 적용할 수 있는가?

주요 결과

• i.i.d. 표준 복소 가우시안 분포를 따르는 입력 행렬에 대해 모든 고유값과 랜덤 입력에 대해 평균을 구한 PEVP의 기대 제곱 조건수는 정확히 (d+1)n²−1/d이다.
• 최대 조건수 μ_max(A)에 대해 E[μ_max(A)²] ≤ ((d+1)n²−1)n 이 성립하므로, 최악의 경우의 고유값조차도 일반적으로 잘 조건화되어 있음을 시사한다.
• 로그형 기대 조건수는 (3/2)log n + (1/2)log(d+1)로 유계이므로, 큰 n과 d에 대해서도 고정밀도 해법이 가능함을 의미한다.
• 조건수의 척도 불변성 덕분에 입력 분포의 분산 σ²이 어떤 값이든 결과가 성립한다.
• 이 방법은 다른 고유값 문제로도 일반화 가능하며, 2n개의 해를 가진 희소 PEVP와 2n개의 해를 가진 삼항 문제에 대해 각각 (3n²+3n)/4와 3n²의 기대 제곱 조건수를 도출한다.
• 해의 집합이 거의 확실히 유한하고, 관련 대수다양체가 정리 2.1에 필요한 기하적 조건을 만족함을 검증함으로써 이론적 프레임워크가 검증되었다.