[发明专利]一种能量分解中功率信号滤波方法及系统

权利要求书

1.一种能量分解中功率信号滤波方法，其特征在于，包括：

获取实测功率信号序列以及多尺度分解层数；所述多尺度分解层数用于对所述实测功率信号序列进行分解；

根据所述多尺度分解层数确定多层低通滤波器截止频率；

根据所述多尺度分解层数以及所述多层低通滤波器截止频率对所述实测功率信号序列逐层进行低通滤波处理，确定每一层滤波后的实测功率信号序列；

获取当前层滤波后的实测功率信号序列；

根据所述当前层滤波后的实测功率信号序列确定当前层时间窗口长度；

对所述当前层滤波后的实测功率信号序列内的功率信号进行逐点处理，构建解析序列；

确定所述解析序列在所述当前层时间窗口长度下的伪Wigner-Ville分布；

按照所述多尺度分解层数从大到小的顺序，根据所述伪Wigner-Ville分布恢复每一层所述滤波后的实测功率信号序列，确定恢复实测功率信号序列；

按照所述多尺度分解层数从大到小的顺序对所述恢复实测功率信号序列进行线性插值处理，确定插值后实测功率信号序列；

从所述多尺度分解层数的最后一层开始，根据当前层的所述恢复实测功率信号序列以及上一层的所述插值后实测功率信号序列滤除所述实测功率信号中的非平稳非高斯噪声；对第C层恢复数据进行线性插值处理，得到插值后的数据序列要求的采样率是的两倍，此时的数据采样率与第C-1层的数据采样率是相同的；其中，C为多尺度分解层数的总层数；

对层序号按照递减的顺序进行下面的处理：假设当前层的序号为c；

当前层有两个数据序列：原来的恢复数据以及第c+1层恢复数据经过插值处理后得到的数据序列；根据这两个序列得到当前层的最终数据恢复序列

对当前层的最终数据恢复序列进行线性插值，得到插值后的数据序列要求的采样率是的两倍；此时的数据采样率与第c-1层的数据采样率是相同的；

重复上述步骤，直至c＝1，结束，得到实测功率信号序列经过滤波之后的数据序列此序列即为所求序列，处理结束，滤除实测功率信号序列中的非平稳非高斯噪声。

2.根据权利要求1所述的滤波方法，其特征在于，所述根据所述多尺度分解层数确定多层低通滤波器截止频率，具体包括：

根据公式确定多层低通滤波器截止频率；其中，为多层低通滤波器截止频率；c为多尺度分解层数，C为多尺度分解层数的总层数，c＝1,2,…,C。

3.根据权利要求2所述的滤波方法，其特征在于，所述根据所述当前层滤波后的实测功率信号序列确定当前层时间窗口长度，具体包括：

获取所述当前层滤波后的实测功率信号序列的采样频率以及主频；

根据公式确定当前层时间窗口长度；其中，为当前层时间窗口长度；为所述当前层滤波后的实测功率信号序列P^c的采样频率；为所述当前层滤波后的实测功率信号序列P^c的主频；T^c为所述当前层滤波后的实测功率信号序列的采样间隔。

4.根据权利要求3所述的滤波方法，其特征在于，所述对所述当前层滤波后的实测功率信号序列内的功率信号进行逐点处理，构建解析序列，具体包括：

根据公式对所述当前层滤波后的实测功率信号序列内的功率信号进行逐点处理，构建解析序列；其中，z^c(n)为构建第c层第n点的解析值；μ_c为调制系数，1＜μ_c≤2；P^c(j)为当前层信号序列P^c中的第j个元素。