一种求解分布控制问题的含参预处理最小残量法(英文)

近年来,提出了一些求解约束优化分布控制问题的方法,其中最常用的方法是先离散偏微分方程,然后求解离散得到的线性方程组.文献中提出了一些Krylov子空间预处理方法用来求解该线性方程组.通过分析张晓莹等提出的块对角预处理矩阵(Zhang X Y,Yan H Y.Huang Y M.On preconditionedMINRES method for solving the distributed control problems.Commun Appl Math Comput,2014,28:128-132.),构造了一个含参数的块对角预处理线性方程组,并运用含参数预处理最小残量方法求解该线性方程组.预处理矩阵的谱分析表明当参数大于1时,含参预处理线性方程组的谱分布更加集中.数值实验结果验证了含参数的预处理最小残量方法对于求解分布式控制问题是有效的.     

关于“求解加权线性最小二乘问题的一类预处理GAOR方法”一文的注记

在"求解加权线性最小二乘问题的一类预处理GAOR方法"一文中,作者提出了求解加权线性最小二乘问题等价$2\times 2$块线性系统的一类预处理GAOR方法,并给出了几个比较定理来说明新提出预处理GAOR方法的优越性.本文我们将指出该文中几个比较定理的不完善之处和证明的错误之处,并给出正确的证明.     

GSOR方法最优迭代参数的简化推导

对于增广线性系统,Bai等研究了广义SOR方法（Bai Z Z,Parlett B,Wang Z Q.On generaliged successive overrelaxation methods for augmented linear systems.Numerische Mathematik,2005,102（1）：1-38）,并得到其最优迭代参数.给出了另外一种推导最优迭代参数的简化方法,这种方法对于求解其他参数加速定常迭代方法的最优迭代参数非常有意义.     

求解一类复对称线性系统的新的PMHSS迭代法（英文）

在科学计算及工程应用中经常遇到复对称线性系统问题,近年来对一种特殊类型的复对称线性系统的研究已成为一个热点.基于白中治等的PMHSS方法(Bai Z Z,Benzi M,Chen F,Wang Z Q.Preconditioned MHSS iteration methods for a class of block twoby-two linear systems with applications to distributed control problems.IMA J Numer Anal,2013,33:343-369),提出一类新的PMHSS迭代法用于求解这种特殊形式的复对称线性系统,给出新方法的收敛性理论以及最优参数的表达式,最后用数值例子展示了新方法的有效性.     

一类Toeplitz线性代数方程组的预处理GMRES方法

本文研究一类来源于分数阶特征值问题的Toeplitz线性代数方程组的求解.构造Strang循环矩阵作为预处理矩阵来求解该Toeplitz线性代数方程组,分析了预处理后系数矩阵的特征值性质.提出求解该线性代数方程组的预处理广义极小残量法（PGMRES）,并给出该算法的计算量.数值算例表明了该方法的有效性.     

求解加权线性最小二乘问题的一类预处理GAOR方法

为了快速求解一类来自加权线性最小二乘问题的2×2块线性系统,本文提出一类新的预处理子用以加速GAOR方法,也就是新的预处理GAOR方法.得到了一些比较结果,这些结果表明当GAOR方法收敛时,新方法比原GAOR方法和之前的一些预处理GAOR方法有更好的收敛性.而且,数值算例也验证了新预处理子的有效性．     

二维Helmholtz方程的联合紧致差分离散方程组的预处理方法

对于二维的Helmholtz方程,本文用联合紧致差分格式(CCD)离散,该差分格式具有六阶精度,三点差分和隐式的特点.本文基于CCD格式离散得到的线性系统和循环矩阵的快速傅里叶变换,提出了一种循环型预处理算子用于广义极小残量迭代算法(GMRES).给出了循环型预处理子的求解算法,证明了该预处理算子能使迭代算法具有较快的收敛速度.本文还与其他算法的预处理算子作比较,数值结果表明本文提出的循环型预处理算子具有更好的稳定性,并且对于较大的波数k,收敛速度也更快.     

刚性目标形状反演的一种非线性最优化方法

发展了从声散射场的远场分布的信息来再现声刚性目标形状反问题的一种非线性最优化方法,它是通过独立地求解一个不适定的线性系统和一个适定的非线性最小化问题来实现的。对反问题的非线性和不适定性的这种分离式数值处理,使所建立方法的数值实现是非常容易和快速的,因为在确定声刚性障碍物形状的非线性最优化步中,只需求解一个只有一个未知函数的小规模的最小平方问题。该方法的另一个特别的性质是,只需要远场分布的一个Fourier系数,即可对未知的刚性目标作物形设别。进而提出了数值实现该方法的一种两步调整迭代算法。对具有各种形状的二维刚性障碍物的数值试验保证了本算法是有效和实用的。     

求解M-张量方程的两种新型算法

多重线性系统在当今的工程计算和数据挖掘等领域有很多实际应用,许多问题可以转化为多重线性系统求解问题.在本文中,我们首先提出了一种新的迭代算法来求解系数张量为M-张量的多重线性系统,在此基础上又提出了一种新的改进算法,并对两种算法的收敛性进行了分析.数值算例的结果表明,本文提出的两种算法是有效的并且改进算法的迭代时间更少.     

求解一类复对称线性系统的改进的SNS和SSS迭代法

本文在Bai的基础上提出改进的斜正规分裂(MSNS)和斜尺度化分裂(MSSS)迭代法,用以求解一类应用广泛的复对称线性系统,并证实MSNS和MSSS迭代法是无条件收敛的.通过利用一些Krylov子空间方法,本文给出相对应的非精确版本的MSNS(MSSS)方法.数值实验说明了所给方法的有效性.     

Robust Iterative Solution of a Class of Time-Dependent Optimal Control Problems

The fast iterative solution of optimal control problems, and in particular PDE-constrained optimization problems, has become an active area of research in applied mathematics and numerical analysis. In this paper, we consider the solution of a class of time-dependent PDE-constrained optimization problems, specifically the distributed control of the heat equation. We develop a strategy to approximate the (1, 1)-block and Schur complement of the saddle point system that results from solving this problem, and therefore derive a block diagonal preconditioner to be used within the MINRES algorithm. We present numerical results to demonstrate that this approach yields a robust solver with respect to step-size and regularization parameter. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

椭圆PDE-约束优化问题的一个预条件子

针对由Galerkin有限元离散椭圆PDE-约束优化问题产生的具有特殊结构的3×3块线性鞍点系统,提出了一个预条件子并给出了预处理矩阵特征值及特征向量的具体表达形式.数值结果表明了该预条件子能够有效地加速Krylov子空间方法的收敛速率,同时也验证了理论结果.     

Additive block diagonal preconditioning for block two-by-two linear systems of skew-Hamiltonian coefficient matrices

For a class of block two-by-two systems of linear equations with certain skew-Hamiltonian coefficient matrices, we construct additive block diagonal preconditioning matrices and discuss the eigen-properties of the corresponding preconditioned matrices. The additive block diagonal preconditioners can be employed to accelerate the convergence rates of Krylov subspace iteration methods such as MINRES and GMRES. Numerical experiments show that MINRES preconditioned by the exact and the inexact additive block diagonal preconditioners are effective, robust and scalable solvers for the block two-by-two linear systems arising from the Galerkin finite-element discretizations of a class of distributed control problems.     

On a type of matrix splitting preconditioners for a class of block two-by-two linear systems

Recently, Bai et al. (2013) proposed an effective and efficient matrix splitting iterative method, called preconditioned modified Hermitian/skew-Hermitian splitting (PMHSS) iteration method, for two-by-two block linear systems of equations. The eigenvalue distribution of the iterative matrix suggests that the splitting matrix could be advantageously used as a preconditioner. In this study, the CGNR method is utilized for solving the PMHSS preconditioned linear systems, and the performance of the method is considered by estimating the condition number of the normal equations. Furthermore, the proposed method is compared with other PMHSS preconditioned Krylov subspace methods by solving linear systems arising in complex partial differential equations and a distributed control problem. The numerical results demonstrate the difference in the performance of the methods under consideration.     

Right preconditioned MINRES for singular systems†

We consider solving large sparse symmetric singular linear systems. We first introduce an algorithm for right preconditioned minimum residual (MINRES) and prove that its iterates converge to the preconditioner weighted least squares solution without breakdown for an arbitrary right‐hand‐side vector and an arbitrary initial vector even if the linear system is singular and inconsistent. For the special case when the system is consistent, we prove that the iterates converge to a min‐norm solution with respect to the preconditioner if the initial vector is in the range space of the right preconditioned coefficient matrix. Furthermore, we propose a right preconditioned MINRES using symmetric successive over‐relaxation (SSOR) with Eisenstat's trick. Some numerical experiments on semidefinite systems in electromagnetic analysis and so forth indicate that the method is efficient and robust. Finally, we show that the residual norm can be further reduced by restarting the iterations.     

A preconditioning technique for a class of PDE-constrained optimization problems

We investigate the use of a preconditioning technique for solving linear systems of saddle point type arising from the application of an inexact Gauss?CNewton scheme to PDE-constrained optimization problems with a hyperbolic constraint. The preconditioner is of block triangular form and involves diagonal perturbations of the (approximate) Hessian to insure nonsingularity and an approximate Schur complement. We establish some properties of the preconditioned saddle point systems and we present the results of numerical experiments illustrating the performance of the preconditioner on a model problem motivated by image registration.     

A generalization of parameterized inexact Uzawa method for generalized saddle point problems

Recently, a class of parameterized inexact Uzawa methods has been proposed for generalized saddle point problems by Bai and Wang [Z.-Z. Bai, Z.-Q. Wang, On parameterized inexact Uzawa methods for generalized saddle point problems, Linear Algebra Appl. 428 (2008) 2900–2932], and a generalization of the inexact parameterized Uzawa method has been studied for augmented linear systems by Chen and Jiang [F. Chen, Y.-L. Jiang, A generalization of the inexact parameterized Uzawa methods for saddle point problems, Appl. Math. Comput. (2008)]. This paper is concerned about a generalization of the parameterized inexact Uzawa method for solving the generalized saddle point problems with nonzero (2, 2) blocks. Some new iterative methods are presented and their convergence are studied in depth. By choosing different parameter matrices, we derive a series of existing and new iterative methods, including the preconditioned Uzawa method, the inexact Uzawa method, the SOR-like method, the GSOR method, the GIAOR method, the PIU method, the APIU method and so on. Numerical experiments are used to demonstrate the feasibility and effectiveness of the generalized parameterized inexact Uzawa methods.     

关于具优势对称部分的不定线性代数方程组的分裂极小残量算法

１．引 言 考虑大型稀疏线性代数方程组 为利用系数矩阵的稀疏结构以尽可能减少存储空间和计算开销，Ｋｒｙｌｏｖ子空间迭代算法［１，１６，２３］及其预处理变型［６，８，１３，１８，１９］通常是求解（１）的有效而实用的方法．当系数矩阵对称正定时，共轭梯度法（ＣＧ（     

Preconditioned iterative methods for diffusion problems with high‐contrast inclusions

This paper is concerned with robust numerical treatment of an elliptic PDE with high‐contrast coefficients, for which classical finite‐element discretizations yield ill‐conditioned linear systems. This paper introduces a procedure by which the discrete system obtained from a linear finite element discretization of the given continuum problem is converted into an equivalent linear system of the saddle‐point type. Three preconditioned iterative procedures—preconditioned Uzawa, preconditioned Lanczos, and preconditioned conjugate gradient for the square of the matrix—are discussed for a special type of the application, namely, highly conducting particles distributed in the domain. Robust preconditioners for solving the derived saddle‐point problem are proposed and investigated. Robustness with respect to the contrast parameter and the discretization scale is also justified. Numerical examples support theoretical results and demonstrate independence of the number of iterations of the proposed iterative schemes on the contrast in parameters of the problem and the mesh size.     

Equivalent operator preconditioning for elliptic problems with nonhomogeneous mixed boundary conditions

The numerical solution of linear elliptic partial differential equations often involves finite element discretization, where the discretized system is usually solved by some conjugate gradient method. The crucial point in the solution of the obtained discretized system is a reliable preconditioning, that is to keep the condition number of the systems under control, no matter how the mesh parameter is chosen. The PCG method is applied to solving convection-diffusion equations with nonhomogeneous mixed boundary conditions. Using the approach of equivalent and compact-equivalent operators in Hilbert space, it is shown that for a wide class of elliptic problems the superlinear convergence of the obtained preconditioned CGM is mesh independent under FEM discretization.