Skip to content

非线性方程求根问题

Root-finding problem · Nonlinear equation solving

定义标量与向量求根问题,并区分根的存在、隔离、重数、残差和局部敏感性。

条目类型
模型

形式陈述

标量求根问题给定函数 f:DRR,寻找 rD 使

f(r)=0.

向量问题给定 F:DRnRn,寻找 F(r)=0。输出规格必须说明是寻找一个根、给定区域内的全部根,还是只判断根是否存在;这三项任务需要的信息和算法并不相同。

若存在邻域使 r 是其中唯一的零点,称 r 为孤立根。标量函数在 r 附近若可写成

f(x)=(xr)mg(x),g(r)0,

r 是重数为 m 的根;m=1 时为简单根。对充分可微的函数,这等价于 f(r)==f(m1)(r)=0f(m)(r)0。向量问题通常用 Jacobian 是否可逆描述非退化根,而不直接照搬标量重数。

区间 [a,b] 若端点函数值异号,在连续性假设下构成括根区间,保证其中至少有一个实根。隔离区间的要求更强:它包含指定集合中的恰好一个根。异号本身不证明唯一,也无法检测两侧同号的偶重根;“括住一个根”与“隔离出唯一根”必须分开表述。

近似值 r^ 的残差是 f(r^)F(r^),前向误差则是 r^r。在标量简单根附近,若函数受到小的加性扰动 Δf,一阶关系为

ΔrΔf(r)f(r).

因此 1/|f(r)| 给出一种绝对局部敏感性尺度。向量情形由 JF(r)1 承担同一角色。相对条件数还取决于输入如何参数化、函数如何缩放,以及允许扰动函数值、系数还是其他数据,不能只给一个无上下文数字。

直觉

求根先是一项关于对象的任务,算法随后才决定怎样搜索。存在性回答“有没有交点”,隔离回答“这段区域内有几个”,条件性回答“曲线稍微移动后交点会走多远”。二分、Newton 和割线法分别利用不同信息,不能用某个方法的更新公式代替这些问题规格。

曲线陡直穿过横轴时,小的竖直扰动通常只让交点水平移动一点;曲线贴着横轴相切时,同样的竖直扰动可能让根分裂、消失或移动很远。导数是否远离零因而既影响局部敏感性,也会影响许多迭代法的收敛行为。

例子与边界

f(x)=cosxx.

[0,1] 上有 f(0)=1>0f(1)=cos11<0,由介值定理至少存在一根。又因

f(x)=sinx1<0

在该区间成立,函数严格递减,所以根唯一;其数值约为 0.7390851332。存在性来自连续与异号,唯一性来自单调性,这两个论证不能合并成“二分找到了唯一根”。

函数 (x1)2x=1 有二重根,根两侧函数值同号,符号变化法无法发现它。扰动为

fε(x)=(x1)2ε

时,ε>0 产生两个根 1±ε;若改成加上正的 ε,实根则消失。根的位移按 ε 而不是 ε 缩放,显示重根附近的病态性。

小残差也不能脱离尺度解释。对 f(x)=α(x1),近似值 r^=0 的根误差恒为 1,残差却是 α,可通过缩小 α 变得任意小。只有结合导数下界、逆 Jacobian 范数或其他条件性信息,残差才能转换成根误差界。

多项式求根只是求根问题的一类。寻找一个已括住的实根可以依靠实区间与符号;寻找多项式的全部复根则要处理复共轭结构、聚集根和系数扰动,不能把同一括区间算法机械推广过去。连续函数还可能有一整段零点,此时“输出一个根”和“描述解集”也完全不同。

推论与应用

求根页沿用数值问题与数值算法的接口区分:问题指定 f、搜索域和根的含义,算法只能通过函数值、导数或 Jacobian 等允许信息产生近似。停止状态和失败状态属于算法合同,不能被“存在一个根”的数学命题替代。

若方程可改写为 x=g(x)不动点迭代xk+1=g(xk) 实现求根;在包含初值的完备区域内,压缩映射条件同时给出存在唯一性和线性收敛。不同的等价改写会产生完全不同的 g 与收敛域,因此“代数上能改写成不动点”并不等于迭代可用。

二分求根法只使用连续性和端点符号来维护全局括区间;Newton 法利用导数或 Jacobian 获得局部速度;割线法用两次函数值近似斜率。混合方法通常保留括区间作为安全外壳,再在条件允许时采用更快的局部步。

残差与误差估计提供停止证书的共同语言,问题条件性解释为何同样大小的残差在简单根和重根附近代表不同误差。任何求根结果都应同时报告搜索区域、残差定义、停止条件和失败状态,而不是只给若干小数位。

参考资料
  • NIST Digital Library of Mathematical Functions, §3.8: Nonlinear Equations, accessed 2026.
  • Richard L. Burden, J. Douglas Faires, and Annette M. Burden, Numerical Analysis, 10th ed., Cengage, 2016, Ch. 2.
  • Nicholas J. Higham, Accuracy and Stability of Numerical Algorithms, 2nd ed., SIAM, 2002, Ch. 1.
关系图谱5 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系