A stochastic exponential Euler scheme for simulation of stiff biochemical reaction systems
A stochastic exponential Euler scheme for simulation of stiff biochemical reaction systems
复制标题
DOI：
10.1007/s10543-014-0485-1
复制
发表时间：
2014-04
影响因子：
1.5
通讯作者：
Y. Komori;K. Burrage
中科院分区：
文献类型：
--
作者：
Y. Komori;K. Burrage
为了模拟刚性生化反应系统，针对具有半线性漂移项的多维非交换随机微分方程导出了显式指数欧拉格式。该方案具有二分之一强阶且均方 A 稳定。将此方法与投影方法相结合，在处理人类以太 a-go-go 相关基因离子通道模型的化学 Langevin 方程的替代公式的数值实验中显示出良好的性能。
In order to simulate stiff biochemical reaction systems, an explicit exponential Euler scheme is derived for multi-dimensional, non-commutative stochastic differential equations with a semilinear drift term. The scheme is of strong order one half and A-stable in mean square. The combination with this and the projection method shows good performance in numerical experiments dealing with an alternative formulation of the chemical Langevin equation for a human ether a-go-go related gene ion channel model.