生物化学反应动力学的随机模拟

来源 :南京农业大学 | 被引量 : 0次 | 上传用户:ntudqliweiwei
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随机问题广泛出现于物理学、化学、生物学等自然科学与工程的许多领域([1,2,13,54]).随机微分方程(stochastic differential equation,SDE)是描写随机现象的一种常用工具([23,27,29,33]).然而大多数随机微分方程由于其复杂性(如非线性性)无法求得其精确解,因此数值解法就成为了一个不可替代的途径.由Gillespie创建的随机模拟算法(stochastic simulation algorithm,SSA)([4,9,12,14,18])是数值模拟化学反应过程的一种精确而有效的方法.但是对复杂化学反应系统(反应物较多或反应速度差异较大),SSA的模拟效率较低,需要改进.本文研究了随机模拟算法的优缺点及一些衍生算法,随机模拟算法与传统的隐式方法的复合算法,以及随机Nystrom方法.本文共分四章:第一章给出了本文的研究背景,介绍了 SSA和随机微分方程的发展及基础知识.第二章主要研究了 SSA及其衍生的r-leap算法、化学Langevin方程(chemical Langevin equation,CLE)、慢尺度随机模拟算法(Slow-scale SSA,ssSSA).本章详细给出了生物化学反应的到下一个反应的时间间隔和下一个反应序号的推导过程以及一般情况下倾向性函数的描述和慢尺度倾向性函数的计算方法,然后进行一些数值实验加以佐证;除此之外,本章还研究了 τ-leap算法和化学Langevin方程导出的τ-leap EM方法,并且受Euler-Heun方法([4,5])的启发提出了 τ-leap Heun方法,并对之做了分析.由于SSA每一步都是精确模拟,而且生物化学反应与微分方程可以相互转化(并不总是能转化),再联系到数值方法中的隐式方法在使用计算机模拟的时候需要提前知道下一步的值,于是第三章将SSA与隐式方法(隐式Euler法,隐式中点法,隐式梯形法)结合,记作隐式SSA,通过大量的数值实验证明了对于化学反应率方程来说,隐式SSA方法比传统的隐式方法更适合长时间大步长的计算.最后,从SSA倾向性函数的角度阐述了为什么隐式中点法和梯形法优于隐式Euler法.对于生物化学反应来说,SSA更多倾向于系统内部噪音,受此启发,第四章主要研究其外部噪音,研究方法就是带加性噪音的随机Nystrom方法,记作SN.本章先给出了 SN的定义,然后将之与随机Runge-Kutta方法(SRK)从阶矩的角度作了比较,说明了 SN比SRK精度更高;然后再讨论了其强弱阶、稳定性和收敛性;最后进行了数值试验,从简单的分子分解模型到基因调控网络模型均说明了随机Nystrom方法能更好地表现细胞行为.最后一章本文总结了做出的主要贡献以及一些尚需解决的问题,然后再展望一下未来可以着手研究的方向.
其他文献
本文以汉语能愿动词句特点为基础研究,借助HSK动态作文语料库,进一步研究面向俄罗斯汉语学习者的能愿动词句教学问题。在研究过程中,运用了对比分析和偏误分析的方法。本文的研究过程为:首先对汉语能愿动词句性质、类型、句子成分和表达形式进行梳理与归纳,以确定能愿动词句的特点;其次探索汉语能愿动词句在俄语中的表达形式,并从HSK动态作文语料库中选取偏误实例进行分析,以了解学习者对能愿动词句的使用和掌握情况;
自我国开始实行失信被执行人制度以来,在各级机关、法院的共同努力下,被执行人的生存空间已日益狭隘,形成了上了“黑名单”,寸步难行的局面。该制度的实施不仅有效缓解了判决难以执行的现象,同时提高了司法效率,使债权人的合法权益得以更好的保护。其在提高我国的司法地位、健全我国的信用制度等方面也具有重大意义。但该制度在实施过程中,也存在着对失信被执行人权利保护认知不足、某些法律依据缺位、限制对象的筛选程序不够
目的:基于大样本病案进行多角度、多层次的数据挖掘分析,全面总结广东省名中医周大桥教授辨治原发性肝癌的独到临床经验及学术思想,为临床医生提供具有参考价值且可行的中医诊疗思路。方法:筛选2017年1月1日至2019年12月31日深圳市中医院肝病科周大桥教授门诊就诊的原发性肝癌患者病案,建立数据库,利用中医传承辅助平台(V2.5版本)进行数据挖掘,分析周教授辨治原发性肝癌的证候规律、药物关联规则以及核心
本文以壮族民间故事题材的动画设计为研究对象,结合情感化理念设计的相关理论,通过深入挖掘和梳理壮族民间故事,尝试将壮族民间故事当中所体现的民族文化、民族精神、民族信仰等融入到动画创作当中。同时,通过对情感化设计理论的探究,试图将情感化设计理念与壮族民间故事关联起来,创作带来民族风格倾向的动画。探索情感化设计在壮族民间故事题材动画创作应用的规律,研究壮族民间故事在动画创作中的设计要素与情感化处理方法,
网络的发展和运用使得网络的交往变得十分的快捷,网络信息呈现了快速的扩散趋势,在网络信息中各种虚假信息和负面信息,可能会对网络信息使用者带来了一定的影响,容易引发社会问题。高中生处于人生成长的关键阶段,由于年龄和知识的增长,高中生对于外界信息获取的渴望越来越高,但高中生心智发展还不成熟,其社会经验和社会辨别力还较弱。在网络信息的充斥下,各种信息出现了交流,影响到了高中生的正常学习和生活。随着网络的发
目的:本文以制何首乌(Radix Po lygoni multiflora praeparata)为研究对象,在制何首乌化学成分的研究基础上,将传统经典方法与现代分析技术相结合,提升制何首乌质量标准,并进行制何首乌饮片质量等级评价标准研究,为制何首乌饮片的质量等级标准制订提供一种新方法、新思路。方法:1.采用正相硅胶、ODS反相硅胶层析法、LH-20凝胶及半制备HPLC等方法对制何首乌70%乙醇提
填筑用红黏土的改良,主要通过添加石灰、水泥等碱性材料来实现,但由于在现场拌和时,石灰和水泥等碱性稳定材料很难均匀拌和入填筑材料内部,只能附着在土团的表面,导致改良效果不甚理想。同时,室内试验采用的填筑用土,根据各规范要求严格限制了土粒尺寸,致使稳定材料在室内试验研究与实际的工程应用中表现出一定的性能差异。并且在长期服役过程中,红黏土呈弱酸性的特点易与碱性改良材料发生互损效应,使其耐久性劣化,对工程
目的:制定紧张型头痛患者报告结局指标评价量表初稿,对量表初稿的性能进行考核,用于在中西医领域中对紧张型头痛患者病情及治疗效果的评估。方法:根据美国食品药品管理局(FDA)、指南及国际药物经济学和成果研究学会(ISPOR)制定的PRO量表设计方法,采用扎根理论进行研究,收集紧张型头痛患者关于在紧张型头痛发病、治疗过程中所关注的症状、疗效等信息进行三级编码,定义所需要测量的概念,设置测量维度与条目,形
我国目前仍是发展中国家,基础设施还不完善,建筑用房、工程建设都需要大量的砂子,河沙资源日益枯竭。各个地区工程多用河沙,而大量工程对河沙的开采导致环境污染,地质地貌被破坏。由于国家对自然资源的管控,对工程来说砂子的成本也将提高。因此大量学者研究了沙漠砂用于建筑细骨料来代替河沙以减少对河沙的需求和依赖。新疆地处西域,新疆沙漠砂资源丰富,便于开采,沙漠砂混凝土由此产生。然而由于沙漠砂比表面积大,搅拌混凝
目的:本研究旨在建立一种可靠的在动脉瘤性蛛网膜下腔出血(Aneurysmal subarachnoid hemorrhage,a SAH)患者行开颅夹闭或介入栓塞术前预测迟发性脑缺血(delayed cerebral ischemia,DCI)发生的风险预测模型,以期为临床实践提供帮助。方法:我们收集了2013年1月1日至2019年7月30日期间在我院神经外科接受开颅夹闭术或介入栓塞术治疗的a S