确定电池荷电状态的方法及装置，存储介质及电子设备与流程

1.本公开涉及电池技术领域，具体地，涉及一种确定电池荷电状态的方法及装置，存储介质及电子设备。


背景技术：

2.电动汽车作为一种新能源汽车，具有降低石油消耗、低污染、低噪声等优点，其被认为是能源危机问题和环境恶化问题的重要解决途径。电池作为电动汽车的动力来源，其荷电状态的准确估计不仅有助于提升电池系统的均衡控制效率和电动汽车的能量管理效率，更攸关动态工况下电动汽车的整车安全。
3.相关技术中，利用电池模型与观测器进行组合，对电池的荷电状态进行估计。在对电池荷电状态进行估计时，由于采样元件校准不合格、采样元件老化等原因，电池数据的测量精度难以得到保证，从而容易导致等效电路模型出现误差，进而导致电池荷电状态的估计出现误差，不利于车辆的安全运行和高效管理。并且，电池模型与观测器组合时，可能加重电池管理系统的运行负担。


技术实现要素：

4.本公开的目的是提供一种确定电池荷电状态的方法及装置，存储介质及电子设备，以解决难以合理确定电池荷电状态的问题。
5.为了实现上述目的，本公开第一方面提供一种确定电池荷电状态的方法，所述方法包括：
6.获取电池的状态数据，所述状态数据包括电流数据、电压数据、温度数据；
7.根据所述电池的等效电路模型的元件参数值、所述状态数据、电池特性数据以及观测器初始参数，确定每类观测器的soc估计值和修正量；
8.根据预设修正边界值以及每类观测器的所述修正量，确定每类观测器的权重信息；
9.基于所述权重信息以及所述每类观测器的soc估计值，确定所述电池的soc目标值。
10.本公开第二方面提供一种确定电池荷电状态的装置，所述装置包括：
11.获取模块，用于获取电池的状态数据，所述状态数据包括电流数据、电压数据、温度数据；
12.第一确定模块，用于根据所述电池的等效电路模型的元件参数值、所述状态数据、电池特性数据以及观测器初始参数，确定每类观测器的soc估计值和修正量；
13.第二确定模块，用于根据预设修正边界值以及每类观测器的所述修正量，确定每类观测器的权重信息；
14.第三确定模块，用于基于所述权重信息以及所述每类观测器的soc估计值，确定所述电池的soc目标值。
15.本公开第三方面提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述中任一项所述方法的步骤。
16.本公开第四方面提供一种电子设备，包括：
17.存储器，其上存储有计算机程序；
18.处理器，用于执行所述存储器中的所述计算机程序，以实现上述任一项所述方法的步骤。
19.通过上述技术方案，至少可以达到以下有益效果：
20.通过最小二乘法rls预估模型确定所述等效电路模型中的元件参数值。并且，在确定所述电池的等效电路模型中的元件参数值时，还考虑了至少包括电压数据的采样误差因子和/或电流数据的采样误差因子的误差信息，从而能够降低采样误差所带来的影响，进而能够提升确定的所述等效电路模型中的元件参数值的准确度，最终起到提升所确定的等效电路模型的精度的效果。进一步地，基于不同类型的观测器确定的电池soc值，确定电池soc目标值，进而提升确定电池荷电状态的合理性，保证车辆的高效管理和可靠运行。
21.本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
22.附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：
23.图1是本公开一示例性实施例所示出的一种确定电池荷电状态的方法的流程图。
24.图2是本公开一示例性实施例所示出的另一种确定电池荷电状态的方法的流程图。
25.图3是本公开一示例性实施例所示出的一种l阶初始等效电路模型的示意图。
26.图4是本公开一示例性实施例所示出的另一种l阶初始等效电路模型的示意图。
27.图5是本公开一示例性实施例所示出的另一种确定电池荷电状态的方法的流程图。
28.图6是本公开一示例性实施例所示出的一种确定等效电路模型的装置的框图。
29.图7是本公开一示例性实施例所示出的一种确定等效电路模型的电子设备的框图。
具体实施方式
30.以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本公开，并不用于限制本公开。
31.需要说明的是，本公开说明书和权利要求书以及附图中的术语“第一”、“第二”等用于区别类似的对象，而不必理解为描述特定的顺序或先后次序。
32.在介绍本公开所提供的确定电池荷电状态的方法、装置、存储介质及电子设备之前，首先对本公开各实施例的应用场景进行介绍。本公开的各实施例可以用于确定电池的荷电状态，所述电池例如可以是三元锂电池、磷酸铁锂电池等等。
33.以电动汽车为例，电池作为电动汽车的动力来源，其状态的准确估计不仅有助于提升电池系统的均衡控制效率和电动汽车的能量管理效率，还攸关动态工况下电动汽车的
整车安全。相关技术中，可以通过确定相应的等效电路模型进而分析所述电池的状态，这些状态例如可以是电池荷电状态(state of charge，简称soc)、电池能量状态(state of energy，简称soe)、电池功率状态(state of power，简称sop)、电池健康状态(state of health，简称soh)等，利用多个电池模型分别与观测器进行组合，或者一个电池模型与一个观测器进行组合的方式，确定电池荷电状态。
34.申请人发现，由于采样元件校准不合格、采样元件老化等原因，电池数据的测量结果中还可能包括对应的误差，从而导致相关技术中辨识出的等效电路模型容易出现误差，进而导致电池状态的估计出现误差，不利于车辆的安全运行和高效管理。举例来讲，随着电动车辆的使用，bms(battery management system，电池管理系统)的采样器件不断老化，其测量偏置将再次出现，造成bms的测量噪声不再为白噪声，而变为有色噪声，最终导致辨识出的等效电路模型出现误差，降低电池状态的估计精度。并且，利用不同模型分别与观测器进行组合，无疑加重了车载bms的运算压力。而利用一种电池模型和一种观测器技术的估计电池soc，其精度和收敛速度无法满足不同工况和环境条件的应用需求。
35.为此，本公开提供一种确定电池荷电状态的方法，参照图1所示出的一种确定电池荷电状态的方法的流程示意图，所述方法包括：
36.s11、获取电池的状态数据，所述状态数据包括电流数据、电压数据、温度数据。
37.其中，所述状态数据还可以包括所述电池的容量数据、开路电压-荷电状态曲线等等。
38.以电动汽车为例，在具体实施时可以通过bms来直接和/或间接地获取电池的所述状态数据。例如，bms可以通过电流传感器来直接获取所述电池的电流数据。又或者，bms可以通过温度传感器来获取所述电池的温度数据。在一些实施例中，bms也可以通过相应的数据接口来间接地获取所述状态数据，例如，通过数据接口来获取存储器中所存储的所述电池的开路电压-荷电状态曲线信息。
39.s12、根据所述电池的等效电路模型的元件参数值、所述状态数据、电池特性数据以及观测器初始参数，确定每类观测器的soc估计值和修正量。
40.在一种可能实现的方式中，在步骤s12中，所述等效电路模型的元件参数值是通过如下方式得到的：
41.基于所述电池的离线测试确定等效电路模型，以及所述等效电路模型的元件参数值与荷电状态和电池温度的函数关系式；
42.根据所述温度数据和所述电池前一时刻的荷电状态数据通过所述函数关系式确定所述元件参数值，其中，所述等效电路模型为rc电路模型或者rv电路模型。
43.在另一种可能实现的方式中，在步骤s12中，所述等效电路模型的元件参数值是通过如下方式得到的：
44.基于所述电池的等效电路模型、误差信息、电池特性数据以及所述状态数据，通过最小二乘法rls预估模型，确定等效电路模型中的元件参数值；
45.其中，所述等效电路模型是基于对所述电池进行离线测试得到的，所述电池特性数据包括不同温度下的开路电压-荷电状态曲线以及滞回电压-荷电状态曲线，其中，所述等效电路模型为rc电路模型。
46.可选地，通过上述两种方式确定等效电路模型的元件参数值时，所述等效电路模
型可以为rc电路模型，所述步骤s11包括：
47.获取电池管理系统bms采集的所述状态数据，所述状态数据还包括所述电池的温度数据、所述电池在不同温度下的开路电压-荷电状态曲线以及滞回电压-荷电状态曲线。
48.例如，bms可以通过电流传感器来直接获取所述电池的电流数据。又或者，bms可以通过温度传感器来获取所述电池的温度数据。在一些实施例中，bms也可以通过相应的数据接口来间接地获取所述状态数据，例如，通过数据接口来获取存储器中所存储的所述电池的开路电压-荷电状态曲线信息以及滞回电压-荷电状态曲线信息。
49.所述步骤s12，包括：
50.根据所述温度数据、当前时刻电池的荷电状态数据、所述开路电压-荷电状态曲线以及滞回电压-荷电状态曲线确定目标开路电压以及目标滞回电压。
51.应当理解，在经过离线测试获取所述电池在不同温度下的开路电压-荷电状态曲线以及滞回电压-荷电状态曲线后，可以根据所述bms获取到的电池当前的soc信息、所述开路电压-荷电状态曲线以及滞回电压-荷电状态曲线确定当前时刻的目标开路电压以及目标滞回电压。
52.这样，在步骤s12中，可以基于所述等效电路模型、所述误差信息、所述状态数据、所述目标开路电压以及所述目标滞回电压通过rls预估模型，确定所述等效电路模型中的元件参数值，从而起到在线辨识所述电池的电池模型参数值的效果。同时，由于辨识过程中还考虑了bms采样器件的采样误差，从而能够进一步地提升确定的所述等效电路模型中的元件参数值的准确度，进而能够提升所述电池等效模型的精度。
53.所述误差信息可以包括电压数据的采样误差因子，用于描述采样元件采集的电池电压数据与所述电池的真实电压数据之间的差异。例如，采样元件采集的电池电压数据可以为：其中，u(k)为第k时刻所述电池的电压真值，ε1为电压数据的采样误差因子。
54.这样，上述技术方案在确定所述电池的等效电路模型中的元件参数值时，还考虑了电压数据的采样误差因子，从而能够降低电压采样误差所带来的影响，进而能够提升确定的所述等效电路模型中的元件参数值的准确度，最终起到提升所确定的等效电路模型的精度的效果。
55.此外，在一些实施例中，所述误差信息也可以包括电流数据的采样误差因子，用于描述采样元件采集的电池电流数据与所述电池的真实电流数据之间的差异。例如，采样元件采集的电池电流数据可以为：其中，i(k)为第k时刻所述电池的电流真值，ε2为电流数据的采样误差因子。
56.这样，上述技术方案在确定所述电池的等效电路模型中的元件参数值时，还考虑了电流数据的采样误差因子，从而能够降低电流采样误差所带来的影响，进而能够提升确定的所述等效电路模型中的元件参数值的准确度，最终起到提升所确定的等效电路模型的精度的效果。
57.值得注意的是，在一些实施例中，所述误差信息也可以同时包括所述电压数据的采样误差因子以及电流数据的采样误差因子。在这种情况下，所述步骤s12包括：
58.基于所述电池的等效电路模型的阶数、所述电压数据的采样误差因子、所述电流
数据的采样误差因子以及所述状态数据，通过最小二乘法rls预估模型，确定所述等效电路模型中的元件参数值。
59.采用上述技术方案，在确定所述电池的等效电路模型中的元件参数值时，还考虑了所述电压数据的采样误差因子以及所述电流数据的采样误差因子，从而能够进一步地降低采样误差所带来的影响，进而能够提升确定的所述等效电路模型中的元件参数值的准确度，最终起到提升所确定的等效电路模型的精度的效果。
60.值得说明的是，在一些实施例中，所述误差信息还可以包括所述电压数据和所述电流数据的采样时间差以及电池开路电压的误差。也就是说，在具体实施时，所述采样误差可以包括电压数据的采样误差因子、电流数据的采样误差因子、所述电压数据和所述电流数据的采样时间差以及电池开路电压的误差中的一种或多种，本公开对此不作限定。
61.在一种可能的实施方式中，参照图2所示出的一种确定等效电路模型的流程图，如图所示：
62.在步骤s21中，基于电池离线测试，获得所述电池的初始属性信息，所述初始属性信息包括所述电池的开路电压-荷电状态曲线、滞回电压-荷电状态曲线。
63.其中，所述初始属性信息还可以包括电池容量信息，电池模型参数的协方差初始值等等，所述离线测试可以包括容量测试，脉冲测试以及典型工况测试。
64.具体的，针对所述容量测试，在一实施例中，所述容量测试包括：
65.(1)调整温度为25℃，将电池以电池生产商建议的容量测试电流值(例如1c)放电至电压下限，并静置30分钟；
66.(2)将电池以电池生产商建议的容量测试电流值(例如1c)充电至电压上限(例如4.25v)后，转为恒压充电(恒压值可以为电池生产商建议数值，例如4.25v)，并静置30分钟；
67.(3)统计(1)和(2)分别累计的容量值，重复上述(1)和(2)直至相邻循环之间的容量值差异小于0.1ah，将此时的容量值记为电池容量q
m
。
68.进一步的，在获得所述电池容量q
m
之后，可以对所述电池进行脉冲测试，所述脉冲测试可以包含充电和放电两个部分，其中充电和放电部分均为20组脉冲组合序列。
69.举例来讲，充电部分的前18组脉冲组合序列中，每组脉冲组合序列可以包括：
70.(1)调整温度为25℃，恒流充电脉冲(幅值为1c)，直至脉冲累积安时变化大于等于q
m
的5％；
71.(2)2小时静置；
72.(3)10秒钟恒流充电脉冲(幅值为0.5c)；
73.(4)40秒钟静置；
74.(5)10秒钟恒流放电脉冲(幅值为0.1c)；
75.(6)将0.1c分别替换为0.5c、1c、2c、3c、4c、5c后，重复(3)至(5)；
76.(7)静置24小时；
77.(8)调整温度至55℃，静置2h；
78.(9)将55℃分别替换为40℃、25℃、10℃、0℃、-10℃、-20℃，-30℃，重复(8)；
79.(10)调整温度为25℃；
80.(11)静置24小时。
81.充电部分的最后2组脉冲组合序列中，每组脉冲包括：
82.(1)恒流-恒压脉冲，直至脉冲累积安时变化大于等于q
m
的5％；
83.(2)2小时静置；
84.(3)10秒钟恒流充电脉冲(幅值为0.5c)；
85.(4)40秒钟静置；
86.(5)10秒钟恒流放电脉冲(幅值为0.5c)；
87.(6)将0.5c分别替换为1c、2c、3c、4c、5c后，重复(3)至(5)；
88.(7)静置24小时；
89.(8)调整温度至55℃，静置2h；
90.(9)将55℃分别替换为40℃、25℃、10℃、0℃、-10℃、-20℃，-30℃，重复(8)；
91.(10)调整温度为25℃；
92.(11)静置24小时。
93.此外，针对放电部分，放电部分的20组脉冲组合序列中，每组脉冲可以包括：
94.(1)恒流放电脉冲(幅值为1c)，直至脉冲累积大于或等于q
m
的5％；
95.(2)2小时静置；
96.(3)10秒钟恒流放电脉冲(幅值为0.5c)；
97.(4)40秒钟静置；
98.(5)10秒钟恒流充电脉冲(幅值为0.5c)；
99.(6)将0.5c分别替换为1c、2c、3c、4c、5c后，重复(3)至(5)；
100.(7)静置24小时；
101.(8)调整温度至55℃，静置2h；
102.(9)将55℃分别替换为40℃、25℃、10℃、0℃、-10℃、-20℃，-30℃，重复(8)；
103.(10)调整温度为25℃；
104.(11)静置24小时。
105.这样，通过上述脉冲测试，可以根据充电部分每组脉冲中步骤(8)中的静置2h后的电池电压，从而获得不同电池soc和不同温度条件下，电池充电ocv(open circuit voltage，开路电压)随电池soc的变化曲线。类似的，可以根据放电部分每组脉冲中步骤(8)中的静置2h后的电池电压，从而获得不同电池soc和不同温度条件下，电池放电ocv随电池soc的变化曲线。其中，相同soc的条件下，电池充电ocv与放电ocv的平均值即为电池ocv，充电ocv与放电ocv的差异的1/2记为滞回电压，电池ocv随电池soc和温度的变化曲线即为电池ocv-soc曲线，滞回电压随电池soc和温度的变化曲线即为滞回电压-荷电状态曲线。
106.这样，在步骤s22中，可以根据所述初始属性信息分别建立不同阶数的多个初始等效电路模型。
107.示例地，根据所述初始属性信息建立的l阶初始等效电路模型如图3所示。其中，其中，u
ocv
和u
hys
分别代表电池开路电压和电池滞回电压，i和u分别代表电池电流(放电为正)和电池电压，r0为电池电路模型中的欧姆内阻，r1～r
l
为对应rc网络1～l的极化电阻，v
a1
～v
al
和v
b1
～v
bl
分别为模型参数，为储能元件的端电压。
108.又一示例，根据所述初始属性信息建立的l阶初始等效电路模型如图4所示。其中，u
ocv
和u
hys
分别代表电池开路电压和电池滞回电压，i和u分别代表电池电流(放电为正)和电池电压，r0为电池电路模型中的欧姆内阻，r1～r
l
为对应rc网络1～l的极化电阻，c1～c
l
为对
应rc网络1～l的极化电容，为储能元件的端电压。
109.其中，可以利用遗传算法、粒子群优化算法等优化算法获得电池模型参数。优化算法以参数初始值作为起点，根据交叉变异等优化规则，对模型参数进行不断的迭代更新，最终达到电池模型电压预测残差均方根最小的目标。
110.在具体实施时，针对每一种阶数的初始等效电路模型，可以通过多目标优化算法来确定该种阶数的初始等效电路模型中各元件参数的取值。以粒子群算法为例，可以通过随机初始化各元件参数的取值，并将该种阶数的等效电路模型在典型工况下的电压预测残差均方根作为适应值，通过不断地迭代优化，从而选择出该种阶数的等效电路模型的元件初始参数值(优化条件可以为该种阶数的等效电路模型在典型工况下的电压预测残差均方根小于预设阈值或是迭代次数达到阈值等等)，最终得到各种阶数的初始等效电路模型。此外，在一些实施例中，为了降低模型的复杂度，还可以对所述等效电路模型的阶数进行限制，例如可以将所述等效电路模型的阶数l限制为l≤3，其中l为正整数。
111.在步骤s23中，针对每一阶数的所述初始等效电路模型，分别测试该初始等效电路模型在目标工况下的计算误差信息和计算时间信息。
112.接上一实施例，在确定各种阶数的初始等效电路模型后，可以分别测试每一阶数的初始等效电路模型在典型工况下的计算误差信息和计算时间信息。
113.这样，在步骤s24中，可以根据每一所述初始等效电路模型的参数数量、所述计算误差信息以及所述计算时间信息计算每一所述初始等效电路模型的匹配度。并，在步骤s25中，将所述匹配度最优的初始等效电路模型的阶数确定为所述等效电路模型的阶数。
114.采用上述技术方案，通过对电池进行离线测试，从而获得了所述电池的初始属性信息。进一步的，根据所述初始属性信息建立并计算了不同阶数的初始等效电路模型的匹配度，从而能够提升所述初始等效电路模型模型的精度。
115.可选地，所述误差信息包括电压数据的采样误差因子和电流数据的采样误差因子，所述rls预估模型的辨识形式为：
[0116][0117]
其中，为所述rls预估模型第k时刻的输出信号的测量值，u
ocv
(k)和u
hys
(k)分别代表所述电池第k时刻的目标开路电压和目标滞回电压，为第k时刻所述bms采集的电池电压值，为所述rls预估模型第k时刻的输入信号，为第k时刻所述bms采集的电池电流值，θ(k)为所述rls预估模型第k时刻的参数矩阵，a1～a
l
，a0，b0～b
l
为所述参数矩阵中的参数，l为所述初始等效电路模型的阶数，t为所述bms的采样周期；
[0118]
其中，u(k)为第k时刻所述电池的电压真值，ε1为电压数据的采样
误差因子，i(k)为第k时刻所述电池的电流真值，ε2为电流数据的采样误差因子。
[0119]
以图3为例进行说明，根据图3所示出的l阶初始等效电路模型可以得到拉普拉斯空间中l阶初始等效电路模型的表达通式：
[0120][0121]
其中，u
s
为bms采集到的电池电压数据，i
s
为与电压同步采样的电池电流真值，u
ocv
为电池开路电压，u
hys
为电池滞回电压，r0代表电池电路模型中的欧姆内阻，r1～r
l
为对应rc网络1～l的极化电阻，c1～c
l
为对应rc网络1～l的极化电容。
[0122]
利用双线性变换并令y(z)＝u
ocv
(z)+u
hys
(z)-u(z)可得：
[0123]
其中，t为bms采样周期。
[0124][0125]
定义：
[0126][0127]
则式(2)可以进一步地简化为式(4)：
[0128][0129]
其中，b
i
和a
j
为简化系数(i＝0～l，j＝1～l
i
)，表达式由电池模型阶数l决定，如式(5)～(7)所示：
[0130]
当l＝1时，表达式如式(5)所示：
[0131][0132]
当l＝2时，表达式如式(6)所示：
[0133][0134]
当l＝3时，表达式如式(7)所示：
[0135][0136]
由式(4)可得电池模型的离散化表达形式，如式(8)所示：
[0137][0138]
进一步的，在一些实施例中，由于车载bms测量误差为有色噪声，因而还可以考虑以下两类误差：
[0139]
(1)、电压数据的测量误差：
[0140]
(2)、电流数据的测量误差：
[0141]
其中，和分别表示bms记录的电池电压测量值和电池电流测量值，ε1为电压数据的采样误差因子，ε2为电流数据的采样误差因子。
[0142]
根据式(9)和(10)，则车载环境下有：应当理解，当电池soc已知时，u
ocv
(k)和u
hys
(k)分别根据开路电压-荷电状态曲线和滞回电压-荷电状态曲线获得，因此也可作为已知量，进而结合式(8)可得到下式：
[0143][0144]
定义：则本公开中rls预估模型的辨识形式为：
[0145][0146]
则有：继而可以采用rls方法对参数θ(k)进行估计。如式(13)所示，其中，p(k)代表k时刻的rls协方差，λ为遗忘因子，取值范围可以是0.95～1，为输出y(k)的测量值。
[0147][0148]
之后利用多元方程组求解方法对式(3)、(5)～(7)进行反解，即可求得电池元件参数值p
parameter
＝[r0，r1～r
l
，c1～c
l
]。举例来讲，联立公式12、13可解得θ(k)，在θ(k)已知的情况下，可对式3、5～7(根据电池模型阶数确定)进行反解，从而得到所述电池元件参数值p
parameter
＝[r0，r1～r
l
，c1～c
l
]，进而可以通过代入所述电池模型参数确定电池的等效电路模型。
[0149]
可选地，对于l阶rc电路模型，如下式所示：
[0150][0151]
u(k)＝u
ocv
(soc(k),temp(k))+u
hys
(soc(k),temp(k))
[0152]
其中，t为采样间隔，c
m
为电池可用容量。
[0153]
可选地，对于l阶rv电路模型，如下式所示：
[0154][0155][0156]
其中，i>0，t为采样间隔，c
m
为电池可用容量。
[0157]
申请人发现，在一些场景中，所述误差信息包括电压数据的采样误差因子、电流数据的采样误差因子、所述电压数据和所述电流数据的采样时间差以及电池开路电压的误差中的至少一种。
[0158]
可选地，所述误差信息包括电压数据的采样误差因子、电流数据的采样误差因子、所述电压数据和所述电流数据的采样时间差以及电池开路电压的误差，所述rls预估模型的辨识形式为：
[0159][0160]
其中，为所述rls预估模型第k时刻的输出信号的测量值，u
ocv
(k)和u
hys
(k)分别代表所述电池第k时刻的目标开路电压和目标滞回电压，为第k时刻所述bms采集的电池电压值，为所述rls预估模型第k时刻的输入信号，为第k时刻所述bms采集的电池电流值，θ(k)为所述rls预估模型第k时刻的参数矩阵，a1～a
l
，a0，c1，d0～d
l
，c2为所述参数矩阵中的参数，l为所述初始等效电路模型的阶数，t为所述bms的采样周期；
[0161]
其中，u(k)为第k时刻所述电池的电压真值，ε1为电压数据的采样误差因子，i为与电压同步采样的电池电流，ε2为电流数据的采样误差因子，ε3为电流数据与电压数据之间的采样时间差值，为电流数据与电压数据之间的采样时间差值，和ε4分别表示电池soc存在误差时的电池ocv值以及对应的电池开路电压的误差。
[0162]
接上一实施例的式(10)进行说明，在一些场景中，还可以考虑如下误差：
[0163]
(3)、电流数据与电压数据之间的采样时间差，即：
[0164]
其中i'(k)＝i(k)+ε2，i为与电压同步采样的电池电流，ε3表示电流数据与电压数据之间的采样时间差值。参照式(10)，i(k)以及ε2的含义本公开在此不做赘述。
[0165]
这样，对式(14)作泰勒展开，可得式(15)：i'(k+ε3)＝i'(k)+ε3·
i'(k)(15)
[0166]
若则式(15)可化为：
[0167]
此外，由于等效电路模型模型的参数值估计过程与电池soc估计过程相互分离，因而在参数值估计过程中还可以考虑电池soc误差的影响。也就是说，当电池soc存在误差时，参数估计需要考虑电池ocv误差：
[0168]
即
[0169][0170]
其中，和ε4分别表示电池soc存在误差时的电池ocv值以及对应的电池开路电压的误差。
[0171]
值得注意的是，对于不同的bms系统来说，电压数据的采样误差因子、电流数据的采样误差因子以及所述电压数据和所述电流数据的采样时间差的数值大小也可能存在着差异。在所述电压数据和所述电流数据的采样时间差的数值大小不可忽略的情况下：
[0172]
考虑式(9)、(10)、(14)～式(17)，则车载环境下有：结合式(8)所示的电池模型的离散化表达形式可得到下式：
[0173][0174]
进一步的，定义进一步的，定义则在此情况下本公开的rls预估模型的辨识形式为：
[0175][0176]
采用上述技术方案，通过最小二乘法rls预估模型对电池模型参数进行反解，从而能够确定所述等效电路模型中的元件参数值。并且，所述最小二乘法rls预估模型还考虑了bms的采样器件所获取到的各类数据的误差，从而能够降低复杂车载条件下，bms的采样器件所获取到的各类数据的误差(例如bms的采样器件不断老化所导致的有色测量噪声，bms测量过程中电流和电压测量过程之间的非同步性所导致的等效电路模型参数值的估计误差等等)对等效电路模型参数值的估计影响，进而能够解决后续电池状态估计精度下降的问题。
[0177]
s13、根据预设修正边界值以及每类观测器的所述修正量，确定每类观测器的权重信息。
[0178]
在一种可能实现的方式中，所述根据预设修正边界值以及每类所述观测器的修正量，确定每类所述观测器的权重信息，包括：
[0179]
根据预设修正边界值以及每类所述观测器的修正量确定每类所述观测器的置信度；
[0180]
基于所述置信度计算得到每类所述观测器的权重信息。
[0181]
可选地，每类观测器计算过程中的预设修正边界值为同一个值，而每类观测器的修正量不同，这样，不同类型的观测器受不同因素的影响大小不同，确定的置信度可能不同。
[0182]
可选地，所述置信度γ
i
是通过如下公式计算得到的：
[0183][0184]
其中，l
i
(k)为第k时刻对应类型观测器的修正量，u
i
(k)为第k时刻所述电压数据，σ为标准差，η为预设修正边界值。
[0185]
基于不同类型观测器的置信度γ
i
，计算当前时刻每类观测器的权重信息f
i
，计算过程的辨析行驶可以如下所示：
[0186][0187]
这样，根据车辆的每一时刻实车状态，可以确定不同类型的观测器的修正量，进而计算对应的置信度，确定每类观测器的权重信息。
[0188]
s14、基于所述权重信息以及所述每类观测器的soc估计值，确定所述电池的soc目标值。
[0189]
在步骤s14中，所述观测器是通过如下方式确定所述电池的soc值的：
[0190]
s1401、基于前一刻的电池状态向量，根据所述元件参数值、所述电流数据以及所述温度数据，利用电池空间方程的输出方程，确定每类所述观测器的状态向量先验值。
[0191]
具体实施时，若车辆处于从下电工况进入上电的工况，则根据各项初始值条件，对每类观测器进行初始化。这种情况下，前一时刻的电池状态向量，根据所述元件参数值、所述电流数据以及所述温度数据即为初始化后的电池状态向量，元件参数值、电流数据以及温度数据。
[0192]
示例地，电池空间方程的输出方程可以是如下辨识形式：
[0193][0194]
其中，为状态向量先验值，x
state
(k-1)为第k-1时刻的状态向量，p
parameter
为所述元件参数值，i(k)为第k时刻的所述电流数据。
[0195]
在一种可能实现的方式中，所述每类所述观测器的修正量是通过以下方式得到的：
[0196]
根据所述状态向量先验值、所述元件参数值、所述电流数据以及所述温度数据确定每类所述观测器的修正量。
[0197]
示例地，通过以下辨识形式确定每类所述观测器的修正量l
i
(k)：
[0198][0199]
其中，i(k)和temp(k)分别为第k时刻的电流数据和温度数据，u(k)为第k时刻的电压数据，p
parameter
为所述元件参数值，为所述状态向量先验值。
[0200]
s1402、基于每类所述观测器的增益矩阵以及所述状态向量先验值确定状态向量后验值，其中，不同类型的观测器的增益矩阵不同。
[0201]
示例地，当所述观测器类型为luenberger观测器时，所述状态向量后验值可以通过如下辨识形式确定：
[0202][0203]
其中，kp
joint
和ki
joint
均为luenberger观测器的预设增益。
[0204]
又一示例，当所述观测器类型为ekf(extended kalman filter扩展卡尔曼滤波)、ukf(unscented kalman filter无迹卡尔曼滤波)或者aukf(adaptive unscented kalman filter自适应无迹卡尔曼滤波)观测器时，所述状态向量后验值可以通过如下辨识形式确定：
[0205][0206]
其中，k
joint
为ekf、ukf、aukf观测器的预设增益。
[0207]
进一步地，当所述观测器为ekf时，其预设增益k
joint
可以通过如下辨识形式确定：k
joint
＝(a
·
p
·
a
t
+q)
·
c
t
·
[c
·
(a
·
p
·
a
t
+q)
·
c
t
+r]
[0208]
其中，p，q，r分别为电池状态协方差初始值、过程噪声方差初始值以及测量噪声方差初始值，
[0209]
当所述观测器为ukf、aukf时，其预设增益k
joint
可以通过如下辨识形式确定：
[0210][0211]
其中，为第j个特征点在计算状态向量协方差时的权重系数，为第j个特征点在计算期望时的权重系数，为第j个特征点输出值与输出的期望之差，r为测量噪声方差初始值。
[0212]
s1403、根据所述状态向量后验值以及电池状态空间方程确定所述电池的soc值。
[0213]
可选地，所述电池状态空间方程为：
[0214][0215]
其中，p
parameter
＝[r0，r1～r
l
，c1～c
l
]，为所述元件参数值的列向量，r0为所述等效电路模型的欧姆内阻，r1～r
l
为所述等效电路模型的极化内阻，c1～c
l
为所述等效电路模型的极化电容，i(k)为第k时刻de所述电流数据，temp(k)为第k时刻所述温度数据，ω(k)为所述过程噪声，γ(k)为所述测量噪声，其方差分别为过程噪声方差和测量噪声方差，f(
·
)和g(
·
)均为非线性函数，soc(k)为所述电池的soc值。
[0216]
最终，在步骤s14中，基于所述电池的soc值，根据如下辨识形式确定所述电池的soc值soc
ec
：
[0217][0218]
其中，soc
i
为不同类型的观测器的电池的soc值。
[0219]
本公开还提供一种确定电池荷电状态的装置，参考图5所示出的一种确定电池荷电状态的装置的框图，所述装置500包括：获取模块510，第一确定模块520，第二确定模块530，第四确定模块540。
[0220]
获取模块510，用于获取电池的状态数据，所述状态数据包括电流数据、电压数据、温度数据；
[0221]
第一确定模块520，用于根据所述电池的等效电路模型的元件参数值、所述状态数据、电池特性数据以及观测器初始参数，确定每类观测器的soc估计值和修正量；
[0222]
第二确定模块530，用于根据预设修正边界值以及每类观测器的所述修正量，确定每类观测器的权重信息；
[0223]
第三确定模块540，用于基于所述权重信息以及所述每类观测器的soc估计值，确定所述电池的soc目标值。
[0224]
上述装置通过最小二乘法rls预估模型确定所述等效电路模型中的元件参数值。并且，在确定所述电池的等效电路模型中的元件参数值时，还考虑了至少包括电压数据的采样误差因子和/或电流数据的采样误差因子的误差信息，从而能够降低采样误差所带来的影响，进而能够提升确定的所述等效电路模型中的元件参数值的准确度，最终起到提升所确定的等效电路模型的精度的效果。进一步地，基于不同类型的观测器确定的电池soc值，确定电池soc目标值，进而提升确定电池荷电状态的合理性，保证车辆的高效管理和可靠运行。
[0225]
可选地，所述第二确定模块包括：第一确定子模块，用于基于所述电池的离线测试确定等效电路模型，以及所述等效电路模型的元件参数值与荷电状态和电池温度的函数关系式；
[0226]
第二确定子模块，用于根据所述温度数据和所述电池前一时刻的荷电状态数据通过所述函数关系式确定所述元件参数值，其中，所述等效电路模型为rc电路模型或者rv电路模型。
[0227]
可选地，所述第二确定模块包括：第九确定子模块，用于基于所述电池的等效电路模型、误差信息、电池特性数据以及所述状态数据，通过最小二乘法rls预估模型，确定等效
电路模型中的元件参数值；
[0228]
其中，所述等效电路模型是基于对所述电池进行离线测试得到的，所述电池特性数据包括不同温度下的开路电压-荷电状态曲线以及滞回电压-荷电状态曲线，其中，所述等效电路模型为rc电路模型。
[0229]
可选地，所述装置还包括第四确定模块，用于确定所述等效电路模型的阶数，所述第四确定模块包括：
[0230]
获取子模块，用于基于电池离线测试，获得所述电池的初始属性信息，所述初始属性信息包括所述电池的开路电压-荷电状态曲线、滞回电压-荷电状态曲线；
[0231]
创建子模块，用于根据所述初始属性信息分别建立不同阶数的多个初始等效电路模型；
[0232]
测试子模块，用于针对每一阶数的所述初始等效电路模型，分别测试该初始等效电路模型在目标工况下的计算误差信息和计算时间信息；
[0233]
计算子模块，用于根据每一所述初始等效电路模型的参数数量、所述计算误差信息以及所述计算时间信息计算每一所述初始等效电路模型的匹配度；
[0234]
确定子模块，用于将所述匹配度最优的初始等效电路模型确定为所述电池的等效电路模型。
[0235]
可选地，所述误差信息包括电压数据的采样误差因子、电流数据的采样误差因子、所述电压数据和所述电流数据的采样时间差以及电池开路电压的误差中的至少一种。
[0236]
可选地，所述误差信息包括电压数据的采样误差因子、电流数据的采样误差因子、所述电压数据和所述电流数据的采样时间差以及电池开路电压的误差，所述rls预估模型的辨识形式为：
[0237][0238]
其中，为所述rls预估模型第k时刻的输出信号的测量值，u
ocv
(k)和u
hys
(k)分别代表所述电池第k时刻的目标开路电压和目标滞回电压，为第k时刻所述bms采集的电池电压值，为所述rls预估模型第k时刻的输入信号，为第k时刻所述bms采集的电池电流值，θ(k)为所述rls预估模型第k时刻的参数矩阵，a1～a
l
，a0，c1，d0～d
l
，c2为所述参数矩阵中的参数，l为所述初始等效电路模型的阶数，t为所述bms的采样周期；
[0239]
其中，u(k)为第k时刻所述电池的电压真值，ε1为电压数据的采样误差因子，i为与电压同步采样的电池电流，ε2为电流数据的采样误差因子，ε3为电流数据与电压数据之间的采样时间差值，为电流数据与电压数据之间的采样时间差值，和ε4分别表示电池soc存在误差时的电池ocv值以及对应的电池开路电压的误差。
[0240]
可选地，所述第四确定模块，包括：
[0241]
第三确定子模块，用于基于前一刻的电池状态向量，根据所述元件参数值、所述电流数据以及所述温度数据，利用电池空间方程的输出方程，确定每类所述观测器的状态向量先验值；
[0242]
第四确定子模块，用于基于每类所述观测器的增益矩阵以及所述状态向量先验值确定状态向量后验值，其中，不同类型的观测器的增益矩阵不同；
[0243]
第五确定子模块，用于根据所述状态向量后验值以及电池状态空间方程确定所述电池的soc值。
[0244]
可选地，所述第四确定模块，包括：
[0245]
第六确定子模块，用于根据所述状态向量先验值、所述元件参数值、所述电流数据以及所述温度数据确定每类所述观测器的修正量。
[0246]
可选地，所述第二确定模块，包括：
[0247]
第七确定子模块，用于根据预设修正边界值以及每类所述观测器的修正量确定每类所述观测器的置信度；
[0248]
第八确定子模块，用于基于所述置信度计算得到每类所述观测器的权重信息。
[0249]
可选地，所述置信度γ
i
是通过如下公式计算得到的：
[0250][0251]
其中，l
i
(k)为第k时刻对应类型观测器的修正量，u
i
(k)为第k时刻所述电压数据，σ为标准差，η为预设修正边界值。
[0252]
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0253]
此外值得说明的是，为描述的方便和简洁，说明书中所描述的实施例均属于优选实施例，其所涉及的部分并不一定是本发明所必须的，例如，第一确定模块和第二确定模块，在具体实施时可以是相互独立的装置也可以是同一个装置，本公开对此不作限定。
[0254]
本公开还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述任一实施例中所述的确定电池荷电状态的方法的步骤。
[0255]
本公开还提供一种电子设备，包括：存储器，其上存储有计算机程序；
[0256]
处理器，用于执行所述存储器中的所述计算机程序，以实现上述任一实施例中所述的确定电池荷电状态的方法的步骤。
[0257]
图6是根据一示例性实施例示出的一种用于确定电池荷电状态的装置600的框图。参照图6，装置600可以包括以下一个或多个组件：处理组件602，存储器604，电力组件606，多媒体组件606，输入/输出(i/o)的接口612，传感器组件614，以及通信组件616。
[0258]
处理组件602通常控制装置600的整体操作，诸如数据的获取、传感器数据的处理，rls算法的求解等等。处理组件602可以包括一个或多个处理器620来执行指令，以完成上述
确定电池荷电状态的方法的全部或部分步骤。此外，处理组件602可以包括一个或多个模块，便于处理组件602和其他组件之间的交互。例如，处理组件602可以包括多媒体模块，以方便多媒体组件606和处理组件602之间的交互。
[0259]
存储器604被配置为存储各种类型的数据以支持在装置600的操作。这些数据的示例包括用于在装置600上操作的任何应用程序或方法的指令，历史电流数据、历史电压数据、电池的开路电压-荷电状态曲线、滞回电压-荷电状态曲线等。存储器604可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(sram)，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，可编程只读存储器(prom)，只读存储器(rom)，磁存储器，快闪存储器，磁盘或光盘。
[0260]
电力组件606为装置600的各种组件提供电力。电力组件606可以包括电源管理系统，一个或多个电源，及其他与为装置600生成、管理和分配电力相关联的组件。
[0261]
多媒体组件606包括在所述装置600和用户之间的提供一个输出接口的屏幕。在一些实施例中，屏幕可以包括液晶显示器(lcd)和触摸面板(tp)。如果屏幕包括触摸面板，屏幕可以被实现为触摸屏，以接收来自用户的输入信号。触摸面板包括一个或多个触摸传感器以感测触摸、滑动和触摸面板上的手势。所述触摸传感器可以不仅感测触摸或滑动动作的边界，而且还检测与所述触摸或滑动操作相关的持续时间和压力。
[0262]
i/o接口612为处理组件602和外围接口模块之间提供接口，上述外围接口模块可以是点击轮，按钮等。
[0263]
传感器组件614包括一个或多个传感器，用于为装置600提供各个方面的状态评估。例如，传感器组件614可以检测电池的温度、电流等等。在一些实施例中，该传感器组件614例如可以包括加温度传感器，速度传感器，电流传感器等等。
[0264]
通信组件616被配置为便于装置600和其他设备之间有线或无线方式的通信。装置600可以接入基于通信标准的无线网络，如wifi，2g或3g，或它们的组合。
[0265]
在示例性实施例中，装置600可以被一个或多个应用专用集成电路(asic)、数字信号处理器(dsp)、数字信号处理设备(dspd)、可编程逻辑器件(pld)、现场可编程门阵列(fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述确定电池荷电状态的方法。
[0266]
在示例性实施例中，还提供了一种包括指令的非临时性计算机可读存储介质，例如包括指令的存储器604，上述指令可由装置600的处理器620执行以完成上述确定电池荷电状态的方法。例如，所述非临时性计算机可读存储介质可以是rom、随机存取存储器(ram)、cd-rom、磁带、软盘和光数据存储设备等。
[0267]
在另一示例性实施例中，还提供一种计算机程序产品，该计算机程序产品包含能够由可编程的装置执行的计算机程序，该计算机程序具有当由该可编程的装置执行时用于执行上述的确定电池荷电状态的方法的代码部分。
[0268]
本公开还提供一种电池管理系统，包括上述任一项确定电池荷电状态的装置。
[0269]
关于上述实施例中的电池管理系统，其中装置执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0270]
以上结合附图详细描述了本公开的优选实施方式，但是，本公开并不限于上述实施方式中的具体细节，在本公开的技术构思范围内，可以对本公开的技术方案进行多种简
单变型，这些简单变型均属于本公开的保护范围。
[0271]
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本公开对各种可能的组合方式不再另行说明。
[0272]
此外，本公开的各种不同的实施方式之间也可以进行任意组合，只要其不违背本公开的思想，其同样应当视为本公开所公开的内容。