锂离子电池SOC估计：EKF与UKF算法实践

1. 项目背景与核心价值

锂离子电池作为新能源汽车的核心部件，其荷电状态（State of Charge, SOC）的准确估计直接关系到整车能量管理效率和电池寿命。传统安时积分法受累积误差影响，而开路电压法需要长时间静置，均难以满足动态工况需求。扩展卡尔曼滤波（EKF）和无迹卡尔曼滤波（UKF）通过融合电池模型与实时测量数据，成为解决这一行业痛点的主流方案。

我在某车企BMS开发项目中实测发现，当SOC估计误差超过5%时，快充效率会下降12%-18%，且电池组温差增大3℃以上。这正是本研究采用Simulink搭建高精度仿真环境的价值所在——在实车测试前验证算法鲁棒性，避免因SOC误判导致的过充/过放风险。

2. 电池建模与参数辨识

2.1 二阶RC等效电路模型构建

选用二阶RC模型（如图1）在复杂度与精度间取得平衡：

code复制Uoc(SOC) ──R0──┬──R1──C1──┬──R2──C2──
                │          │
              电流I      电压U

其中：

• R0：欧姆内阻（主要影响瞬时压降）
• R1/C1：电化学极化环节（时间常数约10-100s）
• R2/C2：浓度极化环节（时间常数约100-1000s）

提示：模型阶数选择需权衡——一阶模型在动态工况下电压误差可达50mV，而三阶模型参数辨识难度指数级增长。

2.2 混合脉冲功率特性（HPPC）测试

通过以下步骤获取模型参数：

1. 在25℃环境舱中以1C倍率完成充放电循环
2. 每10%SOC间隔执行：
   • 静置1h获取开路电压Uoc
   • 施加10s脉冲电流（充/放各一次）
   • 记录弛豫阶段电压响应曲线
3. 采用最小二乘法拟合R0、R1、C1、R2、C2参数矩阵

实测某NMC532电池数据示例：

3. 卡尔曼滤波算法实现

3.1 EKF设计要点

状态方程离散化处理：

code复制X_k = [SOC_k, U1_k, U2_k]^T
U1_k = exp(-Δt/(R1*C1)) * U1_{k-1} + R1*(1-exp(-Δt/(R1*C1)))*I_{k-1}
U2_k = exp(-Δt/(R2*C2)) * U2_{k-1} + R2*(1-exp(-Δt/(R2*C2)))*I_{k-1}

雅可比矩阵计算关键项：

code复制∂f/∂SOC = [1, 0, 0]  // SOC传递只与电流相关
∂h/∂SOC = dUoc/dSOC  // 开路电压-SOC曲线斜率

3.2 UKF改进策略

采用比例修正对称采样（Scaled Unscented Transformation）：

1. 选取sigma点：
   χ_i = x̂ ± (√((n+λ)P_x))_i, i=1,...,2n
   λ = α²(n+κ) - n （通常α=1e-3, κ=0）
2. 权值计算：
   W_m^0 = λ/(n+λ)
   W_c^0 = W_m^0 + (1-α²+β) （β=2最优）
3. 实测对比：在SOC突变工况下，UKF的收敛速度比EKF快约30%

4. Simulink仿真架构

4.1 模块化设计

code复制[电池模型] → [EKF/UKF estimator] → [SOC显示]
    ↑               ↑
[电流传感器]    [电压传感器]

4.2 关键参数配置

1. 求解器选择ode23tb（适合刚性系统）
2. 固定步长设为0.1s（对应BMS实际采样周期）
3. 噪声设置：
   • 过程噪声Q = diag([1e-4, 1e-6, 1e-6])
   • 观测噪声R = 1e-4

5. 充电控制策略优化

5.1 多阶段CC-CV充电改进

1. 恒流阶段：
   • 根据SOC估计值动态调整电流：
     I_charge = min(I_max, Kp*(SOC_target-SOC_est))
   • Kp参数自适应：温度>40℃时降低30%
2. 恒压阶段：
   • 截止条件改为dI/dt < 0.01C/min

5.2 温度补偿策略

建立三维查表：

• 输入：SOC、T、I
• 输出：最大允许充电电流修正系数
  实测效果：在-10℃环境下可避免析锂，充电时间仅延长15%

6. 实测问题与解决方案

6.1 电压测量漂移

现象：持续充电时UKF估计SOC突然跳变5%以上
排查：发现传感器地线阻抗导致共模干扰
解决：增加RC低通滤波（fc=10Hz）并重新校准

6.2 模型失配处理

当电池老化时（容量衰减>20%）：

1. 在线更新R0参数：R0_new = ΔU/ΔI|t=0+
2. 容量重标定：当SOC从100%→5%时，Q_actual = ∫I dt / 0.95

7. 仿真与实车数据对比

UDDS工况测试结果：

在快充场景下（30%-80% SOC）：

• EKF估计误差导致充电电流偏低14%
• UKF方案充电时间缩短至22分钟（比基准快18%）