数学理论与应用 ›› 2023, Vol. 43 ›› Issue (4): 76-92.doi: 10.3969/j.issn.1006-8074.2023.04.005

• • 上一篇    下一篇

系统生物学中的随机微分方程数值仿真算法

牛原玲1,*,陈琳2,陈洛南3   

  1. 1.中南大学数学与统计学院,长沙, 410083; 2. 江西财经大学统计与数据科学学院,南昌, 330013; 3. 中国科学院系统生物学重点实验室,分子细胞科学卓越创新中心(上海生物化学与细胞生物学研究所),上海, 200031
  • 出版日期:2023-12-28 发布日期:2024-01-03
  • 通讯作者: 牛原玲; E-mail: yuanlingniu@csu.edu.cn
  • 基金资助:
    国家自然科学基金项目(Nos. 12071488, 11971488, 12371417, 11961029)资助

Numerical Simulation Algorithms for Stochastic Differential Equations in Systems Biology

Niu Yuanling1,*, Chen Lin2, Chen luonan3   

  1. 1. School of mathematics and Statistics, Central South University, Changsha 410083, China; 2. School of Statistics and Data Science, Jiangxi University of Finance and Economics, Nanchang 330013, China; 3. Key Laboratory of Systems Biology, Center for Excellence in Molecular Cell Science (Shanghai Institute of Biochemistry and Cell Biology), Chinese Academy of Sciences, Shanghai 200031, China
  • Online:2023-12-28 Published:2024-01-03

摘要: 系统生物学中的诸多现象, 如生物化学反应过程、生态系统的演变、传染病的传播等, 都可以用随机微分方程来描述. 由于考虑了随机因素的影响, 随机微分方程模型往往能比确定性的微分方程模型更为准确地刻画变量随时间的演化规律. 但是随机微分方程的真解大多不可得到, 有的即使可以求出真解, 但解的形式极其复杂, 用起来十分不便. 因此, 在计算机上对其进行数值仿真就显得十分必要. 系统生物学中的随机微分方程模型一般呈现出高维、高度非线性、真解位于某些特定的区域等特点, 对它们的数值模拟需要做专门的研究. 本文概述求解几类常见的系统生物学模型(生物化学反应模型、生态系统模型、传染病模型、群体遗传学模型、细胞分化模型)的数值算法及这些数值算法各自的优缺点.

关键词: 系统生物学, 随机微分方程, 数值模拟, 算法

Abstract: Many phenomena in systems biology, such as the biochemical reaction process, the evolution of ecosystems, the spread of infectious diseases, can be described by stochastic differential equations (SDEs). Considering the influence of randomness, stochastic differential equation models can describe the evolution of variables over time more accurately than deterministic differential equation models. However, the analytical solutions of most stochastic differential equations cannot be obtained. Even though some of them can be obtained, the forms of the solutions are usually extremely complex. One therefore requires proper numerical methods to approximate their solutions on computers. These stochastic differential equation models in systems biology usually have the properties of high dimension, high nonlinearity, and the solutions being located in a specified region. It is difficult to simulate them numerically. This paper reviews the numerical simulation algorithms of several typical models in systems biology (biochemical reaction models, ecosystem models, infectious disease models, population genetics models, cell differentiation models), and briefly introduces their advantages and disadvantages.

Key words: Systems biology, Stochastic differential equation, Numerical simulation, Algorithm