Keyboard shortcuts

Press ← or → to navigate between chapters

Press S or / to search in the book

Press ? to show this help

Press Esc to hide this help

综合自测卷 · 第一章(100 分 · 建议 90 分钟)

闭卷、计时、独立完成——查书或看解答做出来的分数没有意义。做完先按文末评分标准自评,再按漏点定位表回查。范围:第一章 1.1–1.3 与真实应用案例。本卷题目全部为新增(与习题集、考场模拟卷不重叠),覆盖同样的知识点。

时间分配建议:第一部分 2020 分钟(每题约 22 分钟,写「判断 + 一句话理由」)→ 第二部分 1515 分钟 → 第三部分 3535 分钟(每题约 1212 分钟,计算题宁慢勿跳步)→ 第四部分 1818 分钟(每题约 99 分钟,写清四件套)→ 留最后 22 分钟回头检查第一部分有没有只写对错的题。

阅卷口径(先读,能省不少分)
  • 判断必须带机制词:写「错」不给分,写「错,因为有效数字由误差限决定、与打印位数无关」才拿满。常用机制词:末位半单位、误差放大因子、次可乘性、相消、条件数、算子范数。
  • 误差限必须写清是谁的限:ε(a)\varepsilon(a)、εr(a)\varepsilon_{\mathrm r}(a)、kk(丢失位数)各是哪个量的,混用即扣分。
  • 算完要给交代:报结果时同时给出「值 + 误差限 + 有效数字位数」(缺一项扣 11 分);计算题按步给分,只写结论最多得 33 分。
  • 开放题四件套:方案 / 理由 / 失败判据 / 验证计划,缺一项扣 22 分(第四部分每题 1010 分)。

第一部分 · 概念判断(10 题 × 3 分 = 30 分)

判断对错,并给出一句话理由或反例。只写「对 / 错」不得分。

1. 用双精度计算,程序打印出的 1515 位数字都可靠。

解答

错。 可靠位数由误差限决定,与打印格式无关:eps\mathrm{eps} 只管单步舍入(u≈1.11×10−16u\approx1.11\times10^{-16}),nn 步运算可累到约 nunu,相近数相减或病态问题放得更多(1.2.5、1.2.8 错误 7)。失分点:只说「浮点有误差」而不落到「位数由误差限判定」。

2. 相对误差限应当用 ε(a)/∣x∣\varepsilon(a)/\lvert x\rvert 计算,因为相对误差的定义式分母是准确值 xx。

解答

错。 定义式用 xx 只是理论写法;xx 未知,实用公式一律用 ∣a∣\lvert a\rvert,两者相差 O(er2)O(e_{\mathrm r}^2)。此外 x=0x=0 处相对误差无定义(定义 1.2.6 后的推导、errata N1)。

3. 若 aa 是准确值 xx 四舍五入到小数点后 33 位得到的近似值,则 ε(a)=0.0005\varepsilon(a)=0.0005,且「从末位到第一位非零数字」之间至少有 33 位有效数字。

解答

错(前半对、后半错)。 误差限确实是末位的半个单位(0.0005=12×10−30.0005=\tfrac12\times10^{-3}),但有效数字要从第一位非零数字数到末位:若 a=0.0042a=0.0042,写成 0.42×10−20.42\times10^{-2}(m=−2m=-2),12×10m−n=12×10−3\tfrac12\times10^{m-n}=\tfrac12\times10^{-3} 给出 n=1n=1——只有 11 位;只有 aa 的量级使首位落在小数点前(如 a=1.234a=1.234)时才是 44 位(定义 1.2.7、例 1.2.3)。

4. 对递推 yn=αn+βyn−1y_n=\alpha_n+\beta y_{n-1}(β\beta 为常数),∣β∣<1\lvert\beta\rvert\lt1 的方向一定数值稳定。

解答

对。 误差满足 δn=β δn−1\delta_n=\beta\,\delta_{n-1},故 δn=βnδ0\delta_n=\beta^n\delta_0,∣β∣<1\lvert\beta\rvert\lt1 时误差按几何级数收缩,且途中新产生的舍入误差同样被压制(1.2.4 原则 (1))。补充:该判别法只对常系数线性递推成立。

5. 用六个数学上完全等价的公式计算 (2−1)6(\sqrt2-1)^6 会差出五个数量级;改用更高精度(如 float128)就能把这个差距压到机器精度。

解答

错。 差距来自公式结构与相近数相减(第 3 个公式 99−70299-70\sqrt2 是两个 ≈102\approx10^2 的量相减,相对误差被放大到 O(1)O(1)),与方法误差/放大因子有关、与字长无关;提高精度只推迟灾难(1.1.3、1.2.8 错误 8)。

6. 由范数等价定理,用 ∥⋅∥1\lVert\cdot\rVert_1 与 ∥⋅∥2\lVert\cdot\rVert_2 分别判断同一个迭代序列是否收敛,可能得到相反结论。

解答

错。 等价定理给出 m∥x∥β≤∥x∥α≤M∥x∥βm\lVert\boldsymbol x\rVert_\beta\le\lVert\boldsymbol x\rVert_\alpha\le M\lVert\boldsymbol x\rVert_\beta(m,Mm,M 与 x\boldsymbol x 无关),一侧趋于零则另一侧必趋于零——收敛性与范数选择无关;会变的是误差界的大小与停机步数(定理 1.3.3、1.3.6)。

7. ∥A∥F\lVert A\rVert_F 与向量 2-范数相容,所以它是从属于 ∥⋅∥2\lVert\cdot\rVert_2 的矩阵范数。

解答

错。 相容 ≠ 从属:任何算子范数都满足 ∥I∥=1\lVert I\rVert=1,而 ∥I∥F=n≠1\lVert I\rVert_F=\sqrt n\neq1(n≥2n\ge2),故 F-范数不从属于任何向量范数(定义 1.3.4 的证明 (c))。

8. 设 ∥A∥=0.5\lVert A\rVert=0.5(算子范数),则用 I+A+A2I+A+A^2 近似 (I−A)−1(I-A)^{-1} 时,误差不超过 0.250.25。

解答

对。 由 Neumann 级数余项估计 ∥(I−A)−1−∑k=0NAk∥≤∥A∥N+11−∥A∥\left\lVert(I-A)^{-1}-\sum_{k=0}^{N}A^k\right\rVert\le\dfrac{\lVert A\rVert^{N+1}}{1-\lVert A\rVert},取 N=2N=2、∥A∥=0.5\lVert A\rVert=0.5 得 0.530.5=0.25\dfrac{0.5^3}{0.5}=0.25 ✓(定理 1.3.7)。

9. 对任何方阵 AA,都有 ∥A∥2=λmax⁡(A)\lVert A\rVert_2=\lambda_{\max}(A)(AA 的最大特征值)。

解答

错。 正确式是 ∥A∥2=λmax⁡(ATA)=\lVert A\rVert_2=\sqrt{\lambda_{\max}(A^{\mathrm T}A)}= 最大奇异值;∥A∥2=λmax⁡(A)\lVert A\rVert_2=\lambda_{\max}(A) 只对对称半正定矩阵成立。反例 A=(0100)A=\begin{pmatrix}0&amp;1\\0&amp;0\end{pmatrix}:λmax⁡(A)=0\lambda_{\max}(A)=0 而 ∥A∥2=1\lVert A\rVert_2=1(1.3.7 错误 3)。

10. 两个相近数相减会大量丢失有效数字,其根本原因是浮点减法本身带有较大的舍入误差。

解答

错。 Sterbenz 引理:y/2≤x≤2yy/2\le x\le2y 时 fl⁡(x−y)=x−y\operatorname{fl}(x-y)=x-y 精确。丢掉的数字来自输入本身(中间结果的有效数字已经很少),所以处方是改写(有理化 / 因式分解 / 和差化积),而不是「换一种减法」(引理 1.2.1、1.2.4)。

第二部分 · 填空与简答(5 题 × 4 分 = 20 分)

写清依据与依据来源(给分:结果正确 22 分 + 关键说明 22 分;数值或单位错一处扣 11 分)。

11. 把下列四个情境归入四类误差(模型 / 观测 / 截断 / 舍入),并写出一句话规则说明「数值分析管哪两类」: (a) 把钢筋混凝土梁简化为 Euler–Bernoulli 一维弹性梁; (b) 用卫星轨道预报(星历)计算的距离与真实距离之差; (c) 用 1−x22+x4241-\dfrac{x^2}{2}+\dfrac{x^4}{24} 代替 cos⁡x\cos x(x=0.1x=0.1); (d) 在 float32 下把 0.10.1 存成 0.10000000149⋯0.10000000149\cdots。

解答

11. (a) 模型误差(建模时丢了信息);(b) 观测误差(输入数据来自测量/预报,带噪声);(c) 截断误差(Taylor 级数被截断,余项约 ∣x∣6/720≈1.4×10−9\lvert x\rvert^6/720\approx1.4\times10^{-9});(d) 舍入误差(0.10.1 的二进制展开是循环小数,存不进有限尾数)。一句话规则:数值分析主要研究截断误差与舍入误差——只有这两类能靠算法设计控制;模型与观测误差归建模、测量学科,但进入计算后仍服从 1.2.3 的传播规律(1.2.1 鸟瞰表)。

12. 某近似值 a=0.0123a=0.0123(作为某准确值的近似): (1) 若 ε(a)=0.0005\varepsilon(a)=0.0005,问 aa 有几位有效数字? (2) 若仪器给出的误差限其实小得多,ε(a)=0.00005\varepsilon(a)=0.00005,问 aa 有几位有效数字?这个结论说明了什么?

解答

12. (1) 2 位。 a=0.0123=0.123×10−1a=0.0123=0.123\times10^{-1},m=−1m=-1;12×10m−n=12×10−1−n=0.0005=12×10−3\tfrac12\times10^{m-n}=\tfrac12\times10^{-1-n}=0.0005=\tfrac12\times10^{-3} 得 −1−n=−3-1-n=-3,即 n=2≤k=3n=2\le k=3 ✓——小数点前那个 00 不算有效数字。 (2) 没有有效数字。 此时要求 12×10−1−n=12×10−5\tfrac12\times10^{-1-n}=\tfrac12\times10^{-5},即 n=4>k=3n=4\gt k=3:误差限要求第 44 位准确,而 aa 只写出 33 位数字——写法与误差限不相称,结论是「没有有效数字」(与 例 1.2.3 的 c=0.86×105c=0.86\times10^5 同类)。含义:标误差限时必须与写出的位数配套;位数不足时要补写数字或放宽误差限。

13. 两个递推格式: yn=13yn−1+bn,zn=3zn−1+cn.y_n=\frac13y_{n-1}+b_n,\qquad z_n=3z_{n-1}+c_n . (1) 各自的误差每步被乘上什么因子?哪个稳定? (2) 若初始误差为 10−810^{-8},2020 步之后两种格式的误差各是多少(保留两位有效数字)? (3) 对不稳定的那个格式,写出你会在工程上采取的补救办法。

解答

13. (1) 误差满足 δn=βδn−1\delta_n=\beta\delta_{n-1}:yy 格式每步乘 β=1/3\beta=1/3(稳定,误差按 3−n3^{-n} 收缩);zz 格式每步乘 β=3\beta=3(不稳定,误差按 3n3^{n} 爆炸)。 (2) 3−20≈2.87×10−103^{-20}\approx2.87\times10^{-10} ⟹ yy 的误差 10−8×2.87×10−10≈2.9×10−1810^{-8}\times2.87\times10^{-10}\approx2.9\times10^{-18}(已在双精度之下);320=3.49×1093^{20}=3.49\times10^{9} ⟹ zz 的误差 10−8×3.49×109≈3.5×10110^{-8}\times3.49\times10^{9}\approx3.5\times10^{1}——荒谬,与真值无关。 (3) 补救:把递推反向走(解出 zn−1z_{n-1} 的表达式,使每步乘 1/31/3),并用积分估计式/边界条件把远端初值先定住(1.2.4 原则 (1)、errata N3)。

14. IEEE 双精度(float64): (1) 写出 eps\mathrm{eps} 与单位舍入误差 uu 的定义式并给出数值; (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\mathrm{eps}=\beta^{1-t}=2^{-52}\approx2.22\times10^{-16};u=12β1−t=2−53≈1.11×10−16u=\tfrac12\beta^{1-t}=2^{-53}\approx1.11\times10^{-16}(定义 1.2.10)。 (2) 11 附近的浮点间距是 2−522^{-52},相邻两数中点是 1+2−531+2^{-53}:增量恰为半个 ulp 时按最近舍入回到 1.01.0;要跨过一个网格间距至少需要 2−522^{-52}。 (3) 1016≈253.1510^{16}\approx2^{53.15} 落在 [253,254)[2^{53},2^{54}) 内,该区间间距为 2k−t=253−52=22^{k-t}=2^{53-52}=2,增量 11 只是半个间距,被舍回 101610^{16}。 共同说明:浮点误差取决于数值的绝对位置(间距 2k−t2^{k-t} 随量级成倍变宽),不是「一律 10−1610^{-16}」——大数吃小数是网格太稀的必然(1.2.5 三条后果)。

15. 取 x=(2,−3,6)T\boldsymbol x=(2,-3,6)^{\mathrm T}:算出 ∥x∥1,∥x∥2,∥x∥∞\lVert\boldsymbol x\rVert_1,\lVert\boldsymbol x\rVert_2,\lVert\boldsymbol x\rVert_\infty,逐条核对 1.3.2 的三条显式不等式,并说明「换范数不改变收敛结论」的依据。

解答

15. ∥x∥1=2+3+6=11\lVert\boldsymbol x\rVert_1=2+3+6=11;∥x∥2=4+9+36=49=7\lVert\boldsymbol x\rVert_2=\sqrt{4+9+36}=\sqrt{49}=7;∥x∥∞=6\lVert\boldsymbol x\rVert_\infty=6。 核对(n=3n=3):∥x∥∞=6≤∥x∥2=7≤3×6≈10.392\lVert\boldsymbol x\rVert_\infty=6\le\lVert\boldsymbol x\rVert_2=7\le\sqrt3\times6\approx10.392 ✓;∥x∥2=7≤∥x∥1=11≤3×7≈12.124\lVert\boldsymbol x\rVert_2=7\le\lVert\boldsymbol x\rVert_1=11\le\sqrt3\times7\approx12.124 ✓;∥x∥∞=6≤∥x∥1=11≤3×6=18\lVert\boldsymbol x\rVert_\infty=6\le\lVert\boldsymbol x\rVert_1=11\le3\times6=18 ✓。 依据:等价定理(定理 1.3.3)保证任意两种范数被与 x\boldsymbol x 无关的常数双向夹住,故「是否趋于零」与范数无关;常数 n,nn,\sqrt n 不可再改进([1,…,1]T[1,\dots,1]^{\mathrm T} 与单分量向量取等)。给分:三个范数 22 分;三条核对 1.51.5 分;依据 0.50.5 分。

第三部分 · 推导计算(3 题 × 10 分 = 30 分)

16.(10 分) 用游标卡尺量得圆柱体的尺寸、用天平称得质量:d=12.50d=12.50 mm、h=38.00h=38.00 mm、m=36.60m=36.60 g,仪器给出的误差限分别为 ε(d)=ε(h)=0.02\varepsilon(d)=\varepsilon(h)=0.02 mm、ε(m)=0.02\varepsilon(m)=0.02 g(各量按四舍五入读到末位)。

(1) 按 ρ=4mπd2h\rho=\dfrac{4m}{\pi d^2h} 计算 ρ\rho(取 π=3.14159\pi=3.14159 为准确值),并给出 εr(ρ)\varepsilon_{\mathrm r}(\rho) 与 ε(ρ)\varepsilon(\rho); (2) 判断 ρ\rho 能有几位有效数字,按规定写法重写结果(带单位:1 mm3=10−3 cm31\ \mathrm{mm^3}=10^{-3}\ \mathrm{cm^3},请用 g/cm3\mathrm{g/cm^3} 报出); (3) 要把 εr(ρ)\varepsilon_{\mathrm r}(\rho) 压到 0.3%0.3% 以下,最该先改进哪个量?给出量化依据。

第 16 题答案与评分要点

(1) 由多元传播公式(式 1.1),三个量在 ρ\rho 里各出现一次的相对贡献是 εr(ρ)≈εr(m)+2εr(d)⏟∂ρ/∂d 带因子 2+εr(h),\varepsilon_{\mathrm r}(\rho)\approx\varepsilon_{\mathrm r}(m)+\underbrace{2\varepsilon_{\mathrm r}(d)}_{\partial\rho/\partial d\ \text{带因子}\ 2}+\varepsilon_{\mathrm r}(h), (ρ∝md−2h−1\rho\propto m d^{-2}h^{-1}:对 dd 的相对放大因子是 ∣∂ln⁡ρ/∂ln⁡d∣=2\lvert\partial\ln\rho/\partial\ln d\rvert=2。)代入: εr(m)=0.0236.60≈0.0546%,εr(d)=0.0212.50=0.16% (×2=0.32%),εr(h)=0.0238.00≈0.0526%.\varepsilon_{\mathrm r}(m)=\frac{0.02}{36.60}\approx0.0546\%,\quad \varepsilon_{\mathrm r}(d)=\frac{0.02}{12.50}=0.16\%\ (\times2=0.32\%),\quad \varepsilon_{\mathrm r}(h)=\frac{0.02}{38.00}\approx0.0526\%. εr(ρ)≈0.0546%+0.32%+0.0526%≈0.427%≈4.3×10−3.\varepsilon_{\mathrm r}(\rho)\approx0.0546\%+0.32\%+0.0526\%\approx0.427\%\approx4.3\times10^{-3}. 数值:πd2h=3.14159×156.25×38.00=18653.2 mm3\pi d^2h=3.14159\times156.25\times38.00=18653.2\ \mathrm{mm^3}, ρ=4×36.6018653.19=7.8485×10−3 g/mm3=7.8485 g/cm3,\rho=\frac{4\times36.60}{18653.19}=7.8485\times10^{-3}\ \mathrm{g/mm^3}=7.8485\ \mathrm{g/cm^3}, ε(ρ)≈4.27×10−3×7.8485≈0.0335 g/cm3.\varepsilon(\rho)\approx4.27\times10^{-3}\times7.8485\approx0.0335\ \mathrm{g/cm^3}.

(2) ρ=7.8485=0.78485×101\rho=7.8485=0.78485\times10^1,m=1m=1。检验 nn:12×101−n≥0.0335\tfrac12\times10^{1-n}\ge0.0335——n=1n=1 时 0.5≥0.03350.5\ge0.0335 ✓;n=2n=2 时 0.05≥0.03350.05\ge0.0335 ✓;n=3n=3 时 0.005<0.03350.005\lt0.0335 ✗。故 22 位有效数字,结果写作 ρ≈7.8 g/cm3(误差限 0.034 g/cm3).\rho\approx7.8\ \mathrm{g/cm^3}\quad(\text{误差限}\ 0.034\ \mathrm{g/cm^3}). 注意:把 7.84857.8485 原样抄进报告是错误用法——可靠位数由误差限决定(1.2.8 错误 7)。

(3) 三项贡献中 2εr(d)=0.32%2\varepsilon_{\mathrm r}(d)=0.32\% 占总误差 0.427%0.427\% 的约 75%75\%,所以先改进直径 dd 的测量:换成误差限 0.0040.004 mm 的千分尺,则 εr(ρ)≈0.0546%+2×0.00412.50%+0.0526%=0.0546%+0.064%+0.0526%≈0.171%<0.3% ✓\varepsilon_{\mathrm r}(\rho)\approx0.0546\%+2\times\frac{0.004}{12.50}\%+0.0526\%=0.0546\%+0.064\%+0.0526\%\approx0.171\%\lt0.3\%\ \checkmark 而只把 ε(m)\varepsilon(m) 减半(0.02→0.010.02\to0.01 g)只能把总误差降到 0.40%0.40\%,达不到目标。结论:改进要挑「相对贡献 × 放大因子」最大的那一项。

评分:(1) 相对误差分解 33 分(三个贡献值各 11 分,dd 的因子 22 漏掉扣 1.51.5)+ 数值 ρ\rho 与 ε(ρ)\varepsilon(\rho) 33 分;(2) 位数判定 22 分(判定过程 11 + 按规定写法 11);(3) 选对量与量化依据 22 分。未写单位或未给误差限,各扣 0.50.5 分。

17.(10 分) 设 A=(2−3104−2113).A=\begin{pmatrix}2&amp;-3&amp;1\0&amp;4&amp;-2\1&amp;1&amp;3\end{pmatrix}.

(1) 求 ∥A∥1\lVert A\rVert_1、∥A∥∞\lVert A\rVert_\infty、∥A∥F\lVert A\rVert_F; (2) 求 ∥A∥2\lVert A\rVert_2:写出 ATAA^{\mathrm T}A、说明特征多项式的三个系数怎样不算展开就得到(迹、全部二阶主子式之和、det⁡(ATA)=(det⁡A)2\det(A^{\mathrm T}A)=(\det A)^2),再给出 λmax⁡\lambda_{\max} 与 ∥A∥2\lVert A\rVert_2(数值保留四位有效数字); (3) 核对 ∥A∥2≤∥A∥F\lVert A\rVert_2\le\lVert A\rVert_F 与 ∥A∥2≤∥A∥1∥A∥∞\lVert A\rVert_2\le\sqrt{\lVert A\rVert_1\lVert A\rVert_\infty},并说明为什么 ∥A∥1\lVert A\rVert_1 与 ∥A∥∞\lVert A\rVert_\infty 一般不相等、转置会交换它们。

第 17 题答案与评分要点

(1) 列绝对值和:第 11 列 2+0+1=32+0+1=3,第 22 列 3+4+1=83+4+1=8,第 33 列 1+2+3=61+2+3=6,故 ∥A∥1=max⁡(3,8,6)=8.\lVert A\rVert_1=\max(3,8,6)=8 . 行绝对值和:第 11 行 66,第 22 行 66,第 33 行 55,故 ∥A∥∞=max⁡(6,6,5)=6.\lVert A\rVert_\infty=\max(6,6,5)=6 . ∥A∥F=4+9+1+0+16+4+1+1+9=45≈6.708.\lVert A\rVert_F=\sqrt{4+9+1+0+16+4+1+1+9}=\sqrt{45}\approx6.708 .

(2) 由各列内积(AA 的列为 (2,0,1)T,(−3,4,1)T,(1,−2,3)T(2,0,1)^{\mathrm T},(-3,4,1)^{\mathrm T},(1,-2,3)^{\mathrm T}): ATA=(5−55−526−85−814).A^{\mathrm T}A=\begin{pmatrix}5&amp;-5&amp;5\\-5&amp;26&amp;-8\\5&amp;-8&amp;14\end{pmatrix}. 特征多项式 det⁡(ATA−λI)=−λ3+c2λ2−c1λ+c0\det(A^{\mathrm T}A-\lambda I)=-\lambda^3+c_2\lambda^2-c_1\lambda+c_0 的三个系数不必硬展开:

系数 来源 值
二次项(迹) tr⁡(ATA)\operatorname{tr}(A^{\mathrm T}A) 5+26+14=455+26+14=45
一次项 全部二阶主子式之和 (5⋅26−25)+(5⋅14−25)+(26⋅14−64)=105+45+300=450(5\cdot26-25)+(5\cdot14-25)+(26\cdot14-64)=105+45+300=450
常数项 det⁡(ATA)=(det⁡A)2\det(A^{\mathrm T}A)=(\det A)^2 det⁡A=30⇒900\det A=30\Rightarrow 900

(det⁡A=2(12+2)+3(0+2)+(0−4)=28+6−4=30\det A=2(12+2)+3(0+2)+(0-4)=28+6-4=30。)于是 λ3−45λ2+450λ−900=0 ⟹ λmax⁡≈31.70,\lambda^3-45\lambda^2+450\lambda-900=0\ \Longrightarrow\ \lambda_{\max}\approx31.70, ∥A∥2=31.70≈5.630.\lVert A\rVert_2=\sqrt{31.70}\approx5.630 . (代入核对:31.703−45×31.702+450×31.70−900≈031.70^3-45\times31.70^2+450\times31.70-900\approx0 ✓。)

(3) 核对:∥A∥2=5.630≤∥A∥F=6.708\lVert A\rVert_2=5.630\le\lVert A\rVert_F=6.708 ✓;∥A∥1∥A∥∞=48≈6.928≥5.630\sqrt{\lVert A\rVert_1\lVert A\rVert_\infty}=\sqrt{48}\approx6.928\ge5.630 ✓,且 ∥A∥2≤∥A∥F\lVert A\rVert_2\le\lVert A\rVert_F 与 ∥A∥2≤∥A∥1∥A∥∞\lVert A\rVert_2\le\sqrt{\lVert A\rVert_1\lVert A\rVert_\infty} 对任意矩阵都成立。 为什么两个范数一般不等:∥A∥1\lVert A\rVert_1 数的是列、∥A∥∞\lVert A\rVert_\infty 数的是行;转置把行列互换,所以 ∥AT∥1=∥A∥∞\lVert A^{\mathrm T}\rVert_1=\lVert A\rVert_\infty、∥AT∥∞=∥A∥1\lVert A^{\mathrm T}\rVert_\infty=\lVert A\rVert_1(1.3.4 记忆法)。

评分:(1) 三个范数各 11 分;(2) ATAA^{\mathrm T}A 1.51.5 分、三个系数各 11 分、λmax⁡\lambda_{\max} 与 ∥A∥2\lVert A\rVert_2 1.51.5 分;(3) 两条核对 1.51.5 分、转置交换关系 1.51.5 分。把 ∥A∥2\lVert A\rVert_2 写成 λmax⁡(A)\lambda_{\max}(A)(本题真正的主特征值是 4.44224.4422,且另两个特征值为复数——ρ(A)=4.4422\rho(A)=4.4422 与 ∥A∥2=5.630\lVert A\rVert_2=5.630 是两个不同的数)本节 00 分。

18.(10 分) 计算积分序列 yn=∫01xnx+3 dx(n=0,1,2,… ).y_n=\int_0^1\frac{x^n}{x+3},dx\qquad(n=0,1,2,\dots).

(1) 推导递推关系,并给出 y0y_0 的表达式与数值; (2) 判断正向递推 yn=1n−3yn−1y_n=\dfrac1n-3y_{n-1} 的稳定性:写出误差递推式与 2020 步的放大倍数;结合 y0y_0 的四位十进制近似说明后果; (3) 写出稳定的反向递推格式,并用积分上下界给出 y50y_{50} 的起点估计区间与取值建议。

第 18 题答案与评分要点

(1) 因 xn+3xn−1=xn−1(x+3)x^n+3x^{n-1}=x^{n-1}(x+3),分子恰被分母整除: yn+3yn−1=∫01xn+3xn−1x+3 dx=∫01xn−1dx=1n(n≥1),y_n+3y_{n-1}=\int_0^1\frac{x^n+3x^{n-1}}{x+3}\,dx=\int_0^1x^{n-1}dx=\frac1n\qquad(n\ge1), y0=∫01dxx+3=ln⁡4−ln⁡3=ln⁡43≈0.28768.y_0=\int_0^1\frac{dx}{x+3}=\ln4-\ln3=\ln\frac43\approx0.28768 .

(2) 设 yn−1y_{n-1} 带误差 δn−1\delta_{n-1},由 yn=1n−3yn−1y_n=\frac1n-3y_{n-1} 得 δn=−3δn−1 ⟹ ∣δn∣=3n∣δ0∣,\delta_n=-3\delta_{n-1}\ \Longrightarrow\ \lvert\delta_n\rvert=3^n\lvert\delta_0\rvert, 每步放大 33 倍;320=3.487×1093^{20}=3.487\times10^{9}(约 3.5×1093.5\times10^{9} 倍)。四位十进制下 y0≈0.2877y_0\approx0.2877,其舍入误差约 ∣0.2877−0.287682∣≈1.8×10−5\lvert0.2877-0.287682\vert\approx1.8\times10^{-5},传到 y20y_{20} 时放大约 6.3×1046.3\times10^{4}——而 y20∼10−2y_{20}\sim10^{-2} 量级,结果完全失真(教材实测 n=30n=30 时前推值已达 −2.0×1016-2.0\times10^{16},真值 5.4×10−35.4\times10^{-3})。结论:∣β∣=3>1\lvert\beta\rvert=3\gt1,正向递推不稳定,提高字长只推迟灾难(@ex 类比见 1.2.4)。

(3) 反向:由 yn+3yn−1=1ny_n+3y_{n-1}=\frac1n 解出 yn−1=13n−13yn,β=−13, ∣β∣<1 (每步误差缩小 3 倍).y_{n-1}=\frac{1}{3n}-\frac13y_n,\qquad\beta=-\frac13,\ \lvert\beta\rvert\lt1\ (\text{每步误差缩小}\ 3\ \text{倍}). 起点:[0,1][0,1] 上 14≤1x+3≤13\frac14\le\frac{1}{x+3}\le\frac13,故 xn4≤xnx+3≤xn3\frac{x^n}{4}\le\frac{x^n}{x+3}\le\frac{x^n}{3},积分得 14(n+1)≤yn≤13(n+1) ⟹ y50∈[1204,1153]=[0.0049020, 0.0065359].\frac{1}{4(n+1)}\le y_n\le\frac{1}{3(n+1)}\ \Longrightarrow\ y_{50}\in\Bigl[\frac{1}{204},\frac{1}{153}\Bigr]=[0.0049020,\ 0.0065359]. 取区间中点 y50≈0.0057190y_{50}\approx0.0057190,起点误差不超过半宽 ≈8.2×10−4\approx8.2\times10^{-4};反向递推每步把它压缩 33 倍,传到 y0y_0 时被乘 (1/3)50≈1.4×10−24(1/3)^{50}\approx1.4\times10^{-24}——起点误差与途中舍入误差都被彻底压掉。

评分:(1) 递推推导 22 分 + y0y_0 11 分;(2) 误差递推式 1.51.5 分、3203^{20} 数值 11 分、后果说明 1.51.5 分;(3) 反向格式 1.51.5 分、上下界 11 分、中点与结论 0.50.5 分。只写「反向更稳」而不给误差因子,扣 22 分。

第四部分 · 开放论述(2 题 × 10 分 = 20 分)

19.(10 分) 某小组用「游标卡尺 + 电子天平」测量一批圆柱体试样的密度(上题第 16 题的公式与仪器),要交一份误差预算报告。要求: (1) 写出这条测量链上四类误差各自的代表来源(每类至少一条),并说明本预算公式只覆盖了其中哪一部分; (2) 在「改进仪器」与「改进公式」之间做选择,给出判据; (3) 写出报告的结果规范(一条结果行里应包含哪些字段); (4) 给出至少两条可测量的失败判据。

第 19 题参考答案与评分要点

(1) 四类误差的来源清单与预算边界

类别 本链上的代表来源 能否进本预算公式
模型误差 把试样当理想正圆柱(实际有圆度 / 锥度误差)、忽略称量时的空气浮力 不进(属于公式本身,V=πd2h/4V=\pi d^2h/4 已假定是圆柱)
观测误差 卡尺读数分辨与对准、天平示值漂移、温湿度导致的尺寸/质量漂移 部分进(以 ε(d),ε(h),ε(m)\varepsilon(d),\varepsilon(h),\varepsilon(m) 的形式进)
截断误差 本次计算中 π\pi 只取 3.141593.14159(若取近似值) 可进(当作一个额外输入量,εr(π)≈1.6×10−6\varepsilon_{\mathrm r}(\pi)\approx1.6\times10^{-6},可忽略)
舍入误差 商与乘积的浮点执行(双精度,相对 u∼10−16u\sim10^{-16}) 可进(量级远小于 10−310^{-3},实际可忽略)

关键结论:预算公式算出的 0.43%0.43\% 只覆盖观测误差(+ 可忽略的截断/舍入);模型误差一旦显著(例如锥度 0.10.1 mm 级),任何减小 ε(d)\varepsilon(d) 的努力都是白费——必须先做模型检查。

(2) 改进仪器还是改进公式:判据 判据写法:「先做残差检验,再看相对贡献占比」。具体:① 对同一试样重复测量 n≥10n\ge10 次,把样本标准差与仪器误差限对比——若实测离散显著大于仪器限(例如 0.050.05 mm vs 0.020.02 mm),说明操作/环境主导,改进仪器无效,应先改测量方法(支架、恒温、多点平均);② 若离散与仪器限一致,则看各项相对贡献:本链中 2εr(d)2\varepsilon_{\mathrm r}(d) 占 75%75\%,故换更高精度的直径测量(千分尺)比改公式有效;③ 只有当模型误差成为主导(锥度、椭圆度使 VV 的真值与 πd2h/4\pi d^2h/4 偏离超过 0.4%0.4\%)时,才改公式——例如沿高度测三个截面、用分段圆柱求和。

(3) 结果规范(一行结果必须含四个字段)

ρ=7.8 g/cm3\rho = 7.8\ \mathrm{g/cm^3}(ε(ρ)=0.03 g/cm3\varepsilon(\rho)=0.03\ \mathrm{g/cm^3},εr=0.43%\varepsilon_{\mathrm r}=0.43\%,22 位有效数字;依据 ρ=4m/(πd2h)\rho=4m/(\pi d^2h),ε(d)=ε(h)=0.02\varepsilon(d)=\varepsilon(h)=0.02 mm,ε(m)=0.02\varepsilon(m)=0.02 g,双精度计算)

四个字段:值 + 绝对(或相对)误差限 + 有效数字位数 + 数据与依据;若引用了范数(向量情形)还必须注明范数,因为它决定误差界的大小。

(4) 失败判据(≥2 条,可测量) ① 重复测量的极差超过仪器误差限的 44 倍(range>4ε\mathrm{range}\gt4\varepsilon)→ 操作或环境主导,预算表失效,先查测量流程; ② 预算预测的相对误差与「高精度参考测量」(如用千分尺 + 精密天平复测同一试样)的实际偏差相差超过一个数量级(实测 ≫\gg 预测)→ 存在未计入的模型误差; ③ 计算结果写出的位数与 ε(ρ)\varepsilon(\rho) 不匹配(如报了 7.84857.8485)→ 报告不合格,直接退回; ④ 换了更高精度仪器后总误差没有下降(改进前后 εr\varepsilon_{\mathrm r} 差异 <10%&lt;10\%)→ 说明主导项判断错了,回到 ①。

评分:(1) 33 分(四类各 0.50.5、预算边界 11);(2) 33 分(判据必须可执行:残差检验 + 贡献占比各 1.51.5);(3) 22 分(四个字段齐全);(4) 22 分(两条各 11 分,须可测量)。只写「提高精度」而不给判据,(2) 得 00 分。

20.(10 分) 一个迭代求解程序用 ∥xk+1−xk∥∞≤10−12\lVert\boldsymbol x_{k+1}-\boldsymbol x_k\rVert_\infty\le10^{-12} 作为唯一停机判据。请: (1) 写出「相邻差」与「真实误差」之间的关系式(设迭代格式为 xk+1=Hxk+d\boldsymbol x_{k+1}=H\boldsymbol x_k+\boldsymbol d,用 ∥H∥\lVert H\rVert 表示),并说明什么时候「相邻差小」推不出「误差小」; (2) 给出三条改进方案,并说明每条治的是什么病; (3) 用第一章的结论解释:这条判据为什么不能充当误差估计,它与「残差 ≠ 误差」是同一件事的两面吗; (4) 给出至少两条可测量的失败判据。

第 20 题参考答案与评分要点

(1) 关系式。 记精确解 x∗=(I−H)−1d\boldsymbol x^*=(I-H)^{-1}\boldsymbol d、误差 ek=xk−x∗\boldsymbol e_k=\boldsymbol x_k-\boldsymbol x^*,则 ek+1=Hek\boldsymbol e_{k+1}=H\boldsymbol e_k、ek=Hke0\boldsymbol e_k=H^k\boldsymbol e_0,于是 xk+1−xk=ek+1−ek=(H−I)ek ⟹ ek=(H−I)−1(xk+1−xk).\boldsymbol x_{k+1}-\boldsymbol x_k=\boldsymbol e_{k+1}-\boldsymbol e_k=(H-I)\boldsymbol e_k\ \Longrightarrow\ \boldsymbol e_k=(H-I)^{-1}(\boldsymbol x_{k+1}-\boldsymbol x_k). 当 ∥H∥=:q<1\lVert H\rVert=:q\lt1 时(定理 1.3.6 保证 I−HI-H 可逆)取范数放缩得 ∥ek∥≤q1−q∥xk+1−xk∥以及∥e0∥≤∥x1−x0∥1−q,\lVert\boldsymbol e_k\rVert\le\frac{q}{1-q}\lVert\boldsymbol x_{k+1}-\boldsymbol x_k\rVert\quad\text{以及}\quad\lVert\boldsymbol e_0\rVert\le\frac{\lVert\boldsymbol x_1-\boldsymbol x_0\rVert}{1-q}, ∥ek∥≤qk∥e0∥≤qk1−q∥x1−x0∥.\lVert\boldsymbol e_k\rVert\le q^k\lVert\boldsymbol e_0\rVert\le\frac{q^k}{1-q}\lVert\boldsymbol x_1-\boldsymbol x_0\rVert . 什么时候推不出误差小:放大因子 q1−q\dfrac{q}{1-q} 在 q→1q\to1 时爆炸(q=0.99q=0.99 时约 9999、q=0.9999q=0.9999 时约 10410^4):相邻差可以小到 10−1210^{-12},而真实误差仍在 10−810^{-8} 量级甚至更大;若 ∥H∥≥1\lVert H\rVert\ge1 则上式根本给不出界(此时才需要补充篇的谱半径与 ρ<1\rho\lt1 的渐近结论)。

(2) 三条改进方案与它治的病

方案 做法 治什么
① 用后验误差估计替代相邻差 先估 q=∥H∥q=\lVert H\rVert(或每步估比值 ∥Δxk+1∥/∥Δxk∥→q\lVert\Delta\boldsymbol x_{k+1}\rVert/\lVert\Delta\boldsymbol x_k\rVert\to q),再用 ∥e∥≲q1−q∥Δx∥\lVert\boldsymbol e\rVert\lesssim\frac{q}{1-q}\lVert\Delta\boldsymbol x\rVert 治「判据与实际误差差一个放大因子」
② 双判据:增广残差 + 条件数 对方程组类问题同时监控 cond(A)⋅∥r∥/∥b∥\mathrm{cond}(A)\cdot\lVert\boldsymbol r\rVert/\lVert\boldsymbol b\rVert(1.3.6) 治「残差/相邻差都小、误差却大」的病态问题
③ 相对化 + 注明范数 改判据为 ∥xk+1−xk∥∥xk+1∥+1≤τ\dfrac{\lVert\boldsymbol x_{k+1}-\boldsymbol x_k\rVert}{\lVert\boldsymbol x_{k+1}\rVert+1}\le\tau,并在报告里写明用的是哪个范数 治「∥⋅∥∞\lVert\cdot\rVert_\infty 只盯最坏分量、尺度大时绝对阈值失效」

(补充:若 qq 接近 11、步数不可预测,正确做法是换更好的预条件 / 换算法,而不是继续调阈值。)

(3) 为什么它不能当误差估计。 停机判据度量的是算法是否稳定下来(增量信息 / 后向信息),误差度量的是解离真值多远(前向信息);两者的桥梁恰好是问题的条件数或迭代矩阵的 q/(1−q)q/(1-q)。这与 1.2.6 的「conditioning 与 stability 必须分开谈」、1.2.6 的「求解器只能保证残差、用户要的是误差」同一句话;矩阵版 e=A−1r\boldsymbol e=A^{-1}\boldsymbol r 的完整推导在补充篇 2.4.4。是同一件事的两面:残差对应 Ax=bA\boldsymbol x=\boldsymbol b 的「相邻/后向」信息,相邻差对应迭代格式的增量信息,二者都要乘上放大因子才能变成误差。

(4) 失败判据(≥2 条,可测量) ① 停机后用一个高精度参考解(或第二种算法/更高精度重算)对照,实测相对误差 >10−6&gt;10^{-6}——说明 q/(1−q)q/(1-q) 被低估,判据过松; ② 停机前 ∥Δxk+1∥/∥Δxk∥\lVert\Delta\boldsymbol x_{k+1}\rVert/\lVert\Delta\boldsymbol x_k\rVert 持续大于 0.90.9(衰减比不趋于常数 q<1q\lt1)——迭代矩阵谱半径过于接近 11,步数不可预测,应换算法或预条件; ③ 相邻差已低于阈值,但残差(或 cond⋅∥r∥/∥b∥\mathrm{cond}\cdot\lVert\boldsymbol r\rVert/\lVert\boldsymbol b\rVert)仍在下降——说明是「增量停滞」而不是「误差已小」,必须继续迭代或改判据; ④ 同一问题换用 ∥⋅∥2\lVert\cdot\rVert_2 后停机步数变化超过 1010 倍——说明阈值 τ\tau 与范数不匹配,报告里没有注明范数。

评分:(1) 44 分(误差递推 11、关系式 1.51.5、放大因子与失效条件 1.51.5);(2) 33 分(三条各 11 分,须写出「治什么病」);(3) 22 分(区分前向/后向 + 指出桥梁是条件数);(4) 11 分(两条可测量判据)。

评分标准

部分题量分值给分要点
一 · 概念判断1030判断 11 分 + 理由 / 反例 22 分;只写对错得 11 分
二 · 填空与简答520结果正确 22 分 + 关键说明 22 分;数值 / 单位错一处扣 11 分
三 · 推导计算330首行或公式 33 分 + 链条完整 44 分 + 位数 / 边界等收口 33 分;只写结论最多 33 分
四 · 开放论述220方案具体 44 分 + 判据可测量 44 分 + 验证计划 22 分

段位换算:≥90\ge90 优秀(能给同学讲清机制);75 ⁣− ⁣8975!-!89 良好(能独立做误差分析);60 ⁣− ⁣7460!-!74 及格(公式记得住、条件记不住);<60&lt;60 回到对应节的「即时自测」与习题集重做。自评纪律:① 判断理由必须包含机制词(末位半单位、放大因子、次可乘、相消、条件数、算子范数),否则按半对给 11 分;② 计算题凡是没给误差限或没给位数的,从收口里扣 11 分;③ 开放题只写方案、不写失败判据的,第四部分上限 66 分。

漏点定位表

错题回去读立刻补做
1、5、12、141.2.2 判定流程、1.2.5、1.1.3习题 A3、习题 C4、回忆卡 A4/A7
2、10定义 1.2.6 后的推导、引理 1.2.1习题 A2、回忆卡 D2/D4
3、11、151.2.1 鸟瞰表、定义 1.2.7、定理 1.3.3习题 B2、习题 A1
4、13、181.2.4 原则 (1) 与处方表习题 C6、习题 C3
6、7、8、9、17、201.3.2、1.3.4、1.3.7、定理 1.3.7、1.3.5 预告框习题 B3、习题 B4、习题 B6
16、191.2.3 四则表与式 (1.1)、真实应用案例习题 C2、习题 D1

最后一句:本卷的目的不是打分,而是把「以为会了」的地方暴露出来。凡「能算出数但说不清依据」的题,回到对应节的 intuition 框与回忆卡当天再练一遍;同一类错误出现两次以上,那不是粗心,而是概念缺口。