[发明专利]无创胎儿心电信号R波的检测方法

权利要求书

1.无创胎儿心电信号R波的检测方法，其特征在于，它包括以下步骤：

步骤一、对孕妇腹部原始信号进行预处理，对电源的工频干扰信号进行数字带通滤波；

步骤二、检测孕妇的心电信号R波，从而获得孕妇心电信号QRS波群；

步骤三、在孕妇的两个QRS波群之间，检测胎儿心电信号的R波；

步骤四、根据心电信号的准周期特性，从已检测的胎儿心电信号的R波，预测与孕妇心电信号QRS波群重合的胎儿心电信号R波的区域位置；

步骤五、根据6个相邻的胎儿心电信号QRS波群，使用奇异值分解方法建立一个胎儿心电信号QRS波群模板；根据6个相邻的不与胎儿信号重合孕妇心电信号QRS波群，使用奇异值分解方法建立一个孕妇心电信号QRS波群模板；

步骤六、根据胎儿心电信号QRS波群模板和孕妇心电信号QRS波群模板，在预测出区域位置融合，其融合后的信号与实际的信号计算相关系数，其最大的相关系数，即为胎儿心电信号R波的准确位置。

2.根据权利要求1所述的无创胎儿心电信号R波的检测方法，其特征在于，在步骤一中，在步骤五中，所述奇异值分解方法具体为：根据已经检测到的6个相邻的不与胎儿信号重合孕妇心电信号，我们可以建立以下矩阵，

其中每行，代表一个已检测孕妇心电信号，X矩阵为6行，即p＝6，根据孕妇心电信号的特点，我们把孕妇心电信号长度设置为0.16s，如果采样为1khz，即q＝160。其中所有信号R波列对齐，位置在信号的中间。

X矩阵可以表示为如下分解：

X＝USV^T

其中，U是m×m阶酉矩阵；S是半正定m×n阶对角矩阵；V是n×n阶酉矩阵。这样的分解就称作X的奇异值分解。S对角线上的元素即为X的奇异值。

3.根据权利要求1所述的无创胎儿心电信号R波的检测方法，其特征在于，在步骤一中，在步骤五中，所述胎儿心电信号的QRS波群模板为：

m＝[t_mQRS]_q

所述心电信号的QRS波群模板为：

f_i＝[0 t_fQRS 0]_q

其中，t_mQRS孕妇心电信号QRS波群的模板，q代表了采样信号的长度；t_fQRS胎儿心电信号QRS波群的模板，q代表了采样信号的长度。

4.根据权利要求1所述的无创胎儿心电信号R波的检测方法，其特征在于，在步骤一中，在步骤六中，根据获得10％refRR区域，进行，定义相关系数计算公式：

代表着M_i向量，t_c代表着真实的信号；

其中，r_i是相关系数，和是向量和t_c的均值；

确定位置如下公式：I＝argmaxr_i，即在获得相关系数，通过比较，获得的最大的相关系数，即是确定的胎儿心电信号的位置。