基于二阶RC等效电路模型的FFRLS + EKF联合SOC估计探究-编程实验室

基于二阶RC等效电路模型的FFRLS+EKF的联合SOC估计具体思路：采用遗忘因子最小二乘法估计电池参数，并将辨识结果导入到扩展卡尔曼滤波EKF算法中，实现FFRLS+EKF的联合估计，基于动态工况能保证运行，simulink模型和仿真结果可见展示图片，估计效果能完全跟随soc的变化内容：simulink模型、电池数据、参考文献，详细遗忘因子最小二乘法、EKF估计电池soc原理

在电池管理系统中，准确估计电池的荷电状态（SOC）至关重要。本文将探讨基于二阶RC等效电路模型，运用遗忘因子最小二乘法（FFRLS）与扩展卡尔曼滤波（EKF）联合估计SOC的方法。

一、二阶RC等效电路模型简介

二阶RC等效电路模型能够较好地模拟电池动态特性。它包含一个电压源、一个内阻以及两个RC并联支路。这种模型可以更精确地反映电池在充放电过程中的暂态响应。

二、遗忘因子最小二乘法（FFRLS）

（一）原理

遗忘因子最小二乘法的核心思想是在传统最小二乘法基础上引入遗忘因子，以更有效地跟踪时变系统参数。对于电池系统，电池参数会随着使用时间、温度等因素发生变化，FFRLS能适应这种变化。

假设我们有一组测量数据\(y(k)\)，并且假设系统的输出可以表示为输入\(u(k)\)与未知参数向量\(\theta\)的线性组合：

\[y(k)=\varphi^T(k)\theta + v(k)\]

其中\(\varphi(k)\)是回归向量，\(v(k)\)是噪声。

传统最小二乘法通过最小化误差平方和来估计参数\(\theta\)：

\[J(\theta)=\sum_{i = 1}^{k}(y(i)-\varphi^T(i)\theta)^2\]

而遗忘因子最小二乘法在计算时给旧数据加上遗忘因子\(\lambda(0<\lambda\leq1)\)，使得新数据具有更大权重。其目标函数变为：

\[J{\lambda}(\theta)=\sum{i = 1}^{k}\lambda^{k - i}(y(i)-\varphi^T(i)\theta)^2\]

（二）代码示例（简单示意，Python）

import numpy as np def ffrls(y, u, lambda_=0.98, P=None, theta=None): n = len(y) if P is None: P = np.eye(len(u[0])) * 10000 if theta is None: theta = np.zeros(len(u[0])) for k in range(n): phi = np.array(u[k]).reshape(-1, 1) K = P.dot(phi) / (lambda_ + phi.T.dot(P).dot(phi)) theta = theta + K.flatten() * (y[k] - phi.T.dot(theta)) P = (P - K.dot(phi.T).dot(P)) / lambda_ return theta

（三）代码分析

在上述代码中，首先对遗忘因子\(\lambda\)、初始协方差矩阵\(P\)和初始参数估计\(\theta\)进行了初始化。然后通过循环，每次根据当前测量数据更新参数估计值\(\theta\)和协方差矩阵\(P\)。K为卡尔曼增益，它通过当前协方差矩阵\(P\)、回归向量\(\phi\)和遗忘因子\(\lambda\)计算得出。新的参数估计\(\theta\)则是在旧估计值基础上，根据测量值与预测值的误差进行修正。

三、扩展卡尔曼滤波（EKF）估计电池SOC原理

（一）原理

EKF主要用于处理非线性系统的状态估计问题。对于电池SOC估计，电池的动态特性是非线性的，因此适合使用EKF。

我们定义系统状态方程和观测方程。以二阶RC等效电路模型为例，状态方程可表示为：

\[\begin{bmatrix}

SOC(k + 1)\\

x_1(k + 1)\\

x_2(k + 1)

\end{bmatrix}=\begin{bmatrix}

1 - \frac{\Delta t}{C_{bat}} & 0 & 0\\

0 & e^{-\frac{\Delta t}{R1C1}} & 0\\

0 & 0 & e^{-\frac{\Delta t}{R2C2}}

\end{bmatrix}\begin{bmatrix}

SOC(k)\\

x_1(k)\\

x_2(k)

\end{bmatrix}+\begin{bmatrix}

\frac{\Delta t}{C_{bat}}\\

0\\

\end{bmatrix}I(k)\]

其中\(SOC\)为荷电状态，\(x1\)、\(x2\)为RC支路的状态变量，\(\Delta t\)为采样时间，\(C{bat}\)为电池容量，\(R1\)、\(R2\)、\(C1\)、\(C_2\)为电路参数，\(I\)为电流。

观测方程为：

\[V(k)=OCV(SOC(k)) - R0I(k)-x1(k)-x_2(k)\]

其中\(V\)为电池端电压，\(OCV\)为开路电压。

EKF通过对状态方程和观测方程进行线性化处理，利用卡尔曼滤波的框架进行状态估计。

（二）代码示例（简单示意，Python）

import numpy as np def ekf(y, u, dt, Q, R, x_hat_prev, P_prev, params): # 解包参数 C_bat, R_0, R_1, C_1, R_2, C_2 = params F = np.array([[1 - dt / C_bat, 0, 0], [0, np.exp(-dt / (R_1 * C_1)), 0], [0, 0, np.exp(-dt / (R_2 * C_2))]]) B = np.array([dt / C_bat, 0, 0]).reshape(-1, 1) H = np.array([1, -1, -1]).reshape(1, -1) # 预测步骤 x_hat_minus = F.dot(x_hat_prev) + B * u P_minus = F.dot(P_prev).dot(F.T) + Q # 更新步骤 K = P_minus.dot(H.T) / (H.dot(P_minus).dot(H.T) + R) x_hat = x_hat_minus + K.flatten() * (y - H.dot(x_hat_minus)) P = (np.eye(3) - K.dot(H)).dot(P_minus) return x_hat, P

（三）代码分析

在这段代码中，首先根据输入的电池参数构建状态转移矩阵\(F\)、控制输入矩阵\(B\)和观测矩阵\(H\)。然后进入预测步骤，利用上一时刻的状态估计值\(xhatprev\)和协方差矩阵\(Pprev\)预测当前时刻的状态估计值\(xhatminus\)和协方差矩阵\(Pminus\)。接着在更新步骤中，根据测量值\(y\)和预测值\(xhatminus\)计算卡尔曼增益\(K\)，进而更新状态估计值\(x_hat\)和协方差矩阵\(P\)。

四、联合估计实现

我们采用FFRLS估计电池参数，将辨识结果导入到EKF算法中，实现联合估计。具体步骤为：先利用FFRLS对电池的内阻、RC参数等进行估计，然后将这些估计参数代入到EKF的状态方程和观测方程中，进行SOC的估计。

五、Simulink模型

通过在Simulink中搭建二阶RC等效电路模型，并结合FFRLS和EKF模块实现联合估计。模型主要包含电池模型模块、FFRLS参数估计模块、EKF估计模块以及数据采集和显示模块。在电池模型模块中设置电池的基本参数，FFRLS模块根据输入的电流、电压数据估计电池参数，EKF模块利用FFRLS估计的参数和输入数据进行SOC估计。