基于周期性快速独立向量分析的胎儿心电信号识别方法技术

本发明专利技术公开了一种基于周期性快速独立向量分析的胎儿心电信号识别方法，包括：1、采集孕妇腹部处的表面心电信号；2、表面心电信号采用快速独立向量分析分解腹部心电信号，估计初始母体心电R波尖峰序列；3、将初始母体心电R波尖峰序列作为相对约束，采用基于周期性约束的快速独立向量分析估计母体心电R波尖峰序列；4、采用奇异值分解估计母体心电信号波形；5、从胎儿多通道心电信号中去除母体心电信号波形，得到估计的剩余信号，并作为胎儿多通道剩余心电信号；6、按照步骤2

【技术实现步骤摘要】
基于周期性快速独立向量分析的胎儿心电信号识别方法


[0001]本专利技术属于生物信号处理
，具体涉及一种基于周期性快速独立向量分析的胎儿心电信号识别方法，主要应用于胎儿心率检测。

技术介绍

[0002]胎儿心电信号(Fetal Electrocardiogram,FECG)能够反映胎儿心脏活动的电生理信号，其中，胎儿的R波代表心脏左右心室进行去极化操作，R波提取是胎儿心脏信号特征识别的重要步骤，根据相邻R波得到R
‑
间隔，能够对胎儿心率进行计算。腹部表面心电信号(Abdominal surface Electrocardiogram,AECG)是一种在孕妇腹部表面采集得到的心电信号，具有无创、安全、便于携带和采集等优点，研究者希望能够从腹部心电信号提取胎儿心电信号，以实现无创胎儿心率监测。然而，由于FECG信号振幅非常微弱，且在腹部采集的心电信号是由母体心电信号(Maternal Electrocardiogram,MECG)、FECG、表面肌电信号(Surface Electromyography,SEMG)，测量噪声本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于周期性快速独立向量分析的胎儿心电信号识别方法，其特征在于，是按如下步骤进行：步骤1、获取孕妇腹部处表面所采集的N通道表面心电信号并进行去噪处理后，得到腹部表面心电信号数据，记为E
′
(t)＝[e1(t)，e2(t)，...，e
l
(t)，...，e
N
(t)]
T
，其中，e
l
(t)表示t时刻去噪后的第l个通道的心电信号；T表示转置；步骤2、采用快速独立向量分析法分解所述腹部表面心电信号数据E
′
(t)，以估计初始母体心电R波尖峰序列r：步骤2.1、定义并初始化母体多通道剩余心电信号e
r
＝E
′
(t)；步骤2.2、对母体多通道剩余心电信号e
r
进行扩展与白化处理，得到母体球化信号步骤2.3、对母体球化信号进行快速独立分量分解，得到母体解混向量W；步骤2.4、计算初始心电信号成分从Y＝[y1，y2，...，y
N
]
T
中选择一列分量作为初始母体心电信号成分y，将y中高于选择阈值的所有局部峰值作为初始母体心电R波尖峰序列r；步骤3、将初始母体心电R波尖峰序列r作为约束，采用基于周期性约束的快速独立向量分析法估计最终的母体心电R波尖峰序列步骤3.1、定义并初始化迭代次数k＝1；定义并初始化第k次迭代的两个拉格朗日乘子α
k
，μ
k
，定义并初始化两个惩罚因子γ，β；定义并随机初始化一个解混向量w在第k次迭代的计算结果为w
k
；定义第k次迭代的相关性函数为定义第k次迭代的周期性函数为定义第k次迭代的周期性函数为其中，为y
k
经过延迟θ的信号，为第k次迭代中的相关性下界，为第k次迭代中的周期性下界；y
k
表示第k次迭代中根据求得的母体心电信号成分，(y
k
)
T
为y
k
的转置；E表示期望；步骤3.2、通过式(1)得到解混向量w在第k次迭代正交化后的计算结果步骤3.2、通过式(1)得到解混向量w在第k次迭代正交化后的计算结果式(1)中，G
′
表示非多项式G的一阶函数；表示的一阶函数，表示的一阶函数；步骤3.3、通过式(2)得到解混向量w在第k+1次迭代的计算结果w
k+1
：步骤3.4、通过式(3)和式(4)分别得到第k+1次迭代的两个拉格朗日乘子α
k+1
，μ
k+1
：：步骤3.5、若k＜θ，则将赋值k+1给k，返回步骤3.2顺序执行，否则，表示得到解混向量w在第θ次迭代的计算结果w
θ
；其中，θ代表的迭代次数上限；步骤3.6、根据w
θ
计算最终的母体心电信号成分将y
*

【专利技术属性】
技术研发人员：张旭，骆轩玉，崔嘉雯，李心慧，陈香，
申请(专利权)人：中国科学技术大学，
类型：发明
国别省市：