CN116953688A - 一种跟踪与探测背景下阵元选择与功率分配联合优化方法 - Google Patents
一种跟踪与探测背景下阵元选择与功率分配联合优化方法 Download PDFInfo
- Publication number
- CN116953688A CN116953688A CN202310506787.8A CN202310506787A CN116953688A CN 116953688 A CN116953688 A CN 116953688A CN 202310506787 A CN202310506787 A CN 202310506787A CN 116953688 A CN116953688 A CN 116953688A
- Authority
- CN
- China
- Prior art keywords
- target
- tracking
- power
- model
- array element
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000001514 detection method Methods 0.000 title claims abstract description 46
- 238000005457 optimization Methods 0.000 title claims abstract description 42
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 54
- 239000013598 vector Substances 0.000 claims abstract description 27
- 239000002245 particle Substances 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims description 34
- 239000011159 matrix material Substances 0.000 claims description 32
- 238000005562 fading Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 description 25
- 238000013468 resource allocation Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 8
- 230000000739 chaotic effect Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 230000035772 mutation Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 230000004927 fusion Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 108090000623 proteins and genes Proteins 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
本发明属于雷达技术领域,具体地涉及一种跟踪与探测背景下阵元选择与功率分配联合优化方法。包括:步骤1:建立MIMO雷达预警空域内目标的回波信号模型;步骤2:建立MIMO雷达预警空域内目标的运动模型,在回波信号模型和运动模型的基础上,建立基于PCRLB的目标状态的BFIM;步骤3:根据回波信号模型,在两种似然假设条件下,建立分布式MIMO雷达的检测模型;步骤4:引入表示天线发射、接收阵元的选择与否,以及天线分配功率的矢量,并将矢量和步骤2的目标状态BFIM结合,建立跟踪与探测联合执行场景下的目标函数;步骤5:建立有限阵元和功率约束下的优化模型;步骤6:采用粒子群算法求解步骤5的优化模型的最优解,该最优解即为阵元与功率的最优分配。
Description
技术领域
本发明属于雷达技术领域,具体地涉及一种跟踪与探测背景下阵元选择与功率分配联合优化方法。
背景技术
相比于传统的单基地相控阵雷达和多基地组网雷达,分布式MIMO雷达具有优越的“四反”潜能,在未来战场中有着广阔的应用前景。理论上而言,分布式MIMO雷达包含的天线数量越多,发射功率越高,其检测、参数估计性能越佳。但在实际应用中,过多的子阵会带来较大的数据传输量,给信息融合中心造成计算负担;并且对于常见的车载、机载雷达而言,能量供应有限。因此,如何合理分配有限的资源,充分发挥MIMO雷达的潜能,是实际应用中的一个重要课题。另一方面,目前大多数研究均关注于某类任务(如目标定位、跟踪或成像)在执行过程中的资源分配,对于多元任务场景很少涉及。而在实际中,尤其是军事应用中,MIMO雷达需要同时执行搜索、跟踪、截获、制导、通信等多重任务。因此,研究多种任务背景下的资源分配,更加符合现实需求。
由于分布式MIMO雷达的检测性能与阵列配置和发射功率紧密相关,故在多元任务场景下,MIMO雷达的阵列选择和发射参数分配必定与跟踪、定位等单类任务场景中的资源分配有着诸多不同。
相比于传统的单基地相控阵雷达和多基地组网雷达,分布式MIMO雷达具有优越的“四反”潜能,在未来战场中有着广阔的应用前景。理论上而言,分布式MIMO雷达包含的天线数量越多,发射功率越高,其检测、参数估计性能越佳。但在实际应用中,过多的子阵会带来较大的数据传输量,给信息融合中心造成计算负担;并且对于常见的车载、机载雷达而言,能量供应有限。因此,如何合理分配有限的资源,充分发挥MIMO雷达的潜能,是实际应用中的一个重要课题。
目前多数文献均关注于某类任务(如目标定位、跟踪或成像)在执行过程中的资源分配,对于多元任务场景很少涉及。而在实际尤其是军事应用中,MIMO雷达需要同时执行搜索、跟踪、截获、制导、通信等多重任务。因此,研究多种任务背景下的资源分配,更加符合现实需求。另一方面,分布式MIMO雷达的检测性能与阵列配置和发射功率紧密相关。故在多元任务场景下,MIMO雷达的阵列选择和发射参数分配必定与跟踪、定位等单类任务场景中的资源分配有着诸多不同。
本文选取具有代表性的两类周期性任务——多目标跟踪和潜在威胁目标探测,以此为背景,对分布式MIMO雷达中的阵元选择与功率分配问题进行研究。
发明内容
针对上述存在的阵元与功率联合优化问题,本发明提出一种改进的PSO算法,采用混沌序列初始化种群,使得粒子群体获得良好的初始解集;根据粒子个体对约束条件的违反程度设计层级惩罚函数,以增强种群的局部搜索能力;对种群最优适应值的更新代数设置门限,当超过门限,最优值还未更新,则引入交叉、变异算子,以增加种群的多样性,促进PSO的全局寻优进程。同时,利用目标函数反比于SNR的特点,对SNR贡献大的天线分配更多的功率,以实现多任务场景下阵元与功率的同时分配。
为了实现上述目的,本发明所采用的技术方案如下:
一种跟踪与探测背景下阵元选择与功率分配联合优化方法,包括:
步骤1:建立MIMO雷达预警空域内目标的回波信号模型;
步骤2:建立MIMO雷达预警空域内目标的跟踪模型,跟踪模型包括运动模型和量测模型,在此基础上,建立基于PCRLB的目标状态的BFIM;
步骤3:根据步骤1的回波信号模型,在两种似然假设条件下,建立分布式MIMO雷达的探测模型;
步骤4:引入表示天线发射、接收阵元的选择与否,以及天线分配功率的矢量,并将矢量和步骤2的目标状态BFIM结合,建立跟踪与探测联合执行场景下的目标函数;
步骤5:基于步骤4,建立有限阵元和功率约束下的优化模型;
步骤6:采用粒子群算法求解步骤5的优化模型的最优解,该最优解即为阵元与功率的最优分配。
优选的,所述步骤4包括:
步骤4.1:引入如下矢量:
式中,tsm∈{0,1},表示第m个发射阵元的选择与否,rsn∈{0,1},表示第n个接收阵元的选择与否,Ptm表示第m个天线分配所得的功率;
步骤4.2:结合步骤4.1的矢量,建立目标状态的BFIM矩阵
F表示状态转移矩阵,Q表示过程噪声矩阵的协方差矩阵,表示目标状态向量,为关于目标状态的贝叶斯Fisher信息矩阵,/>和/>分别表示雅克比矩阵和关于/>的FIM在目标状态预测点/>处的近似值;
步骤4.3:在步骤4.2的基础上,MIMO雷达阵元选择与功率分配的BFIM表示为:
其中,为先验信息的FIM;
步骤4.4:在步骤4.3的基础上,选择最差的PCRLB为目标跟踪的优化标准,则跟踪与探测背景下的双目标函数构建为:
式中,tr(.)为矩阵求迹运算,Pdk为第k个采样间隔检测到新目标的概率;
步骤4.5:引入加权系数λ1,平衡步骤4.4中目标函数的位置和速度的量纲,即则将步骤4.4中的双目标函数转化为单目标优化函数:
式中,λ1为加权系数,且λ1∈(0,+∞)。
优选的,所述步骤5包括:
在限阵元和功率的约束下最优化雷达的跟踪与搜索性能,可选取的阵元总量受到限制:
K为可选取的阵元总数,且发射、接收阵元的选择数量最少为1:
雷达的功率资源是一定的,则:
式中,Ptot为雷达的总功率,为1×M维的全1矢量;
用虚警概率来表示检测性能:
Pfa=Pfa0
则:优化模型为:
tsm∈{0,1},rsn∈{0,1}
Pfa=Pfa0。
优选的,所述步骤6的功率分配公式为:
式中,m′属于已选择的发射阵元;q′代表对于所有选取的发射阵元路径衰落条件最差的目标;为第j个检测点到第m′个阵元的距离,/>为目标反射的复增益;/>为第m′个发射天线到第q′个目标的距离,λ1为加权系数;/>为第j个检测点到第m′个阵元的距离。
与现有技术相比,本发明的有益效果是:
(1)对于阵元选择和功率分配问题,所提出的MPSO优化算法具备良好的性能和较高的求解效率。并且,面对经过权重调整的目标函数以及时变的RCS模型时,MPSO算法展现了较强的鲁棒性和自适应性;
(2)对于已有目标跟踪与潜在威胁探测而言,具有较强路径衰减特性的子阵通常不被选择,算法一般会选取具有较好传播条件的子阵;
(3)在所提模型框架下,距跟踪误差最大的目标和检测中心最近的阵元将被分配给较多的功率资源,因为其对于SNR的贡献最大。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。
在附图中:
图1为本发明交叉操作图示;
图2为交叉操作图示;
图3为目标与MIMO雷达的几何关系;
图4为λ1=10时算法性能对比:(a)跟踪性能对比;(b)检测性能对比;
图5为算法运行时间对比;
图6为λ1=3时算法性能对比:(a)跟踪性能对比;(b)检测性能对比;
图7为λ1=20时算法性能对比:(a)跟踪性能对比;(b)检测性能对比;
图8为无RCS起伏条件下的资源分配结果;
图9为第二种目标RCS模型;
图10为第二种RCS模型条件下的资源分配结果:(a)功率分配结果,(b)接收阵元分配结果;
图11为本发明的方法流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
实施例:
参照附图1-10所示,一种跟踪与探测背景下阵元选择与功率分配联合优化方法,包括:
步骤1:建立MIMO雷达预警空域内目标的回波信号模型。包括:
假定一个分布式MIMO雷达系统有M个发射天线和N个接收天线,发射天线和接收天线的位置分别为(xtm,ytm)和(xrn,yrn),其中,m∈{1,...,M};n∈{1,...,N},Q个匀速运动的目标位于雷达警戒区域,在第k个采样间隔,第q个目标的位置和速度分别表示为和则第m个发射天线发射归一化的正交信号为:
其中,τ表示时延,正交信号的有效带宽和有效时宽可以分别表示为:
βm=(∫f2|Sm(f)|2df)/(∫|Sm(f)|2df) (3)
式中,f表示信号频率,Sm(f)为信号的傅氏变换,基于上述假设,第q个目标经由第(m,n)条路径的回波信号可表示为:
式中,为路径衰落影响;/>为目标反射的复增益,服从于复高斯分布和/>分别为时延和多普勒频移。/>为复高斯白噪声,服从可以具体表示为:
式中,为第m个发射天线到第q个目标的距离;/>为该目标到第n个接收天线的距离;Ptm为第m个天线的发射功率;τc为单次照射时段内的有效脉冲宽度;Gt和Gr分别为发射和接收增益;Ip为接收机处理增益;λ为载波波长;Lc和Lr分别为散射衰落和接收机衰落。和/>之间的关系为:
式中,c为光速。多普勒频移项计算为:
式中,和/>分别为目标与第m个发射天线和第n个接收天线相对于x轴方向的角度,满足:
步骤2:建立MIMO雷达预警空域内目标的跟踪模型,跟踪模型包括运动模型和量测模型,在此基础上,建立基于PCRLB的目标状态的BFIM;包括:
步骤2.1:建立MIMO雷达预警空域内目标的运动模型:
假设MIMO雷达预警空域内存在Q个目标,且目标均在二维平面内做匀速运动,则第q个目标的运动方程为:
目标状态向量用表示;F和/>分别为状态转移矩阵和过程噪声矩阵;/>
F和Q的具体形式为:
式中,Ts为采样间隔;为Kronecker积;nf为噪声系数,决定噪声强度。
步骤2.2:建立量测模型:
通过合适的信号处理机制可以获得丰富的目标信息,目标的时延和多普勒频移可以在第n个接收天线通过最大似然估计得到,假定量测为线性条件下获得的的,则量测函数为:
式中,H为4×4的单位阵;为量测噪声,服从/>此外,/>是关于时延和多普勒频移联合估计的CRLB的函数。
步骤2.3:建立基于PCRLB的目标状态的BFIM:
假设为/>的无偏估计,则总有式(1)成立:
式中,表示对目标状态/>和观测值/>取均值操作,/>为目标状态/>的贝叶斯Fisher信息矩阵(BFIM):
式中,表示对矢量求二阶偏导,目标状态和观测值的/>能够分解为:
式中,为目标状态的PDF(条件概率密度函数);/>为目标状态关于观测的似然函数。通过对BFIM进行分块操作,采用一种迭代计算PCRLB的方法,避免对大规模矩阵的运算操作:
式中,和/>分别为先验信息和数据的FIM:
在中,
将式(10)的目标运动模型代入式(6),并去掉期望符号,得到:
结合矩阵求逆运算,可以表示为:
对于式(5),将接收到的所有信号用向量形式表示:
则观测值关于目标状态/>的条件概率密度函数(PDF)可写作:
对于将(23)代入(18)得:
式中,为观测向量/>关于目标状态/>的雅克比矩阵:
具体形式为:
式中,Aτm,n和分别为Aτ和/>中的第mn列矩阵块。此外,因为经由每一条路径的回波信号相互独立,/>的FIM(Fisher信息矩阵)可写作矩阵分块形式:
式中,diag(.)代表MN×MN维对角矩阵,从而,可以得到BFIM的最终表达式:
从推导过程可以看出,仅与目标的运动模型有关,而阵元选取与功率分配通过影响/>进而影响PCRLB。因此,可以通过改变子阵和功率来获取最佳的PCRLB。但式(13)的第二项是一个求期望的过程,因此,为达到PCRLB,需要利用蒙特卡洛方法求取均值。为满足资源分配中的实时性要求,在此做如下近似:
式中,和/>分别表示雅克比矩阵和关于/>的FIM在目标状态预测点/>处的近似值。
步骤3:根据步骤1的回波信号模型,在两种似然假设条件下,建立分布式MIMO雷达的检测模型;具体为:
根据式(5)的信号模型,在H1和H2两种似然假设条件下,分布式MIMO雷达的检测模型为:
式中,aR和aI分别为目标反射复增益的实部和虚部。H1和H2分别代表目标存在与否。/>为复加性高斯白噪声,满足:
式中,(.)*为共轭运算,将经由全部路径所得的观测值和状态值写为向量的形式,则有:
y=[y1,1,...,yM,1,y1,2,...,yM,N]T (33)
x=[x1,1,...,xM,1,x1,2,...,xM,N]T (34)
在H1假设下,联合条件PDF为:
式中,且假定a服从0均值、方差为/>的高斯分布,IMN为MN×MN维的单位矩阵。由于似然比中常数项系数可最终并入检测门限,为简洁起见,利用常数项k0、k1和η0来代表计算过程中的无关系数。式(35)的推导过程中,用以下两式:
则似然比表示为:
将常数项移除,可以得到充分统计量:
注意的是,为复高斯随机变量,在H1和H2假设下,其协方差分别为/>和/>其中,h0是对所有hm,n求取均值:
即平均路径衰落。从而,可以得到充分统计量的分布:
式中,为自由度为2MN的中心化χ2分布。由此,可以得到虚警概率的表达式:
根据式(42)反解出检测门限η:
将式(43)代入检测概率的表达式,可得:
步骤4:引入表示天线发射、接收阵元的选择与否,以及天线分配功率的矢量,并将矢量和步骤2的目标状态BFIM结合,建立跟踪与探测联合执行场景下的目标函数;具体为:
步骤4.1:为描述天线选择与功率分配问题,首先引入如下矢量:
式中,,tsm∈{0,1},表示第m个发射阵元的选择与否,rsn∈{0,1},表示第n个接收阵元的选择与否,Ptm表示第m个天线分配所得的功率;
在多目标跟踪与探测联合执行场景下,分布式MIMO雷达需要实时改变天线配置和功率分配,以发挥最优性能。对于目标跟踪,PCRLB经常被用作优化标准。由于其不仅为状态估计提供了一个下界,更具有可预测性。因此,步骤2推导出的PCRLB是一个合适的选择。
步骤4.2:根据式(21),结合步骤4.1的矢量,为所有路径和的形式,三个矢量通过影响/>来影响PCRLB,进而建立目标状态的BFIM矩阵/>
F表示状态转移矩阵,Q表示过程噪声矩阵的协方差矩阵,表示目标状态向量,/>为关于目标状态的贝叶斯Fisher信息矩阵,/>和/>分别表示雅克比矩阵和关于/>的FIM在目标状态预测点/>处的近似值;
步骤4.3:在步骤4.2的基础上,结合公式(22),则MIMO雷达阵元选择与功率分配的BFIM表示为:
其中,为先验信息的FIM;
步骤4.4:在多目标跟踪过程中,资源应优先分配给较为紧迫的任务。在此,选择最差的PCRLB为目标跟踪的优化标准,则跟踪与探测背景下的双目标函数构建为:
式中,tr(.)为矩阵求迹运算,由于位置和速度的PCRLB量纲不同,直接相加并无意义,在此,提取位置的PCRLB作为目标函数。Pdk为第k个采样间隔检测到新目标的概率;阵元和功率是雷达优化参数。通过优化阵元和功率,提高目标的跟踪和探测精度。公式中的和1/Pdk(ts,rs,Pt)分别代表了雷达的跟踪和探测性能。
步骤4.5:从式(48)可以看出,阵元选择与功率分配为双目标优化问题,其帕累托前沿存在一系列解集。在求解过程中,着重考虑以下两点:(1)若罗列出任意时刻的所有解集,则需要求解每个解集对应的下一时刻的优化模型。当观测时间周期较长时,待求解的模型数量将过于庞大;
(2)PCRLB和检测概率的大小和单位均不相同,直接相加无任何意义。
因此,引入加权系数λ1,平衡步骤4.4中目标函数的位置和速度的量纲,即则将步骤4.4中的双目标函数转化为单目标优化函数:
式中,λ1为加权系数,且λ1∈(0,+∞)。其不仅平衡了跟踪与检测性能,而且统一了二者的意义。
步骤5:建立有限阵元和功率约束下的优化模型;包括:
对于多目标跟踪与探测联合执行中的阵元选择与功率分配问题,可以描述为如何在有限阵元和功率的约束下最优化雷达的跟踪与搜索性能。首先,减小融合中心的计算负担,可选取的阵元总量受到限制:
K为可选取的阵元总数,并且,为保证系统的有效性,发射、接收阵元的选择数量最少为1:
其次,雷达的功率资源是一定的,则:
式中,Ptot为雷达的总功率,为1×M维的全1矢量,用一定条件下的虚警概率来表示检测性能:
Pfa=Pfa0 (53)
则:优化模型为:
步骤6:采用粒子群算法求解步骤5的优化模型的最优解,该最优解即为阵元与功率的最优分配。
步骤5的式(54)为二元变量与连续型变量的混合优化问题,也是典型的NP难题。尽管可以通过穷举搜索罗列出所有可行解,其计算量将随问题规模呈指数型增长。当阵元总数较多时,计算量相当可观。为解决该问题,本发明提出一种改进的PSO搜索算法,其具备以下特征:
(1)引入混沌序列来优化初值。混沌序列拥有良好的随机性和遍历性,通过混沌优化,可以获得高质量的种群初始值,加速PSO的寻优进程;
(2)根据粒子个体对约束条件的违反程度设计层级惩罚函数,从而增强PSO的搜索精度;
(3)利用预置门限来引入扰动算子:当群体最优值达到门限后还未更新,则引入交叉、变异操作,提高种群多样性,促使种群跳出局部最优,实现全局最优。
整个改进PSO算法框架如图1所示,算法步骤如表1所示。
表1改进PSO算法步骤
PSO优化算法
在搜索过程中,每个粒子对自身寻得的历史最优解集pbest和种群寻得的最优解集gbest进行跟踪,通过不断迭代,最终收敛到全局最优解。在第t次迭代中,第i个粒子的速度和位置更新公式分别为:
/>
式中,和/>分别代表粒子个体的速度和位置矢量;w为惯性权重,其决定了粒子的局部搜索和全局搜索能力;c1和c2为认知参数,其决定了粒子飞向pbest和gbest的步长,通常取c1=c2=2。r1和r2是0~1之间服从正态分布的随机数。PSO算法本用来解决连续变量的优化问题,对于子阵选择问题而言,可以通过对小数取整操作将解空间划分为二元{0,1}区间。
除了优势之外,PSO算法中粒子群体的多样性将随着迭代次数快速下降,使得PSO算法易陷入局部极值。因此,提出混沌初始化优化机制,来增强PSO的全局搜索能力。
采用Logistic映射来产生混沌序列:
η(p+1)=μη(p)[1-η(p)] (57)
式中,η∈[0,1];μ为混沌控制参数。通过pmax次迭代,当μ=4且时,序列将呈现出良好的混沌特性。此时,η的轨迹将规则的分布于整个解空间内。
层级处罚函数
根据个体对约束的违反程度设计不同的惩罚因子,以促进PSO算法的局部搜索能力:
式中,K′为预置值。通过对个体施加层级处罚函数,可以使得位于可行解空间边界的劣解个体得以保留,并通过后续的交叉、变异操作转化为可行解,进而扩展PSO算法的搜索范围,促进搜索进程。
交叉、变异操作
首先,对种群中的每个个体产生随机数pc,且pc∈(0,1)。选择pc≤Pc的个体,作为父代。然后,对于所有选取的个体,随机挑出两个不相同的父代i和j,执行单点交叉操作。产生交叉点cp,且cp∈{2,…,M+N-1},交换父代i和j位于交叉点cp后的所有基因。最后,得到两个子代。交叉操作图示见图2,其中,gk为基因,且k∈{1,…,M+N}。
不同于交叉操作,变异操作基于单个个体。类似于交叉操作中的个体选择,首先随机选取个个体作为父代,对于每个个体,利用下式进行基因更新:
式中,为向下取整运算;⊙为Hardmard积;randn为服从标准正态分布的伪随机序列。
在子阵与功率联合分配过程中,需要注意以下两点:
(1)计算量问题。子阵与功率联合分配算法的计算量将是二者单独优化方法的乘积。在已有的单独优化方案中,各类算法均需要多次迭代才能达到收敛;
(2)耦合性问题。子阵与功率的联合优化是混合变量优化问题,且二者存在耦合,即单独的子阵选择最优与单独的功率分配最优难以保证联合最优。
针对以上两点,提出一种基于贪婪思想的功率分配方案:若发射天线对系统SNR的贡献较高,则对其分配更多的发射功率。这是基于目标函数中BFIM和检测概率均反比于SNR。则第m′个天线获得的功率分配比例可计算为:
式中,m′属于已选择的发射阵元;q′代表对于所有选取的发射阵元路径衰落条件最差的目标;为第j个检测点到第m′个阵元的距离,/>为目标反射的复增益;/>为第m′个发射天线到第q′个目标的距离,λ1为加权系数;/>为第j个检测点到第m′个阵元的距离。
仿真校验:
1.参数设置
选取M=N=10的分布式MIMO雷达系统进行分析。雷达警戒空域内包含Q=4批匀速运动目标,目标与MIMO雷达的几何关系如图3所示,在目标跟踪过程中,MIMO雷达将按照发射子阵1~接收子阵10的方向对警戒区域进行扫描。MIMO雷达的参数配置为:发射功率Ptot=10kW;单次照射有效脉冲宽度τc=0.2ms;载频fc=1GHz;发射、接收增益Gt=Gr=20dB;散射衰落和接收机衰落Lc=Lr=34.7dB;接收机处理增益Ip=20dB;有效带宽βm=1MHz;有效时宽Tm=1ms;采样间隔Ts=5s;虚警概率Pfa0=10-6。各目标的初始状态为:
10帧数据用于每次仿真。最优化模型中,系数λ1=10;选择的子阵数量限制为K=6。改进PSO优化算法中,种群基数Npop=50;惯性权重w=0.8;最大迭代次数tmax=50;混沌优化迭代次数pmax=3000;预置参数K′=3;imax=10;交叉概率Pc=0.6;变异概率Pm=0.3。
仿真中,引入以下基准算法进行对比。(1)穷举阵元选择与平均功率分配算法(exhaustive search selection with uniform power allocation algorithm,ES-UPA);(2)多起点贪婪阵元选择与平均功率分配算法(Multi-start local search selectionwith uniform power allocation algorithm,MSLS-UPA);(3)随机阵元选择与最优功率分配算法(Random selection with optimal power allocation algorithm,RS-OPA)。仿真基于MATLAB2014a软件,电脑配置为3.7GHz处理器与8GB运存。结果取100次蒙特卡洛实验的平均值。
2.仿真结果及分析:
2.1性能对比
对所提算法与三种基准算法的性能进行对比。假定所有目标对于发射阵元和接收阵元的反射系数均为1。
图4(a)和图4(b)分别为四种算法的跟踪和检测性能对比。从中可以看出,本发明所提出的改进PSO算法(modified PSO,MPSO)的性能与ES-UPA性能相近。并且,相比于其它两种算法,优势明显。这主要是因为在MPSO中运用了多种优化机制,如混沌初始化,层级处罚函数以及交叉、变异操作等。因而,MPSO获得了近似最优的性能。与之对比的是,MSLS算法在单次搜索过程中仅对本次的最优阵元进行选择,忽略了其它阵元对于未来搜索过程的影响。因此,极易陷入局部极值点。
2.2计算运行时间对比
为显示本文算法的高效性,我们对ES-UPA、MSLS-UPA和MPSO的运行时间进行对比,结果如图5所示。
从中可以看出,ES-UPA在阵元选取与功率分配过程中有着过高的处理器运行时间,已基本失去意义。MSLS-UPA性能适中,所提的MPSO算法结构轻巧,运行时间最短,更便于实用。
2.3MPSO算法的鲁棒性
对MPSO算法的鲁棒特性进行分析,图6和图7给出了λ1调整情况下跟踪与检测性能的对比。显然,MPSO算法依旧优于MSLS-UPA和RS-OPA,并且与ES-UPA性能相近,展示了较强地鲁棒特性。
2.4目标反射系数的影响
探究目标反射系数对于资源分配结果的影响。图8为所提MPSO算法在无RCS起伏条件下的阵元与功率分配结果。图中,每帧中的深蓝色区域代表tsm=0或rsn=0,其它区域代表tsm=1或rsn=1。同时,图8(a)中的颜色深度也代表各发射阵元的功率分配结果,定义为:
相应的,图8(b)中的颜色深度代表各接收阵元的选择频率,定义为:
式中,NRn代表第n个接收阵元被选取的次数;Nsim为仿真次数。从仿真结果可以看出,MPSO通常选取最近的阵元来对目标进行跟踪和检测;并且,更多的功率资源被分配给距目标和检测点最近的发射子阵,这是因为这些阵元对SNR的贡献较大。
接下来考虑第二种RCS起伏模型。各目标对于发射阵元3的反射系数如图9所示,对于其它发射子阵和接收子阵,各目标的反射系数均为1。该模型显示,所有目标对3号发射阵元有着较低的RCS。
图10给出了第二种RCS模型下的资源分配结果。与图8不同,此时算法不再选择3号发射阵元,而是用其它具有较强反射系数的发射阵元进行代替。与之对应,更多的功率被分配给这些子阵,以利用其距离目标和检测点较近的优势。以上结果也显示了所提算的自适应性。
2.5MPSO算法的最优参数
多目标跟踪与探测背景下的阵元选择与功率分配属于动态性问题,即最优解集是时变的。而PSO算法属于问题依赖性算法,对于某个问题性能优越的PSO算法在面对其它问题时可能性能下降严重。因此,需要对MPSO算法的最优参数进行探究。通过反复试验,MPSO算法的最优参数集合为Npop=50~100;w=0.75~1;tmax=100;pmax≥2000;imax=10~15;
Pc=0.6~0.8;Pm=0.25~0.40。
结论:
(1)对于阵元选择和功率分配问题,所提出的MPSO优化算法具备良好的性能和较高的求解效率。并且,面对经过权重调整的目标函数以及时变的RCS模型时,MPSO算法展现了较强的鲁棒性和自适应性;
(2)对于已有目标跟踪与潜在威胁探测而言,具有较强路径衰减特性的子阵通常不被选择,算法一般会选取具有较好传播条件的子阵;
(3)在所提模型框架下,距跟踪误差最大的目标和检测中心最近的阵元将被分配给较多的功率资源,由于其对于SNR的贡献最大。
仿真验证了MPSO算法的优越性能。此外,仿真结果也表明,具有较好传播条件的子阵有更大的被选择概率,更多功率资源将分配给距离目标和检测点较近的发射阵元。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (4)
1.一种跟踪与探测背景下阵元选择与功率分配联合优化方法,其特征在于:包括:
步骤1:建立MIMO雷达预警空域内目标的回波信号模型;
步骤2:建立MIMO雷达预警空域内目标的跟踪模型,跟踪模型包括运动模型和量测模型,在此基础上,建立基于PCRLB的目标状态的BFIM;
步骤3:根据步骤1的回波信号模型,在两种似然假设条件下,建立分布式MIMO雷达的探测模型;
步骤4:引入表示天线发射、接收阵元的选择与否,以及天线分配功率的矢量,并将矢量和步骤2的目标状态BFIM结合,建立跟踪与探测联合执行场景下的目标函数;
步骤5:基于步骤4,建立有限阵元和功率约束下的优化模型;
步骤6:采用粒子群算法求解步骤5的优化模型的最优解,该最优解即为阵元与功率的最优分配。
2.根据权利要求1所述的一种跟踪与探测背景下阵元选择与功率分配联合优化方法,其特征在于:所述步骤4包括:
步骤4.1:引入如下矢量:
式中,tsm∈{0,1},表示第m个发射阵元的选择与否,rsn∈{0,1},表示第n个接收阵元的选择与否,Ptm表示第m个天线分配所得的功率;
步骤4.2:结合步骤4.1的矢量,建立目标状态的BFIM矩阵
F表示状态转移矩阵,Q表示过程噪声矩阵的协方差矩阵,表示目标状态向量,/>为关于目标状态的贝叶斯Fisher信息矩阵,/>和/>分别表示雅克比矩阵和关于/>的FIM在目标状态预测点/>处的近似值;
步骤4.3:在步骤4.2的基础上,MIMO雷达阵元选择与功率分配的BFIM表示为:
其中,为先验信息的FIM;
步骤4.4:在步骤4.3的基础上,选择最差的PCRLB为目标跟踪的优化标准,则跟踪与探测背景下的双目标函数构建为:
式中,tr(.)为矩阵求迹运算,Pdk为第k个采样间隔检测到新目标的概率;
步骤4.5:引入加权系数λ1,平衡步骤4.4中目标函数的位置和速度的量纲,即则将步骤4.4中的双目标函数转化为单目标优化函数:
式中,λ1为加权系数,且λ1∈(0,+∞)。
3.根据权利要求2所述的一种跟踪与探测背景下阵元选择与功率分配联合优化方法,其特征在于:所述步骤5包括:
在限阵元和功率的约束下最优化雷达的跟踪与搜索性能,可选取的阵元总量受到限制:
K为可选取的阵元总数,且发射、接收阵元的选择数量最少为1:
雷达的功率资源是一定的,则:
式中,Ptot为雷达的总功率,为1×M维的全1矢量;
用虚警概率来表示检测性能:
Pfa=Pfa0
则:优化模型为:
tsm∈{0,1},rsn∈{0,1}
Pfa=Pfa0。
4.根据权利要求3所述的一种跟踪与探测背景下阵元选择与功率分配联合优化方法,其特征在于:所述步骤6的功率分配公式为:
式中,m′属于已选择的发射阵元;q′代表对于所有选取的发射阵元路径衰落条件最差的目标;为第j个检测点到第m′个阵元的距离,/>为目标反射的复增益;/>为第m′个发射天线到第q′个目标的距离,λ1为加权系数;/>为第j个检测点到第m′个阵元的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310506787.8A CN116953688A (zh) | 2023-05-08 | 2023-05-08 | 一种跟踪与探测背景下阵元选择与功率分配联合优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310506787.8A CN116953688A (zh) | 2023-05-08 | 2023-05-08 | 一种跟踪与探测背景下阵元选择与功率分配联合优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116953688A true CN116953688A (zh) | 2023-10-27 |
Family
ID=88450134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310506787.8A Pending CN116953688A (zh) | 2023-05-08 | 2023-05-08 | 一种跟踪与探测背景下阵元选择与功率分配联合优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116953688A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117706541A (zh) * | 2024-02-06 | 2024-03-15 | 四川省华盾防务科技股份有限公司 | 相控阵多目标跟踪方法、装置、设备及存储介质 |
-
2023
- 2023-05-08 CN CN202310506787.8A patent/CN116953688A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117706541A (zh) * | 2024-02-06 | 2024-03-15 | 四川省华盾防务科技股份有限公司 | 相控阵多目标跟踪方法、装置、设备及存储介质 |
CN117706541B (zh) * | 2024-02-06 | 2024-04-16 | 四川省华盾防务科技股份有限公司 | 相控阵多目标跟踪方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Power and bandwidth allocation for multi-target tracking in collocated MIMO radar | |
CN109283522B (zh) | 一种联合时空资源管理的共址mimo雷达目标跟踪方法 | |
Zhang et al. | Sensor scheduling and resource allocation in distributed MIMO radar for joint target tracking and detection | |
CN111025275B (zh) | 基于射频隐身的多基雷达辐射参数多目标联合优化方法 | |
CN111323773A (zh) | 基于射频隐身的组网雷达功率与带宽联合优化分配方法 | |
CN110412515B (zh) | 基于射频隐身的组网雷达多目标跟踪发射功率分配方法 | |
CN111190176B (zh) | 共址mimo雷达组网系统的自适应资源管理方法 | |
CN108896985B (zh) | 基于射频隐身的组网雷达多目标跟踪采样间隔控制方法 | |
CN110673131B (zh) | 多波束集中式mimo雷达时空资源-波形选择管理方法 | |
WO2021258734A1 (zh) | 博弈条件下基于低截获性能的组网雷达最优波形设计方法 | |
CN116953688A (zh) | 一种跟踪与探测背景下阵元选择与功率分配联合优化方法 | |
CN115567353B (zh) | 一种面向雷达组网系统的干扰多波束调度和干扰功率联合优化方法 | |
CN110412534A (zh) | 基于射频隐身的组网雷达多目标跟踪驻留时间优化方法 | |
Zhang et al. | Joint resource optimization for a distributed MIMO radar when tracking multiple targets in the presence of deception jamming | |
CN113671487B (zh) | 一种基于混合相控阵-mimo雷达的目标搜索资源优化方法 | |
KR102665123B1 (ko) | 저피탐지를 위한 통합망 레이더 체류 시간과 방사 전력의 공동 최적화 방법 | |
CN116381641A (zh) | 基于lpi的cmimo目标跟踪的自适应资源管理方法 | |
CN115963451A (zh) | 基于极化特性自适应阵元选择的同时多波束雷达的波束驻留调度方法 | |
CN114325685B (zh) | 一种分布式mimo雷达多目标跟踪中的波束空-时分配算法 | |
Li et al. | Resource optimization strategy in phased array radar network for multiple target tracking when against active oppressive interference | |
Du et al. | Joint power and bandwidth allocation algorithm for radio frequency stealth radar | |
Reddy et al. | An Excavated PSO-SQP Mechanisms for an Effective Radar Transmission in MIMO-OFDM Communication Systems. | |
Wang et al. | LPI time‐based TMS against high‐sensitivity ESM | |
CN114325685A (zh) | 一种分布式mimo雷达多目标跟踪中的波束空-时分配算法 | |
Hoseinzade et al. | Power control and beamforming of phased-MIMO radars to achieve a low probability interception |
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 |