用法:遮住答案,在空白纸上写;写完再展开对照。每张卡末尾标记 ✓ 能完整写出 / △ 有关键词但缺条件 / × 完全想不起,下一轮只重做 △ 与 ×。全页 40 张卡,一轮约 30 分钟。本页负责检索(从空白纸里把知识叫出来);压缩在公式与结论卡,辨析在易混对照卡,考场节奏在综合自测卷;补充篇的对应卡片见补充篇回忆卡。
A1. 四类误差各由谁造成?逐条写出定义(模型误差 / 观测误差 / 截断误差 / 舍入误差)。
解答
模型误差=数学模型与实际问题的差别(写公式的人);观测误差=观测数据与实际数据之差(测数据的人);截断误差(方法误差)=用近似方法代替精确求解所得的近似解与模型精确解之差(选算法的人);舍入误差(计算误差)=字长有限导致数据与运算结果存放产生的误差(机器)。只有后两类能由算法设计控制(定义 1.2.1–定义 1.2.4、1.2.1)。(✓ / △ / ×)
A2. 写出绝对误差 e、绝对误差限 ε 的定义,以及 x,a,ε 三者的关系式。
解答
e=x−a;若 ∣e∣≤ε 则称 ε 为误差限;关系写作 a−ε<x<a+ε,即 x=a±ε——a 已知、x 未知(定义 1.2.5)。(✓ / △ / ×)
A3. 相对误差的真定义式与实用公式分别是什么?x=0 时怎么办?
解答
真定义 er=xx−a;实用公式 er≈ax−a(分母换成已知的 a),二者相差 O(er2)。相对误差限 εr=ε(a)/∣a∣;x=0 处相对误差无定义,近似 0 只能用绝对误差谈论(定义 1.2.6、errata N1)。(✓ / △ / ×)
A4. 「a 有 n 位有效数字」的判据是什么?尾零为什么不能省?
解答
先把 a 写成规格化形式 a=±0.a1a2⋯ak×10m(a1=0):若 ε(a)=21×10m−n 且 n≤k,就说 a 有 n 位有效数字——从第一位非零数字数到「敢担保」的那一位。尾零是「这几位也查过了」的书面声明:2.1400 是 5 位、2.14 是 3 位(定义 1.2.7)。(✓ / △ / ×)
A5. 什么叫「算法」?衡量算法优劣的两个标准是什么?计算复杂性衡量哪两样东西?
解答
算法=从输入数据算出数值问题解的有限的基本运算序列(四则、逻辑、基本函数)。两个标准:① 有可靠的理论基础(正确性、收敛性、数值稳定性、可做误差分析);② 有良好的计算复杂性。复杂性=达到给定精度所需的时间(计算量)与空间(内存)(定义 1.2.8)。(✓ / △ / ×)
A6. 「数值稳定」的算法怎么定义?它与「问题良态 / 病态」是一回事吗?
解答
数值稳定=能控制舍入误差的传播:误差在逐步执行中不放大(放大因子有界)。它与病态不是一回事:病态是问题本身对输入扰动的敏感(条件数大,换算法无效);不稳定是算法的属性(换算法可救)(定义 1.2.9、1.2.6 交叉表)。(✓ / △ / ×)
A7. 写出浮点数集 F(β,t,L,U) 中数的形式,以及机器精度 eps 与单位舍入误差 u。
解答
x=±m×βe,m=(0.b1b2⋯bt)β,L≤e≤U;eps=β1−t、u=21β1−t,且 xfl(x)−x≤u。双精度(β=2,t=53):eps=2−52≈2.22×10−16、u≈1.11×10−16(定义 1.2.10)。(✓ / △ / ×)
A8. 向前误差与向后误差各问什么问题?什么叫向后稳定?
解答
向前误差是「结果离真值多远」(f(x)−y~);向后误差是「使精确模型恰好算出 y~ 的最小输入扰动」——结果是谁家的精确答案。向后误差与机器舍入同量级时称算法向后稳定,此时向前误差≈条件数 × 向后误差(定义 1.2.11、1.2.6)。(✓ / △ / ×)
A9. 向量范数的三条公理、矩阵范数的四条公理各是什么?「相容」又指什么?
解答
向量:正定、齐次、三角不等式。矩阵:前三条 + 次可乘性 ∥AB∥≤∥A∥∥B∥。相容指 ∥Ax∥≤∥A∥∥x∥ 对一切 A 与 x 成立——矩阵范数与向量范数要配对使用(定义 1.3.1、定义 1.3.2、定义 1.3.3)。(✓ / △ / ×)
A10. 什么叫从属范数(算子范数)?为什么 Frobenius 范数不是算子范数?
解答
从属范数=最大拉伸倍数 ∥A∥=max∥x∥=1∥Ax∥。任何算子范数都满足 ∥I∥=1,而 n≥2 时 ∥I∥F=n=1,故 ∥⋅∥F 不从属于任何向量范数——它与 ∥⋅∥2 相容,但相容不等于从属(定理 1.3.4、定义 1.3.4)。(✓ / △ / ×)
A 组的适用条件(背定义时要一起带走的东西)
- 误差限 ε 是上界不是误差本身;x=a±ε 里只有 a 是已知的。
- 相对误差两种写法只在 ∣er∣≪1 时等价;∣er∣ 到 0.1 量级时近似值已接近失效。
- 有效数字判据里的 ε 必须与该数写出的位数相称:只写 2 位数字却要求第 7 位准确,结论只能是「没有有效数字」(例 1.2.3 的 c=0.86×105)。
- eps 与 u 差一个 2:eps 是相邻浮点数的相对间距,u 是一次舍入的相对误差上界,写公式时不要互换。
- 从属范数只对配对的向量范数有意义:∥A∥p 里的 p 说的是向量的 p-范数。
B1. 四类误差里哪两类是数值分析的研究对象?另外两类混进计算后还服从本章规律吗?
解答
研究对象是截断误差与舍入误差(唯一能靠算法设计控制的两类);模型误差与观测误差分别归建模、测量学科,但一旦进入计算,仍按 1.2.3 的传播规律起作用,所以本章工具对它们一视同仁(1.2.1 鸟瞰表)。(✓ / △ / ×)
B2. 为什么由四舍五入得到的近似值,误差限恰是末位的半个单位?
解答
四舍五入到小数点后 d 位时,被舍掉的部分 δ∈[−21×10−d, 21×10−d],故 ∣x−a∣≤21×10−d;而 a 的末位正落在 10−d 位上——于是「从末位一路数到第一位非零数字」全是有效数字(例 1.2.2 后的推导)。(✓ / △ / ×)
B3. 有效数字与相对误差限怎么互相换算?「每多一位有效数字」意味着什么?
解答
n 位有效数字 ⟹ εr(a)≤2a11×101−n≤21×101−n;反之若 εr(a)≤2(a1+1)1×101−n,则至少有 n 位。每多一位有效数字,相对误差缩小十倍(定理 1.2.1)。(✓ / △ / ×)
B4. 函数求值的误差为什么乘一阶导数?一元、多元各写成什么?
解答
由 Taylor 一阶展开:e(u~)≈f′(a)e(a)、ε(u~)≈∣f′(a)∣ε(a);多元则一阶项逐项求和(式 1.1)。前提是 ∣f′′/f′∣ 不很大,否则一阶项消失、须用高阶项(定理 1.2.2)。(✓ / △ / ×)
B5. 四则运算的误差限公式有哪四条?最危险的一步是什么、危险来自哪里?
解答
ε(a±b)=ε(a)+ε(b);乘除只需相对误差限相加;a−b 的相对误差限分母是 ∣a−b∣。危险来自相减(分母趋零时相对误差爆炸,例 1.2.7 一步吃光 4 位有效数字)与小除数(∣b∣2 二次放大)。(✓ / △ / ×)
B6. 递推 yn=αn+βyn−1 的误差如何随步数变化?怎样据此选方向?
解答
误差满足 δn=βδn−1,故 δn=βnδ0:∣β∣<1 收缩(稳定),∣β∣>1 爆炸(不稳定)。遇到 ∣β∣>1 的方向就把递推倒过来走,并用估计式定住远端初值(1.2.4 原则 (1)、errata N3)。(✓ / △ / ×)
B7. Sterbenz 引理说什么?由此看,「避免相近数相减」到底在避免什么?
解答
引理:y/2≤x≤2y 时 fl(x−y)=x−y 精确成立——浮点减法本身不引入新误差。所以丢掉的数字在减法之前就随输入消失了:要避免的是「把有效数字所剩无几的中间结果送去相减」,办法是按处方改写(有理化 / 因式分解 / 和差化积)(引理 1.2.1、1.2.4 处方表)。(✓ / △ / ×)
B8. 浮点算术基本公理是什么?它给出「n 次运算」的什么定量结论?
解答
公理:fl(x op y)=(x op y)(1+δ),∣δ∣≤u——每一步都等于精确结果乘一个微因子。配合相对误差相加得 n 次运算的相对误差约 ≲nu,即「双精度算一万步仍可靠、算一亿步要留神」(1.2.5 公理框)。(✓ / △ / ×)
B9. 病态问题与不稳定算法交叉起来有四种情形,各该怎么办?
解答
良态+稳定 → 万事大吉;病态+稳定 → 向前误差仍大,只能重述问题;良态+不稳定 → 换算法;两者都坏 → 重述+换算法。一句话:病态归问题、不稳定归算法,二者独立(1.2.6 交叉表)。(✓ / △ / ×)
B10. 三种常用向量范数之间有哪些显式不等式?为什么说「换范数不改变收敛结论」?
解答
∥x∥∞≤∥x∥2≤∥x∥1≤n∥x∥∞,细写成 ∥x∥∞≤∥x∥2≤n∥x∥∞、∥x∥2≤∥x∥1≤n∥x∥2;常数 n,n 不可再改进(取等例:全 1 向量或单分量向量)。等价定理保证「是否趋于零」与范数无关,会变的是误差界大小与停机步数(定理 1.3.3)。(✓ / △ / ×)
B11. ∥A∥<1 能换来哪三件东西?
解答
① I±A 非奇异;② ∥(I±A)−1∥≤(1−∥A∥)−1(该式要求算子范数,因为用了 ∥I∥=1);③ Neumann 级数 ∑k≥0Ak=(I−A)−1 收敛,且部分和有余项估计(定理 1.3.6、定理 1.3.7)。(✓ / △ / ×)
B 组的适用条件与常见越界
- 「相减危险」说的是相对误差限;绝对误差限对加减一视同仁,不要写成「减法比加法误差大」。
- 传播公式是一阶近似:∣f′′/f′∣ 不小时必须改用高阶项或直接做区间估计。
- 递推方向判别法只对常系数线性递推成立(β 为常数);变系数时按每步乘积逐段估计。
- 等价定理里的 m,M 依赖 n 与范数对,不是万能常数;报误差必须注明范数。
- ∥A∥<1 只是充分条件:I−A 可逆时 ∥A∥ 完全可以大于 1(如 A=10I)。
C1. 误差三件套与有效数字判据。
解答
e=x−a;∣e∣≤ε;x=a±ε;εr=ε(a)/∣a∣;a=±0.a1⋯ak×10m,ε(a)=21×10m−n(n≤k)⟺ n 位有效数字(定义 1.2.5–定义 1.2.7)。(✓ / △ / ×)
C2. 由误差限反求有效数字位数的操作式(含「没有有效数字」的情形)。
解答
解 21×10m−n=ε(a),即 n=m−log10(2ε);n≤k 成立,n>k 则没有有效数字。快捷核对:误差限是不是某一位的半个单位(1.2.2 判定流程框)。(✓ / △ / ×)
C3. 有效数字 ⟹ 相对误差限(两个档次)。
解答
εr(a)≤2a11×101−n≤21×101−n(定理 1.2.1)。(✓ / △ / ×)
C4. 式 (1.1):多元函数求值的误差估计。
解答
ε(u~)≈∑i=1n∂xi∂f(a1,…,an)ε(ai);带符号版 e(u~)≈∑i∂xi∂fe(ai)(1.2.3)。(✓ / △ / ×)
C5. 四则运算误差限表(绝对 + 相对)。
解答
ε(a±b)=ε(a)+ε(b);ε(ab)≈∣a∣ε(b)+∣b∣ε(a);ε(a/b)≈∣b∣2∣a∣ε(b)+∣b∣ε(a);乘除的相对误差限相加;εr(a±b)=∣a±b∣ε(a)+ε(b)。(✓ / △ / ×)
C6. 递推误差与放大因子。
解答
δn=βnδ0;同一个递推正反两方向的稳定性可以完全相反(β=5:前推放大 5n、反推缩小 5−n)(1.2.4)。(✓ / △ / ×)
C7. 浮点的两组数字:eps / u / 间距。
解答
eps=β1−t、u=21β1−t;double:t=53,eps≈2.22×10−16、u≈1.11×10−16;float32:t=24,eps≈1.19×10−7。区间 [2k,2k+1) 内相邻浮点数间距恒为 2k−t(定义 1.2.10、1.2.5)。(✓ / △ / ×)
C8. 浮点公理与 n 步运算的误差量级。
解答
fl(x op y)=(x op y)(1+δ),∣δ∣≤u;n 步基本运算的相对误差 ≲nu(1.2.5)。(✓ / △ / ×)
C9. 函数条件数与「丢几位」的估计式。
解答
cond(f,x)=f(x)xf′(x);相对误差满足 εr(u)≲cond(f,x)⋅εr(a);当输入相对误差为 u 时,丢失位数 k≈log10cond(定义 1.2.11、1.3.6)。(✓ / △ / ×)
C10. 五个矩阵范数公式与两条收口关系。
解答
∥A∥1=maxj∑i∣aij∣;∥A∥∞=maxi∑j∣aij∣;∥A∥2=λmax(ATA);∥A∥F=(∑i,jaij2)1/2。相容性 ∥Ax∥p≤∥A∥p∥x∥p;等价性 ∥x∥∞≤∥x∥2≤∥x∥1≤n∥x∥∞。(✓ / △ / ×)
C11. Neumann 界与级数(含余项)。
解答
∥A∥<1 ⟹ ∥(I±A)−1∥≤1−∥A∥1,∑k=0∞Ak=(I−A)−1,(I−A)−1−∑k=0NAk≤1−∥A∥∥A∥N+1(余项估计要求算子范数)(定理 1.3.6、定理 1.3.7)。(✓ / △ / ×)
C 组的适用条件
- 误差限公式都是最坏情形:给的是界,不是期望值。
- ∣f′(a)∣ 是放大因子、不是误差;∣f′∣≤1 时函数反而吸收误差。
- 除法公式要求 b=0,且只有 ∣b∣≪∣a∣ 时才危险(分母 ∣b∣2 二次放大)。
- ∥A∥2=λmax(ATA) 的流程是先转置相乘、再取最大特征值、最后开方;∥A∥2=λmax(A) 只对对称半正定矩阵成立。
- k≈log10cond 是量级估计,实测误差常比它小几倍;设计系统时按界取,评估结果时按实测报。
D1. 「a=3.1416 是 π 的 4 位有效数字近似——它保留了 4 位小数。」
解答
两处都错。有效数字位数 ≠ 小数位数:a=0.31416×101,实际误差 ≈7.3×10−6;n=5 时 21×10−4=5×10−5≥7.3×10−6、n=6 时 5×10−6<7.3×10−6,故是 5 位(例 1.2.4)。(✓ / △ / ×)
D2. 「相对误差限应当用 ε(a)/∣x∣ 计算,因为相对误差的定义分母是准确值 x。」
解答
错。定义式用 x 是理论写法,但 x 未知,实用公式一律用 ∣a∣;两者相差 O(er2),只在 ∣er∣ 达 0.1 量级时才需警惕(定义 1.2.6 后的推导、errata N1)。(✓ / △ / ×)
D3. 「ε(ab)≈∣a∣ε(b)+∣b∣ε(a) 是恒等式,可以拿来做严格证明。」
解答
错。它是 Taylor 一阶近似,余项为 O(ε2);精确的带符号关系是 e(ab)=be(a)+ae(b)+e(a)e(b),相对形式含乘积项 er(a)er(b)。只有加减法的带符号版本才是精确的(1.2.3 思考框)。(✓ / △ / ×)
D4. 「相近数相减之所以丢有效数字,是因为浮点减法本身的舍入误差很大。」
解答
错。Sterbenz 引理保证相近数相减在浮点内精确(1.870−1.865 算出 0.005 完全准确);丢掉的数字来自输入本身只剩 4 位有效数字。正确结论:别把有效数字快用完的中间结果送去相减(引理 1.2.1、1.2.4)。(✓ / △ / ×)
D5. 「递推一律正向计算;稳定性只与步长和精度有关,与递推方向无关。」
解答
错。稳定性由递推系数与方向绑定:δn=βnδ0,∣β∣>1 的方向必然爆炸(5100≈7.9×1069),提高精度只推迟灾难、不消除它。应反向递推并用估计式(如 6(n+1)1<yn<5(n+1)1)定住远端初值(1.2.4 原则 (1))。(✓ / △ / ×)
D6. 「方差的两遍法公式 n1∑xi2−xˉ2 与中心化算法只是写法不同,精度一样。」
解答
错。数据共同偏移大时 ∑xi2 与 nxˉ2 是相近数相减,实测可输出 0.0(真值 5.25):∑xi2≈8×1018 的绝对舍入误差已达个位数量级。应先平移或改用 Welford 单遍法(1.2.7 实验二、1.2.8 错误 6)。(✓ / △ / ×)
D7. 「打印出 17 位数字说明这 17 位都可信;判断两个浮点数相等直接写 == 即可。」
解答
错。可靠位数由误差限决定,与打印位数无关;十进制的「整」在二进制里往往不整(0.1+0.2=0.30000000000000004)。判等应写 ∣x−y∣≤εmax(∣x∣,∣y∣) 或用 math.isclose(1.2.5 三条后果、1.2.8 错误 7)。(✓ / △ / ×)
D8. 「∥AB∥=∥A∥∥B∥;F-范数是从属范数;∥A∥2 就是 λmax(A)。」
解答
三处都错。① 公理只给次可乘 ∥AB∥≤∥A∥∥B∥,一般严格小于;② ∥I∥F=n=1,故 F-范数不从属(相容但非算子范数);③ 正确式是 λmax(ATA),只有对称半正定矩阵才有 ∥A∥2=λmax(A)——反例 A=(0010):λmax(A)=0 而 ∥A∥2=1(1.3.7 错误 1–3)。(✓ / △ / ×)
D 组的典型误用清单(考前扫一眼)
① 把有效数字当小数位数;② 相对误差限分母写成 x;③ 把一阶传播公式当恒等式;④ 相近数相减还指望保住有效数字;⑤ 递推不看系数与方向;⑥ 两遍法方差直接上机;⑦ 浮点用 ==、相信打印位数;⑧ 以为升级精度就能救方法性误差;⑨ 把次可乘当等式、把 F-范数当算子范数、把 ∥A∥2 写成 λmax(A);⑩ 以为换范数会改变收敛结论。
- 建卡组:打开 Anki → 新建卡组「数值分析-第 1 章」(与「数值分析-全书」并列,不要混进默认卡组)。
- 逐条录入:正面放「编号 + 题干」,背面放「答案 + 依据」(带上 [@编号] 与节号,便于回正文核对);纯文本批量导入用
Tab 分隔「正面 / 背面」,一行一张卡——本页 40 张按 A1–A10 / B1–B11 / C1–C11 / D1–D8 的顺序录。
- 复习节点(1 / 3 / 7 天):当天 → 隔天(第 1 天)→ 第 3 天 → 第 7 天共四轮;每轮只重做标 △ 与 × 的卡;出现「能背出但说不出适用条件」的情况一律按 △ 处理。
- 加一道保险:C 组(公式卡)背面附一张手写照片;四轮后仍为 × 的卡,抄进勘误与补注第四节的待核实条目单独跟踪。分工记住一句话:公式与结论卡负责压缩、易混对照卡负责辨析、综合自测卷负责考场节奏,本页负责检索。