CN110346736A - 基于改进小树变换的npc三电平逆变器故障诊断方法 - Google Patents
基于改进小树变换的npc三电平逆变器故障诊断方法 Download PDFInfo
- Publication number
- CN110346736A CN110346736A CN201910749114.9A CN201910749114A CN110346736A CN 110346736 A CN110346736 A CN 110346736A CN 201910749114 A CN201910749114 A CN 201910749114A CN 110346736 A CN110346736 A CN 110346736A
- Authority
- CN
- China
- Prior art keywords
- imf
- component
- signal
- fault
- matrix
- 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.)
- Granted
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/40—Testing power supplies
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Inverter Devices (AREA)
- Tests Of Electronic Circuits (AREA)
Abstract
本发明涉及一种基于改进小树变换的NPC三电平逆变器故障诊断方法,包括:构建NPC三电平逆变器电路仿真模型,模拟故障过程,测量桥臂的电压波形作为故障信号;将故障信号分解为多个IMFs分量;对每个IMFs分量进行希尔伯特Hilbert变换,得到其时频分布和振幅,选择前8个IMFs分量;使用包络分析拟合出包络信号,筛选出故障特征参数;用高斯核函数对小树变换进行优化,产生彼此独立的特征向量样本;将样本数据按3∶7分为训练集和测试集;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障。本发明采用EEMD对故障信号进行分解,再结合希尔伯特变换,收集了更多的频域特征,更适合处理三电平逆变器电路故障产生的非线性、非平稳信号。
Description
技术领域
本发明涉及电路器件故障诊断技术领域,尤其是基于改进小树变换的NPC三电平逆变器电路故障诊断方法。
背景技术
相对于传统的两电平逆变器,三电平逆变器具有器件所承受的电压应力低、电磁干扰小、输出电压的谐波少和逆变效率高等优点,但三电平电路使用了更多的功率开关器件,导致电路的可靠性降低,任何一个器件故障都有可能导致整个电路停止工作,甚至会影响到其他电路的安全,造成不可估量的经济损失。因此如何快速地检测出故障时刻和类型,并准确地对故障信号进行分析,实现故障暂态和扰动时刻的准确定位,是现阶段亟待解决的问题。
从20世纪80年代起,相关专家学者们就已经做了大量的研究工作,提出了多种NPC三电平故障诊断技术。汤清泉,颜世超等在《三电平逆变器的功率管开路故障诊断》(中国电机工程学报,2008,28(21):26-32)中提出了以输出侧PWM电压波形和输出电流的极性作为故障特征的诊断方法,该方法速度快、可靠性高,但需要人工查找故障器件,不能够实现准确定位。廖俊勃,帕孜来·马合木提等在《三电平逆变器IGBT的开路故障诊断研究》(电测与仪表,2015,52(20):35-40)中以测量二极管桥臂上两端电压作为故障信息,根据故障信号和小波函数的特点,利用小波变换将故障信号转化为故障特征向量。小波变换在时域和频域都能得到较高的分辨率,但它的效果取决于小波基函数的选取。此外,一旦确定了基函数,就必须用这一种基函数来处理所有信号,不能自适应。陈丹江,叶银忠在《基于多神经网络的三电平逆变器器件开路故障诊断方法》(电工技术学报,2013,28(6):120-126)中采用傅里叶变换的方法进行故障特征提取,傅里叶变换是一种全域变换,能够在频域内获得较高的分辨率,但却不能给出频率随时间的变化情况,不能满足诊断对时域的要求,使得诊断精度较差。
发明内容
本发明的目的在于提供一种更适合处理三电平逆变器电路故障产生的非线性、非平稳信号,拥有更高的准确率的基于改进小树变换的NPC三电平逆变器电路故障诊断方法。
为实现上述目的,本发明采用了以下技术方案:一种基于改进小树变换的NPC三电平逆变器电路故障诊断方法,该方法包括下列顺序的步骤:
(1)构建NPC三电平逆变器电路仿真模型,模拟故障过程,测量桥臂的电压波形作为故障信号;
(2)用集合经验模态分解EEMD将故障信号分解为多个IMFs分量;
(3)对每个IMFs分量进行希尔伯特Hilbert变换,得到其时频分布和振幅,选择前8个IMFs分量;
(4)对经过步骤(3)所选取的IMFs分量,使用包络分析拟合出包络信号,筛选出故障特征参数;
(5)用高斯核函数对小树变换进行优化,优化后的小树变换可以使步骤(4)中的故障特征参数降维、聚类,产生彼此独立的特征向量样本;
(6)将步骤(5)得到的样本数据按3:7分为训练集和测试集;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障。
所述步骤(2)具体包括以下步骤:
(2a)在原始故障信号u(t)中加入正态分布的白噪声序列a(t),得到U(t):
U(t)=u(t)+a(t)
式中,U(t)是加入白噪声序列后的故障信号,u(t)是原始故障信号,a(t)是白噪声序列;
(2b)对U(t)进行EMD分解,得到一组IMFs分量bj(t)和残留量rj(t):
式中,bj(t)是IMFs分量,rj(t)是残留量,j=1,2,...,n,n为IMFs分量的个数;
(2c)每次加入不同的白噪声序列ai(t),重复步骤(2a)和步骤(2b):
式中,Ui(t)为加入第i个白噪声序列ai(t)后的故障信号,bij(t)为加入第i个白噪声序列ai(t)后的IMFs分量,rij(t)为加入第i个白噪声序列ai(t)后的残留量,N为加入高斯白噪声的次数;
(2d)利用白噪声频谱的均匀分布,当信号加在遍布整个时频空间的白噪声背景上时,不同时间尺度的信号会自动分布到合适的参考尺度上;由于零均值噪声的特性,经过多次平均后,噪声将相互抵消,每次得到的IMFs集成均值就可作为最终结果;IMFs分量平均值表示为:
式中,为第j个IMFs分量的平均值;
最后,原始故障信号表示为:
式中,为原始故障信号,r(t)为EEMD分解后的残余量。
所述步骤(3)具体包括以下步骤:
(3a)对每个IMFs分量做Hilbert变换后,第j个IMFs分量的时频分布为:
式中,为第j个IMFs分量的时频分布,ω为角频率,Aj(t)是第j个IMFs分量的振幅,ωj(t)为第j个IMFs分量的角频率;原始故障信号的时频分布为:
式中,Z(ω,t)为原始故障信号的时频分布;
(3b)选择前8个IMFs分量。
在所述步骤(4)中,对所选取的前8个IMFs分量,使用包络分析得到包络信号:
式中,S(ω)为故障信号的包络信号,Zj(ω,t)为第j个IMFs分量的时频分布;通过观察包络信号频谱中不同频率分量的特征,筛选出故障特征参数,记为{x1,x2,...,xf},f为故障特征参数个数。
所述步骤(5)具体包括以下步骤:
引入高斯核函数代替原始小树变换中的相关系数矩阵,其中,xe,xg为样本,σ为高斯核函数的超参数;
(5a)将故障特征参数{x1,x2,...,xf}输入高斯核函数中,记K(xe,xg)=Keg,f为故障特征参数个数;构造f×f阶高斯核矩阵K,视Keg为K的第e行第g列元素;
(5b)初始化分解层数l=1,...,L,其最大分解层数为L=f-1;在l=0层,初始化基矩阵B0为f×f的单位矩阵,初始化和变量的下标集δ={1,2,...,f};
(5c)在高斯核矩阵K(l-1)找出数值最大的元素,记其所在行为α,所在列为β:
(α,β)=argmaxK(l-1)(e,g)
其中:argmax表示在矩阵中寻找最大值,K(l-1)(e,g)表示在第l-1层的高斯核矩阵K(l-1)在坐标(e,g)处的值,且e<g,e和g必须属于和变量的下标集,差变量不做处理;
(5d)在变量树上进行局部PCA变换,求雅可比旋转矩阵J:
其中:c=cosθl,s=sinθl,旋转角θl可以通过 得到;
式中,表示高斯核矩阵K(l)在第α行第α列处的值,表示高斯核矩阵K(l)在第β行第β列处的值,表示高斯核矩阵K(l)在第α行第β列处的值;
(5e)根据雅可比矩阵更新基矩阵B(l)和高斯核矩阵K(l):
B(l)=B(l-1)J,K(l)=JTK(l-1)J
(5f)进行多尺度分析;
(5g)重复步骤(5c)至步骤(5f),直到l=L层,第l层的小树分解为:
式中,R为第l层的小树分解,φl为尺度函数,为细节函数,φl,d是φl与上一层尺度向量φl-1,d的联合,hl,d是原始数据在尺度向量集上的投影;
(5h)将经过小树变换后得到的特征向量矩阵X′进行归一化处理,记为X:
式中,X(o,p)为矩阵X在第o行第p列的值,X′(o,p)为矩阵X′在第o行第p列的值,min(X′)为矩阵X′中的最小值,max(X′)为矩阵X′中的最大值。
所述步骤(6)具体包括以下步骤:
将步骤(5)得到的样本数据X按3:7分为训练集X1和测试集X2;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障:
(6a)将X1输入支持向量机进行训练,构造最优化模型,目标函数为:
约束条件为:yk(WTx+z)≥1-ξk;
其中,W是特征空间中分类超平面的系数向量,C为惩罚系数,ξk为松弛变量,m为松弛变量个数,yk为样本标签,x为输入样本,z是分类面的阈值;
(6b)结合拉格朗日函数:
式中γk,k=1,2,...,m为拉格朗日乘子,w为特征空间中分类超平面的系数,xk为输入样本;
对w,z,ξ求偏导,得:
其中,γk,γv均表示拉格朗日乘子,yk、yv均表示样本标签;
(6c)将高斯核函数代入上述方程,其中xe,xg为样本,σ为高斯核函数的超参数,最优化问题转化为:
其最小值取决于参数C,σ的选择,使用网格搜索法确定C,σ的取值与步骤(5a)中的相同,故障训练模型建立完毕;
(6d)将测试样本X2输入训练模型,完成诊断与分类。
由上述技术方案可知,本发明的优点在于:第一,采用EEMD对故障信号进行分解,克服了传统方式不能自适应不同信号的缺陷以及EMD造成的模态混叠现象;再结合希尔伯特变换,解决了傅里叶变换不能得到瞬时频率的问题;收集了更多的频域特征,更适合处理三电平逆变器电路故障产生的非线性、非平稳信号。第二,在使用改进的小树变换进行特征提取和维数降低之后,正常和故障的数据集在特征空间中很好地聚类并区分开来,处理时间比未改进前减少,相较于K-Means、DBSCAN、GCD方法在相同情况下拥有更高的准确率;利用改进的小树变换使正常和故障的数据集在特征空间中很好地聚类并区分开来,处理时间缩短,准确率更高。第三,SVM提供了一种避开高维空间,直接用此空间的核函数,再利用线性可分情况下的求解方法直接求解对应的高维空间的决策问题的方法,适合处理非线性函数的拟合,它比BP神经网络具有较好的泛化推广能力。
附图说明
图1是本发明的故障诊断流程图;
图2是NPC三电平逆变器电路拓扑图;
图3是NPC三电平逆变器电路A相拓扑图;
图4是无故障时上桥臂电压波形图;
图5是Qa1故障时电压波形图;
图6是Qa2故障时上桥臂电压波形图;
图7是VDa5故障时上桥臂电压波形图;
图8是Qa1和Qa2故障时上桥臂电压波形图;
图9是Qa1和Qa3故障时上桥臂电压波形图;
图10是Qa1和Qa4开路时上桥臂电压波形图;
图11是Qa2和Qa3故障时上桥臂电压波形图;
图12是Qa2开路时下桥臂电压波形图;
图13是Qa1和Qa2开路时下桥臂电压波形图;
图14是无故障时桥臂电压经验模态分解波形图;
图15是单管开路时桥臂电压经验模态分解波形图;
图16是双管开路时桥臂电压经验模态分解波形图;
图17是四种不同算法的ROC曲线图。
具体实施方式
如图1所示,一种基于改进小树变换的NPC三电平逆变器电路故障诊断方法,该方法包括下列顺序的步骤:
(1)构建NPC三电平逆变器电路仿真模型,模拟故障过程,测量桥臂的电压波形作为故障信号;
(2)用集合经验模态分解EEMD将故障信号分解为多个IMFs分量;
(3)对每个IMFs分量进行希尔伯特Hilbert变换,得到其时频分布和振幅,选择前8个IMFs分量;
(4)对经过步骤(3)所选取的IMFs分量,使用包络分析拟合出包络信号,筛选出故障特征参数;
(5)用高斯核函数对小树变换进行优化,优化后的小树变换可以使步骤(4)中的故障特征参数降维、聚类,产生彼此独立的特征向量样本;
(6)将步骤(5)得到的样本数据按3:7分为训练集和测试集;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障。
所述步骤(2)具体包括以下步骤:
(2a)在原始故障信号u(t)中加入正态分布的白噪声序列a(t),得到U(t):
U(t)=u(t)+a(t)
式中,U(t)是加入白噪声序列后的故障信号,u(t)是原始故障信号,a(t)是白噪声序列;
(2b)对U(t)进行EMD分解,得到一组IMFs分量bj(t)和残留量rj(t):
式中,bj(t)是IMFs分量,rj(t)是残留量,j=1,2,...,n,n为IMFs分量的个数;
(2c)每次加入不同的白噪声序列ai(t),重复步骤(2a)和步骤(2b):
式中,Ui(t)为加入第i个白噪声序列ai(t)后的故障信号,bij(t)为加入第i个白噪声序列ai(t)后的IMFs分量,rij(t)为加入第i个白噪声序列ai(t)后的残留量,N为加入高斯白噪声的次数;
(2d)利用白噪声频谱的均匀分布,当信号加在遍布整个时频空间的白噪声背景上时,不同时间尺度的信号会自动分布到合适的参考尺度上;由于零均值噪声的特性,经过多次平均后,噪声将相互抵消,每次得到的IMFs集成均值就可作为最终结果;IMFs分量平均值表示为:
式中,为第j个IMFs分量的平均值;
最后,原始故障信号表示为:
式中,为原始故障信号,r(t)为EEMD分解后的残余量。
所述步骤(3)具体包括以下步骤:
(3a)对每个IMFs分量做Hilbert变换后,第j个IMFs分量的时频分布为:
式中,为第j个IMFs分量的时频分布,ω为角频率,Aj(t)是第j个IMFs分量的振幅,ωj(t)为第j个IMFs分量的角频率;原始故障信号的时频分布为:
式中,Z(ω,t)为原始故障信号的时频分布;
(3b)选择前8个IMFs分量。
在所述步骤(4)中,对所选取的前8个IMFs分量,使用包络分析得到包络信号:
式中,S(ω)为故障信号的包络信号,Zj(ω,t)为第j个IMFs分量的时频分布;通过判断包络信号频谱中不同频率分量的特征,筛选出故障特征参数,记为{x1,x2,...,xf},f为故障特征参数个数。
所述步骤(5)具体包括以下步骤:
引入高斯核函数代替原始小树变换中的相关系数矩阵,其中,xe,xg为样本,σ为高斯核函数的超参数;
(5a)将故障特征参数{x1,x2,...,xf}输入高斯核函数中,记K(xe,xg)=Keg,f为故障特征参数个数;构造f×f阶高斯核矩阵K,视Keg为K的第e行第g列元素;
(5b)初始化分解层数l=1,...,L,其最大分解层数为L=f-1;在l=0层,初始化基矩阵B0为f×f的单位矩阵,初始化和变量的下标集δ={1,2,...,f};
(5c)在高斯核矩阵K(l-1)找出数值最大的元素,记其所在行为α,所在列为β:
(α,β)=argmaxK(l-1)(e,g)
其中:argmax表示在矩阵中寻找最大值,K(l-1)(e,g)表示在第l-1层的高斯核矩阵K(l-1)在坐标(e,g)处的值,且e<g,e和g必须属于和变量的下标集,差变量不做处理;
(5d)在变量树上进行局部PCA变换,求雅可比旋转矩阵J:
其中:c=cosθl,s=sinθl,旋转角θl可以通过 得到;
式中,表示高斯核矩阵K(l)在第α行第α列处的值,表示高斯核矩阵K(l)在第β行第β列处的值,表示高斯核矩阵K(l)在第α行第β列处的值;
(5e)根据雅可比矩阵更新基矩阵B(l)和高斯核矩阵K(l):
B(l)=B(l-1)J,K(l)=JTK(l-1)J
(5f)进行多尺度分析;
(5g)重复步骤(5c)至步骤(5f),直到l=L层,第l层的小树分解为:
式中,R为第l层的小树分解,φl为尺度函数,为细节函数,φl,d是φl与上一层尺度向量φl-1,d的联合,hl,d是原始数据在尺度向量集上的投影;
(5h)将经过小树变换后得到的特征向量矩阵X′进行归一化处理,记为X:
式中,X(o,p)为矩阵X在第o行第p列的值,X′(o,p)为矩阵X′在第o行第p列的值,min(X′)为矩阵X′中的最小值,max(X′)为矩阵X′中的最大值。
所述步骤(6)具体包括以下步骤:
将步骤(5)得到的样本数据X按3:7分为训练集X1和测试集X2;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障:
(6a)将X1输入支持向量机进行训练,构造最优化模型,目标函数为:
约束条件为:yk(WTx+z)≥1-ξk;
其中,W是特征空间中分类超平面的系数向量,C为惩罚系数,ξk为松弛变量,m为松弛变量个数,yk为样本标签,x为输入样本,z是分类面的阈值;
(6b)结合拉格朗日函数:
式中γk,k=1,2,...,m为拉格朗日乘子,w为特征空间中分类超平面的系数,xk为输入样本;
对w,z,ξ求偏导,得:
其中,γk,γv均表示拉格朗日乘子,yk、yv均表示样本标签;
(6c)将高斯核函数代入上述方程,其中xe,xg为样本,σ为高斯核函数的超参数,最优化问题转化为:
其最小值取决于参数C,σ的选择,使用网格搜索法确定C,σ的取值与步骤(5a)中的相同,故障训练模型建立完毕;
(6d)将测试样本X2输入训练模型,完成诊断与分类。
构建NPC三电平逆变器电路实验模型,主电路如图2所示,电路由三相桥臂构成,共有2个钳位电容、12个IGBT功率开关元件、12个续流二极管和6个中性点钳位二极管。由于该电路是三相对称的,故仅以A相为例分析,B、C相同A相。A相电路的拓扑图如图3所示。本发明只针对一个功率器件发生开路故障和两个功率器件同时发生开路故障的情形进行研究。故障分类如表1所示:
表1故障分类
标签 | 故障类型 |
0 | 无故障 |
1 | Q<sub>al</sub> |
2 | Q<sub>a2</sub> |
3 | Q<sub>a3</sub> |
4 | Q<sub>a4</sub> |
5 | VD<sub>a5</sub> |
6 | VD<sub>a6</sub> |
12 | Q<sub>a1</sub>和Q<sub>a2</sub> |
13 | Q<sub>a1</sub>和Q<sub>a3</sub> |
14 | Q<sub>a1</sub>和Q<sub>a4</sub> |
23 | Q<sub>a2</sub>和Q<sub>a3</sub> |
24 | Q<sub>a2</sub>和Q<sub>a4</sub> |
34 | Q<sub>a3</sub>和Q<sub>a4</sub> |
如图2所示,主电路包括A、B、C三相,三相完全对称,包括12个IGBT功率开关器件(Qa1~Qa4、Qb1~Qb4、Qc1~Qc4)、12个续流二极管(VDa1~VDa4、VDb1~VDb4、VDc1~VDc4)、2个钳位电容(C1、C2)和6个中性点钳位二极管(VDa5~VDa6、VDb5~VDb6、VDc5~VDc6)。输出侧有三种电压:+1/2UD、-1/2UD和0。每个IGBT功率开关器件都与一个续流二极管反并联连接,为电流提供反向导通回路。两个钳位电容串联在直流侧,提供相同的直流电压,C1=C2。三相的钳位二极管通过中点相连,并连接到直流侧中性点,使桥臂上与其相连的点的电位等于中性点电位。
如图3所示,A相由6个开关器件:4个功率开关器件(Qa1~Qa4)、4个续流二极管(VDa1~VDa4)和2个中性点钳位二极管(VDa5~VDa6)组成。Qa1、Qa2、Qa3、Qa4串联,VDa1~VDa4分别反向并联在Qa1~Qa4上。Qa1、Qa2组成上半桥,Qa3、Qa4组成下半桥。在VDa5和VDa6的作用下,桥臂输出P,O,N三种电平。
本发明以上桥臂电压为测量对象,得到各种故障状态下的电压波形如图4至图11所示。由于图6(Qa2开路时的电压波形)与图8(Qa1和Qa2同时开路时的电压波形)相同,故对这两种情况引入新的测量量——下桥臂电压,波形分别如图12、图13所示。对各个故障状态下的电压进行EEMD分解,分解结果如图14-图16所示。绘制接收者操作特征(ROC)曲线来衡量改进后的小树变换与其他方法的优劣,如图17所示,AUC为ROC曲线下的面积,介于0和1之间,值越大说明聚类效果越好。由图可见,改进后的小树变换(Improved-Treelet)优于K-Means、DBSCAN及GCD算法。
如图4所示,桥臂电压波形平稳、规则,呈方波状,无畸变。电压值在0和800V间有序交替。
如图5所示,在正半周期,桥臂电压维持400V不变;在负半周期,桥臂电压先持续维持800V,后在周期交替瞬间产生脉冲电压。
如图6所示,在正半周期,桥臂电压在0和800V之间交替;在负半周期,电压维持800V不变。
如图7所示,在正半周期,桥臂电压维持800V;在负半周期,电压维持400V。
如图8所示,在正半周期,桥臂电压在0和800V之间起伏变化;在负半周期,电压维持800V不变。
如图9所示,在正半周期,桥臂电压维持400V不变;在负半周期,电压在0和800V之间起伏变化。
如图10所示,在正半周期和负半周期稳定时,电压维持400V不变;在周期交替瞬间,桥臂产生脉冲电压。
如图11所示,桥臂电压波形为不规则的矩形波,在正半周期和负半周期分别有一次突变。
如图12所示,桥臂电压波形为受到谐波信号干扰的不规则的锯齿波。
如图13所示,在正半周期和负半周期,桥臂电压值呈稳定上升趋势;正半周期的谐波干扰较负半周期多;在周期交替瞬间,有脉冲电压产生。
对比图14-图16可以得出,在不同故障情况下,同一级的IMF分量的频率和幅值都不同,这表明使用经验模态分解之后,故障信息可以被很好的区分、表现出来。
如图17所示,该曲线通常用来衡量分类的结果好坏。AUC表示ROC曲线与横坐标0和1围成的面积,取值范围是(0,1),越接近1说明分类效果越好。由图17可见,改进小树变换(Improved-Treelet)的AUC值为0.958,K均值聚类算法(K-Means)的AUC值为0.807,具有噪声的基于密度的聚类方法(DBSCAN)的AUC值为0.733,图团体检测方法(GCD)的AUC值为0.653。因为0.958﹥0.807﹥0.733﹥0.653,所以改进小树变换的分类效果最好。
表2为不同诊断方法的诊断结果性能评估表,由表可见,本发明提出的方法诊断准确率高达98.59%,优于BP神经网络、极限学***均训练时间和平均测试时间也少于其他三种方法,运算效率有较大提升。
表2不同方法的诊断性能
综上所述,本发明采用EEMD对故障信号进行分解,克服了传统方式不能自适应不同信号的缺陷以及EMD造成的模态混叠现象;再结合希尔伯特变换,解决了傅里叶变换不能得到瞬时频率的问题;收集了更多的频域特征,更适合处理三电平逆变器电路故障产生的非线性、非平稳信号;利用改进的小树变换使正常和故障的数据集在特征空间中很好地聚类并区分开来,处理时间缩短,准确率更高。
Claims (6)
1.一种基于改进小树变换的NPC三电平逆变器故障诊断方法,其特征在于:该方法包括下列顺序的步骤:
(1)构建NPC三电平逆变器电路仿真模型,模拟故障过程,测量桥臂的电压波形作为故障信号;
(2)用集合经验模态分解EEMD将故障信号分解为多个IMFs分量;
(3)对每个IMFs分量进行希尔伯特Hilbert变换,得到其时频分布和振幅,选择前8个IMFs分量;
(4)对经过步骤(3)所选取的IMFs分量,使用包络分析拟合出包络信号,筛选出故障特征参数;
(5)用高斯核函数对小树变换进行优化,优化后的小树变换可以使步骤(4)中的故障特征参数降维、聚类,产生彼此独立的特征向量样本;
(6)将步骤(5)得到的样本数据按3∶7分为训练集和测试集;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障。
2.根据权利要求1所述的基于改进小树变换的NPC三电平逆变器故障诊断方法,其特征在于:所述步骤(2)具体包括以下步骤:
(2a)在原始故障信号u(t)中加入正态分布的白噪声序列a(t),得到U(t):
U(t)=u(t)+a(t)
式中,U(t)是加入白噪声序列后的故障信号,u(t)是原始故障信号,a(t)是白噪声序列;
(2b)对U(t)进行EMD分解,得到一组IMFs分量bj(t)和残留量rj(t):
式中,bj(t)是IMFs分量,rj(t)是残留量,j=1,2,...,n,n为IMFs分量的个数;
(2c)每次加入不同的白噪声序列ai(t),重复步骤(2a)和步骤(2b):
式中,Ui(t)为加入第i个白噪声序列ai(t)后的故障信号,bij(t)为加入第i个白噪声序列ai(t)后的IMFs分量,rij(t)为加入第i个白噪声序列ai(t)后的残留量,N为加入高斯白噪声的次数;
(2d)利用白噪声频谱的均匀分布,当信号加在遍布整个时频空间的白噪声背景上时,不同时间尺度的信号会自动分布到合适的参考尺度上;由于零均值噪声的特性,经过多次平均后,噪声将相互抵消,每次得到的IMFs集成均值就可作为最终结果;IMFs分量平均值表示为:
式中,为第j个IMFs分量的平均值;
最后,原始故障信号表示为:
式中,为原始故障信号,r(t)为EEMD分解后的残余量。
3.根据权利要求1所述的基于改进小树变换的NPC三电平逆变器故障诊断方法,其特征在于:所述步骤(3)具体包括以下步骤:
(3a)对每个IMFs分量做Hilbert变换后,第j个IMFs分量的时频分布为:
式中,为第j个IMFs分量的时频分布,ω为角频率,Aj(t)是第j个IMFs分量的振幅,ωj(t)为第j个IMFs分量的角频率;原始故障信号的时频分布为:
式中,Z(ω,t)为原始故障信号的时频分布;
(3b)选择前8个IMFs分量。
4.根据权利要求1所述的基于改进小树变换的NPC三电平逆变器故障诊断方法,其特征在于:在所述步骤(4)中,对所选取的前8个IMFs分量,使用包络分析得到包络信号:
式中,S(ω)为故障信号的包络信号,Zj(ω,t)为第j个IMFs分量的时频分布;通过观察包络信号频谱中不同频率分量的特征,筛选出故障特征参数,记为{x1,x2,...,xf},f为故障特征参数个数。
5.根据权利要求1所述的基于改进小树变换的NPC三电平逆变器故障诊断方法,其特征在于:所述步骤(5)具体包括以下步骤:
引入高斯核函数代替原始小树变换中的相关系数矩阵,其中,xe,xg为样本,σ为高斯核函数的超参数;
(5a)将故障特征参数{x1,x2,...,xf}输入高斯核函数中,记K(xe,xg)=Keg,f为故障特征参数个数;构造f×f阶高斯核矩阵K,视Keg为K的第e行第g列元素;
(5b)初始化分解层数l=1,...,L,其最大分解层数为L=f-1;在l=0层,初始化基矩阵B0为f×f的单位矩阵,初始化和变量的下标集δ={1,2,...,f};
(5c)在高斯核矩阵K(l-1)找出数值最大的元素,记其所在行为α,所在列为β:
(α,β)=argmaxK(1-1)(e,g)
其中:argmax表示在矩阵中寻找最大值,K(l-1)(e,g)表示在第l-1层的高斯核矩阵K(l-1)在坐标(e,g)处的值,且e<g,e和g必须属于和变量的下标集,差变量不做处理;
(5d)在变量树上进行局部PCA变换,求雅可比旋转矩阵J:
其中:c=cos θl,s=sin θl,旋转角θl可以通过 得到;
式中,表示高斯核矩阵K(l)在第α行第α列处的值,表示高斯核矩阵K(l)在第β行第β列处的值,表示高斯核矩阵K(l)在第α行第β列处的值;
(5e)根据雅可比矩阵更新基矩阵B(l)和高斯核矩阵K(l):
B(l)=B(l-1)J,K(l)=JTK(l-1)J
(5f)进行多尺度分析;
(5g)重复步骤(5c)至步骤(5f),直到l=L层,第l层的小树分解为:
式中,R为第l层的小树分解,φl为尺度函数,为细节函数,φl,d是φl与上一层尺度向量φl-1,d的联合,hl,d是原始数据在尺度向量集上的投影;
(5h)将经过小树变换后得到的特征向量矩阵X′进行归一化处理,记为X:
式中,X(o,p)为矩阵X在第O行第p列的值,X′(o,p)为矩阵X′在第O行第p列的值,min(X′)为矩阵X′中的最小值,max(X′)为矩阵X′中的最大值。
6.根据权利要求1所述的基于改进小树变换的NPC三电平逆变器故障诊断方法,其特征在于:所述步骤(6)具体包括以下步骤:
将步骤(5)得到的样本数据X按3∶7分为训练集X1和测试集X2;训练集用来构造SVM分类器模型,测试集用来实际诊断电路故障:
(6a)将X1输入支持向量机进行训练,构造最优化模型,目标函数为:
约束条件为:yk(WTx+z)≥1-ξk;
其中,W是特征空间中分类超平面的系数向量,C为惩罚系数,ξk为松弛变量,m为松弛变量个数,yk为样本标签,x为输入样本,z是分类面的阈值;
(6b)结合拉格朗日函数:
式中γk,k=1,2,...,m为拉格朗日乘子,w为特征空间中分类超平面的系数,xk为输入样本;
对w,z,ξ求偏导,得:
其中,γk,γv均表示拉格朗日乘子,yk、yv均表示样本标签;
(6c)将高斯核函数代入上述方程,其中xe,xg为样本,σ为高斯核函数的超参数,最优化问题转化为:
其最小值取决于参数C,σ的选择,使用网格搜索法确定C,σ的取值与步骤(5a)中的相同,故障训练模型建立完毕;
(6d)将测试样本X2输入训练模型,完成诊断与分类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910749114.9A CN110346736B (zh) | 2019-08-14 | 2019-08-14 | 基于改进小树变换的npc三电平逆变器故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910749114.9A CN110346736B (zh) | 2019-08-14 | 2019-08-14 | 基于改进小树变换的npc三电平逆变器故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110346736A true CN110346736A (zh) | 2019-10-18 |
CN110346736B CN110346736B (zh) | 2021-07-02 |
Family
ID=68185000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910749114.9A Active CN110346736B (zh) | 2019-08-14 | 2019-08-14 | 基于改进小树变换的npc三电平逆变器故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110346736B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562517A (zh) * | 2020-06-10 | 2020-08-21 | 合肥工业大学 | 一种npc三电平逆变器开关管开路故障诊断方法 |
CN112036301A (zh) * | 2020-08-31 | 2020-12-04 | 中国矿业大学 | 一种基于类内特征迁移学习与多源信息融合的驱动电机故障诊断模型构建方法 |
CN112560328A (zh) * | 2020-11-18 | 2021-03-26 | 电子科技大学 | 基于表面微应变信号的igbt键合引线故障诊断方法 |
CN112597731A (zh) * | 2021-03-04 | 2021-04-02 | 清华四川能源互联网研究院 | 变流器电磁暂态仿真方法、装置、电子设备及存储介质 |
CN112731019A (zh) * | 2020-12-21 | 2021-04-30 | 合肥工业大学 | 一种anpc三电平逆变器故障诊断方法 |
CN113028117A (zh) * | 2021-03-25 | 2021-06-25 | 南京晨光集团有限责任公司 | 基于比例阀轴控器快速响应及其智能控制方法 |
CN113777485A (zh) * | 2021-09-07 | 2021-12-10 | 江南大学 | 基于价值函数误差的永磁同步电机开路故障诊断方法 |
CN117131336A (zh) * | 2023-10-26 | 2023-11-28 | 深圳市欧康精密技术有限公司 | 一种电子设备连接器数据处理方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663686A (zh) * | 2012-03-19 | 2012-09-12 | 西安电子科技大学 | 基于Treelet变换和高斯尺度混合模型的图像去噪方法 |
CN103048137A (zh) * | 2012-12-20 | 2013-04-17 | 北京航空航天大学 | 一种变工况下的滚动轴承故障诊断方法 |
CN104698397A (zh) * | 2015-03-16 | 2015-06-10 | 浙江万里学院 | 一种多电平逆变器故障诊断方法 |
CN105095566A (zh) * | 2015-06-29 | 2015-11-25 | 南京航空航天大学 | 一种基于小波分析和svm的逆变器故障诊断方法 |
US20160011268A1 (en) * | 2014-07-10 | 2016-01-14 | National Cheng Kung University | Motor fault detecting method and motor fault detecting system |
CN106291162A (zh) * | 2016-07-20 | 2017-01-04 | 江南大学 | 一种光伏二极管箝位式三电平逆变器的故障诊断方法 |
CN106443297A (zh) * | 2016-09-07 | 2017-02-22 | 江南大学 | 光伏二极管箝位型三电平逆变器的决策树svm故障诊断方法 |
CN106682303A (zh) * | 2016-12-26 | 2017-05-17 | 江南大学 | 一种基于经验模态分解和决策树rvm的三电平逆变器故障诊断方法 |
CN108229544A (zh) * | 2017-12-22 | 2018-06-29 | 南昌大学 | 一种npc三电平光伏逆变器开路故障诊断方法 |
CN109815934A (zh) * | 2019-02-12 | 2019-05-28 | 西华大学 | 一种应用于逆变器故障的多尺度oglpe特征提取方法 |
-
2019
- 2019-08-14 CN CN201910749114.9A patent/CN110346736B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663686A (zh) * | 2012-03-19 | 2012-09-12 | 西安电子科技大学 | 基于Treelet变换和高斯尺度混合模型的图像去噪方法 |
CN103048137A (zh) * | 2012-12-20 | 2013-04-17 | 北京航空航天大学 | 一种变工况下的滚动轴承故障诊断方法 |
US20160011268A1 (en) * | 2014-07-10 | 2016-01-14 | National Cheng Kung University | Motor fault detecting method and motor fault detecting system |
CN104698397A (zh) * | 2015-03-16 | 2015-06-10 | 浙江万里学院 | 一种多电平逆变器故障诊断方法 |
CN105095566A (zh) * | 2015-06-29 | 2015-11-25 | 南京航空航天大学 | 一种基于小波分析和svm的逆变器故障诊断方法 |
CN106291162A (zh) * | 2016-07-20 | 2017-01-04 | 江南大学 | 一种光伏二极管箝位式三电平逆变器的故障诊断方法 |
CN106443297A (zh) * | 2016-09-07 | 2017-02-22 | 江南大学 | 光伏二极管箝位型三电平逆变器的决策树svm故障诊断方法 |
CN106682303A (zh) * | 2016-12-26 | 2017-05-17 | 江南大学 | 一种基于经验模态分解和决策树rvm的三电平逆变器故障诊断方法 |
CN108229544A (zh) * | 2017-12-22 | 2018-06-29 | 南昌大学 | 一种npc三电平光伏逆变器开路故障诊断方法 |
CN109815934A (zh) * | 2019-02-12 | 2019-05-28 | 西华大学 | 一种应用于逆变器故障的多尺度oglpe特征提取方法 |
Non-Patent Citations (1)
Title |
---|
周进: "基于EEMD样本熵的三电平逆变器故障诊断方法的研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111562517A (zh) * | 2020-06-10 | 2020-08-21 | 合肥工业大学 | 一种npc三电平逆变器开关管开路故障诊断方法 |
CN112036301A (zh) * | 2020-08-31 | 2020-12-04 | 中国矿业大学 | 一种基于类内特征迁移学习与多源信息融合的驱动电机故障诊断模型构建方法 |
CN112036301B (zh) * | 2020-08-31 | 2021-06-22 | 中国矿业大学 | 一种基于类内特征迁移学习与多源信息融合的驱动电机故障诊断模型构建方法 |
CN112560328A (zh) * | 2020-11-18 | 2021-03-26 | 电子科技大学 | 基于表面微应变信号的igbt键合引线故障诊断方法 |
CN112560328B (zh) * | 2020-11-18 | 2022-04-19 | 电子科技大学 | 基于表面微应变信号的igbt键合引线故障诊断方法 |
CN112731019A (zh) * | 2020-12-21 | 2021-04-30 | 合肥工业大学 | 一种anpc三电平逆变器故障诊断方法 |
CN112597731A (zh) * | 2021-03-04 | 2021-04-02 | 清华四川能源互联网研究院 | 变流器电磁暂态仿真方法、装置、电子设备及存储介质 |
CN113028117A (zh) * | 2021-03-25 | 2021-06-25 | 南京晨光集团有限责任公司 | 基于比例阀轴控器快速响应及其智能控制方法 |
CN113777485A (zh) * | 2021-09-07 | 2021-12-10 | 江南大学 | 基于价值函数误差的永磁同步电机开路故障诊断方法 |
CN113777485B (zh) * | 2021-09-07 | 2022-05-17 | 江南大学 | 基于价值函数误差的永磁同步电机开路故障诊断方法 |
CN117131336A (zh) * | 2023-10-26 | 2023-11-28 | 深圳市欧康精密技术有限公司 | 一种电子设备连接器数据处理方法 |
CN117131336B (zh) * | 2023-10-26 | 2024-01-30 | 深圳市欧康精密技术有限公司 | 一种电子设备连接器数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110346736B (zh) | 2021-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110346736A (zh) | 基于改进小树变换的npc三电平逆变器故障诊断方法 | |
US10234495B2 (en) | Decision tree SVM fault diagnosis method of photovoltaic diode-clamped three-level inverter | |
CN112710969B (zh) | 一种单相半桥五电平逆变器开关管开路故障诊断方法 | |
WO2018120077A1 (zh) | 一种基于经验模态分解和决策树rvm的三电平逆变器故障诊断方法 | |
Ahmadipour et al. | Islanding detection method using ridgelet probabilistic neural network in distributed generation | |
CN105095566B (zh) | 一种基于小波分析和svm的逆变器故障诊断方法 | |
CN104698397B (zh) | 一种多电平逆变器故障诊断方法 | |
Saini et al. | Detection and classification of power quality disturbances in wind‐grid integrated system using fast time‐time transform and small residual‐extreme learning machine | |
CN107607303A (zh) | 基于小波包与som网络的高压断路器机械故障识别方法 | |
CN109782091A (zh) | 基于深度小波极限学***逆变器故障诊断方法 | |
CN105807165B (zh) | Mmc换流站运行状态模糊综合评判方法 | |
CN110068776A (zh) | 基于优化支持向量机的三电平逆变器开路故障诊断方法 | |
CN113125992B (zh) | 基于dbn的npc三电平逆变器故障诊断方法及*** | |
Rodriguez et al. | Classification of power quality disturbances using hilbert huang transform and a multilayer perceptron neural network model | |
CN108898182A (zh) | 一种基于核主元分析与支持向量机的mmc故障诊断方法 | |
Sałat et al. | Support vector machine for soft fault location in electrical circuits | |
CN115712871A (zh) | 一种结合重采样和集成学习的电力电子***故障诊断方法 | |
Khan et al. | An efficient islanding classification technique for single phase grid connected photovoltaic system | |
Deng et al. | A high-accuracy-light-AI data-driven diagnosis method for open-circuit faults in single-phase PWM rectifiers | |
Qin et al. | VWM-DCRNN: A method of combining signal processing with deep learning for fault diagnosis in single-phase PWM rectifier | |
Xu et al. | Power quality dectection and classification in active distribution networks based on improved empiricial wavelet transform and disperson entropy | |
CN110968073B (zh) | 一种hvdc***换相失败故障原因双层溯源辨识方法 | |
CN114019310B (zh) | 一种储能变流器的开路故障预判方法、设备及介质 | |
CN106410836B (zh) | Mmc-hvdc暂态扰动交互影响状态的评价方法 | |
CN114814391A (zh) | 充电桩故障识别方法和存储介质 |
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 |