胡 程 廖 鑫 向 寅 曾 涛
(北京理工大学信息与电子学院 北京 100081)
(北京市嵌入式实时信息处理技术重点实验室 北京 100081)
一种生命探测雷达微多普勒测量灵敏度分析新方法
胡 程 廖 鑫 向 寅*曾 涛
(北京理工大学信息与电子学院 北京 100081)
(北京市嵌入式实时信息处理技术重点实验室 北京 100081)
生命探测雷达技术是近些年发展起来的一种新型非接触式生命检测技术,其通过对生命体的振动信息进行测量实现生命体检测和定位。由于生命体的振动幅度一般较小,噪声会对生命探测雷达的测量性能造成严重的干扰,而噪声(热噪声和相位噪声)对此类雷达测量性能的影响目前尚无有效的评价指标。针对这一问题,该文通过建模分析定义了生命探测雷达的微多普勒测量灵敏度并提出了分析该指标的新方法,有效地解决了生命探测雷达测量性能无法定量分析的问题。最后通过仿真给出了给定系统在不同距离的测振性能曲线。
生命探测雷达;微多普勒测量;振动测量;等效平均振幅
生命探测雷达技术是一种综合运用生物医学、电子工程、雷达信号处理等技术的新型非接触式生命检测技术[1],该技术利用电磁波的多普勒效应实现生命体的振动信息探测和提取[2],从而对人体目标进行非接触式的探测、定位、成像等。由于该技术不受环境温度、热物体和声音干扰等优点[3],被广泛应用于地震、坍塌、建筑物倒塌下的废墟救援[4,5],反恐行动中对人员情况的掌握[6,7],生物医疗中对人体呼吸心跳信号的检测[8]等领域。
自从20世纪70年代 Caro C. G 和 Bloice J. A将雷达技术应用于生命体非接触式检测[9],生命探测雷达技术就引起了科学界的广泛关注,相关科研人员对振动测量反演方法[5]、穿墙探测杂波抑制方法[6]、运动模态识别技术[10]和近距离目标成像技术[11]等研究领域开展了研究,取得了大量的研究成果。
理论上取得突破的同时,为了满足实际应用的需求,美国、英国、加拿大、中国等国家的研究者研制了多个基于多普勒原理的生命探测雷达系统[3,12–14],但是相关的文献的内容主要集中在雷达信号处理和生命信号检测方面,并未给出具体测振性能的理论分析。然而,在实际雷达系统的研制和性能评估过程中,这一性能的分析是极其重要的。由于目前并无该性能有效的评价指标,使得生命探测雷达系统的参数设计和性能评估均很难精确地进行。
针对上述问题,本文仿照合成孔径雷达中等效噪声截面积的概念提出了微多普勒测量灵敏度的指标,用于衡量生命探测雷达系统振动测量性能。该指标将相位噪声和热噪声对测振精度的影响转化为等效平均振动幅度,然后定义目标平均振幅与相位噪声、热噪声等效平均振幅之比为相位测量信噪比,可以根据输出信噪比满足最小信噪比要求得到系统的微多普勒测量灵敏度。最后,利用本文提出的系统微多普勒测量灵敏度给出了示例系统并进行了分析,给出了该系统在不同距离处的测振性能。
2.1 线性近似测振模型
生命探测雷达一般可分为脉冲体制和连续波体制,两者均能很好地对微振动进行有效测量。本文提出的分析方法对两种体制均有效,但推导过程略有差异,不失一般性地,本部分先针对脉冲体制雷达进行建模分析,连续波体制雷达的分析将在第5部分给出,第5部分之前,所有的发射信号均认为是脉冲信号。
假设雷达向目标发射脉冲信号,发射信号可以写为:
其中,L(t)为发射信号的波形函数,A0为幅度,fc为载波频率,φ(t)为相位噪声。
记目标的微振动信号为M(t),目标的振动中心与雷达的距离为R0,则t时刻目标与雷达的距离为:
雷达接收到的回波信号可以写为:
其中,σ是目标后向散射系数,c为光速,N(t)是热噪声。
正交解调后的基带信号可以表示为复数形式:
因为雷达的脉冲持续时间很短,并且脉冲重复频率远远高于目标振动信号最高频率的两倍。因此,可以按照走-停-走模型建立微弱振动信号与高分辨1维像之间的调制关系。即可以近似认为在1帧时间内,目标静止不动;在相邻两帧之间目标表面位置发生微变,该微变可以由目标所在分辨单元相位变化同步反映。记tk表示快时间,tm表示慢时间,则式(4)可重新写为:
然后,对回波进行高分辨处理,得到含振动目标的高分辨1维像,如图1所示。高分辨合成以后的目标场景1维像时域表达为:
其中,AF(t)为L(t)的距离维模糊函数。根据匹配滤波理论,s(tk,tm)的峰值点应出现在处,从式(6)中不难看出,的相位中包含了振动的信息,通过对多个PRT的峰值位置相位分析可以提取出目标的振动。然而,式(6)的相位同时受到热噪声和相位噪声的干扰,为了验证上述振动测量方法的有效性和准确性,应对噪声对振动测量精度的影响程度进行建模分析。
图1 目标场景高分辨1维像Fig. 1 High resolution 1-D image of target scene
在式(6)基础上提取目标单元的慢时间信号。不失一般性,设目标在振动时不发生越距离单元走动,并且忽略目标相位中心与距离单元中心的距离差异,则目标单元处的高分辨像复数测量值可以表示为:
其中,S0表示目标1维像幅度,φ0是目标后向散射系数σ所携带的相位,Δφn(tm)=φn(tm)-φn(tm-2R0/c)是收发系统相位噪声在目标回波相位中的综合结果,是随时间变化的非白噪声。E, S, N分别表示目标所在单元总体回波以及其中的目标散射回波、接收机热噪声。目标单元的复矢量E是单元内信号复矢量S和噪声复矢量N的加和关系,如图2所示。
图2 目标距离单元复回波在IQ平面中的表示Fig. 2 Representation of the echo of the target range unit in the IQ plane
假设
即目标振动相位、目标上的相位噪声、热噪声与目标单元静止回波幅度的比值远远小1。针对微弱振动的测量,该假设通常是满足的。对式(7)做泰勒展开得到如下表达式:
将目标在静止状态下的相位作为固定相位从E的相位测量结果中减去,可以得到只与目标振动和噪声相关的相位项,即
其中,Re( ), Im( )分别表示复数的实部和虚部,o是各项高阶小量的统一表示。
在式(8)的小信号条件假设下,式(10)在1阶近似的条件下可以写成:
式(11)即为微振动高分辨相位导出测振模型,它揭示了目标振动、目标回波、相位噪声、热噪声对目标单元相位测量结果的贡献。
2.2 近似模型成立条件分析
相位导出测振的模型在建立过程使用了线性近似,仿真分析近似模型成立的条件。设实际测量的相位信号为通过模型计算的相位信号为则建模误差为
通过仿真分析不同振动幅度,建模误差RMSE和接收机信噪比SNR的关系。图3是在振动幅度分别为10–1λ, 10–2λ, 10–3λ, 10–4λ,建模误差RMSE随信噪比SNR变化的关系图。随着信噪比SNR的增大或者振动幅度的减小,建模误差将逐渐减小,但当振动幅度小于10–2λ时,信噪比将是影响建模误差的主要因素。从仿真结果可以看出当SNR>20 dB,振动幅度小于10–1λ时,都能满足RMSE< –20 dB。
图3 不同振幅建模误差RMSE随SNR变化Fig. 3 RMSE of different amplitude modeling varied with SNR
3.1 相位噪声等效平均振动幅度
式(11)中热噪声和系统相位噪声在测量相位上线性叠加,共同影响测振精度。式(11)中的Δφn(tm)是相位噪声噪声引起的,将相位噪声的平均功率记为PΦ,则:
其中,SΦ(fm)是相位噪声双边带功率谱密度,fm是频率,Bv为信号的带宽。相位噪声的双边带功率谱的形式[15]如图4所示。
图4 相位噪声双边带功率谱密度Fig. 4 Phase noise double sideband power spectral density
系统相位噪声主要取决于雷达电子系统器件水平。相位噪声对测振精度影响衡量指标可以仿造合成孔径雷达中的NEσ0的概念[16],将相位噪声的影响等效平均振动幅度,记为PEA0,则
3.2 热噪声等效平均振动幅度
其中,Pt为发射机发射功率,Gt为发射天线增益,Gr为接收天线增益,σ为目标后向散射系数,λ为波长,R为目标的距离,L为系统损耗,β为信号处理增益。
热噪声的平均功率与接收机的噪声系数,热噪声的带宽有关。将热噪声的平均功率记为PN,则:
其中,k为波尔兹曼常数,T0为室温,Fn为接收机噪声系数,Bn为热噪声的带宽。
同理,将热噪声对测振精度影响衡量指标转化为热噪声等效平均振动幅度,记为NEA0,则:
3.3 微多普勒测量灵敏度
定义式(11)中的目标平均振幅与相位噪声、热噪声等效平均振幅之比为相位测量信噪比,记为SNRM。则相位测量的信噪比SNRM为:
在输出信噪比满足最小信噪比要求SNRM≥SNRmin时,定义满足SNRM=SNRmin时的平均振动幅度为系统微多普勒测量灵敏度Mmin,则
下面利用上述推导结论进行系统测量性能分析的仿真。仿真所用的雷达系统参数如表1所示。
表1 雷达系统设计参数Tab. 1 Radar system design parameters
基于作者所在单位某在设计生命探测雷达在中心频率处的相噪仿真结果,可得到系统在X波段的相位噪声功率谱密度函数如图5所示。考虑到相位噪声功率谱密度函数具有轴对称性,图5只给出了单边的功率谱特性,另一侧的函数与图5的函数关于纵轴对称。
图5 X 波段雷达相位噪声指标Fig. 5 X-band radar phase noise index
将图5给出相位噪声指标代入到式(13),可以计算得到PEA0随距离变化的曲线如图6所示。从图6中不难看出,随着距离的增大,PEA0逐渐增大,且增长趋于平缓,与理论结果相符。
图6 PEA0随距离变化的曲线Fig. 6 PEA0curve varied with distance
基于上述参数设计,通过式(16)计算得到NEA0随距离变化的曲线如图7所示。可以看出,热噪声等效平均振幅的变化特性与相位噪声的变化趋势相似,但振幅要比相位噪声等效平均振幅大两个数量级以上,因此可以得出结论:热噪声对振动测量的影响远远大于相位噪声的影响,在系统分析与设计时,应将热噪声作为重点分析的对象。
图7 NEA0随距离变化的曲线Fig. 7 NEA0curve varied with distance
假设最终输出最小信噪比要求SNRmin=5 dB,根据式(18)计算系统的多普勒测量灵敏度随距离变化的曲线,如图8所示。在20 m的时候能够测量到0.5 mm左右的振动,50 m时能够测量到3 mm左右的振动,100 m时能够测量到12 mm左右的振动。
图8 微多普勒测量灵敏度随距离变化曲线Fig. 8 Micro-Doppler measurement sensitivity curve varied with distance
从该雷达系统的分析结果看出微多普勒测量灵敏度随着距离的增大而变低,其中影响微多普勒测量灵敏度的主要因素是热噪声。
一般而言,连续波体制的雷达往往采用两种波形––单频信号和调频信号。调频连续波信号一般采用去斜处理或匹配滤波处理,其处理结果的峰值相位与相应调频方式的脉冲信号相同,因此本文提出的方法对调频连续波体制的测振雷达同样有效。下面重点讨论单频连续波的情况。
单频连续波雷达发射的信号为单频连续波,其表达式为:
其中,A0为幅度,fc为载波频率,φ(t)为相位噪声。与第2部分的分析类似,可以得出正交解调后的基带信号为:
与脉冲体制的雷达不同,单频连续波体制的雷达回波具有恒定的幅度且无需进行匹配滤波处理,式(20)所示的回波信号的相位直接反映了振动的信息,与式(6)的相位相同。为了便于分析,不妨对其进行等间隔采样,记采样之后的时间轴为tm,则可重新写为:
此时得到的结论与(6)形式完全一致,不同的是由于单频连续波没有进行匹配滤波处理,没有距离分辨能力,在不考虑杂波的前提下其测振原理与调频连续波没有区别。
本文针对生命探测雷达系统设计和性能评估过程中缺少有效测振性能评价指标的问题,通过系统回波建模分析提出了微多普勒测量灵敏度的概念,并利用该概念给出了系统测振性能的具体分析方法,对生命探测雷达的系统参数设计提供了理论基础,并能很好地分析系统在不同场景、不同距离下的测振性能,对于生命探测雷达的后续研究具有重要的指导意义。
[1]路国华. 生物雷达目标信息识别技术的实验研究[D]. [博士论文], 第四军医大学, 2005. Lu Guohua. Experimental study on the technology of recognizing thetarget information based on bioradar[D]. [Ph.D. dissertation], Fourth Military Medical University, 2005.
[2]Skolnik M I. Introduction to Radar System-2nd edition[M]. New York: McGraw Hill, 1980.
[3]张杨, 王健琪, 荆西京, 等. 基于自适应滤波的生物雷达干扰抑制方法[J]. 仪器仪表学报, 2009, 30(4): 807–811. Zhang Yang, Wang Jianqi, Jing Xijing, et al.. Study on interference restraining of bioradar through adaptive filtering[J]. Chinese Journal of Scientific Instrument, 2009, 30(4): 807–811.
[4]Iii E F G. Radar flashlight for through-the-wall detection of humans[J]. Proceedings of SPIE - The International Society for Optical Engineering, 1998, 3375: 280–285.
[5]罗漫江. 生命探测仪的技术研究[D]. [硕士论文], 西安电子科技大学, 2004. Luo Manjiang. Research on technologies of radar bodydetect system[D]. [Master dissertation], Xidian University, 2004.
[6]Partida, Almorox, Mateo, et al.. Through-the-wall surveillance with millimeter-wave LFMCW radars[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(6): 1796–1805.
[7]He F, Zhu G F, Huang X T, et al.. Preliminary results of ultra-wideband through-the-wall life-detecting radar[C]. IEEE National Radar Conference, Washington, DC, 2010: 1327–1330.
[8]Kiriazi J E, Boric-Lubecke O, and Lubecke V M. Dualfrequency technique for assessment of cardiopulmonary effective RCS and displacement[J]. IEEE Sensors Journal, 2012, 12(3): 574–582.
[9]Bloice J A and Caro C G. Contactless apnoea detector using low energy radar[J]. Journal of Physiology, 1972, 223(1): 3–4.
[10]张翼, 邱兆坤, 朱玉鹏. 基于微多普勒信息的人体步态参数估计[J]. 信号处理, 2010, 26(6): 917–922. Zhang Yi, Qiu Zhaokun, and Zhu Yupeng. Human Gait parameter estimation based on micro Doppler feature[J]. Signal Processing, 2010, 26(6): 917–922.
[11]吴世有, 黄琼, 陈洁, 等. 基于超宽带穿墙雷达的目标定位识别算法[J]. 电子与信息学报, 2010, 32(11): 2624–2629. Wu Shiyou, Huang Qiong, Chen Jie, et al.. Target localization and identification algorithm for ultra wideband through-wall radar[J]. Journal of Electronics & Information Technology, 2010, 32(11): 2624–2629.
[12]Wang J Q, Zheng C X, Jin X J, et al.. Study on a noncontact life parameter detection system using millimeterwave[J]. Space Medicine & Medical Engineering, 2004, 17(3): 157–61.
[13]Franks C I, Watson J B, Brown B H, et al.. Respiratory patterns and risk of sudden unexpected death in infancy[J]. Archives of Disease in Childhood, 1980, 55(8): 595–599.
[14]Agrawal A, Fougere B, Oravakandy M, et al.. Design, construction and testing of a microwave radar system for through-wall surveillance[D]. University of Manitoba, Canada, 2004.
[15]Droitcour A D, Boric-Lubecke O, Lubecke V M, et al.. Range correlation and I/Q performance benefits in singlechip silicon Doppler radars for noncontact cardiopulmonary monitoring[J]. IEEE Transactions on Microwave Theory and Techniques, 2004, 52(3): 838–848.
[16]Skolnik M I. RADAR Handbook-2nd edition[M]. New York: McGraw Hill, 1990.
Novel Analytic Method for Determining Micro-Doppler Measurement Sensitivity in Life-detection Radar
Hu Cheng Liao Xin Xiang Yin Zeng Tao
(School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China)
(Beijing Key Laboratory of Embedded Real-Time Information Processing Technology, Beijing 100081, China)
In recent years, a new non-contact life detecting device has been developed, known as life-detection radar, which can measure bodily movement and locate human subjects. Typically, the amplitude of the vibration being measured is quite small, so the measurement is easily contaminated by noise in the radar system. To date, there is no effective index for judging the influence of noise on the vibration measurements in this radar system. To solve this problem, in this paper, we define the micro-Doppler measurement sensitivity to analyze the influence of noise on the measurement. We then perform a simulation to generate a performance curve for the radar system.
Life-detection radar; Micro-Doppler measurement; Vibration measurement; Equivalent average amplitude
TN959
A
2095-283X(2016)05-0455-07
10.12000/JR16090
胡程, 廖鑫, 向寅, 等. 一种生命探测雷达微多普勒测量灵敏度分析新方法[J]. 雷达学报, 2016, 5(5): 455–461.
10.12000/JR16090.
Reference format: Hu Cheng, Liao Xin, Xiang Yin, et al.. Novel analytic method for determining micro-Doppler measurement sensitivity in life-detection radar[J]. Journal of Radars, 2016, 5(5): 455–461. DOI: 10.12000/JR16090.
胡 程(1981–),湖南岳阳人,男,教授,博士生导师,北京理工大学雷达技术研究所副所长,IET Fellow、IEEE Senior Member。主要研究方向为地球同步轨道SAR、双基地/前向散射雷达信号处理及昆虫雷达等。
E-mail: hucheng.bit@gmail.com
廖 鑫(1991–),四川南充人,男,硕士研究生。主要研究方向为微多普勒雷达。
E-mail: liaoxinbit@126.com
向 寅(1981–),男,博士后,2010年于中国科技大学获得博士学位,后在中国科学院电子学研究所任博士后,研究方向为压缩感知、合成孔径雷达成像和信号处理。
E-mail: xy_overlimit@sina.cn
曾 涛(1971–),天津人,男,教授,博士生导师,北京理工大学雷达技术研究所实验室主任、航空学会电子专业委员会委员。主要研究方向为SAR成像技术和实时雷达信号处理。
E-mail: zengtao@bit.edu.cn
2016-08-17;改回日期:2016-10-18;
2016-10-28
向寅 xy_overlimit@sina.cn
国家自然科学基金(61120106004, 61225005),高等学校学科创新引智计划(B14010)
Foundation Items: The National Natural Science Foundation of China (61120106004, 61225005), The 111 Project of China (B14010)