基于改进Thevenin模型的EKF算法模型

**锂电池二阶RC等效电路参数辨识(改进Thevenin模型)***

相比于1阶RC等效电路模型，改进的Thevenin等效电路非线性特性好，一阶RC等效电路能快速便捷模拟锂电池的动态特征，二阶RC等效电路能模拟动力锂电池系统动态和静态特征。因此需要辨识的参数也随之增加，从而增大计算量。

二阶RC等效电路用RΩ表示欧姆压降特征，一个RC网络表示电化学极化特征，另一个RC网络表示浓差极化特征。二阶RC等效电路图如下所示（PS：电池Em测试难度大，一般采用脉冲充放电静置稳定后的OCV为电动势Em，因此OCV ≈ Em）。（参考百度）

二阶RC等效电路主要由元器件RΩ、Rd、Cp、Rp、Ce、Re、Uocv恒压源组成，其中：

（1）RΩ为欧姆内阻

（2）Rd为自放电内阻

（3）Cp为电化学极化电容

（4）Rp为电化学极化内阻

（5）Ce为浓差极化电容

（6）Re为浓差极化内阻

（7）Uocv为开路电压≈ Em电池电动势，（利用查表法OCV-SOC）

（8）Uload为端电压

（9）Up电化学极化内阻端电压

（10）Ue为浓差极化内阻端电压

Simulink二阶rc等效模型的建立：

2RC等效模型的扩展卡尔曼滤波（EKF）的SOC估计模型如下：

部分EKF代码

unction [sys,x0,str,ts]=mdlInitializeSizes
sizes=simsizes;
sizes.NumContStates=0;
sizes.NumDiscStates=0;
sizes.NumOutputs=2;
sizes.NumInputs=3;
sizes.DirFeedthrough=1;
sizes.NumSampleTimes=1;
***********
x0=[0 0 0 0 0]';
str=[];
ts=[-1,0];
%=====================================================
******************
**********
********
******
*****
****
***
**
Ls=L-M;
% J=0.005;
k=0.382;
T=0.002;%T=0&x5-per=pi/3&x5-per=2*pi/3&x5-per=pi&x5-per=4*pi/3&x5-per=5*pi/3&x5-per<2*pi
q=5*pi/3;
a11=-2/3;a12=4/pi;
a21=-2/3;a22=-2/pi;
a31=4/3;a32=-2/pi;

结果分析：

***算法的协防方根误差***

*****SOC估计结果****