[논문 리뷰] Block Kronecker Linearizations of Matrix Polynomials and their Backward Errors

Froilán M. Dopico, Piers W. Lawrence외 2명|MIMS EPrints (University of Southampton)|2017. 07. 16.

Matrix Theory and Algorithms참고 문헌 45인용 수 5

한 줄 요약

이 논문은 행렬 다항식을 위한 새로운 강한 선형화 가족인 블록 크로네커 선형화를 소개하며, 완전한 다항식 고유값 문제의 역으로 안정적인 해법을 가능하게 한다. 유한한 편미분에 대해 엄밀한 유한 편차 후방 오차 분석을 제공하여, 계산된 고유구조가 약간 변형된 다항식에 대해 정확하게 되며, 넓은 선형화 클래스에 대해 유효한 정확한 경계를 제시한다.

ABSTRACT

We introduce a new family of strong linearizations of matrix polynomials---which we call "block Kronecker pencils"---and perform a backward stability analysis of complete polynomial eigenproblems. These problems are solved by applying any backward stable algorithm to a block Kronecker pencil, such as the staircase algorithm for singular pencils or the QZ algorithm for regular pencils. This stability analysis allows us to identify those block Kronecker pencils that yield a computed complete eigenstructure which is exactly that of a slightly perturbed matrix polynomial. The global backward error analysis in this work presents for the first time the following key properties: it is a rigurous analysis valid for finite perturbations (i.e., it is not a first order analysis), it provides precise bounds, it is valid simultaneously for a large class of linearizations, and it establishes a framework that may be generalized to other classes of linearizations. These features are related to the fact that block Kronecker pencils are a particular case of the new family of "strong block minimal bases pencils", which are robust under certain perturbations and, so, include certain perturbations of block Kronecker pencils. We hope that this robustness property will allow us to extend the results in this paper to other contexts.

연구 동기 및 목표

• 완전한 다항식 고유값 문제를 해결하는 데 있어 후방 안정성을 보장하는 행렬 다항식을 위한 새로운 강한 선형화 클래스를 개발하는 것.
• 일차 근사가 아닌 유한한 편미분에 대해 유효한 엄밀한 후방 오차 분석을 수행하는 것.
• 블록 크로네커 펜슬을 포함한 넓은 가족의 선형화에 동시에 적용 가능한 정확한 후방 오차 경계를 설정하는 것.
• 계산된 고유구조가 약간 변형된 행렬 다항식의 고유구조와 정확히 일치하는 조건을 규명하는 것.
• 강한 블록 최소 기저 펜슬의 프레임워크를 통해 다른 선형화 유형으로의 분석 일반화를 위한 기반을 마련하는 것.

제안 방법

• 이중 최소 기저와 블록 최소 기저의 구조에서 유도된 새로운 강한 선형화 가족인 블록 크로네커 펜슬을 제안한다.
• 계수 행렬을 블록 대각 형식으로 변환하기 위해 직교 변환과 행/열 순열을 적용한다.
• 구조적 행렬 $\hat{C}_{\varepsilon-1}(L_{\varepsilon}(\lambda))$ 와 $\hat{C}_{\varepsilon}(L_{\varepsilon}(\lambda))$ 의 특이값 분석을 통해 후방 오차를 경계한다.
• 이 계수 행렬의 최소 특이값이 $\sigma_{\min}(M_\varepsilon) = 2\sin\left(\frac{\pi}{4\varepsilon+2}\right)$ 라고 나타나며, 이는 후방 오차 경계를 결정한다.
• 변환된 행렬의 직교성과 동치성을 활용하여 특이값을 유지하고 정밀한 오차 정량화를 가능하게 한다.
• 강한 블록 최소 기저 펜슬의 프레임워크를 활용하여 편미분에 대한 강건성을 확보하고, 다른 선형화 가족으로의 일반화를 가능하게 한다.

실험 결과

연구 질문

• RQ1 유한한 편미분에 대해 유효하고 일차 근사가 아닌 엄밀한 후방 오차 분석을 어떻게 개발할 수 있는가?
• RQ2 계산된 고유구조가 약간 변형된 행렬 다항식의 고유구조와 정확히 일치하도록 하기 위해 선형화가 가져야 할 조건은 무엇인가?
• RQ3 블록 크로네커 펜슬은 다른 선형화들과 비교해 후방 오차 경계와 안정성 측면에서 어떻게 다를까?
• RQ4 후방 오차 분석을 블록 크로네커 펜슬을 초월한 다른 선형화 유형으로 일반화할 수 있는가?
• RQ5 최소 인덱스와 이중 최소 기저는 강건하고 후방 안정적인 선형화를 구성하는 데 어떤 역할을 하는가?

주요 결과

• 후방 오차 분석은 엄밀하고 유한한 편미분에 대해 유효하며, 일차 근사에 국한되지 않는다.
• 핵심 행렬 $M_\varepsilon$ 의 최소 특이값은 $\sigma_{\min}(M_\varepsilon) = 2\sin\left(\frac{\pi}{4\varepsilon+2}\right)$ 이며, 이는 후방 오차 경계를 결정한다.
• 이 최소 특이값은 $\sigma_{\min}(M_\varepsilon) \geq \frac{3}{2(\varepsilon+1)}$ 를 만족하여 오차 분석에 유용한 하한을 제공한다.
• 분석 결과, 블록 크로네커 펜슬의 계산된 고유구조는 후방 오차 크기로 제한된 약간 변형된 행렬 다항식의 고유구조와 정확히 일치함을 보여준다.
• 강한 블록 최소 기저 펜슬의 프레임워크는 편미분에 대한 강건성을 보장하며, 결과를 다른 선형화 가족으로 확장하는 데 기여한다.
• 계수 행렬 $C_{\varepsilon-1}(L_{\varepsilon}(\lambda))$ 와 $C_{\varepsilon}(L_{\varepsilon}(\lambda))$ 의 특이값은 직교 변환에 의해 유지되어 정밀한 오차 정량화가 가능하다.