[논문 리뷰] Third Order Newton's Method for Zernike Polynomial Zeros

Richard J. Mathar|arXiv (Cornell University)|2007. 05. 09.

Adaptive optics and wavefront sensing참고 문헌 23인용 수 3

한 줄 요약

이 논문은 제르니케 반경 다항식의 20차 이하 양의 실근을 효율적으로 계산하기 위해 제3차 수렴을 달성하는 제3차 뉴턴 방법을 제시한다. 이는 초함수 표현, 도함수 비율을 위한 종결 연속분수, 그리고 재귀적 평가를 활용하여 다항식을 직접 평가하지 않고도 수행된다. 주요 기여는 고정밀도 근 표를 포함한 PARI/GP 구현체이다. 이는 광학적 웨이브플런 분석에 활용된다.

ABSTRACT

The Zernike radial polynomials are a system of orthogonal polynomials over the unit interval with weight x. They are used as basis functions in optics to expand fields over the cross section of circular pupils. To calculate the roots of Zernike polynomials, we optimize the generic iterative numerical Newton's Method that iterates on zeros of functions with third order convergence. The technique is based on rewriting the polynomials as Gauss Hypergeometric Functions, reduction of second order derivatives to first order derivatives, and evaluation of some ratios of derivatives by terminating continued fractions. A PARI program and a short table of zeros complete up to polynomials of 40th order are included.

연구 동기 및 목표

• 광학적 웨이브플런 분석에 필수적인 제르니케 반경 다항식의 양의 실근을 빠르고 고차수의 반복 방법으로 계산하는 것.
• 제르니케 다항식과 그 도함수의 직접 평가가 필요 없도록, f/f' 및 f''/f'의 재귀적 연속분수 표현을 사용하여 이를 제거하는 것.
• 제르니케 다항식이 만족하는 구조 관계와 미분방정식을 활용하여 근 찾기 과정에서 제3차 수렴을 달성하는 것.
• 광학 및 간섭계 응용에 적합한 수치적으로 안정적이고 고정밀도 알고리즘을 제공하는 것.
• 웨이브플런 복원 및 피 pupils 분석에 활용하기 위한, 최대 20차까지의 제르니케 다항식 근에 대한 포괄적인 표를 생성하는 것.

제안 방법

• 제르니케 다항식을 가우스 초함수로 재구성하여 그 미분 구조에 접근한다.
• 연속된 제르니케 다항식 간의 구조 관계를 이용해 f/f'를 종결 연속분수로 표현함으로써, f 또는 그 도함수의 직접 평가를 피한다.
• 두 번째 및 세 번째 도함수는 미분방정식을 통해 f/f'와 관련지어지며, f''/f' 및 f'''/f'의 평가가 명시적 도함수 계산 없이도 가능해진다.
• 유도된 비율을 사용하여 반복적으로 제3차 뉴턴 업데이트 공식 Δx = -(f/f') / (1 - 0.5*(f/f')*(f''/f')) 을 적용한다.
• 초기 근 추정치는 해석적 근사치로 계산되며, n=m+2인 경우 정확한 값, 그 외에는 유리함수 히우리스틱을 사용한다.
• 고차수 근의 초기 추정치를 향상시키기 위해 타일러 전개 기반의 이차 근사가 사용되며, 이후 뉴턴 반복이 적용된다.

실험 결과

연구 질문

• RQ1 제르니케 다항식의 근 찾기 과정에서 다항식이나 도함수를 직접 평가하지 않고도 제3차 수렴을 달성할 수 있는가?
• RQ2 다만 재귀 관계와 연속분수를 사용하여 f/f' 및 f''/f'의 비율을 효율적이고 안정적으로 계산할 수 있는가?
• RQ3 단일 제르니케 다항식의 다수의 근에 대해 반복적 근 찾기 과정을 초기화하는 데 최적의 전략은 무엇인가?
• RQ4 제르니케 다항식이 만족하는 미분방정식의 구조는 계산 비용을 얼마나 줄일 수 있는가?
• RQ5 실제 광학 시스템 설계 및 웨이브플런 센싱 응용에서 사용할 수 있을 정도로 결과 근 표의 정확성과 신뢰성은 어느 정도인가?

주요 결과

• 이 방법은 다항식의 도함수를 직접 평가하지 않고도 재귀 관계와 종결 연속분수를 사용하여 제르니케 다항식의 근 찾기에서 제3차 수렴을 달성한다.
• f/f'의 계산은 종결 연속분수로 단순화되어 호너 방법을 피하고 수치적 취소 현상을 최소화한다.
• 두 번째 도함수 비율 f''/f'는 미분방정식과 f/f' 비율에서 유도되며, 저비용으로 고차수 수렴을 가능하게 한다.
• 알고리즘은 최대 20차까지의 제르니케 다항식 R_n^m(x)의 (0,1) 내 (n−m)/2개의 양의 실근을 30자리 정밀도로 성공적으로 계산한다.
• PARI/GP 구현체가 제공되며, n≤20 및 m≤n, n−m가 짝수인 모든 (n,m) 쌍에 대해 전체 근 표를 출력하는 주 routine이 포함되어 있다.
• 이 방법은 수치적으로 안정적이며 고차수의 직교 다항식을 직접 평가할 때 발생하는 정밀도 손실을 피한다.