“区间Newton法将中心残差与整个盒的导数范围组合,先证明修正不会丢掉原盒内的根,再检查排除或存在唯一性。Krawczyk算子改用固定预条件矩阵,严格内包含可同时认证可逆性和唯一根;仅仅算出…”
形式陈述
先定义集合值修正
设
实际算法不必精确求出这个一般并非盒的集合。它调用区间线性系统的包围,寻找包含全部解
- 保留:
中每个根都属于 ,所以取交可以安全收缩候选盒。 - 排除: 若
,则 内无根。 - 存在与唯一: 若
,则 中恰有一个根,并且它属于 。
最后一项在这里允许非严格包含,因为区间Jacobian的正则性已经单独作为前提验证。若换成别的算子或省去正则性,不能自动保留这个边界结论。
一维可直接计算
对
分子是中心点处的数值,不是整个
直觉
普通Newton在一个点画一条切线。区间Newton同时允许盒内出现的全部斜率或线性化矩阵,再把这些可能线性方程的零点包起来。只要真根位于原盒,它对应的某个平均线性化也在允许范围内,所以真根不能被修正集合丢掉。
这解释了“保留”与“存在”的区别。把一个空候选区域收缩得很小仍可能没有根;只有修正集合整个留在原盒中,配合可逆性,才得到连续自映射,进而获得存在性。
证明:平均矩阵不是同一点的Jacobian
对
逐项区间是凸的,故平均矩阵仍在其中。多维时一般不存在一个共同的
若
剩下的是存在性。平均矩阵随
是连续映射。当
例子与边界
平方根的第一张证书
令
这一步既证明原区间恰一根,也把它包进宽度
读者可分别平方两端,核对它们仍严格夹住
如果将普通Newton的表达式直接按区间代入,算成
没有通过时能说什么
取
它与
盒中每个实际Jacobian都可逆,比整个区间矩阵正则弱。区间会允许从不同点独立拼接各个元素,可能包含实际从未出现的奇异矩阵。此时应缩盒、保留相关性或改用别的证明,不能只抽样几个Jacobian就宣称前提已满足。
推论与应用
二次收缩需要哪些额外控制
在一维简单根
这个估计假定中心值和区间端点按精确实数计算。实际可靠外包围还会有舍入或函数求值尾界的宽度,固定精度下可能到达平台,不能无限宣称宽度继续平方。若只知道
标量一步只需中心函数值、导数区间和常数次区间操作。多维的主要成本是线性解集的可靠包围,取决于所用子算法;不能把集合值
参考资料
- R. Baker Kearfott, “Interval Newton Methods”, Encyclopedia of Optimization, 2001,作者稿pp. 1–2,标量公式、多维线性修正接口与存在/唯一性。本页算例重新使用精确分数计算。
- Siegfried M. Rump, “Verification Methods: Rigorous Results Using Floating-Point Arithmetic”, 2010,作者修订稿§13、Theorems 13.1–13.2;平均矩阵、正则前提及闭包含的Brouwer证明。