[发明专利]一种微弱表面肌电信号中的工频干扰消除方法

权利要求书

1.一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述消除方法包括以下步骤：

获取微弱表面肌电原信号；

将微弱表面肌电原信号按照固定时间间隔采样，得到表面肌电采样信号，所述采样时钟与工频干扰严格同步；

将表面肌电采样信号分组组成多组表面肌电采样分组信号；

将多组表面肌电采样分组信号进行滤波处理，拼合成为工频干扰信号；

从微弱表面肌电原信号中减去工频干扰信号，得到消除工频干扰后的信号。

2.根据权利要求1所述的一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述消除方法还包括表面肌电采样信号过滤步骤：

变化率计算步骤：计算表面肌电采样分组信号的变化率k_n；

离散距离计算步骤：计算变化率k_n的离散距离D_n；

阈值计算步骤：根据离散距离D_n计算阈值G；

过滤步骤：根据阈值G对表面肌电采样分组信号进行过滤。

3.根据权利要求2所述的一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述变化率计算步骤具体方法如下：

对于采样信号：Q_电＝{Q₁,Q₂,Q₃……Q_n},计算采样信号的变化率k为：

采样时间的常数t通过采样频率得到。

4.根据权利要求3所述的一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述离散距离计算步骤的具体方法如下：

变化率k的聚散中心点N为：

利用二阶Minkowski距离法，计算变化率k与聚散中心点的距离：

5.根据权利要求4所述的一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述阈值计算步骤的具体计算方法如下：

G＝g*m(D)

g为根据系统属性设立的安全系数；m(D)为中位数求值函数。

6.根据权利要求5所述的一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述过滤步骤的具体方法如下：

剔除系数为a_n,D_nG时，a_n＝0，D_n≤G时，a_n＝1。

7.根据权利要求2所述的一种微弱表面肌电信号中的工频干扰消除方法，其特征在于，所述消除方法还包括工频干扰同步步骤：

对于过滤后的表面肌电采样分组信号：Q_滤＝{Q₁,Q₂,Q₃……Q_n},计算采样信号的变化率k为：

所述采样时间的常数t通过采样频率得到；

根据变化率k的数据集合k＝{k₂,k₃……k_n}，获计算离散距离D_n，

判断离散距离D_n是否在预设阈值范围之内，若超出预设阈值，则在获取下一组表面肌电采样分组信号时增加或减少采样间隔。