CN103948381A - 一种用于大动态范围多普勒生物信号成像方法 - Google Patents
一种用于大动态范围多普勒生物信号成像方法 Download PDFInfo
- Publication number
- CN103948381A CN103948381A CN201410140529.3A CN201410140529A CN103948381A CN 103948381 A CN103948381 A CN 103948381A CN 201410140529 A CN201410140529 A CN 201410140529A CN 103948381 A CN103948381 A CN 103948381A
- Authority
- CN
- China
- Prior art keywords
- signal
- motion
- amplitude
- dynamic range
- electromagnetic wave
- 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
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种用于大动态范围多普勒生物信号成像方法。射频前端模块通过天线发射电磁波到待测运动生物体表面;下变频模块将反射回来的信号数字正交下变频到基带信号;对基带信号进行动态圆心追踪算法去除无用的直流分量并且保留有用的直流分量,再应用拓展反正切算法解调出的相位信息不受到相位模糊度的限制,从而恢复出被测物体的运动信息。本发明相比于传统的非接触式生物信号探测系统只能探测亚波长级别的小幅度运动,所述的系统具有多功能和应用场景的优点,能突破小幅度运动的限制,应用于数米级别的多个波长大幅度运动的探测,来实现生物医疗等领域的相关非接触大幅度运动探测应用。
Description
技术领域
本发明涉及多普勒生物信号成像方法,尤其是涉及一种用于大动态范围多普勒生物信号成像方法。
背景技术
近几年来,国内外科研工作者在研究使用电磁波进行非接触式心率测量等小幅度运动探测方面,取得了很大的进展。然而,由于正交基带信号中直流偏移不确定的问题,一般反正切运算的值域限制导致相位模糊度问题,探测大动态范围运动,尤其是多个电磁波波长级别的大幅度运动探测受到了极大的限制。
发明内容
本发明目的在于提供一种用于大动态范围多普勒生物信号成像方法。相比于传统的非接触式生物信号探测系统只能探测有限的一个电磁波波长范围内的小幅度运动,本发明具有多功能和多应用场景的优点,能同时应用于数米级别的40倍以上波长的大幅度运动探测,也能与现行的电磁波非接触测量方案应用相同,应用于一个电磁波波长范围内的微小幅度运动的探测。
本发明采用的技术方案是:
射频前端模块通过天线发射电磁波到待测运动生物体表面;下变频模块将反射回来的信号数字正交下变频到基带信号;对基带信号进行动态圆心追踪算法去除无用的直流分量并且保留有用的直流分量,再应用拓展反正切算法解调出的相位信息不受到相位模糊度的限制,从而恢复出被测物体的运动信息。
其具体工作过程如下:
(1)所述射频前端模块对着待测运动生物体表面发射单频载波电磁波信号,经过待测运动生物体表面被反射,射频前端模块接收到反射电磁波信号,并将反射电磁波信号数字正交下变频,生成信号Q和信号I两路,分别满足:
其中,AI(t),AQ(t)分别为信号I和信号Q的幅度;由于数字域的正交解调,忽略幅度不平衡,即AI(t)=AQ(t)=AR(t);θ为一固定相移,与射频前端模块距离待测物体的初始距离有关;x(t)为运动物体的位移信息;λ为载波波长;DCI(t)和DCQ(t)分别为信号I和信号Q的直流偏移;
(2)在小幅度运动测量情况下,信号Q和信号I中的直流偏移变化很小,幅度变化也很小,此时认为不变;因此信号I和信号Q组成了一段以直流偏移(DCI(t),DCQ(t))为圆心、幅度AR(t)为半径的圆弧曲线:
[I(t)-DCI(t)]2+[Q(t)-DCQ(t)]2=AR 2(t)
而当测量大幅度运动时,信号Q和信号I中包含的直流偏移是时变的;此时其构成的圆弧的圆心位置偏移到坐标(DCI(t),DCQ(t));为了消除变化的直流偏移信号对最终成像结果的影响,需要使用动态圆心追踪算法获取坐标(DCI(t),DCQ(t));所述需要使用动态圆心追踪算法获取坐标(DCI(t),DCQ(t));就是将大幅度运动划分为直流偏置和功率不变的小幅度运动,每段小幅度段运动上,对(DCI,DCQ,AR)的方程求最优解:
(3)之后在信号Q和信号I中对其进行补偿消除,获得以下的运动位移信息:
但是由于反三角函数具有(-π/2,π/2)的值域限制,该模块的测量结果中会包含一些不连续点。为提高系统稳定性,使用拓展的反正切算法进行计算。
本发明具有的有益效果是:
相比于传统的非接触式生物信号探测系统只能探测亚波长级别的小幅度运动,所述的系统具有多功能和多应用场景的优点,能突破小幅度运动的限制,应用于数米级别的多个波长大幅度运动的探测,来实现生物医疗等领域的相关非接触大幅度运动探测应用。
附图说明
图1是本发明的结构框图。
图2(a)是15ms内,本发明对音叉震动探测的时域上结果图。
图2(b)是本发明对音叉震动探测的频域上结果图。
图3(a)是20s时间内人体五米来回走动的信号I和信号Q星座图。
图3(b)是20s时间内人体五米来回走动恢复的运动位移图和相应的速度图。
具体实施方式
以下结合附图,具体阐述本发明的工作原理和实施方式:
如图1所示,射频前端模块包括发射天线TX和接收天线RX。首先通过天线TX发射电磁波到待测运动生物体表面;接收天线RX接收到反射信号,然后下变频模块将反射回来的信号数字正交下变频到两路基带信号:信号Q和信号I;对基带信号进行动态圆心追踪算法去除无用的直流分量并且保留有用的直流分量,再应用拓展反正切算法解调出的相位信息不受到相位模糊度的限制,从而恢复出被测物体的运动信息。
其具体工作过程如下:
1)所述射频前端模块对着待测运动生物体表面发射单频载波电磁波信号,经过待测运动生物体表面被反射,射频前端模块接收到反射电磁波信号,并将反射电磁波信号数字正交下变频,生成信号Q和信号I两路,分别满足:
其中,AI(t),AQ(t)分别为信号I和信号Q的幅度;由于数字域的正交解调,忽略幅度不平衡,即AI(t)=AQ(t)=AR(t);θ为一固定相移,与射频前端模块距离待测物体的初始距离有关;x(t)为运动物体的位移信息;λ为载波波长;DCI(t)和DCQ(t)分别为信号I和信号Q的直流偏移;
(2)在小幅度运动测量情况下,信号Q和信号I中的直流偏移变化很小,幅度变化也很小,此时认为不变;因此信号I和信号Q组成了一段以直流偏移(DCI(t),DCQ(t))为圆心、幅度AR(t)为半径的圆弧曲线:
[I(t)-DCI(t)]2+[Q(t)-DCQ(t)]2=AR 2(t)
而当测量大幅度运动时,信号Q和信号I中包含的直流偏移是时变的;此时其构成的圆弧的圆心位置偏移到坐标(DCI(t),DCQ(t));为了消除变化的直流偏移信号对最终成像结果的影响,需要使用动态圆心追踪算法获取坐标(DCI(t),DCQ(t));
(3)之后在信号Q和信号I中对其进行补偿消除,获得以下的运动位移信息:
但是由于反三角函数具有(-π/2,π/2)的值域限制,该模块的测量结果中会包含一些不连续点。为提高系统稳定性,使用拓展的反正切算法进行计算。
所述需要使用动态圆心追踪算法获取坐标(DCI(t),DCQ(t));就是将大幅度运动划分为直流偏置和功率不变的小幅度运动,每段小幅度段运动上,对(DCI,DCQ,AR)的方程求最优解:
图2是在发射单频18GHz电磁波情况下对一个波长范围内小幅度音叉震动的探测结果:
图2(a)是15ms内本发明对音叉震动探测的时域上结果图。上图是应用本发明方法之前,单路信号音叉震动在时域上的变化,下图是应用本发明方法之后,音叉震动在时域上的变化,可以探测到音叉运动的幅度为285μm。
图2(b)是本发明对音叉震动探测的频域上结果图。探测到音叉的震动频率为425.4Hz,与生产说明书上的424Hz基本一致。
图3是在发射2.36GHz电磁波情况下,对20s时间内人体来回走动的探测结果:
图3(a)是人体五米来回走动的信号I和信号Q星座图。图3(b)上图是恢复的运动位移图,下图是相应恢复运动的速度图。从A点开始为正向的匀速直线运动行走,在B时刻折返,C到D时刻人进行了一次跳跃,所以速度忽然变大。可以清楚的看出整个过程的运动都能被恢复辨识。
Claims (3)
1.一种用于大动态范围多普勒生物信号成像方法,其特征在于:射频前端模块通过天线发射电磁波到待测运动生物体表面;下变频模块将反射回来的信号数字正交下变频到基带信号;对基带信号进行动态圆心追踪算法去除无用的直流分量并且保留有用的直流分量,再应用拓展反正切算法解调出的相位信息不受到相位模糊度的限制,从而恢复出被测物体的运动信息。
2.根据权利要求1所述的一种用于大动态范围多普勒生物信号成像方法,其特征在于,其具体工作过程如下:
(1)所述射频前端模块对着待测运动生物体表面发射单频载波电磁波信号,经过待测运动生物体表面被反射,射频前端模块接收到反射电磁波信号,并将反射电磁波信号数字正交下变频,生成信号Q和信号I两路,分别满足:
其中,AI(t),AQ(t)分别为信号I和信号Q的幅度;由于数字域的正交解调,忽略幅度不平衡,即AI(t)=AQ(t)=AR(t);θ为一固定相移,与射频前端模块距离待测物体的初始距离有关;x(t)为运动物体的位移信息;λ为载波波长;DCI(t)和DCQ(t)分别为信号I和信号Q的直流偏移;
(2)在小幅度运动测量情况下,信号Q和信号I中的直流偏移变化很小,幅度变化也很小,此时认为不变;因此信号I和信号Q组成了一段以直流偏移(DCI(t),DCQ(t))为圆心、幅度AR(t)为半径的圆弧曲线:
[I(t)-DCI(t)]2+[Q(t)-DCQ(t)]2=AR 2(t)
而当测量大幅度运动时,信号Q和信号I中包含的直流偏移是时变的;此时其构成的圆弧的圆心位置偏移到坐标(DCI(t),DCQ(t));为了消除变化的直流偏移信号对最终成像结果的影响,需要使用动态圆心追踪算法获取坐标(DCI(t),DCQ(t));
(3)之后在信号Q和信号I中对其进行补偿消除,获得以下的运动物体的位移信息x(t):
但是由于反三角函数具有(-π/2,π/2)的值域限制,该模块的测量结果中会包含一些不连续点。为提高系统稳定性,使用拓展的反正切算法进行计算。
3.根据权利要求2所述的一种用于大动态范围多普勒生物信号成像方法,其特征在于:所述需要使用动态圆心追踪算法获取坐标(DCI(t),DCQ(t));就是将大幅度运动划分为直流偏置和功率不变的小幅度运动,每段小幅度段运动上,对(DCI,DCQ,AR)的方程求最优解:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410140529.3A CN103948381B (zh) | 2014-04-09 | 2014-04-09 | 一种用于大动态范围多普勒生物信号成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410140529.3A CN103948381B (zh) | 2014-04-09 | 2014-04-09 | 一种用于大动态范围多普勒生物信号成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103948381A true CN103948381A (zh) | 2014-07-30 |
CN103948381B CN103948381B (zh) | 2015-11-04 |
Family
ID=51325843
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410140529.3A Expired - Fee Related CN103948381B (zh) | 2014-04-09 | 2014-04-09 | 一种用于大动态范围多普勒生物信号成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103948381B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016045455A1 (zh) * | 2014-09-28 | 2016-03-31 | 成都维客亲源健康科技有限公司 | 适用于可穿戴设备的高可靠低计算量的心律识别电路与方法 |
CN106405520A (zh) * | 2016-09-30 | 2017-02-15 | 浙江大学 | 基于多通道连续波多普勒雷达的物体运动模式识别方法 |
CN108852327A (zh) * | 2018-04-16 | 2018-11-23 | 浙江大学 | 一种从运动干扰中非接触检测微弱生命信号的方法 |
CN110115585A (zh) * | 2019-05-06 | 2019-08-13 | 浙江大学 | 一种心动图的非接触测量方法 |
CN110133642A (zh) * | 2019-05-06 | 2019-08-16 | 浙江大学 | 一种用于雷达传感器的直流偏移补偿方法 |
CN110987150A (zh) * | 2019-12-16 | 2020-04-10 | 上海交通大学 | 一种基于单频连续波雷达的超微幅振动测量方法与系统 |
CN111308465A (zh) * | 2020-03-26 | 2020-06-19 | 中国气象局气象探测中心 | 一种基于多普勒速度方位滤波识别空中生物的方法 |
CN112741618A (zh) * | 2020-12-22 | 2021-05-04 | 浙江大学 | 一种基于fmcw雷达的舌姿检测系统与方法 |
CN113589249A (zh) * | 2021-07-16 | 2021-11-02 | 中山艾朗格科技有限公司 | 校准单频连续波多普勒雷达的直流偏移的信号处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4513748A (en) * | 1983-08-30 | 1985-04-30 | Rca Corporation | Dual frequency heart rate monitor utilizing doppler radar |
WO1999031525A1 (en) * | 1997-12-15 | 1999-06-24 | Milkovich Systems Engineering | Signal processing architecture which improves sonar and pulse doppler radar performance and tracking capability |
KR20120000398A (ko) * | 2010-06-25 | 2012-01-02 | 한국생산기술연구원 | 비접촉식 생체 신호 측정 장치 및 그 장치에서의 비접촉식 생체 신호 측정 방법 |
CN102508256A (zh) * | 2011-11-23 | 2012-06-20 | 中国科学院微电子研究所 | 非接触式实时生理信号监测系统 |
-
2014
- 2014-04-09 CN CN201410140529.3A patent/CN103948381B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4513748A (en) * | 1983-08-30 | 1985-04-30 | Rca Corporation | Dual frequency heart rate monitor utilizing doppler radar |
WO1999031525A1 (en) * | 1997-12-15 | 1999-06-24 | Milkovich Systems Engineering | Signal processing architecture which improves sonar and pulse doppler radar performance and tracking capability |
KR20120000398A (ko) * | 2010-06-25 | 2012-01-02 | 한국생산기술연구원 | 비접촉식 생체 신호 측정 장치 및 그 장치에서의 비접촉식 생체 신호 측정 방법 |
CN102508256A (zh) * | 2011-11-23 | 2012-06-20 | 中国科学院微电子研究所 | 非接触式实时生理信号监测系统 |
Non-Patent Citations (7)
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016045455A1 (zh) * | 2014-09-28 | 2016-03-31 | 成都维客亲源健康科技有限公司 | 适用于可穿戴设备的高可靠低计算量的心律识别电路与方法 |
CN106405520A (zh) * | 2016-09-30 | 2017-02-15 | 浙江大学 | 基于多通道连续波多普勒雷达的物体运动模式识别方法 |
CN108852327B (zh) * | 2018-04-16 | 2020-06-19 | 浙江大学 | 一种从运动干扰中非接触检测微弱生命信号的方法 |
CN108852327A (zh) * | 2018-04-16 | 2018-11-23 | 浙江大学 | 一种从运动干扰中非接触检测微弱生命信号的方法 |
CN110115585A (zh) * | 2019-05-06 | 2019-08-13 | 浙江大学 | 一种心动图的非接触测量方法 |
CN110133642A (zh) * | 2019-05-06 | 2019-08-16 | 浙江大学 | 一种用于雷达传感器的直流偏移补偿方法 |
CN110133642B (zh) * | 2019-05-06 | 2020-07-10 | 浙江大学 | 一种用于雷达传感器的直流偏移补偿方法 |
CN110115585B (zh) * | 2019-05-06 | 2020-07-10 | 浙江大学 | 一种心动图的非接触测量方法 |
WO2020224600A1 (zh) * | 2019-05-06 | 2020-11-12 | 浙江大学 | 一种心动图的非接触测量方法 |
CN110987150A (zh) * | 2019-12-16 | 2020-04-10 | 上海交通大学 | 一种基于单频连续波雷达的超微幅振动测量方法与系统 |
CN110987150B (zh) * | 2019-12-16 | 2021-07-23 | 上海交通大学 | 一种基于单频连续波雷达的超微幅振动测量方法与系统 |
CN111308465A (zh) * | 2020-03-26 | 2020-06-19 | 中国气象局气象探测中心 | 一种基于多普勒速度方位滤波识别空中生物的方法 |
CN112741618A (zh) * | 2020-12-22 | 2021-05-04 | 浙江大学 | 一种基于fmcw雷达的舌姿检测系统与方法 |
CN113589249A (zh) * | 2021-07-16 | 2021-11-02 | 中山艾朗格科技有限公司 | 校准单频连续波多普勒雷达的直流偏移的信号处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103948381B (zh) | 2015-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103948381A (zh) | 一种用于大动态范围多普勒生物信号成像方法 | |
Wang et al. | TagScan: Simultaneous target imaging and material identification with commodity RFID devices | |
US9645234B2 (en) | RFID device, methods and applications | |
CN103793059A (zh) | 一种基于时域多普勒手势恢复识别方法 | |
CN105786185A (zh) | 基于连续波微多普勒雷达的非接触式手势识别系统与方法 | |
Gu et al. | A two-tone radar sensor for concurrent detection of absolute distance and relative movement for gesture sensing | |
Gao et al. | Barcode based hand gesture classification using AC coupled quadrature Doppler radar | |
CN104375145A (zh) | 一种毫米波成像方法及系统 | |
Zheng et al. | Doppler bio-signal detection based time-domain hand gesture recognition | |
CN110133642B (zh) | 一种用于雷达传感器的直流偏移补偿方法 | |
CN114296141A (zh) | 多目标生命征象侦测器及其侦测方法 | |
CN106338716A (zh) | 基于三坐标雷达技术的民用无人机捕获技术及其系统 | |
Hsu et al. | Wearable vital sign sensor using a single-input multiple-output self-injection-locked oscillator tag | |
Fan et al. | Hand gesture recognition based on Wi-Fi chipsets | |
CN106323450A (zh) | 一种基于多普勒雷达的大型柔性结构的振动监测方法 | |
Lai et al. | Finger gesture sensing and recognition using a Wi-Fi-based passive radar | |
CN105699992B (zh) | 高动态gnss载波信号跟踪方法及系统 | |
Wang et al. | Review on microwave/millimeter-wave systems for vital sign detection | |
CN104704386A (zh) | 无线传感器装置 | |
Fairchild et al. | Determining human target facing orientation using bistatic radar micro-Doppler signals | |
Fathy et al. | Overview of human vital signs detection using radar techniques | |
JP2009210372A (ja) | 距離測定装置 | |
Li et al. | Non-contact multi-radar smart probing of body orientation based on micro-Doppler signatures | |
Sippel et al. | Quasi-coherent phase-based localization and tracking of incoherently transmitting radio beacons | |
Zhai et al. | An effective velocity detection method for moving uhf-rfid tags |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151104 |