Iterative Methods for Computing the Moore--Penrose Inverse of Split-Quaternion Matrices with Applications

기술 요약: 분할 사원수 행렬의 무어-펜로즈 역행렬 계산을 위한 반복법

문제 정의
본 논문은 분할 사원수($\mathbb{Q}_s$) 대수 위의 행렬에 대한 무어-펜로즈(Moore–Penrose, MP) 역행렬 계산 문제를 다룬다. 해밀턴 사원수와 달리, 분할 사원수는 부정적 이차 형식을 가지며, 이로 인해 자명하지 않은 제로 디바이저(zero divisors)와 멱영 요소(nilpotent elements)가 존재한다. 이러한 대수적 특징은 역행렬의 정의와 표준 스펙트럼 이론의 적용을 복잡하게 만든다. 특히, 이차 형식의 부정적 특성은 실수 또는 복소수 행렬로부터의 스펙트럼 결과의 직접적인 전이를 방해하며, 제로 디바이저의 존재는 특이(singular) 및 계수 부족(rank-deficient) 사례에 대한 세심한 처리를 요구한다. 저자들은 이 환경에서 MP 역행렬을 정의하기 위한 견고한 프레임워크를 구축하고, 특히 직사각형 및 계수 부족 행렬에 대한 효율적인 계산 알고리즘을 개발하는 것을 목표로 한다.

방법론
핵심적인 방법론적 기여는 $2 \times 2$ 실수 표현($\sigma$)과 $i$-켤레 전치($H$)에 기반한 일관된 프레임워크의 확립이다.

1. 실수 표현: 저자들은 곱셈을 보존하고 $i$-켤레 전차를 표준 실수 전치로 매핑하는 생단사 함수 $\sigma: \mathbb{Q}_s \to \mathbb{R}^{2 \times 2}$를 정의한다. 행렬 $A \in \mathbb{Q}_s^{m \times n}$에 대하여, 실수 표현체 $A^\sigma \in \mathbb{R}^{2m \times 2n}$를 사용하면 MP 역행렬을 실수 MP 역행렬을 통해 직접 정의할 수 있다: $A^\dagger = \sigma^{-1}((A^\sigma)^\dagger)$. 이 정의는 제로 디바이저의 처리를 명확히 하고 역행렬의 존재와 유일성을 보장한다.
2. 뉴턴-슐츠 반복법(Newton–Schulz Iteration): 본 논문은 MP 역행렬을 계산하기 위한 뉴턴-슐츠 반복법을 분석한다. 두 가지 구현 방식이 제시된다:
   • 임베디드(Embedded): 실수 표현체 $A^\sigma$에 표준 실수 뉴턴-슐츠 반복법을 적용한다.
   • 네이티브(Native): 분할 사원수 산술을 직접 사용하여 수행한다 ($X_{k+1} = X_k(2I - AX_k)$).
     저자들은 이 두 수열이 표현 $\sigma$ 하에서 정확히 동일함을 증명한다.
3. 수렴 분석: 실수 표현체의 얇은 특이값 분해(SVD)를 사용하여, 저자들은 반복 항과 잔차(residual)에 대한 정확한 공식을 도출한다. 초기 스케일링 파라미터 $\alpha$가 $0 < \alpha < 2/\|A^\sigma\|_2^2$를 만족할 때, 이 반복법은 임의의 계수(계수 부족 사례 포함)를 가진 행렬에 대해 MP 역행렬로 수렴함을 입증한다. 분석은 잔차와 오차 사이의 정확한 관계를 제공하며, 잔차의 이차 수렴성을 보여준다.
4. 다항식 초기화: 수렴을 가속화하기 위해, 소리우(Souriau)의 역수 재귀에서 영감을 받은 저차 다항식 초기화를 제안한다. 이는 최소제곱(LS) 접근법을 사용하여 정규화된 그람 인자($A^\sigma (A^\sigma)^T$ 또는 $(A^\sigma)^T A^\sigma$)의 역수에 적합한 다항식을 피팅하는 과정을 포함한다. 이 방법은 초기 잔차가 수렴 영역 내에 있는지 확인하기 위한 스펙트럼 수용 테스트를 포함하며, 테스트 실패 시 안전한 폴백(fallback)으로 스케일링된 전치를 사용한다.

주요 기여

• 일관된 프레임워크: 본 논문은 실수 표현과 $i$-켤레 전치를 통해 분할 사원수 행렬의 MP 역행렬을 정의하고, 존재성, 유일성 및 기본 성질(제로 디바이저에 대한 스칼라 공식 포함)을 증명한다.
• 수렴 이론: 직사각형 및 계수 부족 행렬에 대한 뉴턴-슐츠 반복법의 완전한 수렴 분석을 제공한다. 저자들은 이전의 해밀턴 사원수 행렬 분석에서 발견된 풀 랭크(full-rank) 제한을 제거하고, 정확한 오차 재귀 및 잔차-오차 항등식을 도출한다.
• 구현의 동등성: 실수 기반 임베디드 뉴턴-슐츠 반복법과 네이티브 분할 사원수 뉴턴-슐츠 반복법 사이의 정확한 동등성을 입증하여, 이론적 보장 손실 없이 네이티브 산술을 사용할 수 있음을 검증한다.
• LS–Gram 초기화: 그람 행렬의 특정 스펙트럼에 적응하는 새로운 다항식 초기화를 도입한다. 이는 표준 스케일링 전치 방식에 비해 필요한 반복 횟수를 줄여준다는 것을 보여준다.
• CUR 근사: 저자들은 교차(cross) 및 CUR 근사 이론을 분할 사원수 행렬으로 확장한다. 고정된 샘플 행 및 열에 대해 $C^\dagger A R^\dagger$가 프로베니우스 최적(Frobenius-optimal) 중간 인자임을 증명한다. 또한, 저렴한 교차 인자 $W^\dagger$($W$는 교차 행렬)가 정확한 재구성을 생성하는 랭크-리빌링(rank-revealing) 조건을 확립한다.

결과
MATLAB 및 Python으로 수행된 수치 실험은 이론적 발견을 검증한다:

• 정확도: 반복법은 직접적인 SVD 기반 계산과 유사한 정확도(상대 오차 약 $10^{-11}$ ~ $10^{-14}$)를 달ink한다.
• 효율성: 중간 크기의 조밀 행렬(dense matrices)에 대해, 뉴턴-슐츠 반복법은 직접적인 SVD 계산보다 현저히 빠르다 (예: $200 \times 250$ 행렬의 경우 약 3.7배 빠름).
• 초기화 성능: LS–Gram 다항식 초기화는 스케일링된 전치 방식에 비해 뉴턴-슐츠 반복 횟수를 일관되게 감소시킨다. 조건수가 $10^4$에 달하는 실험에서, LS–Gram 방식은 스케일링된 전치 대비 반복 횟수를 약 20% 줄였으며, 정확한 스펙트럼 경계가 필요한 체비쇼프(Chebyshev) 다항식 초기화와 비교해도 경쟁력 있는 성능을 보였다.
• CUR 근사: 실험은 랭크-리빌링 조건이 충족될 때 $W^\dagger$를 사용하는 교차 근사가 저계수(low-rank) 행렬에 대해 정확한 재구성을 달성할 수 있음을 확인하며, 핵심 역행렬의 반복적 계산이 직접적인 방법보다 훨씬 빠름을 보여준다.

의의 및 주장
본 논문은 제로 디바이저와 계수 부족을 실수 표현을 통해 명시적으로 처리하는, 분할 사원수 행렬의 MP 역행렬에 대한 최초의 일관된 프레임워크를 제공한다고 주장한다. 실수 행렬에서 뉴턴-슐츠 반복법의 분석을 분할 사원수 환경으로 전이함으로써, 저자들은 임의의 계수를 가진 행렬에 대한 엄격한 수렴 보장을 확립한다. 제안된 LS–Gram 초기화는 표준 초기화 전략보다 계산 비용을 줄이면서도 (체비쇼프 방법과 달리) 사전 스펙트럼 경계 지식이 필요 없는 실용적인 개선책을 제공한다. 마지막으로, CUR 근사 결과의 확장은 분할 사원수 신호 처리 및 물리 응용 분야에서의 저계수 근사를 위한 이론적 토대를 제공하며, 특정 인자의 최적성과 정확한 재구성을 위한 조건을 규명한다. 저자들은 이러한 방법들이 효율적이긴 하지만, 주요 기여는 이론적 통합과 견고한 반복 솔버의 개발에 있으며, 향후 희소 행렬(sparse matrices) 및 드라질(Drazin) 역행렬에 대한 연구를 제언한다.