CN107356922B - 一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 - Google Patents
一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 Download PDFInfo
- Publication number
- CN107356922B CN107356922B CN201710761011.5A CN201710761011A CN107356922B CN 107356922 B CN107356922 B CN 107356922B CN 201710761011 A CN201710761011 A CN 201710761011A CN 107356922 B CN107356922 B CN 107356922B
- Authority
- CN
- China
- Prior art keywords
- frame
- target
- distance
- frequency
- radar
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005259 measurement Methods 0.000 claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 27
- 238000013139 quantization Methods 0.000 claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims abstract 2
- 230000007704 transition Effects 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 4
- 230000001629 suppression Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 abstract 1
- 238000009825 accumulation Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 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
- 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
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
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
技术领域
本发明属于雷达目标检测与跟踪技术领域,特别涉及多高重频雷达解距离模糊和多帧联合积累技术研究。
背景技术
高PRF雷达具有良好的多普勒分辨率,能够得到在多普勒维上不模糊的量测值,但是观测的距离信息会高度模糊。在雷达系统中常采用多重频工作的方式提高雷达测量的精度,同时对距离模糊问题进行求解。运用中国余数定理原理解模糊的算法时间复杂度高,其不仅要求各重频下进行求解的模糊距离两两互质,且未考虑各重频下所接收回波数据的量测和量化误差。因此,部分量测信息未能被正确解模糊,目标信息丢失严重,导致无法对目标实现正确跟踪;但若考虑各重频下对点迹解模糊所需的容许误差,则会产生虚假点迹,同样无法对目标进行精确跟踪。
检测前跟踪方法可利用目标回波间的相关性在多帧联合积累后对目标进行门限判决,可以实现对目标信号的增强和噪声杂波的抑制,抑制无关虚假点迹,增大信噪比。该技术可以有效避免因单帧处理不当导致的目标点迹丢失问题,改善对目标的检测性能。其中一种应用于雷达目标检测领域的经典方法是动态规划检测前跟踪方法,但该算法的量测模型是基于非模糊量测空间建立的,而高重频雷达体制下所接收到的量测信息距离高度模糊,回波数据无法反映目标的真实位置,因此不满足经典的动态规划检测前跟踪方法的假设,无法直接应用于高重频雷达。
在文献“A Modified DPA for Weak Target Detection via HPRF Radar,Radar(Radar)2011 IEEE CIE International Conference on,vol.2,pp.1737-1740,2011”中通过使用前一时刻预测得到的目标速度对DP积累时的状态转移区域进行修正,以达到对不连续航迹进行值函数积累的目的,但此模型只建立在单重频模式下,且过于依赖对目标初始状态的准确预测,当场景内存在多个模糊目标时,无法有效进行跟踪;此外,由于其需要对目标下一时刻由于距离模糊而可能折叠到的区域进行等概率预判,因此算法复杂度高,不易于工程实现。在文献“A new maneuvering target tracking method for HPRF radars,IET International Radar Conference 2013,Xi'an,2013,pp.1-6.”中考虑了机动目标在高重频雷达体制下的跟踪问题,将脉冲间隔数及其增量作为目标离散状态的一个变量,通过建立交互多模型实现对模糊距离的求解以及目标状态的跟踪;但其需预先得到IMM的模型权重参数,且对脉冲间隔数的准确估计与否将严重影响该算法跟踪性能。
发明内容
本发明的发明目的在于:针对上述存在的问题,提供一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法,解决了传统解模糊算法处理计算代价高、有效真实量测信息少,无法有效恢复目标航迹的问题。
本发明的解决方案是先在雷达最大探测距离内对待解模糊距离单元进行扩展,考虑距离和方位解模糊所需容许误差,对多重频量测数据建立辅助扩展表,根据M/N准则(即对N重频解模糊,解模糊表格中,某一距离单元在距离误差单元Mr内有M个点被标记则认为该距离单元为目标的真实距离)恢复目标真实距离,得到无模糊量测平面;然后联合多帧扫描进行动态规划积累及门限判决,最后恢复出目标真实航迹。该方法有效解决了传统解模糊算法处理后量测信息失真的问题,从而实现了对多高重频雷达体制下目标的准确跟踪。
本发明的扩表辅助法解模糊的多重频雷达目标检测前跟踪方法,包括下列步骤:
基于每帧的量测数据积累动态规划值函数:
首先,对当前帧的量测数据Zk作低门限预处理,保留过门限的点的数据信息,并将未过门限的点的数据值置0,其中第k帧量测 表示第k帧第i个重频接收到的量测数据,i=1,…,N,N为雷达重频个数,且其中表示第k帧第i个重频下的距离单元编号m与方位单元编号n对应量测数据的幅度值;Naz分别表示第i个重频下的距离单元总数、方位单元总数;
然后,对低门限预处理后的量测数据Zk扩表辅助法解模糊:
解距离模糊,输出无模糊量测:在解模糊表格中,若任意一列对应的距离单元在在距离误差单元数Mr(预设值)内至少有M个点被标记,则认为当前距离单元为目标的真实距离,其中M为预设判决阈值,即在解模糊表格的任意列m′的搜索范围内,若至少有M个点被标记,则将对应第m′列的距离单元标记为真实距离,其中第m′列的搜索范围为:第m′列及其Mr个近邻列;
最后,积累动态规划值函数:
用I(xk)表示目标状态xk对应的值函数,τ(xk)表示上一帧(即第k-1帧)所有可能转移到xk的目标状态集合,Ψ(xk)表示上一帧(即第k-1帧)的目标状态(即Ψ(xk)记录帧间目标状态转移关系),其中状态包括目标位置(距离、方位)及速度(距离向速度、方位向速度);
门限判决及航迹恢复,恢复目标航迹:
对观测的尾帧K的任意目标状态xk,如果有maxI(xK)>VT,其中VT为检测门限,则以作为对应尾帧K的观测时刻的估计状态,则利用Ψ(xK)记录的帧间状态转移关系恢复出目标状态xk对应的航迹若maxI(xK)<VT,则表示目标不存在。
通过上述步骤,完成对高重频PD雷达的解模糊以及目标检测跟踪的过程。
本发明的有益效果:本发明通过建立考虑量测和量化误差的扩表辅助法解距离模糊,使得尽可能多的模糊点迹能够成功解模糊,再联合多帧解模糊后的量测数据进行值函数的积累,门限判决后恢复出目标航迹,从而解决了在多高重频下对监测区域实现目标跟踪的问题。本发明的优点是解模糊过程简单、算法复杂度低,且同时利用检测前跟踪算法对解模糊阶段产生的虚假点迹进行有效抑制,而对具有运动特性的点迹进行能量积累,达到对目标航迹进行精确跟踪的目的。
附图说明
图1为本发明的流程图;
图2为本发明的算法示意图;
图3为扩表辅助法示意图,其中横坐标表示重频个数,纵坐标表示雷达最大探测范围内的距离单元数;
图4为各重频下对目标进行观测得到的运动轨迹;
图5为各重频下低门限CFAR(恒虚警处理)处理后点迹平面;
图6为解模糊后的目标跟踪结果与目标真实航迹的对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作进一步地详细描述。
实施例
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2013上验证正确。参见图1、2,本发明的具体实现过程包括下列步骤:
步骤1:系统参数初始化:
初始化系统参数包括:雷达最大探测距离Rmax=2000;雷达重频个数N=3;各重频之比PRF1:PRF2…PRFN=401:431:457;各重频下的最大不模糊距离Ru1=401,Ru2=431,Ru3=457;观测总帧数K=6;雷达距离分辨力Δr=1;雷达方位分辨力Δaz=1;雷达在各个重频下距离维分辨单元个数雷达在最大探测距离下距离维分辨单元个数为Nmax=2000;雷达方位维分辨单元个数Naz=50;利用蒙特卡洛仿真实验计算出检测门限VT=13.3250;当前帧数k=1;
步骤3:对接收到的数据作低门限预处理,保留过门限的点的数据信息,并将未过门限的点的数据值置0。
步骤4:扩表辅助法解模糊:
步骤4.2、扩展待解模糊点迹:
初始化解模糊表格C3×2000为全零平面,C3×2000为一个二维矩阵,行数代表各个重频,列代表距离单元数;初始化距离误差单元数Mr=1;初始化n=1;
步骤4.3、解距离模糊,输出无模糊量测:
取N个重频在相同方位单元数的待解模糊点迹,对各个重频下点迹的距离数在雷达最大探测距离内进行扩展,并标记解模糊表格中扩展后的距离数,其余仍为0;
根据M/N准则,即对N重频解模糊,解模糊表格中,某一距离单元在距离误差单元Mr内有M个点被标记则认为该距离单元为目标的真实距离;即对C3×2000任一距离单元m′(m′=1,…,Nmax),解模糊时对应的搜索范围为:[m′-Mr,m′+Mr],例如Mr=1,每次搜索范围包括3个距离单元,如图3所示,解模糊结果为:若搜索范围内[m′-Mr,m′+Mr]有至少M个点被标记则认为距离单元m′为目标的真实距离,即当搜索范围内的灰色背景个数大于或等于M时,则认为其为真实距离。对于C3×2000的第一个和最后一列,则分别通过补上Mr列全0来构成对应的搜索范围,例如图3中的第1和Nmax列,则分别补上一列全0来搜索,其实质等价于:若m′=1,则对应的搜索范围为:[m′,m′+Mr];若若m′=Nmax,则对应的搜索范围为:[m′+Mr,m′]。
步骤4.4:若n<Waz,令n=n+1,返回步骤4.3;
步骤5:积累动态规划值函数:
I(xk)表示目标目标状态xk对应的值函数,τ(xk)表示第k-1帧所有可能转移到xk的目标目标状态集合,目标状态xk对应的上一帧目标目标状态表示为Ψ(xk);
如果k=1,对每个目标状态对应的值函数赋初值,用各个目标状态对应的第一帧解模糊数据的幅度值初始化该目标状态的值函数;
如果2≤k≤K帧,更新每个目标状态对应的值函数:
步骤6:如果k<K,令k=k+1,返回步骤2;否则,执行步骤7;
步骤7:门限判决及航迹恢复,恢复目标航迹:
对第K帧的任意目标状态xk,如果有maxI(xK)>VT,其中VT为检测门限,则以作为K帧的观测时刻的估计目标状态,则利用Ψ(xK)记录的帧间目标状态转移关系恢复出目标状态xk对应的航迹若maxI(xK)<VT,则宣布目标不存在。得到的点航迹平面如图5所示。
通过上述步骤,可以在得到高重频PD雷达中的目标无模糊航迹,如图6所示。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (1)
1.一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法,其特征在于,包括下列步骤:
基于每帧的量测数据积累动态规划值函数:
对当前帧的量测数据Zk作低门限预处理,保留过门限的点的数据信息,并将未过门限的点的数据值置0,其中第k帧量测 表示第k帧第i个重频接收到的量测数据,i=1,…,N,N为雷达重频个数,且其中表示第k帧第i个重频下的距离单元编号m与方位单元编号n对应量测数据的幅度值;Naz分别表示第i个重频下的距离单元总数、方位单元总数;
对低门限预处理后的量测数据Zk扩表辅助法解模糊:
解距离模糊,输出无模糊量测:在解模糊表格的任意列m′的搜索范围内,若至少有M个点被标记,则将对应第m′列的距离单元标记为真实距离,其中M为预设判决阈值,第m′列的搜索范围为:第m′列及其Mr个近邻列,Mr为预设的距离误差单元数;
若当前帧不是观测首帧,则基于上一帧的值函数I(xk-1)更新当前帧的每个目标状态xk对应的值函数并更新目标Ψ(xk)为:其中,τ(xk)表示上一帧所有可能转移到目标状态xk的目标状态集合,arg max表示输出使目标函数最大时自变量的值;
门限判决及航迹恢复,恢复目标航迹:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710761011.5A CN107356922B (zh) | 2017-08-30 | 2017-08-30 | 一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710761011.5A CN107356922B (zh) | 2017-08-30 | 2017-08-30 | 一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107356922A CN107356922A (zh) | 2017-11-17 |
CN107356922B true CN107356922B (zh) | 2020-05-12 |
Family
ID=60290349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710761011.5A Expired - Fee Related CN107356922B (zh) | 2017-08-30 | 2017-08-30 | 一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107356922B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108536171B (zh) * | 2018-03-21 | 2020-12-29 | 电子科技大学 | 一种多约束下多无人机协同跟踪的路径规划方法 |
CN108519594B (zh) * | 2018-03-26 | 2020-01-21 | 电子科技大学 | 一种重频间帧间联合积累的多重频雷达检测前跟踪方法 |
CN109471091A (zh) * | 2018-11-07 | 2019-03-15 | 中国人民解放军国防科技大学 | 中重频pd雷达同时解模糊与检测跟踪方法 |
CN110221271B (zh) * | 2019-07-02 | 2023-03-03 | 中国航空工业集团公司雷华电子技术研究所 | 一种雷达干涉测角解模糊方法、装置及雷达系统 |
CN110456315B (zh) * | 2019-08-29 | 2023-03-07 | 西安电子工程研究所 | 一种基于位置预测的参差重频目标探测方法 |
CN111289953B (zh) * | 2020-01-14 | 2022-04-26 | 北京理工大学 | 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法 |
CN116755073B (zh) * | 2023-06-21 | 2024-03-26 | 上海雷骥电子科技有限公司 | 查找表解距离模糊的方法及应用 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103513244A (zh) * | 2013-10-17 | 2014-01-15 | 电子科技大学 | 一种基于动态规划的多帧相参积累目标检测前跟踪方法 |
CN104297748A (zh) * | 2014-10-20 | 2015-01-21 | 西安电子科技大学 | 一种基于轨迹增强的雷达目标检测前跟踪方法 |
CN104865570A (zh) * | 2015-05-11 | 2015-08-26 | 电子科技大学 | 一种快速的动态规划检测前跟踪方法 |
CN106772352A (zh) * | 2016-12-01 | 2017-05-31 | 中国人民解放军海军航空工程学院 | 一种基于Hough和粒子滤波的PD雷达扩展微弱目标检测方法 |
-
2017
- 2017-08-30 CN CN201710761011.5A patent/CN107356922B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103513244A (zh) * | 2013-10-17 | 2014-01-15 | 电子科技大学 | 一种基于动态规划的多帧相参积累目标检测前跟踪方法 |
CN104297748A (zh) * | 2014-10-20 | 2015-01-21 | 西安电子科技大学 | 一种基于轨迹增强的雷达目标检测前跟踪方法 |
CN104865570A (zh) * | 2015-05-11 | 2015-08-26 | 电子科技大学 | 一种快速的动态规划检测前跟踪方法 |
CN106772352A (zh) * | 2016-12-01 | 2017-05-31 | 中国人民解放军海军航空工程学院 | 一种基于Hough和粒子滤波的PD雷达扩展微弱目标检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107356922A (zh) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107356922B (zh) | 一种扩表辅助法解模糊的多重频雷达目标检测前跟踪方法 | |
CN107561508B (zh) | 一种用于匀加速运动目标的相参积累检测方法 | |
CN104076355B (zh) | 基于动态规划的强杂波环境中弱小目标检测前跟踪方法 | |
CN102901956B (zh) | 一种雷达微弱目标检测方法 | |
CN107576959B (zh) | 一种基于区域映射解模糊的高重频雷达目标检测前跟踪方法 | |
CN103439697B (zh) | 基于动态规划的目标检测方法 | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN108519594B (zh) | 一种重频间帧间联合积累的多重频雷达检测前跟踪方法 | |
CN112130142B (zh) | 一种复杂运动目标微多普勒特征提取方法及系统 | |
CN107544066B (zh) | 一种基于检测前跟踪分布式异步迭代滤波融合方法 | |
CN110007299B (zh) | 一种基于混合坐标伪谱技术的微弱目标检测跟踪方法 | |
CN107037424B (zh) | 一种基于序贯优化的多普勒雷达多帧相干检测前跟踪方法 | |
CN111505599B (zh) | 一种基于精确距离演化模型的相参积累检测方法 | |
CN105842688A (zh) | 单脉冲雷达对空目标快速捕获方法 | |
CN106772299B (zh) | 一种基于距离匹配的pd雷达微弱目标动态规划检测方法 | |
CN113534115A (zh) | 一种主被动联合处理的声纳目标高精度线谱提取方法 | |
CN116449320A (zh) | 频率捷变雷达体制下的长时间积累与参数估计方法 | |
CN104375128B (zh) | 一种基于互相关函数的高机动目标快速积累检测方法 | |
CN111796266B (zh) | 一种匀加速运动目标rd平面检测前跟踪方法 | |
CN106033120B (zh) | 一种多站雷达异步多帧联合检测方法 | |
JP5036392B2 (ja) | 目標検出装置 | |
CN105652256B (zh) | 一种基于极化信息的高频地波雷达tbd方法 | |
CN108828584B (zh) | 基于航迹折叠因子解模糊的多重频目标检测前跟踪方法 | |
CN116299387B (zh) | 非均匀杂波下干扰正交抑制的目标智能检测方法 | |
CN116520311A (zh) | 一种基于glmb的自适应航迹起始方法 |
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 | ||
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: 20200512 |