矩阵方程AXB+CYD=E的M对称解的迭代算法.pdfVIP

2015年 5月 计 算 数 学 第 37卷第 2期 May,2015 MATHEM ATICA NUMERICA SINICA Vo1．37，No．2 矩阵方程 AXB+CYD= 的 对称解的迭代算法术 周海林 (南京理工大学泰州科技学院，江苏泰州225300) 摘 要 在共轭梯度思想的启发下，结合线性投影算子，给出迭代算法求解了线性矩阵方程 AXB+ CryD ：E的M 对称解 ，y】及其最佳逼近．当矩阵方程AXB+ yD=E有M 对称解 时，应用迭代算法，在有限的误差范围内，对任意初始M 对称矩阵对 [X1，Y1]，经过有限步迭代可 得到矩阵方程的M 对称解；选取合适的初始迭代矩阵，还可得到极小范数 M 对称解．而且，对 任意给定的矩阵对 ，y]，矩阵方程A B+ yD=E的最佳逼近可以通过迭代求解新的矩阵 方程 AXB+ayD =E 的极小范数 M 对称解得到．文中的数值例子证实了该算法的有效性． 关键词：共轭梯度；投影算子；M 对称解；极小范数 M 对称解；最佳逼近 MR (2000)主题分类：65F10 1．引 言 记 R 表示m ×n阶实矩阵集合，厶表示 n阶单位矩阵， 表示 的转置矩阵．对 任意的矩阵 A，BER ，定义矩阵内积 (，B)=tr(BT )，则由它诱导的范数为Frobenius 范数，即IIAII=、／／(A，A)=tr[(A丁’A)]． 定义 1．设矩阵MER ，若矩阵XER 满足 (M XM)T=MTXM，则称矩阵 为n阶M 对称矩阵；若矩阵 XER 满足 (M XM)T=一M M，则称矩阵 为n阶 M 反对称矩阵． 显然，M 对称矩阵 (M 反对称矩阵)是一种广义对称 (反对称)矩阵．记SRF 为全体n 阶M 对称矩阵集合，ASRM~~ 为全体 n阶 M 反对称矩阵集合． 定义 2．对矩阵M，NER ，若 (M，N)=0，则称矩阵M，Ⅳ相互正交． 本文主要考虑下面两个 问题． 问题I．给定矩阵 ，C∈R ，B，D ∈R P，E∈R P，求 ，y】，其中X，YEsR ， 使得 AXB+CYD=E． (1) 问题II．当问题I相容时，记其解集合为 SE，对给定的 [， ， ， ∈R ，求 [， ∈ SE，使得 lI 一 ll。+ll一 l『 m in {ll 一 ll+lly一 ll。)． (2) ，y1s 约束矩阵方程 问题是指求矩阵方程具有一定约束条件的矩阵解．约束条件不同或方程不 同，就得到不同的约束矩阵方程 问题．矩阵方程的求解具有重要的理论意义和很高的应用价 2014年 9月 17日收到 2期 周海林：矩阵方程 A B+CyD=E 的M 对称解的迭代算法 187 值．典型的线性矩阵方程，譬如，a)AXB = C，b)AIXB1+A2XB2+… + = C， C)AXB+C3：D=E等 (参见 [1])，已经受到矩阵论和数值分析学者的关注和重视，取得了一 系列成果．例如，对于矩阵方程 AXB+CYD=E，文 [2-3】通过广义奇异值分解 (GSVD)，给 出了其相容条件和一般解；Xu等 [J应用矩阵对标准相关分解 (CCD)得到了它的最小二乘 解；文 [5】则讨论了AXAT+BYBT=C的一般解和对称解；Liao[0J研究了其相应的最佳逼近 问题．2003年，S．Y．Shim和Y．Chen[j给出了有效算法来求得矩阵方程AXB +CYD =E 的最小二乘