作者中心
专家审稿
责编办公
编辑办公
Analysis and Experimental Validation of Key Technologies for Unmanned Aerial Vehicle-borne Bistatic Interferometric Synthetic Aperture Radar
-
摘要: 双站干涉合成孔径雷达(SAR)技术可突破单站双天线干涉SAR的基线限制，获得更加灵活的基线构型，有利于提升干涉测量精度。无人机载双站干涉SAR机动性好、飞行成本低，具有很高的应用价值，也能为无人机载分布式干涉、三维成像等提供关键支撑，具有重要的研究意义。中国科学院空天信息创新研究院牵头设计研制了国内首套无人机载双站干涉SAR系统，并在内蒙古百灵机场开展了飞行实验。该文简要介绍了该系统方案设计、基本构成和主要性能，并重点介绍了该系统的空间、时间及相位同步和数据处理关键技术，然后介绍了首次飞行实验的方案和实施研究情况，最后给出了实验数据处理结果，验证了关键技术和该无人机载双站干涉SAR系统0.5 m高程测量精度等主要指标，为后续多无人机平台协同开展分布式干涉数据获取及处理研究提供了基础。Abstract: Bistatic interferometric Synthetic Aperture Radar (SAR) overcomes the baseline length limit of the configuration of single-station interferometric SAR with two antennas and has become the primary method of terrain mapping using spaceborne interferometric SAR. To reduce the cost of surveying and mapping while promoting the development and application of Unmanned Aerial Vehicle (UAV)-borne bistatic interferometric SAR, the Aerospace Information Research Institute, Chinese Academy of Sciences took the lead in designing and developing a UAV-borne bistatic interferometric SAR processing system and performed flight experiments at Bailing Airport in Inner Mongolia. Herein, the system design, composition, and performance are introduced, and the scheme and implementation of the first flight experiment, along with the preliminary data processing results, are presented. In addition, the key performance metrics of the system, such as 0.5 m-elevation measurement accuracy, are verified in this study. The system serves as a foundation for future research topics, such as distributed InSAR using a multiaviation platform and tomography data acquisition and processing.
-
1. 引言
当今及未来战场环境已经呈现出如下特点：(1)作战要素、作战样式复杂化；(2)战场资源种类繁多且属性各异；(3)资源配置具有一定时间紧迫性。传统的单基地雷达受到平台局限性的影响，难以完成既定目标探测任务，也无法满足作战要求[1,2]。因此，异构多雷达网络逐渐崭露头角。异构多雷达网络是指由具有不同工作体制的若干种类型雷达网络构成的有机整体[3]，在空间复用、波形分集以及抗干扰能力等方面都具有显著优势。异构多雷达网络系统通过收集不同雷达网络多个雷达节点获取的量测信息来实现远距离多目标的探测，具有较强的分布式并行感知计算能力、优越的鲁棒性和较低的任务容错率，在不同领域都具有广阔的应用前景[4]。
随着信号处理技术的迅速发展，面向多目标跟踪的雷达辐射资源管理研究成为主流，如何高效分配雷达系统的有限辐射资源以提升多目标协同跟踪性能，已经成为国内外学者聚焦的研究课题[5-9]。2019年，Yan等人[10]基于Pareto理论提出了一种相控阵雷达驻留时间优化算法，该算法以提高多目标跟踪精度以及目标信噪比作为双优化目标，并采用线性规划方法结合Min-Max算法进行求解，仿真结果表明，该算法能够同时提升多目标跟踪精度以及目标信噪比。Zhang等人[11]针对分布式多输入多输出(Multiple-Input Multiple-Output, MIMO)雷达多目标检测与跟踪问题，以最坏情况下的目标后验克拉美罗下界(Posterior Cramér-Rao Lower Bound, PCRLB)以及检测概率作为代价函数，采用改进粒子群算法自适应优化设计系统发射机-接收机选择方式、辐射功率和信号带宽的分配，仿真结果表明，该算法具有更优越的实时性以及更好的多目标跟踪性能。2020年，靳标等人[12]基于合作博弈思想，研究了基于多目标跟踪的去中心化组网雷达功率分配问题，仿真结果表明，相较于功率均匀分配算法，所提算法可以有效提升多目标跟踪性能；相较于基于非合作博弈的功率分配算法，所提算法能够分配给对提升目标跟踪性能作用更显著的雷达节点更多的功率，提升了资源利用效率。文献[13]提出了一种大规模MIMO雷达网络认知目标跟踪功率和阵列选择联合优化算法，推导了目标跟踪的预测条件克拉美罗下界(Predicted Conditional-Cramér-Rao Lower Bound, PC-CRLB)并将其作为目标函数进行优化，并提出了基于局部搜索的求解算法对功率和阵列选择进行设计。Lu等人[14]针对频谱共存下组网相控阵雷达系统多目标跟踪场景，提出了一种基于目标跟踪性能的动态资源调度算法，根据通信区域干扰模型推导出目标的PCRLB并将其作为多目标跟踪性能指标，实现了雷达与通信频谱共存下的多目标协同跟踪性能优化。
2021年，Sun等人[15]在传感器位置不确定情况下，针对组网雷达多目标跟踪场景以及出现的雷达量测信息冗余问题，提出了一种量测信息选择与功率联合优化分配算法，该算法可以同时优化辐射功率和量测信息选择传播路径，以提升多目标跟踪精度和降低数据处理负担，采用稀疏性增强的序列凸规划算法进行求解并通过仿真实验进行了有效验证。同年，他们又研究了杂波环境下组网集中式MIMO雷达系统多目标跟踪问题，提出了一种雷达选择和波束图设计联合优化分配算法，旨在优化每部雷达的调度方法和波形相关矩阵来最大限度地提高目标跟踪精度[16]。2022年，文献[17]针对欺骗干扰下的多目标跟踪场景，研究了分布式MIMO雷达的联合资源优化问题，通过基于三步优化循环的算法对系统发射机、接收机联合选择以及辐射功率进行自适应优化设计，仿真结果表明，离干扰机较近且归一化欺骗参数较低的雷达发射机会被分配更多的功率资源，从而有效提升多目标跟踪性能。Su等人[18]提出了一种机动目标跟踪场景下的组网集中式MIMO雷达波形控制与资源联合优化算法，以最小化机动目标跟踪误差作为优化目标，采用改进粒子群算法进行求解，达到了提升机动目标跟踪精度的目的。
雷达的系统性能不仅受自身辐射参数的影响，还与自身系统构型紧密相关，相较于同构雷达，异构雷达具备更灵活的资源分配方式和更卓越的目标探测能力。此外，在实际中，各个雷达节点的时间基准、初始采样时间与采样间隔可能是不同的，因此，雷达系统得到的目标量测信息是异步的。围绕异构雷达系统资源管理以及异步目标跟踪问题，国内外学者进行了大量研究[19,20]。2013年，文献[21]针对异步目标跟踪问题，提出了一种多基地雷达系统功率分配算法，根据量测信息的最优顺序融合推导了目标跟踪精度计算表达式并作为优化目标函数，以辐射功率上下限以及总辐射功率限制为约束条件，建立了多基地雷达系统的异步观测模型，结合凸松弛方法以及梯度投影法进行求解。2019年，文献[22]研究了一种杂波环境下异步雷达网络目标跟踪的功率带宽联合分配算法，在给定系统辐射功率和发射带宽限制的约束条件下，采用分支归约界定算法来进行求解，仿真结果表明，相较于均匀功率带宽分配对比算法，该算法可以明显提升目标跟踪精度。2020年，文献[23]将异构雷达网络与异步多目标跟踪相结合，提出了两种最优资源分配方法，并结合梯度投影法和双上升法对每个融合采样间隔内异构雷达网络的辐射功率和驻留时间进行优化设计，在不同情况下分别有效提升了系统的多目标跟踪精度和射频隐身性能。
然而，上述研究均是在理想检测环境下开展的。实际上，目标雷达散射截面(Radar Cross Section, RCS)起伏以及信号辐射衰减等不确定因素会导致雷达系统的检测概率小于1[24,25]。针对该问题，国内外学者开展研究，取得了一定进展。2020年，Sun等人[26]提出了一种非理想检测下的多基地雷达系统多目标跟踪功率分配算法，该算法推导了非理想检测下的目标状态估计误差PCRLB表达式，并将其作为目标函数，在满足节点选择和功率限制的约束条件下，采用了基于Zoutendijk可行方向法的求解算法对优化模型进行求解，仿真结果表明，该算法可以显著提升系统多目标跟踪精度。2022年，在文献[26]的基础上，文献[27]研究了检测概率变化下的多基地雷达系统波束和功率联合优化分配问题，分析了检测概率与目标RCS、波束分配方式以及辐射功率等参数的数学关系，推导了非理想检测下的目标PCRLB并基于此构建了总体代价函数，以给定的波束和功率资源限制为约束条件，建立了相应的数学优化模型，采用结合贪婪策略和Zoutendijk可行方向法的三步分解算法进行求解，仿真结果验证了所提算法对检测概率变化时多目标跟踪性能的提升。张永平[28]针对雷达系统多目标跟踪背景，提出了一种非理想检测环境下的资源管理算法，通过枚举计算方式得到了非理想检测下的目标PCRLB并作为优化目标函数，对发射波束和功率等参数进行优化分配，从而提升了多目标跟踪精度。
上述研究成果为提升雷达系统的目标跟踪性能奠定了坚实基础，然而，已有研究并未考虑在非理想检测环境下针对异构多雷达网络异步多目标跟踪场景进行资源联合优化。因此，本文围绕非理想检测、异步多目标跟踪以及异构多雷达网络辐射资源管理等领域，研究了非理想检测下面向异步多目标跟踪的异构多雷达网络功率时间联合优化问题。首先，针对异构多雷达网络，在非理想检测下，建立了面向异步多目标跟踪的功率时间联合优化模型，以最小化异步多目标跟踪误差为优化目标，在满足给定的系统辐射资源限制前提下，通过联合优化异构多雷达网络中各雷达节点的选择方式、辐射功率和驻留时间等发射参数，提升系统的异步多目标跟踪精度。仿真结果表明，与现有算法相比，所提算法能够有效降低异构多雷达网络的异步多目标跟踪误差，提升跟踪精度。
2. 系统模型
假设异构多雷达网络系统由M部不同体制和工作模式的雷达子网络组成，系统中第
$m\left( {m = 1,2, \cdots ,M} \right)$ 部雷达网络中包含${N_m}$ 个雷达节点，且第m部雷达网络中的第$n\left( {n = 1,2, \cdots ,{N_m}} \right)$ 个雷达节点的坐标为$ \left( {{x_{n,m}},{y_{n,m}}} \right) $ 。同时，由于不同雷达网络中的雷达节点工作模式不同，各雷达节点发射信号的载波频率可能不同，对于目标q，它们的初始采样时间$ t_{n,m,q}^{{\text{ini}}} $ 和采样间隔$ T_{n,m,q}^0 $ 也可能是不同的。为了后续小节中的进一步推导与计算以及更好地展现多雷达网络的异构特性，定义如下3种类型的雷达网络及发射参数约束：(1) 集中式MIMO雷达网络：在同时多波束工作模式下，雷达节点照射目标的驻留时间固定，但辐射功率会发生变化[29]；
(2) 相控阵雷达网络：该雷达网络中的雷达节点可以快速自适应地旋转波束以照射多个目标，基于此，该雷达网络中的雷达节点辐射功率固定，但驻留时间不同；
(3) 机械雷达扫描网络：该雷达网络中的雷达节点工作参数固定，即所有雷达节点以固定的辐射功率、驻留时间照射目标。
此外，假设异构多雷达网络系统的监视区域中存在Q个做匀速直线运动的独立目标，且目标
$q\left( {q = 1,2, \cdots ,Q} \right)$ 的初始位置为$ \left( {x_{{\text{ini}}}^q,y_{{\text{ini}}}^q} \right) $ ，运动速度为$ \left( {v_{x,{\text{ini}}}^q,v_{y,{\text{ini}}}^q} \right) $ 。2.1 目标运动模型
假设异构多雷达网络系统的融合采样间隔为
$ {T_{{\text{fusion}}}} $ ，且第k个融合采样间隔的区间为$ \left( {{t_{k - 1}},{t_k}} \right] $ ，则$ {T_{{\text{fusion}}}} = {t_k} - {t_{k - 1}} $ 。因此，在第k个融合采样间隔，根据目标初始位置和运动速度，目标q的状态向量可以表示为Xqtk=[xqtk,yqtk,vqx,tk,vqy,tk]T (1) 其中，
$ {\left\{ \cdot \right\}^{\text{T}}} $ 表示矩阵的转置运算，根据状态向量${\boldsymbol{X}}_{{t_k}}^q$ ，将第k个融合采样间隔时目标q的运动状态方程表示为Xqtk=FqXqtk−1+Uqtk−1 (2) 其中，
$ {{\boldsymbol{F}}^q} $ 表示目标q的状态转移矩阵。由于目标作匀速直线运动，可以将$ {{\boldsymbol{F}}^q} $ 计算为Fq=I2⊗[1Tfusion01]=[1Tfusion000100001Tfusion0001] (3) 其中，
$ {{\boldsymbol{I}}_2} $ 表示2阶单位矩阵，$ \otimes $ 表示矩阵的直积运算。式(2)中的$ {\boldsymbol{U}}_{{t_{k - 1}}}^q $ 表示零均值高斯过程白噪声，其协方差矩阵${\boldsymbol{Q}}_{k - 1}^q$ 可以写为Qqk−1=(γqI2)⊗[(Tfusion)33(Tfusion)22(Tfusion)22Tfusion] (4) 其中，
$ {\gamma ^q} $ 表示目标q的过程噪声强度。2.2 异步雷达量测模型
为了简化后续小节中的推导与计算，在本节中做出如下假设：(1)在每个融合采样间隔，异构多雷达网络中的每个雷达节点最多只能被分配去跟踪一个目标；(2)所有雷达节点只能接收来自自身发射信号的回波。
在异构多雷达网络进行异步多目标跟踪时，系统中各个雷达节点的初始采样时间和采样间隔可能不同。因此，需要考虑雷达量测信息的异步特性。假设第k个融合采样间隔内第m部雷达网络中的第n个雷达节点对目标q的量测信息数目为
$ S_{n,m,k}^q $ ，其中，第$s\left( {s = 1,2, \cdots ,S_{n,m,k}^q} \right)$ 个量测信息的到达时刻为$ t_{n,m,k}^{q,s} $ 。基于此，图1给出了目标q在第k个融合采样间隔内的异步量测模型。
图 1 目标q在第k个融合采样间隔内的异步量测模型Figure 1. Asynchronous measurement model of target q in the k-th fusion sampling interval在此基础上，定义一个二元变量来描述异构多雷达网络中各雷达节点与目标q的配对情况，即
μqn,m,k={1, k时刻第m部雷达网络中的第n个雷达节点对目标q进行照射0, 其他 (5) 在第k个融合采样间隔
$ \left( {{t_{k - 1}},{t_k}} \right] $ ，结合雷达节点-目标配对指标$ \mu _{n,m,k}^q $ ，可以得到第m部雷达网络中的第n个雷达节点对目标q第s个量测的方程：zq,sn,m,k={hqn,m,k(Xqtk)+ζq,sn,m,k, μqn,m,k=10, μqn,m,k=0 (6) 其中，
$ {\boldsymbol{z}}_{n,m,k}^{q,s} $ 表示第s个量测的量测矢量，$ {\boldsymbol{h}}_{n,m,k}^q\left( \cdot \right) $ 表示第s个量测的非线性观测函数。$ {\boldsymbol{h}}_{n,m,k}^q\left( {{\boldsymbol{X}}_{{t_k}}^q} \right) $ 的计算式可以写为hqn,m,k(Xqtk)=[Rq,sn,m,kϑq,sn,m,k]=[√(xqtq,sn,m,k−xn,m)2+(yqtq,sn,m,k−yn,m)2arctan(yqtq,sn,m,k−yn,mxqtq,sn,m,k−xn,m)] (7) 其中，
$ R_{n,m,k}^{q,s} $ 和$ \vartheta _{n,m,k}^{q,s} $ 分别表示第m部雷达网络中的第n个雷达节点在第s个量测时与目标q之间的实际距离和方位角。$ x_{t_{n,m,k}^{q,s}}^q $ 和$ y_{t_{n,m,k}^{q,s}}^q $ 表示在到达时刻$ t_{n,m,k}^{q,s} $ 时目标q的位置坐标。此外，式(6)中的$ {\boldsymbol{\zeta }}_{n,m,k}^{q,s} $ 表示第s个量测的误差，其协方差矩阵可以表示为Gq,sn,m,k=[σ2Rq,sn,m,k00σ2ϑq,sn,m,k] (8) 它们与雷达相关发射参数的关系可以表示为[10]
{σ2Rq,sn,m,k∝{Pq,sn,m,kTq,sn,m,k[βq,sk]2}−1σ2ϑq,sn,m,k∝[Pq,sn,m,kTq,sn,m,k/Pq,sn,m,kTq,sn,m,kBNNBNN]−1 (9) 其中，
$ P_{n,m,k}^{q,s} $ 表示第m部雷达网络中的第n个雷达节点在第s个量测时对目标q的辐射功率；$ T_{n,m,k}^{q,s} $ 表示第m部雷达网络中的第n个雷达节点在第s个量测时对目标q的驻留时间；$ \beta _k^{q,s} $ 表示发射信号带宽，$ {B_{{\text{NN}}}} $ 表示雷达接收天线的波束宽度。可以看出，雷达节点的辐射功率与驻留时间等参数与目标的量测误差有关。3. 异构多雷达网络功率时间联合优化算法
3.1 非理想检测下的异步目标BCRLB推导
由于贝叶斯克拉美罗下界(Bayesian Cramér-Rao Lower Bound, BCRLB)可以为参数无偏估计的均方误差(Mean Square Error, MSE)提供一个下界[30]，因此，大多数学者将其作为目标跟踪精度的衡量指标。在第k个融合采样间隔，将异构多雷达网络系统中不同雷达节点关于目标q的复合量测信息按顺序依次列为
Zqk=[zq,11,1,k,⋯,zq,Sqn,m,kn,m,k,⋯,zq,SqNm,M,kNM,M,k]T (10) 结合复合量测信息
$ {\boldsymbol{Z}}_k^q $ 进行后续对目标BCRLB的推导。根据文献[31]，在理想检测下，第k个融合采样间隔时目标q状态估计误差的BCRLB表达式可以写为CqBCRLB,k=[[Qqk−1+FqJ−1(Xqtk−1)(Fq)T]−1+M∑m=1Nm∑n=1μqn,m,k(Hq,sn,m,k)T(Gq,sn,m,k)−1Hq,sn,m,k]−1 (11) 其中，
$ {\boldsymbol{H}}_{n,m,k}^{q,s} $ 表示$ {\boldsymbol{h}}_{n,m,k}^q\left( \cdot \right) $ 的雅可比矩阵。然而，根据文献[27]，在非理想检测环境下，异构多雷达网络中的雷达节点并不都能对所有跟踪目标成功检测，还可能存在漏检。基于此，定义如下二元变量来表征第k个融合采样间隔时异构多雷达网络对目标q的检测情况：
dqn,m,k={1, k时刻第m部雷达网络中的第n个雷达节点检测到目标q 0, 其他 (12) 根据系统模型，可知异构多雷达网络中共存在
$\displaystyle\sum\nolimits_{m = 1}^M {{N_m}}$ 个雷达节点。因此，在每个融合采样间隔，共有${2^{\sum\nolimits_{m = 1}^M {{N_m}} }}$ 种检测情况，将这些情况具体描述为Dqk={Dqi,k|i=1,2,⋯,2M∑m=1Nm} (13) 其中，
${\boldsymbol{D}}_{i,k}^q$ 表示第k个融合采样间隔异构多雷达网络对目标q的第i种检测情况。假设在第i种检测情况中，多雷达网络共有$\psi _{{i,k}}^{q}$ 个雷达节点检测到目标q，则未检测到目标q的雷达节点数目为$\displaystyle\sum\nolimits_{m = {\text{1}}}^M {{N_m}} - \psi _{i,k}^q$ 。基于此，定义第i种检测情况出现的概率为Pr{Dqi,k}=(pqD,n,m,k)ψqi,k⋅(1−pqD,n,m,k)M∑m=1Nm−ψqi,k (14) 其中，
$p_{{\text{D}},n,m,k}^q$ 表示第k个融合采样间隔第m部雷达网络中的第n个雷达节点对目标q的检测概率。实际上，不同雷达节点对目标的检测概率可能是不一样的，它与目标RCS、雷达节点与目标的距离等参数有关。为简化本文后续推导与计算，$p_{{\text{D}},n,m,k}^q$ 在整个跟踪过程中被视为常数。综上所述，结合式(11)—式(14)，在非理想检测下，结合不同雷达节点对目标q的量测信息数目
$ S_{n,m,k}^q $ ，将第k个融合采样间隔目标q的异步跟踪误差BCRLB计算为[27]˜CqBCRLB,k=2M∑m−1Nm∑i=1J−1(Xqtk|Dqi,k)⋅Pr{Dqi,k} (15) 其中，
${\boldsymbol{J}}\left( {{\boldsymbol{X}}_{{t_k}}^q\left| {{\boldsymbol{D}}_{i,k}^q} \right.} \right)$ 表示目标状态的贝叶斯信息矩阵(Bayesian Information Matrix, BIM)，可以写为J(Xqtk|Dqi,k)=[Qqk−1+FqJ−1(Xqtk−1|Dqi,k−1)(Fq)T]−1+M∑m=1Nm∑n=1Sqn,m,k∑s=1μqn,m,kdqn,m,k(Hq,sn,m,k)T⋅(Gq,sn,m,k)−1Hq,sn,m,k (16) 从式(16)中可以看出，
$ {\boldsymbol{J}}\left( {\left. {{\boldsymbol{X}}_{{t_k}}^q} \right|{\boldsymbol{D}}_{i,k}^q} \right) $ 中的第1项，即先验信息的费舍尔信息矩阵(Fisher Information Matrix, FIM)与目标运动模型以及前一时刻目标q状态的BIM有关；$ {\boldsymbol{J}}\left( {\left. {{\boldsymbol{X}}_{{t_k}}^q} \right|{\boldsymbol{D}}_{i,k}^q} \right) $ 中的第2项，即量测信息的FIM与雷达节点与目标的配对指标$ \mu _{n,m,k}^q $ 、雷达节点对目标的检测指标$ d_{n,m,k}^q $ 、第k个融合采样间隔内第m部雷达网络中的第n个雷达节点对目标q的量测信息数目$ S_{n,m,k}^q $ 以及量测误差协方差矩阵$ {\boldsymbol{G}}_{n,m,k}^{q,s} $ 有关。因此，量测信息的FIM受到第k个融合采样间隔内雷达节点选择以及辐射功率和驻留时间优化结果的影响。${{\tilde {\boldsymbol{C}}}}_{{\text{BCRLB}},k}^q$ 矩阵对角线上的第1个元素${{\tilde {\boldsymbol{C}}}}_{{\text{BCRLB}},k}^q \left( {1,1} \right)$ 与第2个元素${{\tilde {\boldsymbol{C}}}}_{{\text{BCRLB}},k}^q\left( {2,2} \right)$ 表征目标预测估计MSE的距离下界。因此，将它们的和作为异步多目标跟踪精度的衡量指标：F(μqk,Pqt,k,Tqk)≜ (17) 其中，
$ {\boldsymbol{\mu }}_k^q $ ,$ {\boldsymbol{P}}_{{\text{t,}}k}^q $ 和$ {\boldsymbol{T}}_k^q $ 分别表示异构多雷达网络关于目标q的雷达节点选择、辐射功率以及驻留时间矩阵，且维度均为$\displaystyle\sum\nolimits_{m = 1}^M {{N_m}} \times 1$ ：\left\{ \begin{gathered} {\boldsymbol{\mu }}_k^q = {\left[ {\mu _{1,1,k}^q, \cdots ,\mu _{n,m,k}^q, \cdots ,\mu _{{N_M},M,k}^q} \right]^{\text{T}}} \\ {\boldsymbol{P}}_{{\text{t,}}k}^q = {\left[ {P_{1,1,k}^q, \cdots ,P_{n,m,k}^q, \cdots ,P_{{N_M},M,k}^q} \right]^{\text{T}}} \\ {\boldsymbol{T}}_k^q = {\left[ {T_{1,1,k}^q, \cdots ,T_{n,m,k}^q, \cdots ,T_{{N_M},M,k}^q} \right]^{\text{T}}} \end{gathered} \right. (18) 3.2 优化模型建立
本文提出了一种非理想检测下面向异步多目标跟踪的异构多雷达网络功率时间联合优化算法，通过在各个融合采样间隔联合优化雷达节点选择方式、辐射功率以及驻留时间，在满足给定系统射频辐射资源限制的约束条件下，最小化异构多雷达网络异步多目标跟踪误差。
然而，在异步多目标跟踪时，在不同雷达节点不同量测信息的到达时刻
$ t_{n,m,k}^{q,s} $ 都进行功率时间联合优化是不合理的。基于此，本文所提算法应用于每个融合采样间隔$ {T_{{\text{fusion}}}} $ 。为了简化后续求解，对于目标q，假设在每个融合采样间隔内，不同雷达网络中的各个雷达节点以相同的发射参数进行照射：\left\{ {\begin{array}{*{20}{c}} {P_{n,m,k}^q = P_{n,m,k}^{q,1} = \cdots = P_{n,m,k}^{q,S_{n,m,k}^q}} \\ {T_{n,m,k}^q = T_{n,m,k}^{q,1} = \cdots = T_{n,m,k}^{q,S_{n,m,k}^q}} \end{array}} \right. (19) 结合式(19)，建立优化模型如下：
\left\{\begin{aligned} & \underset{{{\boldsymbol{\mu}} }_{k}^{q},{{\boldsymbol{P}}}_{\text{t,}k}^{q},{{\boldsymbol{T}}}_{k}^{q},\forall q}{\mathrm{min}}\mathbb{F}\left({{\boldsymbol{\mu}} }_{k}^{q},{{\boldsymbol{P}}}_{\text{t},k}^{q},{{\boldsymbol{T}}}_{k}^{q}\right)\\ & \left\{\begin{aligned} & {{\bar P}_{\mathrm{min}}}\le {P}_{n,m,k}^{q}\le {{\bar P}_{\mathrm{max}}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{CMIMO}} \\ & {P}_{n,m,k}^{q}={P}_{\text{PAR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{PAR}} \\ & {P}_{n,m,k}^{q}={P}_{\text{MSR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{MSR}} \end{aligned}\right. \\ & \left\{\begin{aligned} & {T}_{n,m,k}^{q}={T}_{\text{CMIMO}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{CMIMO}} \\ & {{\bar T}_{\mathrm{min}}}\le {T}_{n,m,k}^{q}\le {{\bar T}_{\mathrm{max}}},\text{ }\, n,m\in {{\boldsymbol{\psi}} }_{\text{PAR}} \\ &{T}_{n,m,k}^{q}={T}_{\text{MSR}},\text{ }\, n,m\in {{\boldsymbol{\psi}} }_{\text{MSR}} \end{aligned}\right. \\ & \displaystyle \sum _{m=1}^{M}{\displaystyle \sum _{n=1}^{{N}_{m}}{P}_{n,m,k}^{q}={P}_{\text{total}}^{q},{\displaystyle \sum _{m=1}^{M}{\displaystyle \sum _{n=1}^{{N}_{m}}{T}_{n,m,k}^{q}={T}_{\text{total}}^{q},\forall q}}} \\ & {\displaystyle \sum _{q=1}^{Q}{\mu }_{n,m,k}^{q}\le 1,\forall n,m} \\ & {\displaystyle \sum _{m=1}^{M}{\displaystyle \sum _{n=1}^{{N}_{m}}{\mu }_{n,m,k}^{q}={L}_{\mathrm{max}},\forall q}} \\ & {\mu }_{n,m,k}^{q}\in \left\{0,1\right\},\forall n,m,q \end{aligned}\right. (20) 其中，
${{\boldsymbol{\psi }}_{{\text{CMIMO}}}}$ 表示集中式MIMO雷达网络节点集合；$ {{\boldsymbol{\psi }}_{{\text{PAR}}}} $ 表示相控阵雷达网络节点集合；$ {{\boldsymbol{\psi }}_{{\text{MSR}}}} $ 表示机械扫描雷达网络节点集合；${{\bar P_{\min }}}$ 和${{\bar P_{\max }}}$ 分别表示集中式MIMO雷达网络中雷达节点辐射功率的下限和上限；$ {P_{{\text{PAR}}}} $ 和$ {P_{{\text{MSR}}}} $ 分别表示相控阵雷达网络和机械扫描雷达网络雷达节点的固定辐射功率；${{\bar T_{\min }}}$ 和${{\bar T_{\max }}}$ 分别表示相控阵雷达网络中雷达节点驻留时间的下限和上限；$ {T_{{\text{CMIMO}}}} $ 和$ {T_{{\text{MSR}}}} $ 分别表示集中式MIMO雷达网络和机械扫描雷达网络的固定驻留时间；$ P_{{\text{total}}}^q $ 和$ T_{{\text{total}}}^q $ 分别表示异构多雷达网络对于目标q的总辐射功率约束和总驻留时间约束；$\displaystyle\sum\nolimits_{q = 1}^Q {\mu _{n,m,k}^q \le 1}$ 表示第k个融合采样间隔内每个雷达节点最多只能跟踪1个目标，$\displaystyle\sum\nolimits_{m = 1}^M {\sum\nolimits_{n = 1}^{{N_m}} {\mu _{n,m,k}^q = {L_{\max }}} }$ 表示第k个融合采样间隔内异构多雷达网络固定分配$ {L_{\max }} $ 个雷达节点去跟踪单个目标。从式(20)中可以看出，第1个约束条件表示不同类型雷达网络中雷达节点辐射功率取值范围；第2个约束条件表示不同类型雷达网络中雷达节点驻留时间取值范围；第3个约束条件表示跟踪单个目标时的异构多雷达网络辐射资源总量限制；第4个约束表示单个雷达节点跟踪目标数量的限制；第5个约束表示跟踪单个目标的雷达节点数量限制。
3.3 优化模型求解
由于式(20)描述的优化问题中存在雷达节点选择指标
$ \mu _{n,m,k}^q $ 以及其他在约束条件中互相耦合的变量$ P_{n,m,k}^q $ 和$ T_{n,m,k}^q $ ，因此，该优化问题是非线性、非凸的。基于此，本文提出了一种非理想检测下面向异步多目标跟踪的异构多雷达网络功率时间联合优化算法，采用基于序列二次规划(Sequential Quadratic Programming, SQP)算法和循环最小法的四步分解算法对式(20)进行求解，具体求解步骤如下：步骤1 求解雷达节点选择问题。
在异构多雷达网络的起始融合采样间隔指定目标q进行跟踪，并设置符合各雷达网络发射参数要求的初始辐射功率矩阵
$ {\boldsymbol{\bar P}}_{{\text{t,}}k}^q $ 和初始驻留时间矩阵$ {\boldsymbol{\bar T}}_k^q $ ，系统根据矩阵$ {\boldsymbol{\bar P}}_{{\text{t,}}k}^q $ 和$ {\boldsymbol{\bar T}}_k^q $ 分配给每个雷达节点相应的辐射资源进行照射。同时，松弛二元变量$ \mu _{n,m,k}^q \in \left\{ {0,1} \right\} $ ，此时的优化模型可以简化为\begin{split} & \mathop {\min }\limits_{{\boldsymbol{\mu }}_k^q,\forall q} \mathbb{F}\left( {{\boldsymbol{\mu }}_k^q,{\boldsymbol{\bar P}}_{{\text{t}},k}^q,{\boldsymbol{\bar T}}_k^q} \right), \\ & {\text{s}}{\text{.t}}{\text{.}}\;\left\{ {\begin{array}{*{20}{c}} {\displaystyle\sum\limits_{q = 1}^Q {\mu _{n,m,k}^q \le 1,\forall n,m} } \\ {0 \le \mu _{n,m,k}^q \le 1,\forall n,m,q} \end{array}} \right. \end{split} (21) 从式(21)可以看出，二元变量
$ \mu _{n,m,k}^q \in \left\{ {0,1} \right\} $ 被松弛为$0 \le \mu _{n,m,k}^q \le 1$ ，且优化模型中的目标函数以及所有约束条件只与$ \mu _{n,m,k}^q $ 相关，根据文献[26]，该优化问题是一个凸问题。对该模型进行求解可以得到雷达节点选择的结果，具体步骤如下：首先，采用SQP算法求解得到雷达节点选择系数矩阵$ {\boldsymbol{\bar \mu }}_k^q $ ；其次，将$ {\boldsymbol{\bar \mu }}_k^q $ 中的选择系数降序排列，依次选取该系数所对应的雷达节点直到选取的数目等于$ {L_{\max }} $ ；最后，将已选取的$ {L_{\max }} $ 个雷达节点对应的二元变量设为$ \mu _{n,m,k}^q = 1 $ ，将未选取的$\displaystyle\sum\nolimits_{m = 1}^M {{N_m}} - {L_{\max }}$ 个雷达节点对应的二元变量设为$ \mu _{n,m,k}^q = 0 $ ，即可得到雷达节点选择结果$ {\boldsymbol{\hat \mu }}_k^q $ 。步骤2 依次求解辐射功率和驻留时间优化分配问题。
根据步骤1得到的雷达节点选择结果
$ {\boldsymbol{\hat \mu }}_k^q $ ，将优化模型重新写为\begin{aligned} & \underset{{{\boldsymbol{P}}}_{\text{t,}k}^{q},{{\boldsymbol{T}}}_{k}^{q},\forall q}{\mathrm{min}}\mathbb{F}\left({\hat{{\boldsymbol{\mu}} }}_{k}^{q},{{\boldsymbol{P}}}_{\text{t},k}^{q},{{\boldsymbol{T}}}_{k}^{q}\right),\\ & \text{s}\text{.t}\text{.:}\left\{\begin{aligned} &\left\{\begin{aligned} & {{\bar P}_{\mathrm{min}}}\le {P}_{n,m,k}^{q}\le {{\bar P}_{\mathrm{max}}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{CMIMO}} \\ & {P}_{n,m,k}^{q}={P}_{\text{PAR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{PAR}} \\ & {P}_{n,m,k}^{q}={P}_{\text{MSR}},\text{ }\, n,m\in {{\boldsymbol{\psi}} }_{\text{MSR}} \end{aligned}\right. \\ &\left\{\begin{aligned} & {T}_{n,m,k}^{q}={T}_{\text{CMIMO}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{CMIMO}} \\ & {{\bar T}_{\mathrm{min}}}\le {T}_{n,m,k}^{q}\le {{\bar T}_{\mathrm{max}}},\text{ }\, n,m\in {{\boldsymbol{\psi}} }_{\text{PAR}} \\ & {T}_{n,m,k}^{q}={T}_{\text{MSR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{MSR}} \end{aligned}\right. \\ & \displaystyle \sum _{m=1}^{M}\displaystyle \sum _{n=1}^{{N}_{m}}{P}_{n,m,k}^{q}={P}_{\text{total}}^{q}\\ & \displaystyle \sum _{m=1}^{M}\displaystyle \sum _{n=1}^{{N}_{m}}{T}_{n,m,k}^{q} ={T}_{\text{total}}^{q},\forall q \end{aligned}\right.\\[-7pt]\end{aligned} (22) 从优化模型(22)可以看出，此时的优化模型包含辐射功率变量
$ P_{n,m,k}^q $ 、驻留时间变量$ T_{n,m,k}^q $ 以及相应的资源约束，因此，它仍然是一个非线性、非凸问题。基于此，对优化模型(22)进行等价转换，可以得到分别优化辐射功率资源和驻留时间资源的两个子优化问题，如下所示：\begin{aligned} & \underset{{{\boldsymbol{P}}}_{\text{t,}k}^{q},\forall q}{\mathrm{min}}\mathbb{F}\left({\hat{{\boldsymbol{\mu}} }}_{k}^{q},{{\boldsymbol{P}}}_{\text{t},k}^{q},{{\boldsymbol{T}}}_{k}^{q}\right)\\ & \text{s}\text{.t}\text{.:}\left\{\begin{aligned} & \left\{\begin{aligned} & {{\bar P}_{\mathrm{min}}}\le {P}_{n,m,k}^{q}\le {{ \bar P}_{\mathrm{max}}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{CMIMO}} \\ & {P}_{n,m,k}^{q}={P}_{\text{PAR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{PAR}} \\ & {P}_{n,m,k}^{q}={P}_{\text{MSR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{MSR}} \end{aligned} \right. \\ & {\displaystyle \sum _{m=1}^{M}{\displaystyle \sum _{n=1}^{{N}_{m}}{P}_{n,m,k}^{q}={P}_{\text{total}}^{q},\forall q}}\\[-7pt] \end{aligned}\right.\end{aligned} (23) \begin{aligned} & \underset{{{\boldsymbol{T}}}_{k}^{q},\forall q}{\mathrm{min}}\mathbb{F}\left({\hat{{\boldsymbol{\mu}} }}_{k}^{q},{{\boldsymbol{P}}}_{\text{t},k}^{q},{{\boldsymbol{T}}}_{k}^{q}\right)\\ & \text{s}\text{.t}\text{.:}\left\{\begin{aligned} &\left\{\begin{aligned} & {T}_{n,m,k}^{q}={T}_{\text{CMIMO}},\text{ }\,\, n,m\in {{\boldsymbol{\psi}} }_{\text{CMIMO}} \\ & {{\bar T}_{\mathrm{min}}}\le {T}_{n,m,k}^{q}\le {{\bar T}_{\mathrm{max}}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{PAR}} \\ & {T}_{n,m,k}^{q}={T}_{\text{MSR}},\text{ } n,m\in {{\boldsymbol{\psi}} }_{\text{MSR}} \end{aligned}\right. \\ & {\displaystyle \sum _{m=1}^{M}{\displaystyle \sum _{n=1}^{{N}_{m}}{T}_{n,m,k}^{q}={T}_{\text{total}}^{q},\forall q}}\\[-7pt] \end{aligned}\right.\end{aligned} (24) 针对优化模型(23)和优化模型(24)进行求解。首先，根据
$ {\boldsymbol{\bar T}}_k^q $ 固定驻留时间，此时，优化模型(23)是一个凸问题[26]，采用SQP算法对其进行求解，可以得到跟踪目标q时的辐射功率最优分配结果$ {\boldsymbol{\hat P}}_{{\text{t,}}k}^q $ ；其次，将得到的$ {\boldsymbol{\hat P}}_{{\text{t,}}k}^q $ 代入优化模型(24)中，同样采用SQP算法进行求解，可以得到驻留时间最优分配结果${\boldsymbol{\hat T}}_k^q$ 。步骤3 指定下一个跟踪目标进行资源优化分配。
移除步骤1中选取的雷达节点以及目标q，指定下一个跟踪目标，重复步骤1和步骤2，直到所有跟踪目标都被分配相应的雷达节点、辐射功率和驻留时间资源。
步骤4 采用循环最小法得到最优资源分配结果。
采用循环最小法，重复步骤1至步骤3，直到连续两次迭代的目标函数值(即异步目标跟踪精度)小于一个预设的固定值，循环终止。基于此，可以得到异构多雷达网络系统异步多目标跟踪的最优雷达节点选择结果
$ {{\boldsymbol{\hat \mu }}_k} $ 、辐射功率分配结果$ {{\boldsymbol{\hat P}}_{{\text{t,}}k}} $ 和驻留时间分配结果$ {{\boldsymbol{\hat T}}_k} $ 。此外，将求解得到的$ {{\boldsymbol{\hat \mu }}_k} $ ,$ {{\boldsymbol{\hat P}}_{{\text{t,}}k}} $ 和$ {{\boldsymbol{\hat T}}_k} $ 作为下一融合采样间隔资源联合优化的初始参数。综上所述，非理想检测下面向异步多目标跟踪的异构多雷达网络功率时间联合优化算法如算法1所示。
1 非理想检测下面向异步多目标跟踪的异构多雷达网络功率时间联合优化算法1. Joint optimization of transmit power and dwell time for asynchronous multi-target tracking in heterogeneous multiple radar networks with imperfect detection步骤1 根据功率和时间初始矩阵分配给各雷达节点相应的辐射
功率$ {\boldsymbol{\bar P}}_{{\text{t,}}k}^q $和驻留时间$ {\boldsymbol{\bar T}}_k^q $，松弛$ \mu _{n,m,k}^q $，采用SQP算法求
解式(21)，得到所有的雷达节点选择系数，降序排列后
依次选取$ {L_{\max }} $个雷达节点；步骤2 将优化模型式(22)分解为两个子优化问题，采用SQP算
法依次迭代求解式(23)和式(24)，得到辐射功率和驻留
时间优化分配结果$ {\boldsymbol{\hat P}}_{{\text{t,}}k}^q $和$ {\boldsymbol{\hat T}}_k^q $；步骤3 异构多雷达网络系统指定下一个目标进行相应的资源优
化分配，重复步骤1和步骤2，直到照射所有目标；步骤4 利用循环最小法得到最优雷达节点选择和功率时间分配
结果。4. 仿真结果及分析
为了验证本文所提非理想检测下面向异步多目标跟踪的异构多雷达网络功率时间联合优化算法的可行性和优越性，本文设计了如下仿真场景。考虑异构多雷达网络系统由
$ M = 3 $ 部独立工作且不同类型的雷达网络组成，其中，雷达网络1为集中式MIMO雷达网络，包含${N_1} = 5$ 个雷达节点；雷达网络2为相控阵雷达网络，包含${N_2} = 5$ 个雷达节点；雷达网络3为机械扫描雷达网络，包含${N_3} = 2$ 个雷达节点，各雷达节点的发射参数如表1所示。此外，异构多雷达网络系统监视区域中的运动目标个数为$ Q = 2 $ ，图2给出了异构多雷达网络分布与目标运动轨迹，目标具体运动轨迹参数如表2所示。同时，表3给出了每个融合采样间隔内各雷达节点的初始采样时间和采样间隔。在仿真场景中，假设异构多雷达网络固定分配不同雷达网络中共计$ {L_{\max }} = 5 $ 个雷达节点去跟踪每个目标，且将目标RCS固定为$ 1\;{{\text{m}}^{\text{2}}} $ 。表 1 异构多雷达网络雷达节点发射参数设置Table 1. Parameter settings of radar node in heterogeneous multiple radar networks参数 数值 参数 数值 ${T_{{\text{fusion}}}}$ $6\;{\text{s} }$ ${ {\bar P_{\min } } }$ $100\;{\text{W}}$ ${ {\bar P_{ {} {\text{max} } } } }$ $1000\;{\text{ W} }$ $ {P_{{\text{PAR}}}} $ $300\;{\text{W}}$ $ {P_{{\text{MSR}}}} $ $300\;{\text{W} }$ ${ {\bar T_{\min } } }$ $0.003\;{\text{s}}$ ${ {\bar T_{\max } } }$ $0.03\;{\text{s} }$ $ {T_{{\text{CMIMO}}}} $ $0.006\;{\text{s}}$ $ {T_{{\text{MSR}}}} $ $0.006\;{\text{s}}$ $ \beta _k^q $ $1\;{\text{MHz}}$ $ P_{{\text{total}}}^q $ $1500\;{\text{W}}$ $ T_{{\text{total}}}^q $ $0.03\;{\text{s}}$ 
图 2 异构多雷达网络分布与目标运动轨迹Figure 2. Deployment of heterogeneous multiple radar networks and trajectories of moving targets表 2 目标具体运动轨迹参数Table 2. Specific trajectory parameters of targets目标编号 初始位置($ {\text{km}} $) 运动速度($ {{\text{m}} \mathord{\left/ {\vphantom {{\text{m}} {\text{s}}}} \right. } {\text{s}}} $) 1 $ \left( { - 10,42.5} \right) $ $ \left( { - 125, - 340} \right) $ 2 $ \left( {30, - 35} \right) $ $ \left( {160,300} \right) $ 表 3 各雷达节点的初始采样时间和采样间隔Table 3. Initial sampling time and sampling interval of each radar node$ \left( {n,m} \right) $ $t_{n,m,q}^{{\text{ini} } }$ (s) $ T_{n,m,q}^0 $ (s) (1,1) 2 1 (2,1) 3 1 (3,1) 4 1 (4,1) 1 1 (5,1) 5 1 (1,2) 2 2 (2,2) 4 2 (3,2) 3 2 (4,2) 5 2 (5,2) 1 2 (1,3) 3 3 (2,3) 5 3 为了更好地说明异构多雷达网络雷达节点选择与功率时间联合优化的规律与内在关系，图3和图4分别给出了
$p_{{\text{D}},n,m,k}^q = 0.9$ 时目标1和目标2的雷达节点选择与功率时间优化结果，黑色部分表示该雷达节点在对应融合采样间隔未被赋予跟踪任务，即$ \mu _{n,m,k}^q = 0 $ 。从中可以看出，异构多雷达网络会优先选择距离目标更近雷达网络中的雷达节点进行目标跟踪，各雷达节点的辐射功率和驻留时间资源随着目标的运动自适应地进行优化分配，从而最大限度地降低异构多雷达网络异步多目标跟踪误差，提升多目标跟踪性能。例如，对于目标1，在跟踪的前25个融合采样间隔，异构多雷达网络都选择了集中式MIMO雷达网络的第1~4个节点进行跟踪，在最后10个融合采样间隔，由于目标1运动到距离相控阵雷达网络更近的位置，系统会逐渐选择相控阵雷达网络中的雷达节点进行跟踪。此外，雷达节点的初始采样时间以及采样间隔也会影响雷达节点选择方式与功率时间优化结果，异构多雷达网络会优先选择在融合采样间隔内次数较多的雷达节点，以降低多目标跟踪误差。例如，对于目标2，异构多雷达网络更多地选择了集中式MIMO雷达网络以及相控阵雷达网络中的雷达节点，而并不选择距离目标2更近的机械扫描雷达网络中的第2个节点进行跟踪。分析其原因为：机械扫描雷达网络中的第2个节点采样间隔为3 s，大于另外两个雷达网络中雷达节点的采样间隔，因此，在同一个融合采样间隔中，其采样次数较少，能够获取的雷达量测信息数目也较少，这会导致异步多目标跟踪精度变差。
图 3 目标1的雷达节点选择与功率时间优化结果Figure 3. Radar node selection, power and dwell time optimization results of target 1
图 4 目标2的雷达节点选择与功率时间优化结果Figure 4. Radar node selection, power and dwell time optimization results of target 2为了更好地展现检测概率变化对异构多雷达网络异步多目标跟踪精度的影响，图5给出了不同检测概率下各目标的BCRLB对比。从图中可以看出，相较于理想检测条件(即检测概率为1.0)，非理想检测下目标状态估计误差的BCRLB更大，这表明目标跟踪误差也更大。此外，随着检测概率的不断降低，异构多雷达网络异步多目标跟踪的检测环境逐渐变得恶劣，当检测概率从0.9降低到0.8时，各目标跟踪的BCRLB均明显增大，这表明在系统可用辐射资源限制相同的情况下，检测概率变化会对异步多目标跟踪精度产生明显影响。
图 5 不同检测概率下的各目标BCRLBFigure 5. BCRLB of each target with different probabilities of detection为了更有效地说明本文所提算法对异构多雷达网络异步多目标跟踪性能的提升，本节采用如下4种算法作为对比，比较分析它们与本文所提算法之间的性能差异：
(1) 节点随机选择算法：在每个融合采样间隔，分配给各个目标的雷达节点是随机选择的，同时通过求解优化问题(20)，将辐射功率和驻留时间的优化结果分配给所选择的雷达节点；
(2) 功率均匀分配算法：固定异构多雷达网络跟踪目标的辐射功率，仅优化驻留时间资源；
(3) 时间均匀分配算法：固定异构多雷达网络跟踪目标的驻留时间，仅优化辐射功率资源；
(4) 功率时间均匀分配算法：同时固定异构多雷达网络跟踪目标的辐射功率资源和驻留时间资源。
定义目标跟踪的平均均方根误差(Average Root Mean Square Error, ARMSE)来表征异步多目标跟踪精度：
{\text{ARMSE}}\left( q \right) \triangleq \\ {\text{ }}\sqrt {\frac{1}{{{N_{{{\rm{MC}}} }}}}\sum\limits_{l = 1}^{{N_{{{\rm{MC}}} }}} {\frac{1}{{N_k^q\left( l \right)}}\sum\limits_{k = 1}^{N_k^q\left( l \right)} {\left\{ {{{\left[ {x_k^q - \hat x_{l,k\left| k \right.}^q} \right]}^2} + {{\left[ {y_k^q - \hat y_{l,k\left| k \right.}^q} \right]}^2}} \right\}} } } (25) 其中，
${N_{{{\rm{MC}}} }} = 100$ 为蒙特卡罗实验次数；$ N_k^q\left( l \right) $ 表示目标q在第l次蒙特卡罗试验时被异构多雷达网络辐射的总次数。图6给出了所有算法的各目标ARMSE对比图。从图中可以看出，本文所提算法能够在满足异构多雷达网络射频辐射资源限制的约束条件下，通过联合优化系统中各雷达节点选择方式、辐射功率和驻留时间，求解得到相比于其他4种对比算法更优越的异步多目标跟踪精度。此外，本文所提算法在检测概率为1.0, 0.9和0.8时的单帧CPU计算耗时分别为0.46 s, 1.03 s和1.08 s，远远小于预设的融合采样间隔。因此，本文所提算法可以在满足作战场景实时性的同时，有效提升异步多目标跟踪性能。
5. 结语
本文考虑了非理想检测、异构多雷达网络和异步多目标跟踪对射频辐射资源消耗以及多目标跟踪性能的影响，提出了一种面向异步多目标跟踪的异构多雷达网络功率时间联合优化算法。该算法以最小化异构多雷达网络异步多目标跟踪误差为优化目标，在满足给定系统射频辐射资源限制的前提下，对雷达节点选择方式、辐射功率和驻留时间等发射参数进行自适应优化设计，以提升异构多雷达网络的异步多目标跟踪精度。仿真结果表明，与节点随机选择算法、功率均匀分配算法、时间均匀分配算法以及功率时间均匀分配算法相比，本文所提算法能够有效降低异构多雷达网络的异步多目标跟踪误差，达到提升其目标跟踪性能的目的。然而，在实际作战场景中，检测概率与目标RCS、密集环境杂波等多种不确定因素息息相关，在多目标跟踪过程中，检测概率会不断发生变化，下一步将围绕检测概率变化下面向异步多目标跟踪的异构多雷达网络多域资源联合优化问题进行研究。
-
图 2 主站L波段全极化SAR和双向同步链集成
Figure 2. Master L-band full-polarization SAR and bidirectional synchronous chain integration
图 3 从站L波段双极化接收系统和双向同步链集成
Figure 3. Slave L-band dual-polarization receiving system and bidirectional synchronous chain integration
图 5 通过UWB技术实现双机协同定位技术
Figure 5. Dual-machine co-positioning technology is realized through UWB technology
图 7 UWB测距和差分GPS定位的融合处理流程图
Figure 7. Flowchart of the fusion processing of UWB ranging and differential GPS positioning
图 8 双站InSAR系统双向同步链示意图
Figure 8. Diagram of the bidirectional synchronous chain of bistatic InSAR system
图 10 从站雷达时间和相位同步后的测量结果
Figure 10. Measurements after time and phase synchronization of slave radar
图 11 无人机载双站InSAR成像处理流程图
Figure 11. Flow chart of UAV-borne bistatic InSAR imaging processing
图 13 无人机载双站InSAR飞行的典型航迹和姿态测量数据
Figure 13. Typical trajectory and attitude measurements of UAV-borne bistatic InSAR
图 18 某三面角成像结果的方位向曲线
Figure 18. Curve of trihedral corner reflector in azimuth in imaging results
图 19 干涉相位误差与距离向像素关系
Figure 19. Relation of interference phase error and pixel location in range direction
图 20 极化校正前后的主站Pauli图
Figure 20. Pauli diagrams of the master station before and after polarization correction
图 24 去平地和滤波后的相位图
Figure 24. Interference phase diagram after flat-earth phase removal and filtering
图 28 去平地及滤波后的干涉相位图
Figure 28. Interference phase diagram after flat-earth phase removal and filtering
表 1 无人机载双站InSAR系统误差分配
Table 1. Systematic error distribution of UAV-borne bistatic InSAR
参数名称 参数值 飞行高度(相对于地面) 2 km 入射角 45° 基线角 0° 有效基线长度 30 m 水平基线长度 42.43 m 航迹控制精度(分配值) 1 m 时间同步引入的包络对齐误差 0.1个采样单元 运补残余误差导致的干涉相位误差 1° 同步链的干涉相位误差(分配值) 5° 去相干引入的相位误差 2.54° 基线测量误差(分配值) 4 mm 基线角误差(分配值) 0.003° 高程误差(理论分析值) 0.46 m 
下载: 导出CSV
表 2 无人机载双站InSAR系统总体构成
Table 2. Overall composition of the UAV-borne bistatic InSAR system
名称 说明 L-SAR(主站) 由L波段全极化SAR和双向同步链组成，其中，L波段全极化SAR由信号产生与发射通道、功放、双极化接收通道、天线及天线支架、开关组合、数据采集存储模块等组成 L-接收系统(从站) 由L波段双极化接收系统和双向同步链组成，其中L波段双极化接收系统由双极化接收通道、天线及天线支架、数据采集存储模块等组成 无人机双机编队 由两架固定翼HC-140无人机和两套协同控制器组成，其中，单架无人机最大作业载荷30 kg、最大翼展尺寸5.4 m，最大飞行速度45 m/s 位姿测量系统 由差分GPS模块和微型惯性测量组成；测量数据更新频率为200 Hz，地面后处理后，航迹测量精度0.05 m，偏航角测量精度0.02°，横滚和俯仰角测量精度0.015°
下载: 导出CSV
表 3 L-SAR载荷参数
Table 3. L-SAR load parameters
名称 参数名称 参数值 系统 波段 L波段 系统同步误差 ≤5° 噪声等效后向散射系数(NEσ0) ≤40 dB 主站 分辨率 ≤0.5 m 最大信号带宽 400 MHz 最大作用距离 ≥6 km 最大测绘带宽 ≥3 km 主站/从站 极化方式 全极化 极化隔离度 ≥25 dBc 极化通道相位不平衡度 ≤8° 幅度不平衡度 ≤0.3 dB
下载: 导出CSV
表 4 双站InSAR构型参数
Table 4. Configuration parameters of bistatic InSAR
模式 飞行高度(m) 基线(m) 说明 模式1 1000 30 双站角0.87o，基线长度探底双机协同飞行安全距离的下限 模式2 2000 50 双站角0.72o，各项参数均比较合适
下载: 导出CSV
表 5 分辨率测试结果
Table 5. Results of resolution test
分辨率 距离向(m) 方位向(m) 主站 0.3747 0.4534 从站 0.3747 0.4860
下载: 导出CSV
表 6 主站数据极化定标前后的三面角极化质量评价
Table 6. Polarization quality evaluation of master data before and after polarization calibration based on trihedral corners
状态 三面角编号 隔离度(dB) 幅度不平衡(dB) 相位不平衡(°) 极化校正前 J2* 29.3217 0.21983 99.8529 J3 36.6858 0.44359 97.3340 J4 28.0391 –0.05674 98.3664 极化校正后 J2* 322.3542 5.786E–15 –7.2788E–15 J3 32.2467 0.15738 –3.2333 J4 37.7835 –0.32067 –1.2498 注：表中带*的J2为定标中使用的三面角，其余定标中未用
下载: 导出CSV
表 7 从站数据极化定标前后的三面角极化质量评价
Table 7. Polarization quality evaluation of slave data before and after polarization calibration based on trihedral corners
状态 三面角编号 隔离度(dB) 幅度不平衡(dB) 相位不平衡(°) 极化校正前 J2* 28.4590 –0.23950 –90.2117 J3 26.8158 –0.08019 –87.2605 J4 27.9540 –0.71433 –86.5173 极化校正后 J2* 311.252 –2.893E–15 3.317E–14 J3 31.8702 0.065833 2.3704 J4 33.8167 –0.493960 3.3900 注：表中带*的J2为定标中使用的三面角，其余定标中未用
下载: 导出CSV
表 8 双站SAR干涉反演高度误差表
Table 8. Table of elevation inversion error of bistatic interferometric SAR
定标点 相干系数 反演高度(m) 实测高度(m) 误差(m) J1 0.97 1385.89 1385.71 0.18 J2 0.89 1383.39 1383.23 0.17 J3 0.93 1385.14 1384.86 0.28 J4 0.97 1382.56 1383.68 –1.12 J5 0.86 1383.23 1384.45 –1.21 J6 0.98 1383.81 1384.38 –0.56 J7 0.99 1382.40 1382.11 0.29 J8 0.96 1387.60 1387.05 0.55 J9 0.92 1383.41 1384.28 –0.87 J10 0.87 1380.68 1383.05 –2.37 J11 0.98 1382.76 1383.38 –0.61 J12 0.96 1383.26 1383.61 –0.35 J13 0.98 1383.29 1382.98 0.31 J14 0.96 1385.93 1386.03 –0.09 J15 0.66 1386.80 1387.79 –1.00 J16 0.89 1388.00 1387.68 0.32 J17 0.48 1385.26 1387.25 –1.99 J18 0.88 1386.49 1386.65 –0.17 J19 0.97 1385.82 1386.09 –0.27 J20 0.93 1385.39 1385.57 –0.18 J21 0.98 1385.53 1385.23 0.30
下载: 导出CSV
-
[1] ROSEN P A, HENSLEY S, JOUGHIN I R, et al. Synthetic aperture radar interferometry[J]. Proceedings of the IEEE, 2000, 88(3): 333–382. doi: 10.1109/5.838084 [2] 仇晓兰, 丁赤飚, 胡东辉. 双站SAR成像处理技术[M]. 北京: 科学出版社, 2010: 8–11.QIU Xiaolan, DING Chibiao, and HU Dong-hui. Bistatic SAR Imaging Algorithms[M]. Beijing: Science Press, 2010: 8–11. [3] 孙亚飞, 江利明, 柳林, 等. TanDEM-X双站SAR干涉测量及研究进展[J]. 国土资源遥感, 2015, 27(1): 16–22. doi: 10.6046/gtzyyg.2015.01.03SUN Yafei, JIANG Liming, LIU Lin, et al. TanDEM-X bistatic SAR interferometry and its research progress[J]. Remote Sensing for Natural Resources, 2015, 27(1): 16–22. doi: 10.6046/gtzyyg.2015.01.03 [4] 章皖秋, 岳彩荣, 颜培东. TanDEM-X极化干涉SAR森林冠层高度反演[J]. 东北林业大学学报, 2017, 45(1): 47–54. doi: 10.3969/j.issn.1000-5382.2017.01.011ZHANG Wanqiu, YUE Cairong, and YAN Peidong. Forest canopy height retrieval by PolInSAR with TanDEM-X data[J]. Journal of Northeast Forestry University, 2017, 45(1): 47–54. doi: 10.3969/j.issn.1000-5382.2017.01.011 [5] LIANG Da, LIU Kaiyu, ZHANG Heng, et al. The processing framework and experimental verification for the noninterrupted synchronization scheme of LuTan-1[J]. IEEE Transactions on Geoscience and Remote Sensing, 2021, 59(7): 5740–5750. doi: 10.1109/TGRS.2020.3024561 [6] LIN Haoyu, DENG Yunkai, ZHANG Heng, et al. On the processing of dual-channel receiving signals of the LuTan-1 SAR System[J]. Remote Sensing, 2022, 14(3): 515. doi: 10.3390/rs14030515 [7] HENDRIKS I, META A, TRAMPUZ C, et al. MetaSensing’s novel L-band airborne SAR sensors for the BelSAR project: First Bistatic results[EB/OL]. https://www.researchgate.net/profile/Izzy-Hendriks/publication/326317566_MetaSensing%27s_Novel_L-Band_Airborne_SAR_Sensors_for_the_BelSAR_Project_First_Bistatic_Results/links/5b45af4b0f7e9b1c72236322/MetaSensings-Novel-L-Band-Airborne-SAR-Sensors-for-the-BelSAR-Project-First-Bistatic-Results.pdf, 2017. [8] META A, TRAMPUZ C, COCCIA A, et al. First results of the BelSAR L band airborne bistatic fully polarimetric Synthetic aperture radar campaign[C]. 2017 IEEE International Geoscience and Remote Sensing Symposium, Fort Worth, USA, 2017: 1040–1042. [9] YANG Jianyu, HUANG Yulin, YANG Haiguang, et al. A first experiment of airborne bistatic forward-looking SAR - Preliminary results[C]. 2013 IEEE International Geoscience and Remote Sensing Symposium, Melbourne, Australia, 2013: 4202–4204. [10] LI Zhongyu, WU Junjie, YI Qingying, et al. Bistatic forward-looking SAR ground moving target detection and imaging[J]. IEEE Transactions on Aerospace and Electronic Systems, 2015, 51(2): 1000–1016. doi: 10.1109/TAES.2014.130539 [11] LI Zhongyu, YE Hongda, LIU Zhutian, et al. Bistatic SAR clutter-ridge matched STAP method for nonstationary clutter suppression[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 5216914. doi: 10.1109/TGRS.2021.3125043 [12] 傅志豪. 无人机载干涉SAR系统性能分析及应用研究[D]. [硕士论文], 应急管理部国家自然灾害防治研究院, 2021.FU Zhihao. Performance analysis and application of UAV SAR system[D]. [Master dissertation], National Institute of Natural Hazards, 2021. [13] 仇晓兰, 焦泽坤, 杨振礼, 等. 微波视觉三维SAR关键技术及实验系统初步进展[J]. 雷达学报, 2022, 11(1): 1–19. doi: 10.12000/JR22027QIU Xiaolan, JIAO Zekun, YANG Zhenli, et al. Key technology and preliminary progress of microwave vision 3D SAR experimental system[J]. Journal of Radars, 2022, 11(1): 1–19. doi: 10.12000/JR22027 [14] LV Zexin, QIU Xiaolan, CHENG Yao, et al. Multi-rotor UAV-borne PolInSAR data processing and preliminary analysis of height inversion in urban area[J]. Remote Sensing, 2022, 14(9): 2161. doi: 10.3390/rs14092161 [15] 林春辉. 单基/双基SAR成像若干关键问题研究[D]. [博士论文], 西安电子科技大学, 2019.LIN Chunhui. Study on some imaging issues of monostatic and bistatic SAR[D]. [Ph.D. dissertation], Xidian University, 2019. [16] 刘松林. 双基SAR时频同步系统研究[D]. [硕士论文], 南京航空航天大学, 2016.LIU Songlin. Research on time and frequency synchronization system of bistatic SAR[D]. [Master dissertation], Nanjing University of Aeronautics and Astronautics, 2016. [17] 汤晓涛, 楼良盛, 刘志铭. 编队卫星InSAR系统的相位同步分析[J]. 地球信息科学, 2008, 10(6): 798–801. doi: 10.3969/j.issn.1560-8999.2008.06.020TANG Xiaotao, LOU Liangsheng, and LIU Zhiming. Analyses of phase synchronization on InSAR system based on formation- flying satellites[J]. Journal of Geo-information Science, 2008, 10(6): 798–801. doi: 10.3969/j.issn.1560-8999.2008.06.020 [18] 丁赤飚, 李芳芳, 胡东辉, 等. 机载干涉合成孔径雷达数据处理技术[M]. 北京: 科学出版社, 2017: 17–36.DING Chibiao, LI Fangfang, HU Donghui, et al. Data Processing Technology of Airborne Interferometric Synthetic Aperture Radar[M]. Beijing: Science Press, 2017: 17–36. [19] 李芳芳, 仇晓兰, 孟大地, 等. 机载双天线InSAR运动补偿误差的影响分析[J]. 电子与信息学报, 2013, 35(3): 559–567. doi: 10.3724/SP.J.1146.2012.00850LI Fangfang, QIU Xiaolan, MENG Dadi, et al. Effects of motion compensation errors on performance of airborne dual-antenna InSAR[J]. Journal of Electronics &Information Technology, 2013, 35(3): 559–567. doi: 10.3724/SP.J.1146.2012.00850 [20] 刘琦, 岳彩荣, 章皖秋, 等. 极化干涉SAR森林冠层高反演的地形坡度改正[J]. 东北林业大学学报, 2017, 45(1): 55–60, 70. doi: 10.3969/j.issn.1000-5382.2017.01.012LIU Qi, YUE Cairong, ZHANG Wanqiu, et al. Terrain slope correction on PolInSAR forest canopy height inversion[J]. Journal of Northeast Forestry University, 2017, 45(1): 55–60, 70. doi: 10.3969/j.issn.1000-5382.2017.01.012 [21] 朱刚. 超宽带(UWB)原理与干扰[M]. 北京: 清华大学出版社, 2009: 2–5.ZHU Gang. Ultra Wideband (UWB) Principle and Interference[M]. Beijing: Tsinghua University Press, 2009: 2–5. [22] LUECKEN H, STEINER C, and WITTNEBEN A. Location-aware UWB communication with generalized energy detection receivers[J]. IEEE Transactions on Wireless Communications, 2012, 11(9): 3068–3078. doi: 10.1109/TWC.2012.070912.110101 [23] 潘莉娟. 星间高精度时间同步和测距系统的研究[D]. [硕士论文], 中国科学技术大学, 2008.PAN Lijuan. Research on high-precision time synchronization and ranging systems between satellites[D]. [Master dissertation], University of Science and Technology of China, 2008. [24] 张方辉, 梁兴东, 周良将. 双站SAR时间同步误差建模及分析[J]. 国外电子测量技术, 2010, 29(8): 36–40. doi: 10.3969/j.issn.1002-8978.2010.08.013ZHANG Fanghui, LIANG Xingdong, and ZHOU Liangjiang. Modeling and analyzing of time synchronization errors in bistatic SAR[J]. Foreign Electronic Measurement Technology, 2010, 29(8): 36–40. doi: 10.3969/j.issn.1002-8978.2010.08.013 [25] YOUNIS M, METZIG R, and KRIEGER G. Performance prediction of a phase synchronization link for Bistatic SAR[J]. IEEE Geoscience and Remote Sensing Letters, 2006, 3(3): 429–433. doi: 10.1109/LGRS.2006.874163 [26] 雷科. 机载双基地合成孔径雷达系统同步问题研究[D]. [硕士论文], 电子科技大学, 2008.LEI Ke. Research on synchronization of airborne bistatic synthetic-aperture radar system[D]. [Master dissertation], University of Electronic Science and Technology of China, 2008. [27] WEIB M. Synchronisation of bistatic radar systems[C]. 2004 IEEE International Geoscience and Remote Sensing Symposium, Anchorage, USA, 2004: 1750–1753. [28] 向建冰, 吕孝雷, 付希凯, 等. 天绘二号双星InSAR成像与DSM生成技术[J]. 测绘学报, 2022, 51(12): 2493–2500. doi: 10.11947/j.AGCS.2022.20210373XIANG Jianbing, LÜ Xiaolei, FU Xikai, et al. Bistatic InSAR interferometry imaging and DSM generation for TH-2[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(12): 2493–2500. doi: 10.11947/j.AGCS.2022.20210373 [29] 孟大地. 机载合成孔径雷达运动补偿算法研究[D]. [博士论文], 中国科学院电子学研究所, 2006.MENG Dadi. Research on motion compensation algorithm for airborne SAR[D]. [Ph. D. dissertation], Institute of Electronics, Chinese Academy of Sciences, 2006. [30] QIU Xiaolan, HAN Bin, MENG Dadi, et al. An azimuth resample method for bistatic SAR motion compensation[C]. 8th European Conference on Synthetic Aperture Radar, Aachen, Germany, 2010: 1–4. [31] DALL J, GRINDER-PEDERSEN J, and MADSEN S N. Calibration of a high resolution airborne 3D SAR[C]. IEEE International Geoscience and Remote Sensing Symposium Proceedings. Remote Sensing - A Scientific Vision for Sustainable Development, Singapore, 1997: 1018–1021. [32] 张薇. 机载双天线干涉SAR定标方法研究[D]. [博士学位论文], 中国科学院电子学研究所, 2009.ZHANG Wei. Airborne dual-antenna InSAR’s interferometric calibration method research[D]. [Ph.D. dissertation], Institute of Electrics, Chinese Academy of Sciences, 2009. [33] 吴迪, 李焱磊, 周良将, 等. 一种基于Whitt算法的SAR极化定标改进方法[J]. 雷达科学与技术, 2018, 16(2): 125–132. doi: 10.3969/j.issn.1672-2337.2018.02.002WU Di, LI Yanlei, ZHOU Liangjiang, et al. An improved method for SAR polarimetric calibration based on Whitt algorithm[J]. Radar Science and Technology, 2018, 16(2): 125–132. doi: 10.3969/j.issn.1672-2337.2018.02.002 [34] LI Fangfang, DING Chibiao, ZHANG Yueting, et al. Airborne InSAR interferometric phase analysis, unwrapping method, and fast implementation in low coherence areas[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2020, 13: 5241–5250. doi: 10.1109/JSTARS.2020.3020148 [35] 黄海风, 张永胜, 董臻. 星载合成孔径雷达干涉新技术[M]. 北京: 科学出版社, 2015: 52–55.HUANG Haifeng, ZHANG Yongsheng, and DONG Zhen. New Interferometric Technology of Spaceborne Synthetic Aperture Radar[M]. Beijing: Science Press, 2015: 52–55. [36] GOLDSTEIN R M and WERNER C L. Radar interferogram filtering for geophysical applications[J]. Geophysical Research Letters, 1998, 25(21): 4035–4038. doi: 10.1029/1998GL900033 [37] BARAN I, STEWART M P, KAMPES B M, et al. A modification to the Goldstein radar interferogram filter[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(9): 2114–2118. doi: 10.1109/TGRS.2003.817212 [38] 张俊娜, 邓喀中, 范洪冬, 等. InSAR相位解缠方法应用比较[J]. 现代测绘, 2011, 34(4): 4. doi: 10.3969/j.issn.1672-4097.2011.04.004ZHANG Junna, DENG Kazhong, FAN Hongdong, et al. Comparison on application of InSAR phase unwrapping methords[J]. Modern Surveying and Mapping, 2011, 34(4): 4. doi: 10.3969/j.issn.1672-4097.2011.04.004 -
下载:
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 0
