拓冰网站定制 
简介锂电池荷电状态SOC估算一直是电池管理系统BMS中的核心难题电池本身是强非线性、时变且与温度强相关的电化学系统传统安时积分法存在累积误差扩展卡尔曼滤波在大非线性场景下线性化误差明显。容积卡尔曼滤波CKF基于三阶球面-径向容积规则无需复杂调参且数值稳定性好适合强非线性状态递推Transformer则凭借自注意力机制捕获长程时序依赖有效补偿等效电路模型在低温、动态工况下的失配误差。将二者融合既保留物理模型的可解释性又借助数据驱动能力提升复杂场景的泛化精度。该方案在BJDST、DST、FUDS、US06四种动态工况及0°C、25°C、45°C三档温度下完成验证SOC估计RMSE大幅降低尤其在低温场景提升显著。这一技术路径可为BMS算法设计与工业时序建模提供参考。 锂离子电池的SOC荷电状态估算是BMS里最让人头疼的模块之一。这个项目做了一套基于Transformer与CKF容积卡尔曼滤波融合的SOC高精度估计系统在BJDST、DST、FUDS、US06四种动态工况和0°C、25°C、45°C三档温度下做了完整验证。整套方案的核心思路是让CKF负责电池模型动力学的状态递推让Transformer去补偿模型失配和复杂工况带来的非线性误差两者各干各擅长的活。对做BMS算法、电池状态估计或者想了解Transformer在工业时序任务里怎么落地的同学来说这套架构值得花几分钟看完——尤其是融合策略和温度泛化这两块我踩了不少坑后面都会讲到。1. 项目背景与核心问题剖析1.1 SOC估计为什么是BMS里最难啃的部分SOCState of Charge荷电状态说白了就是电池还剩多少电通常用百分比表示0%到100%。它不像电压电流那样拿万用表直接量出来必须通过模型或算法去“猜”。这个“猜”字背后藏着整个BMS算法的核心矛盾电池是一个强非线性、时变、且和温度强相关的电化学系统OCV开路电压-SOC曲线在中段非常平缓电压稍微有点测量噪声SOC就会被放大得很离谱。早期BMS里最常用的是安时积分法也就是把电流对时间做积分初始SOC加上电流累积量再除以额定容量。这个方法在实验室里精度还不错但实际用车场景里问题很大一是电流传感器存在零漂和噪声积分误差会越积越大二是电池可用容量会随温度、老化、放电倍率变化你拿一个固定容量去除算出来的SOC自然不准。跑一个完整的充放电循环安时积分误差到5%以上很常见而且这个误差不会自己消失。后来大家开始用开路电压法做校正但OCV必须等电池静置足够长时间才能测准车总不能停在路边等两个小时再去算电量。再后来就是基于等效电路模型ECM加卡尔曼滤波的一套方法从EKF扩展卡尔曼滤波到UKF无迹卡尔曼滤波再到UKF和神经网络结合这条路走了很多年。核心问题是模型再精细也是对真实电化学过程的近似模型误差在剧烈动态工况下会被滤波算法部分“吸收”成状态误差最终反映在SOC估计偏差上。1.2 方案选型Transformer与CKF各承担什么角色这个项目最终选了Transformer和CKF做融合不是拍脑袋决定的而是把SOC估计任务拆解成两个子问题之后自然得出的结论。第一个子问题是“状态递推”。SOC本身是一个状态量它符合连续时间动态系统的规律——上一时刻的SOC和电流、温度基本决定了当前时刻的SOC。卡尔曼滤波这类方法天生就是干这个的它能把模型预测和量测校正结合起来输出一个带方差信息的估计值。在卡尔曼滤波家族里我选了CKF而不是更普及的EKF或UKF原因后面会详细说简单讲就是EKF在大非线性场景下雅可比矩阵线性化误差大UKF的sigma点分布需要调参、数值稳定性在状态维数变高时不太好保证而CKF用三阶球面-径向容积规则没有额外超参数数值稳定性好很适合电池这种中等维数3到5维的强非线性系统。第二个子问题是“模型误差补偿”。不管等效电路模型建得多细温度和老化导致的参数漂移、电池内部的迟滞效应、工况突变带来的动力学变化都会让模型和真实系统之间存在偏差。这种偏差是高度非线性的很难用显式公式描述。这时候就需要一个表达能力强的数据驱动模型来兜底。Transformer的优势在于自注意力机制能直接建立任意两个时间步之间的依赖关系对SOC估计这种长序列任务它比LSTM能更自然地抓住“几分钟前一个大电流脉冲对当前电压的影响”这类长程特征。所以架构就很清晰了CKF做主递推Transformer做误差补偿两者不是竞争关系而是互补关系。后续第四章的实验数据也证明融合后的RMSE比单独用CKF提升明显尤其在0°C低温这种模型参数严重失配的场景下Transformer补偿的效果格外突出。2. 总体架构与方案设计思路2.1 系统架构概览滤波主骨架加数据驱动补偿支路整套系统的数据流是这样的电池在特定工况下放电测试台以1Hz采样率记录端电压、负载电流和表面温度这三个时间序列作为输入经过滑动窗口切片窗口长度256个采样点约4.3分钟分别送入两个支路。主支路是CKF滤波器它使用一个二阶RC等效电路模型后续简称2RC模型作为状态方程。状态向量是三维的SOC、电化学极化电压U1、浓度极化电压U2。输入是电流I量测方程是端电压Ut的预测值。每来一个新采样点CKF执行一次“预测-校正”流程输出当前SOC的滤波估计值和协方差。补偿支路是Transformer编码器。窗口内的电流、电压、温度先做归一化加上位置编码后进入4层编码器。Transformer输出的是两个东西一个是端电压的预测残差ΔU用来修正量测方程让后续的卡尔曼更新更准另一个是SOC的修正量ΔSOC直接叠加到滤波输出上。这两个输出分别接到CKF的两个不同接口等于给滤波算法加了一个“环境感知”的外挂。最终融合后的SOC估计值再通过一个低通滤波器做平滑处理防止卡尔曼增益突变导致输出跳变。整个系统在Python里完成训练和离线验证导出ONNX模型后可以嵌入到BMS原型控制器里做实时推理。实测在单核ARM Cortex-A53上一次完整的模型推理加滤波更新耗时约18毫秒满足实时性要求。2.2 为什么是“TransformerCKF”而不是“LSTMEKF”我在前面已经简单提到了选型方向这里再展开讲讲对比过程因为这是很多人会问的问题。先看数据驱动部分。早期序列建模基本是LSTM和GRU的天下它们在SOC估计里也有不少应用效果比纯EKF好。但LSTM有个天然问题按时间步顺序计算当前时刻的隐状态依赖于上一时刻的隐状态无法并行训练速度慢而且长程信息依赖门控结构一步步传递链长了之后早期信息容易被稀释。Transformer的自注意力机制绕开了这个问题——它把整个序列一次性输入每个位置可以直接和其他所有位置计算注意力权重想用哪一段历史信息用哪一段不依赖循环递推。在SOC估计里这特别适合捕捉“几十秒前那个大电流放电对当前电压极化的影响”。再看滤波部分。EKF需要计算非线性函数的雅可比矩阵在OCV-SOC曲线平台区SOC在30%到70%之间导数很小EKF容易对这个区间的量测校正不敏感SOC估计误差会增大。UKF用无迹变换避免了雅可比计算但需要设定alpha、beta、kappa三个参数这些参数对滤波效果影响大调起来非常玄学。CKF用的容积点都是从协方差矩阵的平方根变换得到的权重全部为正不存在UKF在状态维数较高时sigma点权重为负导致的协方差非正定问题而且不需要调参。从论文复现角度的检验也表明CKF在高维强非线性系统中的数值稳定性最好。一句话总结LSTMEKF是“上一个时代”的成熟方案TransformerCKF是“当前在准确率和工程可行性之间平衡得比较好的方案”。如果你手里有充足的历史数据、计算资源充裕这个选型值得试试。2.3 数据链路与多工况设计为什么要覆盖四种工况和三档温度算法模型不能只在实验室最理想的条件里跑得准工程应用讲究泛化能力。这个项目设计了四组标准动态工况——BJDST、DST、FUDS、US06——来覆盖不同驾驶风格和路况对电池的激励模式。四个工况的差异本质上是电流激励的动态频率和幅值不同。DSTDynamic Stress Test是美国先进电池联合会USABC制定的简化动态工况电流曲线由周期性的放电、充电和静置段组成适合做算法初调。FUDSFederal Urban Driving Schedule模拟的是城市道路上走走停停的驾驶状态电流变化比较频繁但整体幅值不大。US06则是美国环保署定义的高速高功率工况电流变化快、峰值大、动态最剧烈对模型和滤波算法的动态响应能力要求最高。BJDST可以理解为考虑了北京路况特点的改进版DST低速蠕行和频繁启停的比例更高。三档温度的选择也对应真实场景0°C是冬季户外低温环境电池内阻显著升高、极化效应增强最容易出现SOC估计偏差25°C是常温算法表现最稳的状态45°C是高温工况不仅电化学参数变化还涉及热失控风险SOC估计必须保持稳定。在实际验证中我遵循的原则是“低温场景优先保不跑偏高温场景重点看稳定性”。模型训练时把三档温度的数据都喂进去测试时每个温度单独评估这样可以暴露泛化短板。3. 核心模块原理与工程实现细节3.1 Transformer模块设计与实现要点Transformer部分并不是直接拿NLP领域的标准架构照搬而是针对电池时序数据做了几个关键修改。输入特征有三列端电压Ut、负载电流I、温度T。由于SOC动态范围很大直接喂原始值会让模型训练不稳定所以先对每列做独立的Min-Max归一化映射到0到1之间。为了增强模型对局部动态变化的感知我还额外算了两个差分特征dUt当前电压和上一步电压的差值和dI当前电流和上一步电流的差值。这两个差分特征在实验里对精度的提升很有帮助相当于把“变化趋势”显式告诉了模型。位置编码是Transformer处理时序数据时绕不开的一步。由于自注意力机制本身是位置无关的它不关心输入顺序必须把时间步信息嵌入到输入里。这里用标准的正弦位置编码因为它的相对位置信息可以通过旋转矩阵表达模型比较容易学到时间步之间的相对距离概念。模型结构采用4层编码器Encoderd_model设为128注意力头数为4FFN中间维度是512dropout统一用0.1。解码器没有用——原因是SOC估计本质上是回归映射任务我们只需要把输入序列编码成一个特征向量然后接一个全连接输出层直接回归当前时刻的SOC修正量和端电压残差不需要像机器翻译那样逐词生成。这个决策能省不少参数量和推理时间。训练时的损失函数是Huber Loss同时关注SOC修正量和端电压残差两部分用加权方式组合。SOC修正量权重给0.8端电压残差给0.2因为端口电压只是个中间量SOC才是最终目标。优化器用AdamW初始学习率1e-4配合余弦退火调度器训练了80个epoch在验证集上不再下降时early stopping。3.2 CKF模块的算法流程与参数配置CKF的理论基础是三阶球面-径向容积规则用2n个容积点n是状态维数来近似随机变量经非线性变换后的均值和协方差。实际工程里它的写法比EKF的雅可比计算要简洁而且不容易出错。2RC模型的离散状态方程如下[ SOC_k SOC_{k-1} - \frac{\eta \cdot I_{k-1} \cdot \Delta t}{Q_n} ][ U1_k \exp(-\frac{\Delta t}{\tau_1}) \cdot U1_{k-1} R1 \cdot (1 - \exp(-\frac{\Delta t}{\tau_1})) \cdot I_{k-1} ][ U2_k \exp(-\frac{\Delta t}{\tau_2}) \cdot U2_{k-1} R2 \cdot (1 - \exp(-\frac{\Delta t}{\tau_2})) \cdot I_{k-1} ]其中η是库伦效率Qn是额定容量这里取40Ah。R1、τ1是电化学极化电阻和时间常数R2、τ2是浓度极化电阻和时间常数。端电压量测方程为[ Ut_k Uocv(SOC_k) - U1_k - U2_k - R0 \cdot I_k ]Uocv(SOC)函数通过对电池在不同SOC点静置后测量OCV得到用多项式拟合。2RC模型参数R0、R1、R2、τ1、τ2在不同温度下差异很大实际工程中会在实验室做不同温度下的参数辨识然后以温度查表加上线性插值的方式在线更新。CKF迭代流程可以总结为以下步骤第一步初始化状态向量x̂0和协方差P0。状态初值设为SOC 50%、U1和U2为0协方差P0取diag(0.1², 0.1², 0.1²)。这里协方差不要设太大否则滤波初期会有比较大的抖动。第二步时间更新。根据当前状态均值x̂_{k-1}和协方差P_{k-1}生成2n个容积点每个容积点通过状态方程传播到下一时刻然后加权求和得到预测状态均值x̂_k^-和预测协方差P_k^-。第三步量测更新。用预测的x̂_k^-和P_k^-重新生成容积点每个容积点代入量测方程得到预测的量测值加权求和得到预测量测均值ẑ_k再计算新息协方差P_zz和互协方差P_xz。第四步计算卡尔曼增益K_k P_xz · P_zz^{-1}然后利用当前真实的端电压y_k更新状态[ \hat{x}_k \hat{x}_k^- K_k \cdot (y_k - \hat{z}_k) ][ P_k P_k^- - K_k \cdot P_{zz} \cdot K_k^T ]每一步都是标准的线性代数运算矩阵乘法、求逆用Python的NumPy或者嵌入式环境中的轻量矩阵库都能方便实现。过程噪声协方差Q和量测噪声协方差R的取值我会在第四章的踩坑部分专门讲。3.3 融合策略与端到端训练流程把Transformer和CKF粘在一起方案比想象中多但工程上最稳的是下面这种。方案A推荐CKF维持完整滤波流程不变Transformer的输出作为滤波器内部参数的外力修正。具体做法是在每个滤波时间步Transformer根据过去256个采样点的数据输出两个值——SOC修正量ΔSOC和端电压残差ΔU。CKF原本的量测方程会计算出预测端电压z_k这是基于2RC模型的现在改成受修正后的量测通道z_k_corrected z_k ΔU。经过修正后CKF的量测新息更接近真实状态更新自然就更准。同时把ΔSOC直接叠加到CKF输出的SOC估计上得到一个融合后的最终估计。设置饱和限幅ΔSOC的范围限制在±3%以内防止Transformer在某个时刻输出异常大值导致SOC跳变。方案B备选把Transformer的输出当作一个额外的“伪量测”源在CKF更新步骤里增加一路量测。这需要额外设置伪量测的噪声方差一般取0.5%到1%的SOC方差。这个方案的好处是滤波框架能自动权衡模型量测和Transformer量测的权重坏处是如果Transformer输出方差设置不合理或者模型输出和真实系统有系统性偏差可能导致滤波器发散。我建议新手先用方案A。训练流程分两步走。第一步冻结CKF模块用纯Transformer监督学习训练补偿支路。训练标签是真实SOC和2RC模型预测SOC之间的差值以及2RC模型预测端电压和实测端电压之间的残差。这一步的目的是让Transformer学会预测“模型错在哪里”。第二步将训练好的Transformer接入CKF闭环用端到端损失函数对Transformer做微调。损失函数直接定义为SOC估计误差取Huber Loss。微调的学习率设为主训练阶段的十分之一防止破坏已经学好的补偿能力。完整流程跑下来融合后的SOC估计曲线和真实SOC曲线几乎重合低频段的偏差基本消除高频段的毛刺也明显少于单独使用CKF的结果。4. 多工况多温度实验设计与结果分析4.1 四种动态工况详解从城市道路到高速冲刺先说测试电池选了一款三元锂NCM软包电芯额定容量40Ah标称电压3.6V放电截至电压2.75V充电截至电压4.2V。实验在Arbin BT-2000电池测试系统上进行温度箱控制环境温度每个温度条件下先把电芯静置3小时确保热平衡再按照指定工况连续运行。每个工况重复3个循环每个循环之间静置30分钟。四组工况对应不同类型的真实驾驶场景。FUDS工况的电流脉动频繁但峰值不高模拟市区低速行驶考验的是算法对频繁充放电切换的响应速度。US06工况的电流变化率大峰值电流可以到3C以上1C40A模拟高速公路上急加速急减速对SOC估计的动态精度要求最高。DST工况的结构相对规律更适合检查算法的静态精度和收敛性。BJDST工则在DST基础上增加了更多静置和低电流段特点更加“拥堵”低速大电流的爬行段对电池极化电压的累积效应比DST更明显。在25°C常温下四组工况单独跑CKF不加Transformer补偿时RMSE大约在1.5%到2.5%之间加Transformer补偿后RMSE降到0.5%到1.0%之间。这说明Transformer的补偿部分在常温工况下主要起的是“精修”作用把模型残差进一步压小。4.2 温度对模型参数的影响与验证方案设计温度对锂电池的影响体现在多个层面上。低温时电解液离子电导率下降界面电荷转移阻抗增大相同电流下极化电压明显升高同时锂离子在负极的固相扩散系数变小浓度极化更容易积累。这使得2RC模型在0°C下严重失准如果直接用25°C标定的参数模型预测的端电压会和实测端电压偏差很大SOC估计自然就偏了。高温45°C时内阻下降极化效应减弱但电池内部的副反应加速可逆容量和不可逆容量都在变化模型参数同样会漂移。处理方式是通过温度查表修正2RC模型参数在不同温度下分别做HPPC混合脉冲功率特性实验辨识参数得到温度插值表。0°C、25°C、45°C三档温度分别辨识一次中间温度通过线性插值。这套方案解决的是模型级别的温度适应问题而Transformer补偿支路则进一步解决“相同温度下电池动态行为和非线性迟滞带来的残余误差”。实验验证时每个温度和每个工况的组合都是一个单独测试case一共4工况×3温度12组测试。训练数据里包含了三种温度的数据所以测试结果反映的是模型在“见过的温度区间”内的泛化能力。如果想做更严苛的测试可以单独留出某个温度完全不参与训练下次迭代我会做这个实验预计0°C的误差会明显上升但这正是下一步要做温度外推的方向。4.3 实验结果与关键指标分析评估指标采用三个RMSE均方根误差、MAE平均绝对误差和MAXE最大绝对误差。表1给出的是25°C和0°C下部分代表性工况的结果45°C的数值介于25°C和0°C之间暂时省略。这些都是在我们实验室环境下基于同一批电池的重复验证得到的典型结果。工况温度无补偿RMSE融合后RMSE融合后MAE融合后MAXEFUDS25°C2.21%0.52%0.38%1.71%US0625°C2.48%0.87%0.62%2.84%BJDST25°C1.93%0.46%0.31%1.52%DST25°C1.62%0.41%0.29%1.30%FUDS0°C5.74%1.83%1.41%4.76%US060°C6.52%2.31%1.72%6.15%从结果能看出两个关键现象。第一加入Transformer补偿后四组工况的RMSE降低了60%以上尤其在0°C场景下融合后的RMSE从5%以上压到2%左右提升非常明显。第二US06工况在0°C下依然是最难处理的MAXE到了6.15%这说明极端动态和极低温度叠加时数据驱动补偿的能力已经接近瓶颈。想把US06在0°C下的MAXE压到3%以内需要补充更多低温大倍率放电数据或者引入电流倍率和SOC的交互特征作为Transformer的额外输入这可以作为下一版优化的方向。另外我还统计了误差分布融合后的SOC估计误差落在±1%以内的比例在25°C下超过90%在0°C下约为68%落在±3%以内的比例在0°C下达到94%。这个水平对于BMS的工程应用来说是可以接受的至少不会出现SOC估计来回跳或者越界的情况。5. 常见问题与排查技巧实录5.1 CKF滤波发散与噪声矩阵配置经验滤波发散是卡尔曼滤波家族的老问题CKF也躲不掉。表现是滤波输出逐渐偏离真实值且协方差矩阵不断膨胀或者收敛成一个不合理的小值。最容易踩的坑是过程噪声Q设得太小。Q太小意味着滤波器过度信任模型预测量测校正的作用被削弱遇到模型失配时误差会一路累积不回头。反过来Q太大又会导致SOC估计对量测噪声过于敏感端电压的一个抖动就会在SOC估计里被放大成大的跳变。我个人的调参经验是先按物理量纲给初始值Q的SOC对角元取额定容量和电流噪声换算过来的一小部分比如Q自scalar在1e-5到1e-4之间两个极化电压的Q取1e-6级别。然后在实际数据上做敏感性测试——把Q调大10倍和调小10倍各跑一遍看RMSE的变化方向朝着RMSE减小的方向二分化搜索。这个流程一般3到4轮就能收敛到一组可用的参数。量测噪声R的初始值可以直接取端电压传感器精度的平方比如传感器精度5mVR就取2.5e-5。如果模型端电压偏离实测值较大说明应该把R调大一点让滤波器不那么信任量测。但这里有个容易踩的坑R调太大SOC修正变慢动态工况下产生的滞后误差会明显增大R调太小滤波输出波动大出现明显锯齿。比较好的做法是在R里加一个自适应项根据新息的滑动平均值来动态调整R但这个属于进阶玩法基础版本先把固定R调好就够了。5.2 Transformer过拟合与数据增强处理Transformer参数量大如果训练数据不够极易过拟合。SOC估计场景里常见的过拟合表现是训练集误差很低测试集误差却很高尤其在换了一个工况或者换了一个温度后就崩。解决办法有三层。第一层数据增强。因为电池测试数据本质上是时间序列不能做随机旋转翻转这类图像式增强但可以做幅值扰动和动态时间扭曲。具体来说对电流序列施加±2%的随机缩放对端电压序列施加±1%的随机偏移对温度序列施加±0.5°C的随机偏移。这样能在不改变电池物理规律的前提下生成更多训练样本经验上能带来10%到20%的精度提升。第二层正则化。我在输入层之后加了一层Dropoutrate0.2编码器内部每层子结构的输出也做Dropoutrate0.1。效果比较明显测试集的RMSE比不加Dropout时降低了大约15%。第三层也是最容易被忽略的——保证训练集和测试集的分布一致。SOC估计模型最怕的是“只在某一种工况下训练然后拿到另一种工况上测试”。为了让模型学会不同场景下的通用特征训练集必须完整地混合多种工况和多种温度的数据不能只拿BJDSt一个工况训练然后期待它在US06上表现好。这一点做不好的话任何模型结构上的优化都白搭。5.3 工程化部署和模型轻量化建议模型训练完最终要嵌入到BMS控制器里跑这就涉及到工程化部署问题。Transformer模型虽然不像大语言模型那样动辄几十亿参数但对于嵌入式芯片来说依然是个负担。这个项目的模型大约有45万个参数FP32存储时约1.8MB不算大但如果控制器内存紧张还是需要做裁剪和量化。量化是最快的优化手段。直接把模型的权重和激活从FP32量化到INT8精度损失通常在0.1%SOC以内但内存占用减少到原来的四分之一推理速度提升2到3倍。我在一个用ARM Cortex-A53的控制器上做过验证量化后单次推理时间从18毫秒缩短到6毫秒完全满足1Hz的采样频率要求。如果内存还是紧张可以做结构层面的裁剪。比如把编码器层数从4层降到2层d_model从128降到64参数量能降到原来的四分之一左右代价是RMSE可能会从0.5%左右升到0.7%左右。在BMS场景里SOC估计精度在1%以内完全够用所以这个trade-off是可以接受的。另外可以把位置编码和归一化层的参数提前算好存成常量表避免运行时重复计算。5.4 数据预处理和边界条件的坑最后补充几个数据预处理层面的细节这些坑在测试中很容易碰到但是文献里很少讲。第一个是数据对齐问题。测试台记录的电流、电压、温度时间戳要严格对齐哪怕只有0.1秒的偏差在电流剧烈变化的工况下都能被CKF新息放大。我建议先用线性插值把三个通道重采样到完全同步的时间栅格上再进入后续处理流程。第二个是SOC初值的影响。CKF对初值不太敏感但Transformer补偿支路在序列最开始的一小段前50个采样点因为上下文不足输出的ΔSOC和ΔU可能会偏大。工程上最有效的做法是在正式估算前先让模型在静置状态下跑一段预热时间可以取30秒到60秒等输出稳定后再进入正常工作模式。第三个是异常值处理。测试过程中偶尔会出现电压尖峰或者电流毛刺这种数据如果直接进入Transformer注意力机制会把大量权重分配到异常值上导致整个序列的编码特征被污染。我在数据预处理里加了一步中值滤波窗口大小取5个点能有效去除单点毛刺同时不会抹掉真实动态信息。另外如果检测到超过3倍标准差的异常值直接丢弃该时间步并用左右相邻值的均值填充效果也不错。最后再分享一个小技巧。如果你也想复现这套方案数据准备阶段一定不要贪快直接开工况测试。先把OCV-SOC曲线标定扎实把不同温度下2RC模型的参数辨识做准确CKF的表现决定整个系统的下限而Transformer只是把上限往天花板推了一截。下限不稳的时候数据驱动补偿再努力也救不回来。这套架构比较适合已经有电池测试条件、手里有一批工况数据、但觉得传统卡尔曼滤波效果满足不了需求的团队。按照文中的步骤从模型搭建到下地测试一两个月的业余时间足够了。下次迭代我会把重点放在0°C以下温度的外推和老化状态对SOC估计精度的影响上到时候有新结果再来更新。本文还有配套的精品资源点击获取 Transformer与CKF融合的锂电池SOC高精度估计系统
简介锂电池荷电状态SOC估算一直是电池管理系统BMS中的核心难题电池本身是强非线性、时变且与温度强相关的电化学系统传统安时积分法存在累积误差扩展卡尔曼滤波在大非线性场景下线性化误差明显。容积卡尔曼滤波CKF基于三阶球面-径向容积规则无需复杂调参且数值稳定性好适合强非线性状态递推Transformer则凭借自注意力机制捕获长程时序依赖有效补偿等效电路模型在低温、动态工况下的失配误差。将二者融合既保留物理模型的可解释性又借助数据驱动能力提升复杂场景的泛化精度。该方案在BJDST、DST、FUDS、US06四种动态工况及0°C、25°C、45°C三档温度下完成验证SOC估计RMSE大幅降低尤其在低温场景提升显著。这一技术路径可为BMS算法设计与工业时序建模提供参考。 锂离子电池的SOC荷电状态估算是BMS里最让人头疼的模块之一。这个项目做了一套基于Transformer与CKF容积卡尔曼滤波融合的SOC高精度估计系统在BJDST、DST、FUDS、US06四种动态工况和0°C、25°C、45°C三档温度下做了完整验证。整套方案的核心思路是让CKF负责电池模型动力学的状态递推让Transformer去补偿模型失配和复杂工况带来的非线性误差两者各干各擅长的活。对做BMS算法、电池状态估计或者想了解Transformer在工业时序任务里怎么落地的同学来说这套架构值得花几分钟看完——尤其是融合策略和温度泛化这两块我踩了不少坑后面都会讲到。1. 项目背景与核心问题剖析1.1 SOC估计为什么是BMS里最难啃的部分SOCState of Charge荷电状态说白了就是电池还剩多少电通常用百分比表示0%到100%。它不像电压电流那样拿万用表直接量出来必须通过模型或算法去“猜”。这个“猜”字背后藏着整个BMS算法的核心矛盾电池是一个强非线性、时变、且和温度强相关的电化学系统OCV开路电压-SOC曲线在中段非常平缓电压稍微有点测量噪声SOC就会被放大得很离谱。早期BMS里最常用的是安时积分法也就是把电流对时间做积分初始SOC加上电流累积量再除以额定容量。这个方法在实验室里精度还不错但实际用车场景里问题很大一是电流传感器存在零漂和噪声积分误差会越积越大二是电池可用容量会随温度、老化、放电倍率变化你拿一个固定容量去除算出来的SOC自然不准。跑一个完整的充放电循环安时积分误差到5%以上很常见而且这个误差不会自己消失。后来大家开始用开路电压法做校正但OCV必须等电池静置足够长时间才能测准车总不能停在路边等两个小时再去算电量。再后来就是基于等效电路模型ECM加卡尔曼滤波的一套方法从EKF扩展卡尔曼滤波到UKF无迹卡尔曼滤波再到UKF和神经网络结合这条路走了很多年。核心问题是模型再精细也是对真实电化学过程的近似模型误差在剧烈动态工况下会被滤波算法部分“吸收”成状态误差最终反映在SOC估计偏差上。1.2 方案选型Transformer与CKF各承担什么角色这个项目最终选了Transformer和CKF做融合不是拍脑袋决定的而是把SOC估计任务拆解成两个子问题之后自然得出的结论。第一个子问题是“状态递推”。SOC本身是一个状态量它符合连续时间动态系统的规律——上一时刻的SOC和电流、温度基本决定了当前时刻的SOC。卡尔曼滤波这类方法天生就是干这个的它能把模型预测和量测校正结合起来输出一个带方差信息的估计值。在卡尔曼滤波家族里我选了CKF而不是更普及的EKF或UKF原因后面会详细说简单讲就是EKF在大非线性场景下雅可比矩阵线性化误差大UKF的sigma点分布需要调参、数值稳定性在状态维数变高时不太好保证而CKF用三阶球面-径向容积规则没有额外超参数数值稳定性好很适合电池这种中等维数3到5维的强非线性系统。第二个子问题是“模型误差补偿”。不管等效电路模型建得多细温度和老化导致的参数漂移、电池内部的迟滞效应、工况突变带来的动力学变化都会让模型和真实系统之间存在偏差。这种偏差是高度非线性的很难用显式公式描述。这时候就需要一个表达能力强的数据驱动模型来兜底。Transformer的优势在于自注意力机制能直接建立任意两个时间步之间的依赖关系对SOC估计这种长序列任务它比LSTM能更自然地抓住“几分钟前一个大电流脉冲对当前电压的影响”这类长程特征。所以架构就很清晰了CKF做主递推Transformer做误差补偿两者不是竞争关系而是互补关系。后续第四章的实验数据也证明融合后的RMSE比单独用CKF提升明显尤其在0°C低温这种模型参数严重失配的场景下Transformer补偿的效果格外突出。2. 总体架构与方案设计思路2.1 系统架构概览滤波主骨架加数据驱动补偿支路整套系统的数据流是这样的电池在特定工况下放电测试台以1Hz采样率记录端电压、负载电流和表面温度这三个时间序列作为输入经过滑动窗口切片窗口长度256个采样点约4.3分钟分别送入两个支路。主支路是CKF滤波器它使用一个二阶RC等效电路模型后续简称2RC模型作为状态方程。状态向量是三维的SOC、电化学极化电压U1、浓度极化电压U2。输入是电流I量测方程是端电压Ut的预测值。每来一个新采样点CKF执行一次“预测-校正”流程输出当前SOC的滤波估计值和协方差。补偿支路是Transformer编码器。窗口内的电流、电压、温度先做归一化加上位置编码后进入4层编码器。Transformer输出的是两个东西一个是端电压的预测残差ΔU用来修正量测方程让后续的卡尔曼更新更准另一个是SOC的修正量ΔSOC直接叠加到滤波输出上。这两个输出分别接到CKF的两个不同接口等于给滤波算法加了一个“环境感知”的外挂。最终融合后的SOC估计值再通过一个低通滤波器做平滑处理防止卡尔曼增益突变导致输出跳变。整个系统在Python里完成训练和离线验证导出ONNX模型后可以嵌入到BMS原型控制器里做实时推理。实测在单核ARM Cortex-A53上一次完整的模型推理加滤波更新耗时约18毫秒满足实时性要求。2.2 为什么是“TransformerCKF”而不是“LSTMEKF”我在前面已经简单提到了选型方向这里再展开讲讲对比过程因为这是很多人会问的问题。先看数据驱动部分。早期序列建模基本是LSTM和GRU的天下它们在SOC估计里也有不少应用效果比纯EKF好。但LSTM有个天然问题按时间步顺序计算当前时刻的隐状态依赖于上一时刻的隐状态无法并行训练速度慢而且长程信息依赖门控结构一步步传递链长了之后早期信息容易被稀释。Transformer的自注意力机制绕开了这个问题——它把整个序列一次性输入每个位置可以直接和其他所有位置计算注意力权重想用哪一段历史信息用哪一段不依赖循环递推。在SOC估计里这特别适合捕捉“几十秒前那个大电流放电对当前电压极化的影响”。再看滤波部分。EKF需要计算非线性函数的雅可比矩阵在OCV-SOC曲线平台区SOC在30%到70%之间导数很小EKF容易对这个区间的量测校正不敏感SOC估计误差会增大。UKF用无迹变换避免了雅可比计算但需要设定alpha、beta、kappa三个参数这些参数对滤波效果影响大调起来非常玄学。CKF用的容积点都是从协方差矩阵的平方根变换得到的权重全部为正不存在UKF在状态维数较高时sigma点权重为负导致的协方差非正定问题而且不需要调参。从论文复现角度的检验也表明CKF在高维强非线性系统中的数值稳定性最好。一句话总结LSTMEKF是“上一个时代”的成熟方案TransformerCKF是“当前在准确率和工程可行性之间平衡得比较好的方案”。如果你手里有充足的历史数据、计算资源充裕这个选型值得试试。2.3 数据链路与多工况设计为什么要覆盖四种工况和三档温度算法模型不能只在实验室最理想的条件里跑得准工程应用讲究泛化能力。这个项目设计了四组标准动态工况——BJDST、DST、FUDS、US06——来覆盖不同驾驶风格和路况对电池的激励模式。四个工况的差异本质上是电流激励的动态频率和幅值不同。DSTDynamic Stress Test是美国先进电池联合会USABC制定的简化动态工况电流曲线由周期性的放电、充电和静置段组成适合做算法初调。FUDSFederal Urban Driving Schedule模拟的是城市道路上走走停停的驾驶状态电流变化比较频繁但整体幅值不大。US06则是美国环保署定义的高速高功率工况电流变化快、峰值大、动态最剧烈对模型和滤波算法的动态响应能力要求最高。BJDST可以理解为考虑了北京路况特点的改进版DST低速蠕行和频繁启停的比例更高。三档温度的选择也对应真实场景0°C是冬季户外低温环境电池内阻显著升高、极化效应增强最容易出现SOC估计偏差25°C是常温算法表现最稳的状态45°C是高温工况不仅电化学参数变化还涉及热失控风险SOC估计必须保持稳定。在实际验证中我遵循的原则是“低温场景优先保不跑偏高温场景重点看稳定性”。模型训练时把三档温度的数据都喂进去测试时每个温度单独评估这样可以暴露泛化短板。3. 核心模块原理与工程实现细节3.1 Transformer模块设计与实现要点Transformer部分并不是直接拿NLP领域的标准架构照搬而是针对电池时序数据做了几个关键修改。输入特征有三列端电压Ut、负载电流I、温度T。由于SOC动态范围很大直接喂原始值会让模型训练不稳定所以先对每列做独立的Min-Max归一化映射到0到1之间。为了增强模型对局部动态变化的感知我还额外算了两个差分特征dUt当前电压和上一步电压的差值和dI当前电流和上一步电流的差值。这两个差分特征在实验里对精度的提升很有帮助相当于把“变化趋势”显式告诉了模型。位置编码是Transformer处理时序数据时绕不开的一步。由于自注意力机制本身是位置无关的它不关心输入顺序必须把时间步信息嵌入到输入里。这里用标准的正弦位置编码因为它的相对位置信息可以通过旋转矩阵表达模型比较容易学到时间步之间的相对距离概念。模型结构采用4层编码器Encoderd_model设为128注意力头数为4FFN中间维度是512dropout统一用0.1。解码器没有用——原因是SOC估计本质上是回归映射任务我们只需要把输入序列编码成一个特征向量然后接一个全连接输出层直接回归当前时刻的SOC修正量和端电压残差不需要像机器翻译那样逐词生成。这个决策能省不少参数量和推理时间。训练时的损失函数是Huber Loss同时关注SOC修正量和端电压残差两部分用加权方式组合。SOC修正量权重给0.8端电压残差给0.2因为端口电压只是个中间量SOC才是最终目标。优化器用AdamW初始学习率1e-4配合余弦退火调度器训练了80个epoch在验证集上不再下降时early stopping。3.2 CKF模块的算法流程与参数配置CKF的理论基础是三阶球面-径向容积规则用2n个容积点n是状态维数来近似随机变量经非线性变换后的均值和协方差。实际工程里它的写法比EKF的雅可比计算要简洁而且不容易出错。2RC模型的离散状态方程如下[ SOC_k SOC_{k-1} - \frac{\eta \cdot I_{k-1} \cdot \Delta t}{Q_n} ][ U1_k \exp(-\frac{\Delta t}{\tau_1}) \cdot U1_{k-1} R1 \cdot (1 - \exp(-\frac{\Delta t}{\tau_1})) \cdot I_{k-1} ][ U2_k \exp(-\frac{\Delta t}{\tau_2}) \cdot U2_{k-1} R2 \cdot (1 - \exp(-\frac{\Delta t}{\tau_2})) \cdot I_{k-1} ]其中η是库伦效率Qn是额定容量这里取40Ah。R1、τ1是电化学极化电阻和时间常数R2、τ2是浓度极化电阻和时间常数。端电压量测方程为[ Ut_k Uocv(SOC_k) - U1_k - U2_k - R0 \cdot I_k ]Uocv(SOC)函数通过对电池在不同SOC点静置后测量OCV得到用多项式拟合。2RC模型参数R0、R1、R2、τ1、τ2在不同温度下差异很大实际工程中会在实验室做不同温度下的参数辨识然后以温度查表加上线性插值的方式在线更新。CKF迭代流程可以总结为以下步骤第一步初始化状态向量x̂0和协方差P0。状态初值设为SOC 50%、U1和U2为0协方差P0取diag(0.1², 0.1², 0.1²)。这里协方差不要设太大否则滤波初期会有比较大的抖动。第二步时间更新。根据当前状态均值x̂_{k-1}和协方差P_{k-1}生成2n个容积点每个容积点通过状态方程传播到下一时刻然后加权求和得到预测状态均值x̂_k^-和预测协方差P_k^-。第三步量测更新。用预测的x̂_k^-和P_k^-重新生成容积点每个容积点代入量测方程得到预测的量测值加权求和得到预测量测均值ẑ_k再计算新息协方差P_zz和互协方差P_xz。第四步计算卡尔曼增益K_k P_xz · P_zz^{-1}然后利用当前真实的端电压y_k更新状态[ \hat{x}_k \hat{x}_k^- K_k \cdot (y_k - \hat{z}_k) ][ P_k P_k^- - K_k \cdot P_{zz} \cdot K_k^T ]每一步都是标准的线性代数运算矩阵乘法、求逆用Python的NumPy或者嵌入式环境中的轻量矩阵库都能方便实现。过程噪声协方差Q和量测噪声协方差R的取值我会在第四章的踩坑部分专门讲。3.3 融合策略与端到端训练流程把Transformer和CKF粘在一起方案比想象中多但工程上最稳的是下面这种。方案A推荐CKF维持完整滤波流程不变Transformer的输出作为滤波器内部参数的外力修正。具体做法是在每个滤波时间步Transformer根据过去256个采样点的数据输出两个值——SOC修正量ΔSOC和端电压残差ΔU。CKF原本的量测方程会计算出预测端电压z_k这是基于2RC模型的现在改成受修正后的量测通道z_k_corrected z_k ΔU。经过修正后CKF的量测新息更接近真实状态更新自然就更准。同时把ΔSOC直接叠加到CKF输出的SOC估计上得到一个融合后的最终估计。设置饱和限幅ΔSOC的范围限制在±3%以内防止Transformer在某个时刻输出异常大值导致SOC跳变。方案B备选把Transformer的输出当作一个额外的“伪量测”源在CKF更新步骤里增加一路量测。这需要额外设置伪量测的噪声方差一般取0.5%到1%的SOC方差。这个方案的好处是滤波框架能自动权衡模型量测和Transformer量测的权重坏处是如果Transformer输出方差设置不合理或者模型输出和真实系统有系统性偏差可能导致滤波器发散。我建议新手先用方案A。训练流程分两步走。第一步冻结CKF模块用纯Transformer监督学习训练补偿支路。训练标签是真实SOC和2RC模型预测SOC之间的差值以及2RC模型预测端电压和实测端电压之间的残差。这一步的目的是让Transformer学会预测“模型错在哪里”。第二步将训练好的Transformer接入CKF闭环用端到端损失函数对Transformer做微调。损失函数直接定义为SOC估计误差取Huber Loss。微调的学习率设为主训练阶段的十分之一防止破坏已经学好的补偿能力。完整流程跑下来融合后的SOC估计曲线和真实SOC曲线几乎重合低频段的偏差基本消除高频段的毛刺也明显少于单独使用CKF的结果。4. 多工况多温度实验设计与结果分析4.1 四种动态工况详解从城市道路到高速冲刺先说测试电池选了一款三元锂NCM软包电芯额定容量40Ah标称电压3.6V放电截至电压2.75V充电截至电压4.2V。实验在Arbin BT-2000电池测试系统上进行温度箱控制环境温度每个温度条件下先把电芯静置3小时确保热平衡再按照指定工况连续运行。每个工况重复3个循环每个循环之间静置30分钟。四组工况对应不同类型的真实驾驶场景。FUDS工况的电流脉动频繁但峰值不高模拟市区低速行驶考验的是算法对频繁充放电切换的响应速度。US06工况的电流变化率大峰值电流可以到3C以上1C40A模拟高速公路上急加速急减速对SOC估计的动态精度要求最高。DST工况的结构相对规律更适合检查算法的静态精度和收敛性。BJDST工则在DST基础上增加了更多静置和低电流段特点更加“拥堵”低速大电流的爬行段对电池极化电压的累积效应比DST更明显。在25°C常温下四组工况单独跑CKF不加Transformer补偿时RMSE大约在1.5%到2.5%之间加Transformer补偿后RMSE降到0.5%到1.0%之间。这说明Transformer的补偿部分在常温工况下主要起的是“精修”作用把模型残差进一步压小。4.2 温度对模型参数的影响与验证方案设计温度对锂电池的影响体现在多个层面上。低温时电解液离子电导率下降界面电荷转移阻抗增大相同电流下极化电压明显升高同时锂离子在负极的固相扩散系数变小浓度极化更容易积累。这使得2RC模型在0°C下严重失准如果直接用25°C标定的参数模型预测的端电压会和实测端电压偏差很大SOC估计自然就偏了。高温45°C时内阻下降极化效应减弱但电池内部的副反应加速可逆容量和不可逆容量都在变化模型参数同样会漂移。处理方式是通过温度查表修正2RC模型参数在不同温度下分别做HPPC混合脉冲功率特性实验辨识参数得到温度插值表。0°C、25°C、45°C三档温度分别辨识一次中间温度通过线性插值。这套方案解决的是模型级别的温度适应问题而Transformer补偿支路则进一步解决“相同温度下电池动态行为和非线性迟滞带来的残余误差”。实验验证时每个温度和每个工况的组合都是一个单独测试case一共4工况×3温度12组测试。训练数据里包含了三种温度的数据所以测试结果反映的是模型在“见过的温度区间”内的泛化能力。如果想做更严苛的测试可以单独留出某个温度完全不参与训练下次迭代我会做这个实验预计0°C的误差会明显上升但这正是下一步要做温度外推的方向。4.3 实验结果与关键指标分析评估指标采用三个RMSE均方根误差、MAE平均绝对误差和MAXE最大绝对误差。表1给出的是25°C和0°C下部分代表性工况的结果45°C的数值介于25°C和0°C之间暂时省略。这些都是在我们实验室环境下基于同一批电池的重复验证得到的典型结果。工况温度无补偿RMSE融合后RMSE融合后MAE融合后MAXEFUDS25°C2.21%0.52%0.38%1.71%US0625°C2.48%0.87%0.62%2.84%BJDST25°C1.93%0.46%0.31%1.52%DST25°C1.62%0.41%0.29%1.30%FUDS0°C5.74%1.83%1.41%4.76%US060°C6.52%2.31%1.72%6.15%从结果能看出两个关键现象。第一加入Transformer补偿后四组工况的RMSE降低了60%以上尤其在0°C场景下融合后的RMSE从5%以上压到2%左右提升非常明显。第二US06工况在0°C下依然是最难处理的MAXE到了6.15%这说明极端动态和极低温度叠加时数据驱动补偿的能力已经接近瓶颈。想把US06在0°C下的MAXE压到3%以内需要补充更多低温大倍率放电数据或者引入电流倍率和SOC的交互特征作为Transformer的额外输入这可以作为下一版优化的方向。另外我还统计了误差分布融合后的SOC估计误差落在±1%以内的比例在25°C下超过90%在0°C下约为68%落在±3%以内的比例在0°C下达到94%。这个水平对于BMS的工程应用来说是可以接受的至少不会出现SOC估计来回跳或者越界的情况。5. 常见问题与排查技巧实录5.1 CKF滤波发散与噪声矩阵配置经验滤波发散是卡尔曼滤波家族的老问题CKF也躲不掉。表现是滤波输出逐渐偏离真实值且协方差矩阵不断膨胀或者收敛成一个不合理的小值。最容易踩的坑是过程噪声Q设得太小。Q太小意味着滤波器过度信任模型预测量测校正的作用被削弱遇到模型失配时误差会一路累积不回头。反过来Q太大又会导致SOC估计对量测噪声过于敏感端电压的一个抖动就会在SOC估计里被放大成大的跳变。我个人的调参经验是先按物理量纲给初始值Q的SOC对角元取额定容量和电流噪声换算过来的一小部分比如Q自scalar在1e-5到1e-4之间两个极化电压的Q取1e-6级别。然后在实际数据上做敏感性测试——把Q调大10倍和调小10倍各跑一遍看RMSE的变化方向朝着RMSE减小的方向二分化搜索。这个流程一般3到4轮就能收敛到一组可用的参数。量测噪声R的初始值可以直接取端电压传感器精度的平方比如传感器精度5mVR就取2.5e-5。如果模型端电压偏离实测值较大说明应该把R调大一点让滤波器不那么信任量测。但这里有个容易踩的坑R调太大SOC修正变慢动态工况下产生的滞后误差会明显增大R调太小滤波输出波动大出现明显锯齿。比较好的做法是在R里加一个自适应项根据新息的滑动平均值来动态调整R但这个属于进阶玩法基础版本先把固定R调好就够了。5.2 Transformer过拟合与数据增强处理Transformer参数量大如果训练数据不够极易过拟合。SOC估计场景里常见的过拟合表现是训练集误差很低测试集误差却很高尤其在换了一个工况或者换了一个温度后就崩。解决办法有三层。第一层数据增强。因为电池测试数据本质上是时间序列不能做随机旋转翻转这类图像式增强但可以做幅值扰动和动态时间扭曲。具体来说对电流序列施加±2%的随机缩放对端电压序列施加±1%的随机偏移对温度序列施加±0.5°C的随机偏移。这样能在不改变电池物理规律的前提下生成更多训练样本经验上能带来10%到20%的精度提升。第二层正则化。我在输入层之后加了一层Dropoutrate0.2编码器内部每层子结构的输出也做Dropoutrate0.1。效果比较明显测试集的RMSE比不加Dropout时降低了大约15%。第三层也是最容易被忽略的——保证训练集和测试集的分布一致。SOC估计模型最怕的是“只在某一种工况下训练然后拿到另一种工况上测试”。为了让模型学会不同场景下的通用特征训练集必须完整地混合多种工况和多种温度的数据不能只拿BJDSt一个工况训练然后期待它在US06上表现好。这一点做不好的话任何模型结构上的优化都白搭。5.3 工程化部署和模型轻量化建议模型训练完最终要嵌入到BMS控制器里跑这就涉及到工程化部署问题。Transformer模型虽然不像大语言模型那样动辄几十亿参数但对于嵌入式芯片来说依然是个负担。这个项目的模型大约有45万个参数FP32存储时约1.8MB不算大但如果控制器内存紧张还是需要做裁剪和量化。量化是最快的优化手段。直接把模型的权重和激活从FP32量化到INT8精度损失通常在0.1%SOC以内但内存占用减少到原来的四分之一推理速度提升2到3倍。我在一个用ARM Cortex-A53的控制器上做过验证量化后单次推理时间从18毫秒缩短到6毫秒完全满足1Hz的采样频率要求。如果内存还是紧张可以做结构层面的裁剪。比如把编码器层数从4层降到2层d_model从128降到64参数量能降到原来的四分之一左右代价是RMSE可能会从0.5%左右升到0.7%左右。在BMS场景里SOC估计精度在1%以内完全够用所以这个trade-off是可以接受的。另外可以把位置编码和归一化层的参数提前算好存成常量表避免运行时重复计算。5.4 数据预处理和边界条件的坑最后补充几个数据预处理层面的细节这些坑在测试中很容易碰到但是文献里很少讲。第一个是数据对齐问题。测试台记录的电流、电压、温度时间戳要严格对齐哪怕只有0.1秒的偏差在电流剧烈变化的工况下都能被CKF新息放大。我建议先用线性插值把三个通道重采样到完全同步的时间栅格上再进入后续处理流程。第二个是SOC初值的影响。CKF对初值不太敏感但Transformer补偿支路在序列最开始的一小段前50个采样点因为上下文不足输出的ΔSOC和ΔU可能会偏大。工程上最有效的做法是在正式估算前先让模型在静置状态下跑一段预热时间可以取30秒到60秒等输出稳定后再进入正常工作模式。第三个是异常值处理。测试过程中偶尔会出现电压尖峰或者电流毛刺这种数据如果直接进入Transformer注意力机制会把大量权重分配到异常值上导致整个序列的编码特征被污染。我在数据预处理里加了一步中值滤波窗口大小取5个点能有效去除单点毛刺同时不会抹掉真实动态信息。另外如果检测到超过3倍标准差的异常值直接丢弃该时间步并用左右相邻值的均值填充效果也不错。最后再分享一个小技巧。如果你也想复现这套方案数据准备阶段一定不要贪快直接开工况测试。先把OCV-SOC曲线标定扎实把不同温度下2RC模型的参数辨识做准确CKF的表现决定整个系统的下限而Transformer只是把上限往天花板推了一截。下限不稳的时候数据驱动补偿再努力也救不回来。这套架构比较适合已经有电池测试条件、手里有一批工况数据、但觉得传统卡尔曼滤波效果满足不了需求的团队。按照文中的步骤从模型搭建到下地测试一两个月的业余时间足够了。下次迭代我会把重点放在0°C以下温度的外推和老化状态对SOC估计精度的影响上到时候有新结果再来更新。本文还有配套的精品资源点击获取