本页属“补充篇 · 数值线性代数工具箱”,不是课堂第 2 章;课堂第 2 章见线性方程组的解法。
用法:遮住答案,在空白纸上写;写完再展开对照。每张卡末尾标记 ✓ 能完整写出 / △ 有关键词但缺条件 / × 完全想不起,下一轮只重做 △ 与 ×。全页 46 张卡,一轮约 35 分钟。本页负责检索;压缩看公式与结论卡,辨析看易混对照卡,考场节奏看综合自测卷;第一章的对应卡片见第一章回忆卡。
A1. 写出内积的三条公理与诱导范数;除 xTy 外再举一族合法的内积。
解答
对称性 ⟨x,y⟩=⟨y,x⟩;对第一变元线性(对第二变元由对称性同理);正定性 ⟨x,x⟩≥0 且为零当且仅当 x=0。诱导范数 ∥x∥=⟨x,x⟩、距离 d(x,y)=∥x−y∥。族:加权内积 ⟨x,y⟩D=xTDy(D 对称正定;对角正权只是特例)(定义 2.1.1、例 2.1.1)。(✓ / △ / ×)
A2. 「正交组」「标准正交组」「标准正交基」三者的差别是什么?零向量算不算正交?
解答
正交组=非零向量组且两两内积为零;标准正交组=正交组中每个向量范数为 1;标准正交基=含 dimV 个向量的标准正交组。零向量与任何向量约定为正交(定义 2.1.2)。(✓ / △ / ×)
A3. 正交矩阵与「列正交矩阵」怎么定义?二者差在哪?
解答
方阵 Q 满足 QTQ=I 称正交矩阵;长方阵 Q∈Rm×n(m≥n)满足 QTQ=In 只称列正交矩阵——此时 QQT=Im,它是一个投影矩阵(QR 里的 Q 就是这一类)(定义 2.1.3)。(✓ / △ / ×)
A4. 写出四个基本子空间的定义,并标明各自住在哪个空间里。
解答
列空间 C(A)={Ax}⊆Rm;零空间 N(A)={x:Ax=0}⊆Rn;行空间 C(AT)⊆Rn;左零空间 N(AT)⊆Rm(定义 2.2.1)。(✓ / △ / ×)
A5. 正交补 S⊥ 的定义是什么?检验 x∈S⊥ 时要不要对 S 的每个向量都查一遍?
解答
S⊥={x:⟨x,s⟩=0 ∀s∈S}——与 S 中每个向量都正交的向量全体。实际只需查 S 的一组基:内积对第二变元线性,基上为零则全体为零(定义 2.2.2)。(✓ / △ / ×)
A6. 什么叫最小二乘解?什么叫残差?二者分别可不可以算出来?
解答
x^ 最小化 ∥Ax−b∥2 时称最小二乘解;r=b−Ax^ 称残差。残差可以算(b 与 x^ 都在手上),真解与误差 x−x^ 一般算不出来(定义 2.2.3)。(✓ / △ / ×)
A7. Householder 变换的公式与几何意义是什么?反射向量与镜面分别是谁?
解答
H=I−vTv2vvT:关于超平面 v⊥ 的镜像——v 方向翻转(Hv=−v),v⊥ 上不动(定义 2.3.1、定理 2.3.3)。(✓ / △ / ×)
A8. Givens 旋转怎么定义?它与 Householder 的关键区别(消元粒度、行列式)是什么?
解答
单位阵基础上把 (i,i),(i,j),(j,i),(j,j) 四个位置换成 c,s,−s,c(c2+s2=1)。它只作用在两个坐标上、一次消一个元素、detG=+1(保定向);Householder 一次改动整个子块、detH=−1(反定向)(定义 2.3.2、2.3.5)。(✓ / △ / ×)
A9. 矩阵条件数怎么定义?谱条件数等于什么?它衡量的是「元素多大」还是别的?
解答
cond(A)=∥A∥∥A−1∥(依赖所选范数,报告必须注明);cond2(A)=σmax/σmin=σ1/σr。它量的是离奇异有多远(行列向量之间的接近平行程度),对整体缩放免疫:cond(cA)=cond(A)(定义 2.4.1、定理 2.4.1)。(✓ / △ / ×)
A10. 特征值、谱、谱半径分别怎么定义?实矩阵的特征值一定是实数吗?
解答
Av=λv(v=0)时 λ 为特征值,λ⟺det(A−λI)=0;全部特征值(含重数)构成谱 σ(A);ρ(A)=max∣λi∣。实矩阵也可以有复特征值(旋转矩阵即例);只有对称矩阵保证全实(定义 2.5.1、定理 2.5.1)。(✓ / △ / ×)
A 组的适用条件
- 内积公理只保证「对称 + 线性 + 正定」,任何满足它们的二元函数都是内积;换内积就换正交、换投影、换最优解。
- 「正交」永远相对某个内积而言:写 xTy=0 时默认 Euclid 内积,切换内积后这句话不再是判据。
- 长方阵的 Q 只有列正交,「逆等于转置」只对方阵成立。
- 条件数依赖范数(cond1,cond2,cond∞ 可差 n 倍量级),但「丢几位」的结论一致——这是范数等价性在起作用。
- 谱半径 ρ 与范数 ∥A∥ 是两个不同的数,只在正规矩阵上相等。
B1. 内积诱导的范数自动满足哪些性质?怎么判断一个范数不来自内积?
解答
自动满足 Cauchy–Schwarz、三角不等式、正定与齐次——即确实是范数。判据是平行四边形恒等式 ∥x+y∥2+∥x−y∥2=2∥x∥2+2∥y∥2:满足它才来自内积(内积由极化恒等式反算回来)。∥⋅∥1 与 ∥⋅∥∞ 在 x=(1,0),y=(0,1) 上立刻失败(定理 2.1.1、定理 2.1.2、例 2.1.2)。(✓ / △ / ×)
B2. 正交组为什么必线性无关?标准正交基下「求坐标」与「求长度」各退化成什么?
解答
设 ∑icivi=0,两边与 vj 作内积,其余项因正交消失,得 cj∥vj∥2=0,故 cj=0。标准正交基下坐标就是内积 αi=⟨x,qi⟩,长度按 Parseval ∥x∥2=∑i∣⟨x,qi⟩∣2,正交向量之和满足勾股定理(定理 2.1.3、定理 2.1.4)。(✓ / △ / ×)
B3. 正交矩阵有哪些等价刻画(至少写六条)?正交变换为什么「不放大误差」?
解答
QTQ=I;QQT=I;Q−1=QT;列(行)向量组标准正交;保内积、保长度、保距离;∣detQ∣=1;奇异值全为 1(∥Q∥2=cond2(Q)=1);复特征值模全为 1。误差层面:∥QA∥2=∥A∥2、∥QA∥F=∥AP∥F=∥A∥F——放大因子恒为 1,误差与矩阵大小同比例、不随变换次数累积(定理 2.1.5、定理 2.1.6)。(✓ / △ / ×)
B4. 写出行秩=列秩、两个维数公式与两对正交补,并说出它们是哪种分解。
解答
rankA=dimC(A)=dimC(AT)=r;dimN(A)=n−r、dimN(AT)=m−r;C(AT)=N(A)⊥、C(A)=N(AT)⊥,即 Rn=C(AT)⊕N(A)、Rm=C(A)⊕N(AT) 两个正交直和。由此任何 b 唯一分解为 b=p+e(p∈C(A)、e∈N(AT)),且 e⊥p(定理 2.2.1、定理 2.2.2)。(✓ / △ / ×)
B5. 正交分解定理说了几件事?「投影」与「最近点」是什么关系?
解答
① 每个 x 唯一写成 x=p+e(p∈S、e∈S⊥);② 勾股 ∥x∥2=∥p∥2+∥e∥2;③ 有标准正交基时 p=∑i⟨x,qi⟩qi。最近性:∥x−s∥>∥x−p∥ 对任何 s=p 成立——投影就是最近点,最近 ⟺ 余项垂直(定理 2.2.3、定理 2.2.4)。(✓ / △ / ×)
B6. 投影矩阵 P=A(ATA)−1AT 的三条性质是什么?反过来如何判定一个矩阵是正交投影?
解答
对称、幂等、特征值只取 0/1 且 trP=rankP=rankA(P=O 时 ∥P∥2=1)。反之「对称 + 幂等」就是某个子空间上的正交投影——这是最省事的指纹判据;只幂等而不对称的是斜投影,余项不再垂直于子空间(定理 2.2.5、2.2.3 指纹框)。(✓ / △ / ×)
B7. 最小二乘解的存在性与唯一性各靠什么?正态方程的三件套是什么?
解答
总存在(ATb∈C(AT)=C(ATA));唯一 ⟺ A 列满秩(此时 ATA 非奇异)。三件套:① ATAx^=ATb⟺ATr=0⟺ 最小二乘解;② x^=(ATA)−1ATb、Ax^=Pb、r∈N(AT);③ ∥r∥2=∥b∥2−∥Ax^∥2(勾股)(定理 2.2.6)。(✓ / △ / ×)
B8. 加权最小二乘的解是什么?它与「换内积」是什么关系?与普通最小二乘有什么关系?
解答
x^=(ATWA)−1ATWb,满足加权正态方程 ATW(Ax−b)=0;几何上就是 b 在 C(A) 上、按加权内积 ⟨⋅,⋅⟩W 作正交投影的结果——公式结构不变,只换内积。又若 W=LTL,则加权最小二乘 ⟺ 方程两边同乘 L 后做普通最小二乘(定理 2.2.7)。(✓ / △ / ×)
B9. Gram–Schmidt 给出什么?QR 分解什么时候唯一?「分层张成」为什么重要?
解答
给出 A=QR:Q=[q1,…,qn] 列正交、R 上三角且 rjj=∥v∥2>0;在「R 对角元全为正」的规范化下分解唯一(否则可按列成对变号)。分层张成 span{q1..qk}=span{a1..ak} 是 QR 区别于 SVD 的关键性质,也是 QR 迭代求特征值的基础;秩亏时改用列主元 QR(AP=QR,∣r11∣≥⋯)(定理 2.3.1、定理 2.3.2)。(✓ / △ / ×)
B10. 四种正交化方法(CGS / MGS / Householder / Givens)的正交性损失各是什么量级?怎么选?
解答
CGS:cond2u(cond≳u−1/2≈108 后彻底失效,误差饱和在 O(n)——最坏情形);MGS:condu;Householder 与 Givens:u(与 cond 无关)。选法:稠密用 Householder(LAPACK geqrf、NumPy 默认),稀疏 / 带状 / 局部修改用 Givens(一次只动两行),需要显式 Q 且条件数温和时可用 MGS,CGS 只用于教学(2.3.3 对比表、定理 2.3.3、定理 2.3.4)。(✓ / △ / ×)
B11. 用 QR 解最小二乘要解什么方程?为什么它比正态方程好,好在哪个量级?
解答
把 A 的全 QR 代入:∥Ax−b∥2=∥Rx−c∥2+∥d∥2(c 为 QTb 前 n 行),故解上三角方程 Rx^=c(回代 O(n2)),最小残差范数 =∥d∥2。全程不形成 ATA,误差量级从正态方程的 cond2u 降到 condu——cond≳108 时必须弃用正态方程(定理 2.3.5、2.3.6)。(✓ / △ / ×)
B12. 条件数的五条基本性质 + 两条扰动界 + 「丢几位」的估计。
解答
cond(A)≥1(算子范数下);cond(A−1)=cond(A);cond(cA)=cond(A);cond2(Q)=1;cond(AB)≤cond(A)cond(B)(一般不是等式)。右端扰动 ∥x∥∥δx∥≤cond∥b∥∥δb∥(紧界);系数扰动 ≤1−condρAcond(ρA+ρb),要求 cond⋅ρA<1。丢失位数 k≈log10cond;cond≳1016=u−1 时双精度形同奇异(定理 2.4.1–定理 2.4.3、2.4.2)。(✓ / △ / ×)
B13. 「残差小」为什么推不出「误差小」?正确的判据是什么?
解答
因为 e=x^−x=−A−1r,故 ∥x∥∥e∥≤cond(A)∥b∥∥r∥:两者之间隔着条件数。正确判据是 cond⋅∥r∥/∥b∥,不是 ∥r∥(实测:残差相对量 2.5×10−12 时误差可达 50%)。求解器能保证的是残差(后向),用户想要的是误差(前向)(定理 2.4.4)。(✓ / △ / ×)
B14. 谱半径与范数是什么关系?Gerschgorin 圆盘定理说了什么、怎么用?
解答
ρ(A)≤∥A∥ 对任何相容范数成立,但可远小于范数(Jordan 块:ρ=1、∥⋅∥2≈1000)。Gerschgorin:∣λ−aii∣≤ri=∑j=i∣aij∣(行版本,列版本对 AT 用),全部特征值落在圆盘并集内;若 k 个圆盘构成与其他圆盘不相交的连通区域,则其中恰含 k 个特征值;行、列两版本取交集可收紧围栏(注意:逐盘配对取交只是启发式)。严格对角占优 ∣aii∣>ri ∀i ⟹ 0∈/∪Di ⟹ 非奇异(定理 2.5.2、2.5.2)。(✓ / △ / ×)
B15. Rayleigh 商的两个极值性质是什么?Courant–Fischer 说了什么?幂法的收敛速度由谁决定?
解答
R(x)=xTAx/xTx 是特征值的加权平均,故 λmin=minR、λmax=maxR,极值点恰是对应特征向量。Courant–Fischer:λk=mindimS=kmax0=x∈SR(x)(对 −A 用一次得第二等式)——特征值问题由此变成优化问题,Lanczos / 子空间迭代都从它出发。幂法:向量的方向误差 O(∣λ2/λ1∣k),Rayleigh 商误差 O(∣λ2/λ1∣2k)——同样精度所需步数减半(定理 2.5.3–定理 2.5.5)。(✓ / △ / ×)
B16. 特征值对扰动有多敏感?对称情形、一般情形、重特征值情形各是什么结论?
解答
对称(Weyl):∣λi(A+E)−λi(A)∣≤∥E∥2,1-Lipschitz,一点都不放大;一般可对角化(Bauer–Fike):≤cond2(V)∥E∥2,放大因子就是特征向量矩阵的条件数(对称时可取 V 正交,退化为情形一);非正规 k 重特征值:位移 O(∥E∥1/k)(Jordan 块 k=2 时精确等于 t)——固有病态,任何算法都救不了(定理 2.5.6)。(✓ / △ / ×)
B17. SVD 一次性交出哪些东西?(四个子空间的基 / 三个量 / 最优低秩逼近)
解答
A=UΣVT,σ1≥⋯≥σr>0:C(A)=span{u1..ur}、N(AT)=span{ur+1..um}、C(AT)=span{v1..vr}、N(A)=span{vr+1..vn}(都是标准正交基,不必再补 Gram–Schmidt);∥A∥2=σ1、∥A∥F=∑σi2、cond2=σ1/σr、rankA=#{σi>0};Eckart–Young:Ak=∑i≤kσiuiviT 是最优秩 k 逼近,∥A−Ak∥2=σk+1 精确成立(定理 2.6.1–定理 2.6.4)。(✓ / △ / ×)
B18. 伪逆怎么定义、有哪些关键性质?为什么「不截断的伪逆」很危险?数值秩怎么定?
解答
A+=VΣ+UT(非零 σi 取倒数、零仍记零):满足 Penrose 四条件且唯一;AA+、A+A 分别是到 C(A)、C(AT) 的正交投影;x^=A+b 是唯一的最小范数最小二乘解(秩亏时最小的那一点,正态方程只能给出「直线上任意一点」);列满秩时 A+=(ATA)−1AT。危险:1/σr 会把相应方向放大 cond(A) 倍,必须按噪声水平截断。数值秩 #{i:σi>tσ1},t 按数据噪声水平取(matrix_rank 的默认阈值 =max(m,n)u 只适合「精确矩阵数值降秩」)(定理 2.6.3、2.6.4)。(✓ / △ / ×)
B 组的适用条件与常见越界
- 平行四边形恒等式是「范数来自内积」的充要判据(有限维版 Jordan–von Neumann),不能只查一半。
- P=A(ATA)−1AT 的前提是 A 列满秩;秩亏时只能用伪逆或列主元 QR。
- 「QR 唯一」只在 R 对角元全为正的规范形式下成立;比较两套 Q,R 前先统一符号。
- ρ<1 只保证渐近收敛(ρ=limk∥Ak∥1/k),非正规矩阵可能先「瞬态增长」;稳妥判据是 ∥A∥<1。
- Eckart–Young 保证误差最小,不保证最优矩阵唯一(σk=σk+1 时最优子空间可旋转)。
C1. 诱导范数、平行四边形恒等式、极化恒等式。
解答
∥x∥=⟨x,x⟩;∥x+y∥2+∥x−y∥2=2∥x∥2+2∥y∥2;⟨x,y⟩=41(∥x+y∥2−∥x−y∥2)(定义 2.1.1、定理 2.1.2)。(✓ / △ / ×)
C2. 标准正交基下:坐标、Parseval、投影公式。
解答
αi=⟨x,qi⟩;∥x∥2=∑i∣αi∣2;p=∑i⟨x,qi⟩qi(定理 2.1.4、定理 2.2.3)。(✓ / △ / ×)
C3. 四个子空间的维数与两个正交直和。
解答
dimC(A)=dimC(AT)=r,dimN(A)=n−r,dimN(AT)=m−r;Rn=C(AT)⊕N(A)、Rm=C(A)⊕N(AT)(定理 2.2.1、定理 2.2.2)。(✓ / △ / ×)
C4. 投影、投影矩阵、正态方程、加权版——四式一行收口。
解答
p=A(ATA)−1ATb(记 P=A(ATA)−1AT);ATAx=ATb;x^=(ATA)−1ATb;加权 x^=(ATWA)−1ATWb(定理 2.2.5–定理 2.2.7)。(✓ / △ / ×)
C5. 三种正交化的算法式要点(CGS / MGS 系数、Householder 消元向量、Givens 的 c,s,r)。
解答
CGS:rij=qiTaj(原始列);MGS:rij=qiTv(当前列,就地更新);Householder:α=−sgn(x1)∥x∥2、v=x−αe1、H=I−2vvT/vTv(首分量取相加,避免相消);Givens:r=hypot(x1,x2)、c=x1/r、s=x2/r(别手写 x12+x22)(定理 2.3.1、定理 2.3.3、定理 2.3.4)。(✓ / △ / ×)
C6. QR 解最小二乘的两式 + 四种方法的正交性损失与复杂度。
解答
Rx^=(QTb)1:n、∥r∥2=∥d∥2。损失:CGS cond2u、MGS condu、Householder ≈u、Givens ≈u;复杂度:GS 2mn2、Householder 2mn2−32n3、Givens 约 6mn2(定理 2.3.5、2.3.5 对比表)。(✓ / △ / ×)
C7. 条件数全量速记(定义 + 两界 + 位数 + 残差关系)。
解答
cond=∥A∥∥A−1∥、cond2=σ1/σr;∥x∥∥δx∥≤cond∥b∥∥δb∥、∥x∥∥δx∥≤1−condρAcond(ρA+ρb);k≈log10cond;e=A−1r、∥x∥∥e∥≤cond∥b∥∥r∥(定理 2.4.1–定理 2.4.4)。(✓ / △ / ×)
C8. 谱侧四式:ρ≤∥A∥、Gerschgorin、幂法/Rayleigh 收敛率、扰动界。
解答
ρ(A)≤∥A∥;∣λ−aii∣≤∑j=i∣aij∣;向量误差 O(∣λ2/λ1∣k)、Rayleigh 商误差 O(∣λ2/λ1∣2k);对称 ∣Δλ∣≤∥E∥2、一般 ≤cond2(V)∥E∥2、非正规 k 重 O(∥E∥1/k)(定理 2.5.2、定理 2.5.3、定理 2.5.6)。(✓ / △ / ×)
C9. R(x) 的定义、Courant–Fischer 两个等式、反幂法要点。
解答
R(x)=xTAx/xTx;λk=mindimS=kmaxx∈SR=maxdimT=n−k+1minx∈TR;反幂法把 A 换成 (A−μI)−1,收敛到最接近 μ 的特征值,每步解一个方程组而不显式求逆(定义 2.5.2、定理 2.5.5、2.5.3)。(✓ / △ / ×)
C10. SVD 全量:分解、外积、范数、伪逆、Eckart–Young、与 PCA 的关系。
解答
A=UΣVT;A=∑i≤rσiuiviT;∥A∥2=σ1、∥A∥F=∑iσi2、cond2=σ1/σr;A+=VΣ+UT;minrankB≤k∥A−B∥2=σk+1、F 版 =∑i>kσi2;PCA:主方向 =V 的列、主成分方差 =σi2/(n−1)、σi(Xc)=(n−1)λi(C)(中心化后)(定理 2.6.1–定理 2.6.4、2.6.6)。(✓ / △ / ×)
C 组的适用条件
- 所有「ATA 型」公式(正态方程、cond2=σ1/σr 的 ATA 算法)都受 cond2 影响:能用 SVD / QR 就别显式形成 ATA。
- Householder 的符号选择、Givens 的
hypot、MGS 的「就地更新」——三处都是防相消/防溢出的约定,改一行就换稳定性量级。
- Eckart–Young 的等式 ∥A−Ak∥2=σk+1 对任何 k<r 都成立;用于压缩时还要看能量占比 ∑i≤kσi2/∑iσi2。
- PCA 的方差解释比、低秩去噪的阈值,都必须先中心化 / 定噪声水平,否则结论无意义。
D1. 「一组向量范数全为 1,就是标准正交基。」
解答
错。范数为 1 只保证「标准」;「正交」必须逐个验证两两内积。教材 例 2.1.3 的 r 组范数全对、r1Tr2=0 也对,但 r2Tr3=−2/18=0——套用坐标公式与 Parseval 会同时静默失效(拼回去 =(20/3,10/3,−1/3)=x)。(✓ / △ / ×)
D2. 「把 ∥⋅∥1 换成内积范数来定义正交,就能对 1-范数做投影与 SVD。」
解答
错。∥⋅∥1 与 ∥⋅∥∞ 不满足平行四边形恒等式(取 x=(1,0),y=(0,1):左端 8 或 2,右端都是 4),不来自任何内积,因此没有正交、投影、SVD 这套结构。见到「正交」二字,默认内积是 Euclid 内积(例 2.1.2、2.1.4 错误 2)。(✓ / △ / ×)
D3. 「P2=P 说明 P 是正交投影矩阵」;「长方阵满足 QTQ=In,所以 QQT=Im。」
解答
两处都错。① 正交投影要求对称 + 幂等;只幂等的是斜投影(例 S=(1010)),其余项不垂直于子空间。② QQT 是到 C(Q) 的投影矩阵,tr(QQT)=n<m,故 QQT=Im——「逆等于转置」只对方阵成立(2.2.6 错误 5、2.1.4 错误 4)。(✓ / △ / ×)
D4. 「Gram–Schmidt 就用经典版(CGS),它与修正版数学上等价」;「用 QR 解最小二乘,先算 ATA 更省事。」
解答
两处都错。① 精确算术下 CGS 与 MGS 结果相同,浮点下差一个 cond:cond≳108 时 CGS 的 ∥QTQ−I∥2 饱和到 O(n)(实测 n=16 时 6.55),MGS 为 4.7×10−4、Householder 为 3.0×10−15。② 形成 ATA 就把条件数平方(cond(ATA)=cond(A)2),误差由 condu 退化为 cond2u——实测 cond=6×1010 时正态方程误差 28(真值 1)、QR 仍为 2.2×10−6(2.3.6 实测表)。(✓ / △ / ×)
D5. 「最小二乘解总是存在且唯一,直接解正态方程就行。」
解答
错两处。① 存在性没问题(ATb∈C(ATA)),但唯一性要求 A 列满秩;秩亏时解是一整条仿射集,要用伪逆取最小范数解。② 病态时正态方程把 cond 平方;另外算完必须做 ATr=0 自检——它是最省事的错误探测器(残差各分量之和为零不代表 ATr=0,必须按 xi 加权求和)(2.2.6 错误 1、4)。(✓ / △ / ×)
D6. 「求解器报的残差是 10−15,所以解准到 15 位」;「同一个病态矩阵,换 Cholesky / QR / SVD 就能算准。」
解答
两处都错。① 残差与误差之间隔着条件数:∥e∥/∥x∥≤cond⋅∥r∥/∥b∥;实测残差 2.5×10−12 而误差 50%。② 病态是问题的属性:Hilbert 12 阶用任何解法误差都是 0.22,换算法无效——只能重述问题(正则化)、提高精度或迭代精化;另外先做行均衡化再看条件数(示例中 cond∞ 从 1020 降到 4×1010),且别手写 inv 算条件数(定理 2.4.4、2.4.6 错误 2、3、4、7)。(✓ / △ / ×)
D7. 「用 ρ(A) 代替 ∥A∥ 估收敛;ρ<1 就每步收缩」;「重特征值也能算准」;「Gerschgorin 圆盘给的是特征值近似值。」
解答
三处都错。① ρ≤∥A∥ 是不等式:N=(1010001) 的 ρ=1 而 ∥N∥2≈1000;非正规矩阵还会瞬态增长(∥Nk∥≈1000k),故稳妥判据是 ∥A∥<1。② 非正规 k 重特征值对 O(ε) 扰动移动 O(ε1/k):k=2 时双精度只剩约 8 位有效数字,是固有病态。③ 圆盘给的是包含区间,只有对角占优时圆心才接近特征值;别忘了列版本与取交集(定理 2.5.2、2.5.6 错误 1–4)。(✓ / △ / ×)
D8. 「σi(A)=∣λi(A)∣」;「先算 ATA 的特征值再开方就是 SVD」;「伪逆不用截断、数值秩用库函数默认阈值即可」;「PCA 不中心化也一样」。
解答
四处都错。① 奇异值总非负、共 min(m,n) 个;特征值可负可复。反例 A=(0010):λi=0 而 σ1=1。② cond(ATA)=cond(A)2——白丢一半有效数字;实现里用 np.linalg.svd(A)(Golub–Kahan 直接作用在 A 上)。③ 未截断的伪逆会把 1/σr 放大 cond 倍,必须按噪声水平设 rcond;matrix_rank 默认阈值 max(m,n)u 只适合精确矩阵(实测同一矩阵默认判秩 20、按噪声阈值只有 3)。④ 不中心化做 SVD,第一主方向会指向均值方向,所谓「解释方差」根本不是方差(2.6.7 错误 1–7、例 2.6.5)。(✓ / △ / ×)
D 组的典型误用清单(考前扫一眼)
① 范数全 1 就当标准正交;② 拿 1-/∞-范数做「正交」;③ 幂等就当正交投影、长方阵当成方阵用 Q−1=QT;④ 生产代码用 CGS、MGS 写成 CGS;⑤ 病态最小二乘用正态方程、算完不做 ATr=0 自检;⑥ 拿残差当误差、以为换算法能救病态、忘行均衡化;⑦ 拿 ρ 当 ∥A∥、指望算准重特征值、把 Gerschgorin 当精确值;⑧ 把奇异值当特征值、对 ATA 做特征值分解求 SVD、伪逆不截断、PCA 不中心化。
- 建卡组:打开 Anki → 新建卡组「数值分析-补充篇」(与「数值分析-第 1 章」并列,不要混在默认卡组里)。
- 逐条录入:正面放「编号 + 题干」,背面放「答案 + 依据」(带 [@编号] 与节号);纯文本批量导入用
Tab 分隔「正面 / 背面」,一行一张卡——本页 46 张按 A1–A10 / B1–B18 / C1–C10 / D1–D8 的顺序录。
- 复习节点(1 / 3 / 7 天):当天 → 隔天 → 第 3 天 → 第 7 天共四轮;每轮只重做标 △ 与 × 的卡。本页的错误几乎全部发生在条件上(列满秩、R 对角元为正、cond⋅ρA<1、噪声水平定阈值),凡「能背出结论但说不出条件」一律按 △ 处理。
- 加一道保险:C 组(公式卡)背面附一张手写照片;四轮后仍为 × 的卡抄进第二章的小结与陷阱对应的那一条,作为个人薄弱点单独跟踪。全章公式的两页压缩见公式与结论卡与易混对照卡。