一类Sylvester矩阵方程的迭代解法-JournalofNortheasternUniversity.PDF

第３８卷第６期 东 北 大 学 学 报 ( 自 然 科 学 版 ) Ｖｏｌ.３８ꎬＮｏ.６ ２０１７年 ６ 月 Ｊｏｕｒｎａｌ ｏｆ Ｎｏｒｔｈｅａｓｔｅｒｎ Ｕｎｉｖｅｒｓｉｔｙ(Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ) Ｊｕｎ. ２０ １７ 　 ｄｏｉ:１０.３９６９/ｊ.ｉｓｓｎ.１００５－３０２６.２０１７.０６.０２９ 一类Ｓｙｌｖｅｓｔｅｒ矩阵方程的迭代解法 邵新慧ꎬ彭　 程 (东北大学 理学院ꎬ辽宁沈阳　 １１０８１９) 摘　 　 　 要: 针对Ｓｙｌｖｅｓｔｅｒ矩阵方程给出了一种基于梯度的迭代解法.通过引入一个松弛参数和应用层次 识别原理ꎬ构建了一种新型的迭代方法求解一类 Ｓｙｌｖｅｓｔｅｒ矩阵方程.收敛分析表明ꎬ在一定的假设条件下对 于任意初始值ꎬ迭代解都收敛到精确解.数值算例也表明了所给方法的有效性和优越性. 关　 键　 词: Ｓｙｌｖｅｓｔｅｒ矩阵方程ꎻ迭代解法ꎻ梯度ꎻ收敛ꎻ松弛参数 中图分类号: Ｏ ２４１ꎻＴＰ ３０１６　 　 　 文献标志码: Ａ　 　 　 文章编号: １００５－３０２６(２０１７)０６－０９０９－０４ Ｉｔｅｒａｔｉｖｅ Ｓｏｌｕｔｉｏｎｓ ｔｏ Ｓｙｌｖｅｓｔｅｒ Ｍａｔｒｉｘ Ｅｑｕａｔｉｏｎｓ ＳＨＡＯＸｉｎ￣ｈｕｉꎬＰＥＮＧ Ｃｈｅｎｇ (Ｓｃｈｏｏｌ ｏｆ ＳｃｉｅｎｃｅｓꎬＮｏｒｔｈｅａｓｔｅｒｎ ＵｎｉｖｅｒｓｉｔｙꎬＳｈｅｎｙａｎｇ １１０８１９ꎬＣｈｉｎａ. Ｃｏｒｒｅｓｐｏｎｄｉｎｇ ａｕｔｈｏｒ:ＳＨＡＯ Ｘｉｎ￣ｈｕｉꎬ Ｅ￣ｍａｉｌ:ｘｉｎｈｕｉ１００２＠１２６.ｃｏｍ) Ａｂｓｔｒａｃｔ: Ｔｈｅ ｇｒａｄｉｅｎｔ￣ｂａｓｅｄ ｉｔｅｒａｔｉｖｅ ｓｏｌｕｔｉｏｎｓ ｔｏ Ｓｙｌｖｅｓｔｅｒ ｍａｔｒｉｘ ｅｑｕａｔｉｏｎｓ ａｒｅ ｇｉｖｅｎ. Ｂｙ ｉｎｔｒｏｄｕｃｉｎｇ ａ ｒｅｌａｘａｔｉｏｎ ｐａｒａｍｅｔｅｒ ａｎｄ ａｐｐｌｙｉｎｇ ｔｈｅ ｈｉｅｒａｒｃｈｉｃａｌ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｐｒｉｎｃｉｐｌｅꎬ ａｎ ｉｔｅｒａｔｉｖｅ ａｌｇｏｒｉｔｈｍ ｉｓ ｃｏｎｓｔｒｕｃｔｅｄ ｔｏ ｓｏｌｖｅ Ｓｙｌｖｅｓｔｅｒ ｍａｔｒｉｘ ｅｑｕａｔｉｏｎｓ. Ｃｏｎｖｅｒｇｅｎｃｅ ａｎａｌｙｓｉｓ ｉｎｄｉｃａｔｅｓｔｈａｔ ｔｈｅ ｉｔｅｒａｔｉｖｅ ｓｏｌｕｔｉｏｎｓ ｃｏｎｖｅｒｇｅ ｗｉｔｈ ｔｈｅ ｅｘａｃｔ ｓｏｌｕｔｉｏｎｓ ｔｏ ａｎｙ ｉｎｉｔｉａｌ ｖａｌｕｅ ｕｎｄｅｒ ｃｅｒｔａｉｎ ａｓｓｕｍｐｔｉｏｎｓ. Ｎｕｍｅｒｉｃａｌ ｅｘａｍｐｌｅｓ ａｒｅ ｇｉｖｅｎ ｔｏ ｔｅｓｔｉｆｙ ｔｈｅ ｅｆｆｉｃｉｅｎｃｙ ｏｆ ｔｈｅ ｐｒｏｐｏｓｅｄ ｍｅｔｈｏｄ. Ｋｅｙ ｗｏｒｄｓ: Ｓｙｌｖｅｓｔｅｒ ｍａｔｒｉｘ ｅｑｕａｔｉｏｎꎻ ｉｔｅｒａｔｉｖｅ ｓｏｌｕｔｉｏｎꎻ ｇｒａｄｉｅｎｔꎻ ｃｏｎｖｅｒｇｅｎｃｅꎻ ｒｅｌａｘａｔｉｏｎ ｐａｒａｍｅｔｅｒ 　 　 Ｓｙｌｖｅｓｔｅｒ 矩阵方程在很多领域中都有涉及ꎬ 可能得到精确解ꎬ此时数值解的研究就非常必要. 例如矩形域上的椭圆边值问题离散后会得到一个 近些年ꎬ国内外很多学者对各种类型的Ｓｙｌｖｅｓｔｅｒ Ｓｙｌｖｅｓｔｅｒ形式的矩形方程ꎻ在常微分方程定性理 矩阵方程的数值解进行研究ꎬ得到了一些有效的 论研究及数值求解的隐式Ｒｕｎｇｅ－Ｋｗｕｔｔａ方法与 数值求解算法[１－１５]. 块方法中常会碰到这类方程的求解问题ꎻ在线性 本文从Ｓｙｌｖｅｓｔｅｒ 矩阵方程ＡＸＢ ＝Ｆ 的迭代 统计领域中也会有Ｓｙｌｖｅｓｔｅｒ矩阵方程ꎻ在控制理 解法入手ꎬ先给出了它的迭代解的收敛条件ꎬ然后 论及应用中ꎬ极点配置、观测器及构造 Ｓｙｌｖｅｓｔｅｒ 介绍了计算 Ｓｙｌｖｅｓｔｅｒ 矩阵方程ＡＸＢ ＋ＣＸＤ ＝Ｆ 函数中都涉及上述方程的求解.除此之外ꎬ特殊形