基于DP-TBD的协同组网雷达方位标定误差估计与目标检测.pdf

第4 4卷 第3期2 0 2 3年9月制 导 与 引 信GU I D AN C E&F U Z E V o l.4 4 N o.3S e p.2 0 2 3文章编号:1 6 7 1-0 5 7 6(2 0 2 3)0 3-0 0 0 8-0 8基于D P-T B D的协同组网雷达方位标定误差估计与目标检测吴熙芃1,孙希平1,赵泽亚2,吴建新1,金 雪2(1.中山大学深圳电子与通信工程学院,广东 深圳 5 1 8 1 0 7;2.北京跟踪与通信技术研究所,北京 1 0 0 0 9 4)摘 要:提出一种基于动态规划-检测前跟踪(d y n a m i c p r o g r a mm i n g t r a c k-b e f o r e-d e t e c t,D P-T B D)的协同组网雷达方位角标定误差估计与目标检测算法。该算法将D P-T B D方法和多雷达协同处理相结合,将各雷达的坐标系统一到基准雷达坐标系上,利用多部雷达进行D P-T B D处理,得到各雷达的多帧累积值函数,直接利用量测数据估计组网雷达的方位角标定误差,再进行多部雷达的协同积累处理,实现目标检测。该算法可显著提高低信噪比条件下的目标检测性能。仿真结果验证了该算法的有效性。关键词:动态规划-检测前跟踪;方位角标定误差;协同组网雷达中图分类号:T N 9 5 7.5 2文献标志码:AD O I:1 0.3 9 6 9/j.i s s n.1 6 7 1-0 5 7 6.2 0 2 3.0 3.0 0 2收稿日期:2 0 2 2-0 6-2 2基金项目:国家自然科学基金(6 2 1 0 1 6 0 3);航天科学基金(2 0 1 9 2 0 0 M 1 0 0 1);广东省自然科学基金(2 0 2 1 A 1 5 1 5 0 1 1 9 7 9);中山大学中央高校基本科研业务费专项资金(2 2 q n t d 0 4 0 1)作者简介:吴熙芃(1 9 9 5),男,博士研究生,主要从事雷达目标检测与跟踪研究。通信作者:吴建新(1 9 8 2),男,博士,副教授,主要从事雷达目标检测与抗干扰研究。A z i m u t h C a l i b r a t i o n E r r o r E s t i m a t i o n a n d T a r g e t D e t e c t i o nf o r C o o p e r a t i v e N e t w o r k e d R a d a r s B a s e d o n D P-T B DWU X i p e n g1,S UN X i p i n g1,ZHA O Z e y a2,WU J i a n x i n1,J I N X u e2(1.S c h o o l o f E l e c t r o n i c s a n d C o mm u n i c a t i o n E n g i n e e r i n g,S h e n z h e n C a m p u s o f S u nY a t-s e n U n i v e r s i t y,S h e n z h e n 5 1 8 1 0 7,G u a n g d o n g,C h i n a;2.B e i j i n g I n s t i t u t e o fT r a c k i n g a n d T e l e c o mm u n i c a t i o n s T e c h n o l o g y,B e i j i n g 1 0 0 0 9 4,C h i n a)A b s t r a c t:A z i m u t h c a l i b r a t i o n e r r o r e s t i m a t i o n a n d t a r g e t d e t e c t i o n a l g o r i t h m b a s e d o n d y n a m i c p r o g r a mm i n g t r a c k-b e f o r e-d e t e c t(D P-T B D)f o r c o o p e r a t i v e n e t w o r k e d r a d a r s w a s p r o p o s e d.T h e a l g o r i t h m c o m b i n e s t h e D P-T B D m e t h o d a n d m u l t i-r a d a r c o o p e r a t i v e p r o c e s s i n g,u n i f i e s t h e c o o r d i n a t e s y s t e m s o f e a c h r a d a r t o t h e r e f e r e n c e r a d a r c o o r d i n a t e s y s t e m,a n d o b t a i n s t h e m u l t i-f r a m e m e r i t f u n c t i o n o f e a c h r a d a r b y u s i n g D P-T B D p r o c e s s i n g o f m u l t i p l e r a d a r s,w h i c h c a n e s t i m a t e t h e a z i m u t h c a l i b r a t i o n e r r o r o f t h e n e t w o r k e d r a d a r s d i r e c t l y f r o m t h e m e a s u r e m e n t d a t a,a n d t h e n p e r f o r m s c o o p e r a t i v e 第3期吴熙芃,等:基于D P-T B D的协同组网雷达方位标定误差估计与目标检测a c c u m u l a t i o n p r o c e s s i n g o f m u l t i p l e r a d a r s t o a c h i e v e t a r g e t d e t e c t i o n.T h e a l g o r i t h m c a n s i g n i f i c a n t l y i m p r o v e t h e t a r g e t d e t e c t i o n p e r f o r m a n c e a t l o w s i g n a l-t o-n o i s e r a t i o.T h e s i m u l a t i o n r e s u l t s v e r i f y t h e e f f e c t i v e n e s s o f t h e a l g o r i t h m.K e y w o r d s:D P-T B D;a z i m u t h c a l i b r a t i o n e r r o r;c o o p e r a t i v e n e t w o r k e d r a d a r s0 引言 在多雷达组网的目标检测与跟踪系统中,雷达系统标定误差的存在将直接影响组网雷达数据融合后的目标检测与跟踪。其方位角标定误差的存在将导致雷达量测数据不准确、量测结果可信度降低。因此寻找有效的雷达方位角标定方法具有重要的现实意义。组网雷达系统方位角标定误差的估计与校正,是指组网雷达系统误差配准中的空间配准。首先对各雷达的方位角标定误差进行估计,然后再利 用 该 估 计 值 对 雷 达 系 统 量 测 数 据 进 行校正1。在传统组网雷达信号处理中,估计雷达系统标定误差的方法主要分为两类2:第一类是将雷达量测数据与理想的目标位置数据进行比较,估计系统的标定误差;第二类是利用多部雷达对同一目标的量测数据进行处理并对比,估计系统的标定误差。第二类方法主要考虑的是固定平台雷达的空间配准3,包括实时质量控制误差配准算法4、最小二乘/广义最小二乘误差配准算法5-6、精确极大似然配准算法7、基于大地坐标系的误差配准算法8等。上述方法均需要对各雷达的量测数据进行处理,并对目标状态和系统标定误差进行估计。这些方法在低信噪比环境下估计得到的系统标定误差较大,会影响目标的检测与跟踪。动 态 规 划-检 测 前 跟 踪(d y n a m i c p r o g r a mm i n g t r a c k-b e f o r e-d e t e c t,D P-T B D)是一种多帧回波数据联合处理的非相参积累方法,能够在提高雷达对微弱目标的检测性能的同时,对目标的状态信息(距离、方位角等)进行估计。相比于 传 统 的 跟 踪 前 检 测(d e t e c t b e f o r e t r a c k,D B T)方法,D P-T B D方法能够最大程度地保留目标的状态信息9-1 0。在检测前跟踪的众多实现方法中,D P-T B D方法由于其广泛的适用性及易于实现等特点,成为学者们研究的热点1 1-1 4。针对传统方位角标定误差校正方法在低信噪比环境下估计得到的系统标定误差较大的问题,本文采用多部固定平台雷达对同一目标的量测数据进行处理的方法,对各雷达的方位角标定误差进行估计。该算法利用多部雷达协同量测数据,将各雷达的坐标系统一到基准雷达(假设该雷达不存在方位角标定误差)坐标系上,利用多部雷达进行D P-T B D处理,得到各雷达的多帧累积值函数,估计各雷达的方位角标定误差,再进行多部雷达的协同积累处理,对目标进行检测。1 算法原理1.1 方位角误差建模及多雷达坐标统一 假设 组 网 雷 达 系 统 包 含M部 地 基 雷 达,M部雷达对同一个目标进行观测,各雷达坐标已知,但由伺服标定引起的方位角标定误差未知。建立各雷达和目标的坐标系,雷达与目标位置及方位角标定误差如图1所示。各雷达采用数字波束形成(d i g i t a l b e a m f o r m i n g,D B F)的方式测量目标方位角。在各雷达坐标系中,雷达m的目标观测方位角和距离分别表示为 m和Rm,目标的真实方位角为m,目标的方位角标定误差为 m,其中m=1,2,M;雷达间的基线长度为dh,其中h=1,2,M-1。以雷达1作为标定基准雷达,假设其方位角标定误差为0,即 1=0。误差标定后雷达m的目标观测方位角 m=m+m。各雷达接收数据处理均在距离-方位角坐标系下完成,现采用M部雷达协同的方式对 m进行估计。为了进行雷达间数据配准与融合处理,需要将各雷达的距离-方位角坐标系统一变换到雷达1的距离-方位角坐标系中1 5。变换过程中,雷达1为融合中心雷达,其余为参与融合的雷9 制 导 与 引 信 第4 4卷图1 雷达与目标位置及方位角标定误差示意图达。根据余弦定理和正弦定理,参与融合的雷达的目标观测距离Rm和观测方位角 m的量测值可表示为Rm=m-1h=1dh 2+R21-2m-1h=1dhR1c o s 1 m=-a r c s i nR1s i n 1Rm ,m=2,3,M(1)统一各雷达数据处理的坐标系可使用二维插值1 6方法实现。在低信噪比情况下,采用D P-T B D方法估计目标在各雷达坐标系下的方位角。统一坐标系后,对多部雷达进行协同非相参积累,对目标进行检测,同时可求得方位角标定误差。另外需要注意的是,在本文中,各雷达是针对单目标进行检测和跟踪并估计各雷达的方位角标定误差的,如图1所示。各雷达方位角量测误差是随机的,而各雷达方位角标定误差固定但未知,此时无法区分方位角量测误差和方位角标定误差。从式(1)可以看出,组网雷达间基线长度的量测误差也将影响到统一坐标系后各雷达的方位角量测值和距离量测值。同时,由于各雷达方位角标定误差是固定的,可在方位角标定误差的估计过程中考虑距离量测误差和方位角量测误差,而不影响方位角标定误差的估计。为了突出重点,本文假设每部雷达均已知其他雷达在各自坐标系中的准确位置。1.2 D P-T B D处理 D P-T B D是一种多帧帧间非相参积累方法,能够最大程度地保留目标的状态信息。D P-T B D能够对目标的方位角进行估计,因此D P-T B D方法能够直接应用到方位角标定误差的估计中,同时能对目标进行检测9-1 0。在建立的方位角误差模型和观测模型的基础上,D P-T B D处理通过初始化、迭代递推累积和航迹回溯等步骤,得到目标的距离-方位角航迹估计序列,最后选择最接近目标真实轨迹的序列1 7。设雷达共录取到K帧量测数据,则第k帧量测数据中,各量测单元的量测值zk=w(k),无目标A(k)+w(k),有目标,k=1,2,K(2)式中:w(k)为量测噪声幅度;A(k)为目标信号幅度。本文根据雷达方程计算各雷达接收的目标信号幅度。D P-T B D处理具体包括4个步骤。步骤1:初始化。首先建立距离-方位角状态空间,状态空间离散化后的第k帧目标状态Sk=rk k T,其中rk,k分别表示第k帧目标的距离和方位角,T为矩阵转置运算符。雷达第1次扫描(k=1)的目标状态S1=r1 1 T,有IS1 =z1S1 (S1)=0(3)式中:IS1 表示目标状态S1对应的阶段值;z1S1 表 示 目 标 状 态S1对 应 的 量 测 值;(S1)表示转移到当前帧状态S1的上一帧的状态值,因当前帧为初始帧,因此将其设为0。本文采用幅度作为D P-T B D处理的值函数。步骤2:迭代递推累积。当2kK时,对目标状态Sk=rk k T,有ISk =m a xSk-1 Sk (ISk-1 )+zkSk (Sk)=a r g m a xSk-1(ISk-1 )(4)式中:ISk 表示目标状态Sk对应的阶段值;Sk 表示第k-1帧所有的可能转移到状态Sk的离散状态的集合;zkSk 表示目标状态Sk对应的量测值;(Sk)表示转移到当前帧目标状态Sk的上一帧状态值;m a x()为求最大值函数;a r g m a x()为 求 最 大 值 对 应 的 自 变 量 的01第3期吴熙芃,等:基于D P-T B D的协同组网雷达方位标定误差估计与目标检测函数。步骤3:迭代结束。假设当迭代累积到第K帧时,能量累积的最大值超过检测门限,有Sk=a r g m a xSK(ISk ),s.t.m a x(ISk )VD T,k=K-1,K-2,1(5)式中:Sk表示第k帧中能量累积的最大值超过检测门限的状态;VD T表示D P-T B D的检测门限。步骤4:航迹回溯。对每个状态Sk,有Sk=(Sk)(6)则估计目标的航迹矩阵S1:K=S1 S2 SK。1.3 协同雷达方位角标定误差估计与目标检测算法 设 雷 达m得 到 的 全 部K帧 量 测 数 据 为Zm,1:K,D P-T B D处理后的全部K帧值函 数为Im,1:K。各雷达对同一个目标进行观测,各雷达除方位角标定误差和排布位置外,所有参数均相同。为了将各雷达量测的初始距离门对准,在探测空域已知的情况下,各雷达可在单目标所处的大致距离范围内搜索,并对单目标进行检测,确定单目标在各雷达坐标系下所处的初始距离门,最后将初始距离门对准。由于雷达n(n=2,3,M)与雷达1的差异仅体现在方位角标定误差 n上,而且各雷达的位置已知。可预先设置雷达n与雷达1之间的方位角偏移网格数,组成向量Gn(l),l=1,2,L,L为向量长度。因此,估计 n等价于从Gn(l)中找到最优的方位角偏移网格数,并且此时已进行了准确的坐标转换。Gn(l)的设置可采用如下方法进行。首先提取雷达1全部K帧量测数据Z1,1:K 的最大值所在的 方 位 角 网 格 位 置,组 成 向 量1=1,11,21,k1,K ;然后将雷达n的全部K帧量测数据Zn,1:K按照1.1节方法统一到雷达1的坐标系中得到Z-n,1:K。再提取Z-n,1:K的最大值所 在 的 方 位 角 网 格 位 置,组 成 向 量-n=-n,1-n,2-n,k-n,K ;最后计算雷达n与雷达1量测数据最大值所在方位角的偏移网格数向量-=-n-1,提取-中元素的最大值-m a x。则Gn(l)=1,-m a x+C L,其中:1,-m a x+C L表示生成1-m a x+C 之间 的L个点,这些点的间距为-m a x+C-1 /(L-1);表示向上取整;C表示可调整的网格参数,设为整数。为了找到最优的方位角偏移网格数,可采用搜索的方法进行。将Zn,1:K偏移Gn(l)中的一个元素值,得到ZGn(l)n,1:K,并将其统一到雷达1的坐标系中得到ZGn(l)n,1:K。对ZGn(l)n,1:K按1.2节的D P-T B D方法进行处理,得到值函数向量IGn(l)n,1:K。最后将IGn(l)n,1:K累加 到I1,1:K上,得 到 第l次 搜 索 的 值函数IGn(l)n,a l l=IGn(l)n,1:K+I1,1:K(7)在遍历搜索L个点之后,将得到的值函数组成向量IGn(1:L)n,a l l=IGn(1)n,a l lIGn(2)n,a l lIGn(l)n,a l lIGn(L)n,a l l (8)通过获取最优的方位角偏移网格数l,判定此时已准确进行了坐标转换。采用准则为l=a r g m a xl(IGn(1:L)n,a l l)(9)在搜索得到l后,可得到各雷达方位角标定误差 n=Gnl (1 0)式中:为方位角分辨率。协同雷达方位角标定误差估计与目标检测算法流程如下。步骤1:初始化。采集雷达n的全部K帧量测数据Zn,1:K,令l=1。步骤2:量测数据偏移方位角网格。将雷达n的全部K帧量测数据Zn,1:K偏移Gn(l),得到ZGn(l)n,1:K。步骤3:统一坐标系。将雷达n的全部K帧量测数据ZGn(l)n,1:K统一到雷达1的坐标系中,得到ZGn(l)n,1:K。11 制 导 与 引 信 第4 4卷步骤4:D P-T B D处理。对ZGn(l)n,1:K进行D P-T B D处理,得到值函数IGn(l)n,1:K,并按式(7)累加到雷达1的D P-T B D处理后的值函数上。步骤5:循环。令l=l+1,重复步骤2步骤4,直至l=L。步骤6:提取最优的偏移方位角网格数,判定已准确进行坐标转换。在循环终止之 后,得到IGn(1:L)n,a l l,并按式(9)取Gnl ,认为此时雷达n的值函数与雷达1重合,判定已准确进行坐标转换,再按式(1 0)进行各雷达方位角标定误差的计算。经由上述处理之后,可将各雷达的值 函数IGn(ln)n,1:K,其中ln表示雷达n的最优方位角偏移网格数,全部累加到雷达1的D P-T B D处理后的值函数上,并结合式(5),利用该值函数对目标进行检测。协同雷达方位角标定误差估计的处理流程如图2所示。图2 协同雷达方位角标定误差估计处理流程图2 实验验证 采用仿真数据验证算法的有效性。实验设置4部地基雷达,其中雷达1为基准雷达。考虑到工程应用中,L过大将导致搜索量的急剧增加,设置Gn(l)的向量长度L的网格数为1 0个,可调整参数C的网格数为2 0个。各雷达除方位角标定误差和排布位置外,其余参数均相同,各雷达波门开启时刻对应的距离均为3 0k m,录取数据格式为距离-方位角,采用D B F方法测角。主要仿真参数:信噪比为7 d B(高斯噪声背景下),各雷达均积累4帧数据(帧间间隔 时 间 为1s),过 程 噪 声 功 率 谱 密 度 为0.0 0 1m2/s4;雷达中心频率为3 0 GH z,带宽为1 0 0MH z,脉宽为4 s,距离采样率为1 2 0MH z,方位角分辨率为0.0 0 1 r a d。雷达在x-y平面上排布,设定各雷达的位置坐标和方位角标定误差。雷达1的位置坐标为(0 m,0 m),方位角标定误差为0 r a d;雷达2的位置坐标为(4 0 m,0 m),方位角标定误差为0.0 1 75 r a d;雷达3的位置坐标为(8 0 m,0 m),方位角标定误差为0.0 3 49 r a d;雷达4的位置坐标为(1 2 0 m,0 m),方位角标定误差为0.0 5 24 r a d。设定目标的初始状态,即目标在x-y平面上的初始位置坐标x0,y0和速度分量vx,vy分别为3 0 0,3 0 0m和1 4.3,1 6.5m/s。目标在x-y平面上匀速运动。雷达2雷达4的距离-方位角坐标系将均统一到基准雷达1的距离-方位角坐标系中。2.1 方位角标定误差估计 根据1.3节所提算法,估计雷达2雷达4的方位角标定误差。以雷达4为例,将雷达4的D P-T B D处理后的值函数累加到雷达1的D P-T B D处理后的值函数上,对比直接进行坐标转换和方位角误差校准后进行坐标转换(已取得最优的偏移方位角网格数,第8个网格)两种情况下的方位角标定误差,如图3所示。当直接进行坐标转换时,雷达1和雷达4的目标方位角量测值之间存在偏差,如图3(a)所示。为了能够利用D P-T B D同时对目标进行检测,采用搜索的方法继续进行坐标转换,直到找到最优的方位角偏移网格数,使雷达1和雷达4的目标观测方位角重合。在搜索到第8个网格时,雷达4找到最优的方位角偏移网格数,认为此时已准确进行坐标转换,如图3(b)所示。在准确进行坐标转换之后,本文算法估计所得的方位角标定误差如表1所示。可以看到,本文算法估计所得的方位角标定误差的精度与所取的方位角偏移网格数向量Gn(l)有关。由于网格数只能取整数,最后雷达2估计所得的方位角21第3期吴熙芃,等:基于D P-T B D的协同组网雷达方位标定误差估计与目标检测图3 方位角误差校准前后坐标转换的对比图 标定误差为0.0 1 7 r a d(对应网格数1 7,第5个网格),雷 达3估 计 所 得 的 方 位 角 标 定 误 差 为0.0 3 5 r a d(对应网格数3 5,第7个网格),雷达4估计所得的方位角标定误差为0.0 5 0 r a d(对应网格数5 0,第8个网格)。为了进一步对本文算法的估计精度进行分析,对雷达2雷达4在不同信噪比下各进行1 0 0次蒙特卡罗实验。设定目标在x-y平面上的初始位置坐标x0,y0和速度分量vx,vy分别为3 0 0,3 0 0m 和1 4.3,1 6.5m/s。雷达2雷达4的方位角标定误差估计的均方根误差(RM S E)与信噪比的关系如图4所示。可以看出,采用本文算法 估 计 得 到 的 各 雷 达 方 位 角 标 定 误 差 的RM S E随着信噪比的提升变化并不明显。在目标运动过程中,各雷达量测方位角随着雷达-目标距离的变化而变化。为了分析不同的目标运动参数对方位角标定误差估计的影响,设 表1 雷达方位角标定误差的估计雷达真实方位角标定误差/r a d本文算法估计方位角标定误差/r a d方位角偏移网格数向量绝对误差/r a d相对误差/%雷达20.0 1 7 50.0 1 71,5,9,1 3,1 7,2 1,2 5,2 9,3 3,3 70.0 0 0 52.6 0雷达30.0 3 4 90.0 3 51,7,1 3,1 8,2 4,2 9,3 5,4 0,4 6,5 10.0 0 0 10.2 7雷达40.0 5 2 40.0 5 01,8,1 5,2 2,2 9,3 6,4 3,5 0,5 7,6 40.0 0 2 44.5 1图4 雷达方位角标定误差估计的均方根误差1定目标在x-y平面上的初始位置坐标x0,y0和速度分量vx,vy分别为3 0 0,3 0m,和1,1m/s,各雷达参数 及 处 理 流 程 与 前 文 相 同。对 雷 达2雷达4在不同信噪比下各进行1 0 0次蒙特卡罗实验。雷达2雷 达4方 位 角 标 定 误 差 估 计 的RM S E与信噪比的关系如图5所示。从图5可以看出,在该目标运动参数下,通过本文算法 估计 得 到 的 各 雷 达 方 位 角 标 定 误差的RMS E随 着 信 噪 比 的 提 升,变 化 并 不明显。对比图4和图5可以看出,在两种不同的目图5 雷达方位角标定误差估计的均方根误差2标运动参数下,各雷达方位角标定误差的RM S E并不相同,说明目标运动参数的改变将导致各雷达方位角标定误差RM S E的改变。此外,在图4中,雷达4的RM S E明显大于雷达2和雷达3的;在图5中,雷达3的RM S E明显大于雷达2和雷达4的。这是由于受到噪声影响,各雷达依据本文方法所获取的方位角偏移网格数在每次蒙特卡罗实验中不尽相同,使得各雷达估计的方位角标定误差在每次蒙特卡罗实验中不尽相同,最终导致各雷达的方位角标定误差的RM S E不尽相同。31 制 导 与 引 信 第4 4卷2.2 检测与跟踪性能分析 D P-T B D是一种帧间非相参积累方法。其性能分析采用蒙特卡罗实验进行:在高斯白噪声背景下,当虚警概率为0.0 0 0 1时,进行2 0 0 0 0次蒙特卡罗实验,对1个目标进行观测。检测性能分析共分3种情况进行。情况1:单雷达节点多帧非相参积累;情况2:不存在方位角标定误差时多雷达节点非相参积累;情况3:存在方位角标定误差时多雷达节点非相参积累。不同情况下的检测概率如图6所示。图6 不同情况下的检测概率从图6可以看出,相比于单雷达多帧积累的情况,对校正标定误差之后的多节点组网雷达进行多帧积累,可显著提高低信噪比条件下的目标检测性能,而当存在标定误差时,组网雷达多帧积累的检测性能几乎没有改善。跟踪性能分析在信噪比为7 d B的条件下进行。首先对比方位角标定误差校正前各雷达恢复的目标航迹和误差校正后多雷达融合的目标航迹,如图7所示。情况1:理想目标航迹;情况2:基准雷达1恢复航迹;情况3:校正前雷达2恢复航迹;情况4:校正前雷达3恢复航迹;情况5:校正前雷达4恢复航迹;情况6:多雷达融合目标航迹。可 以 看 到,在 方 位 角 标 定 误 差 校 正 前,雷达2雷达4各自恢复的3条航迹均与理想目标航迹存在差距,主要是在方位角上存在差距;在方位角标定误差校正后,将各雷达航迹融合,与理想目标航迹进行对比,可以看到融合后的航迹与理想目标航迹具有良好的一致性,尤其是在方位角上并没有明显的差距。图7中的目标跟踪误差(绝对误差)均值,如表2所示。可以看到,多雷达融合后的目标跟踪误差均值比雷达2雷达4校正前的目标跟踪误图7 误差校正前后恢复的航迹差均值小。表2 目标跟踪误差均值雷达距离/m方位角/r a d雷达1(基准)0.6 3 5 30.0 0 0 6雷达2(校正前)0.4 9 1 10.0 1 7 2雷达3(校正前)1.5 5 4 90.0 3 1 4雷达4(校正前)4.0 5 4 90.0 4 3 9多雷达(融合)0.4 0 4 20.0 0 1 63 结论 本文提出了一种基于D P-T B D的协同组网雷达方位角标定误差估计方法。该算法将D P-T B D和多雷达协同处理结合起来,能直接利用量测数据估计组网雷达的方位角标定误差,并进行多雷达协同积累,实现目标检测,可显著提高低信噪比条件下的目标检测性能。仿真结果验证了算法的有效性。参考文献1 欧阳志宏.雷达组网误差配准算法研究D.长沙:国防科学技术大学,2 0 0 7.2 杨宏文,郁文贤,胡卫东,等.基于数据补偿的雷达系 统 误 差 估 计 J.火 力 与 指 挥 控 制,2 0 0 0,2 5(2):2 3-2 7.3 何友,修建娟,张晶炜,等.雷达数据处理及应用M.3版.北京:电子工业出版社,2 0 1 3:3 3 5-3 4 6.4 王波,王灿林,董云龙.R T Q C误差配准算法性能分析J.系统 仿真 学 报,2 0 0 6,1 8(1 1):3 0 6 7-3 0 8 1.5 李鸿艳,冯新喜.一种基于E C E F坐标转换的最小二乘 配 准 算 法 J.系 统 工 程 与 电 子 技 术,2 0 0 2,2 4(1):9 2-9 5.6 拜兵虎,白树林,张玉梅.广义最小二乘法配准组41第3期吴熙芃,等:基于D P-T B D的协同组网雷达方位标定误差估计与目标检测网雷达 误 差 的 技 术 研 究 J.无 线 互 联 科 技,2 0 1 4(7):1 5 8-1 5 9.7 Z HOU Y F,L E UN G H,Y I P P C.A n e x a c t m a x i m u m l i k e l i h o o d r e g i s t r a t i o n a l g o r i t h m f o r d a t a f u s i o nJ.I E E E T r a n s a c t i o n s o n S i g n a l P r o c e s s i n g,1 9 9 7,4 5(6):1 5 6 0-1 5 7 3.8 Z HOU Y F,L E UNG H,B L AN CHE T T E M.S e n s o r a l i g n m e n t w i t h e a r t h-c e n t e r e d e a r t h-f i x e d(E C E F)c o o r d i n a t e s y s t e mJ.I E E E T r a n s a c t i o n s o n A e r o s p a c e a n d E l e c t r o n i c S y s t e m s,1 9 9 9,3 5(2):4 1 0-4 1 7.9 易伟.基于检测前跟踪技术的多目标跟踪算法研究D.成都:电子科技大学,2 0 1 3.1 0 B A R N I V Y.D y n a m i c p r o g r a mm i n g s o l u t i o n f o r d e t e c t i n g d i m m o v i n g t a r g e t sJ.I E E E T r a n s-a c t i o n s o n A e r o s p a c e a n d E l e c t r o n i c S y s t e m s,1 9 8 5,A E S-2 1(1):1 4 4-1 5 6.1 1 L E D,HUANG Q,L I U Z L,e t a l.A d a p t i v e-t r a n s i t i o n-s e t t r a c k b e f o r e d e t e c t a l g o r i t h m b a s e d o n d y n a m i c p r o g r a mm i n gC/2 0 1 9 I n t e r n a t i o n a l C o n f e r e n c e o n C o n t r o l,A u t o m a t i o n a n d I n f o r m a t i o n S c i e n c e s(I C C A I S),C h e n g d u,C h i n a.P i s c a t a w a y,N J:I E E E P r e s s,2 0 2 0:1-5.1 2 YU H B,TAN S C,C AO Q,e t a l.A d y n a m i c p r o g r a mm i n g b a s e d T B D a l g o r i t h m f o r n e a r s p a c e t a r g e t s u n d e r r a n g e a m b i g u i t yC/C o mm u n i-c a t i o n s,S i g n a l P r o c e s s i n g a n d S y s t e m,S i n g a p o r e.B e r l i n:S p r i n g e r-V e r l a g,2 0 2 0:1 0 2 9-1 0 3 7.1 3 L I Y,L I C P,T I AN M,e t a l.T w o-s t e p t h r e s h o l d s T B D a l g o r i t h m f o r t i m e s e n s i t i v e t a r g e t b a s e d o n d y n a m i c p r o g r a mm i n gJ.I E E E A c c e s s,2 0 2 0,8:2 0 9 2 6 7-2 0 9 2 7 7.1 4 李昕哲,王首勇,郑岱堃.一种杂波自适应动态规划检测 前 跟 踪 算 法 J.空 军 预 警 学 院 学 报,2 0 1 6,3 0(3):1 6 2-1 6 6,1 7 2.1 5 U L AN D E R L,HE L L S T E N H,S T E N S T R OM G.S y n t h e t i c-a p e r t u r e r a d a r p r o c e s s i n g u s i n g f a s t f a c t o r i z e d b a c k-p r o j e c t i o nJ.I E E E T r a n s a c t i o n s o n A e r o s p a c e a n d E l e c t r o n i c S y s t e m s,2 0 0 3,3 9(3):7 6 0-7 7 6.1 6 F R O L I N D P O,U L AN D E R L.E v a l u a t i o n o f a n g u l a r i n t e r p o l a t i o n k e r n e l s i n f a s t b a c k-p r o j e c t i o n S A R p r o c e s s i n gJ.I E E P r o c e e d i n g s-R a d a r,S o n a r a n d N a v i g a t i o n,2 0 0 6,1 5 3(3):2 4 3-2 4 9.1 7 文鸣.中重频雷 达 的检 测前 跟踪 技 术研 究D.成都:电子科技大学,2 0 1 9.(上接第7页)3 X I AO B,L I U Y X,HUO K,e t a l.O F DM i n t e g r a t i o n w a v e f o r m d e s i g n f o r r a d a r a n d c o mm u n i c a t i o n a p p l i c a t i o nC/I n t e r n a t i o n a l C o n f e r e n c e o n R a d a r S y s t e m s(R a d a r 2 0 1 7),B e l f a s t.P i s c a t a w a y,N J:I E E E P r e s s,2 0 1 7:1-5.4 胡苏,高 原,王 军.通 信 雷 达 一 体 化 波 形 设 计M.北京:国防工业出版社,2 0 1 8:5-1 3.5 Z HAN G W X,Z HANG H.T h e d e s i g n o f i n t e g r a t e d w a v e f o r m b a s e d o n M S K-L FM s i g n a lC/2 0 2 0 1 5 t h I E E E I n t e r n a t i o n a l C o n f e r e n c e o n S i g n a l P r o c e s s i n g(I C S P),B e i j i n