CN103197313A - 子阵列时间反演镜探测方法 - Google Patents
子阵列时间反演镜探测方法 Download PDFInfo
- Publication number
- CN103197313A CN103197313A CN2013101464808A CN201310146480A CN103197313A CN 103197313 A CN103197313 A CN 103197313A CN 2013101464808 A CN2013101464808 A CN 2013101464808A CN 201310146480 A CN201310146480 A CN 201310146480A CN 103197313 A CN103197313 A CN 103197313A
- Authority
- CN
- China
- Prior art keywords
- omega
- signal
- time reversal
- subarray
- time
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种子阵列时间反演镜探测方法,属于时间反演技术领域,针对在探测多个目标时,传统时间反演镜方法无法得到正确的目标图像的问题,本发明的方法具体将时间反演探测阵列单元分成若干个子阵列,然后对每个子阵列的数据进行单独处理,通过寻找各子阵列的共同聚焦点来发现目标。与传统时间反演镜探测方法使用信号强度来标示目标位置不同,本发明的方法使用信号聚焦强度来寻找目标,可以实现多目标探测。
Description
技术领域
本发明属于时间反演技术领域,涉及目标探测领域,具体涉及一种多目标探测方法。
背景技术
首先介绍时间反演的数学描述过程。考虑两个传感器,分别定义为A和B,向传感器A发送激励信号v(t),其傅里叶变换为V(ω)。设通过传感器后,空间中所得到的激励信号的傅里叶变换为P(ω)。由于传感器自身具有非线性特性Ka(ω)(若传感器为天线,Ka(ω)则为天线的时域特性;若传感器为激励超声波信号的压电传感器,则Ka(ω)为压电传感器的电机械转化效率),P(ω)可以表示成:
P(ω)=Ka(ω)V(ω) (1-1)
设传感器B的非线性特性为Kb(ω),传感器B的时域输出电信号vB(t)可以表示为:
其中,Gab(ω)为传感器A到传感器B的传递函数。
出于方便的目的,定义Kab(ω)=Ka(ω)Kb(ω),在输出信号vB(t)中截取部分信号,设截取信号的时间长度为T,则公式(1-2)的时间反演可以表示为:
其中,“*”表示频域上的共轭,由此公式可以看出,时域上的时间反演等同于频域上的相位共轭。
在传感器B,重新激励上述时间反演信号,此时空间中的场与传感器A发送原始激励信号所形成的场,在时序上是相反的。传感器A所接收到的时间反演信号可以表示为
将传感器A接收到的时间反演信号再次进行时间反转,得到:
从公式(1-5)可以看到,传感器A接收到的信号所对应的时间反演版本是近似于原始发射信号的,所以,时间反演技术可以在时域和空域上恢复发射信号,即实现时空上的聚焦。
然后推导一般时间反演镜成像过程。出于技术简化,忽略阵列单元的时域非线性响应特性,设空间中有S个散射目标分别位于rs(1≤s≤S),探测阵列的信号接收单元共有N个,分别位于rn(1≤n≤N),散射目标和各信号接收单元之间的传递函数为G(rs,rn,ω),设目标信号的傅里叶变换为X(ω),则位于rn的第n个接收单元接收到的信号的傅里叶变换为:
以时间零点为起点取时间长度为T的一段上述信号,并且进行时间反演处理。根据公式(1-3),则该段信号的时间反演版本为:
其中,“*”表示共轭运算。
假设在成像阶段,在第n个接收点回传时间反演信号FTR(ω,rn),那么,在成像空间rk处的信号的傅里叶变换可以表示为:
如果所有时间反演信号被同时回传辐射,由公式(1-8)得,空间rk处的信号可以表示为:
对函数Y(ω,rk)进行反傅里叶变换,得到y(rk,t),
定义信号幅值函数c(rk,t):
于是,y(rk,t)可以重写为
考虑单个目标的情况,即S=1,那么位于成像空间rk处的总信号的傅里叶变换可以表示为:
对于rk=rs处,公式(1-14)可以表示为:
根据公式(1-12)、公式(1-15)和时间反演信号的聚焦效应,可以得到点rk=rs的时间反演信号的最大值为:
当N足够大的时候,目标点处的时间反演信号将成为整个成像空间的最大信号。因此,只要找出空间中最大信号的位置就能找到目标,这就是传统时间反演镜探测方法用来发现目标的标准。但是当N较小的时候,目标处的时间反演聚焦信号可能不再是空间最大信号。对此,可以进行如下解释。
由时间反演基本原理可知,所有时间反演信号不会聚焦在目标以外的空间点,所以在这些点处,各接收单元的时间反演信号不相互叠加。为了方便分析,考虑这样一些空间点,该点的最大信号将由某一接收单元的时间反演信号构成,例如成像空间中时间反演信号的回传辐射点。这里,假设该点位于它的最大时间反演信号为:
由于传统时间反演镜探测方法根据成像空间中的最大信号值来发现目标,因此要求但是由于自由空间损耗和媒质对信号的衰减,G(rn,rk,ω)随着rk的增加而迅速变小。所以对于较远的目标(即rs较大),如果不对信号进行补偿,根据公式(1-16)和公式(1-17),只有当N足够大时,即接收单元的数量足够多时,才能满足而且,所需的接收单元的数量会随着rs的增加而增加。因此当N较小的时候,目标处的时间反演聚焦信号可能不再是空间最大信号,此时若采用传统时间反演镜探测方法进行探测,将获得错误的目标图像。
接下来考虑两个目标,分别为位于的散射体s1和位于的散射体s2,这里假设,散射体s1与阵列距离较近,其散射信号较大;散射体s2与阵列距离较远,其散射信号较小,即:按传统的时间反演镜探测方法,将所有的接收信号进行时间反转处理后回传辐射。那么,根据公式(1-15)和公式(1-16),处和处的聚焦信号幅值分别为:
由于 因此,
发明内容
本发明的目的是为了解决在探测多个目标时,传统时间反演镜方法无法得到正确的目标图像的问题,提出一种子阵列时间反演镜探测方法。
本发明的技术方案为:一种子阵列时间反演镜探测方法,具体包括如下步骤:
S1.在成像阶段将探测阵列分为若干个子阵列;
S2.各子阵列分别回传各自的时间反演信号;
S3.记录成像区域中各子阵列在每一时刻的时间反演信号分布Ep(rk,t),并且进行归一化,得到E′p(rk,t);
S4.将各子阵列在每个像素点的归一化信号值按时间先后进行叠加,得到u(rk,t),
S5.将每个像素点的最大叠加信号值作为该点的像素值:I(rk)=umax(rk,t)。
本发明的有益效果是:本发明的方法是将时间反演探测阵列单元分成若干个子阵列,然后对每个子阵列的数据进行单独处理,通过寻找各子阵列的共同聚焦点来发现目标。与传统时间反演镜探测方法使用信号强度来标示目标位置不同,本发明的方法使用信号聚焦强度来寻找目标,可以实现多目标探测。
附图说明
图1是本发明方法的原理结构图。
图2是本发明方法的流程示意图。
附图标记说明:目标(rs,s=1,2,…,S)、信号发送单元(rm)、信号接收单元(rn,n=1,2,…,N)。
具体实施方式
下面结合附图对本发明的实施例做进一步的说明。
在设计新的时间反演多目标探测方法的时候,总希望利用时间反演时域聚焦特性。因此首先需要获取在多个不同目标处的时间反演信号所具有时域物理特性,接着根据这些物理特性才能构建新的探测方法。在推导这些物理性质之前,首先建立如附图1所示的探测模型。
设空间中有S个散射目标分别位于rs(1≤s≤S)。为了推导简单,这里设信号发送阵列只由一个信号发送单元构成,并且其位于rm;并且忽略阵列单元的时域非线性特性,该阵列信号发送单元和各散射目标之间的传递函数可以表示为G(rs,rm,ω),接收阵列的信号接收单元共有N个,分别位于rn(1≤n≤N),散射目标体和各信号接收单元之间的传递函数为G(rs,rn,ω)。
设探测信号为X(ω),则位于rn的第n个接收单元接收到的信号的傅里叶变换为:
以时间零点为起点取时间长度为T的一段上述信号,并且进行时间反演处理。
根据公式(1-3),则该段信号的时间反演公式为:
其中,“*”表示共轭运算。
假设在成像阶段,在第n个接收点发送反演信号YTR(ω,rn),那么位于探测空间的rk处信号的傅里叶变换可以表示为:
如果探测阵列中的Ni个探测点同时发送时间反演波,由公式(2-3)得,空间点rk处的信号可以表示为:
对函数F(ω,rk)进行反傅里叶变换,得到f(rk,t),
当rk=rs时,即在第s个散射体处,由正交性知,公式(2-4)可以表示为:
因此
式中,α(ω)是X(ω)的相位,β(rs,rm,ω)是G(rs,rm,ω)的相位。
定义信号幅值函数d(rk,t):
于是,f(rk,t)可以重写为:
所以,该Ni个接收点构成的阵列所发送的时间反演信号,在rk=rs的位置,获得最大值dmax(rk)的时间,即聚焦时间为:
再根据公式(2-8)和(2-11),可得
ts=T-β′(rs,rm,ω0)-α′(ω0) (2-12)
这里,假设在成像阶段,将整个目标信号接收阵列分成P个子阵列,每个子阵列由Np个接收单元构成,1≤p≤P。那么,根据公式(2-7)、(2-8)和公式(2-12),各子阵列在空间中S个目标散射体处的聚焦时间为:
可见,每个子阵列发送的时间反演信号在各目标点处获得最大值的时间与子阵列接收单元的位置以及接收单元的数量无关,所得的各子阵列在目标处的聚焦时间是一样的。
基于上述的时间反演信号在目标处的时域特性,提出一种子阵列时间反演镜探测方法,具体流程如图2所示,包括如下步骤:
(1)在成像阶段,将探测阵列分为若干个子阵列。(实验表明,阵列的排列方式对结果影响不大。)
(2)各子阵列分别回传各自的时间反演信号。(因为各探测点采用的是独立的传感器,所以不用考虑信号间的相互影响。)
(3)记录成像区域中各子阵列在每一时刻的时间反演信号分布Ep(rk,t),并且进行归一化,得到E′p(rk,t)。(这里的归一化,具体进行计算Ep(rk,t)/|Ep(rk,t)|。)
(4)将各子阵列在每个像素点的归一化信号值按时间先后进行叠加,得到u(rk,t)。
(5)将每个像素点的最大叠加信号值作为该点的像素值:
I(rk)=umax(rk,t) (2-15)
根据公式(2-13)可以知道各目标处的聚焦时间是不一样的,所以当按时间顺序对归一化信号进行叠加后,各目标点出的最大信号出现的时间是不同的,所以需要在整个时间范围内作求最大值操作,从而得到各像素点的最大叠加信号值。
由于时间反演技术的时空聚焦效应,各散射体的时间反演信号在各自的聚焦时刻是最大信号,因此在步骤(3)中,各目标处的聚焦信号将被归一化为1,在步骤(4)中,由于不同的子阵列的时间反演信号在同一目标处的聚焦时间相同,所以各目标点处的信号将会叠加,并且由于归一化的原因,所有目标处信号强度在数值上等于子阵列的个数。这里需要另外说明的是,归一化后,单一子阵列的最大信号值为1,并且信号叠加次数上限为子阵列的个数,所以在叠加后的成像区域中,各空间点的信号值上限在数值上也等于子阵列的个数,即叠加后的各目标点处的最大信号为成像区域中的最大信号。但是根据公式(2-13)可以知道各目标处的聚焦时间是不一样的,并且按时间先后对归一化信号进行叠加,这就意味着叠加后,各目标点处的最大信号出现的时间是不相同的,因而不能通过画出某一时刻的场分布得到所有目标的图像。因此在步骤(5)中,将叠加后的各像素点所对应的时域最大信号值作为该像素点的像素值,以便将所有目标处的最大信号值在同一成像区域中表示出来,由于经过归一化和叠加后各目标点处的最大信号相等且都为成像区域中的最大信号,所以在所得的最后图像中,所有的目标将以最大像素值的形式在成像区域中被标示出来,可见本发明的方法可以发现多个目标。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (1)
1.一种子阵列时间反演镜探测方法,具体包括如下步骤:
S1.在成像阶段将探测阵列分为若干个子阵列;
S2.各子阵列分别回传各自的时间反演信号;
S3.记录成像区域中各子阵列在每一时刻的时间反演信号分布Ep(rk,t),并且进行归一化,得到E′p(rk,t);
S4.将各子阵列在每个像素点的归一化信号值按时间先后进行叠加,得到u(rk,t),
S5.将每个像素点的最大叠加信号值作为该点的像素值:I(rk)=umax(rk,t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310146480.8A CN103197313B (zh) | 2013-04-25 | 2013-04-25 | 子阵列时间反演镜探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310146480.8A CN103197313B (zh) | 2013-04-25 | 2013-04-25 | 子阵列时间反演镜探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103197313A true CN103197313A (zh) | 2013-07-10 |
CN103197313B CN103197313B (zh) | 2014-11-05 |
Family
ID=48719989
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310146480.8A Expired - Fee Related CN103197313B (zh) | 2013-04-25 | 2013-04-25 | 子阵列时间反演镜探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103197313B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777519A (zh) * | 2015-04-21 | 2015-07-15 | 电子科技大学 | 一种电磁信号时间反演变换的方法及系统 |
CN108119764A (zh) * | 2017-12-26 | 2018-06-05 | 东莞理工学院 | 时间反演自适应网格输气管道检漏方法 |
CN108167657A (zh) * | 2017-12-26 | 2018-06-15 | 东莞理工学院 | 输气管道泄漏区域时间反演自动高分辨率定位方法 |
CN110456342A (zh) * | 2019-07-11 | 2019-11-15 | 西安电子科技大学 | 单发射天线雷达的远场多运动目标探测方法 |
CN112804633A (zh) * | 2020-12-28 | 2021-05-14 | 华南理工大学 | 一种时间反演高精度室内定位方法与系统 |
CN114895273A (zh) * | 2022-03-29 | 2022-08-12 | 中国人民解放军国防科技大学 | 基于时间反演的分布式侦干探一体化系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101873692A (zh) * | 2010-06-23 | 2010-10-27 | 电子科技大学 | 基于时间反演的无线传感器网络节点定位方法 |
US20120127020A1 (en) * | 2010-11-18 | 2012-05-24 | Eung Gi Paek | Distributed time-reversal mirror array |
-
2013
- 2013-04-25 CN CN201310146480.8A patent/CN103197313B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101873692A (zh) * | 2010-06-23 | 2010-10-27 | 电子科技大学 | 基于时间反演的无线传感器网络节点定位方法 |
US20120127020A1 (en) * | 2010-11-18 | 2012-05-24 | Eung Gi Paek | Distributed time-reversal mirror array |
Non-Patent Citations (3)
Title |
---|
ANTHONY J.DEVANEY ET AL: "Time-reversal-based imaging and inverse scattering of multiply scattering point targets", 《ACOUSTICAL SOCIETY OF AMERICA》, vol. 118, no. 5, 30 November 2005 (2005-11-30), XP012073430, DOI: doi:10.1121/1.2042987 * |
丁帅 等: "电磁信号时间反演的实现方法研究进展", 《电子科技大学学报》, vol. 40, no. 5, 30 September 2011 (2011-09-30) * |
张光旻: "新型时间反演探测方法和信息传输方法的研究", 《中国博士学位论文全文数据库 基础科学辑》, no. 05, 15 May 2013 (2013-05-15) * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104777519A (zh) * | 2015-04-21 | 2015-07-15 | 电子科技大学 | 一种电磁信号时间反演变换的方法及系统 |
CN104777519B (zh) * | 2015-04-21 | 2017-06-09 | 电子科技大学 | 一种电磁信号时间反演变换的方法及系统 |
CN108119764A (zh) * | 2017-12-26 | 2018-06-05 | 东莞理工学院 | 时间反演自适应网格输气管道检漏方法 |
CN108167657A (zh) * | 2017-12-26 | 2018-06-15 | 东莞理工学院 | 输气管道泄漏区域时间反演自动高分辨率定位方法 |
CN108119764B (zh) * | 2017-12-26 | 2020-02-07 | 东莞理工学院 | 时间反演自适应网格输气管道检漏方法 |
CN110456342A (zh) * | 2019-07-11 | 2019-11-15 | 西安电子科技大学 | 单发射天线雷达的远场多运动目标探测方法 |
CN112804633A (zh) * | 2020-12-28 | 2021-05-14 | 华南理工大学 | 一种时间反演高精度室内定位方法与系统 |
CN112804633B (zh) * | 2020-12-28 | 2021-11-02 | 华南理工大学 | 一种时间反演高精度室内定位方法与系统 |
CN114895273A (zh) * | 2022-03-29 | 2022-08-12 | 中国人民解放军国防科技大学 | 基于时间反演的分布式侦干探一体化系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103197313B (zh) | 2014-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103197313B (zh) | 子阵列时间反演镜探测方法 | |
Wen et al. | Joint 2D-DOD, 2D-DOA, and polarization angles estimation for bistatic EMVS-MIMO radar via PARAFAC analysis | |
Wen et al. | Fast direction finding for bistatic EMVS-MIMO radar without pairing | |
CN103760547B (zh) | 基于互相关矩阵的双基mimo雷达角度估计方法 | |
CN104749552A (zh) | 基于稀疏重构的互质阵列波达方向角估计方法 | |
CN107957574B (zh) | 基于ifft和混合匹配追踪的时分地基mimo滑坡雷达成像方法 | |
CN103650366B (zh) | 波束形成方法和设备 | |
Capistrano–Filho et al. | Internal controllability of the Korteweg–de Vries equation on a bounded domain | |
CN104535998B (zh) | 旋转镜像综合孔径辐射计及测量方法 | |
CN104635221B (zh) | 一种基于内定标数据的子带拼接方法 | |
CN102809746B (zh) | 一种基于mimo技术的高分辨扇扫成像方法 | |
CN104155648A (zh) | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 | |
CN103245943A (zh) | 一种使用线性调频信号的mimo阵列扇扫成像方法 | |
CN103605112B (zh) | 多发多收干涉合成孔径雷达时频二维信号波形设计方法 | |
CN103605113A (zh) | 多发多收干涉合成孔径雷达空时二维信号波形设计方法 | |
Li et al. | Reduced-dimensional ESPRIT for direction finding in monostatic MIMO radar with double parallel uniform linear arrays | |
WO2016005819A3 (en) | Multi-cell transducer | |
CN103163527B (zh) | 基于ccd辅助定位-双光楔光束调节的激光监听装置及采用该装置实现监听的方法 | |
CN106054121A (zh) | 一种确定雷达目标角度的方法及装置 | |
De Marchi et al. | Spiral-shaped piezoelectric sensors for Lamb waves direction of arrival (DoA) estimation | |
CN104034801A (zh) | 基于合成时反的结构损伤迭代聚焦成像监测方法 | |
CN102394686B (zh) | 高精度阵列天线接收系统角度估计的装置及其方法 | |
CN103217671A (zh) | 色噪声环境下的多输入多输出雷达收发角度快速估计方法 | |
CN103308877A (zh) | 基于分离式偶极子对阵列的多参数估计方法 | |
Mao et al. | Transmit subaperturing for MIMO radars with nested arrays |
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 |
Granted publication date: 20141105 Termination date: 20150425 |
|
EXPY | Termination of patent right or utility model |