一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置

1.本发明涉及一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置，属于雷达技术领域。


背景技术：

2.在实际的工程应用中，由于存在系统误差、量测误差和时间误差，很难满足空间完全对齐的条件，特别是机载雷达等动平台处于快速机动或机动转弯的场景下。针对非随机的静态对齐误差的消除，前人已经提出了大量的对齐算法，主要分为量测级对齐算法、跟踪级对齐算法。对于随机对齐误差，上述方法并不能有效消除。因此，即使完成了误差校正，在实际工程应用中仍会有较小的对齐误差存在。
3.在对齐误差存在的情况下，基于真实目标空间位置相关，而虚假目标空间位置不相关的特点，提出了数据级融合抗欺骗式干扰算法。该算法在单个真实目标情况下，可较好实现真实目标的鉴别。但是，多站雷达系统被欺骗概率受欺骗距离影响严重，且随着欺骗距离减小，被欺骗概率不断增大。


技术实现要素：

4.本发明的目的在于克服现有技术中的不足，提供一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置，可以在多站雷达系统中存在对齐误差的场景下完成有源假目标鉴别。
5.为达到上述目的，本发明是采用下述技术方案实现的：
6.第一方面，本发明提供了一种对齐误差下多站雷达协同抗欺骗式干扰方法，包括：
7.根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；
8.根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；
9.根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；
10.根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量；
11.根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果。
12.进一步地，所述根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数，包括：
13.获取多站雷达的站址误差向量，计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量求得目标对齐误差向量的误差协方差矩阵；
14.根据误差协方差矩阵，得到测量目标距离的标准差和测量目标角度的标准差，根
据测量目标距离的标准差和测量目标角度的标准差求得误差椭圆的距离维长度和角度维长度；
15.根据距离维长度和角度维长度，计算得到所有空间分集通道的距离分辨单元数和角度分辨单元数。
16.进一步地，所述获取多站雷达的站址误差向量，计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量求得目标对齐误差向量的误差协方差矩阵，包括：
17.设置一个二维坐标系，第m个发射机的坐标位置为(x
t
,y
t
)
t
，第n个接收机的坐标位置为(xr,yr)
t
，目标在接收机中的量测值为z＝(ρr,θr)
t
，则根据其量测值可得到目标的位置坐标为：
[0018][0019]
其中，rr为从目标到接收机的距离，(
·
)
t
表示矩阵转置运算，ρr表示发射机到目标再到接收机的距离，θr表示目标相对于接收机的角度，下标t表示发射站，下标r表示接收站；
[0020]
对公式(1)两侧微分，由站址误差向量(dxr,dyr,dx
t
,dy
t
)
t
得到目标对齐误差向量(dρr,dθr)
t
的公式为：
[0021][0022]
其中，dx
t
和dy
t
分别表示发射站在坐标系上沿x轴和y轴的误差向量，dxr和dyr分别表示接收站在坐标系上沿x轴和y轴的误差向量，dρr表示接收机得到的发射机到目标再到接收机的距离和的误差向量，dθr表示目标相对于接收机的角度的误差向量，c
r1
＝cosθr，c
r2
＝sinθr，c
t1
＝cosθ
t
，c
t2
＝sinθ
t
，θ
t
表示目标相对于发射机的角度，rr和r
t
是从目标分别到接收机和发射机的距离，矩阵b是双基探测下由站址误差到目标对齐误差的转换矩阵；
[0023]
根据矩阵b，求目标对齐误差向量(dρr,dθr)
t
的误差协方差矩阵p，公式为：
[0024]
p＝bλb
t
(3)
[0025]
其中，λ是发射站和接收站的站址误差矢量(dxr,dyr,dx
t
,dy
t
)
t
的误差协方差矩阵，阵，和分别表示接收站在x轴和y轴上的量测方差，和分别表示发射站在x轴和y轴上的量测方差，diag(
·
)表示以向量每一个元素为对角元素构成的对角矩阵。
[0026]
进一步地，所述根据误差协方差矩阵，得到测量目标距离的标准差和测量目标角度的标准差，根据测量目标距离的标准差和测量目标角度的标准差求得误差椭圆的距离维
长度和角度维长度，包括：
[0027]
根据误差协方差矩阵p，得到测量目标距离的标准差和测量目标角度的标准差
[0028][0029][0030]
其中，[
·
]
ij
表示取矩阵的第i行第j列数据；
[0031]
设目标存在于目标对齐误差向量的误差椭圆的概率为α，计算出该误差椭圆的距离维长度a和角度维长度b，离维长度a和角度维长度b，
[0032]
进一步地，所述根据距离维长度和角度维长度，计算得到所有空间分集通道的距离分辨单元数和角度分辨单元数，包括：
[0033]
根据距离维长度a和角度维长度b，得到第mn个空间分集通道中目标存在的误差椭圆中的距离分辨单元数q
ρ
(m,n)和角度分辨单元数q
θ
(m,n)，分别为：
[0034][0035][0036]
其中，表示舍入到大于或等于它的最接近整数，δr及δθ分别是接收站的距离和角度分辨率，第mn个空间分集通道表示由第m个发射站发射、第n个接收站接收构成的空间分集通道；
[0037]
根据第mn个空间分集通道中目标存在的误差椭圆中的距离分辨单元数q
ρ
(m,n)和角度分辨单元数q
θ
(m,n)，计算得到多站雷达中所有mn个空间分集通道对应的对齐误差椭圆的距离分辨单元数q
ρ
(m,n)和角度分辨单元数q
θ
(m,n)，其中，m、n分别为多站雷达系统中发射站和接收站的数量，1≤m≤m，1≤n≤n。
[0038]
进一步地，所述根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量，包括：
[0039]
根据每个空间分集通道的q
ρ
(m,n)和q
θ
(m,n)，得到目标在第mn个空间分集通道中的距离存在范围d(m,n)和角度存在范围分别为：
[0040][0041][0042]
根据所述距离存在范围d(m,n)和角度存在范围得到目标在第mn个空间分
集通道中的距离-角度二维存在范围为：
[0043][0044]
取出目标在第mn个空间分集通道的距离-角度二维存在范围内的所有分辨率单元的回波信号构成信号矢量y
mn
，公式如下：
[0045]ymn
＝[y
mn
(1),y
mn
(2),
…
,y
mn
(q
mn
)]，1≤m≤m，1≤n≤n
ꢀꢀꢀꢀꢀꢀ
(11)
[0046]
其中，y
mn
(p)是第p个距离-角度分辨率单元的回波信号，p＝1，2
……qmn
,q
mn
＝q
ρ
(m,n)q
θ
(m,n)是所有距离-角度分辨率单元的个数；
[0047]
根据最大似然算法，计算在第mn个空间分集通道内目标回波的估计值y
mn
为y
mn
＝max(y
mn
)，1≤m≤m，1≤n≤n；
[0048]
根据得到的各空间分集通道内目标回波的估计值y
mn
，得到目标在多站雷达下的接收信号矢量x，
[0049]
x＝[y
11
,y
12
,
…
,y
mn
]
ꢀꢀꢀꢀꢀꢀꢀ
(12)
[0050]
其中，y
mn
第mn个空间分集通道内目标回波的估计值。
[0051]
进一步地，所述根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数，包括：
[0052]
检测区域中检测到k个目标，接收信号矢量记为xk,k＝1,2,...,k，k为检测到目标个数，计算不同目标之间的相关系数，形成相关系数矩阵ω，其第i行第j列的元素[ω]
ij
表示为：
[0053][0054]
式中，(
·
)h表示共轭转置，||
·
||表示向量的欧几里得范数，xi表示第i行的接收信号矢量，xj表示第j行的接收信号矢量。
[0055]
进一步地，所述根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量，包括：
[0056]
设置门限η，对相关系数矩阵ω进行二进制量化，得到量化后的相关系数矩阵
[0057][0058]
其中，是对角元素全为1的对称矩阵，所有元素非零即一，表示：第i个目标和第j个目标的接收信号矢量是不相关的，两个目标记为真实目标1次；表示：第i个目标和第j个目标的接收信号矢量是相关的，两个目标记为虚假目标1次；
[0059]
根据量化后相关系数矩阵计算目标被记录为假目标的次数，形成目标相关次数矢量φ，其第i个元素[φ]i为：
[0060][0061]
其中，k为检测到目标个数。
[0062]
进一步地，所述根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果，包括：
[0063]
设置目标判别阈值κ，根据目标相关次数矢量φ，得到第i个目标的真假鉴别结果[φ]i，目标的真假鉴别结果表示为：
[0064][0065]
第二方面，本发明提供一种对齐误差下多站雷达协同抗欺骗式干扰装置，包括：
[0066]
分辨单元数确定模块，用于根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；
[0067]
接收信号矢量估计模块，用于根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；
[0068]
相关系数计算模块，用于根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；
[0069]
目标相关次数矢量计算模块，用于根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量；
[0070]
鉴别结果获取模块，用于根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果。
[0071]
第三方面，本发明提供一种电子设备，包括处理器及存储介质；
[0072]
所述存储介质用于存储指令；
[0073]
所述处理器用于根据所述指令进行操作以执行根据前述任一项所述方法的步骤。
[0074]
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现前述任一项所述方法的步骤。
[0075]
与现有技术相比，本发明所达到的有益效果：
[0076]
1、本发明提供一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置，适用于随机对齐误差无法有效消除的场景下，可以在多站雷达系统中对齐误差无法有效的场景下消除完成有源假目标鉴别。
[0077]
2、本发明提供一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置，利用了真假目标接收信号矢量相关性的差异进行目标鉴别，与欺骗干扰类型无关，可适用于任意调制产生的有源假目标。
附图说明
[0078]
图1是本发明的实现流程图；
[0079]
图2是对齐误差随站址误差的变化示意图；
[0080]
图3在q
mn
(m,n)＝3时有源目标鉴别算法的真实目标鉴别概率示意图；
[0081]
图4在q
mn
(m,n)＝3时有源目标鉴别算法的有源假目标误判概率示意图；
[0082]
图5在q
mn
(m,n)＝5时有源目标鉴别算法的真实目标鉴别概率示意图；
[0083]
图6在q
mn
(m,n)＝5时有源目标鉴别算法的有源假目标误判概率示意图；
[0084]
图7在q
mn
(m,n)＝7时有源目标鉴别算法的真实目标鉴别概率示意图；
[0085]
图8在q
mn
(m,n)＝7时有源目标鉴别算法的有源假目标误判概率示意图。
[0086]
实施方式
[0087]
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0088]
本实施例介绍一种对齐误差下多站雷达协同抗欺骗式干扰方法，包括：
[0089]
根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；
[0090]
根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；
[0091]
根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；
[0092]
根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量；
[0093]
根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果。
[0094]
如图1所示，本实施例提供的对齐误差下多站雷达协同抗欺骗式干扰方法，其应用过程具体涉及如下步骤：
[0095]
步骤1、根据各雷达的站址误差计算得到多站雷达的对齐误差向量，确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数。
[0096]
设多站雷达系统由m个发射站和n个接收站组成，共可构成mn个空间分集通道。以第m个发射机和第n个接收机组成的空间分集通道(记为第mn个空间分集通道)为例说明该通道目标存在的误差椭圆的距离维分辨单元数q
ρ
(m,n)和角度维分辨单元数q
θ
(m,n)的计算方法。
[0097]
1a)假设有一个二维坐标系，第m个发射机坐标位置为(x
t
,y
t
)
t
，第n个接收机坐标位置为(xr,yr)
t
，目标在接收机中的量测值为z＝(ρr,θr)
t
，则根据其量测值可得到目标的位置坐标为：
[0098][0099]
其中，rr为从目标到接收机的距离，(
·
)
t
表示矩阵转置运算，ρr表示发射机到目标再到接收机的距离，θr表示目标相对于接收机的角度，下标t表示发射站，下标r表示接收站；
[0100]
1b)对公式两侧微分，由站址误差向量(dxr,dyr,dx
t
,dy
t
)
t
得到目标对齐误差向量(dρr,dθr)
t
的公式为：
[0101][0102]
其中，dx
t
和dy
t
分别表示发射站在坐标系上沿x轴和y轴的误差向量，dxr和dyr分别表示接收站在坐标系上沿x轴和y轴的误差向量，dρr表示接收机得到的发射机到目标再到
接收机的距离和的误差向量，dθr表示目标相对于接收机的角度的误差向量，c
r1
＝cosθr，c
r2
＝sinθr，c
t1
＝cosθ
t
，c
t2
＝sinθ
t
，θ
t
表示目标相对于发射机的角度，rr和r
t
是从目标分别到接收机和发射机的距离，矩阵b是双基探测下由站址误差到目标对齐误差的转换矩阵；
[0103]
1c)根据1b)中的矩阵b，求目标对齐误差向量(dρr,dθr)
t
的误差协方差矩阵p，公式为：
[0104]
p＝bλb
t
(3)
[0105]
其中，λ是发射站和接收站的站址误差矢量(dxr,dyr,dx
t
,dy
t
)
t
的误差协方差矩阵，阵，和分别表示接收站在x轴和y轴上的量测方差，和分别表示发射站在x轴和y轴上的量测方差，diag(
·
)表示以向量每一个元素为对角元素构成的对角矩阵。
[0106]
1d)根据1c)中的误差协方差矩阵p，得到测量目标距离的标准差和测量目标角度的标准差
[0107][0108][0109]
其中，[
·
]
ij
表示取矩阵的第i行第j列数据。
[0110]
1e)设目标存在于目标对齐误差向量的误差椭圆的概率为α，计算出该误差椭圆的两轴长度，距离维长度角度维长度
[0111]
1f)根据步骤1e)中计算得到的距离维长度a和角度维长度b，可以得到第mn个空间分集通道中目标存在的误差椭圆中的距离分辨单元数q
ρ
(m,n)和角度分辨单元数q
θ
(m,n)分别为：
[0112][0113][0114]
其中，表示舍入到大于或等于它的最接近整数，δr及δθ分别是接收站的距离和角度分辨率，第mn个空间分集通道表示由第m个发射站发射、第n个接收站接收构成的空间分集通道。
[0115]
重复上述步骤可得到多站雷达中所有mn个空间分集通道对应的对齐误差椭圆的
距离分辨单元数q
ρ
(m,n)和角度分辨单元数q
θ
(m,n)，其中，1≤m≤m，1≤n≤n。
[0116]
步骤2、根据目标各通道误差椭圆内的回波接收数据，采用最大似然(ml)算法从误差椭圆中得到所有检测到目标的接收信号矢量。
[0117]
具体实现步骤如下：
[0118]
2a)根据每个空间分集通道的q
ρ
(m,n)和q
θ
(m,n)，可得到目标在第mn个空间分集通道中的距离存在范围d(m,n)和角度存在范围分别为
[0119][0120][0121]
2b)根据步骤2a)的距离存在范围d(m,n)和角度存在范围得到目标在第mn个空间分集通道中的距离-角度二维存在范围为
[0122][0123]
2c)取出目标在第mn个空间分集通道的距离-角度二维存在范围内的所有分辨率单元的回波信号构成信号矢量y
mn
，
[0124]ymn
＝[y
mn
(1),y
mn
(2),
…
,y
mn
(q
mn
)]，1≤m≤m，1≤n≤n
ꢀꢀ
(11)
[0125]
其中，y
mn
(p)是第p个距离-角度分辨率单元的回波信号，p＝1，2
……qmn
,q
mn
＝q
ρ
(m,n)q
θ
(m,n)是所有距离-角度分辨率单元的个数。
[0126]
2d)根据最大似然算法，计算在第mn个空间分集通道内目标回波的估计值y
mn
为y
mn
＝max(y
mn
)，1≤m≤m，1≤n≤n。
[0127]
2e)根据2d)得到的各空间分集通道内目标回波的估计值y
mn
，得到目标在多站雷达下的接收信号矢量x，
[0128]
x＝[y
11
,y
12
,
…
,y
mn
]
ꢀꢀ
(12)
[0129]
其中，y
mn
第mn个空间分集通道内目标回波的估计值。
[0130]
步骤3、根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数。
[0131]
具体实现步骤如下：
[0132]
检测区域中检测到k个目标，接收信号矢量记为xk,k＝1,2,...,k，k为检测到目标个数，计算不同目标之间的相关系数，形成相关系数矩阵ω(维数为k
×
k矩阵)，其第i行第j列的元素为
[0133][0134]
式中，(
·
)h表示共轭转置，||
·
||表示向量的欧几里得范数，xi表示第i行的接收信号矢量，xj表示第j行的接收信号矢量。
[0135]
步骤4、根据目标接收信号矢量之间的相关系数，统计各目标在另一雷达中与之相关的目标个数，设计阈值完成假目标鉴别。
[0136]
具体实现步骤如下：
[0137]
4a)设置门限η，对相关系数矩阵ω进行二进制量化，得到量化后的相关系数矩阵
[0138][0139]
在没有任何先验知识的情况下，门限η＝0.5。是对角元素全为1的对称矩阵，所有元素非零即一，表示第i个目标和第j个目标的接收信号矢量是不相关的，两个目标记为真实目标1次；表示第i个目标和第j个目标的接收信号矢量是相关的，两个目标记为虚假目标1次。
[0140]
4b)根据量化后相关系数矩阵计算目标被记录为假目标的次数，形成目标相关次数矢量φ，其第i个元素[φ]i为
[0141][0142]
其中，k为检测到目标个数。
[0143]
如果第i个目标是真实目标，[φ]i是一个小整数，特别是在独立的情况下，[φ]i等于零，如果第i个目标是假目标，[φ]i是一个更大的整数。
[0144]
4c)设置目标判别阈值κ，根据目标相关次数矢量φ，得到第i个目标的真假鉴别结果：
[0145][0146]
通常，可以设置目标判别阈值κ＝1或2。
[0147]
实验场景：
[0148]
多站雷达系统由m＝2个发射站和n＝5个接收站组成，发射站位置坐标为[0,0]m、[300,0]m。接收站位置坐标分别为[-500,0]m、[-250,0]m、[0,0]m、[250,0]m、[500,0]m，对某飞机编队所在区域进行探测。该飞机编队中有四架飞机作为雷达真实目标分别为[30,31]km、[31,30]km、[29,30]km、[30,29]km，其中两架飞机携带有源干扰机，分别为[29,30]km、[30,29]km，对多站雷达系统实施欺骗式干扰。
[0149]
假设每个通道的距离分辨率为50m，角度分辨率为2
°
。以由[0,0]m处的发射机和[500,0]m处的接收机以及[30,31]km处的目标组成的信道为例。所有发射机和接收机站的天线增益相同，系统波长为λ＝0.1m。目标存在于目标对齐误差向量的误差椭圆的概率为α＝0.99。
[0150]
为了保护飞机编队，每个干扰机一次生成8个有源假目标，所有真实目标的snr相同，在第一信道中设置为8db，所有主动假目标的jnr相同，在第一个信道中从7.5db变为20db。可以根据雷达方程计算其他分集信道中的snr或jnr。阈值η设置为η＝0.5,以及目标鉴别阈值κ设置为κ＝2。在对齐误差是3、5和7的情况下，对所提出的目标鉴别方法进行了仿真，以评估其鉴别性能。
[0151]
实验内容：
[0152]
(1)根据步骤1，得到目标所在误差椭圆的长轴短轴的距离和，得到对齐误差随站
址误差的变化图，如图2所示。
[0153]
(2)通过105次蒙特卡罗仿真实验，在不同对齐误差q
mn
(m,n)下，仿真了真实目标的鉴别概率p
t
和有源假目标的误判概率pf随jnr的变化情况，如图3-5所示，每个图中的曲线对应于目标尺寸d的不同值，其中d＝0m、15m和30m。
[0154]
实验结果：
[0155]
(1)从图2以及仿真结果可以看出，对齐误差与站址误差成比例，误差椭圆中的角度分辨单元数q
θ
＝1，在站址误差小于100m的情况下，对齐误差不大于9，误差椭圆包含9个距离-角度分辨单元数。
[0156]
(2)如图3-8所示，在所有仿真条件下，真实目标的鉴别概率始终大于95％，虚假目标的误判概率保持在0.9％以下。换言之，所提出的方法可以在对齐误差下获得预期的鉴别性能，在保留真实目标的情况下有效地鉴别假目标，这表明了该算法的可行性。如图3-图7所示，真实目标鉴别概率p
t
的曲线在10db附近有一个拐点，当jnr《10db时，鉴别概率随jnr减小。当jnr》10db时，随着jnr增大。如图4-图8所示，当jnr《10db时，对齐误差增加了虚假目标的误判概率。当jnr大于10db时，对齐误差对误判概率的影响较小，几乎为零。
[0157]
基于同一设计构思，本发明提供一种对齐误差下多站雷达协同抗欺骗式干扰装置，包括：
[0158]
分辨单元数确定模块，用于根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；
[0159]
接收信号矢量估计模块，用于根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；
[0160]
相关系数计算模块，用于根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；
[0161]
目标相关次数矢量计算模块，用于根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量；
[0162]
鉴别结果获取模块，用于根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果。
[0163]
基于同一设计构思，本发明提供一种电子设备，包括处理器及存储介质；
[0164]
所述存储介质用于存储指令；
[0165]
所述处理器用于根据所述指令进行操作以执行根据实施例1中任一项所述方法的步骤。
[0166]
基于同一设计构思，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现实施例1中任一项所述方法的步骤。
[0167]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。

技术特征：
1.一种对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，包括：根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量；根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果。2.根据权利要求1所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数，包括：获取多站雷达的站址误差向量，计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量求得目标对齐误差向量的误差协方差矩阵；根据误差协方差矩阵，得到测量目标距离的标准差和测量目标角度的标准差，根据测量目标距离的标准差和测量目标角度的标准差求得误差椭圆的距离维长度和角度维长度；根据距离维长度和角度维长度，计算得到所有空间分集通道的距离分辨单元数和角度分辨单元数。3.根据权利要求2所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述获取多站雷达的站址误差向量，计算得到多站雷达的对齐误差向量，根据多站雷达的对齐误差向量求得目标对齐误差向量的误差协方差矩阵，包括：设置一个二维坐标系，第m个发射机的坐标位置为，第n个接收机的坐标位置为，目标在接收机中的量测值为，则根据其量测值可得到目标的位置坐标为：其中，为从目标到接收机的距离，，表示矩阵转置运算，表示发射机到目标再到接收机的距离， 表示目标相对于接收机的角度，下标t表示发射站，下标r表示接收站；对公式（1）两侧微分，由站址误差向量得到目标对齐误差向量的公式为：其中，和分别表示发射站在坐标系上沿x轴和y轴的误差向量，和分别表示接收站在坐标系上沿x
轴和y轴的误差向量，表示接收机得到的发射机到目标再到接收机的距离和的误差向量，表示目标相对于接收机的角度的误差向量，，，，，表示目标相对于发射机的角度，，和是从目标分别到接收机和发射机的距离，，，矩阵b是双基探测下由站址误差到目标对齐误差的转换矩阵；根据矩阵b，求目标对齐误差向量的误差协方差矩阵p，公式为：其中，是发射站和接收站的站址误差矢量的误差协方差矩阵，，和分别表示接收站在x轴和y轴上的量测方差，和分别表示发射站在x轴和y轴上的量测方差，表示以向量每一个元素为对角元素构成的对角矩阵。4.根据权利要求3所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据误差协方差矩阵，得到测量目标距离的标准差和测量目标角度的标准差，根据测量目标距离的标准差和测量目标角度的标准差求得误差椭圆的距离维长度和角度维长度，包括：根据误差协方差矩阵p，得到测量目标距离的标准差和测量目标角度的标准差：其中，表示取矩阵的第i行第j列数据；设目标存在于目标对齐误差向量的误差椭圆的概率为，计算出该误差椭圆的距离维长度a和角度维长度b，，。5.根据权利要求4所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据距离维长度和角度维长度，计算得到所有空间分集通道的距离分辨单元数和角度分辨单元数，包括：根据距离维长度a和角度维长度b，得到第mn个空间分集通道中目标存在的误差椭圆中的距离分辨单元数和角度分辨单元数，分别为：
其中，表示舍入到大于或等于它的最接近整数，及分别是接收站的距离和角度分辨率，第mn个空间分集通道表示由第m个发射站发射、第n个接收站接收构成的空间分集通道；根据第mn个空间分集通道中目标存在的误差椭圆中的距离分辨单元数和角度分辨单元数，计算得到多站雷达中所有mn个空间分集通道对应的对齐误差椭圆的距离分辨单元数和角度分辨单元数，其中，m、n分别为多站雷达系统中发射站和接收站的数量，，。6.根据权利要求5所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量，包括：根据每个空间分集通道的和，得到目标在第mn个空间分集通道中的距离存在范围和角度存在范围分别为：根据所述距离存在范围和角度存在范围，得到目标在第mn个空间分集通道中的距离-角度二维存在范围为：取出目标在第mn个空间分集通道的距离-角度二维存在范围内的所有分辨率单元的回波信号构成信号矢量，公式如下：，其中，是第p个距离-角度分辨率单元的回波信号，,是所有距离-角度分辨率单元的个数；根据最大似然算法，计算在第mn个空间分集通道内目标回波的估计值为，，；根据得到的各空间分集通道内目标回波的估计值，得到目标在多站雷达下的接收信号矢量，其中，为第mn个空间分集通道内目标回波的估计值。7.根据权利要求6所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数，包括：检测区域中检测到k个目标，接收信号矢量记为, ，k为检测到目标个数，计算不同目标之间的相关系数，形成相关系数矩阵，其第i行第j列的元素表示为：
式中，表示共轭转置，表示向量的欧几里得范数，表示第i行的接收信号矢量，表示第j行的接收信号矢量。8.根据权利要求7所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据目标接收信号矢量之间的相关系数，计算目标被记录为假目标的次数，形成目标相关次数矢量，包括：设置门限，对相关系数矩阵进行二进制量化，得到量化后的相关系数矩阵：其中，是对角元素全为1的对称矩阵，所有元素非零即一，表示：第i个目标和第j个目标的接收信号矢量是不相关的，两个目标记为真实目标1次；表示：第i个目标和第j个目标的接收信号矢量是相关的，两个目标记为虚假目标1次；根据量化后相关系数矩阵，计算目标被记录为假目标的次数，形成目标相关次数矢量，其第i个元素为：其中，k为检测到目标的个数。9.根据权利要求8所述的对齐误差下多站雷达协同抗欺骗式干扰方法，其特征在于，所述根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果，包括：设置目标判别阈值k，根据目标相关次数矢量，得到第i个目标的真假鉴别结果，目标的真假鉴别结果表示为：。10.一种对齐误差下多站雷达协同抗欺骗式干扰装置，其特征在于，包括：分辨单元数确定模块，用于根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；接收信号矢量估计模块，用于根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；相关系数计算模块，用于根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；目标相关次数矢量计算模块，用于根据目标接收信号矢量之间的相关系数，计算目标
被记录为假目标的次数，形成目标相关次数矢量；鉴别结果获取模块，用于根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果。

技术总结
本发明公开了一种对齐误差下多站雷达协同抗欺骗式干扰方法及装置，所述方法包括根据多站雷达的站址误差计算得到多站雷达的对齐误差向量，确定目标存在的误差椭圆的距离维和角度维的两轴分辨单元数；根据误差椭圆的距离维和角度维的两轴分辨单元数,采用最大似然算法从误差椭圆中得到所有检测到目标的接收信号矢量；根据所有目标的接收信号矢量，计算不同目标接收信号矢量之间的相关系数；根据目标接收信号矢量之间的相关系数，计算目标相关次数矢量；根据目标相关次数矢量，与预设的目标判别阈值进行比较，得到目标的真假鉴别结果，本发明可以在多站雷达系统中存在对齐误差的场景下完成有源假目标鉴别。场景下完成有源假目标鉴别。场景下完成有源假目标鉴别。


技术研发人员：赵珊珊 谢飚 刘子威 易民举
受保护的技术使用者：南京邮电大学
技术研发日：2023.05.15
技术公布日：2023/7/12