[发明专利]一种高精度卫星导航装置干扰测向方法

权利要求书

1.一种高精度卫星导航装置干扰测向方法，其特征在于，该方法适用于对阵列天线的测向干扰进行计算，所述阵列天线包括多个阵元，所述方法包括：

步骤1，根据每一个所述阵元接收到的阵列信号，采用干扰网格遍历测向的方式，依次计算每一个测向干扰对应的测向结果，其中，所述测向结果包括网格测向结果、第一网格测向误差以及网格偏差估计量以及第二网格测向误差；

步骤2，当判定所述测向结果与测向平均值之间的偏差小于或等于预设阈值时，将所述测向结果记作测向稳定值，并计算第一预设数量个稳定值的平均值，将所述平均值记作干扰测向值，

所述步骤1中所述网格测向结果的计算方法，包括：

步骤11，对所述阵列天线的接收区域进行第一网格划分，并依次计算每一个第一网格的第一导向矢量，其中，所述第一导向矢量按照卫星导航信号频点序号、俯仰网格序号、水平网格序号进行标记；

步骤12，采用遍历的方式，依次增加所述俯仰网格序号、所述水平网格序号，并根据所述第一导向矢量，计算每一个所述第一网格的第一网格投影，所述第一网格投影的计算公式为：

式中，为所述第一网格投影，β为干扰阵列信号初始化状态为零中频数字阵列信号，为所述第一导向矢量，为所述第一导向矢量的转置共轭，(*)^-1为矩阵求逆运算，θ为所述俯仰网格序号，9≥θ≥0，为所述水平网格序号，f_i为第i个卫星导航信号频点对应的频率，i＝1,2,…,I，I为卫星导航信号频点最大序号；

步骤13，将所述第一网格投影的最大值，记作所述网格测向结果。

2.如权利要求1所述的高精度卫星导航装置干扰测向方法，其特征在于，所述步骤2中所述测向平均值的计算方法，其特征在于，所述计算方法包括：

统计当前所述测向结果对应时刻后的第二预设数量个测向和值；

计算统计出的测向结果的平均值，记作所述测向平均值，其中，所述测向和值为所述网格测向结果与对应的网格偏差估计量的和值。

3.如权利要求2所述的高精度卫星导航装置干扰测向方法，其特征在于，所述预设阈值的取值为3°。

4.如权利要求1所述的高精度卫星导航装置干扰测向方法，其特征在于，所述步骤1中所述第一网格测向误差的计算公式为：

式中，δ(m)为第m个测向干扰的所述第一网格测向误差，m＝1,2,…,M，M为测向干扰的个数，δ(0)为第一初始误差，DOA_g(m)为第m个测向干扰的所述网格测向结果。

5.如权利要求4所述的高精度卫星导航装置干扰测向方法，其特征在于，所述步骤1中所述网格偏差估计量的计算方法，包括：

步骤14，对所述网格测向结果DOA_g(m)对应的第一网格进行第二网格划分，并根据第二网格的第二导向矢量和所述第一网格的第一导向矢量，计算每一个所述第二网格的网格偏差；

步骤15，根据所述网格偏差，计算每一个所述第二网格的小波投影，并根据所述小波投影和所述第一网格测向误差，计算所述网格偏差估计量，所述网格偏差估计量的计算公式为：

式中，DOA_c(m)为第m个测向干扰的网格偏差估计量，为第个第二网格的小波投影，θ′为第二网格的俯仰网格序号，为第二网格的水平网格序号，为网格偏差，t为时间序号，t＝1,2,…,T，T的取值由第二网格的数量确定，sin^*(/10)为小波基函数sin(/10)的转置共轭，为小波投影的转置共轭。

6.如权利要求5所述的高精度卫星导航装置干扰测向方法，其特征在于，所述步骤1中所述第二网格测向误差的计算公式为：

式中，E(m)为第m个测向干扰的所述第二网格测向误差，m＝1,2,…,M，M为测向干扰的个数，E(0)为第二初始误差。

7.如权利要求1至6中任一项所述的高精度卫星导航装置干扰测向方法，其特征在于，所述阵列信号为多路零中频数字阵列信号。