1.一种地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，包括：S1.确定GNSS干扰源数目以及各GNSS干扰源的位置坐标以及配置参数；S2.根据各GNSS干扰源的配置参数，确定各GNSS干扰源的待分析区域；对于任一GNSS干扰源，通过对比GNSS干扰源在自由空间中传播的最大作用距离和GNSS干扰源在仅考虑地球曲率时的最大作用距离，取其中较小的值作为GNSS干扰源的最大作用距离，GNSS干扰源的待分析区域即以GNSS干扰源为圆心、GNSS干扰源的最大作用距离为半径的区域；将各GNSS干扰源对应的待分析区域存储在数组中；S3.在各GNSS干扰源的待分析区域内，采用基于参考面可视域分析方法对各GNSS干扰源的待分析区域内的每个目标点进行受干扰情况判断，GNSS干扰源以及与GNSS干扰源邻近网格的目标点表示受干扰；得到各GNSS干扰源作用区域；所述邻近网格包括以GNSS干扰源为原心，相邻一个单元的网格；S4.按照经纬度拼接所有GNSS干扰源作用区域得到受地形遮挡影响的多GNSS干扰源作用区域。
2.根据权利要求1所述的地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，S1中所述GNSS干扰源的配置参数包括GNSS干扰源的自身高度、干扰源频率、干扰信号波长、发射天线增益、接收天线增益、接收机抗干扰能力干信比、接收信号功率和干扰源发射功率。
3.根据权利要求2所述的地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，对于任一GNSS干扰源，通过下式计算GNSS干扰源在自由空间中传播的最大作用距离：其中：d si 表示第i个GNSS干扰源在自由空间中传播的最大作用距离；λ i 表示第i个GNSS干扰源的干扰信号波长；G Ti 表示第i个GNSS干扰源的发射天线增益；G Ri 表示接收天线增益；M Ji 表示接收机抗干扰能力干信比；P Si 表示接收信号功率P Si ；P JTi 表示第i个GNSS干扰源的发射功率。
4.根据权利要求2所述的地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，对于任一GNSS干扰源，通过下式计算GNSS干扰源在仅考虑地球曲率时的最大作用距离；其中设定地球为均匀球体其半径为R；d earthi 表示第i个GNSS干扰源在仅考虑地球曲率时的最大作用距离；H IFi ＝h i +H Si ，H IFi 表示第i个GNSS干扰源的高度，H Si 为第i个GNSS干扰源所在位置的实际地形高程值，h i 为第i个GNSS干扰源的自身高度。
5.根据权利要求1至4中任一权利要求所述的地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，S3中，对于任一GNSS干扰源，通过以下步骤获得其GNSS干扰源作用区域：(1)采用规则DEM格网将GNSS干扰源的待分析区域进行二维行列划分，GNSS干扰源的待分析区域划分为多个单元网格，GNSS干扰源对应DEM格网的中心网格点，DEM格网的每个网格点对应待分析区域内的一个目标点；(2)以GNSS干扰源为原点，将待分析区域用正北、西北、正西、西南、正南、东南、正东及东北这八个方向线划分为八个区域块；(3)对于位于八个方向线上的目标点，根据目标点与GNSS干扰源之间的距离，按照由近及远的判断顺序对各目标点进行受干扰情况判断；(4)对位于八个区域块内的目标点，根据目标点与GNSS干扰源之间的距离，按照由近及远的判断顺序对各目标点进行受干扰情况判断。
6.根据权利要求5所述的地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，步骤(3)中，GNSS干扰源以及与距离GNSS干扰源仅一个单元网格的八个方向线上的目标点均赋值为1，代表受干扰；对距离GNSS干扰源2个以及2个以上单元网格的八个方向线上的目标点进行受干扰情况判断，方法如下：在待分析区域objDistrict maxi 对应的二维DEM格网中，对于方向线上当前待进行受干扰情况判断的目标点O m ，该目标点的实际高程值为 与目标点O m 相邻的八个网格点中离GNSS干扰源最近的一个网格点即为目标点O m 的辅助网格点rp，辅助网格点rp对应的实际高程值是考虑地形遮挡后此点可受干扰的最小高程，即为该网格点的可视高程值visual_H rp ；在待分析区域objDistrict maxi 的三维空间中，连接GNSS干扰源与目标点O m 的辅助网格点rp并延长，计算目标点O m 在该直线上投影的垂直高程值 若 则目标点O m 受到干扰，赋值为1，且目标点O m 的可视高程值 若 则目标点O m 不会受到干扰，赋值为0，且该点的可视高程值
7.根据权利要求6所述的地形遮挡影响的GNSS干扰源作用区域计算方法，其特征在于，步骤(4)中，距离GNSS干扰源仅一个单元网格的八个区域块内的目标点均赋值为1，代表受干扰；对距离GNSS干扰源2个以及2个以上单元网格的八个区域块内的目标点进行受干扰情况判断，方法如下：在待分析区域objDistrict maxi 对应的二维DEM格网中，对于区域块内当前待进行受干扰情况判断的目标点O n ，该目标点的实际高程值为 与目标点O n 相邻的八个网格点中离GNSS干扰源最近的两个网格点即为目标点O n 的两个辅助网格点rp1与rp2，辅助网格点rp1与rp2对应的实际高程值即为辅助网格点rp1与rp2的可视高程值visual_H rp1 和visual_H rp2 ；在待分析区域objDistrict maxi 的三维空间中，GNSS干扰源与目标点O n 的两个辅助网格点在空间中确定一个平面，计算目标点O n 在该平面上投影的垂直高程值 若 则目标点O n 受到干扰，且目标点的可视高程值 若 则目标点O n 不会受到干扰，且目标点O n 的可视高程值
8.一种存储介质，其上存储有计算机程序，其特征在于：所述计算机程序在被处理器运行时执行权利要求1所述的地形遮挡影响的GNSS干扰源作用区域计算方法的步骤。
9.一种计算机系统，包括机体与设在机体内的机载电路板，所述机载电路板上设有处理器与存储器，所述存储器存储有计算机程序，其特征在于：所述处理器执行所述计算机程序时实现权利要求1所述的地形遮挡影响的GNSS干扰源作用区域计算方法的步骤。
