CN106405520B - 基于多通道连续波多普勒雷达的物体运动模式识别方法 - Google Patents
基于多通道连续波多普勒雷达的物体运动模式识别方法 Download PDFInfo
- Publication number
- CN106405520B CN106405520B CN201610874009.4A CN201610874009A CN106405520B CN 106405520 B CN106405520 B CN 106405520B CN 201610874009 A CN201610874009 A CN 201610874009A CN 106405520 B CN106405520 B CN 106405520B
- Authority
- CN
- China
- Prior art keywords
- movement
- antenna
- arcsine
- under test
- phase
- 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
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
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
本发明公开了一种基于多通道连续波多普勒雷达的物体运动模式识别方法。通过1个天线发射电磁波,照射待测运动物体表面,被反射的信号经2‑3个接收天线接收并下变频到基带,用反正弦算法,分别解调出不受相位模糊度限制的相位信息,得到每个接收天线与运动物体间的距离变化信息,再通过追踪公式得到物体的运动轨迹,在运动轨迹的基础上进行模式识别。本发明避免反正切类函数的相位模糊问题,适用高采样率和低采样率的工作条件,并在低采样率的条件下保持一定的准确度。实现在对待测物体定位和跟踪的基础上进行模式识别。识别准确率高,抗干扰能力强,架构简单、成本低,解调出来的运动大多为线性关系,无需复杂数据处理过程,节省硬件资源。
Description
技术领域
本发明涉及非接触式对物体运动模式识别的方法,尤其是涉及一种基于多通道连续波多普勒雷达的物体运动模式识别方法。
背景技术
使用连续波多普勒雷达或雷达传感器技术进行运动检测近年来取得了很大的进展。连续波雷达传感器采用单一频率工作,具有结构简单、成本低、可集成性好、抗干扰能力强等突出优点。此前其所用的基带解调算法一般采用小角近似原理或反正切运算,无法实现低采样率场景下的大动态范围运动的测量。另外,此前的连续波多普勒雷达或雷达传感器多用于测量一维运动,并且很少用于运动模式识别方面的应用。
发明内容
为了解决背景技术中存在的问题,本发明的目的在于提供一种基于多通道连续波多普勒雷达的物体运动模式识别方法。上述方法针对解调得到的基带信号,因此适用于所有采用零中频、副载波调制、低中频、超外差或数字中频架构的雷达或雷达传感器。
本发明采用的技术方案是:
本发明通过一个天线发射电磁波,照射待测运动物体表面,被反射的信号经2-3个接收天线接收并下变频到基带,再应用反正弦算法,即Arcsine算法,分别解调出不受相位模糊度限制的相位信息,得到每个接收天线与待测物体间的距离变化信息,再通过提出的追踪公式得到物体的运动轨迹,在运动轨迹的基础上进行模式识别。
所述Arcsine算法是:
对于每个接收天线的任意正交下变频接收机,其输出的正交I、Q信号经采样后,待测物体不受相位模糊度限制的运动相位信息Φ[n]表示为:
式中n表示对I、Q信号的采样点数,当采样率足够高时,即采样率大于20v/λ,其中v为待测物体的最大运动速度,λ为发射电磁波的波长,上式可以近似等于下式:
Arcsine算法避免了反正切类函数的相位模糊问题,同时在高采样率即采样率大于20v/λ和低采样率即采样率小于20v/λ但大于4v/λ的条件下保持准确。
当待测物体运动在三维空间中时,使用所述的1个发射天线和3个接收天线,使用如下追踪公式对待测物体在三维空间中的运动轨迹进行记录
其中(xT,yT,zT)为发射天线坐标,(x(t),y(t),z(t))为目标物体的实时坐标,(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)分别为3个接收天线的坐标,d1、d2、d3分别是电磁波从发射天线发射经待测物体反射后再到各个接收天线走过的距离。
当待测物体运动在二维平面上时,使用所述的1个发射天线和2个接收天线对待测物体的运动轨迹进行记录;以物体运动平面作为xoy平面,将1个发射天线放置于坐标(0,0,zT),2个接收天线分别放置于坐标(x1,0,zT)、(x2,0,zT),则上面公式简化为:
其另一表达方式为:
本发明具有的有益效果是:
本发明的反正弦(Arcsine)算法避免了反正切类函数的相位模糊问题,同时适用于高采样率和低采样率的工作条件,并在低采样率的条件下保持准确。综合来说,本发明实现了在对待测物体定位和跟踪的基础上进行模式识别。具有识别准确率高,抗干扰能力强,架构简单、成本低的优点,解调出来的运动大多为线性关系,无需复杂数据处理过程,节省硬件资源。
附图说明
图1是本发明Arcsine算法的图解。
图2是1个发射天线和3个接收天线追踪物体运动轨迹的立体图。
图3是实施例Arcsine算法和直接反正切算法以及反正切加相位解模糊算法的对比结果。
图4是实施例使用Arcsine解调待测物体单击和移动两种运动模式下的距离变化信息的实验结果。
图5是实施例在轨迹追踪基础上识别待测物体单击和移动两种运动模式的实验结果。
具体实施方式
下面结合附图和实施例对本发明做进一步的说明。
本发明的工作原理和实施方式:
射频前端通过1个发射天线发射电磁波,照射待测运动物体表面,反射的信号经两个接收天线接收并下变频到基带。对于任意正交下变频接收机,其输出的正交I、Q信号经采样后,待测物体的距离变化信息x(t)可以表示为:
其中DCI(t)和DCQ(t)是直流偏移,AI(t)和AQ(t)是I、Q信号的幅度(正交解调时可认为两者相等),λ、θ0分别表示电磁波的波长、收发机的残余相位噪声(在相关解调中可认为等于0)和收发信号的额外相差。所以x(t)可由下式得到
由于反三角函数具有(-π/2,π/2)的值域限制,得到的x(t)会包含不连续点,极大地影响测量精度和范围。为提高系统的鲁棒性,将接收到的I、Q信号在去掉直流偏移后,在I-Q星座图上表示,如图1所示,记向量 和分别表示第k次和第k-1次I、Q信号采样结果。所以两个向量夹角可以写为下式,其中是两个向量所构成平面的法向量。
对ΔΦk累积求和,得:
上式即为无相位模糊度的解调结果,其中n表示对I、Q信号的采样点数。当和足够接近,即采样率足够高时,上式可以近似等于下式:
通过反正弦(Arcsine)算法得到的相位信息再乘以系数得到物体的运动引起的d1、d2的变化量,结合初始时刻d1、d2的值d01、d02,即可得到每个时刻d1、d2的值。测量d01、d02的技术方案为:发射双频副载波到待测物体,接收机分别对两个频率的信号进行解调。以d01为例,双频副载波两个频率的信号经过反正弦(Arcsine)算法解调后各自的相位分别为Φ1、Φ2,并满足以下关系:
当将上面两式的相位之差控制在一个模糊度范围内时,即k1=k2,将上两式相减得到:
对于空间中任意放置的1个发射天线和3个接收天线,如图2所示,使用如下公式(或其扩展、化简或等效的公式)对目标在三维空间中的运动模式进行记录。
其中(xT,yT,zT)为发射天线坐标,(x(t),y(t),z(t))为目标物体的实时坐标,(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)分别为3个接收天线的坐标,d1、d2、d3分别是电磁波从发射天线发射经目标反射后再到各个接收天线走过的距离。
当待测物体运动在二维平面上时,使用1个发射天线和2个接收天线对物体的运动模式进行记录和进行运动模式识别。以物体运动平面作为xoy平面,将1个发射天线和2个接收天线放置在一个平行于xoy平面的直线上,例如发射天线坐标为(0,0,zT),2个接收天线坐标分别为(x1,0,zT)、(x2,0,zT),则上面的公式简化为:
其另一表达方式为:
得到轨迹后将其与已构建的模型进行特征匹配来进行识别。特征匹配方法有直接匹配法、动态时间规整法、隐性马尔科夫模型(HHM)法、神经网络模型法等。本发明实施例选用的是基于概率统计的隐性马尔科夫模型(HHM)法,它尤其适用时间序列的建模,对复杂度高的动作也具有很高的识别精度,易于添加或修改手势库。使用此种方法首先根据运动轨迹进行模式分类,然后开始训练,为每一种模式建立一个HHM模型,识别时取概率最大的一个HHM即可。
图3是反正弦(Arcsine)算法和直接反正切算法以及反正切加相位解模糊算法在有噪声情况下的解调效果对比图,可以看出,Arcsine算法相位连续性最好,其余两种算法均出现了不同程度的相位跳变。
图4是使用Arcsine解调待测物体单击和移动两种运动模式下的距离变化信息的实验结果。
图5是在轨迹追踪基础上识别待测物体单击和移动两种运动模式的实验结果,从图4和图5中可以清楚地分辨待测物体的两种运动模式。
Claims (1)
1.基于多通道连续波多普勒雷达的物体运动模式识别方法,其特征在于:通过1个发射天线发射电磁波,照射待测运动物体表面,被反射的信号经2-3个接收天线接收并下变频到基带,再应用反正弦算法,即Arcsine算法,分别解调出不受相位模糊度限制的相位信息,得到每个接收天线与待测物体间的距离变化信息,再通过提出的追踪公式得到物体的运动轨迹,在运动轨迹的基础上进行模式识别;
所述Arcsine算法是:
对于每个接收天线的任意正交下变频接收机,其输出的正交I、Q信号经采样后,待测物体不受相位模糊度限制的运动相位信息Φ[n]表示为
式中n表示对I、Q信号的采样点数,当采样率足够高时,即采样率大于20v/λ,其中v为待测物体的最大运动速度,λ为发射电磁波的波长,上式可以近似等于下式
Arcsine算法避免了反正切类函数的相位模糊问题,同时在高采样率即采样率大于20v/λ和低采样率即采样率小于20v/λ但大于4v/λ的条件下保持准确;
当待测物体运动在三维空间中时,使用所述的1个发射天线和3个接收天线,使用如下追踪公式对待测物体在三维空间中的运动轨迹进行记录
其中(xT,yT,zT)为发射天线坐标,(x(t),y(t),z(t))为目标物体的实时坐标,(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)分别为3个接收天线的坐标,d1、d2、d3分别是电磁波从发射天线发射经待测物体反射后再到各个接收天线走过的距离;
当待测物体运动在二维平面上时,使用所述的1个发射天线和2个接收天线对待测物体的运动轨迹进行记录;以待测物体的运动平面作为xoy平面,将1个发射天线放置于坐标(0,0,zT),2个接收天线分别放置于坐标(x1,0,zT)、(x2,0,zT),则公式简化为
其另一表达方式为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610874009.4A CN106405520B (zh) | 2016-09-30 | 2016-09-30 | 基于多通道连续波多普勒雷达的物体运动模式识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610874009.4A CN106405520B (zh) | 2016-09-30 | 2016-09-30 | 基于多通道连续波多普勒雷达的物体运动模式识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106405520A CN106405520A (zh) | 2017-02-15 |
CN106405520B true CN106405520B (zh) | 2018-11-13 |
Family
ID=59229597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610874009.4A Active CN106405520B (zh) | 2016-09-30 | 2016-09-30 | 基于多通道连续波多普勒雷达的物体运动模式识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106405520B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108363043A (zh) * | 2018-01-26 | 2018-08-03 | 浙江大学 | 分布放置连续波多普勒雷达传感器及多运动目标探测方法 |
CN108490427A (zh) * | 2018-02-07 | 2018-09-04 | 浙江大学 | 一种运动目标室内定位及实时追踪方法 |
CN108896984B (zh) * | 2018-04-20 | 2020-11-24 | 浙江大学 | 一种基于射频单频连续波的运动物体定位跟踪方法 |
CN110163130B (zh) * | 2019-05-08 | 2021-05-28 | 清华大学 | 一种用于手势识别的特征预对齐的随机森林分类系统及方法 |
CN110133584B (zh) * | 2019-05-16 | 2023-04-21 | 中国科学技术大学 | 一种基于无源反射信号高精度实时手写轨迹追踪方法 |
TWI730436B (zh) * | 2019-10-02 | 2021-06-11 | 立積電子股份有限公司 | 辨識物體的方法、物體辨識系統及更新物體辨識系統中主要分類器的方法 |
CN113567947B (zh) * | 2021-07-21 | 2023-07-28 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 一种多普勒相移信号的解调方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9335825B2 (en) * | 2010-01-26 | 2016-05-10 | Nokia Technologies Oy | Gesture control |
CN103793059A (zh) * | 2014-02-14 | 2014-05-14 | 浙江大学 | 一种基于时域多普勒手势恢复识别方法 |
CN103948381B (zh) * | 2014-04-09 | 2015-11-04 | 浙江大学 | 一种用于大动态范围多普勒生物信号成像方法 |
US9817109B2 (en) * | 2015-02-27 | 2017-11-14 | Texas Instruments Incorporated | Gesture recognition using frequency modulated continuous wave (FMCW) radar with low angle resolution |
CN105824020B (zh) * | 2016-03-12 | 2018-06-12 | 浙江大学 | 副载波调制的连续波多普勒雷达传感器及运动解调方法 |
CN105786185B (zh) * | 2016-03-12 | 2019-01-18 | 浙江大学 | 基于连续波微多普勒雷达的非接触式手势识别系统与方法 |
-
2016
- 2016-09-30 CN CN201610874009.4A patent/CN106405520B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106405520A (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106405520B (zh) | 基于多通道连续波多普勒雷达的物体运动模式识别方法 | |
CN111656217B (zh) | 用于虚拟孔径雷达跟踪的系统和方法 | |
Miesen et al. | UHF RFID localization based on synthetic apertures | |
Molchanov et al. | Short-range FMCW monopulse radar for hand-gesture sensing | |
CN106918807B (zh) | 一种雷达回波数据的目标点迹凝聚方法 | |
Falcone et al. | Two‐dimensional location of moving targets within local areas using WiFi‐based multistatic passive radar | |
US7928896B2 (en) | Application of time reversal to synthetic aperture imaging | |
CN105005040B (zh) | 一种雷达测角的方法 | |
CN105487060B (zh) | 一种双通道四斜率调制的多目标提取方法 | |
WO2012094825A1 (zh) | 采用射频识别阅读器定位标签的方法以及射频识别阅读器 | |
CN103353597A (zh) | 一种用于超高频rfid定位的相位式测距方法 | |
CN103901425B (zh) | 一种测量运动物体横向速度的雷达及方法 | |
CN108562879A (zh) | 基于fpga的舰载雷达恒虚警检测方法 | |
JP6192151B2 (ja) | 信号選別装置、信号選別方法、及びレーダ装置。 | |
Hyun et al. | Human-vehicle classification scheme using doppler spectrum distribution based on 2D range-doppler FMCW radar | |
CN110857979B (zh) | 一种基于相参信号的极化频选特性材料抗干扰辨识系统及方法 | |
CN109541579B (zh) | 基于Bezier模型的霍夫变换的多普勒穿墙雷达定位方法 | |
CN109507654B (zh) | 一种基于ls的复杂环境下相位信息计算方法 | |
Sim et al. | Road environment recognition for automotive FMCW radar systems through convolutional neural network | |
Lee et al. | Digit recognition in air-writing using single millimeter-wave band radar system | |
Liu et al. | Suppressing coupling and stationary clutters in FMCW radars with temporal filtering | |
Endo et al. | Multi-person position estimation based on correlation between received signals using mimo fmcw radar | |
KR20200114467A (ko) | 레이다 표적 검출 연산량 감소 방법 및 장치 | |
CN114545384A (zh) | 一种应用于跑道异物检测的整合方法及系统 | |
WO2012114816A1 (ja) | 信号識別装置、信号識別方法、およびレーダ装置 |
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 |