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

权利要求书

1.一种基于对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处理，得到时频分布。

2.根据权利要求1所述的方法，其特征在于，所述步骤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频率数组中取出相同行数中频率值作为主要频率的值。

3.根据权利要求2所述的方法，其特征在于，所述k1等于0.05，所述k2等于0.01。

4.一种基于对ENPEMF信号进行BSWT-DDTFA时频分析的系统，用于对地震进行预测，其特征在于，包含：

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

各时频分布均采用下述模块获取得到：

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

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

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

时频分布获取模块，用于将得到的IMF分量的瞬时频率特性中各IMF分量分别进行SWT处理，得到时频分布。