一般耦合矩阵方程组的迭代算法研究.pdfVIP

第 30卷 哈尔滨师范大学 自然科学学报 Vo1．30，No．52014 第 5期 NATURALSCIENCESJOURNALOFHARBIN NORMALUNIVERSITY 一 般耦合矩阵方程组的迭代算法研究 陶金钱 (哈尔滨工业大学) 【摘 要】通过推广共轭梯度法思想给出一种迭代算法去求解一般耦合矩阵 方程组的广义双对称解，并对算法性质给予介绍说明，将证明若一般耦合矩阵方程 组关于广义双对称解相容，那么在不考虑误差的情况下，对于任意给定的初始广义 双对称矩阵组，利用所构造出的迭代算法，都能在有限步之 内迭代得到其广义双对 称解．若取定特殊的初始矩阵，则可获得其极小Frobenius范数约束解，进一步解决 最佳逼近 问题． 【关键词】一般耦合矩阵方程组；广义双对称解；Frobenius范数；最佳逼近解 在文献[7—11]中提出了求解几类耦合矩阵方 0 引言 程组的广义中心对称解及反对称解的有效迭代 矩阵方程像 Sylvester和Lyapunov矩阵方程 算法．该文研究下述含多个矩阵变量方程组的广 在系统和控制理论方面有着重要应用，国内外许 义双对称解及最佳逼近问题，下面将文中所用术 多专家都致力于寻找它们的解析解及数值解．一 语做一简要说明：用尺 和SOR 分别表示n阶 般来说，矩阵方程 问题可通过直接法和迭代法进 实方阵和实对合阵；对任意的A，B×R ，定义内 行求解．直接法的一种思路是利用 Kronecker积 积为 (A，B)= r(rA)，由此 内积生成的矩阵 及 向量映射，通过矩阵按行拉直运算将方程组转 Frobenius范数为 ll ll= (A，B)；用Ao 表 化为线性方程组求解，这方面的突 出工作是 示它们的Kronecker积；vec(A)表示矩阵A按行 Lancaster#]1用Kronecker积和拉直映射算子得到 拉直的列向量． 了线性矩阵方程的一般解存在的充要条件和通 定义1 设P和 Q是任意的实对称正交阵， 解的表达式…．而在实际问题中，更多时候只需 若矩阵 满足PXQ=X且X=X ，则称 为关 要快速的得到满足一定条件的特解，因此迭代法 于P，Q的广义双对称矩阵，记为X ∈BSR 。． 的研究应运而生 ，并且取得了一定 的研究成果． 问题 I 给定矩阵A ∈R ，B ∈ ， 这些成果大都利用共轭梯度法、参数迭代法等算 法，解决了多种类型的矩阵方程组的一般解、对 C ∈R ，P，Q ∈SORi ，i，=1，2，…，z，求矩 称解及反对称解等 问题．陈国 良和 Xingping 阵组 ( ， ，…，Xf)，f∈BsR 使得 Sheng给出了矩阵方程组AXB =C和A：XB：= fAl1X1Bl14-A12X2B12+… +A1zXfB “=C1 C：的一种迭代解法 J．周彬和段广仁等在文献 A21X,BE1+A22X2B22+… +A21XIB2z C2 (1) [3—6]中提出了求解几类Sylvester矩阵方程组 的迭代解法，并且研究了一