[发明专利]心磁信号噪声自适应滤波消除设计方法

权利要求书

1.心磁信号噪声自适应滤波消除设计方法，其特征在于，包括以下步骤： (1) 对实时采集到的心磁信号进行S变换，获得信号时频域内特征矩阵S变换时频矩阵； (2) 根据S变换时频矩阵时频分辨率的实时要求高低，在S变换中引入时频调节因子                                               、构建广义S变换，利用广义S变换时频调节因子调节信号时频特征的时频分辨率；(3) 采用奇异值分解方法分解广义S变换二维时频特征矩阵，得到左特征时域矩阵、右特征频域矩阵及特征值组成的对角阵，将单个奇异值对角阵分别与左右特征时域矩阵相乘得到单奇异值S变换时频矩阵，单奇异值S变换时频矩阵经反S变换即可得到单个主要奇异值所对应的信号时域成份；(4) 计算单个奇异值时域信号与噪声的相关系数，取相关系数r<0.5且单奇异值占总奇异值比例大于10%的奇异值所对应奇异矩阵区域为有效心磁信号奇异值分布区域，对有效心磁信号时频域进行反S变换重构信号从而实现心磁信号自适应时频滤波。

2.根据权利要求1所述的心磁信号噪声自适应滤波消除设计方法，其特征在于，所述步骤（1）中，所述心磁信号是指工作在无磁屏蔽室条件下利用高温射频超导量子干涉仪系统所测量到的心脏电心理过程在躯干上面产生的时变磁场强度信号；

所述步骤（1）中，所述S变换是指以高斯函数为窗函数的信号短时傅里叶变换，其表达式为，式中，为时间，为频率，为时移因子，为虚数，为表示实时连续心磁信号，为高斯窗函数，其中。

3.根据权利要求2所述的心磁信号噪声自适应滤波消除设计方法，其特征在于，所述步骤（2）中，所述广义S变换是指在S变换中引入时频调节因子、，其表达式为；其中，为高斯窗幅度拉伸因子，为频率尺度拉伸因子，为时间，为频率，为时移因子，为虚数，为表示实时连续心磁信号。

4.根据权利要求3所述的心磁信号噪声自适应滤波消除设计方法，其特征在于，所述步骤（3）中，所述奇异值分解方法是指一种非退化正交矩阵分解法，具体是指心磁信号经广义S变换后的二维时频特征矩阵正交分解方法。