6.2 正定矩阵
6.2.1-3
略。
6.2.4 设A对称正定,B是实矩阵。
1. 证明,对任意整数,Ak也正定。
2. 若存在正整数r,使得ArB=BAr,证明AB=BA.
1. 证明:
显然
A有谱分解
A=QΛQT,则
Ak=QΛkQT。而
Λ>0⟹Λk>0,所以
Ak也正定。
2. 证明:
因为
A对称正定,所以存在正交矩阵
Q,使
QAQT=Λ=λ1I⋱λnI,其中
λ1,λ2,⋯,λn均为正数且互不相等。于是
ArB=BAr⟹QTΛrQB=BQTΛrQ⟹ΛrQBQT=QBQTΛr
这说明
QBQT与
Λr可交换。所以
QBQT=B11B22⋱Bnn,显然其与
Λ也可交换。于是
AB=QTΛQB=QTΛQBQTQ=QTQBQTΛQ=BQTΛQ=BA.
6.2.5 对n阶实对称矩阵A,证明,当实数t充分大时,tIn+A正定。
证明:
xT(tIn+A)x=txTx+xTAx=t∥x∥+xTAx显然当
t充分大时,
t∥x∥+xTAx>0对任意
x=0都成立,所以
tIn+A正定。
6.2.6 对n阶实对称矩阵A,证明,存在正实数c,使得对任意x∈Rn,都有∣xTAx∣≤cxTx.
证明:
∣xTAx∣≤cxTx⟹xTx∣xTAx∣≤c。显然只需要
c≥∣λ∣max即可。
6.2.7 (Hadamard不等式) 给定n阶对称正定矩阵A,求证:
1.对任意y,det([AyTy0])≤0;
证明:
显然
A可逆,于是
AyTy0=AyT0−yTA−1y=−(yTA−1y)det(A),而
A−1显然也正定,即对任意
y,
yTA−1y≥0,等号成立当且仅当
y=0。而
det(A)>0,所以对任意
y,
det([AyTy0])≤0
2.记A=[aij],则det(A)≤anndet(An−1),其中An−1是A的n−1阶顺序主子阵;
证明:
显然
An−1可逆,于是
det(A)=An−1aTaann=AaT0ann−aTAn−1−1a=(ann−aTAn−1−1a)det(An−1)
显然
aTAn−1−1a≥0对任意
a成立,所以
det(A)≤anndet(An−1)。
3.det(A)≤a11a22⋯ann.
证明:
由上面结论立得。
利用上述结论证明:如果实矩阵T=[tij]可逆,那么det(T)2≤i=1∏n(t1i2+t2i2+⋯+tni2).
证明:
det(T)2=det(TTT),而矩阵
TTT显然正定,且其第
i行第
i列的对角元为
t1i2+t2i2+⋯+tni2。由上述结论得
det(T)2≤i=1∏n(t1i2+t2i2+⋯+tni2).
6.2.8 证明A=[i+j1]n×n正定。
证明:
法一:
注意到
i+j1=∫01ti+j−1dt
设
x=x1x2⋮xn,则
xTAx=1≤i≤n1≤j≤n∑i+jxixj=∫01t(1≤i≤n1≤j≤n∑xixjti+j−2)=∫01t(x1+x2t+⋯+xntn−1)2dt
显然对任意
x=0,
xTAx>0恒成立
⟹ A正定。
法二:
设
An=[i+j1]n×n,于是
det(An)=1+112+11⋮n+111+212+21⋮n+21⋯⋯⋯1+n12+n1⋮n+n1=(1+1)(n+1)n−1(2+1)(n+1)n−2⋮n+11(1+2)(n+2)n−1(2+2)(n+2)n−2⋮n+21⋯⋯⋯(1+n)(n+n)n−1(2+n)(n+n)n−2⋮n+n1=1≤i≤n∏n+i11≤i<n∏(n−i)1+112+11⋮11+212+21⋮1⋯⋯⋯1+n12+n1⋮1=1≤i≤n∏n+i11≤i<n∏(n−i)(1+1)(1+n)n−1(2+1)(2+n)n−1⋮0(1+2)(1+n)n−2(2+2)(2+n)n−2⋮0⋯⋯⋯1+n12+n1⋮1=1≤i≤n∏(n+i)1≤i<n∏(n−i)(1+1)(1+n)n−1(2+1)(2+n)n−1⋮(n−1+1)(n−1+n)n−1(1+2)(1+n)n−2(2+2)(2+n)n−2⋮(n−1+2)(n−1+n)n−2⋯⋯⋯(1+n−1)(1+n)n−(n−1)(2+n−1)(2+n)n−(n−1)⋮(n−1+n−1)(n−1+n)n−(n−1)=1≤i≤n∏(n+i)1≤i<n∏(i+n)1≤i<n∏(n−i)21+112+11⋮n−1+111+212+21⋮n−1+21⋯⋯⋯1+n−112+n−11⋮n−1+n−11=(n+n)1≤i<n∏(n+i)1≤i<n∏(i+n)1≤i<n∏(n−i)2An−1
而
A1=21,于是
det(An)=1≤i<j≤n∏(i−j)21≤i≤n1≤j≤n∏(i+j)1,显然任意顺序主子式都为正,所以
A正定。
6.2.9 证明Hilbert矩阵Hn=[i+j−11]n×n正定。
证明:
与6.2.8类似,也有两种方法。
法一:
注意到
i+j−11=∫01ti+j−2dt,
设
x=x1x2⋮xn,则
xTAx=1≤i≤n1≤j≤n∑i+j−1xixj=∫01(1≤i≤n1≤j≤n∑xixjti+j−2)=∫01(x1+x2t+⋯+xntn−1)2dt
显然对任意
x=0,
xTAx>0恒成立
⟹ A正定。
法二:
法二与6.2.8中法二几乎完全一致,故从略。
6.2.10 设A,B是实对称矩阵,如果A−B正定,则记作A≻B。求证:
1.A≻B,B≻C可以推出 A≻C.
2.A≻B和B≻A不可能同时成立。
3.对任意实对称矩阵A,都存在实数k1,k2使得k1In≻A≻k2In.
1. 证明:
A−B正定
⟹ 对任意
x=0,xT(A−B)x。
B−C正定
⟹ 对任意
x=0,xT(B−C)x。
两式相加立得
xT(A−C)x,所以
A−C正定,即
A≻C.
2. 证明:
假设同时成立,则
{xT(A−B)x>0xT(B−A)x>0,两式相加得
xT(A−B+B−A)x=0>0,矛盾。所以
A≻B和
B≻A不可能同时成立。
3. 证明:
先证存在实数
k1使得
k1In≻A,即
(k1In−A)正定。
对任意非零向量
x,
xT(k1In−A)x=k1xTx−xTAx≥k1xTx−λmaxxTx。显然当
k1足够大时,上式大于0对任意非零向量
x均成立。于是存在实数
k1使得
k1In≻A。
同理可得存在足够小的实数
k2使得
A≻k2In。
6.2.11 举例说明,实对称矩阵A的所有顺序主子式都非负,但A并不半正定.
解:
[000−1]所有顺序主子式都为0,但不半正定。
6.2.12 证明命题 6.2.4. 对实对称矩阵A,以下叙述等价: 1. A半正定。2. A的特征值都是非负数。3. 存在列满秩矩阵T,使得A=TTT。4. A存在LDLT分解,且D的对角元都是非负数。
证明:
采用轮转证法。
1⟹2 任取
A的一个特征向量
x,则
xTAx=λxTx≥0,于是
λ≥0.
2⟹3设
A=QΛQT,其中
Λ=diag(λ1,⋯,λr,0,⋯,0)。因为
Λ对角元非负,所以设
Λ=D2,则
A=(QD)(QD)T=[λ1q1⋯λrqr0⋯0]λ1q1T⋮λrqrT0T⋮0T,设
T=[λ1q1⋯λrqr],显然
T列满秩且
A=TTT.
3⟹4设
TT的QR分解为
QLT,其中
Q是列正交矩阵,
L是下三角矩阵。于是
A=TTT=LQTQLT=LLT,令
L=LD,其中
L是单位下三角阵,
D是对角阵则
A=LLT=LD2LT,令
D=D2立得。
4⟹1 xTAx=(LTx)TD(LTx)≥0.
6.2.13 证明,实对称矩阵半正定,当且仅当它的所有主子式都非负,其中主子式的定义见练习 4.3.22.
证明:
法一:
先证:实对称矩阵
A半正定
⟹ A所有主子式都非负。
显然,
A的某个
k阶主子阵可表示为
PTAP,P=[ei1ei2⋯eik],其中
i1<i2<⋯<ik,
ei是标准基。于是对任意
x=0,都有
xTPTAPx=(Px)TA(Px)≥0⟹所有主子阵都半正定
⟹A所有主子式都非负。
再证:实对称矩阵
A所有主子式都非负
⟹ A半正定。
用数学归纳法,按照阶数归纳。
阶数
n=1时显然成立。现假设阶数小于
n时命题成立,考虑
n阶实对称矩阵
A。分两种情况讨论:
(a)若
rank(A)<n。
设
x1,⋯,xn−r是
N(A)的一组基。逐一考虑标准基
e1,⋯,en,若
ei与
x1,⋯,xn−r线性无关,则将其加入基中。这样可以得到
Rn的一组基
(ei1,⋯,eir,x1,⋯,xn−r),其中
i1<i2<⋯<ik。设
X=[ei1⋯eirx1⋯xn−r],则
XTAX=[Cr×rOOO].
因为
C是
A的一个主子阵,且
A所有主子式都非负,所以
C所有主子式也都非负。根据归纳假设,
C半正定,于是
XTAX也半正定,而合同变换不改变正定性,所以
A也半正定。
这说明若阶数小于
n时成立,矩阵阶数为
n时也成立,于是
n≥1时命题均成立。
(b)若
rank(A)=n。
显然此时
det(A)>0,这说明
A只可能有0个或偶数个负数特征值。
根据归纳假设,
A的
n−1阶顺序主子阵
An−1半正定。假设
A有偶数个负数特征值,任意取其中两个负特征值
λi,λj,对应的特征向量是
x,y,注意
xTy=yTx=0.
设
u=αx+βy,存在一组
α,β使得
u的最后一个元素是为0,但
u=0。即
u=[u10],u1=0。于是
uTAu=u1TAn−1u1=(αx+βy)TA(αx+βy)=α2λixTx+αβλjxTy+αβλiyTx+β2λjyTy=α2λixTx+β2λjyTy<0
即存在
u1使得
u1TAn−1u1<0,这与归纳假设相矛盾。所以
A没有负数特征值
⟹ A半正定。
这说明若阶数小于
n时成立,矩阵阶数为
n时也成立,于是
n≥1时命题均成立。
法二:
先证:实对称矩阵
A半正定
⟹ A所有主子式都非负。
考虑矩阵
A+ϵI,其中
ϵ是任意小的正实数。于是对任意
x=0,有
xT(A+ϵI)x=xTAx+ϵxTx>0⟹A+ϵI正定
⟹A+ϵI的任意
k阶主子式
det(Ak+ϵI)>0。而
det(Ak+ϵI)是关于
ϵ的连续函数,所以
det(Ak)≥0,即
A所有主子式都非负。
再证:实对称矩阵
A所有主子式都非负
⟹ A半正定。
设
Ak是
A的任一
k阶主子阵,
det(Ak)≥0。
考虑矩阵
A+ϵI,其中
ϵ是任意小的正实数。其任一
k阶主子式为
det(Ak+ϵIk)。设
λ是
Ak的特征值,则
det(Ak−λIk)=0⟹det(Ak+ϵIk−(λ+ϵ)Ik)=0.这说明若
λ是
Ak的特征值,则
λ+ϵ是
Ak+ϵIk的特征值。因为
ϵ是任意小的正实数,所以若
λ=0,
λ+ϵ不改变符号,若
λ=0,
λ+ϵ>0,于是
det(Ak+ϵIk)>0⟹A+ϵI任意顺序主子阵为正
⟹A+ϵI正定。
所以对任意
x=0,有
xT(A+ϵI)x=xTAx+ϵxTx>0,而
xTAx+ϵxTx是关于
ϵ的连续函数,于是
xTAx≥0⟹A半正定。
6.2.14 设A是实对称矩阵,证明:
1.A半正定,当且仅当存在实对称矩阵B,使得A=B2.
2.若A半正定,则存在唯一的半正定实对称矩阵B,使得A=B2.
1.证明:
A是实对称矩阵且半正定
⟺A=QTΛQ,其中
Λ是对角元均为非负数的对角阵
⟺存在对角阵
D,使得
A=QTD2Q=QTDQQTDQ=B2.
2.证明:
A=QTΛQ,其中
Λ是对角元均为非负数的对角阵,
Λ=diag(λ1,λ2,⋯,λn),显然存在唯一的分解
Λ=D2使得
D也是对角元均为非负数的对角阵,即
D=diag(λ1,λ2,⋯,λn)。此时
A=QTD2Q=QTDQQTDQ=B2,且
B=QTDQ是半正定矩阵。
6.2.15 证明,若实对称矩阵对角占优,且对角元素全为正数,则该矩阵正定。
证明:
法一:
设
A对角占优,则
det(A)=0.
考虑矩阵
A+tIn,其中
t是充分大的正实数。显然
A+tIn也对角占优,且
xT(A+tIn)x=xTAx+txTx>0,于是
A+tIn正定。
设
Ak是
A的
k阶顺序主子阵,则对任意
k,都有
det(Ak+tIk)>0。又因为
det(Ak+tIk)是关于
t的连续函数,且对任意
t≥0,
det(Ak+tIk)=0,而
t充分大时,
det(Ak+tIk)>0。所以
det(Ak)>0对任意
k均成立
⟹A正定。
法二:
设
x=x1x2⋮xn,则
xTAx=1≤i≤n1≤j≤n∑aijxixj=1≤i<j≤n∑2aijxixj+1≤i≤n∑aiixi2(由对称性)>1≤i<j≤n∑2aijxixj+1≤i≤n∑xi2i=j∑∣aij∣(由对角占优)=1≤i<j≤n∑2aijxixj+1≤i<j≤n∑∣aij∣(xi2+xj2)(由对称性)≥1≤i<j≤n∑2aijxixj+1≤i<j≤n∑2∣aij∣xixj(由Cauchy不等式)≥0
即
xTAx>0对任意
x=0恒成立
⟹A正定。
6.2.16 证明或举出反例:
1.如果A对称正定,则A−1正定。
证明:
A对称正定则其特征值全是正数,而
A−1的特征值是
A特征值的倒数,显然也全为正数,所以
A−1正定。
2.如果A,B对称正定,则A+B正定。
显然正确。
3.如果A,B对称半正定,则A+B半正定。
显然正确。
4.如果A,B对称不定,则A+B不定。
不正确。
设
A=[−1111],其特征值是
±2,不定。
B=[1−1−1−2],其特征值是
2−1±13,不定。而
A+B=[000−1],其特征值是
0或
−1,半负定。
5.如果A列满秩,B对称正定,则ATBA正定。
证明:
A列满秩,所以当
x=0时,
Ax=0。而
B对称正定,即对任意
x=0,xTBx>0。所以
xTATBAx=(Ax)TB(Ax)>0⟹ATBA正定。
6.如果S=ATA且A有简化 QR 分解A=QR,则S=RTR是S的 Cholesky 分解。
证明:
S=ATA=(QR)T(QR)=RTQTQR=RTR,其中
R是上三角阵。
7.如果A,B对称正定,则AB的特征值都是正数。
设
AB的特征值是
λ,对应的特征向量是
x,则
ABx=λx。两边左乘
xTBT得
xTBTABx=λxTBTx⟹(Bx)TA(Bx)=λxTBTx,而
A,B都正定,
(Bx)TA(Bx)>0,xTBTx>0,所以
λ>0.
6.2.17 证明命题 6.2.7. 方阵的合同关系是等价关系。
证明:
反身性:
A=InTAIn,因此
A和
A合同。
对称性:若
A,B合同,则存在可逆矩阵
X,使得
XTAX=B,显然同时也有
X−TBX−1=A。所以
B和
A也合同。
传递性:若
A,B合同,
B,C合同,则存在可逆矩阵
X,Y,使得
XTAX=B,YTBY=C。即
(XY)TA(XY)=C,而
XY是可逆矩阵,因此
A,C也合同。
所以方阵的合同关系是等价关系。
6.2.18 考虑实对称矩阵S=[ABTBC].证明S正定当且仅当A及其 Schur 补都正定。
证明:
注意到
[ABTBC]=[IBTA−1I][AC−BTA−1B][IA−1BI]
于是
S正定
⟺S有
LDLT分解
⟺存在下三角阵
L1,L2使得
A=L1DL1T,C−BTA−1B=L2DL2T⟺ A及其Schur补都正定。
6.2.19 设A,B是n阶实对称矩阵,A正定。证明,存在可逆矩阵T,使得TTAT和TTBT同时是对角矩阵。
证明:
A是正定的实对称矩阵,所以存在正交矩阵
QA使得
QATAQA=ΛA。设
ΛA=D2,则
(QAD−1)TA(QAD−1)=In.
注意到
(QAD−1)TB(QAD−1)也是实对称矩阵,存在谱分解
QT(QAD−1)TB(QAD−1)Q=Λ,而
QT(QAD−1)TA(QAD−1)Q=QTQ=In。所以存在可逆矩阵
T=QAD−1Q使得
TTAT和
TTBT同时是对角矩阵。
6.2.20 设A,B是n阶实对称矩阵,A,B半正定。证明,存在可逆矩阵T,使得TTAT和TTBT同时是对角矩阵。
证明:
A,B半正定
⟹任意
x=0,xTAx≥0,xTBx≥0。分两种情况讨论:
(a)若
N(A)∩N(B)={0}.
对任意
x=0,xT(A+B)x=xTAx+xTBx>0⟹A+B正定。显然
A+B也是对称矩阵,由6.2.19,存在可逆矩阵
T,使得
TT(A+B)T=Λ1和
TTAT=Λ2同时是对角矩阵。而
TTBT=TT(A+B)T−TTAT=Λ1−Λ2也是对角矩阵。所以存在可逆矩阵
T,使得
TTAT和
TTBT同时是对角矩阵。
(b)若
N(A)∩N(B)={0}.
取
N(A)∩N(B)的一组正交基
(q1,⋯,qs),并扩张为正交矩阵
Q,则
QTAQ=[OA1],QTBQ=[OB1].
显然此时
N(A1)∩N(B1)={0},问题即转化为(a)中情况。