专利名称：：一种电池剩余电量的混合快速估计方法
技术领域：
：本发明属于电池
技术领域：
，涉及一种电池剩余电量的混合快速估计方法。技术背景电池作为备用电源已在通讯、电力系统、军事装备等领域得到了广泛的应用。同传统燃油汽车相比，电动汽车可实现零排放，因此是未来汽车的主要发展方向。在电动汽车中电池直接作为主动能量供给部件，因此其工作状态的好坏直接关系到整个汽车的行驶安全性和运行可靠性。为确保电动汽车中的电池组性能良好，延长电池组使用寿命，须及时、准确地了解电池的运行状态，对电池进行合理有效的管理和控制。电池荷电状态(StateofCharge,以下简称S0C)的精确估算是电池能量管理系统中最核心的技术。电池的SOC无法用一种传感器直接测得，它必须通过对一些其他物理量的测量，并采用一定的数学模型和算法来估计得到。目前常用的电池SOC估计方法有开路电压法、安时法等。开路电压法进行电池S0C估计时电池必须静置较长时间以达到稳定状态，而且只适用于电动汽车在停车状态下的SOC估计，不能满足在线检测要求。安时法易受到电流测量精度的影响，在高温或电流波动剧烈情况下，精度很差。
发明内容本发明的目的就是克服现有技术的不足，提出一种电池剩余电量的混合估计方法，可以适用于所有电池，且估计速度快，估计精度较高。本发明的电池剩余电量混合快速估计方法，具体步骤是步骤(l).用状态方程和观测方程表示电池的各个时刻的荷电状态:状态方程zi+1=/(z""十w=z广2尺观测方程&=g(A,4)+、=《0-—~-《2zt+尺3ln&+《4ln(l—z4)+vt、其中z为电池的荷电状态(StateofCharge,S0C)，即剩余电量；；/,.为电池的放电比例系数，反映的是放电速率、温度、自放电、老化等因素对电池SOC的影响程度；A,为测量时间间隔，^为处理噪声，4为放电电流，2是电池在室温25。c条件下、以1/30倍额定电流的放电速率放电时所能得到的额定总电量，《。、《、《2、&、^为常数，i&《2《3I4f，/>为电池观测模型的参数，是一个列向量，对同类型的电池它们是不变的；i为电池的内阻，^为观测噪声；k为测量时刻。设定处理噪声^的方差、和观测噪声、的方差&分别为2w=10-5，&=10-2;放电比例系数7,的确定方法为(a)将完全充满电的电池以不同放电速率C,.(0<C,、C，C为电池的额定放电电流)恒流放电iV(A^10)次，计算相应放电速率下的电池总电量Q，1S/S7Vo(b)根据最小二乘方法拟合出2,.与C,.间的二次曲线关系，即在最小均方误差准则下求出同时满足Q.=《2+6C,.+c，a》，c为最优系数。(C)在放电电流为iA时，对应的放电比例系数仏为++c此处，最优系数a力，c对于同一类型的电池只需确定一次，确定后可作为已知常数直接用于所有同类型电池的剩余电量估计。电池的内阻i以及常数^、&、A、A、^的确定方法为(d)在室温25。C条件下、以1/30倍额定电流对充满电的电池进行恒定电流放电直至电量耗尽；(e)在放电过程中以时间间隔"测量电池在s时刻的端电压凡，s=0，l,2，...M，其中5=0对应电池充满后的起始放电时刻，^=^对应电池电量耗尽的终止时刻。(f)计算^时刻的剩余电量&=1-WM。(g)记7=y。11_—一z。lnz0ln(l_z0)i丄30130zAlnzjln(l-zj—zMlnzMln(l—zM)《2则;^(H广iTy,也就得到内阻i以及常数K。、&、&、K3、&。对同一类型的电池，这些参数只需确定一次，确定后可作为已知常数直接用于所有同类型电池的剩余电量估计。步骤(2).将电池静置较长时间后，测量初始时刻*=0时电池两端的开路电压t/。，计算电池的初始剩余电量z。z。=A(t/。)。其中A为电池的开路电压到电池剩余电量的映射函数。设定、的方差P。为1(T2。映射函数A的确定方法为(a)将电池充满电，在室温25。C条件下、以%。(C为定电流)的放电速率对电池进行恒流放电r.Af(r为放电次数，r=l,2,...)时间后，将电池开路。在电池静置较长时间后，测量电池两端的开路电压R。然后在室温25。c条件下、以%。(c为额定电流)的放电速率对电池进行恒流放电直至完全放电，测得此时的放电时间为;，计算得到相应的电池剩余总电量^为。(b)选取多项式函数A(x)-"。+fl,;c+a2;c2+…+Z，(1《《^5)，根据所测得的t/,及计算得到的^，采用最小二乘法拟合得出多项式函数的最优系数K^，…,)，即得到映射函数/。对同一类型的电池，映射函数A只需确定一次。步骤(3).电池实际工作过程中，在时刻A"，2,3，…，测量电池两端的端电压^及电池的供电电流4。步骤(4).根据所测得的^及4，交替使用扩展卡尔曼滤波方法及安时法进行电池剩余电量的估计，直至电池剩余电量少于严/。，其中15Si^20;其中扩展卡尔曼滤波方法所连续利用的测量数据点数为S点(S>200)，安时法所连续利用的测量数据点数为i点(i>S)。在时刻/o利用扩展卡尔曼滤波方法进行电池剩余电量估计的过程为(a)计算状态方程对状态变量的偏导数^_1:(b)根据状态方程计算状态的估计值z::=Vl-;7,V^/"(C)计算状态估计值的协方差/T:g、dK+、-g—1+SW(d)根据观测方程计算观测的估计值j^A=K。_及4_《/z:_《2z:+《3Inz:+《4ln(l—《)(e)计算观测方程对状态变量的偏导数&:<formula>formulaseeoriginaldocumentpage10</formula>(f)计算卡尔曼滤波增益i^(g)计算荷电状态^及其方差g:在时刻it，利用安时法进行电池剩余电量估计的方法为步骤(5)当电池剩余电量少于尸％时，一直使用上述步骤(4)中的扩展卡尔曼滤波方法进行电池剩余电量的估计。根据本发明的第一方面，公开了一种电池剩余电量的混合快速估计方法中扩展卡尔曼滤波方法所采用的状态方程和观测方程。其中观测方程中电池模型参数通过最小二乘法确定。根据本发明的第二方面，公开了一种电池剩余电量的混合快速估计方法中初始剩余电量的估计方法。初始剩余电量采用开路电压法确定。根据本发明的第三方面，公开了一种电池剩余电量的混合快速估计方法所依赖的测量物理量。主要包括电池的端电压和供电电流。根据本发明的第四和第五方面，公开了一种电池剩余电量的混合快速估计方法的具体流程。当电池剩余电量大于尸％时，采用扩展卡尔曼滤波方法和安时法交替进行的方式进行电池剩余电量的估计；当电池剩余电量小于尸％时，采用扩展卡尔曼滤波方法进行电池剩余电量的估计。本发明可以方便地进行电池SOC的快速估计，该方法收敛速度快，估计精度高，而且适用于各种电池soc的快速估计。具体实施方式电池剩余电量的混合快速估计方法的具体歩骤是步骤(l).用状态方程和观测方程表示电池的各个时刻的荷电状态状态方程Zw=,4)+w=&4,观测方程=g(^,4)+、=K。一-~^—《2&+《3ln&+《4ln(l_2J+、其中z为电池的荷电状态(StateofCharge,S0C)，即剩余电量；77,.为电池的放电比例系数，反映的是放电速率、温度、自放电、老化等因素对电池SOC的影响程度；"是测量时间间隔，^为处理噪声，2是电池在室温25°c条件下、以1/30倍额定电流的放电速率放电时所能得到的额定总电量。I。、《、&、&、《4为常数，p=[《。w&尺2&x4;T，p为电池观测模型的参数，是一个列向量，对同类型的电池它们是不变的；i为电池的内阻，、为观测噪声。设定处理噪声^的方差、和观测噪声、的方差Sv分别为Sw=10-52V=10-2。放电比例系数7/i的确定方法为(a)将完全充满电的电池以不同放电速率C,.(0<C,、C，C为电池的额定放电电流)恒流放电7V(iV〉10)次，计算相应放电速率下的电池总电量Q，W《iV。(b)根据最小二乘方法拟合出Q与C,.间的二次曲线关系，即在最小均方误差准则下求出同时满足Q,^C,2+M:,.+c，a,&,c为最优系数。(C)在放电电流为/4时，对应的放电比例系数7,.为7,=.25".-flZt++c此处，最优系数",6,c对于同一类型的电池只需确定一次，确定后可作为已知常数直接用于所有同类型电池的剩余电量估计。电池的内阻i以及常数K。、《、A、&、A的确定方法为.-(d)在室温25。c条件下、以1/30倍额定电流对充满电的电池进行恒定电流放电直至电量耗尽；(e)在放电过程中以时间间隔A,测量电池在s时刻的端电压凡，s=0,l,2,...M，其中^=0对应电池充满后的起始放电时刻，^M对应电池电量耗尽的终止时刻。(f)计算s时刻的剩余电量z,l-WM。1(g)记r少o1:一一—z0lnz。ln(l_z0)i—二30i丄30-ZjIn^ln(l-5)zMlnzMln(l—zM)《《2则^=(^//)-W7，也就得到内阻及以及常数尺。、&、《2、&、《4。对同一类型的电池，这些参数只需确定一次，确定后可作为已知常数直接用于所有同类型电池的剩余电量估计。步骤(2).将电池静置较长时间后，测量初始时刻it=0时电池两端的开路电压t/。，计算电池的初始剩余电量z。0)=A(t/。)。其中A为电池的开路电压到电池剩余电量的映射函数。设定z。的方差A为10-2。映射函数A的确定方法为(h)将电池充满电，在室温25。c条件下、以%0(c为额定电流)的放电速率对电池进行恒流放电r"(r为放电次数，r=l,2,...)时间后，将电池开路。在电池静置较长时间后，测量此时电池两端的开路电压R。然后在室温25°c条件下、以%。(c为额定电流)的放电速率对电池进行恒流放电直至完全放电，测得此时的放电时间为？；，计算得到相应的电池剩余总电量^为'30""(0选取合适的多项式函数/^)=+^+"2义2+."+浐，(B《S5)，根据所测得的R及^，采用最小二乘法拟合得出多项式函数的最优系数K^p…,aj，即得到映射函数L对同一类型的电池，映射函数A只需确定一次。步骤(3).电池实际工作过程中，在时刻&=1,2,3,."，测量电池两端的端电压A及电池的供电电流"步骤(4).根据所测得的A及^，交替使用扩展卡尔曼滤波方法及安时法进行电池剩余电量的估计，直至电池剩余电量少于尸％;其中扩展卡尔曼滤波方法所连续利用的测量数据点数为S点(^>200)，安时法所连续利用的测量数据点数为/点。在时刻^，利用扩展卡尔曼滤波方法进行电池剩余电量估计的过程为(a)计算状态方程对状态变量的偏导数^=1(b)根据状态方程计算状态的估计值^:z:=Zt—(c)计算状态估计值的协方差P。&"w^《+、^—1+SV(d)根据观测方程计算观测的估计值义<formula>formulaseeoriginaldocumentpage14</formula>(e)计算观测方程对状态变量的偏导数C,:(f)计算卡尔曼滤波增益i^(g)计算状态^及其方差g:在时刻A，利用安时法进行电池剩余电量估计的方法为步骤(5)当电池剩余电量少于尸％时，一直使用扩展卡尔曼滤波方法进行电池剩余电量的估计。整个估计过程是在线完成的，即在电池实际工作过程中同步完成各时刻电池剩余电量的估计。本发明方法混合应用了开路电压法、扩展卡尔曼滤波方法及安时法进行电池剩余电量的估计。在整个电池剩余电量估计过程中，较长时间内使用安时法进行电池剩余电量的估计，主要利用的是安时法快速估计的特点；另一方面，扩展卡尔曼滤波方法和安时法交替进行，可以利用扩展卡尔曼滤波方法的高精确性对安时法的估计结果进行修正，从而保证了电池剩余电量的估计精度。在剩余电量较低时直接采用扩展卡尔曼滤波方法进行剩余电量的估计进一步保证了估计精度。权利要求1、一种电池剩余电量的混合快速估计方法，其特征在于该方法的具体步骤是步骤(1)用状态方程和观测方程表示电池的各个时刻的荷电状态状态方程<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><msub><mi>i</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mi>k</mi></msub><mo>=</mo><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><mrow><mo>(</mo><mfrac><mrow><msub><mi>η</mi><mi>i</mi></msub><mi>Δt</mi></mrow><msub><mi>Q</mi><mi>n</mi></msub></mfrac><mo>)</mo></mrow><msub><mi>i</mi><mi>k</mi></msub><mo>+</mo><msub><mi>w</mi><mi>k</mi></msub></mrow>]]></math>id="icf0001"file="A2009101009780002C1.tif"wi="67"he="11"top="56"left="58"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>观测方程<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>y</mi><mi>k</mi></msub><mo>=</mo><mi>g</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><msub><mi>i</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>v</mi><mi>k</mi></msub><mo>=</mo><msub><mi>K</mi><mn>0</mn></msub><mo>-</mo><mi>R</mi><msub><mi>i</mi><mi>k</mi></msub><mo>-</mo><mfrac><msub><mi>K</mi><mn>1</mn></msub><msub><mi>z</mi><mi>k</mi></msub></mfrac><mo>-</mo><msub><mi>K</mi><mn>2</mn></msub><msub><mi>z</mi><mi>k</mi></msub><mo>+</mo><msub><mi>K</mi><mn>3</mn></msub><mi>ln</mi><msub><mi>z</mi><mi>k</mi></msub><mo>+</mo><msub><mi>K</mi><mn>4</mn></msub><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>z</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>v</mi><mi>k</mi></msub></mrow>]]></math>id="icf0002"file="A2009101009780002C2.tif"wi="117"he="10"top="74"left="58"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中z为电池的荷电状态，即剩余电量；ηi为电池的放电比例系数；Δt是测量时间间隔，wk为处理噪声；K0、K1、K2、K3、K4为常数，p＝[K0RK1K2K3K4]T，p为电池观测模型的参数，是一个列向量；R为电池的内阻，vk为观测噪声；设定处理噪声wk的方差∑w和观测噪声vk的方差∑v分别为∑w＝10-5∑v＝10-2.放电比例系数ηi的确定方法为(a)将完全充满电的电池以不同放电速率Ci恒流放电N(N＞10)次，，其中0＜Ci≤C，C为电池的额定放电电流；计算相应放电速率下的电池总电量Qi，1≤i≤N；(b)根据最小二乘方法拟合出Qi与Ci间的二次曲线关系，即在最小均方误差准则下求出同时满足<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><mi>a</mi><msubsup><mi>C</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mi>b</mi><msub><mi>C</mi><mi>i</mi></msub><mo>+</mo><mi>c</mi><mo>,</mo></mrow>]]></math>id="icf0003"file="A2009101009780002C3.tif"wi="33"he="4"top="190"left="78"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>a，b，c为最优系数；(c)在放电电流为ik时，对应的放电比例系数ηi为<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msub><mi>η</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>Q</mi><mi>n</mi></msub><mrow><msubsup><mi>ai</mi><mi>k</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>bi</mi><mi>k</mi></msub><mo>+</mo><mi>c</mi></mrow></mfrac></mrow>]]></math></maths>此处，最优系数a，b，c对于同一类型的电池只需确定一次，确定后可作为已知常数直接用于所有同类型电池的剩余电量估计；电池的内阻R以及常数K0、K1、K2、K3、K4的确定方法为(d)在室温25℃条件下、以1/30倍额定电流对充满电的电池进行恒定电流放电直至电量耗尽；(e)在放电过程中以时间间隔Δt测量电池在s时刻的端电压ys，s＝0，1，2，...M，其中s＝0对应电池充满后的起始放电时刻，s＝M对应电池电量耗尽的终止时刻；(f)计算s时刻的剩余电量zs＝1-s/M；(g)记<mathsid="math0005"num="0005"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>y</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>M</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math>id="icf0005"file="A2009101009780003C1.tif"wi="18"he="24"top="89"left="46"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0006"num="0006"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mi>C</mi><mn>30</mn></mfrac></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><msub><mi>z</mi><mn>0</mn></msub></mfrac></mtd><mtd><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub></mtd><mtd><mi>ln</mi><msub><mi>z</mi><mn>0</mn></msub></mtd><mtd><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mi>C</mi><mn>30</mn></mfrac></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><msub><mi>z</mi><mn>1</mn></msub></mfrac></mtd><mtd><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub></mtd><mtd><mi>ln</mi><msub><mi>z</mi><mn>1</mn></msub></mtd><mtd><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mi>C</mi><mn>30</mn></mfrac></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><msub><mi>z</mi><mi>M</mi></msub></mfrac></mtd><mtd><mo>-</mo><msub><mi>z</mi><mi>M</mi></msub></mtd><mtd><mi>ln</mi><msub><mi>z</mi><mi>M</mi></msub></mtd><mtd><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>z</mi><mi>M</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math>id="icf0006"file="A2009101009780003C2.tif"wi="84"he="42"top="80"left="68"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths><mathsid="math0007"num="0007"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>K</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mi>R</mi></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>4</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math>id="icf0007"file="A2009101009780003C3.tif"wi="15"he="37"top="82"left="159"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>则p＝(HTH)-1HTY，也就得到内阻R以及常数K0、K1、K2、K3、K4；对同一类型的电池，这些参数只需确定一次，确定后可作为已知常数直接用于所有同类型电池的剩余电量估计；步骤(2)将电池静置较长时间后，测量初始时刻k＝0时电池两端的开路电压U0，计算电池的初始剩余电量z0Q0＝h(U0)；其中h为电池的开路电压到电池剩余电量的映射函数；设定z0的方差P0为10-2；映射函数h的确定方法为(h)将电池充满电，在室温25℃条件下、以C/30放电速率对电池进行恒流放电r·Δt时间后，将电池开路，其中C为额定电流，r＝1，2，...L；在电池静置较长时间后，测量此时电池两端的开路电压Ur；然后在室温25℃条件下、以C/20的放电速率对电池进行恒流放电直至完全放电，测得此时的放电时间为Tr，计算得到相应的电池剩余总电量zr为<mathsid="math0008"num="0008"><math><![CDATA[<mrow><msub><mi>z</mi><mi>r</mi></msub><mo>=</mo><mfrac><mi>C</mi><mn>30</mn></mfrac><msub><mi>T</mi><mi>r</mi></msub><mo>/</mo><msub><mi>Q</mi><mi>n</mi></msub><mo>,</mo></mrow>]]></math>id="icf0008"file="A2009101009780003C4.tif"wi="25"he="9"top="237"left="134"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其中Qn是电池在室温25℃条件下、以1/30倍额定电流的放电速率放电时所能得到的额定总电量；(i)选取合适的多项式函数h(x)＝a0+a1x+a2x2+…+aKxK，1≤K≤5，根据所测得的Ur及zr，采用最小二乘法拟合得出多项式函数的最优系数{a0，a1，…，aK}，即得到映射函数h；对同一类型的电池，映射函数h只需确定一次；步骤(3)电池实际工作过程中，在时刻k＝1，2，3，…，测量电池两端的端电压yk及电池的供电电流ik；步骤(4)根据所测得的yk及ik，交替使用扩展卡尔曼滤波方法及安时法进行电池剩余电量的估计，直至电池剩余电量少于P％；其中扩展卡尔曼滤波方法所连续利用的测量数据点数为S点，S＞200，安时法所连续利用的测量数据点数为R点；在时刻k，利用扩展卡尔曼滤波方法进行电池剩余电量估计的过程为(a)计算状态方程对状态变量的偏导数Ak-1<mathsid="math0009"num="0009"><math><![CDATA[<mrow><msub><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><msub><mi>i</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow><msub><mrow><mo>∂</mo><mi>z</mi></mrow><mi>k</mi></msub></mfrac><msub><mo>|</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>=</mo><msubsup><mover><mi>z</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>+</mo></msubsup></mrow></msub><mo>=</mo><mn>1</mn></mrow>]]></math></maths>(b)根据状态方程计算状态的估计值zk-<mathsid="math0010"num="0010"><math><![CDATA[<mrow><msubsup><mi>z</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>z</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>η</mi><mi>i</mi></msub><msub><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mi>Δt</mi><mo>/</mo><msub><mi>Q</mi><mi>n</mi></msub></mrow>]]></math>id="icf0010"file="A2009101009780004C2.tif"wi="37"he="4"top="159"left="121"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>(c)计算状态估计值的协方差Pk-<mathsid="math0011"num="0011"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msubsup><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><mo>+</mo><msub><mi>Σ</mi><mi>w</mi></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>Σ</mi><mi>w</mi></msub></mrow>]]></math>id="icf0011"file="A2009101009780004C3.tif"wi="56"he="4"top="170"left="106"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>(d)根据观测方程计算观测的估计值id="icf0012"file="A2009101009780004C4.tif"wi="5"he="4"top="181"left="111"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/><mathsid="math0012"num="0012"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msub><mi>K</mi><mn>0</mn></msub><mo>-</mo><mi>R</mi><msub><mi>i</mi><mi>k</mi></msub><mo>-</mo><msub><mi>K</mi><mn>1</mn></msub><mo>/</mo><msubsup><mi>z</mi><mi>k</mi><mo>-</mo></msubsup><mo>-</mo><msub><mi>K</mi><mn>2</mn></msub><msubsup><mi>z</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>K</mi><mn>3</mn></msub><mi>ln</mi><msubsup><mi>z</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>K</mi><mn>4</mn></msub><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>z</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow></mrow>]]></math></maths>(e)计算观测方程对状态变量的偏导数Ck<mathsid="math0013"num="0013"><math><![CDATA[<mrow><msub><mi>C</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>h</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><msub><mi>i</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>z</mi><mi>k</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>=</mo><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup></mrow></msub><mo>=</mo><mfrac><msub><mi>K</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>-</mo><msub><mi>K</mi><mn>2</mn></msub><mo>+</mo><mfrac><msub><mi>K</mi><mn>3</mn></msub><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup></mfrac><mo>-</mo><mfrac><msub><mi>K</mi><mn>4</mn></msub><mrow><mn>1</mn><mo>-</mo><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>k</mi><mo>-</mo></msubsup></mrow></mfrac></mrow>]]></math></maths>(f)计算卡尔曼滤波增益Lk<mathsid="math0014"num="0014"><math><![CDATA[<mrow><msub><mi>L</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>C</mi><mi>k</mi><mi>T</mi></msubsup></mrow><mrow><msub><mi>C</mi><mi>k</mi></msub><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msubsup><mi>C</mi><mi>k</mi><mi>T</mi></msubsup><mo>+</mo><msub><mi>Σ</mi><mi>v</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msub><mi>C</mi><mi>k</mi></msub></mrow><mrow><msup><mrow><mo>(</mo><msub><mi>C</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>Σ</mi><mi>v</mi></msub></mrow></mfrac></mrow>]]></math></maths>(g)计算状态zk及其方差Pk<mathsid="math0015"num="0015"><math><![CDATA[<mrow><msub><mi>z</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>z</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>L</mi><mi>k</mi></msub><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></mrow>]]></math></maths><mathsid="math0016"num="0016"><math><![CDATA[<mrow><msub><mi>P</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>L</mi><mi>k</mi></msub><msub><mi>C</mi><mi>k</mi></msub><mo>)</mo></mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup></mrow>]]></math></maths>在时刻k，利用安时法进行电池剩余电量估计的方法为zk＝zk-1-ηiik-1Δt/Qn；步骤(5)当电池剩余电量少于P％时，一直使用上述步骤(4)中的扩展卡尔曼滤波方法进行电池剩余电量的估计。全文摘要本发明涉及一种电池剩余电量的混合快速估计方法。现有方法不能满足在线检测要求，并且精度差。本发明方法首先开路电压法估计出电池的初始剩余电量，然后交替使用扩展卡尔曼滤波方法和安时法进行电池剩余电量的估计，当剩余电量较低时直接使用扩展卡尔曼滤波方法进行剩余电量的估计。本发明方法可以方便地进行电池剩余电量的快速估计，收敛速度快、估计精度较高。文档编号G01R31/36GK101625397SQ200910100978公开日2010年1月13日申请日期2009年8月6日优先权日2009年8月6日发明者何志伟,杰徐,毓曾,高明煜,黄继业申请人:杭州电子科技大学