一种基于变分模态分解的短时异常电池识别方法与流程

1.本发明属于电池技术领域，具体涉及一种根据电池放电过程中的电压数据识别故障电池单体的方法。


背景技术：

2.随着时代的发展，能源短缺和环境污染问题日益突出，燃油汽车被电动汽车代替的趋势不可避免。近几年电动汽车市场占有率逐年递增，锂电池作为电动汽车的关键部件，被广泛应用在纯电动车上。
3.锂离子电池系统是由许多的锂离子电池单体组成的一个非常复杂系统，具有高度非线性和温度敏感特性，而且每个电池单体不可能做到完全一致，所以不一致性也是锂离子电池系统的特性之一。这些特性会导致长期使用的锂离子电池系统存在一定的安全隐患。
4.及时发现电动汽车的电池系统故障，把电池系统的健康状态实时反馈给驾驶员，不仅可以提高电动汽车的质量，促进电动汽车的发展，也可以保障驾驶员的人身安全。
5.引入变分模态分解(vmd)算法，利用粒子群优化算法搜索变分模态分解算法的最佳影响参数组合，根据搜索结果设定变分模态分解算法的惩罚参数和分量个数，并从参数优化变分模态分解算法得到的最佳分量中计算特征参数。
6.支持向量机(support vector machine)是一种应用比较广泛算法，无论是在学术研究和工业设计都比较常用，与机器学习算法相比较，支持向量机(svm)在处理非线性问题上的效果会比较突出。
7.变分模态分解(vmd)和支持向量机(support vector machine)分类算法的应用，可改善电池故障诊断精度。


技术实现要素：

8.鉴于此，有必要提供一种能提高电池故障诊断准确性的方法。
9.一种基于变分模态分解的短时异常电池识别方法，其特征在于，包括以下步骤：
10.s1，采集一组包含多个电池单体的电池组放电电压数据f(t)；
11.s2，基于粒子群算法得到最优组合[α,k]，并对f(t)进行变分模态分解，在k个分量中选择熵值最小对应的模态分量y(t)；
[0012]
s3，从模态分量y(t)中求取电池单体信号对应的自相关系数和相邻电池单体信号之间的互相关系数；
[0013]
s4，选取部分自相关系数和互相关系数对支持向量机进行训练，选择另一部分自相关系数和互相关系作为测试样本，建立svm模型，最后利用svm模型实现电池故障诊断。
[0014]
在上述的基于变分模态分解的短时异常电池识别方法，所述步骤s2包括：
[0015]
s2.1，粒子群算法用熵函数的值作为选取[α,k]的指标，当粒子群算法搜寻到最小熵值时寻优结束，粒子群算法的极小熵值对应的[α,k]即为最优组合；
[0016]
s2.2，对电池组放电电压数据f(t)进行变分模态分解；
[0017]
s2.3，构造变分问题，将k个本征模态函数的和等于输入的电池放电数据作为约束条件，寻求k个本征模态函数h
k
(t)；
[0018]
s2.4，确保k个分量为具有中心频率的有限带宽的模态分量，同时各模态的估计带宽之和最小，将乘法算子λ(t)引入lagrange，得到一个无约束的变分问题；
[0019]
s2.5，为了找到扩展lagrange表达式的鞍点，利用交替方向乘子法对分量函数中心频率和乘法算子λ
n 1
进行更新；
[0020]
s2.6，给定精度ε，若满足式(1)，迭代停止，输出k个满足条件的模态分量；
[0021]
s2.7，
[0022]
3.如权利要求1所述的基于变分模态分解的短时异常电池识别方法，其特征在于，所
[0023]
述步骤s3包括：
[0024]
s3.1，从被选模态分量y(t)中，求取电池单体信号对应的自相关系数；
[0025]
s3.2，从被选模态分量y(t)中，求取相邻电池单体信号之间的互相关系数。在上述的基于变分模态分解的短时异常电池识别方法，所述步骤s4包括：
[0026]
s4.1，选取一部分带标签的自相关系数和互相关系数样本，对支持向量机进行训练；
[0027]
s4.2，选择另一部分自相关系数和互相关系数作为测试样本，对支持向量机进行测试，
[0028]
通过反复的训练和测试，建立基于支持向量机的故障诊断模型；
[0029]
s5.3，以自相关系数和互相关系数作为特征参量，送入建立的故障诊断模型，完成故障电池的识别。
[0030]
本发明可以有效地提高电池故障诊断的准确性，对驾驶员人身安全和车辆的行驶安全有了更好的保障。
附图说明
[0031]
图1基于变分模态分解的短时异常电池故障诊断流程图.
具体实施方式
[0032]
下面将对本发明提供的故障电池单体检测方法作进一步的详细说明。
[0033]
本次发明提供的对输入电池信号进行故障检测和诊断的方法，包括以下具体步骤：
[0034]
s1，提供一组包含多个电池单体的电池组的放电电压数据f(t)；
[0035]
s2，粒子群算法预设惩罚因子α、分量个数k的组合；具体包括以下步骤；
[0036]
s2.1，初始化粒子群算法的各项参数并引入熵函数作为适应度函数；
[0037]
s2.2，选择一个t维度的空间，选取x个粒子，其中x个粒子是参数α和k的组合，组合[α,k]作为粒子的位置；
[0038]
s2.3，随机产生一定数量数[α,k]组合作为粒子的初始位置，随机初始化每个粒子的移动速度，第i个粒子在t维空间的位置为l
i
，速度为v
i
；
[0039]
s2.4，将包络熵函数作为粒子群算法适应度函数，零均值信号x(j)(j＝1,2,...,n)的包络熵函数为：
[0040][0041]
s2.5，a(j)归一化后得到u
j
，x(j)经过hilbert解调后得到包络信号a(j)；
[0042]
s3，粒子群预设一个组合[α,k]后，开始对电池放电数据f(t)做vmd运算；具体包括以下步骤；
[0043]
s3.1，构造变分问题，以k个本征模态函数的和等于输入的电池放电数据为约束条件，寻求k个本征模态函数h
k
(t)；
[0044]
s3.2，对分解出来的每个模态函数进行hilbert变换，求出边际谱，得到单侧频谱的表达式：
[0045][0046]
s3.3，将每个模态函数单侧频谱调制到基频带，加入指数项，表达式为：
[0047][0048]
s3.4，确定每个模态分量带宽，将目标问题转变成求解带约束的变分问题：
[0049][0050]
s3.5，式(4)中：h
k
＝{h1,h1,...,h
k
}为分解得到的k个模态函数分量，ω
k
＝{ω1,ω1,
…
,ω
k
}为各个模态函数的中心频率；
[0051]
s4，为求取约束变分模型的最优解，将乘法算子λ(t)引入lagrange，得到一个无约束的变分问题；具体包括以下步骤；
[0052]
s4.1，粒子群算法已经预设惩罚因子α，将乘法算子λ(t)引入lagrange，增广lagrange的表达式：
[0053][0054]
s4.2，式(6)中α也称为平衡约束参数，λ为lagrange乘法算子；
[0055]
s5，为了找到扩展lagrange表达式的鞍点，利用交替方向乘子法对分量函数中心频率和乘法算子λ
n 1
进行更新，具体包括以下步骤；
[0056]
s5.1，其中取值问题可表述为：
[0057][0058]
s5.2，式中ω
k
等同于
[0059]
s5.3，采用等距变换，将取值问题转换到频域，得到待求解的二次优化问题表达式：
[0060][0061]
s5.4，根据同样的过程，将中心频率的取值问题转换到频域，表达式为：
[0062][0063]
s5.5，得到中心频域的更新方法为；
[0064][0065]
s5.6，式中，相当于当前剩余量的维纳滤波；
[0066][0067]
s5.6，初始化{λ1}，并且令n＝0；
[0068]
s5.7，开始循环迭代，根据表达式(9)(10)(11)更新h
k
、ω
k
、λ，直至k＝k；
[0069]
s5.8，k为粒子群算法的预设值；
[0070]
s5.9，满足迭代停止的条件是满足式(1)；最后得到k个分量；
[0071]
s6，粒子群算法预设一次k值，就会对就会通过适应度函数计算一次熵值，直到搜寻到最小熵值对应的k值；具体包括以下步骤；
[0072]
s6.1，粒子群算法选取第一个位置，预设k值，开始对电池放电数据信号f(t)变分模态分解，并记录一个熵值；
[0073]
s6.2，根据式(12)更新粒子的速度和位置；
[0074][0075]
s6.3，式中，ω为惯性权重；d＝1，2，...d；i＝1,2,...m，k为当前的迭代次数，c1和c2为加速度因子，η为介于[0,1]间的随机数；
[0076]
s6.4，将计算得到的k个模态分量解调后的信号处理成一个概率分布序列u
j
，再由式(1)计算熵值；
[0077]
s6.5，直至迭代次数达到最大设定的k值后输出最佳适应度值及粒子的位置；
[0078]
s6.6，找到最佳适应度值及粒子的位置后，粒子群寻优结束，所得的[α,k]即为最优组合。
[0079]
s7，通过计算k个分量的熵值筛选最优分量，取最小熵值对应的模态分量y(t)；
[0080]
s8，得到最优模态分量y(t)，计算出相应的电池单体信号对应的自相关系数zxg
i
和连续两个电池单体信号之间的相关系数dcjxg
i,i 1
；具体包括以下步骤：
[0081]
s8.1，计算时间窗口内第i个电池单体信号的平均值y
i
(t)与第i个电池单体信号的方差fc
i
；
[0082]
s8.2，根据式(13)计算时间窗口内第i个电池单体信号的自相关系数zxg
i
；
[0083][0084]
s8.3，计算时间窗口内第i个和第i 1个电池单体信号的平均值y
i
(t)、y
i 1
(t)，第i个电池单体信号的方差fc
i
和第i 1个电池单体信号方差fc
i 1
；
[0085]
s8.4，根据式(14)计算时间窗口内第i个单体电池和第i个电池单体信号的相关系数dcjxg
i,i 1
；
[0086][0087]
s9，根据所计算出来的特征参数的值，搭建支持向量机模型，对电池进行故障诊断，具体包括以下步骤：
[0088]
s9.1，选取一部分自相关系数zxg
i
和互相关系数dcjxg
i,i 1
作为训练样本，对支持向量机进行训练；
[0089]
s9.2，选取另一部分自相关系数zxg
i
和互相关系数dcjxg
i,i 1
测试样本，对支持向量机进行测试，通过反复的训练和测试，建立基于支持向量机的故障诊断模型；
[0090]
s9.3，svm模型采用的核函数为高斯核函数，表达式为：
[0091]
k(x1,x2)＝exp(
‑
||x1‑
x2||2/2σ2)
ꢀꢀꢀꢀ
(15)
[0092]
s9.4，式(15)中σ为带宽；
[0093]
s9.5，以自相关系数zxg
i
和互相关系数dcjxg
i,i 1
作为特征参量，送入建立的故障诊断模型，完成故障电池的识别。
[0094]
与现有技术相比，本发明所给出的故障电池单体检测方法通过对放电过程中采集到的电压数据进行更加详细的分析，将电池放电信号用vmd算法进行处理，为了避免人为主观因素的干预，利用粒子群算法搜索影响vmd算法的参数组合[α,k]，用最后得到的模态分量计算特征参数，最后搭建svm模型，对电池进行故障诊断。本发明提供的电池单体故障诊断方法对于提高使用电池组产品的安全性有重要作用。
[0095]
另外，电池技术领域的相关技术人员还可以在本发明精神内做其它改变，当然本发明所要求保护范围也包含这些依据本发明精神所做的改变。