[发明专利]一种资源充足情形下最小观测成本的卫星任务规划方法

权利要求书

1.一种资源充足情形下最小观测成本的卫星任务规划方法，其特征是应用于由一个待观测矩形区域R和n个成像卫星的覆盖机会集合S＝{s₁,s₂,...,s_i,...s_n}构成的任务规划场景中；其中，s_i表示第i个覆盖机会，1≤i≤n；所述卫星任务规划方法是按如下步骤进行：

步骤1.参数定义和初始化：

以所述待观测矩形区域R的任意一个顶点为原点o，以与所述原点相邻的两条边分别为x轴和y轴，从而建立坐标系o-xy；

所述任务规划场景中，每个覆盖机会对应一个卫星到地面的垂直投影点直线轨迹，简称为星下点轨迹，记第i个覆盖机会s_i对应的星下点轨迹为o_i；

每个覆盖机会对应一个卫星到地面高度，记第i个覆盖机会s_i对应的卫星到地面高度为h_i；

每个覆盖机会对应一个最大观测长度，记第i个覆盖机会s_i对应的最大观测长度为d_i；

每个覆盖机会对应一个最大偏转角度，记第i个覆盖机会s_i对应的最大偏转角度为p_i；

每个覆盖机会对应一个相机视场角度，记第i个覆盖机会s_i对应的相机视场角度为w_i；

每个覆盖机会对应一个单位距离观测成本，记第i个覆盖机会s_i对应的单位距离观测成本为m_i；

由所述第i个覆盖机会s_i的星下点轨迹o_i、卫星到地面高度h_i、最大观测长度d_i、最大偏转角度p_i、相机视场角度w_i和单位距离观测成本m_i共同构成所述第i个覆盖机会s_i的属性；

步骤2.定义当前迭代次数为k，并初始化k＝1，定义最大迭代次数为K；

步骤3.将所述待观测矩形区域R进行第k次划分，得到若干个大小相等的正方形网格所组成的第k个网格R_k，将第k个网格R_k中的每个小正方形称为单元格；并对每个单元格进行编号为1,2,…,j,…,Q_k，且记录每个单元格的四个顶点的坐标位置；

步骤4.对每个覆盖机会得到左单元格集合、上单元格集合和下单元格集合：

步骤4.1.初始化i＝1；

步骤4.2.从第k个网格R_k中遍历得到第i个覆盖机会s_i的所有左单元格，并构成左单元格集合其中，表示第i个覆盖机会s_i的第m个左单元格；M_i表示第i个覆盖机会s_i的左单元格总数,m＝1,2,…,M_i；

步骤4.3.从第k个网格R_k中遍历得到第m个左单元格的上单元格，并构成上单元格集合其中，表示第i个覆盖机会s_i的第m个左单元格的第e个上单元格,e＝1,2,…,E_m,从而得到第i个覆盖机会s_i的所有左单元格的上单元格集合

步骤4.4.从第k个网格R_k中遍历得到第m个左单元格的第e个上单元格的第y个下单元格，并构成下单元格集合其中，表示第i个覆盖机会s_i的第m个左单元格的第e个上单元格的第y个下单元格，y＝1,2,…,Y_e,从而得到第i个覆盖机会s_i的所有左单元格的上单元格的下单元格集合其中

步骤4.5.将i+1赋值给i，判断in是否成立，若成立则执行步骤5，否则转步骤4.2；

步骤5.对每个覆盖机会生成若干个覆盖模式，即观测条带：

步骤5.1.令i＝1，初始化第k个覆盖模式总集合C_k为空集；

步骤5.2.初始化第k个覆盖模式总集合C_k中第i个覆盖机会s_i的覆盖模式集合为空集，令m＝1；

步骤5.3.令e＝1；

步骤5.4.令y＝1；

步骤5.5.根据第i个覆盖机会s_i及其左单元格集合U_i中第m个左单元格以及其第e个上单元格以及其第y个下单元格利用覆盖模式生成法生成覆盖模式并存入覆盖模式集合中；

步骤5.6.将y+1赋值给y，并判断yY_e是否成立，若成立，则执行步骤5.7，否则转步骤5.5；

步骤5.7.将e+1赋值给e，并判断eE_m是否成立，若成立，则执行步骤5.8，否则转步骤5.4；

步骤5.8.将m+1赋值给m，并判断mM_i是否成立，若成立，则将覆盖模式集合存入第k个覆盖模式总集合C_k后，执行步骤5.9，否则转步骤5.3；

步骤5.9.将i+1赋值给i，判断in是否成立，若成立，则获得第k个覆盖模式总集合并执行步骤6，否则转步骤5.2；

步骤6.使用基于动态贪婪的启发式算法从第k个覆盖模式总集合C_k中选出部分覆盖模式构成第k个可行解P_k；

步骤7.将k+1赋值给k，判断kK是否成立，若成立，输出第K个可行解P_K，否则执行步骤8；

步骤8.将第k-1个网格R_k-1作为父网格，将父网格中的每个单元格等分为小正方形，得到若干个大小相等的小正方形所组成的第k个网格R_k，称R_k为子网格，将所述第k个网格R_k中的小正方形作为新的单元格，记录每个新单元格的四个顶点坐标，从而完成一次嵌套父子网格的构造；

步骤9.将第k-1个可行解P_k-1中各个覆盖模式映射到第k个网格R_k下，即将各个覆盖模式的左单元格和上单元格更新为第k个网格R_k中的左单元格和上单元格，并对映射后的第k个网格R_k上的每个覆盖模式构造若干个临近覆盖模式后，一起存入第k个覆盖模式总集合C_k，再返回步骤6。