[논문 리뷰] Approximations for the Caputo derivative (II)

Yuri Dimitrov, Venelin Todorov외 1명|arXiv (Cornell University)|2016. 04. 25.

Fractional Differential Equations Solutions참고 문헌 20인용 수 6

한 줄 요약

이 논문은 다항로그함수 전개와 중점 구적법을 사용하여 캡투오 분수도수 도함수의 고차수 수치 근사법을 개발하며, 베타 함수와 거듭제곱 함수의 캡투오 도함수에 대한 리만 합 근사의 점근적 오차 전개를 유도한다. 주요 기여는 수렴 속도가 향상된 새로운 유형의 이동형 그룬발트 유형 근사법을 제안한 것으로, 분수차수 2−α 및 2−2α를 가진 분수적 미분방정식의 수치적 해를 통해 검증되었으며, 테스트 케이스에서 두 번째 차수 정확도를 달성하였다.

ABSTRACT

In the present paper we use the expansion formula of the polylogarithm function to construct approximations of the Caputo derivative which are related to the midpoint approximation of the integral in the definition of the Caputo derivative. The asymptotic expansion formula of the Riemann sum approximation of the beta function and the first terms of the expansion formulas of the approximations of the Caputo derivative of the power function are obtained in the paper. The induced shifted approximations of the Grünwald formula and the approximations of the Caputo derivative studied in the first part of the paper are constructed and applied for numerical solution of fractional differential equations.

연구 동기 및 목표

• 표준 L1 스킴을 초월하여 캡투오 분수도수 도함수에 대한 고차수 수치 근사법을 개발하기 위해.
• 이동형 그룬발트 유형 공식을 구성하여 분수적 미분방정식을 해결할 때 수렴 속도를 향상시키기 위해.
• 베타 함수와 거듭제곱 함수의 캡투오 도함수에 대한 리만 합 근사의 점근적 전개 공식을 도출하기 위해.
• 알려진 해석적 해를 가진 두 항 및 세 항 분수적 미분방정식의 수치적 해를 통해 제안된 근사법을 검증하기 위해.
• 최대 l∞ 노름을 사용하여 신규 근사법의 수렴 행동과 오차 차수를 분석하기 위해.

제안 방법

• 적분 정의에 기반한 중점 구적법을 활용하여 다항로그함수 함수 전개를 사용해 캡투오 도함수의 근사법을 구성한다.
• 베타 함수의 리만 합 근사에 대한 점근적 전개 공식을 도출하여 수치적 방법의 오차 분석을 가능하게 한다.
• 기본 가중치를 수정하여 그룬발트 공식의 유도된 이동형 근사법을 구성함으로써 수렴 차수를 향상시킨다.
• 유도된 근사법을 두 항 및 세 항 분수적 미분방정식의 수치적 해법에 적용하며, 안정성과 수렴성 분석을 포함한 반복적 방법을 사용한다.
• 오일러-마클라우린 공식과 제타 함수 항등식을 활용하여 해의 도함수와 베르누이 수를 포함하는 전개 항을 유도한다.
• 알려진 해석적 형태를 가진 분수적 미분방정식의 정확한 해를 사용하여 방법을 검증하며, 최대 l∞ 오차와 수렴 차수를 계산한다.

실험 결과

연구 질문

• RQ1 다항로그함수 기반 전개는 표준 L1 방법에 비해 캡투오 도함수에 대해 더 고차수 근사법을 제공할 수 있는가?
• RQ2 베타 함수와 거듭제곱 함수의 캡투오 도함수에 대한 리만 합 근사의 점근적 오차 구조는 무엇인가?
• RQ3 이동형 그룬발트 유형 근사법은 분수적 미분방정식을 해결할 때 수렴 속도를 어떻게 향상시키는가?
• RQ4 두 항 및 세 항 분수적 미분방정식에 적용했을 때 신규 수치적 방법의 관측된 수렴 차수는 무엇인가?
• RQ5 유도된 전개 공식은 분수도수 도함수 근사의 맥락에서 제타 함수와 베르누이 수와 어떻게 관련되는가?

주요 결과

• 제안된 이동형 근사법은 두 항 분수적 미분방정식에서 두 번째 차수 수렴(2−α)을 달성하며, 수치 실험에서 관측된 수렴 차수는 1.99에 가까웠다.
• α=0.25인 두 항 방정식의 경우, h=0.003125에서 최대 오차는 0.7875×10⁻⁶였으며, 수렴 차수는 1.9998였다.
• α=0.4인 세 항 방정식의 수치적 해는 h=0.003125에서 수렴 차수 1.1727을 기록하여 개선된 안정성과 정확도를 보였다.
• α=0.6일 때 두 번째 차수 방법은 수렴 차수 1.9292를 달성하여 이론적 2−α 수렴 속도를 확인하였다.
• 오차 분석 결과, 리만 제타 함수와 해의 이阶 도함수를 포함하는 점근적 전개 항이 주요 오차 항에 크게 기여하는 것으로 나타났다.
• 다항로그함수 및 제타 함수 항등식에서 유도된 고차수 보정 항을 포함함으로써, 이 방법은 L1 근사 프레임워크를 성공적으로 확장하였다.