非对称鞍点问题的松弛型预处理子

来源 :上海大学 | 被引量 : 2次 | 上传用户:David_storm
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在流体力学、带有限制条件的二次优化问题及电磁学等应用领域中,对所研究问题进行线性化及有限元(有限差分)离散处理后,通常都会归结为同一类大型稀疏线性方程组—鞍点问题的求解.对于鞍点问题的求解方法主要包括直接法和迭代方法两种.但是直接法在求解大规模稀疏系统时,会产生“填充”现象,运算量比较大.而迭代方法虽然在鞍点问题的求解中发挥了举足轻重的作用,但是对于像鞍点问题这样具有大规模系数矩阵的方程组而言,在实际计算过程中要经过很多步迭代才能收敛.因此,近年来对鞍点问题的求解,都是考虑采用预处理技术.本论文在潘建瑜等人2006年提出的DPSS预处理子的基础上,提出了求解非对称鞍点问题的三种松弛型预处理子,并给出了相应的理论成果及数值实验,具体研究工作可概括如下:1. DPSS预处理子是求解非对称鞍点问题的有效预处理子.由它进行预处理的鞍点问题,当选取合适迭代参数时,系数矩阵的特征值将趋于原点和(2,0).本文通过删除DPSS预处理子差矩阵的第一个分解矩阵中的一个对角位移项,而得到求解非对称鞍点问题的VDPSS预处理子.相比较DPSS预处理子,VDPSS预处理子更加贴近原鞍点问题的系数矩阵.且由该预处理子预处理的系数矩阵一定会有n个特征值为1.本论文不仅从理论上分析了预处理方程组的Krylov子空间情况,还通过数值实例来验证VDPSS预处理子在特征值分布、迭代步数以及CPU时间方面的优越性.2.通过直接删除DPSS预处理子和原鞍点问题系数矩阵的差矩阵中一个位移项,本文提出了求解非对称鞍点问题的第二种松弛型预处理子,即RDPSS预处理子.文中不仅证明了预处理后系数矩阵也具有n个特征值为1,还分析了在迭代参数趋于0和+∞时预处理后系数矩阵的剩余特征值均趋于0.更重要的是本文还证明了RDPSS预处理子对应的迭代方法无条件收敛,且对预处理后方程组的Krylov子空间情况也进行了分析.数值实例表明,RDPSS预处理子无论是从特征值分布方面,还是迭代步数,亦或是CPU时间方面,都能更加有效求解鞍点问题,且该预处理子对迭代参数的敏感性不强.3.由于DPSS预处理子与原鞍点问题系数矩阵的差矩阵中存在平衡α和a-1的问题,本文在前两种松弛型预处理子的基础上将差矩阵中对角线上的位移项全部换为零矩阵,提出了求解非对称鞍点问题的一种新的松弛型预处理子.该种新松弛型预处理子避免了需要在α和α-1之间平衡选择a.并且本文还证明了对应迭代方法的无条件收敛性,分析了预处理后系数矩阵特征值的分布情况,且理论上得出了保证迭代收敛最快的最优参数的选取.本文还对预处理后系数矩阵的特征向量的结构进行了分析,因为迭代方法的收敛速度不仅和谱性质有关,还与特征向量结构有关.针对计算过程中需要求解系数矩阵中含有A-1的方程组,文中提出了一种有效的非精确近似计算方法.数值实验表明该新松弛型预处理子相比其他松弛型预处理子在求解非对称鞍点问题时更加有效.其非精确近似计算方法虽然在迭代步数上增加了,但是CPU时间却大大减少了,因此是新松弛型预处理子的一个有效近似.
其他文献
离散可积系统的变换和约化对于探讨方程间的联系以及构造精确解等具有重要意义.本文分为以下三部分进行讨论.第一,利用ABS链方程的分解构造Backlund变换.讨论了 ABS链方程的分解性质,并按照两个方向的平移将其分解,通过选取不同的函数h来构造BT.当h为仿线性时,给出了所有满足多维相容性的相容三重组.当h非仿线性时发现该BT可以用来构造仿线性方程与多二次方程的联系.作为BT的一些应用,以Q1方程
我国猕猴桃面积和产量均居于世界第一位。根据中国农业农村部最新数据统计,贫困地区猕猴桃栽培面积达到11.14万hm2,占全国猕猴桃总规模的46.4%,已成为助推脱贫攻坚和产业兴旺的重要产业之一。十堰市处于秦岭南麓,已选育出‘汉美’‘武当1号’[1]等猕猴桃品种,是猕猴桃的适宜生长区。随着"十三五"脱贫攻坚的完成,在"十四五"期间如何实现乡村振兴,
期刊
作为典型的拓扑密堆相(TCP),σ相具有较宽的溶解度范围,并存在大量的实验测量数据,如晶格常数、占位分数等,可以作为TCP的研究原型。本文通过使用第一性原理计算、CALPHAD方法并结合文献中的实验数据,对二元系σ相的物理性质进行了系统的研究(包括原子分布、摩尔体积、形成焓和体弹模量)。首先,我们发现σ相的原子分布(即组成原子在不同晶体结构位置的分布规律或占位优先顺序)受到尺寸因素和组成元素电子排
基于光正交频分复用(OOFDM)技术的无源光网络以其卓越的色散容忍能力、频谱效率高等特点被认为是下一代光接入网的候选技术,受到广泛关注。如何有效降低OFDM-PON ONU端独立光源成本、提高带宽资源利用率,是该技术实用化的关键难点之一。本文围绕低成本RSOA强度调制器在OFDM-PON系统中的应用进行研究,提出了两种提升系统性能的方案,进行了数值仿真分析和实验验证。最后针对OFDM-PON系统中
随着技术的发展,人们可以用3D设备快速获取真实场景的彩色点云数据,并在建模、仿真、漫游等方面获得广泛应用。但面对高精度,数量庞大的点云数据,如何寻求有效的方法进行压缩,以节省存贮空间,实现点云有效传输和实时的渲染,是目前亟待解决的问题;而对采集、传输、压缩及复原过程中引起的误差和失真,也需要给出一个客观质量评价,以衡量各种数据处理算法性能的优劣。因此,研究彩色点云压缩和质量评价算法,有着重要的理论
由于网络资源和负载的快速增长及分布不均衡,互联网“尽力而为”的服务模式决定了网络拥塞无法避免,其不会因为网络处理能力的提高及网络链路带宽的增加而消除,故有效的网络拥塞控制机制对保障网络的QoS(Quality of Service,服务质量)至关重要。因此,在过去二十年,网络拥塞控制问题研究一直是非常重要且极具挑战性的研究热点,吸引了全世界众多科研人员的持续关注并取得一些研究成果,但尚有许多方面需
在有限群论的研究中,子群的正规性与子群间的某种交换性是人们研究的基本出发点,而子群的正规化子与中心化子是子群的正规性和交换性的一种度量.于是利用子群的正规化子和中心化子来研究有限群的结构成为人们非常感兴趣的研究课题,并且获得了大量的研究成果.本文也将从子群的正规化子和中心化子出发来研究有限p-群的结构,同时我们对秩为2的有限2-群也做了进一步的研究.第三章研究了2平衡p-群的结构.称有限群G为n平
本文采用定向凝固手段,以Fe-Ni、Pb-Bi、Cu-Sn和Zn-Cu四种包晶合金为研究对象,系统研究了稳恒磁场(以下简称磁场)对包晶合金定向凝固组织的影响。模拟了磁场下包晶合金定向凝固过程中不同尺度上热电磁流动和热电磁力,考察了横向弱磁场对包晶合金带状组织、岛状组织、偏析、枝晶生长和取向的影响,考察了纵向强磁场对包晶合金宏观偏析、固态相变和枝晶生长的影响。构建了磁场下热电磁效应三维数值模型,研究
近年来,稀疏性问题在科学研究和工程实践等诸多领域里引起了广泛的重视.特别是在压缩感知理论被提出之后,使得稀疏性问题在信号恢复,图像处理以及统计推断等领域得到了大量的应用.因此,稀疏性问题的研究具有重要的理论和现实意义.本文主要是基于(?)q正则化的稀疏优化问题研究.当0
Kirchhoff型微分方程是德国物理学家G. Kirchhoff于1883年研究弦振动时提出的一种模型,它修正了经典的达朗贝尔波动方程,从而更加精确地描述了弦振动的过程.这个模型在非牛顿流体力学、天体物理、图像处理、血浆问题和弹性理论等诸多领域都有广泛应用Orlicz-Sobolev (Musielak-Orlicz-Sobolev)空间克服了Sobolev空间不能处理非齐次算子的缺陷,为上述非