小波域变步长ICA算法提取胎儿心电信号.pdfVIP

第 34卷第 1期 杭 州 电 子 科 技 大 学 学 报 Vo1．34。No．1 2014年1月 Journal of ／／an~hou Dianzi University Jan．2014 doi：10．3969／j．issn．1001—9146．2014．Ol-005 小波域变步长 ICA算法提取胎儿心电信号 占海龙 ，赵治栋 (1．杭州电子科技大学通信工程学院，浙江杭州I310018； 2．杭州电子科技大学电子信息学院，浙江杭州310018) 摘要：胎儿心电信号是反映胎儿宫内生理活动的一种最主要的客观指标，从母体体表混合信号 中 分离出胎儿心电信号是近年来信号处理领域研究的热点问题。该文利用小波变换将采集到的各 导联混合信号变换到小波域 ；在小波域结合模拟退火法变步长进行独立分量的提取；将提取到的 各独立分量进行小波逆变换 ，从中选择胎儿心电信号成分。通过对临床实际数据进行实验，实验 结果提取出了清晰的胎儿心电信号。 关键词：胎儿心电信号提取；小波变换；独立分离提取 中图分类号 ：RS18 文献标识码 ：A 文章编号 ：1001—9146(2014)01—0025—05 0 引 言 胎儿在母亲宫内的健康状况近年来 日益引起人们的重视，多普勒胎心听诊、胎心电子监护仪等都被 用来检测胎儿的状况，但这些监护方法存在缺点，影响着传统胎儿监护技术的稳定性 J。通过置电极 于母体腹部体表可采集到包含胎儿心电信号的体表混合信号，然而从母体腹部提取到的信号中包含了 很多干扰成分，用常规的滤波方法很难提取到清晰的胎儿心电信号(Fetalelectrocardiogram，FECG)J。 近年来国内外众多学者在胎儿心电信号提取领域进行了广泛的研究，并取得了一定的研究成果。自适 应消噪法应用在胎儿心电信号提取中 ；利用奇异值分解方法分离出较为清晰的胎儿心电信号 ；利 用奇异值分解法结合多项式神经网络方法提取出清晰的胎儿心电信号 ；利用改进独立分量分析 (IndependentComponentAnalysis，ICA)算法分离出较为清晰的胎儿心电信号 J。但这些方法都存在一 定的缺陷，如 自适应消噪法对参考电极的要求高，奇异值分解法对电极放置位置要求较高等。本文主要 基于小波变换和独立分量分析，利用系统辨识数据库中的8导孕妇体表信号进行分析方法，结合模拟退 火法中的指数衰减型变步长算法来优化 ICA算法中的步长选择问题，方法能够在兼顾收敛速度和稳态 误差的前提下准确地提取出胎儿心电信号。 1 基于小波域的ICA算法 设 S=[s，s，…，s]为n个相互统计独立的源信号，A为一非奇异混合矩阵，观测信号为X=[x。， X，X，，…，X]，则有X=AS。ICA算法的目的是通过寻找一解混矩阵w，使得输出信号Y中的各分量 相互统计独立，且为源信号s中各分量的估计 ，允许幅度和顺序存在差异，即：Y=WX=WAS，其中Y= [Y，Y：，Y，，…，Y]。解混矩阵w 的更新过程如下： 收稿 日期：2013—06—20 基金项 目：浙江省科技计划公益技术研究资助项目(2012C21005) 作者简介：占海龙(1988一)，男，江西上饶人，在读研究生，信号与信息处理 杭 州 电子 科 技 大 学 学 报 2014年 W(k+1)=W(k)+u(I一_厂(Y())Y (k))W(k) (1) w(k)为k时刻矩阵w 的值，I为单位矩阵 ·)为激活函数，U为步长。对于离散信号x(k)，对 其进行小波变换后的小波系数为： cJ，q=X(k)，Oj，(k) (2) 由于内积具有线性特征，即：(x(k)+Z(k))，Y(k)=X(k)，y(k)+Z(k)，Y(k)， a】I【(k)，Y(k)=ax(k)，Y(k)，因此，对于混合信号