Sylvester-Kac矩阵的奇异值所在区间估计

2024-01-29白明月秦建国

宁夏师范学院学报 2024年1期

白明月,秦建国

(郑州商学院通识教育中心,河南巩义 451200)

设P是如下形式的(n+1)×(n+1)的三对角矩阵

(1)

其中n是一个非负整数.这个矩阵在统计力学和量子物理学中应用比较多,比如,布朗运动、在超大立体上随机漫步和氢原子的斯塔克效应等许多问题中.SYLVESTER JJ[1]猜想出计算P的特征值的公式λi=n-2i,i=0,1,2,…,n.MUIR T[2]第一个给出了此公式的证明.矩阵P的对称形式出现在SCHRODINGER E[3]的论文中,这些结果包含了波动力学的基础.后来,KAC M[4]应用生成函数的办法,证明Sylvester公式,还给出了相应特征向量的分量的多项式表示.VINCZE I[5]和近些年的EDELMAN A[6]通过把P表示成适当定义的函数空间上的微分算子的矩阵,给出了Sylvester公式简洁的证明.长期以来,关于矩阵P的谱分解的研究可谓如火如荼,但关于矩阵P的奇异值分解的研究却并不充分.BOROS T[7]和CARLOS M[8]分别给出了在n=2m,即P的阶数是奇数时,P至少有m+1个奇异值的平方是正整数这一猜想的证明.近年,人们对于矩阵奇异值的研究一直没有间断[9].本文给出矩阵P的其余m个奇异值存在的区间估计.

1 基本结果

(2)

证明先考虑m为奇数的情形.

下面利用置换矩阵将P化为四分块矩阵.

令Π=[e2,e4,…,e2m,e1,e3,e5,…,e2m+1],这里

ei=(0,…,0,1,0,…,0),

第i个位置元素是1,其余元素全为零,长度为2m+1的第i个单位向量,i=1,2,…,2m+1.不难验证

(3)

(3)式中Q∈R(m+1)×m和S∈Rm×(m+1)分别是如下形式的矩阵

又因为P、Q、S均为实矩阵,有

(4)

所以,矩阵QQ*=QQT和SS*=SST的所有特征值的算术平方根构成矩阵P的奇异值集合.而矩阵QQT的特征值或Q的奇异值已由BOROS T[7]给出.现在只需要讨论矩阵SST的特征值问题.由于目前尚未找到计算P的这些奇异值的公式,下面给出其存在的区间估计,为此计算SST,

下面利用关于矩阵特征值的Gerschgorin圆盘定理[10],估计SST的特征值所属范围.由于SST是一个m阶的三对角矩阵,应用Gerschgorin圆盘定理,加之SST的非零元素仅位于主对角线和紧邻该主对角线的上和下2条对角线上,只需考虑SST的这3条对角线上的元素.由于SST的元素又满足条件

(SST)i,j=(SST)m-i+1,m-j+1,i,j=1,2,…,m,

(5)

Dk={z||z-[k2+(2m-k)2]|≤k[2m-(k-2)]+(k+2)(2m-k)},k=3,5,…,m.

而第一行的元素不能用上式,把第一行所对应的圆盘单列,于是可知SST的特征值必然位于下述圆盘的并集之内

∪{z||z-[1+(2m-1)2]|≤3(2m-1)}.

(6)

但(SST)*=SST,表示矩阵SST是一个Hermitian矩阵,因此它的特征值是实数,即它的特征值全在复平面的实数轴上,这时(6)式中表示复变量的z可以换为表示实变量的x.进一步,矩阵SST还是一个半正定矩阵,从而是一个Hermitian半正定矩阵,它的特征值全部非负.于是,上述全部闭圆盘可以全部改为闭区间.

将(6)式中的圆盘进一步化简,发现上述m个圆盘(闭区间)具有逐个包含关系.

事实上,

即

先比较除第一行外所对应的区间的包含关系,由此取∀k1,k2∈{3,5,…,m},k1

由于这2个区间

D1=[4[(m-k1)2-m],4m(m+1)],D2=[4[(m-k2)2-m],4m(m+1)]

的右端点相同,而左端点之差是

4[(m-k1)2-m]-4[(m-k2)2-m]=4[2m-(k1+k2)](k2-k1),

由于k1,k2∈{3,5,…,m},k10,k2-k1>0,4[(m-k1)2-m]-4[(m-k2)2-m]>0.这意味着,D1⊂D2.

[0,4m(m+1)]⊂R.

(7)

再将此区间与第一个(圆盘)区间[4m2-10m+5,4m2+2m-1]比较,当m≥3时,