[논문 리뷰] Jordan Normal and Rational Normal Form Algorithms

Bernard Parisse, Morgane Vaughan|arXiv (Cornell University)|2004. 12. 02.

Polynomial and algebraic computation참고 문헌 7인용 수 3

한 줄 요약

이 논문은 팔레프 공식과 (λI − A)의 코마트릭스를 사용하여 행렬의 조르당 표준형과 유리 표준형을 계산하는 결정론적 알고리즘을 제시한다. 비영인 B(λ₀)의 열이 고유값 λ₀에 대한 고유벡터임을 이용하여, 특성다항식이 인수분해된 경우 O(n⁴)의 체(field) 연산으로 효율적인 계산이 가능하다.

ABSTRACT

In this paper, we present a determinist Jordan normal form algorithms based on the Fadeev formula: \[(λ\cdot I-A) \cdot B(λ)=P(λ) \cdot I\] where $B(λ)$ is $(λ\cdot I-A)$'s comatrix and $P(λ)$ is $A$'s characteristic polynomial. This rational Jordan normal form algorithm differs from usual algorithms since it is not based on the Frobenius/Smith normal form but rather on the idea already remarked in Gantmacher that the non-zero column vectors of $B(λ_0)$ are eigenvectors of $A$ associated to $λ_0$ for any root $λ_0$ of the characteristical polynomial. The complexity of the algorithm is $O(n^4)$ field operations if we know the factorization of the characteristic polynomial (or $O(n^5 \ln(n))$ operations for a matrix of integers of fixed size). This algorithm has been implemented using the Maple and Giac/Xcas computer algebra systems.

연구 동기 및 목표

• 대수적으로 닫힌 체 위에서 행렬의 조르당 표준형을 계산하는 결정론적 알고리즘을 개발하는 것.
• 체 확장을 필요로 하지 않고 임의의 체 위에서 유리 표준형을 계산하는 방법으로 이 방법을 확장하는 것.
• 팔레프/스미스 표준형 기반 접근법의 대안으로, (λI − A)의 코마트릭스를 사용하는 것.
• 특성다항식의 인수분해가 알려진 조건 하에서 O(n⁴) 복잡도를 달성하여 병렬 처리에 적합하게 만드는 것.
• 컴퓨터 대수 시스템인 Maple과 Giac/Xcas에 알고리즘을 구현하고 통합하여 실용적 사용을 가능하게 하는 것.

제안 방법

• 특성다항식 P(λ)와 (λI − A)의 코마트릭스인 행렬다항식 B(λ)를 계산하기 위해 팔레프 알고리즘을 활용한다.
• 항등식 (λI − A)B(λ) = P(λ)I를 이용하여 추적 기반 공식을 통해 B(λ) 계수의 재귀 관계를 유도한다.
• 고유값 λ₀에서 B(λ₀)를 효율적으로 평가하기 위해 호른어의 방법(Hörner's method)을 적용하여 고유벡터를 추출한다.
• P(λ)의 기약 인수 Q에 대해 Q(A)를 포함하는 선형계를 풀어 일반화된 고유벡터를 구성한다.
• 이항계수 항등식을 활용하여 재귀식 (9)을 통해 기저 벡터에서 고차원 순환벡터를 재구성한다.
• 복소수 조르당 표준형 계산과 1차 이상의 기약 인수에 대해 이를 유리 표준형으로 변환하는 하이브리드 방법을 Maple과 Giac/Xcas에 구현한다.

실험 결과

연구 질문

• RQ1 코마트릭스 (λI − A)와 고유값 평가만을 사용하여 조르당 표준형을 결정론적으로 계산할 수 있는가?
• RQ2 체 확장을 피하면서도 행렬다항식 연산만으로 유리 표준형을 어떻게 구성할 수 있는가?
• RQ3 팔레프/스미스 표준형 방법과 비교해 볼 때, 코마트릭스 기반 알고리즘의 계산 복잡도는 어떻게 되는가?
• RQ4 이 알고리즘은 효율적으로 병렬화하거나 기호적 오버헤드를 최소화한 컴퓨터 대수 시스템에 어떻게 통합할 수 있는가?
• RQ5 다항식 행렬 연산을 통해 기저 고유벡터에서 일반화된 고유벡터를 어떻게系통적으로 재구성할 수 있는가?

주요 결과

• 특성다항식이 완전히 인수분해된 경우 알고리즘이 O(n⁴)의 체 연산으로 조르당 표준형을 계산한다.
• 고정된 크기의 정수 행렬에 대해서는 다항식 인수분해 비용으로 인해 복잡도가 O(n⁵ ln n)으로 증가한다.
• 유리 표준형은 먼저 의사-유리형을 계산한 후 구조적 변환을 통해 변환하여 구성된다.
• 팔레프 또는 스미스 표준형에 대한 의존도를 피함으로써, 유리 표준형 계산을 위한 새로운 길을 제시한다.
• 알고리즘은 Maple과 Giac/Xcas에 성공적으로 구현되었으며, jordan 및 rat_jordan와 같은 함수를 통해 직접 사용이 가능하다.
• 핵심 통찰은 B(λ₀)의 비영인 열이 고유값 λ₀에 대한 고유벡터임을 보여주며, 이는 (A − λ₀I)v = 0을 풀지 않고도 고유벡터를 직접 추출할 수 있음을 의미한다.