CN108919818A - 基于混沌种群变异pio的航天器姿态轨道协同规划方法 - Google Patents
基于混沌种群变异pio的航天器姿态轨道协同规划方法 Download PDFInfo
- Publication number
- CN108919818A CN108919818A CN201810366021.3A CN201810366021A CN108919818A CN 108919818 A CN108919818 A CN 108919818A CN 201810366021 A CN201810366021 A CN 201810366021A CN 108919818 A CN108919818 A CN 108919818A
- Authority
- CN
- China
- Prior art keywords
- population
- individual
- iteration
- chaos
- operator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 87
- 230000003044 adaptive effect Effects 0.000 claims abstract description 24
- 230000035772 mutation Effects 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000011218 segmentation Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 3
- 230000008929 regeneration Effects 0.000 claims description 3
- 238000011069 regeneration method Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 10
- 230000006866 deterioration Effects 0.000 abstract description 3
- 230000008602 contraction Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 15
- 230000008901 benefit Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000007423 decrease Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了基于混沌种群变异PIO的航天器姿态轨道协同规划方法，属于卫星姿态轨道控制的技术领域。本发明在算法的演化阶段采取了“探索”，“搜寻”，“变异”，“归巢”的鸽群动态寻优策略。在地图指南针阶段针对种群的初始化问题加入混沌算子进行初始化操作，在种群完成初始化后加入了自适应算子实现了种群可依据当前种群演化状态自适应演化，同时加入了变异算子来解决种群陷入局部最优解的问题；在地标算子阶段针对种群收缩问题加入收缩算子，解决了优秀个体流失过快、种群退化问题，使规划结果更加平滑，种群演化更加深入，局部最优解问题以及算法发散问题得到了解决，进一步减少了算法的计算量。
Description
技术领域
本发明公开了基于混沌种群变异PIO的航天器姿态轨道协同规划方法，属于卫星姿态轨道控制的技术领域。
背景技术
PIO(Pigeon-Inspired Optimization，鸽群启发优化)算法是受鸽群在归巢过程中的导航过程的启发发展而来的，鸽群在远离目的地的时候主要依靠太阳以及地磁进行导航，而在抵近目的地后则会使用地标进行导航，并且群体中不熟悉地标的个体将会跟随熟悉地标的群体飞行。
现有PIO算法能够实现无人机以及机器人的避障以及路径规划、图像边缘识别，但是尚未解决航天器姿态轨道在复杂约束条件以及强耦合关系下的协调规划问题。现有PIO算法存在以下四类缺陷：(1)研究发现种群状态的初始化会直接影响前期的种群演化结果，并且当前的PIO算法在种群演化时无法根据种群当前的演化状态自适应地对后面的演化进行调整，这使得算法存在易陷入局部最优解、演化停滞、提前收敛、演化不深入等问题；(2)现有PIO算法在地标算法阶段存在种群数量下降过快的问题，这使得种群在最后阶段未充分演化便将优秀个体淘汰，算法有发散甚至退化的可能；(3)在适应度函数方面，当前PIO算法无法针对规划结果的平滑程度进行筛选。(4)研究发现现有PIO算法的结果的方差较大，这使得航天器的控制机构频繁启动，消耗航天器宝贵能源的同时对机械结构也有一定的磨损。
而当前针对航天器的姿态轨道的规划算法普遍是将姿态和轨道分开进行规划的，但是航天器的姿态轨道实际上存在直接的耦合关系，在进行制导时应考虑到航天器的位置，在航天器轨道规划结果的基础上对姿态进行规划。目前，针对航天器姿态规划问题，McInness通过势函数构建了姿态限制性区域的模型并且使用了李雅普诺夫第二法则得到控制量的输入，该方法虽然计算量较低，可有效减少星载计算机的资源占用，但是并没有将航天器姿态机动的时间以及耗能考虑在内，同时，在计算过程中使用的欧拉角在运动学描述中产生奇点；Melton从最优控制的角度使用数值解法，但是其计算量较大并且可能无法适用于无障碍运动；苏抗等人从卫星低RCS特性的角度出发对卫星的姿态进行规划，但是其并未对滚转角进行规划，且俯仰角仅可在[0-180°]内规划；Kjellberg等人采用二十面体离散化技术并使用A*算法实现最优路径的规划，但是算法中并未考虑边界约束和动态约束；Kim Y等人将路径规划问题转换成半定规划问题，但该方法的复杂度会随着约束的增加而急剧增大。
本发明旨在对现有PIO算法的核心迭代算法以及适应度函数进行深度改进以实现航天器集群的姿态轨道协同规划。
发明内容
本发明的发明目的是针对上述背景技术的不足，提供了基于混沌种群变异PIO的航天器姿态轨道协同规划方法，采用改进型PIO算法实现了强耦合关系下的航天器姿态轨道的协同规划，解决了现有航天器姿态轨道规划方案尚未实现姿态与轨道的协同规划以及现有姿态规划方案或轨道规划方案未考虑全部动态约束条件和影响规划结果的全部因素且复杂度大的技术问题。
本发明为实现上述发明目的采用如下技术方案：
基于混沌种群变异PIO的航天器姿态轨道协同规划方法，采用基于混沌种群变异的自适应鸽群算法规划航天器路径，依据路径节点与航天器的姿态性强制约束以及限制性约束的相对关系建立姿态规划模型，将姿态规划模型映射到R参数空间，采用基于混沌种群变异的自适应鸽群算法规划航天器姿态得到耦合的姿态轨道规划结果。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的进一步优化方案，基于混沌种群变异的自适应鸽群算法包含以下两个阶段：
指南针算子迭代演化阶段：在每次迭代过程中分别选出当前迭代过程中的局部最优位置以及所有个体在当前迭代过程中产生的全局最优位置，并引入表示个体对鸽群整体在之前迭代中所产生的精英个体学习的全局更新算子和表示前两次迭代演化趋势的自适应算子以及对种群进行变异操作的变异算子，迭代更新个体的位置和速度，变异算子仅在种群适应度变化率达判定门限值时启动，所述变异算子启动后生成一个以个体当前位置为均值且以当前种群适应度变化率的倒数为方差的随机数，在当前迭代次数达到设定值时进入地标算子迭代演化阶段；
地标算子迭代演化阶段：依据种群适应度对种群进行筛选进而迭代更新局部最优个体和全局最优个体。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的进一步优化方案，指南针算子描述为：
Vi(t+1)分别为第i只个体在第t次迭代、第t+1次迭代中的速度，Xi(t)、Xi(t+1)分别为第i只个体在第t次迭代、第t+1次迭代中的位置，w为惯性因子，wmax和wmin为惯性因子最大值和最小值，Tmax为最大迭代次数，为第i只个体在第t次迭代中的局部最优位置，Cp为自适应算子，fitness(t-3)、fitness(t-2)、fitness(t-1)分别为种群在第t-3次、第t-2次、第t-1次迭代中的适应度值，Xg为所有个体在第t次迭代中的全局最优位置，Cg为全局更新算子，Cr为激活因子，Rd(t+1)为变异算子在第t+1次迭代中生成的随机数，fitness(t)为种群在第t次迭代中的适应度值。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的进一步优化方案，所述种群适应度由各个体的适应度确定，第i只个体的适应度fitness(i)的表达式为：
threat_inj(i)、w1为第i只个体进入威胁区域j时的威胁度及其权值，threat_outj(i)、w2为第i只个体在威胁区域j外时的威胁度及其权值，distance(i)、w3为第i只个体规划路径的整体距离消耗对整体规划距离的影响程度及其权值，angle(i)、w4为第i只个体规划路径的机动角度及其权值，w1+w2+w1+w4＝1。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的更进一步优化方案，第i只个体进入威胁区域j时的威胁度threat_inj(i)为：Li,k为第i只个体所规划路径第k个分段的长度，tj为威胁区域j的威胁等级，K为第i只个体所规划路径的分段数目，d0.1,k、d0.3,k、d0.5,k、d0.7,k、d0.9,k分别为第i只个体所规划路径第k个分段的十分之一处、十分之三处、十分之五处、十分之七处、十分之九处与威胁区域j中心的距离，Tj为威胁区域j的中心位置，Rj为威胁区域j的半径，Xi(t)为在t次迭代中第i只个***置的三维坐标。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的更进一步优化方案，第i只个体在威胁区域j外时的威胁度threat_outj(i)为：
tj为威胁区域j的威胁等级，K为第i只个体所规划路径的分段数目，d0.1,k、d0.3,k、d0.5,k、d0.7,k、d0.9,k分别为第i只个体所规划路径第k个分段的十分之一处、十分之三处、十分之五处、十分之七处、十分之九处与威胁区域j中心的距离，Tj为威胁区域j的中心位置，Rj为威胁区域j的半径，Xi(t)为在t次迭代中第i只个***置的三维坐标。
再进一步的，基于混沌种群变异PIO的航天器姿态轨道协同规划方法中，威胁区域j的威胁等级tj为： 表示从威胁区域j的中心位置指向第i-1只个***置的向量，表示从第i-1只个***置指向第i只个***置的向量。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的更进一步优化方案，第i只个体规划路径的机动角度angle(i)为：M为个体数目，表示从第i-2只个***置指向第i-1只个***置的向量，表示从第i-1只个***置指向第i只个***置的向量。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的进一步优化方案，地标算子迭代演化阶段依据表达式：更新种群数量，NP(t)为第t次迭代的种群数量，N为种群数量，Tmax为最大迭代次数。
作为基于混沌种群变异PIO的航天器姿态轨道协同规划方法的再进一步优化方案，在指南针算子迭代演化阶段开始前，采用Tent Map映射混沌初始化种群。
本发明采用上述技术方案，具有以下有益效果：
(1)本文针对航天器在复杂约束下的姿态轨道协同规划问题设计了一种基于混沌种群变异的自适应鸽群改进型算法，并且将航天器姿态规划空间映射到Rodrigues参数空间中进行规划，以满足欧拉角规划解耦的问题，使得卫星在空间中的每一个位置对应一个姿态，得到了姿态限制性的规划结果；
(2)引入的自适应算子依据种群当前状态来动态调整种群搜索的范围，在种群进化缓慢时加大搜索范围，在种群速度过快时减小搜搜范围，可提升演化深度并在一定程度上避免局部最优解，引入的变异算子在种群陷入局部最优解时使种群脱离局部最优解而在种群得到全局最优解时保留全局最优解，在地标算子阶段针对种群收缩问题加入收缩算子以更新种群数量，解决了优秀个体流失过快、种群退化问题，使规划结果更加平滑，种群演化更加深入，局部最优解问题以及算法发散问题得到了解决，算法具有更快的收敛速度，大大降低了计算量，实现了航天器姿态轨道的最优规划；
(3)适应度函数引入平滑程度评价因子以对种群演化结果的平滑程度进行筛选，使优化结果尽量平滑，减少抖动，对威胁等级的改进是针对航天器编队与威胁区域距离的进行的，航天器进入威胁区域后威胁等级十分高且上升趋势随距离的靠近而加快，可在筛选时避免选择进入碰撞范围的解，而威胁等级在航天器离开威胁区域后下降平缓，可筛选离开碰撞范围的解，对航天器的位置相要求相对松散。
附图说明
图1为使用Tent Map映射种群初始化的示意图。
图2为种群数量随迭代次数变化趋势的示意图。
图3为姿态轨道协同规划的示意图。
图4为本申请利用基于混沌种群变异的自适应鸽群算法规划航天器姿态轨道的流程图。
图5为经典款PIO算法轨道规划的仿真图。
图6为粒子群算法轨道规划的仿真图。
图7为本申请改进型混沌种群变异自适应鸽群算法轨道规划的仿真图。
图8为本申请提出的CPVAPIO算法和PIO算法及PSO算法计算成本的对比图。
图9为采用本申请CPVAPIO算法得到的姿态限制性规划结果。
图10为采用本申请CPVAPIO算法得到的航天器姿态角规划结果。
具体实施方式
下面结合附图对发明的技术方案进行详细说明。
本文针对鸽群算法的种群初始化，种群迭代演化，以及适应度函数三个方面进行了重点研究，并提针对现有遗传算法的不足，提出了基于“探索”，“搜寻”，“变异”，“归巢”的鸽群动态寻优策略以解决现有算法在收敛速度，局部最优，演化深度方面存在的问题。
(1)探索：
在“探索”阶段，种群进行初始化，作为非线性***中一种典型的现象——混沌现象，具有随机性、遍历性以及规律性等特点。在该阶段利用混沌现象的遍历性可使初始化种群尽量分散地分布在整个需要探索的空间内，同时，种群中的个体有更高的概率接近搜索空间中的全局最优解，在一定程度上改善种群初始化过于集中所导致的迭代逐步陷入局部最优解的缺陷。
于是，在种群的初始化中使用Tent Map映射混沌初始化种群：
Sn+1＝ξ*(1-2*|Sn-0.5|),n＝0,1,2...,N (1)，
其中，0＜S0＜1，随机数ξ取1时则***处于完全混沌状态，Sn为混沌空间内对应规划空间内第i只个***置的混沌个体，Sn+1为对应规划空间内第i+1只个***置的混沌个体，同时所需要初始化的空间的范围是[Xmin,Xmax]，规划的维度为D，于是，个体初始化位置为：
如图1所示，Tent Map映射可以很好的覆盖0-1的空间。
(2)搜寻：
在种群完成初始化后开始迭代对整个空间进行搜索寻优，本部分针对经典PIO算法的核心迭代算子进行了研究，种群在演化的过程中主要参考全局最优个体的位置信息来更新，而迭代初期的种群进化并不深入，全局最优个体相对普通个体并不具有十分明显的优势，从而对整个种群的参考价值并不高，与此同时，经典PIO算法在更新速度时，全局最优解的权值为一个随机数。于是，本文针对全局最优个体加入了全局更新算子Cg，使全局最优个体只有在优势明显时才具有更高的权值，而迭代初期的权值较小。
研究发现，在迭代初期，单个个体在已有迭代中的演化经验相比初期的全局最优个体具有更高的参考价值，可以达到帮助收敛的目的，而经典PIO算法并没有对该信息加以利用，于是，本文创新性地引入了自适应算子Cp，使之可以依据种群中个体当前迭代状态与前两次迭代的演化趋势进行自适应调整，在鸽群进化缓慢时增大搜索的范围，以达到动态搜索最优个体的目的。
自适应改进型地图指南针的优化算子更新公式如下：
Xi(t+1)＝Xi(t)+Vi(t+1) (4)，
其中，Vi(t)、Vi(t+1)分别为第i只个体在第t次迭代、第t+1次迭代中的速度，Xi(t)为第i只个体在第t次迭代中的位置，为第i只个体在t次迭代中的局部最优位置，Xg为所有个体在第t次迭代中的全局最优位置，fitness(t-3)、fitness(t-2)、fitness(t-1)分别为第i只个体在第t-3次、第t-2次、第t-1次迭代中的适应度值，w为惯性因子，wmax为惯性因子最大值，此处取0.9，wmin惯性因子最小值，此处取0.4，Tmax为最大迭代次数，t为当前迭代次数。
(3)变异：
局部最优解是经典PIO算法以及同类遗传算法共同存在的问题，如何抑制算法陷入局部最优解以及在算法陷入局部最优解后如何激活种群脱离局部最优解是解决此类问题的关键，于是，本文在上述优化算子的基础上进一步加入了种群变异算子。当种群适应度值的变化率连续30次低于0.01时，则种群有可能陷入了局部最优解，此时，对鸽群中距离当前局部最优解较远的一半(即按适应度值排序的后一半个体)进行变异，变异算子取一个以个体当前位置为均值且以当前适应度变化率的倒数为方差的随机数。可使后一半种群在下次迭代中被激活，进而实现大范围搜索，将种群带出局部最优解。同时，由于仅对后一半种群进行变异，如果算法已找到全局最优解，大范围搜索的结果不会优于当前已找到的全局最优解，故不会使算法发散。
最终地图指南针的优化算子如下：
其中，Cr(Rd(t+1)-Xi(t))为变异算子，Cr为激活因子，Rd(t+1)为以Xi(t)为均值且以为方差的高斯随机变量。
(4)归巢：
在结束地图指南针算子迭代演化阶段后，算法进入地标算子迭代演化阶段，地标算子在该阶段快速减少种群数量以保留全局最优解、抑制种群发散。但在经典PIO算法中，该阶段种群下降速度为每次减少一半，减少速度过快以至于超过了种群收敛的速度，从而导致种群退化。于是，本申请对种群降低策略进行了改进，使种群数量在前期快速下降筛选优势个体而在后期缓慢下降保留优势个体。
种群数量NP(t)按下式进行更新：
其中，N为种群数量，此处，取N＝30，Tmax＝50，种群数量随迭代次数的变化趋势如图2所示。
在经典PIO算法的适应度函数中，仅对规划路径的长短以及与障碍区域的距离做了定义，而在实际问题中，非威胁区域的路径规划依然对航天器的整体机动最优化程度有着重要的影响。
在经典PIO算法中，用于定义障碍区域威胁程度的tk为一变量，可针对威胁区域的不同设定不同的威胁等级，而航天器避障与飞机穿过雷达区域不同，任何与威胁区域的接触均会引起碰撞，碰撞是航天器路径规划需要规避的问题。于是，本申请将tk改进为与航天器同障碍航天器距离相关的函数，距离越近表示威胁程度越高，而在航天器远离障碍时，威胁程度迅速降低，以避免航天器不必要的机动。
同时，对于航天器机动问题，优化的轨迹应尽量平滑以减少大角度大范围的机动进而延长航天器的使用寿命，于是在适应度函数中应加入机动角度因子。
于是，在原适应度函数的基础上进行改进如下式：
其中，fitness(i)为第i只个体的适应度，threat_inj(i)、w1为第i只个体进入威胁区域j时的威胁度及其权值，threat_outj(i)、w2为第i只个体在威胁区域j外时的威胁度及其权值，distance(i)、w3为第i只个体规划路径的整体距离消耗对整体规划距离的影响程度及其权值，angle(i)、w4为第i只个体规划路径的机动角度及其权值，w1+w2+w1+w4＝1。该适应度函数，在规划空间中针对路径进入威胁区域内与在威胁区域外的情形建立不同的适应度函数，并且加入评价路径平滑程度的angle(i)。
当规划路径进入威胁区域时，其适应度函数为：
其中，Li,k为第i只个体所规划路径第k个分段的长度，tj为威胁区域j的威胁等级，K为第i只个体所规划路径的分段数目，d0.1,k、d0.3,k、d0.5,k、d0.7,k、d0.9,k分别为第i只个体所规划路径第k个分段的十分之一处、十分之三处、十分之五处、十分之七处、十分之九处与威胁区域j中心的距离，Tj为威胁区域j的中心位置，Rj为威胁区域j的半径，实际中，Rj为第j个障碍航天器的外轮廓半径，Xi(t)为在第t次迭代中第i只个***置的三维坐标。
当规划路径未进入威胁区域时，其适应度函数为：
在以上两式中，tj表征当前的威胁等级：
其中，表示从威胁区域j的中心位置指向第i-1只个***置的向量，表示从第i-1只个***置指向第i只个***置的向量。
路径平滑适应度函数为：
其中，M为个体数目，表示从第i-2只个***置指向第i-1只个***置的向量。
航天器姿态轨道协同规划方案
航天器的姿态轨道规划问题是航天领域的热点问题，但是目前已有的研究主要将姿态与轨道单独进行研究，而在实际中由于航天器轨道的变化，观测任务目标以及其它天体的相对位置也发生了变化。从而轨道的规划问题与姿态的规划问题是存在耦合关系的，应将二者协同考虑。
于是，本文首先依据航天器机动变轨时需要规避的障碍建立路径规划空间的模型，初始化障碍区域以及航天器的当前位置以及目标位置，利用本文中的算法进行航天器的路径规划。随后在已完成路径规划的基础上，依据路径点与航天器的姿态强制性约束以及限制性约束的相对关系建立姿态规划模型，并将模型映射到R参数空间中，从而再次利用本文提出的基于混沌种群变异的自适应鸽群优化方法(下文简称CPAVPIO算法)进行规划完成姿态轨道的协同规划，输出耦合的姿态轨道规划结果，协同规划的流程如图3所示。
本发明采用混沌种群变异的自适应鸽群算法实现航天器的姿态轨道协同规划如图4所示。
(1)初始化鸽群位置以及算法参数
在鸽群的位置初始化中应用Tent Map混沌理论，使鸽群的初始化位置尽量均匀且无规律的分散在整个规划空间中。并且设置鸽群数量Q、搜索维度D、鸽群的初始化位置为X0＝[x0 y0 z0]，鸽群的目标位置为XD＝[xD yD zD]，地图指南针算子迭代次数T1、最大迭代次数TMAX、第j个威胁区域的中心坐标Tj＝[xT yT zT]T、第j个威胁区域的半径Rj。
(2)启动地图指南针算子进行种群演化
为了使增加种群的演化深度，在种群演化过程中尽量避免陷入局部最优解，相比传统的算法中对全局最优解单一参考，应该使种群依据当前演化的状态来决定下一次迭代演化中对全局最优解的参考。于是，创新性地设计了如式(3)至式(7)所示的自适应改进型地图指南针的优化算子。
(3)局部最优解脱离
当前同类型算法共同存在的一个问题即为种群有可能会陷入局部最优解。而在陷入局部最优解后如何使种群脱离局部最优解，同时，使得种群不至于太活跃而使算法发散，以及在将全局最优解误判为局部最优解后使算法脱离全局最优解的问题是一个难点问题。为解决此问题，在算法中设计了变异算子。当种群适应度值的变化率连续30次低于0.01时则种群有可能陷入了局部最优解，此时，对鸽群中距离当前局部最优解较远的一半(即按适应度值排序的后一半个体)进行变异，变异算子将产生以当前值为均值并以当前适应度变化率的倒数为方差的一个随机数，以这个随机数作为参考进行后续演化，这样使得当前的解得到了保留。如果当前解为局部最优解，距离此解较远的后一半种群在变异后可将种群带离此解；如果当前解已为全局最优解，则后一半种群的演化不会产生优于此解的结果，因而保留了全局最优解。最终得到如式(8)至式(10)所示的地图指南针的优化算子。
(4)启动地标算子进行种群演化
原始算法的地标算子部分，种群减少速率过快甚至超过了种群收敛的速度，于是，本申请对种群的下降速率进行设计，使种群数量在前期快速下降以筛选优势个体而在后期缓慢下降保留优势个体，种群数量按式(11)进行更新。
(5)适应度函数
在演化的过程中依据适应度函数对种群进行筛选排序，产生局部最优解与全局最优解。
仿真实验与结果验证
实验设定当前机动航天器在运动坐标系下开始机动时位于坐标系原点(0,0,0,)KM，目标位置为(65,80,30)KM，设定变轨航天器前方障碍航天器的相对位置为(45,50,5)KM，(8,25,15)KM，(40,68,7)KM，航天器的最小安全距离为5KM。
仿真实验条件为主核四核2.8Hz、16GB内存、64位计算机，在此基础上设定算法内部参数分别为：
●经典款PIO算法与改进型PIO算法：个体数量M＝30，规划维度D＝20，迭代次数为200次(地图指南针算子150次，地标算子50次)。
●粒子群算法：个体数量M＝30，规划维度D＝20，迭代次数为200次。
●混沌种群变异的自适应鸽群算法(CPAVPIO算法)：个体数量M＝30，规划维度D＝20，迭代次数为200次。
1.1航天器轨道规划实验结果及算法性能对比
经典款PIO算法轨道规划结果如图5所示，粒子群(PSO)算法轨道规划结果如图6所示，本申请基于混沌种群变异自适应鸽群(CPVAPIO)改进算法轨道规划的结果如图7所示，CPVAPIO算法和PIO算法及PSO算法计算成本如图8所示，CPVAPIO算法和PIO算法及PSO算法计算量对比如表1所示。
表1CPVAPIO算法和PIO算法及PSO算法计算量对比
采用本申请CPVAPIO算法得到的姿态限制性规划结果如图9所示，采用本申请CPVAPIO算法得到的航天器姿态角规划结果如图10所示。
可见，本文针对航天器在复杂约束下的姿态轨道协同规划问题设计了一种基于混沌种群变异的自适应鸽群改进型算法，并且将航天器姿态规划空间映射到Rodrigues参数空间中进行规划，以满足欧拉角规划解耦的问题。通过对比实验可以看出本文的算法相比于经典PIO算法以及PSO算法，在规划结果上更加平滑，可减少不必要的姿轨调整；仿真结果表明，混沌种群变异的自适应鸽群算法(CPAVPIO算法)相比PIO以及PSO算法，种群进化更加深入，适应度下降值在10左右，而其余两种算法均在5左右，同时，算法具有更快的收敛速度；并且通过计算量的对比可以看出本文的改进型算法大大降低了计算量，实现了航天器姿态轨道的最优规划。
Claims (10)
1.基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，采用基于混沌种群变异的自适应鸽群算法规划航天器路径，依据路径节点与航天器的姿态性强制约束以及限制性约束的相对关系建立姿态规划模型，将姿态规划模型映射到R参数空间，采用基于混沌种群变异的自适应鸽群算法规划航天器姿态得到耦合的姿态轨道规划结果。
2.根据权利要求1所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，所述基于混沌种群变异的自适应鸽群算法包含以下两个阶段：
指南针算子迭代演化阶段：在每次迭代过程中分别选出当前迭代过程中的局部最优位置以及所有个体在当前迭代过程中产生的全局最优位置，并引入表示个体对鸽群整体在之前迭代中所产生的精英个体学习的全局更新算子和表示前两次迭代演化趋势的自适应算子以及对种群进行变异操作的变异算子，迭代更新个体的位置和速度，变异算子仅在种群适应度变化率达判定门限值时启动，所述变异算子启动后生成一个以个体当前位置为均值且以当前种群适应度变化率的倒数为方差的随机数，在当前迭代次数达到设定值时进入地标算子迭代演化阶段；
地标算子迭代演化阶段：依据种群适应度对种群进行筛选进而迭代更新局部最优个体和全局最优个体。
3.根据权利要求2所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，所述指南针算子描述为：Vi(t)、Vi(t+1)分别为第i只个体在第t次迭代、第t+1次迭代中的速度，Xi(t)、Xi(t+1)分别为第i只个体在第t次迭代、第t+1次迭代中的位置，w为惯性因子，wmax和wmin为惯性因子最大值和最小值，Tmax为最大迭代次数，为第i只个体在第t次迭代中的局部最优位置，Cp为自适应算子，fitness(t-3)、fitness(t-2)、fitness(t-1)分别为种群在第t-3次、第t-2次、第t-1次迭代中的适应度值，Xg为所有个体在第t次迭代中的全局最优位置，Cg为全局更新算子，Cr为激活因子，Rd(t+1)为变异算子在第t+1次迭代中生成的随机数，fitness(t)为种群在第t次迭代中的适应度值。
4.根据权利要求2所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，所述种群适应度由各个体的适应度确定，第i只个体的适应度fitness(i)的表达式为：threat_inj(i)、w1为第i只个体进入威胁区域j时的威胁度及其权值，threat_outj(i)、w2为第i只个体在威胁区域j外时的威胁度及其权值，distance(i)、w3为第i只个体规划路径的整体距离消耗对整体规划距离的影响程度及其权值，angle(i)、w4为第i只个体规划路径的机动角度及其权值，w1+w2+w1+w4＝1。
5.根据权利要求4所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，第i只个体进入威胁区域j时的威胁度threat_inj(i)为：Li,k为第i只个体所规划路径第k个分段的长度，tj为威胁区域j的威胁等级，K为第i只个体所规划路径的分段数目，d0.1,k、d0.3,k、d0.5,k、d0.7,k、d0.9,k分别为第i只个体所规划路径第k个分段的十分之一处、十分之三处、十分之五处、十分之七处、十分之九处与威胁区域j中心的距离，Tj为威胁区域j的中心位置，Rj为威胁区域j的半径，Xi(t)为在t次迭代中第i只个***置的三维坐标。
6.根据权利要求4所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，第i只个体在威胁区域j外时的威胁度threat_outj(i)为：|Xi(t)-Tj|＞Rj，tj为威胁区域j的威胁等级，K为第i只个体所规划路径的分段数目，d0.1,k、d0.3,k、d0.5,k、d0.7,k、d0.9,k分别为第i只个体所规划路径第k个分段的十分之一处、十分之三处、十分之五处、十分之七处、十分之九处与威胁区域j中心的距离，Tj为威胁区域j的中心位置，Rj为威胁区域j的半径，Xi(t)为在t次迭代中第i只个***置的三维坐标。
7.根据权利要求5或6所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，威胁区域j的威胁等级tj为： 表示从威胁区域j的中心位置指向第i-1只个***置的向量，表示从第i-1只个***置指向第i只个***置的向量。
8.根据权利要求4所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，第i只个体规划路径的机动角度angle(i)为：M为个体数目，表示从第i-2只个***置指向第i-1只个***置的向量，表示从第i-1只个***置指向第i只个***置的向量。
9.根据权利要求2所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，地标算子迭代演化阶段依据表达式：更新种群数量，NP(t)为第t次迭代的种群数量，N为种群数量，Tmax为最大迭代次数。
10.根据权利要求2所述基于混沌种群变异PIO的航天器姿态轨道协同规划方法，其特征在于，在指南针算子迭代演化阶段开始前，采用Tent Map映射混沌初始化种群。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810366021.3A CN108919818B (zh) | 2018-04-23 | 2018-04-23 | 基于混沌种群变异pio的航天器姿态轨道协同规划方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810366021.3A CN108919818B (zh) | 2018-04-23 | 2018-04-23 | 基于混沌种群变异pio的航天器姿态轨道协同规划方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108919818A true CN108919818A (zh) | 2018-11-30 |
| CN108919818B CN108919818B (zh) | 2020-08-04 |
Family
ID=64403124
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810366021.3A Active CN108919818B (zh) | 2018-04-23 | 2018-04-23 | 基于混沌种群变异pio的航天器姿态轨道协同规划方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108919818B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109813304A (zh) * | 2019-03-25 | 2019-05-28 | 西北工业大学 | 一种基于进化策略的分段搜索地磁仿生导航方法 |
| CN110562491A (zh) * | 2019-08-21 | 2019-12-13 | 中国地质大学（武汉） | 基于种群分布状态对太空发电站做姿态控制的方法及*** |
| CN112596390A (zh) * | 2020-12-18 | 2021-04-02 | 哈尔滨工业大学 | 一种航天器姿控***pwpf调制器参数确定方法 |
| CN113095538A (zh) * | 2020-06-08 | 2021-07-09 | 华北电力大学 | 面向灵活性运行的热电联产机组宽负荷运行动态特性建模方法 |
| CN113619815A (zh) * | 2021-08-04 | 2021-11-09 | 西北工业大学 | 一种航天器集群动态路径规划方法 |
| CN114397818A (zh) * | 2022-01-04 | 2022-04-26 | 南京航空航天大学 | 基于mapio的航天器集群轨道重构路径规划方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101226561A (zh) * | 2007-12-28 | 2008-07-23 | 南京航空航天大学 | 用于微型航天器姿态轨道控制***的微型仿真支持***及工作方法 |
| US7410130B2 (en) * | 2003-12-03 | 2008-08-12 | The Boeing Company | Star-tracker-based attitude determination for spinning spacecraft |
| CN101554926A (zh) * | 2009-05-20 | 2009-10-14 | 上海微小卫星工程中心 | 航天器的姿态控制***及方法 |
| CN106950833A (zh) * | 2017-03-28 | 2017-07-14 | 南京航空航天大学 | 基于改进pio算法的卫星姿态动态规划方法 |
-
2018
- 2018-04-23 CN CN201810366021.3A patent/CN108919818B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7410130B2 (en) * | 2003-12-03 | 2008-08-12 | The Boeing Company | Star-tracker-based attitude determination for spinning spacecraft |
| CN101226561A (zh) * | 2007-12-28 | 2008-07-23 | 南京航空航天大学 | 用于微型航天器姿态轨道控制***的微型仿真支持***及工作方法 |
| CN101554926A (zh) * | 2009-05-20 | 2009-10-14 | 上海微小卫星工程中心 | 航天器的姿态控制***及方法 |
| CN106950833A (zh) * | 2017-03-28 | 2017-07-14 | 南京航空航天大学 | 基于改进pio算法的卫星姿态动态规划方法 |
Non-Patent Citations (2)
| Title |
|---|
| BING HUA ,ETC: "Intelligent attitude planning algorithm based on the characteristics of low radar cross section characteristics of microsatellites under complex constraints", 《PROCEEDINGS OF THE INSTITUTION OF MECHANICAL ENGINEERS》 * |
| 张翀， 王彦恺: "航天器编队姿态与轨道耦合协同控制器设计", 《新型工业化》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109813304A (zh) * | 2019-03-25 | 2019-05-28 | 西北工业大学 | 一种基于进化策略的分段搜索地磁仿生导航方法 |
| CN110562491A (zh) * | 2019-08-21 | 2019-12-13 | 中国地质大学（武汉） | 基于种群分布状态对太空发电站做姿态控制的方法及*** |
| CN113095538A (zh) * | 2020-06-08 | 2021-07-09 | 华北电力大学 | 面向灵活性运行的热电联产机组宽负荷运行动态特性建模方法 |
| CN112596390A (zh) * | 2020-12-18 | 2021-04-02 | 哈尔滨工业大学 | 一种航天器姿控***pwpf调制器参数确定方法 |
| CN112596390B (zh) * | 2020-12-18 | 2021-10-19 | 哈尔滨工业大学 | 一种航天器姿控***pwpf调制器参数确定方法 |
| CN113619815A (zh) * | 2021-08-04 | 2021-11-09 | 西北工业大学 | 一种航天器集群动态路径规划方法 |
| CN113619815B (zh) * | 2021-08-04 | 2022-07-12 | 西北工业大学 | 一种航天器集群动态路径规划方法 |
| CN114397818A (zh) * | 2022-01-04 | 2022-04-26 | 南京航空航天大学 | 基于mapio的航天器集群轨道重构路径规划方法 |
| CN114397818B (zh) * | 2022-01-04 | 2023-11-07 | 南京航空航天大学 | 基于mapio的航天器集群轨道重构路径规划方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108919818B (zh) | 2020-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108919818A (zh) | 基于混沌种群变异pio的航天器姿态轨道协同规划方法 | |
| CN106969778B (zh) | 一种多无人机协同施药的路径规划方法 | |
| Flint et al. | Cooperative control for multiple autonomous UAV's searching for targets | |
| CN112016812B (zh) | 多无人机任务调度方法、***及存储介质 | |
| Edison et al. | Integrated task assignment and path optimization for cooperating uninhabited aerial vehicles using genetic algorithms | |
| Wen et al. | UAV online path planning algorithm in a low altitude dangerous environment | |
| CN109032168A (zh) | 一种基于dqn的多无人机协同区域监视的航路规划方法 | |
| CN108388250B (zh) | 一种基于自适应布谷鸟搜索算法的水面无人艇路径规划方法 | |
| CN111121784B (zh) | 一种无人侦察机航路规划方法 | |
| CN113238232B (zh) | 面向海洋静态目标的自主水下航行器***目标搜索方法 | |
| CN111664852A (zh) | 一种无人机路径规划方法及装置 | |
| CN115202394A (zh) | 一种基于改进遗传算法的无人机全覆盖路径规划方法 | |
| CN114840020A (zh) | 一种基于改进鲸鱼算法的无人机飞行轨迹规划方法 | |
| CN115903879A (zh) | 基于地形数据插值技术的无人机航迹规划方法 | |
| CN115903888A (zh) | 一种基于天牛群算法的旋翼无人机自主路径规划方法 | |
| CN111797966B (zh) | 一种基于改进羊群算法的多机协同全局目标分配方法 | |
| CN116088576A (zh) | 一种基于改进鲸鱼算法的无人机三维路径规划方法 | |
| Polycarpou et al. | Cooperative control of distributed multi-agent systems | |
| CN112733251A (zh) | 一种多无人飞行器协同航迹规划方法 | |
| CN116772880A (zh) | 一种基于无人机视觉的无人车路径规划方法 | |
| CN113220008B (zh) | 多火星飞行器的协同动态路径规划方法 | |
| Carolina Jara Ten Kathen et al. | A comparison of pso-based informative path planners for autonomous surface vehicles for water resource monitoring | |
| CN117806346A (zh) | 一种基于自适应蜣螂优化器的无人机路径规划方法 | |
| Zhu et al. | A cooperative task assignment method of multi-UAV based on self organizing map | |
| CN116501089A (zh) | 基于改进蛇优化算法的无人机三维路径规划方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |