[논문 리뷰] Algebraic properties of robust Pad\\'e approximants

Beckermann, Bernhard, Ana C. Matos|arXiv (Cornell University)|2013. 10. 09.

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

한 줄 요약

이 논문은 구조적 행렬—특히 직사각형 토플리츠 행렬 $ C $, 무늬가 있는 토플리츠 행렬 $ T $, 삼각형 토플리츠 행렬 $ Q $, 그리고 실베스터 유사 행렬 $ S $—를 분석하여 잘 조건화된 패드 근사에서 전방 안정성과 임의의 가짜 극점의 부재를 확립한다. 주요 기여는 조건수 $ \kappa(S) $ 가 적당한 수준인 잘 조건화된 근사에서 프로이사르 듀블릿과 작은 잔여항에 대한 강건성을 보장함으로써 단위 원판 내 수렴의 신뢰성을 확보한다는 점이다.

ABSTRACT

For a recent new numerical method for computing so-called robust Pad\\'e approximants through SVD techniques, the authors gave numerical evidence that such approximants are insensitive to perturbations in the data, and do not have so-called spurious poles, that is, poles with a close-by zero or poles with small residuals. A black box procedure for eliminating spurious poles would have a major impact on the convergence theory of Pad\\'e approximants since it is known that convergence in capacity plus absence of poles in some domain D implies locally uniform convergence in D. In the present paper we provide a proof for forward stability (or robustness), and show absence of spurious poles for the subclass of so-called well-conditioned Pad\\'e approximants. We also give a numerical example of some robust Pad\\'e approximant which has spurious poles, and discuss related questions. It turns out that it is not sufficient to discuss only linear algebra properties of the underlying rectangular Toeplitz matrix, since in our results other matrices like Sylvester matrices also occur. These types of matrices have been used before in numerical greatest common divisor computations.

연구 동기 및 목표

• 해석적 계속성과 함수 근사에서 수렴을 방해하는 오랜 동안의 문제인 가짜 극점과 프로이사르 듀블릿을 해결하기 위해.
• 특히 데이터 편향이나 노이즈 존재 시에도 관찰되는 SVD 기반 패드 근사 방법의 수치적 강건성에 대한 이론적 기반을 제공하기 위해.
• 강건한 패드 근사와 잘 조건화된 패드 근사 간의 차이를 명확히 하여, 오직 잘 조건화된 근사만이 가짜 극점과 작은 잔여항의 부재를 보장한다는 것을 보여주기 위해.
• 특히 실베스터 유사 행렬의 조건수 $ \kappa(S) $ 와 패드 근사의 수치적 안정성 및 수렴 성질 간의 연결 고리를 설정하기 위해.
• 강건한 패드 근사가 가짜 극점이 없는가라는 열린 질문을 다루며, 강건성만으로는 부족하지만, 잘 조건화됨으로써 그 부재가 보장된다는 것을 보여주기 위해.

제안 방법

• 계수 계산을 위한 선형계를 표현하는 데 사용되는 네 가지 구조적 행렬—$ C $ (직사각형 토플리츠), $ T $ (무늬가 있는 토플리츠), $ Q $ (삼각형 토플리츠), $ S = QT $ (실베스터 유사)—을 통한 패드 근사 문제 분석.
• 특이값 분해(SVD)와 의사역행렬 분석을 활용하여, 특히 $ \kappa(S) $ 를 수치적 안정성의 척도로 삼아 기저 선형계의 조건수를 연구.
• 적당한 조건수 $ \kappa(S) $ 를 가진 잘 조건화된 패드 근사의 개념을 도입하여, 역방향 안정성과 프로이사르 듀블릿, 작은 잔여항에 대한 제어를 보장.
• 균일한 현수선 거리(metric)와 정규화된 계수 벡터 차이를 사용하여 근사치를 비교함으로써, 단위 원판 내 진짜 함수에의 가까이움을 측정하는 데 등가성을 입증.
• SVD 기반 방법의 대체로 재귀적 방법인 카바이-멜레시코 알고리즘을 적용하여 강건한 근사를 계산할 수 있음을 보여주지만, 일부 잘 조건화된 근사를 놓칠 수 있음을 시사.
• 수치 예제를 통해 강건한 근사에서도 여전히 가짜 극점이 존재할 수 있음을 보여주며, 강건성만으로는 부족함을 입증함—즉, $ S $ 의 잘 조건화됨이 가짜 잔여항의 부재를 보장하기 위해 필수적임.

실험 결과

연구 질문

• RQ1 수치적 증거가 시사하는 것처럼, 강건한 패드 근사가 진정으로 가짜 극점이 없는가?
• RQ2 실베스터 유사 행렬 $ S $ 의 조건수와 패드 근사에서 프로이사르 듀블릿 또는 작은 잔여항의 부재 사이의 정확한 관계는 무엇인가?
• RQ3 패드 근사에서 전방 안정성과 역방향 안정성이 동시에 달성될 수 있는가? 어떤 행렬 조건수 가정 하에?
• RQ4 조건수 $ \kappa(S) \lesssim 1 $ 으로 정의되는 잘 조건화된 패드 근사의 집합은 강건한 근사의 진부분집합이며, 가짜 극점 없이 용량 수렴을 보장하는가?
• RQ5 구조적 조건수 $ \kappa(S) $ 와 비구조적 조건수 $ \kappa(T) $ 는 어떤 관계가 있으며, 구조적 조건수는 패드 근사 오차에 대해 더 날카로운 경계를 제공할 수 있는가?

주요 결과

• 적당한 $ \kappa(S) $ 를 가진 잘 조건화된 패드 근사는 전방 잘 조건화되어 있으며, 가짜 극점이나 프로이사르 듀블릿이 존재하지 않는다.
• 잘 조건화된 근사에 대해 실수 계수 패드 사상은 역방향 잘 조건화되어 있으며, $ \|Q^{-1}T\| \lesssim \kappa(S) $ 를 만족하여 수치적 안정성을 보장한다.
• 균일한 현수선 거리와 정규화된 계수 벡터 차이 간의 거리는 단위 원판 내 진짜 함수에의 가까이움을 측정하는 데 등가적인 메트릭이다.
• 잘 조건화된 근사에서는 분자와 분모 다항식의 수치적 서로소성이 보장되며, 이는 [4, 7] 의 기준과 일치한다.
• 강건한 근사에서도 여전히 가짜 극점이 존재할 수 있으며, 이는 2.3 예제에서 입증되었으며, 강건성만으로는 프로이사르 듀블릿을 제거할 수 없음을 보여준다.
• 코로나리에 6.3과 정리 1.4에 따르면, 고정된 대각선을 따라 진짜 함수 $ f $ 와 가까운 잘 조건화된 패드 근사는 유한 개 이하로 제한된다.