一类特殊矩阵方程的并行预处理变形共轭梯度算法

研究了求解一类矩阵方程AXB=C,提出了一种并行预处理变形共轭梯度法．该方法给出一种迭代法的预处理模式．首先给出的预处理矩阵是严格对角占优矩阵,构造并行迭代求解预处理矩阵方程的迭代格式,进而使用变形共轭梯度法并行求解．通过数值试验,预处理变形共轭梯度法与直接使用变形共轭梯度法相比较,该算法不仅有效提高了收敛速度,而且具有很高的并行性．     

一种求解鞍点问题的预处理并行算法

研究了一种求解鞍点问题的并行预处理变形共轭梯度算法.通过应用迭代法进行预处理后,再采用变形共轭梯度求解的模式.首先构造系数矩阵近似逆的多项式表达式,以此作为预处理矩阵的逆矩阵,对方程组进行预处理;然后采用变形共轭梯度法并行求解预处理后的线性方程组.为减少运算量,采用迭代方式并行计算多项式与向量的乘法运算.通过调整迭代次数,即调整多项式次数,检验各种次数的多项式进行预处理后的求解方程的效果.数值试验结果表明,该算法明显优于未预处理的变形共轭梯度法,且当预处理迭代次数取4时效果最好.     

求解Sylvester矩阵方程的一种改进的梯度方法

提出了一种改进的梯度迭代算法来求解Sylvester矩阵方程和Lyapunov矩阵方程.该梯度算法是通过构造一种特殊的矩阵分裂,综合利用Jaucobi迭代算法和梯度迭代算法的求解思路.与已知的梯度算法相比,提高了算法的迭代效率.同时研究了该算法在满足初始条件下的收敛性.数值算例验证了该算法的有效性.     

求解混合型Lyapunov矩阵方程的参数迭代法

采用参数迭代法求一类混合型Lyapunov矩阵方程A~TX XA B~TXB=C的对称解.在方程相容的条件下,给出了迭代法收敛的充要条件和一些充分条件,以及参数的选取方法.最后,利用数值算例对有关结果进行了验证.     

一类广义Lyapunov矩阵方程的正定解

研究了双线性系统中的一类广义Lyapunov矩阵方程的正定解.基于混合单调算子不动点定理,给出新的存在正定解的充分条件,构造了求其正定解的不动点迭代方法,并给出了迭代误差估计公式.数值实验表明新方法是可行的.     

关于部分变元Lyapunov矩阵方程几个问题的讨论

关于Lyapunov矩阵方程A^TB+BA=-C的解与线性定常系统x=Ax之零解的部分变元渐近稳定性的关系,本文就最近发表的一些结果讨论了如下几个问题。一、由于全变元正定函数也满足部分正定性的条件,有必要引进严格部分正定函数的定义。严格部分正定函数与全变元正定函数是互不相包的;二、求解矩阵方程即意味着对于给定的矩连A(它使系统x=Ax之零解对部分变元渐近稳定),矩阵C应满足什么条件使矩阵方程有解B,此即Lyapunov函数的存在与构造问题;本文还指出C的秩不必为m,当rank C>m时     

求解带Toeplitz矩阵的线性互补问题的一类预处理模系矩阵分裂迭代法

针对系数矩阵为对称正定Toeplitz矩阵的线性互补问题,本文提出了一类预处理模系矩阵分裂迭代方法.先通过变量替换将线性互补问题转化为一类非线性方程组,然后选取Strang或T.Chan循环矩阵作为预优矩阵,利用共轭梯度法进行求解.我们分析了该方法的收敛性.数值实验表明,该方法是高效可行的.     

摄动离散矩阵Lyapunov方程向后误差分析

研究摄动离散矩阵Lyapunov方程解的向后误差,利用矩阵Kronecker积的性质以及矩阵范数的性质,给出方程近似解的向后误差界,最后通过数值例子说明解的向后稳定性.     

关于H-矩阵的H-预处理子（英文）

设A为一实对称正定的严格对角占优矩阵.设A=D-B为A的Jacobi分裂.为了求解线性方程组Ax=b,在新提出的预处理子的基础上,我们采用预处理共轭梯度方法(PCG)来求解该问题.新提出的预处理子Pv=D+νvv~T,其中v=|B|e,e=(1,...,1)~T,ν=v~TBv/||v||_2~4,且ν使||cvv~T-B||_F达到极小.我们得到了预处理矩阵P_v~(-1)A特征值的上下界,它的界比JIN提出的预处理子的界简单紧凑.数值结果表明我们的预处理子的有效性.     

几类典型矩阵方程的梯度矩阵的计算

应用矩阵微分思想,计算了几类典型的矩阵方程的梯度矩阵并给予了证明.     

一类Lyapunov矩阵方程的双反对称解的迭代算法

基于求线性代数方程组的共轭梯度法的思想,建立一种求Lyapunov矩阵方程的双反对称解的迭代算法,对任意给定的初始双反对称矩阵,算法能够在有限步迭代计算后得到矩阵方程的极小范数双反对称解,同时在上述解集中也可得出指定矩阵的最佳逼近双反称矩阵.数值算例表明,迭代算法是有效的.     

用矩阵方程法求解循环卷积的改进算法

从循环卷积的定义出发,描述了用以m为变量的方法求解循环卷积的步骤.特别详述了其中的矩阵方程法,指出了该方法的不足——含有冗余项,提出了矩阵方程法的改进算法,消除了冗余项,简化了矩阵方程,给出了矩阵方程内部各列的构成规律.实例表明该方法简便、有效.     

Some Applications of the Lyapunov Matrix Equation

The matrix equation A'P+PA = –Q originates in the stabilityanalysis of the system of linear differential equations x =Ax by Lyapunov's direct method. Many other aspects of the stabilityof such systems and of related problems in matrix theory canalso be examined by this matrix equation. Some of these arediscussed in this paper and new applications to the stabilityof second order damped dynamic systems and to stable quasi-Jacobimatrices are given.     

利用初等变换求解线性矩阵方程

本文给出了一般线性矩阵方程ＡｍｎＸｎｓ＝Ｂｍｓ，ＸｍｎＡｎｓ＝Ｂｍｓ，ＡｍｎＸｎｓＢｓｔ＝Ｃｍｔ的解的结构定理，并介绍了一种利用初等变换求解上述三类线性矩阵方程的方法．     

求解线性规划的一种单纯形矩阵法

线性规划的单纯形法一直是运筹学教学中的难点,是求解线性规划的一种重要方法.通过实例从代数角度探讨了单纯形法的迭代思想,提出了用单纯形矩阵求解线性规划的方法.同传统的单纯形表计算比较而言,此方法操作简单,不易出错,为线性规划的求解提供了一种行之有效的方法。     

解大系统稳定性的积分方程法

本文通过分解积分方程组，并建立积分方程法比较原理，化大系统为低维系统，进而讨论了带时滞的时变大系统的稳定性，给出了新的结果，这一方法也可以用来讨论其它类型大系统的稳定性问题．     

闭凸集约束下线性矩阵方程求解的松弛交替投影算法

研究线性矩阵方程AXB=C在闭凸集合R约束下的数值迭代解法.所考虑的闭凸集合R为(1)有界矩阵集合,(2)Q-正定矩阵集合和(3)矩阵不等式解集合.构造松弛交替投影算法求解上述问题,并用算子理论证明了由该算法生成的序列具有弱收敛性.给出了矩阵方程AXB=C求对称非负解和对称半正定解的数值算例,大量数值实验验证了该算法的可行性和高效性,并说明该算法与交替投影算法和谱投影梯度算法比较在迭代效率上的明显优势.     

AOR Iterative Method for Coupled Lyapunov Matrix Equations

An AOR(Accelerated Over-Relaxation)iterative method is suggested by introducing one more parameter than SOR(Successive Over-Relaxation)method for solving couple...     

新的辅助方程法构造KdV方程的行波解

应用一种新的辅助方程法成功地获得了(1+1)维KdV方程的多个含有参数的精确行波解,所得的解涵盖了已有结果．与其它方法相比,所给出的方法具有简单高效、计算量小、速度快、易于求解等特点．另外,所给的方法还可以用来求解其它的一大类非线性发展方程的精确行波解．