研究生

研究生数值分析方程的数值解


原始发现


寻找根的问题是数值分析中的常见挑战,其目标是确定一个给定函数何时变为零的值。对于数学表达为 ( f(x) = 0 ) 的函数求根是一个基本任务,在科学计算、工程、物理和数学中有许多应用。在这次全面讨论中,我们将探讨用于寻找根的各种数值技术和算法,并提供可视化示例和详细解释。

求根介绍

在数学分析中,一个函数的根或零点是满足 ( f(x) = 0 ) 的数 ( x )。在处理低次多项式时求根有时是简单的,因为存在用于直接解决线性、二次、三次和四次方程的公式。然而,对于高次多项式或非多项式函数,通常不存在解析解,必须使用数值方法。

求根重要性的故事

求根在各个领域中扮演重要角色。在物理学中,根可以表示平衡点。在工程学中,求根可以意味着评估系统在特定条件下的行为。在计算机图形学中,当计算交点时及在优化问题中用方程表示约束时,它显得尤为重要。因此,学习各种求根技术成为必要的知识。

包络法

包络法是一类用于求根的流行数值方法。这些方法依赖于中值定理,该定理指出如果一个连续函数在一个区间上改变符号,则该区间内至少有一个根。这些方法从两个初始猜测 (a, b) 开始,满足 ( f(a) cdot f(b) < 0 )。

二分法

二分法是最简单且最可靠的包络法之一。它通过将包含根的区间反复分成两半,选择符号变化发生的子区间来缩小包含根的区间。让我们通过一个示例和可视化来探索这个方法。

考虑函数 ( f(x) = x^2 - 4 )。我们想要找到根,开始时的初始猜测为 ( a = 0 ) 和 ( b = 3 ),其中 ( f(a) = -4 ) 和 ( f(b) = 5 )。

初始区间: [0, 3]

迭代1:
  中点: c = (0 + 3) / 2 = 1.5
  计算: f(c) = 1.5^2 - 4 = -1.75
  新区间: [1.5, 3] (因为 f(1.5) < 0)

迭代2:
  中点: c = (1.5 + 3) / 2 = 2.25
  计算: f(c) = 2.25^2 - 4 = 1.0625
  新区间: [1.5, 2.25] 

继续这一过程,直到间隙足够小。
0 1.5 (中点) 3

二分法的优势在于其简单性和在根附近保证收敛的能力,尽管收敛速度相对较慢且为线性。

开放法

与包络法不同,开放法不需要初始点来包络根,这使得它们在某些情况下更快,但它们可能不太可靠且并不总是收敛。

牛顿–拉弗森法

牛顿-拉弗森法是一种开放法,利用导数来递次逼近实值函数的根。该方法使用迭代点的切线更好地逼近根。

对于一个函数 ( f(x) ),如果 ( x_n ) 是一个根的估计值,则下一个估计值 ( x_{n+1} ) 由以下公式给出:

x_{n+1} = x_n - frac{f(x_n)}{f'(x_n)}

示例

假设我们从初始猜测 ( x_0 = 1 ) 开始,找出 ( f(x) = x^2 - 2 ) 的根。

f(x) = x^2 - 2
f'(x) = 2x

迭代1:
  x_1 = 1 - frac{1^2 - 2}{2(1)} = 1.5

迭代2:
  x_2 = 1.5 - frac{1.5^2 - 2}{2(1.5)} = 1.41666667

继续,直到 x 接近真实根。

这种方法在根附近以二次速度收敛,比二分法更快,因为它是以二次速度收敛的。

割线法

割线法可以被视为牛顿–拉弗森法的有限差分近似,不需要计算导数。这种方法使用两个初始猜测,通过以下方式迭代计算割线:

x_{n+1} = x_n - f(x_n) cdot frac{x_n - x_{n-1}}{f(x_n) - f(x_{n-1})}

示例

假设我们想从 ( x_0 = 1 ) 和 ( x_1 = 2 ) 开始,找到 ( f(x) = x^2 - 3 ) 的根。

迭代1:
  x_2 = 2 - frac{4 - 3}{2 - 1} = 1.5

迭代2:
  x_3 = 1.5 - frac{2.25 - 3}{1.5 - 2} = 1.75

这个模式继续并到达根。
x_0 = 1 x_1 = 2

不动点迭代

不动点迭代使用形式为 ( x = g(x) ) 的重排函数,并且迭代进行 ( x_{n+1} = g(x_n) )。这方法包括重新编写函数,使得对该函数的迭代收敛到根。

示例:考虑找到 ( f(x) = x^3 - x - 2 = 0 ) 的根。

重排: x = g(x) = (x + 2)^{1/3}

要重复:
  x_0 = 初始猜测,假设为 1.5
  x_1 = g(x_0) = (1.5 + 2)^{1/3} ≈ 1.650
  x_2 = g(x_1) = (1.650 + 2)^{1/3} ≈ 1.553

继续迭代,直到收敛。

图形示例

g(x) y=x x_0 x_1 x_2

在实践中,收敛依赖于 ( g(x) ) 的选择和初始猜测 ( x_0 )。

结论

求根是数值分析的重要方面,涉及找到表达为 ( f(x) = 0 ) 函数的根。根据函数的特性、领域知识和所需精度,可能有更适合特定问题的方法。包络法和开放法甚至图形解释之间的选择很大程度上取决于导数的可用性以及初始猜测与真实根的接近程度。理解和正确应用这些数值方法是解决所有科学和工程学科中实际数学模型的关键部分。

进一步的学习和通过在计算软件或编程环境中实现这些方法的实践,有助于增强理论理解与实际技能。了解每个方法的优点、局限性和收敛特性对于它们的有效应用很重要。


研究生 → 6.1.1


U
username
0%
完成于 研究生


评论