[논문 리뷰] A new algorithm for computing Pade approximants

В. М. Адуков, Ibryaeva O. L|arXiv (Cornell University)|2011. 12. 24.

Statistical and numerical algorithms참고 문헌 7인용 수 6

한 줄 요약

이 논문은 토플리츠 행렬의 커널 구조에서 유도된 최소 차수 분모를 기반으로 한 축소 패드 근사법을 정의하여 패드 근사의 새로운 알고리즘을 제안한다. 이 방법은 수치 반올림 오차로 인한 가짜 다항식 근(폴-제로 쌍)인 프로이산 이중극을 피하기 위해 분모가 차수에서 유일하게 최소가 되도록 보장함으로써, 비정상적인 경우에도 안정적이고 정확한 유리근사가 가능하다.

ABSTRACT

In paper a new definition of reduced Pade approximant and algorithm for its computing is proposed. Our approach is based on the investigation of the kernel structure of the Toeplitz matrix. It is shown that the reduced Pade approximant always has nice properties which classical Pade approximant possesses only in the normal case. The new algorithm allows us to avoid Froissart doublets appearance induced by computer roundoff in the non-normal Pade table.

연구 동기 및 목표

• 신호 처리 및 기타 응용 분야에서 수치 반올림 오차로 인해 발생하는 가짜 다항식 근(폴-제로 쌍)인 프로이산 이중극 문제를 해결하기 위해.
• 비정상적인 경우에 비유일성과 불안정성으로 악화되는 고전적 패드 근사의 대안으로 강건한 방법을 개발하기 위해.
• 결과로 얻어진 유리근사가 정상적인 경우에만 보장되던 정확한 차수와 결함 지수를 유지하도록 보장하기 위해.
• 패워 시리즈 계수와 관련된 토플리츠 행렬의 커널에서 직접 최소 차수 분모를 계산할 수 있는 수치적으로 안정적인 알고리즘을 제공하기 위해.

제안 방법

• 패워 시리즈 계수와 관련된 토플리츠 행렬의 커널에서 유도된 가능한 모든 분모 중에서 차수를 최소화한 유리함수를 축소 패드 근사로 정의한다.
• 토플리츠 행렬의 커널 구조를 이용해 가능한 모든 분모를 매개변수화하고, 랭크 분석을 통해 최소 차수의 분모를 식별한다.
• 단계 3.2 정리를 적용하여, 잘라낸 토플리츠 부분행렬의 랭크에서 유도된 본질적 지수 μ₁에 기반해 첫 번째 본질적 다항식(최소 차수 분모)을 계산한다.
• 단계 4.2 정리와 단계 4.3 정리를 활용해 분모와 분자의 가짜 영계수를 탐지하고 제거하여, 최종 유리함수가 기약 형태임을 보장한다.
• d가 계수 삭제를 제어하는 프로시저 `ReducedPade(f, m, n, a, d)`를 사용해 맵라에서 알고리즘을 구현하고, 근 찾기 및 정확한 유리함수와의 비교를 통해 결과를 검증한다.
• 수치 실험을 통해 축소 근사가 원래 함수와 정확히 일치함(기계 정밀도 범위 내)을 보이며, 가짜 근이 발생하지 않음을 입증한다.

실험 결과

연구 질문

• RQ1 비정상적인 경우에도 분모가 차수에서 유일하게 최소가 되도록 패드 근사를 정의할 수 있는가?
• RQ2 어떻게 토플리츠 행렬의 커널 구조를 활용하여 수치적 오염 요소 없이 최소 차수 분모를 계산할 수 있는가?
• RQ3 실제 계산에서 단일 부동소수점 반올림 오차로 인한 프로이산 이중극의 출현을 얼마나 효과적으로 방지할 수 있는가?
• RQ4 고전적 패드 근사가 실패할 경우에도 축소 패드 근사가 유지하는 성질(예: 차수 유일성, 공통 근이 없는 성질)은 무엇인가?
• RQ5 맵라와 같은 기호 계산 시스템에서 정확성과 가짜 다항식 근-제로 쌍을 피하면서도 효율적으로 알고리즘을 구현할 수 있는가?

주요 결과

• 축소 패드 근사는 상수 인자 외에는 유일하게 정의되며, 분자와 분모의 정확한 차수와 알려진 결함 지수를 가진다.
• 알고리즘은 토플리츠 행렬의 커널을 분석하여 최소 차수 분모를 성공적으로 계산함으로써, 비유일성 또는 가짜 해를 유도하지 않는다.
• 예제 5.1에서 $ f(z) = \frac{(z+1)(z-2)}{(z+2.1)(z-1)} $에 대한 (4,4) 근사에서 정확한 극(폴) $ z = -2.1, 1 $과 정확한 영(제로) $ z = -1, 2 $를 복원하며, 가짜 근이 없음을 확인하였다.
• 예제 5.2에서 $ f(z) = \frac{z+1.01}{(z+2)(z-2.01)} $에 대한 (2,3) 근사에서 정확한 영 $ z = -1.01 $과 극 $ z = -2, 2.01 $을 복원하였으며, 프로이산 이중극이 없음을 확인하였다.
• 예제 5.3에서 $ d = 1 $을 설정해 영계수를 유지하면 가짜 근(예: $ \sim 10^5 $)가 발생하지만, $ d = 0 $으로 설정해 영계수를 삭제하면 정확한 극과 영이 복원됨을 통해 계수 삭제의 필요성을 입증하였다.
• d = 0일 경우 알고리즘이 원래 유리함수의 근 구조와 정확히 일치하는 근사치를 일관되게 생성함으로써, 반올림 오차로 인한 이중극에 대한 강건성을 입증하였다.