CN107144793B - 一种动力电池的分数阶KiBaM模型参数辨识方法及系统 - Google Patents

一种动力电池的分数阶KiBaM模型参数辨识方法及系统 Download PDF

Info

Publication number
CN107144793B
CN107144793B CN201710304516.9A CN201710304516A CN107144793B CN 107144793 B CN107144793 B CN 107144793B CN 201710304516 A CN201710304516 A CN 201710304516A CN 107144793 B CN107144793 B CN 107144793B
Authority
CN
China
Prior art keywords
power battery
fractional order
model
identification
linear
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710304516.9A
Other languages
English (en)
Other versions
CN107144793A (zh
Inventor
李岩
赵阳
张承慧
段彬
周忠凯
李泽元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN201710304516.9A priority Critical patent/CN107144793B/zh
Publication of CN107144793A publication Critical patent/CN107144793A/zh
Application granted granted Critical
Publication of CN107144793B publication Critical patent/CN107144793B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367—Software therefor, e.g. for battery testing using modelling or look-up tables

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Secondary Cells (AREA)

Abstract

本发明公开一种动力电池的分数阶KiBaM模型参数辨识方法及系统，其中，该辨识方法包括步骤1：初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；步骤2：利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；步骤3：采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数；否则，返回步骤1。

Description

一种动力电池的分数阶KiBaM模型参数辨识方法及系统
技术领域
本发明属于动力电池建模领域，尤其涉及一种动力电池的分数阶KiBaM模型参数辨识方法及系统。
背景技术
随着世界技术的革新，面临石油消耗的危机和大气环境污染的加剧，增加了对混合动力汽车、纯电动汽车、燃料电池汽车等电动汽车研发的投入与政策的扶持力度成为各国能源发展战略的重要举措。在新能源发展体系中，动力电池作为电动汽车动力系统、光伏电站和储能电站能源存储的关键部分，如何更好降低动力电池的生产成本，延长动力电池寿命和充分使用电池可用容量便成为了当前的研究热点。电池荷电状态估计是电动汽车车用动力电池系统研究的难点之一，电池剩余电量估算的精度直接关系到电动汽车的使用性能，而电池模型的建立是电池荷电状态估计的基础。电池模型描述了电池内部状态变量与外部特性之间的特定关系，准确的电池荷电状态估计是需要建立一个有效的电池模型作为保障的。
现有常用的动力电池模型有两种：无量纲的电化学模型和有量纲的等效电路模型。电化学模型对应电池内部微观尺度的动力学信息，即使应用首性原理(firstprinciple)仍很难与电压、电流等外特性量化关联。
但是，现有的分数阶模型辨识要么在频域内采用阻抗谱分析、伯德图、奈奎斯特曲线等方法，利用电化学工作站等动力电池测试仪器获得无损激励下动力电池的精确模型；要么通过获取电压、电流等外部数据应用整数阶建模的方法在时域内实现分数阶模型参数的辨识。二者各有利弊且均存在一定的缺陷：首先，阻抗谱分析针对的是半电池，其输入必须是对电池无损的持续激励，不适用于实际工况，其应用局限在电池的研发阶段。第二，时域内的模型辨识可以获得非常精确的模型，但是分数阶模型的时域辨识缺乏相应的理论方法因此存在一定困难，如利用Grunwald-Letnikov等方法获得的离散模型存在计算量巨大从而无法应用等问题)。第三，分数阶阶次的辨识是分数阶模型辨识的关键环节，然而现有的分数阶阶次辨识方法多采用区间分割法容易陷入伪阶次的情况，极大地影响了模型辨识的精度。
发明内容
为了解决现有技术的缺点，本发明的第一目的是提供一种动力电池的分数阶KiBaM模型参数辨识方法。该方法能够保留全部模型参数的物理意义，利用分数阶阶次的全参数自适应辨识方法，在其他参数未精确已知时能够通过迭代学习获得精确的阶次，不仅能够提高模型的精度，同时满足实际工况下动力电池测试、模拟和BMS的实际需求，更进一步通过结合其他技术能够更好的实现动力电池内部多状态的精确综合估计。
本发明的一种动力电池的分数阶KiBaM模型参数辨识方法，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力的电池时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；所述动力电池的分数阶KiBaM模型参数辨识方法包括：
步骤1：初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
步骤2：利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
步骤3：采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数；否则，返回步骤1。
进一步的，在所述步骤1中，利用最小二乘辨识方法辨识非线性参数向量。
进一步的，在所述步骤2中，利用最小二乘辨识方法辨识线性参数向量。
在辨识和参数估计领域中，最小二乘辨识法即可用于动态系统也可以用于静态系统，即可用于线性系统也可用于非线性系统，即可用于离线估计又可用于在线估计；而且在随机的环境下，利用最小二乘法时，并不要求观测数据提供出有关它的概率统计方面的信息，而用这种方法所获得的估计结果，有相当好的统计特性，而且利用最小二乘原理所拟定的辨识算法在实施上比较简单。
本发明的第二目的是提供一种动力电池的分数阶KiBaM模型参数辨识系统。该系统能够保留全部模型参数的物理意义，利用分数阶阶次的全参数自适应辨识，在其他参数未精确已知时能够通过迭代学习获得精确的阶次，不仅能够提高模型的精度，同时满足实际工况下动力电池测试、模拟和BMS的实际需求，更进一步通过结合其他技术能够更好的实现动力电池内部多状态的精确综合估计。
本发明的一种动力电池的分数阶KiBaM模型参数辨识系统，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力的电池时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；所述动力电池的分数阶KiBaM模型参数辨识系统包括：
非线性参数向量辨识模块，其用于初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
线性参数向量辨识模块，其用于利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
模型最优参数获取模块，其用于采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数。
进一步的，在所述非线性参数向量辨识模块中，利用最小二乘辨识方法辨识非线性参数向量。
进一步的，在所述线性参数向量辨识模块中，利用最小二乘辨识方法辨识线性参数向量。
在辨识和参数估计领域中，最小二乘辨识法即可用于动态系统也可以用于静态系统，即可用于线性系统也可用于非线性系统，即可用于离线估计又可用于在线估计；而且在随机的环境下，利用最小二乘法时，并不要求观测数据提供出有关它的概率统计方面的信息，而用这种方法所获得的估计结果，有相当好的统计特性，而且利用最小二乘原理所拟定的辨识算法在实施上比较简单。
本发明还提供了另一种动力电池的分数阶KiBaM模型参数辨识系统。
本发明的一种动力电池的分数阶KiBaM模型参数辨识系统，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力的电池时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；所述动力电池的分数阶KiBaM模型参数辨识系统包括：
数据采集装置，其被配置为：采集动力电池的实际输出并传送至服务器；
所述服务器，其被配置为：
初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数。
进一步的，所述服务器，还被配置为：利用最小二乘辨识方法辨识非线性参数向量。
进一步的，所述服务器，还被配置为：利用最小二乘辨识方法辨识线性参数向量。
在辨识和参数估计领域中，最小二乘辨识法即可用于动态系统也可以用于静态系统，即可用于线性系统也可用于非线性系统，即可用于离线估计又可用于在线估计；而且在随机的环境下，利用最小二乘法时，并不要求观测数据提供出有关它的概率统计方面的信息，而用这种方法所获得的估计结果，有相当好的统计特性，而且利用最小二乘原理所拟定的辨识算法在实施上比较简单。
进一步的，所述服务器还与显示装置相连，所述显示装置用于展示动力电池的分数阶KiBaM模型的最优参数。
本发明的有益效果为：
(1)本发明的该方法保留了全部模型参数的物理意义，如KiBaM非线性项对应着动力电池的欧姆内阻，W和CPE分别对应动力电池内部不同位置的扩散性质等；该方法适用于Caputo，Riemann-Liouville，Grunwald-Letnikov等多种分数阶导数的定义，且兼容各式分数阶算子的高精度计算方法；该方法首次提出了分数阶阶次的全参数自适应辨识方法，在其他参数未精确已知时能够通过迭代学习获得精确的阶次。
(2)本方法不仅能够提高模型的精度，同时满足实际工况下动力电池测试、模拟和BMS的实际需求，更进一步通过结合其他技术能够更好的实现动力电池内部多状态的精确综合估计。
(3)本方法的迭代学习分数阶阶次辨识方法是全参数自适应的，避免了区间分割法可能会造成的伪阶次问题，从而有效减少了模型的标定次数。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。
图1是动力电池的分数阶KiBaM模型；
图2是动力电池的分数阶KiBaM模型参数辨识方法流程图；
图3是锂电池等效电路结构；
图4是辨识得到的动力电池模型输出与实际输出对比效果图；
图5是一种动力电池的分数阶KiBaM模型参数辨识系统结构示意图；
图6是另一种动力电池的分数阶KiBaM模型参数辨识系统结构示意图。
具体实施方式
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
鉴于充放电的需求，其负极的多孔(分形)性必然存在，从而引出了瓦尔堡阻抗(W元件)和CPE元件的概念，他们已成为动力电池等效电路模型的必备元件，而这两个元件均为分数阶的，因此分数阶等效电路模型蕴含着更多的物理意义，其参数与动力电池内部特性紧密关联，在实际应用中可有效避免模型的频繁标定。
1.1电池模型建立
动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成，如图1所示。
动态线性环节为如下等比例阶次分数阶传递函数形式：
其中，正整数r和h均假设已知，piα、pjα分别表示iα、jα阶次分数阶微分，bi、aj分别表示iα、jα阶次分数阶微分相应的系数且均为常数。
其中，动态线性环节对应动力电池的时间常数。分数阶阶次α对应电池内部电子的扩散。
式(1)中含有h+r+1个线性参数构成线性参数向量θl：
θl＝[ah，…，a1，br，…，b0]T。
非线性静态部分f(u)由多项式表示：
f(u)＝β0+β1u+β2u2+…+βmum。
其中，u表示动力电池输入量；βk表示非线性静态部分的系统，k＝0,1,…,m，m为自然数。
非线性参数向量θn由m+1个参数构成：
θn＝[β0 β1 β2 … βm]T
其中，非线性参数向量对应电池的非线性容量特征。
非线性参数向量的物理意义为电池不同部分间的扩散过程，即对应不同工作负载和电极材料的缓冲区电子。
1.2参数辨识问题的定义
电池模型建立好后，模型中参数辨识问题定义为：给定采集到的系统输入数据u(i)和输出数据y(i)，其中，i＝1,2,…,N；N为正整数。寻找动力电池的分数阶KiBaM模型的分数阶以及相关参数使得动力电池系统的成本函数达到最小值。
动力电池系统的成本函数为：
其中，v＝y-G(p)f(u)。
图2是动力电池的分数阶KiBaM模型参数辨识方法流程图。
如图2所示，本发明的动力电池的分数阶KiBaM模型参数辨识方法包括：
步骤1：初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量。
以图3所示的锂电池等效电路结构为例，N元件表示电池的非线性容量特征；RC环节(RC值)对应电池的时间常数；W元件对应电子在负极的扩散。
在步骤1中，利用最小二乘辨识方法辨识非线性参数向量。
对该结构建立分数阶KiBaM模型，其中Rs表示非线性环节，后部分表示线性环节。根据DST(动态压力测试)标准，采集得到的I/V数据作为待辨识模型的输入输出数据。测试根据锂电池性质中选取Warburg阻抗阶次为β＝0.5，等效电路中的其它参数分别为Rs＝19.1mΩ,Rp＝15.6mΩ,Y0＝1.24Mho,Y1＝370Mho,α＝0.665。由此得到交流阻抗为：
其中，非线性项Rs随误差变化。
步骤2：利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型。
在步骤2中，利用最小二乘辨识方法辨识线性参数向量。
在辨识和参数估计领域中，最小二乘辨识法即可用于动态系统也可以用于静态系统，即可用于线性系统也可用于非线性系统，即可用于离线估计又可用于在线估计；而且在随机的环境下，利用最小二乘法时，并不要求观测数据提供出有关它的概率统计方面的信息，而用这种方法所获得的估计结果，有相当好的统计特性，而且利用最小二乘原理所拟定的辨识算法在实施上比较简单。
步骤3：采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数；否则，返回步骤1。
在本实施例中，动力电池的电流信号作为整个动力电池系统的输入信号，动力电池的电压信号作为整个动力电池系统的输出信号。
实验辨识得到的分数阶模型与实际输出的对比效果见图4,其中Rs∈[34,37.7]。
本发明的该方法保留了全部模型参数的物理意义，如KiBaM非线性项对应着动力电池的欧姆内阻，W和CPE分别对应动力电池内部不同位置的扩散性质等；该方法适用于Caputo，Riemann-Liouville，Grunwald-Letnikov等多种分数阶导数的定义，且兼容各式分数阶算子的高精度计算方法；该方法首次提出了分数阶阶次的全参数自适应辨识方法，在其他参数未精确已知时能够通过迭代学习获得精确的阶次。
本方法不仅能够提高模型的精度，同时满足实际工况下动力电池测试、模拟和BMS的实际需求，更进一步通过结合其他技术能够更好的实现动力电池内部多状态的精确综合估计。
图5是一种动力电池的分数阶KiBaM模型参数辨识系统结构示意图；
如图5所示，动力电池的分数阶KiBaM模型参数辨识系统，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力的电池时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；所述动力电池的分数阶KiBaM模型参数辨识系统包括：
(1)非线性参数向量辨识模块，其用于初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
在所述非线性参数向量辨识模块中，利用最小二乘辨识方法辨识非线性参数向量。
(2)线性参数向量辨识模块，其用于利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
在所述线性参数向量辨识模块中，利用最小二乘辨识方法辨识线性参数向量。
在辨识和参数估计领域中，最小二乘辨识法即可用于动态系统也可以用于静态系统，即可用于线性系统也可用于非线性系统，即可用于离线估计又可用于在线估计；而且在随机的环境下，利用最小二乘法时，并不要求观测数据提供出有关它的概率统计方面的信息，而用这种方法所获得的估计结果，有相当好的统计特性，而且利用最小二乘原理所拟定的辨识算法在实施上比较简单。
(3)模型最优参数获取模块，其用于采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数。
在本实施例中，动力电池的电流信号作为整个动力电池系统的输入信号，动力电池的电压信号作为整个动力电池系统的输出信号。
本发明的动力电池的分数阶KiBaM模型参数辨识系统。该系统能够保留全部模型参数的物理意义，利用分数阶阶次的全参数自适应辨识，在其他参数未精确已知时能够通过迭代学习获得精确的阶次，不仅能够提高模型的精度，同时满足实际工况下动力电池测试、模拟和BMS的实际需求，更进一步通过结合其他技术能够更好的实现动力电池内部多状态的精确综合估计。
图6是另一种动力电池的分数阶KiBaM模型参数辨识系统结构示意图。
如图6所示，本发明的一种动力电池的分数阶KiBaM模型参数辨识系统，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力的电池时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；所述动力电池的分数阶KiBaM模型参数辨识系统包括：
(1)数据采集装置，其被配置为：采集动力电池的实际输出并传送至服务器；
(2)所述服务器，其被配置为：
初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数。
在本实施例中，动力电池的电流信号作为整个动力电池系统的输入信号，动力电池的电压信号作为整个动力电池系统的输出信号。
其中，所述服务器，还被配置为：利用最小二乘辨识方法辨识非线性参数向量。
所述服务器，还被配置为：利用最小二乘辨识方法辨识线性参数向量。
在辨识和参数估计领域中，最小二乘辨识法即可用于动态系统也可以用于静态系统，即可用于线性系统也可用于非线性系统，即可用于离线估计又可用于在线估计；而且在随机的环境下，利用最小二乘法时，并不要求观测数据提供出有关它的概率统计方面的信息，而用这种方法所获得的估计结果，有相当好的统计特性，而且利用最小二乘原理所拟定的辨识算法在实施上比较简单。
在另一实施例，服务器还与显示装置相连，所述显示装置用于展示动力电池的分数阶KiBaM模型的最优参数。
上述典例虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形(如分数阶高阶模型和分数阶Hammerstein-Wiener模型等)仍在本发明的保护范围以内。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (10)

1.一种动力电池的分数阶KiBaM模型参数辨识方法，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力电池的时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；其特征在于，所述动力电池的分数阶KiBaM模型参数辨识方法包括：
步骤1：初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
步骤2：利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
步骤3：采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数；否则，返回步骤1。
2.如权利要求1所述的一种动力电池的分数阶KiBaM模型参数辨识方法，其特征在于，在所述步骤1中，利用最小二乘辨识方法辨识非线性参数向量。
3.如权利要求1所述的一种动力电池的分数阶KiBaM模型参数辨识方法，其特征在于，在所述步骤2中，利用最小二乘辨识方法辨识线性参数向量。
4.一种动力电池的分数阶KiBaM模型参数辨识系统，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力电池的时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；其特征在于，所述动力电池的分数阶KiBaM模型参数辨识系统包括：
非线性参数向量辨识模块，其用于初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
线性参数向量辨识模块，其用于利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
模型最优参数获取模块，其用于采集动力电池的实际输出，计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数。
5.如权利要求4所述的一种动力电池的分数阶KiBaM模型参数辨识系统，其特征在于，在所述非线性参数向量辨识模块中，利用最小二乘辨识方法辨识非线性参数向量。
6.如权利要求4所述的一种动力电池的分数阶KiBaM模型参数辨识系统，其特征在于，在所述线性参数向量辨识模块中，利用最小二乘辨识方法辨识线性参数向量。
7.一种动力电池的分数阶KiBaM模型参数辨识系统，所述动力电池的分数阶KiBaM模型由一个非线性静态环节与线性动态环节串联构成；线性动态环节为动力电池的时间常数构成的线性参数向量，非线性静态环节为动力电池的非线性放电特性相关的非线性参数向量；其特征在于，所述动力电池的分数阶KiBaM模型参数辨识系统包括：
数据采集装置，其被配置为：采集动力电池的实际输出并传送至服务器；
所述服务器，其被配置为：
初始化线性参数向量，利用迭代学习辨识方法来辨识所述分数阶KiBaM模型的分数阶阶次，得到分数阶阶次初值，由此辨识非线性参数向量；
利用辨识的非线性参数向量，再结合分数阶阶次，辨识线性参数向量，进而确定出动力电池的分数阶KiBaM模型；
计算并判断所确定的动力电池的分数阶KiBaM模型输出与采集的动力电池实际输出偏差的范数是否达到最小值，若是，则结束辨识，得到动力电池的分数阶KiBaM模型的最优参数。
8.如权利要求7所述的一种动力电池的分数阶KiBaM模型参数辨识系统，其特征在于，所述服务器，还被配置为：利用最小二乘辨识方法辨识非线性参数向量。
9.如权利要求7所述的一种动力电池的分数阶KiBaM模型参数辨识系统，其特征在于，所述服务器，还被配置为：利用最小二乘辨识方法辨识线性参数向量。
10.如权利要求7所述的一种动力电池的分数阶KiBaM模型参数辨识系统，其特征在于，所述服务器还与显示装置相连，所述显示装置用于展示动力电池的分数阶KiBaM模型的最优参数。
CN201710304516.9A 2017-05-03 2017-05-03 一种动力电池的分数阶KiBaM模型参数辨识方法及系统 Active CN107144793B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710304516.9A CN107144793B (zh) 2017-05-03 2017-05-03 一种动力电池的分数阶KiBaM模型参数辨识方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710304516.9A CN107144793B (zh) 2017-05-03 2017-05-03 一种动力电池的分数阶KiBaM模型参数辨识方法及系统

Publications (2)

Publication Number Publication Date
CN107144793A CN107144793A (zh) 2017-09-08
CN107144793B true CN107144793B (zh) 2019-05-21

Family

ID=59774376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710304516.9A Active CN107144793B (zh) 2017-05-03 2017-05-03 一种动力电池的分数阶KiBaM模型参数辨识方法及系统

Country Status (1)

Country Link
CN (1) CN107144793B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108037459B (zh) * 2017-11-24 2022-02-11 北京机械设备研究所 动力电池模型最优阶次确定方法及系统
CN108519555A (zh) * 2018-04-11 2018-09-11 北京理工大学 一种锂离子电池的改进分数阶模型和参数辨识方法
CN108647434A (zh) * 2018-05-10 2018-10-12 燕山大学 一种基于改进粒子滤波算法的二元荷电状态估算方法
CN110333667A (zh) * 2019-08-14 2019-10-15 大连海事大学 一种基于非线性新息辨识模型的船舶运动仿真及控制方法
CN110750912B (zh) * 2019-10-28 2020-12-29 清华大学 一种基于无量纲化模型的电池热参数辨识方法
CN112486020A (zh) * 2020-12-01 2021-03-12 中国科学技术大学 基于自适应模型预测控制的混合储能系统功率分配系统及方法
CN113158606A (zh) * 2021-03-30 2021-07-23 闽南理工学院 一种考虑非线性容量效应的锂电池模型

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102937704A (zh) * 2012-11-27 2013-02-20 山东省科学院自动化研究所 一种动力电池rc等效模型的辨识方法
CN103941196A (zh) * 2014-05-07 2014-07-23 吉林大学 锂离子电池荷电状态估计方法
CN104112036A (zh) * 2014-06-12 2014-10-22 湖南文理学院 混联式混合动力电动汽车的仿真方法
CN104392080A (zh) * 2014-12-19 2015-03-04 山东大学 一种锂电池分数阶变阶等效电路模型及其辨识方法
CN105045971A (zh) * 2015-07-01 2015-11-11 西安交通大学 一种锂电池分数阶离散化阻抗模型

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102937704A (zh) * 2012-11-27 2013-02-20 山东省科学院自动化研究所 一种动力电池rc等效模型的辨识方法
CN103941196A (zh) * 2014-05-07 2014-07-23 吉林大学 锂离子电池荷电状态估计方法
CN104112036A (zh) * 2014-06-12 2014-10-22 湖南文理学院 混联式混合动力电动汽车的仿真方法
CN104392080A (zh) * 2014-12-19 2015-03-04 山东大学 一种锂电池分数阶变阶等效电路模型及其辨识方法
CN105045971A (zh) * 2015-07-01 2015-11-11 西安交通大学 一种锂电池分数阶离散化阻抗模型

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种估计锂电池充电状态的分数阶阻抗模型;杨晴霞 等;《西安交通大学学报》;20150831;第49卷(第8期);第128-132页
基于分数阶理论的车用锂离子电池建模及荷电状态估计;刘树林 等;《电工技术学报》;20170228;第32卷(第4期);第189-195页

Also Published As

Publication number Publication date
CN107144793A (zh) 2017-09-08

Similar Documents

Publication Publication Date Title
CN107144793B (zh) 一种动力电池的分数阶KiBaM模型参数辨识方法及系统
CN106855612B (zh) 计及非线性容量特性的分数阶KiBaM电池模型及参数辨识方法
CN104991980B (zh) 锂离子电池的电化学机理建模方法
Moura et al. PDE estimation techniques for advanced battery management systems—Part I: SOC estimation
CN105093114B (zh) 一种电池在线建模与荷电状态的联合估计方法及系统
CN104965179B (zh) 一种锂离子蓄电池的温度组合电路模型及其参数识别方法
CN109765496A (zh) 一种基于在线电化学阻抗谱测量的电池健康状态估计方法
CN108846227B (zh) 一种基于随机森林及容量自恢复效应分析的锂离子电池容量退化预测评估方法
WO2023274194A1 (zh) 一种适用于富锂锰基电池的高阶模型参数辨识方法和系统
CN107255787B (zh) 基于信息熵的电池组不一致性综合评价方法及系统
CN106250576A (zh) 一种基于动态阻抗的锂电池模型的建模方法
CN108761341A (zh) 一种锂离子电池电化学模型参数获取方法
CN109143083A (zh) 一种数据驱动的电动车辆锂离子电池析锂诊断方法
Nikolian et al. Classification of Electric modelling and Characterization methods of Lithium-ion Batteries for Vehicle Applications
CN103439666A (zh) 一种锂离子电池容量衰退评估的几何方法
CN106872901B (zh) KiBaM-分数阶等效电路综合特征电池模型及参数辨识方法
CN106777786A (zh) 一种锂离子电池soc估算方法
CN107064809A (zh) 计及电池电化学特性的动态等效电路模型及其工作方法
CN109613446A (zh) 一种基于时间序列分析的锂电池老化检测方法
CN103983920A (zh) 一种建立电动车辆的动力电池的模型的方法
CN106951605A (zh) 一种锂离子电池组等效模型构建方法
CN106093615A (zh) 超级电容储能模块的健康状态估计方法
CN105911480A (zh) 一种动力电池soc估计方法
Li et al. A framework for states co-estimation of hybrid energy storage systems based on fractional-order theory
Wang et al. Battery pack topology structure on state-of-charge estimation accuracy in electric vehicles

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant