劉桂敏,張美黎,呂洪斌
(北華大學數學與統計學院,吉林 吉林 132013)
張量的特征值問題有重要的應用背景,如在盲源分離[1]、磁共振成像[2-3]、分子構象[4]等方面都有重要應用.其中,非負張量的特征值和特征向量有許多研究結果[5-8].本文給出一個具有一般形式的非負張量譜半徑(最大特征值)的估計不等式.
如果ai1i2…im≥0,ij=1,2,…,n,j=1,2,…,m,則稱為非負張量.我們記所有m階n維非負張量的集合為.
一個m階n維張量=(δi1…im)稱為單位張量,如果
定義1[5-6]對于m階n維張量和一個向量x=(x1,x2,…,xn)T,xm-1是一個向量,其第i個分量為
(
如果有λ∈和一個非零向量x=(x1,x2,…,xn)T,使得
類似于非負矩陣理論,我們稱ρ()=sup{|λ|:λ∈σ()}是張量的譜半徑,其中σ()是張量的特征值集合.
文獻[7]將不可約矩陣的概念推廣到張量.
定義2[7]對于m階n維張量=(ai1i2…im),如果存在一個非空的真子集I?〈n〉,使得
ai1i2…im=0, ?i1∈I, ?i2,…,im?I,
關于非負張量的特征值和特征向量文獻[7]給出了如下結果:
定理1[7]設=(ai1i2…im)∈,則存在λ0≥0和一個非負向量x0≠0,使得
(1)
定理2[7]設=(ai1i2…im)∈不可約,則方程(1)中的(λ0,x0)滿足:
(ⅰ)λ0>0是的一個特征值.
(ⅱ)x0>0,即x0的所有分量是正的.
(ⅲ)如果λ是的一個具有非負特征向量的特征值,則λ=λ0.此外,在不計倍數的情況下非負特征向量是唯一的.
(ⅳ)如果λ是的特征值,則|λ|≤λ0.
由定理2(ⅱ)和(ⅳ)可知,非負張量的譜半徑即是其按模最大特征值.
定義
引理1[9]設=(ai1i2…im)∈,=(ci1i2…im)∈,0≤≤,則ρ()≤ρ().
由引理1我們有
引理2設=(ai1i2…im)∈不可約,則ρ()>ai…i,i∈〈n〉.
文獻[9]給出了非負張量譜半徑的上下界估計結果:
定理3[9]設=(ai1i2…im)∈,則
文獻[10]給出了非負張量譜半徑上下界估計的改進結果:
定理4[10]設=(ai1i2…im)∈,且n≥2,則
δ()≤ρ(),
其中
δ(),
δi, j(.
本文給出非負張量譜半徑上下界估計的一個具有一般形式的不等式,在特別情況下改進了定理3和定理4.
定理5設=(ai1i2…im)∈,則
κ()≤ρ(),
其中
κ(),
(2)
不失一般性,假設
xt1≥xt2≥…≥xtn-1≥xtn>0.
(ⅰ)首先證明
ρ().
由方程組(2)有
因此,對l∈〈n-1〉,有
(ρ(
r′t1(
即
(ρ()-at1…t1-r′t1(.
(3)
同理,由方程組(2)有
因此
(ρ(.
(4)
(ρ()-at1…t1-r′t1(
因此
ρ(at1…t1+at2…t2+r′t1()+
從而有
ρ().
(ⅱ)現在證明
ρ()=κ().
由方程組(2)有
和
與(ⅰ)的證明相似,對l∈〈n-1〉容易得到
ρ().
因此,有
ρ()=κ().
綜合(ⅰ)、(ⅱ)有
κ()≤ρ().
其中,
推論1設=(ai1i2…im)∈,且n≥2,則
Δ()≤ρ(),
其中
Δ(),
Δi, j(ai…i+aj…j+r′i(.
例1設其中
計算知ρ()=10.688 1.由定理3有
10≤ρ()≤13.
由定理4計算有
δ1,2=10.831 0,δ2,1=10.099 0,δ1,3=11.164 4,δ3,1=11.844 3,δ1,4=11.446 2,δ4,1=12.708 2,
δ2,3=10.352 3,δ3,2=11.831 0,δ2,4=10.520 8,δ4,2=12.765 0,δ3,4=12.211 1,δ4,3=12.855 7.
因此
δ(
即有
10.099 0=δ()≤ρ()=12.855 7.
由推論1計算有
Δ1,2=10.577 7,Δ2,1=10.348 5,Δ1,3=11.402 2,Δ3,1=11.639 4,Δ1,4=11.787 1,Δ4,1=12.389 9,
Δ2,3=10.681 1,Δ3,2=11.262 1,Δ2,4=11.000 0,Δ4,2=12.068 4,Δ3,4=12.345 9,Δ4,3=12.700 6.
因此
Δ()}=12.700 6,
即有
10.348 5=Δ()≤ρ()=12.700 6.
因此,定理5的推論1是定理4和定理3的一個很好的改進.