1.一种基于分区最佳空间构型遍历与迭代优化的选星方法,其特征在于,包括:S100,获取当前历元下所有卫星组成第一卫星集合;S200,将俯仰角最大的卫星作为选定卫星,并从所述第一卫星集合中去除俯仰角低于截止俯仰角的卫星以及所述选定卫星得到第二卫星集合;S300,以所述选定卫星为基准,在所述第二卫星集合中寻找使得空间构型体积最大对应的三个目标卫星;将所述选定卫星与所述三个目标卫星组成基本卫星子集;S400,从所述第二卫星集合中去除所述基本卫星子集中的卫星得到第三卫星集合;以无放回方式在所述第三卫星集合中每次选择一个待定卫星,并将每个待定卫星与所述基本卫星子集构成一个卫星组合得到多个卫星组合;S500,选定对卫星组合影响程度系数最大的待定卫星所在的卫星组合,并确定该卫星组合的GDOP值;S600,如果S500中的所述GDOP值满足GDOP门限值要求或所述卫星组合满足卫星数目要求,则确定S500中选定的卫星组合为最优卫星组合;S200包括:S210,获取所述第一卫星集合中的所有卫星在地心地固直角坐标系中的三维坐标,并将所述三维坐标转换至站心坐标系下得到站心坐标系下的卫星坐标;S220,利用S210中的所述卫星坐标确定所有卫星的方位角和俯仰角;S230,选择俯仰角最大的卫星作为选定卫星,并从所述第一卫星集合中去除俯仰角低于截止俯仰角的卫星以及所述选定卫星得到第二卫星集合;S210包括:S211,获取所述第一卫星集合中所有卫星在地心地固直角坐标系下的三维坐标、观测站在所述地心地固直角坐标系中的三维坐标以及所述观测站的大地经纬高坐标;S212,利用S211中的观测站在所述地心地固直角坐标系中的三维坐标以及所述观测站的大地经纬高坐标,将S211中的所有卫星在地心地固直角坐标系下的三维坐标变换到站心坐标系下,得到所有卫星在所述站心坐标系下的卫星坐标;S300包括:S310,将所述第二卫星集合中所有卫星投影到二维圆面得到卫星投影点;S320,通过分层聚类对所有卫星投影点进行聚类处理,得到在二维平面内包含至少一个卫星投影点的三个卫星投影点区域;S330,以所述选定卫星为基准,在所述三个卫星投影点区域中寻找使得空间构型体积最大对应的三个目标卫星;S340,将所述三个目标卫星与所述选定卫星组成基本卫星子集;S320包括:S321,将所有卫星投影点分别作为初始卫星投影点区域,并根据每对卫星投影点之间的距离计算初始的距离矩阵;S322,使用平均链接合并策略以及所述初始的距离矩阵,寻找距离较近的卫星投影点区域进行迭代合并,直到得到在二维平面内的三个卫星投影点区域,每个卫星投影点区域至少包括一个卫星投影点;S322包括:S3221,通过所述初始的距离矩阵寻找最接近的两个卫星投影点区域,并将其合并成一个卫星投影点区域;S3222,更新卫星投影点区域;S3223,重新计算各个卫星投影点区域之间的距离,并利用所述距离更新距离矩阵;S3224,重复S3221-S3224直到得到在二维平面内的三个卫星投影点区域。
2.根据权利要求1所述的基于分区最佳空间构型遍历与迭代优化的选星方法,其特征在于,S500包括:S510,利用所述基本卫星子集的观测矩阵,计算每个卫星组合的观测矩阵;S520,针对每个卫星组合,利用该卫星组合的观测矩阵计算对应的待定卫星对该卫星组合的影响程度系数;S530,选定影响程度系数最大的待定卫星所在的卫星组合,并确定该卫星组合对应的GDOP值。
3.根据权利要求1所述的基于分区最佳空间构型遍历与迭代优化的选星方法,其特征在于,在S600之后,所述基于分区最佳空间构型遍历与迭代优化的选星方法还包括:如果S500中的所述GDOP值不满足GDOP门限值要求且卫星组合满足卫星数目要求,则用所述卫星组合替换所述基本卫星子集,并重复S400至S500直至得到最优卫星组合。