第五章 特征值的估计及对称矩阵的极性本章主要讨论数值代数中的三个特殊理论,
即特征值的估计广义特征值问题以及实对称矩阵(一般是Hermite矩阵)特征值的极小极大原理,其次也涉及到一些特征值和奇异值的扰动问题,最后简要地介绍矩阵直积的一些性质及其在线性矩阵方程求解方面的应用。这几方面的内容,在矩阵的理论研究与实际应用当中都有着相当重要的作用。
5.1特征值的估计特征值的界首先给出直接估计矩阵特征值模的上界的一些方法定理5.1 设A=(ars)(Rn×n,令
M=
表示A任一特征值,则的虚部Im()
满足不等式
(5.1.1)
推论 实对称矩阵的特征值都是实数。
事实上,当A这实对称矩阵时,M=0.
由定理5.1可得Im()=0,即为实数。
引理1 设B(Cn×n,列向量y(Cn满足||y||2=1,
则|yHBy|.
定理5.2 设A(Cn×n,则A的任一特征值满足
||||A||
 (5.1.3)
 (5.1.4)
推论,Hermite矩阵的特征值都是实数,
反Hermite矩阵的特征值为零或纯虚数。
事实上,当A为Hermite矩阵时,由式(5.1.4)
知Im()=0,即为实数;
当A为反Hermite矩阵时,由式(5.1.3)知
Re()=0,即为为零或纯虚数。
定义.5.1设

,
则称矩阵A按行(弱)对角占优。
定义5.2 设A(Cn×n。如果AT按行严格对角占优,则称A按列严格对角占优;
如果AT按行(弱)对角占优,
则称A按列(弱)对角占优。
对直接估计矩阵特征之乘积的模的界,
再给出以下两个方法。
定理5.3 设A=(ars)(Cn×n,
令Mr=|arr|+
如果A按行严格对角占优,则
(5.1.5)
且当ars=0(s>r)时,式(5.1.5)中等号成立。
定理5.4 (Hadamard’s inequality)
设A=(ars)(Cn×n,则有
(5.1.7)
且式(5.1.7)中等号成立的充分条件是某
as0=0或者(ar,as)=0(r(s),
这里a1,…,an表示A的n个列向量。
定理5.5 (Schur’s inequality)
设A=(ars)(Cn×n的特征值为(1,…,(n,则有
 (5.1.9)
定义5.3 设A=(aij)(Cn×n,称由不等式
|z-aii|≤Ri (5.1.10)
在复平面上确定的区域为矩阵A的第i个
Gerschgorin圆(盖尔圆)Gi的半径(i=1,…,n).
定理5.6 (Gerschgorin theorem1)
矩阵A=(aij)(Cn×n的一切特征值都在它的n个盖尔圆的并集之内。
定理5.7 (Gerschgorin theorem2)
由矩阵A的所有盖尔圆组成的连通部分中有且仅有A的k个特征值(盖尔圆相重时重复计数,特征值相同时也重复计数).
推论 若将式(5.1.10)中的Ri改作
 (5.1.13)
则定理5.6与定理5.7的结论仍然成立。
其中 (i为任意的正实数。
利用推论,有时能够得到更精确的特征值的包含区域。
例5,严格对角占优矩阵非奇异。
定理5.8 设不可约矩阵*A=(aij)n×n有一个特征值(在其n个盖尔圆|z-aii|≤Ri(I=1,…,n)
并集的边界上,则所有 n个圆周
|z-aii|=Ri(I=1,…,n) (5.1.16)
都通过点(。
谱与范数的关系:((A)(||A||
利用定理5.8,加强式(5.1.15) 的结果如下定理所述。
定理5.9如果A=(aij)n×n不可约,且存在i0使得

则有p(A)<||A||
定理5.10 (Ky Fan)设A=(aij)(Cn×n,B=(bij)(Rn×n,
如果bij(|aij|(I,j=1,…,n)则对A的任一特征值(,
必有I使
 (5.1.17)
引理2 设(和(是两个非负实数,0≤(≤1,
则有 (((1-(≤ ((+(1-()(
定理5.11 (Ostrowski theorem1)
设A=(aij)(Cn×n,0≤(≤1,(是A的任一特征值。
则存在i使得
 (5.1.19)
推论1,在定理5.11的条件下,存在i使得
|(-aii|≤(Ri(A)+(1-()Ri(AT) (5.1.24)
关于定理5.11,还有以下的推论推论2如果A奇异,取0≤(≤1,则存在i使得
(1) |aii|≤[Ri(A)]([Ri(AT)]1-(;
|aii|≤(Ri(A)+(1-()Ri(AT).
推论3 对于0≤(≤1,恒成立
(1) ((A) ≤max{|aii|+[Ri(A)]([Ri(AT)]1-(};
(2) ((A) ≤max{|aii|+(Ri(A)+(1-()Ri(AT)}.
推论4记取0≤(≤1,则有
(5.1.25)
推论5 (Farnell A B)

推论6 (BrauerA)

推论7 (Parker W V)

推论8 (Browne E T)

定理5.12 (Ostrowski theorem 2)
设A=(aij)(Cn×n(n(2),则对A的任一特征值(,
存在I,j(i(j),使(属于
(ij(A)={z|z(C,|z-aii| |z-ajj|Ri(A)Rj(A)}
推论 设A=(aij)(Cn×n(n(2),如果对于所有的i(j,恒有|aii| |ajj|>Ri(A)Rj(A),则detA(0.
广义特征值问题在振动理论中,常常会碰到形式如下的特征值问题,求数(使方程
Ax=(Bx (5.2.1)
有非零解x,这里A为n阶对称矩阵,
B为n阶实对称正定矩阵x为n维列向量。
当B=I时,式(5.2.1)就成为普通的特征值问题,
因此式(5.2.1)可以看作是对普通特征值问题的推广。
定义5.5 称形如式(5.2.1)的特征值问题为矩阵A相对于矩阵B的广义特征值问题,
简称为广义特征值问题;称满足式(5.2.1)要求的数(为矩阵A相对于矩阵B的特征值;
而与(相对应的非零解x称之为属于(的特征向量定义5.6 满足式(5.2.5) 向量系x1,…,xn称为按B标准正交化向量系;式(5.2.5)的第一式称作B正交(共轭)条件,
按B标准正交化向量x1,…,xn具有以下的性质。
性质1 xi(0(I=1,…,n)
性质2 x1,…,xn线性无关。
5.3对称矩阵特征值的极性在许多实际问题中,所产生的矩阵往往都具有对称性。如用等距的差分格式求解调和方程的第一类边值问题时所出现的矩阵,以及用有限元法求解某些结构问题时所产生的刚度矩阵,
一般都是对称的,特别是,实对称矩阵在理论研究与实际应用当中占有比较重要的地位。
因此,本节将着重讨论实对称矩阵的一些性质。
实对称矩阵的Rayleigh商的极性先引入如下定义定义5.7 设A是n阶实对称矩阵,x(Rn。称
 (5.3.1)
为矩阵A的Rayleigh商。
Rayleigh商(5.3.1)具有以下的特殊性。
性质1 R(x)是x的连续函数。
在0点不连续。
性质2 R(x)是x的零次齐次函数。
事实上,对任意的实数((0,有


性质3 x(L(x0)(x0(0)时,R(x)是一常数。
性质4 R(x)的最大值和最小值存在,且能够在单位球面S={x|x(Rn,||x||2=1}上达到。
补充性质:(Toeplitz定理)设A为任意复数矩阵,
定义
S={ z,z=xHAx,||x||=1},
则S为复平面上的闭凸集。
定理5.16 设A为实对称矩阵,则
 (5.3.3)
推论1 在||x||2=1上,p1和pn分别是R(x)的一个极小点和极大点,即有
R(p1)=(1,R(pn)=(n (5.3.5)
推论2 如果(1=…=(n(1kn),则在||x||2=1上,
R(x)的所有极小点为(1p1+…+(kpk
其中(i(R(i=1,…,k),且满足(21+…+(2k=1
定理5.17 设x(L(pr,pr+1,…,ps),1rsn,
则有
 (5.3.7)
定理5.18 (Courant-Fischer)设实对称矩阵A的特征值按式(5.3.2)的次序排列,则A的k个特征值
}(5.3.8)
其中Vk是Rn的任意一个k维子空间,1kn.
定理5.19 设实对称矩阵A和A+Q的特征值分别为
(1(2…(n和(1(2…(n,则有
|(I-(I|||Q||2 (I=1,…,n) (5.3.11)
定理5.20 (Hoffman-Wielandt)设实对称矩阵A,A+Q
和Q的特征值分别是(1(2…(n
和(1(2…(n,和(1(2…,( n
并定义向量u=((1,…,(n)T,v=((1,…,(n)T,w=((1,…,( n)T,
则 ||u-v||2||w||2.
定理5.21 (Lidskii-Wielandt)在定理5.20的条件下,
u落在形为v+pw向量集的凸包(即包含该向量的最小凸集)中,其中p取遍所有可能的排列矩阵。
定义5.8 设A,B为n阶实对称矩阵,且B正定,
x(R n.称为矩阵A相对于矩阵B的广义Rayleigh商。
定理5.22 非零向量x0是R(x)的驻点的充要条件是x0为Ax=(Bx的属于特征值(的特征向量,
这里A和B的意义同式(5.3.12).
推论 若是 Ax=(Bx的特征向量,则R()是与之对应的特征值。
定理5.23 设Vk为Rn的任意一个k维子空间,
则广义特征问题(5.2.1)的第k个特征值和第n-k+1
个特征值具有下列的极小极大性质
 (5.3.14)
 (5.3.15)
推论2 设Vn-k+1是Rn的任意一个n-k+1维子空间,
则定理5.23或推论1的结论可写成如下形式
 (5.3.18)
 (5.3.19)
对称正定矩阵的补充性质:
设A,B为对称(半)正定矩阵,矩阵C为A和B
的点积,即 cij=aij(bij,则C为对称(半)正定矩阵。
定理5.24 设A(Rrm×n的奇异值排列为
0=(1=…(n-r<(n-r+1…(n (5.3.21)
则A的第三k个奇异值和第n-k+1个奇异值具有下列的极性质
 (5.3.22)
 (5.3.23)
其中Vk是Rn的任一k维子空间。
定理5.25 设A(Rrm×n的奇异值排列同式(5.3.21),
(A+Q)(  的奇异排列为
0=(1=…(n-r<(n-r+1…(n (5.3.25)
则有
|(I-(I|||Q||2 (I=1,…,n) (5.3.26)
定理5.26 设A(Rrm×n和(A+Q)(的奇异值排列分别同式(5.3.21)和(5.3.25),Q(的排列值为
0=(1=… <…(n (5.3.29)
定义向量
u=((1,…,(n)T,v=((1,…,(n)T,w=((1,…,(n)T
则有 ||u-v||2||w||2
定义5.9 设A=(aij)(Cm×n,B=(bij)(Cp×q,则称如下的分块矩阵
 (5.4.1)
为A与B的直积(Kronecker积)。
定理5.27 设Am×m的特征值为(1,(2,…,(m,
Bn×n的特征值为(1,(2,…,(n,则f(A,B) 的全体特征值为f((1,(j)(I=1,2,…,m;j=1,2,…,n)
推论1设Am×m的特征值为(1,(2,…,(m,
Bn×n的特征值为(1,(2,…,(n,则的全体特征值为(i(j,(I=1,2,…,m;j=1,2,…,n).
推论2 设A(Cm×m,B=Cn×n,则有
det()=(detA)n(detB)m
推论3 设A(Cm×m,B=Cn×n,则有
tr()=(trA)(trB)
定理5.28 方程(5.4.6)有解的充要条件是
.这里,R(A)
表示矩阵A的列空间。
定理5.29 设Am×m的特征值为(1,(2,…,(m,
Bn×n的特征值为(1,(2,…,(n,则方程(5.4.5)有惟一解的充要条件是(I+(j(0(I=1,2,…,n).
推论1设Am×m的特征值为(1,(2,…,(m,Bn×n的特征值为(1,(2,…,(n,则齐次方程AX+XB=O
有非零解的充要条件是存在i0与j0,使=0。
推论2 设A是m阶矩阵,则齐次方程AX-XA=O一定有非零解。
定理5.30 设Am×m的特征值为(1,(2,…,(m,Bn×n的特征值为(1,(2,…,(n,则方程 有惟一解的充要条件是存在1+((I(j)+…+((I(j)l(0
(I=1,2,…,m;j=1,2,…,n).
齐次方程非零解的充要条件是存在i0与j0,使1+()+…+()l=0
引理3 设A(Cm×m,B=Cn×n,F(Cm×n如果A与B的特征值的实部都小于零,则积分存在。
定理5.31 设A(Cm×m,B=Cn×n,F(Cm×n且A与B的特征值之和不等于零,那么,如果积分
存在,则方程(5.4.5)的惟一解X=-。
推论1 设A(Cm×m与B=Cn×n,的特征值满足
Re((I)<0(I=1,…,m),Re((I)<0(j=1,…,n)
则方程(5.4.5)的惟一解为X=-
推论2 设A(Cm×m的特征值满足Re((I)<0(I=1,…,m),,则方程AHX+XA=-F的惟一解为
X=如果F(Cm×m是正定矩阵,则解矩阵X也是正定的。