对流扩散方程最优控制问题中线性方程组的预处理算法

来源 :华东师范大学 | 被引量 : 0次 | 上传用户:alimamaai
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本论文主要讨论对流扩散方程最优控制问题中线性方程组的数值求解.对此类问题进行求解时,需要将解优化过程和解状态方程统一结合起来.目前主要的处理方法有先离散后优化和先优化后离散两种途径.无论采用哪种方法,最后都需要求解一个大规模的线性方程组.对于这样的大规模线性方程组,我们通常采用Krylov子空间迭代方法进行求解.当系数矩阵的条件数或谱分布比较差时,迭代法(如:共轭梯度法,广义最小残量法等)可能收敛会很慢,有时甚至会不收敛.为了改善迭代法的收敛性质,我们可以采用预处理技术.在科学计算领域中,预处理是当今开发有效求解器的关键技术之一,好的预处理子不仅能够提高迭代算法的收敛速度,还可以改善算法的稳定性.对流扩散方程最优控制问题通过先离散后优化和先优化后离散方法处理后,将分别产生两类不同的线性方程组,其中一类是对称的,而另一类是非对称的.本文主要讨论这两类线性方程组的预处理方法.具体工作如下:(1)对于先离散后优化的情形,我们分别讨论了当正则化参数较小时和正则化参数较大时的预处理子选取问题,包括块对角,块三角和约束型预处理子的构造.在此基础上,我们提出了两类新的预处理子.(2)对于先优化后离散的情形,我们给出了相应的块对角和块三角预处理子,并提出了一类基于不完全分解的新预处理子.(3)利用系数矩阵的特殊结构,我们对原线性方程组进行了降阶处理,将其转化为一个等价的但规模较小的线性方程组,然后讨论了其预处理方法.(4)最后,我们对提出的预处理子进行数值试验,验证了这些预处理方法的有效性.
其他文献
病毒传播和信息传播是复杂网络中十分重要的研究课题。目前,很多工作解释了信息传播迅速的原因以及如何加速传播。在实际的网络中,有一类节点对于传播的过程有着十分重要的作
本文主要考虑了边传递图直积的独立数也满足Tardif问题等式和一般图的直积与其对应线图的直积独立数分别同时都满足Tardif问题等式的条件.文章正文由三章组成.第一章主要介绍
本课题由国家自然科学基金项目“小功率微波微等离子体的研究”(批准号:61072007)资助。小功率平面微带微波微等离子体是一项基于微波和微电子以及微等离子体的高新技术。在18
本论文研究了基于微悬臂梁的微光腔中的光热耦合机制中光热冷却效率的温度依赖性。我们的实验研究以磁共振力显微镜(MRFM)的实验平台作为基础,搭建了适合于研究光力学中耦合
特定mRNA的输出、转录和翻译是基因表达的必要前提。然而,这些连续步骤之间的功能联系尚不清楚。在酵母及人细胞中,mRNA输出因子Glel对mRNA从细胞核到细胞质的输出是必需的,
一个次数d≥2的代数整数a,若α>max2≤i≤d|αi|其中αi(2≤i≤d)为α的除它自身外的所有共轭元,那么称这个代数整数α为Perron数.关于最小Perron数,有著名的Lind-Boyd猜想[1
这是一篇关于最小包含球问题及其相关问题的综述.最小包含球问题早在19世纪由Sylvester根据这样的背景提出:(1)社区医院问题:我们将社区中的每户人家看作平面中的一个点,为了社
设H是复可分Hilbert空间,B(H)是H上所有的有界线性算子组成的代数.在本文中,我们研究具有闭的数值域的算子的紧扰动问题,证明了具有闭的数值域的算子在B(H)中是范数拓扑下稠
本文研究具有测量误差的自回归模型统计推断问题.我们主要考虑可观测的变量(即带有误差的变量,也称替代变量)与准确真实变量之间的函数关系是未知的,在这种情况下对自回归模
在基础研究领域内,例如离子在加速器内的运输和寿命问题,离子与物质相互作用机制,原子光谱,天体物理及等离子体物理等许多方面,离子原子碰撞过程机制起到非常重要的作用。其