4.3.1
略.
4.3.2 利用按照行(列)展开求下列行列式;按照哪一行(列)展开,使得计算最简单?
按照0最多的行(列)展开即可.
4.3.3 计算a1b1c1d1e1a2b2c2d2e2a3b3000a4b4000a5b5000.
解:3-5行一定线性相关,行列式为0.
4.3.4
略.
4.3.5 计算λ−1λ⋱⋱−1λ−1anan−1⋮a2λ+a1.
解:化为行列式立得,原式
=λn+λn−1a1+λn−2a2+⋯+an.
4.3.6 计算λn1λ⋱2⋱⋱⋱⋱2⋱λ1nλ.
解:将第
n行加到第
n−1行,再将第
n−1行加到第
n−2行,以此类推。原式=
λ+nnλ+nλ+n−1n−1⋯λ+nλ+n−2⋱⋯⋯λ+n⋱2λ+nλ+n⋯⋱λ+11λ+nλ+nλ+nλ+nλ
=(λ+n)1n1λ+n−1n−1⋯λ+nλ+n−2⋱⋯⋯λ+n⋱21λ+n⋯⋱λ+111λ+nλ+nλ+nλ
用第
n列减第
n−1列,再用第
n−1列减第
n−2列,以此类推。
=(λ+n)1n0λ−1n−1⋯1λ−1⋱⋯2⋱20⋱λ−110n−1λ−1=(λ+n)λ−1n−11λ−1⋱2⋱2⋱λ−11n−1λ−1
设原式=
Sn(λ).上面的计算表明
Sn(λ)=(λ+n)Sn−1(λ−1).所以原式=
(λ+n)(λ+n−1)⋯(λ+2)S1(λ−(n−1))=(λ+n)(λ+n−1)⋯(λ+2)λ−(n−1)11λ−(n−1)=(λ+n)(λ+n−1)⋯(λ+2)(λ−n+2)(λ−n)
4.3.7 回顾例1.7.4中的对称、上三角和下三角Pascal矩阵。在四阶的情形,这三种矩阵分别为S4=1111123413610141020,U4=1000110012101331,L4=1111012300130001.另外,存在LU分解Sn=LnUn.
1. 求
det(Ln),det(Un),det(Sn).
2. 求
Sn右下角元素的代数余子式。
3. 将
Sn右下角的元素减
1得到矩阵
An,求
det(An).
1. 解:
det(Ln)=det(Un)=det(Sn)=1.
2. 解:
Cn,n=(−1)n+ndet(Sn−1)=det(Sn−1)=1.
3. 解:
令
Sn=s11s21⋮sn1s12s22⋮sn2⋯⋯⋯s1ns2n⋮snn。按照最后一行展开,则
det(Sn)=sn1Cn1+sn2Cn2+⋯+snnCnn,同理
det(An)=sn1Cn1+sn2Cn2+⋯+(snn−1)Cnn。所以
det(Sn)−det(An)=Cnn=det(Sn−1)=1,即
det(An)=0.
4.3.8 给定An=2−1−12⋱⋱⋱−1−12,Bn=1−1−12⋱⋱⋱−1−12.
1. 利用展开式得到
det(Bn)关于
n的递推关系,并计算
det(Bn).
2. 利用
det(An)与
det(Bn)的关系计算
det(An).
1. 解:
按照第一行展开后根据行列式的列线性,有
det(Bn)=1−1−12−1−12⋱⋱⋱−1−12=2−1−12⋱⋱⋱−1−12+−1−12⋱⋱⋱−1−12=1−1−12⋱⋱⋱−1−12+10−12⋱⋱⋱−1−12+−1−12⋱⋱⋱−1−12=1−1−12⋱⋱⋱−1−12+0=det(Bn−1)
而
det(B2)=1,于是
det(Bn)=1.
2. 解:
同样利用列线性,
det(An)=2−1−12−1−12⋱⋱⋱−1−12=1−1−12−1−12⋱⋱⋱−1−12+1−12−1−12⋱⋱⋱−1−12=det(Bn)+det(An−1)=1+det(An−1)
而
det(A2)=3,于是
det(An)=n+1.
4.3.9 (行列式中的Fibonacci数列) 如果一个矩阵比上(下)三角矩阵仅仅多一排非零对角元,则称之为上(下)Hessenberg矩阵。例如,H4=2100121011211112就是上Hessenberg矩阵。上Hessenberg矩阵在数值分析中很有用。
1. 令
Hn为
n阶上Hessenberg矩阵,其对角元素都是
2,其他非零元素都是
1。证明
det(Hn+2)=det(Hn+1)+det(Hn),即这些行列式组成了Finonacci数列。
2. 令
Sn是对角元素为
3,与对角元相邻的元素为
1的
n阶三对角矩阵(见练习1.3.12),例如,
S4=3100131001310013。它既是上Hessenburg矩阵也是下Hessenburg矩阵。求
Sn的递归公式,并分析与Fibonacci数列的关系。
3. 设
n阶三对角矩阵的行列式的完全展开式中,最多有
tn项非零,求
tn的递归公式。
1. 证明:
det(Hn)=2112111⋱⋱⋯⋯⋱⋱111⋮12n×n
按照第一列展开,
2112111⋱⋱⋯⋯⋱⋱111⋮12n×n=2211⋱⋱⋯⋱⋱11⋮12(n−1)×(n−1)−1112⋱⋯⋱⋱11⋮12(n−1)×(n−1)=2det(Hn−1)−1112⋱⋯⋱⋱11⋮12(n−1)×(n−1)
由列线性,
1112⋱⋯⋱⋱11⋮12(n−1)×(n−1)=2112⋱⋯⋱⋱11⋮12(n−1)×(n−1)−10121⋯⋱⋱⋱⋯⋱⋱11⋮⋮12(n−1)×(n−1)=2112⋱⋯⋱⋱11⋮12(n−1)×(n−1)−2112⋱⋯⋱⋱11⋮12(n−2)×(n−2)=det(Hn−1)−det(Hn−2)
于是
det(Hn)=2det(Hn−1)−(det(Hn−1)−det(Hn−2))=det(Hn−1)+det(Hn−2)
而
det(H1)=2,det(H2)=3, 所以
det(Hn) 是 Fibonacci 数列。
2. 解:
det(Sn)=3113⋱⋱⋱113
按照第一列展开,
3113⋱⋱⋱113n×n=33113⋱⋱⋱113(n−1)×(n−1)−113⋱1⋱⋱⋱⋱113(n−1)×(n−1)=3det(Sn−1)−311⋱⋱⋱⋱113(n−2)×(n−2)=3det(Sn−1)−det(Sn−2)
由
det(S1)=3,det(S2)=8得
det(Sn) 是 Fibonacci 数列的第
2n+2项。
3. 解:
设
n−1阶三对角矩阵
Sn−1的完全展开式有
tn−1项非零。现考虑
Sn=Sn−111snn 的完全展开式。显然
snn和
Sn−1的完全展开式中所有非零项相乘也非零,这里有
tn−1个非零项。最右一列的
1和最下一行的
1和
Sn−2的完全展开式中所有非零项相乘也非零,这里有
tn−2个非零项。于是
tn=tn−1+tn−2.
4.3.10 求下列推广的Vandermonde行列式:
11⋮1x1x2⋮xn⋯⋯⋯x1n−2x2n−2⋮xnn−2x1nx2n⋮xnn.
解:
可以看出原式和Vandermonde行列式之间差一列
n−1次项,我们不妨把
n−1次项补上去,同时为了保持方阵的形式还需要另外加一行。设得到的Vandermonde矩阵是
D,即,
D=11⋮11x1x2⋮xny⋯⋯⋯⋯x1n−2x2n−2⋮xnn−2yn−2x1n−1x2n−1⋮xnn−1yn−1x1nx2nxnnyn.
且
det(D)=(y−x1)(y−x2)⋯(y−xn)(xn−x1)(xn−x2)⋯(x2−x1).
将
D按照最后一行展开可以发现,原式即为
det(D)中
yn−1系数的相反数。所以原式=
(x1+x2+⋯+xn)1≤j<i≤n∏(xi−xj).
4.3.11 对n阶方阵A,证明,det(λIn−A)是λ首项系数为1的n次多项式.
证明是显然的,略.
4.3.12
1.设
A是正交矩阵,且
det(A)<0,证明
In+A不可逆.由此可得,存在非零向量
x,使得
Ax=−x.
2.设
A是奇数阶正交矩阵,且
det(A)>0,证明
In−A不可逆.由此可得,存在非零向量
x,使得
Ax=x.偶数阶的情形,结论是否成立?
1.证明:
det(AT(A+I))=det(A)det(A+I)=det(I+AT)=det(I+A).
即
det(A)det(A+I)=det(A+I).
又因为
det(A)<0,所以
det(A+I)=0,即
In+A不可逆.
2.证明:
A是正交矩阵且
det(A)>0⇒det(A)=1.
所以
det(AT(In−A))=det(In−A)
而
det(AT(In−A))=det(AT−In)=(−1)ndet(In−A).
即
det(In−A)=(−1)ndet(In−A).
当
n是奇数时,
det(In−A)=−det(In−A)⇒det(In−A)=0.而当
n是偶数时不成立.
4.3.13 设n阶方阵A=[aij]对角占优,且对角元素都是正数,证明det(A)>0.
证明:
注意到
det(A+tI)是关于
t的多项式函数,处处连续可导.
A+tI(t>0)也是对角占优矩阵,对任意
t>0,
det(A+tI)=0.又因为
limt→∞det(A+tI)=+∞. 所以对任意
t>0,det(A+tI)>0.于是当
t=0时,
det(A+tI)=det(A)≥0.
又因为
A是对角占优矩阵
⇒det(A)=0,所以
det(A)>0.
4.3.14 证明若A不可逆,则其补矩阵的秩是0或1.
证明:
令
A=[a1a2⋯an],补矩阵的第
i列为
ci=ci1ci2⋮cin.
若
rank(A)≤n−2,去掉一列后的矩阵仍然不可逆,代数余子式全部为0.即
rank(C)=0.
若
rank(A)=n−1,因为
A不可逆
⇒det(A)=0.所以当
i=j时,
ajTci=0.当
i=j时
ajTci=det(A)=0.
这表明
ci⊥R(A),所以
rank(C)≤1.
即
rank(C)≤1.
4.3.15 给定所有元素全为整数的可逆矩阵A,证明,A−1的所有元素全为整数,当且仅当∣detA∣=1.
由4.3.7逆矩阵公式立得.
4.3.16
略.
4.3.18 设n阶方阵A,B满足AB=BA,证明:
1.
∣det([A−BBA])∣=∣det(A2+BB)∣
2.
det([A−BBA])=det(A2+B2)
1.证明:
∣det([A−BBA])∣=det([A−BBA][AB−BA])=det([A2+B200A2+B2])=∣det(A2+B2)∣.
2.证明:
法一:
若
A可逆,
det([A−BBA])=det([AOBA+BA−1B])=det(A2+ABA−1B)=det(A2+B2).
若
A不可逆,
det(A)=0,则
det(A+ϵI) 一定不等于0.其中
ϵ是任意小的正实数.(这是因为此时行列式是关于
ϵ的多项式函数).
显然有
det([A+ϵI−BBA+ϵI])=det((A+ϵI)2+B2).
因为行列式是关于
ϵ的多项式函数,处处连续可导,所以:
det([A−BBA])=ϵ→0limdet([A+ϵI−BBA+ϵI])=ϵ→0limdet((A+ϵI)2+B2)=det(A2+B2).
法二:
利用行变换和列变换
A−BBA=行变换A−B−iABA−iB=列变换A+iBOBA−iB=det((A+iB)(A−iB))=det(A2+B2)
4.3.19-4.3.20
略.
4.3.21 在空间R3中,证明由向量a1,a2围出的平行四边形的面积是det(ATA),其中A=[a1a2.]
证明:
det(ATA)=det((QR)TQR)=det(RTR)=∣det(R)∣.
其中
R是上三角阵,且
R的对角元分别是
a1的长度和
a2垂直于
a1分量的长度.于是
S=∣det(R)∣=det(ATA).
4.3.22 对n阶方阵A在1,2,⋯,n中任取k个数i1<i2<⋯<ik,A的第i1,i2,⋯,ik行与i1,i2,⋯,ik列的交叉点处元素组成的k阶方阵,称为A的一个k阶主子阵,其行列式称为A的一个k阶主子式。显然A有Cnk个主子阵,也有Cnk个主子式。证明:
1.
det(A+λIn)是关于
λ的
n次多项式.
2.记
det(A+λIn)=a0λn+a1λn−1+⋯+an,则
a0=1,而对
k=1,2,⋯,n,
ak是
A的所有
k阶主子式的和.
4.3.23 (行列式游戏)甲和乙两个人构造一个n阶方阵,轮流填写矩阵中的元素,直到填满为止。如果矩阵的行列式非零,则甲胜,否则,乙胜。
1.如果乙先开始,且
n=2,则谁有必胜策略。
解:显然无论是谁先手,
n=2时乙必胜。乙只需要在任意一行(列)全填0即可。当甲先手时,若甲下(1,1),乙下(2,2)就能保证胜利。
2.如果乙先开始,且
n=3,则谁有必胜策略。
解:乙有必胜策略,乙通过填0即可获胜。乙先在任意一个位置填0,接下来无论甲在任何地方下另一个数,乙只需要在与自己刚刚填写的数字位于同一列的另一个空格处填0,问题就转变为双方争夺
n=2的余子阵的问题,而在
n=2的情况下乙必胜。例子:乙在(2,1)填0,甲有两种选择(a)在(2,2)或(2,3)填1;(b)在其他列填1.若甲选(a),乙就在同一列另外一个空格处填0,问题即转变为
n=2情况。若甲选(b),乙就在(2,2)处填0,接下来甲只能在(2,3)填1,否则立输。此时问题显然也变为
n=2的情况。乙必胜。
3.如果甲先开始,且
n=3,则谁有必胜策略。
解:乙有必胜策略。和2中讨论的一样,乙通过填0即可获胜。
4.如果甲先开始,且
n是偶数,则谁有必胜策略。
解:乙有必胜策略,乙通过在对应行填和甲一样的数构造两个相同的行即可获胜。例子:甲在(2,1)填3,乙也在(2,2)填3,甲在(3,2)填5,乙也在(3,1)填5...以此类推,就能得到两个相同的行。乙必胜。
感兴趣的读者可自行思考
n是偶数且乙先手或
n是奇数的情况。相关讨论见
https://mathoverflow.net/questions/2193/variation-on-a-matrix-game