一种适用于载荷凝视成像的多目标观测轨迹智能感知方法与流程

本发明涉及航天控制领域，特别是一种适用于载荷凝视成像的多目标观测轨迹智能感知方法。

背景技术：

随着地面观测任务规模的大幅提升，传统非敏捷航天器难以满足对地观测需求。近年来敏捷航天器发展迅速，由于其具备三轴快速姿态机动能力，使得对地观测具有更高灵活性和更多选择性。为了在特定时间窗口内完成更多的点目标观测任务，如何进行合理、高效的观测轨迹感知优化已成为首要问题。一方面，直接设计更高性能的智能优化算法，来满足对地多目标点观测需求；另一方面，则对目标进行聚类，间接提高观测效率。

多目标观测轨迹感知优化是整个敏捷航天器平台观测任务的顶层决策，是完成对地多目标点观测的基础。多目标观测轨迹感知的质量直接影响观测任务的效率和完成度，进而影响观测收益。近年来，在理论研究与工程领域，多目标观测轨迹感知已引起了广泛关注。针对点目标的轨迹感知优化通常分为多目标点的聚类和观测序列决策。通过多目标点聚类模型，将能够被同时观测的目标点进行合并，减少航天器平台姿态机动的次数，为敏捷航天器观测任务的高效执行奠定基础。而观测序列决策是在多目标点聚类的基础上，通过对聚类后目标的观测时间、观测序列进行智能感知优化，最大程度完成对尽可能多目标点的观测任务。

现阶段针对敏捷航天器平台对地多目标观测轨迹感知的研究已经取得了一定的成果。但是，已有基于遗传算法、免疫算法、粒子群算法、禁忌搜索算法等智能优化算法的研究成果，主要考虑的载荷观测模式是推扫式的，尚未针对载荷凝视观测任务的特殊性开展深入的研究工作。实际上，随着面阵传感器的逐步成熟，其具有视场大、无需扫描机构等比较优势，基于载荷面阵凝视的点目标观测受到了越来越多的关注。综上所述，针对载荷凝视成像任务的特点，开展敏捷航天器平台的多目标观测轨迹智能感知方法研究，具有显著的现实意义和前瞻性。

技术实现要素：

本发明所要解决的技术问题是：克服现有技术的不足，提供一种适用于载荷凝视成像的多目标观测轨迹智能感知方法，通过多目标点的聚类模型，建立观测轨迹的感知优化模型，并利用改进的智能蚁群算法求解得到最优观测轨迹。

本发明所采用的技术方案是：一种适用于载荷凝视成像的多目标观测轨迹智能感知方法，包括步骤如下：

(1)建立所需观测多目标点的聚类模型，聚类模型以聚类点集合c的形式进行描述，在聚类点集合c中包含了多个聚类点子集合，在每个聚类点子集合中包含了个数不等的目标点；

具体步骤如下：

步骤1.1：对参数进行初始化设置，参数包括：任务点集合v、聚类点集合c、联通点集合l，进入步骤1.2；

步骤1.2：判断任务点集合v是否为空集；若任务点集合v为空集，则完成聚类模型建立，进入步骤(2)；若任务点集合v不为空集，则令聚类点集合c为空集，将任务点集合v中的任意点k加入到聚类点集合c中，并从任务点集合v中删除点k，进入步骤1.3；

步骤1.3：搜索任务点集合v中剩余的点中与聚类点集合c中所有节点都连通的点，建立连通点集合l，进入步骤1.4；

若点k与点l满足如下公式，则判定点k与点l之间是连通的：

其中：和分别表示点k与点l在地心惯性坐标系中的位置，rs、re和αp分别表示航天器平台飞行的轨道半径、地球半径和航天器载荷的视场角；

步骤1.4：判断连通点集合l是否为空集；若连通点集合l为空集，则进入步骤1.5；若连通点集合l不为空集，则将点l加入到聚类点集合c中，并将点l从任务点集合v中删除，重复步骤1.4直至任务点集合v中不再存在与点k连通的点，进入步骤1.5；

步骤1.5：得到聚类点集合c，返回步骤1.2；

(2)建立观测轨迹的感知优化模型，考虑约束条件的限制并设计优化目标，从聚类点集合c中选取满足约束条件的聚类点子集合，所建立的观测轨迹感知优化模型，具体形式如下：

其中：si和fij为智能感知优化模型的待优化变量，其计算公式分别为

n表示聚类点集合c中所包含聚类点子集合ci的个数，为正整数；i＝1,…,n；j＝1,…,n；

wi表示聚类点子集合ci的观测收益、k1表示权重系数、trij表示从聚类点子集合ci到聚类点子集合cj航天器平台姿态机动的时间、tmax表示航天器平台中执行机构能够提供的最大输出力矩；

tei表示聚类点子集合ci最早可观测的时刻、tsi表示聚类点子集合ci的实际开始观测时刻、tli表示聚类点子集合ci最晚可观测的时刻、di表示聚类点子集合ci需要持续观测的时长、tsj表示聚类点子集合cj的实际开始观测时刻；

a表示航天器平台姿态机动的角加速度、amax表示航天器平台所能达到的最大角加速度、rs表示航天器平台的位置矢量、σij表示从聚类点子集合ci到聚类点子集合cj所需机动的姿态角度；μ为地球引力常数；

聚类点子集合ci和cj之间没有交集；

(3)通过改进的智能蚁群算法进行感知优化模型的求解，得到航天器平台多目标观测的最优轨迹，具体步骤如下：

步骤3.1：设置最大迭代次数及蚂蚁数量；

步骤3.2：根据最大迭代次数，判断是否完成所有迭代，若完成则进入步骤3.10结束感知优化模型求解；否则，进入步骤3.3；

步骤3.3：将所有蚂蚁放置在感知优化开始时刻的星下点位置，并清空禁忌表、时间记录表和路径记录表，进入步骤3.4；

步骤3.4：判断是否还有未参与感知优化模型求解的蚂蚁，若有，则进入步骤3.5，否则，进入步骤3.9；

步骤3.5：令每一只蚂蚁进行路径搜索，并记录其已经访问过的聚类子集合和尚未访问的聚类子集合；对于尚未访问过的聚类子集合，计算从当前位置及时刻转移至该子集合的时间，并进入步骤3.6；

步骤3.6：根据步骤3.5所计算的时间，判断是否存在满足时间窗口约束要求的集合，若不存在则该蚂蚁完成一次完整的路径搜索，返回步骤3.4；否则进入步骤3.7；

步骤3.7：基于所设计的转移概率以聚类点子集合ci为出发点，来选择下一个需要观测的聚类点子集合cj，并进入步骤3.8；

其中：表示第k2只蚂蚁从聚类点子集合ci到聚类点子集合cj的转移概率；τij表示聚类点子集合ci到聚类点子集合cj的信息素浓度；tej表示聚类点子集合cj最早可观测的时刻；tanh(·)表示双曲正切函数；tlj表示聚类点子集合cj的最晚可观测时刻；δ表示防止奇异运算的常值项；表示从聚类点子集合ci到下一次能够访问的聚类点子集合的集合；qim表示从聚类点子集合ci到下一次能够访问的集合中任意可能子集合cm的转移概率；α、γ、β、λ、χ为权重系数；

步骤3.8：记录蚂蚁在搜索路径上留下的信息素、观测序列、以及观测各个聚类点子集合的时刻，并更新当前时刻返回步骤3.5；

步骤3.9：计算每只蚂蚁的收益，找到收益最大的蚂蚁，并将其与当前全局最大收益进行比较，获得此次迭代的全局最优观测轨迹，更新两个聚类子集合之间的信息素，返回步骤3.2；

步骤3.10：结束感知优化模型求解。

本发明与现有技术相比的有益效果是：

(1)本发明方法建立了一种载荷凝视模式下多目标点快速聚类模型，该模型无需考虑目标可见时间窗口、航天器平台位置等约束，仅需要知道目标点位置和航天器载荷的视场范围，给出聚类条件并生成连通图，基于图论理论完成了任务聚类，计算量小，更加符合实际工程需求。

(2)本发明方法与现有的多目标观测轨迹感知方法相比，在多目标点快速聚类的基础上，对传统蚁群算法进行了改进，引入了启发因子提高了智能感知优化算法的寻优能力，降低了陷入局部极小值的概率。

(3)本发明的方法根据目标点之间的距离和航天器载荷视场，基于图论理论建立了多目标点聚类模型，在此基础上，建立了观测轨迹感知优化模型，并基于改进的智能蚁群算法得到了最优观测轨迹。

附图说明

图1为本发明方法中多目标点的聚类建模流程；

图2为本发明方法中观测轨迹的智能感知优化算法流程；

图3为本发明方法中聚类后多目标点的智能感知观测轨迹；

图4为本发明方法中聚类后多目标点的智能感知轨迹综合观测收益；

图5为本发明方法中未聚类的多目标点观测轨迹。

具体实施方式

结合附图和实施例对本发明进行进一步说明。

如图1、图2所示，一种适用于载荷凝视成像的多目标观测轨迹智能感知方法，其特征在于，包括步骤如下：

(1)建立所需观测多目标点的聚类模型，聚类模型以聚类点集合c的形式进行描述，在聚类点集合c中包含了多个聚类点子集合，在每个聚类点子集合中包含了个数不等的目标点；

具体步骤如下：

步骤1.1：对参数进行初始化设置，参数包括：任务点集合v、聚类点集合c、联通点集合l，进入步骤1.2；

步骤1.2：判断任务点集合v是否为空集；若任务点集合v为空集，则完成聚类模型建立，进入步骤(2)；若任务点集合v不为空集，则令聚类点集合c为空集，将任务点集合v中的任意点k加入到聚类点集合c中，并从任务点集合v中删除点k，进入步骤1.3；

步骤1.3：搜索任务点集合v中剩余的点中与聚类点集合c中所有节点都连通的点，建立连通点集合l，进入步骤1.4；

若点k与点l满足如下公式，则判定点k与点l之间是连通的：

其中：和分别表示点k与点l在地心惯性坐标系中的位置，rs、re和αp分别表示航天器平台飞行的轨道半径、地球半径和航天器载荷的视场角；

步骤1.4：判断连通点集合l是否为空集；若连通点集合l为空集，则进入步骤1.5；若连通点集合l不为空集，则将点l加入到聚类点集合c中，并将点l从任务点集合v中删除，重复步骤1.4直至任务点集合v中不再存在与点k连通的点，进入步骤1.5；

步骤1.5：得到聚类点集合c，返回步骤1.2；

(2)建立观测轨迹的感知优化模型，考虑约束条件的限制并设计优化目标，从聚类点集合c中选取满足条件的聚类点子集合，所建立的观测轨迹感知优化模型，具体形式如下：

其中：si和fij为智能感知优化模型的待优化变量，其计算公式分别为

n表示聚类点集合c中所包含聚类点子集合ci的个数，为正整数；i＝1,…,n；j＝1,…,n；

wi表示聚类点子集合ci的观测收益、k1表示权重系数、trij表示从聚类点子集合ci到聚类点子集合cj航天器平台姿态机动的时间、tmax表示航天器平台中执行机构能够提供的最大输出力矩；

tei表示聚类点子集合ci最早可观测的时刻、tsi表示聚类点子集合ci的实际开始观测时刻、tli表示聚类点子集合ci最晚可观测的时刻、di表示聚类点子集合ci需要持续观测的时长、tsj表示聚类点子集合cj的实际开始观测时刻；

a表示航天器平台姿态机动的角加速度、amax表示航天器平台所能达到的最大角加速度、rs表示航天器平台的位置矢量、σij表示从聚类点子集合ci到聚类点子集合cj所需机动的姿态角度；μ为地球引力常数；

聚类点子集合ci和cj之间没有交集；

(3)通过改进的智能蚁群算法进行感知优化模型的求解，得到适用于载荷凝视成像模式的敏捷航天器多目标观测的最优轨迹，具体步骤如下：

步骤3.1：设置最大迭代次数及蚂蚁数量；

步骤3.2：根据最大迭代次数，判断是否完成所有迭代，若完成则进入步骤3.10结束感知优化模型求解；否则，进入步骤3.3；

步骤3.3：将所有蚂蚁放置在感知优化开始时刻的星下点位置，并清空禁忌表、时间记录表和路径记录表，进入步骤3.4；

步骤3.4：判断是否还有未参与规划的蚂蚁，若有，则进入步骤3.5，否则，进入步骤3.9；

步骤3.5：令每一只蚂蚁进行路径搜索，并记录其已经访问过的聚类子集合和尚未访问的聚类子集合；对于尚未访问过的聚类子集合，计算从当前位置及时刻转移至该子集合的时间，并进入步骤3.6；

步骤3.6：根据步骤3.5所计算的时间，判断是否存在满足时间窗口约束要求的集合，若不存在则该蚂蚁完成一次完整的路径搜索，返回步骤3.4；否则进入步骤3.7；

步骤3.7：基于所设计的转移概率以聚类点子集合ci为出发点，来选择下一个需要观测的聚类点子集合cj，并进入步骤3.8；

其中：表示第k2只蚂蚁从聚类点子集合ci到聚类点子集合cj的转移概率；τij表示聚类点子集合ci到聚类点子集合cj的信息素浓度；tej表示聚类点子集合cj最早可观测的时刻；tanh(·)表示双曲正切函数；tlj表示聚类点子集合cj的最晚可观测时刻；δ表示防止奇异运算的常值项；表示从聚类点子集合ci到下一次能够访问的聚类点子集合的集合；qim表示从聚类点子集合ci到下一次能够访问的集合中任意可能子集合cm的转移概率；α、γ、β、λ、χ为权重系数；

步骤3.8：记录蚂蚁在搜索路径上留下的信息素、观测序列、以及观测各个聚类点子集合的时刻，并更新当前时刻返回步骤3.5；

步骤3.9：计算每只蚂蚁的收益，找到收益最大的蚂蚁，并将其与当前全局最大收益进行比较，获得此次迭代的全局最优观测轨迹，更新两个聚类子集合之间的信息素，返回步骤3.2；

步骤3.10：结束感知优化模型求解。

实施例：

所需观测目标范围为：东经115°至东经125°(e115°～e125°)，北纬25°至北纬35°(n25°～n35°)的50个密集目标点。航天器平台的轨道参数、姿态机动和载荷性能具体如下表所示。

表1相关参数设置

改进智能蚁群算法的相关参数设置为：最大迭代次数40次，蚂蚁数量30只，转移概率各项权值α＝2，β＝1，γ＝1.5，λ＝1.5，χ＝3，优化指标权值k＝10^-6。目标点的观测收益wi从0～1之间随机选取。各目标点需要观测时间为10s。图3～图5分别给出了聚类后多目标点的智能感知观测轨迹、综合观测收益、以及未聚类的多目标点观测轨迹。从实施例的结果中可以看出，本发明方法获得的观测收益为19.8704，未采用本发明方法获得的观测收益为10.9948；这说明本发明方法在提高观测效率、减少航天器平台姿态机动角度，降低在轨运行成本等方面具有比较优势。

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。