一种航天器交会与接近路径规划算法.pdf

第44卷第8 期2023年8 月宇航学报Journal of AstronauticsVol.44No.8August2023一种航天器交会与接近路径规划算法黄成，王涛，许家忠（哈尔滨理工大学自动化学院，哈尔滨150 0 8 0）摘要：针对非合作交会与接近任务中近程导引段航天器的三维路径规划问题，提高其在多约束条件下的寻优性能，提出一种融合混沌映射、自适应随机动态扰动、随机探索机制以及反向学习策略的混合蝙蝠算法。首先，利用Circle混沌映射进行种群初始化，提高种群多样性；在全局搜索中引入自适应随机动态扰动系数，并在位置更新时进行位置限定，提高算法全局寻优的性能；在局部搜索中融入柯西分布的逆累积分布函数，避免算法进入局部最优。同时，引人分段随机的反向学习策略,进一步提高种群多样性，平衡全局寻优与局部探索的性能。然后，以基于升序排列的路径节点进行三次样条插值编码，构造躲避障碍物求解航天器三维路径的适应度函数，结合混合蝙蝠算法求得平滑路径曲线。最后，通过对比仿真表明了所提方法的最优路径次数较对比算法增加，最差路径次数减少，路径平均解减少。关键词：航天器；非合作交会；蝙蝠算法；混沌映射；反向学习；自适应随机探索中图分类号：V249D0I:10.3873/j.issn.1000-1328.2023.08.011A Path Planning Algorithm for Spacecraft Rendezvous and Approach文献标识码：A文章编号：10 0 0-132 8(2 0 2 3)0 8-12 2 5-13HUANG Cheng,WANG Tao,XU Jiazhong(School of Automation,Harbin University of Science and Technology,Harbin 150080,China)Abstract:Aiming at the three-dimensional path planning problem of the spacecraft in the short range guidance phase ofthe non-cooperative rendezvous and approach mission,to improve its optimization performance under multiple constraints,ahybrid bat algorithm integrating chaotic mapping,adaptive random dynamic disturbance,random exploration mechanismand reverse learning strategy is proposed.Firstly,Circle chaotic mapping is used for population initialization to improvepopulation diversity.The adaptive stochastic dynamic perturbation coefficient is introduced into the global search,and theposition is limited when the position is updated to improve the performance of the global search algorithm.The inversecumulative distribution function of Cauchy distribution is incorporated into local search to avoid local optimization.At thesame time,the segmented random reverse learning strategy is introduced to further improve the population diversity andbalance the performance of global optimization and local exploration,Then,the ascending path nodes are coded with cubicspline interpolation to construct the fitness function of obstacle avoidance to solve the three-dimensional path of thespacecraft,and the smooth path curve is obtained by combining the hybrid bat algorithm.Finally,comparison simulationshows that the optimal path time of the proposed method is increased,the worst path time is reduced and the average pathsolutions are reduced.Key words:Spacecraft;Non-cooperative trade fair;Bat algorithm;Chaotic mapping;Reverse learning;Adaptiverandom exploration收稿日期：2 0 2 3-0 1-0 1；修回日期：2 0 2 3-0 4-16基金项目：国家自然科学基金（52 10 2 455）12260 引 言空间交会与接近任务是指追踪星抵近飞行至另一个空间目标附近,并完成诸如绕飞、在轨操作、对接等操作的整个过程，根据目标星能否提供导航信息，可以将其分为合作与非合作两类。相比合作式情形,非合作目标的交会与接近在空间监测与服务、深空探测等领域具有更广阔的应用前景，但由于缺乏信息沟通，存在更大的挑战性。为了执行复杂多样的非合作目标近距离操作任务，航天器路径规划成为该研究领域的一项关键技术。目前约有19 0 0颗在轨航天器与150 0 0 枚空间碎片存在于外太空1,采用有效的航天器路径规划技术可以保证在不发生碰撞的情况下找到一条从当前位置到达目标位置的最优路径，同时满足消耗燃料最少、飞行时间最短等2 要求。否则,航天器在飞行过程中可能发生碰撞事件或出现燃料不足等问题，导致操作任务不能正常进行。目前路径规划方法大致分为两类，一类是以可视图法、单元分解法3、随机采样法4、人工势场法等为代表的传统算法,另一类则是以A算法5-6、遗传算法7-8、粒子群算法”】、蚁群算法等为代表的启发式算法。文献10 为解决无人驾驶飞机在避障条件下的路径规划问题，提出了一种实时避碰算法,该算法所规划出的路径能够避免碰撞,并且能够保证最小安全距离。文献11针对传统人工势场法存在局部极小值以及目标不可达的问题，提出一种增加逃逸力因数的人工势场法。文献12 基于粒子群优化方法提出了大角度姿态机动时航天器在多种约束下的路径规划算法。文献13针对空间在轨操控机器人的路径规划问题，提出一种基于拓邻域搜索的蚁群算法，该算法可以有效地规划出机器人全局路径。然而，随着环境系统复杂性及任务难度的增加,传统算法难以取得理想的效果,而常规启发式算法也存在局部最优、收敛速度慢以及效率低等问题。Yangl141于2 0 10 年基于群体智能提出蝙蝠算法（Ba t a l g o r i t h m，BA），这是一种可以搜索全局最优解的启发式优化算法，与其他群智能优化算法相比结构相对简单、全局搜索能力较强、参数较少，因此更适用于解决要求运算速度更快、环境适应性更强、宇航学报性能更优的路径规划问题。文献15提出一种在蝙蝠算法中融入黄金正弦算法解决移动机器人路径规划问题的方法，该方法具有较快的收敛能力以及全局搜索能力。文献16 针对水下机器人避障的路径优化问题,通过线性渐变、高斯柯西变异以及界限随机重置机制等策略来改进蝙蝠算法，该方法使得水下机器人规划的路径更短、更平滑。文献17 提出了一种基于改进蝙蝠算法的移动机器人路径规划方法，该方法平衡了算法的全局寻优与局部探索能力,提高了算法的收敛速度与寻优精度。蝙蝠算法也常用于求解连续型的优化问题，在调度问题18、函数优化19、故障诊断2 0 图像识别等方面有着广泛的应用。然而,标准蝙蝠算法存在易陷入局部极值、后期收敛速度慢等缺点，限制了其在复杂、快速收敛性任务中的实际应用，为了改善上述缺点，拓展其应用领域，需要对标准蝙蝠算法进行优化。普遍情况下,研究者常采用传统直线或者圆弧分段插值的方法来解决路径规划的优化问题，但这些方法构造出来的移动路径会产生较多的转折点，因此能够导致路径曲线出现连续性以及平滑性较差的问题。而采用三次样条插值方法对目标路径进行优化，可以避免路径产生较多的转折点，从而拟合出一条连续且平滑的曲线。因此本文以文献17 中的方法为理论基础,提出一种新的改进蝙蝠算法,进一步提高算法的寻优能力，并结合改进编码方式的三次样条插值算法对航天器路径进行三维规划。综合考虑近程导引段任务特点及其对路径规划算法寻优性能的需求,本文将优化后的蝙蝠算法和三次样条插值方法相结合，提出一种收敛速度更快、路径更短的航天器路径规划方法。主要创新点为：1)引人Circle混沌映射进行种群初始化,有效地解决了初始解聚集等问题；2)位置更新时进行位置限定，同时在全局搜索阶段引人自适应随机动态扰动系数,在局部搜索阶段融入柯西分布的逆累积分布函数,实现了算法权重的动态变化,并增加了种群变异机制;3)融合分段随机的反向学习策略，调节了算法的搜索范围,进一步扩大了蝙蝠种群的多样性。1混合蝙蝠算法 CPTDBA1.1相对运动参考坐标系定义地心惯性坐标系OXYZ、轨道坐标系0 1x1y131第44卷第8 期以及目标星本体坐标系0 2 x2y232之间的关系如图1所示。其中地心作为地心惯性坐标系的原点0，X轴在赤道平面内，并且指向春分点；Z轴沿着地球自转轴并指向北极；Y轴满足右手定则。轨道坐标系原点0 1为参考星的质心，轴沿地心指向参考星质心；y轴在参考星所在的轨道平面内,与轴相垂直,并且指向运动方向；z轴通过右手定则进行确定。目标星本体坐标系原点0 2 为目标星的质心，3个正交惯性的目标星主轴分别为2，y 2，z 2。Z地球X图1坐标系定义Fig.1Coordinate system definition为了方便计算数据，本文以轨道坐标系作为航天器的相对运动参考坐标系。该坐标系中的参考星可以是真实的航天器，也可以是人为虚构出来的航天器，因此本文以路径规划开始时的航天器作为该坐标系的参考星。1.2蝙蝠算法数学建模算法中的声波频率更新式为：f;=fmin+(fmx-fmn)速度更新式为：VtI=Vi+(Xt+I-X,)f.:全局位置更新式为：XiI=X+V*式中：f.为第i只蝙蝠的频率；fmin与fmax分别为频率的最小值与最大值；为随机数,在0,1之间均匀分布；V+为第t+1代第i只蝙蝠的速度；X+为第t+1代第i只蝙蝠的位置；X，为蝙蝠当前最优的位置。局部搜索的位置更新式为：Xnew=Xold+sA黄成等：一种航天器交会与接近路径规划算法式中：Xnew为蝙蝠个体的新位置；Xold为当前蝙蝠个体的位置；E-1,1,且为一个随机数；A是当前蝙蝠种群所有个体的平均响度。声波响度更新式为：A+=A;脉冲发射频率更新式为：+l=r(1-e-)式中：E（0,1）是响度衰减系数；r*是t+1代第i只蝙蝠的脉冲发射频率；r是第i只蝙蝠的初始脉冲发射频率；0 是脉冲发射频率增强系数。目标星1.3Circle混沌映射初始化种群在标准蝙蝠算法与文献17 提出的PTRBA算法中,初始蝙蝠种群的选取存在随机且分布相对不均匀的问题，因此容易产生全局搜索能力不强，收敛速度较慢2 1以及陷人局部最优解的问题。为了解参考星决这一问题,提高算法的搜索能力,本文将Circle混Y沌映射算子引人到蝙蝠算法中。原Circle混沌映射的表达式为Zn+1=mod(zm+a2式中：=0.2，b=0.5,但这种Circle混沌映射值的分布依然不均匀,在0.15,0.6 之间的取值较为集中,如图2(b)所示。因此本文对原Circle混沌映射进行改进，得到带有改进Circle混沌映射的种群初始化表达式：X,=mod(3.85X,-1+0.4.式中：X，代表第n只蝙蝠个体的初始位置，n=1，(1)2,Npop-1,，Np o p 为蝙蝠的种群数量。由图2 对比可知，（a）（b）图中原Circle 混沌(2)值在0.15,0.6 之间分布的较为集中,而（c）（d）图中改进Circle混沌值分布的相对均匀,因此,引人(3)改进Circle混沌映射对种群进行初始化有利于扩大蝙蝠种群的多样性,从而提高蝙蝠算法的寻优性能，避免算法过度早熟的问题。1.4自适应随机动态扰动系数在标准蝙蝠算法中,利用全局位置更新式（3）：X=Xi+V*，在全局搜索阶段进行蝙蝠个体位置更新时，容易使得蝙蝠算法陷入局部最优。针对(4)这一问题,本文提出一种将自适应随机动态扰动系1227(5)(6)6sin(2Tzm0.7(7)(8)1228数引人到全局位置更新式的方法。该方法进一步扩大了蝙蝠种群捕捉目标的范围,增强了算法的全局搜索能力。自适应随机动态扰动系数入计算式为：W=exp(9)max入=W+$B(bi,b2)式中：t是当前蝙蝠种群的迭代次数；Tmax是蝙蝠种群的最大迭代次数；W为自适应权重策略协调收敛因子，能够平衡全局寻优和局部探索的能力，提高算法的搜索精度；=0.1,0.5上均匀分布的随机数；B(b1,b2）表示服从贝塔分布的随机数，本文中取b1=1，b=2。加人服从贝塔分布的随机数后，蝙蝠算法有机会取得较大的权值，能够使得权重进行动态的变化,从而使蝙蝠算法的收敛速度得到提高。因此在全局位置更新式中引入自适应随机动态扰动系数，可以提高算法全局寻优的能力。全局位置更新式修改为宇航学报1.01601400.81201000.6800.460400.2200500(a)原Circle混沌映射值1.00.80.60.40.20第44卷10001500维度5001000维度(c)改进Circle混沌映射值Fig.2Circle chaos map value distribution2.00015002.000图2 Circle混沌映射值分布图1.5柯西分布的逆累积分布函数标准蝙蝠算法在局部位置更新时，由于位置更新式(4):Xnew=Xold+A中的=-1,1为随机数,因此蝙蝠种群在局部搜索时具有随机性,不容(10)易寻找到局部最优解。本文采用柯西分布的逆累积分布函数对蝙蝠种群进行变异,利用柯西分布具有“尾巴”长的特点，扩大蝙蝠个体的变异范围。同时,当蝙蝠种群进行柯西逆累积分布函数变异时,蝙蝠个体将采用螺旋行走的方式进行局部寻优，因此能够避免蝙蝠种群进行盲目变异,从而提高局部搜索的能力。柯西逆累积分布函数式（12)所示,蝙蝠局部搜索的位置更新式修改为式（13）。F-l(p;xo;)=xo+tan(Xnew=Xold+A tan00.10.20.30.40.50.60.7混沌值(b)原Circle混沌映射直方图1401201008060402000.10.20.30.40.50.60.70.80.91.0混沌值(d)改进Circle混沌映射直方图XtI=XI+V+(TX0.80.91.02(12)(13)2(11)第8 期式中：Xnew为蝙蝠个体新位置；Xold为当前蝙蝠个体位置；当前蝙蝠种群所有个体的平均响度为A；K E0,1均匀分布的随机数,所以tan（k-1/2）具有正切随机探索的能力，可以使蝙蝠个体在当前局部最优解附近进行搜索，并且使得算法跳出局部最优,从而避免蝙蝠算法盲目搜索，使局部寻优探索机制更加准确。1.6分段随机反向学习策略如果蝙蝠算法在解空间中陷人了局部最优，那么算法在后期进行寻优时精度就会降低,因此,本文提出一种分段随机的反向学习策略。文献2 2 中定义了在n维空间中，点P(x1,x2，,x,）的反向解为：x=a;+b;-x;式中：x;a i,b,，i=1,2,3,n。结合式（14），得到如式（15）、（16）所示的分段随机反向解：X=(a;+b;-X),X;=ai,c;(15)X=a;+b-X,X,=c,b,(16)式中：X为蝙蝠个体最优解位置；X,为当前蝙蝠个体位置；=0.1,0.5上均匀分布的随机数；c;=（a;+b,）/2。通过采用反向学习策略对蝙蝠种群位置进行调整，进一步扩大种群的多样性，提高算法寻优精度与收敛速度，平衡全局寻优与局部探索的能力，使算法更快地跳出局部最优。2三次样条插值函数2.1平面三次样条插值设区间,b内的n+1个点,把区间分为n个形如（o,x），（x 2,x）,（x n-1,x n)的小区间,其中=，=b。每个小区间内的曲线都是1个三次多项式，而三次样条插值函数S（x）满足以下条件：1)S(x)=yi,i=0,1,2,.,n;2)在每个小区间x;,i+1内，S(）=S(）都是1个三次方程；3)S(x),S(),S(x)连续。三次样条插值函数在每段小区间内的形式可构造成如下形式：S,(x)=a,x+b,x*+c,x+d;,i=0,1,2,n-1从S（）形式可以看出,每个小区间内有4个未知数（ai,bi,C;,d,）。而区间,b内共有n个黄成等：一种航天器交会与接近路径规划算法小区间,所以整个三次样条插值函数（17)共有4n个待定系数。因此求解这些待定系数，就需要构造4n个方程。S(x)=ax3+bx?+cx+d首先，区间内所有的点；均满足插值条件，即式（18）,共n+1个方程。除端点外,其余n-1个点均满足式（19）,共2（n1）个方程,而加上两个端点分别满足第1个和第n个三次方程的条件,共2 n个方程。其次，n1个插值点的S(）是一阶可导且连续的，即满足式(2 0）,共n1个方程。另外，S（）也是二阶可导且连续的,即满足式（2 1）,同样地,共n-1 个方程。S(x,)=yi,i=0,1,2,n(14)S_(x,)=S(x,),i=1,2,n-1(19)S(x,)=S(x),i=1,2,n-1(20)S(x;)=S(x),i=1,2,n-1(21)此时共有4n2 个方程，再补充两个方程就可以得到S（），这两个方程常用以下边界条件来补充：1）自然边界：S（o）=S（n）=0,即指定端点二阶导数等于0；2）固定边界：S（o）=A，S（n）=B,即指定端点一阶导数；3)周期边界条件：端点处的函数值或导数值满足周期条件。通过式（18）（2 1）以及上述任意两个边界条件就可以求解出式（17）中的待定系数（a;,b;,Ci，d;）,因此便能得到三次样条插值函数的表达式S()。2.2空间三次样条插值由于本文求解的是航天器的三维路径，属于空间规划问题，需要进行空间的三次样条插值。根据投影法中,投影到两个平面上的两条曲线可以确定出唯一一条空间曲线的原理,本文采用投影法,将空间问题转化成平面问题进行解决。即利用投影法，结合求解平面三次样条插值的方法，计算形如式（2 2）的空间三次样条插值函数。2=S(x,y)=ao+aix+ax2+a+bo+biy+b2y2+b,y3(22)即把空间曲线S(x,y）投影到xoy平面和xoz平面上分别进行平面三次样条插值。本文根据式1229(17)(18)1230(17）构造出xoy平面和xoz平面上的三次样条插值函数式（2 3）、（2 4），并分别求解出对应的平面三次样条插值函数与z。y=ajx+bix*+Cix+d,2=ax+bax+cax+d2将式（2 2）（2 4)联立求解，就能得到空间三次样条插值函数S（，y），然后通过结合混合蝙蝠算法就能够规划出航天器的三维路径曲线。3航天器三维路径规划3.1粒子编码设计在三次样条插值中，每个大区间j,b,之间的连接处可称为路径节点，而每个小区间x;,xi+1之间的连接处则称为插值点。路径节点就相当于整条路径上最大的转向次数，通常情况下,3 5次转向就可以避开所有障碍物到达目标位置。粒子的编码就是路径节点的坐标位置,可以根据当前的环境进行选择。本文提出一种基于升序排列的路径节点进行三次样条插值编码设计的方法。假设通过蝙蝠算法得到j个坐标为（1,y1,z），（x j,j,z）的路径节点,并且xj，y j，z,均进行对应的升序排列。已知路径起点与终点的坐标分别为（x o,y o,z o）、（x m,m，z m），通过三次样条插值，在区间xo,xi,xj,xm、y o,y i,y j,y m 以及zo,21，,z j,m 上分别构造n个插值点，就可以得到n个形如（Xi,Yi,Z）,，（X,Y,Z）的三维插值点坐标，将起点、插值点、路径节点以及终点进行连线，便可以得到航天器最优的三维路径。通过将路径节点进行升序排列，避免了连线过程中曲线弯曲交织、弯曲程度较大等问题,减少了不必要路径的产生,从而提高了算法的运行效率。3.2适应度函数构建航天器在规划三维路径时需要满足以下条件：1)能够躲避障碍物,并且不能产生碰撞；2)航天器消耗的燃料最少、飞行时间最短，即规划的三维路径最短。本文就上述两个条件进行适应度函数的构建。适应度函数表达式如下：F=L(1+w)式中：L为将起点、插值点、路径节点以及终点进行连线的曲线长度，计算式为（2 6）。为避障系数，选宇航学报取合适的值可以剔除穿越障碍物的路径，本文通过实验对比,决定取=10 0 0。为标志变量,初始值为0,计算结果与障碍物相关。(23)L=(24)Z/(Xi.-X,)+(Y-Y,)+(Z.-Z.)(26)式中：（X,Y,Z,）为第i个插值点的坐标。为了方便计算适应度函数，本文将障碍物设置成圆形，个数为N,因此的计算形式可如下所示：n=n+Average(0.)dki=0 ki=max(d,=d=V(X,-X,)+(Y,-Yh)?+(Z,-Zh)2)(29)式中k=1,2,N,i=1,2,.,n;dh为当前路径上第i个插值点与第k个障碍圆心的距离，而d则表示n个d组成的集合，（Xrk，Yr k，Zh）为第个障碍的圆心坐标；.为标志变量数值集合；r为第k个圆形障碍的半径。当路径穿过第k个障碍物时，k0,反之，Oki=0;若整条路经均未穿过第k个圆形障碍,则集合,内的数均为0。为累加值，当整条路经均未穿过障碍时，=0,否则，0。3.3CPTDBA算法实现步骤步骤1.初始化种群规模Npop、声波频率最大值fmax和最小值fmin、速度V、声波响度A；、脉冲发射频率r；以及迭代总次数Tmax。确定路径起始点（o，y o，z o）和终止点（xm，y m，z m），设定圆形障碍物个数N、圆心（X,Y,Z.）以及半径r,并根据当前航天器所处环境中障碍物的个数确定路径节点个数以及插值点个数n。步骤2.初始化蝙蝠种群位置,第1只蝙蝠位置（Xo,Yo,Zo）进行随机初始化,其余 Npop1只蝙蝠位置利用式（8）进行初始化。步骤3.根据蝙蝠所在位置，参考分段随机反向学习计算式（15）、（16）,求出对应的反向解。步骤4将上述求得的x，y j，坐标分别进行(25)升序排列,根据式(2 2)计算出空间三次样条插值函数S,y）,并利用S(,y）求出n个插值点的坐标(Xi,Yi,Z,),(Xn,Yn,Zn)。第44卷(27)0),i=1,2,n(28)第8 期步骤5.利用式（2 6)求出路径长度L，并利用式(2 7）（2 9)计算出标志变量的值，从而判断路径是否穿过圆形障碍。步骤6.利用式(2 5)计算出适应度函数的值以及步骤3中反向解的适应度值。步骤7.比较当前蝙蝠个体位置与其对应反向解的适应度值大小,若前者大于后者,则用反向解代替当前蝙蝠个体的位置。步骤8.利用式（1）、（2）、（11）、（13）进行蝙蝠个体频率f、速度V.以及位置X，的更新,并对更新后的X,进行位置限定，将其固定在Xmin，Xma x 范围内；并利用式（5）、（6)进行蝙蝠个体响度A；、脉冲发射频率T；的更新。步骤9.排列适应度值，找出位置最优解。步骤10.在迭代次数Tmax内循环执行步骤3 8。步骤11.输出航天器最短的三维路径。开始初始化参数初始化种群位置下计算种群位置反向解计算当前位置以及反向解插值点计算路径长度以及标志变量计算当前位置以及反向解适应度值N当前位置解大于反向解适应度值？TY当前位置不变反向解代替当前位置图3CPTDBA算法流程图Fig.3Flow chart of CPTDBA algorithm4实验仿真为了保证某些参数设置的合理性，保证后续实验的正常进行，本文通过三种对比试验，确定了标志变量的取值、限定位置Xmin，Xmx 的范围以及路径节点升序排列的可行性。同时，为保证提出的黄成等：一种航天器交会与接近路径规划算法CPTDBA算法在求解航天器三维路径问题上的有效性和准确性,本文将 CPTDBA 算法与文献17 中提出的具有反向学习和正切随机探索的蝙蝠算法PTRBA、标准蝙蝠算法BA以及标准粒子群算法PSO,在简单和复杂两种环境下进行对比实验。4.1实验环境及参数设置为了减少其他因素对本文实验的干扰,CPT-DBA、PT RBA、BA、PSO 4种算法均在Python3.9软件平台进行仿真运行。本文分别在简单环境与复杂环境中对航天器三维路径进行对比实验,其中部分参数值的设置保持一致，如4种算法的种群数量Npop=150，迭代总次数Tmax=100,3种蝙蝠算法的声波频率最大值fmax=5kHz、最小值fmin=0kHz,声波响度A,=0.5dB,脉冲发射频率r;=0.25kHz,响度衰减系数=0.3,脉冲发射频率增强系数=1。为了确定对航天器三维路径规划影响较大的参数,即确定标志变量w的取值、限定位置Xmin，Xmax的范围以及路径节点升序排列的可行性,提高算法寻优能力，找到符合条件的最短路径，本文分别对其进行了相应的对比实验。而为了使对比结果更具有代表性以及说服力，本文将在复杂环境下以CPTDBA算法为例进行实验，实验结合路径曲折程计算频率、速度并度、种群收敛速度以及运行30 次后路径平均解、最进行位置更新短路径解、最长路径解等因素进行综合分析,从而确计算响度、发射频率定出相关参数的取值。找出当前最优解1)标志变量Q令=10，=10 0，=10 0 0,这3种取值都N送代次数最大？IY得到位置最优解输出三维路径结束1231能规划出航天器最短路径,且路径长度差距不大,但=10 0 0 时路径不存在碰撞，路径曲折程度最小，收敛速度最快，最优解和陷入局部最优的次数也最小。因此,本文设定=10 0 0。2)限定位置Xmin,Xmax分1，7，0.5，7.5，0.1，7.9 3种情况进行对比,其中3种位置限定均能规划出无碰撞的航天器三维路径,但0.1，7.9迭代速度较快,陷人局部最优的次数较少，最优解分布较为均匀，路径曲折程度、路径平均解以及最短路径解也均优于其他两种情况。因此，本文设置限定位置的范围为0.1,7.9。3）路径节点升序排列的可行性分为路径节点排序和路径节点未排序两种路径1232情况,其中两种路径情况都能规划出符合条件的航天器三维路径,但路径节点排序时的路径曲折程度小,迭代速度快，且路径平均解、最短路径解以及最长路径解都要优于未排序时的结果。因此,本文对插值的路径节点进行升序排列。4.2单目标系统下的三维路径规划单目标系统分别在简单环境和复杂环境两种情况下，对4种算法进行三维路径规划的对比实验，设CPTDBA.PSO宇航学报定路径起始点为（0，0，0）,终止点为（8，8，8）。在简单环境下设置空间障碍物个数N=6，路径节点个数j=3,在复杂环境下设置空间障碍物个数N=11，路径节点个数j=4，且两种情况下三次样条插值点个数n均为10 0。由上述确定参数的对比实验可知,标志变量=10 0 0，限定范围为0.1，7.9,且路径节点进行升序排列。两种路情况下的对比实验结果如图4及表1所示。PTRBABA起始点终止点第44卷CPTDBA-PTRBA.-PSO起始点终止点BA80765432100765432-08U/zU/z81234362346y/m788034/m567801633.452y/m(a)简单环境下三维路径15.014.814.614.4(b)复杂环境下三维路径CPTDBA15.2PTRBA.BAPSOCPTDBAPTRBABA15.0PSO14.814.614.4上14.20102030405060708090100送代次数(c）简单环境下送代曲线图4两种路径情况下4种算法的路径对比图Fig.4Path comparison of four algorithms under two path conditions表1两种路径情况下4种算法的最优路径情况对比Table 1 Comparison of optimal path conditions of fouralgorithms under two path conditions路径最优路径最路径情况算法平均解/mCPTDBA14.296 6PTRBA14.340 3简单环境BAPSOCPTDBAPTRBA复杂环境BAPSO-01020304050607080 90100送代次数(d)复杂环境下送代曲线由图4（a）（b)可知,4种算法在两种路径情况下都能够规划出航天器最短三维路径，但CPTDBA算法规划出的路径曲线更为光滑，路径的路径曲折程度相对较小，特别是在复杂环境下该结果尤次数差次数1416614.401 4414.349 0614.393 11414.455 5614.497 6414.524 461为明显。从图4(c）（d)中可以明显的看出,在两种环17境下,本文提出的混合蝙蝠算法CPTDBA收敛速度6均要优于另外3种算法，并且最优解也均小于另外13种算法，这是因为在全局搜索阶段加入的自适应6随机动态扰动系数可以使算法的权重发生动态变9化,从而提高算法收敛速度以及寻优能力。在图414（c）（d)中也可以看出CPTDBA算法能够快速逼-第8 期近最优解，并且陷人局部最优解的时间较短，能够快速跳出局部最优，这是因为在局部搜索阶段引入的柯西分布的逆累积分布函数避免了算法盲目搜索，从而提高了局部搜索能力，使得算法能够及时跳出局部最优。在两种路径情况下，每种算法独立运行30 次后,CPTDBA算法的最优解分布较为均匀,在简单和复杂两种环境下的变动范围分别为（0，0.15）与（0，0.35），而另外3种算法的最优解分布变动范围较大，最小变动范围也在(0，0.5)之间,在复杂环境下该结果更加明显。这是因为CPTDBA算法利用Circle混沌映射的方法对种群进行了初始化，避免了初始种群分布不均匀的问题,扩大了初始种群的多样性，使得算法较为稳定,进一步提高了算法的寻优性能。如表1所示，每种算法独立运行30 次后，两种路径情况下CPTDBA算法的路径平均解分别为14.2966m、14.3931m，要小于另外3种算法，并且路径最优（最优解小于其他算法)次数以及路径最差（最优解大于其他算法）次数均要优于另外3种算法。这是因为采用的随机分段的反向学习策略进一步扩大了种群的多样性,提高了寻优精度与收敛速度,并且平衡了全局寻优与局部探索的能力,因此能够让算法更加稳定。4.3多目标系统下的三维路径规划多目标系统的三维路径在规划时不仅要求航天器不能与空间障碍物发生碰撞,并且航天器与航天一机器人1-机器人2-机器人3机器人1-机器人2-机器人3起始点终止点起始点终止点黄成等：一种航天器交会与接近路径规划算法器之间也不能发生碰撞。针对上述要求，本文基于虚拟障碍物2 3思想,提出一种以插值点构造虚拟障碍物进行航天器多目标三维路径规划的方法，该方法具体流程如下：步骤1.首先利用上述算法规划出第1个航天器的三维路径，然后利用第1条路径的插值点,从第1个插值点开始,以m个插值点之间的距离为直径，并以个插值点间的中心坐标为圆心，人工虚拟出第1个空间圆形障碍,其余nm 个插值点以此类推，可虚拟出M个障碍。步骤2.此时环境中共有N+M个障碍,在此环境下再规划出第2 个航天器的三维路径。由于第1个航天器路径已经人工虚拟成障碍物,所以规划出的第2 条路径不会发生碰撞。步骤3.基于上述思想，以此类推,规划出多目标系统中所有航天器的三维路径。多目标系统同样也是在简单和复杂两种环境下对4种算法进行对比实验，不同的是多目标系统在两种环境下又根据起止点的异同分为四种情况，因此每种算法都有8 种运算环境。同理,算法中标志变量为=10 0 0，限定范围为0.1，7.9,路径节点进行升序排列。在简单环境下设置空间障碍物个数N=4,在复杂环境下设置空间障碍物个数N=9，其余参数与单目标系统相同。本文只给出复杂环境情况下的实验结果图，简单环境情况下只给出表格结果。4种算法的实验结果如下所示：机器人1-机器人2-机器人3起始点终止点1233一机器人1-机器人2-机器人3起始点终止点786428644/my/m68(1)CPTDBA一机器人1-机器人2-机器人3起始点点终止点864208264X/m6(2)PTRBA机器人1-机器人2-机器人3起始点终止点864N200826442X/m6J/m80y/m(1)CPTDBA机器人1-机器人2-机器人3起始点终止点/864200242X/m680(2)PTRBA机器人1-机器人2-机器人3起始点终止点y/mX64三N282644x/m68y/m(3)BA(a)相同起点相同终点864282644x/m268y/m(4)PSO864u/z2082644x/m268y/m(3)BA(b)相同起点不同终点8642008242X/m68y/m(4)PSO21234机器人1-机器人2-机器人3起始点终止点宇航学报机器人1-机器人2-机器人3机器人1-机器人2-机器人3起始点终止点第44卷机器人1-机器人2-机器人3起始点终止点起始点终止点864U/z20024x/m68864200862442x/m6y/m80y/m864200862442X/m686200862442/m6y/m8y/m642(1)CPTDBA一机器人1-机器人2-机器人3起始点终止点(2)PTRBA机器人1-机器人2-机器人3一起始点终止点(1)CPTDBA机器人1-机器人2-机器人3起始点终止点(2)PTRBA机器人1-机器人2-机器人3一起始点终止点86420082644x/m680y/m864200826442X/m680y/m864200826442X/m686200242x/m6y/m642J/m(3)BATable 2Comparison of optimal path conditions of the four algorithms under complex paths起止点机器人1算法状态CPTDBA相同起点PTRBA相同终点BAPSOCPTDBA相同起点PTRBA不同终点BAPSOCPTDBA不同起点PTRBA相同终点BAPSOCPTDBA不同起点PTRBA不同终点BAPSO(4)PSO(c)不同起点相同终点图5复杂路径下多目标三维路径图Fig.5Multi-target 3D path diagram under complex paths表2 复杂路径下4种算法的最优路径情况对比机器人2平均解/m平均解/m14.315 614.495 714.354 414.575 314.413 714.547 314.463 914.679 214.292.012.666 914.351 212.732.314.419 512.712.914.416 312.852514.307 613.172.714.354 613.188 514.470 613.090 314.451