[发明专利]一种用于含噪ENPEMF信号的BGabor-NSPWVD三维时频分析方法

说明书

本发明提供了一种用于含噪ENPEMF信号的BGabor‑NSPWVD三维时频分析方法，所述方法包括步骤：获取输入的ENPEMF信号；对所述ENPEMF信号进行Gabor变换得到Gabor数组；将所述Gabor数组二值化得到BGabor数组；对所述ENPEMF信号进行SPWVD变换得到SPWVD数组；将所述SPWVD数组进行归一化处理得到NSPWVD数组；所述NSPWVD数组点除以所述BGabor数组得到临时数组；设置合适的上限值与下限值调整所述临时数组得到调整后的临时数组；修正SPWVD数组得到修正后的SPWVD数组；所述修正后的SPWVD数组点除以调整后的临时数组得到BGabor‑NSPWVD数组。本发明能够适应噪声信号处理，实现自适应频率赋值，降低含噪信号在扫频时产生的虚假分量，在ENPEMF数据分析中得到良好的应用。

技术领域

本发明涉及ENPEMF分析领域，具体涉及一种用于含噪ENPEMF信号的BGabor-NSPWVD三维时频分析方法。

背景技术

非平稳信号的时频分析研究主要应用于ENPEMF信号分析领域，包括：地表磁异常监测、烃类油气勘探、滑坡地表电磁异常监测、大坝电磁安全监测等等多个领域磁异常监测领域。同时也可应用于其他非平稳时序信号领域，如地球物理探测领域、信息通信领域、自然灾害和环境监测领域、非平稳故障信号诊断领域。但是，ENPEMF信号分析存在三维时频分布不准确的问题。因此，如何降低ENPEMF含噪信号在扫频时产生的虚假分量就成为亟待解决的问题。

发明内容

本发明提供了一种用于含噪ENPEMF信号的BGabor-NSPWVD三维时频分析方法，可以有效解决上述问题。

本发明提供的技术方案是：一种用于含噪ENPEMF信号的BGabor-NSPWVD三维时频分析方法，所述方法包括步骤：获取输入的ENPEMF信号；对所述ENPEMF信号进行Gabor变换得到Gabor数组；将所述Gabor数组二值化得到BGabor数组；对所述ENPEMF信号进行SPWVD变换得到SPWVD数组；将所述SPWVD数组进行归一化处理得到NSPWVD数组；所述NSPWVD数组点除以所述BGabor数组得到临时数组；设置合适的上限值与下限值调整所述临时数组得到调整后的临时数组；修正SPWVD数组得到修正后的SPWVD数组；所述修正后的SPWVD数组点除以调整后的临时数组得到BGabor-NSPWVD数组。

本发明的有益效果是：本发明提供了一种用于含噪ENPEMF信号的BGabor-NSPWVD三维时频分析方法，通过综合Gabor变换与和SPWVD时频分辨率高的特点，不仅克服了Gabor变换时频分辨率固定且不够的缺点，还克服了WVD会产生交叉项的问题。同时，对Gabor数组进行二值化处理，对SPWVD数组进行归一化处理，使得最终方法具有高度的鲁棒性。本发明能够适应噪声信号处理，实现自适应频率赋值，降低含噪信号在扫频时产生的虚假分量，在ENPEMF数据分析中得到良好的应用。

附图说明

图1是本发明实施例中一种用于含噪ENPEMF信号的BGabor-NSPWVD三维时频分析方法的整体流程图；

图2是本发明实施例中使用不同算法处理ENPEMF信号的时频二维效果示意图；

图3是本发明实施例中使用不同算法处理ENPEMF信号的时频-幅值三维效果示意图；

图4是本发明实施例中使用BGabor-NSPWVD方法对2013年芦山地震期间的ENPEMF信号进行时频分析的时频-幅值效果示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述，下文中提到的具体技术细节，如：方法等，仅为使读者更好的理解技术方案，并不代表本发明仅局限于以下技术细节。