[论文解读] An efficient method to solve the mathematical model of HIV infection for CD8+ T-cells

Samad Noeiaghdam, Emran Khoshrouye Ghiasi|arXiv (Cornell University)|Jun 28, 2019

Fractional Differential Equations Solutions参考文献 52被引用 7

一句话总结

本文提出同伦分析变换法（HATM）求解涉及CD8+ T细胞的HIV感染非线性数学模型。通过结合同伦分析法与拉普拉斯变换，作者通过收敛性分析、$\bar{h}$-曲线及$N=5$和$N=10$的残差误差图，展示了高精度解，验证了该方法在复杂生物系统中的准确性和鲁棒性。

ABSTRACT

In this paper, the mathematical model of HIV infection for CD8+ T-cells is illustrated. The homotopy analysis method and the Laplace transformations are combined for solving this model. Also, the convergence theorem is proved to demonstrate the abilities of presented method for solving non-linear mathematical models. The numerical results for N = 5, 10 are presented. Several h-curves are plotted to show the convergence regions of solutions. The plots of residual error functions indicate the precision of presented method.

研究动机与目标

• 开发一种高效半解析方法，求解涉及CD8+ T细胞的HIV感染动力学非线性常微分方程组模型。
• 将同伦分析法（HAM）与拉普拉斯变换结合，构建新颖的HATM框架，以提升收敛性与精度。
• 通过严格的收敛性定理，证明给定HIV-CD8+ T细胞模型的HATM解的收敛性。
• 通过不同近似阶数（$N=5$，$N=10$）的残差误差函数与$\bar{h}$-曲线分析，验证方法的精度。

提出的方法

• 首先对HIV-CD8+ T细胞模型（1）中的每个方程应用拉普拉斯变换，将系统转化为代数方程组。
• 利用辅助参数$\bar{h}$构造同伦，解表示为嵌入参数$q$的幂级数形式，并推导出$n$阶变形方程。
• 对$n$阶变形方程应用反拉普拉斯变换，得到$T_n(t), I_n(t), V_n(t), Z_n(t), Z_{a,n}(t)$的递推表达式。
• 近似解通过截断级数获得：$T_N(t) = \sum_{n=0}^{N} T_n(t)$，其他变量形式类似。
• 通过绘制$\bar{h}$-曲线优化收敛控制参数$\bar{h}$，以识别级数收敛的区域。
• 计算并绘制残差误差函数$E_{N,T}(t), E_{N,I}(t), \dots$，以量化解的精度。

实验结果

研究问题

• RQ1 HATM能否有效求解涉及CD8+ T细胞的HIV感染动力学非线性常微分方程组模型？
• RQ2 HATM解在HIV-CD8+ T细胞模型中的收敛区域是什么？其随近似阶数$N$如何变化？
• RQ3 残差误差函数衡量下，HATM解与精确解相比的精度如何？
• RQ4 收敛控制参数$\bar{h}$的不同取值如何影响HATM解的收敛性与稳定性？
• RQ5 HATM是否为求解HIV-CD8+ T细胞系统等复杂非线性生物数学模型提供了可靠且灵活的框架？

主要发现

• 当$N=5$时，HATM解被显式推导为$t$、$\bar{h}$及其幂的级数形式，系数包含模型参数。
• 当$N=5$时，解的收敛区域为$-0.9 \leq \bar{h}_T \leq -0.2$，$-0.8 \leq \bar{h}_I \leq -0.2$，$-0.8 \leq \bar{h}_V \leq -0.2$，$-1.2 \leq \bar{h}_Z \leq -0.6$，以及$-0.8 \leq \bar{h}_{Z_a} \leq -0.4$。
• 当$N=10$时，收敛区域扩展为$-0.9 \leq \bar{h}_T \leq -0.4$，$-1 \leq \bar{h}_I \leq -0.2$，$-1 \leq \bar{h}_V \leq -0.3$，$-1.2 \leq \bar{h}_Z \leq -0.4$，以及$-1 \leq \bar{h}_{Z_a} \leq -0.4$，表明高阶近似提升了稳定性。
• 在$\bar{h} = -0.8$时，$N=5$和$N=10$的残差误差函数在时间域内均接近零，证实了解的高精度。
• 收敛性定理得到严格证明，确立HATM作为求解HIV-CD8+ T细胞模型等非线性系统可靠方法的地位。
• 该方法展现出强大的灵活性与精度，当$N=10$时，残差误差始终低于$10^{-3}$，表明具有高度计算可靠性。