CN111142084B - 一种微动太赫兹谱识别检测算法 - Google Patents
一种微动太赫兹谱识别检测算法 Download PDFInfo
- Publication number
- CN111142084B CN111142084B CN201911291687.8A CN201911291687A CN111142084B CN 111142084 B CN111142084 B CN 111142084B CN 201911291687 A CN201911291687 A CN 201911291687A CN 111142084 B CN111142084 B CN 111142084B
- Authority
- CN
- China
- Prior art keywords
- terahertz
- detection algorithm
- time
- frequency
- micromotion
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 30
- 238000001228 spectrum Methods 0.000 title claims abstract description 26
- 230000001629 suppression Effects 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 238000001035 drying Methods 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 230000001678 irradiating effect Effects 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 claims description 9
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000035508 accumulation Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 206010053615 Thermal burn Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
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/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/418—Theoretical aspects
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
本发明属于太赫兹探测技术领域,涉及一种微动太赫兹谱识别检测算法。本发明的检测算法包括:采用MIMO太赫兹雷达照射目标;对得到的两路I/Q信号进行静态分量抑制;对每一路I/Q信号进行反正切解调;对解调后的信号进行互相干处理;对互相干处理后的信号进行时‑频分析;基于时‑频分析结果,重构回波信号;对重构的回波信号进行时域滤波、变时窗分解;估计目标的太赫兹多普勒普。本发明的检测算法,可以在非视距环境下实现对目标太赫兹多普勒谱的检测,进一步拓宽了太赫兹探测技术的应用范围和领域。
Description
技术领域
本发明属于太赫兹探测技术领域,涉及一种微动太赫兹谱识别检测算法。
背景技术
太赫兹波在电磁波谱中的特殊位置及物理特殊性,使其在公共安全、生物医学、通信、军用雷达等领域具有独特的应用价值和广阔的市场前景。太赫兹雷达是太赫兹波诸多应用中最重要的研究方向之一。太赫兹雷达具有频率高、方向性强等特点,上述特点赋予了太赫兹雷达巨大的应用潜力,成为世界各主要国家的研究热点。相比于激光红外雷达,太赫兹雷达对非极性和弱极性介质材料具有更强的透射能力,烟尘、雾霾等恶劣环境下的目标检测性能基本不受影响,可实现全天候工作,同时回波具有更好的相参性,具备更高的抗干扰能力。相比于微波雷达,太赫兹雷达波长较小,低径向速度的目标在太赫兹雷达中可产生更大的多普勒频移,能够实现慢速运动或蠕动目标的检测,且更易实现大信号带宽和窄天线波束,非常有利于获取目标特征结构细节,实现目标的高精度识别。
基于MIMO雷达的太赫兹谱检测技术具有广泛应用,如在军事目标预警监测应用中,能够服务于高超声速目标的探测、反导拦截等领域,实现隐身军事装备的反制。在安检中,能够服务于人员分布和活动情况的侦查,实现对无驻留人体的高精度实时成像。在灾后救援应用中,能够服务于抗震救灾中废墟下人体的寻找,实现伤员的及时搜救。在生物医学应用中,能够实现非接触式生命信号监测,实现对大面积烧烫伤病人、婴幼儿、防区监控范围内的流动人群呼吸和心跳进行监控,排查安全隐患。然而,目前大部分太赫兹谱检测技术多是视距环境下,非视距环境下很难实现目标太赫兹多普勒谱的检测。
发明内容
本发明的目的是解决现有技术中存在的问题,提供一种非视距环境下太赫兹多普勒普的检测方法。
本发明解决其技术问题采用的技术方案是:一种微动太赫兹谱识别检测算法,包括:
采用MIMO太赫兹雷达照射目标;
对得到的两路I/Q信号进行静态分量抑制;
对每一路I/Q信号进行反正切解调;
对解调后的信号进行互相干处理;
对互相干处理后的信号进行时-频分析;
基于时-频分析结果,重构回波信号;
对重构的回波信号进行时域滤波、变时窗分解;
估计目标的太赫兹多普勒普。
进一步地,所述的静态分量抑制采用加权距离轮廓方法分别对两路I/Q信号进行静态分量抑制,利用空间背景估计方法实现背景估计,表示为:
Q1(i+1)=λQ1(i+1)+(1-λ)Q1(i),i=1,…,N (2)
Q2(i+1)=λQ2(i+1)+(1-λ)Q2(i),i=1,…,N (4)
其中,N表示数字化回波信号的样本点数,λ表示加权因子,取值范围为0~1。
进一步地,对每一路I/Q信号进行反正切解调,表示为:
进一步地,互相干处理表示为:
W=Ω1*Ω2 (7)
其中,*表示互相干运算。
进一步地,利用短时傅里叶变换对互相干处理结果进行时-频分析,表示为:
其中,p表示3-9Hz范围内的频率分量;
Ξ表示汉明窗函数,表示为:
其中,α+β=1。
进一步地,重构的回波信号表示为:
进一步地,利用均值滤波器实现信号增强与噪声抑制,表示为:
进一步地,变时窗分解包括:
将信号分为ψ组,第一组点数r,第二组点数r+o,第三组点数r+2o,…,以此类推,第ψ组包括点数N;其中:r=2n,且r的采样时长至少为3~5倍的呼吸或心跳周期,n为正整数;o为正整数;
在每组数据后加0,直至每组数据具有N个值;
对变时窗得到的每组数据进行傅里叶变换,并对其进行加窗处理,表示为:
Tj[n]=χ[n]{FFT{Sj[n]}}j=1,....,Ψ;n∈K*;K*=k*,k*+1,...,k*+κ-1 (12)
其中,FFT{Sj[n]}表示第j个时窗内的信号的傅里叶变换,k*表示加窗频率下限;
对于得到的具有相同频谱分辨率的ψ组频域特征做如下处理:
对信号进行多重频率累积,表示为:
X[n]=(Π[n]+jΠ[n])*(Π[n]+jΠ[n]) (14)
进一步地,所述MIMO太赫兹雷达的发射信号功率范围为1mW-1W,带宽93.2-96.8GHz,信号模式可为有效带宽内的任意的点频连续波、线性调频连续波或步进频连续波。
本发明的检测算法,可以在非视距环境下实现对目标太赫兹多普勒谱的检测,进一步拓宽了太赫兹探测技术的应用范围和领域。
附图说明
图1是本发明实施例提供的方法的示意图;
图2是本发明实施例的算法流程图;
图3是解调后的回波信号;
图4是互相干信号示意图;
图5是时频分析示意图;
图6是重构回波信号示意图;
图7是均值滤波结果示意图;
图8是变时窗分解方法示意图;
图9是多普勒谱结果示意图。
具体实施方式
为了便于理解本发明,下面结合附图和具体实施例,对本发明进行更详细的说明。附图中给出了本发明的较佳的实施例。但是,本发明可以以许多不同的形式来实现,并不限于本说明书所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
本实施例的一种微动太赫兹谱识别检测算法,如图1所示,MIMO太赫兹雷达发射96GHz单频连续波,以人体为目标,人体与雷达间的距离为6m,两者之间以木门为间隔,检测算法流程如图2所示,具体包括:
MIMO雷达同时得到两路I/Q信号,人体目标回波信号中通常包括静态分量,信号较强,导致多普勒信息被淹没,因此,本实施例先利用加权距离轮廓算法分别对两路I/Q信号进行静态分量抑制,可表示为:
利用空间背景估计算法实现背景估计,表示为:
Q1(i+1)=λQ1(i+1)+(1-λ)Q1(i),i=1,…,N (2)
Q2(i+1)=λQ2(i+1)+(1-λ)Q2(i),i=1,…,N (4)
其中,N表示数字化回波信号的样本点数,λ表示加权因子,取值范围为0-1,本实施例中取值为0.75。
对每一路I/Q信号进行反正切解调,可表示为:
解调后的回波信号,以Ω1为例,如图3所示。
对解调后的信号Ω1与Ω2,进行互相干处理,可表示为:
W=Ω1*Ω2 (7)
其中,*表示互相干运算,结果如图4所示。
利用短时傅里叶变换对W进行时-频分析,可表示为:
其中,p表示3-9Hz范围内的频率分量。
Ξ表示汉明窗函数,其可表示为
其中,α+β=1;本发明中α=0.42,β=0.58。
时-频分析结果如图5所示。
基于时-频分析结果K,重构回波信号可表示为:
重构回波信号如图6所示。
利用均值滤波器实现信号增强与噪声抑制,可表示为:
为了实现人体生命体征的提取,本实施例公开的变时窗分解方法如图8所示,将信号分为ψ组,第一组包含点数r,第二组点数r+o,第三组点数r+2o,…以此类推,第ψ组包括点数N。为了确保频率估计精度,要求r=2n,且r的采样时长至少为3-5倍的呼吸或心跳周期,n为正整数;o为正整数。本实施例中取r=512,o=100。
在对每组数据进行傅立叶变化之前,为了确保数据的频谱具有相同的频率分辨率,在每组数据后加0,直至每组数据具有N个值。
对变时窗得到的每组数据进行傅里叶变换,同时为了抑制带外噪声,对其进行加窗处理,可表示为:
Tj[n]=χ[n]{FFT{Sj[n]}}j=1,....,Ψ;n∈K*;K*=k*,k*+1,...,k*+κ-1 (12)
其中,FFT{Sj[n]}表示第j个时窗内的信号的傅里叶变换,k*表示加窗频率下限,本实施例中时窗范围为0.18-2.4Hz。
对于得到的具有相同频谱分辨率的ψ组频域特征做如下处理:
同时为了抑制带内噪声,对信号进行多重频率累积,可表示为:
X[n]=(Π[n]+jΠ[n])*(Π[n]+jΠ[n]) (14)
由此可得到生命体征微动信号的频率估计值,如图9所示。
Claims (8)
1.一种微动太赫兹谱识别检测算法,其特征在于,包括:
采用MIMO太赫兹雷达照射目标;
对得到的两路I/Q信号进行静态分量抑制;
对每一路I/Q信号进行反正切解调;
对解调后的信号进行互相干处理;
对互相干处理后的信号进行时-频分析;
基于时-频分析结果,重构回波信号;
对重构的回波信号进行时域滤波、变时窗分解;其中,变时窗分解包括:
将信号分为ψ组,第一组点数r,第二组点数r+o,第三组点数r+2o,…,以此类推,第ψ组包括点数N;其中:r=2n,且r的采样时长为3-5倍的呼吸或心跳周期,n为正整数;o为正整数;所有组的起点位置相同;
在每组数据后加0,直至每组数据具有N个值;
对变时窗得到的每组数据进行傅里叶变换,并对其进行加窗处理,表示为:
Τj[n]=χ[n]{FFT{Sj[n]}}j=1,....,Ψ;n∈K*;K*=k*,k*+1,...,k*+κ-1;
其中,FFT{Sj[n]}表示第j个时窗内的信号的傅里叶变换,k*表示加窗频率下限;
对于得到的具有相同频谱分辨率的ψ组频域特征做如下处理:
对信号进行多重频率累积,表示为:
X[n]=(Π[n]+jΠ[n])*(Π[n]+jΠ[n])
估计目标的太赫兹多普勒普。
4.根据权利要求3所述的微动太赫兹谱识别检测算法,其特征在于,互相干处理表示为:
W=Ω1*Ω2;
其中,*表示互相干运算。
8.根据权利要求1-7任一项所述的微动太赫兹谱识别检测算法,其特征在于,所述MIMO太赫兹雷达的发射信号功率范围为1mW-1W,带宽93.2-96.8GHz,信号模式可为有效带宽内的任意的点频连续波、线性调频连续波或步进频连续波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911291687.8A CN111142084B (zh) | 2019-12-11 | 2019-12-11 | 一种微动太赫兹谱识别检测算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911291687.8A CN111142084B (zh) | 2019-12-11 | 2019-12-11 | 一种微动太赫兹谱识别检测算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111142084A CN111142084A (zh) | 2020-05-12 |
CN111142084B true CN111142084B (zh) | 2023-04-07 |
Family
ID=70518406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911291687.8A Active CN111142084B (zh) | 2019-12-11 | 2019-12-11 | 一种微动太赫兹谱识别检测算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111142084B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111982856B (zh) * | 2020-08-27 | 2023-08-01 | 中电科思仪科技股份有限公司 | 一种基于太赫兹波的物质无标志检测识别方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976566A (zh) * | 2010-07-09 | 2011-02-16 | 瑞声声学科技(深圳)有限公司 | 语音增强方法及应用该方法的装置 |
CN103220440A (zh) * | 2012-01-19 | 2013-07-24 | 索尼公司 | 噪声抑制装置、噪声抑制方法以及程序 |
CN103674001A (zh) * | 2013-11-19 | 2014-03-26 | 南京航空航天大学 | 一种基于增强自适应时频峰值滤波的光纤陀螺去噪方法 |
CN105232026A (zh) * | 2015-10-29 | 2016-01-13 | 无锡南理工科技发展有限公司 | 一种非接触式生命体征检测***的心跳频率检测算法 |
CN108564965A (zh) * | 2018-04-09 | 2018-09-21 | 太原理工大学 | 一种抗噪语音识别*** |
CN109805931A (zh) * | 2019-02-18 | 2019-05-28 | 中电科仪器仪表有限公司 | 基于太赫兹多普勒雷达的远距离生命微动信号检测方法 |
CN110448290A (zh) * | 2019-08-19 | 2019-11-15 | 中电科仪器仪表有限公司 | 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及*** |
-
2019
- 2019-12-11 CN CN201911291687.8A patent/CN111142084B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976566A (zh) * | 2010-07-09 | 2011-02-16 | 瑞声声学科技(深圳)有限公司 | 语音增强方法及应用该方法的装置 |
CN103220440A (zh) * | 2012-01-19 | 2013-07-24 | 索尼公司 | 噪声抑制装置、噪声抑制方法以及程序 |
CN103674001A (zh) * | 2013-11-19 | 2014-03-26 | 南京航空航天大学 | 一种基于增强自适应时频峰值滤波的光纤陀螺去噪方法 |
CN105232026A (zh) * | 2015-10-29 | 2016-01-13 | 无锡南理工科技发展有限公司 | 一种非接触式生命体征检测***的心跳频率检测算法 |
CN108564965A (zh) * | 2018-04-09 | 2018-09-21 | 太原理工大学 | 一种抗噪语音识别*** |
CN109805931A (zh) * | 2019-02-18 | 2019-05-28 | 中电科仪器仪表有限公司 | 基于太赫兹多普勒雷达的远距离生命微动信号检测方法 |
CN110448290A (zh) * | 2019-08-19 | 2019-11-15 | 中电科仪器仪表有限公司 | 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及*** |
Non-Patent Citations (2)
Title |
---|
李媛 ; 铁勇 ; 那顺乌力吉 ; 李树华 ; 王强 ; .基于扩展谱减法的语音增强算法及实现.内蒙古大学学报(自然科学版).2008,(第01期),全文. * |
陈学华 ; 贺振华 ; 黄德济 ; .广义S变换及其时频滤波.信号处理.2008,(第01期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111142084A (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Vital sign detection method based on multiple higher order cumulant for ultrawideband radar | |
Liang et al. | An improved algorithm for through-wall target detection using ultra-wideband impulse radar | |
Narayanan | Through-wall radar imaging using UWB noise waveforms | |
Li et al. | Advanced signal processing for vital sign extraction with applications in UWB radar detection of trapped victims in complex environments | |
Ding et al. | Noncontact multiple targets vital sign detection based on VMD algorithm | |
Wu et al. | Improved human respiration detection method via ultra‐wideband radar in through‐wall or other similar conditions | |
Li et al. | Research on feature extraction of ship-radiated noise based on multi-scale reverse dispersion entropy | |
Lv et al. | Characterization and identification of IR-UWB respiratory-motion response of trapped victims | |
Li et al. | A new method for non-line-of-sight vital sign monitoring based on developed adaptive line enhancer using low centre frequency UWB radar | |
CN110161491B (zh) | 一种针对微弱生命体的测距和呼吸频率估计方法 | |
Liang et al. | Ultra-wide band impulse radar for life detection using wavelet packet decomposition | |
CN111142084B (zh) | 一种微动太赫兹谱识别检测算法 | |
Nejadgholi et al. | Time-frequency based contactless estimation of vital signs of human while walking using PMCW radar | |
Bi et al. | Adaptive target extraction method in sea clutter based on fractional fourier filtering | |
Choi et al. | Robust cardiac rate estimation of an individual | |
Yang et al. | Applications of TVF‐EMD in Vital Signal Detection for UWB Radar | |
Chen et al. | Ultrawideband synthetic aperture radar for respiratory motion detection | |
Ram et al. | Sparsity‐based autoencoders for denoising cluttered radar signatures | |
Li et al. | A new algorithm for through wall human respiration monioring using GPR | |
Li et al. | Macro-motion detection using ultra-wideband impulse radar | |
Xue et al. | A dynamic clutter interference suppression method for multiple static human targets detection using ultra‐wideband radar | |
Zhu et al. | Non-contact detection of vital signs with a hybrid feature extraction method using a uwb radar sensor | |
Lu et al. | Detection of human's motions through a wall using UWB radar | |
Xian et al. | Improved Denoising Method for UWB Vital Signs Detection and Extraction | |
Qian | Radar clutter suppression solution based on ICA |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |