一种适用于中高轨卫星的电磁干扰定位方法与流程

技术特征：
1.一种适用于中高轨卫星的电磁干扰定位方法，其特征在于，该电磁干扰定位方法包括：步骤1，根据卫星系统中的电磁监测结果，确定受干扰的卫星；步骤2，根据所述受干扰的卫星的覆盖范围，计算全球地图覆盖矩阵，并对所述全球地图覆盖矩阵求和，选取求和后全球地图覆盖矩阵中最大元素对应区域作为初始定位区域；步骤3，采用遍历的方式，依次判断选取出的所述初始定位区域内是否满足三颗或三颗以上卫星共视，若是，执行步骤4，否则，重新获取所述卫星系统中的电磁监测结果，执行步骤1；步骤4，根据共视卫星的电磁监测结果，进行多星时间差定位，计算干扰源定位坐标，其中，所述共视卫星为所述初始定位区域内共视的卫星。2.如权利要求1所述的适用于中高轨卫星的电磁干扰定位方法，其特征在于，所述步骤2中计算全球地图覆盖矩阵，具体包括：步骤21，根据经纬度和预设间隔，将全球的地图展为二维图，并进行网格化，根据网格化后的二维图，生成全球地图覆盖初始矩阵，其中，所述全球地图覆盖初始矩阵中的元素取值为0；步骤22，根据所述受干扰的卫星的覆盖范围，计算所述受干扰的卫星与所述网格化后的二维图中各网格中心位置之间的距离；步骤23，当判定所述距离小于或等于星地极限距离时，将所述距离对应的所述全球地图覆盖初始矩阵中的元素取值赋值为1，以计算出所述全球地图覆盖矩阵。3.如权利要求1所述的适用于中高轨卫星的电磁干扰定位方法，其特征在于，所述步骤4中多星时间差定位的计算公式为：式中，δt
i
为干扰源的干扰源信号到达主星与第i个副星的时间差，主星为所述共视卫星中的任一个，副星为除所述主星外的其余共视卫星，x＝[x y z]
t
为所述干扰源定位坐标，x0＝[x
0 y
0 z0]
t
为所述主星的位置，x
i
＝[x
i y
i z
i
]
t
为第i个所述副星的位置，p为所述副星的数量，c为电磁波传播速度，||
·
||为范数函数。4.如权利要求1或3所述的适用于中高轨卫星的电磁干扰定位方法，其特征在于，所述步骤4中还包括：计算干扰源信号到达主星与副星的时间差，具体包括：步骤41，对所述主星接收到的干扰信号和所述副星接收到的干扰信号进行离散傅里叶变换，并计算离散傅里叶变换后的所述主星与所述副星接收到的干扰信号的互谱参数，所述互谱参数的计算公式为：
式中，k为离散点，[
·
]为取整函数，n为总离散点数；步骤42，对所述互谱参数进行傅里叶反变换，计算r
12
(τ)，对应的计算公式为：式中，r
12
(τ)为互谱参数的傅里叶反变参数，f
k
为dft离散傅里叶变换之后的第k个离散点的数字频率所对应的射频频率，τ为参考时间差；步骤43，选取所述r
12
(τ)最大值对应的参考时间差，记作所述干扰源信号到达主星与副星的时间差。5.如权利要求4所述的适用于中高轨卫星的电磁干扰定位方法，其特征在于，所述计算干扰源信号到达主星与副星的时间差，具体还包括：当判定信噪比小于预设阈值时，根据预设功率阈值，选取互谱参数大于所述预设功率阈值的离散点组成较强频点集合k
i
，对应的计算公式为：式中，γ0为所述预设功率阈值；根据所述较强频点集合k
i
，计算所述互谱参数的傅里叶反变参数r
12
(τ)，对应的计算公式为：6.如权利要求1所述的适用于中高轨卫星的电磁干扰定位方法，其特征在于，所述电磁干扰定位方法中还包括：步骤5：根据预设距离误差，将计算出的干扰源定位坐标与干扰知识库中的已知干扰源位置坐标进行比对，当判定所述计算出的干扰源定位坐标与所述干扰知识库中的已知干扰源位置坐标之间的距离误差小于或等于所述预设距离误差时，将所述已知干扰源位置坐标作为匹配结果，否则，将所述计算出的干扰源定位坐标作为匹配结果，并将所述计算出的干扰源定位坐标作为新的已知干扰源位置坐标，更新所述干扰知识库。

技术总结
本申请公开了一种适用于中高轨卫星的电磁干扰定位方法，该定位方法包括：步骤1，根据卫星系统中的电磁监测结果，确定受干扰的卫星；步骤2，根据受干扰的卫星的覆盖范围，计算全球地图覆盖矩阵，并对全球地图覆盖矩阵求和，选取求和后全球地图覆盖矩阵中最大元素对应区域作为初始定位区域；步骤3，采用遍历的方式，依次判断选取出的初始定位区域内是否满足三颗或三颗以上卫星共视，若是，执行步骤4，否则，重新获取卫星系统中的电磁监测结果，执行步骤1；步骤4，根据共视卫星的电磁监测结果，进行多星时间差定位，计算干扰源定位坐标。通过本申请中的技术方案，实现了对地表大功率干扰源监测与定位，进而提升卫星系统上行链路的安全性。全性。全性。


技术研发人员：陈金平 王梦丽 赵金贤 卫恒 胡进 王志印 洪源 杨洋 刘学伟 范令志 程宇峰
受保护的技术使用者：南京电子设备研究所
技术研发日：2022.04.01
技术公布日：2022/7/29