[发明专利]一种心磁信号背景噪声S变换域去除方法及系统

权利要求书

1.一种心磁信号背景噪声S变换域去除方法，其特征在于，包括以下步骤：

S1：采集噪声数据点

实时采集仪器所在环境下的背景噪声N个数据点；

S2：对噪声进行S变换

对采集的N个数据点背景噪声采用S变换方法变换成二维时频矩阵；

S3：绘制最大值二维时频曲面

连续采集多次，计算每次采集信号对应的二维时频矩阵，得到多组二维时频矩阵，计算统计过程中多组每一个时频点时频矩阵的最大值，得到统计情况下的背景噪声S变换二维时频矩阵，绘制矩阵得到最大值二维时频曲面；

S4：确定区分临界界面

取最大值曲面的设定置信度曲面作为噪声与心磁信号的区分临界界面；

S5：比较心磁信号二维时频矩阵与区分临界界面的关系

对采集的心磁信号采用S变换方法变换成二维时频矩阵，当采集的心磁信号二维时频矩阵大于区分临界界面时，保留心磁信号的二维时频矩阵；当心磁信号的二维时频矩阵小于或等于区分临界界面时，将心磁信号的二维时频矩阵置为0，对心磁信号的二维时频矩阵进行修改；

S6：进行反S变换

对修改后的心磁信号二维时频矩阵进行反S变换，得到去除背景噪声后的心磁信号。

2.根据权利要求1所述的一种心磁信号背景噪声S变换域去除方法，其特征在于：在所述步骤S1中，所述仪器所在环境为心磁信号测量前的仪器所在环境。

3.根据权利要求1所述的一种心磁信号背景噪声S变换域去除方法，其特征在于：在所述步骤S2中，S变换为采用高斯窗函数的短时傅里叶变换，其表达式为：其中，t为时间，f为频率，τ为时移因子，i为虚数，x(t)为表示实时连续噪声。

4.根据权利要求1所述的一种心磁信号背景噪声S变换域去除方法，其特征在于：在所述步骤S4中，设定置信度为95％。

5.根据权利要求1所述的一种心磁信号背景噪声S变换域去除方法，其特征在于：在所述步骤S6中，反S变换为采用高斯窗函数的短时傅里叶逆变换，其表达式为：其中，u，t均为时间，f为频率，τ为时移因子，i为虚数，x(u)为表示去除背景噪声后的心磁信号。

6.一种心磁信号背景噪声S变换域去除系统，采用如权利要求1～5任一项所述的背景噪声去除方法对心磁信号中的背景噪声进行去除，包括：

噪声采集模块，用于实时采集仪器所在环境下的背景噪声N个数据点；

噪声变换模块，用于对采集的N个数据点背景噪声采用S变换方法变换成二维时频矩阵；

绘制模块，用于连续采集多次，计算每次采集信号对应的二维时频矩阵，得到多组二维时频矩阵，计算统计过程中多组每一个时频点时频矩阵的最大值，得到统计情况下的背景噪声S变换二维时频矩阵，绘制矩阵得到最大值二维时频曲面；

临界界面确定模块，用于取最大值曲面的设定置信度曲面作为噪声与心磁信号的区分临界界面；

比较模块，用于对采集的心磁信号采用S变换方法变换成二维时频矩阵，当采集的心磁信号二维时频矩阵大于区分临界界面时，保留心磁信号的二维时频矩阵；当心磁信号的二维时频矩阵小于区分临界界面时，将心磁信号的二维时频矩阵置为0，对心磁信号的二维时频矩阵进行修改；

反变换模块，用于对修改后的心磁信号二维时频矩阵进行反S变换，得到去除背景噪声后的心磁信号；

中央处理模块，用于向其他模块发出指令，完成相关动作；

所述噪声采集模块、噪声变换模块、绘制模块、临界界面确定模块、比较模块、反变换模块均与中央处理模块电连接。