一种评估星群变轨机动的多重可达覆盖性方法

本发明涉及一种评估星群变轨机动的多重可达覆盖性方法，适用于卫星星座多星可达区域的快速评估，属于航空航天领域。

背景技术：

1、卫星轨道机动可达域可表征其在未来一段时间内可能到达的空间范围，对维护航天器在轨安全、改善空间态势感知能力具有重要意义。针对单个卫星的可达域计算已经较为全面，包括脉冲可达域和连续推力可达域。这类可达域采用三维网格表示，对于点目标的可达判定较为简便，等效于判断点是否位于多面体内部。对于三维空间目标需要将其表面进行离散取点再进行点的可达判定。针对多星可达域计算的研究较少，现有方法多是将每颗卫星对于目标的可达特性进行统计分析。对于分别由n个面描述可达包络的m颗卫星，m个离散点表示的目标区域的多星可达性需要进行mmn次计算判定，所需计算量难以接受。同时，仅表面离散取点无法量化计算多星可达目标区域的体积，而目标整体区域的离散取点所需m的数目极大，进一步增大了计算量。因此，多星可达域的量化计算和几何可视化通常借助于现有三维网格模型的布尔运算算法。

2、先技术[1](参见三维网格模型的空间布尔运算[j]《华中科技大学学报(自然科学版)》,毕林,王李管,陈建宏,冯兴隆,2008,294(5):82-85)进行三维网格模型的空间布尔运算。该方法通过相交测试求出两两相交三角形面片之间的交线，再由相交三角形与它的交线得到多边形，并对该多边形进行重新三角网格化，最后根据重组后模型判断其他网格的取舍。然而，若简单地进行相交测试，对于由n个网格三角面片描述的两个三维区域进行一次布尔运算需要进行n(n-1)/2次相交判断，这种做法的耗时将不可接受。

3、因此通常利用特定技术进行快速相交部位的检测。先技术[2](参见基于八叉树的三维网格模型体素化方法[j]《工程图学学报》,吴晓军,刘伟军,王天然,2005(4):1-7)递归地将空间立方体进行划分，划分为8个小空间立方体，只要对位于相同子空间内的三角面片进行相交检测，即可快速定位相交的三角形，提高布尔运算算法的整体运行效率。

4、另一方面，先技术[3](参见obbtree:a hierarchical structure for rapidinterference detection[c]《proceedings of the 23rd annual conference oncomputer graphics and interactive techniques》,gottschalk s,lin m c,manocha d,1996)提出obb包围盒算法，提前预处理筛选排除掉部分不可能相交的三角面片，减少需要求交的面片对的数量，提高布尔运算操作的运行效率。

5、然而，这类算法是基于相交环提取最终模型网格，因此需要考虑各类特殊、复杂的情形，算法的鲁棒性有待提闻。同时，这类算法始终只能对于两两模型依次操作，对于m颗卫星进行n重可达域计算需要进行m！/(n-1)！(m-n)！次布尔运算，所需计算量依然难以接受。因此，本发明提出一种基于线覆盖的多星可达域计算方法。

技术实现思路

1、本发明主要目的是提供一种评估星群变轨机动的多重可达覆盖性方法，在地心固连球坐标系下将目标区域在距离和天顶角两个维度进行等体积离散，并以离散单元中心线的可达方位角区间作为评价指标；通过计算中心线与可达域包络的交点并进行详细判断，得到单星可达方位角区间；采用一种节点排序和一次性顺序读取方法，得到对应目标中心线的多星可达方位角区间，同时结合目标方位角区域对于多星可达方位角区间进行截取；整合所有目标中心线的多星可达区间得到星座对于目标区域的多星可达区域，即实现星群变轨机动的多重可达覆盖性评估。

2、本发明的目的是通过下述技术方案实现的。

3、本发明公开的一种评估星群变轨机动的多重可达覆盖性方法，包括如下步骤：

4、步骤一：在地心固连坐标系si或地心惯性坐标系sj下，将目标区域s用部分三维球层表示，将部分三维球层在距离r和天顶角θ维度进行等体积离散，并将方位角维度进行解析，得到i×j个等体积环rij，并以等体积环中心线cij的可达情况作为等体积环rij的可达情况。

5、首先定义采用的坐标系，考虑地球惯性坐标系、地球固连坐标系和卫星轨道坐标系：

6、地心惯性坐标系sj，其坐标原点位于地球引力中心e，基准平面定义为地球平均赤道平面。x轴指向平春分点，z轴垂直于赤道面指向北极，y轴由右手法则确定：y＝z×x。

7、地心固连坐标系si，原点位于地球引力中心e，基准平面定义为地球平均赤道平面，xi轴沿格林尼治子午面和地球赤道平面的交线，zi轴垂直于赤道面指向北极，yi轴由右手法则确定：yi＝zi×xi。

8、卫星轨道坐标系s0，以卫星质心o为坐标原点，以航天器当前位置矢量方向rsat为x0轴，z0轴垂直轨道平面的指向正法线方向，y0由右手法则确定，即y0＝z0×x0。

9、对于单星可达域计算，可达域包络表示在轨道坐标系s0下。对于多星可达域计算，可达域包络需要统一表示在惯性坐标系sj下。对于由轨道六根数表示为[a,e,i,ω,ω,f]的具体卫星，通过如公式(1)所示转换矩阵将可达域包络需要统一表示在惯性坐标系sj下：

10、

11、其中u＝ω+f，m为绕对应轴的旋转矩阵。

12、地心惯性系与地心固连系的转换中仅考虑地球匀速自转因素，忽略岁差、章动和极移的影响。坐标系仅围绕z轴旋转，对应转换矩阵为：

13、

14、其中gast为格林尼治平恒星时，为计算卫星位置时刻平春分点到格林尼治子午线的地球旋转角。

15、根据目标特性不同在地心固连坐标系si或地心惯性坐标系sj下，目标区域s由经纬度和距离给定，或者由三维网格包络给定。根据s和任务需要设定距离区间[rl,ru]和天顶角区间[θl,θu]定义将三维球层表示为：

16、

17、则部分三维球层的体积为：

18、

19、对于三维球层在距离r和天顶角θ维度进行离散，并在方位角维度进行解析，得到的i×j个等体积离散圆环rij表示为：

20、

21、其中ri和θj为进行离散的节点满足：

22、r0＝rl,ri＝ru,θ0＝θl,θj＝θu

23、给定三维信息的离散环体积vij为：

24、

25、随着天顶角和距离的增大，等间距(同δθ＝θj-θj-1和δr＝ri-ri-1)网格法所表示的离散单元体积差距越来越大，不利于星座的可达特性分析。采用等体积离散法确定节点ri和θj。部分三维球层的体积由式(3)给出，则每个离散环的体积vij为

26、

27、结合式(4)和式(5)，节点ri和θj表示为：

28、

29、

30、将离散单元rij的中心线定义为cij，表示为：

31、

32、其中cij是由节点ri和θj决定的圆，中心线的可达情况由星座中卫星可达域对应的方位角区间表示。

33、当离散单元的参数i和j足够大时，所划分的离散环体积vij较小，对于单个离散单元的可达情况rrd(ij)使用其中心线可达情况φrd(ij)的近似性能上升，表示为：

34、

35、以距离和天顶角二维离散计算量i×j实现三维空间区域目标的等体积离散，以解析的方位角中心环可达区间表示空间区域目标，适用于以经纬度或者网格包络给定的任意目标区域。

36、步骤二：针对于步骤一等体积离散得到的等体积环中心线cij和转换到惯性坐标系sj下的单星可达包络，通过将单星可达三维包络拆分为多个三角面片，在球坐标系下计算中心线cij和三角面片的面内交点，得到对应的单星可达方位角矩阵和目标方位角矩阵φtar(ij)。

37、单星可达包络在轨道坐标系s0下表示为r(x,y,z)0，其中三轴位置均由m×n的网格矩阵给出。将其转换到惯性坐标系sj下表示为r(x,y,z)j，或在固连坐标系si下表示为r(x,y,z)i。

38、计算卫星可达域对应的方位角区间即求中心圆cij与构成包络的所有网格面交点。为简化计算将包络表示在球坐标下，同样采用m×n的网格矩阵表示，问题等效为求r＝(ri-1+ri)/2且θ＝(θj-1+θj)/2的直线与构成包络的所有网格面交点。

39、将m×n网格矩阵中相邻的2×2矩阵等效为一个非平面的网格面单元，共k＝(n-1)(m-1)个。分别计算距离和天顶角的极值，判断cij是否位于网格面的r和θ范围内，筛选出的网格面表示为集合k。

40、对于初步筛选出的网格面集合k，判断不重复点数，点数为3时构成1个三角面片，点数为4时构成两个三角面片。对于所有三角面片求与cij的交点，并判断交点是否位于三角形内部。

41、计算集合k中所有满足要求的交点p的方位角。由于包络为闭合曲面，交点个数必然为偶数2km。将第m颗卫星的可达域对于cij的交点按方位角从小到大顺序排列得到序列为：

42、

43、其中+1表示可达区间起点，-1表示可达区间终点，共km个区间

44、由于将所有网格点转换到球坐标系下时，方位角取值范围为[0,2π]。当时，存在两种情况：1)z轴与包络相交；2)包络与x轴正向和z轴正向构成的平面相交，z轴与包络不相交。将方位角取值范围为变更为[-π,π]，若仍确定为情况1；反之，确定为情况2。

45、对于情况1，当cij与球坐标系下的包络存在交点时，进一步判断交点对应的是可达区域还是不可达区域。在笛卡尔固连坐标系下，取两个方位角节点的中点沿z轴正向作射线与可达域包络求交点，方式与球坐标下相同。若交点为奇数，则对应方位角区间为可达区域；反之，为不可达区域，变化对应两个节点的正负性，即

46、若cij与球坐标系下的包络不存在交点时，需要进一步判断cij是否位于包络内部。同样笛卡尔惯性坐标系下，任取cij上一点沿z轴正向作射线与可达域包络求交点。若交点为奇数，则cij位于包络内部，可达方位角区间为[0,2π]。

47、对于情况1和情况2，求得交点的方位角可能位于[-π,0]区间，转换到[0,2π]内以便下步骤进行计算。将方位角位于[-π,0]区间的节点等效转换到[0,2π]区间内，再次按方位角从小到大顺序排列。若第一列第二行为-1，最后一列第二行必为+1，添加一个起始节点和一个终止节点

48、通过将单星可达三维包络拆分为多个三角面片，在球坐标系下计算中心线cij和三角面片的面内交点，根据公式(5)得到任意条件下单星可达域与目标中心线cij的方位角序列目标方位角序列φtar(ij)采用同样的方式得到。

49、步骤三：对于步骤二得到的采用起始和终止节点表示的第m颗卫星的可达域对于cij的可达方位角序列通过节点排序和一次性顺序读取得到星座对于中心线cij的多星可达矩阵根据φtar(ij)中目标方位角区间对进行截取，得到给定目标区域的多星可达矩阵

50、将步骤二得到的所有m颗卫星的可达方位角进行合并得到：

51、

52、此时，矩阵中的节点数目为：

53、

54、接着，将矩阵按照第一行升序排列得到：

55、

56、将可达星数定义为n并设定其初始值n0＝0。然后，顺序读取矩阵的第二行元素。当第s个元素为+1(s∈{1,2,...,s})，对应的第一行角度表示一个新可达区间的起始点，因此可达星数增加ns＝ns-1+1。相反，当第s个元素为-1(s∈{1,2,...,s})，对应的第一行角度表示一个已有可达区间的终止点，因此可达星数减小ns＝ns-1-1。当读取完所有元素后，将对应的ns表示在矩阵第三行得到一个全新的矩阵为：

57、

58、其中ns为区间的可达星数。

59、随后，需要根据目标方位角区间φtar(ij)对于多星可达矩阵进行变换。将φtar(ij)中的2p个方位角节点表示为并与的节点一同按照方位角升序排列，得到新矩阵为：

60、

61、其中所在列的第三行为前一列第三行的可达星数n，向中第一列添加全0列，使得时当多个点具有与相同的天顶角时，需要确保位于排序最后，位于排序最前。

62、定位和所在列并对进行截取，得到p个对应目标多星可达矩阵表示为：

63、

64、通过节点排序和一次性顺序读取得到星座对于中心线cij给定目标区域的多星可达矩阵

65、步骤四：根据步骤三得到的中心线cij给定目标区域的多星可达矩阵计算单个离散单元rij的n星可达方位角γrdnx(ijp)(n)，整合所有目标中心线的多星可达区间得到星座对于目标区域的多星可达区域，即实现星群变轨机动的多重可达覆盖性评估。

66、最大可达星数定义为nmax＝max(nps)，初始多星可达角度定义为γrdnx(ijp)(n)＝0，其中n∈{1,...,nmax}。然后，遍历第三行中ns＝n对应的方位角节点，则n星可达方位角度为：

67、

68、在得到每个离散单元中心线cij对应目标区间的n星可达角度之后，可达角度γcovnx(ijp)(n)与目标角度的占比被近似为整个离散单元中目标区域的可达率。则中心线cij所有目标的多星可达体积vrdnx(ij)(n)满足：

69、

70、当单个离散单元的可达区间给定后，其可达体积可直接得到。所有离散单元的可达体积之和即是整个空间目标区域的可达体积。对于i×j个等体积离散圆环，其可达体积表示为：

71、

72、其中vrdnx(n)是整个空间目标区域的n星可达体积。

73、还包括步骤五：根据步骤四评估得到的星群变轨机动的多重可达覆盖性，实现多星可达区域的可视化表征，提高对多星协同中多星可达区域评估的精度和效率，进而提高多星协同任务执行的精度和效率。所述多星协同任务包括观测、制导、拦截。

74、有益效果：

75、1、本发明公开的一种评估星群变轨机动的多重可达覆盖性方法，在地心固连球坐标系下，以i×j个离散圆环中心线的方位角区间表示可达区域和目标区域，以离散单元中心线的可达方位角区间作为评价指标，通过将单星可达三维包络拆分为多个三角面片并解析计算面内交点，可达方位角精度和计算量与三维网格包络相近。同时，以起始和终止节点表示的n个方位角区间能够进行基于节点排序和顺序读取的多星可达区间快速计算，时间复杂度为o(n2+n)。而三维布尔运算计算n颗卫星对应的m重覆盖区域的时间复杂度为o(n！/(m-1)！(n-m)！)。

76、2、本发明公开的一种评估星群变轨机动的多重可达覆盖性方法，通过分常规有效区域相交和z轴与包络相交包络、x轴正向和z轴正向构成的平面相交且z轴与包络不相交两种特殊情况，得到空间任意位置和任意形状的三维包络与离散圆环中心线的相交方位角区间，适用于以经纬度或者网格包络给定的任意卫星可达区域和目标区域。

77、3、本发明公开的一种评估星群变轨机动的多重可达覆盖性方法，在实现上述有益效果1、2基础上，实现多星可达域的量化评估和几何可视化，根据多星可达区域的可视化表征结果，提高对多星协同中多星可达区域评估的精度和效率，进而提高多星协同任务执行的精度和效率。所述多星协同任务包括观测、制导、拦截等。