CN106291524A - 基于人体模型的lfmcw雷达探测运动人体目标的方法 - Google Patents
基于人体模型的lfmcw雷达探测运动人体目标的方法 Download PDFInfo
- Publication number
- CN106291524A CN106291524A CN201610642463.7A CN201610642463A CN106291524A CN 106291524 A CN106291524 A CN 106291524A CN 201610642463 A CN201610642463 A CN 201610642463A CN 106291524 A CN106291524 A CN 106291524A
- Authority
- CN
- China
- Prior art keywords
- gait
- frequency
- data
- lfmcw
- lfmcw 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.)
- Granted
Links
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于人体模型的LFMCW雷达探测运动人体目标的方法,步骤如下:将LFMCW雷达回波差拍信号进行模数转换,并对每个PRI回波差拍数据进行FFT处理,获得人体慢速运动时LFMCW雷达差拍信号的频域数据并在慢时间‑距离域二维平面重新排列;构建Boulic人体步行模型,并根据人体运动步态频率和初始相位范围,在步态频率‑初始相位二维平面上构建多个不同参数的人体躯干运动模型,构造出LFMCW雷达回波非线性相位补偿信号并对重排数据进行相位补偿,然后进行FFT处理得到四维空间的处理结果;通过恒虚警检测处理判决出某距离单元是否存在目标,并在存在目标的情况下得到目标的距离、速度与步态频率信息。本发明提升了LFMCW雷达检测的信噪比,从而改善了其在强地杂波中对人体目标的探测性能。
Description
技术领域
本发明属于雷达信息探测技术领域,特别是一种基于人体模型的LFMCW雷达探测运动人体目标的方法。
技术背景
在雷达系统中,目标探测是最重要的任务之一,线性调频连续波(LinearFrequency Modulation Continuous Wave,LFMCW)雷达具有体积小、重量轻、距离分辨力高、无距离盲区和抗截获能力强等独特优势,在地面侦察场景中可以用于探测敌方的人员、车和坦克等运动目标,并获得目标的距离、速度和方位等信息。
在实际环境中,由于车辆、坦克等目标的雷达散射截面积(Radar Cross Section,RCS)较大并且运动速度快,LFMCW雷达能够较为容易地实现探测。而对于人这种地面慢速弱目标的探测,技术难度却很大,其瓶颈主要在于:1)这类目标多处于复杂的强地杂波中,具有反射系数小、RCS起伏大的特点,导致雷达对此类目标探测能力较弱。2)人体慢速目标在多普勒域与地物杂波谱(尤其是运动杂波)存在严重的交叠,当目标位于风吹动的草地或树林时,草地和树林杂波能量扩展到慢速目标的多普勒通道,使得目标回波与地物杂波混叠,难以从多普勒维上分辨出目标和杂波。3)人体行走属于典型的非刚体运动。各个关节同时协调摆动导致人体行走时的雷达回波多普勒成分复杂,使得雷达回波在多普勒域上产生较大的多普勒扩展,导致在积累时目标的能量分散,加大目标探测难度。
LFMCW雷达对地面慢速人体目标的探测主要依赖于接收到的回波信噪比。根据雷达方程,我们可以通过提高雷达发射功率、增大天线孔径、降低接收机噪声系数、优选发射频率等措施来提高雷达接收机的灵敏度,但是以上方法都需要增加硬件成本。相对而言,利用信号处理的方法来提升LFMCW雷达系统的信噪比从而实现对人体目标的探测,成本较低,手段灵活,能够尽最大可能的发挥雷达的潜力。
信号处理中提高检测信噪比的方法有增加相参积累时间、多帧数据联合处理等方法。然而人体步行并非匀速运动,而是一个类似于钟摆的运动。在前进过程中,人体四肢周期性的旋转带动躯干重复摆动,从而产生微多普勒回波。这些微多普勒分布在人体平均径向速度两侧,导致长时间相参积累后多普勒的扩散,使得目标信号的能量不能得到充分积累。
发明内容
本发明的目的在于提供一种LFMCW雷达探测运动人体目标的方法,将Boulic人体运动模型与LFMCW雷达信号处理相结合,通过人体运动非线性相位补偿将多普勒扩展能量聚焦,提升LFMCW雷达检测信噪比,从而改善其在强地杂波中对人体目标的探测性能。
实现本发明目的技术解决方案为:一种基于人体模型的LFMCW雷达探测运动人体目标的方法,包括以下步骤:
步骤1,将LFMCW雷达回波差拍信号进行模数转换,然后依次对每个脉冲重复周期回波差拍数据进行快速傅里叶变换处理,获得人体慢速运动时LFMCW雷达差拍信号的频域数据,并将该数据在慢时间-距离域二维平面重新排列;
步骤2,根据LFMCW雷达接收回波时长构建Boulic人体步行模型,并根据Boulic模型中的人体运动步态频率和初始相位范围,在步态频率-初始相位二维平面上构建多个不同参数的人体躯干运动模型,将躯干模型与发射波形参数结合,构造出LFMCW雷达回波非线性相位补偿信号;
步骤3,使用构建的相位补偿信号对步骤1中重排数据进行四维空间相位补偿,然后对补偿后的结果进行FFT处理得到多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果;
步骤4,对每个距离单元下不同初始相位、步态频率对应的多普勒维数据进行恒虚警检测处理判决出该距离单元是否存在目标,并在存在目标的情况下得到目标的距离、速度与步态频率信息。
进一步地,步骤1所述将LFMCW雷达回波差拍信号进行模数转换,然后依次对每个脉冲重复周期回波差拍数据进行快速傅里叶变换处理,获得人体慢速运动时LFMCW雷达差拍信号的频域数据,并将该数据在慢时间-距离域二维平面重新排列,具体如下:
(1.1)对每个脉冲重复周期内雷达回波差拍信号Δf(t)进行A/D采样处理,将第m个PRI的采样结果记为xm:
xm=[x(m,1),x(m,2),...,x(m,n),...,x(m,N)]T (1)
其中,x(m,n)表示在第m个PRI内模数转换后的第n点回波时域数据;
(1.2)对xm进行N点快速傅里叶变换处理,获得第m个PRI内N点频域结果Xm:
Xm=[X(m,1),X(m,2),...,X(m,n),...,X(m,N)]T (2)
其中,X(m,n)表示在第m个PRI内FFT后的第n点频域数据;
(1.3)重复步骤(1.1)~(1.2),采集M个PRI内的差拍频域数据作为一个相参处理间隔,记为X:
X=[X1,X2,...,Xm,...,XM]T (3)
将X按照慢时间-距离域二维平面排列,完成频域数据重排。
进一步地,步骤2所述根据LFMCW雷达接收回波时长构建Boulic人体步行模型,并根据Boulic模型中的人体运动步态频率和初始相位范围,在步态频率-初始相位二维平面上构建多个不同参数的人体躯干运动模型,将躯干模型与发射波形参数结合,构造出LFMCW雷达回波非线性相位补偿信号,具体如下:
(2.1)根据Boulic人体运动模型,进行步态频率-初始相位平面人体躯干运动模型的构造,在构造躯干位移模型时,躯干的运动时长与X时长对应,数据X的总时长Tm为:
Tm=M·Tr (4)
其中Tr是LFMCW雷达发射信号重复周期;
(2.2)在Boulic模型中,人体相对步行速度vr范围是0m/s到3m/s,根据相对步行速度vr和步态频率fgait间的关系,得fgait范围是0Hz到1.3Hz,关系式如下:
vr=(fgait·1.346)2 (5)
将fgait离散为P点,则其中第p个步态频率fgait(p)为
同理初始相位范围是0到2π,将离散成Q点,则第q个初始相位为:
(2.3)在fgait(p)和条件下第m个脉冲重复周期对应的人体躯干运动位移sp,q(m)如下:
依次改变m取值获得单个相参处理间隔下人体躯干运动模型sp,q:
sp,q=[sp,q(1),sp,q(2),...,sp,q(m),...,sp,q(M)]T (9)
(2.4)将sp,q与LFMCW雷达发射波形参数相结合,构建fgait(p)和条件下单个相参处理间隔对应的LFMCW雷达回波非线性相位补偿信号Hp,q:
其中,fc是LFMCW雷达发射信号载频,c是光速;
(2.5)分别遍历所有p、q的值,完成步态频率-初始相位平面LFMCW回波相位补偿信号的构造。
进一步地,步骤3所述使用构建的相位补偿信号对步骤1中重排数据进行四维空间相位补偿,然后对补偿后的结果进行FFT处理得到多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果,具体如下:
(3.1)对于步骤1中重排数据X,用表示X中第n个距离单元下M点慢时间维数据,如下:
(3.2)用Hp,q对进行相位补偿,即
其中表示克罗内克积,所得结果Yn,p,q表示第n个距离单元中fgait(p)和条件下相位补偿后的M点慢时间维数据;
(3.3)针对依次改变n的取值,同时遍历步态频率-初始相位二维平面内所有p、q对应的Hp,q,重复式(12)的相位补偿处理,完成慢时间-距离单元-步态频率-初始相位四维空间相位补偿;
(3.4)对四维空间补偿结果中慢时间维M点数据Yn,p,q进行快速傅里叶变换处理,得到第n个距离单元下fgait(p)和对应的多普勒维数据Fn,p,q,即
Fn,p,q=FFT(Yn,p,q) (13)
依次对不同的Yn,p,q进行式(13)操作,完成四维空间慢时间维M点快速傅里叶变换处理,从而获得多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果。
进一步地,步骤4所述对每个距离单元下不同初始相位、步态频率对应的多普勒维数据进行恒虚警检测处理判决出该距离单元是否存在目标,并在存在目标的情况下得到目标的距离、速度与步态频率信息,具体如下:
依次对四维空间中多普勒维数据Fn,p,q进行恒虚警检测,判断出Fn,p,q对应的距离单元有无目标的存在,若在该距离单元发现目标,则根据相应距离单元下幅值最大值对应的多普勒单元、fgait(p)和计算出目标的速度、最佳运动参数和从而得到目标运动信息。
本发明与现有技术相比,其显著优点为:(1)将人体运动的特性和LFMCW雷达信号处理相结合,在步态频率-初始相位二维平面上构建的多普勒相位补偿信号和雷达人体回波中的非线性相位相互抵消,可以使人体回波信号同相相加,改善了人体目标回波能量积累的效果;(2)简化了Boulic人体运动模型,只考虑人体躯干运动引起的多普勒扩展,原理简单,实现容易;(3)针对慢速运动的人体目标,利用Boulic模型中运动速度和步态频率间的关系来限定步态频率变化范围,减小了在步态频率-初始相位二维平面补偿处理的计算量;(4)抑制了动目标检测(Moving Target Detection,MTD)处理后多普勒扩散现象,有效地提高了LFMCW雷达检测信噪比(Signal-to-Noise Ratio,SNR),增强了其对慢速运动人体的探测能力。
附图说明
图1为相位补偿法LFMCW雷达信号处理流程图。
图2为LFMCW雷达收发信号、回波差拍信号与采样时刻示意图。
图3为差拍频域数据重排示意图。
图4为四维空间中单个距离单元内相位补偿示意图。
图5为四维空间中单个距离单元内FFT处理示意图。
具体实施方式
本发明为一种基于人体模型的LFMCW雷达探测运动人体目标的方法,思路如下:人体步行时躯干和四肢摆动产生的微多普勒随时间是有规律变化的,根据Boulic模型可知,人体步行是一种关节连接的运动,人体的步态是大脑、肌肉、神经、关节和骨骼高度协调的周期性运动,而LFMCW回波相位中的多普勒调制是由人体各个部位在雷达径向上运动位移引起的。由此可得,人体步行对LFMCW雷达回波相位的调制也具有周期性的变化。利用人体运动的模型可以推导出LFMCW雷达回波相位变化的规律,根据回波相位变化规律可以在长时间相参积累前对雷达回波进行相位补偿,这样在雷达回波进行相参积累时可以使目标信号同相相加,从而提升LFMCW雷达的信噪比,达到增强LFMCW雷达在强地杂波下探测运动人体目标能力的目的。
结合图1,本发明LFMCW雷达探测运动人体目标的方法,具体步骤如下:
步骤1,将LFMCW雷达回波差拍信号进行模数转换,然后依次对每个脉冲重复周期(Pulse Repetition Interval,PRI)回波差拍数据进行快速傅里叶变换(Fast FourierTransform,FFT)处理,获得人体慢速运动时LFMCW雷达差拍信号的频域数据,并将该数据在慢时间-距离域二维平面重新排列,具体如下:
(1.1)对每个PRI内雷达回波差拍信号Δf(t)进行图1中“A/D采样”处理。其中LFMCW雷达收发信号、回波差拍信号Δf(t)及采样时刻如图2所示。将第m个PRI的采样结果记为xm:
xm=[x(m,1),x(m,2),...,x(m,n),...,x(m,N)]T (14)
x(m,n)表示在第m个PRI内模数转换后的第n点回波时域数据。在A/D采样时,需要将最大作用距离对应回波时延之前的数据舍弃。
(1.2)接着对xm进行图1中“N点FFT处理”,获得第m个PRI内N点频域结果Xm:
Xm=[X(m,1),X(m,2),...,X(m,n),...,X(m,N)]T (15)
其中X(m,n)表示在第m个PRI内FFT后的第n点频域数据。在Xm中,点与点之间的间隔对应LFMCW雷达的距离分辨率,这N点频域信息代表了N个距离单元。
(1.3)重复步骤(1.1)~(1.2),采集M个PRI内的差拍频域数据作为一个相参处理间隔(Coherent Processing Interval,CPI),记为X:
X=[X1,X2,...,Xm,...,XM]T (16)
将X按照图3所示的慢时间-距离域二维平面排列,完成图1中“频域数据重排”。
步骤2,根据LFMCW雷达接收回波时长构建Boulic人体步行模型。在Boulic模型中,人体步行是由头部、躯干、四肢等12个关键部位相互协调运动组成。但与人体其它部位相比,躯干的RCS最大,回波中躯干能量最强,因此只考虑人体躯干运动引起的多普勒扩展。在构建躯干运动模型时,根据Boulic模型中给出的人体运动步态频率和初始相位范围,在步态频率-初始相位二维平面上构建多个不同参数的人体躯干运动模型。然后将躯干模型与发射波形参数结合,构造出LFMCW雷达回波非线性相位补偿信号,具体如下:
(2.1)首先根据Boulic人体运动模型,进行图1的“步态频率-初始相位平面人体躯干运动模型的构造”。在构造躯干位移模型时,躯干的运动时长需要与X时长对应。数据X的总时长Tm为
Tm=M·Tr (17)
其中Tr是LFMCW雷达发射信号重复周期。由于人体运动缓慢,因此躯干运动位移在单个PRI内近似不变,只在PRI之间变化。
(2.2)在Boulic模型中,人体相对步行速度vr范围是0m/s到3m/s。根据相对步行速度vr和步态频率fgait间的关系
vr=(fgait·1.346)2 (18)
可得fgait范围是0Hz到1.3Hz。将fgait离散为P点,则其中第p个步态频率fgait(p)为
同理初始相位范围是0到2π,将离散成Q点,则第q个初始相位为
(2.3)这样,在fgait(p)和条件下第m个PRI对应的人体躯干运动位移:
依次改变m取值获得单个CPI下人体躯干运动模型:
sp,q=[sp,q(1),sp,q(2),...,sp,q(m),...,sp,q(M)]T (22)
(2.4)然后将sp,q与LFMCW雷达发射波形参数相结合,忽略次要因素的影响,构建fgait(p)和条件下单个CPI对应的LFMCW雷达回波非线性相位补偿信号Hp,q
其中fc是LFMCW雷达发射信号载频,c是光速。
(2.5)分别遍历所有p、q的值,重复上述处理,完成图1中“步态频率-初始相位平面LFMCW回波相位补偿信号的构造”。
步骤3,使用构建的相位补偿信号对(1)中重排数据进行四维空间相位补偿,然后对补偿后的结果进行FFT处理得到多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果,具体如下:
(3.1)对于步骤1中重排数据X,用表示X中第n个距离单元下M点慢时间维数据,如下:
(3.2)首先用Hp,q按照图4所示对进行相位补偿,即
其中表示克罗内克积,所得结果Yn,p,q表示第n个距离单元中fgait(p)和条件下相位补偿后的M点慢时间维数据。
(3.3)针对依次改变n的取值,同时遍历步态频率-初始相位二维平面内所有p、q对应的Hp,q,重复式(12)的相位补偿处理,完成图1中“慢时间-距离单元-步态频率-初始相位四维空间相位补偿”。
(3.4)然后对四维空间补偿结果中慢时间维M点数据Yn,p,q进行如图5所示FFT处理,得到第n个距离单元下fgait(p)和对应的多普勒维数据Fn,p,q,即
Fn,p,q=FFT(Yn,p,q) (26)
依次对不同的Yn,p,q进行式(13)操作,完成图1中“四维空间慢时间维M点FFT处理”,从而获得多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果。
步骤4,对每个距离单元下不同初始相位、步态频率对应的多普勒维数据进行恒虚警(Constant False Alarm Rate,CFAR)检测处理判决出该距离单元是否存在目标,并在存在目标的情况下得到目标的距离、速度与步态频率等信息。
依次对四维空间中多普勒维数据Fn,p,q进行图1中“CFAR检测”,判断出Fn,p,q对应的距离单元有无目标的存在。若在该距离单元发现目标,则根据相应距离单元下幅值最大值对应的多普勒单元、fgait(p)和计算出目标的速度、最佳运动参数和从而得到目标运动信息。
实施例1
本专利测试采用的LFMCW雷达系统由收发天线、发射机、接收机、频率综合器、信号处理器和目标显示器等组成。该雷达整机性能参数如表1:
表1雷达整机性能参数
用上述雷达采集人在不同距离情况下缓慢前进的频域数据。其中人面向雷达径向运动,速度范围为0.8m/s到1.8m/s。通过对采集的数据进行直接MTD和相位补偿处理,对比得出两种方法下目标峰值能量和噪声平均功率的变化以及SNR的提升情况,其中在相位补偿处理时将步态频率-初始相位二维平面分为P×Q点(P为76,Q为72)。
表2、表3是人距离雷达27m、45m、63m和81m时,采集实测数据并分别用MTD和相位补偿两种方法处理的结果。从表2、表3可以看出,雷达发射相同功率信号时,目标距离雷达越远,采用MTD处理和相位补偿处理后的回波中目标峰值能量都变小,而系统的噪声平均功率几乎不变。比较表2、表3中目标位于相同距离单元下两种方法处理的结果可以发现:相位补偿处理后的目标峰值能量比MTD后的能量有所增加,而两者的噪声基底平均功率变化不大。
表2不同距离时MTD处理结果
表3不同距离时相位补偿处理结果
表4不同距离时两者处理SNR变化
表4是不同距离时两者处理后SNR变化情况。从表4能够得出,相位补偿处理后的SNR和MTD处理相比有较好的改善,提升大约4dB到5dB。
相位补偿之所以可以增加SNR是因为该方法通过构建多个不同运动参数下的多普勒相位补偿信号,采用遍历的方式与LFMCW雷达回波信号相位匹配,估计出目标的最佳运动参数(和),使得回波中的非线性相位相互抵消,让目标信号可以同相累加。这样MTD处理后的多普勒扩散能够得到充分利用,可以使目标峰值能量增加,而噪声平均功率变化很小,从而达到提升LFMCW雷达检测SNR的目的。由此可得,与MTD处理相比,对LFMCW雷达强地杂波下的慢速人体回波信号进行相位补偿处理,目标峰值能量增加,噪声平均功率不变,检测SNR有较好的提升,LFMCW雷达的人体慢速运动探测性能增强。
Claims (5)
1.一种基于人体模型的LFMCW雷达探测运动人体目标的方法,其特征在于,包括以下步骤:
步骤1,将LFMCW雷达回波差拍信号进行模数转换,然后依次对每个脉冲重复周期回波差拍数据进行快速傅里叶变换处理,获得人体慢速运动时LFMCW雷达差拍信号的频域数据,并将该数据在慢时间-距离域二维平面重新排列;
步骤2,根据LFMCW雷达接收回波时长构建Boulic人体步行模型,并根据Boulic模型中的人体运动步态频率和初始相位范围,在步态频率-初始相位二维平面上构建多个不同参数的人体躯干运动模型,将躯干模型与发射波形参数结合,构造出LFMCW雷达回波非线性相位补偿信号;
步骤3,使用构建的相位补偿信号对步骤1中重排数据进行四维空间相位补偿,然后对补偿后的结果进行FFT处理得到多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果;
步骤4,对每个距离单元下不同初始相位、步态频率对应的多普勒维数据进行恒虚警检测处理判决出该距离单元是否存在目标,并在存在目标的情况下得到目标的距离、速度与步态频率信息。
2.根据权利要求1所述的基于人体模型的LFMCW雷达探测运动人体目标的方法,其特征在于,步骤1所述将LFMCW雷达回波差拍信号进行模数转换,然后依次对每个脉冲重复周期回波差拍数据进行快速傅里叶变换处理,获得人体慢速运动时LFMCW雷达差拍信号的频域数据,并将该数据在慢时间-距离域二维平面重新排列,具体如下:
(1.1)对每个脉冲重复周期内雷达回波差拍信号Δf(t)进行A/D采样处理,将第m个PRI的采样结果记为xm:
xm=[x(m,1),x(m,2),...,x(m,n),...,x(m,N)]T (1)
其中,x(m,n)表示在第m个PRI内模数转换后的第n点回波时域数据;
(1.2)对xm进行N点快速傅里叶变换处理,获得第m个PRI内N点频域结果Xm:
Xm=[X(m,1),X(m,2),...,X(m,n),...,X(m,N)]T (2)
其中,X(m,n)表示在第m个PRI内FFT后的第n点频域数据;
(1.3)重复步骤(1.1)~(1.2),采集M个PRI内的差拍频域数据作为一个相参处理间隔,记为X:
X=[X1,X2,...,Xm,...,XM]T (3)
将X按照慢时间-距离域二维平面排列,完成频域数据重排。
3.根据权利要求1所述的基于人体模型的LFMCW雷达探测运动人体目标的方法,其特征在于,步骤2所述根据LFMCW雷达接收回波时长构建Boulic人体步行模型,并根据Boulic模型中的人体运动步态频率和初始相位范围,在步态频率-初始相位二维平面上构建多个不同参数的人体躯干运动模型,将躯干模型与发射波形参数结合,构造出LFMCW雷达回波非线性相位补偿信号,具体如下:
(2.1)根据Boulic人体运动模型,进行步态频率-初始相位平面人体躯干运动模型的构造,在构造躯干位移模型时,躯干的运动时长与X时长对应,数据X的总时长Tm为:
Tm=M·Tr (4)
其中Tr是LFMCW雷达发射信号重复周期;
(2.2)在Boulic模型中,人体相对步行速度vr范围是0m/s到3m/s,根据相对步行速度vr和步态频率fgait间的关系,得fgait范围是0Hz到1.3Hz,关系式如下:
vr=(fgait·1.346)2 (5)
将fgait离散为P点,则其中第p个步态频率fgait(p)为
同理初始相位范围是0到2π,将离散成Q点,则第q个初始相位为:
(2.3)在fgait(p)和条件下第m个脉冲重复周期对应的人体躯干运动位移sp,q(m)如下:
依次改变m取值获得单个相参处理间隔下人体躯干运动模型sp,q:
sp,q=[sp,q(1),sp,q(2),...,sp,q(m),...,sp,q(M)]T (9)
(2.4)将sp,q与LFMCW雷达发射波形参数相结合,构建fgait(p)和条件下单个相参处理间隔对应的LFMCW雷达回波非线性相位补偿信号Hp,q:
其中,fc是LFMCW雷达发射信号载频,c是光速;
(2.5)分别遍历所有p、q的值,完成步态频率-初始相位平面LFMCW回波相位补偿信号的构造。
4.根据权利要求1所述的基于人体模型的LFMCW雷达探测运动人体目标的方法,其特征在于,步骤3所述使用构建的相位补偿信号对步骤1中重排数据进行四维空间相位补偿,然后对补偿后的结果进行FFT处理得到多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果,具体如下:
(3.1)对于步骤1中重排数据X,用表示X中第n个距离单元下M点慢时间维数据,如下:
(3.2)用Hp,q对进行相位补偿,即
其中表示克罗内克积,所得结果Yn,p,q表示第n个距离单元中fgait(p)和条件下相位补偿后的M点慢时间维数据;
(3.3)针对依次改变n的取值,同时遍历步态频率-初始相位二维平面内所有p、q对应的Hp,q,重复式(12)的相位补偿处理,完成慢时间-距离单元-步态频率-初始相位四维空间相位补偿;
(3.4)对四维空间补偿结果中慢时间维M点数据Yn,p,q进行快速傅里叶变换处理,得到第n个距离单元下fgait(p)和对应的多普勒维数据Fn,p,q,即
Fn,p,q=FFT(Yn,p,q) (13)
依次对不同的Yn,p,q进行式(13)操作,完成四维空间慢时间维M点快速傅里叶变换处理,从而获得多普勒频率-距离单元-初始相位-步态频率四维空间的处理结果。
5.根据权利要求1所述的基于人体模型的LFMCW雷达探测运动人体目标的方法,其特征在于,步骤4所述对每个距离单元下不同初始相位、步态频率对应的多普勒维数据进行恒虚警检测处理判决出该距离单元是否存在目标,并在存在目标的情况下得到目标的距离、速度与步态频率信息,具体如下:
依次对四维空间中多普勒维数据Fn,p,q进行恒虚警检测,判断出Fn,p,q对应的距离单元有无目标的存在,若在该距离单元发现目标,则根据相应距离单元下幅值最大值对应的多普勒单元、fgait(p)和计算出目标的速度、最佳运动参数和从而得到目标运动信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610642463.7A CN106291524B (zh) | 2016-08-05 | 2016-08-05 | 基于人体模型的lfmcw雷达探测运动人体目标的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610642463.7A CN106291524B (zh) | 2016-08-05 | 2016-08-05 | 基于人体模型的lfmcw雷达探测运动人体目标的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291524A true CN106291524A (zh) | 2017-01-04 |
CN106291524B CN106291524B (zh) | 2019-01-22 |
Family
ID=57666461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610642463.7A Active CN106291524B (zh) | 2016-08-05 | 2016-08-05 | 基于人体模型的lfmcw雷达探测运动人体目标的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291524B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107064927A (zh) * | 2017-03-21 | 2017-08-18 | 西北大学 | 一种高速弹载雷达目标相参积累方法 |
CN107907878A (zh) * | 2017-11-08 | 2018-04-13 | 零八电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN108537144A (zh) * | 2018-03-21 | 2018-09-14 | 特斯联(北京)科技有限公司 | 一种多维人体步态识别方法与设备 |
CN109917347A (zh) * | 2019-04-10 | 2019-06-21 | 电子科技大学 | 一种基于时频域稀疏重构的雷达行人检测方法 |
CN110488279A (zh) * | 2019-08-30 | 2019-11-22 | 成都纳雷科技有限公司 | 一种基于雷达的行人与树木区分方法与装置 |
CN110501685A (zh) * | 2019-08-23 | 2019-11-26 | 北京电子工程总体研究所 | 一种基于雷达信号相位补偿的多帧相参积累方法 |
CN110584631A (zh) * | 2019-10-10 | 2019-12-20 | 重庆邮电大学 | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 |
CN110673135A (zh) * | 2018-07-03 | 2020-01-10 | 松下知识产权经营株式会社 | 传感器、推定装置、推定方法和程序记录介质 |
CN111160376A (zh) * | 2019-12-31 | 2020-05-15 | 联想(北京)有限公司 | 一种信息处理方法、装置、电子设备和存储介质 |
CN111337897A (zh) * | 2020-04-21 | 2020-06-26 | 湖南红船科技有限公司 | 一种lfmcw雷达快速目标识别方法 |
CN111693981A (zh) * | 2020-05-27 | 2020-09-22 | 清华大学 | 变重复周期的二次相关处理方法 |
CN112327286A (zh) * | 2020-10-13 | 2021-02-05 | 深圳大学 | 低复杂度下日常活动分类方法、装置、设备及存储介质 |
CN112639523A (zh) * | 2020-06-30 | 2021-04-09 | 华为技术有限公司 | 一种雷达探测方法及相关装置 |
CN113017590A (zh) * | 2021-02-26 | 2021-06-25 | 清华大学 | 生理数据监测方法、装置、计算机设备和存储介质 |
CN113064131A (zh) * | 2021-03-18 | 2021-07-02 | 江苏蛮酷科技有限公司 | 雷达波处理方法、装置和计算机存储介质 |
CN113096791A (zh) * | 2020-01-08 | 2021-07-09 | 富士通株式会社 | 基于无线信号的人员识别装置、方法及系统 |
CN116453227A (zh) * | 2023-06-19 | 2023-07-18 | 武汉理工大学 | 一种基于双毫米波雷达的船舶环境下步态识别方法 |
CN117930213A (zh) * | 2024-03-25 | 2024-04-26 | 湘江实验室 | 一种基于毫米波雷达探测人体存在的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009062773A1 (de) * | 2007-11-15 | 2009-05-22 | Robert Bosch Gmbh | Verfahren zum erkennen von lebewesen |
CN102928835A (zh) * | 2012-10-09 | 2013-02-13 | 北京航空航天大学 | 一种基于改进广义s变换的人体目标运动状态识别方法 |
WO2014120289A1 (en) * | 2012-10-10 | 2014-08-07 | Raytheon Company | Detection of concealed object on a body using radio frequency signatures on frequencies and polarizations |
CN104237865A (zh) * | 2014-10-07 | 2014-12-24 | 电子科技大学 | 一种基于太赫兹雷达回波的人体目标微动信号时频分析方法 |
CN104267394A (zh) * | 2014-10-07 | 2015-01-07 | 电子科技大学 | 高分辨率人体目标运动特征检测方法 |
CN105137423A (zh) * | 2015-09-30 | 2015-12-09 | 武汉大学 | 一种穿墙雷达对多运动目标实时探测、分离的方法 |
-
2016
- 2016-08-05 CN CN201610642463.7A patent/CN106291524B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009062773A1 (de) * | 2007-11-15 | 2009-05-22 | Robert Bosch Gmbh | Verfahren zum erkennen von lebewesen |
CN102928835A (zh) * | 2012-10-09 | 2013-02-13 | 北京航空航天大学 | 一种基于改进广义s变换的人体目标运动状态识别方法 |
WO2014120289A1 (en) * | 2012-10-10 | 2014-08-07 | Raytheon Company | Detection of concealed object on a body using radio frequency signatures on frequencies and polarizations |
CN104237865A (zh) * | 2014-10-07 | 2014-12-24 | 电子科技大学 | 一种基于太赫兹雷达回波的人体目标微动信号时频分析方法 |
CN104267394A (zh) * | 2014-10-07 | 2015-01-07 | 电子科技大学 | 高分辨率人体目标运动特征检测方法 |
CN105137423A (zh) * | 2015-09-30 | 2015-12-09 | 武汉大学 | 一种穿墙雷达对多运动目标实时探测、分离的方法 |
Non-Patent Citations (4)
Title |
---|
SEVGI Z. GÜRBÜZ,ET AL: "A Nonlinear-Phase Model-Based Human Detector for Radar", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
ZHILI ZHANG,ET AL: "Feature Extraction and Classification of Human Motions with LFMCW Radar", 《IEEE INTERNATIONAL WORKSHOP ON ELECTROMAGNETICS (IWEM)》 * |
杨杉: "基于雷达回波的人体检测和步态特征提取", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
颜松: "基于雷达回波的人体探测算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107064927A (zh) * | 2017-03-21 | 2017-08-18 | 西北大学 | 一种高速弹载雷达目标相参积累方法 |
CN107907878A (zh) * | 2017-11-08 | 2018-04-13 | 零八电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN107907878B (zh) * | 2017-11-08 | 2021-05-11 | 零八一电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN108537144A (zh) * | 2018-03-21 | 2018-09-14 | 特斯联(北京)科技有限公司 | 一种多维人体步态识别方法与设备 |
CN108537144B (zh) * | 2018-03-21 | 2019-01-29 | 特斯联(北京)科技有限公司 | 一种多维人体步态识别方法与设备 |
CN110673135A (zh) * | 2018-07-03 | 2020-01-10 | 松下知识产权经营株式会社 | 传感器、推定装置、推定方法和程序记录介质 |
CN109917347A (zh) * | 2019-04-10 | 2019-06-21 | 电子科技大学 | 一种基于时频域稀疏重构的雷达行人检测方法 |
CN109917347B (zh) * | 2019-04-10 | 2021-09-24 | 电子科技大学 | 一种基于时频域稀疏重构的雷达行人检测方法 |
CN110501685A (zh) * | 2019-08-23 | 2019-11-26 | 北京电子工程总体研究所 | 一种基于雷达信号相位补偿的多帧相参积累方法 |
CN110488279A (zh) * | 2019-08-30 | 2019-11-22 | 成都纳雷科技有限公司 | 一种基于雷达的行人与树木区分方法与装置 |
CN110488279B (zh) * | 2019-08-30 | 2022-12-09 | 成都纳雷科技有限公司 | 一种基于雷达的行人与树木区分方法与装置 |
CN110584631B (zh) * | 2019-10-10 | 2022-03-25 | 重庆邮电大学 | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 |
CN110584631A (zh) * | 2019-10-10 | 2019-12-20 | 重庆邮电大学 | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 |
CN111160376B (zh) * | 2019-12-31 | 2023-11-24 | 联想(北京)有限公司 | 一种信息处理方法、装置、电子设备和存储介质 |
CN111160376A (zh) * | 2019-12-31 | 2020-05-15 | 联想(北京)有限公司 | 一种信息处理方法、装置、电子设备和存储介质 |
CN113096791B (zh) * | 2020-01-08 | 2024-06-04 | 富士通株式会社 | 基于无线信号的人员识别装置、方法及系统 |
CN113096791A (zh) * | 2020-01-08 | 2021-07-09 | 富士通株式会社 | 基于无线信号的人员识别装置、方法及系统 |
CN111337897A (zh) * | 2020-04-21 | 2020-06-26 | 湖南红船科技有限公司 | 一种lfmcw雷达快速目标识别方法 |
CN111693981A (zh) * | 2020-05-27 | 2020-09-22 | 清华大学 | 变重复周期的二次相关处理方法 |
CN111693981B (zh) * | 2020-05-27 | 2023-10-31 | 清华大学 | 变重复周期的二次相关处理方法 |
CN112639523A (zh) * | 2020-06-30 | 2021-04-09 | 华为技术有限公司 | 一种雷达探测方法及相关装置 |
CN112639523B (zh) * | 2020-06-30 | 2022-04-29 | 华为技术有限公司 | 一种雷达探测方法及相关装置 |
WO2022000333A1 (zh) * | 2020-06-30 | 2022-01-06 | 华为技术有限公司 | 一种雷达探测方法及相关装置 |
CN112327286A (zh) * | 2020-10-13 | 2021-02-05 | 深圳大学 | 低复杂度下日常活动分类方法、装置、设备及存储介质 |
CN112327286B (zh) * | 2020-10-13 | 2023-10-20 | 深圳大学 | 低复杂度下日常活动分类方法、装置、设备及存储介质 |
CN113017590B (zh) * | 2021-02-26 | 2022-04-01 | 清华大学 | 生理数据监测方法、装置、计算机设备和存储介质 |
CN113017590A (zh) * | 2021-02-26 | 2021-06-25 | 清华大学 | 生理数据监测方法、装置、计算机设备和存储介质 |
CN113064131A (zh) * | 2021-03-18 | 2021-07-02 | 江苏蛮酷科技有限公司 | 雷达波处理方法、装置和计算机存储介质 |
CN113064131B (zh) * | 2021-03-18 | 2023-12-26 | 江苏蛮酷科技有限公司 | 雷达波处理方法、装置和计算机存储介质 |
CN116453227B (zh) * | 2023-06-19 | 2023-09-19 | 武汉理工大学 | 一种基于双毫米波雷达的船舶环境下步态识别方法 |
CN116453227A (zh) * | 2023-06-19 | 2023-07-18 | 武汉理工大学 | 一种基于双毫米波雷达的船舶环境下步态识别方法 |
CN117930213A (zh) * | 2024-03-25 | 2024-04-26 | 湘江实验室 | 一种基于毫米波雷达探测人体存在的方法 |
CN117930213B (zh) * | 2024-03-25 | 2024-06-04 | 湘江实验室 | 一种基于毫米波雷达探测人体存在的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106291524B (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291524A (zh) | 基于人体模型的lfmcw雷达探测运动人体目标的方法 | |
CN106646446B (zh) | 脉冲压缩捷变频雷达的动目标检测方法 | |
CN101738606B (zh) | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 | |
CN109143179A (zh) | 一种基于变重频技术的捷变频雷达信号处理方法及装置 | |
CN105158748A (zh) | 一种高速目标多通道补偿聚焦与tbd混合积累检测方法 | |
CN109613506B (zh) | 一种随机跳频重频捷变雷达目标回波信号检测方法 | |
CN104977571B (zh) | 基于俯仰频率分集stap的距离模糊杂波抑制方法 | |
CN104502898B (zh) | 将修正rft和修正mdcft相结合的机动目标参数估计方法 | |
CN106093870A (zh) | 高超声速飞行器下降段的sar‑gmti杂波抑制方法 | |
CN106842166A (zh) | 一种适用于线性调频连续波雷达系统的解速度模糊方法 | |
CN103399310A (zh) | 基于相位差分Radon-Lv分布的雷达微弱动目标检测方法 | |
CN103278820A (zh) | 临近空间慢速平台sar动目标检测方法及成像方法 | |
CN106646447A (zh) | 基于线性调频连续波的雷达目标长时间积累检测方法 | |
CN106338723A (zh) | 一种基于互质脉冲重复间隔的空时自适应处理方法及装置 | |
CN106199549B (zh) | 一种采用谱减法提升lfmcw雷达信噪比的方法 | |
CN108535730A (zh) | 一种多普勒气象雷达解速度模糊方法和系统 | |
CN110133654A (zh) | 一种高轨卫星sar动目标检测方法 | |
CN105116395B (zh) | 一种天基双基地雷达的动目标长时间相参积累方法 | |
CN109613507B (zh) | 一种针对高阶机动目标雷达回波的检测方法 | |
CN105738897A (zh) | 基于联合稀疏特性的信号距离像重构方法 | |
CN106199602A (zh) | 一种Staggered‑SAR回波信号重建方法 | |
CN104375128B (zh) | 一种基于互相关函数的高机动目标快速积累检测方法 | |
CN109031299A (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
CN106597445A (zh) | 一种基于自适应Chirp分解的SAR动目标检测方法 | |
CN107942307B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |