CN107450070A - 基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法 - Google Patents

基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法 Download PDF

Info

Publication number
CN107450070A
CN107450070A CN201710242213.9A CN201710242213A CN107450070A CN 107450070 A CN107450070 A CN 107450070A CN 201710242213 A CN201710242213 A CN 201710242213A CN 107450070 A CN107450070 A CN 107450070A
Authority
CN
China
Prior art keywords
mrow
target
msub
msubsup
mtd
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
Application number
CN201710242213.9A
Other languages
English (en)
Other versions
CN107450070B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710242213.9A priority Critical patent/CN107450070B/zh
Publication of CN107450070A publication Critical patent/CN107450070A/zh
Application granted granted Critical
Publication of CN107450070B publication Critical patent/CN107450070B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • G01S13/723Radar-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/726Multiple 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

基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法
技术领域
本发明属于雷达跟踪技术领域,涉及多目标跟踪的相控阵雷达波束和驻留时间资源联合管理技术研究。
背景技术
相控阵雷达作为一种先进的有源电扫阵列多功能雷达,它能同时形成一系列独立的波束,因而可以实现同时跟踪空域中的上百个目标,并且可以完成目标的搜索、辨识、跟踪、无源探测等多种功能。相控阵雷达天线通过计算机对移位器的控制来实现对雷达波束的控制,从而可以任意改变波束在目标上的驻留时间。由于这些优点,相控阵雷达在多目标跟踪上有着广泛的应用。
雷达资源管理就是根据不同的任务、目标特性和周围的环境对有限的系统资源进行分配,在保证系统性能的前提下,最有效地利用系统资源。然而,相控阵雷达的波束资源和驻留时间是有限的,如何根据周围环境和目标的特点自适应的实现雷达资源管理,成为现在我们研究的主要问题之一。在多目标跟踪环境下,就波束调度和驻留时间分配而言,不但要解决单目标跟踪时遇到的驻留时间分配问题,还有最重要的一点,波束指向问题,简单的说就是如何使雷达根据不同目标的跟踪需求来调整其波束指向和驻留时间。在文献“Adaptive Beam Scheduling Algorithm for an Agile Beam Radar in Multi-TargetTracking,International Conference on.IEEE,2006”中,作者提出了一种基于预测协方差来分配雷达波束的办法,但该方法只考虑了如何有效调度波束的问题,未考虑波束的驻留时间分配,解决的问题比较单一。文献“Simultaneous multi-beam resourceallocation scheme for multiple target tracking,63(12):3110-3122,2015”基于性能驱动联合分配了波束指向和发射功率,但该方法假设驻留时间平均,且希望目标跟踪性能越高越好,但实际跟踪中,只要目标跟踪精度满足一定要求即可,该方法并未考虑资源节约问题。因此,以节约资源为目的的波束和驻留时间联合分配具有一定应用价值意义。
发明内容
本发明为解决上述技术问题,提出了一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,解决相控阵雷达在执行跟踪任务时,由于波束数量少于目标数时波束调度不合理及平均分配波束驻留时间造成的资源浪费问题。
本发明采用的技术方案是:一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,包括:
S1、建立系统模型;具体为:设定相控阵雷达的位置为(x,y),能同时发射M个波束,在监控区域广泛分布Q个目标,且M≤Q;确定目标的初始位置以及目标的运动模式,确定驻留时间ΔTq,k,并计算目标的回波信号及回波信噪比;
其中,q表示目标序号,q=1,…,Q,k表示第k个跟踪时刻;
所述目标运动形式为:匀速直线运动或匀加速直线运动或协同转弯运动;
S2、根据步骤S1确定的目标运动形式得到目标运动模型,并根据步骤S1计算得到的目标的回波信号,得到目标观测模型;根据目标运动模型与目标观测模型,得到包含驻留时间的目标量测噪声协方差;
S3、根据步骤S2得到的目标运动模型和目标的量测噪声协方差,得到各目标的预测贝叶斯克拉美罗界;
根据得到贝叶斯克拉美罗界,确定各目标跟踪精度的准则;
建立以总驻留时间为目标函数,目标跟踪精度和波束为约束的资源管控优化函数,对资源管控优化函数进行求解,得到波束和驻留时间的分配结果;
S4、根据步骤S3得到的波束和驻留时间的分配结果,采用粒子滤波实现雷达对多目标的跟踪。
进一步地,步骤S1还包括:引入二元变量Γk为变量,表示资源分配前相控阵雷达波束分配结果;
其中,q=1,…,Q;为变量,表示资源分配前目标q上的波束指向。
进一步地,步骤S3所述的根据步骤S2得到的目标运动模型和目标的量测噪声协方差,得到各目标的预测贝叶斯克拉美罗界;具体为:
假设目标估计状态为则其满足:
其中,表示对目标状态和观测求数学期望;表示目标状态的贝叶斯信息矩阵,zq,k表示目标观测模型,表示目标运动模型,(·)T表示转置运算,表示对求倒数;
所述表达式为:
其中,表示目标的先验FIM,表示目标的数据FIM;
若在k时刻,目标q不被波束照射,则驻留时间为零且目标的数据FIM不存在,所以贝叶斯信息矩阵转化为:
进行求逆,得到目标的预测贝叶斯克拉美罗界为:
更进一步地,步骤S3所述根据得到贝叶斯克拉美罗界,确定各目标跟踪精度的准则;具体为:根据计算得到目标跟踪精度的指标;
其中,CCRLB(1,1)和CCRLB(3,3)分别表示贝叶斯克拉美罗界对角线上的第一个和第三个分量。
进一步地,步骤S3所述资源管控优化函数表达式为:
其中,表示Q个1组成的行向量,ΔTk表示各个目标上驻留时间组成的列向量,ΔTk=[ΔT1,k,ΔT2,k,…,ΔTQ,k]T;ΔTmin和ΔTmax分别为驻留时间ΔTq,k的上界和下界;Sk为k时刻用于跟踪的波束数量;Ttrack为k时刻用于跟踪的总时间;η表示目标跟踪门限。
进一步地,步骤S3所述的根据最优化算法,得到波束和驻留时间的分配结果,具体为:
首先对k时刻每个目标的先验求逆得:
其中,为目标的先验的倒数,表示预测跟踪性能的克拉美罗界;
然后令并将其作为衡量先验信息的准则;接着将和跟踪门限η进行比较,得到各目标的波束分配结果:
其中,表示预测跟踪精度的指标,分别表示预测跟踪性能克拉美罗界对角线上的第一个和第三个分量;为定值,表示目标q上的波束分配指向;根据各目标的波束分配指向,得到相控阵雷达波束分配结果
最后根据波束分配结果将资源管控优化函数转化为凸优化问题;
采用梯度投影法来求解,得到驻留时间分配最终得到波束和驻留时间分配结果
本发明的有益效果:本发明的方法,首先建立信号模型、目标动态模型和量测模型,初始化系统参数,然后,推导出带有被管理资源参数的贝叶斯克拉美罗界并将其作为目标跟踪性能的准则,将资源分配问题转化为一个最小化驻留时间,以目标跟踪性能为约束的数学优化模型,然后根据最优化算法理论,提出一个先解决波束指向再分配驻留时间的两步分解算法,最后根据得到的资源分配结果,采用粒子滤波实现雷达对多目标的跟踪;有效解决了波束数量有限和平均分配驻留时间带来的相控阵雷达资源浪费问题,从而实现了既能保证所有目标满足预定跟踪性能又使用于跟踪任务的资源消耗量有所减少,更合理的分配了系统资源;并且本发明可以应用于多目标跟踪的资源优化领域。
附图说明
图1是基于目标跟踪的相控阵雷达波束和驻留时间联合管理流程图。
图2是相控阵雷达多波束工作模式示意图。
图3是多目标航迹与雷达位置分布图。
图4是本方法与平均分配的目标5的BCRLB和RMSE对比图。
图5是用于跟踪的时间和总时间的对比图。
图6是波束分布结果图。
图7是脉冲数目分布结果图。
图8是波束分布结果图。
图9是脉冲数目分布结果图。
图10是资源分配后检测概率。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
图1给出了基于目标跟踪的相控阵雷达波束和驻留时间联合管理流程图,本发明的技术方案为:基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,相控阵雷达多波束工作模式如图2所示,本申请的技术方案包括以下步骤:
S1、建立系统模型,初始化系统参数
假设相控阵雷达的位置为(x,y),能同时发射M个波束,在监控区域广泛分布Q个目标,且M≤Q,即每个时刻最多有M个目标能被波束照射。该雷达系统对这些目标进行跟踪,目标运动形式可为:匀速直线运动或匀加速直线运动或协同转弯运动,但后两种运动形式比较复杂,为了便于处理,本实施例中假设每个目标匀速运动,则目标q(q=1,…,Q)的初始位置和速度分别为那么在k时刻,目标q的位置和速度分别为考虑到k时刻并不是所都有目标被照射,引入二元变量Γk为变量,表示资源分配前相控阵雷达波束分配结果;
其中,q=1,…,Q;为变量,表示资源分配前目标q上的波束指向。
假设雷达波束发射脉冲信号,脉冲重复周期为Tpri,若k时刻有Nq,k个脉冲照射到某一目标,则雷达波束在该目标的驻留时间为ΔTq,k=Nq,k·Tpri,若雷达平均发射功率为Pav,则每个脉冲的能量为En=PavTpri;在k时刻目标q上回波信号表达式可写为:
其中,sp(t)为发射信号的归一化复包络,αq,k为信道损失,wq,k(t,n)为噪声样本,hq,k为目标雷达散射截面积RCS:Radar-Cross Section,τq,k为回波信号的时延,fq,k为多普勒频移。
为了改善和提高雷达对目标的检测和跟踪性能,需要对回波样本进行相参积累来提高信噪比,本发明采用动目标检测(MTD,moving targets detection)技术来实现脉冲的相参积累。不考虑杂波,因此当目标位于某一特定距离单元内,对照射到目标上回波信号rq,k(t,n),n=1,…,Nq,k在慢时间维上进行离散傅里叶变换即可,得到其频谱:
理论上,当f=fq,k时,|Sq,k(t,f)|2的峰值为假设噪声样本wq,k(t,n)为彼此独立,且均值为零,方差为的高斯白噪声,则其总的噪声功率为得到回波信噪比SNRq,k(ΔTq,k):
从上式可以看出,回波信噪比和驻留时间成一定数值关系。
S2、根据步骤S1确定的目标运动模式得到目标运动模型,本申请中假设目标匀速运动,目标q在k时刻的状态为:则其动态方程可写为:
其中,Fk表示状态转移矩阵,表示过程噪声,且本申请中为均值为零,方差为Qq,k-1的高斯白噪声;其中Qq,k-1的表达式如下:
其中,I2为二阶单位矩阵,T0表示采样间隔,ρq表示目标q的过程噪声强度。
从步骤S1计算得到的目标的回波信号,得到目标观测模型;具体为:从回波信号中提取距离和角度信息,在第k个跟踪时刻,目标量测方程可表示为:
其中,量测距离和角度分别为:
根据目标运动模型与目标观测模型,得到包括驻留时间的各目标量测噪声协方差;具体如下:量测噪声vq,k为均值为零,方差为Rq,k的高斯白噪声,量测噪声与回波信噪比有关,其关系式可表示为:
如此,便得到了包含驻留时间ΔTq,k的目标量测噪声协方差Rq,k(ΔTq,k)表达式。
S3、根据步骤S2得到的目标运动模型和目标的量测噪声协方差,得到各目标的预测贝叶斯克拉美罗界;
由于贝叶斯克拉美罗界为目标状态估计最小均方误差MSE提供了一个下界,且具有一定的预测性。因此,本申请采用它作为跟踪性能的准则。假设目标估计状态为则其满足:
其中,表示对目标状态和观测求数学期望;表示目标状态的贝叶斯信息矩阵(BFIM),其可以分为两个部分:
其中,分别表示目标的先验FIM(Fisher Information Matrix,信息矩阵)和目标的数据FIM,它们具体表达式如下:
其中,的转置表示量测函数zq,k的雅克比行列式,
在k时刻,对于目标q,它可能不被波束照射,此种情况下,驻留时间为零且数据FIM不存在,所以BFIM可以写为
对上式求逆,求得预测贝叶斯克拉美罗界为:
根据得到贝叶斯克拉美罗界,确定各目标跟踪精度的准则;因为目标贝叶斯克拉美罗界的对角线元素可反映目标状态向量各个分量估计方差的下界,所以可将下式作为各个目标跟踪精度的指标:
其中,CCRLB(1,1)和CCRLB(3,3)分别表示贝叶斯克拉美罗界对角线上的第一个和第三个分量,因此为目标位置估计的最小均方根误差(RMSE),即跟踪精度,提供了一个下界,并将它作为跟踪精度的准则。
建立以总驻留时间为目标函数,目标跟踪精度和波束为约束的资源管控优化函数,根据最优化算法,得到波束和驻留时间的分配结果;
在保证所有目标满足一定跟踪精度,即给定跟踪门限η的情况下,要使总的波束驻留时间消耗最少,确定目标函数为波束在目标上总的驻留时间结合波束和驻留时间约束,此优化问题可以建立为:
其中,表示Q个1组成的行向量,ΔTk表示各个目标上驻留时间组成的列向量,ΔTk=[ΔT1,k,ΔT2,k,…,ΔTQ,k]T;Sk为k时刻用于跟踪的波束数量;ΔTmin和ΔTmax分别为驻留时间ΔTq,k的上界和下界;Ttrack为k时刻用于跟踪的总时间;为变量,表示资源分配前目标q上的波束指向;η表示目标跟踪门限。
式(16)即为带有资源变量与ΔTq,k的优化问题模型。
尽管目标函数和其它约束是线性、凸的,但由于是个二元变量,因此该优化问题是非凸的。针对此非线性、非凸的优化问题,本申请提出一个两步分解算法来实现波束和驻留时间的分配,具体包括以下步骤:
1、基于目标先验FIM分配雷达波束指向,实际上,分配波束就是选择哪些目标被照射。
首先对k时刻每个目标的先验求逆得:然后令并将作为衡量先验信息的准则;接着将和跟踪门限η进行比较:
其中,为目标的先验的倒数,表示预测跟踪性能的克拉美罗界;表示预测跟踪精度的指标,分别表示预测跟踪性能克拉美罗界对角线上的第一个和第三个分量;为定值,表示目标q上的波束分配指向;根据各目标的波束分配指向,得到相控阵雷达波束分配结果
最终,目标跟踪精度的准则可以写为:
从而可得到相控阵雷达波束的分配结果
2、基于第一步得到的波束分配结果可将式(16)转成一个如下的凸优化问题:
对于此凸优化问题,本申请采用梯度投影法来求解,得到驻留时间分配最终通过提出的两部分解算法得到波束和驻留时间分配结果
S4、根据资源优化结果,采用粒子滤波实现多目标的跟踪;由于观测模型是非线性的,本发明采用粒子滤波算法来实现对目标的状态估计。在第k时刻,如果目标q被波束照射,则该目标的后验概率密度分布为
其中,为目标状态的预测概率密度函数,为似然函数,其表达式如下:
如果目标未被照射,则其后验状态为预测状态。由此可见只要得到雷达波束分布情况和波束驻留时间就可以得出每个目标的状态估计,实现多目标跟踪。
下面根据一个MATLAB仿真例子,对本发明的内容作进一步说明。
S1、假设雷达位置为(0,0),考虑到计算复杂度,本实施例用三个波束去跟踪五个目标,各目标初始状态及目标距离雷达的距离如表1所示,并估计出目标第一帧状态。假设脉冲重复周期Tpri=1ms,每个跟踪时刻用于跟踪的总时间为Ttrack=0.5s,波束驻留时间的约束为0.005Ttrack≤ΔTq,k≤0.9Ttrack,发射信号的载频和有效边带分别为fq,k=1GHz和βq,k=1MHz。根据初始参数利用式(4)计算出回波信噪比。
表1目标初始状态及距离雷达距离
目标标号 1 2 3 4 5
目标位(km) (150,50) (-150,100) (-30,200) (-130,200) (195,190)
目标速(km/s) (0.2,0.38) (0.5,-0.2) (0.6,0.05) (-0.35,-0.4) (-0.35,0.3)
距离(km) 158 178 200 238 272
S2、根据式(5)计算出目标每一帧的预测状态,根据回波信噪比,式(7)和式(9)计算出每个目标的量测噪声协方差Rq,k(ΔTq,k)。目标的航迹与雷达位置分布图如图3所示。
S3、将目标的预测状态和量测噪声协方差Rq,k(ΔTq,k)代入式(16)求出每个目标的预测贝叶斯克拉美罗界,根据式(15)得出目标跟踪精度准则的计算表达式,将计算得到的目标跟踪精度的准则代入式(16),形成需要的资源优化模型,根据提出的两步分解算法步骤,最终得出波束和驻留时间分配结果。具体仿真分析如下:
本实施例将目标跟踪门限设为η=200m,目标的过程噪声统一设为S1:[ρ1,...,ρ5]=[0.022,0.022,0.022,0.022,0.022],目标RCS统一设为均匀分布H:hk=[h1,k,h2,k,...,h5,k]T=[1,1,...,1]T,即未考虑过程噪声和目标RCS对资源分配的影响,所以每个目标的资源分布只跟目标到雷达的距离有关。
图4给出了本方法与平均分配的目标5的BCRLB和RMSE对比图,可看出在资源均匀分配的操作模式下,平均分配的目标5未能达到预定跟踪精度,可见本发明在维持远目标的跟踪性能方面更有效。
图5给出了本发明方法用于跟踪的时间和总时间的对比图,说明本方法在保证所有目标跟踪精度的情况下还可以节省资源。
从图6和图7可以看出,目标5一直被波束照射,且照射到目标5上的脉冲数一直高于其他目标,这是因为他距离雷达位置最远。其他目标上的脉冲数目根据它与雷达的位置依次变化,距离雷达越近,脉冲数越少。
为了探讨过程噪声资源分配的影响,本实施例设置了第二种过程噪声模型S2:[ρ1,...,ρ5]=[0.12,0.022,0.22,0.0052,0.0052]。同时在目标跟踪门限为η=200m,RCS为H1分布的情况下,波束和脉冲分布结果分别为图8和图9所示。对比图6和图7可以看出,由于目标1和目标3的过程噪声变大,机动性变强,它们一直被波束照射且脉冲数目明显增多。目标4和5由于过程噪声的变小,被波束照射的次数有所减少。
图10为分别将跟踪门限设为η1=180m,η2=200m,η3=240m时的资源分配后,雷达波束检测概率的变化,可以看出,和将总的跟踪资源用尽相比,本申请提出的办法更加节约资源,且跟踪门限越低,即对跟踪性能的要求越低,本申请提出的办法节约的资源越多,波束检测概率越高。
S4、将通过步骤S3得到的资源结果代入式(21)得到似然函数,将目标预测状态和似然函数代入式(20),从而得到每个目标的估计状态。最后的目标估计航迹如图3所示。
通过本发明具体实施方式可以看出,和传统的固定资源分配相比,本发明可以在保证所有目标跟踪精度的前提下,使相控阵雷达系统用于跟踪任务的资源消耗量有所减少,大概节约了15%~20%的资源,而且能保证较远目标的跟踪性能。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。

Claims (6)

1.一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,其特征在于,包括:
S1、建立系统模型;具体为:设定相控阵雷达的位置为(x,y),能同时发射M个波束,在监控区域广泛分布Q个目标,且M≤Q;确定目标的初始位置以及目标的运动模式,确定驻留时间ΔTq,k,并计算目标的回波信号及回波信噪比;
其中,q表示目标序号,q=1,…,Q,k表示第k个跟踪时刻;
所述目标运动形式为:匀速直线运动或匀加速直线运动或协同转弯运动;
S2、根据步骤S1确定的目标运动形式得到目标运动模型,并根据步骤S1计算得到的目标的回波信号,得到目标观测模型;根据目标运动模型与目标观测模型,得到包含驻留时间的目标量测噪声协方差;
S3、根据步骤S2得到的目标运动模型和目标的量测噪声协方差,得到各目标的预测贝叶斯克拉美罗界;
根据得到贝叶斯克拉美罗界,确定各目标跟踪精度的准则;
建立以总驻留时间为目标函数,目标跟踪精度和波束为约束的资源管控优化函数,对资源管控优化函数进行求解,得到波束和驻留时间的分配结果;
S4、根据步骤S3得到的波束和驻留时间的分配结果,采用粒子滤波算法实现雷达对多目标的跟踪。
2.根据权利要求1所述的一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,其特征在于,步骤S1还包括:引入二元变量Γk为变量,表示资源分配前相控阵雷达波束分配结果;
其中,q=1,…,Q;为变量,表示资源分配前目标q上的波束指向。
3.根据权利要求1所述的一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,其特征在于,步骤S3所述的根据步骤S2得到的目标运动模型和目标的量测噪声协方差,得到各目标的预测贝叶斯克拉美罗界;具体为:
假设目标估计状态为则其满足:
其中,表示对目标状态和观测求数学期望;表示目标状态的贝叶斯信息矩阵,zq,k表示目标观测模型,表示目标运动模型,(·)T表示转置运算,表示对求倒数;
所述表达式为:
其中,表示目标的先验FIM,表示目标的数据FIM;
若在k时刻,目标q不被波束照射,则驻留时间为零且目标的数据FIM不存在,所以贝叶斯信息矩阵转化为:
进行求逆,得到目标的预测贝叶斯克拉美罗界为:
<mrow> <msub> <mi>C</mi> <mrow> <mi>C</mi> <mi>R</mi> <mi>L</mi> <mi>B</mi> </mrow> </msub> <mrow> <mo>(</mo> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mi>q</mi> </msubsup> <mo>,</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mo>|</mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>q</mi> </msubsup> </msub> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mi>J</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mi>q</mi> </msubsup> <mo>,</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <msub> <mo>|</mo> <msubsup> <mi>x</mi> <mi>k</mi> <mi>q</mi> </msubsup> </msub> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>.</mo> </mrow>
4.根据权利要求3所述的一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,其特征在于,步骤S3所述根据得到贝叶斯克拉美罗界,确定各目标跟踪精度的准则;具体为:根据计算得到目标跟踪精度的指标;
其中,CCRLB(1,1)和CCRLB(3,3)分别表示贝叶斯克拉美罗界对角线上的第一个和第三个分量。
5.根据权利要求4所述的一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,其特征在于,步骤S3所述资源管控优化函数表达式为:
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mi>min</mi> </mtd> <mtd> <mrow> <msubsup> <mn>1</mn> <mi>Q</mi> <mi>T</mi> </msubsup> <msub> <mi>&amp;Delta;T</mi> <mi>k</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mi>q</mi> </msubsup> <mo>=</mo> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>}</mo> <mo>,</mo> <mi>q</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mi>Q</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mn>1</mn> <mi>Q</mi> <mi>T</mi> </msubsup> <msub> <mi>&amp;Gamma;</mi> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>S</mi> <mi>k</mi> </msub> <mo>&amp;le;</mo> <mi>M</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>,</mo> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mi>q</mi> </msubsup> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <mn>0</mn> <mo>,</mo> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mi>q</mi> </msubsup> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>F</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mi>q</mi> </msubsup> <mo>,</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mi>&amp;eta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mn>1</mn> <mi>Q</mi> <mi>T</mi> </msubsup> <msub> <mi>&amp;Delta;T</mi> <mi>k</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>T</mi> <mrow> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>c</mi> <mi>k</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,表示Q个1组成的行向量,ΔTk表示各个目标上驻留时间组成的列向量,ΔTk=[ΔT1,k,ΔT2,k,…,ΔTQ,k]T;ΔTmin和ΔTmax分别为驻留时间ΔTq,k的上界和下界;Sk为k时刻用于跟踪的波束数量;Ttrack为k时刻用于跟踪的总时间;η表示目标跟踪门限。
6.根据权利要求5所述的一种基于目标跟踪的相控阵雷达波束和驻留时间联合分配办法,其特征在于,步骤S3所述的根据最优化算法,得到波束和驻留时间的分配结果,具体为:
首先对k时刻每个目标的先验求逆得:
其中,为目标的先验的倒数,表示预测跟踪性能的克拉美罗界;
然后令并将其作为衡量先验信息的准则;接着将和跟踪门限η进行比较,得到各目标的波束分配结果:
其中,表示预测跟踪精度的指标,分别表示预测跟踪性能克拉美罗界对角线上的第一个和第三个分量;为定值,表示目标q上的波束分配指向;根据各目标的波束分配指向,得到相控阵雷达波束分配结果
最后根据波束分配结果将资源管控优化函数转化为凸优化问题;
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <munder> <mi>min</mi> <mrow> <msub> <mi>&amp;Delta;T</mi> <mi>k</mi> </msub> </mrow> </munder> </mtd> <mtd> <mrow> <msubsup> <mn>1</mn> <mi>Q</mi> <mi>T</mi> </msubsup> <msub> <mi>&amp;Delta;T</mi> <mi>k</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>&amp;le;</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>,</mo> <msubsup> <mi>&amp;Gamma;</mi> <mi>k</mi> <mrow> <mi>o</mi> <mi>p</mi> <mi>t</mi> </mrow> </msubsup> <mo>=</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>F</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>q</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mi>&amp;eta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mn>1</mn> <mi>Q</mi> <mi>T</mi> </msubsup> <msub> <mi>&amp;Delta;T</mi> <mi>k</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>T</mi> <mrow> <mi>t</mi> <mi>r</mi> <mi>a</mi> <mi>c</mi> <mi>k</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
采用梯度投影法来求解,得到驻留时间分配最终得到波束和驻留时间分配结果
CN201710242213.9A 2017-04-14 2017-04-14 基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法 Active CN107450070B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710242213.9A CN107450070B (zh) 2017-04-14 2017-04-14 基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710242213.9A CN107450070B (zh) 2017-04-14 2017-04-14 基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法

Publications (2)

Publication Number Publication Date
CN107450070A true CN107450070A (zh) 2017-12-08
CN107450070B CN107450070B (zh) 2020-02-21

Family

ID=60486272

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710242213.9A Active CN107450070B (zh) 2017-04-14 2017-04-14 基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法

Country Status (1)

Country Link
CN (1) CN107450070B (zh)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107656264A (zh) * 2017-08-02 2018-02-02 南京航空航天大学 杂波环境下的机会阵雷达多目标跟踪的功率资源管理方法
CN108288107A (zh) * 2017-07-03 2018-07-17 上海机电工程研究所 多传感器协同探测任务规划方法
CN108333583A (zh) * 2018-01-22 2018-07-27 西安电子科技大学 基于相控阵雷达搜索和跟踪双目标优化的资源分配方法
CN108802734A (zh) * 2018-06-20 2018-11-13 清华大学 一种控制雷达系统时序同步的方法及装置
CN108830007A (zh) * 2018-06-28 2018-11-16 电子科技大学 一种基于二次代价函数的雷达发射功率分配方法
CN109270512A (zh) * 2018-10-31 2019-01-25 北京无线电测量研究所 分布式相参合成雷达的调度方法及系统
CN109587802A (zh) * 2019-01-21 2019-04-05 中国电子科技集团公司第二十九研究所 一种基于相控阵捷变波束的上行信道资源分配方法
CN110673131A (zh) * 2019-11-25 2020-01-10 电子科技大学 多波束集中式mimo雷达时空资源-波形选择管理方法
CN111025275A (zh) * 2019-11-21 2020-04-17 南京航空航天大学 基于射频隐身的多基雷达辐射参数多目标联合优化方法
CN111090079A (zh) * 2019-12-24 2020-05-01 中国航天科工集团八五一一研究所 基于无源传感器协同的雷达组网辐射间隔优化控制方法
CN111323773A (zh) * 2020-02-20 2020-06-23 南京航空航天大学 基于射频隐身的组网雷达功率与带宽联合优化分配方法
CN112068124A (zh) * 2020-08-20 2020-12-11 南京航空航天大学 面向低截获的组网雷达驻留时间与辐射功率联合优化方法
CN112859064A (zh) * 2021-01-18 2021-05-28 中国船舶重工集团公司第七二四研究所 一种被动相控阵雷达自适应辐射源跟踪调度方法
CN113093172A (zh) * 2021-03-15 2021-07-09 南京航空航天大学 面向目标跟踪的多平台雷达路径与波形选择联合优化方法
CN114442082A (zh) * 2022-01-25 2022-05-06 中国船舶重工集团公司第七二四研究所 一种基于辐射源截获参数的变驻留时长目标跟踪方法
CN114609589A (zh) * 2022-03-09 2022-06-10 电子科技大学 一种基于启发式回溯的实时相控阵雷达波束驻留调度方法
CN115052295A (zh) * 2022-05-12 2022-09-13 电子科技大学 一种实时的相控阵雷达组网系统的资源管理方法
CN115993596A (zh) * 2023-03-27 2023-04-21 中国人民解放军63921部队 一种特性参数测量雷达资源调配方法、装置及计算机存储介质
CN117687018A (zh) * 2024-01-31 2024-03-12 四川省华盾防务科技股份有限公司 一种相控阵波束跟踪方法、装置、系统及存储介质
CN114442082B (zh) * 2022-01-25 2024-10-25 中国船舶集团有限公司第七二四研究所 一种基于辐射源截获参数的变驻留时长目标跟踪方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103913730A (zh) * 2014-04-10 2014-07-09 四川九洲电器集团有限责任公司 利用异步滑动式mtd提升相控阵雷达的探测能力的方法
CN104007419A (zh) * 2014-06-03 2014-08-27 西安电子科技大学 关于驻留时间与重访间隔的雷达时间资源联合分配方法
CN105116381A (zh) * 2015-08-12 2015-12-02 西安电子科技大学 多波束雷达时间功率资源联合分配方法
CN105182317A (zh) * 2015-08-20 2015-12-23 电子科技大学 一种基于集中式mimo雷达搜索模式下的资源管理方法
CN105785354A (zh) * 2016-03-11 2016-07-20 中国人民解放军国防科学技术大学 一种相控阵雷达探测高速运动目标的搜索参数优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103913730A (zh) * 2014-04-10 2014-07-09 四川九洲电器集团有限责任公司 利用异步滑动式mtd提升相控阵雷达的探测能力的方法
CN104007419A (zh) * 2014-06-03 2014-08-27 西安电子科技大学 关于驻留时间与重访间隔的雷达时间资源联合分配方法
CN105116381A (zh) * 2015-08-12 2015-12-02 西安电子科技大学 多波束雷达时间功率资源联合分配方法
CN105182317A (zh) * 2015-08-20 2015-12-23 电子科技大学 一种基于集中式mimo雷达搜索模式下的资源管理方法
CN105785354A (zh) * 2016-03-11 2016-07-20 中国人民解放军国防科学技术大学 一种相控阵雷达探测高速运动目标的搜索参数优化方法

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108288107A (zh) * 2017-07-03 2018-07-17 上海机电工程研究所 多传感器协同探测任务规划方法
CN107656264A (zh) * 2017-08-02 2018-02-02 南京航空航天大学 杂波环境下的机会阵雷达多目标跟踪的功率资源管理方法
CN107656264B (zh) * 2017-08-02 2021-04-27 南京航空航天大学 杂波环境下的机会阵雷达多目标跟踪的功率资源管理方法
CN108333583A (zh) * 2018-01-22 2018-07-27 西安电子科技大学 基于相控阵雷达搜索和跟踪双目标优化的资源分配方法
CN108333583B (zh) * 2018-01-22 2021-08-03 西安电子科技大学 基于相控阵雷达搜索和跟踪双目标优化的资源分配方法
CN108802734A (zh) * 2018-06-20 2018-11-13 清华大学 一种控制雷达系统时序同步的方法及装置
CN108830007A (zh) * 2018-06-28 2018-11-16 电子科技大学 一种基于二次代价函数的雷达发射功率分配方法
CN109270512A (zh) * 2018-10-31 2019-01-25 北京无线电测量研究所 分布式相参合成雷达的调度方法及系统
CN109587802B (zh) * 2019-01-21 2020-11-24 中国电子科技集团公司第二十九研究所 一种基于相控阵捷变波束的上行信道资源分配方法
CN109587802A (zh) * 2019-01-21 2019-04-05 中国电子科技集团公司第二十九研究所 一种基于相控阵捷变波束的上行信道资源分配方法
CN111025275A (zh) * 2019-11-21 2020-04-17 南京航空航天大学 基于射频隐身的多基雷达辐射参数多目标联合优化方法
CN111025275B (zh) * 2019-11-21 2021-10-08 南京航空航天大学 基于射频隐身的多基雷达辐射参数多目标联合优化方法
CN110673131A (zh) * 2019-11-25 2020-01-10 电子科技大学 多波束集中式mimo雷达时空资源-波形选择管理方法
CN111090079A (zh) * 2019-12-24 2020-05-01 中国航天科工集团八五一一研究所 基于无源传感器协同的雷达组网辐射间隔优化控制方法
CN111090079B (zh) * 2019-12-24 2023-10-13 中国航天科工集团八五一一研究所 基于无源传感器协同的雷达组网辐射间隔优化控制方法
CN111323773A (zh) * 2020-02-20 2020-06-23 南京航空航天大学 基于射频隐身的组网雷达功率与带宽联合优化分配方法
CN112068124A (zh) * 2020-08-20 2020-12-11 南京航空航天大学 面向低截获的组网雷达驻留时间与辐射功率联合优化方法
CN112859064A (zh) * 2021-01-18 2021-05-28 中国船舶重工集团公司第七二四研究所 一种被动相控阵雷达自适应辐射源跟踪调度方法
CN112859064B (zh) * 2021-01-18 2024-08-16 中国船舶集团有限公司第七二四研究所 一种被动相控阵雷达自适应辐射源跟踪调度方法
CN113093172A (zh) * 2021-03-15 2021-07-09 南京航空航天大学 面向目标跟踪的多平台雷达路径与波形选择联合优化方法
CN113093172B (zh) * 2021-03-15 2023-12-05 南京航空航天大学 面向目标跟踪的多平台雷达路径与波形选择联合优化方法
CN114442082A (zh) * 2022-01-25 2022-05-06 中国船舶重工集团公司第七二四研究所 一种基于辐射源截获参数的变驻留时长目标跟踪方法
CN114442082B (zh) * 2022-01-25 2024-10-25 中国船舶集团有限公司第七二四研究所 一种基于辐射源截获参数的变驻留时长目标跟踪方法
CN114609589A (zh) * 2022-03-09 2022-06-10 电子科技大学 一种基于启发式回溯的实时相控阵雷达波束驻留调度方法
CN114609589B (zh) * 2022-03-09 2023-08-11 电子科技大学 一种基于启发式回溯的实时相控阵雷达波束驻留调度方法
CN115052295A (zh) * 2022-05-12 2022-09-13 电子科技大学 一种实时的相控阵雷达组网系统的资源管理方法
CN115052295B (zh) * 2022-05-12 2023-10-24 电子科技大学 一种实时的相控阵雷达组网系统的资源管理方法
CN115993596A (zh) * 2023-03-27 2023-04-21 中国人民解放军63921部队 一种特性参数测量雷达资源调配方法、装置及计算机存储介质
CN117687018A (zh) * 2024-01-31 2024-03-12 四川省华盾防务科技股份有限公司 一种相控阵波束跟踪方法、装置、系统及存储介质
CN117687018B (zh) * 2024-01-31 2024-04-19 四川省华盾防务科技股份有限公司 一种相控阵波束跟踪方法、装置、系统及存储介质

Also Published As

Publication number Publication date
CN107450070B (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN107450070B (zh) 基于目标跟踪的相控阵雷达波束和驻留时间联合分配方法
CN107728139B (zh) 一种基于多目标跟踪的相控阵雷达组网系统资源管理方法
Yuan et al. Scaled accuracy based power allocation for multi-target tracking with colocated MIMO radars
Zhang et al. Antenna selection for target tracking in collocated MIMO radar
Yan et al. Cooperative target assignment and dwell allocation for multiple target tracking in phased array radar network
CN103760556B (zh) 基于集中式mimo雷达的多目标认知跟踪方法
CN106199579A (zh) 分布式mimo雷达目标跟踪精度联合资源优化方法
CN106777679B (zh) 一种基于脉冲交错的isar成像雷达资源自适应调度方法
CN114442088A (zh) 基于雷达的目标组生成
CN110780290B (zh) 基于lstm网络的多机动目标跟踪方法
CN111027677B (zh) 基于深度确定性策略梯度ddpg的多机动目标跟踪方法
CN108562894B (zh) 雷达波束指向与发射功率的分配方法
CN107942310B (zh) 分布式mimo雷达系统多目标位置估计的资源联合优化方法
Narykov et al. Algorithm for resource management of multiple phased array radars for target tracking
CN110376580B (zh) 一种异步多目标跟踪性能驱动的异构雷达网资源分配方法
Sun et al. JPBA of ARN for target tracking in clutter
CN107576945B (zh) 基于预测贝叶斯克拉美罗界的相控阵雷达重返和驻留时间分配方法
Schöpe et al. Multi-task sensor resource balancing using lagrangian relaxation and policy rollout
Yuan et al. Joint tracking sequence and dwell time allocation for multi-target tracking with phased array radar
CN108830007A (zh) 一种基于二次代价函数的雷达发射功率分配方法
Lu et al. Adaptive beam burden management for multi-target tracking in distributed radar network
Junkun et al. Joint power and bandwidth allocation for centeralized target tracking in multiple radar system
de Boer et al. Radar resource management for multi-target tracking using model predictive control
Zhang et al. Trajectory Planning for Airborne Radar in Extended Target Tracking Based on Deep Reinforcement Learning
Charlish et al. Information theoretic measures for MFR tracking control

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