使用自定义函数拟合辨识HPPC工况下的电池数据(适用于一阶RC、二阶RC等电池模型)

本文主要是介绍使用自定义函数拟合辨识HPPC工况下的电池数据(适用于一阶RC、二阶RC等电池模型),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

该程序可以离线辨识HPPC工况下的电池数据,只需要批量导入不同SOC所对应的脉冲电流电压数据,就可以瞬间获得SOC为[100% 90% 80% 70% 60% 50% 40% 30% 20% 10% 0%]的所有电池参数,迅速得到参数辨识的结果并具有更高的精度,可以很大程度上降低参数辨识的时间,提高参数辨识效率。

1.首先是代码中要更改的参数,我们要根据我们不同的模型选择不同的初始值。

如果是一阶模型iniPar = [0.06;0.03;1000]; 0.06是R0的初始值,0.03是R1的初始值,1000是C1的初始值。

如果是二阶模型iniPar = [0.05;0.05;2000;0.005;800]; 0.05是R0的初始值,0.05是R1的初始值,2000是C1的初始值,0.005是R2的初始值,800是C2的初始值。

如果是1RC+H模型,iniPar = [0.05;0.01;1000;0.001;0.1];

如果是2RC+H模型,iniPar = [0.05;0.01;1500;0.003;200;0.000005;50]

%% 代码中要更改的参数
currentLimit = 10; % 高于此限值的电流将被视为1C电流  
iniPar = [0.06;0.03;1000]; % 1RC
%iniPar = [0.05;0.05;2000;0.005;800]; % 2RC
%iniPar = [0.05;0.01;1000;0.001;0.1]; % 1RC+H
%iniPar = [0.05;0.01;1500;0.003;200;0.000005;50]; % 2RC+H

2.模型阶数的选择 

%% Parameters used in algorithm 
xIni = 0; % 1RC
%xIni = [0; 0]; % 2RC and 1RC+H
%xIni = [0; 0; 0]; % 2RC+H
[m,n] = size(dataHPPC{1});
outputPar = zeros(3,n); % %3 - 1RC, 5 - 1RC+H and 2RC, 7 - 2RC+H
modelError = zeros(m,n); 
vModel = zeros(m,n); 
model = 1; %1 - c1RC, 2 - 2RC, 3 - 1RC+H, 4 - 2RC+H

3.运行function函数,求解参数

所有的function函数如下图所示:

二阶RC模型的function函数如下:

function [vModel,XTimUp1,XTimUp2] = Model_2RC (current, tSample, xPrev1,xPrev2, ocvVoltage, ECM_Parameters)% Model Parameters
r0 = ECM_Parameters(1);
r1 = ECM_Parameters(2);
c1 = ECM_Parameters(3);
r2 = ECM_Parameters(4);
c2 = ECM_Parameters(5);tau1 = r1*c1;
tau2 = r2*c2;%% State Time Update
XTimUp1 = exp(-tSample/(tau1))*xPrev1 + r1*(1-exp(-tSample/tau1))*current;
XTimUp2 = exp(-tSample/(tau2))*xPrev2 + r2*(1-exp(-tSample/tau2))*current;%% The voltage response from the model
vModel = ocvVoltage-XTimUp1-XTimUp2-r0*current;end
for i = 1:ntime = dataHPPC{1}(:,i);current = dataHPPC{2}(:,i);vExp = dataHPPC{3}(:,i);% Find the ocvCurve value right before the HPPC test starts index = find(current>currentLimit);ocvVoltage = vExp(index(1)-0);% Obtaining the parameter estimatesfun = @(beta,x)ObjectiveFunction(beta,x,ocvVoltage,xIni,model);outputPar(:,i) = nlinfit([time,current],vExp,fun,iniPar);% Running the voltage model with the given battery parametersvModel(:,i) = ObjectiveFunction(outputPar(:,i),[time,current],ocvVoltage,xIni,model);modelError(:,i) =  abs((vModel(:,i) - vExp)./vExp)*100;

主要使用了MATLAB的nlinfit函数。

nlinfit函数的用法为:beta = nlinfit(X, Y, modelfun, beta0)

X与Y分别是时间与电压

beta:估计出的非线性函数的系数,即我们辨识的参数值

modelfun:指定的要拟合的非线性函数,即y=a-b*exp(-c*t)-d*exp(-f*t)

beta0:回归系数的初始值,即待辨识参数的初始值。

辨识得到的参数结果保存在变量outputPar里面,从左到右依次对应的是不同SOC下所对应的参数,从上到下的五行分别是参数R0、R1、R2、C1、C2。

 绘制不同SOC的参数拟合结果如下所示:

4. 绘图,包括根据辨识的参数值仿真出的端电压与真实测量端电压的对比图、端电压的误差

    figureplot(vModel(:,i))hold onplot(vExp)polt(vModel(:,i)-vExp)

二阶RC模型的端电压拟合情况对比

一阶RC模型的端电压拟合情况对比:

 

这篇关于使用自定义函数拟合辨识HPPC工况下的电池数据(适用于一阶RC、二阶RC等电池模型)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/350166

相关文章

使用Python从PPT文档中提取图片和图片信息(如坐标、宽度和高度等)

《使用Python从PPT文档中提取图片和图片信息(如坐标、宽度和高度等)》PPT是一种高效的信息展示工具,广泛应用于教育、商务和设计等多个领域,PPT文档中常常包含丰富的图片内容,这些图片不仅提升了... 目录一、引言二、环境与工具三、python 提取PPT背景图片3.1 提取幻灯片背景图片3.2 提取

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

使用Python实现图像LBP特征提取的操作方法

《使用Python实现图像LBP特征提取的操作方法》LBP特征叫做局部二值模式,常用于纹理特征提取,并在纹理分类中具有较强的区分能力,本文给大家介绍了如何使用Python实现图像LBP特征提取的操作方... 目录一、LBP特征介绍二、LBP特征描述三、一些改进版本的LBP1.圆形LBP算子2.旋转不变的LB

Maven的使用和配置国内源的保姆级教程

《Maven的使用和配置国内源的保姆级教程》Maven是⼀个项目管理工具,基于POM(ProjectObjectModel,项目对象模型)的概念,Maven可以通过一小段描述信息来管理项目的构建,报告... 目录1. 什么是Maven?2.创建⼀个Maven项目3.Maven 核心功能4.使用Maven H

Python中__init__方法使用的深度解析

《Python中__init__方法使用的深度解析》在Python的面向对象编程(OOP)体系中,__init__方法如同建造房屋时的奠基仪式——它定义了对象诞生时的初始状态,下面我们就来深入了解下_... 目录一、__init__的基因图谱二、初始化过程的魔法时刻继承链中的初始化顺序self参数的奥秘默认

SpringBoot使用GZIP压缩反回数据问题

《SpringBoot使用GZIP压缩反回数据问题》:本文主要介绍SpringBoot使用GZIP压缩反回数据问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录SpringBoot使用GZIP压缩反回数据1、初识gzip2、gzip是什么,可以干什么?3、Spr

Spring Boot 集成 Quartz并使用Cron 表达式实现定时任务

《SpringBoot集成Quartz并使用Cron表达式实现定时任务》本篇文章介绍了如何在SpringBoot中集成Quartz进行定时任务调度,并通过Cron表达式控制任务... 目录前言1. 添加 Quartz 依赖2. 创建 Quartz 任务3. 配置 Quartz 任务调度4. 启动 Sprin

Linux下如何使用C++获取硬件信息

《Linux下如何使用C++获取硬件信息》这篇文章主要为大家详细介绍了如何使用C++实现获取CPU,主板,磁盘,BIOS信息等硬件信息,文中的示例代码讲解详细,感兴趣的小伙伴可以了解下... 目录方法获取CPU信息:读取"/proc/cpuinfo"文件获取磁盘信息:读取"/proc/diskstats"文

Java使用SLF4J记录不同级别日志的示例详解

《Java使用SLF4J记录不同级别日志的示例详解》SLF4J是一个简单的日志门面,它允许在运行时选择不同的日志实现,这篇文章主要为大家详细介绍了如何使用SLF4J记录不同级别日志,感兴趣的可以了解下... 目录一、SLF4J简介二、添加依赖三、配置Logback四、记录不同级别的日志五、总结一、SLF4J

使用Python实现一个优雅的异步定时器

《使用Python实现一个优雅的异步定时器》在Python中实现定时器功能是一个常见需求,尤其是在需要周期性执行任务的场景下,本文给大家介绍了基于asyncio和threading模块,可扩展的异步定... 目录需求背景代码1. 单例事件循环的实现2. 事件循环的运行与关闭3. 定时器核心逻辑4. 启动与停