[发明专利]一种利用数据高维嵌入的微震信号滤波方法和系统

权利要求书

1.所述一种利用数据高维嵌入的微震信号滤波方法，其特征在于，包括：

步骤101获取按时间顺序采集的信号序列S；

步骤102求取困惑度P_i和联合自相似值p_ij，求取公式为：和其中，为第i个Shannon熵；i＝1,2,···,N为所述Shannon熵序号；j＝1,2,···,N为所述联合自相似值p_ij的列序号；N为所述信号序列S的长度；为第一自相似值；为第二自相似值；σ₀为所述信号序列S的均方差；s_i为所述信号序列S的第i个元素；s_j为所述信号序列S的第j个元素；k为求和序号，取值范围为k＝1,2,···,N；

步骤103创建迭代控制参数o，并令其值为1；求取参数矢量d的初始化值d⁰；所述初始化值d⁰的第l个元素为其中l为元素序号，取值范围为l＝1,2,···,N；表示均值为0均方差为10^-4的高斯随机分布函数；令所述参数矢量d的第一步值d¹＝d⁰；

步骤104所述迭代控制参数o的值加1并更新所述参数矢量d的值，具体为：d^o＝d^o-1+α(o)[d^o-1-d^o-2]；其中，d^o为所述参数矢量d的第o步值；d^o-1为所述参数矢量d的第o-1步值；d^o-2为所述参数矢量d的第o-2步值；为高维嵌入因子；

步骤105求取逼近误差ε，求取公式为：ε＝|d^o-d^o-1|；

步骤106判断所述逼近误差ε是否大于或等于预设阈值0.001；如果所述逼近误差ε大于或等于所述预设阈值0.001，则返回所述步骤104、所述步骤105和所述步骤106；直至所述逼近误差ε小于所述预设阈值0.001；

步骤107求取滤除了噪声的信号序列S_NEW，具体为S_NEW＝d^o+α(o)[d⁰-d^o-1]。

2.所述一种利用数据高维嵌入的微震信号滤波系统，其特征在于，包括：

模块201获取按时间顺序采集的信号序列S；

模块202求取困惑度P_i和联合自相似值p_ij，求取公式为：和其中，为第i个Shannon熵；i＝1,2,···,N为所述Shannon熵序号；j＝1,2,···,N为所述联合自相似值p_ij的列序号；N为所述信号序列S的长度；为第一自相似值；为第二自相似值；σ₀为所述信号序列S的均方差；s_i为所述信号序列S的第i个元素；s_j为所述信号序列S的第j个元素；k为求和序号，取值范围为k＝1,2,···,N；

模块203创建迭代控制参数o，并令其值为1；求取参数矢量d的初始化值d⁰；所述初始化值d⁰的第l个元素为其中l为元素序号，取值范围为l＝1,2,···,N；表示均值为0均方差为10^-4的高斯随机分布函数；令所述参数矢量d的第一步值d¹＝d⁰；

模块204所述迭代控制参数o的值加1并更新所述参数矢量d的值，具体为：d^o＝d^o-1+α(o)[d^o-1-d^o-2]；其中，d^o为所述参数矢量d的第o步值；d^o-1为所述参数矢量d的第o-1步值；d^o-2为所述参数矢量d的第o-2步值；为高维嵌入因子；

模块205求取逼近误差ε，求取公式为：ε＝|d^o-d^o-1|；

模块206判断所述逼近误差ε是否大于或等于预设阈值0.001；如果所述逼近误差ε大于或等于所述预设阈值0.001，则返回所述模块204、所述模块205和所述模块206；直至所述逼近误差ε小于所述预设阈值0.001；

模块207求取滤除了噪声的信号序列S_NEW，具体为S_NEW＝d^o+α(o)[d⁰-d^o-1]。