马氏链计算软件开发及其在生物信息与复杂网络中的应用

来源 :上海大学 | 被引量 : 0次 | 上传用户:saxthon
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本论文是以作者攻读硕士学位期间所承担研究课题的工作为基础,着重马氏链的相关数值计算软件开发及其在复杂网络和生物信息研究中的应用。 首先,典型的复杂网络是一种具有增长和择优机制的网络。在Barabasi和Albert使用动力学方法(平均场)给出了BA(Barabasi-Albert)模型的度分布的解析解后,人们在处理更为复杂的模型时,并不能给出度分布的解析解。因此,史定华,陈庆华和刘黎明提出了应用非齐次马氏链计算复杂网络度分布的方法,本人的工作则从马氏链的转移概率矩阵出发,根据矩形迭代算法,进行具体的数值计算软件的开发,分别得到基本BA模型,m可变BA模型和两个演化网络模型的网络度分布,并与动力学方法得到的结果相比较,说明利用马氏链计算复杂网络度分布的有效性。 其次,多序列联配是生物信息学中重要的生物序列分析手段,而剖面隐马氏模型是获得多序列联配的一种有效方法,一般需要经过模型的初始化、训练、联配三个过程。然而,目前广泛采用的Baum-Welch训练算法假设各条可观察序列相互独立,这与实际情况有所不符。本文对剖面隐马氏模型,给出可观察序列在不独立情况下的改进Baum-Welch算法,对两种特殊情况下(相互独立和一致依赖),得到了改进算法的具体表达式,然后进行了具体的数值软件开发,并通过二个具体的蛋白质家族的多序列联配来说明改进算法的有效性。
其他文献
以不同倍性的亚洲百合品种为亲本进行了21组杂交,观察杂交亲本和杂交子代的染色体数目,旨在研究百合杂交亲和性与染色体倍性之间的相关性。结果发现:Renoir和Gironde为二倍体
本文主要研究人寿保险中的精算模型,对人寿保险的各个模型分别进行了分析和研究。根据人寿保险模型的几个类别详细讨论寿险品种的发展历程和寿险精算模型的发展过程。 首先
学位
本文利用非线性最优规划问题的Karush-Kuhn-Tucker(KKT)条件,研究了非线性切换系统的指数稳定性和非线性时滞切换系统的渐近稳定性.对于线性切换系统的研究主要基于不变集理论
本文讨论了交错环链补中的不可压缩、分段不可压缩曲面的性质。设L是S3中的一个交错环链,将L投影到S2上,L的每个交叉点都对应一个“bubble”,用来体现L的交叉点性质.如果L有n个
矩阵不等式作为一个广阔的数学领域,从某种意义上说不等式比等式有更大的用处。本文研究了一种矩阵块Kronecker积—Khatri-Rao积,建立了若干关于这种矩阵乘积的矩阵不等式、特