CN103728605A - 一种新型基于uwb雷达的非接触式生命特征信号提取方法 - Google Patents
一种新型基于uwb雷达的非接触式生命特征信号提取方法 Download PDFInfo
- Publication number
- CN103728605A CN103728605A CN201310749757.6A CN201310749757A CN103728605A CN 103728605 A CN103728605 A CN 103728605A CN 201310749757 A CN201310749757 A CN 201310749757A CN 103728605 A CN103728605 A CN 103728605A
- Authority
- CN
- China
- Prior art keywords
- record section
- vital sign
- data
- extracting method
- radar
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种新型基于UWB雷达的非接触式生命特征信号提取方法,包括以下步骤:(1)对原始雷达探测数据记录剖面进行道间以及道内数据压缩(2)对第二记录剖面在距离向上进行滑动加窗积分处理,(3)对第三记录剖面进行背景消除以及数据无量纲化处理,(4)对第四记录剖面在距离向上对时间向上的雷达数据进行频谱分析,(5)对第一时间记录序列采用自适应滤波以滤除环境噪声和随机噪声得到第二时间记录序列,即为最终所需要的人体生命特征信号。本发明所公开的新型基于UWB雷达的非接触式生命特征信号提取方法,实现简单,计算量小,便于DSP实时化片上处理,可实现人体生命特征信号的快速探测。
Description
技术领域
本发明涉及非接触式特征信号和微弱信号提取领域,尤其涉及基于UWB雷达的非接触式生命特征信号提取方法。
背景技术
利用UWB雷达的穿透能力强、抗干扰等特性进行人体生命体探测和搜救是当前的一个热点和新点,该新兴技术融合了雷达技术和生物医学工程技术,是国际科技界公认的一个非常重要的前沿技术领域。由于该技术可以非接触式远距离穿透非金属介质(墙体、废墟、砖块等)对人体目标进行有效的探测和识别,所以被广泛的应用于地震、泥石流、雪崩、火灾等灾害后被压埋幸存者或被困人员的搜救、反恐斗争中穿墙侦察等领域,除此之外还可用于战场伤员的搜救以及医疗非接触危重病人的远程监护等领域。
目前现行的基于UWB雷达的非接触式生命特征信号提取方法多采用构造高阶带通滤波器对经过预处理后的携带有人体生命特征信号的雷达回波进行滤波以实现噪声与生命特征信号分离,此方法一大缺点是生命特征信号为频率低,且频带窄,不便于通过构造带通滤波器来有效去除噪声,另外一点当预处理后的信号信噪比较低时无法有效提取人体生命特征信号。对基于UWB雷达的非接触式生命特征信号提取方法上还有一些研究者采用小波分析法、时频分析法、BP神经网络法以及相关法来进行分析处理,但这些方法计算复杂,不利于实时化处理。
发明内容
本发明所要解决的技术问题就是提供一种简单便于实时化的新型基于UWB雷达的非接触式生命特征信号提取方法。
本发明采用如下技术方案:
一种新型基于UWB雷达的非接触式生命特征信号提取方法,包括以下步骤:
(1)对原始雷达探测数据记录剖面也即第一记录剖面进行道间以及道内数据压缩,以减小雷达数据量形成第二记录剖面;
(2)对第二记录剖面在距离向上进行滑动加窗积分处理,用以增强生命特征信号形成第三记录剖面;
(3)对第三记录剖面进行背景消除以及数据无量纲化处理,用以消除强背景信号和距离因素的影响形成第四记录剖面;
(4)对第四记录剖面在距离向上对时间向上的雷达数据进行频谱分析,用以确定可能存在人体生命特征信号的距离位置,提取出第四记录剖面对应距离位置时间向上的雷达数据记为第一时间记录序列;
(5)对第一时间记录序列采用自适应滤波以滤除环境噪声和随机噪声得到第二时间记录序列,即为最终所需要的人体生命特征信号。
进一步的,步骤(1)的道内和道间数据压缩方法为:首先将第一记录剖面的相邻N道相加平均,其中N取值与UWB雷达的扫描速度相关,UWB雷达扫描速度越快则N取值越大,反之亦然;其次将第一记录剖面的每道内相邻M个数据点相加平均,其中M取值4或8。
进一步的,步骤(2)的距离向上进行滑动加窗处理为:依次对第二记录剖面中的每道数据中的每个数据点进行滑动加窗积分。
进一步的,步骤(3)的背景消除以及数据无量纲化处理为:将记录每道减去全部记录道的相加平均值,同时再除以相应时间序列的标准差。
进一步的,步骤(4)的频谱分析采用FFT快速完成,在32道/秒的扫描速度下,选择512点FFT为最佳。
进一步的,步骤(5)自适应滤波算法为LMS,自适应滤波器的长度选择为32,滤波器系数的初始值选择为滤波器长度的倒数即为1/32,收敛因子在10-8以下。
本发明的有益效果为:
本发明所公开的新型基于UWB雷达的非接触式生命特征信号提取方法,实现简单,计算量小,便于DSP实时化片上处理,可实现人体生命特征信号的快速探测。
附图说明
图1是本发明实施例1所公开的方法的雷达信号处理流程图;
图2是本发明实施例1中原始雷达探测数据记录剖面;
图3是本发明实施例1中预处理后雷达探测数据记录剖面;
图4是本发明实施例1中第一时间序列信号;
图5是本发明实施例1中自适应滤波算法原理框图;
图6是本发明实施例1中第二时间序列信号。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1,本实施例公开了一种新型基于UWB雷达的非接触式生命特征信号提取方法,包括以下步骤:
(1)对原始雷达探测数据记录剖面也即第一记录剖面进行道间以及道内数据压缩,以减小雷达数据量形成第二记录剖面;
(2)对第二记录剖面在距离向上进行滑动加窗积分处理,用以增强生命特征信号形成第三记录剖面;
(3)对第三记录剖面进行背景消除以及数据无量纲化处理,用以消除强背景信号和距离因素的影响形成第四记录剖面;
(4)对第四记录剖面在距离向上对时间向上的雷达数据进行频谱分析,用以确定可能存在人体生命特征信号的距离位置,提取出第四记录剖面对应距离位置时间向上的雷达数据记为第一时间记录序列;
(5)对第一时间记录序列采用自适应滤波以滤除环境噪声和随机噪声得到第二时间记录序列,即为最终所需要的人体生命特征信号。
进一步的,步骤(1)的道内和道间数据压缩方法为:首先将第一记录剖面的相邻N道相加平均,其中N取值与UWB雷达的扫描速度相关,UWB雷达扫描速度越快则N取值越大,反之亦然;其次将第一记录剖面的每道内相邻M个数据点相加平均,其中M取值4或8。
进一步的,步骤(2)的距离向上进行滑动加窗处理为:依次对第二记录剖面中的每道数据中的每个数据点进行滑动加窗积分。
进一步的,步骤(3)的背景消除以及数据无量纲化处理为:将记录每道减去全部记录道的相加平均值,同时再除以相应时间序列的标准差。
进一步的,步骤(4)的频谱分析采用FFT快速完成,在32道/秒的扫描速度下,选择512点FFT为最佳。
进一步的,步骤(5)自适应滤波算法为LMS,自适应滤波器的长度选择为32,滤波器系数的初始值选择为滤波器长度的倒数即为1/32,收敛因子在10-8以下。
具体的说,如图1至图6,图1为本发明整个雷达信号处理流程图,图2为原始雷达探测数据记录剖面,横坐标为时间序列,纵坐标为雷达时窗内采样点数,在实际UWB雷达的非接触式生命探测试验中获得的探测数据记录剖面,即为原始雷达探测数据记录剖面记为第一记录剖面。其中雷达的扫描速度为32道/秒,时窗内采样点数为4096点,第一记录剖面总共由1248个记录道组成。由雷达扫描速度和时窗内采样点数可以看出雷达数据量相对DSP处理器存储芯片而言相当大,要实现实时处理困难较大。
图3为由图2所示原始雷达探测数据记录首先经过道内雷达数据4点压缩,用以减小雷达探测记录数据量;再进行16点雷达距离向上滑动加窗积分处理,用以提高信噪比;最后进行背景消除以及数据无量纲化处理后所得记录剖面,用以消除雷达记录道数据的衰减及增益差异。其中背景消除和无量纲化处理具体公式为:
其中 X(i)为雷达记录剖面横向时间序列。其中S值的计算可以借助Matlab软件中的自带函数std来完成。
对图3中的所获得第四记录剖面进行频谱分析,可以采用FFT来快速完成,FFT点数的选择需要满足频域分辨率的要求,根据雷达的扫描速度而定。在32道/秒的扫描速度下,选择512点FFT为最佳。将FFT转换至频域信号在0.1Hz至0.8Hz范围内的最大数值作为比对对象,选出在距离向上的最大者并记下相应雷达距离向上的采样点值,提取出对应上述采样点值所对应的雷达横向时间序列记为第一时间记录序列。
图5为本发明所采用的自适应算法原理框图,采用的算法为LMS(最小均方误差),自适应滤波器的长度选择为32,滤波器系数的初始值选择为滤波器长度的倒数即为1/32,收敛因子的选择较为关键要足够小,一般在10-8以下较佳。在这里可以直接调用Matlab中的LMS函数来完成。
对图4采用图5中的自适应滤波方法可最终提取获得人体生命特征信号,如图6所示。
本实施例所公开的新型基于UWB雷达的非接触式生命特征信号提取方法,实现简单,计算量小,便于DSP实时化片上处理,可实现人体生命特征信号的快速探测。
Claims (6)
1.一种新型基于UWB雷达的非接触式生命特征信号提取方法,其特征在于:包括以下步骤:
(1)对原始雷达探测数据记录剖面也即第一记录剖面进行道间以及道内数据压缩,以减小雷达数据量形成第二记录剖面;
(2)对第二记录剖面在距离向上进行滑动加窗积分处理,用以增强生命特征信号形成第三记录剖面;
(3)对第三记录剖面进行背景消除以及数据无量纲化处理,用以消除强背景信号和距离因素的影响形成第四记录剖面;
(4)对第四记录剖面在距离向上对时间向上的雷达数据进行频谱分析,用以确定可能存在人体生命特征信号的距离位置,提取出第四记录剖面对应距离位置时间向上的雷达数据记为第一时间记录序列;
(5)对第一时间记录序列采用自适应滤波以滤除环境噪声和随机噪声得到第二时间记录序列,即为最终所需要的人体生命特征信号。
2.根据权利要求1所述的生命特征信号提取方法,其特征在于:步骤(1)的道内和道间数据压缩方法为:首先将第一记录剖面的相邻N道相加平均,其中N取值与UWB雷达的扫描速度相关,UWB雷达扫描速度越快则N取值越大,反之亦然;其次将第一记录剖面的每道内相邻M个数据点相加平均,其中M取值4或8。
3.根据权利要求1所述的生命特征信号提取方法,其特征在于:步骤(2)的距离向上进行滑动加窗处理为:依次对第二记录剖面中的每道数据中的每个数据点进行滑动加窗积分。
4.根据权利要求1所述的生命特征信号提取方法,其特征在于:步骤(3)的背景消除以及数据无量纲化处理为:将记录每道减去全部记录道的相加平均值,同时再除以相应时间序列的标准差。
5.根据权利要求1所述的生命特征信号提取方法,其特征在于:步骤(4)的频谱分析采用FFT快速完成,在32道/秒的扫描速度下,选择512点FFT为最佳。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310749757.6A CN103728605A (zh) | 2013-12-31 | 2013-12-31 | 一种新型基于uwb雷达的非接触式生命特征信号提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310749757.6A CN103728605A (zh) | 2013-12-31 | 2013-12-31 | 一种新型基于uwb雷达的非接触式生命特征信号提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103728605A true CN103728605A (zh) | 2014-04-16 |
Family
ID=50452760
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310749757.6A Pending CN103728605A (zh) | 2013-12-31 | 2013-12-31 | 一种新型基于uwb雷达的非接触式生命特征信号提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103728605A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110554381A (zh) * | 2019-08-30 | 2019-12-10 | 湖南正申科技有限公司 | 一种用于冲激脉冲式穿墙雷达的人体静止目标加速检测方法 |
CN113050060A (zh) * | 2021-03-25 | 2021-06-29 | 南京敏智达科技有限公司 | 基于稀疏重构与神经网络的生命探测方法及探测装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102058411A (zh) * | 2010-10-27 | 2011-05-18 | 中国人民解放军第四军医大学 | 一种多通道基于uwb雷达式生命探测仪 |
-
2013
- 2013-12-31 CN CN201310749757.6A patent/CN103728605A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102058411A (zh) * | 2010-10-27 | 2011-05-18 | 中国人民解放军第四军医大学 | 一种多通道基于uwb雷达式生命探测仪 |
Non-Patent Citations (3)
Title |
---|
何永波: "超宽带雷达信号微动特征识别研究", 《中国优秀硕士学位论文全文数据库基础科学辑》, no. 2, 15 February 2010 (2010-02-15), pages 1 - 30 * |
费翔宇、王君超、王春和: "穿墙定位雷达静止人体目标的信号处理", 《信号处理》, vol. 27, no. 5, 31 May 2011 (2011-05-31), pages 1 - 4 * |
费翔宇、王春和、黄素贞: "雷达式非接触生命参数信号的提取及判别", 《军事通信》, no. 9, 31 December 2009 (2009-12-31), pages 1 - 2 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110554381A (zh) * | 2019-08-30 | 2019-12-10 | 湖南正申科技有限公司 | 一种用于冲激脉冲式穿墙雷达的人体静止目标加速检测方法 |
CN113050060A (zh) * | 2021-03-25 | 2021-06-29 | 南京敏智达科技有限公司 | 基于稀疏重构与神经网络的生命探测方法及探测装置 |
CN113050060B (zh) * | 2021-03-25 | 2024-04-26 | 南京敏智达科技有限公司 | 基于稀疏重构与神经网络的生命探测方法及探测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105919584B (zh) | 用于可穿戴心率监测设备的心率估计方法及装置 | |
CN102058411B (zh) | 一种多通道基于uwb雷达式生命探测仪 | |
CN106127110B (zh) | 一种基于uwb雷达和最优svm的人体细粒度运动识别方法 | |
CN102008291A (zh) | 一种可用于多目标探测的单通道基于uwb的雷达式生命探测仪 | |
CN104266894B (zh) | 一种基于相关性分析的矿山微震信号初至波时刻提取方法 | |
CN106842164B (zh) | 基于小波去噪的非合作脉冲压缩雷达微弱目标检测方法 | |
CN103824302B (zh) | 基于方向波域图像融合的sar图像变化检测方法 | |
CN103245976B (zh) | 基于uwb生物雷达的人体目标与周围环境结构兼容探测方法 | |
CN107688553B (zh) | 基于小波变换和逻辑回归算法检测心电波形特征的方法 | |
CN106019254B (zh) | 一种uwb冲击生物雷达多人体目标距离向分离辨识方法 | |
CN101828918A (zh) | 基于波形特征匹配的心电信号r波峰检测方法 | |
CN107167802A (zh) | 一种基于超宽带雷达的呼吸信号检测算法 | |
CN104268883A (zh) | 一种基于边缘检测的时频谱曲线提取方法 | |
CN105223614B (zh) | 一种基于dwt_sta/lta的含噪信号p波初至峰度拾取方法 | |
CN101632587A (zh) | 一种基于小波变换的步态信号提取方法 | |
CN105426814A (zh) | 一种基于手机的老人跌倒检测方法 | |
CN101251445B (zh) | 旋转机械碰摩声发射信号的分形特征分析方法 | |
CN102636775B (zh) | 基于模糊逻辑识别的风廓线雷达回波谱重构方法 | |
CN109805936B (zh) | 基于地面振动信号的人体摔倒检测系统 | |
CN109633633A (zh) | 一种基于分段分类增强处理的生命信号增强方法 | |
CN104714215A (zh) | 一种基于杂波强度分区的恒虚警检测方法 | |
CN102018515B (zh) | 单通道基于uwb雷达式生命探测仪的多个人体目标识别和一维定位方法 | |
CN104036455A (zh) | 基于二代小波的红外图像细节增强方法 | |
CN102217932B (zh) | 一种abr信号波峰检测的算法 | |
CN102258368A (zh) | 胎心电检测的时域稀疏性线性混叠盲分离模型的判别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140416 |