所谓部分矩阵是给定部分元素而其余元素是自由变量的矩阵。矩阵的填充是通过选取部分矩阵自由变量的值,使填充后的矩阵具有给定的某一性质,例如具有给定的逆、特征值、范数、行列式等。矩阵填充在很多科学领域和工程领域都具有重要的应用价值,比如运用矩阵填充原理使缺失的数据恢复
[1 ] 和其在推荐系统中的应用
[2 ] 。
国内外的许多学者对于部分矩阵逆的填充问题进行了大量研究,并取得了重要的研究成果。1986年,Fiedler等
[3 ] 在任意数域上对填充问题
A 11 A 12 A 21 ? - 1 = ? ? ? B 22 进行了研究。1996年,Hua
[4 ] 又进一步研究了上述矩阵的对称填充解。2013年,Lin等
[5 ] 给出了填充问题
? A 12 A 21 A 22 - 1 = B 11 ? ? ? 的解。2016年,赵琳琳等
[6 ] 给出了两类填充问题
A 11 A 12 A 21 ? - 1 = ? B 12 ? ? 和
A 11 X X * A 22 - 1 = ? B 12 B 12 * ? 有解的充分必要条件及其解的表示。2019年,白晓丽等
[7 ] 研究了一类Hamilton矩阵逆的填充,给出了解的表示并举例进行了验证。当然,有关部分矩阵的其他一些填充问题也取得了许多研究成果,例如矩阵最大最小秩,舒尔补的填充,见文献[
8 ]和文献[
9 ]。
现讨论两类部分矩阵逆的填充问题,给出了它们存在解的充分必要条件和解的表示,并举例进行了验证。令所有的m × n 阶复矩阵用C m × n 表示,用符号N * ,r a n k N ,n u l l N ,I m N 和K e r N 分别表示矩阵N 的共轭转置,秩,零空间的维数,值域和零空间。符号I 表示具有对应阶数的单位矩阵。
1 研究的问题
主要研究了以下两个填充问题。
问题1 已知A ∈ C n × m ,B ∈ C n × q ,C ∈ C p × m ,D ∈ C m × p ,E ∈ C q × n ,m + q = n + p ,求矩阵X ∈ C p × q ,使
问题2 已知A = A * ∈ C m × m , B ∈ C m × p , C ∈ C m × p ,求矩阵X = X * ∈ C p × p ,使
引理1[10 ] 令
α ⊆ 1 , 2 , ⋯ , n ,
α 的补
1 , 2 , ⋯ , n \ α 记为
α c ,如果
α , β ⊆ 1 , 2 , ⋯ , n ,那么
A α | β 表示由
A 中位于
α 行,
β 列的元素构成的子矩阵。若
A 可逆,则
引理2 A 是Hermite矩阵,若A 可逆,则A - 1 也是Hermite矩阵。
证明 设A 是Hermite矩阵,由其定义可知A * = A 。若A 可逆,则有A * - 1 = A - 1 ,即A - 1 * = A - 1 ,所以A - 1 也是Hermite矩阵。
事实上,问题1可描述为:已知A ∈ C n × m ,B ∈ C n × q ,C ∈ C p × m ,D ∈ C m × p ,E ∈ C q × n ,m + q = n + p ,找矩阵X ,Y ,Z 使
即找矩阵X ,Y ,Z 使
由引理2,问题2可描述为:已知A = A * ∈ C m × m ,B ∈ C m × p ,C ∈ C m × p ,找Hermite矩阵X ,Y ,Z 使
即找Hermite矩阵X ,Y ,Z 使
A B B * X Y C C * Z = Y C C * Z A B B * X = I (7)
2 主要结果及证明
定理1 问题1有解当且仅当r a n k A C = m ,r a n k A B = n ,并且存在矩阵Z ∈ C q × p 满足:
(a) n u l l A = n u l l Z ;
(b) A D + B Z = 0 ;
(c) E A + Z C = 0 ;
(d) K e r Z ⊆ K e r I - C D ;
(e) I m I - E B ⊆ I m Z 。
证明 显然,存在矩阵X 使A B C X 非奇异必须也一定有r a n k A C = m ,r a n k A B = n 。
必要性:问题1有解即存在矩阵X ,Y ,Z 使(5)式成立,则:
A Y + B E = I (8)
A D + B Z = 0 (9)
C Y + X E = 0 (10)
C D + X Z = I (11)
Y A + D C = I (12)
Y B + D X = 0 (13)
E A + Z C = 0 (14)
E B + Z X = I (15)
条件(a)可由引理1得出,条件(b)、(c)、(d)和(e)分别可由式(9 )、(14 )、(11 )和式(15) 得到。
充分性:假设存在矩阵Z 满足条件(a)~ (e)。下须证,存在矩阵X ,Y 满足式(4) 。式(4) 在块对角等价下是不改变的,所以设P ,Q ,R ,S 可逆,有
P 0 0 R A B C X Q 0 0 S = P A Q P B S R C Q R X S ,
Q - 1 0 0 S - 1 Y D E Z P - 1 0 0 R - 1 = Q - 1 Y P - 1 Q - 1 D R - 1 S - 1 E P - 1 S - 1 Z R - 1 ,则式(4) 等价于:
P A Q P B S R C Q R X S - 1 = Q - 1 Y P - 1 Q - 1 D R - 1 S - 1 E P - 1 S - 1 Z R - 1
其中,若将A ,B ,C ,D ,E 和Z 分别用P A Q ,P B S ,R C Q ,Q - 1 D R - 1 ,S - 1 E P - 1 和S - 1 Z R - 1 代替,条件(a)~ (e)仍成立。所以,可以选取P ,Q ,R ,S 使P A Q ,S - 1 Z R - 1 具有简单形式。不失一般性,可假设矩阵A ,Z 有下面的形式:
由条件(a)得A 与Z 中(2,2)零块相等,A 中的单位矩阵与Z 中的单位矩阵不一定同阶,不影响计算情况下,没有做标识。并假设B ,C ,D ,E ,X ,Y 是根据A ,Z 进行的分块。根据A 和Z 特殊的形式以及条件(b),得到
I 0 0 0 D 11 D 12 D 21 D 22 + B 11 B 12 B 21 B 22 I 0 0 0 = 0 0 0 0 ,
可得出
根据A 和Z 特殊的形式以及条件(c),得到
E 11 E 12 E 21 E 22 I 0 0 0 + I 0 0 0 C 11 C 12 C 21 C 22 = 0 0 0 0 ,
可得出
此外,由条件(a)可得D 22 ,E 22 都是方阵。由条件(d)可知存在矩阵X 1 = X 1 11 X 1 12 X 1 21 X 1 22 ,使X 1 Z = I - C D 成立,即
X 1 11 0 X 1 21 0 = I 0 0 I - C 11 0 C 21 C 22 D 11 0 D 21 D 22 。
比较上式(2 ,2 )处的元素得到C 22 D 22 = I ,因此D 22 可逆。上式可化简得
X 1 11 = I - C 11 D 11 ,
X 1 21 = - C 21 D 11 + C 22 D 21 。
由条件(e)可知存在矩阵X 2 = X 2 11 X 2 12 X 2 21 X 2 22 ,使Z X 2 = I - E B 成立,即
X 2 11 X 2 12 0 0 = I 0 0 I - E 11 E 12 0 E 22 B 11 B 12 0 B 22 ,
比较上式(2 ,2 )处的元素得到E 22 B 22 = I ,因此B 22 可逆,上式可化简得
X 2 11 = I - E 11 B 11 = I - C 11 D 11 ,
X 2 12 = - E 11 B 12 + E 12 B 22 。
可见X 1 11 = X 2 11 ,现选取满足上述两条件的X
X = I - C 11 D 11 - E 11 B 12 + E 12 B 22 - C 21 D 11 + C 22 D 21 X 22 (16)
其中X 22 任意。
接下来根据X 的取值,来选取Y 值以符合(4)式。由必要条件A Y + B E = I 和Y A + D C = I ,即
Y 11 Y 12 0 0 + B 11 B 12 0 B 22 E 11 E 12 0 E 22 = I 0 0 I ,
Y 11 0 Y 21 0 + D 11 0 D 21 D 22 C 11 0 C 21 C 22 = I 0 0 I ,
得
Y 11 = I - B 11 E 11 = I - D 11 C 11 , Y 12 = - B 11 E 12 + B 12 E 22 , Y 21 = - D 21 C 11 + D 22 C 21 。
由式(13) Y B + D X = 0 ,即
Y 11 Y 12 Y 21 Y 22 B 11 B 12 0 B 22 + D 11 0 D 21 D 22 X 11 X 12 X 21 X 22
= 0 0 0 0 , 得
Y 22 = - Y 21 B 12 + D 21 X 12 + D 22 X 22 B 22 - 1
从而,可选
Y = I - B 11 E 11 - B 11 E 12 + B 12 E 22 - D 21 C 11 + D 22 C 21 - Y 21 B 12 + D 21 X 12 + D 22 X 22 B 22 - 1 (17)
容易验证由式(16 )、(17 )所表示的矩阵X ,Y 满足式(4) 。
定理2 问题2有解当且仅当存在矩阵Z = Z * ∈ C p × p 满足:
(Ⅰ) n u l l A = n u l l Z ;
(Ⅱ) A C + B Z = 0 ;
(Ⅲ) K e r Z ⊆ K e r I - B * C ;
(Ⅳ) I m I - C * B ⊆ I m Z 。
证明 必要性:问题2有解即存在Hermite矩阵X ,Y ,Z 使式(7) 成立,则:
A Y + B C * = I (18)
A C + B Z = 0 (19)
B * Y + X C * = 0 (20)
B * C + X Z = I (21)
Y A + C B * = I (22)
Y B + C X = 0 (23)
C * A + Z B * = 0 (24)
C * B + Z X = I (25)
条件(Ⅰ)可由引理1得出,条件(Ⅱ),(Ⅲ)和(Ⅳ)分别可由(19),(21)和(25)式得到。
充分性:假设存在Hermite矩阵Z 满足条件(Ⅰ)~(Ⅳ)。下须证,存在Hermite矩阵X ,Y 满足式(6) 。A ,Z 是Hermite矩阵,则A ,Z 可以进行酉对角化。式(6) 在酉块对角相似下等价不变,则存在P ,R 为酉矩阵,有
P 0 0 R A B B * X P * 0 0 R * = P A P * P B R * R B * P * R X R * ,
P 0 0 R Y C C * Z P * 0 0 R * = P Y P * P C R * R C * P * R Z R * ,
使式(6) 等价于:
P A P * P B R * R B * P * R X R * - 1 = P Y P * P C R * R C * P * R Z R * 。
其中,若将A ,B ,B * ,C 和C * 分别用P A P * ,P B R * ,R B * P * ,P C R * 和R C * P * 代替,条件(Ⅰ)~(Ⅳ)仍成立。所以,可选取P ,R 使P A P * ,R Z R * 具有简单形式。不失一般性,可假设矩阵A ,Z 有下面的形式:
其中Λ 、Δ 为实对角矩阵并且可逆,并假设B ,B * ,C ,C * ,X ,Y 是根据A ,Z 进行的分块。根据A 和Z 特殊的形式以及条件(Ⅱ),得到
Λ C 11 Λ C 12 0 0 + B 11 Δ 0 B 21 Δ 0 = 0 0 0 0 ,
可得出
Λ C 11 + B 11 Δ = 0 ,Λ C 12 = 0 ,B 21 Δ = 0 ,即Λ C 11 = - B 11 Δ ,C 12 = 0 ,B 21 = 0 。
此外,由条件(Ⅰ)可得C 22 和B 22 是方阵。由条件(Ⅲ)可知存在矩阵X 1 = X 1 11 X 1 12 X 1 21 X 1 22 ,使X 1 Z = I - B * C 成立,即
X 1 11 Δ 0 X 1 21 Δ 0 = I 0 0 I - B 11 * 0 B 12 * B 22 * C 11 0 C 21 C 22 。
比较上式(2 ,2 )处的元素得到B 22 * C 22 = I ,因此C 22 可逆,上式可化简得
X 1 11 = I - B 11 * C 11 Δ - 1 = Δ - 1 + B 11 * Λ - 1 B 11 ,
X 1 21 = - B 12 * C 11 + B 22 * C 21 Δ - 1 。
由条件(Ⅳ)可知存在矩阵X 2 = X 2 11 X 2 12 X 2 21 X 2 22 ,使Z X 2 = I - C * B 成立,即
Δ X 2 11 Δ X 2 12 0 0 = I 0 0 I - C 11 * C 21 * 0 C 22 * B 11 B 12 0 B 22 ,
可化简得
X 2 11 = Δ - 1 I - C 11 * B 11 = Δ - 1 + B 11 * Λ - 1 B 11 ,
X 2 12 = - Δ - 1 C 11 * B 12 + C 21 * B 22 。
可见X 1 11 = X 2 11 ,X 1 21 * = X 2 12 ,现选取满足上述两条件的X ,即
X = Δ - 1 + B 11 * Λ - 1 B 11 - Δ - 1 C 11 * B 12 + C 21 * B 22 - B 12 * C 11 + B 22 * C 21 Δ - 1 X 22 (26)
其中X 22 为任意Hermite矩阵,所以X 为Hermite矩阵。
接下来根据X 的取值,选取Y 值以符合式(6) ,由必要条件A Y + B C * = I 和Y A + C B * = I ,即
Λ Y 11 Λ Y 12 0 0 + B 11 B 12 0 B 22 C 11 * C 21 * 0 C 22 * = I 0 0 I ,
Y 11 Λ 0 Y 21 Λ 0 + C 11 0 C 21 C 22 B 11 * 0 B 12 * B 22 * = I 0 0 I ,
得
Y 11 = Λ - 1 I - B 11 C 11 * = Λ - 1 + C 11 Δ - 1 C 11 * ,
Y 12 = - Λ - 1 B 11 C 21 * + B 12 C 22 * ,
Y 21 = - C 21 B 11 * + C 22 B 12 * Λ - 1 。
由必要条件Y B + C X = 0 ,即
Y 11 Y 12 Y 21 Y 22 B 11 B 12 0 B 22 + C 11 0 C 21 C 22 X 11 X 12 X 21 X 22 = 0 0 0 0 ,得
Y 22 = - Y 21 B 12 + C 21 X 12 + C 22 X 22 B 22 - 1 = B 22 - 1 * B 12 * Λ - 1 B 12 B 22 - 1 + C 21 Δ - 1 C 21 * - B 22 - 1 * X 22 B 22 - 1 (27)
从而,可选
Y = Λ - 1 + C 11 Δ - 1 C 11 * - Λ - 1 B 11 C 21 * + B 12 C 22 * - C 21 B 11 * + C 22 B 12 * Λ - 1 Y 22 (28)
其中Y 22 满足式(27) ,容易验证Y 是Hermite矩阵。且由式(26 )、(28 )所表示的矩阵X ,Y 满足式(6) 。
3 算例
例1 给定部分矩阵
A B C X = 1 2 0 1 - 1 2 i 2 0 0 0 2 - 1 - 1 2 i 0 ? ? 2 2 3 ? ? ,
求矩阵X ,Y ,Z 使下式成立,
A B C X - 1 = Y D E Z = ? ? - 4 + 2 i 0 ? ? 3 3 2 4 + 2 i 8 ? ? 0 1 2 ? ? 。
根据定理1,当取
P = 1 0 0 2 , Q = 2 0 0 3 , R = 1 0 0 1 2 , S = 2 0 0 1 2
时,上述两个矩阵等价于
A ' B ' C ' X ' = 1 0 2 - i 1 0 0 0 2 - 2 - i 0 ? ? 2 1 ? ? 和
Y ' D ' E ' Z ' = ? ? - 2 + i 0 ? ? 1 1 2 + i 2 ? ? 0 1 2 ? ? ,
方程解不唯一。若取
根据式(16 )、(17 )可取
X ' = - 4 - 6 - i 3 - 2 i 1 - i , Y ' = - 4 - 9 2 + 2 i i 5 2 + 1 2 i
经计算
A B C X = 1 2 0 1 - 1 2 i 2 0 0 0 2 - 1 - 1 2 i 0 - 2 - 12 - 2 i 2 2 3 3 - 2 i 4 - 4 i ,
Y D E Z = - 8 - 18 + 8 i - 4 + 2 i 0 3 i 15 + 3 i 3 3 2 4 + 2 i 8 2 0 0 1 2 0 0 ,
经验证A B C X - 1 = Y D E Z 。
例2 给定部分矩阵
A B B * X = 2 0 - 4 + 2 i - 2 + 2 i 2 + i 0 0 0 0 i - 4 - 2 i 0 ? ? ? - 2 - 2 i 0 ? ? ? 2 - i - i ? ? ? ,
求Hermite矩阵X ,Y ,Z 使下式
A B B * X - 1 = Y C C * Z = ? ? 2 - i 1 - i 0 ? ? 1 + i 2 i 2 + i 1 - i ? ? ? 1 + i 2 ? ? ? 0 - i ? ? ? 成立。
根据定理2,方程解不唯一,现取Z = 1 0 0 0 1 0 0 0 0 ,X 22 = 2 ,根据式(26 )、(28 ),可取
X = 11 6 - 2 i - 4 - 5 i 6 + 2 i 5 - 1 - 5 i - 4 + 5 i - 1 + 5 i 2 ,
Y = 15 2 5 2 - 4 i 5 2 + 4 i 13 2 ,
经验证A B B * X - 1 = Y C C * Z ,
即
2 0 - 4 + 2 i - 2 + 2 i 2 + i 0 0 0 0 i - 4 - 2 i 0 11 6 - 2 i - 4 - 5 i - 2 - 2 i 0 6 + 2 i 5 - 1 - 5 i 2 - i - i - 4 + 5 i - 1 + 5 i 2 - 1
= 15 2 5 2 - 4 i 2 - i 1 - i 0 5 2 + 4 i 13 2 1 + i 2 i 2 + i 1 - i 1 0 0 1 + i 2 0 1 0 0 - i 0 0 0 。
4 结论
基于矩阵与其逆矩阵之间秩的相等关系,通过块对角等价对复数域上的一类部分矩阵进行了逆的填充,得到了填充解存在的充分必要条件及相关解的表达式。此外,在充分利用Hermite矩阵性质和结构特点的基础上,又结合矩阵的秩理论,通过酉块对角相似,研究了一类Hermite矩阵逆的填充,同样得到了填充解存在的充分必要条件及相关解的表达式。这些结果为进一步研究部分矩阵的填充问题提供了有效的理论价值,也在矩阵填充的应用领域具有相应的实用价值。