[논문 리뷰] Generalized Jacobi Functions and Their Applications to Fractional Differential Equations

Sheng Chen, Jie Shen외 1명|arXiv (Cornell University)|2014. 07. 31.

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

한 줄 요약

이 논문은 분수계 미분방정식(FDEs)의 특이성과 비국소성 문제를 해결하기 위해 조정 가능한 매개변수를 가진 일반화된 자코비 함수(GJFs)를 도입한다. GJF-페트로프-갈레르킨 스펙트럴 방법을 구축함으로써 저밀도이고 잘 조절된 선형 시스템과 데이터의 부드러움에만 의존하는 스펙트럴 수렴 속도를 달성하여, 경계 특이성이 있는 해에 대해서도 높은 정확도를 확보한다.

ABSTRACT

In this paper, we consider spectral approximation of fractional differential equations (FDEs). A main ingredient of our approach is to define a new class of generalized Jacobi functions (GJFs), which is intrinsically related to fractional calculus, and can serve as natural basis functions for properly designed spectral methods for FDEs. We establish spectral approximation results for these GJFs in weighted Sobolev spaces involving fractional derivatives. We construct efficient GJF-Petrov-Galerkin methods for a class of prototypical fractional initial value problems (FIVPs) and fractional boundary value problems (FBVPs) of general order, and show that with an appropriate choice of the parameters in GJFs, the resulted linear systems can be sparse and well-conditioned. Moreover, we derive error estimates with convergence rate only depending on the smoothness of data, so truly spectral accuracy can be attained if the data are smooth enough. The idea and results presented in this paper will be useful to deal with more general FDEs associated with Riemann-Liouville or Caputo fractional derivatives.

연구 동기 및 목표

• 기존 수치적 방법들이 비효율적인 분수계 미분방정식(FDEs)의 특이해와 비국소 분수도함수 문제를 해결한다.
• 비국소성으로 인해 밀도 높고 잘 조절되지 않는 행렬을 생성하는 유한차분법 및 유한요소법의 한계를 극복한다.
• 분수계 미분방정식 해의 본질적 특이성을 반영하는 기저 함수를 활용해 스펙트럴 프레임워크를 개발한다.
• 해의 정규성에 관계없이 최적의 수렴 속도를 확보하여, 데이터가 부드러울 경우 스펙트럴 정확도를 달성한다.
• 기존에 수치적으로 다루지 못했던 고차수의 초기값 및 경계값 문제를 포함한 일반차수 FDEs에까지 스펙트럴 방법의 적용 범위를 확장한다.

제안 방법

• FDE 해에 내재된 경계 특이성을 반영하기 위해 두 개의 실수 매개변수를 가진 새로운 일반화된 자코비 함수(GJFs) 클래스를 정의한다.
• 리만-리우빌 도함수 및 카프토 도함수의 도함수 공식을 포함한 GJFs의 분수계 미분법 성질을 수립한다.
• 분수도함수를 포함한 가중 치비셰프 공간을 활용해 GJF-페트로프-갈레르킨 스펙트럴 방법을 구성한다.
• 결과로 얻어지는 선형 시스템이 저밀도이고 잘 조절되도록 GJFs의 매개변수를 선택하여 효율적인 해법을 가능하게 한다.
• 해의 정규성에 영향을 받지 않고 데이터의 부드러움에만 의존하는 가중 치비셰프 노름에서의 오차 추정을 유도한다.
• 실수 매개변수를 가진 자코비 다항식의 성질, 특히 도함수 규칙과 재귀관계를 활용해 안정적인 계산을 가능하게 한다.

실험 결과

연구 질문

• RQ1 분수계 미분방정식 해의 특이 행동을 자연스럽게 반영할 수 있는 새로운 기저 함수 클래스를 구성할 수 있는가?
• RQ2 비국소성으로 인해 전통적인 스펙트럴 방법이나 유한차분법보다 저밀도이고 잘 조절된 선형 시스템을 생성할 수 있는 스펙트럴 방법을 어떻게 설계할 수 있는가?
• RQ3 가중 치비셰프 공간에서 분수도함수를 포함한 GJF 기반 스펙트럴 근사의 오차 추정은 무엇인가?
• RQ4 적절한 기저 함수를 사용함으로써 데이터가 부드러운 경우 해가 특이하더라도 스펙트럴 수렴을 달성할 수 있는가?
• RQ5 제안된 GJF 프레임워크는 고차수의 초기값 및 경계값 문제를 포함한 일반차수 FDEs까지 얼마나 넓게 확장할 수 있는가?

주요 결과

• 제안된 일반화된 자코비 함수(GJFs)는 분수계 미분학과 본질적으로 연결되어 있으며, 매개변수 조절을 통해 해의 특이성을 정확히 반영할 수 있다.
• GJF-페트로프-갈레르킨 방법은 저밀도이고 잘 조절된 선형 시스템을 유도하여 기존 스펙트럴 또는 유한차분 방법 대비 계산 비용을 크게 감소시킨다.
• 오차 추정은 해의 정규성에 영향을 받지 않고 데이터의 부드러움에만 의존함을 보여주며, 데이터가 해석적일 경우 스펙트럴 정확도를 달성한다.
• 데이터가 부드러운 경우, 해가 특이하더라도 지수 수렴을 달성하여 방법의 강건성을 입증한다.
• 이 프레임워크는 일반차수의 분수계 초기 및 경계값 문제, 특히 이전에 수치적으로 다루지 못했던 고차수 문제까지도 성공적으로 처리한다.
• 가중 치비셰프 공간에서의 근사 성질을 통해 이론적 결과가 검증되었으며, 엄밀한 수렴 분석이 보장된다.