矩阵方程AXA^T＋BYB^T=C的亚半正定解.pdfVIP

第23卷第6期 大 学 数 学 Vo1．23。№．6 2007年 12月 C()I LEGE MATHEMATICS Dec．2007 矩阵方程AxAT+BYB丁--C的亚半正定解 潘秋华 (江苏科技大学 数理系，江苏镇江 212003) [摘 要]讨论了矩阵方程AXA +BⅧ 一c关于亚半正定矩阵x，y有解的充要条件，并在有解时给出 了解的通式． [关键词]矩阵方程；矩阵分解；亚半正定矩阵；通解 [中图分类号]O151．21 [文献标识码]C [文章编号]1672—1454(2007)06—0150—04 1 引言与引理 文中 ，0 ，s ”分别表示m×rt实矩阵，rt阶实正交矩阵和对称半正定矩阵的全体； 曼n≥X。”：{A∈ IxrAX≥O，VxE )表示n阶亚半正定矩阵集合；A≥O表示A是对称半正定矩阵， R(A)，N(A)分别表示矩阵A的值域和零空问；A ，A 分虽表示矩阵A的转置及Moore—Penrose广义 逆．1I·lI是Frobenius范数，J为单位矩阵．记s(A)一A+A ，S (A)表示s(A)的MP逆． 关于X，l，的矩阵方程 AXA +BYB 一C (1) 的研究已有很多工作，文[1][2]研究了矩阵方程(1)的对称解，[3]研究了(1)的对称正定解．本文考虑 矩阵方程(1)的亚半正定解，讨论了方程有解的充要条件，并在有解时给出了解的一般表达式． 先给出几个有用的引理． 引理1l_3 矩阵A∈ ”的充要条件是S(A)∈S墨 ”． 引理2嘲 A∈S n≥X。 的充要条件是对适当阶数的任意矩阵P，都有P APES ． 引理3_4] A：fA ll A． l≥o(其中A 一A A~2=A22)的充要条件为lA 2 A22 j A11A+A12一A12，A11≥O， A22--A~zA+A1 2≥O． (或相应地A A去A =A『2， A ≥O， A 一A A去A『2≥O)、 由引理3及亚半正定矩阵的定义容易证明以下引理4． fA11 A1 2 f 引理4 设A—l“ 1∈ ”，则A∈ ”的充要条件为 【Azl Azz J i)S(A11)S (A11)(A12+A )一A12+A ； ii)S(A11)≥O； iii)S(A22)一(A21+A )S (A11)(A12+A )≥O (或相应地，S(A22)S (A22)(A21+A )一A21+A ，A22≥O，S(A11)一(A12+A五)S (A22)(A21+A ) ≥O)． 引理5Es 设D∈ ，则矩阵方程x+x 一D有解的充分必要条件为D—D ，其通解为 X=ID+H ， [收稿日期]2oo6—03—25 第6期 潘秋华：矩阵方程 AXA +B 一C的亚半正定解 151 其中H是任意反对称矩阵． 2 主要结论 定理1 已知A∈ mX n， ∈ ，C∈ ，矩阵对(A，B)的广义奇异值分解为 A—GQ1 H1， B—GQ2H2，