[논문 리뷰] Solving generalized Abel's integral equations of the first and second kinds via Taylor-collocation method

Eisa Zarei, Samad Noeiaghdam|arXiv (Cornell University)|2018. 04. 23.

Fractional Differential Equations Solutions참고 문헌 28인용 수 11

한 줄 요약

이 논문은 일반화된 아벨 적분방정식의 제1종 및 제2종을 해결하기 위해 타일러 다항식과 콜로케이션을 조합한 새로운 수치적 방법을 제안한다. 미지 함수를 타일러 전개로 근사하고, 이산적인 점들에서 콜로케이션 조건을 적용함으로써, 이 적분방정식을 선형 시스템으로 변환한다. 수치적 결과들은 높은 정확도와 빠른 수렴 속도를 보이며, 테스트된 예제들에서 n=9일 때 최대 오차가 1.4×10⁻¹³ 이하로 감소함을 보여준다.

ABSTRACT

The aim of this paper is to present an efficient numerical procedure to approximate the generalized Abel's integral equations of the first and second kinds. For this reason, the Taylor polynomials and the collocation method are applied. Also, the error analysis of presented method is illustrated. Several examples are approximated and the numerical results show the accuracy and efficiency of this method.

연구 동기 및 목표

• 일반화된 아벨 적분방정식의 제1종 및 제2종을 효율적으로 해결하기 위한 수치적 체계를 개발하기 위해.
• 분수계수 특이핵을 가진 볼테라 유형의 적분방정식에서 발생하는 특이핵 문제를 다루기 위해.
• 하이브리드 타일러-콜로케이션 접근법을 통해 기존 방법에 비해 정확도와 수렴 속도를 향상시키기 위해.
• 제안된 방법에 대한 엄밀한 오차 한계와 수렴 분석을 제공하기 위해.
• 정확한 해가 알려진 기준 문제들에 대해 방법을 검증하기 위해.

제안 방법

• 지정된 점 z를 중심으로 차수 n인 타일러 다항식을 사용하여 미지 함수 Φ(x)를 근사하기 위해.
• n+1개의 등간격 점 xi = a + (b−a)i/n (i=0에서 n까지)을 선택하여 콜로케이션 방법을 적용하기 위해.
• 타일러 전개를 적분방정식에 대입하고, 각 콜로케이션 점에서 방정식을 강제로 만족시키기 위해.
• 선형 시스템 AX = G를 구성하기 위해, 여기서 A는 (t−z)^j / (ϕ(xi)−ϕ(t))^α의 적분을 포함하고, X는 타일러 계수 Φ^(j)(z)를 담고, G는 g(xi)를 포함한다.
• 선형 시스템을 풀어 타일러 계수를 결정하고, 근사해를 재구성하기 위해.
• 최종 다항식을 사용하여 다양한 x 및 n 값에서 수치적 해와 절대 오차를 계산하기 위해.

실험 결과

연구 질문

• RQ1 타일러-콜로케이션 방법은 일반화된 아벨 적분방정식의 제1종 및 제2종을 해석할 때 높은 정확도를 달성할 수 있는가?
• RQ2 제안된 방법의 수렴 속도는 타일러 다항식의 차수 n에 따라 어떻게 달라지는가?
• RQ3 다양한 n 값과 테스트 문제들에 대해 최대 절대 오차의 크기는 어느 정도인가?
• RQ4 특히 수치적 일致성 측면에서 정확한 해가 알려진 문제들에 대해 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 0 < α < 1 인 특이핵을 효과적으로 다룰 수 있으며, 안정성과 정확도를 유지할 수 있는가?

주요 결과

• 타일러-콜로케이션 방법은 높은 정확도를 달성하였으며, n=9일 때 예제 1에서 최대 절대 오차가 1.37683×10⁻¹³ 으로 감소하였다.
• 예제 2(Φ(x) = cos x)의 경우, n=9일 때 최대 오차가 3.06992×10⁻¹⁰ 으로 감소하여 빠른 수렴을 보였다.
• 예제 3(Φ(x) = e^x)에서는 n=9일 때 최대 오차가 1.28571×10⁻⁸ 로 나타나, n이 증가함에 따라 일관되게 향상됨을 보였다.
• 예제 4(Φ(x) = x²)에서는 n=5일 때, 예제 5(Φ(x) = 1−2x)에서는 n=3일 때 정확한 해를 도출하여, 특정 다항식 해에 대해 정확성이 입증되었다.
• 오차 분석 결과, 방법의 안정성과 오차 전파의 유한성 확인되었으며, n이 증가함에 따라 오차가 크게 감소함을 확인하였다.
• 모든 예제에서의 수치적 결과는 이 방법의 강건성, 효율성 및 기존 표준 방법에 비해 뛰어난 성능을 보임을 보여주었다.