一种复合锂离子电池健康状况估计方法

1.本发明属于复合锂离子电池健康状况估计领域，尤其涉及一种复合锂离子电池健康状况估计方法。背景技术：2.电池是一个复杂的电化学系统，存在各种老化机制导致容量衰减和功率衰退，因此，基于可用容量的锂离子电池健康状况(soh)估计方法不能在电动车上直接测量，车载soh估计算法给出电池老化机制直接相关的替代参数。根据锂离子电池不同的工作条件，单体电池内部产生不同的老化机制，对于基于石墨阳极的锂离子电池，老化机制包括但不限于固体电解质相(sei)形成、锂电镀、颗粒开裂、活性物质(am)溶解和am隔离等。常用的电池soh估计方法是利用电化学阻抗谱(eis)建立等效电路模型(ecm)参数，另一种方法称为数据驱动方法，该方法缺乏对模型参数的物理解释，需要采集大量数据以涵盖所有可能的电池退化机制。目前，有学者提出增量容量分析法(ica)，该方法是一种强大的电池老化检测算法。weng等人利用ica分析部分充电数据，dubarry等人将容量分析结果编制成n维查找表，将其输入电池管理系统中(bms)，用于车载soh估计。n维查找表能够利用部分重要特征来识别电池的不同老化机制，例如可循环锂损失、负极am损失和正极活性材料损失。改进的soh估计方法是基于恒流电流计算经验ecm参数或电化学模型参数。bartlett等人提出一种基于降阶增强型单粒子模型(espm)参数的双soc/soh估计器，该方法利用扩展卡尔曼滤波(ekf)估计高电压门限值的负极锂浓度，确定电池容量。但是，该模型存在不可观测性问题，并且模型收敛性仅适用于低电压的工作条件，不适合车载soh估计。技术实现要素：3.本发明的目的是克服现有技术中的不足，提供一种复合锂离子电池健康状况估计方法；4.这种复合锂离子电池健康状况估计方法，包括以下步骤：5.步骤1、基于增强型单粒子模型(espm)，建立复合锂离子电池电化学模型，并给出动态质量和控制方程；6.步骤2、在循环老化实验的实验周期中，基于nelder-mead算法计算放电容量，并利用循环老化实验中每个微循环(micro-cycle)充电周期对复合锂离子电池电化学模型进行校准；7.步骤3、建立安时吞吐量ah与锂离子电池放电容量q的关系方程，建立安时吞吐量ah与可循环锂离子损耗摩尔数nli的关系方程；8.步骤4、建立复合锂离子电池soh估计方法；9.步骤5、设计循环老化实验条件，根据步骤1至步骤4得出每次老化实验周期中，每个单体锂离子电池循环充电吞吐量。10.作为优选，步骤5中循环老化实验条件包括实验温度、充电水平、老化实验周期和充电吞吐量。11.作为优选，步骤1中控制方程具体为：12.正负电极活性材料的质量守恒方程为：[0013][0014]上式中，cam为活性物质的摩尔浓度，单位为mol·m-3；t为采样时间；ds,am表示活性固体物质的扩散系数；r为离子半径变量；jam(t)为插层电流密度，f为法拉第常数；rp,am表示活性物质粒子直径，单位为m；上式(1)表示锂离子在固体颗粒中的扩散过程；[0015]用电解质的质量守恒方程表示锂离子在电解液中的扩散过程，由菲克第二定律控制：[0016][0017]上式中，εe为电解液粒子体积分数；ce为电解质的摩尔浓度；de为电解质扩散系数，单位为m2s-1；x为负极化学计量值；aam为活性物质的常数系数；为正粒子转移数；[0018]电解质电极电位计算公式为：[0019][0020]上式中，为电解质电位，ie(t)为t时刻电解质产生的电流值，σe表示电解质电导率，单位为ω·m-1；r为热力学常数，t为热力学温度，单位为k；为负粒子转移数；γ为活性物质中的电流分数；电解质电极电位由电解质浓度、电解质电导率和活性化学物质确定；[0021]活性物质的过电位由开路电压uam减去活性物质电极电位以及接触电阻分压rc,amβami(t)得到：[0022][0023]上式中，ηam(t)为t时刻活性物质的活化过电位，单位为v；uam为活性物质的开路电压值，单位为v；为活性物质的电位；rc,am为活性物质的接触电阻，单位为ω；βam为活性物质的电流分数；i(t)为t时刻的电流值，单位为a；[0024]通过总电流除以电活性表面积计算插层电流密度，并决定公式(1)和公式(2)的最终解：[0025][0026]上式中，i0,am表示活性物质在充电状态下产生的电流值；a为电极横截面积；l为电极厚度；εam为活性物质体积分数；[0027]活性物质在充电状态下产生的电流值由活性物质的摩尔浓度和活性物质粒子半径确定[0028][0029]上式中，kam为活性物质的常数系数；cmax,am为活性物质最大摩尔浓度；rp,am为正电极的接触电阻；[0030]电解液中欧姆过电位的欧姆定律为：[0031][0032]上式中，ve为电解质的端电压值，单位为v；[0033]由正负电极过电位差值、欧姆电阻过电位值和电池内阻产生的电压值确定单体电池端电压：[0034][0035]上式中，为t时刻的正极电极电位，为t时刻的负极电极电位，ve(t)为t时刻的电解质的端电压值；rcell为电池内阻，单位为ω。[0036]作为优选，步骤2具体为：[0037]步骤2-1、建立循环锂离子损失摩尔数nli的计算方程：[0038]nli＝nli,n+nli,pꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(9)[0039]上式中，nli,n和nli,p分别表示负极和正极中的可循环锂离子摩尔数，计算公式如(10)和(11)：[0040]nli,n＝x·cmax,n·a·ln·εam,nꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(10)[0041]nli,p＝y·cmax,p·a·lp·εam,pꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(11)[0042]上式中，a表示当前集电极面积；x表示正极的荷电状态，y表示负极的荷电状态；cmax,n表示锂离子在负极中的最大饱和浓度，cmax,p表示锂离子在正极中的最大饱和浓度，单位均为mol·m-3，εam,n表示负极中活性物质体积分数，εam,p表示正极中活性物质体积分数；l为电极厚度，单位为m，ln为负极厚度，lp为正极厚度；[0043]步骤2-2、建立锂离子电池放电容量的计算公式：[0044][0045]上式中，y0表示充电正极荷电状态，yf表示放电正极荷电状态，x0表示充电负极荷电状态，xf表示放电负极荷电状态；根据公式(8)，电极电位由下式得出：[0046][0047]上式中，cs,i(t)为t时刻的正负极固体锂离子浓度，其中负极固体锂离子浓度为cs,n；正极固体锂离子浓度为cs,p；ηi为活化过电位，ηn为负极活化过电位，ηp为正极活化过电位；可循环锂离子损失摩尔数nli与cs,p(t)/cs,n(t)之间的关系方程为：[0048][0049][0050]步骤2-3、利用循环老化实验中每个微循环(micro-cycle)充电周期对复合锂离子电池电化学模型进行校准。[0051]作为优选，步骤3具体包括如下步骤：[0052]步骤3-1、建立m个单体电池放电容量值和安时吞吐量ah之间的函数关系式：[0053]ql＝f1,lah0.5+f2,l l＝1,2,...,mꢀꢀꢀꢀꢀ(16)[0054]上式中，f1,l和f2,l为校准参数；[0055]步骤3-2、建立安时吞吐量ah与可循环锂离子损耗摩尔数关公式：[0056]nli,l＝gl(ah)and rincrease＝hl(ah)ꢀꢀꢀꢀꢀꢀ(17)[0057]上式中，nli,l为第l个电池的可循环锂离子损耗摩尔数，rincrease为内阻增加值；gl(ah)和hl(ah)的表达式为：[0058]gl(ah)＝k1,lah0.5+k2,l l＝1,2,...,mꢀꢀꢀꢀꢀꢀꢀꢀ(18)[0059]hl(ah)＝k3,lah+k4,l l＝1,2,...,mꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(19)[0060]上式中，k1，k2，k3，k4均为方程系数；[0061]步骤3-3、基于仿真电压vsim和实验电压vexp的均方根误差(rmse)建立成本函数，基于espm模型对复合锂离子电池电化学模型进行近似降阶处理，成本函数为：[0062][0063]上式中，n为电压采样个数。[0064]作为优选，步骤4中soh的定义为：[0065][0066]上式中，下标eol表示终止寿命，qactual表示单体电池实际容量值；qnew表示全新电池容量值；qeol表示达到终止寿命时单体电池容量值；根据函数gl(ah)和hl(ah)确定安时吞吐量ah与参数nli、rincrease之间的关系，根据步骤1中的公式(1)至公式(8)，建立soh映射公式：soh＝h(nli,rincrease)，基于nli、rincrease参数对估计当前soh值。[0067]作为优选，步骤4具体包括如下步骤：[0068]步骤4-1、假设tsample表示与终止寿命(eol)相对应的采样时间，下标k为当前采样时间，建立的soh估计模型为：[0069]sohk＝μk·ah0.5+αkꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(22)[0070]上式中，μk和αk为方程系数；[0071]步骤4-2、建立soh估计方法状态方程和量测方程；[0072]方法状态方程为：[0073]xk＝[μk；αk]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(23)[0074][0075]上式中，xk为状态空间矩阵，包含方程系数向量μk和αk；χ为过程噪声，δ为高斯白噪声；[0076]量测方程为，[0077]zk＝sohk＝μk·ah0.5+αk+νkꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(25)[0078]上式中，zk为量测方程，νk为量测误差；利用单个粒子的后验分布p(xk|z1:k)预测锂离子电池剩余使用寿命；[0079]步骤4-3、建立状态估计的密度函数：[0080][0081]上式中，ns为粒子数量，ω为每个粒子的权重系数；[0082]步骤4-4、计算未来每个采样时间j＝1:tsample-k的soh预测值如公式(27)：[0083][0084]由公式(28)得出[0085][0086]本发明的有益效果是：[0087]本发明属于一种基于增强型单粒子模型(espm)的参数估计算法，以及基于循环老化实验数据和电化学模型的复合锂离子电池健康状况估计方法；本发明引入增强型单粒子模型(espm)，提取与单体电池soh直接相关的参数特征，并利用实时充电工况(包括快速充电)获取模型参数特征；本发明选取lmo-nmc-石墨复合锂离子电池进行循环老化实验，验证新型电池模型和soh估计方法的可行性和有效性。[0088]本发明利用车载锂离子电池循环充放电实验数据，建立基于电化学模型的lmo-nmc-石墨复合锂离子电池soh估计方法；引入增强单粒子模型(espm)，基于老化实验周期中部分充电过程数据辨识可循环锂离子损失摩尔数nli。附图说明[0089]图1-1为基于老化实验的nli变化曲线图，图1-2为基于老化实验的r变化曲线图；[0090]图2-1为基于可循环锂离子损失的正极电化学计量变化图，图2-2为基于可循环锂离子损失的负极电化学计量变化图；[0091]图3-1为实施例中电池1的老化实验周期中soh与ah吞吐量关系曲线图，图3-2为实施例中电池2的老化实验周期中soh与ah吞吐量关系曲线图；图3-3为实施例中电池3的老化实验周期中soh与ah吞吐量关系曲线图；图3-4为实施例中电池4的老化实验周期中soh与ah吞吐量关系曲线图；图3-5为实施例中电池5的老化实验周期中soh与ah吞吐量关系曲线图。具体实施方式[0092]下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干修饰，这些改进和修饰也落入本发明权利要求的保护范围内。[0093]实施例一[0094]本技术实施例一提供了一种复合锂离子电池健康状况估计方法：[0095]步骤1、基于增强型单粒子模型(espm)，建立复合锂离子电池电化学模型，并给出动态质量和控制方程；[0096]步骤2、在循环老化实验的实验周期中，基于nelder-mead算法计算放电容量，并利用循环老化实验中每个微循环(micro-cycle)充电周期对复合锂离子电池电化学模型进行校准；[0097]步骤3、建立安时吞吐量ah与锂离子电池放电容量q的关系方程，建立安时吞吐量ah与可循环锂离子损耗摩尔数nli的关系方程；[0098]步骤4、建立复合锂离子电池soh估计方法；[0099]步骤5、设计循环老化实验条件，根据步骤1至步骤4得出每次老化实验周期中，每个单体锂离子电池循环充电吞吐量。[0100]实施例二[0101]在实施例一的基础上，本技术实施例二提供了实施例一中估计方法在现实中的应用：[0102]步骤1、基于增强型单粒子模型(espm)，建立复合锂离子电池电化学模型，并给出动态质量和控制方程；控制方程具体为：[0103]正负电极活性材料的质量守恒方程为：[0104][0105]上式中，cam为活性物质的摩尔浓度，单位为mol·m-3；t为采样时间；ds,am表示活性固体物质的扩散系数；r为离子半径变量；jam(t)为插层电流密度，f为法拉第常数；rp,am表示活性物质粒子直径，单位为m；上式(1)表示锂离子在固体颗粒中的扩散过程；[0106]用电解质的质量守恒方程表示锂离子在电解液中的扩散过程，由菲克第二定律控制：[0107][0108]上式中，εe为电解液粒子体积分数；ce为电解质的摩尔浓度；de为电解质扩散系数，单位为m2s-1；x为负极化学计量值；aam为活性物质的常数系数；为正粒子转移数；[0109]电解质电极电位计算公式为：[0110][0111]上式中，为电解质电位，ie(t)为t时刻电解质产生的电流值，σe表示电解质电导率，单位为ω·m-1；r为热力学常数，t为热力学温度，单位为k；为负粒子转移数；γ为活性物质中的电流分数；电解质电极电位由电解质浓度、电解质电导率和活性化学物质确定；[0112]活性物质的过电位由开路电压uam减去活性物质电极电位以及接触电阻分压rc,amβami(t)得到：[0113][0114]上式中，ηam(t)为t时刻活性物质的活化过电位，单位为v；uam为活性物质的开路电压值，单位为v；为活性物质的电位；rc,am为活性物质的接触电阻，单位为ω；βam为活性物质的电流分数；i(t)为t时刻的电流值，单位为a；[0115]通过总电流除以电活性表面积计算插层电流密度，并决定公式(1)和公式(2)的最终解：[0116][0117]上式中，i0,am表示活性物质在充电状态下产生的电流值；a为电极横截面积；l为电极厚度；εam为活性物质体积分数；[0118]活性物质在充电状态下产生的电流值由活性物质的摩尔浓度和活性物质粒子半径确定[0119][0120]上式中，kam为活性物质的常数系数；cmax,am为活性物质最大摩尔浓度；rp,am为正电极的接触电阻；[0121]电解液中欧姆过电位的欧姆定律为：[0122][0123]上式中，ve为电解质的端电压值，单位为v；[0124]由正负电极过电位差值、欧姆电阻过电位值和电池内阻产生的电压值确定单体电池端电压：[0125][0126]上式中，为t时刻的正极电极电位，为t时刻的负极电极电位，ve(t)为t时刻的电解质的端电压值；rcell为电池内阻，单位为ω。[0127]上式(1)至(8)中，参数下标的定义如表1所示，电池动态质量和控制方程如下表2所示；[0128]表1参数下标定义[0129]下标表示下标定义0充电状态e电解质exp实验数据sim仿真数据cell单体电池am活性物质p正电极n负电极s固体物质[0130]表2电池动态质量和控制方程[0131][0132]步骤2、在循环老化实验的实验周期中，基于nelder-mead算法计算放电容量，并利用循环老化实验中每个微循环(micro-cycle)充电周期对复合锂离子电池电化学模型进行校准；[0133]步骤2-1、建立循环锂离子损失摩尔数nli的计算方程：[0134]nli＝nli,n+nli,pꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(9)[0135]上式中，nli,n和nli,p分别表示负极和正极中的可循环锂离子摩尔数，计算公式如(10)和(11)：[0136]nli,n＝x·cmax,n·a·ln·εam,nꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(10)[0137]nli,p＝y·cmax,p·a·lp·εam,pꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(11)[0138]上式中，a表示当前集电极面积；x表示正极的荷电状态，y表示负极的荷电状态；cmax,n表示锂离子在负极中的最大饱和浓度，cmax,p表示锂离子在正极中的最大饱和浓度，单位均为mol·m-3，εam,n表示负极中活性物质体积分数，εam,p表示正极中活性物质体积分数；l为电极厚度，单位为m，ln为负极厚度，lp为正极厚度；[0139]步骤2-2、建立锂离子电池放电容量的计算公式：[0140][0141]上式中，y0表示充电正极荷电状态，yf表示放电正极荷电状态，x0表示充电负极荷电状态，xf表示放电负极荷电状态；根据公式(8)，电极电位由下式得出：[0142][0143]上式中，cs,i(t)为t时刻的正负极固体锂离子浓度，其中负极固体锂离子浓度为cs,n；正极固体锂离子浓度为cs,p；ηi为活化过电位，ηn为负极活化过电位，ηp为正极活化过电位，对应的正极为liymn2o4，负极为lixc6；可循环锂离子损失摩尔数nli与cs,p(t)/cs,n(t)之间的关系方程为：[0144][0145][0146]步骤2-3、利用循环老化实验中每个微循环(micro-cycle)充电周期对复合锂离子电池电化学模型进行校准；[0147]步骤3、建立安时吞吐量ah与锂离子电池放电容量q的关系方程，建立安时吞吐量ah与可循环锂离子损耗摩尔数nli的关系方程；[0148]步骤3-1、建立m个单体电池放电容量值和安时吞吐量ah之间的函数关系式：[0149]ql＝f1,lah0.5+f2,l l＝1,2,...,mꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(16)[0150]上式中，f1,l和f2,l为校准参数；[0151]步骤3-2、建立安时吞吐量ah与可循环锂离子损耗摩尔数关公式：[0152]nli,l＝gl(ah)and rincrease＝hl(ah)ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(17)[0153]上式中，nli,l为第l个电池的可循环锂离子损耗摩尔数，rincrease为内阻增加值；gl(ah)和hl(ah)的表达式为：[0154]gl(ah)＝k1,lah0.5+k2,l l＝1,2,...,mꢀꢀꢀꢀꢀꢀꢀꢀ(18)[0155]hl(ah)＝k3,lah+k4,l l＝1,2,...,mꢀꢀꢀꢀꢀꢀꢀꢀꢀ(19)[0156]上式中，k1，k2，k3，k4均为方程系数；[0157]步骤3-3、基于仿真电压vsim和实验电压vexp的均方根误差(rmse)建立成本函数，基于espm模型对复合锂离子电池电化学模型进行近似降阶处理，成本函数为：[0158][0159]上式中，n为电压采样个数；[0160]步骤4、建立复合锂离子电池soh估计方法；soh的定义为：[0161][0162]上式中，下标eol表示终止寿命，qactual表示单体电池实际容量值；qnew表示全新电池容量值；qeol表示达到终止寿命时单体电池容量值；根据函数gl(ah)和hl(ah)确定安时吞吐量ah与参数nli、rincrease之间的关系，根据步骤1中的公式(1)至公式(8)，建立soh映射公式：soh＝h(nli,rincrease)，基于nli、rincrease参数对估计当前soh值；[0163]步骤4-1、假设tsample表示与终止寿命(eol)相对应的采样时间，下标k为当前采样时间，建立的soh估计模型为：[0164]sohk＝μk·ah0.5+αkꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(22)[0165]上式中，μk和αk为方程系数；[0166]步骤4-2、建立soh估计方法状态方程和量测方程；[0167]方法状态方程为：[0168]xk＝[μk；αk]ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(23)[0169][0170]上式中，xk为状态空间矩阵，包含方程系数向量μk和αk；χ为过程噪声，δ为高斯白噪声；[0171]量测方程为，[0172]zk＝sohk＝μk·ah0.5+αk+νkꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ(25)[0173]上式中，zk为量测方程，νk为量测误差；利用单个粒子的后验分布p(xk|z1:k)预测锂离子电池剩余使用寿命；[0174]步骤4-3、建立状态估计的密度函数：[0175][0176]上式中，ns为粒子数量，ω为每个粒子的权重系数；[0177]步骤4-4、计算未来每个采样时间j＝1:tsample-k的soh预测值如公式(27)：[0178][0179]由公式(28)得出[0180][0181]步骤5、设计循环老化实验条件，根据步骤1至步骤4得出每次老化实验周期中，每个单体锂离子电池循环充电吞吐量。本实施例设计循环老化实验条件如下表3所示，实验温度为30℃。“2级”充电方案表示恒流恒压充电方案(cccv)，“3级”充电方案表示直流快速充电方案(dc fast charging)，老化实验周期约为三个月，每隔一周进行一次容量测试。表3设置soc最小值(socmin)，例如socmin＝45％、socmin＝35％和socmin＝25％，当soc值达到设置的socmin后，使用恒流恒压充电方案(cccv)对锂离子电池充电。恒流充电阶段的充电倍率为3c/2和5c。充电倍率分为3c/2和5c，表3给出每次老化实验周期中，每个单体锂离子电池循环充电吞吐量(ah)。[0182]表3循环老化实验条件表[0183]老化实验socmin(％)充电水平充电倍率温度(℃)充电吞吐量(ah)1452级3c/23024,3612352级3c/23025,3533252级3c/23027,5974453级5c3024,2925353级5c3020,141[0184]实验结果1：[0185]图1-1至图1-2为基于老化实验充电数据的可循环锂离子摩尔数nli和内阻增加值rincrease的变化曲线：(1)nli初始值为0.75mol，容量衰减率为10％之后，nli减少到0.68mol；[0186](2)容量衰减率为0-6％之间，nli减少量较大(0.05mol)，容量衰减率为6％-10％之间，nli较少量为一半；[0187](3)rincrease初始值约为0.65mω，容量衰减率为10％之后，rincrease增加到1.4mω；[0188](4)容量衰减率为0-6％之间，rincrease增加量较大(0.65mω)，容量衰减率为6％-10％之间，rincrease增加量较小(0.15mω)；[0189](5)两条变化曲线均符合线性关系。[0190]实验结果2：[0191]基于可循环锂离子损失的电极计量变化曲线如图2-1至图2-2所示。y0和yf分别表示0-1之间数值(正极liymn2o4)，x0和xf分别表示0-1之间数值(负极lixc6)。纵坐标up和un分别表示正极电压和负极电压。由图1可知：(1)0附近的初始y0,new和循环老化y0,aged数值相近，1附近的初始yf,new和循环老化yf,aged数值偏差较大；(2)0附近的初始x0,new和循环老化x0,aged数值相近，1附近的初始xf,new和循环老化xf,aged数值偏差较大。[0192]实验结果3[0193]图3-1至图3-5为基于5个老化实验周期的soh估计值和拟合曲线。(1)如图3-1至图3-3，当充放电倍率为2/3c时，单体电池的soh估计值与实测值误差较小，初始soc值为35％时的误差值最小；(2)如图3-4和图3-5，当充放电倍率为5c时，吞吐量在15kah-20kah之间的soh估计值与实测值误差较大；(3)soc初始值影响初始soh估计精度，充放电倍率影响过程soh估计精度。