CN114469025A - 一种用于fmcw雷达的多目标呼吸和心跳信号检测方法 - Google Patents
一种用于fmcw雷达的多目标呼吸和心跳信号检测方法 Download PDFInfo
- Publication number
- CN114469025A CN114469025A CN202111631545.9A CN202111631545A CN114469025A CN 114469025 A CN114469025 A CN 114469025A CN 202111631545 A CN202111631545 A CN 202111631545A CN 114469025 A CN114469025 A CN 114469025A
- Authority
- CN
- China
- Prior art keywords
- target
- distance
- equal
- targets
- matrix
- 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.)
- Granted
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/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
-
- 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/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
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/746—Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- 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/024—Detecting, measuring or recording pulse rate or heart rate
-
- 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
-
- 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
本发明公开一种用于FMCW雷达的多目标呼吸和心跳信号检测方法,包括以下步骤:对雷达中频信号,按照帧周期构建矩阵R(m,n),对R(m,n)中每帧数据分别进行FFT,得到时间距离矩阵D(p,n);对D(p,n)的前N1列数据利用恒虚警率进行检测,得到多个目标的潜在距离;利用K‑means算法对N1列中的多个目标潜在距离聚类,得到K个目标的粗略距离,并确定K个目标距离范围;对于D(p,n)中每一列数据,在K个距离范围内搜索峰值,得到K个目标距离值,并计算K个目标位置的相位;对N列数据中K个目标相位进行解绕和差分,进行滤波分离,得到K个目标的呼吸和心跳信号。本发明可克服杂波影响,快速得到多个目标的粗略位置;改善呼吸和心跳信号的检测精度,适用于复杂环境下多目标生命体征的检测。
Description
技术领域
本发明涉及雷达目标检测与信号处理领域,尤其涉及一种用于FMCW雷达的多目标呼吸和心跳信号检测方法。
背景技术
呼吸和心跳是人体重要的生命体征信号,目前常用的检测方式都为接触式检测方式,在对重症感染病人和重度烧伤患者体征信号检测时有很大的局限性。基于雷达的非接触式生命体征检测技术具有远距离、穿透性和非接触等特点,在医疗监护领域有广阔的应用前景。应用于生命体征信号检测领域的雷达体制主要有三种,分别为连续波(CW)、超宽带(UWB)以及线性调频连续波(FMCW)。FMCW雷达具有UWB雷达的测距能力和CW雷达高灵敏度的特点,而且具有体积小、重量轻、功耗低的优势,因此,FMCW毫米波雷达是非接触式生命体征检测系统的优先选择。
FMCW雷达利用人体胸腔运动产生的微多普勒效应,通过目标距离单元信号相位的变化来检测呼吸和心跳信号。在进行多目标的呼吸和心跳检测时,由于不同距离的目标回波强度存在差异,而且目标回波还受到相近距离单元内的其他物体回波和人体自身抖动信号的干扰,因此,如何准确检测目标距离,对于提高多目标呼吸和心跳信号检测准确性至关重要。
发明内容
针对现有技术存在的不足,本发明提供了一种用于FMCW雷达的多目标呼吸和心跳信号检测方法。该方法对部分时间距离帧数据进行恒虚警率检测和K-means距离聚类,可克服干扰杂波影响,快速得到多个目标的粗略距离;在每个时间距离帧数据的目标距离范围内搜索可进一步提高目标距离估计精度,有利于改善呼吸和心跳信号检测精度。
为了解决上述技术问题,本发明采用的技术方案步骤如下:
S101:利用FMCW雷达接收雷达中频信号,并对中频信号按照帧周期采样构建信号矩阵R(m,n),其中1≤m≤M、1≤n≤N,M为每帧采样的点数,N为帧数;
S102:对信号矩阵R(m,n)的每帧或每列数据分别进行FFT变换,得到时间距离矩阵D(p,n),其中1≤p≤P,P为FFT点数;
S103:选取时间距离矩阵D(p,n)的前N1列数据,利用恒虚警率方法进行检测,得到目标潜在距离矩阵S(i,n1),其中1≤i≤max(L(n1)),1≤n1≤N1,L(n1)为第n1列中检测出的目标数量;
S104:利用K-means算法对目标潜在距离矩阵S(i,n1)中目标距离值进行聚类,得到K个目标的粗略距离S1(k),其中1≤k≤K,K为先验目标数量,并确定K个目标距离范围为Sd(k)=[S1(k)-ΔR/2,S1(k)+ΔR/2],ΔR为距离范围大小;
进一步地,步骤S103中,利用恒虚警率方法进行检测,具体如下:
S201:对时间距离矩阵D(p,n)的前N1列数据取其模值的平方,得到D1(p,n1)=|D(p,n1)|2;
S202:选择CA-CFAR检测器作为恒虚警率检测器,将D1(p,n1)中第1列数据作为检测器的输入,得到该列中目标数量L(1),并记录下目标距离CS1(i)(1≤i≤L(1)),遍历N1列数据,得到目标潜在距离矩阵:
进一步地,步骤S104中利用K-means算法对目标潜在距离矩阵S(i,n1)中目标距离值进行聚类,具体如下:
S302:选择欧式距离作为类间距离度量标准,利用K-means算法对ST(ki)进行聚类,得到K个目标的粗略距离S1(k)。
进一步地,步骤S105的具体过程如下:
S401:对D(p,n)每列数据取模,得到D2(p,n)=|D(p,n)|;
步骤S106中,得到K个目标的呼吸和心跳信号的具体步骤如下:
与现有技术相比,本发明的有益效果包括:
1.利用恒虚警率检测技术和K-means聚类算法相结合,对少量时间距离帧数据进行处理,可有效克服目标附近物体和目标其它部位产生的干扰杂波,快速获取多个目标的粗略距离;
2.利用多个目标粗略距离构建多个距离范围,在多个距离范围内搜索可提高多目标距离估计精度,改善多目标呼吸和心跳信号的检测精度。
附图说明
图1是本发明提供的一种用于FMCW雷达的多目标呼吸和心跳信号检测方法流程示意图;
图2为原始雷达中频信号矩阵;
图3为雷达中频信号的时间距离矩阵;
图4为利用恒虚警率技术对时间距离矩阵的前100列数据进行检测的结果;
图5为检测出的两个目标距离处的相位差分信号;
图6为目标1的呼吸和心跳信号;
图7为目标2的呼吸和心跳信号。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供了一种用于FMCW雷达的多目标呼吸和心跳信号检测方法。请参考图1,图1是本发明方法的流程示意图;方法包括以下步骤:
S101:利用FMCW雷达接收雷达中频信号,并对中频信号按照帧周期采样构建信号矩阵R(m,n),其中1≤m≤M、1≤n≤N,M为每帧采样的点数,N为帧数;
请参考图2,图2是原始雷达中频信号矩阵R(m,n);以该矩阵来说具体说明本发明方案。
两个测试目标分别位于FMCW雷达正前方0.7m和1.3m处,利用FMCW雷达进行测试,其中每帧采样点数为100,采样率为2MHz,对应的快时间长度为50us,帧周期为0.05s,帧数为400,对应的慢时间长度为20s,由此得到100×400的原始中频信号矩阵R(m,n),如图2所示。
S102:对信号矩阵R(m,n)的每帧或每列数据分别进行FFT变换,得到时间距离矩阵D(p,n),其中1≤p≤P,P为FFT点数;
具体的,在本实施例中,对矩阵R(m,n)中每列分别进行128点FFT,得到128×400的时间距离矩阵D(m,n),请参考图3,图3是雷达中频信号的时间距离矩阵;
S103:选取时间距离矩阵D(p,n)的前N1列数据,利用恒虚警率方法进行检测,得到目标潜在距离矩阵S(i,n1),其中1≤i≤max(L(n1)),1≤n1≤N1,L(n1)为第n1列中检测出的目标数量;
步骤S103中,利用恒虚警率方法进行检测,具体如下:
S201:对时间距离矩阵D(p,n)的前N1列数据取其模值的平方,得到D1(p,n1)=|D(p,n1)|2;
S202:选择CA-CFAR检测器作为恒虚警率检测器,将D1(p,n1)中第1列数据作为检测器的输入,得到该列中目标数量L(1),并记录下目标距离CS1(i)(1≤i≤L(1)),遍历N1列数据,得到目标潜在距离矩阵:
具体的,在本实施例中,选取D(p,n)的前100列数据,分别利用(恒虚警率检测器)CA-CFAR技术进行检测,结果如图4所示,其中第35列数据的检测示意图如图4(a)所示,该列中检测的两个目标距离分别为0.74m和1.31m;100列数据的目标潜在距离值如图4(b)所示;
S104:利用K-means算法对目标潜在距离矩阵S(i,n1)中目标距离值进行聚类,得到K个目标的粗略距离S1(k),其中1≤k≤K,K为先验目标数量,并确定K个目标距离范围为Sd(k)=[S1(k)-ΔR/2,S1(k)+ΔR/2],ΔR为距离范围大小;
步骤S104中利用K-means算法对目标潜在距离矩阵S(i,n1)中目标距离值进行聚类,具体如下:
S302:选择欧式距离作为类间距离度量标准,利用K-means算法对ST(ki)进行聚类,得到K个目标的粗略距离S1(k)。
具体的,在本实施例中,利用K-means算法对目标潜在距离值进行聚类,得到两个目标粗略距离分别为0.74m和1.34m,并确定两个目标距离范围分别为[0.57m 0.91m]和[1.17m 1.51m],其中距离范围大小为0.34m。
步骤S105的具体过程如下:
S401:对D(p,n)每列数据取模,得到D2(p,n)=|D(p,n)|;
具体的,在本实施例中,对D(p,n)的每列数据,在两个目标距离范围内搜索峰值,并计算两个峰值处的相位值;
步骤S106中,得到K个目标的呼吸和心跳信号的具体步骤如下:
具体的,在本实施例中,对两个峰值处的相位进行解卷绕和差分,得到的相位差分信号如图5所示,其中图5(a)为目标1(0.74m附近)的相位差分信号,图5(b)为目标2(1.34m附近)的相位差分信号;
设计两个IIR带通滤波器,通带分别[0.1Hz 0.7Hz]和[0.9Hz 2.0Hz],分别对两个目标的相位差分信号进行滤波,得到目标1的呼吸信号和心跳信号分别如图6(a)和6(b)所示,目标2的呼吸信号和心跳信号分别如图7(a)和7(b)所示。
在利用FMCW雷达检测过程中,同时利用接触式脉搏传感器测量两个目标的心跳频率,并通过人工计数方式测量两个目标的呼吸频率,作为参考呼吸和心跳频率。
参考值和雷达检测值如表1所示。从表1中可看出,本发明方法测量结果与参考值之间误差较小,具有较高的检测精度。
表1雷达检测值与参考值
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本申请实施例的实施过程构成任何限定。
本发明的有益效果是:
1.利用恒虚警率检测技术和K-means聚类算法相结合,对少量时间距离帧数据进行处理,可有效克服目标附近物体和目标其它部位产生的干扰杂波,快速获取多个目标的粗略距离;
2.利用多个目标粗略距离构建多个距离范围,在多个距离范围内搜索可提高多目标距离估计精度,改善多目标呼吸和心跳信号的检测精度。
以上所述本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保护范围内。
Claims (5)
1.一种用于FMCW雷达的多目标呼吸和心跳信号检测方法,其特征在于:包括以下步骤:
S101:利用FMCW雷达接收雷达中频信号,并对中频信号按照帧周期采样构建信号矩阵R(m,n),其中1≤m≤M、1≤n≤N,M为每帧采样的点数,N为帧数;
S102:对信号矩阵R(m,n)的每帧或每列数据分别进行FFT变换,得到时间距离矩阵D(p,n),其中1≤p≤P,P为FFT点数;
S103:选取时间距离矩阵D(p,n)的前N1列数据,利用恒虚警率方法进行检测,得到目标潜在距离矩阵S(i,n1),其中1≤i≤max(L(n1)),1≤n1≤N1,L(n1)为第n1列中检测出的目标数量;
S104:利用K-means算法对目标潜在距离矩阵S(i,n1)中目标距离值进行聚类,得到K个目标的粗略距离S1(k),其中1≤k≤K,K为先验目标数量,并确定K个目标距离范围为Sd(k)=[S1(k)-ΔR/2,S1(k)+ΔR/2],ΔR为距离范围大小;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111631545.9A CN114469025B (zh) | 2021-12-28 | 2021-12-28 | 一种用于fmcw雷达的多目标呼吸和心跳信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111631545.9A CN114469025B (zh) | 2021-12-28 | 2021-12-28 | 一种用于fmcw雷达的多目标呼吸和心跳信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114469025A true CN114469025A (zh) | 2022-05-13 |
CN114469025B CN114469025B (zh) | 2023-05-26 |
Family
ID=81496622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111631545.9A Active CN114469025B (zh) | 2021-12-28 | 2021-12-28 | 一种用于fmcw雷达的多目标呼吸和心跳信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114469025B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106175723A (zh) * | 2016-06-27 | 2016-12-07 | 中国人民解放军第三军医大学第附属医院 | 一种基于fmcw宽带雷达的多生命监护系统 |
CN108919247A (zh) * | 2018-03-30 | 2018-11-30 | 中国科学院电子学研究所 | 一种基于恒虚警率检测的多目标体检测及定位方法 |
CN110584631A (zh) * | 2019-10-10 | 2019-12-20 | 重庆邮电大学 | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 |
CN113261925A (zh) * | 2021-03-11 | 2021-08-17 | 山东师范大学 | 一种基于毫米波雷达的人体感知方法及系统 |
-
2021
- 2021-12-28 CN CN202111631545.9A patent/CN114469025B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106175723A (zh) * | 2016-06-27 | 2016-12-07 | 中国人民解放军第三军医大学第附属医院 | 一种基于fmcw宽带雷达的多生命监护系统 |
CN108919247A (zh) * | 2018-03-30 | 2018-11-30 | 中国科学院电子学研究所 | 一种基于恒虚警率检测的多目标体检测及定位方法 |
CN110584631A (zh) * | 2019-10-10 | 2019-12-20 | 重庆邮电大学 | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 |
CN113261925A (zh) * | 2021-03-11 | 2021-08-17 | 山东师范大学 | 一种基于毫米波雷达的人体感知方法及系统 |
Non-Patent Citations (3)
Title |
---|
YANYUN XU,ET AL: "A Novel Method for Automatic Detection of Trapped Victims by Ultrawideband Radar" * |
唐良勇;赵恒;张亚菊;: "基于WA-EMD算法的脉冲式超宽带雷达多目标生命体征检测" * |
荣霞等: "关于探地雷达对地下目标信号优化检测研究" * |
Also Published As
Publication number | Publication date |
---|---|
CN114469025B (zh) | 2023-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110584631B (zh) | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 | |
CN110531336A (zh) | 一种物体检测识别方法及系统 | |
CN106814358A (zh) | 一种用于超宽带穿墙雷达中多人体目标检测的系统及方法 | |
CN104076404B (zh) | 运用多通道相干抑制地磁背景噪声的磁异常探测方法 | |
CN115644840A (zh) | 基于毫米波雷达的生命体征检测方法 | |
CN116077044A (zh) | 基于毫米波雷达的生命体征检测方法 | |
CN112098970A (zh) | 一种用于交通微波检测的速度解模糊算法及相关设备 | |
CN114217301B (zh) | 一种高精度边坡监测雷达目标检测及测距方法 | |
CN105891795A (zh) | 一种双站rcs测量定标方法 | |
CN113812932A (zh) | 一种基于毫米波雷达的车内生命目标探测及定位技术 | |
CN114469025B (zh) | 一种用于fmcw雷达的多目标呼吸和心跳信号检测方法 | |
CN104758005A (zh) | 一种单通道乘法器解调的超声多普勒血流速度测量方法 | |
CN105824017B (zh) | 一种基于汽车补盲雷达的径向速度获取方法及装置 | |
CN112014833A (zh) | 一种高速目标时频域探测方法 | |
CN114983373B (zh) | 一种检测人体心率的方法 | |
CN114690143B (zh) | 一种抑制雷达杂波的方法、装置、雷达和介质 | |
CN114355329A (zh) | 一种调频连续波雷达生命体征距离检测的方法 | |
CN114859297A (zh) | 桥梁振动目标检测方法及装置 | |
CN110726988B (zh) | Pd雷达探测高超声速目标的距离和速度模糊互解方法 | |
CN114002670A (zh) | 一种微波微动传感器的信号处理方法、系统及相关组件 | |
CN109507662B (zh) | 一种低成本高精度的雷达方位测量方法 | |
Jagadesh | Development of Time Delay Estimation Algorithm Using Fuzzy Based Optimized Iterative Unscented Kalman Filter | |
Zhang et al. | An efficient track-before-detect algorithm based on complex likelihood ratio in radar systems | |
Wang et al. | An improved CFAR for life detection of UWB radar | |
CN116509365A (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 |