(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111178293.9 (22)申请日 2021.10.09 (71)申请人 武汉科技大 学 地址 430081 湖北省武汉市青山区和平大 道947号 (72)发明人 张永　涂雷　薛志伟　郑英　赵敏　 (74)专利代理 机构 上海精晟知识产权代理有限 公司 3125 3 代理人 肖爱华 (51)Int.Cl. G06F 30/27(2020.01) G06N 3/04(2006.01) G06N 3/08(2006.01) G06F 119/02(2020.01) G06F 119/04(2020.01) (54)发明名称 一种基于W OUKF算法的锂离子电池剩余寿命 预测方法 (57)摘要 本发明提供一种基于WOUKF算法的锂离子电 池剩余寿命预测方法， 包括如下步骤： 步骤1： 对 电池数据集进行电池容量数据提取， 以电池容量 作为锂电池寿命 预测指标， 根据锂电池的循环周 期选择预测起始点， 利用预测起始点之前的电池 容量数据拟合模型， 并建立状态空间模型； 步骤 2： 利用WOUKF算法追踪已知历史数据， 直至得到 第T时刻状态变量， 获取在预测起始点之前每个 时刻的电池容量预测值与电池容量实测值之间 的残差数据； 步骤3： 残差 数据和电池容量实测值 训练LSTM模型， 获得未来时刻的残差数据预测值 及电池容量预测量测值； 步骤4： 更新状态向量， 预测出电池容量， 重复步骤3直到预测的电池容 量达到失效阈值时停止， 计算出RUL。 权利要求书5页 说明书14页 附图6页 CN 113901716 A 2022.01.07 CN 113901716 A 1.一种基于WOUK F算法的锂离 子电池剩余寿命预测方法， 其特 征在于， 包括以下步骤： 步骤1： 获得电池数据集并对该电池数据集进行电池容量数据提取， 以电池容量作为锂 电池寿命 预测指标， 并根据锂电池的循环周期 选择预测起始 点T， 并利用预测起始 点之前的 电池容量数据在基于傅里叶形式的锂离子电池退化模型上进 行拟合并确定模型初始 参数， 建立状态空间模型； 步骤2： 利用WOUKF算法对已知历史数据进行状态追踪， 直至得到第T时刻的状态变量， 获取在预测起始点之前的每 个时刻的电池容 量预测值与电池容 量实测值之间的残差数据； 步骤3： 将步骤2得到的残差数据和电池容量实测值作为训练集训练LSTM模型， 获得未 来时刻的残差数据预测值以及电池容 量预测量测值； 步骤4： 基于WOUKF算法和残差数据预测值来更新状态向量， 结合傅里叶模型预测出未 来时刻的电池容量， 重复步骤3直到预测的电池容量达到失效阈值时停止预测， 并根据达到 失效阈值的时间计算出锂电池的剩余寿命值RUL。 2.如权利要求1所述的基于WOUK F算法的锂离 子电池剩余寿命预测方法， 其特 征在于： 其中， 步骤1中的所述基于傅里叶形式的锂离 子电池退化模型的定义如下： y＝a+b*cos(d*k)+c *sin(d*k)       (1) 式(1)中， y表示电池容量； k表示循环次数； a,b,c,d表示模型初始参数； *表示乘法运 算； sin表示正弦函数； cos表述 余弦函数； 依据式(1)建立相应的状态空间模型为： 式(2)中， xk表示k时刻的状态向量； xk+1表示k+1时刻的状态向量； ωk表示协方差大小为 Qk的过程噪声； ak,bk,ck,dk表示k时刻的状态方程参数； yk表示k时刻的电池容量； υk表示协 方差大小为Rk的观测噪声。 3.如权利要求1所述的基于WOUK F算法的锂离 子电池剩余寿命预测方法， 其特 征在于： 其中， 步骤2中所述 WOUKF算法包括如下步骤： 步骤2‑1， 在k时刻用采样策略选取了N＝2n+1个sigma点集 计算此N个sigma点 对应的权值 和协方差 权值 步骤2‑2， 对步骤2 ‑1的每个sigma点， 用采样策略继续选取N＝2n+ 1个sigma点 潜在的sigma点为N*N＝(2n+1)2个， 计算此N*N个sigma点对应的权值为 和协方差权值 步骤2‑3， 对均值权 重进行更新： 权　利　要　求　书 1/5 页 2 CN 113901716 A 2式(16)中， 表示k时刻第i个采样点的权值； 表示k+1时刻观测值的后 验概率密度函数； 表示k+1时刻状态向量的后延概率密度函数； 表 示归一化常数； e表示对数函数； 表示k+1时刻状态向量方程； σ 表示比例参数， 以调 整采样点权 重之间的距离； 步骤2‑4， 归一化权值： 步骤2‑5， 更新k+1时刻的状态估计值及协方差： 式(18)(19)中， 表示k+1时刻的状态估计； 表示k+1时刻状态向量的加权平均 值； 表示k时刻第i个采样点的权值； 表示k时刻第i个采样点的协方差； 表示k +1时刻协方差加权平均值； Q表示过程噪声协方差； 步骤2‑6， 状态更新以及协方差更新： 取采样点实测值， 加权计算平均值： 式(20)至(24)中， 表示采样点观测值； 表示观测方程； 表示k时刻第 i个采样点的权值； 表示k+1时刻 采样点观测加权平均值； Pzz表示观测协方差加权平均 值； 表示k时刻第i个采样点的协方差； R表示大小为R的协方差矩阵； Pxz表示交叉协方差 加权平均值； 表示k+1时刻第i个采样点的状态估计加权平均值； 表示k+1时刻权　利　要　求　书 2/5 页 3 CN 113901716 A 3