4. OCV-SOC曲线分段线性插值实现与精度评估

任务拆解与数据准备

本案例要解决的核心问题是：给定一张离散的 OCV-SOC 关系表（例如每 5% SOC 一个节点），如何用分段线性插值把任意电压测量值映射为 SOC，并量化“节点数量”对估算精度的影响。船舶电池的 OCV-SOC 曲线在中间平台区（约 30%~70% SOC）非常平坦，电压微小偏差会放大为较大的 SOC 误差，因此插值节点密度是精度控制的关键参数。

我们先用一个可复现的解析模型生成“高密度真值表”，再从中抽取不同密度的节点做插值，最后用同一批测试电压对比误差。为便于教学，采用如下单调 OCV-SOC 关系（单位：V，SOC 取 0~1）：

OCV(s) = 3.20 + 0.55·s + 0.12·sin(π·s)   // 平台区被正弦项压平，模拟磷酸铁锂特征

该函数在 s∈[0,1] 上单调递增，可作为“真值”。高密度查表用 1001 个节点（步长 0.001）作为精度基准。

分段线性插值算法实现

插值函数需要完成三件事：定位电压落在哪个区间、在该区间内做线性内插、处理边界与越界。下面是完整可运行的 Python 实现，包含真值生成、节点抽取、插值与误差评估：

import numpy as np

def ocv_true(s):
    return 3.20 + 0.55 * s + 0.12 * np.sin(np.pi * s)

def build_table(n_nodes):
    """在 SOC 0~1 上等间隔取 n_nodes 个节点，返回 (soc_nodes, ocv_nodes)"""
    soc = np.linspace(0.0, 1.0, n_nodes)
    return soc, ocv_true(soc)

def interp_soc(v_meas, soc_nodes, ocv_nodes):
    """分段线性插值：输入电压，输出 SOC；越界时钳位到端点"""
    v = np.asarray(v_meas, dtype=float)
    soc = np.interp(v, ocv_nodes, soc_nodes)  # ocv_nodes 单调递增
    return soc

# 高密度基准表
soc_ref, ocv_ref = build_table(1001)

# 测试电压：覆盖平台区与两端，共 401 个点
v_test = np.linspace(ocv_ref[0], ocv_ref[-1], 401)
soc_exact = interp_soc(v_test, soc_ref, ocv_ref)  # 基准 SOC

for n in [6, 11, 21, 41, 81, 161]:
    s_n, v_n = build_table(n)
    soc_est = interp_soc(v_test, s_n, v_n)
    err = soc_est - soc_exact
    print(f"节点数={n:4d}  最大|误差|={np.max(np.abs(err))*100:6.3f}%  "
          f"RMS误差={np.sqrt(np.mean(err**2))*100:6.3f}%")

关键点说明：np.interp 要求自变量（此处为 OCV）单调递增，本曲线满足；越界电压会被钳位到首末 SOC，符合“电压低于最低 OCV 视为 0% SOC”的工程约定。测试电压用 401 点均匀覆盖全量程，能同时暴露平台区与陡峭区的误差差异。

节点数量对精度的影响

运行上述代码可得到一组典型结果（示例数据，由给定解析模型直接算出，非实测）：

误差随节点数增加呈近似二次收敛（分段线性插值误差与区间长度平方成正比），但节点数翻倍带来的精度提升迅速衰减。对船舶 BMS 而言，5% 步长（21 节点）通常已能把插值误差压到 0.3% 以内，远小于电压采样噪声（典型 ±5 mV 对应平台区约 1%~2% SOC 误差）的影响，因此继续加密节点的边际价值有限。

插值误差的空间分布

下图展示了同一组节点（21 节点）下，插值 SOC 与基准 SOC 的偏差沿 SOC 轴的分布形态，以及节点加密对最大误差的收敛趋势。可以直观看到误差在平台区中部最大、在两端趋近于零，这正是分段线性插值“弦线逼近曲线”的固有特征。

结果判读与节点数选择

判读要点有三条：

• 看最大误差而非平均误差：平台区中部的峰值误差决定了最坏情况下的 SOC 偏差，是选型依据。
• 与电压噪声对比：若电压采样噪声在平台区折算的 SOC 不确定度已大于插值误差，则加密节点无意义。可用 dSOC/dV 估算：平台区斜率约 0.3 V/100% SOC，5 mV 噪声对应约 1.7% SOC。
• 看收敛拐点：从 21 节点到 41 节点误差降约 4 倍，从 41 到 81 降约 4 倍，但绝对量已很小；工程上取 21~41 节点即可。