CN111685760B - 一种基于雷达测量的人体呼吸频率计算方法 - Google Patents
一种基于雷达测量的人体呼吸频率计算方法 Download PDFInfo
- Publication number
- CN111685760B CN111685760B CN202010706217.XA CN202010706217A CN111685760B CN 111685760 B CN111685760 B CN 111685760B CN 202010706217 A CN202010706217 A CN 202010706217A CN 111685760 B CN111685760 B CN 111685760B
- Authority
- CN
- China
- Prior art keywords
- human body
- calculating
- frequency
- echo sequence
- respiratory frequency
- 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.)
- Active
Links
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/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- 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/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明涉及雷达信号处理技术,其公开了一种基于雷达测量的人体呼吸频率计算方法,实时、准确地计算呼吸频率,计算复杂度低。该方法包括以下步骤:a.接收使用超宽带雷达向目标人体所处自由空间发射电磁波而反射回来的回波信号;b.对回波信号的距离维信息进行傅里叶变换,然后进行极大值搜索,提取以极大值为中心的预设距离内的数据,得到人体的回波序列;c.通过人体回波序列的频谱图计算参考呼吸频率;d.计算人体回波序列平均幅度差,并结合参考呼吸频率计算准确的呼吸频率。
Description
技术领域
本发明涉及雷达信号处理技术,具体涉及一种基于雷达测量的人体呼吸频率计算方法。
背景技术
呼吸是人体维持生存最基础的生命活动。呼吸频率包含有重要的生理信息,呼吸频率的测量在救援活体探测、医疗病患监护、心肺功能观察、运动效果评估、睡眠质量监测等领域都有广泛应用。市面上的家用产品只能测量心率,不能同时检测呼吸,而医用产品测量呼吸需佩戴检测装置,舒适度不高。
而采用非接触式测量技术来测量呼吸频率可以提高用户的舒适度,超宽带雷达是一种很好的选择。超宽带雷达以其非接触、远距离、可穿透的特点,其在呼吸监测领域受到越来越多的关注。
人体在呼吸时,胸部表面会随着产生周期性的起伏,超宽带雷达能够捕获这些起伏的回波信号。但是超宽带雷达测量呼吸容易受到人体体动、房间其它物体等多种干扰,因此需要一种具有实时性、准确性、低计算复杂度的算法来计算呼吸频率。
发明内容
本发明所要解决的技术问题是:提出一种基于雷达测量的人体呼吸频率计算方法,实时、准确地计算呼吸频率,计算复杂度低。
本发明解决上述技术问题采用的技术方案是:
一种基于雷达测量的人体呼吸频率计算方法,包括以下步骤:
a.接收使用超宽带雷达向目标人体所处自由空间发射电磁波而反射回来的回波信号;
b.对回波信号的距离维信息进行傅里叶变换,然后进行极大值搜索,提取以极大值为中心的预设距离内的数据,得到人体的回波序列;
c.通过人体回波序列的频谱图计算参考呼吸频率;
d.计算人体回波序列平均幅度差,并结合参考呼吸频率计算准确的呼吸频率。
作为进一步优化,步骤c具体包括:
c1.对人体回波序列进行平滑处理;
c2.计算平滑处理后的人体回波序列的频谱图,在频谱图中找出频率范围0.1-0.6Hz的所有峰值,取最大峰对应的频率值,记为fa。
作为进一步优化,步骤d具体包括:
d1.对平滑处理后的人体回波序列做趋势消除处理;
d2.计算趋势消除后人体回波序列的平均幅度差,找出频率范围0.1-0.6Hz的所有谷点;
d3.计算准确的呼吸频率值:
若找出的谷点只有1个,记该点对应的频率值为fb,则呼吸频率=60*fb bpm;
若找出的谷点数大于1个,记这些谷点对应的频率值分别为fb1、fb2、…fbn,分别计算与fa的差的绝对值,取令该绝对值最小时的频率值fbi,则呼吸频率=60*fbi bpm。
本发明的有益效果是:
大幅降低了对数据长度的要求,最少包含两个完整的呼吸周期即可计算呼吸频率,从而保证了计算的实时性;此外,本发明在频谱分析提供的呼吸频率基础上,采用平均幅度差函数得到更高精度的呼吸频率,避免了倍频或分频错误,从而提升了计算的准确性;该方法计算复杂度低,存储空间需求小。
附图说明
图1为本发明中的基于雷达测量的人体呼吸频率计算方法流程图;
图2为接收的人体回波序列示意图;
图3为人体回波序列频谱图;
图4为人体回波序列的平均幅度差示意图。
具体实施方式
本发明旨在提出一种基于雷达测量的人体呼吸频率计算方法,实时、准确地计算呼吸频率,计算复杂度低。如图1所示,其实现步骤包括:a.接收使用超宽带雷达向目标人体所处自由空间发射电磁波而反射回来的回波信号;b.对回波信号的距离维信息进行傅里叶变换,然后进行极大值搜索,提取以极大值为中心的预设距离内的数据,得到人体的回波序列;c.通过人体回波序列的频谱图计算参考呼吸频率;d.计算人体回波序列平均幅度差,并结合参考呼吸频率计算准确的呼吸频率。
实施例:
本实施例中,被试人员选取舒适姿势坐着,避免剧烈的身体运动而影响检测,保持正常呼吸,在距离被试约1米距离安放超宽带雷达,使用发射机向目标人体所处自由空间发射电磁波,通过接收机接收目标人体反射回来的回波信号,其中采样频率为80Hz;
接收到回波信号后,对回波信号的距离维信息进行傅里叶变换,然后进行极大值搜索,提取以极大值为中心的预设距离内的数据,得到人体的回波序列,如图2所示;
接着,我们对人体回波序列进行窗长为80个点的平滑,计算平滑后人体回波序列的频谱图,如图3所示;在频谱图中找出频率范围0.1-0.6Hz的所有峰值,取最大峰对应的频率值,记为fa;
接下来,我们通过计算人体回波序列平均幅度差,并结合参考呼吸频率计算准确的呼吸频率,具体为:
对平滑后的人体回波序列做趋势消除;
计算趋势消除后人体回波序列的平均幅度差,如图4所示;然后,找出频率范围0.1-0.6Hz的所有谷点:
若找出的谷点只有1个,记该点对应的频率值为fb,则呼吸频率=60*fb bpm;
若找出的谷点数大于1个,记这些谷点对应的频率值分别为fb1、fb2、…fbn,分别计算与fa的差的绝对值,取令该绝对值最小时的fbi,则呼吸频率=60*fbi bpm。
Claims (1)
1.一种基于雷达测量的人体呼吸频率计算方法,其特征在于,包括以下步骤:
a.接收使用超宽带雷达向目标人体所处自由空间发射电磁波而反射回来的回波信号;
b.对回波信号的距离维信息进行傅里叶变换,然后进行极大值搜索,提取以极大值为中心的预设距离内的数据,得到人体的回波序列;
c.通过人体回波序列的频谱图计算参考呼吸频率;
d.计算人体回波序列平均幅度差,并结合参考呼吸频率计算准确的呼吸频率;
步骤c具体包括:
c1.对人体回波序列进行平滑处理;
c2.计算平滑处理后的人体回波序列的频谱图,在频谱图中找出频率范围0.1-0.6Hz的所有峰值,取最大峰对应的频率值,记为fa;
步骤d具体包括:
d1.对平滑处理后的人体回波序列做趋势消除处理;
d2.计算趋势消除后人体回波序列的平均幅度差,找出频率范围0.1-0.6Hz的所有谷点;
d3.计算准确的呼吸频率值:
若找出的谷点只有1个,记该点对应的频率值为fb,则呼吸频率=60*fb bpm;
若找出的谷点数大于1个,记这些谷点对应的频率值分别为fb1、fb2、…fbn,分别计算与fa的差的绝对值,取令该绝对值最小时的频率值fbi,则呼吸频率=60*fbi bpm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010706217.XA CN111685760B (zh) | 2020-07-21 | 2020-07-21 | 一种基于雷达测量的人体呼吸频率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010706217.XA CN111685760B (zh) | 2020-07-21 | 2020-07-21 | 一种基于雷达测量的人体呼吸频率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111685760A CN111685760A (zh) | 2020-09-22 |
CN111685760B true CN111685760B (zh) | 2022-06-03 |
Family
ID=72486483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010706217.XA Active CN111685760B (zh) | 2020-07-21 | 2020-07-21 | 一种基于雷达测量的人体呼吸频率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111685760B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113440120B (zh) * | 2021-06-17 | 2022-10-28 | 北京航空航天大学 | 一种基于毫米波雷达的人员呼吸心跳检测方法 |
CN113721213B (zh) * | 2021-07-26 | 2023-10-10 | 森思泰克河北科技有限公司 | 生命体检测方法、终端及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183526A (zh) * | 2006-11-14 | 2008-05-21 | 中兴通讯股份有限公司 | 一种检测语音信号基音周期的方法 |
KR100905102B1 (ko) * | 2008-01-12 | 2009-06-30 | 주식회사 엑스엘 | 비접촉식 생체신호 측정 장치 |
CN101572089A (zh) * | 2009-05-21 | 2009-11-04 | 华为技术有限公司 | 信号周期检测方法及其装置 |
CN108577815A (zh) * | 2018-04-03 | 2018-09-28 | 电子科技大学 | 一种基于超宽带雷达的人体呼吸速率和心率的测定方法 |
CN110464320A (zh) * | 2019-08-22 | 2019-11-19 | 四川长虹电器股份有限公司 | 多目标人体心率和呼吸频率测量系统及方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2796189B1 (fr) * | 1999-07-05 | 2001-10-05 | Matra Nortel Communications | Procedes et dispositifs de codage et de decodage audio |
-
2020
- 2020-07-21 CN CN202010706217.XA patent/CN111685760B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183526A (zh) * | 2006-11-14 | 2008-05-21 | 中兴通讯股份有限公司 | 一种检测语音信号基音周期的方法 |
KR100905102B1 (ko) * | 2008-01-12 | 2009-06-30 | 주식회사 엑스엘 | 비접촉식 생체신호 측정 장치 |
CN101572089A (zh) * | 2009-05-21 | 2009-11-04 | 华为技术有限公司 | 信号周期检测方法及其装置 |
CN108577815A (zh) * | 2018-04-03 | 2018-09-28 | 电子科技大学 | 一种基于超宽带雷达的人体呼吸速率和心率的测定方法 |
CN110464320A (zh) * | 2019-08-22 | 2019-11-19 | 四川长虹电器股份有限公司 | 多目标人体心率和呼吸频率测量系统及方法 |
Non-Patent Citations (2)
Title |
---|
Contact-free measurement of respiratory rate using infrared and vibration sensors;Fatih Erden等;《Infrared Physics&Technology》;20150914;第73卷;第90-91页 * |
一种改进的AMDF求取语音基音的方法;成新民等;《微电子学与计算机》;20051130;第22卷(第11期);第162-167页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111685760A (zh) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5596670B2 (ja) | 患者の非接触呼吸モニタリング | |
CN109875529A (zh) | 一种基于超宽带雷达的生命体征检测方法及系统 | |
Yamamoto et al. | ECG signal reconstruction via Doppler sensor by hybrid deep learning model with CNN and LSTM | |
US20200121207A1 (en) | Method of processing a signal representing a physiological rhythm | |
JP2015027550A (ja) | 心肺パラメータモニタ用装置、心肺パラメータモニタ用システムおよび非一時的な記録媒体 | |
CN103385702A (zh) | 一种无创血压连续检测装置及方法 | |
CN111685760B (zh) | 一种基于雷达测量的人体呼吸频率计算方法 | |
CN108577815A (zh) | 一种基于超宽带雷达的人体呼吸速率和心率的测定方法 | |
CN112244794B (zh) | 基于周期性特征的生命体征检测方法、装置和存储介质 | |
CN110520935A (zh) | 从无线电信号学习睡眠阶段 | |
CN113854992A (zh) | 基于77GHz毫米雷达的非接触式精确心率检测方法 | |
Shi et al. | Neural network based real-time heart sound monitor using a wireless wearable wrist sensor | |
Liu et al. | Human biometric signals monitoring based on wifi channel state information using deep learning | |
CN112363139A (zh) | 基于幅度特征的人体呼吸时长检测方法、装置和存储介质 | |
Okano et al. | Multimodal cardiovascular information monitor using piezoelectric transducers for wearable healthcare | |
Hu et al. | Remote vital signs measurement of indoor walking persons using mm-Wave FMCW radar | |
JP2015217143A (ja) | 心拍計測装置 | |
CN112450900B (zh) | 一种基于智能音箱的非接触式心跳检测方法 | |
Taheri et al. | Non-invasive breathing rate detection using a very low power ultra-wide-band radar | |
CN116172539A (zh) | 基于机器学习的生命体征检测方法、系统、设备及介质 | |
CN114098679B (zh) | 基于深度学习和射频感知的生命体征监测波形的恢复方法 | |
Zhao et al. | T-HSER: Transformer Network Enabling Heart Sound Envelope Signal Reconstruction Based on Low Sampling Rate Millimeter Wave Radar | |
Kawasaki et al. | Continuous blood pressure estimation using millimeter wave radar | |
Kawasaki et al. | Continuous blood pressure monitoring with MMW radar sensor | |
KR20230077607A (ko) | 시간적 위상 코히어런시에 기반한 심박 측정 거리 결정 장치 및 방법 |
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 |