两类部分矩阵逆的填充

张盼盼 ,  张澜

内蒙古工业大学学报(自然科学版) ›› 2024, Vol. 43 ›› Issue (01) : 12 -17.

PDF (564KB)
内蒙古工业大学学报(自然科学版) ›› 2024, Vol. 43 ›› Issue (01) : 12 -17. DOI: 10.13785/j.cnki.nmggydxxbzrkxb.2024.01.003
数理科学

两类部分矩阵逆的填充

作者信息 +

Completion of inverses for two kinds of partial matrices

Author information +
文章历史 +
PDF (577K)

摘要

根据秩理论和矩阵的可逆条件,研究了两类部分矩阵逆的填充问题,证明了两类部分矩阵逆的填充问题有解的充分必要条件,并给出了解的表示,最后分别举例验证结论。

Abstract

The filling problems for the inverses of two classes of partial matrices are studied on the basis of rank theory and invertibility conditions for matrices. The sufficient necessary conditions for the filling problems to have solutions are proved, and the representations of the solution are given, and examples are given to verify the conclusions, respectively.

关键词

矩阵填充 / 埃尔米特矩阵 / 值域 /

Key words

matrix completion / hermitian matrix / range / kernel

引用本文

引用格式 ▾
张盼盼,张澜. 两类部分矩阵逆的填充[J]. 内蒙古工业大学学报(自然科学版), 2024, 43(01): 12-17 DOI:10.13785/j.cnki.nmggydxxbzrkxb.2024.01.003

登录浏览全文

4963

注册一个新账户 忘记密码

所谓部分矩阵是给定部分元素而其余元素是自由变量的矩阵。矩阵的填充是通过选取部分矩阵自由变量的值,使填充后的矩阵具有给定的某一性质,例如具有给定的逆、特征值、范数、行列式等。矩阵填充在很多科学领域和工程领域都具有重要的应用价值,比如运用矩阵填充原理使缺失的数据恢复[1]和其在推荐系统中的应用[2]
国内外的许多学者对于部分矩阵逆的填充问题进行了大量研究,并取得了重要的研究成果。1986年,Fiedler等[3]在任意数域上对填充问题A11A12A21?-1=???B22进行了研究。1996年,Hua[4]又进一步研究了上述矩阵的对称填充解。2013年,Lin等[5]给出了填充问题?A12A21A22-1=B11???的解。2016年,赵琳琳等[6]给出了两类填充问题A11A12A21?-1=?B12??A11XX*A22-1=?B12B12*?有解的充分必要条件及其解的表示。2019年,白晓丽等[7]研究了一类Hamilton矩阵逆的填充,给出了解的表示并举例进行了验证。当然,有关部分矩阵的其他一些填充问题也取得了许多研究成果,例如矩阵最大最小秩,舒尔补的填充,见文献[8]和文献[9]。
现讨论两类部分矩阵逆的填充问题,给出了它们存在解的充分必要条件和解的表示,并举例进行了验证。令所有的m×n阶复矩阵用Cm×n表示,用符号N*rankNnullNImNKerN分别表示矩阵N的共轭转置,秩,零空间的维数,值域和零空间。符号I表示具有对应阶数的单位矩阵。

1 研究的问题

主要研究了以下两个填充问题。

问题1 已知ACn×mBCn×qCCp×mDCm×pECq×nm+q=n+p,求矩阵XCp×q,使

ABCX-1=?DE?

问题2 已知A=A*Cm×m, BCm×p, CCm×p,求矩阵X=X*Cp×p,使

ABB*X-1=?CC*?

引理1[10]α1, 2, , nα的补1, 2, , n\α记为αc,如果α, β1, 2, , n,那么Aα|β表示由A中位于α行,β列的元素构成的子矩阵。若A可逆,则

nullAα|β=nullA-1βc|αc

引理2A是Hermite矩阵,若A可逆,则A-1也是Hermite矩阵。

证明A是Hermite矩阵,由其定义可知A*=A。若A可逆,则有A*-1=A-1,即A-1*=A-1,所以A-1也是Hermite矩阵。

事实上,问题1可描述为:已知ACn×mBCn×qCCp×mDCm×pECq×nm+q=n+p,找矩阵XYZ使

ABCX-1=YDEZ

即找矩阵XYZ使

ABCXYDEZ=YDEZABCX=I

由引理2,问题2可描述为:已知A=A*Cm×mBCm×pCCm×p,找Hermite矩阵XYZ使

ABB*X-1=YCC*Z

即找Hermite矩阵XYZ使

ABB*XYCC*Z=YCC*ZABB*X=I

2 主要结果及证明

定理1 问题1有解当且仅当rankAC=mrankAB=n,并且存在矩阵ZCq×p满足:

(a)  nullA=nullZ

(b)  AD+BZ=0

(c)  EA+ZC=0

(d)  KerZKerI-CD

(e)  ImI-EBImZ

证明 显然,存在矩阵X使ABCX非奇异必须也一定有rankAC=mrankAB=n

必要性:问题1有解即存在矩阵XYZ使(5)式成立,则:

AY+BE=I
AD+BZ=0
CY+XE=0
CD+XZ=I
YA+DC=I
YB+DX=0
EA+ZC=0
EB+ZX=I

条件(a)可由引理1得出,条件(b)、(c)、(d)和(e)分别可由式(9)、(14)、(11)和式(15)得到。

充分性:假设存在矩阵Z满足条件(a)~(e)。下须证,存在矩阵XY满足式(4)式(4)在块对角等价下是不改变的,所以设PQRS可逆,有

P00RABCXQ00S=PAQPBSRCQRXS

Q-100S-1YDEZP-100R-1=Q-1YP-1Q-1DR-1S-1EP-1S-1ZR-1,则式(4)等价于:

PAQPBSRCQRXS-1=Q-1YP-1Q-1DR-1S-1EP-1S-1ZR-1

其中,若将ABCDEZ分别用PAQPBSRCQQ-1DR-1S-1EP-1S-1ZR-1代替,条件(a)~(e)仍成立。所以,可以选取PQRS使PAQS-1ZR-1具有简单形式。不失一般性,可假设矩阵AZ有下面的形式:

A=I000Z=I000

由条件(a)得AZ中(2,2)零块相等,A中的单位矩阵与Z中的单位矩阵不一定同阶,不影响计算情况下,没有做标识。并假设BCDEXY是根据AZ进行的分块。根据AZ特殊的形式以及条件(b),得到

I000D11D12D21D22+B11B12B21B22I000=0000

可得出

B11=-D11B21=0D12=0

根据AZ特殊的形式以及条件(c),得到

E11E12E21E22I000+I000C11C12C21C22=0000

可得出

C11=-E11C12=0E21=0

此外,由条件(a)可得D22E22都是方阵。由条件(d)可知存在矩阵X1=X111X112X121X122,使X1Z=I-CD成立,即

X1110X1210=I00I-C110C21C22D110D21D22

比较上式(22)处的元素得到C22D22=I,因此D22可逆。上式可化简得

X111=I-C11D11
X121=-C21D11+C22D21

由条件(e)可知存在矩阵X2=X211X212X221X222,使ZX2=I-EB成立,即

X211X21200=I00I-E11E120E22B11B120B22

比较上式(22)处的元素得到E22B22=I,因此B22可逆,上式可化简得

X211=I-E11B11=I-C11D11
X212=-E11B12+E12B22

可见X111=X211,现选取满足上述两条件的X

X=I-C11D11-E11B12+E12B22-C21D11+C22D21X22

其中X22任意。

接下来根据X的取值,来选取Y值以符合(4)式。由必要条件AY+BE=IYA+DC=I,即

Y11Y1200+B11B120B22E11E120E22=I00I
Y110Y210+D110D21D22C110C21C22=I00I

Y11=I-B11E11=I-D11C11,Y12=-B11E12+B12E22,Y21=-D21C11+D22C21

式(13)YB+DX=0,即

Y11Y12Y21Y22B11B120B22+D110D21D22X11X12X21X22

=0000,

Y22=-Y21B12+D21X12+D22X22B22-1

从而,可选

Y=I-B11E11-B11E12+B12E22-D21C11+D22C21-Y21B12+D21X12+D22X22B22-1

容易验证由式(16)、(17)所表示的矩阵XY满足式(4)

定理2 问题2有解当且仅当存在矩阵Z=Z*Cp×p满足:

(Ⅰ)  nullA=nullZ

(Ⅱ)  AC+BZ=0

(Ⅲ)  KerZKerI-B*C

(Ⅳ)  ImI-C*BImZ

证明 必要性:问题2有解即存在Hermite矩阵XYZ使式(7)成立,则:

AY+BC*=I
AC+BZ=0
B*Y+XC*=0
B*C+XZ=I
YA+CB*=I
YB+CX=0
C*A+ZB*=0
C*B+ZX=I

条件(Ⅰ)可由引理1得出,条件(Ⅱ),(Ⅲ)和(Ⅳ)分别可由(19),(21)和(25)式得到。

充分性:假设存在Hermite矩阵Z满足条件(Ⅰ)~(Ⅳ)。下须证,存在Hermite矩阵XY满足式(6)AZ是Hermite矩阵,则AZ可以进行酉对角化。式(6)在酉块对角相似下等价不变,则存在PR为酉矩阵,有

P00RABB*XP*00R*=PAP*PBR*RB*P*RXR*
P00RYCC*ZP*00R*=PYP*PCR*RC*P*RZR*

使式(6)等价于:

PAP*PBR*RB*P*RXR*-1=PYP*PCR*RC*P*RZR*

其中,若将ABB*CC*分别用PAP*PBR*RB*P*PCR*RC*P*代替,条件(Ⅰ)~(Ⅳ)仍成立。所以,可选取PR使PAP*RZR*具有简单形式。不失一般性,可假设矩阵AZ有下面的形式:

A=Λ000Z=Δ000

其中ΛΔ为实对角矩阵并且可逆,并假设BB*CC*XY是根据AZ进行的分块。根据AZ特殊的形式以及条件(Ⅱ),得到

ΛC11ΛC1200+B11Δ0B21Δ0=0000

可得出

ΛC11+B11Δ=0ΛC12=0B21Δ=0,即ΛC11=-B11ΔC12=0B21=0

此外,由条件(Ⅰ)可得C22B22是方阵。由条件(Ⅲ)可知存在矩阵X1=X111X112X121X122,使X1Z=I-B*C成立,即

X111Δ0X121Δ0=I00I-B11*0B12*B22*C110C21C22

比较上式(22)处的元素得到B22*C22=I,因此C22可逆,上式可化简得

X111=I-B11*C11Δ-1=Δ-1+B11*Λ-1B11
X121=-B12*C11+B22*C21Δ-1

由条件(Ⅳ)可知存在矩阵X2=X211X212X221X222,使ZX2=I-C*B成立,即

ΔX211ΔX21200=I00I-C11*C21*0C22*B11B120B22

可化简得

X211=Δ-1I-C11*B11=Δ-1+B11*Λ-1B11
X212=-Δ-1C11*B12+C21*B22

可见X111=X211X121*=X212,现选取满足上述两条件的X,即

X=Δ-1+B11*Λ-1B11-Δ-1C11*B12+C21*B22-B12*C11+B22*C21Δ-1X22

其中X22为任意Hermite矩阵,所以X为Hermite矩阵。

接下来根据X的取值,选取Y值以符合式(6),由必要条件AY+BC*=IYA+CB*=I,即

ΛY11ΛY1200+B11B120B22C11*C21*0C22*=I00I
Y11Λ0Y21Λ0+C110C21C22B11*0B12*B22*=I00I

Y11=Λ-1I-B11C11*=Λ-1+C11Δ-1C11*
Y12=-Λ-1B11C21*+B12C22*
Y21=-C21B11*+C22B12*Λ-1

由必要条件YB+CX=0,即

Y11Y12Y21Y22B11B120B22+C110C21C22X11X12X21X22=0000,得
Y22=-Y21B12+C21X12+C22X22B22-1=B22-1*B12*Λ-1B12B22-1+C21Δ-1C21*-B22-1*X22B22-1

从而,可选

Y=Λ-1+C11Δ-1C11*-Λ-1B11C21*+B12C22*-C21B11*+C22B12*Λ-1Y22

其中Y22满足式(27),容易验证Y是Hermite矩阵。且由式(26)、(28)所表示的矩阵XY满足式(6)

3 算例

例1 给定部分矩阵

ABCX=1201-12i20002-1-12i0??223??

求矩阵XYZ使下式成立,

ABCX-1=YDEZ=??-4+2i0??3324+2i8??012??

根据定理1,当取

P=1002Q=2003R=10012S=20012

时,上述两个矩阵等价于

A'B'C'X'=102-i10002-2-i0??21??

Y'D'E'Z'=??-2+i0??112+i2??012??

方程解不唯一。若取

Z'=1000X22'=1-i

根据式(16)、(17)可取

X'=-4-6-i3-2i1-iY'=-4-92+2ii52+12i

经计算

ABCX=1201-12i20002-1-12i0-2-12-2i2233-2i4-4i
YDEZ=-8-18+8i-4+2i03i15+3i3324+2i82001200

经验证ABCX-1=YDEZ

例2 给定部分矩阵

ABB*X=20-4+2i-2+2i2+i0000i-4-2i0???-2-2i0???2-i-i???

求Hermite矩阵XYZ使下式

ABB*X-1=YCC*Z=??2-i1-i0??1+i2i2+i1-i???1+i2???0-i???成立。

根据定理2,方程解不唯一,现取Z=100010000X22=2,根据式(26)、(28),可取

X=116-2i-4-5i6+2i5-1-5i-4+5i-1+5i2
Y=15252-4i52+4i132

经验证ABB*X-1=YCC*Z

20-4+2i-2+2i2+i0000i-4-2i0116-2i-4-5i-2-2i06+2i5-1-5i2-i-i-4+5i-1+5i2-1
=15252-4i2-i1-i052+4i1321+i2i2+i1-i1001+i20100-i000

4 结论

基于矩阵与其逆矩阵之间秩的相等关系,通过块对角等价对复数域上的一类部分矩阵进行了逆的填充,得到了填充解存在的充分必要条件及相关解的表达式。此外,在充分利用Hermite矩阵性质和结构特点的基础上,又结合矩阵的秩理论,通过酉块对角相似,研究了一类Hermite矩阵逆的填充,同样得到了填充解存在的充分必要条件及相关解的表达式。这些结果为进一步研究部分矩阵的填充问题提供了有效的理论价值,也在矩阵填充的应用领域具有相应的实用价值。

参考文献

[1]

方启林. 基于矩阵填充理论的阵元缺失数据的恢复[D]. 合肥:中国科学技术大学, 2020.

[2]

潘伟. 矩阵填充在推荐系统中的应用研究[D]. 赣州:江西理工大学, 2021.

[3]

FIEDLER M, MARKHAM T L. Completing a matrix when certain entries of its inverse are specified[J]. Linear Algebra and Its Applications, 1986, 74:225-237.

[4]

HUA D. Completing a symmetric 2×2 block matrix and its inverse[J]. Linear Algebra and Its Applications, 1996, 235:235-245.

[5]

LIN Y, WANG Q W. Completing a 2×2 block matrix of real quaternions with a partial specified inverse[J]. Journal of Applied Mathematics, 2013, 2013:271978.

[6]

赵琳琳, 张立华, 闫立梅. 部分矩阵逆的填充[J]. 华东师范大学学报(自然科学版), 2016(6):88-93.

[7]

白晓丽, 张澜. 一类Hamilton矩阵逆的填充[J]. 华东师范大学学报(自然科学版), 2019(2):56-62, 68.

[8]

WOERDEMAN H J. Minimal rank completions of partial banded matrices[J]. Linear and Multilinear Algebra, 1993, 36(1):59-68.

[9]

BRUALDI R A, HUANG Z J, ZHAN X Z. Singular, nonsingular and bounded rank completions of ACI-matrices[J]. Linear Algebra and Its Applications, 2010, 433(7):1452-1462.

[10]

BARRETT W W, LUNDQUIST M E, JOHNSON C R, et al. Completing a block diagonal matrix with a partially prescribed inverse[J]. Linear Algebra and Its Applications, 1995, 223/224:73-87.

基金资助

内蒙古自治区自然科学基金项目(2021LHMS01004)

内蒙古自治区自然科学基金项目(2023LHMS01016)

AI Summary AI Mindmap
PDF (564KB)

69

访问

0

被引

详细

导航
相关文章

AI思维导图

/