CN109782249B - 一种两目标相关时延估计算法 - Google Patents
一种两目标相关时延估计算法 Download PDFInfo
- Publication number
- CN109782249B CN109782249B CN201910131733.1A CN201910131733A CN109782249B CN 109782249 B CN109782249 B CN 109782249B CN 201910131733 A CN201910131733 A CN 201910131733A CN 109782249 B CN109782249 B CN 109782249B
- Authority
- CN
- China
- Prior art keywords
- target
- targets
- time delay
- signal
- estimation algorithm
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种两目标相关时延估计算法,针对较为接近的两个目标,利用两目标相关匹配的方法,得到两目标时延的概率分布,在复加性高斯白噪声CAWGN信道下,针对复高斯散射目标,利用两目标相关匹配后进行峰值搜索的方法,对两目标时延,即两目标的位置进行估计。本发明的数值仿真结果验证了算法的可行性,本发明的结论对实际多目标雷达探测系统设计具有重要的理论指导意义。
Description
技术领域
本发明涉及信息传输与处理技术领域,尤其是一种两目标相关时延估计算法。
背景技术
雷达基本功能是分析目标的回波信号来找到目标并检测目标的空间位置。随着雷达信号和数据处理能里的快速提升,现代雷达的作用已不局限在对目标的检测和定位,还包含动目标测速、目标特征的识别和提取、实现目标的成像以及对运动目标的跟踪和实时预测等。
雷达最基本的用途是对目标进行探测。目标探测一直是雷达技术的重要领域,雷达目标探测技术是通过分析回波信号,获取目标的距离-状态-方位的多维特征,可以将其划分为一种算法研究。传统的雷达目标探测技术研究中,往往以发现目标和测量目标为基础,迄今为止,单目标探测技术已经相对完善,然而在实际的探测场景中,目标并不总是单独出现的,多目标探测的研究就显得至关重要。而多目标探测问题相较单目标复杂得多,除了单目标探测遇到的问题外,还需要研究目标间的干扰。因此,在多目标探测中多研究的是目标较远,能够清楚分辨情况下目标的位置,当目标距离较近时,由于探测环境的复杂多变性与干扰的多重性,以及目标相互之间的强干扰,使得临近目标探测到目前为止仍未有令人十分满意的结果,没有形成稳定的技术和方法。因此,寻求简便可行又有效的高分辨算法来估计临近目标位置很有必要。
发明内容
本发明所要解决的技术问题在于,提供一种两目标相关时延估计算法,解决了在目标较接近时单个匹配滤波器难以区分目标的问题。
为解决上述技术问题,本发明提供一种两目标相关时延估计算法,包括如下步骤:
(1)雷达系统接收端接收信号,将接收信号变频到基带,并通过带宽为B/2的理想低通滤波器,则接收信号为:
其中,s(t)表示实际的基带信号,yl表示第l个目标的反射系数,τl表示第l个目标的时延,L为目标总数,w(t)表示总功率为N0的复高斯白噪声;将接收信号r(t)以采样率B进行离散化,假设参考点是观测区间中点,观测区间为[-D/2,D/2),v表示信号传播速度,则由T=2D/v得时延区间为[-T/2,T/2),因此时间带宽积N=TB为归一化探测区间,归一化时延xl=Bτl,得到离散序列:
其中n=-N/2,...,N/2-1,
矢量形式为
z=U(x)y+w;
(2)在给定X条件下,Z的概率密度函数服从复高斯分布,可通过对Y和W求期望求得协方差矩阵R,特别计算当目标跟踪数L=2时,R的表达式;由协方差矩阵R可得概率密度p(z|x);
(3)假设位置X的先验概率密度为p(x)=1/N,由协方差矩阵的逆求得后验概率p(x|z);
(4)假设z0为两目标位置分别位于x1=x10,x2=x20处的接收信号,则可用协方差矩阵R0产生复高斯随机矢量z0,其中
R0=N0(ρ1 2u(x10)u(x10)H+ρ2 2u(x20)u(x20)H+I)
设定相关匹配滤波器两目标之间的相对距离为Δx,将归一化时延x1在探测区间内遍历,那么,x2=x1+Δx也随之遍历,并逐次与产生的接收信号z0进行相关匹配,距离差Δx需在设定的范围内遍历;然后进行峰值搜索,找到最大值时的x1,x2值即为目标的位置的一次估计,多次估计得到目标位置的概率分布。
优选的,步骤(1)中,基带信号s(t)为:s(t)=sinc(Bt)。
优选的,步骤(2)中,当目标跟踪数L=2时,可以求得协方差矩阵R为
R=EY,W[zzH]=N0(ρ1 2u(x1)u(x1)H+ρ2 2u(x2)u(x2)H+I)
其中ρ1 2,ρ2 2分别为目标1,目标2的信噪比,I为单位矩阵。
优选的,步骤(3)条件下,协方差的逆
后验概率为
优选的,步骤(4)中,采用matlab工具产生复高斯随机矢量z0,仿真两目标相关匹配后的位置分布。
本发明的有益效果为:在雷达探测目标中,利用两目标相关匹配方法,能够更准确地找到两个目标的位置,解决了目标较接近时单个匹配滤波器难以区分目标的问题,本发明提出的算法在多目标探测中有很好的理论指导意义。
附图说明
图1为本发明的方法流程示意图。
图2为信噪比为10dB,两目标距离为1的算法仿真示意图。
图3为信噪比为10dB,两目标距离为0.5的算法仿真示意图。
图4为信噪比为10dB,两目标距离为0.3的算法仿真示意图。
具体实施方式
如图1所示,一种两目标相关时延估计算法,包括如下步骤:
步骤1、将接收到的信号变频到基带,并通过带宽为B/2的理想低通滤波器,则接收信号为:
其中,s(t)表示实际的基带信号,yl表示第l个目标的反射系数,τl表示第l个目标的时延,L为目标总数,w(t)表示总功率为N0的复高斯白噪声,本例中s(t)选择sinc函数。将接收信号r(t)以采样率B进行离散化,假设参考点是观测区间中点,观测区间为[-D/2,D/2),v表示信号传播速度,则由T=2D/v得时延区间为[-T/2,T/2),因此时间带宽积N=TB为归一化探测区间,归一化时延xl=Bτl,得到离散序列:
其中n=-N/2,...,N/2-1。
为了简洁表示,可写成矢量形式
z=U(x)y+w
其中z=[z(-N/2),…,z(N/2-1)]T表示离散接收信号;U(x)=[u1(x)…uL(x)]表示时延信号,ul(x)=[sinc(-N/2-xl)…sinc(N/2-1-xl)]T;y=[y1…yL]T表示目标反射系数;w=[w(-N/2),…,w(N/2-1)]T表示噪声矢量。
步骤2、在给定X条件下,Z的概率密度函数服从复高斯分布。可通过对Y和W求期望求得协方差矩阵
R=EY,W[zzH]
当目标跟踪数L=2时,可以求得R为
R=E[(U(x)y+w)(U(x)y+w)H]
=U(x)E[yyH]UH(x)+E[wwH]
=N0(ρ1 2u(x1)u(x1)H+ρ2 2u(x2)u(x2)H+I)
其中ρ1 2,ρ2 2分别为目标1,目标2的信噪比,I为单位矩阵。
则给定X条件下,Z的概率密度为
步骤3、根据协方差矩阵求逆公式,得
假设目标位置在观测区间内均匀分布,即距离X的先验概率密度为p(x)=1/N,那么,后验概率为
步骤4、进一步的,本发明的一种两目标相关时延估计算法,假设z0为两目标位置分别位于x1=x10,x2=x20处的接收信号,则可用协方差矩阵R0产生复高斯随机矢量z0,其中
R0=N0(ρ1 2u(x10)u(x10)H+ρ2 2u(x20)u(x20)H+I)
设定相关匹配滤波器两目标之间的相对距离为Δx,将归一化时延x1在探测区间内遍历,那么,x2=x1+Δx也随之遍历,并逐次与产生的接收信号z0进行相关匹配,即
|(u(x1)+u(x1+Δx))z0|2
距离差Δx需在设定的范围内遍历。然后进行峰值搜索,找到最大值时的x1,x2值即为目标的位置的一次估计。多次估计得到目标位置的概率分布。此方法可以应用于信噪比较高,已追踪到存在两个临近目标时,对两目标距离的估计。
本发明中采用matlab工具仿真,仿真时设定两目标信噪比相同,数值仿真验证了本发明的正确性、可行性和有效性。图2为信噪比为10dB,两目标距离为1的算法仿真图;图3为信噪比为10dB,两目标距离为0.5的算法仿真图;图4为信噪比为10dB,两目标距离为0.3的算法仿真图。从图中可以看出,目标始终可以分辨,且可以根据峰值估计两目标距离差。
Claims (4)
1.一种两目标相关时延估计算法,其特征在于,包括如下步骤:
(1)雷达系统接收端接收信号,将接收信号变频到基带,并通过带宽为B/2的理想低通滤波器,则接收信号为:
其中,s(t)表示实际的基带信号,yl表示第l个目标的反射系数,τl表示第l个目标的时延,L为目标总数,w(t)表示总功率为N0的复高斯白噪声;将接收信号r(t)以采样率B进行离散化,假设参考点是观测区间中点,观测区间为[-D/2,D/2),v表示信号传播速度,则由T=2D/v得时延区间为[-T/2,T/2),因此时间带宽积N=TB为归一化探测区间,归一化时延xl=Bτl,得到离散序列:
其中n=-N/2,...,N/2-1,
矢量形式为
z=U(x)y+w
其中z=[z(-N/2),…,z(N/2-1)]T表示离散接收信号,U(x)=[u1(x)…uL(x)]表示时延信号,y=[y1…yL]T表示目标反射系数,w=[w(-N/2),…,w(N/2-1)]T表示噪声矢量;
(2)在给定X条件下,Z的概率密度函数服从复高斯分布,通过对Y和W求期望求得协方差矩阵R,计算当目标跟踪数L=2时,R的表达式;由协方差矩阵R得概率密度p(z|x);在给定X条件下,Z的概率密度函数服从复高斯分布;通过对Y和W求期望求得协方差矩阵
R=EY,W[zzH]
当目标跟踪数L=2时,求得R为
R=E[(U(x)y+w)(U(x)y+w)H]
=U(x)E[yyH]UH(x)+E[wwH]
=N0(ρ1 2u(x1)u(x1)H+ρ2 2u(x2)u(x2)H+I)
其中ρ1 2,ρ2 2分别为目标1,目标2的信噪比,I为单位矩阵;
(3)假设位置X的先验概率密度为p(x)=1/N,由协方差矩阵的逆求得后验概率p(x|z);则给定X条件下,Z的概率密度为
假设目标位置在观测区间内均匀分布,即距离X的先验概率密度为p(x)=1/N,那么,后验概率为
(4)假设z0为两目标位置分别位于x1=x10,x2=x20处的接收信号,则用协方差矩阵R0产生复高斯随机矢量z0,其中
R0=N0(ρ1 2u(x10)u(x10)H+ρ2 2u(x20)u(x20)H+I)
设定相关匹配滤波器两目标之间的相对距离为Δx,将归一化时延x1在探测区间内遍历,那么,x2=x1+Δx也随之遍历,并逐次与产生的接收信号z0进行相关匹配,距离差Δx需在设定的范围内遍历;然后进行峰值搜索,找到最大值时的x1,x2值即为目标的位置的一次估计,多次估计得到目标位置的概率分布。
2.如权利要求1所述的两目标相关时延估计算法,其特征在于,步骤(1)中,基带信号s(t)为:s(t)=sinc(Bt)。
4.如权利要求1所述的两目标相关时延估计算法,其特征在于,步骤(4)中,采用matlab工具产生复高斯随机矢量z0,仿真两目标相关匹配后的位置分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910131733.1A CN109782249B (zh) | 2019-02-22 | 2019-02-22 | 一种两目标相关时延估计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910131733.1A CN109782249B (zh) | 2019-02-22 | 2019-02-22 | 一种两目标相关时延估计算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109782249A CN109782249A (zh) | 2019-05-21 |
CN109782249B true CN109782249B (zh) | 2022-09-27 |
Family
ID=66485776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910131733.1A Active CN109782249B (zh) | 2019-02-22 | 2019-02-22 | 一种两目标相关时延估计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109782249B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110289922B (zh) * | 2019-07-09 | 2021-08-10 | 南京航空航天大学 | 一种时间调制及解调方法 |
CN111812591B (zh) * | 2020-06-16 | 2022-10-21 | 南京航空航天大学 | 一种基于贝叶斯估计的目标检测方法 |
CN117351042B (zh) * | 2023-10-08 | 2024-05-14 | 哈尔滨工业大学 | 一种基于高斯和滤波的任意步时延目标跟踪方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6044336A (en) * | 1998-07-13 | 2000-03-28 | Multispec Corporation | Method and apparatus for situationally adaptive processing in echo-location systems operating in non-Gaussian environments |
CN108562883B (zh) * | 2017-12-29 | 2022-06-10 | 南京航空航天大学 | 一种多载波雷达系统的最大似然距离估计算法 |
CN108957422B (zh) * | 2018-06-01 | 2022-07-29 | 电子科技大学 | 一种基于量化数据的云mimo雷达的根克拉美罗下界计算方法 |
CN109239683B (zh) * | 2018-08-28 | 2022-08-02 | 天津大学 | 宽带被动mimo雷达的克拉美罗界分析方法 |
KR102289882B1 (ko) * | 2018-12-31 | 2021-08-13 | 한국항공대학교산학협력단 | 잡음 기반 위상 변조 파형의 부엽 제거를 통한 위치 추정 성능 향상 기법 및 다중 경로 환경에서의 충돌 방지를 위한 광대역 레이다 시스템 |
-
2019
- 2019-02-22 CN CN201910131733.1A patent/CN109782249B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109782249A (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108761404B (zh) | 一种基于二次相位函数参数估计及补偿的改进算法 | |
CN108761418B (zh) | 一种基于微多普勒特征的雷达多径目标剔除方法 | |
CN109782249B (zh) | 一种两目标相关时延估计算法 | |
CN107561508B (zh) | 一种用于匀加速运动目标的相参积累检测方法 | |
CN107132534B (zh) | 一种高速雷达目标频域检测的优化方法 | |
US10222472B2 (en) | System and method for detecting heading and velocity of a target object | |
CN109324322A (zh) | 一种基于被动相控阵天线的测向与目标识别方法 | |
CN109164438A (zh) | 一种基于lfm相干脉冲串的到达时间和到达频率的联合估计方法 | |
CN106646447A (zh) | 基于线性调频连续波的雷达目标长时间积累检测方法 | |
CN110988834B (zh) | 一种基于自适应阈值双门限的脉冲到达时间测量方法 | |
CN112835005B (zh) | 一种基于超分辨目标跟踪的微多普勒特征提取方法 | |
CN104483668A (zh) | 一种高精度雷达信号检测和跟踪系统及其方法 | |
CN109655819A (zh) | 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 | |
CN109001671B (zh) | 一种跳频信号的目标检测和参数估计方法及装置 | |
CN109507654B (zh) | 一种基于ls的复杂环境下相位信息计算方法 | |
Hyun et al. | Method to improve range and velocity error using de-interleaving and frequency interpolation for automotive FMCW radars | |
CN108508433A (zh) | 基于模板匹配的无源探测系统快速时频同步方法 | |
CN110967677B (zh) | 时频域分级多目标分辨与测量方法 | |
Wei et al. | Detection and localization of high speed moving targets using a short-range UWB impulse radar | |
JP2003121544A (ja) | ドップラシフト周波数測定装置およびその利用装置 | |
CN109655804B (zh) | 一种基于奇异值分解的临近目标相对距离估计方法 | |
CN114280532A (zh) | 一种基于带内共轭点乘的雷达目标角度估计方法及系统 | |
CN108983192B (zh) | 基于gps辐射源的雷达运动目标参数估计方法 | |
CN115267721B (zh) | 一种基于双频sar的地面动目标径向速度估计方法 | |
CN108318857A (zh) | 基于分数阶傅里叶变换的多个非合作发射源被动定位方法 |
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 |