[논문 리뷰] Backward errors and linearizations for palindromic matrix polynomials

Bibhas Adhikari|ArXiv.org|2008. 12. 22.

Matrix Theory and Algorithms참고 문헌 26인용 수 3

한 줄 요약

이 논문은 $*$-팔린드롬 및 $*$-반팔린드롬 행렬 다항식의 근사 고유요소에 대한 계산 가능한 구조적 후행 오차 표현식을 유도하고, 근사 고유쌍을 정확하게 만들기 위한 최소 구조 유지 편집을 특성화한다. 고유값이 단위원 위에 있을 경우, 구조적 선형화가 후행 오차에 거의 영향을 주지 않음을 보이며, 프로베니우스 노름과 스펙트럴 노름에 대해 각각 $\sqrt{2}$ 및 $\sqrt{(m+1)/(2m)}$의 상한을 제공한다.

ABSTRACT

We derive computable expressions of structured backward errors of approximate eigenelements of *-palindromic and *-anti-palindromic matrix polynomials. We also characterize minimal structured perturbations such that approximate eigenelements are exact eigenelements of the perturbed polynomials. We detect structure preserving linearizations which have almost no adverse effect on the structured backward errors of approximate eigenelements of the *-palindromic and *-anti-palindromic polynomials.

연구 동기 및 목표

• 모든 $*$-팔린드롬 및 $*$-반팔린드롬 행렬 다항식에서 근사 고유요소에 대한 계산 가능한 구조적 후행 오차 표현식을 도출한다.
• 편집된 다항식에서 근사 고유요소를 정확하게 만들기 위한 최소 구조 유지 편집을 특성화한다.
• 근사 고유요소의 구조적 후행 오차에 대해 최소한의 영향을 미치는 구조 유지 선형화를 규명한다.
• 특히 고유값이 단위원 위에 있을 경우 선형화가 후행 오차에 미치는 영향을 분석한다.
• 이러한 구조적 행렬 다항식에 대해 구조 유지 알고리즘을 통해 계산된 고유요소의 정확도 평가를 지원한다.

제안 방법

• 구조 유지 편집 $\triangle\mathrm{P}$ 의 최소 노름을 통해 구조적 후행 오차 $\eta_M^\mathbb{S}(\lambda, x, \mathrm{P})$ 를 유도하며, 이는 $\mathrm{P}(\lambda)x + \triangle\mathrm{P}(\lambda)x = 0$ 을 만족한다.
• 편집의 크기를 측정하기 위해 $\mathbb{C}^{n \times n}$ 상에서 프로베니우스 노름과 스펙트럴 노름을 사용하여 $|\!|\!|\cdot|\!|\!|_M$ 을 정의한다.
• 최소 후행 오차를 달성하는 $\triangle\mathrm{P}$ 의 명시적 표현을 구성하며, 프로베니우스 노름의 경우 유일성, 스펙트럴 노름의 경우 비유일성을 보여준다.
• 후행 오차가 원래 다항식 $\mathrm{P}$ 와 그 선형화 $\mathrm{L}$ 간의 관계를 분석하며, 안사 벡터 $v$ 와 $\Lambda_{m-1} = [\lambda^{m-1}, \dots, \lambda, 1]^T$ 를 사용한다.
• 예를 들어 [3] 에서 제시된 구조적 후행 오차 분석 이론을 $*$-팔린드롬 및 $*$-반팔린드롬 구조에 적용한다.
• 잔차 $r = -\mathrm{P}(\lambda)x$ 와 그 노름을 사용하여 선형화된 형태와 원래 형태 간의 구조적 후행 오차 비율에 대한 상한을 도출한다.

실험 결과

연구 질문

• RQ1 모든 $*$-팔린드롬 또는 $*$-반팔린드롬 행렬 다항식에서 근사 고유요소에 대한 계산 가능한 구조적 후행 오차 표현식은 무엇인가?
• RQ2 편집된 다항식에서 근사 고유요소를 정확하게 만들기 위한 최소 구조 유지 편집은 무엇인가?
• RQ3 원래 다항식의 구조적 후행 오차와 그 구조 유지 선형화의 오차는 어떻게 비교되는가?
• RQ4 고유값이 단위원 위에 있을 경우와 아닐 경우 선형화가 후행 오차에 미치는 영향은 어떠한가?
• RQ5 $*$-팔린드롬/반팔린드롬 다항식에 대해 후행 오차에 대한 악영향을 최소화하는 구조적 선형화는 무엇인가?

주요 결과

• 프로베니우스 노름의 경우 최소 편집 $\triangle\mathrm{P}$ 는 유일하며, 스펙트럴 노름의 경우 이러한 편집은 무수히 많다.
• 만약 $|\lambda| = 1$ 이면, 선형화의 구조적 후행 오차와 원래 다항식의 오차 비율은 프로베니우스 노름과 스펙트럴 노름 모두에서 $\sqrt{(m+1)/(2m)}$ 와 $\sqrt{2}$ 사이로 유계이다.
• $|\lambda| \neq 1$ 인 경우, 후행 오차 비율은 잔차 $r = -\mathrm{P}(\lambda)x$ 와 선형화의 구조에 따라 달라지며, $\|\widehat{r}\|_2$ 와 $\|(1, \lambda)^T\|_2$ 를 포함하는 명시적 상한이 존재한다.
• 고유값이 단위원 위에 있을 경우, 선형화된 시스템 $\mathrm{L}$ 의 후행 오차는 $\eta_F^\mathbb{S}(\lambda, \Lambda_{m-1} \otimes x; \mathrm{L}, v) \leq \sqrt{2} \cdot \eta_F(\lambda, x, \mathrm{P})$ 를 만족한다.
• 단위원 외부에 있는 고유값의 경우, 각각 $\|r_p\|_2$ 또는 $\|r_{ap}\|_2$ 를 최소화함으로써 $H$-팔린드롬 또는 $H$-반팔린드롬 선형화를 선택하는 것이 바람직하다.
• 분석 결과, 고유값이 단위원 위에 있을 경우 구조 유지 선형화가 후행 오차에 거의 영향을 주지 않음을 확인하였으며, 이는 정확도 유지 계산에 이상적임을 시사한다.