CN110221116A - 基于加窗插值和解析模式分解的电压闪变包络检测方法 - Google Patents

Abstract

本发明公开了一种基于加窗插值和解析模式分解的电压闪变包络检测方法，它包括：通过模数转换器将被测的电网电压电流信号从模拟信号转化为数字信号，得到N点长度的采样数据；对N点长度的采样数据进行加窗修正处理，得到加窗修正后的离散电压信号y(n)；利用双谱线插值FFT对离散电压信号y(n)进行频谱分析，求取修正后的基波频率和相位，得到电网电压同步信号；利用解析模式分解方法，提取并得到电压闪变包络信号的分量v(n)；对电压闪变包络信号的分量v(n)进行频谱分析，得到电压闪变包络信号的幅值修正函数和频率修正函数；本发明具有准确提取，检测精度高，抗干扰性强，可用于电压闪变的实时检测。

Description

基于加窗插值和解析模式分解的电压闪变包络检测方法

技术领域

本发明涉及电力系统领域，具体涉及一种基于加窗插值和解析模式分解的电压闪变包络检测方法。

背景技术

波动性负荷的应用导致电力系统中电压产生波动，照明设备发生闪变，对工业生产日常生活造成严重影响，是导致供、用电设备故障与失效的重要原因。闪变是当前电能质量的重要参数之一，对电压闪变进行准确测量，可为抑制和消除电压波动和闪变的影响提供科学依据，以便为电力系统采取措施提供预警。对于分离出电压波动信号的方法，现有的大部分检测算法所使用的都是基于IEC提出来的平方解调检波法。这种获取波动信号的方法简单加上其计算量较小，受到了广泛的应用，但是它的缺点是计算精度较低，如果需要提高计算精度则需要在算法中加入校正步骤来实现，这又导致整个算法结构复杂，计算速度慢。

发明内容

本发明要解决的技术问题是：提供一种基于加窗插值和解析模式分解的电压闪变包络检测方法，以解决现有技术对于分离出电压波动信号采用基于IEC提出来的平方解调检波法，存在计算精度较低，如果需要提高计算精度则需要在算法中加入校正步骤来实现，这又导致整个算法结构复杂，计算速度慢等技术问题。

本发明的技术方案是：

一种基于加窗插值和解析模式分解的电压闪变包络检测方法，它包括：

步骤S101，通过模数转换器将被测的电网电压电流信号从模拟信号转化为数字信号，得到N点长度的采样数据；

步骤S102，对N点长度的采样数据采用四项三阶Nuttall窗函数进行加窗修正处理，得到加窗修正后的离散电压信号y(n)；

步骤S103，利用双谱线插值FFT对修正后的离散电压信号y(n)进行频谱分析，求取修正后的基波频率和相位，得到电网电压同步信号；

步骤S104，利用解析模式分解方法，针对建立的电压闪变信号的调制模型，提取并得到电压闪变包络信号的分量v(n)；

步骤S105，利用四项三阶Nuttall窗和双谱线插值FFT对电压闪变包络信号的分量v(n)进行频谱分析，得到电压闪变包络信号的幅值修正函数和频率修正函数。

所述的一种基于加窗插值和解析模式分解的电压闪变包络检测方法，

所述的电压闪变信号的调制模型为：

式中：U₀为基波电压幅值，ω₀、θ₀为基波电压角频率和初相角；A(t)表示幅度值变化的信号；v(t)为闪变包络信号；h为波动分量的项数；m_j、ω_j、分别为第j项波动分量的幅值系数、角频率和相角；离散表达式为

式中：ω₀＝2πf₀/f_s；ω_j＝2πf_j/f_s；f₀为基波频率，f_s为信号采样频率，f_j为第j项闪变波动分量的频率。

所述四项三阶Nuttall窗函数是一种余弦组合窗，其时域表示为

式中，M为窗函数的项数，N为窗的长度，b_m为窗的系数，b₀＝0.338946，b₁＝0.481976，b₂＝0.161054，b₃＝0.018027。

用窗函数w(n)对电压信号进行处理，可得到加窗后信号的离散傅里叶变换的表达式

式中：W(f)是w(n)的连续频谱；离散抽样间隔为Δf＝f_s/N；N为数据截断长度；h为波动分量的项数；f_j为第j项闪变包络分量的频率；m_j为第j项闪变包络分量的幅值。

设k₀为峰值点，k₂和k₁分别为附近最大谱线和次最大的谱线，k₁≤k₀≤k₂＝k₁+1，α＝k₀-k₁-0.5。这两条谱线的幅值分别为y₁＝|X(k₁Δf)|和y₂＝|X(k₂Δf)|，设参数β＝(y₂-y₁)/(y₂+y₁)，α＝k₀-k₁-0.5，可知α取值范围为[0.5,0.5]，

记式(5)的反函数为α＝h^-1(β)，基于4项3阶Nuttall窗双谱线插值FFT求取基波频率f₀和相位的修正表达式为：

f₀＝k₀Δf＝(α+k₁+0.5)Δf (公式6)

所述式中i＝1,2。

所述基于4项3阶Nuttall窗的5阶逼近式为：

α＝H(β)＝2.95494514β+0.17671943β³+0.09230694β⁵ (公式8)

g(α)＝3.20976143+0.9187393α²+0.14734229α⁴ (公式9)

所述解析模式分解模型为：

S₁＝sin(δt)H[cos(δt)G(t)]-cos(δt)H[sin(δt)G(t)] (公式10)

式中：连续的时间信号G(t)可表示为G(t)＝S₁+S₂形式，且频率较低的成分S₁频率和较高的成分S₂频率成分不重叠，δ为截断频率δ，截断频率δ大于频率较低成分S₁的最大频率，且小于频率较高成分S₂的最小频率，H[X]代表对X进行希尔伯特变换。

所述电网电压同步信号为将提取的同步信号与电网信号相乘；再做解析模式分解，从而提取出电压闪变包络信号。

截断频率δ选取为50Hz。

利用双谱线插值FFT对电压闪变包络信号的分量y(n)进行频谱分析得到的电压闪变包络信号的幅值修正函数和频率修正函数为：

f_j＝k_jΔf＝(α+k_j1+0.5)Δf (公式11)

式中：f_j为第j项电压闪变包络信号的频率修正函数，m_j为第j项电压闪变包络信号的幅值修正，为第j项电压闪变包络信号的相位修正函数。

本发明有益效果：

本发明提出了一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，该方法的提取包络过程避免了复杂滤波器的设计，且不受倍频分量的影响，经仿真实验，结果表明该电压闪变检测算法的检测精度受电网基波频率波动、谐波、间谐波、噪声的影响很小，且对含有单一频率和多频率波动分量的闪变信号都可以实现闪变参数的准确提取，检测精度高，抗干扰性强，可用于电压闪变的实时检测；解决了现有技术对于分离出电压波动信号采用基于IEC提出来的平方解调检波法，存在计算精度较低，如果需要提高计算精度则需要在算法中加入校正步骤来实现，这又导致整个算法结构复杂，计算速度慢等技术问题。

附图说明

图1为本发明的流程示意图；

图2为利用本发明对只含有一种频率波动分量时的电压闪变包络信号参数的仿真检测结果示意图；

图3为利用本发明对存在基波频率偏移时的电压闪变包络信号参数的仿真检测结果示意图；

图4为利用本发明对没有谐波、叠加3次和5次谐波、叠加间谐波、叠加整数次谐波同时叠加间谐波四中情况的电压闪变包络信号参数的仿真检测结果示意图；

图5为利用本发明对含噪声的电压闪变包络信号参数的仿真检测结果示意图。

具体实施方式

一种基于加窗插值和解析模式分解的电压闪变包络检测方法，它包括：

步骤S101，通过模数转换器将被测的电网电压电流信号从模拟信号转化为数字信号，得到N点长度的采样数据；

步骤S102，对N点长度的采样数据采用四项三阶Nuttall窗函数进行加窗修正处理，得到加窗修正后的离散电压信号y(n)；

步骤S103，利用双谱线插值FFT对修正后的离散电压信号y(n)进行频谱分析，求取修正后的基波频率和相位，得到电网电压同步信号；

步骤S104，利用解析模式分解方法，针对建立的电压闪变信号的调制模型，提取并得到电压闪变包络信号的分量v(n)；

步骤S105，利用四项三阶Nuttall窗和双谱线插值FFT对电压闪变包络信号的分量v(n)进行频谱分析，得到电压闪变包络信号的幅值修正函数和频率修正函数。

所述的电压闪变信号的调制模型为：

式中：U₀为基波电压幅值，ω₀、θ₀为基波电压角频率和初相角；A(t)表示幅度值变化的信号；v(t)为闪变包络信号；h为波动分量的项数；m_j、ω_j、分别为第j项波动分量的幅值系数、角频率和相角；离散表达式为

式中：ω₀＝2πf₀/f_s；ω_j＝2πf_j/f_s；f₀为基波频率，f_s为信号采样频率，f_j为第j项闪变波动分量的频率。

所述四项三阶Nuttall窗函数是一种余弦组合窗，其时域表示为

式中，M＝4为窗函数的项数，N为窗的长度，b_m为窗的系数，b₀＝0.338946，b₁＝0.481976，b₂＝0.161054，b₃＝0.018027。

用窗函数w(n)对电压信号进行处理，可得到加窗后信号的离散傅里叶变换的表达式

式中：W(f)是w(n)的连续频谱；离散抽样间隔为Δf＝f_s/N；N为数据截断长度；h为波动分量的项数；f_j为第j项闪变包络分量的频率；m_i为第j项闪变包络分量的幅值。

设k₀为峰值点，k₂和k₁分别为附近最大谱线和次最大的谱线，k₁≤k₀≤k₂＝k₁+1，α＝k₀-k₁-0.5。这两条谱线的幅值分别为y₁＝|X(k₁Δf)|和y₂＝|X(k₂Δf)|，设参数β＝(y₂-y₁)/(y₂+y₁)，α＝k₀-k₁-0.5，可知α取值范围为[0.5,0.5]，

记式(5)的反函数为α＝h^-1(β)，基于4项3阶Nuttall窗双谱线插值FFT求取基波频率f₀和相位的修正表达式为：

f₀＝k₀Δf＝(α+k₁+0.5)Δf (公式6)

所述式中i＝1,2。

所述基于4项3阶Nuttall窗的5阶逼近式为：

α＝H(β)＝2.95494514β+0.17671943β³+0.09230694β⁵ (公式8)

g(α)＝3.20976143+0.9187393α²+0.14734229α⁴ (公式9)

所述解析模式分解模型为：

S₁＝sin(δt)H[cos(δt)G(t)]-cos(δt)H[sin(δt)G(t)] (公式10)

式中：连续的时间信号G(t)可表示为G(t)＝S₁+S₂形式，且频率较低的成分S₁频率和较高的成分S₂频率成分不重叠，δ为截断频率δ，截断频率δ大于频率较低成分S₁的最大频率，且小于频率较高成分S₂的最小频率，H[X]代表对X进行希尔伯特变换。

所述电网电压同步信号为将提取的同步信号与电网信号相乘；再做解析模式分解，从而提取出电压闪变包络信号。

截断频率δ选取为50Hz。

利用双谱线插值FFT对电压闪变包络信号的分量y(n)进行频谱分析得到的电压闪变包络信号的幅值修正函数和频率修正函数为：

f_j＝k_jΔf＝(α+k_j1+0.5)Δf (公式11)

式中：f_j为第j项电压闪变包络信号的频率修正函数，m_j为第j项电压闪变包络信号的幅值修正，为第j项电压闪变包络信号的相位修正函数。

以下将结合图1-图5对本发明的技术方案做进一步详细说明。

步骤S101，通过模数转换器将被测的电网电压电流信号从模拟信号转化为数字信号，得到N点长度的采样数据；

主要参数选择如下：

(1)采样速率：f_s＝800Hz；

(2)电网基波频率：f₀＝50Hz；

(3)电网基波幅值：

(4)采样数据长度：N＝4096；

步骤S102，以电压信号为例，对N点长度的采样数据采用四项三阶Nuttall窗函数进行加窗修正处理，得到加窗修正后的离散电压信号y(n)；

电压闪变信号的调制模型为：

式中：U₀为基波电压幅值，ω₀、θ₀为基波电压角频率和初相角；A(t)表示幅度值变化的信号；v(t)为闪变包络信号；h为波动分量的项数；m_j、ω_j、分别为第j项波动分量的幅值系数、角频率和相角。其离散表达式为

式中：ω₀＝2πf₀/f_s；ω_j＝2πf_j/f_s；f₀为基波频率，f_s为信号采样频率，f_j为第j项闪变波动分量的频率。

本发明选取的四项三阶Nuttall窗函数是一种余弦组合窗，其时域表示为

式中，M＝4为窗函数的项数，N为窗的长度，b_m为窗的系数，b₀＝0.338946，b₁＝0.481976，b₂＝0.161054，b₃＝0.018027。

用窗函数w(n)对电压信号进行处理，可得到加窗后信号的离散傅里叶变换的表达式

式中：W(f)是w(n)的连续频谱；离散抽样间隔为Δf＝f_s/N；N为数据截断长度；h为波动分量的项数；f_j为第j项闪变包络分量的频率；m_i为第j项闪变包络分量的幅值。

步骤S103，利用双谱线插值FFT对修正后的离散电压信号y(n)进行频谱分析，求取修正后的基波频率和相位，得到电网电压同步信号。

设峰值点k₀附近幅值最大和次最大的谱线分别为k₁和k₂，显然，k₁≤k₀≤k_j2＝k_j1+1，这两条谱线的幅值分别为y₁＝|X(k₁Δf)|和y₂＝|X(k₂Δf)|，设β＝(y₂-y₁)/(y₂+y₁)，α＝k₀-k₁-0.5，可知α取值范围为[0.5，0.5]，由式(4)可得

记式(5)的反函数为α＝h^-1(β)，由β可求参数α。

因此，基波频率修正式为：

f₀＝k₀Δf＝(α+k₁+0.5)Δf (6)

基波幅值修正是对k_j1和k_j2两根谱线进行加权平均，其计算式为

电网信号基波相位修正式为：

调用Matlab程序拟合可得基于4项3阶Nuttall窗的5阶逼近式为：

α＝H(β)＝2.95494514β+0.17671943β³+0.09230694β⁵ (9)

g(α)＝3.20976143+0.9187393α²+0.14734229α⁴ (10)

经该步骤可以求取电网信号基波频率与相位，从而得到电网的同步信号。

步骤S104，利用解析模式分解方法，针对建立的电压闪变信号的调制模型，提取并得到电压闪变包络信号的分量v(n)。

设电网同步信号用u_a(t)表示，即u_a(t)＝cos(ω₀t+θ₀)，引入参数g(t)，令g(t)＝u_a(t)u(t)则

式中：第一项是直流分量，第二项是闪变包络，第三项是二倍工频分量，第四项为二倍工频分量经调制的高频分量。滤除直流分量后，积化和差后表达式为：

式中：第一项是0～35Hz闪变包络分量，第二项为100Hz二倍频分量，第三项为100～135Hz的高频分量，第四项为65～100Hz高频分量，其每一项都可以写为一个系数与余弦函数相乘的形式，展开后总共有3k+1项即可变换为如下形式：

式中：b_j为第j项的系数ω_j为第j项的角速度θ_j为第j项的初相。

为得到闪变包络分量，需要将大于35Hz的高频分量去除得到0～35Hz的包络分量。对g₁(t)作解析模式分解，并设ω_c＝50Hz可得

设δ_c＝ω_ct，δ_j＝ω_jt+θ_j，当ω_j>ω_c时，利用希尔伯特性质H[cosδ]＝sinδ，H[sinδ]＝-cosδ,以及三角变换可得

即大于50Hz的高频分量被滤除。

当ω_j<ω_c时，可得

当ω_j<ω_c时，0～35Hz闪变包络信号保留下来，g₂(t)即是要求电压闪变包络分量v(t)。

步骤S105，利用四项三阶Nuttall窗和双谱线插值FFT对电压闪变包络信号的分量v(n)进行频谱分析，得到电压闪变包络信号的幅值修正函数和频率修正函数。

设峰值点k_j附近幅值最大和次最大的谱线分别为k_j1和k_j2，显然，k_j1≤k_j≤k_j2＝k_j1+1，这两条谱线的幅值分别为y₁＝|X(k_j1Δf)|和y₂＝|X(k_j2Δf)|，设β＝(y₂-y₁)/(y₂+y₁)，α＝k_j-k_j1-0.5，可知α取值范围为[0.5,0.5]，

推导原理同步骤S103，则第j项分量的频率修正式为

f_j＝k_jΔf＝(α+k_j1+0.5)Δf (17)

第j项分量的幅值修正是对k_j1和k_j2两根谱线进行加权平均，其计算式为

当N值较大时，式(18)可简化为

A_j＝N^-1(y₁+y₂)g(α) (19)

第j项分量的初相位修正式为

上述实施例中一种基于加窗插值和解析模式分解的电压闪变包络检测新方法是以四项三阶Nuttall窗和双谱线插值FFT为例进行说明的，但是本发明并不局限于此，其还可以是其他优化窗函数和其他多谱线插值FFT相结合。

本发明的效果通过如下仿真实验可以获得，具体如下：

为了对本文提出的算法的可行性进行验证，利用MATLAB软件进行仿真实验分析。

(1)单一频率的电压闪变包络信号

首先对只含有一种频率波动分量的闪变信号进行仿真研究。电网频率f₀为50Hz,基波的幅值U₀为采样频率为500Hz，闪变波动分量的幅值系数m_i在0～0.1pu的范围内变化，闪变波动分量的频率f_i在0.05～35Hz的范围内变化，基波电压的初相角设为任意值，在以上条件下，以本发明提出的方法进行仿真实验得到的闪变波动分量幅值相对误差如图(2)所示。

由图(2)仿真结果可知，当闪变调制信号只有一种频率的波动分量时，本文提出方法测量的闪变调幅波幅值的相对误差小于5×10^-4％，因此在闪变信号只含有单一频率波动分量情况下本文提出的闪变检测算法可以有效实现闪变包络参数的准确提取。

(2)多频率调制电压闪变

为检测本文提出的方法对含有多个频率波动分量的闪变信号电压闪变参数提取效果，依据式(1)设定闪变信号v(t)中含有三个频率分别为5、11、18Hz的信号进行叠加，并设置相应的调幅波系数分别为0.04(pu)、0.07(pu)、0.09(pu)，其余参数同上，仿真结果如表(1)所示。

表1多调制频率下调幅波系数检测误差

通过表(1)分析可知，三个频率的幅值相对误差都小于2×10^-3％。因此在闪变信号中同时含有多个不同频率波动分量情况下，本发明提出的闪变检测算法可以有效的实现对闪变包络参数的准确提取。

(3)系统频率偏差的影响

实际电网中，由于用电负荷和发电功率的不平衡，造成了电网频率波动，为了检测本文提出的闪变检测算法在电网频率波动时的检测效果，设定电网频率在49.5～50.5Hz的范围内变化，单一频率的调幅波系数为0.1，剩余参数设置同上，进行仿真实验，得到的闪变包络信号的调幅波幅值相对误差的仿真结果如图(3)所示。

由图(3)分析可知电网基波频率在49.5～50.5Hz波动时，本发明提出的算法测得的调幅波幅值的相对误差小于8×10^-4％，因此本文提出的方法可以在电网频率波动的情况下实现对闪变参数的准确提取。

(4)谐波和间谐波的影响

由于电网中存在大量的非线性负载，会造成电网中大量的谐波和间谐波的存在，为了验证本文的方法在谐波和间谐波存在情况下电压闪变的检测效果，分别设置没有谐波、叠加3次和5次谐波、叠加间谐波、叠加整数次谐波同时叠加间谐波四中情况进行仿真实验，参数设置如下，3次谐波幅值U₃＝5％；五次谐波幅值U₅＝10％；间谐波频率f_ih＝168Hz；幅值U_ih＝2％。其余仿真参数同上，仿真结果如图(4)所示。

由图(4)分析可知，在有谐波和间谐波影响的情况下，电压闪变参数相对误差小与4×10^-4％，可知本发明的闪变检测方法可以在谐波和间谐波的影响下实现对闪变参数的高精度检测。

(5)白噪声影响分析

为了验证本文提出闪变检测算法在含有噪声情况下的闪变检测的效果，对含有单一频率波动分量的闪变信号中加入白噪声的情况进行仿真实验。设定白噪声的信噪比在20～100dB的范围内变化，单一频率的波动分量的幅值为0.1pu，频率在0～35Hz的范围内变化，其余参数设置同上，仿真结果如图(5)所示。

由图(5)分析可知：当信噪比从20dB上升到40dB过程，调幅波幅值提取误差为从10％降低到0.5％之后随信噪比的增大，相对误差均小于0.5％。电网电压中存在噪声一般大于40dB，因此，采用本文提出的闪变参数检测方法在噪声情况下可实现对闪变参数的准确检测。

本发明提出了种基于加窗插值和解析模式分解的电压闪变包络检测新方法，该方法的提取包络过程避免了复杂滤波器的设计，且不受倍频分量的影响，经仿真实验，结果表明该电压闪变检测算法的检测精度受电网基波频率波动、谐波、间谐波、噪声的影响很小，且对含有单一频率和多频率波动分量的闪变信号都可以实现闪变参数的准确提取，检测精度高，抗干扰性强，可用于电压闪变的实时检测。

Claims (9)

1.一种基于加窗插值和解析模式分解的电压闪变包络检测方法，它包括：

步骤S101，通过模数转换器将被测的电网电压电流信号从模拟信号转化为数字信号，得到N点长度的采样数据；

步骤S102，对N点长度的采样数据采用四项三阶Nuttall窗函数进行加窗修正处理，得到加窗修正后的离散电压信号y(n)；

步骤S103，利用双谱线插值FFT对修正后的离散电压信号y(n)进行频谱分析，求取修正后的基波频率和相位，得到电网电压同步信号；

步骤S104，建立电压闪变信号的调制模型，利用解析模式分解方法，提取并得到电压闪变包络信号的分量v(n)；

步骤S105，利用四项三阶Nuttall窗和双谱线插值FFT对电压闪变包络信号的分量v(n)进行频谱分析，得到电压闪变包络信号的幅值修正函数和频率修正函数。

2.根据权利要求1所述的一种基于加窗插值和解析模式分解的电压闪变包络检测方法，其特征在于：

所述的电压闪变信号的调制模型为：

式中：U₀为基波电压幅值，ω₀、θ₀为基波电压角频率和初相角；A(t)表示幅度值变化的信号；v(t)为闪变包络信号；h为波动分量的项数；m_j、ω_j、分别为第j项闪变波动分量的幅值系数、角频率和相角；离散表达式为

式中：ω₀＝2πf₀/f_s；ω_j＝2πf_j/f_s；f₀为基波频率，f_s为信号采样频率，f_j为第j项闪变波动分量的频率。

3.根据权利要求1所述的一种基于加窗插值和解析模式分解的电压闪变包络检测方法，其特征在于：

所述四项三阶Nuttall窗函数是一种余弦组合窗，其时域表示为

式中，M＝4为窗函数的项数，N为窗的长度，b_m为窗的系数，b0＝0.338946，b₁＝0.481976，b₂＝0.161054，b₃＝0.018027。

用窗函数w(n)对电压信号进行处理，可得到加窗后信号的离散傅里叶变换的表达式

式中：W(f)是w(n)的连续频谱；离散抽样间隔为Δf＝f_s/N；N为数据截断长度；h为波动分量的项数；f_j为第j项闪变包络分量的频率；m_j为第j项闪变包络分量的幅值。

4.根据权利要求4所述的一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，其特征在于：设k₀为峰值点，k₂和k₁为附近幅值最大谱线和次最大的谱线，k₁≤k₀≤k₂＝k₁+1，α＝k₀-k₁-0.5。这两条谱线的幅值分别为y₁＝|X(k₁Δf)|和y₂＝|X(k₂Δf)|，设参数β＝(y₂-y₁)/(y₂+y₁)，α＝k₀-k₁-0.5，可知α取值范围为[0.5，0.5]，

记式(5)的反函数为α＝h^-1(β)，基于4项3阶Nuttall窗双谱线插值FFT求取基波频率f₀和相位的修正表达式为：

f₀＝k₀Δf＝(α+k₁+0.5)Δf (公式6)

所述式中i＝1,2。

5.根据权利要求4所述的一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，其特征在于，所述基于4项3阶Nuttall窗的5阶逼近式为：

α＝H(β)＝2.95494514β+0.17671943β³+0.09230694β⁵ (公式8)

g(α)＝3.20976143+0.9187393α²+0.14734229α⁴ (公式9)。

6.根据权利要求4所述的一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，其特征在于：

所述解析模式分解模型为：

S₁＝sin(δt)H[cos(δt)G(t)]-cos(δt)H[sin(δt)G(t)] (公式10)

式中：连续的时间信号G(t)可表示为G(t)＝S₁+S₂形式，且频率较低的成分S₁频率和较高的成分S₂频率成分不重叠，δ为截断频率δ，截断频率δ大于频率较低成分S₁的最大频率，且小于频率较高成分S₂的最小频率，H[X]代表对X进行希尔伯特变换。

7.根据权利要求1所述的一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，其特征在于：所述电网电压同步信号为将提取的同步信号与电网信号相乘；再做解析模式分解，从而提取出电压闪变包络信号。

8.根据权利要求6所述的一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，其特征在于，截断频率δ选取为50 Hz。

9.根据权利要求1所述的一种基于加窗插值和解析模式分解的电压闪变包络检测新方法，其特征在于，利用双谱线插值FFT对电压闪变包络信号的分量y(n)进行频谱分析得到的电压闪变包络信号的幅值修正函数和频率修正函数为：

f_j＝k_jΔf＝(α+k_j1+0.5)Δf (公式11)

式中：f_j为第j项电压闪变包络信号的频率修正函数，m_j为第j项电压闪变包络信号的幅值修正，为第j项电压闪变包络信号的相位修正函数。