生化反应系统中加速随机模拟算法的研究

来源 :安徽师范大学 | 被引量 : 0次 | 上传用户:shengli1011
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
一切生命活动的基础是细胞,细胞的活动构成一套精密而复杂的化学反应网络。因此研究这些相互作用的化学反应网络已成为系统生物学研究中的一个重要课题。为了研究复杂化学反应的动力学特征,需要对生反应系统进行建模和仿真。由于细胞内的化学反应都具有随机性和离散性,这就使得我们在建模过程中注重“微观离散随机”的模型。  具有代表性反应的算法(Representative Reaction Approach,RRA)的提出为研究生化反应系统开拓出了一个新的研究方向,它主要的思想是以反应2A?B代表整个生化反应系统中所有的反应。本文在RRA算法基础之上,提出“最后所有可能步进的RRA算法(Final All Possible Step RRA,FAPS RRA)”。该算法改变了RRA算法的数据收集方式,目的是在不失RRA算法精度的前提下,减少模拟的运行次数,结果表明在相同精度的前提下,FAPS RRA算法比RRA算法的模拟运行速率有显著的提高。  时滞是生化反应系统中普遍存在的现象,如基因转录和mRNA翻译等都不是立即发生而是经过一段时间后才发生反应。Barrio等人与 Cai提出的精确时滞随机模拟算法(DSSA)均用来解决含时滞的生化反应系统中的动力学问题。由于DSSA每次只能模拟一个反应事件,模拟生化反应进程时速率过低。本文在现有RRA算法的基础上进行拓展,将其应用到含时滞的生化反应系统中,提出了时滞生化反应系统中RRA算法(DRRA)。仿真实验表明:与DSSA算法相比较,我们提出的DRRA算法在模拟含时滞的生化反应系统时更为有效与实用。  考虑到生化反应系统中时滞反应与非时滞反应可能同时发生反应,Bayati提出了含时滞的加速随机模拟算法(D-leaping)。该算法是通过对时滞生化反应系统建立一个以时滞的等待序列为主体的四元结构体,来有效地模拟时滞生化反应系统。在实际运行中,为了得到可靠的统计特性,往往需要对大量的模拟做平均。文中我们提出模拟时滞生化反应系统的加速随机模拟算法,称为“最后所有可能步进的加速算法(DFAPS-leaping)”。该算法没有改变D-leaping算法的模拟过程,却减少了需要运行的次数,比D-leaping算法更有效。数值结果表明,提出的算法能广泛应用到各种时滞生化反应系统中,并与已有算法相比其效率有显著提高。
其他文献
复信道信号恢复技术采用硬件设计电路网络达到使失真信号被恢复的目的,但是用硬件的方法进行信号恢复存在着元件参数需反复调试、对不同类型的通信体制和信道特性的适应能力
本文主要利用匹配渐近展开法、合成展开法、界定函数法等摄动方法和微分不等式理论研究奇摄动内层问题.  第一章引言部分综述了摄动理论与方法的历史发展及有关应用背景,并
全文共分五章:第一章,综述光正交码的研究背景和当前领域的研究状况,并且给出了一些基本的名词和事实.第二章,主要运用重要的组合结构Skew starter,证明了存在一个最优(v,5,1
由于图的带宽和问题在VLSI布局设计和纠错代码设计中的应用背景,引起了一些数学工作者的兴趣.该文讨论了图的圈上带宽和的一些问题.全文分三章.第1章给出了有关圈上带宽与圈
经典的概率极限理论研究的对象主要是随机变量的部分和的弱收敛性或强收敛性,[18][36]就是这方面的经典文献.随着实践的深入,人们越来越发现,许许多多的实际问题都可以转化为
该文运用变分方法和拓扑度方法讨论两类二阶Hamilton系统同宿轨道的存在性和两类二阶常微分方程正同宿轨道的存在性.它由两章组成.第一章借助于Hamilton系统的拉格朗日泛函在
该文研究将平面点集划分为不变凸多边形的有关问题.文献[1]中Kiyoshi Hosono与Masatsugu Urabe讨论了一个平面点集的所有不交分划中不交凸四边形的最大个数.设P为平面中n个点
该文分成两个部分:第一部分:保险公司如何根据多个不同的目标来选择个人代理人的薪酬方案?最常见的做法是由管理层参照他人的做法与已有的经验来进行方案选择.该文在分别对保
由于矩阵广义逆在许多领域中有着广泛的应用,如在微分和积分方程、算子理论、统计学、控制论、Markov链、最优化等.因此,自上个世纪中期以来矩阵广义逆就成为一个非常重要的
二十世纪三十年代P.Jordan最早给出了约当代数的概念.约当代数是一类十分重要的非结合代数,在很多领域都有广泛的应用.本文主要研究三维约当代数上的Rota-Baxter算子和Hom-预