1.一种电池等效电路参数辨识方法,所述等效电路包括依次串联的恒压源、电阻R0、第一RC回路和第二RC回路,其中第一RC回路由电阻R1和电容C1并联构成,第二RC回路由电阻R2和电容C2并联构成,其特征在于,包括以下步骤:S1.对静置后的电池依次进行恒流脉冲放电和搁置测试,记录测试过程的时间序列、端电压序列和电流序列,获得放电电流Ip和脉冲时长tp,以静置结束时刻为零时刻,并以零时刻的端电压作为恒压源电压Voc;
S2.确定待辨识参数向量θ=[θ1,θ2,θ3,θ4,θ5]=[R0,R1,C1,R2,C2]及其下界LBj=L×(prior) (prior)θj 、上界UBj=U×θj ,其中L、U为预设系数且L
S3.生成满足拉丁超立方设计的N×5维均匀分布矩阵X,其中每列X{:,j}在区间[0,1]上均匀分布且分层间隔为1/N;
(i)
在对数空间构造参数样本集,其各参数样本θj 满足下式:(i)
log10(θj )=log10(LBj)+Xi,j×[log10(UBj)‑log10(LBj)]以上i取值为1至N,j取值为1至5;
(i) (i) (i) (i)
S4.将各参数样本θj 空间变换为φj :φj =ln(θj );
按下式构造电路模型fb:
(i) (i) (i) (i) (i) (i)当0
(i) (i) (i) (i)
当t>tp时,Vp(t,θj )=Voc‑V1cexp[‑(t‑tp)/(θ2 θ3 )]‑V2cexp[‑(t‑tp)/(θ4(i)θ5 )],其中:
(i) (i) (i) (i) (i) (i)V1c=θ2 Ip[1‑exp(‑tp/(θ2 θ3 ))],V2c=θ4 Ip[1‑exp(‑tp/(θ4 θ5 ))];
以上t为时间,I(t)为时间t对应的电流测试值,Vp为端电压预测值;
按下式基于fb模型构建电路变换模型:
(i) (i)
Vp(t,φj )=fb(t,exp(φj ),I(t),Voc,Ip,tp)(i)
S5.通过调用多个核心进行并行计算,每个核心针对不同样本φj 执行优化:生成最优解集
其中测试过程时间序列为[t1,t2,…,tk,…,tK],Vm(tk)为时间tk对应的端电压测试值,(i)Vm(tk,φj )为由电路变换模型获得的端电压预测值;
S6.将最优解集 按下式转化回物理空间:(i) (i)
θj =exp(φj );
(i)
S7.从步骤S6中获得的各解θj 中找出最适解作为最终辨识结果。
2.根据权利要求1所述的电池等效电路参数辨识方法,其特征在于,所述电池为锂离子电池或钠离子电池。
3.根据权利要求1所述的电池等效电路参数辨识方法,其特征在于,所述测试过程对时间、端电压和电流进行同步采样且采样时间间隔保持不变。
4.根据权利要求1或3所述的电池等效电路参数辨识方法,其特征在于,所述测试过程的采样时间间隔在0.01s至1s之间。
5.根据权利要求1所述的电池等效电路参数辨识方法,其特征在于,所述步骤S2中的预设系数L在0.0001至0.1之间,预设系数U在10至10000之间。
6.根据权利要求1所述的电池等效电路参数辨识方法,其特征在于,所述步骤S5中用于并行计算的核心数目大于4。
7.根据权利要求1所述的电池等效电路参数辨识方法,其特征在于,所述步骤S5中的并(i)行计算过程,将各样本φj 构成的样本集平均分配至各计算核心,各计算核心依据其分配得到的样本执行独立的优化计算,并将结果汇总成最优解集。
8.根据权利要求1所述的电池等效电路参数辨识方法,其特征在于,所述步骤S7中,将(i)各解θj 分别代入电路模型fb中获得其对应的端电压预测值序列,并与端电压实测值序列比较以评价各解的优劣,从而找出最适解作为最终辨识结果。
9.根据权利要求8所述的电池等效电路参数辨识方法,其特征在于,评价各解的优劣的方式为,通过端电压预测值序列和端电压实测值序列计算得到残差均方根值RMSE或拟合优2
度R,以RMSE最小者对应的解或R2最大者对应的解作为最适解。
10.一种计算机可读取存储介质,存储有多条计算机指令,其特征在于,所述计算机指令适于处理器进行加载,以执行权利要求1‑9任一项所述的电池等效电路参数辨识方法。