CN105137419B - 一种利用栅瓣增益的粒子滤波检测前跟踪方法 - Google Patents
一种利用栅瓣增益的粒子滤波检测前跟踪方法 Download PDFInfo
- Publication number
- CN105137419B CN105137419B CN201510569949.8A CN201510569949A CN105137419B CN 105137419 B CN105137419 B CN 105137419B CN 201510569949 A CN201510569949 A CN 201510569949A CN 105137419 B CN105137419 B CN 105137419B
- Authority
- CN
- China
- Prior art keywords
- target
- graing lobe
- signal
- particle
- tracking
- 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.)
- Expired - Fee Related
Links
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
- G01S13/00—Systems 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/66—Radar-tracking systems; Analogous systems
-
- 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
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
技术领域
本发明属于雷达栅瓣抑制和目标检测跟踪技术领域,特别涉及栅瓣抑制和低信噪比下的目标检测前跟踪技术领域。
背景技术
检测前跟踪是一种适用于低信噪比的情况下对目标进行检测和跟踪的信号处理技术。检测前跟踪技术在单帧内并不宣布检测结果,而是对数据进行多帧处理后同时宣布检测结果与目标的航迹。雷达系统中常用的检测前跟踪技术主要有基于动态规划的检测前跟踪方法、基于粒子滤波的检测前跟踪方法以及基于霍夫变换的检测前跟踪方法等。其中粒子滤波是递推贝叶斯估计的一种蒙特卡洛实现方法。它的主要思想是用大量采样的粒子来近似目标状态。
阵列天线是由多个天线单元按照一定规律进行排列的天线系统。天线阵列的形式多样,包括线阵、圆阵等。均匀线阵是一种实现简单,各阵列单元均匀排列在一条直线上的天线阵列。采用均匀线阵可以增加天线等效孔径,提高天线增益,从而提升方位分辨率。但是,随着阵元间距的增大,均匀线阵系统在获得更窄波束和更高方位分辨率等好处的同时,也带来了严重的栅瓣问题。栅瓣会造成目标遮挡、虚警率高以及距离模糊等影响,因此必须对栅瓣进行有效抑制。
目前,国内外学者已经对栅瓣的抑制方法做了大量的研究。对于均匀线阵雷达,一种常见的栅瓣抑制方法是采用调频信号抑制栅瓣。对于同一目标,由于不同频率的信号产生的栅瓣位置不同,同时所有信号在目标所在位置强度叠加,从而使栅瓣强度相对目标强度降低,达到抑制栅瓣的目的。但是该方法实现复杂,而且不能完全消除栅瓣的影响,同时浪费了分散到栅瓣位置的信号能量。
发明内容
本发明的目的是针对现有均匀线阵雷达中栅瓣抑制存在的问题,从目标的检测和跟踪角度,提供一种利用栅瓣增益的粒子滤波检测前跟踪方法。在系统初始化后,该方法首先对栅瓣与目标的位置关系进行计算,得到栅瓣位置和目标位置的对应关系;然后采用粒子滤波方法,对目标进行检测跟踪;在粒子滤波方法的具体实现中,利用栅瓣增益重新计算粒子权重,等效提高目标的信噪比,从而达到大幅度地提高检测跟踪精度的目的。
本发明提供了一种利用栅瓣增益的粒子滤波检测前跟踪方法,它包括以下步骤:
步骤1、初始化系统参数:
初始化系统参数包括:初始化监视范围对应的数据平面D、发射机个数M、接收机个数N、收发天线位置、发射调频信号起始频率f0、信号步进间隔△f、信号频点数Kf,粒子数目Np、观测总帧数K;
步骤2、计算目标与栅瓣的方位关系;
步骤3、接收第k帧回波信号,经脉冲压缩和BP成像处理后,形成当前帧高分辨率图像;
步骤4、粒子位置采样:其中为目标状态,p(·|·)表示求条件概率,n=1,2,…,Np是粒子标号;
步骤5、计算目标的栅瓣影响区域:
令表示目标位置在时,对应的栅瓣位置的空间集合;
其中表示向上取整,△x和△y分别表示X维和Y维的分辨率,表示到空间原点的距离,θg(xk)是目标xk的方位为θ0时此系统参数条件下对应栅瓣方位的集合;
步骤6、利用栅瓣对目标的影响,对采样粒子的权重进行更新;
其中,(xi,yj)表示量测单元(i,j)对应的空间坐标,表示对应的空间坐标;
步骤7、对粒子进行重采样;
步骤8、状态估计:判断出目标位置;
步骤9、如果k<K,令k:=k+1,返回步骤3。
进一步的,步骤2中计算目标与栅瓣的方位关系为:
其中d是等效接收天线阵列间隔,kf=0,1,...,Kf-1是步进频信号频点标号,是发射信号波长,θ0是目标相对于天线阵列中心的方位,是目标相对于天线阵列的中心;对于所有的{kf},记满足上式的集合为θg(xk),即θg(xk)是目标xk的方位为θ0时此系统参数条件下对应栅瓣方位的集合;
进一步的,所述步骤3中接收第k帧回波信号,经脉冲压缩和BP成像处理后,形成当前帧高分辨率图像Z(k),
Z(k)={zk(i,j)}
其中,1≤i≤NX,1≤j≤NY,i和j分别为图像X维和Y维的量化状态;NX是X维量化的单元个数,NY是Y维量化的单元个数;zk(i,j)表示第k帧图像数据的量测单元(i,j)中的量测值,为图像数据的幅度;
通过上面的步骤,就可以在有栅瓣影响的雷达回波数据平面中对目标进行高精度的检测跟踪。
本发明的创新点在于消除栅瓣对目标检测跟踪的不利影响;利用栅瓣增益提高目标检测跟踪精度;提供在信号处理阶段有效利用栅瓣信息的新思路。
本发明的优点在于利用了稀疏均匀线阵雷达中出现的栅瓣,实现对微弱目标的检测和跟踪,具有精度高、检测跟踪性能好的优点。
附图说明
图1为目标高分辨率图像中栅瓣位置示意图。
图2为本发明的流程图。
图3为没有执行上述利用栅瓣增益的步骤5与执行了该步骤的仿真结果对比。
其中:(a)为两种方法随信噪比变化的概率对比曲线;(b)为两种方法随帧数变化的均方根误差对比曲线。
具体实施方式
本发明主要采用计算机仿真的方法进行验证,所有步骤、结论都在MATLAB-R2012b上验证正确。具体实施步骤如下:
步骤1、初始化系统参数:
初始化系统参数包括:初始化监视范围对应的数据平面D、发射机个数M、接收机个数N、收发天线位置、发射调频信号起始频率f0、信号步进间隔△f、信号频点数Kf,粒子数目Np、观测总帧数K。
步骤2、计算目标与栅瓣的方位关系:
其中d是等效接收天线阵列间隔,kf=0,1,...,Kf-1是步进频信号频点标号,是发射信号波长,θ0是目标相对于天线阵列中心的方位,是目标相对于天线阵列的中心。对于所有的{kf},记满足上式的集合为θg(xk),即θg(xk)是目标xk的方位为θ0时,此系统参数条件下对应栅瓣方位的集合。
步骤3、生成高分辨率图像:
接收第k帧回波信号,经脉冲压缩和BP成像处理后,形成当前帧高分辨率图像Z(k),
Z(k)={zk(i,j)}
其中,1≤i≤NX,1≤j≤NY,i和j分别为图像X维和Y维的量化状态。NX是X维量化的单元个数,NY是Y维量化的单元个数。zk(i,j)表示第k帧图像数据的量测单元(i,j)中的量测值,为图像数据的幅度。
步骤4、粒子位置采样:其中p(·|·)表示求条件概率,n=1,2,…,Np是粒子标号;
步骤5、计算目标的栅瓣影响区域:
令表示目标位置在时,对应的栅瓣位置的空间集合。
其中表示向上取整,△x和△y分别表示X维和Y维的分辨率,表示到空间原点的距离。
步骤6、粒子权重更新:
其中,(xi,yj)表示量测单元(i,j)对应的空间坐标,表示对应的空间坐标。
步骤7、粒子重采样:
7.1.对于所有粒子n=1,2,…,Np,
7.2.重采样Np个粒子;
步骤8、状态估计:
步骤9、如果k<K,令k:=k+1,返回步骤3。
通过上面的步骤,就可以在有栅瓣影响的雷达回波数据平面中对目标进行高精度的检测跟踪。
在上述仿真中,没有执行计算栅瓣影响区域的步骤5与执行了该步骤的仿真对比结果如图3所示。由图3可知,在完成该方法的利用栅瓣增益之后,可以有效地提高目标的检测概率,降低目标跟踪的均方根误差。
通过本发明的具体实施可以看出,本发明在稀疏均匀线阵雷达有栅瓣影响的条件下,利用栅瓣增益,对目标能量进行积累,可以在信号处理阶段消除栅瓣不利影响,实现对微弱目标的检测和跟踪。
Claims (3)
1.一种利用栅瓣增益的粒子滤波检测前跟踪方法,包括以下步骤:
步骤1、初始化系统参数:
初始化系统参数包括:初始化监视范围对应的数据平面D、发射机个数M、接收机个数N、收发天线位置、发射调频信号起始频率f0、信号步进间隔Δf、信号频点数Kf,粒子数目Np、观测总帧数K;
步骤2、计算目标与栅瓣的方位关系;
步骤3、接收第k帧回波信号,经脉冲压缩和BP成像处理后,形成当前帧高分辨率图像;
步骤4、粒子位置采样:其中为目标状态,p(·|·)表示求条件概率,n=1,2,...,Np是粒子标号;
步骤5、计算目标状态的栅瓣影响区域:
令表示目标状态在时,对应的栅瓣位置的空间集合;
其中表示向上取整,Δx和Δy分别表示X维和Y维的分辨率,表示到空间原点的距离,θg(xk)是目标xk的方位为θ0时此系统参数条件下对应栅瓣方位的集合;
步骤6、利用栅瓣对目标的影响,对采样粒子的权重进行更新;
其中,Zk表示第k帧图像数据量测值,(xi,yj)表示量测单元(i,j)对应的空间坐标,表示对应的空间坐标;
步骤7、对粒子进行重采样;
步骤8、状态估计:判断出目标位置;
步骤9、如果k<K,令k:=k+1,返回步骤3。
2.如权利要求1所述的一种利用栅瓣增益的粒子滤波检测前跟踪方法,其特征在于步骤2中计算目标与栅瓣的方位关系为:
其中d是等效接收天线阵列间隔,kf=0,1,...,Kf-1是步进频信号频点标号,是发射信号波长,θ0是目标相对于天线阵列中心的方位,是目标相对于天线阵列的中心;对于所有的{kf},记满足上式的集合为θg(xk),即θg(xk)是目标xk的方位为θ0时此系统参数条件下对应栅瓣方位的集合。
3.如权利要求1所述的一种利用栅瓣增益的粒子滤波检测前跟踪方法,其特征在于所述步骤3中接收第k帧回波信号,经脉冲压缩和BP成像处理后,形成当前帧高分辨率图像Z(k),
Z(k)={zk(i,j)}
其中,1≤i≤NX,1≤j≤NY,i和j分别为图像X维和Y维的量化状态;NX是X维量化的单元个数,NY是Y维量化的单元个数;zk(i,j)表示第k帧图像数据的量测单元(i,j)中的量测值,为图像数据的幅度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510569949.8A CN105137419B (zh) | 2015-09-09 | 2015-09-09 | 一种利用栅瓣增益的粒子滤波检测前跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510569949.8A CN105137419B (zh) | 2015-09-09 | 2015-09-09 | 一种利用栅瓣增益的粒子滤波检测前跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105137419A CN105137419A (zh) | 2015-12-09 |
CN105137419B true CN105137419B (zh) | 2017-08-11 |
Family
ID=54722825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510569949.8A Expired - Fee Related CN105137419B (zh) | 2015-09-09 | 2015-09-09 | 一种利用栅瓣增益的粒子滤波检测前跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137419B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11169254B2 (en) * | 2017-10-06 | 2021-11-09 | National Chung-Shan Institute Of Science And Technology | Operation mode control method and radar system |
CN110376556B (zh) * | 2019-06-11 | 2021-05-11 | 杭州电子科技大学 | 一种基于锦标赛选择的双层粒子滤波检测前跟踪方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL1020287C2 (nl) * | 2002-04-02 | 2003-10-03 | Thales Nederland Bv | Werkwijze voor meerdoelendetectie, met name voor toepassing in rondzoekradars met meerbundelvorming in elevatie. |
CN101566689B (zh) * | 2009-06-10 | 2011-04-13 | 电子科技大学 | 一种消除三维合成孔径雷达图像栅瓣的方法 |
CN102043150A (zh) * | 2010-12-06 | 2011-05-04 | 电子科技大学 | 用于微弱目标检测的改进型粒子滤波检测前跟踪方法 |
CN104237853B (zh) * | 2014-09-16 | 2016-09-28 | 电子科技大学 | 一种针对多帧检测前跟踪点迹序列的粒子滤波方法 |
-
2015
- 2015-09-09 CN CN201510569949.8A patent/CN105137419B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105137419A (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10429488B1 (en) | System and method for geo-locating and detecting source of electromagnetic emissions | |
US8947294B1 (en) | Method and system for adaptively cancelling clutter from the sidelobes of a ground-based radar | |
CN108089167B (zh) | 一种合成孔径雷达跨脉冲干扰信号检测方法 | |
CN112051568B (zh) | 一种两坐标雷达的俯仰测角方法 | |
CN104569947A (zh) | 一种探地雷达的多通道数据实时融合方法 | |
CN109212500A (zh) | 一种基于稀疏重构的ka-stap杂噪协方差矩阵高精度估计方法 | |
US10001550B2 (en) | Method for automatic classification of radar objects | |
CN115598593B (zh) | 一种等长短基线高精度测向定位方法、系统、设备及终端 | |
CN108872977B (zh) | 基于单通道超宽带雷达的生命体双站协同探测方法 | |
Rao et al. | Long‐time coherent integration detection of weak manoeuvring target via integration algorithm, improved axis rotation discrete chirp‐Fourier transform | |
CN110146873B (zh) | 分布式非相参雷达的目标位置和速度的估计方法 | |
CN108414992B (zh) | 一种基于相位信息杂波图的目标检测方法 | |
CN111025256A (zh) | 一种机载雷达的微弱生命体征信号的检测方法及系统 | |
CN105137419B (zh) | 一种利用栅瓣增益的粒子滤波检测前跟踪方法 | |
CN111580099A (zh) | 基于联合熵值的穿墙成像雷达墙体杂波抑制方法 | |
CN111142096A (zh) | 一种基于网格划分的多基地雷达目标定位方法 | |
CN103983975A (zh) | 基于两种雷达的大气运动垂直速度检测方法及系统 | |
RU2665032C2 (ru) | Устройство распознавания воздушно-космических объектов в двухдиапазонных радиолокационных комплексах с активными фазированными антенными решетками (афар) | |
CN113534077B (zh) | 一种雷达辐射源威力反演方法、装置及电子设备 | |
CN113189562B (zh) | 一种基于高程的地形探测算法 | |
Jagadesh et al. | Modeling Target Detection and Performance Analysis of Electronic Countermeasures for Phased Radar. | |
CN105652256B (zh) | 一种基于极化信息的高频地波雷达tbd方法 | |
CN111859255A (zh) | 一种地形遮蔽影响下雷达探测范围计算方法 | |
US20110274217A1 (en) | Method and system for deinterleaving signals | |
RU2304789C1 (ru) | Способ радиолокационного сопровождения траектории объекта |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
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: 20170811 Termination date: 20210909 |