[论文解读] Fractional Order Optimal Control Model For Malaria Infection

Eric Okyere, Francis T. Oduro+2 more|arXiv (Cornell University)|Jul 6, 2016

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

一句话总结

本文提出了一种使用Caputo导数的分数阶最优控制模型，以在传播动力学中引入记忆效应。通过应用广义欧拉格式的前向-后向扫掠法，研究结果表明，与整数阶模型相比，具有时变控制（蚊帐、治疗、杀虫剂喷洒）的分数阶模型能显著减少感染的人类和蚊子数量。

ABSTRACT

We propose and study an optimal control model for malaria infection described by system of fractional differential equations. The model is formulated in terms of the left and right Caputo fractional derivatives. We determine the necessary conditions for the optimality of the controlled dynamical system.The forward-backward sweep method with generalized euler scheme is applied to numerically compute the solutions of the optimality system.

研究动机与目标

• 开发一种用于疟疾传播的分数阶最优控制模型，通过Caputo导数捕捉记忆效应。
• 利用变分法和分数阶分部积分公式，推导受控系统的必要最优性条件。
• 使用前向-后向扫掠法结合广义欧拉格式，数值求解最优性系统。
• 比较不同控制策略（经处理蚊帐、治疗、杀虫剂喷洒）在整数阶与分数阶模型下的有效性。
• 证明具有时变控制的分数阶模型相比整数阶模型能实现更优的疾病控制效果。

提出的方法

• 采用左、右Caputo分数阶导数构建疟疾传播模型，以表征人类和蚊子种群的记忆依赖动力学。
• 利用变分法和分数阶分部积分公式，推导分数阶最优控制问题的必要最优性条件。
• 应用前向-后向扫掠法求解涉及状态方程与伴随方程的两点边值问题。
• 采用广义欧拉方法作为数值格式，近似最优性系统中分数阶微分方程的解。
• 实施多种控制策略：经处理蚊帐（u₁）、治疗（u₂）和杀虫剂喷洒（u₃），每种均为时变函数。
• 采用现有文献中的完整初始条件与参数值集合，模拟并比较不同控制情景及分数阶参数（α = 1, 0.99, 0.95, 0.90）下的结果。

实验结果

研究问题

• RQ1 与整数阶模型相比，分数阶导数在疟疾传播动力学建模中有哪些改进？
• RQ2 在使用Caputo导数的疟疾传播模型中，分数阶最优控制问题的必要最优性条件是什么？
• RQ3 在减少感染人类和蚊子数量方面，单个与组合使用的时变控制（蚊帐、治疗、杀虫剂喷洒）的有效性如何？
• RQ4 整数阶模型（α=1）与分数阶模型（α<1）在疾病控制结果上的相对影响是什么？
• RQ5 在不同控制策略与分数阶参数下，最优控制函数（u₁, u₂, u₃）如何随时间演变？

主要发现

• 具有时变控制的分数阶模型（α < 1）相比整数阶模型（α = 1），能显著减少感染人类与蚊子的数量。
• 同时使用三种控制手段——经处理蚊帐、治疗与杀虫剂喷洒——可使感染人类与蚊子数量的减少幅度最大。
• 当仅使用单一控制时，经处理蚊帐（u₁）对减少感染人类的效果最显著，其次为治疗（u₂）与杀虫剂喷洒（u₃）。
• 最优控制函数（u₁, u₂, u₃）表现出随时间变化的行为，初始值较高并逐渐下降，表明早期高强度干预为最优策略。
• 数值模拟证实，广义欧拉方法能以高精度有效逼近分数阶最优性系统的解。
• 分数阶导数固有的记忆效应增强了模型对现实世界疾病动力学的反映能力，尤其在长期控制策略中表现更优。