像考试那样练,才不会在考场上手生。 这一页是一份闭卷、计时、100 分制的综合自测卷,
覆盖面从误差理论到插值(与课程进度一致),题量与难度贴近期末。它和必考点竞技场
的分工:竞技场练单点反射(快、碎、连击),这里练整卷节奏(慢、整、按点给分)。
做完先自己按文末评分标准打分,再按错题定位回正文与讲座页。
应试要点(老师课堂反复强调的批改口径)
- 必要 / 充分必须分清:写「充要条件」之前先问一句「反方向成立吗」——ρ(G)<1 是充要,∥G∥<1 只是充分;SOR 的 ω∈(0,2) 是必要。
- 前提条件缺不得:幂法的「严格按模最大」、牛顿法的「初值靠近根 / 四条件」——不写前提,结论就错。
- 一句话理由也要写:概念判断题只写「对/错」不给分,必须给理由或反例。
- 方向与范数:递推稳定性看放大因子;算特征值必须 2-范数;算特征向量范数随取。
- 公式要能自己推:证明题按步给分,写结论不写过程最多拿 2 分。
判断对错,并给出一句话理由或反例。只写「对 / 错」不得分。
1. (a) 有限维线性空间上任意两个向量范数都等价;(b) 这一结论对无穷维空间也成立。
2. 设 ∥A∥<1(某个从属范数),则 I−A 一定可逆。反之 I−A 可逆是否必有 ∥A∥<1?
3. 迭代式 x(k+1)=Gx(k)+d 收敛的充要条件是「存在某个范数使 ∥G∥<1」。
4. 对逐次超松弛(SOR)法,只要松弛因子落在 0<ω<2,迭代就一定收敛。
第一部分答案(每题 5 分:判断 2 分 + 理由 3 分)
1. (a) 对、(b) 错。 有限维空间上范数等价性来自单位球紧性;无穷维空间上不同范数可以不等价(典型:函数空间上的 ∥⋅∥∞ 与 ∥⋅∥1)。
2. 前半对,后半错。 ∥A∥<1⇒I−A 可逆,且 (I−A)−1=∑k=0∞Ak(Neumann 级数,因 ∥Ak∥≤∥A∥k 收敛)。反之不真:I−A 可逆只要求 1∈/σ(A),与范数大小无关——例如 A=10I,I−A=−9I 可逆,但 ∥A∥≫1。
3. 错。 收敛的充要条件是谱半径 ρ(G)<1。∥G∥<1 只是充分条件(因 ρ(G)≤∥G∥)——可以 ρ(G)<1 而所有常见范数都 ≥1(非正规矩阵)。
4. 错。 0<ω<2 是 SOR 收敛的必要条件,不是充分条件:区间外一定不收敛,区间内不保证收敛(还需 A 的性质,如对称正定或对角占优)。
5. a=2.71828 是 e 的近似值,a 有几位有效数字?(e=2.718281828…)
6. 设 A=(2−1−12)。求 ∥A∥1、∥A∥∞ 与谱半径 ρ(A)。
7. 写出牛顿迭代格式,并说明它为什么(至少)二阶收敛。
8. 用 n+1 个互异节点做多项式插值,插值多项式的次数最高是多少?「互异」保证了什么?
第二部分答案(每题 5 分)
5. 误差 ∣a−e∣≈1.83×10−6。写 a=0.271828×101,故 m=1;要求 21×101−n≤1.83×10−6:n=6 时 21×10−5=5×10−6≥1.83×10−6 ✓,n=7 时 5×10−7<1.83×10−6 ✗。故 6 位(从首位 2 数起,末尾的 8 也有效)。
6. ∥A∥1=max 列和 =max(2+1, 1+2)=3;∥A∥∞=max 行和 =3;特征值 λ=2±1,故 ρ(A)=3。这里 A 对称,∥A∥2=ρ(A)=3。
7. xk+1=xk−f′(xk)f(xk)。它等价于不动点迭代,迭代函数 φ(x)=x−f′(x)f(x),在根 s 处
φ′(s)=[f′(s)]2f(s)f′′(s)=0(f(s)=0, f′(s)=0),
由收敛阶定理(φ′(s)=0 而一般 φ′′(s)=0)得至少二阶。
8. 次数最高为 n 次。节点互异保证范德蒙德矩阵 detA=∏i<j(xj−xi)=0,从而方程组有唯一解,插值多项式存在且唯一。(注意:唯一性 = 数值可稳定求出——范德蒙德矩阵病态。)
9. 解线性方程组
⎩⎨⎧2x1+x2−x3=1,x1+3x2+2x3=13,x1−x2+3x3=8.
(a) 写出 A 的 LU 分解(或说明为何不必选主元);(b) 用前代 / 回代求出解;(c) 求 detA。
10. 设 A=(2001),取 u0=(1,1)T。按规范化幂法做一步(先单位化、再乘 A):写出 y0,u1,∥u1∥2,y1,并用 y1TAy1 估计最大特征值。
11. 判断下列不动点迭代 xk+1=φ(xk) 在给定区间上是否收敛,并说明依据:
(a) φ(x)=cosx 在 [0,1] 上;(b) φ(x)=x2−2 在 [1,2] 上。
第三部分答案(每题 12 分)
9. (a) 各阶顺序主子式 2>0、2113=5>0,均非零,不必选主元即可 LU 分解。记 L 单位下三角、U 上三角:
A=10.50.501−0.600120012.50−12.55.
(消元:r2←r2−0.5r1、r3←r3−0.5r1 得第 2 列为 2.5,−1.5;再 r3←r3−(−0.6)r2,得 U33=5。)
(b) 解 Ly=b 得 y=(1, 12.5, 15)T;回代 Ux=y:
x3=515=3,x2=2.512.5−2.5×3=2,x1=21−x2+x3=1.
(c) detA=detU=2×2.5×5=25。
✅ 解为 (1,2,3)。 回代验证:2+2−3=1 ✓、1+6+6=13 ✓、1−2+9=8 ✓。考场务必养成回代自查的习惯——本题数值设计得“刚好是整数”,一旦算出分数,多半是消元出了错。
10. ∥u0∥2=2,y0=21(1,1)T。
u1=Ay0=21(2,1)T,∥u1∥2=25,y1=∥u1∥2u1=51(2,1)T.
y1TAy1=51(2,1)(2001)(21)=51(2,1)⋅(4,1)=59=1.8 ⟶ 2=λ1 ✓
11. (a) 收敛。 φ([0,1])=[cos1,1]⊆[0,1](自映射,cos1≈0.5403),且 ∣φ′(x)∣=∣sinx∣≤sin1≈0.8415<1(压缩)。两条件满足,迭代收敛到 [0,1] 内唯一不动点 x∗=cosx∗≈0.7391。
(b) 不满足。 φ([1,2])=[−1,2]⊈[1,2](自映射失败:φ(1)=−1);且 ∣φ′(x)∣=∣2x∣≥2>1,压缩条件也失败。故定理不适用,迭代不保证收敛。
12. 设 ∥A∥<1(从属范数)。证明 I−A 可逆,且
(I−A)−1=k=0∑∞Ak=I+A+A2+⋯.
13. 设 f 二次连续可导、f(s)=0、f′(s)=0。证明牛顿迭代在 s 附近至少二阶收敛,即误差 ek=xk−s 满足 ek+1=O(ek2)。
第四部分答案(每题 12 分)
12. 由次乘性 ∥Ak∥≤∥A∥k,且 ∥A∥<1,故级数
Sn=I+A+⋯+An
满足 ∥Sn+1−Sn∥=∥A∥n+1→0,即 {Sn} 是 Cauchy 列,在(完备的矩阵空间里)收敛,记极限为 S=∑k≥0Ak。
计算 (I−A)Sn=I−An+1。令 n→∞,An+1→O(因 ∥A∥n+1→0),故 (I−A)S=I,即 I−A 可逆、S=(I−A)−1。∎
13. 把 f 在 xk 处展开(或用 φ 的 Taylor):牛顿迭代 xk+1=xk−f′(xk)f(xk)。取 φ(x)=x−f′(x)f(x),则
φ′(x)=[f′(x)]2f(x)f′′(x),φ′′(x)=f′(x)f′′(x)−[f′(x)]2f(x)f′′′(x)+⋯
(只要 f′(s)=0、f 足够光滑)。在 s 处 f(s)=0,故 φ′(s)=0;而一般 φ′′(s)=f′(s)f′′(s)=0。由收敛阶定理
ek+1=2φ′′(ξ)ek2=O(ek2),
即至少二阶收敛(当 f′′(s)=0 时恰为二阶)。∎
这一节考「把公式用到工程上」——期末应用题与平时大作业的原型。 每题先想清楚「用哪个方法、为什么、注意什么」,再对答案。
S1 · 海底温度剖面:采样点插值
潜水器沿深度 xi 每隔一段测一次温度 Ti,得到 6 个点 (xi,Ti),i=0,…,5。你只知道这 6 个数,却要给出任意深度的温度估计。
问:用什么方法?多项式次数取几?为什么不能直接把系数解出来?给出更稳的做法。
参考答案
用多项式插值:6 个点 ⇒ 至多 5 次多项式 P5,使 P5(xi)=Ti。理论上 (xi,Ti) 互异 ⇒ 插值多项式存在且唯一(范德蒙德矩阵非奇异)。
但不能直接解系数:系数矩阵是范德蒙德矩阵,节点一多或一靠近就严重病态,解出的系数会丢掉大量有效数字(「存在唯一 ≠ 能算出来」)。
更稳的做法:换基到拉格朗日基 li(x)=∏j=ixi−xjx−xj,直接写出 P5(x)=∑i=05Tili(x)——系数就是温度值,不用解方程组。工程上还会进一步用分段低次插值(样条)避免高次多项式的 Runge 振荡。
S2 · 量子力学能级:位移反幂法精修
自伴算子的离散化矩阵 A 是对称的,其特征值就是要找的能级。用渐近方法已得到第 k 个能级的近似 p,但精度不够。
问:如何把 p 精修到第 k 个能级的机器精度?写出算法步骤,并说明为什么不显式求逆。
参考答案
- 构造位移矩阵 A~=A−pI,其特征值为 λi−p;由 p 的选取,∣λk−p∣ 最小。
- 对 A~ 用反幂法(即对 A~−1 用幂法)求出模最小的特征值 μ=λk−p 及对应特征向量 xk。
- 加回位移:λk=μ+p。
不显式求逆:每步解线性方程组 A~uj+1=yj(用 LU 分解,分解可预先算一次并复用),因为「我们并不知道 A−1,只知道 A」。反幂法每迭代一步解一个方程组——这正是直接法/迭代法的用武之地。
S3 · 单摆大摆角:非线性方程求根
单摆大摆动周期含第一类完全椭圆积分,令 F(θ)=π2K(sin2θ)−c=0(c 由实测周期定出),要解出摆角 θ。F 在目标区间上严格单调(F′=0)、凸性不变,且已知 F(a)F(b)<0。
问:(a) 用牛顿法还是割线法?初值怎么取?(b) 若 θ 的解析表达式不好求导,怎么办?(c) 收敛有多快?
参考答案
(a) 本题满足牛顿法全局收敛的四条:F(a)F(b)<0、F′′ 不变号、F′ 在区间上不为零、初值取在 F(x0)F′′(x0)>0 的一侧(即曲线的“外侧”)。故用牛顿法,且初值不必很靠近根,取在条件 (iv) 指定的一端即可。
(b) 若 F′ 不好算(如 F 是椭圆积分),改用割线法,用差商 xk−xk−1F(xk)−F(xk−1) 代替导数——代价是需要两个初值。再省一步可把 xk−1 钉死为 x0(单点割线法),但收敛阶降为一阶。
(c) 牛顿法二阶(有效数字约翻倍);割线法超线性、阶为黄金比 1.618(解 r2−r−1=0);单点割线法一阶。
S4 · 电网潮流 / 结构分析:迭代法选型
大规模稀疏线性系统 Ax=b,A 对称正定。手头有 Jacobi、Gauss–Seidel、SOR 三种迭代法,朴素高斯消去 O(n3) 不可接受。
问:为何不用朴素高斯消去?如何在三者中选?若用 SOR,ω 怎么定?收敛性怎么判?
参考答案
为何不用:n 很大时 O(n3) 不可行;A 稀疏,直接法会产生大量填充(fill-in),破坏稀疏性。稀疏 / 迭代法是工程主流。
选型:三者都由矩阵分裂 A=N−P 决定(拆法即名字)。A 对称正定或严格对角占优时,Jacobi、G–S 都收敛;SOR 通常更快(介于 ω=1 与最优 ω 之间),但不保证任意 ω 都收敛。
ω:必要条件是 0<ω<2(区间外必发散)。ω=1 退化为 G–S;ω>1 为超松弛。最优 ω 需 A 的谱信息(如 Jacobi 谱半径 μ,ωopt=1+1−μ22)。
收敛判据:根本是 ρ(G)<1(充要);大规模问题算不出特征根,退而用充分条件「找范数使 ∥G∥<1」。注意:∥G∥ 接近 1 时,相邻迭代值差小不能断定接近真解。
| 部分 | 题量 | 分值 | 给分要点 |
| 一 · 概念判断 | 4 | 20 | 判断 2 分 + 理由/反例 3 分;只写判断不给理由得 0 |
| 二 · 填空简答 | 4 | 20 | 结论正确 3 分 + 过程/依据 2 分 |
| 三 · 计算 | 3 | 36 | 每题:方法正确 5 分 + 过程 4 分 + 结果 3 分;结果错但有回代自查意识加 1 分 |
| 四 · 证明 | 2 | 24 | 每题:关键步骤 8 分 + 逻辑完整 4 分;只写结论 ≤2 分 |
| 得分 | 水平 | 下一步 |
| 90–100 | 可以上考场 | 扫一遍必考点冲刺清单,保持手感即可 |
| 75–89 | 过线有把握 | 错题定位回对应小节精读,再去竞技场把薄弱关刷到 100% |
| 60–74 | 有漏洞 | 按错题分块回讲座页「老师强调」一节 + 重做习题集同型题 |
| <60 | 需系统补 | 回最小必读路径按一周路线重走,重点补公式卡自测 |
真实场景题怎么算分
「真实场景特化训练」不计入 100 分,但应用题是拉开差距的地方。评分只看三点:
①方法选对(该用插值还是迭代?牛顿还是割线?);
②说得出为什么(为什么不用直接法、为什么不显式求逆、ω 为什么有限制);
③知道注意事项(病态、必须两个初值、收敛阶的代价)。
只看这里,每条都是老师课堂点名强调、且竞技场与自测卷反复考的:
- 误差:递推稳定性看放大因子(>1 就不稳定);相近数相减改写(a2−b2=(a−b)(a+b));有效数字从首位非零数字数起。
- 范数:∥A∥1 列和最大、∥A∥∞ 行和最大、∥A∥F 不从属;相容性 ∥AB∥≤∥A∥∥B∥;判 I±A 可逆找一个范数 ∥A∥<1 即可。
- 条件数:cond(A)≥1、cond(I)=1;残差小 = 误差小;cond 大 ⇒ 丢有效数字。
- 线性方程组:列主元防小主元放大误差;LU 需顺序主子式非零;三对角用 O(n) 追赶法。
- 迭代法:收敛充要 ρ(G)<1,充分 ∥G∥<1;对角占优 ⇒ Jacobi/G–S 都收敛;SOR 必要 0<ω<2(必要非充分)。
- 特征值与幂法:先单位化再乘 A;特征向量范数随取、特征值必须 2-范数;「严格按模最大」且必为实数;反幂法每步解方程组;位移 A−pI 求最近 p 的特征值。
- 非线性方程:牛顿法局部(初值须近根)vs 四条件给全局(核心:凸性 + 单调性不变);不动点两条件各管不越界 / 收敛;割线法两初值、阶 1.618;单点割线法一阶。
- 插值:n+1 点 → 唯一 n 次多项式;范德蒙德存在唯一但病态;拉格朗日基系数直接读出 f(xi)。
做完这一页,去必考点竞技场把五关刷到全绿,就基本稳了。