CDF样条插值边界延续异常排查：x>110段未保持水平1反而下降的问题

我完全懂你现在的头疼点——明明想让CDF样条在x>110之后稳稳保持1，结果它反倒掉回0了，这肯定是代码里的拼接逻辑出了问题，咱们来一步步捋清楚：

核心问题分析 #

你当前的代码里有两个关键错误，导致了结尾段的异常：

1. 样条阶数不匹配：中间用splrep拟合的是三次样条（默认k=3），但你手动构造的结尾段是线性的（k=1），直接拼接会破坏三次样条的节点约束规则，导致BSpline无法正确识别边界的固定值要求。
2. 节点与系数拼接逻辑错误：三次样条要实现“某区间恒为常数”的效果，需要在区间边界处设置k+1个重复节点（对于k=3来说就是4个），这样才能让样条的一阶、二阶导数都为0，保证水平延续。你当前的节点拼接没有满足这个要求。

修正后的代码方案 #

下面是调整后的代码，我会逐段解释逻辑：

import numpy as np
from scipy.interpolate import splrep, BSpline

def _find_spline(x_orig: np.ndarray, cdf: np.ndarray) -> BSpline:
    # 保留你原有的平滑因子计算逻辑
    smoothing_factor, iters = _find_highest_smoothing_factor(x_orig, cdf, 1e-5)

    # 中间段：拟合三次样条，保留默认的k=3
    t_middle, c_middle, k_middle = splrep(x_orig, cdf, s=smoothing_factor)
    k = k_middle  # 固定为三次样条的阶数3

    # 处理开头段：x <= x_orig[0]（即x<=40）时恒为0
    # 三次样条需要k+1个重复节点来固定值，在x_orig[0]处添加k个重复节点
    t_begin = np.concatenate([[x_orig[0]] * k, t_middle[:k+1]])
    c_begin = np.concatenate([[0.0] * k, c_middle[:1]])

    # 处理结尾段：x >= x_orig[-1]（即x>=110）时恒为1
    # 同样在x_orig[-1]处添加k个重复节点，确保样条水平延续
    t_end = np.concatenate([t_middle[-k-1:], [x_orig[-1]] * k])
    c_end = np.concatenate([c_middle[-1:], [1.0] * k])

    # 合并节点和系数：去掉中间段与前后段重复的节点，避免冗余
    t_combined = np.concatenate([t_begin[:-k], t_middle[k:-k], t_end[k:]])
    c_combined = np.concatenate([c_begin[:-1], c_middle[1:-1], c_end[1:]])

    # 创建BSpline，关闭自动外插（因为已手动处理边界）
    return BSpline(t_combined, c_combined, k, extrapolate=False)

关键逻辑说明 #

• 开头段处理：在x=40处设置3个重复节点，加上中间段的起始节点，确保三次样条在x<=40时导数全为0，恒为0。
• 结尾段处理：在x=110处设置3个重复节点，加上中间段的结束节点，让三次样条在x>=110时保持水平1。
• 节点拼接：去掉中间段和前后段重复的节点，避免BSpline出现节点冲突，保证样条的连续性。
• 关闭自动外插：因为我们已经手动处理了边界外的情况，不需要依赖BSpline的默认外插逻辑，防止再次出现异常。

用这个代码测试的话，x<40时会恒为0，x>110时恒为1，中间40-110段保持你原来的拟合效果，完全符合你的需求。

备注：内容来源于stack exchange，提问作者lqbweb