心磁信号广义s变换域奇异值分解滤波方法-物理学报.pdf

物理学报 Acta Phys. Sin. Vol. 62, No. 14 (2013) 148702 心磁信号广义S 变换域奇异值分解滤波方法* 尹柏强1† 何怡刚12 吴先明1 1) ( 湖南大学电气与信息工程学院, 长沙 410082 ) 2) ( 合肥工业大学电气与自动化工程学院, 合肥 230009 ) ( 2012年12月27 日收到; 2013年3月26 日收到修改稿) 针对心磁信号工频及背景噪声干扰问题, 提出了广义S 变换奇异值分解(singular value decomposition, SVD) 滤 波方法. 在离散S 变换基础上, 导出了广义矩阵S 变换和逆变换公式. 通过对采样信号进行广义S 变换, 调节时频分 辨率, 利用SVD 分解方法确定有效心磁信号区域, 实现自适应时频滤波. 实验结果表明, 该方法能有效滤除工频及背 景噪声干扰, 且在较少奇异值个数情况下可获得更好的滤波性能. 关键词: 心磁信号, S 变换, 奇异值分解, 时频滤波 PACS: 87.80.y, 87.85.d, 05.40.Ca DOI: 10.7498/aps.62.148702 征, 在消除50 Hz 工频干扰时, 仍需采用自适应滤 1 引言 波方法. 文献[7] 利用小波变换来处理心磁信号, 采 用symlet8 小波函数进行消噪, 取得了良好的滤波 1 心脏疾病是危害人类健康的主要疾病之一 . 效果, 但是寻找合适的小波基函数较为复杂, 且在 目前, 心脏功能检测和诊断方法主要采用心电图来 高频段需要强制消噪处理. S 变换具有与小波变换 2 检测和分析心脏电活动 . 随着磁传感器件的发 8 类似的多分辨特性 , 利用S 变换多分辨二维时频 展, 高温射频超导量子干涉仪作为灵敏度极高的磁 特性可以构建反S 变换时频滤波方法913 . 目前, 3 传感器件广泛地应用于心磁信号的测量中 . 高灵 有关S 变换滤波方法主要采用全局阈值滤波和局 敏度的高温射频超导量子干涉仪采集到的心磁信 部阈值滤波1114 . 局部阈值滤波在有先验知识的 号中存在大量的背景噪声, 其主要成分为50 Hz 工 4 前提下, 可以有效实现时频滤波. 然而, 对于实测环 频干扰及其他环境噪声 . 心磁信号特征的有效提 境下的心磁信号, 除工频干扰外, 背景噪声干扰也 取是心脏疾病诊断的主要关键. 文献[5] 提出了针对背景噪声消除问题的自适 是影响信号质量的主要因素, 如何在背景噪声中快 应理论, 所提出自适应算法适用于心磁信号及环境 速有效地提取心磁信号仍是一个值得深入研究的 噪声可同步采集情形, 但不能解决无参考噪声数据 问题. 综上所述, 本文在S 变换