一种深空探测器自主任务规划时间约束几何处理方法技术

本发明专利技术涉及一种深空探测器自主任务规划时间约束几何处理方法，属于深空探测技术领域。本发明专利技术设计了一种验证时间约束一致性(活动变量值域能够满足所有约束)和约束几何处理方法，能够快速验证深空探测器规划过程中时间约束的一致性并处理时间约束，得到活动变量最终值域，弥补采用时间约束网处理大量活动变量时计算时间长的缺点。与基于时间约束网的时间处理方法相比较，解决同样的时间约束问题，计算时间短、效率高，更加适合实时性要求高的深空探测器。

【技术实现步骤摘要】

本专利技术涉及，属于深空探测

技术介绍
在深空环境下，深空探测器与目标距离远、飞行时间长、所处环境动态多变，传统 地面站-航天器这样的控制方法已经无法满足操作费用、实时性、通讯网络等各项要求。为 达到较高的实时性要求，在探测器中应用自主技术成为重要途径。而自主任务规划技术是 自主技术的关键技术之一。 在深空探测领域中，采用自主任务规划，需要对时间进行表示、并对时间约束进行 处理。目前，时间约束网方法是探测器自主任务规划中采用的时间表示和处理方法，并已经 应用到了探索太空的项目中，如深空一号中自治远程代理系统、美国规划、调度和约束推理 平台EUR0PA等。时间约束网方法的一个缺点是如果改变其中一个点的约束或者新加入变量 点，就要对所有点的值进行计算，若在规划过程中引入一个活动，就会引入2个变量点，2*(n +1)个约束弧(η为原来的变量点），导致计算量急剧增加，不适合快速处理时间约束。 由于时间在二维情况下能够明确表示出时间的开始点、结束点和持续时间，因此 提出利用二维坐标系对时间进行表示并用几何方法处理时间约束。该方法首先由Rit提出， 后来Pu jari、Kumari和Sattar进一步讨论了一下，其他只是简单进行了介绍，而最近 Ullberg使用了该方法表示少量的定量约束解决环境识别问题中的区间推理问题。上述研 究中虽然做了一定研究，但是没有详细给出利用二维坐标系表示时间和处理时间约束的方 法。
技术实现思路
本专利技术的目的是针对目前深空探测器自主任务规划中时间约束处理问题，为了克 服时间约束网计算量大且计算时间长的缺点，提出一种深空探测器自主任务规划时间约束 几何处理方法，是一种验证时间约束一致性(活动变量值域能够满足所有约束)和约束几何 处理方法，在任务规划活动变量众多的情况下，快速实现时间约束处理，从而满足探测器实 时性要求。 本专利技术通过设计时间约束几何处理方法实现，具体实现步骤如下： 步骤1，建立深空探测器时间规划问题模型。 深空探测器时间规划问题由探测器系统状态集合、探测器可执行活动变量集合、 活动间的约束集合、探测器初始状态以及目标状态组成，即V= {V1，v2，. . .，vm}。 其中，V={V1，V2, . . .，vm}为状态变量集合，且VeDv，Dv为探测器能达到的状态值 域，m为状态变量的数量;0= {〇1，〇2,...，on}为探测器可执行的活动变量集合，η为活动变量 的数量，任意一个探测器活动〇χ具有开始点s x、结束点ex、持续时间dx，即ox = {sx，ex，dx}; C为 探测器活动间的约束集合，此处约束为两个活动间的时间约束，即C(〇x，〇y) = "〇x Cons 〇y"，〇x与〇y为活动间约束对，Cons为活动间约束关系；s。为探测器初始状态，即当时间为Ο (相对零点）时，探测器状态变量都对应相应的值;g为目标状态，即探测器各个变量需要达 到的值。 步骤2，对步骤1中探测器活动变量及活动间的时间约束进行表示。 在处理时间约束之前，需要对探测器活动变量及时间约束进行表示。 1.对探测器活动变量及活动自身约束进行表示。 设探测器活动变量集合0= {〇1，02, . . .，Οη}中任意一个探测器活动Ox，且Ox= {sx, ex,dx}。将活动变量ox在二维坐标系下进行表示:横坐标x为开始点，纵坐标y为结束点，y = x +dx与纵坐标交点的值为持续时间。 对于探测器活动变量ox= {sx，ex，dx}，自身约束为sxe 、exE 和 dxe 。其中，slx，s2x为活动o x开始点值域的下界和上界;elx，e2x为活动ox结束点值 域的下界和上界;dl x，d2x为活动ox持续时间值域的下界和上界。如果有相对零点，则slx 2 s2x2 0;slx = s2x时，表示开始点的值域只有一个时刻；同理eld e2x 2 0，dlx2 d2x2 0。 2.对深空探测器活动间约束进行表示。 设探测器两个活动变量ox= {sx，ex，dx}和oy= {sy，ey，dy}。活动间约束包括ox开始 点对〇y开始点的约束，ss 1，ss2为约束值域的下界和上界;ox开始点对Oy结束点的 约束，sel，se2为约束值域的下界和上界;o x结束点对oy开始点的约束， esl，es2为约束值域的下界和上界;ox结束点点对oy结束点的约束，eel，ee2为约 束值域的下界和上界。（ssl，ss2，sel，se2，esl，es2，eel，ee2为任意实数）。 步骤3，根据探测器活动变量自身约束，对所有探测器活动变量值域进行自我削 减。 自我削减方法为：活动变量~={&，&，心}初始值域为^已和dxe 。利用几何表示方法，二维坐标系下开始点值域范围界限为x = slx与X =82。结束点值域范围界限为7 = 6]^与7 = 62\，持续时间值域范围界限为7 = 1+(1]^与7 =叉+ d2x。六条直线相交于12个交点。从12个交点中挑选满足条件的交点：Xe，ye ，y￡ 。剩余η个点(η < 12)，其中任意一点表示为(xn，yn)。 经自我削减后，探测器活动变量0x值域为Sxe、exe[ el'x，e2'xWPdxe ,sl ,x = min(xn),s2 ,x = max(xn) ,el ,x = min(yn),e2 ,x=max(yn) ,(11^ = min(yn_xn)，d2 'X=max(yn-xn)。其中min表示最小值，max表示最大值。 步骤4，从深空探测器活动间约束对中任意选择一对以〇\，<^) = 〃<^〇118(^〃根据(^ 值域和约束关系Cons推导〇y的值域。分别推导〇y的开始点值域、结束点值域和持续时间值 域。步骤4.1，对于Oy的开始点值域，ox和约束、对有影响。在坐标系中，为了两个约束同时满足，〇y的开始点可行区域 Availablex= Π (1) 如果交集为空，则表示没有值同时满足两个约束，整个时间约束处理过程结束，结 论为时间约束不一致，即无法得到满足所有时间约束的活动值域。如果交集不为空，则进行 步骤4.2。 步骤4.2,对于oy的结束点值域，ox和约束和对其有 影响。在坐标系中，为了两个约束同时满足，0y的结束点可行区域为 Availabley= Π (2) 如果交集为空，则表示没有值同时满足两个约束，整个时间约束处理过程结束，结 论为时间约束不一致，即无法得到满足所有时间约束的活动值域。如果交集不为空，则进行 步骤4.3。步骤4.3，对〇y的持续时间值域，由式(3)获得。 Availabled = dl Π d2 Π d3 Π d4 (3) 其中，dl = ，d2=，d3 = ，d4 = ;且 d31=min(dlx-ssl+eel，dlx-ssl+ee2，dlx-ss2+eel，dl x-ss2+ee2)， d32=max(d2x-ssl+eel，d2x-ssl+ee2，d2x-ss2+eel，d2 x-ss2+ee2)， d41=min(-dlx-esl+sel，-dlx-esl+se2，-dlx-es2+sel，-dl x-es2+se2)， d42=ma本文档来自技高网...

【技术保护点】
一种深空探测器自主任务规划时间约束几何处理方法，其特征在于：具体实现步骤如下：步骤1，建立深空探测器时间规划问题模型；深空探测器时间规划问题由探测器系统状态集合、探测器可执行活动变量集合、活动间的约束集合、探测器初始状态以及目标状态组成，即V＝{v1,v2,...,vm}；其中，V＝{v1,v2,...,vm}为状态变量集合，且v∈Dv，Dv为探测器能达到的状态值域，m为状态变量的数量；O＝{o1,o2,...,on}为探测器可执行的活动变量集合，n为活动变量的数量，任意一个探测器活动ox具有开始点sx、结束点ex、持续时间dx，即ox＝{sx,ex,dx}；C为探测器活动间的约束集合，此处约束为两个活动间的时间约束，即C(ox,oy)＝ox Cons oy，ox与oy为活动间约束对，Cons为活动间约束关系；so为探测器初始状态，即当相对零点时，探测器状态变量都对应相应的值；g为目标状态，即探测器各个变量需要达到的值；步骤2，对步骤1中探测器活动变量及活动间的时间约束进行表示；在处理时间约束之前，需要对探测器活动变量及时间约束进行表示；1.对探测器活动变量及活动自身约束进行表示；设探测器活动变量集合O＝{o1,o2,...,on}中任意一个探测器活动ox，且ox＝{sx,ex,dx}；将活动变量ox在二维坐标系下进行表示：横坐标x为开始点，纵坐标y为结束点，y＝x+dx与纵坐标交点的值为持续时间；对于探测器活动变量ox＝{sx,ex,dx}，自身约束为sx∈[s1x,s2x]、ex∈[e1x,e2x]和dx∈[d1x,d2x]；其中，s1x，s2x为活动ox开始点值域的下界和上界；e1x，e2x为活动ox结束点值域的下界和上界；d1x，d2x为活动ox持续时间值域的下界和上界；如果有相对零点，则s1x≥s2x≥0；s1x＝s2x时，表示开始点的值域只有一个时刻；同理e1x≥e2x≥0，d1x≥d2x≥0；2.对深空探测器活动间约束进行表示；设探测器两个活动变量ox＝{sx,ex,dx}和oy＝{sy,ey,dy}；活动间约束包括ox开始点对oy开始点的约束[ss1,ss2]，ss1，ss2为约束值域的下界和上界；ox开始点对oy结束点的约束[se1,se2]，se1，se2为约束值域的下界和上界；ox结束点对oy开始点的约束[es1,es2]，es1，es2为约束值域的下界和上界；ox结束点点对oy结束点的约束[ee1,ee2]，ee1，ee2为约束值域的下界和上界；ss1，ss2，se1，se2，es1，es2，ee1，ee2为任意实数；步骤3，根据探测器活动变量自身约束，对所有探测器活动变量值域进行自我削减；自我削减方法为：活动变量ox＝{sx,ex,dx}初始值域为sx∈[s1x,s2x]、ex∈[e1x,e2x]和dx∈[d1x,d2x]；利用几何表示方法，二维坐标系下开始点值域范围界限为x＝s1x与x＝s2x，结束点值域范围界限为y＝e1x与y＝e2x，持续时间值域范围界限为y＝x+d1x与y＝x+d2x；六条直线相交于12个交点；从12个交点中挑选满足条件的交点：x∈[e1x,e2x]，y∈[e1x,e2x]，y∈[x+d1x,x+d2x]；剩余n个点，n≤12，其中任意一点表示为(xn,yn)；经自我削减后，探测器活动变量ox值域为sx∈[s1'x,s2'x]、ex∈[e1'x,e2'x]和dx∈[d1'x,d2'x]，其中，s1'x＝min(xn)，s2'x＝max(xn)，e1'x＝min(yn)，e2'x＝max(yn)，d1'x＝min(yn‑xn)，d2'x＝max(yn‑xn)；其中min表示最小值，max表示最大值；步骤4，从深空探测器活动间约束对中任意选择一对C(ox,oy)＝oxCons oy根据ox值域和约束关系Cons推导oy的值域；分别推导oy的开始点值域、结束点值域和持续时间值域；步骤4.1，对于oy的开始点值域[s1y,s2y]，ox和约束[ss1,ss2]、[es1,es2]对[s1y,s2y]有影响；在坐标系中，为了两个约束同时满足，oy的开始点可行区域Availablex＝[s1x+ss1,s2x+ss2]∩[e1x+es1,e2x+es2]   (1)如果交集为空，则表示没有值同时满足两个约束，整个时间约束处理过程结束，结论为时间约束不一致，即无法得到满足所有时间约束的活动值域；如果交集不为空，则进行步骤4.2；步骤4.2，对于oy的结束点值域[e1y,e2y]，ox和约束[se1,se2]和[ee1,ee2]对其有影响；在坐标系中，为了两个约束同时满足，oy的结束点可行区域为Availabley＝[s1x+se1,s2x+se2]∩[e1x+ee1,e2x+ee2]   (...

【技术特征摘要】

【专利技术属性】
技术研发人员：徐瑞，李朝玉，崔平远，朱圣英，高艾，乔栋，尚海滨，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11