发明名称 基于相关向量回归的在线预测锂离子电池剩余寿命的方法
摘要 基于相关向量回归的在线预测锂离子电池剩余寿命的方法,属于锂离子电池寿命预测技术领域。它解决了现有锂离子电池采用离线方法预测剩余寿命,预测精度低的问题。它首先选取原始样本,然后进行相空间重构构造训练样本集;再初始化相关向量机RVM模型参数;RVM训练,得到RVM预测模型;得到预测值<img file="DDA00002590473300011.GIF" wi="88" he="46" />将<img file="DDA00002590473300012.GIF" wi="92" he="64" />与y<sub>new</sub>进行比较,若<img file="DDA00002590473300013.GIF" wi="382" he="61" />则构造新的训练集WS=WS∪INS,重新训练RVM,更新RVM预测模型;否则保持RVM预测模型不变;进行递推预测,直到预测值小于失效阈值U时预测完成,从而实现待预测锂离子电池剩余寿命的在线预测。本发明适用于锂离子电池剩余寿命的预测。
申请公布号 CN102968573A 申请公布日期 2013.03.13
申请号 CN201210543701.0 申请日期 2012.12.14
申请人 哈尔滨工业大学 发明人 周建宝;刘大同;马云彤;彭宇;彭喜元
分类号 G06F19/00(2006.01)I 主分类号 G06F19/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 1.一种基于相关同量回归的在线预测锂离于电池剩余寿命的万法,其特征在于:它包括以下步骤:步骤一:选取待预测锂离子电池容量数据IS=(C<sub>1</sub>,C<sub>2</sub>,…C<sub>n</sub>)作为原始样本,C<sub>i</sub>为电池容量,单位为Ah,i=1,2,…,n,n为正整数;进行相空间重构构造训练样本集:设定嵌入维数l=5,迟延d=1,得到训练样本集{(x<sub>1</sub>,y<sub>1</sub>),(x<sub>2</sub>,y<sub>2</sub>),…,(x<sub>n-l</sub>,y<sub>n-l</sub>)},其中x<sub>j</sub>=(C<sub>j</sub>,C<sub>j+1</sub>,…,C<sub>j+l-1</sub>),y<sub>j</sub>=C<sub>j+l</sub>,j=1,2,…n-l,其中x=(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>n-l</sub>)为相关向量机RVM模型输入数据,y=(y<sub>1</sub>,y<sub>2</sub>,…,y<sub>n-l</sub>)为相关向量机RVM模型的输出数据;步骤二:初始化相关向量机RVM模型参数:相关向量机RVM模型的数学表达式为y=Φω+ε,其中ω=(ω<sub>0</sub>,…,ω<sub>n-l</sub>)<sup>T</sup>为模型的权值,ε=(ε<sub>1</sub>,ε<sub>2</sub>,…ε<sub>n-l</sub>)为高斯噪声,且ε<sub>j</sub>~N(0,σ<sup>2</sup>),σ<sup>2</sup>为RVM模型输出数据y的噪声方差,Φ为n×(n+1)的矩阵,且Φ=[φ<sub>1</sub>,φ<sub>2</sub>…φ<sub>n-l</sub>]<sup>T</sup>,φ<sub>j</sub>=[1,K(x<sub>j</sub>,x<sub>1</sub>),…,K(x<sub>j</sub>,x<sub>j</sub>)…K(x<sub>j</sub>,x<sub>n-l</sub>)]K(x<sub>j</sub>,x<sub>n-l</sub>)为核函数:<maths num="0001"><![CDATA[<math><mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>,</mo><msub><mi>x</mi><mrow><mi>n</mi><mo>-</mo><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>n</mi><mo>-</mo><mi>l</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><msup><mi>&eta;</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>η为核参数;设定核参数η=3,最大迭代次数iter=1000,{α<sub>k</sub>}=0.1,k=0,1…n-l,α<sub>k</sub>为权值ω<sub>k</sub>的超参数,σ<sup>2</sup>=var(y)*0.1;步骤三:RVM训练:步骤三一:计算ω的协方差∑和均值μ:∑=(σ<sup>-2</sup>Φ<sup>T</sup>Φ+A)<sup>-1</sup>,μ=σ<sup>-2</sup>∑Φ<sup>T</sup>y<sup>T</sup>其中,∑为(n-l+1)×(n-l+1)的矩阵,μ为(n-l+1)×1的列向量,μ=(μ<sub>0</sub>,μ<sub>1</sub>,…,μ<sub>n-l</sub>)<sup>T</sup>,A=diag(α<sub>0</sub>,α<sub>1</sub>,…,α<sub>n-l</sub>);步骤三二:使用迭代估计法计算得到新的α<sub>k</sub>和σ<sup>2</sup>,记为<img file="FDA00002590473000021.GIF" wi="86" he="54" />和(σ<sup>2</sup>)<sup>new</sup>:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&alpha;</mi><mi>k</mi><mi>new</mi></msubsup><mo>=</mo><mfrac><msub><mi>&gamma;</mi><mi>k</mi></msub><msubsup><mi>&mu;</mi><mi>k</mi><mn>2</mn></msubsup></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msup><mrow><mo>(</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow><mi>new</mi></msup><mo>=</mo><mfrac><mrow><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><mi>&Phi;&mu;</mi><mo>|</mo><mo>|</mo></mrow><mrow><mi>n</mi><mo>-</mo><mi>l</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mi>l</mi></mrow></munderover><msub><mi>&gamma;</mi><mi>k</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中γ<sub>k</sub>=1-α<sub>k</sub>∑<sub>kk</sub>,μ<sub>k</sub>为ω<sub>k</sub>的均值,∑<sub>kk</sub>为协方差∑的第k个对角线元素;设定迭代次数参数L=1;步骤三三:将<img file="FDA00002590473000024.GIF" wi="98" he="61" />和(σ<sup>2</sup>)<sup>new</sup>带入步骤三一的公式中,重复步骤三一和步骤三二,更新μ和∑,令L=L+1,重复该步骤三三,直到L>iter时迭代结束,迭代结束时的(σ<sup>2</sup>)<sup>new</sup>记为<img file="FDA00002590473000025.GIF" wi="123" he="72" />步骤四:步骤三中迭代结束后,删除与ω<sub>k</sub>=0所对应x中的x<sub>j</sub>,剩余x<sub>j</sub>称为相关向量RVs,所有相关向量RVs组成相关向量集IS<sub>RV</sub>,由此获得RVM预测模型<img file="FDA00002590473000026.GIF" wi="281" he="54" />h为预测步长,方差为<img file="FDA00002590473000027.GIF" wi="653" he="66" />并得到工作集WS,WS=IS<sub>RV</sub>;步骤五:将新增样本集INS=[x<sub>new</sub>,y<sub>new</sub>]中x<sub>new</sub>输入到RVM预测模型,得到预测值<img file="FDA00002590473000028.GIF" wi="89" he="46" />将<img file="FDA00002590473000029.GIF" wi="108" he="64" />与y<sub>new</sub>进行比较,若<img file="FDA000025904730000210.GIF" wi="382" he="62" />则构造新的训练集WS=WS∪INS,重新训练RVM,更新RVM预测模型;否则保持RVM预测模型不变;PEB为预测误差限,<maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>new</mi></msub><mo>=</mo><msup><mi>&mu;</mi><mi>T</mi></msup><msub><mi>&phi;</mi><mi>new</mi></msub><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>new</mi><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>&sigma;</mi><mi>MP</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&phi;</mi><mi>new</mi><mi>T</mi></msubsup><mi>&Sigma;</mi><msub><mi>&phi;</mi><mi>new</mi></msub><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><msubsup><mi>&sigma;</mi><mi>new</mi><mn>2</mn></msubsup></math>]]></maths>为预测值的方差;步骤六:重复步骤三至步骤五,进行递推预测,直到预测值小于失效阈值U时预测完成,从而实现待预测锂离子电池剩余寿命的在线预测。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号