发明名称 一种锂电池SOC估算方法
摘要 本发明公开了一种锂电池SOC估算方法,包括:产生sigama点,确定加权系数,UT变换,Kalman滤波框架。本发明所述锂电池SOC估算方法,可以克服现有技术中使用寿命短、安全性差和可靠性低等缺陷,以实现使用寿命长、安全性好和可靠性高的优点。
申请公布号 CN104502852A 申请公布日期 2015.04.08
申请号 CN201410769169.3 申请日期 2014.12.12
申请人 广西科技大学 发明人 刘胜永;李昊
分类号 G01R31/36(2006.01)I 主分类号 G01R31/36(2006.01)I
代理机构 北京中恒高博知识产权代理有限公司 11249 代理人 姜万林
主权项 一种锂电池SOC估算方法,其特征在于,包括:⑴产生sigama点<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&chi;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000011.GIF" wi="1724" he="102" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&chi;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mrow><mo>(</mo><mi>a</mi><msqrt><mi>L</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></msqrt><mo>)</mo></mrow><mi>i</mi></msub><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>L</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000012.GIF" wi="1693" he="126" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&chi;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mrow><mo>(</mo><mi>a</mi><msqrt><mi>L</mi><msub><mover><mi>P</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></msqrt><mo>)</mo></mrow><mi>i</mi></msub><mi>i</mi><mo>=</mo><mi>L</mi><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>2</mn><mi>L</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>18</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000013.GIF" wi="1654" he="122" /></maths>L=n+λ      (3‑19)其中,a是一个较小的正数,通常取a=1;n为系统状态变量个数,针对单体电池,n=2,λ是一比例系数;⑵确定加权系数<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mn>0</mn></msub><mo>=</mo><mfrac><mi>&lambda;</mi><mrow><mi>n</mi><mo>+</mo><mi>&lambda;</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>a</mi><mn>2</mn></msup><mo>+</mo><mi>&beta;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>20</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000014.GIF" wi="1599" he="151" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msup><mi>La</mi><mn>2</mn></msup></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>2</mn><mi>L</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>21</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000015.GIF" wi="1632" he="150" /></maths>λ=a<sup>2</sup>(n+κ)‑n      (3‑22)其中,κ是二阶比例系数,通常取κ=1;β用来合并先验信息,对于高斯白噪声系统,取β=2;⑶UT变换χ<sub>i,k|k‑1</sub>=f(χ<sub>i,k‑1</sub>)      (3‑23)<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>L</mi></mrow></munderover><msub><mi>w</mi><mi>i</mi></msub><msub><mi>&chi;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>24</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000016.GIF" wi="1763" he="159" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mn>2</mn><mi>L</mi></mrow></munderover><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>&chi;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><msub><mi>&chi;</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>+</mo><msub><mi>Q</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>25</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000017.GIF" wi="1721" he="164" /></maths>其中<img file="FDA0000632784250000018.GIF" wi="197" he="103" />是矩阵平方根的第i列,w<sub>i</sub>为权重,Q<sub>k</sub>为过程噪声方差;⑷Kalman滤波框架将UT变换所得的值代入经典的Kalman滤波框架之中便得到UKF算法,递推过程如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>C</mi><mi>k</mi></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>26</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000019.GIF" wi="1752" he="104" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mover><mi>P</mi><mo>^</mo></mover><mi>yy</mi></msub><mo>=</mo><msub><mi>C</mi><mi>k</mi></msub><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><msub><mi>C</mi><mi>k</mi></msub><mi>T</mi></msup><mo>+</mo><msub><mi>R</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>27</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000021.GIF" wi="1729" he="109" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mover><mi>P</mi><mo>^</mo></mover><mi>xy</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><msub><mi>C</mi><mi>k</mi></msub><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>28</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000022.GIF" wi="1751" he="113" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mover><mi>P</mi><mo>^</mo></mover><mi>xy</mi></msub><msubsup><mover><mi>P</mi><mo>^</mo></mover><mi>yy</mi><mo>-</mo></msubsup><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mover><mi>y</mi><mo>^</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>29</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000023.GIF" wi="1739" he="113" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mover><mi>P</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>P</mi><mo>^</mo></mover><mi>xy</mi></msub><msubsup><mover><mi>P</mi><mo>^</mo></mover><mi>yy</mi><mo>-</mo></msubsup><msup><msub><mover><mi>P</mi><mo>^</mo></mover><mi>xy</mi></msub><mo>-</mo></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>30</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632784250000024.GIF" wi="1769" he="115" /></maths>其中<img file="FDA0000632784250000025.GIF" wi="186" he="102" />为状态协方差,C<sub>k</sub>为系统输出矩阵,R<sub>k</sub>为测量噪声方差。
地址 545006 广西壮族自治区柳州市城中区东环大道268号