CN107942322B - 多径利用雷达时间反演目标距离估计方法 - Google Patents

多径利用雷达时间反演目标距离估计方法 Download PDF

Info

Publication number
CN107942322B
CN107942322B CN201711068989.XA CN201711068989A CN107942322B CN 107942322 B CN107942322 B CN 107942322B CN 201711068989 A CN201711068989 A CN 201711068989A CN 107942322 B CN107942322 B CN 107942322B
Authority
CN
China
Prior art keywords
signal
time
target
multipath
echo
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
Application number
CN201711068989.XA
Other languages
English (en)
Other versions
CN107942322A (zh
Inventor
张娟
王晨红
郭玉梅
张林让
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201711068989.XA priority Critical patent/CN107942322B/zh
Publication of CN107942322A publication Critical patent/CN107942322A/zh
Application granted granted Critical
Publication of CN107942322B publication Critical patent/CN107942322B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/08Systems for measuring distance only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

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

本发明公开了一种多径利用雷达时间反演目标距离估计方法,解决了多径环境下雷达目标准确测距的问题。实现步骤是:产生雷达发射信号并得到回波数据;将回波数据进行能量归一化并时间反演,得到时反信号;将时反信号重新发射到环境中,得到时反回波数据;构造匹配滤波器,对时反回波进行匹配滤波,选择幅值门限进行峰值检测,峰值中首径时延对应的距离即为目标距离。本发明将回波数据能量归一化时间反演并重新发射到环境中,通过将多径能量和直达波能量在时间上聚焦,利用多径能量增加了时反回波的信多噪比,实现了对目标的距离估计,且在低信噪比情况下测距精度更高,成本较低,适用范围较广,可用于海面低空雷达点目标距离估计。

Description

多径利用雷达时间反演目标距离估计方法
技术领域
本发明属于雷达技术领域,特别涉及目标距离估计,具体是一种多径利用雷达时间反演目标距离估计方法,可用于对海面低空点目标的测距。
背景技术
雷达在海面低空环境测量目标距离时会受到多径效应及噪声的影响,即雷达收到的回波信号包括来自目标的直达波信号、来自海面反射的多径波信号及噪声信号。受到多径效应和噪声的影响,参照说明书附图6,当信噪比处于[-11dB,-13dB]范围内,若使用传统的匹配滤波方法测距,距离估计均方根误差极大,无法测出目标距离。
针对多径效应对雷达在海面低空环境目标距离估计的影响,现有方法主要通过对多径信号在空域和时域进行抑制来提高测距精度。在空域方面,设置天线位置或设计地平面天线,使多径信号无法进入天线,但无法适用于所有情况,实现成本较高。在时域方面,使用窄相关技术或前边缘技术,但会带来距离估计误差。
受到多径效应的影响,在低信噪比情况下,现有技术中的匹配滤波方法距离估计精度较低,为了提高精度,通过空域和时域的方法对多径信号进行抑制,存在着成本较高、适用范围较低等问题。
发明内容
本发明的目的是针对上述方法的不足,提出一种更加准确的多径利用雷达时间反演目标距离估计方法,利用多径信息对目标进行距离估计。
本发明是一种多径利用雷达时间反演目标距离估计方法,其特征在于,利用多径信息对目标进行测距,包括有如下步骤:
(1)获取雷达回波信号:雷达产生发射信号st(t)并发射到环境中,得到雷达回波信号sr(t);
(2)获取雷达时反信号:将雷达回波数据sr(t)进行能量归一化并进行时间反演处理,得到雷达时反信号
Figure GDA0002802811680000011
(3)获取时反回波信号:将时反信号
Figure GDA0002802811680000012
重新发射到环境中,得到时反回波信号为
Figure GDA0002802811680000021
(4)对目标进行距离估计:首先构造时反匹配滤波器,用该滤波器对时反回波信号
Figure GDA0002802811680000022
进行匹配滤波处理,实现目标距离估计。
本发明通过将回波进行能量归一化并时间反演重发,对再次接收到的时反回波进行匹配滤波,提高了时反回波的信多噪比,减小了目标距离估计误差。
本发明与现有技术相比,具有以下优点:
一、在多径环境下,受到噪声的影响,若信噪比过低,来自目标的回波将被噪声信号淹没,传统的匹配滤波方法无法准确估计目标的距离。本发明将回波信号进行时间反演并重新发送到环境中,通过时间反演使时反信号对信道产生自适应聚焦效果,实现了对多径信号的利用,增强了时反回波的信多噪比。
二、本发明通过对回波信号进行时间反演,利用多径信号的能量增加了时反回波信号的信多噪比,再对时反回波信号进行匹配滤波,可以得到更高的测距精度及更小的距离估计均方根误差。
三、为减弱多径效应对测距的影响,传统方法采取抑制多径效应的方法。本发明将回波进行时反重发,有效利用多径信息进行距离估计,无需考虑对多径信号的抑制。本发明方法成本较低,且适用范围较广。
附图说明
图1是本发明的实现流程图;
图2是在设定的目标参数和环境参数下,传统匹配滤波测距结果图;
图3是在设定的目标参数和环境参数下,本发明时反匹配滤波测距结果图。
图4是在设定的目标参数和环境参数下,本发明时反回波信多噪比增益和信噪比的关系曲线图。
图5是在设定的目标参数和环境参数下,本发明时反回波信多噪比增益和散射系数模值的关系曲线图。
图6是在设定的目标参数和环境参数下,传统匹配滤波距离估计均方根误差和信噪比的关系曲线图。
图7是在设定的目标参数和环境参数下,本发明时反匹配滤波距离估计均方根误差和信噪比的关系曲线图。
具体实施方式
下面结合附图,对本发明详细说明:
实施例1
雷达在海面低空环境中进行目标测距时,传统方法对回波信号进行匹配滤波,但在低信噪比情况下,受到噪声影响,无法准确估计目标距离。针对此现状,本发明展开了研究,针对性地提出了一种成本较低,适用于各种延迟多径,能改善回波信多噪比且精度更高且的海面低空点目标距离估计方法。
本发明是一种多径利用雷达时间反演目标距离估计方法,参照图1,利用多径信息对目标进行测距,包括有如下步骤:
(1)雷达产生发射信号st(t),得到雷达回波数据sr(t)。
(2)将雷达回波数据sr(t)进行能量归一化并进行时间反演处理,得到时反信号
Figure GDA0002802811680000031
受到多径效应和接收噪声等因素的影响,对雷达回波数据sr(t)进行能量归一化,即可保证雷达两次发射的信号功率相等。将能量归一化后的数据进行时间反演,即在时间上将信号数据反序,以便于实现信号对信道的自适应。
(3)将时反信号
Figure GDA0002802811680000032
重新发射到环境中,得到时反回波为
Figure GDA0002802811680000033
将时反信号
Figure GDA0002802811680000034
重新发射到环境中,即可实现信号对信道的自适应效果,且时反回波中的分量将在st *(-t)处实现能量聚焦,时反回波的信多噪比将实现提高。
(4)构造时反匹配滤波器并对时反回波
Figure GDA0002802811680000035
进行匹配滤波处理,实现距离估计。
通过构造时反匹配滤波器对时反回波进行匹配滤波,即可提取目标直达波时延,也即实现了目标距离的估计。同样由于时间反演给时反回波带来了能量聚焦效果,提高了信多噪比,从而距离估计精度也将得到提升。
在海面低空环境下,受到多径信号的影响,若信噪比过低,传统的匹配滤波方法将无法准确估计目标的真实距离。本发明将回波信号能量归一化并时间反演得到时反回波信号,对时反回波进行匹配滤波,设置幅值门限对时反回波进行峰值检测,通过提取峰值中的首径时延可以准确地估计出目标的真实距离。由于对回波信号进行了时间反演重发,时反回波将实现能量聚焦效果,且信多噪比将增大,从而提高了距离估计精度。
实施例2
多径利用雷达时间反演目标距离估计方法同实施例1,本发明步骤(2)中所述将雷达回波数据sr(t)进行能量归一化并进行时间反演处理,得到时反信号
Figure GDA0002802811680000036
具体包括:
2a)根据发射信号st(t)和回波数据sr(t)计算能量归一化因子k:
Figure GDA0002802811680000041
在时间上积分计算发射信号st(t)和回波数据sr(t)的能量比值,即为能量归一化因子,其中T为发射信号时宽,*为对信号取共轭,η1为一阶多径波信道复电磁散射系数,η2为二阶多径波信道复电磁散射系数,
Figure GDA0002802811680000042
为接收机噪声信号方差。
2b)将回波sr(t)进行能量归一化并时间反演:
Figure GDA0002802811680000043
其中*为共轭操作,τ为直达波信道延时,Δτ为一阶多径波信道和直达波信道的延时差,n(t)为接收机噪声信号,其方差为
Figure GDA0002802811680000044
将回波数据sr(t)乘以归一化因子k,即可保证再次发射的信号和发射信号st(t)能量相等。将回波数据进行时间反演,即对信号进行时间上的反序并取共轭,将t变为-t,得到时反信号,且时反信号将对环境信道产生自适应效果。
在多径环境下,若噪声能量较大,信噪比过低,传统的匹配滤波方法将无法准确估计目标的距离。本发明对回波信号进行时间反演并重新发送到环境中,使时反信号具备了对信道的自适应能力,导致接收到的时反回波的信多噪比得到增强。
实施例3
多径利用雷达时间反演目标距离估计方法同实施例1-2,步骤(4)中的对目标进行距离估计,具体包括有如下步骤:
4a)构造时反匹配滤波器并进行匹配滤波处理:构造时反匹配滤波器冲激响应函数为hMF(t),用该滤波器对时反回波信号
Figure GDA0002802811680000045
进行匹配滤波处理,对时反回波信号
Figure GDA0002802811680000046
进行匹配滤波,得到时反匹配滤波输出信号sMF
4b)目标距离估计:选择幅值门限Γm,对时反匹配滤波输出信号进行峰值检测,得到检测结果信号s′MF,提取检测结果信号s′MF中的首径时延τ1,其对应的距离即为目标的距离。
实施例4
多径利用雷达时间反演目标距离估计方法同实施例1-3,步骤4a)中的构造时反匹配滤波器并进行匹配滤波处理,具体包括有如下步骤:
4a1)构造时反匹配滤波器,滤波器的冲激响应函数为hMF(t)=ksr(t);
用回波信号sr(t)构造时反匹配滤波器,其中包含了目标的直达波距离信息,便于进行匹配滤波后估计出目标的距离。
4a2)用该滤波器对时反回波
Figure GDA0002802811680000051
进行匹配滤波,得到时反匹配滤波结果信号sMF
Figure GDA0002802811680000052
其中
Figure GDA0002802811680000053
为卷积符号,k为能量归一化因子,sinc(·)为辛克函数。通过对时反回波
Figure GDA0002802811680000054
进行匹配滤波,sMF中包括了多个峰值函数sinc(·),估计其中的sinc[πB(t-τ)]峰值项对应延时即可估计出目标距离。
实施例5
多径利用雷达时间反演目标距离估计方法同实施例1-4,步骤4b)中的对目标进行距离估计,具体包括有如下步骤:
4b1)选择幅值门限Γm,对时反匹配滤波结果信号进行峰值检测,得到检测结果信号s′MF
Figure GDA0002802811680000055
测出时反匹配滤波输出信号sMF中幅值最大值为Vmax,门限Γm取值范围为:[0.3Vmax,0.5Vmax],本例中门限Γm=0.4Vmax,sMF中的幅值小于门限的部分将被置零,大于等于门限的部分将保留不变,从而实现峰值项sinc[πB(t-τ+Δτ)]和sinc[πB(t-τ+2Δτ)]被滤除。
4b2)提取检测结果信号s′MF中的首径时延τ1,其对应的距离即为目标的距离。
选择门限对sMF进行峰值检测,选择s′MF中的首径信号即第一个峰值对应的信号sinc[πB(t-τ)],估计首径信号中的时延τ,即可估计出雷达到目标距离为R=1.5×108×τ,单位为米。
本发明通过对回波信号进行时间反演并重新发射到环境中,利用多径信号的能量增强了时反回波信号的信多噪比,在低信噪比环境下实现了更高的测距精度及更小的距离估计均方根误差。
下边给出一个更为详尽的例子,对本发明进一步说明。
实施例6
多径利用雷达时间反演目标距离估计方法同实施例1-5,参照图1,本发明的具体实现步骤如下:
步骤1:产生雷达发射信号st(t),信号形式为线性调频信号,其带宽为B,时宽为T,得到雷达回波数据sr(t):
Figure GDA0002802811680000061
其中h(t)为环境总信道传输响应函数,n(t)为接收机噪声信号,看作高斯白噪声信号,方差为
Figure GDA0002802811680000062
η1为一阶多径波信道复电磁散射系数,η2为二阶多径波信道复电磁散射系数,且η2=η1η1
步骤2:将雷达回波数据sr(t)进行能量归一化并进行时间反演处理,得到时反信号
Figure GDA0002802811680000063
根据发射信号st(t)和回波数据sr(t)计算能量归一化因子k:
Figure GDA0002802811680000064
其中T为发射信号时宽,*为共轭操作,η1为一阶多径波信道复电磁散射系数,η2为二阶多径波信道复电磁散射系数,
Figure GDA0002802811680000065
为接收机噪声信号方差。
将回波sr(t)进行能量归一化并时间反演得到时反信号
Figure GDA0002802811680000066
Figure GDA0002802811680000071
本发明对回波sr(t)进行能量归一化并时间反演,归一化保证了两次发射信号的能量相等,且时间反演可以使信号在重发过程中通过利用多径实现能量聚焦。
步骤3:将时反信号
Figure GDA0002802811680000072
重新发射到环境中,得到时反回波为
Figure GDA0002802811680000073
Figure GDA0002802811680000074
其中n(t)为接收机噪声信号,看作高斯白噪声信号,方差为
Figure GDA0002802811680000075
由于对回波信号进行了时间反演,再次接收到的时反回波的信多噪比将得到提高,信多噪比增益SMNRgain为:
Figure GDA0002802811680000076
其中|η1|为一阶多径波信道复电磁散射系数模值,|η2|为二阶多径波信道复电磁散射系数模值,n(t)为接收机噪声信号,看作高斯白噪声信号,方差为
Figure GDA0002802811680000077
k为能量归一化因子。在海面低空环境中,参照图5,图5给出了信噪比取值[-20dB,-7dB,0dB,10dB]情况下信多噪比增益和散射系数模值的关系,散射系数模值取值范围为[0,1],从图中可以看出,在同一信噪比情况下,信多噪比增益随着散射系数模值变大而增加,在同一散射系数模值情况下,信噪比越低,信多噪比增益越高。在低空海面环境下,海面散射系数模值|η1|>0.7,信多噪比增益SMNRgain>0,本发明时间反演回波的信多噪比得到提升,距离估计均方根误差变小,测距精度变高。
步骤4:构造时反匹配滤波器并对时反回波
Figure GDA0002802811680000081
进行匹配滤波处理,实现距离估计。
4a)构造时反匹配滤波器冲激响应函数为hMF(t)=ksr(t)。用该滤波器对时反回波
Figure GDA0002802811680000082
进行匹配滤波,忽略其中的噪声项,得到时反匹配滤波结果sMF
Figure GDA0002802811680000083
从sMF中可以看出,峰值项sinc[πB(t-τ)]对应的延时τ即为所估计的目标距离,滤除sMF中的sinc[πB(t-τ+2Δτ)]和sinc[πB(t-τ+Δτ)]项并取首径信号,即可估计出目标距离。
4b)选择幅值门限Γm,对时反匹配滤波结果sMF进行峰值检测,得到检测结果信号s′MF
Figure GDA0002802811680000084
对sMF过门限峰值检测,本例中门限取值Γm=0.3Vmax,则sMF中小于门限Γm的部分将被滤除,再提取s′MF中的首径时延τ1,即第一个峰值对应的时延,τ1对应的距离即为所估计的目标的距离。
在多径环境下,若信噪比过低,传统的匹配滤波测距方法无法准确测出距离,产生较大的测距误差,而本发明通过对回波信号进行能量归一化并时间反演重新发射到环境中,时反信号可以对环境信道进行自适应聚焦,使时反回波的信多噪比得到提高,在低信噪比情况下估计出目标距离信息,且估计精度更高。为减弱多径效应对测距的影响,传统方法采取抑制多径效应的方法,而本发明则是通过利用多径信息实现时反信号的能量聚焦效果,并对时反回波进行匹配滤波,实现距离估计。本发明无需对多径信号进行抑制,成本较低,适用范围更广、距离估计误差更小。
下面通过仿真实验和数据对本发明的技术效果再做说明。
实施例7
多径利用雷达时间反演目标距离估计方法同实施例1-6,本例对多径环境下的目标距离估计方法通过以下仿真进行验证。
(1)实验场景:
在相同的环境参数、雷达参数情况下,分别采取传统的匹配滤波测距方法和本发明多径利用雷达时间反演目标距离估计方法进行仿真实验。
(2)实验内容和实验结果分析:
实验条件:在计算机中模拟对海面低空点目标的测距,设置环境参数、雷达参数如表1。
表1.海面低空点目标测距参数表
Figure GDA0002802811680000091
在上述参数情况下,考虑一阶和二阶多径对单目标距离估计的影响,发射信号为线性调频信号,采用传统的匹配滤波方法进行距离估计,得到结果图如图2所示,其中横坐标为距离值,纵坐标为归一化匹配滤波幅值。从图2中可以看出,受到多径效应的影响,匹配滤波测距结果图中出现3个尖峰,距离为30000m处对应的峰值为直达波,距离为30020m处对应的峰值为一阶多径波,距离为30040m处对应的峰值为二阶多径波,直达波距离为目标真实距离,一阶多径波和二阶多径波视为由于多径效应而产生的干扰信号,另外在图中所有距离单元都分布着随机噪声带来的信号毛刺。
在同样仿真条件下,采用本发明多径利用雷达时间反演目标距离估计方法进行距离估计,得到结果图如图3所示,其中横坐标为距离值,纵坐标归一化匹配滤波幅值。本发明通过对回波信号进行时间反演,时反回波的信多噪比得到提高,从而测距精度得到提升。从图3中可以看出,提取匹配滤波结果信号中的首径时延,其对应的距离为30000m,与参数表中目标预设距离相同,实现了对目标距离的准确测量。对比图2和图3,本发明获得的图3中测距结果峰值主瓣更窄,随着信噪比的降低,本发明方法距离估计误差更小。
实施例8
多径利用雷达时间反演目标距离估计方法同实施例1-7,本发明对多径环境下的目标距离估计性能通过以下仿真进一步验证。
(1)实验场景:
在相同的环境参数、雷达参数情况下,分别采取传统的匹配滤波测距方法和本发明多径利用雷达时间反演目标距离估计方法进行仿真实验。
(2)实验内容和实验结果分析:
信噪比变化范围为[-20dB,20dB],散射系数模值变化范围为[0,1],采用本发明多径利用雷达时间反演目标距离估计方法进行距离估计,时反回波信多噪比增益和信噪比的关系如图4所示,其中横坐标为信噪比,纵坐标为信多噪比增益。
从图4中可以看出,在同一散射系数模值情况下,随着信噪比下降,信多噪比增益变大,信多噪比增益和信噪比的关系比较稳定。如散射系数|η1|=0.9时,信噪比从20dB变化到-20dB,信多噪比增益随之变大。在同一信噪比情况下,散射系数模值变大,信多噪比增益也随之变大。如信噪比为-10dB时,散射系数从|η1|=0.1变化到|η1|=1,信多噪比增益随之变大。
时反回波信多噪比增益和散射系数模值的关系如图5所示,其中横坐标为散射系数模值,纵坐标为信多噪比增益。从图5中可以看出,在同一信噪比情况下,若散射系数模值越大,时反回波的信多噪比增益越大,在同一散射系数情况下,信噪比越小,时反回波的信多噪比增益越大。
结合图4和图5,对于不同的信噪比,时反回波信多噪比增益和散射系数的变化规律较为稳定,对于不同的散射系数,时反回波信多噪比增益和信噪比的变化规律同样较为稳定。在海面低空环境下,如信噪比取值较小不超过0dB,散射系数模值|η1|>0.7时,可以保证信多噪比增益SMNRgain>0,也证明了本发明在不同信噪比和散射系数情况下的信多噪比增益的稳定性和本发明距离估计方法的可靠性。
实施例9
多径利用雷达时间反演目标距离估计方法同实施例1-8,本发明对多径环境下的目标距离估计性能通过以下仿真进一步验证。
(1)实验场景:
在相同的环境参数、雷达参数情况下,分别采取传统的匹配滤波测距方法和本发明多径利用雷达时间反演目标距离估计方法进行仿真实验,设置环境参数、雷达参数如表2,选择不同的信噪比和散射系数模值,信噪比范围选择为[-13dB,-6.5dB],以0.5dB为间隔,散射系数模值分别选择为[0.75,0.8,0.85],进行300次蒙特卡洛仿真。
表2. 300次蒙特卡洛仿真测距参数表
Figure GDA0002802811680000111
采用传统的匹配滤波方法进行300此蒙特卡洛仿真,计算不同信噪比和散射系数模值情况下距离估计结果和理论值的均方根误差,得到结果图如图6所示,其中横坐标为信噪比,纵坐标为距离估计均方根误差。图6中小图为信噪比从-6.5dB到-9.5dB时距离估计均方根误差和信噪比的关系曲线放大图。从图中可以看出,当信噪比高于-9.5dB,传统匹配滤波距离估计均方根误差较小,其范围处于1m和2m之间。当信噪比低于-9.5dB,传统匹配滤波距离估计均方根误差较大,距离估计均方根误差为数十米,数百米量级。
在同样仿真条件下,采用本发明多径利用雷达时间反演目标距离估计方法进行距离估计,得到结果图如图7所示,其中横坐标为距离值,纵坐标为距离估计均方根误差。图7中小图为信噪比从-6.5dB到-12.5dB时距离估计均方根误差和信噪比的关系曲线放大图。从图7中可以看出,当信噪比高于-12.5dB,时反匹配滤波距离估计误差较小,其范围处于0和4m之间,当信噪比低于-12.5dB,时反匹配滤波距离估计均方根误差较大,量级为数百米。
对比图6和图7,当信噪比在-6.5dB和-9.5dB之间,传统匹配滤波距离估计均方根误差范围在1m和2m之间,本发明时间反演匹配滤波距离估计均方根误差范围在0和1m之间。当信噪比在-10dB之间-13dB时,传统匹配滤波距离估计均方根误差较大,其量级在数十米到数百米之间,无法准确估计目标距离,而本发明时间反演匹配滤波仍然可以保证较低的测距均方根误差,处于在0到4m之间。对比不同情况下的距离估计均方根误差,本发明的距离估计误差更小,证明了本发明可以更准确地实现目标的距离估计。
简而言之,本发明公开的多径利用雷达时间反演目标距离估计方法,解决了多径环境下雷达目标准确测距的问题。实现步骤是:产生雷达发射信号并得到回波数据;将回波数据进行能量归一化并时间反演,得到时反信号;将时反信号重新发射到环境中,得到时反回波数据;构造匹配滤波器,对时反回波进行匹配滤波,选择幅值门限进行峰值检测,峰值中首径时延对应的距离即为目标距离。本发明将回波数据能量归一化时间反演并重新发射到环境中,通过将多径能量和直达波能量在时间上聚焦,利用多径能量增加了时反回波的信多噪比,实现了对目标的距离估计,且在低信噪比情况下测距精度更高,成本较低,适用范围较广,可用于海面低空雷达点目标距离估计。
以上描述仅是本发明的一个具体实例,并未构成对本发明的任何限制,显然对于本领域的专业人员来说,在了解了本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修改和改变,但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。

Claims (2)

1.一种多径利用雷达时间反演目标距离估计方法,其特征在于,利用多径信息对目标进行测距,包括有如下步骤:
(1)获取雷达回波信号:雷达产生发射信号st(t)并发射到环境中,得到雷达回波信号sr(t);
(2)获取雷达时反信号:将雷达回波数据sr(t)进行能量归一化并进行时间反演处理,得到雷达时反信号
Figure FDA0002802811670000011
(3)获取时反回波信号:将时反信号
Figure FDA0002802811670000012
重新发射到环境中,得到时反回波信号为
Figure FDA0002802811670000013
(4)对目标进行距离估计:首先构造时反匹配滤波器,用该滤波器对时反回波信号
Figure FDA0002802811670000014
进行匹配滤波处理,实现目标距离估计,对目标进行距离估计,具体包括有如下步骤:
4a)构造时反匹配滤波器并进行匹配滤波处理:构造时反匹配滤波器冲激响应函数为hMF(t),用该滤波器对时反回波信号
Figure FDA0002802811670000015
进行匹配滤波处理,对时反回波信号
Figure FDA0002802811670000016
进行匹配滤波,得到时反匹配滤波输出信号sMF;构造时反匹配滤波器并进行匹配滤波处理,具体包括有如下步骤:
4a1)构造时反匹配滤波器,其冲激响应函数为hMF(t)=ksr(t);
4a2)用该滤波器对时反回波
Figure FDA0002802811670000017
进行匹配滤波,得到时反匹配滤波结果信号sMF
Figure FDA0002802811670000021
其中,B为带宽,
Figure FDA0002802811670000022
为卷积符号,*为对信号取共轭,k为能量归一化因子,τ为直达波信道延时,Δτ为一阶多径波信道和直达波信道的延时差,η1为一阶多径波信道复电磁散射系数,η2为二阶多径波信道复电磁散射系数,sinc(·)为辛克函数;
4b)目标距离估计:选择幅值门限Γm,对时反匹配滤波输出信号进行峰值检测,得到检测结果信号s′MF,提取检测结果信号s′MF中的首径时延τ1,其对应的距离即为目标的距离;对目标进行距离估计,具体包括有如下步骤:
4b1)选择幅值门限Γm,对时反匹配滤波结果信号进行峰值检测,得到检测结果信号s′MF
Figure FDA0002802811670000023
根据sMF中的sinc[πB(t-τ-Δτ)]和sinc[πB(t-τ+Δτ)]峰值项幅值比值计算门限Γm,以使sMF中sinc[πB(t-τ+2Δτ)]和sinc[πB(t-τ+Δτ)]峰值项被滤除,从而得到s′MF,进行距离估计,具体门限Γm取值范围为:[0.3Vmax,0.5Vmax],其中Vmax为sMF中峰值幅值的最大值;
4b2)提取检测结果信号s′MF中的首径时延τ1,其对应的距离即为目标的距离。
2.根据权利要求1所述的多径利用雷达时间反演目标距离估计方法,其特征在于,步骤(2)中所述的将雷达回波数据sr(t)进行能量归一化并进行时间反演处理,得到时反信号
Figure FDA0002802811670000031
具体包括:
2a)根据发射信号st(t)和回波数据sr(t)计算能量归一化因子k:
Figure FDA0002802811670000032
其中T为发射信号时宽,*为对信号取共轭,η1为一阶多径波信道复电磁散射系数,η2为二阶多径波信道复电磁散射系数,
Figure FDA0002802811670000033
为接收机噪声信号方差;
2b)将回波sr(t)进行能量归一化并时间反演:
Figure FDA0002802811670000034
其中*为共轭操作,τ为直达波信道延时,Δτ为一阶多径波信道和直达波信道的延时差,η1为一阶多径波信道复电磁散射系数,η2为二阶多径波信道复电磁散射系数,n(t)为接收机噪声信号。
CN201711068989.XA 2017-11-03 2017-11-03 多径利用雷达时间反演目标距离估计方法 Active CN107942322B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711068989.XA CN107942322B (zh) 2017-11-03 2017-11-03 多径利用雷达时间反演目标距离估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711068989.XA CN107942322B (zh) 2017-11-03 2017-11-03 多径利用雷达时间反演目标距离估计方法

Publications (2)

Publication Number Publication Date
CN107942322A CN107942322A (zh) 2018-04-20
CN107942322B true CN107942322B (zh) 2021-03-23

Family

ID=61933281

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711068989.XA Active CN107942322B (zh) 2017-11-03 2017-11-03 多径利用雷达时间反演目标距离估计方法

Country Status (1)

Country Link
CN (1) CN107942322B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031231B (zh) * 2018-08-03 2023-02-10 西安电子科技大学 雷达低空目标时间反演相干角度估计方法
CN109782228B (zh) * 2019-02-21 2022-08-23 哈尔滨工程大学 瞬态信号的多定位节点联合波形估计方法
CN110109092B (zh) * 2019-04-24 2023-03-21 西安电子科技大学 多径环境下基于时间反演的雷达测速方法
CN110109077B (zh) * 2019-04-26 2023-03-28 西安电子科技大学 基于时间反演的mimo雷达相干角度估计方法
CN110515052B (zh) * 2019-06-18 2021-07-16 西安电子科技大学 一种基于时间反演的超宽带频域非等间隔采样目标检测方法
CN110596686A (zh) * 2019-09-04 2019-12-20 电子科技大学 一种基于时间反转技术的频控阵时间聚焦方法
CN111273294A (zh) * 2020-02-19 2020-06-12 西北工业大学 一种基于时间反转的水声脉冲回波测距方法
CN113009448B (zh) * 2021-03-09 2022-12-06 森思泰克河北科技有限公司 多径目标的检测方法、装置、设备和存储介质
CN113447983B (zh) * 2021-06-23 2022-10-04 湖南国天电子科技有限公司 一种浅地层剖面仪的数据采集及信号处理方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102005828A (zh) * 2010-11-30 2011-04-06 电子科技大学 基于时间反演的无线传感器网络节点无线充电系统及方法
CN104717173A (zh) * 2015-03-16 2015-06-17 电子科技大学 基于信道去耦的子载波复数均衡tmo无线通信方法
CN106226754A (zh) * 2016-08-22 2016-12-14 西安电子科技大学 基于时间反演的低仰角波达方向估计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102005828A (zh) * 2010-11-30 2011-04-06 电子科技大学 基于时间反演的无线传感器网络节点无线充电系统及方法
CN104717173A (zh) * 2015-03-16 2015-06-17 电子科技大学 基于信道去耦的子载波复数均衡tmo无线通信方法
CN106226754A (zh) * 2016-08-22 2016-12-14 西安电子科技大学 基于时间反演的低仰角波达方向估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Time-Reversal-Based Multipath Mitigation Technique for ISAR Images;Imanol de Arriba-Ruiz等;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20130531;第51卷(第5期);第3119~3138页 *
基于远场时间反演的亚波长天线阵列研究;梁木生 等;《物理学报》;20130315;第62卷(第5期);第058401-1~058401-7页 *

Also Published As

Publication number Publication date
CN107942322A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN107942322B (zh) 多径利用雷达时间反演目标距离估计方法
CN105785324B (zh) 基于mgcstft的线性调频信号参数估计方法
CN107607937B (zh) 基于时间反演的雷达目标测距方法
CN104142496B (zh) 基于连通域划分的统计mimo雷达多目标定位方法
CN104111449B (zh) 一种改进的基于广义内积的空时二维自适应处理方法
CN104730498A (zh) 基于Keystone和加权旋转FFT的目标检测方法
CN108333568B (zh) 冲击噪声环境下基于Sigmoid变换的宽带回波Doppler和时延估计方法
CN110109092B (zh) 多径环境下基于时间反演的雷达测速方法
CN104360251A (zh) 一种变压器局部放电的超声波信号时延估计方法
CN103116154A (zh) 基于杂波环境下的发射与接收联合优化自适应滤波方法
CN111580110A (zh) 一种基于浅海多径时延的复合码水声测距方法
CN112130138A (zh) 一种超视距雷达的电离层相位污染校正方法
US8358233B2 (en) Radar target detection process
CN108957416A (zh) 脉冲噪声环境下基于分数阶功率谱密度的线性调频信号参数估计方法
CN112748429B (zh) 一种快速噪声对消滤波方法
CN102508265A (zh) 基于信号分离估计理论的卫星导航信号多径干扰抑制方法
Nguyen A low complexity parameter estimation technique for LFMCW signals
AU2004248609A1 (en) Estimation of background noise and its effect on sonar range estimation
CN102707157B (zh) 一种基于功率谱的单频脉冲信号参数估值方法
CN104101868A (zh) 基于干扰子空间重构的雷达多假目标干扰抑制方法
CN106597390B (zh) 一种海杂波抑制方法
CN111025258B (zh) 一种针对雷达波形分集的联合失配滤波器及其设计方法
CN104793197A (zh) 基于ifft频谱相除法和梯度自适应格型滤波的直达波抑制方法
CN111600665B (zh) 基于arima滤波器的sage信道参数估计方法
CN112731292A (zh) 局部imf能量加权的低空飞行目标信号时延估计方法

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