《数值分析》期末试卷 (精选01)

大学生期末考试复习平台前往平台

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

一、填空题(每小题 4 分,共 40 分)

  1. 为了提高数值计算精度,当正数 x x x 充分大时,应将 ln ⁡ ( x − x 2 − 1 ) \ln(x - \sqrt{x^2 - 1}) ln(xx21 ) 改写为 ______

::::: details 查看答案与解析
答案: − ln ⁡ ( x + x 2 − 1 ) -\ln(x + \sqrt{x^2 - 1}) ln(x+x21 )

详细解析:
x x x 充分大时, x x x x 2 − 1 \sqrt{x^2 - 1} x21 是两个非常接近的大数,直接相减会产生相近数相减导致的有效数字严重丢失。
利用分子有理化技巧:

x − x 2 − 1 = ( x − x 2 − 1 ) ( x + x 2 − 1 ) x + x 2 − 1 = x 2 − ( x 2 − 1 ) x + x 2 − 1 = 1 x + x 2 − 1 x - \sqrt{x^2 - 1} = \frac{(x - \sqrt{x^2 - 1})(x + \sqrt{x^2 - 1})}{x + \sqrt{x^2 - 1}} = \frac{x^2 - (x^2 - 1)}{x + \sqrt{x^2 - 1}} = \frac{1}{x + \sqrt{x^2 - 1}} xx21 =x+x21 (xx21 )(x+x21 )=x+x21 x2(x21)=x+x21 1

代入原式:

ln ⁡ ( x − x 2 − 1 ) = ln ⁡ ( 1 x + x 2 − 1 ) = − ln ⁡ ( x + x 2 − 1 ) \ln(x - \sqrt{x^2 - 1}) = \ln\left(\frac{1}{x + \sqrt{x^2 - 1}}\right) = -\ln(x + \sqrt{x^2 - 1}) ln(xx21 )=ln(x+x21 1)=ln(x+x21 )

这样改写后,运算变为了加法和求对数,避免了精度损失。

难度:
考点: #数值稳定性 #相近数相减
:::: details 💡 学习锦囊
📖 核心原则: 在数值计算中,应尽量避免两个相近的数相减。常见的处理方法包括分子有理化、三角公式变换或使用泰勒展开。
::::
:::: details 🔄 举一反三
x x x 很小时,如何改写 1 − cos ⁡ x x 2 \frac{1 - \cos x}{x^2} x21cosx 以提高精度?
::: details 查看答案与解析
答案: 2 sin ⁡ 2 ( x / 2 ) x 2 \frac{2\sin^2(x/2)}{x^2} x22sin2(x/2)。利用二倍角公式避免了 1 − cos ⁡ x 1 - \cos x 1cosx 的精度损失。
:::
::::
:::::

  1. 已知 f ( x 2 , x 3 ) = f ( x 3 ) − f ( x 2 ) x 3 − x 2 = 6 − 1 4 − 2 = 5 2 f(x_2, x_3) = \frac{f(x_3) - f(x_2)}{x_3 - x_2} = \frac{6 - 1}{4 - 2} = \frac{5}{2} f(x2,x3)=x3x2f(x3)f(x2)=4261=25,且 f ( x 1 ) = 3 , x 1 = 0 f(x_1) = 3, x_1 = 0 f(x1)=3,x1=0,则二阶差商 f [ x 1 , x 2 , x 3 ] = f[x_1, x_2, x_3] = f[x1,x2,x3]= ______

::::: details 查看答案与解析
答案: 0.875 0.875 0.875 (或 7 / 8 7/8 7/8

详细解析:

  1. 计算一阶差商 f [ x 1 , x 2 ] f[x_1, x_2] f[x1,x2]
    已知 x 1 = 0 , f ( x 1 ) = 3 x_1 = 0, f(x_1) = 3 x1=0,f(x1)=3;从题干推知 x 2 = 2 , f ( x 2 ) = 1 x_2 = 2, f(x_2) = 1 x2=2,f(x2)=1
    f [ x 1 , x 2 ] = f ( x 2 ) − f ( x 1 ) x 2 − x 1 = 1 − 3 2 − 0 = − 1 f[x_1, x_2] = \frac{f(x_2) - f(x_1)}{x_2 - x_1} = \frac{1 - 3}{2 - 0} = -1 f[x1,x2]=x2x1f(x2)f(x1)=2013=1
  2. 计算二阶差商 f [ x 1 , x 2 , x 3 ] f[x_1, x_2, x_3] f[x1,x2,x3]
    已知 f [ x 2 , x 3 ] = 2.5 f[x_2, x_3] = 2.5 f[x2,x3]=2.5 x 1 = 0 , x 3 = 4 x_1 = 0, x_3 = 4 x1=0,x3=4
    f [ x 1 , x 2 , x 3 ] = f [ x 2 , x 3 ] − f [ x 1 , x 2 ] x 3 − x 1 = 2.5 − ( − 1 ) 4 − 0 = 3.5 4 = 0.875 f[x_1, x_2, x_3] = \frac{f[x_2, x_3] - f[x_1, x_2]}{x_3 - x_1} = \frac{2.5 - (-1)}{4 - 0} = \frac{3.5}{4} = 0.875 f[x1,x2,x3]=x3x1f[x2,x3]f[x1,x2]=402.5(1)=43.5=0.875

难度: ⭐⭐
考点: #差商定义 #均差表
:::: details 💡 学习锦囊
📖 相关公式: n n n 阶差商公式为 f [ x 0 , x 1 , … , x n ] = f [ x 1 , … , x n ] − f [ x 0 , … , x n − 1 ] x n − x 0 f[x_0, x_1, \dots, x_n] = \frac{f[x_1, \dots, x_n] - f[x_0, \dots, x_{n-1}]}{x_n - x_0} f[x0,x1,,xn]=xnx0f[x1,,xn]f[x0,,xn1]
::::
:::: details 🔄 举一反三
f ( x ) f(x) f(x) 是一个 2 次多项式,已知其二阶差商为 1 1 1,则其最高次项系数为?
::: details 查看答案与解析
答案: 1。根据性质: n n n 次多项式的 n n n 阶差商等于其最高次项系数。
:::
::::
:::::

  1. X = ( 2 , 3 , − 4 ) T X = (2, 3, -4)^T X=(2,3,4)T,则其向量范数 ∥ X ∥ 2 = \| X\|_2 = X2= ______

::::: details 查看答案与解析
答案: 29 \sqrt{29} 29

详细解析:
向量的 2-范数(欧几里得范数)定义为各分量平方和的开方:

∥ X ∥ 2 = ∑ i = 1 n x i 2 \|X\|_2 = \sqrt{\sum_{i=1}^n x_i^2} X2=i=1nxi2

代入分量:

∥ X ∥ 2 = 2 2 + 3 2 + ( − 4 ) 2 = 4 + 9 + 16 = 29 \|X\|_2 = \sqrt{2^2 + 3^2 + (-4)^2} = \sqrt{4 + 9 + 16} = \sqrt{29} X2=22+32+(4)2 =4+9+16 =29

难度:
考点: #向量范数 #2-范数
:::: details 💡 学习锦囊
📖 常用范数:

  • 1-范数 ∥ X ∥ 1 = ∑ ∣ x i ∣ \|X\|_1 = \sum |x_i| X1=xi(分量绝对值之和)
  • ∞ \infty -范数 ∥ X ∥ ∞ = max ⁡ ∣ x i ∣ \|X\|_\infty = \max |x_i| X=maxxi(分量绝对值最大值)
  • 2-范数 ∥ X ∥ 2 = X T X \|X\|_2 = \sqrt{X^T X} X2=XTX
    ::::
    :::: details 🔄 举一反三
    求向量 X = ( 1 , − 5 , 2 ) T X = (1, -5, 2)^T X=(1,5,2)T 的 1-范数和 ∞ \infty -范数。
    ::: details 查看答案与解析
    答案: ∥ X ∥ 1 = 8 \|X\|_1 = 8 X1=8 ∥ X ∥ ∞ = 5 \|X\|_\infty = 5 X=5
    解析: ∥ X ∥ 1 = ∣ 1 ∣ + ∣ − 5 ∣ + ∣ 2 ∣ = 8 \|X\|_1 = |1| + |-5| + |2| = 8 X1=∣1∣+5∣+∣2∣=8 ∥ X ∥ ∞ = max ⁡ ( 1 , 5 , 2 ) = 5 \|X\|_\infty = \max(1, 5, 2) = 5 X=max(1,5,2)=5
    :::
    ::::
    :::::
  1. 取步长 h = 0.1 h = 0.1 h=0.1,用欧拉法解初值问题 { y ′ = y x 2 + y y ( 1 ) = 1 \begin{cases} y' = \frac{y}{x^2} + y \\ y(1) = 1 \end{cases} {y=x2y+yy(1)=1 的计算公式为 ______

::::: details 查看答案与解析
答案: y n + 1 = y n ( 1.1 + 0.1 x n 2 ) y_{n+1} = y_n(1.1 + \frac{0.1}{x_n^2}) yn+1=yn(1.1+xn20.1)

详细解析:
欧拉法(Euler Method)的迭代公式为:

y n + 1 = y n + h f ( x n , y n ) y_{n+1} = y_n + h f(x_n, y_n) yn+1=yn+hf(xn,yn)

本题中 f ( x , y ) = y x 2 + y f(x, y) = \frac{y}{x^2} + y f(x,y)=x2y+y h = 0.1 h = 0.1 h=0.1
代入公式:

y n + 1 = y n + 0.1 ( y n x n 2 + y n ) y_{n+1} = y_n + 0.1 \left( \frac{y_n}{x_n^2} + y_n \right) yn+1=yn+0.1(xn2yn+yn)

提取公因子 y n y_n yn

y n + 1 = y n ( 1 + 0.1 ( 1 x n 2 + 1 ) ) = y n ( 1.1 + 0.1 x n 2 ) y_{n+1} = y_n \left( 1 + 0.1 \left( \frac{1}{x_n^2} + 1 \right) \right) = y_n \left( 1.1 + \frac{0.1}{x_n^2} \right) yn+1=yn(1+0.1(xn21+1))=yn(1.1+xn20.1)

难度: ⭐⭐
考点: #欧拉法 #常微分方程数值解
:::: details 💡 学习锦囊
📖 相关公式:

  • 显式欧拉法 y n + 1 = y n + h f ( x n , y n ) y_{n+1} = y_n + h f(x_n, y_n) yn+1=yn+hf(xn,yn)(一阶精度,条件收敛)
  • 隐式欧拉法 y n + 1 = y n + h f ( x n + 1 , y n + 1 ) y_{n+1} = y_n + h f(x_{n+1}, y_{n+1}) yn+1=yn+hf(xn+1,yn+1)(一阶精度,绝对收敛)
    ::::
    :::: details 🔄 举一反三
    若改用后退欧拉法(隐式)解上述问题,其迭代公式是什么?
    ::: details 查看答案与解析
    答案: y n + 1 = y n 1 − 0.1 ( 1 / x n + 1 2 + 1 ) y_{n+1} = \frac{y_n}{1 - 0.1(1/x_{n+1}^2 + 1)} yn+1=10.1(1/xn+12+1)yn
    解析: 后退欧拉公式为 y n + 1 = y n + h ( y n + 1 x n + 1 2 + y n + 1 ) y_{n+1} = y_n + h(\frac{y_{n+1}}{x_{n+1}^2} + y_{n+1}) yn+1=yn+h(xn+12yn+1+yn+1),移项得 y n + 1 ( 1 − h ( 1 / x n + 1 2 + 1 ) ) = y n y_{n+1}(1 - h(1/x_{n+1}^2 + 1)) = y_n yn+1(1h(1/xn+12+1))=yn
    :::
    ::::
    :::::
  1. 已知 G = ( 10 12 1 1 ) G = \begin{pmatrix} 10 & 12 \\ 1 & 1 \end{pmatrix} G=(101121),则其无穷范数下的条件数 cond ⁡ ∞ ( G ) = \operatorname{cond}_{\infty}(G) = cond(G)= ______

::::: details 查看答案与解析
答案: 143 143 143

详细解析:

  1. 计算 ∥ G ∥ ∞ \|G\|_\infty G
    ∥ G ∥ ∞ = max ⁡ ( 10 + 12 , 1 + 1 ) = 22 \|G\|_\infty = \max(10+12, 1+1) = 22 G=max(10+12,1+1)=22
  2. 计算 G − 1 G^{-1} G1
    行列式 det ⁡ ( G ) = 10 × 1 − 12 × 1 = − 2 \det(G) = 10 \times 1 - 12 \times 1 = -2 det(G)=10×112×1=2
    G − 1 = 1 − 2 ( 1 − 12 − 1 10 ) = ( − 0.5 6 0.5 − 5 ) G^{-1} = \frac{1}{-2} \begin{pmatrix} 1 & -12 \\ -1 & 10 \end{pmatrix} = \begin{pmatrix} -0.5 & 6 \\ 0.5 & -5 \end{pmatrix} G1=21(111210)=(0.50.565)
  3. 计算 ∥ G − 1 ∥ ∞ \|G^{-1}\|_\infty G1
    ∥ G − 1 ∥ ∞ = max ⁡ ( 6.5 , 5.5 ) = 6.5 \|G^{-1}\|_\infty = \max(6.5, 5.5) = 6.5 G1=max(6.5,5.5)=6.5
  4. 计算条件数
    cond ⁡ ∞ ( G ) = ∥ G ∥ ∞ ⋅ ∥ G − 1 ∥ ∞ = 22 × 6.5 = 143 \operatorname{cond}_{\infty}(G) = \|G\|_\infty \cdot \|G^{-1}\|_\infty = 22 \times 6.5 = 143 cond(G)=GG1=22×6.5=143

难度: ⭐⭐
考点: #矩阵条件数 #无穷范数
:::: details 💡 学习锦囊
📖 核心概念: 条件数 cond ⁡ ( A ) = ∥ A ∥ ⋅ ∥ A − 1 ∥ \operatorname{cond}(A) = \|A\| \cdot \|A^{-1}\| cond(A)=AA1

  • 它衡量了方程组解对输入数据变化的敏感程度。
  • 条件数越大,矩阵越病态(Ill-conditioned)。
    ::::
    :::: details 🔄 举一反三
    计算该矩阵 G G G 的 1-范数条件数 cond ⁡ 1 ( G ) \operatorname{cond}_1(G) cond1(G)
    ::: details 查看答案与解析
    答案: 143 143 143
    解析: ∥ G ∥ 1 = max ⁡ ( 11 , 13 ) = 13 \|G\|_1 = \max(11, 13) = 13 G1=max(11,13)=13 ∥ G − 1 ∥ 1 = max ⁡ ( 1 , 11 ) = 11 \|G^{-1}\|_1 = \max(1, 11) = 11 G11=max(1,11)=11;故 cond ⁡ 1 ( G ) = 13 × 11 = 143 \operatorname{cond}_1(G) = 13 \times 11 = 143 cond1(G)=13×11=143
    :::
    ::::
    :::::
  1. 用牛顿下山法求解方程 x 5 3 − 2 x = 0 \frac{x^5}{3} - 2x = 0 3x52x=0 根的迭代公式是 ______

::::: details 查看答案与解析
答案: x k + 1 = x k − λ x k 5 − 6 x k 5 x k 4 − 6 x_{k+1} = x_k - \lambda \frac{x_k^5 - 6x_k}{5x_k^4 - 6} xk+1=xkλ5xk46xk56xk

详细解析:

  1. 确定函数及其导数
    f ( x ) = 1 3 x 5 − 2 x    ⟹    f ′ ( x ) = 5 3 x 4 − 2 f(x) = \frac{1}{3}x^5 - 2x \implies f'(x) = \frac{5}{3}x^4 - 2 f(x)=31x52xf(x)=35x42
  2. 写出牛顿迭代部分
    f ( x k ) f ′ ( x k ) = 1 3 x k 5 − 2 x k 5 3 x k 4 − 2 = x k 5 − 6 x k 5 x k 4 − 6 \frac{f(x_k)}{f'(x_k)} = \frac{\frac{1}{3}x_k^5 - 2x_k}{\frac{5}{3}x_k^4 - 2} = \frac{x_k^5 - 6x_k}{5x_k^4 - 6} f(xk)f(xk)=35xk4231xk52xk=5xk46xk56xk
  3. 引入下山因子 λ \lambda λ
    牛顿下山法公式为: x k + 1 = x k − λ f ( x k ) f ′ ( x k ) x_{k+1} = x_k - \lambda \frac{f(x_k)}{f'(x_k)} xk+1=xkλf(xk)f(xk)
    代入得: x k + 1 = x k − λ x k 5 − 6 x k 5 x k 4 − 6 x_{k+1} = x_k - \lambda \frac{x_k^5 - 6x_k}{5x_k^4 - 6} xk+1=xkλ5xk46xk56xk

难度: ⭐⭐
考点: #牛顿下山法 #非线性方程求根
:::: details 💡 学习锦囊
📖 为什么需要下山法? 标准牛顿法虽然具有二阶收敛速度,但其收敛性高度依赖于初值的选取。当下山条件 ∣ f ( x k + 1 ) ∣ < ∣ f ( x k ) ∣ |f(x_{k+1})| < |f(x_k)| f(xk+1)<f(xk) 不满足时,通过减小 λ \lambda λ 来调整步长,可以保证收敛。
::::
:::: details 🔄 举一反三
若使用弦截法求解此方程,其迭代公式(双点形式)是什么?
::: details 查看答案与解析
答案: x k + 1 = x k − f ( x k ) ( x k − x k − 1 ) f ( x k ) − f ( x k − 1 ) x_{k+1} = x_k - \frac{f(x_k)(x_k - x_{k-1})}{f(x_k) - f(x_{k-1})} xk+1=xkf(xk)f(xk1)f(xk)(xkxk1)
解析: 弦截法用两点连线的斜率代替导数,避免了计算 f ′ ( x ) f'(x) f(x)
:::
::::
:::::

  1. l j ( x ) l_j(x) lj(x) ( j = 0 , 1 , … , n ) (j=0,1,\dots,n) (j=0,1,,n) 是区间 [ a , b ] [a, b] [a,b] 上的一组 n n n 次插值基函数,则插值型求积公式中的求积系数 A j = A_j = Aj= ______

::::: details 查看答案与解析
答案: ∫ a b l j ( x ) d x \int_a^b l_j(x) dx ablj(x)dx

详细解析:
在插值型求积公式中,我们将被积函数 f ( x ) f(x) f(x) 用其 n n n 次插值多项式 L n ( x ) = ∑ j = 0 n f ( x j ) l j ( x ) L_n(x) = \sum_{j=0}^n f(x_j)l_j(x) Ln(x)=j=0nf(xj)lj(x) 代替。

∫ a b f ( x ) d x ≈ ∫ a b L n ( x ) d x = ∑ j = 0 n f ( x j ) ( ∫ a b l j ( x ) d x ) \int_a^b f(x) dx \approx \int_a^b L_n(x) dx = \sum_{j=0}^n f(x_j) \left( \int_a^b l_j(x) dx \right) abf(x)dxabLn(x)dx=j=0nf(xj)(ablj(x)dx)

A j = ∫ a b l j ( x ) d x A_j = \int_a^b l_j(x) dx Aj=ablj(x)dx,即得到求积公式 ∑ j = 0 n A j f ( x j ) \sum_{j=0}^n A_j f(x_j) j=0nAjf(xj)

难度:
考点: #插值型求积公式 #求积系数
:::: details 💡 学习锦囊
📖 核心知识点:

  • 求积系数 A j A_j Aj 只与节点 x j x_j xj 的选取有关,与被积函数 f ( x ) f(x) f(x) 无关。
  • 牛顿-科特斯公式:节点等距分布时的插值型求积公式。
    ::::
    :::: details 🔄 举一反三
    对于区间 [ − 1 , 1 ] [-1, 1] [1,1] 上的 2 节点插值型求积公式(节点为 − 1 -1 1 1 1 1),其求积系数 A 0 , A 1 A_0, A_1 A0,A1 是多少?
    ::: details 查看答案与解析
    答案: A 0 = 1 , A 1 = 1 A_0 = 1, A_1 = 1 A0=1,A1=1
    解析: 区间长度为 2,等权分配。
    :::
    ::::
    :::::
  1. 当插值节点为等距分布时,如果要估计结果的舍入误差,应该选用插值公式中的 ______

::::: details 查看答案与解析
答案: 差分形式(或牛顿插值公式)

详细解析:
在等距节点下,使用差分形式(如牛顿前向/后向插值公式)可以更直观地观察到各阶差分的量级。当某阶差分由于舍入误差开始出现不规则波动(噪声)时,可以据此判定插值阶数的选取上限,从而有效地估计和控制舍入误差的影响。

难度: ⭐⭐
考点: #舍入误差 #等距节点插值
:::: details 💡 学习锦囊
📖 计算稳定性提示: 高次插值(如 Lagrange 插值)容易出现 Runge 现象,导致边缘误差极大。在实际应用中,通常采用分段低次插值来平衡计算复杂度和误差控制。
::::
:::: details 🔄 举一反三
为了抑制 Runge 现象,应该如何选择插值节点?
::: details 查看答案与解析
答案: 使用 Chebyshev 节点(切比雪夫节点)。
解析: Chebyshev 节点在区间两端分布较密,中间较稀,能使插值余项的最大值达到最小。
:::
::::
:::::

  1. 解方程 f ( x ) = 0 f(x) = 0 f(x)=0 的简单迭代法 x k + 1 = φ ( x k ) x_{k+1} = \varphi(x_k) xk+1=φ(xk) 满足在有根区间内 ______,则在有根区间内任意取一点作为初始值,迭代解都收敛。

::::: details 查看答案与解析
答案: ∣ φ ′ ( x ) ∣ ⩽ L < 1 |\varphi'(x)| \leqslant L < 1 φ(x)L<1

详细解析:
根据局部收敛性定理(或压缩映像原理),若迭代函数 φ ( x ) \varphi(x) φ(x) 在有根区间 [ a , b ] [a, b] [a,b] 内满足:

  1. φ ( x ) ∈ [ a , b ] \varphi(x) \in [a, b] φ(x)[a,b](映射自身);
  2. φ ( x ) \varphi(x) φ(x) 可导且满足 ∣ φ ′ ( x ) ∣ ⩽ L < 1 |\varphi'(x)| \leqslant L < 1 φ(x)L<1(压缩性)。
    则该迭代过程对区间内任意初值均收敛。

难度:
考点: #简单迭代法 #收敛性条件
:::: details 💡 学习锦囊
📖 收敛速度:

  • 0 < ∣ φ ′ ( α ) ∣ < 1 0 < |\varphi'(\alpha)| < 1 0<φ(α)<1,则为线性收敛
  • φ ′ ( α ) = 0 \varphi'(\alpha) = 0 φ(α)=0,则至少为平方收敛(如牛顿法)。
    ::::
    :::: details 🔄 举一反三
    如果迭代序列收敛较慢,可以采用哪种方法进行加速?
    ::: details 查看答案与解析
    答案: Aitken 加速法(埃特金加速)。
    解析: 通过利用序列的前几项预测极限位置,大幅提高收敛阶。
    :::
    ::::
    :::::
  1. n n n 个求积节点的高斯求积公式的代数精度为 ______

::::: details 查看答案与解析
答案: 2 n − 1 2n-1 2n1

详细解析:
高斯求积公式通过最优地选择求积节点 x j x_j xj,使得对于 n n n 个节点,其代数精度达到最高。理论证明,选取 n n n 次正交多项式的零点作为节点时,公式对所有次数不超过 2 n − 1 2n-1 2n1 的多项式均能精确成立。

难度:
考点: #高斯求积公式 #代数精度
:::: details 💡 学习锦囊
📖 常见高斯求积:

  • Gauss-Legendre 公式:最常用,节点是 Legendre 多项式的根。
  • Gauss-Chebyshev 公式:用于处理带权函数的积分。
    ::::
    :::: details 🔄 举一反三
    3 个节点的高斯求积公式代数精度是多少?
    ::: details 查看答案与解析
    答案: 5。利用公式 2 ( 3 ) − 1 = 5 2(3) - 1 = 5 2(3)1=5
    :::
    ::::
    :::::
Logo

openEuler 是由开放原子开源基金会孵化的全场景开源操作系统项目,面向数字基础设施四大核心场景(服务器、云计算、边缘计算、嵌入式),全面支持 ARM、x86、RISC-V、loongArch、PowerPC、SW-64 等多样性计算架构

更多推荐