一种基于S变换与谱峭度的局部放电白噪声抑制方法技术

本发明专利技术公开了一种基于S变换与谱峭度的局部放电白噪声抑制方法，步骤如下：输入染噪信号h(j)；对h(j)进行S变换，得到时频矩阵A；对A进行奇异值分解，求取该奇异值的差分谱并初步确定局部放电信号的有效阶次p；取前p个奇异值进行信号重构得到时频矩阵A*并进行S逆变换，得到疑似PD信号y(z)发生的各个时间段范围[t

【技术实现步骤摘要】
一种基于S变换与谱峭度的局部放电白噪声抑制方法


[0001]本专利技术涉及电力
，具体涉及一种基于S变换与谱峭度的局部放电白噪声抑制方法。

技术介绍

[0002]局部放电检测是评判设备的绝缘状态、及时发现设备内部缺陷的重要手段。在局部放电检测现场时，电气设备周围通常存在大量电磁干扰，绝缘缺陷产生的微弱的局部放电信号容易淹没在严重的背景噪声中。因此研究如何有效降噪，对电力设备的状态检测具有重大的意义。
[0003]白噪声是由于电气设备的热噪声引起的，大量存在于局部放电检测现场。在现有技术中，目前局部放电中的白噪声去除主要有小波分析法，经验模态分析法等。小波方法去噪的主要思路是通过采用阈值修正小波系数，再利用修正的小波系数进行信号重构，即可得到去噪后的信号，但其去噪效果受基小波类型、阈值选择方法等影响较大。EMD存在模式混合和端点效应等缺陷。在去除白噪声方面，上述方法仍然需要进行改进。近来利用奇异值分解法用于抑制局部放电噪声得到了越来越多的关注，使用奇异值分解进行去噪的关键在于轨迹矩阵与有效信号的奇异值阶次的确定。有文献利用基于特征分解的奇异值分解来抑制局部放电信号中的窄带干扰和白噪声，取得了较好的效果，但该方法不能自动选取局部放电信号的有效奇异值阶次，并且文中所加的窄带干扰只有一个，不具有一般性。S变换是由Stockwell提出的，它结合了短时傅里叶变换单频率独立分析和小波变换的多分辨分析的能力，是一个可以完全无损地实现逆变换的时频分析工具，通过S变换得到的时频矩阵从时间和频率上分解了信号，含有丰富的时频信息，故可利用S变换得到的时频矩阵作为轨迹矩阵进行奇异值分解，与其他方法结合有效地确定局部放电信号的奇异值的有效阶次，继而完成白噪声的抑制。

技术实现思路

[0004]本专利技术的目的是为了解决现有技术中的上述缺陷，提供一种基于S变换与谱峭度的局部放电白噪声抑制方法，该去噪方法能有效地去除局部放电信号中混杂的白噪声信号。
[0005]本专利技术的目的可以通过采取如下技术方案达到：
[0006]一种基于S变换与谱峭度的局部放电白噪声抑制方法，所述局部放电白噪声抑制方法包括下列步骤：
[0007]S1、输入染噪信号h(j),j＝1,2,...,N，其中j为信号采样点；
[0008]S2、对染噪信号h(j)进行S变换，得到时频矩阵A；
[0009]S3、对时频矩阵A进行奇异值分解，求取该奇异值的差分谱，基于奇异值差分谱初步确定局部放电信号的有效阶次p；
[0010]S4、取前p个奇异值进行信号重构得到时频矩阵A*并进行S逆变换，得到疑似局部
放电(Partial Discharge,PD)PD信号y(z)发生的各个时间段范围[t
l1
,t
l2
]，l＝1,2,...,p，利用谱峭度算法进一步确定信号的真实有效阶次*
[0011]p；
[0012]S5、根据确定信号的真实有效阶次p
*
进行奇异值重构得到时频矩阵C，对时频矩阵C进行S逆变换，得到去噪后的局部放电信号g(j)。
[0013]进一步地，所述步骤S1中，采用单指数衰减振荡函数和双指数衰减振荡函数两种数学模型模拟局部放电信号和服从w～N(0,0.42)的白噪声信号共同构成染噪信号h(j)。
[0014][0015]进一步地，所述步骤S2过程如下：
[0016]对染噪信号h(j)进行S变换，由卷积定理可得S变换离散形式为：
[0017]式中，N为采样总点数，T为采样间间隔，k,m＝0,1,...,N-1；n＝1,...,N-1，是h(j)的快速傅里叶变换，为经过S变换后得到的时频矩阵，记为时频矩阵A。
[0018]进一步地，所述步骤S3过程如下：
[0019]S31、对时频矩阵A进行奇异值分解，奇异值分解公式如下：
[0020][0021]其中，矩阵U、V分别为w
×
w、q
×
q的正交矩阵，D＝diag(λ1,λ2,...,λ
i
,...,λ
r
)，r＝min(w,q)，奇异值λ
i
满足λ1≥λ2≥...≥λi...≥λr＞0，u
i
和v
i
表示w、q维列向量；
[0022]S32、求取该奇异值的差分谱：
[0023]令b
i
＝λ
i-λ
i+1
i＝1,2,......,r-1
[0024]选取b
i
组成奇异值差分谱B＝(b1,b2,...,b
i
,...,b
r-1
)，其中b
i
为相邻奇异值的差值；
[0025]S33、通过奇异值差分谱初步确定疑似局部放电的有效阶次p，具体如下：
[0026]选择差分谱的最大峰值b
p
对应的奇异值序号作为疑似局部放电的有效阶次p，前p个奇异值记为[λ1,λ2,...,λ
p
]。
[0027]进一步地，所述步骤S4过程如下：
[0028]S41、根据步骤S31中的奇异值分解公式，利用前p个奇异值[λ1,λ2,...,λ
p
]进行重构得到时频矩阵A*；
[0029]S42、对时频矩阵A*进行S逆变换，得到疑似局部放电信号y(z)，并确定疑似局部放电y(z)发生的时间段范围[t
l1
,t
l2
]，其中S逆变换公式为：
[0030][0031]式中：N为采样总点数，T为采时间间隔，n,k＝0,1,...,N-1,为重构得
到的时频矩阵A*；
[0032]S43、在染噪信号h(j)上截取时间段[t
l1
,t
l2
]的信号，记为离散信号x(o)，采用短时傅里叶变换STFT计算离散信号x(o)的谱峭度，其中离散信号x(o)的短时傅里叶变换：
[0033][0034]谱峭度公式如下：
[0035][0036]其中，f为频率，w(o)为窗函数，Sn(f)为信号n阶谱矩，E_为均值计算，|.|为取模计算；
[0037]S44、根据白噪声的谱峭度为0及局部放电信号的谱峭度大于0的特点，检查所有的疑似局部放电信号y(z)是否为局部放电信号，若所有的疑似局部放电信号y(z)经过谱峭度分析均判断为局部放电信号，则继续令p＝p+1，返回步骤S41继续执行，直到出现第一个非局部放电信号，转到步骤S45，若不均为PD信号，则p＝p-1，再返回步骤S41执行，直到找到所有的信号均为PD信号，转到步骤S45；
[0038]S45、将当前的有效阶次P*＝P作为最佳奇异值阶次。
[0039]进一步地，所述步骤S5过程如下：
[0040]S51、利用步骤S31中的奇异值分解公式使用前p
*
个奇异值进行重本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于S变换与谱峭度的局部放电白噪声抑制方法，其特征在于，所述局部放电白噪声抑制方法包括下列步骤：S1、输入染噪信号h(j),j＝1,2,...,N，其中j为信号采样点；S2、对染噪信号h(j)进行S变换，得到时频矩阵A；S3、对时频矩阵A进行奇异值分解，求取该奇异值的差分谱，基于奇异值差分谱初步确定局部放电信号的有效阶次p；S4、取前p个奇异值进行信号重构得到时频矩阵A*并进行S逆变换，得到疑似局部放电PD信号y(z)发生的各个时间段范围[t
l1
,t
l2
]，l＝1,2,...,p，利用谱峭度算法进一步确定信号的真实有效阶次p
*
；S5、根据确定信号的真实有效阶次p
*
进行奇异值重构得到时频矩阵C，对时频矩阵C进行S逆变换，得到去噪后的局部放电信号g(j)。2.根据权利要求1所述的一种基于S变换与谱峭度的局部放电白噪声抑制方法，其特征在于，所述步骤S1中，采用单指数衰减振荡函数和双指数衰减振荡函数两种数学模型模拟局部放电信号和服从w～N(0,0.42)的白噪声信号共同构成染噪信号h(j)。3.根据权利要求1所述的一种基于S变换与谱峭度的局部放电白噪声抑制方法，其特征在于，所述步骤S2过程如下：对染噪信号h(j)进行S变换，由卷积定理可得S变换离散形式为：式中，N为采样总点数，T为采样间间隔，k,m＝0,1,...,N-1；n＝1,...,N-1，是h(j)的快速傅里叶变换，为经过S变换后得到的时频矩阵，记为时频矩阵A。4.根据权利要求1所述的一种基于S变换与谱峭度的局部放电白噪声抑制方法，其特征在于，所述步骤S3过程如下：S31、对时频矩阵A进行奇异值分解，奇异值分解公式如下：其中，矩阵U、V分别为w
×
w、q
×
q的正交矩阵，D＝diag(λ1,λ2,...,λ
i
,...,λ
r
)，r＝min(w,q)，奇异值λ
i
满足λ1≥λ2≥...≥λ
i
...≥λ
r
＞0，u
i
和v
i
表示w、q维列向量；S32、求取该奇异值的差分谱：令b
i
＝λ
i-λ
i+1
，i＝1,2,......,r-1选取b

【专利技术属性】
技术研发人员：罗新，艾嘉伟，牛海清，宋廷汉，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：