CN107157474B - 一种便携式多胞胎心电监测仪及其监测方法 - Google Patents
一种便携式多胞胎心电监测仪及其监测方法 Download PDFInfo
- Publication number
- CN107157474B CN107157474B CN201710500220.4A CN201710500220A CN107157474B CN 107157474 B CN107157474 B CN 107157474B CN 201710500220 A CN201710500220 A CN 201710500220A CN 107157474 B CN107157474 B CN 107157474B
- Authority
- CN
- China
- Prior art keywords
- signals
- signal
- matrix
- imf
- fetal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000001605 fetal effect Effects 0.000 title claims description 37
- 210000002458 fetal heart Anatomy 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 210000003754 fetus Anatomy 0.000 claims abstract description 22
- 230000008774 maternal effect Effects 0.000 claims abstract description 20
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 210000001015 abdomen Anatomy 0.000 claims abstract description 6
- 230000003321 amplification Effects 0.000 claims abstract description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 5
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 6
- 230000003187 abdominal effect Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000000747 cardiac effect Effects 0.000 claims description 3
- 238000004891 communication Methods 0.000 claims description 3
- 238000000718 qrs complex Methods 0.000 claims description 3
- 239000004973 liquid crystal related substance Substances 0.000 claims 1
- 230000003183 myoelectrical effect Effects 0.000 claims 1
- 230000035935 pregnancy Effects 0.000 abstract description 4
- 238000000926 separation method Methods 0.000 abstract description 3
- 210000003815 abdominal wall Anatomy 0.000 abstract description 2
- 230000003044 adaptive effect Effects 0.000 abstract 1
- 230000036541 health Effects 0.000 description 6
- 238000002565 electrocardiography Methods 0.000 description 4
- 210000003954 umbilical cord Anatomy 0.000 description 4
- 238000004804 winding Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 208000006893 Fetal Hypoxia Diseases 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 208000019622 heart disease Diseases 0.000 description 2
- 230000035764 nutrition Effects 0.000 description 2
- 235000016709 nutrition Nutrition 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 206010003497 Asphyxia Diseases 0.000 description 1
- 208000032943 Fetal Distress Diseases 0.000 description 1
- 206010016855 Foetal distress syndrome Diseases 0.000 description 1
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000005246 left atrium Anatomy 0.000 description 1
- 229910001416 lithium ion Inorganic materials 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000009984 peri-natal effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000034225 regulation of ventricular cardiomyocyte membrane depolarization Effects 0.000 description 1
- 210000005245 right atrium Anatomy 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 210000004291 uterus Anatomy 0.000 description 1
- 210000001215 vagina Anatomy 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/366—Detecting abnormal QRS complex, e.g. widening
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4343—Pregnancy and labour monitoring, e.g. for labour onset detection
- A61B5/4362—Assessing foetal parameters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/04—Constructional details of apparatus
- A61B2560/0431—Portable apparatus, e.g. comprising a handle or case
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Signal Processing (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Pediatric Medicine (AREA)
- Pregnancy & Childbirth (AREA)
- Gynecology & Obstetrics (AREA)
- Reproductive Health (AREA)
- Power Engineering (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种便携式胎心电监测仪及其监测方法，属于医疗设备领域。该装置包括生物电传感器、信号放大与滤波器、A/D转换器、嵌入式微控制器、存储器、LCD显示器、充电电池、时钟频率模块、蓝牙和NFC模块；该装置检测方法包括步骤：从母体腹部采集生物电信号；采用自适应滤波算法滤除掉肌电信号、噪声干扰；采用小波算法提取出MECG信号；对FECGErr信号，加入高斯白噪声，进行集合经验模态分解；采用X表示分解后的IMF信号，对X进行奇异值分解，求出矩阵X的k个奇异值；并对该k个奇异值排序；对胎儿个数预测。本发明只采用了一路母体腹壁信号实现母胎心电信号的分离，可在孕妇怀孕5个孕周后初次使用，并可以对胎儿个数进行预估。
Description
技术领域
本发明属于医疗设备监测领域，涉及一种便携式多胞胎心电监测仪及其监测方法。
背景技术
全世界每年约有300多万胎儿在孕期死亡，约400多万婴儿在出生后前四周内死亡。仅中国每年就有数千万的孕妇和约2000万的新生儿出生，2002年一年仅北京的婴儿死亡率就有约0.59％，2007年浙江省围产儿死亡率约为1.174％。而在贫困地区，新生儿平均死亡率就更高，约为3.976％。由此可见，为了提高新生儿的质量和生存率，早发现和早预防就显得十分重要。产前检查的主要内容是胎儿健康监护，胎儿健康监护能提前检测到某些胎儿的心脏疾病、胎儿被脐带缠绕引起的窒息和胎儿窘迫等异常情况，实现优生和胎儿健康安全预防，能大大减少胎儿和新生儿的死亡率。
胎儿健康监测的主要指标就是胎心率(fetalheartrate,FHR)，它简单方便的通过数值直观地显示了胎儿的心脏功能，并且与子宫内的营养、供氧状况和脐带缠绕等危险信息息息相关，当胎心率出现明显变化时，医生不仅会考虑胎儿可能的心脏疾病和危险，还会考虑孕期妇女本身的健康状况、营养情况等，并给出诊断建议，实现孕期妇女安全怀孕、顺利生产等。所以，在妇女孕期实施胎心监测，对提高围产期胎儿质量和降低新生儿的死亡率具有十分重要的现实意义。所以，胎儿健康监测在某种意义上就等同于胎心率监测。
目前常用的家用胎心监测仪是基于超声波多普勒效应实现的胎心音监测，胎心音监测仪虽然能够实时计算出胎心率，并能够提供一种直观的胎心跳动的声音。但是胎心音监测仪需要在孕妇腹部涂上耦合剂，并且要移动超声探头直到寻找到最佳监听部位为止。因此胎心音监测仪不适合长时间在线监测，并且超声波对胎儿的影响并不十分明确。而胎心电(Fetal Electrocardiogram,FECG)监测仪可以更方便、更简单地进行胎心监测，胎心电信号中除了可以计算出胎心率之外还蕴含了更多的信息，而且胎心电监测仪能够实现长时间在线监测，以便将来可以实现远程医疗等。
最早的胎心电监测是有创的，需要将电极通过阴道插入子宫，贴到胎儿头皮上。这种方法操作复杂，可能会对孕妇和胎儿造成一定损伤，且不适于早期的胎心电监测。因此，目前的胎心电监测法都是无创的，是在孕妇体表采集信号，从而采用一定的信号处理方法提取出所需的胎心电信号。
使用得最多的胎心电提取方法是基于盲信号处理的胎心电信号提取法，这种方法需要在孕妇体表采集多路信号。并且目前的胎心电信号提取算法都没有考虑多胞胎的情况。所以，目前的胎心电监测仪都只适合医院使用，并不太适合家用的便携式操作环境。
发明内容
有鉴于此，本发明的目的在于提供一种基于时间多尺度实现了母体心电、胎儿心电的分离，并且基于奇异值分解的方法提出了一种多胞胎便携式胎心电监测仪的监测方法，为达到上述目的，本发明提供如下技术方案：
一种便携式多胞胎心电监测仪，包括用于采集母体肌电信号的生物电传感器、用于放大信号与滤掉噪声干扰的信号放大与滤波器、用于把模拟信号转换成数字信号的A/D转换器、用于控制和处理信号的嵌入式微控制器、用于存储信息的存储器、用于显示监测数据的LCD显示器、用于提供电源的充电电池、时钟频率模块、用于非接触式传输数据的蓝牙和近距离无线通信技术(Near Field Communication,NFC)模块；
所述生物电传感器将采集到的母体肌电信号通过所述信号放大与滤波器进行滤波和除噪，再通过所述A/D转换器转换成数字信号输入到所述嵌入式微控制器进行信号处理，将监测结果传送到所述LCD显示器，并通过蓝牙和NFC模块进行无线数据传输；所述嵌入式微控制器将信息存储到所述存储器，所述时钟频率模块控制所述嵌入式微控制器的时钟频率。
一种便携式多胞胎心电监测仪的检测方法，该方法包括以下步骤：
S1：采集母体腹部的信号；
S2：采用自适应滤波算法滤除掉肌电信号、噪声干扰；
S3：采用小波算法提取出母体心电信号MECG；
S4：对去除母体心电信号后的残余信号FECGErr，加入高斯白噪声，进行集合经验模态分解(Ensemble Empirical Mode Decomposition,EEMD)，从而把残余信号分解为若干个本征模态信号(Intrinsic Mode Function,IMF)；采用X表示分解后的IMF信号，即X＝(c1,c2,…,cn,r)，其中ci表示第i个本征模态信号，共n个本征模态信号，r表示分解后的残差；
S5：对X进行奇异值分解(Singular Value Decomposition,SVD)，求出矩阵X*XT的k个特征值，设为λi，其中1≤i≤k，该特征值的算术平方根σi即为矩阵X的k个奇异值；并对该k个奇异值按从大到小的顺序排序为：σ1＞σ2＞…＞σk；
S6：设μi＝σi/σi+1，其中i＝1,2,…,k-1；
求出μi中的最大值，设当i＝j时，μi取得k个奇异值的降低梯度中的最大值μj，即相邻奇异值比值中的最大值μj；
则对胎儿的预测个数为j个，且j个胎心电信号为：
Fecg＝(c1，c2，…，cj-1，cj+cj+1+…cn+r)。
进一步，所述步骤S1具体为：将便携式胎心电监测仪的传感器贴于母体肚脐眼正下方三指宽处，采集到一路母体腹部肌电信号。
进一步，步骤S4中所述EEMD具体包括以下步骤：
S401：在FecgErr中添加一定量的高斯白噪声；
S402：进行初始化操作：l0(t)＝FecgErr，i＝1；
S403：得到第i个IMF；
S403a：h0(t)＝li-1(t)，j＝1；
S403b：找出hj-1(t)的局部极值点；
S403c：对hj-1(t)的极大和极小值点分别进行三次样条函数插值，分别形成hj-1(t)的上下包络线；
S403d：计算上下包络线的平均值mj-1(t)；
S403e：计算hj(t)＝hj-1(t)-mj-1(t)，若hj(t)是IMF函数，则第i个本征函数imfi(t)＝hj(t)，否则j＝j+1，转至步骤S403b；
S404：计算li(t)＝li-1(t)-imfi(t)；
S405：如果li(t)的极值点数仍多于2个，则i＝i+1，转至步骤S403；否则，分解结束；最终，li(t)是残余分量。
进一步，步骤S5中所述的SVD具体包括以下步骤：
S501：求出矩阵X的转置矩阵XT、XT*X或X*XT；
S502：选择矩阵XT*X和X*XT中尺度较小的矩阵，进入步骤S503；
S503：求出矩阵XT*X或X*XT的k个特征值为λi，其中1≤i≤k；
S504：求出λi的算数平方根σi即为矩阵X的k个奇异值；
进一步，根据步骤S6中所述的胎心电信号Fecg和QRS波群的个数计算出胎儿的瞬时心率。
进一步，根据步骤S6中所述的胎心电信号Fecg该胎心电信号分析脐带缠绕、胎儿缺氧的情况。
本发明的有益效果在于：本发明只采用了一路传感器从母体腹壁采集信号，然后滤除肌电等干扰和噪声，从而实现了母胎心电信号的分离。由于胎心电信号是在强噪声中的一种弱信号，所以本发明必须采用高信噪比和灵敏度的生物电传感器。本发明可在孕妇怀孕5个孕周后初次使用，并可以对胎儿个数进行预估。在医院用B超确诊胎儿个数后，可进一步对算法进行修正。实现对孕妇的远程医疗、远程监护和延续护理。本发明推广之后将具有较大的经济价值和社会效益，对远程医疗和移动监护的意义较为重大。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：
图1为本发明示意图；
图2为母体腹部传感器粘贴位置示意图；
图3为便携式胎心电监测仪结构原理图；
图4为EEMD分解出的本征模态信号IMF；
图5为母胎心电信号分离结果。
具体实施方式
下面将结合附图，对本发明的优选实施例进行详细的描述。
如图3所示，一种便携式多胞胎心电监测仪，包括用于采集母体肌电信号的生物电传感器、用于放大信号与滤掉噪声干扰的信号放大与滤波器、用于把模拟信号转换成数字信号的A/D转换器、用于控制和处理信号的嵌入式微控制器、用于存储信息的存储器、用于显示监测数据的LCD显示器、用于提供电源的充电电池、时钟频率模块、用于非接触式传输数据的蓝牙和近距离无线通信技术(Near Field Communication,NFC)模块；
所述生物电传感器将采集到的母体肌电信号通过所述信号放大与滤波器进行滤波和除噪，再通过所述A/D转换器转换成数字信号输入到所述嵌入式微控制器进行信号处理，将监测结果传送到所述LCD显示器，并通过蓝牙和NFC模块进行无线数据传输；所述嵌入式微控制器将信息存储到所述存储器，所述时钟频率模块控制所述嵌入式微控制器的时钟频率。
如图1所示，一种便携式多胞胎心电监测仪的检测方法，该方法包括以下步骤：
S1：采集母体腹部的信号；
S2：采用自适应滤波算法滤除掉肌电信号、噪声干扰；
S3：采用小波算法提取出母体心电信号MECG；
S4：对去除母体心电信号后的残余信号FECGErr，加入高斯白噪声，进行集合经验模态分解(EnsembleEmpiricalModeDecomposition,EEMD)，从而把残余信号分解为若干个本征模态信号(IntrinsicModeFunction,IMF)；采用X表示分解后的IMF信号，即X＝(c1,c2,…,cn,r)，其中ci表示第i个本征模态信号，共n个本征模态信号，r表示分解后的残差；
S5：对X进行奇异值分解(SingularValueDecomposition,SVD)，求出矩阵X*XT的k个特征值，设为λi，其中1≤i≤k，该特征值的算术平方根σi即为矩阵X的k个奇异值；并对该k个奇异值按从大到小的顺序排序为：σ1＞σ2＞…＞σk；
S6：设μi＝σi/σi+1，其中i＝1,2,…,k-1；
求出μi中的最大值，设当i＝j时，μi取得k个奇异值的降低梯度中的最大值μj，即相邻奇异值比值中的最大值μj；
则对胎儿的预测个数为j个，且j个胎心电信号为：
Fecg＝(c1,c2,…,cj-1,cj+cj+1+…cn+r)。
所述步骤S1具体为：将便携式胎心电监测仪的传感器贴于母体肚脐眼正下方三指宽处，采集到一路母体腹部肌电信号。
步骤S4中所述EEMD具体包括以下步骤：
S401：在FecgErr中添加一定量的高斯白噪声；
S402：进行初始化操作：l0(t)＝FecgErr，i＝1；
S403：得到第i个IMF；
S403a：h0(t)＝li-1(t)，j＝1；
S403b：找出hj-1(t)的局部极值点；
S403c：对hj-1(t)的极大和极小值点分别进行三次样条函数插值，分别形成hj-1(t)的上下包络线；
S403d：计算上下包络线的平均值mj-1(t)；
S403e：计算hj(t)＝hj-1(t)-mj-1(t)，若hj(t)是IMF函数，则第i个本征函数imfi(t)＝hj(t)，否则j＝j+1，转至步骤S403b；
S404：计算li(t)＝li-1(t)-imfi(t)；
S405：如果li(t)的极值点数仍多于2个，则i＝i+1，转至步骤S403；否则，分解结束；最终，li(t)是残余分量。
步骤S5中所述的SVD具体包括以下步骤：
S501：求出矩阵X的转置矩阵XT、XT*X或X*XT；
S502：选择矩阵XT*X和X*XT尺度较小的矩阵，进入步骤S503；
S503：求出矩阵XT*X或X*XT的k个特征值为λi，其中1≤i≤k；
S504：求出λi的算数平方根σi即为矩阵X的k个奇异值；
根据步骤S6中所述的胎心电信号Fecg和QRS波群的个数计算出胎儿的瞬时心率。
根据步骤S6中所述的胎心电信号Fecg该胎心电信号分析脐带缠绕、胎儿缺氧的情况。
具体地，首先，将高灵敏度的传感器贴在孕妇腹部如图2所示的④号位置，经过专门的信号放大、调理电路等将采集到的信号输入到嵌入式胎心电监护仪中，并经过A/D模块将其转换为数字信号，以便后续处理。
先进行滤波处理：采用截止频率为10Hz至1kHz的带通滤波器。以便滤出高频噪声和低频肌电信号等干扰信号，由于本发明采用了高品质的锂离子充电电池，所以没有设计50Hz的陷波器，即本发明未考虑工频干扰。
采用小波多尺度分析算法提取出母体心电信号：需要选取合适的小波包及分解层数等。
如图4所示，对剔除母体心电信号之后的残余信号进行EEMD分解，把残余信号分解为若干个本振模态信号IMF，再对IMF进行SVD分解，求出IMF的奇异值，并将这些奇异值按从大到小顺序排序。
再求出奇异值的变化梯度序列，即相邻两个奇异值的比值。找出最大变化梯度值及其在序列中的位置，该位置数即为预估胎儿个数。
再根据预估胎儿个数及IMF信号重构出胎心电信号。
本发明最终提取的胎儿心电信号和母体心电信号如图5所示，并得出结论为该案例中胎儿个数为1，及胎儿心率和母体心率分别为每分钟80次左右和100次左右。并且在胎儿心电信号的基础上，能够进行P波、QRS波(P波：左右心房除极时的电位变化；QRS波：代表心室除极的电位变化)等波型分析，以便得出更多的结论。
最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。
Claims (6)
1.一种便携式胎心电监测仪的检测方法，其特征在于：该方法包括以下步骤：
S1：采集母体腹部的信号；
S2：采用自适应滤波算法滤除掉肌电信号、噪声干扰；
S3：采用小波算法提取出母体心电信号MECG；
S4：对去除母体心电信号后的残余信号FECGErr，加入高斯白噪声，进行集合经验模态分解EEMD，从而把残余信号分解为若干个本征模态信号IMF；采用X表示分解后的IMF信号，即X＝(c1,c2,…,cn,r)，其中ci表示第i个本征模态信号，共n个本征模态信号，r表示分解后的残差；
S5：对X进行奇异值分解SVD，求出矩阵X*XT的k个特征值，设为λi，其中1≤i≤k，该特征值的算术平方根σi即为矩阵X的k个奇异值；并对该k个奇异值按从大到小的顺序排序为：σ1＞σ2＞…＞σk；
S6：设μi＝σi/σi+1，其中i＝1,2,…,k-1；
求出μi中的最大值，设当i＝j时，μi取得k个奇异值的降低梯度中的最大值μj，即相邻奇异值比值中的最大值μj；
则对胎儿的预测个数为j个，且j个胎心电信号为：
Fecg＝(c1,c2,…,cj-1,cj+cj+1+…cn+r)。
2.如权利要求1所述的一种便携式胎心电监测仪的检测方法，其特征在于：所述步骤S1具体为：将便携式胎心电监测仪的传感器贴于母体肚脐眼正下方三指宽处，采集到一路母体腹部肌电信号。
3.如权利要求1所述的一种便携式胎心电监测仪的检测方法，其特征在于：步骤S4中所述EEMD具体包括以下步骤：
S401：在FecgErr中添加一定量的高斯白噪声；
S402：进行初始化操作：l0(t)＝FecgErr，i＝1；
S403：得到第i个IMF；
S403a：h0(t)＝li-1(t)，j＝1；
S403b：找出hj-1(t)的局部极值点；
S403c：对hj-1(t)的极大和极小值点分别进行三次样条函数插值，分别形成hj-1(t)的上下包络线；
S403d：计算上下包络线的平均值mj-1(t)；
S403e：计算hj(t)＝hj-1(t)-mj-1(t)，若hj(t)是IMF函数，则第i个本征函数imfi(t)＝hj(t)，否则j＝j+1，转至步骤S403b；
S404：计算li(t)＝li-1(t)-imfi(t)；
S405：如果li(t)的极值点数仍多于2个，则i＝i+1，转至步骤S403；否则，分解结束；最终，li(t)是残余分量。
4.如权利要求1所述的一种便携式胎心电监测仪的检测方法，其特征在于：步骤S5中所述的SVD具体包括以下步骤：
S501：求出矩阵X的转置矩阵XT、XT*X或X*XT；
S502：选择矩阵XT*X和X*XT中尺度较小的矩阵，进入步骤S503；
S503：求出矩阵XT*X或X*XT的k个特征值为λi，其中1≤i≤k；
S504：求出λi的算数平方根σi即为矩阵X的k个奇异值；
5.如权利要求1所述的一种便携式胎心电监测仪的检测方法，其特征在于：根据步骤S6中所述的胎心电信号Fecg和QRS波群的个数计算出胎儿的瞬时心率。
6.基于权利要求1所述方法的便携式胎心电监测仪，其特征在于：包括用于采集母体肌电信号的生物电传感器、用于放大信号与滤掉噪声干扰的信号放大与滤波器、用于把模拟信号转换成数字信号的A/D转换器、用于控制和处理信号的嵌入式微控制器、用于存储信息的存储器、用于显示监测数据的LCD显示器、用于提供电源的充电电池、时钟频率模块、用于非接触式传输数据的蓝牙和近距离无线通信技术NFC模块；
所述生物电传感器将采集到的母体腹部体表肌电信号通过所述信号放大与滤波器进行滤波和除噪，再通过所述A/D转换器转换成数字信号输入到所述嵌入式微控制器进行信号处理，将监测结果传送到所述LCD显示器，并通过蓝牙和NFC模块进行无线数据传输；所述嵌入式微控制器将信息存储到所述存储器，所述时钟频率模块控制所述嵌入式微控制器的时钟频率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710500220.4A CN107157474B (zh) | 2017-06-27 | 2017-06-27 | 一种便携式多胞胎心电监测仪及其监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710500220.4A CN107157474B (zh) | 2017-06-27 | 2017-06-27 | 一种便携式多胞胎心电监测仪及其监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107157474A CN107157474A (zh) | 2017-09-15 |
| CN107157474B true CN107157474B (zh) | 2020-04-28 |
Family
ID=59826264
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710500220.4A Expired - Fee Related CN107157474B (zh) | 2017-06-27 | 2017-06-27 | 一种便携式多胞胎心电监测仪及其监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107157474B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111557661B (zh) * | 2020-05-15 | 2024-03-08 | 深圳清华大学研究院 | 心电信号处理方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6149602A (en) * | 1996-09-10 | 2000-11-21 | Arcelus; Almudena | User-worn electrocardiogram viewer device |
| US7892178B1 (en) * | 2009-09-28 | 2011-02-22 | Impact Sports Technologies, Inc. | Monitoring device for an interactive game |
| CN103222862A (zh) * | 2013-04-23 | 2013-07-31 | 江苏矽望电子科技有限公司 | 基于声学传感器的被动式胎心监护系统 |
| CN203493632U (zh) * | 2013-04-23 | 2014-03-26 | 江苏矽望电子科技有限公司 | 被动式胎心率测量系统 |
| CN103705270A (zh) * | 2012-09-28 | 2014-04-09 | 清华大学 | 胎心监测设备 |
| CN103876731A (zh) * | 2014-03-25 | 2014-06-25 | 电子科技大学 | 一种胎儿心电信号提取装置及方法 |
-
2017
- 2017-06-27 CN CN201710500220.4A patent/CN107157474B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6149602A (en) * | 1996-09-10 | 2000-11-21 | Arcelus; Almudena | User-worn electrocardiogram viewer device |
| US7892178B1 (en) * | 2009-09-28 | 2011-02-22 | Impact Sports Technologies, Inc. | Monitoring device for an interactive game |
| CN103705270A (zh) * | 2012-09-28 | 2014-04-09 | 清华大学 | 胎心监测设备 |
| CN103222862A (zh) * | 2013-04-23 | 2013-07-31 | 江苏矽望电子科技有限公司 | 基于声学传感器的被动式胎心监护系统 |
| CN203493632U (zh) * | 2013-04-23 | 2014-03-26 | 江苏矽望电子科技有限公司 | 被动式胎心率测量系统 |
| CN103876731A (zh) * | 2014-03-25 | 2014-06-25 | 电子科技大学 | 一种胎儿心电信号提取装置及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107157474A (zh) | 2017-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7333850B2 (en) | Maternal-fetal monitoring system | |
| Abdulhay et al. | Monitoring techniques | |
| US7831302B2 (en) | Monitoring electrical muscular activity | |
| CN104523263B (zh) | 基于移动互联网的孕产妇健康监护系统 | |
| US20190125246A1 (en) | System And Method For Analyzing Progress of Labor And Preterm Labor | |
| EP2451345B1 (en) | Monitoring uterine activity | |
| Alberola-Rubio et al. | Comparison of non-invasive electrohysterographic recording techniques for monitoring uterine dynamics | |
| Fanelli et al. | Telefetalcare: A first prototype of a wearable fetal electrocardiograph | |
| CN109199375A (zh) | 一种无创胎儿心电检测装置及心电信号数据处理方法 | |
| CN107411736A (zh) | 胎儿心电信号检测系统 | |
| CN102160787A (zh) | 一种基于时频变换的胎儿心电盲提取方法 | |
| KR101890072B1 (ko) | 자궁 근전도 신호를 기반으로 임산부의 조산예측기능을 갖는 임산부 모니터링 시스템 | |
| Queyam et al. | An IoT based multi-parameter data acquisition system for efficient bio-telemonitoring of pregnant women at home | |
| Viunytskyi et al. | Fetal ECG and heart rhythm analyzing using BabyCard | |
| CN107157474B (zh) | 一种便携式多胞胎心电监测仪及其监测方法 | |
| CN204500673U (zh) | 基于移动互联网的孕产妇健康监护系统 | |
| Hu et al. | A wireless and wearable system for fetal heart rate monitoring | |
| Fanelli | TeleFetal Care: development of a wearable system for fetal monitoring during pregnancy | |
| Fanelli et al. | A smart wearable prototype for fetal monitoring | |
| Desai et al. | A comparison and quantification of fetal heart rate variability using Doppler ultrasound and direct electrocardiography acquisition techniques | |
| Tiwari et al. | Labour monitoring in pregnant women using electrocardiography and electromyography | |
| Preethi et al. | Design of H-IoT based pregnancy monitoring device in free-living environment | |
| JP7284548B1 (ja) | 信号出力装置 | |
| Ungureanu et al. | Prenatal Telemedicine-Advances in Fetal Monitoring | |
| CN116898418A (zh) | 一种孕期母婴生理指标动态监测方法及监测系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200428 Termination date: 20210627 |