闭卷、计时、独立完成——查书或看解答做出来的分数没有意义。做完先按文末评分标准自评,再按漏点定位表回查。范围:第一章 1.1–1.3 与真实应用案例。本卷题目全部为新增(与习题集、考场模拟卷不重叠),覆盖同样的知识点。
时间分配建议:第一部分 20 分钟(每题约 2 分钟,写「判断 + 一句话理由」)→ 第二部分 15 分钟 → 第三部分 35 分钟(每题约 12 分钟,计算题宁慢勿跳步)→ 第四部分 18 分钟(每题约 9 分钟,写清四件套)→ 留最后 2 分钟回头检查第一部分有没有只写对错的题。
阅卷口径(先读,能省不少分)
- 判断必须带机制词:写「错」不给分,写「错,因为有效数字由误差限决定、与打印位数无关」才拿满。常用机制词:末位半单位、误差放大因子、次可乘性、相消、条件数、算子范数。
- 误差限必须写清是谁的限:ε(a)、εr(a)、k(丢失位数)各是哪个量的,混用即扣分。
- 算完要给交代:报结果时同时给出「值 + 误差限 + 有效数字位数」(缺一项扣 1 分);计算题按步给分,只写结论最多得 3 分。
- 开放题四件套:方案 / 理由 / 失败判据 / 验证计划,缺一项扣 2 分(第四部分每题 10 分)。
判断对错,并给出一句话理由或反例。只写「对 / 错」不得分。
1. 用双精度计算,程序打印出的 15 位数字都可靠。
解答
错。 可靠位数由误差限决定,与打印格式无关:eps 只管单步舍入(u≈1.11×10−16),n 步运算可累到约 nu,相近数相减或病态问题放得更多(1.2.5、1.2.8 错误 7)。失分点:只说「浮点有误差」而不落到「位数由误差限判定」。
2. 相对误差限应当用 ε(a)/∣x∣ 计算,因为相对误差的定义式分母是准确值 x。
解答
错。 定义式用 x 只是理论写法;x 未知,实用公式一律用 ∣a∣,两者相差 O(er2)。此外 x=0 处相对误差无定义(定义 1.2.6 后的推导、errata N1)。
3. 若 a 是准确值 x 四舍五入到小数点后 3 位得到的近似值,则 ε(a)=0.0005,且「从末位到第一位非零数字」之间至少有 3 位有效数字。
解答
错(前半对、后半错)。 误差限确实是末位的半个单位(0.0005=21×10−3),但有效数字要从第一位非零数字数到末位:若 a=0.0042,写成 0.42×10−2(m=−2),21×10m−n=21×10−3 给出 n=1——只有 1 位;只有 a 的量级使首位落在小数点前(如 a=1.234)时才是 4 位(定义 1.2.7、例 1.2.3)。
4. 对递推 yn=αn+βyn−1(β 为常数),∣β∣<1 的方向一定数值稳定。
解答
对。 误差满足 δn=βδn−1,故 δn=βnδ0,∣β∣<1 时误差按几何级数收缩,且途中新产生的舍入误差同样被压制(1.2.4 原则 (1))。补充:该判别法只对常系数线性递推成立。
5. 用六个数学上完全等价的公式计算 (2−1)6 会差出五个数量级;改用更高精度(如 float128)就能把这个差距压到机器精度。
解答
错。 差距来自公式结构与相近数相减(第 3 个公式 99−702 是两个 ≈102 的量相减,相对误差被放大到 O(1)),与方法误差/放大因子有关、与字长无关;提高精度只推迟灾难(1.1.3、1.2.8 错误 8)。
6. 由范数等价定理,用 ∥⋅∥1 与 ∥⋅∥2 分别判断同一个迭代序列是否收敛,可能得到相反结论。
解答
错。 等价定理给出 m∥x∥β≤∥x∥α≤M∥x∥β(m,M 与 x 无关),一侧趋于零则另一侧必趋于零——收敛性与范数选择无关;会变的是误差界的大小与停机步数(定理 1.3.3、1.3.6)。
7. ∥A∥F 与向量 2-范数相容,所以它是从属于 ∥⋅∥2 的矩阵范数。
解答
错。 相容 ≠ 从属:任何算子范数都满足 ∥I∥=1,而 ∥I∥F=n=1(n≥2),故 F-范数不从属于任何向量范数(定义 1.3.4 的证明 (c))。
8. 设 ∥A∥=0.5(算子范数),则用 I+A+A2 近似 (I−A)−1 时,误差不超过 0.25。
解答
对。 由 Neumann 级数余项估计 (I−A)−1−∑k=0NAk≤1−∥A∥∥A∥N+1,取 N=2、∥A∥=0.5 得 0.50.53=0.25 ✓(定理 1.3.7)。
9. 对任何方阵 A,都有 ∥A∥2=λmax(A)(A 的最大特征值)。
解答
错。 正确式是 ∥A∥2=λmax(ATA)= 最大奇异值;∥A∥2=λmax(A) 只对对称半正定矩阵成立。反例 A=(0010):λmax(A)=0 而 ∥A∥2=1(1.3.7 错误 3)。
10. 两个相近数相减会大量丢失有效数字,其根本原因是浮点减法本身带有较大的舍入误差。
解答
错。 Sterbenz 引理:y/2≤x≤2y 时 fl(x−y)=x−y 精确。丢掉的数字来自输入本身(中间结果的有效数字已经很少),所以处方是改写(有理化 / 因式分解 / 和差化积),而不是「换一种减法」(引理 1.2.1、1.2.4)。
写清依据与依据来源(给分:结果正确 2 分 + 关键说明 2 分;数值或单位错一处扣 1 分)。
11. 把下列四个情境归入四类误差(模型 / 观测 / 截断 / 舍入),并写出一句话规则说明「数值分析管哪两类」:
(a) 把钢筋混凝土梁简化为 Euler–Bernoulli 一维弹性梁;
(b) 用卫星轨道预报(星历)计算的距离与真实距离之差;
(c) 用 1−2x2+24x4 代替 cosx(x=0.1);
(d) 在 float32 下把 0.1 存成 0.10000000149⋯。
解答
11. (a) 模型误差(建模时丢了信息);(b) 观测误差(输入数据来自测量/预报,带噪声);(c) 截断误差(Taylor 级数被截断,余项约 ∣x∣6/720≈1.4×10−9);(d) 舍入误差(0.1 的二进制展开是循环小数,存不进有限尾数)。一句话规则:数值分析主要研究截断误差与舍入误差——只有这两类能靠算法设计控制;模型与观测误差归建模、测量学科,但进入计算后仍服从 1.2.3 的传播规律(1.2.1 鸟瞰表)。
12. 某近似值 a=0.0123(作为某准确值的近似):
(1) 若 ε(a)=0.0005,问 a 有几位有效数字?
(2) 若仪器给出的误差限其实小得多,ε(a)=0.00005,问 a 有几位有效数字?这个结论说明了什么?
解答
12. (1) 2 位。 a=0.0123=0.123×10−1,m=−1;21×10m−n=21×10−1−n=0.0005=21×10−3 得 −1−n=−3,即 n=2≤k=3 ✓——小数点前那个 0 不算有效数字。
(2) 没有有效数字。 此时要求 21×10−1−n=21×10−5,即 n=4>k=3:误差限要求第 4 位准确,而 a 只写出 3 位数字——写法与误差限不相称,结论是「没有有效数字」(与 例 1.2.3 的 c=0.86×105 同类)。含义:标误差限时必须与写出的位数配套;位数不足时要补写数字或放宽误差限。
13. 两个递推格式:
yn=31yn−1+bn,zn=3zn−1+cn.
(1) 各自的误差每步被乘上什么因子?哪个稳定?
(2) 若初始误差为 10−8,20 步之后两种格式的误差各是多少(保留两位有效数字)?
(3) 对不稳定的那个格式,写出你会在工程上采取的补救办法。
解答
13. (1) 误差满足 δn=βδn−1:y 格式每步乘 β=1/3(稳定,误差按 3−n 收缩);z 格式每步乘 β=3(不稳定,误差按 3n 爆炸)。
(2) 3−20≈2.87×10−10 ⟹ y 的误差 10−8×2.87×10−10≈2.9×10−18(已在双精度之下);320=3.49×109 ⟹ z 的误差 10−8×3.49×109≈3.5×101——荒谬,与真值无关。
(3) 补救:把递推反向走(解出 zn−1 的表达式,使每步乘 1/3),并用积分估计式/边界条件把远端初值先定住(1.2.4 原则 (1)、errata N3)。
14. IEEE 双精度(float64):
(1) 写出 eps 与单位舍入误差 u 的定义式并给出数值;
(2) 解释为什么 1.0 + 2**-53 == 1.0 为真、而 1.0 + 2**-52 > 1.0 为真;
(3) 解释 1e16 + 1 - 1e16 == 0.0 为真;这三个现象共同说明「浮点误差与什么有关」。
解答
14. (1) eps=β1−t=2−52≈2.22×10−16;u=21β1−t=2−53≈1.11×10−16(定义 1.2.10)。
(2) 1 附近的浮点间距是 2−52,相邻两数中点是 1+2−53:增量恰为半个 ulp 时按最近舍入回到 1.0;要跨过一个网格间距至少需要 2−52。
(3) 1016≈253.15 落在 [253,254) 内,该区间间距为 2k−t=253−52=2,增量 1 只是半个间距,被舍回 1016。
共同说明:浮点误差取决于数值的绝对位置(间距 2k−t 随量级成倍变宽),不是「一律 10−16」——大数吃小数是网格太稀的必然(1.2.5 三条后果)。
15. 取 x=(2,−3,6)T:算出 ∥x∥1,∥x∥2,∥x∥∞,逐条核对 1.3.2 的三条显式不等式,并说明「换范数不改变收敛结论」的依据。
解答
15. ∥x∥1=2+3+6=11;∥x∥2=4+9+36=49=7;∥x∥∞=6。
核对(n=3):∥x∥∞=6≤∥x∥2=7≤3×6≈10.392 ✓;∥x∥2=7≤∥x∥1=11≤3×7≈12.124 ✓;∥x∥∞=6≤∥x∥1=11≤3×6=18 ✓。
依据:等价定理(定理 1.3.3)保证任意两种范数被与 x 无关的常数双向夹住,故「是否趋于零」与范数无关;常数 n,n 不可再改进([1,…,1]T 与单分量向量取等)。给分:三个范数 2 分;三条核对 1.5 分;依据 0.5 分。
16.(10 分) 用游标卡尺量得圆柱体的尺寸、用天平称得质量:d=12.50 mm、h=38.00 mm、m=36.60 g,仪器给出的误差限分别为 ε(d)=ε(h)=0.02 mm、ε(m)=0.02 g(各量按四舍五入读到末位)。
(1) 按 ρ=πd2h4m 计算 ρ(取 π=3.14159 为准确值),并给出 εr(ρ) 与 ε(ρ);
(2) 判断 ρ 能有几位有效数字,按规定写法重写结果(带单位:1 mm3=10−3 cm3,请用 g/cm3 报出);
(3) 要把 εr(ρ) 压到 0.3% 以下,最该先改进哪个量?给出量化依据。
第 16 题答案与评分要点
(1) 由多元传播公式(式 1.1),三个量在 ρ 里各出现一次的相对贡献是
εr(ρ)≈εr(m)+∂ρ/∂d 带因子 22εr(d)+εr(h),
(ρ∝md−2h−1:对 d 的相对放大因子是 ∣∂lnρ/∂lnd∣=2。)代入:
εr(m)=36.600.02≈0.0546%,εr(d)=12.500.02=0.16% (×2=0.32%),εr(h)=38.000.02≈0.0526%.
εr(ρ)≈0.0546%+0.32%+0.0526%≈0.427%≈4.3×10−3.
数值:πd2h=3.14159×156.25×38.00=18653.2 mm3,
ρ=18653.194×36.60=7.8485×10−3 g/mm3=7.8485 g/cm3,
ε(ρ)≈4.27×10−3×7.8485≈0.0335 g/cm3.
(2) ρ=7.8485=0.78485×101,m=1。检验 n:21×101−n≥0.0335——n=1 时 0.5≥0.0335 ✓;n=2 时 0.05≥0.0335 ✓;n=3 时 0.005<0.0335 ✗。故 2 位有效数字,结果写作
ρ≈7.8 g/cm3(误差限 0.034 g/cm3).
注意:把 7.8485 原样抄进报告是错误用法——可靠位数由误差限决定(1.2.8 错误 7)。
(3) 三项贡献中 2εr(d)=0.32% 占总误差 0.427% 的约 75%,所以先改进直径 d 的测量:换成误差限 0.004 mm 的千分尺,则
εr(ρ)≈0.0546%+2×12.500.004%+0.0526%=0.0546%+0.064%+0.0526%≈0.171%<0.3% ✓
而只把 ε(m) 减半(0.02→0.01 g)只能把总误差降到 0.40%,达不到目标。结论:改进要挑「相对贡献 × 放大因子」最大的那一项。
评分:(1) 相对误差分解 3 分(三个贡献值各 1 分,d 的因子 2 漏掉扣 1.5)+ 数值 ρ 与 ε(ρ) 3 分;(2) 位数判定 2 分(判定过程 1 + 按规定写法 1);(3) 选对量与量化依据 2 分。未写单位或未给误差限,各扣 0.5 分。
17.(10 分) 设
A=201−3411−23.
(1) 求 ∥A∥1、∥A∥∞、∥A∥F;
(2) 求 ∥A∥2:写出 ATA、说明特征多项式的三个系数怎样不算展开就得到(迹、全部二阶主子式之和、det(ATA)=(detA)2),再给出 λmax 与 ∥A∥2(数值保留四位有效数字);
(3) 核对 ∥A∥2≤∥A∥F 与 ∥A∥2≤∥A∥1∥A∥∞,并说明为什么 ∥A∥1 与 ∥A∥∞ 一般不相等、转置会交换它们。
第 17 题答案与评分要点
(1) 列绝对值和:第 1 列 2+0+1=3,第 2 列 3+4+1=8,第 3 列 1+2+3=6,故
∥A∥1=max(3,8,6)=8.
行绝对值和:第 1 行 6,第 2 行 6,第 3 行 5,故
∥A∥∞=max(6,6,5)=6.
∥A∥F=4+9+1+0+16+4+1+1+9=45≈6.708.
(2) 由各列内积(A 的列为 (2,0,1)T,(−3,4,1)T,(1,−2,3)T):
ATA=5−55−526−85−814.
特征多项式 det(ATA−λI)=−λ3+c2λ2−c1λ+c0 的三个系数不必硬展开:
| 系数 |
来源 |
值 |
| 二次项(迹) |
tr(ATA) |
5+26+14=45 |
| 一次项 |
全部二阶主子式之和 |
(5⋅26−25)+(5⋅14−25)+(26⋅14−64)=105+45+300=450 |
| 常数项 |
det(ATA)=(detA)2 |
detA=30⇒900 |
(detA=2(12+2)+3(0+2)+(0−4)=28+6−4=30。)于是
λ3−45λ2+450λ−900=0 ⟹ λmax≈31.70,
∥A∥2=31.70≈5.630.
(代入核对:31.703−45×31.702+450×31.70−900≈0 ✓。)
(3) 核对:∥A∥2=5.630≤∥A∥F=6.708 ✓;∥A∥1∥A∥∞=48≈6.928≥5.630 ✓,且 ∥A∥2≤∥A∥F 与 ∥A∥2≤∥A∥1∥A∥∞ 对任意矩阵都成立。
为什么两个范数一般不等:∥A∥1 数的是列、∥A∥∞ 数的是行;转置把行列互换,所以 ∥AT∥1=∥A∥∞、∥AT∥∞=∥A∥1(1.3.4 记忆法)。
评分:(1) 三个范数各 1 分;(2) ATA 1.5 分、三个系数各 1 分、λmax 与 ∥A∥2 1.5 分;(3) 两条核对 1.5 分、转置交换关系 1.5 分。把 ∥A∥2 写成 λmax(A)(本题真正的主特征值是 4.4422,且另两个特征值为复数——ρ(A)=4.4422 与 ∥A∥2=5.630 是两个不同的数)本节 0 分。
18.(10 分) 计算积分序列
yn=∫01x+3xndx(n=0,1,2,…).
(1) 推导递推关系,并给出 y0 的表达式与数值;
(2) 判断正向递推 yn=n1−3yn−1 的稳定性:写出误差递推式与 20 步的放大倍数;结合 y0 的四位十进制近似说明后果;
(3) 写出稳定的反向递推格式,并用积分上下界给出 y50 的起点估计区间与取值建议。
第 18 题答案与评分要点
(1) 因 xn+3xn−1=xn−1(x+3),分子恰被分母整除:
yn+3yn−1=∫01x+3xn+3xn−1dx=∫01xn−1dx=n1(n≥1),
y0=∫01x+3dx=ln4−ln3=ln34≈0.28768.
(2) 设 yn−1 带误差 δn−1,由 yn=n1−3yn−1 得
δn=−3δn−1 ⟹ ∣δn∣=3n∣δ0∣,
每步放大 3 倍;320=3.487×109(约 3.5×109 倍)。四位十进制下 y0≈0.2877,其舍入误差约 ∣0.2877−0.287682∣≈1.8×10−5,传到 y20 时放大约 6.3×104——而 y20∼10−2 量级,结果完全失真(教材实测 n=30 时前推值已达 −2.0×1016,真值 5.4×10−3)。结论:∣β∣=3>1,正向递推不稳定,提高字长只推迟灾难(@ex 类比见 1.2.4)。
(3) 反向:由 yn+3yn−1=n1 解出
yn−1=3n1−31yn,β=−31, ∣β∣<1 (每步误差缩小 3 倍).
起点:[0,1] 上 41≤x+31≤31,故 4xn≤x+3xn≤3xn,积分得
4(n+1)1≤yn≤3(n+1)1 ⟹ y50∈[2041,1531]=[0.0049020, 0.0065359].
取区间中点 y50≈0.0057190,起点误差不超过半宽 ≈8.2×10−4;反向递推每步把它压缩 3 倍,传到 y0 时被乘 (1/3)50≈1.4×10−24——起点误差与途中舍入误差都被彻底压掉。
评分:(1) 递推推导 2 分 + y0 1 分;(2) 误差递推式 1.5 分、320 数值 1 分、后果说明 1.5 分;(3) 反向格式 1.5 分、上下界 1 分、中点与结论 0.5 分。只写「反向更稳」而不给误差因子,扣 2 分。
19.(10 分) 某小组用「游标卡尺 + 电子天平」测量一批圆柱体试样的密度(上题第 16 题的公式与仪器),要交一份误差预算报告。要求:
(1) 写出这条测量链上四类误差各自的代表来源(每类至少一条),并说明本预算公式只覆盖了其中哪一部分;
(2) 在「改进仪器」与「改进公式」之间做选择,给出判据;
(3) 写出报告的结果规范(一条结果行里应包含哪些字段);
(4) 给出至少两条可测量的失败判据。
第 19 题参考答案与评分要点
(1) 四类误差的来源清单与预算边界
| 类别 |
本链上的代表来源 |
能否进本预算公式 |
| 模型误差 |
把试样当理想正圆柱(实际有圆度 / 锥度误差)、忽略称量时的空气浮力 |
不进(属于公式本身,V=πd2h/4 已假定是圆柱) |
| 观测误差 |
卡尺读数分辨与对准、天平示值漂移、温湿度导致的尺寸/质量漂移 |
部分进(以 ε(d),ε(h),ε(m) 的形式进) |
| 截断误差 |
本次计算中 π 只取 3.14159(若取近似值) |
可进(当作一个额外输入量,εr(π)≈1.6×10−6,可忽略) |
| 舍入误差 |
商与乘积的浮点执行(双精度,相对 u∼10−16) |
可进(量级远小于 10−3,实际可忽略) |
关键结论:预算公式算出的 0.43% 只覆盖观测误差(+ 可忽略的截断/舍入);模型误差一旦显著(例如锥度 0.1 mm 级),任何减小 ε(d) 的努力都是白费——必须先做模型检查。
(2) 改进仪器还是改进公式:判据
判据写法:「先做残差检验,再看相对贡献占比」。具体:① 对同一试样重复测量 n≥10 次,把样本标准差与仪器误差限对比——若实测离散显著大于仪器限(例如 0.05 mm vs 0.02 mm),说明操作/环境主导,改进仪器无效,应先改测量方法(支架、恒温、多点平均);② 若离散与仪器限一致,则看各项相对贡献:本链中 2εr(d) 占 75%,故换更高精度的直径测量(千分尺)比改公式有效;③ 只有当模型误差成为主导(锥度、椭圆度使 V 的真值与 πd2h/4 偏离超过 0.4%)时,才改公式——例如沿高度测三个截面、用分段圆柱求和。
(3) 结果规范(一行结果必须含四个字段)
ρ=7.8 g/cm3(ε(ρ)=0.03 g/cm3,εr=0.43%,2 位有效数字;依据 ρ=4m/(πd2h),ε(d)=ε(h)=0.02 mm,ε(m)=0.02 g,双精度计算)
四个字段:值 + 绝对(或相对)误差限 + 有效数字位数 + 数据与依据;若引用了范数(向量情形)还必须注明范数,因为它决定误差界的大小。
(4) 失败判据(≥2 条,可测量)
① 重复测量的极差超过仪器误差限的 4 倍(range>4ε)→ 操作或环境主导,预算表失效,先查测量流程;
② 预算预测的相对误差与「高精度参考测量」(如用千分尺 + 精密天平复测同一试样)的实际偏差相差超过一个数量级(实测 ≫ 预测)→ 存在未计入的模型误差;
③ 计算结果写出的位数与 ε(ρ) 不匹配(如报了 7.8485)→ 报告不合格,直接退回;
④ 换了更高精度仪器后总误差没有下降(改进前后 εr 差异 <10%)→ 说明主导项判断错了,回到 ①。
评分:(1) 3 分(四类各 0.5、预算边界 1);(2) 3 分(判据必须可执行:残差检验 + 贡献占比各 1.5);(3) 2 分(四个字段齐全);(4) 2 分(两条各 1 分,须可测量)。只写「提高精度」而不给判据,(2) 得 0 分。
20.(10 分) 一个迭代求解程序用
∥xk+1−xk∥∞≤10−12
作为唯一停机判据。请:
(1) 写出「相邻差」与「真实误差」之间的关系式(设迭代格式为 xk+1=Hxk+d,用 ∥H∥ 表示),并说明什么时候「相邻差小」推不出「误差小」;
(2) 给出三条改进方案,并说明每条治的是什么病;
(3) 用第一章的结论解释:这条判据为什么不能充当误差估计,它与「残差 ≠ 误差」是同一件事的两面吗;
(4) 给出至少两条可测量的失败判据。
第 20 题参考答案与评分要点
(1) 关系式。 记精确解 x∗=(I−H)−1d、误差 ek=xk−x∗,则 ek+1=Hek、ek=Hke0,于是
xk+1−xk=ek+1−ek=(H−I)ek ⟹ ek=(H−I)−1(xk+1−xk).
当 ∥H∥=:q<1 时(定理 1.3.6 保证 I−H 可逆)取范数放缩得
∥ek∥≤1−qq∥xk+1−xk∥以及∥e0∥≤1−q∥x1−x0∥,
∥ek∥≤qk∥e0∥≤1−qqk∥x1−x0∥.
什么时候推不出误差小:放大因子 1−qq 在 q→1 时爆炸(q=0.99 时约 99、q=0.9999 时约 104):相邻差可以小到 10−12,而真实误差仍在 10−8 量级甚至更大;若 ∥H∥≥1 则上式根本给不出界(此时才需要补充篇的谱半径与 ρ<1 的渐近结论)。
(2) 三条改进方案与它治的病
| 方案 |
做法 |
治什么 |
| ① 用后验误差估计替代相邻差 |
先估 q=∥H∥(或每步估比值 ∥Δxk+1∥/∥Δxk∥→q),再用 ∥e∥≲1−qq∥Δx∥ |
治「判据与实际误差差一个放大因子」 |
| ② 双判据:增广残差 + 条件数 |
对方程组类问题同时监控 cond(A)⋅∥r∥/∥b∥(1.3.6) |
治「残差/相邻差都小、误差却大」的病态问题 |
| ③ 相对化 + 注明范数 |
改判据为 ∥xk+1∥+1∥xk+1−xk∥≤τ,并在报告里写明用的是哪个范数 |
治「∥⋅∥∞ 只盯最坏分量、尺度大时绝对阈值失效」 |
(补充:若 q 接近 1、步数不可预测,正确做法是换更好的预条件 / 换算法,而不是继续调阈值。)
(3) 为什么它不能当误差估计。 停机判据度量的是算法是否稳定下来(增量信息 / 后向信息),误差度量的是解离真值多远(前向信息);两者的桥梁恰好是问题的条件数或迭代矩阵的 q/(1−q)。这与 1.2.6 的「conditioning 与 stability 必须分开谈」、1.2.6 的「求解器只能保证残差、用户要的是误差」同一句话;矩阵版 e=A−1r 的完整推导在补充篇 2.4.4。是同一件事的两面:残差对应 Ax=b 的「相邻/后向」信息,相邻差对应迭代格式的增量信息,二者都要乘上放大因子才能变成误差。
(4) 失败判据(≥2 条,可测量)
① 停机后用一个高精度参考解(或第二种算法/更高精度重算)对照,实测相对误差 >10−6——说明 q/(1−q) 被低估,判据过松;
② 停机前 ∥Δxk+1∥/∥Δxk∥ 持续大于 0.9(衰减比不趋于常数 q<1)——迭代矩阵谱半径过于接近 1,步数不可预测,应换算法或预条件;
③ 相邻差已低于阈值,但残差(或 cond⋅∥r∥/∥b∥)仍在下降——说明是「增量停滞」而不是「误差已小」,必须继续迭代或改判据;
④ 同一问题换用 ∥⋅∥2 后停机步数变化超过 10 倍——说明阈值 τ 与范数不匹配,报告里没有注明范数。
评分:(1) 4 分(误差递推 1、关系式 1.5、放大因子与失效条件 1.5);(2) 3 分(三条各 1 分,须写出「治什么病」);(3) 2 分(区分前向/后向 + 指出桥梁是条件数);(4) 1 分(两条可测量判据)。
| 部分 | 题量 | 分值 | 给分要点 |
| 一 · 概念判断 | 10 | 30 | 判断 1 分 + 理由 / 反例 2 分;只写对错得 1 分 |
| 二 · 填空与简答 | 5 | 20 | 结果正确 2 分 + 关键说明 2 分;数值 / 单位错一处扣 1 分 |
| 三 · 推导计算 | 3 | 30 | 首行或公式 3 分 + 链条完整 4 分 + 位数 / 边界等收口 3 分;只写结论最多 3 分 |
| 四 · 开放论述 | 2 | 20 | 方案具体 4 分 + 判据可测量 4 分 + 验证计划 2 分 |
段位换算:≥90 优秀(能给同学讲清机制);75−89 良好(能独立做误差分析);60−74 及格(公式记得住、条件记不住);<60 回到对应节的「即时自测」与习题集重做。自评纪律:① 判断理由必须包含机制词(末位半单位、放大因子、次可乘、相消、条件数、算子范数),否则按半对给 1 分;② 计算题凡是没给误差限或没给位数的,从收口里扣 1 分;③ 开放题只写方案、不写失败判据的,第四部分上限 6 分。
最后一句:本卷的目的不是打分,而是把「以为会了」的地方暴露出来。凡「能算出数但说不清依据」的题,回到对应节的 intuition 框与回忆卡当天再练一遍;同一类错误出现两次以上,那不是粗心,而是概念缺口。