形式陈述
在实数系 理路 实数系 Real number system · Ordered complete field 满足序域公理与上确界完备性的数系。 上,标量求根问题给定函数 理路 函数 Function · Map · Mapping 由定义域、陪域和单值图共同组成,并把每个输入送到唯一输出的映射。 f : D ⊆ R → R ,寻找 r ∈ D 使
f ( r ) = 0. 向量问题给定 F : D ⊆ R n → R n ,寻找 F ( r ) = 0 。输出规格必须说明是寻找一个根、给定区域内的全部根,还是只判断根是否存在;这三项任务需要的信息和算法并不相同。
若存在邻域使 r 是其中唯一的零点,称 r 为孤立根。对整数 m ≥ 1 ,标量函数在 r 的某邻域内若可写成
f ( x ) = ( x − r ) m g ( x ) , g ( r ) ≠ 0 , 且 g 在 r 连续,则 r 是重数为 m 的根;m = 1 时为简单根。对在根附近属于 C m 的函数,这等价于 f ( r ) = ⋯ = f ( m − 1 ) ( r ) = 0 且 f ( m ) ( r ) ≠ 0 。g 的连续性不可省略:若允许 g ( x ) = 1 / x (x ≠ 0 )并另设 g ( 0 ) = 1 ,就会把 f ( x ) = x = x 2 g ( x ) 的简单根错误地叫作二重根。向量问题通常用 Jacobian 是否可逆描述非退化根,而不直接照搬标量重数。
区间 [ a , b ] 若端点函数值异号,在连续性假设下构成括根区间,保证其中至少有一个实根。隔离区间的要求更强:它包含指定集合中的恰好一个根。异号本身不证明唯一,也无法检测两侧同号的偶重根;“括住一个根”与“隔离出唯一根”必须分开表述。
近似值 r ^ 的残差是 f ( r ^ ) 或 F ( r ^ ) ,前向误差则是 r ^ − r 。在标量情形,设 f ( r ) = 0 、f 在 r 处可微且 f ′ ( r ) ≠ 0 ,并固定在 r 附近连续的 g 。对扰动族 f ε = f + ε g ,若选取满足 f ε ( r ε ) = 0 且 r ε → r 的局部根分支,则
Δ r := r ε − r = − ε g ( r ) f ′ ( r ) + o ( | ε | ) = − Δ f ( r ) f ′ ( r ) + o ( | ε | ) . 这里 Δ f = ε g ,估计描述这条趋向 r 的根分支,不另行保证任意扰动下根的存在或唯一。
因此 1 / | f ′ ( r ) | 给出一种绝对局部敏感性尺度。向量情形由 J F ( r ) − 1 承担同一角色。相对条件数还取决于输入如何参数化、函数如何缩放,以及允许扰动函数值、系数还是其他数据,不能只给一个无上下文数字。
直觉
求根先是一项关于对象的任务,算法随后才决定怎样搜索。存在性回答“有没有交点”,隔离回答“这段区域内有几个”,条件性回答“曲线稍微移动后交点会走多远”。二分、Newton 和割线法 理路 割线法 Secant method 用最近两点的割线斜率近似导数,在无导数的一维求根中取得局部超线性收敛。 分别利用不同信息,不能用某个方法的更新公式代替这些问题规格。
曲线陡直穿过横轴时,小的竖直扰动通常只让交点水平移动一点;曲线贴着横轴相切时,同样的竖直扰动可能让根分裂、消失或移动很远。导数是否远离零因而既影响局部敏感性,也会影响许多迭代法的收敛行为。
例子与边界
令
f ( x ) = cos x − x . 在 [ 0 , 1 ] 上有 f ( 0 ) = 1 > 0 、f ( 1 ) = cos 1 − 1 < 0 ,由介值定理 理路 介值定理 Intermediate value theorem 连续函数在实区间上不能跳过中间高度;由完备性证明存在性,并厘清二分法、唯一性和不动点迭代。 至少存在一根。又因
f ′ ( x ) = − sin x − 1 < 0 在该区间成立,函数严格递减,所以根唯一;其数值约为 0.7390851332 。存在性来自连续与异号,唯一性来自单调性,这两个论证不能合并成“二分找到了唯一根”。
函数 ( x − 1 ) 2 在 x = 1 有二重根,根两侧函数值同号,符号变化法无法发现它。扰动为
f ε ( x ) = ( x − 1 ) 2 − ε 时,ε > 0 产生两个根 1 ± ε ;若改成加上正的 ε ,实根则消失。根的位移按 ε 而不是 ε 缩放,显示重根附近的病态性。
小残差也不能脱离尺度解释。对 f ( x ) = α ( x − 1 ) ,近似值 r ^ = 0 的根误差恒为 1 ,残差却是 − α ,可通过缩小 α 变得任意小。只有结合导数下界、逆 Jacobian 范数或其他条件性信息,残差才能转换成根误差界。
多项式求根只是求根问题的一类。寻找一个已括住的实根可以依靠实区间与符号;寻找多项式的全部复根则要处理复共轭结构、聚集根和系数扰动,不能把同一括区间算法机械推广过去。连续函数还可能有一整段零点,此时“输出一个根”和“描述解集”也完全不同。
两种复根证书把这个区别落实为交付物:完整边界计数 理路 整条复边界的绕数证书 Certified complex winding · Validated polygon winding · 复零点围道计数证书 以避零凸圆盘覆盖真实像路径,认证采样折线的整数绕数;给出精确射线交叉规则、完整边分割复核、有限停止条件与预算未决出口。 提交全部源边小段与准确绕数,系数族根簇 理路 Pellet 圆盘计数与根簇证书 Pellet root cluster certificate · Pellet Tk test · Coefficient-family root clustering 通过平移系数的单项优势认证圆盘根总重数,传递系数误差预算,并以不相交圆盘和次数闭合交付全部复根簇而非伪称不同根隔离。 提交局部优势、系数误差余量、不相交性和次数闭合。它们可以认证总重数与覆盖;若边界碰根或误差使单根身份不可识别,仍应明确返回未决或根簇。
推论与应用
求根页沿用数值问题与数值算法 理路 数值问题与数值算法 Numerical problem and algorithm 区分数学问题、有限数据、求解算法与实际执行,并据此追踪误差和计算成本。 的接口区分:问题指定 f、搜索域和根的含义,算法只能通过函数值、导数或 Jacobian 等允许信息产生近似。停止状态和失败状态属于算法合同,不能被“存在一个根”的数学命题替代。
若方程可改写为 x = g ( x ) ,不动点迭代 理路 不动点迭代 Fixed-point iteration · Picard iteration 把方程改写为不动点问题后反复应用同一映射,并用压缩性控制收敛、误差与停机。 用 x k + 1 = g ( x k ) 实现求根;在包含初值的完备区域内,压缩映射条件同时给出存在唯一性和线性收敛。不同的等价改写会产生完全不同的 g 与收敛域,因此“代数上能改写成不动点”并不等于迭代可用。
二分求根法 理路 二分求根法 Bisection method 以端点异号区间为不变量,给出可验证误差界、对数成本和有限精度停止条件。 只使用连续性和端点符号来维护全局括区间;Newton 法 理路 Newton 非线性方程法 Newton's method · Newton-Raphson method · Newton method for nonlinear systems 在当前点解一阶线性化方程来修正非线性方程的近似解,并分析其局部二次收敛与失效边界。 利用导数或 Jacobian 获得局部速度;割线法用两次函数值近似斜率。混合方法通常保留括区间作为安全外壳,再在条件允许时采用更快的局部步。
残差与误差估计 理路 残差、误差估计与停止准则 Residual and error estimation · Stopping criterion 区分可计算残差与未知真误差,并说明把缺陷转成误差界和停止证书所需的条件。 提供停止证书的共同语言,问题条件性 理路 问题条件性与条件数 Conditioning of a problem · Condition number 度量问题真解对输入微扰的局部敏感性,并把这种固有敏感性与算法误差分开。 解释为何同样大小的残差在简单根和重根附近代表不同误差。任何求根结果都应同时报告搜索区域、残差定义、停止条件和失败状态,而不是只给若干小数位。
参考资料
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.