一种开关电源剩余使用寿命预测方法及系统与流程

1.本发明设计电源系统寿命预测技术领域，特别涉及一种开关电源剩余使用寿命预测方法及系统。


背景技术：

2.开关电源具有体积小、功耗低、效率高、稳压范围宽等特点。因此，它被广泛应用于车辆、船舶、飞机等行业。在开关电源单元中，由于电压高、温度高、振动大等原因，电容器和场效应晶体管的退化率和故障率较高。对于控制系统来说，如果开关电源单元出现故障是非常危险的，因为它会导致错误控制或控制故障。因此，有必要对开关电源的剩余使用寿命预测进行研究。
3.剩余使用寿命预测方法可分为两大类：基于模型的方法和数据驱动的方法。基于模型的方法通常从对象的失效机理出发，建立性能退化量与应力(高温、高压、强振动等)之间的关系，并获得性能退化的物理模型，从而分析对象的剩余使用寿命。但通常情况下，由于环境和工作条件的多样性，很难建立准确的物理模型。而且对于开关电源，由于大量电子元件之间的复杂相关性，使用基于模型的方法来分析退化更加困难。如何实现开关电源的寿命的预测成为一个亟待解决的技术问题。


技术实现要素：

4.本发明的目的是提供一种开关电源剩余使用寿命预测方法及系统，以基于数据驱动的方法实现开关电源剩余使用寿命预测。
5.为实现上述目的，本发明提供了如下方案：
6.一种开关电源剩余使用寿命预测方法，所述方法包括如下步骤：
7.获取开关电源的输出电压的当前时频谱；
8.根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态；
9.根据开关电源的当前性能退化状态和当前时频谱，采用sir与hsmm结合的算法，对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间；
10.根据当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，确定所述开关电源的剩余寿命。
11.可选的，所述根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态，之前还包括：
12.获取已失效开关电源样本的全寿命周期的输出电压数据；
13.根据已失效开关电源样本的全寿命周期的输出电压数据，利用公式
计算已失效开关电源样本的全寿命周期的波纹电压数据；其中，v
ω
表示已失效开关电源样本的波纹电压，v0表示已失效开关电源样本的输出电压，v
i
表示已失效开关电源样本的输入电压，f
s
表示已失效开关电源样本的开关频率，l是已失效开关电源样本的电感值，c是已失效开关电源样本的电容值，esr是已失效开关电源样本的电容的等效电阻；
14.根据已失效开关电源样本的全寿命周期的波纹电压数据，采用模糊c均值聚类算法，确定开关电源的性能退化状态的数量，及每个所述波纹电压数据对应的性能退化状态；
15.将每个所述波纹电压数据对应性能退化状态作为所述波纹电压数据对应的输出电压数据的标签，分别将每个所述输出电压数据的时频谱数据作为输入，每个所述输出电压数据的标签作为输出构建训练集和测试集；
16.利用训练集和测试集对所述resnet模型进行训练和测试，获得训练好的resnet模型。
17.可选的，所述根据开关电源的当前性能退化状态和当前时频谱，采用sir与hsmm结合的算法，对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，具体包括：
18.根据开关电源的当前性能退化状态，确定hsmm模型的初始性能退化状态概率分布和初始性能退化状态转移概率矩阵；
19.利用所述初始性能退化状态概率分布构建sir算法的初始粒子集；
20.以所述初始性能退化状态转移概率矩阵和初始粒子集为基础，采用sir算法对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间。
21.可选的，所述以所述初始性能退化状态转移概率矩阵和初始粒子集为基础，采用sir算法对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，具体包括：
22.初始化性能退化次数m的数值为0；
23.根据k
‑
1时刻得到的粒子集通过k
‑
1时刻得到的性能退化状态转移概率矩阵a
k
‑1，预测k时刻的粒子集为：其中，表示根据预测得到的概率，表示k时刻的粒子i；
24.根据hsmm模型的观测状态概率矩阵，利用公式和计算k时刻的每个粒子的权值；其中，和分别表示k时刻的粒子i和粒子j的归一化后的权值，i＝1,2,3，
…
，n
k
，j＝1,2,3，
…
，n
k
，n
k
表示k时刻的粒子的数量，
表示k时刻的粒子i的归一化前的权值，y
k
表示已失效开关电源样本的全寿命周期的时频谱数据中对应的k时刻的时频谱，b表示hsmm模型的观测状态概率矩阵，表示在粒子集中y
k
的概率(即权值)；
25.根据k时刻的粒子集和k时刻的每个粒子的权值，计算k时刻性能退化的估计值其中，表示k时刻的粒子i，n
k
表示k时刻的粒子的数量；
26.根据k时刻性能退化的估计值判断开关电源是否退化至下一个性能退化状态，获得第一判断结果；
27.若所述第一判断结果表示是，则将k的数值赋值给性能退化时刻k
m
，令m的数值增加1，并判断r m的数值是否大于性能退化状态数量总值，获得第二判断结果；其中，r表示当前性能退化状态；
28.若所述第一判断结果表示否或所述第二判断结果表示否，则利用公式计算k时刻的累计时间，利用公式a
k
＝p(d
k
) (i
‑
p(d
k
))
·
a0，计算k时刻的性能退化状态转移概率矩阵a
k
；其中，为k
‑
1时刻性能退化状态的估计值，a0为初始性能退化状态转移概率矩阵，d
k
‑1和d
k
分别为k
‑
1时刻和k时刻的累计时间，p(d
k
)为驻留时间为d
k
的概率矩阵，i为单位矩阵；并令k的数值为1，返回步骤“根据k
‑
1时刻得到的粒子集通过k
‑
1时刻得到的性能退化状态转移概率矩阵a
k
‑1，预测k时刻的粒子集为：
29.若所述第二判断结果表示是，则输出当前性能退化状态的剩余时间当前性能退化状态之后的各个性能退化状态的驻留时间
30.可选的，所述根据当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，确定所述开关电源的剩余寿命，具体包括：
31.当前性能退化状态的剩余时间利用公式计算当前性能退化状态的剩余寿命；其中，表示当前性能退化状态的剩余寿命，表示所有已失效开关电源样本在性能退化状态n的驻留时间的均值；表示当前性能退化状态r所占的比重，其中，y1，y2，分别表示已失效开关电源样本的全寿命周期的时频谱数据中对应的1时刻、2时刻和k0时刻的时频谱，λ表示hsmm模型，表示在hsmm模型下，性能退化状态为当前性能退化
状态r时的概率；
32.根据当前性能退化状态之后的各个性能退化状态的驻留时间利用公式计算所述当前性能退化状态之后的各个性能退化状态的剩余寿命；
33.其中，表示当前性能退化状态之后的第m个性能退化状态的剩余寿命，表示当前性能退化状态r之后的第m个性能退化状态所占的比重；
34.计算当前性能退化状态的剩余时间与当前性能退化状态之后的各个性能退化状态的剩余时间的和，作为所述开关电源的剩余寿命。
35.可选的，根据k
‑
1时刻得到的粒子集通过k
‑
1时刻得到的性能退化状态转移概率矩阵a
k
‑1，预测k时刻的粒子集为：之前还包括：
36.利用公式计算k
‑
1时刻得到的粒子集中有效粒子的数量；其中，表示k
‑
1时刻的粒子i的归一化后的权值，n
k
‑1表示k
‑
1时刻得到的粒子集中粒子的数量，r
eff
表示k
‑
1时刻得到的粒子集中有效粒子的数量；
37.根据k
‑
1时刻得到的粒子集中有效粒子的数量，对k
‑
1时刻得到的粒子集进行重采样更新。
38.可选的，以所述初始性能退化状态转移概率矩阵和初始粒子集为基础，采用sir算法对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，之前还包括：
39.通过对已失效开关电源样本的全寿命周期的输出电压数据的时频谱数据及输出电压数据的时频谱数据的标签进行统计学分析，构建hsmm模型的观测状态概率矩阵为：
40.b
gq
＝p(w
q
|s
t
＝g)，1≤g≤n，1≤q≤m
41.其中，b
gq
为观测状态概率矩阵中第(g，q)个元素，表示已失效开关电源样本在第g个性能退化状态时，观测到第q个输出电压数据的时频谱数据w
q
的概率，m为已失效开关电源样本的输出电压数据的时频谱数据的个数，n表示性能退化状态的数量。
42.一种开关电源剩余使用寿命预测系统，所述系统包括：
43.当前时频谱采集模块，用于获取开关电源的输出电压的当前时频谱；
44.当前性能退化状态确定模块，用于根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态；
45.状态估计模块，用于根据开关电源的当前性能退化状态和当前时频谱，采用sir与hsmm结合的算法，对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间；
46.寿命预测模块，用于根据当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，确定所述开关电源的剩余寿命。
47.可选的，所述系统，还包括：
48.数据获取模块，用于获取已失效开关电源样本的全寿命周期的输出电压数据；
49.波纹电压数据计算模块，用于根据已失效开关电源样本的全寿命周期的输出电压数据，利用公式计算已失效开关电源样本的全寿命周期的波纹电压数据；其中，v
ω
表示已失效开关电源样本的波纹电压，v0表示已失效开关电源样本的输出电压，v
i
表示已失效开关电源样本的输入电压，f
s
表示已失效开关电源样本的开关频率，l是已失效开关电源样本的电感值，c是已失效开关电源样本的电容值，esr是已失效开关电源样本的电容的等效电阻；
50.模糊c均值聚类模块，用于根据已失效开关电源样本的全寿命周期的波纹电压数据，采用模糊c均值聚类算法，确定开关电源的性能退化状态的数量，及每个所述波纹电压数据对应的性能退化状态；
51.样本数据集构建模块，用于将每个所述波纹电压数据对应性能退化状态作为所述波纹电压数据对应的输出电压数据的标签，分别将每个所述输出电压数据的时频谱数据作为输入，每个所述输出电压数据的标签作为输出构建训练集和测试集；
52.resnet模型训练模块，用于利用训练集和测试集对所述resnet模型进行训练和测试，获得训练好的resnet模型。
53.可选的，所述状态估计模块，具体包括：
54.hsmm模型初始参数求解子模块，用于根据开关电源的当前性能退化状态，确定hsmm模型的初始性能退化状态概率分布和初始性能退化状态转移概率矩阵；
55.粒子集初始化子模块，用于利用所述初始性能退化状态概率分布构建sir算法的初始粒子集；
56.状态估计模块子模块，用于以所述初始性能退化状态转移概率矩阵和初始粒子集为基础，采用sir算法对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间。
57.根据本发明提供的具体实施例，本发明公开了以下技术效果：
58.本发明所提供的一种开关电源剩余使用寿命预测方法，所述方法包括如下步骤：获取开关电源的输出电压的当前时频谱；根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态；根据开关电源的当前性能退化状态和当前时频谱，采用sir与hsmm结合的算法，对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间；根据当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，确定所述开关电源的剩余寿命。本发明基于开关电源的输出电压的时频谱，对开关电源进行状态识别，并进一步的基于sir与hsmm结合的算法对其性能退化状态进行递推估计，即基于数据驱动的方法实现开关电源剩余使用寿命预测，无需建立开关电源的物理模型，克服了由于开关电源的结构组成的复杂性造成基于物理模型的预测方法中的物理模型建立困难的技术缺陷。
59.本发明还考虑到纹波电压与内部元件的关系，电容、mosfet(金属
‑
氧化物半导体场效应晶体管)等内部元件的退化会反映在纹波电压中。因此，可以选择纹波电压作为开关电源单元的性能指标，进行开关电源性能退化状态的评估和分类。
附图说明
60.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
61.图1为本发明提供的一种开关电源剩余使用寿命预测方法的流程图；
62.图2为本发明提供的一种开关电源剩余使用寿命预测方法的原理图；
63.图3为本发明提供的开关电源性能退化状态划分的结果示意图；
64.图4为本发明提供的基于小波变换提取的纹波电压的特征示意图；
65.图5为本发明提供的resnet模型的结构图；
66.图6为本发明提供的resnet模型训练流程图；
67.图7为本发明提供的resnet模型的识别结果图；图7a为训练得到的resnet模型对训练集中的不同性能退化状态下数据的识别结果图，图7b为训练得到的resnet模型对测试集中的不同性能退化状态下数据的识别结果图；
68.图8为本发明提供的采用sir与hsmm结合的算法进行性能退化状态递推估计的流程图；
69.图9为本发明提供的采用sir与hsmm结合的算法进行性能退化状态递推估计的结果图；
70.图10为本发明提供的基于sir与hsmm结合的算法获得的剩余使用寿命估计结果图。
具体实施方式
71.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
72.本发明的目的是提供一种开关电源剩余使用寿命预测方法及系统，以基于数据驱动的方法实现开关电源剩余使用寿命预测。
73.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
74.数据驱动法是一种更普遍的方法，因为它直接从数据中提取特征来估计剩余使用寿命，而不需要任何参数简化和模型假设。对于开关电源，由于大量电子元件之间的复杂相关性，使用基于模型的方法来分析退化更加困难，所以主要采用数据驱动的方法来研究开关电源单元的退化。考虑到纹波电压与内部元件的关系，电容、mosfet(金属
‑
氧化物半导体场效应晶体管)等内部元件的退化会反映在纹波电压中。因此，可以选择纹波电压作为开关
电源单元的性能指标，采用数据驱动的方法来分析剩余使用寿命。具体方案如下：
75.如图1和2所示，本发明提供一种开关电源剩余使用寿命预测方法，所述开关电源系统主要指电容、mosfet等内部元件和测量输出电压的传感器。
76.所述方法包括如下步骤：
77.步骤101，获取开关电源的输出电压的当前时频谱。
78.步骤102，根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态。
79.步骤102所述根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态，之前还包括：
80.获取已失效开关电源样本的全寿命周期的输出电压数据。根据已失效开关电源样本的全寿命周期的输出电压数据，利用公式计算已失效开关电源样本的全寿命周期的波纹电压数据；其中，v
ω
表示已失效开关电源样本的波纹电压，v0表示已失效开关电源样本的输出电压，v
i
表示已失效开关电源样本的输入电压，f
s
表示已失效开关电源样本的开关频率，l是已失效开关电源样本的电感值，c是已失效开关电源样本的电容值，esr是已失效开关电源样本的电容的等效电阻。根据已失效开关电源样本的全寿命周期的波纹电压数据，采用模糊c均值聚类算法，确定开关电源的性能退化状态的数量，及每个所述波纹电压数据对应的性能退化状态。将划分好的5类性能退化状态的开关电源纹波电压分别采用离散小波变换做特征处理，得到各个性能退化状态的输出可观测信号的时频谱。将每个所述波纹电压数据对应性能退化状态作为所述波纹电压数据对应的输出电压数据的标签，分别将每个所述输出电压数据的时频谱数据作为输入，每个所述输出电压数据的标签作为输出构建训练集和测试集；训练集用于模型训练，测试集用于验证算法的可行性。利用训练集和测试集对所述resnet模型进行训练和测试，获得训练好的resnet模型。
81.其中，据已失效开关电源样本的全寿命周期的波纹电压数据，采用模糊c均值聚类算法，确定开关电源的性能退化状态的数量，及每个所述波纹电压数据对应的性能退化状态，即，获取开关电源系统工作状态下的输出电压数据v0，计算得到纹波电压v
ω
，并将其分为训练集x(t)和测试集y(t)，其中训练集x(t)用于模型训练，测试集y(t)用于验证算法的可行性。采用模糊c均值聚类算法确定开关电源性能退化状态的数量n，具体包括如下步骤：
82.通常，开关电源的性能下降是由内部元件的退化引起的。根据纹波电压计算公式，可以得出纹波电压随电容、电感和mosfet压降的变化而变化的结论。当开关电源性能变差时，纹波电压也会随之变化。因此根据纹波电压的变化对开关电源的性能状态进行分类。
[0083][0084]
其中，v
ω
是纹波电压，v
i
是输入电压，v0是输出电压，f
s
是开关频率，l是电感值，c是
电容值，esr是电容的等效电阻。
[0085]
一般在描述一个元器件的性能情况时，通常分为健康状态、轻度退化状态、中度退化状态和重度退化状态。设备达到重度退化状态后，再经历一段时间便会引发失效。将训练集x(t)＝{x1,x2,x3,...,x
u
}分为n类，记为s＝{s
c
},c＝1,2,...,n。m
c
(c＝1,2,3,...,n)是聚类中心，μ
c
(x
i
)是第i个数据对第c个聚类的隶属度函数，b为设定的模糊度指数。隶属度函数的聚类损失函数j
f
为：
[0086][0087]
令j
f
对m
c
和μ
c
(x
i
)求偏导，并使求导结果为零，得到聚类损失函数的极小值条件为：
[0088][0089][0090]
通过对上两式的不断迭代，直到满足收敛的条件其中γ是迭代步数，ε是一个很小的常数表示误差阈值，得到最优解n＝5。图3表示了开关电源性能退化状态与观测状态的划分。
[0091]
将划分好的5类性能退化状态的开关电源纹波电压s＝{s
c
},c＝1,2,...,5分别采用离散小波变换做特征处理，得到各个性能退化状态的输出可观测信号的时频谱w：
[0092][0093]
其中α是尺度参数，控制小波函数长度的伸缩变化。β是平移量，控制小波函数平移的量级大小。ψ(t)是小波母函数。m为频率范围指数，n为时间步长变化指数。提取的不同性能状态的输出纹波电压特征如图4所示，将一维纹波电压信号转换成时频谱，其中横轴表示时域特征的变化，纵轴表示频率特征的变化，颜色的深浅变化表示幅值能量大小的变化。原始的纹波电压信号并没有明显的变化特征，但从时频谱上可以发现，随着性能的退化，频率特征逐渐发生变化。
[0094]
训练resnet(残差卷积神经网络)模型得到一个开关电源性能退化的可观测状态的识别模型，对时频谱w进行识别分类，实现对可观测状态的识别。
[0095]
对于纹波电压，其特征变化相对较弱，resnet在cnn(卷积神经网络)中加入直接连通的信道，有较高的模型训练速度和精度。建立50层深度的resnet模型，整体结构如图5所示，具体结构组成如表1所示。
[0096]
表1 resnet模型结构组成
[0097][0098]
搭建的resnet模型由输入层、隐藏层、输出层和多个残差学习单元组成，其中隐藏层包括卷积层、池化层、激活层和全连接层。
[0099]
(1)卷积层的输出为：
[0100][0101]
(2)最大池化函数：
[0102][0103]
(3)选择relu作为激活函数：
[0104][0105]
其中a
l
是第l个卷积层的结果，σ(
·
)是激活函数，是第l个池化层的结果，ω是卷积核，b
l
是第l个卷积层的偏差，h是池化尺寸，q是池化层输出的尺寸，
[0106]
(4)残差学习单元的输入s和输出f(s)直连有两种情况，一种是特征图的维度相同，直接将两个通道的输出相加；另一种是特征图的维度不同，需要对输入w做卷积处理，改变特征图维度之后再相加，得到残差学习单元的输出y，如下式所示，其中γ是卷积操作：
[0107]
y＝f(s) s或y＝f(s) γ(s)
[0108]
(5)全连接层的输出值被传送到输出层，以执行状态估计。
[0109]
建立好resnet模型后，对其进行训练。将每个时频谱特征图与其相对应的状态进行标签。用每类性能退化状态的时频谱数据分别训练resnet模型。resnet残差学习网络训练过程如图6所示，模型的准确率随着损失函数的下降逐渐提升，最终模型的准确率稳定在80％以上。训练得到的resnet模型对不同状态之间的识别结果如图7a所示，能够实现对微小变化的观测特征量之间的区分。然后对测试集y(t)的数据进行分类，如图7b所示。结果表明可以对不同的时频谱数据进行观测状态分类。
[0110]
步骤103，根据开关电源的当前性能退化状态和当前时频谱，采用sir与hsmm结合的算法，对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间。
[0111]
利用时频谱求解hsmm各个参数，包括初始状态概率分布、状态转移概率矩阵、观测状态概率矩阵和各状态的驻留时间分布；
[0112]
步骤103采用如下步骤实现：根据初始时刻的状态概率分布生成初始粒子集；根据
指定前一时刻得到的粒子集，通过状态转移概率矩阵得到指定时刻状态的粒子集；根据指定时刻的观测值以及观测状态概率矩阵得到各预测粒子的权值；将计算得到的各粒子的权重值进行归一化处理；根据指定时刻的粒子集和各粒子权重，计算当前状态的估计值；根据各粒子的权重计算有效粒子数，并对粒子集进行重采样更新，作为下一个时刻状态估计的粒子集；根据各状态的驻留时间计算新的状态转移概率矩阵；重复上述过程，实现对当前状态的递推估计。
[0113]
步骤103使用时频谱w作为性能退化的可观测特征，采用sir与hsmm结合的方法确定开关电源系统的性能退化模型，实现对开关电源当前性能退化状态的递推估计。状态估计过程如图8所示，具体步骤为：
[0114]
开关电源系统的状态方程和观测方程分别为：
[0115]
x
k
＝f
k
(w
k
‑1,v
k
‑1)
[0116]
y
k
＝h
k
(w
k
,n
k
)
[0117]
其中w
k
表示时刻k系统的状态，即时频谱。v
k
表示时刻k的随机噪声，f
k
(
·
)表示动态系统的状态方程，与hsmm中的状态转移概率矩阵a相对应，y
k
表示时刻k动态系统的观测值，n
k
表示时刻k的观测随机噪声，h
k
(
·
)表示动态系统从状态到观测的观测方程，与hsmm中的观测状态概率矩阵b相对应。开关电源性能退化状态估计过程如下：
[0118]
(1)利用全寿命周期的时频谱数据w求解hsmm各参数λ＝(π,a0,b,θ)。
[0119]
初始状态概率分布
[0120]
π＝{π
i
},π
i
＝p(s1＝g),1≤g≤n
[0121]
s1是开关电源在初始时刻的性能状态，即健康状态；
[0122]
状态转移概率矩阵，表示开关电源在运行的过程中状态之间发生转移的概率，计算公式为
[0123]
a0＝{a
ij
}，a
ij
＝p(s
t 1
＝j|s
t
＝i)，1≤i,j≤n
[0124]
其中a
ij
表示开关电源在运行的过程中从状态i转移到状态j的概率。
[0125]
观测状态概率矩阵
[0126]
b
gq
＝p(w
q
|s
t
＝g)，1≤g≤n，1≤q≤m
[0127]
m为开关电源可观测状态(即时频谱w)的个数。b
gq
表示性能状态为s
t
＝g时，观测到第q个可观测状态w
q
的概率；
[0128]
各状态的驻留时间分布
[0129]
θ＝{θ1,
…
θ
n
}
[0130]
θ
i
是概率密度函数的参数。
[0131]
(2)生成初始粒子集：根据初始时刻的状态概率分布π生成初始粒子集
[0132]
(3)状态转移(预测)：根据k
‑
1时刻得到的粒子集通过状态转移概率矩阵a
k
‑1得到k时刻状态的粒子集
[0133][0134]
(4)计算粒子权值(更新)：根据步骤三resnet模型识别得到的k时刻的观测值y
k
以
及观测状态概率矩阵b得到各预测粒子的权值
[0135][0136]
(5)权值归一化：将计算得到的各粒子的权重值进行归一化处理：
[0137][0138]
(6)状态估计：根据k时刻的粒子集和各粒子权重计算当前性能退化状态的估计值：
[0139][0140]
(7)重采样：根据各粒子的归一化权值计算有效粒子数，并对粒子集进行重采样更新，作为下一个时刻状态估计的粒子集。有效粒子数r
eff
计算公式如下：
[0141][0142]
(8)状态转移概率矩阵更新：根据开关电源各性能退化状态的驻留时间d
k
计算新的状态转移概率矩阵a
k
：
[0143]
a
k
＝p(d
k
) (i
‑
p(d
k
))
·
a0[0144]
其中，
[0145][0146]
上述过程中，(2)～(7)为sir方法，(1)和(8)为hsmm方法。将这两种方法的计算过程进行结合，并重复该过程，实现对开关电源当前性能退化状态的递推估计，状态序列递推估计结果如图9所示。
[0147]
从结果对比可以看出，通过sir和hsmm结合的方法可以很好地通过前一时刻的状态和当前时刻的观测递推实现对隐藏状态序列的估计，解决了hsmm需要整个观测序列才能估计隐藏状态序列的缺点。
[0148]
步骤104，根据当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，确定所述开关电源的剩余寿命。
[0149]
确定了开关电源当前的性能退化状态后，剩余使用寿命的计算分为两部分：
[0150]
(1)计算开关电源当前状态的剩余时间。
[0151][0152]
其中是hsmm模型λ＝(π,a0,b,θ)中参数θ
i
中的均值，δ
k
(i)表示不同状态所占的比重，通过加权求和得到该状态下开关电源剩余时间的估计。
[0153]
(2)计算开关电源后续状态的剩余时间。
[0154]
根据初始状态转移概率矩阵a0计算下一个状态，直至失效状态。定义开关电源下一个状态可能出现的概率为：
[0155][0156]
其中概率最大的为下一个时刻可能出现的状态：
[0157][0158]
如果达到失效状态，则开关电源会在该状态下达到驻留时间时发生失效。计算得到开关电源各状态下的剩余时间：
[0159][0160]
按上述过程迭代计算，最终可以得到开关电源剩余使用寿命的估计值：
[0161]
rul＝∑d
[0162]
图10为sir与hsmm结合的剩余使用寿命估计结果，能够看出该方法可以很好地实现对开关电源剩余使用寿命的估计。
[0163]
本发明还提供一种开关电源剩余使用寿命预测系统，所述系统包括：
[0164]
当前时频谱采集模块，用于获取开关电源的输出电压的当前时频谱；
[0165]
当前性能退化状态确定模块，用于根据所述当前时频谱，利用训练好的resnet模型，确定所述开关电源的当前性能退化状态；
[0166]
状态估计模块，用于根据开关电源的当前性能退化状态和当前时频谱，采用sir与hsmm结合的算法，对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间。
[0167]
所述状态估计模块，具体包括：hsmm模型初始参数求解子模块，用于根据开关电源的当前性能退化状态，确定hsmm模型的初始性能退化状态概率分布和初始性能退化状态转移概率矩阵；粒子集初始化子模块，用于利用所述初始性能退化状态概率分布构建sir算法的初始粒子集；状态估计模块子模块，用于以所述初始性能退化状态转移概率矩阵和初始粒子集为基础，采用sir算法对开关电源的性能退化状态进行递推估计，确定所述开关电源的当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间。
[0168]
寿命预测模块，用于根据当前性能退化状态的剩余时间和所述当前性能退化状态之后的各个性能退化状态的驻留时间，确定所述开关电源的剩余寿命。
[0169]
作为一种优选的实施方式，所述系统，还包括：数据获取模块，用于获取已失效开关电源样本的全寿命周期的输出电压数据；波纹电压数据计算模块，用于根据已失效开关电源样本的全寿命周期的输出电压数据，利用公式计算已失效开关电源样本的全寿命周期的波纹电压数据；其中，v
ω
表示已失效开关电源样本的波纹电压，v0表示已失效开关电源样本的输出电压，v
i
表示已失效开关电源样本的输入电
压，f
s
表示已失效开关电源样本的开关频率，l是已失效开关电源样本的电感值，c是已失效开关电源样本的电容值，esr是已失效开关电源样本的电容的等效电阻；模糊c均值聚类模块，用于根据已失效开关电源样本的全寿命周期的波纹电压数据，采用模糊c均值聚类算法，确定开关电源的性能退化状态的数量，及每个所述波纹电压数据对应的性能退化状态；样本数据集构建模块，用于将每个所述波纹电压数据对应性能退化状态作为所述波纹电压数据对应的输出电压数据的标签，分别将每个所述输出电压数据的时频谱数据作为输入，每个所述输出电压数据的标签作为输出构建训练集和测试集；resnet模型训练模块，用于利用训练集和测试集对所述resnet模型进行训练和测试，获得训练好的resnet模型。
[0170]
本发明首先获取开关电源系统工作状态下的输出电压数据，计算得到纹波电压，采用模糊c均值聚类算法确定开关电源性能退化状态的数量；将划分好的5类性能退化状态的开关电源纹波电压分别采用离散小波变换做特征处理，得到各个性能状态的输出可观测信号的时频谱；用得到的时频谱作为性能退化的可观测特征，训练resnet模型得到一个开关电源性能退化的可观测状态的识别模型，从而实现对可观测状态的识别；然后采用sir与hsmm结合的方法确定开关电源系统的性能退化模型；依据当前所处状态和各状态驻留时间分布实现开关电源剩余使用寿命预测。本发明能够对开关电源的剩余使用寿命进行较准确的预测。
[0171]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0172]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。