两类线性方程组的随机迭代算法及化学主方程的反位移Arnoldi算法

来源 :上海大学 | 被引量 : 0次 | 上传用户:willian1019
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大型稀疏线性方程组的求解和大型矩阵指数函数的计算一直是很多科学计算和人工智能领域的核心问题,构造这两类问题的高效算法也一直是数值代数领域的研究热点之一.本文一方面将贪婪的随机Kaczmarz(GRK)算法和随机Gauss-Seidel(GRGS)算法应用于岭回归和分解线性系统的求解当中,构造了这两类线性系统的松驰型随机迭代算法;另一方面,利用反位移技术和重正交化的Arnoldi过程,提出了一种新的逼近化学主方程精确解(矩阵指数函数与向量的乘积)的数值算法(SIRA).主要创新工作包括:1.对岭回归问题,考虑到它的正规方程组的系数矩阵是对称正定的,本文首先利用对称正定矩阵的性质修改GRK算法的概率准则和迭代式,构造了求解岭回归问题的变式GRK算法并分析了其收敛性.其次,通过在变式GRK算法的迭代式中引入位于区间(0,2)内的松弛参数,得到了松驰型的变式GRK算法并证明了它的收敛性.最后,为了最大限度利用迭代过程中所计算出的正规方程组的残差信息,本文又提出了松弛型的变式GRK算法的加速迭代格式.数值实验表明,本文提出的三种算法的收敛速率明显比文献[51]中方法快,并且含参数的变式GRK算法及其加速格式所用CPU时间也要少得多.2.RK-RK算法和REK-RK算法是目前求分解线性系统的两种最新的随机迭代算法.考虑到GRK算法和GRGS算法的快速收敛性,本文分别构造了针对相容和不相容分解线性系统的松弛型的GRK-GRK算法和GRGS-GRK算法并给出了这两种算法的收敛性分析.数值实验表明,对于分解线性系统,本文所提出的算法在迭代步和CPU时间方面明显优于RK-RK算法和REK-RK算法.3.FSP算法和Krylov FSP算法是目前求解化学主方程(CME)的两种经典的降阶算法.本文利用反位移技术和重正交化的Arnoldi过程,构造了逼近化学主方程精确解的SIRA算法,该算法不需要确定初始有限状态投影集合及状态投影展开方法,计算过程简单.通过对具体的生化反应系统模型进行实验,结果表明对中等规模的化学主方程来说,SIRA算法比FSP算法和Krylov FSP算法的精度都高.
其他文献
冶金熔体物化性质直接影响着冶炼工艺的优化及产品质量的控制。由于高温实验困难,因此在有限实验基础上的模型预报是获取此类性质数据的重要方法。目前,多元熔体物化性质的计算主要依靠唯象模型,其中,几何模型因其简洁有效的形式而成为应用广泛的一类唯象模型。几何模型主要是根据多元系边界性质对其内部性质进行预测的一类方法,其发展至今主要面临三大难题:a)如何发展边界拟合方法及权重分配方式?b)如何发展局部互溶几何
现代设计在不断渗透人们生活和工作的过程中,展现出其强大的影响力,而中国平面设计深受现代设计思潮的影响,无论是设计硬软件,还是设计理念或者模式,都不同程度地受到西方设计思维和机制的影响。显然在满足人们平面审美诉求的过程中,如果仅仅停留于模仿,往往难以达到理想的状态,此时就需要切实地思考这样的问题:如何将优秀传统文化元素融入到中国平面设计中去。从这个角度入手,首先思考现代设计思潮对于中国平面设计的影响
热电材料是一种能够实现热能与电能直接转换的材料,可以利用其Seebeck效应进行温差发电,也可以利用其Peltier效应实现制冷和制热,在废热发电、电子器件制冷领域具有重要的应用价值。Bi2Te3基化合物,在室温区具有优异的热电性能,已经广泛应用于商用制冷器件。通过掺杂、固溶、纳米化等手段,可以大幅提升Bi2Te3基块体材料的热电性能。作为一种成熟的商用热电材料,需要具备优异的热电性能、较短的制作
热电材料是利用塞贝克效应与帕尔帖效应实现电能与热能相互转换的功能材料,在废热回收再利用和固态制冷领域有着广泛的应用。热电材料的能量转换效率正比于无量纲热电优值z T,而z T由材料内部的电子和声子的输运过程所共同决定。因此,提高热电材料的z T值需要电声输运的协同优化。然而,决定电声输运的三个主要性能参数(电导率,塞贝克系数,热导率)相互耦合作用强烈,利用传统优化手段针对某一个参数进行优化通常会导
针对低品位、复合共伴生矿产资源利用率低的现状,应用“低温还原-磁选或熔分”思路,将铁与其他有价元素进行分离。为了提高低温还原的反应效率,本文首次提出将稳恒磁场作用于铁矿还原过程,以期实现磁场对铁矿还原的强化作用,为开发低品位、复合共伴生铁矿低温磁场强化还原技术提供理论基础与技术支持。基于经典热力学模型、固态相变理论及磁能理论,将磁能引入含铁矿物固态还原过程,提出了磁能反应热力学模型、磁能形核热力学
水资源危机是本世纪人类面临的最大危机之一。为了应对水资源危机,大力发展水处理技术成为必然选择之一。电容去离子技术是近些年来新兴起的水处理技术,具有低能耗、高效率、无污染、可循环使用的优点,因此有很大的发展前景。电容去离子技术的核心是电极材料,现有的电极材料主要包括活性炭、碳气凝胶等,但是它们的容量有限,性能不能满足大规模工业化生产的要求。石墨烯作为一种新型碳材料,具有高导电性和高理论比表面积的优点
微分动力系统分析方法的理论研究及其在其它领域中的应用,是应用数学中的两个重要内容.本博士学位论文就在动力系统的理论和应用两个方面进行了研究,主要做了三方面的工作:建立且分析了一个有向关联网络上的传染病模型,推广了平面动力系统中的旋转数理论,并且利用Conley指标理论研究了双吸引子间的动力学结构.文中的主要内容概述如下:首先,越来越多的迹象表明,许多真实的网络之间是相互作用的,并且有时网络上的疾病
社区结构(Community Structure)是复杂网络的一个重要特征,它具有社区内部节点连接相对紧密而社区之间节点连接相对稀疏的特点。发现网络中的社区结构对理解网络拓扑结构和挖掘网络底层隐含信息起着重要作用。因此,在复杂网络中进行社区发现算法研究成为近年来的研究热点之一。传统的社区发现算法主要研究非重叠社区结构的网络,即网络中的节点能且能归属于一个社区。然而,在现实世界的网络中,大多数实体都
本论文的研究分为两部分.论文第一部分主要研究求解大型稀疏线性方程组的随机Kaczmarz方法的相关数值算法问题,并将其应用到压缩感知、信号重构等问题中.所做工作的相关内容概括如下:1.通过对Kaczmarz算法的研究,利用随机抽样的方法来计算部分残差,作者得到了一种新的求解大型线性系统的数值迭代方法.根据每一次随机选取部分行计算其残差,并选取残差最大的一行来进行Kaczmarz迭代,作者构造了随机
输液管道广泛应用于航空、航天、机械以及能源等各个领域。由于流体流动和外界环境引发的管道振动会对整个结构产生严重的影响,因此对输液管道振动的研究有重要的工程应用价值。同时,输液管道作为一种典型的陀螺连续系统,对其振动特性的分析也有重要的理论意义。当流速处于亚临界范围,管道在零平衡位形附近振动。当流速进入超临界领域,管道在非平凡静平衡位形附近振动。本文基于Timoshenko梁理论和广义Hamilto