[论文解读] An approximation formula for the Katugampola integral

Ricardo Almeida, Nuno R. O. Bastos|arXiv (Cornell University)|Nov 24, 2015

Fractional Differential Equations Solutions参考文献 10被引用 8

一句话总结

本文提出了一种新颖的近似公式，用于Katugampola分数阶积分，仅使用一阶导数，将分数阶问题转化为标准常微分方程。该方法在近似参数N增大时表现出高精度，并收敛于精确解，通过数值例子和对分数阶积分方程的应用得到验证。

ABSTRACT

The objective of this paper is to present an approximation formula for the Katugampola fractional integral, that allows us to solve fractional problems with dependence on this type of fractional operator. The formula only depends on first-order derivatives, and thus we convert the fractional problem into a standard one. With some examples we show the accuracy of the method, and then we present the utility of the method by solving a fractional integral equation.

研究动机与目标

• 开发一种用于求解涉及Katugampola分数阶积分算子的分数阶微分与积分方程的实用数值方法。
• 通过将直接分数阶积分简化为常微分方程系统，克服直接分数阶积分的计算复杂性。
• 将现有针对Riemann–Liouville和Hadamard分数阶积分的近似技术推广至更一般的Katugampola框架。
• 提供一种计算高效且收敛的数值方法，适用于物理与工程领域中的实际分数阶模型。
• 通过数值实验和对具有已知解析解的测试分数阶积分方程的应用，验证该方法。

提出的方法

• 通过依赖于分数阶阶数α和参数ρ的系数的级数展开，推导出左Katugampola分数阶积分的近似公式。
• 引入辅助函数W_k(t)，定义为x(τ)加权于(t^ρ - τ^ρ)^{k-1}的积分，通过初值问题求解。
• 将分数阶积分表示为包含A、B_k系数和W_k(t)的有限和，余项E_N(t)在N → ∞时趋于零。
• 将相同方法应用于右Katugampola分数阶积分，得到对称的近似公式。
• 通过将近似公式代入原分数阶积分方程，将原始方程转化为常微分方程组。
• 使用数值积分和ODE求解器（例如MATLAB）计算给定N和参数下变换后系统的解。

实验结果

研究问题

• RQ1 能否构建一种基于一阶导数的近似方法，用于Katugampola分数阶积分，从而避免直接进行分数阶积分？
• RQ2 近似精度如何依赖于截断参数N？当N增大时，误差是否收敛至零？
• RQ3 所提出的近似方法能否有效将分数阶积分方程转化为可求解的常微分方程组？
• RQ4 该方法在不同分数阶阶数α和形状参数ρ下表现如何？
• RQ5 与具有已知解析解的测试函数相比，该近似的数值行为如何？

主要发现

• 近似误差E_N(t)在N → ∞时收敛至零，其上界与1/N^α成正比，确保理论收敛性。
• 对于测试函数x(t) = t^{2ρ}，数值近似与精确分数阶积分高度吻合，且随着N增大，误差持续减小。
• 该方法成功求解了具有已知解x(t) = t^{2ρ}的分数阶积分方程，在中等N值下以高精度恢复精确解。
• 图2中的数值结果表明，该方法在不同α和ρ值下均保持高精度，包括非整数和较大的分数阶阶数。
• 收敛速率被定量界定为O(1/N^α)，且余项E_N(t)通过Gamma和Beta函数显式估计。
• 该方法可将复杂的分数阶问题转化为标准ODE系统，显著简化数值求解过程。