发明名称 一种估计电动汽车的动力电池的荷电状态的方法
摘要 本发明涉及电动汽车的动力电池管理系统。为提高估计精度,本发明提出一种估计电动汽车的动力电池的荷电状态的方法,采集动力电池的端电压、充放电电流及电池表面温度;根据其开路电压与荷电状态的对应关系解析得出开路电压模型,根据动力电池的动态电压系统输入矩阵、参数矩阵以及输出矩阵,建立其动态电压模型;对状态观测器进行时间更新和测量更新,得到系统状态的和系统状态估计误差的协方差矩阵在k时刻的先验估计值和后验估计值,循环上述更新操作,直至估计完成,得到动力电池的荷电状态SoC。采用该估计方法可将荷电状态SoC的估计精度提高到5%以内,并在动力电池的荷电状态SoC的设定值偏离其参考值时,将荷电状态SoC迅速收敛至其参考值。
申请公布号 CN104569835A 申请公布日期 2015.04.29
申请号 CN201410784128.1 申请日期 2014.12.16
申请人 北京理工大学 发明人 孙逢春;熊瑞;何洪文
分类号 G01R31/36(2006.01)I 主分类号 G01R31/36(2006.01)I
代理机构 北京市邦道律师事务所 11437 代理人 段君峰;薛艳
主权项 一种估计电动汽车的动力电池的荷电状态的方法,其特征在于,该方法包括如下步骤:步骤1、在所述动力电池充放电的过程中,对该动力电池的端电压U<sub>t</sub>、充放电电流i<sub>L</sub>及电池表面温度T<sub>em</sub>进行采样,采样时间间隔为Δt;步骤2、建立所述动力电池的开路电压模型和动态电压模型:根据所述动力电池的开路电压U<sub>oc</sub>与荷电状态SoC、电池表面温度T<sub>em</sub>三者之间的对应关系解析得出该动力电池的开路电压模型U<sub>oc</sub>(z),且U<sub>oc</sub>(z)=α<sub>0</sub>+α<sub>1</sub>z+α<sub>2</sub>z<sup>2</sup>+α<sub>3</sub>z<sup>3</sup>+α<sub>4</sub>/z+α<sub>5</sub>ln z+α<sub>6</sub>ln(1‑z)其中,z为所述动力电池的荷电状态SoC的表示符号,α<sub>j</sub>(j=0,1,2,……,6)为所述动力电池的开路电压模型的拟合系数,且该拟合系数与电池表面温度T<sub>em</sub>相关;根据所述动力电池的动态电压系统输入矩阵<img file="FDA0000635141440000011.GIF" wi="117" he="67" />参数矩阵θ<sub>n,k</sub>以及输出矩阵y<sub>n,k</sub>,建立所述动力电池的动态电压模型,<img file="FDA0000635141440000012.GIF" wi="1602" he="385" />其中,E<sub>L,k</sub>为所述动力电池的动态电压E<sub>L</sub>在k时刻的值,E<sub>L,k‑1</sub>为所述动力电池的动态电压E<sub>L</sub>在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,且E<sub>L,k</sub>=U<sub>t,k</sub>‑U<sub>oc</sub>(z<sub>k</sub>),i<sub>L,k</sub>为所述动力电池的充放电电流i<sub>L</sub>在k时刻的值,i<sub>L,k‑1</sub>为所述动力电池的充放电电流i<sub>L</sub>在k时刻前一个时刻的值,二者之间的时间间隔为一个采样时间间隔Δt,<img file="FDA0000635141440000022.GIF" wi="477" he="89" />为所述动态电压模型的模型参数,n为所述动力电池的等效电路中RC网络的阶数,且n为正整数,并采用最小二乘法辨识出所述动态电压模型的模型参数,得出所述动力电池的动态电压模型;步骤3、估计所述动力电池的荷电状态SoC:建立所述动力电池的状态方程和观测方程,<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mi>k</mi></msub><mo>=</mo><msub><mi>Ax</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>Bu</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>&omega;</mi><mi>k</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Y</mi><mi>k</mi></msub><mo>=</mo><msub><mi>U</mi><mrow><mi>oc</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>U</mi><mrow><mi>D</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>i</mi><mrow><mi>L</mi><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>R</mi><mi>i</mi></msub><mo>+</mo><msub><mi>&rho;</mi><mi>k</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000635141440000021.GIF" wi="714" he="187" /></maths>其中,x<sub>k</sub>为所述动力电池在k时刻的系统估计状态,且x<sub>k</sub>=[U<sub>D,k</sub>,z<sub>k</sub>]<sup>T</sup>,U<sub>D,k</sub>为所述动力电池在k时刻的极化电压,z<sub>k</sub>为所述动力电池在k时刻的荷电状态,u<sub>k‑1</sub>为所述动力电池在k‑1时刻的系统输入信息,ω<sub>k</sub>为所述动力电池在k时刻的系统状态白噪声,其均值为零,协方差为Q<sub>k</sub>,Y<sub>k</sub>为所述动力电池在k时刻的系统观测矩阵,U<sub>oc,k</sub>为所述动力电池在k时刻的开路电压,ρ<sub>k</sub>为所述动力电池在k时刻的系统参数白噪声,其均值为零,协方差为R<sub>k</sub>,<img file="FDA0000635141440000031.GIF" wi="1357" he="564" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>&Delta;t</mi><mrow><msub><mi>R</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><msub><mi>C</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>R</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>&Delta;t</mi><mrow><msub><mi>R</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub><msub><mi>C</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>R</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>&Delta;t</mi><mrow><msub><mi>R</mi><mi>Dn</mi></msub><msub><mi>C</mi><mi>Dn</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><msub><mi>R</mi><mi>Dn</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mi>&Delta;t</mi><msub><mi>C</mi><mi>n</mi></msub></mfrac></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000635141440000032.GIF" wi="675" he="662" /></maths>对状态观测器中的X,P,Q,R进行初始化设置得到X<sub>0</sub>、P<sub>0</sub>、Q<sub>0</sub>、R<sub>0</sub>,其中,X<sub>0</sub>为所述动力电池的系统状态X的初始值,P<sub>0</sub>为所述动力电池的系统状态估计误差的协方差矩阵P的初始值,Q<sub>0</sub>为所述动力电池的过程噪声Q的初始值,R<sub>0</sub>为所述动力电池的观测噪声的协方差R的初始值;对所述状态观测器进行时间更新,且更新的时间长度为一个采样时间间隔Δt,得到所述系统状态X和所述系统状态估计误差的协方差矩阵P在k时刻的先验估计值<img file="FDA0000635141440000033.GIF" wi="63" he="91" />和<img file="FDA0000635141440000034.GIF" wi="89" he="67" />对所述状态观测器进行测量更新,得到所述系统状态X和所述系统状态估计误差的协方差矩阵P在k时刻的后验估计值<img file="FDA0000635141440000035.GIF" wi="62" he="100" />和<img file="FDA0000635141440000036.GIF" wi="99" he="73" />循环上述更新操作,得到所述动力电池的荷电状态SoC。
地址 100081 北京市海淀区中关村南大街5号