CN104605831A - 非接触式生命体征监护系统的呼吸与心跳信号分离算法 - Google Patents
非接触式生命体征监护系统的呼吸与心跳信号分离算法 Download PDFInfo
- Publication number
- CN104605831A CN104605831A CN201510056983.5A CN201510056983A CN104605831A CN 104605831 A CN104605831 A CN 104605831A CN 201510056983 A CN201510056983 A CN 201510056983A CN 104605831 A CN104605831 A CN 104605831A
- Authority
- CN
- China
- Prior art keywords
- signal
- waveform
- heartbeat
- real
- obtains
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Pulmonology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明提供一种基于非接触式生命体征监护系统的呼吸与心跳信号分离算法,包括:基于反正切处理原理,将模拟接收端接收到的I/Q两路信号进行反正切处理得到正交解调信号;采用基于插值运算的波形重构算法,选取内节点对呼吸信号进行拟合;将正交解调信号中的拟合呼吸信号波形进行剔除从而获得心跳信号波形;对拟合呼吸信号波形和心跳信号波形进行信号平滑以及频率提取操作,得到呼吸和心跳信号的实时频率。本发明方法可有效对正交解调信号中的呼吸信号和心跳信号进行分离提取,并准确获取两者频率信息,方便医务人员准确、实时、连续掌握病人的呼吸和心跳状况。
Description
技术领域
本发明属于雷达领域,特别涉及一种基于非接触式生命体征监护系统的呼吸与心跳信号分离算法。
背景技术
心跳及呼吸是人体重要的身体特征状况信息。目前常用的处理生命体征监护系统的信号处理方法为采用数字滤波器组对呼吸和心跳信号进行分离提取,运用该种方法时,需要极高阶数的滤波器才能满足分离的性能要求,并且滤波器本身具有的相移和衰减也会导致生命体征信号产生形变,因此该法测量精度存在较大的误差。
回波信号中的呼吸成分较心跳成分来说,其幅度大、频率低,即主要变化趋势是由呼吸信号来决定的。因此,采用基于插值运算的重构算法对呼吸信号进行拟合,得到包含实际呼吸大部分信息的拟合呼吸信号波形是可行的。
插值运算并不是试图找到适合于所有自变量数组x的全局最优拟合函数y=f(x),而是希望能够找到一个解析函数连接自变量相邻的两个点(xi,xi+1),并由此找到两点之间的数值。在许多工程问题上,由于只能获得无规律的离散点,插值运算可以帮我们获得近似的连续过程,便于数学的解析方法对已有数据进行处理和运算。
本发明采用插值运算的原理,在复杂的正交解调信号当中选择恰当的内节点,并用曲线将内节点光滑链接起来得到拟合呼吸信号波形,再从正交解调信号中将拟合呼吸信号波形进行剔除,从而得到心跳信号波形,可以有效解决滤波器处理时面临的众多问题
发明内容
本发明依托非接触式生命体征监护系统进行实现。
该系统采用现代雷达技术和生物医学技术相融合的方式,基于多普勒检测原理,运用无线测量的方式实时获取被测目标的体征信息,包括呼吸、心跳等;主要包括雷达信号发射部分、雷达信号模拟接收部分和数字信号处理部分,其中,雷达信号模拟接收部分分别与雷达信号发射部分和数字信号处理部分相连。
发射部分用于通过发射天线向被测目标发出连续波信号并照射在目标的胸腔,连续波信号经过由呼吸和心跳等周期信号引起的胸腔运动对其进行调制后得到回波信号;回波信号经过模拟接收部分解调成为两路包含呼吸和心跳信息的正交数字信号,一路为I路信号,另一路为Q路信号;I路和Q路信号经过数字信号处理部分进行生命体征信号的提取及频率检测。
实现本发明目的的技术解决方案为:
步骤1,系统模拟端送来两路包含呼吸和心跳信息的正交数字信号,一路为I路信号,另一路为Q路信号,对I路和Q路进行反正切处理得到正交解调信号。
步骤2,利用基于插值运算的波形重构算法将步骤1获得的正交解调信号进行波形重构处理得到拟合呼吸信号波形,具体步骤如下:
步骤3,将步骤1中的正交解调信号中剔除步骤2获得的拟合呼吸信号得到心跳信号波形。
步骤4,对步骤2中获得的拟合呼吸信号波形进行傅里叶变换处理得到实时呼吸频率,对步骤3中获得的心跳信号波形进行平滑处理后得到实时心跳波形,再经过傅立叶变换处理得到实时心跳频率。
本发明进一步采用以下技术手段实现发明目的:
(1)步骤2的具体过程为:
步骤2.1,选取内节点间隔n;
步骤2.2,根据步骤2.1中的内节点间隔n,采用插值运算中的三次样条插值算法对正交解调信号进行波形重构,得到拟合呼吸信号波形。
内节点间隔n应满足:其中,n0=1,2,3...,TR与TH分别为呼吸信号周期和心跳信号周期,fs为系统采样频率。
(2)步骤4的具体步骤如下:
步骤4.1,对拟合呼吸信号波形进行FFT运算,得到实时呼吸频率;
步骤4.2,对心跳信号波形进行平滑滤波处理,得到实时心跳波形;
步骤4.3,对实时心跳波形进行FFT运算,得到实时心跳频率。
本发明与现有技术相比,其显著优点:(1)可有效对实际测得呼吸信号和心跳信号进行分离提取,所得计算结果符合实际情况,同时基于本发明实施的双支路正交解调可以有效的解决空点带来的信号频率无法检测问题;(2)本发明分离效果好,计算精度高,实际应用性强,方便医务人员准确、实时、连续掌握病人的心跳及呼吸状况。
下面结合说明书附图对本发明做进一步描述。
附图说明
图1为本发明的非接触式生命体征监护系统的信号处理流程框图;
图2为本发明的实际I、Q两路信号和正交解调信号波形;
图3为本发明的内节点和拟合呼吸信号波形以及心跳信号波形;
图4为本发明的拟合呼吸信号波形及其频谱图;
图5为本发明的实时心跳波形及其频谱图。
具体实施方式
结合图1,步骤1,系统模拟端送来两路包含呼吸和心跳信息的正交数字信号,一路为I路信号(同相支路),另一路为Q路信号(正交支路),对I路和Q路进行反正切处理得到正交解调信号。
步骤2,利用基于插值运算的波形重构算法将步骤1获得的正交解调信号进行波形重构处理得到拟合呼吸信号波形。
步骤3,将步骤1中的正交解调信号中剔除步骤2获得的拟合呼吸信号得到心跳信号波形。
步骤4,对步骤2中获得的拟合呼吸信号波形进行傅里叶变换处理得到实时呼吸频率,对步骤3中获得的心跳信号波形进行平滑处理后得到实时心跳波形,再经过傅立叶变换处理得到实时心跳频率。
步骤2的具体过程为:
步骤2.1,选取内节点间隔n;
步骤2.2,根据步骤2.1中的内节点间隔n,采用插值运算中的三次样条插值算法对正交解调信号进行波形重构,得到拟合呼吸信号波形。
内节点间隔n应满足:其中,n0=1,2,3...,TR与TH分别为呼吸信号周期和心跳信号周期,fs为系统采样频率。
步骤4的具体步骤如下:
步骤4.1,对拟合呼吸信号波形进行FFT运算,得到实时呼吸频率;
步骤4.2,对心跳信号波形进行平滑滤波处理,得到实时心跳波形;
步骤4.3,对实时心跳波形进行FFT运算,得到实时心跳频率。
结合图2,系统模拟端送来两路包含呼吸和心跳信息的正交数字信号,一路为I路信号,另一路为Q路信号,对I路和Q路进行反正切处理得到正交解调信号,经过步骤1后得到正交解调信号如图2中最下方图形所示。
结合图3,利用基于插值运算的波形重构算法将步骤1获得的正交解调信号进行波形重构处理得到拟合呼吸信号波形。经过步骤2.1后得到内节点如图3中星号所示;经过步骤2.2后得到拟合呼吸信号波形如图3中虚线所示。将步骤1中的正交解调信号中剔除步骤2获得的拟合呼吸信号得到心跳信号波形。经过步骤3后得到心跳信号波形如图3中加粗实线所示。
结合图4、图5,对步骤2中获得的拟合呼吸信号波形进行傅里叶变换处理得到实时呼吸频率,对步骤3中获得的心跳信号波形进行平滑处理后得到实时心跳波形,再经过傅立叶变换处理得到实时心跳频率。经过步骤4.1后得到呼吸信号频谱图如图4中最下方图形所示,实时呼吸频率为17.58次/分钟(即0.293Hz);经过步骤4.2后得到心跳信号波形如图5中最上方图形所示;经过步骤4.3后得到心跳信号频谱图如图5中最下方图形所示,实时心跳频率为60.06次/分钟(即1.001Hz)。
Claims (5)
1.一种非接触式生命体征监护系统的呼吸与心跳信号分离算法,其特征在于,包括以下步骤:
步骤1,非接触式生命体征监护系统送来正交的两路心跳和呼吸数字混合信号,一路为I路信号,另一路为Q路信号,对I路和Q路进行反正切处理得到正交解调信号;
步骤2,利用插值运算,将步骤1获得的正交解调信号进行波形重构后得到拟合呼吸信号波形;
步骤3,将步骤1中的正交解调信号中剔除步骤2获得的拟合呼吸信号得到心跳信号波形;
步骤4,对步骤2获得的拟合呼吸信号波形进行傅里叶变换处理得到实时呼吸频率,对步骤3获得的心跳信号波形进行平滑处理后得到实时心跳波形,在经过傅立叶变换处理得到实时心跳频率。
2.根据权利要求1所述的非接触式生命体征监护系统的呼吸与心跳信号分离算法,其特征在于,步骤2的具体过程如下:
步骤2.1,选取内节点间隔n;
步骤2.2,根据步骤2.1中得到的内节点间隔n,采用三次样条插值算法对正交解调信号进行波形重构,得到拟合呼吸信号波形。
3.根据权利要求2所述的非接触式生命体征监护系统的呼吸与心跳信号分离算法,其特征在于:根据奈奎斯特采样定理,内节点间隔n应满足
其中,n0=1,2,3...,TR与TH分别为呼吸信号周期和心跳信号周期,fs为系统采样频率。
4.根据权利要求1所述的非接触式生命体征监护系统的呼吸与心跳信号分离算法,其特征在于,步骤3中心跳信号获取包括如下步骤:将正交解调信号与步骤2中得到的拟合呼吸信号波形做差,得到心跳信号波形。
5.根据权利要求1所述的非接触式生命体征监护系统的呼吸与心跳信号分离算法,其特征在于,步骤4实时频率获取包括如下步骤:
步骤4.1,对拟合呼吸信号波形进行FFT运算,得到实时呼吸频率;
步骤4.2,对心跳信号波形进行平滑滤波处理,得到实时心跳波形;
步骤4.3,对实时心跳波形进行FFT运算,得到实时心跳频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510056983.5A CN104605831A (zh) | 2015-02-03 | 2015-02-03 | 非接触式生命体征监护系统的呼吸与心跳信号分离算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510056983.5A CN104605831A (zh) | 2015-02-03 | 2015-02-03 | 非接触式生命体征监护系统的呼吸与心跳信号分离算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104605831A true CN104605831A (zh) | 2015-05-13 |
Family
ID=53140692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510056983.5A Pending CN104605831A (zh) | 2015-02-03 | 2015-02-03 | 非接触式生命体征监护系统的呼吸与心跳信号分离算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104605831A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017169870A (ja) * | 2016-03-24 | 2017-09-28 | 新日本無線株式会社 | 心肺機能測定装置 |
CN107981841A (zh) * | 2017-10-27 | 2018-05-04 | 深圳和而泰智能控制股份有限公司 | 一种信号处理方法、装置、设备及介质 |
CN109199348A (zh) * | 2018-09-13 | 2019-01-15 | 芜湖博高光电科技股份有限公司 | 一种智能生命体征探测装置 |
CN109316173A (zh) * | 2018-10-10 | 2019-02-12 | 大连理工大学 | 一种非接触式生命体征无线监测装置及方法 |
CN109875529A (zh) * | 2019-01-23 | 2019-06-14 | 北京邮电大学 | 一种基于超宽带雷达的生命体征检测方法及系统 |
CN110596705A (zh) * | 2019-08-22 | 2019-12-20 | 南京理工大学 | 基于生命体征sar成像的人体目标身份识别方法及系统 |
CN111685741A (zh) * | 2020-06-11 | 2020-09-22 | 中山大学 | 基于正交解调脉冲超宽带雷达探测人体呼吸率心率的方法 |
CN112716462A (zh) * | 2020-12-09 | 2021-04-30 | 北京航空航天大学 | 一种可控照射指向的窄波束毫米波人体心跳/呼吸体征监测装置 |
CN114098679A (zh) * | 2021-12-30 | 2022-03-01 | 中新国际联合研究院 | 基于深度学习和射频感知的生命体征监测波形的恢复方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4958638A (en) * | 1988-06-30 | 1990-09-25 | Georgia Tech Research Corporation | Non-contact vital signs monitor |
US20110208060A1 (en) * | 2010-02-24 | 2011-08-25 | Haase Wayne C | Non-contact Biometric Monitor |
-
2015
- 2015-02-03 CN CN201510056983.5A patent/CN104605831A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4958638A (en) * | 1988-06-30 | 1990-09-25 | Georgia Tech Research Corporation | Non-contact vital signs monitor |
US20110208060A1 (en) * | 2010-02-24 | 2011-08-25 | Haase Wayne C | Non-contact Biometric Monitor |
Non-Patent Citations (1)
Title |
---|
岳宇: "生物雷达检测技术中心跳和呼吸信号分离技术的研究", 《 中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017169870A (ja) * | 2016-03-24 | 2017-09-28 | 新日本無線株式会社 | 心肺機能測定装置 |
CN107981841A (zh) * | 2017-10-27 | 2018-05-04 | 深圳和而泰智能控制股份有限公司 | 一种信号处理方法、装置、设备及介质 |
CN109199348A (zh) * | 2018-09-13 | 2019-01-15 | 芜湖博高光电科技股份有限公司 | 一种智能生命体征探测装置 |
CN109316173A (zh) * | 2018-10-10 | 2019-02-12 | 大连理工大学 | 一种非接触式生命体征无线监测装置及方法 |
CN109875529A (zh) * | 2019-01-23 | 2019-06-14 | 北京邮电大学 | 一种基于超宽带雷达的生命体征检测方法及系统 |
CN110596705A (zh) * | 2019-08-22 | 2019-12-20 | 南京理工大学 | 基于生命体征sar成像的人体目标身份识别方法及系统 |
CN111685741A (zh) * | 2020-06-11 | 2020-09-22 | 中山大学 | 基于正交解调脉冲超宽带雷达探测人体呼吸率心率的方法 |
CN111685741B (zh) * | 2020-06-11 | 2021-06-08 | 中山大学 | 基于正交解调脉冲超宽带雷达探测人体呼吸率心率的方法 |
CN112716462A (zh) * | 2020-12-09 | 2021-04-30 | 北京航空航天大学 | 一种可控照射指向的窄波束毫米波人体心跳/呼吸体征监测装置 |
CN112716462B (zh) * | 2020-12-09 | 2022-04-05 | 北京航空航天大学 | 一种可控照射指向的窄波束毫米波人体心跳/呼吸体征监测装置 |
CN114098679A (zh) * | 2021-12-30 | 2022-03-01 | 中新国际联合研究院 | 基于深度学习和射频感知的生命体征监测波形的恢复方法 |
CN114098679B (zh) * | 2021-12-30 | 2024-03-29 | 中新国际联合研究院 | 基于深度学习和射频感知的生命体征监测波形的恢复方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104605831A (zh) | 非接触式生命体征监护系统的呼吸与心跳信号分离算法 | |
CN104644142B (zh) | 一种非接触式生命体征监护的信号处理算法 | |
CN110584631A (zh) | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 | |
CN105962914B (zh) | 基于盲源分离的呼吸与心跳信号的分离方法及装置 | |
CN106264502A (zh) | 一种非接触式生理信号检测方法 | |
CN103070728A (zh) | 一种非接触式生命体征监护设备 | |
WO2016057781A1 (en) | Method and apparatus for non-contact fast vital sign acquisition based on radar signal | |
WO2020072297A1 (en) | Direct rf signal processing for heart-rate monitoring using uwb impulse radar | |
CN105232026A (zh) | 一种非接触式生命体征检测系统的心跳频率检测算法 | |
CN106805940A (zh) | 一种连续波生物雷达体征检测装置 | |
CN110464320A (zh) | 多目标人体心率和呼吸频率测量系统及方法 | |
CN103070687A (zh) | 一种非接触式生命体征监护系统的信号处理算法 | |
CN109330597B (zh) | 一种基于信道状态信息的人体呼吸追踪方法 | |
CN108852327B (zh) | 一种从运动干扰中非接触检测微弱生命信号的方法 | |
CN106175731A (zh) | 非接触式生命体征监测的信号处理系统 | |
CN111513706B (zh) | 一种针对含有异常r波的心电信号的检测方法和装置 | |
CN106805941A (zh) | 一种连续波生物雷达体征检测装置 | |
CN111358464A (zh) | 一种针对卧床病人的非接触式生命体征监测方法 | |
CN206239402U (zh) | 一种基于dsp的汽车疲劳驾驶实时检测系统 | |
CN115399747B (zh) | 一种基于毫米波雷达的多人生命体征监测方法 | |
CN115644840A (zh) | 基于毫米波雷达的生命体征检测方法 | |
Lee et al. | Doppler radar in respiratory monitoring: Detection and analysis | |
CN105395214B (zh) | 基于rf数据超声成像处理方法及系统 | |
KR20210066332A (ko) | 타겟의 생체 정보 결정 방법 및 장치 | |
Rong et al. | Novel Respiration-Free Heartbeat Detection Algorithm Using Millimeter-Wave Radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150513 |
|
WD01 | Invention patent application deemed withdrawn after publication |