关于Dekker数值根查找法区间更新逻辑的技术问询

一、对跖点选择逻辑的解释 #

Dekker方法的核心约束是始终维护区间[a, b]，保证两端点的函数值符号相反（即根必在该区间内）。维基百科里的这段描述，本质是在确定新的区间左端点：

新对跖点的选择需满足f(a_k+1)与f(b_k+1)符号相反。若f(a_k)与f(b_k+1)符号相反，则对跖点保持不变：a_k+1 = a_k；否则，f(b_k+1)与f(b_k)符号相反，新对跖点为a_k+1 = b_k。

• 当f(a_k) * f(b_{k+1}) < 0时：说明新生成的点b_{k+1}和原左端点a_k的函数值异号，根必然在[a_k, b_{k+1}]区间内。所以不需要改动左端点，直接保留a_{k+1}=a_k，只需要把右端点更新为b_{k+1}即可。
• 当f(a_k) * f(b_{k+1}) >= 0时：说明b_{k+1}和a_k同号，那根据原区间[a_k, b_k]的约束（f(a_k)*f(b_k)<0），b_{k+1}必然和b_k异号，根就在[b_k, b_{k+1}]区间内。所以要把左端点更新为原右端点b_k，右端点同样更新为b_{k+1}。

你的代码逻辑是对的，但ak = ak属于无意义的冗余代码，可以简化为：

# 先保存旧的bk值，避免更新后被覆盖
bk_old = bk
# 更新右端点为新生成的bk_next
bk = bk_next
# 判断左端点是否需要更新
if fbk_next * func(ak) >= 0:
    ak = bk_old

二、交换步骤的执行时机 #

维基百科提到的交换步骤，是为了让右端点b始终是区间内函数值绝对值更大的那个点，这样后续用割线法生成下一个迭代点时，能更高效地向根收敛。这个步骤的执行时机是：

若|f(a_k+1)| < |f(b_k+1)|，则a_k+1可能是比b_k+1更优的解猜测值，因此需交换a_k+1与b_k+1的值。

具体来说，要在完成对跖点选择（即确定了新的a_{k+1}和b_{k+1}）之后执行交换：

1. 先完成区间端点的更新（确定新的a和b）
2. 计算abs(func(a))和abs(func(b))
3. 如果abs(func(a)) < abs(func(b))，交换a和b的值

代码片段示例：

# 完成区间端点更新后
fa = func(ak)
fb = func(bk)
if abs(fa) < abs(fb):
    ak, bk = bk, ak
    fa, fb = fb, fa

收敛偏差的可能原因 #

你的代码出现轻微收敛偏差，大概率是以下两个问题：

• 没有正确保存原bk值，导致左端点更新时用了已经修改后的bk（即bk_next），破坏了区间的符号约束
• 交换步骤的时机错误（比如在更新区间端点前执行，或者完全漏掉），导致后续迭代的点生成逻辑偏离了Dekker方法的要求

内容的提问来源于stack exchange，提问作者blov