基于随机Petri网的随机生物过程的模拟算法研究

来源 :哈尔滨工业大学 | 被引量 : 0次 | 上传用户:yanglsm
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
系统生物学的发展使得生物建模和计算机模拟成为生物过程模拟的重要手段,可以有力地支持生物学家及其他研究人员在医学和制药等领域的研究。为建模和模拟像代谢路径、遗传调控网络和细胞信号传导路径这样的复杂生物系统,在过去几年已经发展了多种成功的方法。而且,对于建模活体内的反应,随机建模已经成为一种物理上更实际的选择。  为了正确而有效地建模和模拟生物学系统,我们构建了一个基于随机Petri网的生物模拟软件系统CELLMATLAB,在这个系统中实现了多种随机算法,用于模拟随机Petri网表示的随机生物过程。  本文主要实现了三类不同的模拟算法,即精确的随机模拟算法(SSA)、随机跳步算法和混合模拟算法,并指出其优点和存在的缺陷。SSA具有很高的计算精度,但是计算代价很大;随机跳步算法的模拟效率很高,但是存在一定的精度损失,而且应用范围有限;混合模拟算法介于前两者之间,适于刚性系统的模拟。  这些算法彼此之间虽然存在一定的差异,但都是以分子物种的个数为变量。对于大型的复杂系统,尤其是刚性系统,反应速度和分子数量均在几个数量级上变化,那么这种表示方式显得不够实际。于是,本文提出了一种改进的动态混合模拟算法,根据反应的速度和分子物种的浓度将系统分为连续反应体系和离散反应体系,分别使用一组描述分子物种浓度的瞬时变化的微分方程和精确的随机模拟算法来求解。对于数量很高的分子物种,这种浓度的表示方法更为实际,而且可以近似为连续的过程;对于数量很低的分子,仍旧使用分子数量来表示,看作是离散的个体。由于确定性反应常数和随机性反应常数之间存在一定的转换关系,那么根据反应物的状态来调整其使用的模拟方法是完全可行的。这种动态的混合模拟,既能反映大数量物种的连续性,又能及时捕获少量分子的波动性,在性能上获得了很大的提高。  最后,本文通过简单的四个反应模型、细胞内病毒动力学的模型和E.Coli的热休克模型这三个示例,不仅验证了各种算法的效率和精确度,而且证明本文提出的改进算法是正确的和可用的,尤其适于刚性系统的模拟。
其他文献
近年来,随着无线应用需求和市场的不断增长,特别是随着3G技术的到来,无线互联网将得到前所未有的发展。伴随着无线互联网的发展,人们的学习方式也在不断的改变,从传统的课堂
命名实体识别是目前自然语言处理研究的热点问题。MUC(Message Understanding Conferences)对命名实体的定义是:人们感兴趣的专有名词和特定的数量词,它一般可分为:人名、地
云计算的发展受到全球各国政府和企业高度重视,成为全球信息通信技术产业界公认的发展重点。同时,云计算带来了第三次信息改革浪潮,彻底改变了人们使用IT服务的方式,降低了社
数据复制是分布式数据库中提高系统可用性和可靠性的一项重要技术。长期以来,数据对象均使用静态复制方式,数据复制计划由分布式数据库管理者指定,该计划是固定的,直到管理者
Ontology是近年来知识工程领域的一个研究热点,鉴于其在知识表示和异构数据处理方面的独到之处,而在企业信息化中得到了广泛的应用。随着信息技术的飞速发展,要求构建的Ontol
对于企业来讲,客户关系己经成为企业生存发展所必须关注的问题,企业以客户为中心,不仅是要单纯的发挥主观能动性为客户提供服务,更要在客户的挖掘上下功夫,只有努力赢取客户,维持客
随着云计算的快速发展,云服务提供商将会为越来越多的用户提供服务,由于单个云服务提供商的能力是有限的,它们为用户提供按需的服务时,可能无法满足用户的服务要求,或者所提
RNA(ribonucleic acid)二级结构预测是分子生物学中的关键问题。利用X射线晶体衍射和核磁共振(Nuclear Magnetic Resonance)等实验方法测定RNA分子的立体结构,可以得到比较精确、
同步距离是对两组事件间依赖程度的定量描述,也是刻画系统动态行为的工具。由于同步距离的计算涉及到网系统的结构和网系统的初始标识,这无疑给同步距离的求解带来了很大的难度
嵌入式系统是由硬件和软件组成的,用来完成特定任务的专用系统。嵌入式系统复杂程度不断提高,特别是集成电路设计进入SoC开发时代之后,传统开发方法已难以适应系统设计的需要。