[论文解读] Analysis of some generalized ABC- fractional logistic models

Thabet Abdeljawad, Mohamed Ali Hajji+2 more|arXiv (Cornell University)|Dec 15, 2019

Fractional Differential Equations Solutions参考文献 20被引用 5

一句话总结

本文引入并分析了使用多参数Mittag-Leffler核的广义ABC-分数阶逻辑斯蒂模型。建立了存在性与唯一性定理，通过平衡点扰动法进行稳定性分析，并通过数值模拟表明：当γ或α增大时，收敛至平衡态的速度加快；当μ增大时，收敛速度减慢；随着γ或α增加，早期阶段会出现振荡行为。

ABSTRACT

In this article, some logistic models in the settings of Caputo fractional operators with multi-parametered Mittag-Leffer kernels (ABC) are studied. This study mainly focuses on modified quadratic and cubic logistic models in the presence of a Caputo type fractional derivative. Existence and uniqueness theorems are proved and stability analysis is discussed by perturbing the equilibrium points. Numerical illustrative examples are discussed for the studied models.

研究动机与目标

• 将经典逻辑斯蒂模型扩展至广义ABC-分数阶导数框架，采用多参数Mittag-Leffler核。
• 研究在新型分数阶算子下，修改后的二次与三次逻辑斯蒂方程解的存在性与唯一性。
• 通过扰动技术分析平衡点的稳定性。
• 数值展示分数阶参数（α, μ, γ）对解轨迹与收敛动力学的影响。

提出的方法

• 采用基于积分表达式 $ E^{ ho}_{ heta,eta}( heta, t) $ 的广义左侧行ABC-分数阶导数，其核为多参数Mittag-Leffler函数。
• 应用广义AB分数阶积分算子以逆向ABC导数，实现通过积分方程表示解。
• 采用基于离散卷积与递归系数计算的数值格式，其中 $ ilde{c}^{m}_{j} = sum_{k=0}^{ inf} (h^{ heta k - u + 1} A_{k}^{ heta, u, ho}) ilde{r}^{m}_{k,j} $。
• 推导隐式数值格式：$ x_m = x_0 + sum_{j=0}^{m} ilde{c}^{m}_{j} f_j $，其中 $ f_j $ 为在第j步计算的非线性函数。
• 利用广义ABC导数的性质：$ {}^{AB}I^{ heta, u, ho} {}^{ABC}D^{ heta, u, ho} u(t) = u(t) - u(t_0) $，推导解的形式。
• 通过固定参数（r=0.5, K=2）进行数值模拟，改变α、μ、γ及初始条件 $ x_0 $，研究收敛性与振荡行为。

实验结果

研究问题

• RQ1 广义ABC-分数阶导数中的参数α、μ与γ如何影响逻辑斯蒂模型解的存在性与唯一性？
• RQ2 广义ABC-分数阶逻辑斯蒂模型中平衡点（x=0与x=2）的稳定性行为如何？
• RQ3 初始条件 $ x_0 $ 如何影响收敛至稳定平衡态的速度？
• RQ4 参数γ与α在解轨迹的早期演化中如何引发振荡行为？
• RQ5 参数μ如何影响解达到平衡态所需的时间？

主要发现

• 对所有初始值 $ x_0 $，解轨迹均渐近收敛至平衡点 $ x_2 = 2 $，证实 $ x_2 = 2 $ 为渐近稳定，而 $ x_1 = 0 $ 为不稳定。
• 当参数 $ \nu $ 较小时，收敛至平衡态的速度更快，表明较低的 $ \nu $ 值可加速趋近稳态。
• 当参数 $ \rho $ 从1.0增至3.0时，达到平衡态的时间缩短，并在解轨迹的早期阶段引入振荡行为。
• 当分数阶阶数 $ \alpha $ 从0.7增至0.99时，解达到平衡态所需时间减少，表明更高的 $ \alpha $ 值可实现更快收敛。
• 当初始条件 $ x_0 $ 更接近稳定平衡态 $ x_2 = 2 $ 时，收敛速率更高，表明收敛对初始值具有强依赖性。
• 当 $ \rho $ 或 $ \alpha $ 增大时，早期阶段观察到振荡行为，表明这些参数的较高值会在收敛前引入瞬态振荡。