模拟2 解答
第1题(10分)
设 P2(t) 为次数 ⩽2 的多项式全体构成的线性空间。对任意 f(t),g(t)∈P2(t),定义内积
(f(t),g(t))=∫01f(t)g(t)dt.
(1)求多项式 1−t 的长度
解:
多项式 1−t 的长度即其在上述内积下的范数:
∥1−t∥=(1−t,1−t)=∫01(1−t)2dt=∫01(1−2t+t2)dt=[t−t2+3t3]01=31
(2)记 W=span{1−t},求 W⊥
解:
W⊥={p(t)∈P2(t)∣(p(t),1−t)=0}。设 p(t)=a+bt+ct2,则
(p,1−t)=∫01(a+bt+ct2)(1−t)dt=2a+6b+12c=0,
即 6a+2b+c=0。该方程有 2 个自由未知量,故 dimW⊥=2。
取两组无关整数解:
- a=1,b=0,c=−6 ⇒ p1(t)=1−6t2;
- a=0,b=1,c=−2 ⇒ p2(t)=t−2t2。
可验证 (p1,1−t)=0,(p2,1−t)=0,且 p1,p2 线性无关。
W⊥=span{1−6t2,t−2t2}
第2题(10分)
在线性空间 R2×2 中有下列 5 个向量:
α1=(1120),α2=(−2311),β1=(1300),β2=(1122),β3=(2111).
令 V=span{α1,α2},W=span{β1,β2,β3}。
(1)证明:α2∈W
解:
将 2×2 矩阵按列展开为 R4 中的向量。设
α2=x1β1+x2β2+x3β3,
等价于线性方程组(按列优先向量化):
130011222111x1x2x3=−2311.
对增广矩阵作行化简,得
130011222111−2311⟹1000010000104547−250,
有解。取特解 x1=45,x2=47,x3=−25,则
α2=45β1+47β2−25β3∈W.
验算:
45(1300)+47(1122)−25(2111)=(−2311)=α2.
故 α2∈W。
(2)判别 V+W 是否是直和,并说明理由
解:
V+W 为直和当且仅当 V∩W={0}(等价地,V+W 中每个向量分解 v+w 唯一)。
由(1),α2∈W;又 α2∈V=span{α1,α2},故 α2=0 且 α2∈V∩W。
因此 V∩W={0},V+W 不是直和。
第3题(10分)
已知 4×5 矩阵
A=110−2110−2422020248440,
经初等行变换,得行阶梯形
10001000420022008400.
(1)求 A 的满秩分解和 A+
解:
由行阶梯形,主元在第 1、3 列,故 rank(A)=2。
继续化为行最简形:
100010000100−21000200.
读出列向量关系(记 ci 为 A 的第 i 列):
c2=c1,c4=c3−2c1,c5=2c3.
取两个线性无关列作 F,系数矩阵作 G,得满秩分解 A=FG(F∈R4×2,G∈R2×5):
F=110−24220,G=(101001−2102).
验证:FG=A。
由满秩分解 A=FG,Moore–Penrose 广义逆为
A+=GT(GGT)−1(FTF)−1FT.
先算
GGT=(6−2−26),(GGT)−1=321(6226),
FTF=(66624),(FTF)−1=1081(24−6−66).
代入计算,得
A+=961961321961161288728872885−3211445−721−721721241361−14411−14411−1441487−721.
(2)求 Ax=b 的极小范数最小二乘解,其中 b=(1,1,1,1)T
解:
因 rank(A)=2<4,而 b∈R4,方程组一般不相容。极小范数最小二乘解为
x=A+b.
计算:
x=(−181,−181,181,61,91)T
此时 Ax=AA+b,为 b 在 R(A) 上的正交投影。
(3)写出 A+A 的 Smith 标准型
解:
先算
A+A=163163161−16581163163161−1658116116116316183−165−1651611611818181838143.
A+A 为对称幂等阵,rank(A+A)=2,特征值为 1,1,0,0,0。
在域 R 上:秩为 2,Smith 标准型为
1000001000000000000000000
在 Z 上:A+A 含分母,先同乘 16 化为整数矩阵
B=16⋅A+A=331−52331−5211316−5−51112226212.
rank(B)=2。对 B 施以初等行、列变换(只做整数加减与换行换列),化为对角形:
第 1 步:将最小非零元 1 换至 (1,1) 位——R1↔R3:
BR1↔R3133−52133−52311161−5−5112622212.
第 2 步:以 a11=1 为 pivot,消去第 1 列其余元素:
R2←R2−3R1,R3←R3−3R1,R4←R4+5R1,R5←R5−2R1,
再作列变换 C2←C2−C1,C4←C4−C1,C5←C5−6C1,得
10000000003−8−81600−8−81600−16−16320.
第 3 步:在右下角子块中取最小元 8,R2←−R2,再 R4←R4+2R2 消去,得第二 pivot 8;继续列倍加消去同行其余元,并删去零行零列。
第 4 步:整理得
1000008000000000000000000
对角元 d1=1,d2=8,满足 1∣8,即为 B(亦即 16⋅A+A)在 Z 上的 Smith 标准型。
(4)求 ∥eA+A∥F
解:
A+A 为对称幂等阵(投影阵),特征值为 1,1,0,0,0。
矩阵指数与特征值同步:eA+A 的特征值为 e,e,1,1,1。
A+A 为正规阵(对称),Frobenius 范数等于特征值模长平方和的开方:
∥eA+A∥F=e2+e2+12+12+12=2e2+3≈4.2164.
(亦可利用幂等性:eA+A=I+(e−1)A+A,直接按定义求 Frobenius 范数,结果相同。)
第4题(10分)
已知
A=21112513961543594103268,b=21−1−2.
(1)求 A 的 LU 分解(L 为单位下三角阵,U 为上三角阵)
解:
用 Doolittle 方法,主元存入 U,乘数 lik=aik(k)/akk(k) 存入 L。
消第 1 列(u11=2):
l21=21,l31=21,l41=21.
更新后第 2 列主元 u22=4。
消第 2 列(u22=4):
l32=412=3,l42=48=2.
更新后第 3 列主元 u33=4。
消第 3 列:
l43=432=8.
得
L=1212121013200180001,U=20002400612404862
验算:LU=A。
(2)利用 LU 分解解 Ax=b
解:
令 Ly=b(前代),再令 Ux=y(回代)。
前代 Ly=b:
⎩⎨⎧y1=2,y2=1−21y1=0,y3=−1−21y1−3y2=−2,y4=−2−21y1−2y2−8y3=13.
回代 Ux=y:
⎩⎨⎧x4=2y4=213,x3=4y3−6x4=4−2−39=−441,x2=4y2−12x3−8x4=40+123−52=471,x1=2y1−2x2−6x3−4x4=22−271+2123−26=1.
验算 Ax=b 成立。
x=(1,471,−441,213)T