[发明专利]一种基于对ENPEMF信号进行BSWT-DDTFA时频分析的方法及系统

说明书

本发明公开了一种基于对ENPEMF信号进行BSWT‑DDTFA时频分析的方法及系统，用于对地震进行预测，分别采用BSWT‑DDTFA时频分析方法对实时的地球天然脉冲电磁场数据与地震发生前地球天然脉冲电磁场数据进行分析，分别得到时频分布，将时频分布进行比较后，判断地震发生的几率。本发明的得到时频分布时，能够有效抑制背景噪声和其他高斯白噪声源，得到的时频分布时频聚集特性优异，对地震发生几率的判断准确性高。

技术领域

本发明涉及地震前兆研究与地震预测领域，更具体地说，涉及一种基于对ENPEMF信号进行BSWT-DDTFA时频分析的方法及系统。

背景技术

地震给人类的生活带来了巨大的灾难，据统计，全球的自然灾害之中，地震造成的死亡人数占全部自然灾害死亡人数的54％，堪称自然灾害之最。如何预测地震一直以来都是一个热门敏感的话题。然而，因地震预测有着地球内部的“不可入性”、大地震的“非频发性”、“地震物理过程的复杂性”三大困难，地震预测成为了公认的世界性难题，对于地震预测前兆方面，目前没有一个精确的手段可以预测地震发生。

发明内容

本发明要解决的技术问题在于，针对上述的现有技术中目前没有一个精确的手段可以预测地震发生的技术缺陷，提供了一种基于对ENPEMF信号进行BSWT-DDTFA时频分析的方法及系统，以对地震进行预测。

根据本发明的其中一方面，本发明为解决其技术问题，提供了一种基于对ENPEMF信号进行BSWT-DDTFA时频分析的方法，用于对地震进行预测，包含下述步骤：

将根据实时采集的ENPEMF信号得到的时频分布与根据预先采集的地震前的ENPEMF信号得到的时频分布进行对比，根据二者的一致性判断地震发生的几率，其中二者一致性越好，地震发生几率越高；

各时频分布均采用下述步骤获取得到：

S1、对ENPEMF信号进行SWT处理，得到ENPEMF信号的SWT时频分布矩阵Tf(t,f)以及SWT频率数组；

S2、根据得到的SWT时频分布矩阵Tf(t,f)，从SWT频率数组提取信号的主要频率的值；

S3、将提取出来的频率值乘以2πt分别作为DDTFA方法的初始相位值，得到IMF分量的瞬时频率特性；

S4、将得到的IMF分量的瞬时频率特性中各IMF分量分别进行SWT处理，得到时频分布。

进一步的，在本发明的基于对ENPEMF信号进行BSWT-DDTFA时频分析的方法中，步骤S2包含下述子步骤：

S21、找到SWT时频分布矩阵Tf(t,f)中的最大值max；

S22、将SWT时频分布矩阵Tf(t,f)中的大于k₁*max的值全部置1，小于k₁*max的值全部置0，k₁为预设的系数，0<k₁<1；

S23、对SWT时频分布矩阵Tf(t,f)的每一行的数据进行求和，找到和最大的数hangmax，依次获取和大于k2*hangmax的各行坐标形成数组hangshu()，k₂为预设的系数，0<k₂<1；

S24、将数组hangshu()中每一类数分别求平均，将得到的平均值分别作为各主要频率所在的行坐标，其中数组hangshu()中数的分类按照下述方法实现：判断数组hangshu()中相邻两个数之间的差是否小于预设值，若是，则归为一类，否则归为不同类；

S25、通过各主要频率所在的行坐标，从SWT频率数组中取出相同行数中频率值作为主要频率的值。