-
矩形张量奇异值问题是张量谱理论研究的主要课题之一,对其进行的研究主要集中在3个方面:一是对奇异值的性质进行研究[1-3];二是对某些特殊奇异值(如按模最大奇异值)进行估计或计算[4-8];三是对所有奇异值进行定位,即在复平面上给出所有奇异值的包含集[9-12].最近,张量奇异值定位问题引起了广泛关注并获得了一些初步结果.文献[10]利用对左、右特征向量的按模最大分量进行分类讨论的思想首次给出了张量奇异值包含集.随后,文献[11]利用对矩形张量指标集的划分给出了奇异值的一个S-型包含集.文献[12]利用图的弱连通性给出了奇异值的一个新包含集.本文继续考虑张量奇异值定位问题,拟综合利用文献[10-11]中的技巧和方法给出张量奇异值的更精确的S-型包含集.新包含集的优势是在不增加额外计算量的情形下,仅对文献[10-11]中的某些包含集取交集就可得到比文献[10-12]中包含集更精确的包含集.
全文HTML
-
用
$\mathbb{R}$ ($\mathbb{C}$ )表示实(复)数集,p,q,m,n为正整数,m,n≥2,且N={1,2,…,n}.记A=(ai1…ipj1…jq),若则称A为一个(p,q)阶m×n维实矩形张量,记作A∈
$\mathbb{R}$ [p,q;m,n].若存在λ∈$\mathbb{C}$ 和非零向量x=(x1,…,xm)T,y=(y1,…,yn)T满足方程组则称λ为A的奇异值,x,y为相应于λ的左、右特征向量,其中l=p+q,Axp-1yq和x[l-1]为m维向量,它们的第i个分量为
Axpyq-1和y[l-1]为n维向量,它们的第j个分量为
弹性张量是满足p=q=2且m=n=2,3的矩形张量,其在非线性弹性材料学中有着重要的应用[1].文献[10]就m=n这种情形对矩形张量的奇异值进行了定位,给出了如下包含集定理:
定理1[10] 设A∈
$\mathbb{R}$ [p,q;n,n],则其中σ(A)表示A的所有奇异值所成的集合,
为了减少计算量,文献[11]通过划分N为非空真子集S及其补集S给出了如下S-型奇异值包含集定理:
定理2[11] 设A∈
$\mathbb{R}$ [p,q;n,n],S为N的非空真子集,S为S在N中的补集,则
-
定理3 设A∈
$\mathbb{R}$ [p,q;n,n],S为N的非空真子集,S为S在N中的补集,则证 设λ∈σ(A),x=(x1,x2,…,xn)T和y=(y1,y2,…,yn)T分别为λ对应的左、右特征向量,
则|xs|和|xt|中至少有一个是正数,|yg|和|yh|中至少有一个是正数.下面分4种情形来证明.
情形1 假设Vmax=|xs|,则|xs|>0.任取j∈S,(1)式的第s个方程可写为
情形1.1 若|xj|≥|yj|,则对(3)式取模,并应用三角不等式,得
即
若(4)式中|xj|>0,则由(1)式的第j个方程
得
将(4)式和(5)式相乘,并消去|xs|l-1|xj|l-1>0,可得
若(4)式中|xj|=0,则由|xs|>0得
此时(6)式仍成立.
情形1.2 若|yj|>|xj|,类似地可得
若(7)式中|yj|>0,由(2)式的第j个方程
得
将(7)式和(8)式相乘,并消去|xs|l-1|yj|l-1>0,可得
若(7)式中|yj|=0,则由|xs|>0得
此时(9)式仍成立.
由(6)式和(9)式得
此时
$\lambda \in {{\overset{\wedge }{\mathop{\gamma }}\, }_{s, j}}$ (A).由j∈S的任意性得$\lambda \in \bigcap\limits_{j\in \bar{S}}{{{\overset{\wedge }{\mathop{\gamma }}\, }_{s, j}}}$ (A).由s∈S得$\lambda \in \bigcup\limits_{i\in S}{\bigcap\limits_{j\in \bar{S}}{{{\overset{\wedge }{\mathop{\gamma }}\, }_{i, j}}}}$ (A).类似于情形1的证明,可完成其它3种情形的证明:
情形2 假设Vmax=|yg|,此时
$\lambda \in \bigcup\limits_{i\in S}{\bigcap\limits_{j\in \bar{S}}{{{{\tilde{\gamma }}}_{i, j}}}}$ (A).情形3 假设Vmax=|xt|,此时
$\lambda \in \bigcup\limits_{i\in \bar{S}}{\bigcap\limits_{j\in S}{{{\overset{\wedge }{\mathop{\gamma }}\, }_{i, j}}}}$ (A).情形4 假设Vmax=|yh|,此时
$\lambda \in \bigcup\limits_{i\in \bar{S}}{\bigcap\limits_{j\in S}{{{\overset{\wedge }{\mathop{\gamma }}\, }_{i, j}}}}$ (A).由定理1、定理2和定理3易得如下比较定理:
定理4 设A∈
$\mathbb{R}$ [p,q;n,n],S为N的非空真子集,S为S在N中的补集,则
-
例1 设A∈
$\mathbb{R}$ [2,2;3,3],其中a1111=a1122=a1133=a1222=a1231=a1233=a1313=a1322=a1323=a1132=a2113=a2122=a2123=a2213=a2323=a2331=a2333=a3113=a3121=a3123=a3131=a3211=a3223=a3313=a3321=a3323=a3332=1,a2111=a2133=a2222=a2312=a2332=a3112=a3122=a3132=a3133=a3212=a3213=a3221=a3222=a3231=a3232=a3233=a3311=a3322=a3331=2,a3333=3,a1131=9,a2131=10,a2121=14,其余元素均为0.下面对A的奇异值进行定位.取S={3},S={1,2},由定理1和定理2均得
由文献[12]中定理3.3得
由定理3得
显然γS(A)⊆T(A)⊆ΩS(A)⊆Ω(A).
例1表明:由定理3得到的张量奇异值包含集比由文献[10]中定理2.2、文献[11]中定理1和文献[12]中定理3.3得到的包含集精确.
例2 设A∈
$\mathbb{R}$ [2,2;2,2],其中a1111=a1112=a1222=a2112=a2121=a2221=1,其余元素均为0.经计算得所有奇异值为±3,±1.677 4±0.672 2i,±1.245 2±0.632 2i,±1.141 7±0.201 8i,±1.068 2±1.217 5i,±1,±0.859 9±0.507 2i,±0.822 6,±0.337 3±1.812 5i,±0.209 0±1.037 2i,±0.137 8±1.253 0i,0.取S={1},S={2},由定理3得γS(A)={z∈$\mathbb{C}$ :|z|≤3}.显然σ(A)⊆γS(A).例2表明:由定理3得到的张量奇异值包含集可以恰好包含所有的奇异值.