[논문 리뷰] A Tunably-Accurate Laguerre Petrov-Galerkin Spectral Method for Multi-Term Fractional Differential Equations on the Half Line

Anna Lischke, Mohsen Zayernouri외 1명|arXiv (Cornell University)|2016. 07. 28.

Fractional Differential Equations Solutions참고 문헌 15인용 수 4

한 줄 요약

이 논문은 반직선에서 최대 1인 도함수 차수를 가진 다항수 분수미분방정식을 해결하기 위한 조절 가능한 정확도를 가진 라거르 페트로프-갈레르킨 스펙트럴 방법을 제안한다. 일반화된 관련 라거르 다항식(GALFs)을 기저 함수로 사용하고 재귀관계를 활용함으로써, 스펙트럴 수렴을 달성하며 잘 조절된 희소 행렬을 사용해 O(N log N) 연산으로 선형 시스템을 해소한다.

ABSTRACT

We present a new tunably-accurate Laguerre Petrov-Galerkin spectral method for solving linear multi-term fractional initial value problems with derivative orders at most one and constant coefficients on the half line. Our method results in a matrix equation of special structure which can be solved in $\mathcal{O}(N \log N)$ operations. We also take advantage of recurrence relations for the generalized associated Laguerre functions (GALFs) in order to derive explicit expressions for the entries of the stiffness and mass matrices, which can be factored into the product of a diagonal matrix and a lower-triangular Toeplitz matrix. The resulting spectral method is efficient for solving multi-term fractional differential equations with arbitrarily many terms. We apply this method to a distributed order differential equation, which is approximated by linear multi-term equations through the Gauss-Legendre quadrature rule. We provide numerical examples demonstrating the spectral convergence and linear complexity of the method.

연구 동기 및 목표

• 분수도함수의 국소적이지 않은 성질로 인해 발생하는 다항수 분수미분방정식 해법의 높은 계산 비용을 해결하기 위해.
• 반직선에서 도함수 차수 ≤1인 선형 다항수 분수초기값문제(FIVPs)를 위한 효율적인 스펙트럴 방법을 개발하기 위해.
• 가우스-레지앙드르 적분을 통해 분포형 도함수 미분방정식을 근사하고, 이를 바탕으로 다항수 FIVPs를 해결함으로써 분포형 도함수 미분방정식에 대해 고정밀도와 고효율을 달성하기 위해.
• 수렴 속도를 제어할 수 있는 조절 가능한 매개변수를 도입하여 방법의 정확도를 조절 가능하게 하기 위해.
• GALFs를 기저 함수로 사용하는 페트로프-갈레르킨 설정을 통해 잘 조절된 희소 선형 시스템을 달성하기 위해.

제안 방법

• 반직선에서의 분수 스투름-리우빌 문제의 고유함수로 유도된 일반화된 관련 라거르 다항식(GALFs)을 시험함수와 시험함수로 사용하기 위해.
• 분수 적분 부분적 통합을 통해 시험함수로의 분수도함수 전달이 가능한 영리한 페트로프-갈레르킨 변분법을 설정하기 위해.
• 라거르 다항식의 재귀관계를 활용해 강성행렬과 질량행렬의 항들에 대한 명시적 표현을 유도하여 수치적 적분을 피하기 위해.
• 행렬 항들을 대각행렬과 하삼각 토플리츠 행렬의 곱으로 분해하여 빠른 행렬-벡터 곱셈을 가능하게 하기 위해.
• 구조화된 행렬 성질을 활용하여 빠른 변환을 적용해 O(N log N) 연산으로 최종 선형 시스템을 해소하기 위해.
• 분포형 도함수 미분방정식을 다항수 FIVPs로 근사하기 위해 가우스-레지앙드르 적분을 적용하여, 트라페즈 방식 대비 수렴성을 향상시키기 위해.

실험 결과

연구 질문

• RQ1 도함수 차수 ≤1인 반직선에서의 다항수 분수초기값문제에 대해, 높은 정확도와 계산 효율성을 동시에 확보하는 페트로프-갈레르킨 스펙트럴 방법을 구성할 수 있는가?
• RQ2 기저 함수는 어떻게 선택해야 하며, t^α 형태의 해의 특이성을 자연스럽게 포착하고 스펙트럴 수렴을 보장할 수 있는가?
• RQ3 수렴 속도를 제어할 수 있는 매개변수를 도입함으로써 방법의 정확도를 조절 가능한가?
• RQ4 결과로 생기는 선형 시스템의 조립 및 해소의 계산 복잡도는 얼마이며, 행렬의 구조적 특성을 활용해 O(N²) 이하로 줄일 수 있는가?
• RQ5 분포형 도함수 방정식 근사를 위한 적분 규칙 선택(트라페즈 대비 가우스-레지앙드르)이 전체 정확도에 어떤 영향을 미치는가?

주요 결과

• 방법은 반직선에서 가중 L²-노름에서 스펙트럴 수렴을 달성하며, 이 수렴 속도는 조절 매개변수 α₁에 민감하게 영향을 받는다.
• 강성행렬과 질량행렬은 대각행렬과 하삼각 토플리츠 행렬의 곱으로 표현되며, 이는 O(N log N) 해소 시간을 가능하게 한다.
• 수치적 적분을 피하여 명시적 행렬 항들을 도출함으로써 적분 오차를 방지하고 조립 비용을 감소시킨다.
• 수치 실험 결과, 가우스-레지앙드르 적분은 K=10개의 적분 점으로 기계 정밀도에 도달하는 반면, 트라페즈 법칙은 K=50일 때조차도 약 O(10⁻³) 오차로 정체된다.
• t^λ 특이성을 가진 비연속적인 해에 대해서도, 조절 매개변수 α₁를 특이성 지수와 일치시킬 경우(예: λ = 2 + 1/3일 때 α₁ = 1/3) 높은 정확도를 유지한다.
• 분포형 도함수 방정식은 다항수 FIVPs로 근사함으로써 성공적으로 해소되며, 오차는 K가 충분히 크지 않을 경우에만 적분 오차에 의해 지배된다.