应用科学学报 ›› 2009, Vol. 27 ›› Issue (3): 266-0271.

• 信号与信息处理 • 上一篇    下一篇

改进的τ-leap算法在生化反应系统随机模拟中的应用

刘焕1, 彭新俊2;3, 周文4, 王翼飞1   

  1. 1. 上海大学数学系,上海200444
    2. 上海师范大学数学系,上海200234
    3. 上海市高校科学计算重点实验室,上海200234
    4. 安徽师范大学数学与计算机学院,安徽芜湖241000
  • 收稿日期:2008-09-19 修回日期:2008-12-28 出版日期:2009-05-28 发布日期:2009-05-28
  • 通信作者: 通信作者: 王翼飞,教授,博导,研究方向:计算分子生物学 E-mail:yifei_wang@staff.shu.edu.cn
  • 基金资助:
    国家自然科学基金(No.30871341);国家“863”高技术研究发展计划基金(No.2006AA02Z190)资助项目

Application of Improved τ-leap Algorithm in Stochastic Simulation of Biochemical Reaction Systems

  1.  1. Department of Mathematics, Shanghai University, Shanghai 200444, China
    2. Department of Mathematics, Shanghai Normal University,Shanghai 200234, China
    3. Scientific Computing Key Laboratory of Shanghai Universities, Shanghai 200234, China
    4. Mathematics and Computation College, Anhui Normal University, Wuhu 241000, Anhui Province, China
  • Received:2008-09-19 Revised:2008-12-28 Online:2009-05-28 Published:2009-05-28

摘要:

     提出了一种t -选择策略,有效地反映了生化反应系统中分子数目的改变. 并由此提出了改进的t-leap(improved (τ-leap)算法,该算法对生化反应系统的随机模拟更为有效和实用. 并以两个生化反应系统模型为例,分别用精确的SSA算法、改进的τ-leaping算法以及已有的修正的τ-leap(modified tau-leap)算法进行了模拟计算. 仿真实验结果表明:在具有同等计算复杂度的情况下,改进的τ-leap算法较修正的τ-leap明显地提高了模拟精度.

关键词:  随机模拟算法 , τ-选择策略 , τ-leap算法 , 生化反应系统

Abstract:



    A τ-selected strategy that ascertain the τ-value is proposed in this paper. It is developed as an improved τ-leap algorithm. The τ-selected strategy can effectively reflect change numbers of the species in biochemical systems, while the improved τ-leaping algorithm is more effective and practical to stochastic simulation of biochemical reaction systems than other algorithms. We take two models of biochemical reaction systems as examples, and carry out simulation using accurate SSA algorithm, improved τ--leaping algorithm, and the present modified τ-leap respectively. Numerical results demonstrate that, compared with modified τ-leap algorithm, the improvedτ-leap algorithm can significantly increase simulation precision with the same computation complexity under the same condition.

Key words: stochastic simulation algorithm , τ-selected strategy , τ-leap algorithm , biochemical reaction system

中图分类号: