一种山基的掩星预报方法

1.一种山基的掩星预报方法，其特征在于，包括以下步骤：(1)获取所有GNSS导航星位置通过掩星探测仪的定位设备获取各GNSS导航卫星的广播电文，根据各星座公布ICD接口文件获得所有导航星在时间T的位置PGNSS；(2)获取山基观测站的位置以及经度、纬度和高度利用定位设备获取观测站的位置P，并将位置转化为大地经纬高；(3)获取所有GNSS导航星相对于观测站的俯仰角Elv、方位角Azm及切点高度h方位角和俯仰角Elv，计算公式如下：式中x RELnew ,y RELnew ,z RELnew 是GNSS导航星相对于山基观测站转化为东北天坐标下的相对矢量的各轴分量，由以下转化公式计算而来：RRELnew＝RREL*C式中C为转换矩阵，可以通过步骤(2)获得的观测站经纬高计算取得，用来将WGS-84坐标系下山基观测站的位置计算转换到东北天坐标系；RREL为GNSS导航星相对于山基观测站的相对矢量，由如下公式计算得到：RREL＝PGNSS-P式中P为得到的观测站的位置；PGNSS为步骤(1)中得到的各GNSS导航星在时间T的位置；(4)判断是否符合门限利用步骤(3)中得到的俯仰角Elv、方位角Azm和山基观测站中天线的安装方位，来判断是否满足掩星事件的门限；若满足，则判定该导航星为掩星；(5)如果步骤(4)中的天线方位满足约束，则时间外推20min时间，分别计算步骤(1)～(3)是否满足步骤(4)中的掩星事件门限，如果两次均满足，则判定该导航星为掩星。 2.按照权利要求1所述的山基的掩星预报方法，其特征在于，所述步骤(4)中的掩星事件的门限为：Elv min ≤Elv≤Elv max ，Azm min ≤Azm≤Azm max ；式中：Elv max 为掩星的俯仰角门限上限，Elv min 为掩星的俯仰角门限下限，Azm max 为掩星的方位角门限上限，Azm min 为掩星的方位角门限下限：所述Elv max 、Elv min 、Azm max 及Azm min 由如下公式计算：Elv max ＝5Elv min ＝Elv tanget -2Azm max x ＝w/2+azAzm min x ＝-w/2+az式中：Elv max 根据经验值选取，w为天线的幅角，az为天线安装方位角，Elv tanget 为掩星事件切地球的最低俯仰角；所述最低俯仰角Elv tanget 由如下公式计算：式中：切点最低高度a为地球半径，H为山基观测站的高度。