[发明专利]一种新能源汽车动力电池SOH在线估算的方法

权利要求书

1.一种电动汽车动力电池SOH在线估算的方法，其特征是：包括以下步骤：

S1电池SOH预测模型的建立及优化

采用具有输出-输入反馈机制的改进Elman网络算法：OIF Elman网络，OIF Elman分为3层：输入层、隐层和输出层，其中输入层神经元数学模型为

x_u(k)＝u(k-1)  (1)；

x_u(k)表示k时刻输入层节点状态，u(k)表示k时刻系统输入量；

隐层神经元的数学模型为：

x(k)＝f(W¹¹x_c(k)+W¹²x_u(k)+W¹⁴y_c(k))  (2)

x_c(k)＝αx_c(k-1)+x(k-1)  (3)

y_c(k)＝γy_c(k-1)+y(k-1)  (4)；

其中，x(k)表示k时刻隐层节点的状态；x_c(k)表示k时刻联系单元1节点的状态；y_c(k)表示k时刻联系单元2节点的状态；W¹¹表示隐含节点与联系单元的连接权矩阵；W¹²表示隐层节点与输入节点之间的连接权矩阵；W¹⁴表示隐层节点与联系单元2之间的连接权矩阵；

输出层神经元的数学模型为：

y(k)＝W¹³x(k)  (5)；

其中，y(k)表示k时刻输出层节点的状态；W¹³表示隐层节点与输出节点之间的连接权矩阵；

所述的OIF Elman网络的隐层节点的输入为：输入层节点状态、联系单元1状态、联系单元2状态，设第k步系统的实际输出为y_d(k)，定义误差函数为:

E(k)=12(yd(k)-y(k))T(yd(k)-y(k))---(6);]]>

将E对连接权W¹¹、W¹²、W¹³、W¹⁴分别求偏导，由梯度下降法可得OIF Elman网络的学习算法：

ΔW_ij¹³＝η₃δ_i⁰x_j(k)  i＝1,2,…,m；j＝1,2,…,n  (7)

ΔW_jq¹²＝η₂δ_j^hu_q(k-1)  j＝1,2,…,n；q＝1,2,…,r  (8)

ΔWjl11=ηlΣi=1m(δi0Wij13)∂xi(k)∂Wjl11j=1,2,...,n;l=1,2,...,n---(9)]]>

ΔWjs14=η4Σi=1m(δi0Wij13)∂xi(k)∂Wjs14j=1,2,...,n;s=1,2,...,m---(10)]]>

∂xi(k)∂Wjs14=fj′(·)ys(k-1)+γ∂xi(k-1)∂Wjs14j=1,2,...,n;s=1,2,...,m---(11)]]>

其中，η₁，η₂，η₃，η₄分别为W¹¹、W¹²、W¹³、W¹⁴的学习步长，0≤γ＜1为联系单元2的自连接反馈增益因子，δ_i⁰和δ_j^h为误差信号；

δ_i⁰＝(y_di(k)-y_i(k))  (12)

δjh=Σi=1m(δi0Wij13)fj′(·)---(13)]]>

∂xi(k)∂Wj111=fj′(·)xl(k-1)+α∂xi(k-1)∂Wjl11---(14);]]>

所述的OIF Elman网络学习算法采用梯度下降法，即自适应学习速率动量梯度下降反向传播算法，它既能提高网络的训练速率，又能有效抑制网络陷入局部极小点，学习的目的是用网络的实际输出值与输出样本值的差值来修改权值和阀值，使网络输出层的误差平方和最小；

S2电池SOH模型输入参数的选择

选择电池内阻、电流和温度作为模型的输入参数；

确定隐层神经元的个数n₂和输入层神经元个数n₁：

隐层神经元的个数n₂和输入层神经元个数n₁之间有以下近似关系：

n₂＝2n₁+1  (14)

由样本数据可知，网络输入层有3个神经元，输出层有1个神经元，故分别将中间层神经元的数目设置为7；

S3在Matlab7.1环境下运行

对网络输入层和输出层的各单元值进行标准化处理，并使他们处于[0,1]或[-1,1]中，采用的输入参数的标准化方法为：

x^=x-xminxmax-xmin---(15)]]>

式中，x_max为该组变量的最大值，x_min为该组变量的最小值；

S4验证电池SOH预测模型的仿真分析选用6组相同型号的电池，采用相同充电制度充满电，在不同电流、不同温度的条件下，对3组电池进行放电试验，获得实测数据100组，选择其中80组数据对该模型进行训练，网络训练过程在Matlab7.1环境下运行，仿真中用改进Elman和Elman两种神经网络建模进行预测，分别得到预测结果，对比分析两者与实测SOH的误差。

2.根据权利要求1所述的电动汽车动力电池SOH在线估算的方法，其特征是：所述的步骤S3包括以下子步骤：

S3-1参数初始化；

S3-2批量输入学习样本，并对输入和输出量进行归一化处理；

S3-3计算各层以及联系单元输出；

S3-4计算输出误差；

S3-5判断输出误差满足要求否：若满足要求则存储训练完的权值后结束；若不满足要求则依公式12-13计算各层误差信号、依公式7-10调整各层权值、返回步骤S3-3。