(1.Department of Mathmatics,Minjiang University,Fuzhou 350108,China; 2.School of Mathmatical Sciences,Xiamen University,Xiamen 361005,China)
tensor" target="_blank">loop tensor; tensor" target="_blank"> inverse tensor; tensor
DOI: 10.6043/j.issn.0438-0479.201608007
备注
根据一种张量乘积、循环张量以及张量左右逆的基本定义和性质,给出可求逆的循环张量应该具有的结构特点以及相应逆的性质.通过一些例子给出了详细的求解循环张量逆的过程.
Basic definitions and characteristics of tensor products,loop tensors,and tensors of the inverse were introduced.The character of the structure and the inverse of the loop tensor,whose inverse could be solved,were provided.The solution procedure for the inverse of loop tensor was described along with some examples.
引言
1 预备知识
矩阵的逆在矩阵理论和计算中都有着重要的作用[1-3].尽管张量是矩阵的一种延伸,但对于张量的逆,至今没有一个非常理想的定义.在参考文献[4]中,给出了一种张量逆的定义,而在文献[2]中,给出了循环张量的定义.本文中,将依据这种张量逆和循环张量的定义,求解循环张量的逆,同时给出相应逆的性质(主要是结构特点).
本文中考虑的张量都是复数域C的,张量的表达形式和一些相关记号都来自文献[4].
定义1[4] 令A∈Cn1×n2×…×n2 和B∈Cn2×n3×…×nk+1 分别是阶数 m≥2和k≥1 的张量,那么它们的乘积AB是一个阶数为(m-1)(k-1)+1的张量C,它的元素按下式定义:
ciα1α2…αm-1=∑i1,i2…,im-1∈[n2]aii1i2…im-1bi1α1bi2α2…bim-1αm-1.(1)
这里 i∈[n1],α1,α2…,αm-1∈[n3]×[n4]…×[nk+1].
在上述定义中,如果 n1=n2=…=nk+1=n,那么 AB与文献[6-7]中所介绍的张量乘积一致.定义1 的张量乘积有如下的一些性质[4]:
1)(A1+A2)B=A1B+A2B,这里 A1,A2∈Cn1×n2×…×n22,B∈Cn2×n3×…×nk+1.
2)A(B1+B2)=AB1+AB2,这里A∈Cn1×n2,B1,B2∈Cn2×n3×…×nk+1.
3)AIn2=A,In2B=B,这里A∈Cn<1×n2×…×n2,B∈Cn2×n3×…×nk+1,In2 是维数为 n2 的单位矩阵.
4)A(BC)=(AB)C,这里 A∈Cn1×n2×…×n2,B∈Cn2×n3×…×n3,C∈Cn3×n4×…×nr.
如果一个 m 阶 n 维的张量I=(δi1i2…im)满足 δi1i2…im=1,当 i1=i2=…=im,而其他元素均为 0,则称这样的张量为单位张量.
定义2[4] 令A是m 阶 n 维的张量,B 是 k 阶 n 维的张量.如果 AB=I,那么 A 称作 B 的左逆 m 阶张量,B 称作 A 的右逆 k 阶张量.
一组向量a1∈Cn1,a2∈Cn2…,ak∈Cnk的Segre外积,记作a1a2…ak,其实是一个 k 阶张量 A∈Cn1×n2…×nk,其中的元素满足 ai1…ik=(a1)i1(a2)i2…(ak)ik.一个张量 A∈Cn1×n2…×nk 称为是秩 1 的,如果存在非零向量ai∈Cni(i=1,2,…,k),使得 A=a1a2…ak 成立.
在文献[5]中,Chen等给出了如下的循环张量定义.
定义3[5] 令A=(aj1j2…jm)是 m 阶 n 维张量.如果对任意 jl,kl∈[n],kl=jlmod(n)+1,l∈[m],满足
aj1j2…jm=ak1k2…km,(2)
那么称A是 m 阶 n 维循环张量.
下面将m阶n维张量的集合记作Tm,n,m阶n维循环张量的集合记作Cm,n.
对于循环矩阵而言,第一个行向量可以生成得到其他的行向量.同样的,循环张量也有类似的特点.本文中先介绍一个概念——行张量,对于张量 A=(aj1j2…jm)∈Tm,n,令Ak=(a(k)j1j2…jm-1)∈Tm-1,n 中的元素满足 a(k)j1j2…jm-1 ≡akj1j2…jm-1,称 Ak 为张量 A 的第 k 个行张量.若A是循环张量,则不难发现行张量 Ak,k=2,3,…,n 可以由第一个行张量 A1=(αj1j2…jm-1)得到,这里 αj1j2…jm-1≡a(1)j1/sub>j2…jm-1.本文中称A1 为循环张量 A的根张量.
现在记P=(pjk)∈T2,n 为置换矩阵,且满足 pjj+1=1 当 j∈[n-1],pn1=1,其他元素为 0.即
P=(0 1 … 0 0
0 0 1 … 0
0 0
0 1
1 0 … 0 0).(3)
那么,根据循环张量的定义 3,可以得到以下的一些性质.
引理1[4] 假设 A∈Cm,n 以及矩阵P定义为式(3),那么对于任意的 k∈[n],有
Ak+1=PTAkP,(4)
这里
An+1=A1.
引理2[4] 假设A∈Tm,n 以及矩阵P定义为式(3),那么下面3种情况是等价的:
a)A∈Cm,n.
b)PTAP=A.
c)对任意的循环矩阵C∈C2,n,有CTAC∈Cm,n.
2 循环张量的逆
2.1 循环矩阵的左右逆张量如果一个循环矩阵是非奇异的,那么它的逆矩阵也同样是循环矩阵.由定义 2 可知,循环矩阵不仅有阶数为 2 的左右逆张量,还拥有阶数更高的左右逆张量.它们的结构如何?是否也是循环张量?接下来本文中会通过理论证明和实验数据来详细说明.
定理1(循环矩阵的右逆张量)设A∈C2,n 是非奇异的,那么存在唯一的循环张量B∈Tm,n,使得AB=I∈Tm,n.
证明 存在唯一性:
通过张量乘积的性质,可得:
AB=I,
A-1AB=A-1I,
IB=A-1I,
B=A-1I.(5)
设A-1=(aij),由式(1),可以计算出B中的元素:
Bii1i2…im-1=∑nj=1aijiji1i2…im-1=
{ail, i1=i2…=im-1=l∈[n],
0, 其他.(6)
循环结构:
PTBP=PTA-1IP=
(PTA-1P)(PTIP)=
A-1I=B.(7)
这里矩阵 P 的定义为式(3).
从式(5)可见右逆张量中所有的行张量都是对角的.
例1 设A∈C2,3,它的根张量是(1,2,3).求张量B∈T3,3,使得
AB=I∈T3,3.
求解:
A=(1 2 3
3 1 2
2 3 1)
A-1=(-5/18 7/18 1/18
1/18 -5/18 7/18
7/18 1/18 -5/18).
根据式(5),容易计算出张量 B:
B1=(-5/18 0 0
0 7/18 0
0 0 1/18),
B2=(1/18 0 0
0 -5/18 0
0 0 7/18),
B3=(7/18 0 0
0 1/18 0
0 0 -5/18).(8)
容易验证,张量 B 是循环张量,且每个行张量都是对角的.
定理2(循环矩阵的左逆张量)令A∈C2,n 是非奇异的,那么存在唯一的循环张量 B∈Tm,n,使得BA=I∈Tm,n.
证明 存在唯一性:
通过张量乘积的性质,可得:
BA=I,
BAA-1=IA-1,
BI=IA-1,
B=IA-1.(9)
循环结构:
PTBP=PTIA-1P=
(PTIP)(PTA-1P)=
IA-1=B.(10)
这里矩阵P定义为式(3).
因此,只需要得到循环张量B的根张量B1 即可,设A-1=(aij).
B1=A-TI1A-1=a1…a1_}m-1,(11)
这里
I1=(δi1i2…im-1)∈Tm-1,n
满足
δi1i2…im-1={1, i1=i2…=im-1=1,
0, 其他.
而
a1=(a11,a12,…,a1n)T.
由上面的证明可见,左逆张量所有的行张量都是秩 1 的.
例2 设A∈C2,3,它的根张量是(1,2,3),求张量 B∈T3,3,使得BA=I∈T3,3.
求解:
B1=a1a1=(-5/18,7/18,1/18)T
(-5/18,7/18,1/18)T=
(25/324 -35/324 -5/324
-35/324 49/324 7/324
-5/324 7/324 1/324).(12)
同样可以得到:
B2=(1/324 -5/324 7/324
-5/324 25/324 -35/324
7/324 -35/324 49/324),
B3=(49/324 7/324 -35/324
7/324 1/324 -5/324
-35/324 -5/324 25/324).(13)
容易验证,张量B是循环张量,且每个行张量都是秩1的.
2.2 循环张量的左右逆矩阵对于循环矩阵的左右逆张量,已经清楚知道了其存在性、唯一性以及结构特点.现在,考虑这些结论是否适用于阶数大于等于 3 的循环张量.本文中首先讨论下循环张量的左右逆矩阵.
定理3(循环张量的右逆矩阵)令A∈Cm,n,当它的根张量 A1=a…a_}m-1 且以a为根向量的循环矩阵非奇异时,存在 B∈T2,n,使得 AB=I∈Tm,n.
证明 因为AB=I,可得到
bTj Albk={1, j=l=k∈[n],
0, 其他.(14)
这里bj 是矩阵B的第 j 列向量,Al 是张量 A 的第 l 个行张量.
进而,令
C=(A1b1,Ab2…,Anbn)∈T2,n,
可得到
BTC=I.(15)
这里 I 是维数 n 的单位矩阵.
注意矩阵 B 是非奇异的.由张量乘积的性质可得 A=IB-1.那么
A1=B-TI1B-1=b1…b1_}m-1,
这里bT1是B-1的第一行向量.
满足上述条件时,右逆矩阵虽然存在,但不一定是唯一的,也不一定有循环结构.
例3 设A∈C3,3,它的根张量是A1=(1,2,3)T(1,2,3)T∈T2,3.求矩阵 B∈T2,3,使得 AB=I∈T3,3.
求解:因为根张量
A1=(1,2,3)T(1,2,3)T∈T2,3,
所以
b1=(1,2,3)T或(-1,-2,-3)T.
同理,b2,b3 也都有 2 种情况.综合在一起,矩阵 B 总共有 8 种情况,可以求得 2 种:
B=(-5/18 7/18 1/18
1/18 -5/18 7/18
7/18 1/18 -5/18),
B'=(5/18 -7/18 -1/18
-1/18 5/18 -7/18
-7/18 -1/18 5/18),(16)
是右逆循环矩阵.而其他的6种:
B(1)=(5/18 7/18 1/18
-1/18 -5/18 7/18
-7/18 1/18 -5/18),
B(2)=(-5/18 -7/18 -1/18
1/18 5/18 -7/18
7/18 -1/18 5/18),
B(3)=(-5/18 -7/18 1/18
1/18 5/18 7/18
7/18 -1/18 -5/18),
B(4)=(5/18 7/18 -1/18
-1/18 -5/18 -7/18
-7/18 1/18 5/18),
B(5)=(-5/18 7/18 -1/18
1/18 -5/18 -7/18
7/18 1/18 5/18),
B(6)=(5/18 -7/18 1/18
-1/18 5/18 7/18
-7/18 -1/18 -5/18),(17)
仅是右逆矩阵,并不具有循环结构.
定理4 (循环张量的左逆矩阵)令A∈Cm,n,当它的根张量 A1 是对角的,并且这些对角元素组成的循环矩阵是非奇异时,存在唯一循环矩阵 B∈T2,n,使得BA=I∈Tm,n.
证明 存在唯一性:
由定义1可以得到
(BA)ii1i2…im-1=∑njbijaji1i2…im-1=
{1,i=i1=i2…=im-1=l∈[n];
0,其他.(18)
令
al=(a1l…l,a2l…l,…,anl…l)T,
可以得到方程
Bal=el,(19)
这里 el=(0,…,1,…,0)T 是只有第 l 个元素为 1,其他元素均为 0 的单位向量.
进一步,令A=(a1,a2,…,an),得到方程
BA=I,(20)
这里 I是维数 n 的单位矩阵.不难发现矩阵B是非奇异的.由张量乘积的性质可得A=B-1I.
循环结构:
A=B-1I,
PTAP=PTB-1IP,
A=(PTB-1P)(PTIP),
A=(PTB-1P)I,
B-1I=PTB-1PI,
(B-1-PTB-1P)I=0,
B-1=PTB-1P.(21)
最后一步的等价关系由文献[4]中的引理 2.1 保证.
例4 设A∈C3,3,它的根张量是A1=(1 0 0
0 2 0
0 0 3).求矩阵B∈T2,3,使得 BA=I∈T3,3.
求解:对角元素为(1,2,3),可得
A=(1 3 2
2 1 3
3 2 1).
根据BA=I,求得
B=(-5/18 1/18 7/18
7/18 -5/18 1/18
1/18 7/18 -5/18).
容易验证矩阵 B是左逆循环矩阵.
2.3 循环张量的左右逆张量现在,讨论下高阶循环张量的左右逆张量,研究它们的存在性、唯一性和结构性.
定理5 (循环张量的右逆张量)令A∈Cm,n,当它的根张量 A1=a…a_}m-1,且以a为根向量的循环矩阵非奇异时,存在张量B∈Tk,n,使得AB=I∈T(m-1)(k-1)+1,n,且其行张量都是对角的.
证明 根据定理3,存在非奇异矩阵B,使得
AB=I.(22)
同样,非奇异矩阵 B-1 存在右逆张量 B∈Tk,n,使得 B-1B=I,且其行张量都是对角的.
因此
ABB-1B=AB=I.(23)
根据定理3,循环张量的右逆矩阵不唯一且不一定具有循环结构,因此循环张量的右逆张量同样不唯一且不一定具有循环结构.
定理6 (循环张量的左逆张量)令A∈Cm,n,当它的根张量A1 是对角的,并且这些对角元素组成的循环矩阵是非奇异时,存在唯一循环张量B∈Tk,n,使得 BA=I∈T(k-1)(m-1)+1,n
证明 根据定理4,存在唯一非奇异循环矩阵B,使得
BA=I.(24)
同样,非奇异循环矩阵 B-1 存在唯一左逆循环张量B∈Tk,n,使得BB-1=I,且其行张量都是秩 1 的.
因此
BB-1BA=BA=I.(25)
3 结 论
通过上述讨论,可以发现以参考文献[4]中的张量积定义为基础,能求逆的循环张量是不够全面的,只有满足一定的结构时,才存在左右逆,且左右逆的唯一性和结构性还不一定能得到保证.后续工作希望通过不断的努力,能找到一个更加合适的张量积定义.
- [1] COX D,LITTLE J,O'SHEA D.Using algebraic geometry[M].New York:Springer,1998:1-572.
- [2] GELFAND I,KAPRANOV M,ZELEVINSKY A.Discriminants,resultants and multidimensional determinants[M].Boston:Birkhäuser,1994:1189-1214.
- [3] LIM L H.Tensors and hypermatrices[M]∥Handbook of Linear Algebra.Boca Raton,FL:CRC Press,2013:1-30.
- [4] BU C J,ZHANG X,ZHOU J,et al.The inverse,rank and product of tensors[J].Linear Algebra and Its Applications,2014,446:269-280.
- [5] CHEN Z M,QI L Q.Circulant tensors:native eigenvalues and symmetrization[EB/OL].[2014-06-05].https:∥arxiv.org/abs/1312.2752v6.
- [6] HU S L,HUANG Z H,LING C,et al.On determinants and eigenvalue theory of tensors[J].Journal of Symbolic Computation,2013,50:508-531.
- [7] SHAO J Y.A general product of tensors with applications[J].Linear Algebra Appl,2013,439:2350-2366.