发明名称 一种新能源汽车动力电池SOH在线估算的方法
摘要 一种电动汽车动力电池SOH在线估算的方法:S1.电池SOH预测模型的建立及优化,采用具有输出-输入反馈机制的改进Elman网络算法:OIF Elman网络;S2.电池SOH模型输入参数的选择:电池内阻、电流和温度作为模型的输入参数;S3.在Matlab7.1环境下运行;S4.验证。本发明在样本点较少时,无论在训练速度上,还是在预测精度上,OIF Elman网络明显优于Elman网络。OIF Elman网络提高了网络的泛化能力,既降低了对训练样本个数的需求,又能提高预测精度,在动力电池剩余寿命(SOH)预测中可以成功应用。
申请公布号 CN104035037A 申请公布日期 2014.09.10
申请号 CN201410199230.5 申请日期 2014.05.12
申请人 广东电网公司电力科学研究院 发明人 罗敏;孙卫明;肖勇;赵伟;黄默涵;孟金岭
分类号 G01R31/36(2006.01)I 主分类号 G01R31/36(2006.01)I
代理机构 广州知友专利商标代理有限公司 44104 代理人 周克佑
主权项 一种电动汽车动力电池SOH在线估算的方法,其特征是:包括以下步骤:S1电池SOH预测模型的建立及优化采用具有输出‑输入反馈机制的改进Elman网络算法:OIF Elman网络,OIF Elman分为3层:输入层、隐层和输出层,其中输入层神经元数学模型为x<sub>u</sub>(k)=u(k‑1) (1);x<sub>u</sub>(k)表示k时刻输入层节点状态,u(k)表示k时刻系统输入量;隐层神经元的数学模型为:x(k)=f(W<sup>11</sup>x<sub>c</sub>(k)+W<sup>12</sup>x<sub>u</sub>(k)+W<sup>14</sup>y<sub>c</sub>(k)) (2)x<sub>c</sub>(k)=αx<sub>c</sub>(k‑1)+x(k‑1) (3)y<sub>c</sub>(k)=γy<sub>c</sub>(k‑1)+y(k‑1) (4);其中,x(k)表示k时刻隐层节点的状态;x<sub>c</sub>(k)表示k时刻联系单元1节点的状态;y<sub>c</sub>(k)表示k时刻联系单元2节点的状态;W<sup>11</sup>表示隐含节点与联系单元的连接权矩阵;W<sup>12</sup>表示隐层节点与输入节点之间的连接权矩阵;W<sup>14</sup>表示隐层节点与联系单元2之间的连接权矩阵;输出层神经元的数学模型为:y(k)=W<sup>13</sup>x(k) (5);其中,y(k)表示k时刻输出层节点的状态;W<sup>13</sup>表示隐层节点与输出节点之间的连接权矩阵;所述的OIF Elman网络的隐层节点的输入为:输入层节点状态、联系单元1状态、联系单元2状态,设第k步系统的实际输出为y<sub>d</sub>(k),定义误差函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><msub><mi>y</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000503821400000011.GIF" wi="1320" he="127" /></maths>将E对连接权W<sup>11</sup>、W<sup>12</sup>、W<sup>13</sup>、W<sup>14</sup>分别求偏导,由梯度下降法可得OIF Elman网络的学习算法:ΔW<sub>ij</sub><sup>13</sup>=η<sub>3</sub>δ<sub>i</sub><sup>0</sup>x<sub>j</sub>(k) i=1,2,…,m;j=1,2,…,n (7)ΔW<sub>jq</sub><sup>12</sup>=η<sub>2</sub>δ<sub>j</sub><sup>h</sup>u<sub>q</sub>(k‑1) j=1,2,…,n;q=1,2,…,r (8)<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;</mi><msup><msub><mi>W</mi><mi>jl</mi></msub><mn>11</mn></msup><mo>=</mo><msub><mi>&eta;</mi><mi>l</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msup><msub><mi>&delta;</mi><mi>i</mi></msub><mn>0</mn></msup><msup><msub><mi>W</mi><mi>ij</mi></msub><mn>13</mn></msup><mo>)</mo></mrow><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msup><msub><mi>W</mi><mi>jl</mi></msub><mn>11</mn></msup></mrow></mfrac><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo><mi>l</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000503821400000021.GIF" wi="1368" he="136" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Delta;</mi><msup><msub><mi>W</mi><mi>js</mi></msub><mn>14</mn></msup><mo>=</mo><msub><mi>&eta;</mi><mn>4</mn></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msup><msub><mi>&delta;</mi><mi>i</mi></msub><mn>0</mn></msup><msup><msub><mi>W</mi><mi>ij</mi></msub><mn>13</mn></msup><mo>)</mo></mrow><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msup><msub><mi>W</mi><mi>js</mi></msub><mn>14</mn></msup></mrow></mfrac><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo><mi>s</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000503821400000022.GIF" wi="1403" he="150" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msup><msub><mi>W</mi><mi>js</mi></msub><mn>14</mn></msup></mrow></mfrac><mo>=</mo><msubsup><mi>f</mi><mi>j</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mo>&CenterDot;</mo><mo>)</mo></mrow><msub><mi>y</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>&gamma;</mi><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msup><msub><mrow><mo>&PartialD;</mo><mi>W</mi></mrow><mi>js</mi></msub><mn>14</mn></msup></mfrac><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo><mi>s</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000503821400000023.GIF" wi="1485" he="149" /></maths>其中,η<sub>1</sub>,η<sub>2</sub>,η<sub>3</sub>,η<sub>4</sub>分别为W<sup>11</sup>、W<sup>12</sup>、W<sup>13</sup>、W<sup>14</sup>的学习步长,0≤γ<1为联系单元2的自连接反馈增益因子,δ<sub>i</sub><sup>0</sup>和δ<sub>j</sub><sup>h</sup>为误差信号;δ<sub>i</sub><sup>0</sup>=(y<sub>di</sub>(k)‑y<sub>i</sub>(k)) (12)<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><msub><mi>&delta;</mi><mi>j</mi></msub><mi>h</mi></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msup><msub><mi>&delta;</mi><mi>i</mi></msub><mn>0</mn></msup><msup><msub><mi>W</mi><mi>ij</mi></msub><mn>13</mn></msup><mo>)</mo></mrow><msubsup><mi>f</mi><mi>j</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mo>&CenterDot;</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000503821400000024.GIF" wi="1358" he="131" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><msup><msub><mrow><mo>&PartialD;</mo><mi>W</mi></mrow><mrow><mi>j</mi><mn>1</mn></mrow></msub><mn>11</mn></msup></mfrac><mo>=</mo><msubsup><mi>f</mi><mi>j</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mo>&CenterDot;</mo><mo>)</mo></mrow><msub><mi>x</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>&alpha;</mi><mfrac><mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msup><msub><mrow><mo>&PartialD;</mo><mi>W</mi></mrow><mi>jl</mi></msub><mn>11</mn></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000503821400000025.GIF" wi="1501" he="149" /></maths>所述的OIF Elman网络学习算法采用梯度下降法,即自适应学习速率动量梯度下降反向传播算法,它既能提高网络的训练速率,又能有效抑制网络陷入局部极小点,学习的目的是用网络的实际输出值与输出样本值的差值来修改权值和阀值,使网络输出层的误差平方和最小;S2电池SOH模型输入参数的选择选择电池内阻、电流和温度作为模型的输入参数;确定隐层神经元的个数n<sub>2</sub>和输入层神经元个数n<sub>1</sub>:隐层神经元的个数n<sub>2</sub>和输入层神经元个数n<sub>1</sub>之间有以下近似关系:n<sub>2</sub>=2n<sub>1</sub>+1 (14)由样本数据可知,网络输入层有3个神经元,输出层有1个神经元,故分别将中间层神经元的数目设置为7;S3在Matlab7.1环境下运行对网络输入层和输出层的各单元值进行标准化处理,并使他们处于[0,1]或[‑1,1]中,采用的输入参数的标准化方法为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>x</mi><mo>^</mo></mover><mo>=</mo><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>min</mi></msub></mrow><mrow><msub><mi>x</mi><mi>max</mi></msub><mo>-</mo><msub><mi>x</mi><mi>min</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000503821400000031.GIF" wi="1042" he="128" /></maths>式中,x<sub>max</sub>为该组变量的最大值,x<sub>min</sub>为该组变量的最小值;S4验证电池SOH预测模型的仿真分析选用6组相同型号的电池,采用相同充电制度充满电,在不同电流、不同温度的条件下,对3组电池进行放电试验,获得实测数据100组,选择其中80组数据对该模型进行训练,网络训练过程在Matlab7.1环境下运行,仿真中用改进Elman和Elman两种神经网络建模进行预测,分别得到预测结果,对比分析两者与实测SOH的误差。
地址 510080 广东省广州市越秀区东风东路水均岗8号