当前位置：首页 > news >正文

二阶RC电池模型参数在线辨识：基于BMS电池管理系统的实践

news 2026/7/3 18:13:24

二阶RC电池模型参数在线辨识（BMS电池管理系统）使用遗忘因子最小二乘法 FFRLS 对电池模型进行参数辨识，并利用辨识的参数进行端电压的实时验证，基于动态工况，电压误差不超过20mv，也可以用来与离线辨识做对比，效果见图内容包含做电池Simulink模型、电芯数据、推导公式、参考论文程序已经调试好，可直接运行，也可以替换成自己的数据

在电池管理系统（BMS）领域，对电池模型参数进行精确辨识至关重要。本文将深入探讨如何使用遗忘因子最小二乘法（FFRLS）对二阶RC电池模型进行参数在线辨识，并通过实时验证端电压来评估辨识效果。

一、二阶RC电池模型与Simulink模型搭建

二阶RC电池模型是一种常用的等效电路模型，它能较好地模拟电池的动态特性。在Simulink中搭建此模型时，需要考虑电池的欧姆内阻、极化电阻和极化电容等元件。

例如，我们可以使用Simulink的电气元件库来构建电路结构，代码实现（以MATLAB脚本辅助生成部分参数为例）：

% 定义部分电池模型参数 R0 = 0.01; % 欧姆内阻 R1 = 0.05; % 极化电阻1 C1 = 1000; % 极化电容1 R2 = 0.03; % 极化电阻2 C2 = 800; % 极化电容2

在Simulink模型中，将这些参数连接到对应的电阻、电容模块，构建出完整的二阶RC电路结构。这样就完成了电池Simulink模型的初步搭建。

二、电芯数据的获取与使用

电芯数据是进行参数辨识的基础。这些数据通常包含电池在不同工况下的充放电电流、端电压以及对应的时间戳等信息。

假设我们已经从实验中获取了电芯数据，存储在一个MATLAB矩阵中，每一行代表一个时间点的数据，第一列是电流，第二列是电压，代码如下：

% 假设data矩阵已从外部导入，包含电流和电压数据 data = load('cell_data.txt'); current = data(:, 1); voltage = data(:, 2);

这些数据将作为后续参数辨识算法的输入，用于寻找最佳的模型参数。

三、遗忘因子最小二乘法（FFRLS）推导公式

遗忘因子最小二乘法是一种递推算法，它在传统最小二乘法的基础上引入遗忘因子，以更好地跟踪时变系统的参数变化。

其基本推导公式如下：

设系统的输出模型为 \(y(k) = \boldsymbol{\theta}^T\boldsymbol{\varphi}(k) + v(k)\)，其中 \(y(k)\) 是系统在 \(k\) 时刻的输出，\(\boldsymbol{\theta}\) 是待辨识参数向量，\(\boldsymbol{\varphi}(k)\) 是数据向量，\(v(k)\) 是噪声。

递推公式为：

\(\boldsymbol{K}(k) = \frac{\boldsymbol{P}(k - 1)\boldsymbol{\varphi}(k)}{\lambda + \boldsymbol{\varphi}^T(k)\boldsymbol{P}(k - 1)\boldsymbol{\varphi}(k)}\)

\(\hat{\boldsymbol{\theta}}(k) = \hat{\boldsymbol{\theta}}(k - 1) + \boldsymbol{K}(k)[y(k) - \boldsymbol{\varphi}^T(k)\hat{\boldsymbol{\theta}}(k - 1)]\)

\(\boldsymbol{P}(k) = \frac{1}{\lambda}[\boldsymbol{P}(k - 1) - \boldsymbol{K}(k)\boldsymbol{\varphi}^T(k)\boldsymbol{P}(k - 1)]\)

其中，\(\lambda\) 是遗忘因子，\(\boldsymbol{K}(k)\) 是增益矩阵，\(\hat{\boldsymbol{\theta}}(k)\) 是 \(k\) 时刻的参数估计值，\(\boldsymbol{P}(k)\) 是协方差矩阵。

在我们的二阶RC电池模型参数辨识中，将电池的端电压作为输出 \(y(k)\)，电流等相关量构成数据向量 \(\boldsymbol{\varphi}(k)\)，通过不断迭代上述公式，即可得到模型参数的估计值。

四、参数辨识与端电压实时验证

基于上述公式，编写FFRLS算法的MATLAB代码进行参数辨识：

% 初始化参数 lambda = 0.98; % 遗忘因子 theta_hat = zeros(5, 1); % 初始参数估计值 P = 1000 * eye(5); % 初始协方差矩阵 for k = 1:length(current) phi = [1, current(k), -exp(-1/(R1*C1)), -exp(-1/(R2*C2)), current(k)*(1 - exp(-1/(R1*C1)))/(R1*C1) + current(k)*(1 - exp(-1/(R2*C2)))/(R2*C2)]; K = P * phi / (lambda + phi' * P * phi); theta_hat = theta_hat + K * (voltage(k) - phi' * theta_hat); P = (1/lambda) * (P - K * phi' * P); end

得到辨识参数后，利用这些参数在Simulink模型中计算端电压，并与实际测量的端电压进行对比验证。在动态工况下，我们期望电压误差不超过20mV。