8.1 非负矩阵——不等式及其推广
设 B=[bij]∈Mn,r 和 A=[aij]∈Mn,r ,记
B⩾0 ,如果所有 bij⩾0
B>0 ,如果所有 bij>0
A⩾B ,如果 A−B⩾0
A>B ,如果 A−B>0
类似地可定义反向关系 ⩽ 和 < . 我们定义 ∣A∣=[∣aij∣] . 如果 A⩾0 : 则称 A 是非负矩阵, 又如果 A>0 , 则称 A 是正矩阵. 下列简单事实可直接由定义推出.
练习 设 A,B∈Mn,r . 证明:
(8.1.1) ∣A∣⩾0 对每个 A 成立; ∣A∣=0 当且仅当 A=0 .
(8.1.2) ∣aA∣=∣a∣∣A∣ 对所有 a∈C 成立.
(8.1.3) ∣A+B∣⩽∣A∣+∣B∣ .
(8.1.4) 如果 A⩾0 且 A=0 , 则当 n 或 r 大于 1 时 A>0 不一定成立.
(8.1.5) 如果 A⩾0,B⩾0 ,且 a,b⩾0 ,则 aA+bB⩾0 .
(8.1.6) 如果 A⩾B 且 C⩾D , 则 A+C⩾B+D .
(8.1.7) 如果 A⩾B 且 B⩾C , 则 A⩾C .
练习 现在假定 A,B,C,D∈Mn , 且 x,y∈Cn . 证明:
(8.1.8) ∣Ax∣⩽∣A∣∣x∣ .
(8.1.9) ∣AB∣⩽∣A∣∣B∣ .
(8.1.10) ∣Am∣⩽∣A∣m 对所有 m=1,2,⋯ 成立.
(8.1.11) 如果 0⩽A⩽B 且 0⩽C⩽D , 则 0⩽AC⩽BD .
(8.1.12) 如果 0⩽A⩽B ,则 0⩽Am⩽Bm 对所有 m=1,2,⋯ 成立。
(8.1.13) 如果 A⩾0 ,则 Am⩾0 ;如果 A>0 ,则 Am>0 对所有 m=1,2,⋯ 成立。
(8.1.14) 如果 A>0,x⩾0 ,且 x=0 ,则 Ax>0 .
(8.1.15) 如果 A⩾0 , x>0 且 Ax=0 , 则 A=0 .
(8.1.16) 如果 ∣A∣⩽∣B∣ , 则 ∥A∥2⩽∥B∥2 .
(8.1.17) ∥A∥2=∥∣A∣∥2 .
显然,最后两个论断对任意绝对范数都成立,Frobenius 范数 (l∘ 范数)只是其中一个例子.这些简单关系首先可应用于关于谱半径的不等式.
8.1.18 定理 设 A,B∈Mn . 如果 ∣A∣⩽B , 则 ρ(A)⩽ρ(∣A∣)⩽ρ(B) .
证明:对每个 m=1,2,… ,根据(8.1.10)和(8.1.12),有 ∣Am∣⩽∣A∣m⩽Bm ,因此,根据(8.1.16)和(8.1.17),对所有 m−1,2,… ,不等式
∥Am∥2⩽∥∣A∣m∥2⩽∥Bm∥2 和 ∥Am∥2mm⩽∥∣A∣m∥2mm⩽∥Bm∥2mm 成立,如果现在设 m→∞ 且运用(5.6.14),便导出 ρ(A)⩽ρ(∣A∣)⩽ρ(B).
8.1.19 推论 设 A,B∈Mn . 如果 0⩽A⩽B , 则 ρ(A)⩽ρ(B) .
8.1.20 推论 设 A∈Mn . 如果 A⩾0 , 义如果 A 是 A 的任一主子矩阵, 则 ρ(A)⩽ρ(A) . 特别是, max1,…,naij⩽ρ(A) .
证明:设 1⩽r⩽n ,且设 A 是 A 的 r×r 主子方阵。设 A 表示把 A 的各元排放在它们原来的位置(作为 A 的元)而把 0 元放在其余位置所构成的 n×n 矩阵。于是 ρ(A)=ρ(A) 且 0⩽A⩽A ,所以根据推论(8.1.9)可知, ρ(A)=ρ(A)⩽P(A) □
491
对于一个不一定是Hermite矩阵的谱半径来说,上述推论中的下界 aν⩽ρ(A) 是得到的第一个非平凡下界,不过 A 为非负的假定是本质的.
练习 试构造一个相似于 [0010] 的矩阵,且不含零元。它的谱半径是什么?它是非负矩阵吗?关于推论(8.1.20)的后一部分,这说明什么问题?
练习 证明,如果 A,B∈Mn ,且 0⩽A<B ,则 ρ(A)<ρ(B) 。提示:存在某个 α>1 使得 0⩽Λ⩽αΛ<B 。如果 ρ(A)=0 ,则从推论(8.1.19)可得出结论,如果 ρ(A)=0 ,则把推论(8.1.20)应用于 B 可得出结论。
因为我们马上就会有关于非负矩阵的谱半径的较好上界,定理(8.18)将用来求任意矩阵的谱半径的上界。
8.1.21 引理 设 A∈Mn ,且假定 A⩾0 。如果 A 的各个行和是常数,则 ρ(A)=∣A∣ 。如果 A 的各个列和是常数,则 ρ(A)=∣A∣ 。
证明:我们知道, ρ(A)⩽∥A∥ 对任意矩阵范数 ∥⋅∥ 成立,但是,如果各行和是常数,则 x=[1,⋯,1]t 是关于特征值 ∥A∥ 的特征向量。把同样的论证应用于 AΓ 便得到关于列和的论断。
8.1.22 定理 设 A∈Mn ,且假定 A⩾0 ,则
1…nminj=1∑naij⩽ρ(A)⩽1…nmaxj=1∑naij(8.1.23) 且
1⩽i⩽nmini=1∑nait⩽ρ(A)⩽1⩽j⩽nmaxj=1∑nait(8.1.24) 证明:设 α=min{i1,i2,i3,i4}∑j=1naij ,然后构造新矩阵 B ,使 A⩾B⩾0 ,且 ∑j=1nbij≡α 对所有 i=1,2,⋯ 成立。例如,如果 α=0 ,令 B=0 ,又如果 α>0 ,可以令 bij=αaij(∑j=1naij)21 。根据引理(8.1.21), ρ(B)=α ,又根据推论(8.1.19), ρ(B)⩽ρ(A) 。相应的上界易用类似的方法来证明。把行和界应用于 AT 便得到列和界。
练习 证明上述结果中关于上界的论断.
8.1.25 推论 设 A∈Mn . 如果 A⩾0 , 且对所有 i=1,2,…,n 有 ∑j=1naij>0 , 则 ρ(A)>0 . 特别是, 如果 A>0 , 或者如果 A 是不可约非负矩阵, 则 ρ(A)>0 .
练习 证明不可约矩阵不可能有零行或零列.
因为只要 S 可逆就有 ρ(S−1AS)=ρ(A) ,所以可以引进某些自由参数来推广上述定理。如果 S=diag(x1,⋯,xn) ,且所有 xi>0 ,则当 A⩾0 时 S−1AS⩾0 。把定理(8.1.22)应用于 S−1AS=(aijxjxi)−1 ,便得到下述更一般的结果。
8.1.26 定理 设 A∈Mn ,且假定 A⩾0 。则对任意正向量 x∈Cn ,有
1,r=nminxt1j=1∑naijxj⩽ρ(A)⩽1,r=nmaxxt1j=1∑naijxj(8.1.27) 和
1⩽j⩽nminxji=1∑nxiaij⩽ρ(A)⩽1⩽j⩽nmaxxji=1∑nxiaij.(8.1.28) 8.1.29 推论 设 A∈Mn , x∈Rn ,且假定 A⩾0 和 x>0 。如果 α,β⩾0 使 αx⩽Ax⩽βx ,则 α⩽ρ(A)⩽β 。如果 αx<Ax ,则 α<ρ(A) ;如果 Ax<βx ,则 ρ(A)<β 。
证明:如果 αx⩽Ax ,则 α⩽min1⩽i⩽nxi−1∑j=1naijxj 。由定理(8.1.26)得知 α⩽ρ(A) 。如果 αx<Ax ,则有某个 α′>α ,使得 α′x⩽Ax 。在这种情形 ρ(A)⩾α′>α ,所以 ρ(A)>α 。类似地可证明上界。
练习 完成推论(8.1.29)的证明.
8.1.30 推论 设 A∈Mn ,且假定 A 非负。如果 A 有正特征向量,则相应的特征值是 ρ(A) ;即如果 Ax=λx 且 x>0 和 A⩾0 ,则 λ=ρ(A) 。
证明:如果 x>0 且 Ax=λx ,则 λ⩾0 且 λx⩽Ax⩽λx 。另一方面,由推论(8.1.29)可知, λ⩽ρ(A)⩽λ □
8.1.31 推论 设 A∈Mn ,且假定 A 是非负矩阵。如果 A 有正特征向量,则
ρ(A)=x→0max1⩽i⩽nminxi1j=1∑naijxj=r→0min1⩽i⩽nmaxxi1j=1∑naijxj.(8.1.32) 练习 证明上述推论. 提示:在(8.1.27)中采用正特征向量.
8.1.33 推论 设 A∈Mn ,且假定 A 是非负矩阵。如果 A 有正特征向量 x ,则对所有 m=1
492
2,…和所有 i=1 ,…, n ,有
j=1∑naij(m)⩽[minxkmaxxk]ρ(A)m和[maxxkminxk]ρ(A)m⩽j=1∑naij(m),(8.1.34) 其中 Am≡[aij(m)] 。特别是,如果 ρ(A)>0 ,则对 m=1,2,…,[ρ(A)−1A]m 各元一致有界。
证明:如果 Ax=ρ(A)x ,则 A′′x=ρ(A)′′x 。如果 A≥0 ,则 A′′≥0 ,且有
ρ(A)mmax1⩽k⩽nxk⩾ρ(Am)x1=[Amx]1−∑j=1najimjx,⩾(min1:k:nxk)∑j=1naij(m) 对任意 i=1,2,…,n 成立,因为 x>0 ,所确定的上界可由除法得出.类似地有
ρ(A)mmin1⩽k⩽nxk⩽ρ(A)mxi=[Amx]i=∑j=1naj(m)xj⩽(max1⩽k⩽nxk)∑j=1naij(m) 对任意 i=1,⋯,n 成立. 因为 x>0 , 所确定的下界可由除法得出.
习题
如果 A⩾0 ,且对某个 k , Ak>0 ,证明 ρ(A)>0 .
给出一个 2×2 矩阵 A ,使得 A⩾0 , A 不是正的,且 A2>0 .
假定 A⩾0 ,且 A=0 。如果 A 有正特征向量,证明 ρ(A)>0 。
如果 ρ(A)<1 ,我们知道,当 m→∞ 时 Am→0 。如果 A⩾0 ,且 A 有正特征向量,试用推论(8.1.33)证明, ∣Am∣⩽ρ(A)mC(A) 对所有 m=1,2,⋯ 成立,其中 C(A) 是常数矩阵。试考察 A=[1011] ,说明 A 有正特征向量的假定不能省略。试从推论(5.6.13)来讨论这两个结果。
如果 A⩾0 有特征向量,证明 A 相似于其各行和为常数的非负矩阵。这个常数是什么?提示:利用定理(8.1.26)前面的注释。
我们将在(8.4)节证明,非负不可约矩阵必定有正特征向量。说明非负矩阵可以有正特征向量而且是可约的。
设 A=[aij]∈Mn 是非负矩阵且有正特征向量 x=[xi] . 试用(8.1.33)证明,
(a) n1ρ(A)mminxkmaxxk1−k−n1−k−n⩽max1⩽p⩽naip(m) 对每个 i=1,2,…,n 成立;
(b) limm→∞[∑p=1natp(m)]1−m=ρ(A) 对每个 i=1,2,…,n 成立.
对所有 m=1,2,… ,记 Am=[aij(m)]