[論文レビュー] 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

ひとこと要約

本稿では、CD8+ T細胞を含む非線形なHIV感染の数学的モデルを解くために、ホモトピー解析変換法（HATM）を提案する。ホモトピー解析法とラプラス変換を組み合わせることで、$\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}$-曲線解析を通じて、手法の精度を検証すること。

提案手法

• HATMは、まず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 HIV-CD8+ T細胞モデルにおけるHATM解の収束領域は何か？また、近似次数$ 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} $未満に保たれ、高い計算信頼性を示した。