CN102419597B - 一种限定相对姿态的大规模编队航天器姿态一致控制方法 - Google Patents
一种限定相对姿态的大规模编队航天器姿态一致控制方法 Download PDFInfo
- Publication number
- CN102419597B CN102419597B CN 201110397491 CN201110397491A CN102419597B CN 102419597 B CN102419597 B CN 102419597B CN 201110397491 CN201110397491 CN 201110397491 CN 201110397491 A CN201110397491 A CN 201110397491A CN 102419597 B CN102419597 B CN 102419597B
- Authority
- CN
- China
- Prior art keywords
- sigma
- spacecraft
- attitude
- control
- gesture
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种限定相对姿态的大规模编队航天器姿态一致控制方法,涉及大规模编队航天器的姿态协同控制方法,为了解决目前在有限敏感器、星载处理器能力较弱、通信带宽有限的情况下,难以完成大规模编队飞行任务的问题,它包括具体步骤如下:在控制过程中的每个调整姿态周期中,通过下述方法获得每颗航天器的控制力矩实现姿态调整:步骤一、基于近邻原则设计星间通信拓扑结构;决定敏感范围δ,将距第i颗航天器δ范围内所有的其它航天器都视为与所述第i颗航天器发生信息交互、协同;步骤二、实现两两邻近航天器间的相对姿态不超出预期的上限σ;第三步、获得第i颗航天器的控制输入力矩τi,τi作为输入力矩控制其姿态。用于航天器大规模编队飞行。
Description
技术领域
本发明涉及大规模编队航天器的姿态协同控制方法。
背景技术
随着高分辨率对地观测、深空探测等一系列航天活动的蓬勃开展,近年来的各种空间任务对航天器功能和性能的需求在不断提高,其研制和发射成本大幅增加,任务周期越来越长,风险也随之升高。另外,当前及未来的众多空间任务要求载荷系统具备长基线或多点同步工作的能力,如分布式SAR系统。为此,在上个世纪90年代,国外学者提出了“航天器编队飞行”的概念。所谓编队飞行,就是由若干颗微小型航天器在周期相同的轨道上构成某个特定形状(构型),所有航天器一方面保持这个构型,同时绕地球中心运转。整个编队系统构成一个超级“虚拟航天器”,航天任务是由整个编队系统来完成的,编队中的航天器之间存在紧密的信息交互和协同。
由若干颗航天器构成的编队系统已经得到较多的研究和初级试验。但随着给定控制精度要求的低燃耗需求、有限敏感器、传感器和通信资源的系统、高精度高稳定度快速机动等航天技术的不断发展和空间长基线干涉测量、大范围立体成像、间断式定位导航等空间任务需求的推动,针对由数十颗、上百颗航天器构成的大规模编队飞行系统的研究已经提上日程。
由几颗航天器构成的编队系统可以利用全局信息来设计系统方案。而对于大规模航天器编队飞行系统,在获取和处理全局信息时受到极大限制,面临通信链路复杂、通信网络脆弱和通信带宽有限等诸多问题。因此,目前在有限敏感器、星载处理器能力较弱、通信带宽有限的情况下,难以完成大规模编队飞行任务。
发明内容
本发明的目的是为了解决目前在有限敏感器、星载处理器能力较弱、通信带宽有限的情况下,难以完成大规模编队飞行任务的问题,提供一种限定相对姿态的大规模编队航天器姿态一致控制方法。
一种限定相对姿态的大规模编队航天器姿态一致控制方法,它包括具体步骤如下:
选定各航天器姿态运动学及动力学的描述方式:采用修正的罗德里格斯参数来描述各航天器的姿态,得到拉格朗日方程形式的第i颗航天器的姿态动力学模型:
式中:σi即描述第i颗航天器姿态的罗德里格斯参数,τi为第i颗航天器的控制输入力矩;
在控制过程中的每个调整姿态周期中,通过下述方法获得每颗航天器的控制力矩实现姿态调整:
步骤一、基于近邻原则设计星间通信拓扑结构;任务需求决定敏感器载荷的参数,进而决定敏感范围δ,位于第i颗航天器敏感范围δ内所有的其它航天器的集合为Ni,该集合Ni内的每一个航天器都视为与所述第i颗航天器发生信息交互、协同,根据所述集合Ni得到通信拓扑结构图,采用Laplacian矩阵描述所述通信拓扑结构图;
步骤二、实现两两邻近航天器间的相对姿态不超出预期的上限σ,采用人工势场法设计比例控制系数pσi,使pσi随着||σi-σj||的变化而变化以达到自行调节的控制作用,||σi-σj||是第i颗航天器与第j颗航天器的相对姿态,pσi为>0的比例控制系数;
步骤三、将步骤二获得的比例控制系数pσi带入公式:
计算获得第i颗航天器的控制输入力矩τi,将该控制输入力矩τi作为第i颗航天器的输入力矩控制其姿态。
本发明基于局部(相对)信息,摒弃了传统的基于全局信息的集中方式。作用于各航天器的控制输入构建于其自身与邻近航天器间的相对信息之上,不需要中央通信与数据处理机制,有效地减少了数据处理量、节省了系统资源、提高了系统的可靠性。解决了协同自身姿态的同时保持相对姿态受限的问题,当相对姿态在预期范围内时,基于人工势函数的变比例控制系数的调节作用较小,可忽略;当相对姿态临近期望范围的上限时,控制系数的调节作用将得以快速有效地发挥。这种方法对节省系统资源是有益的,经验证调节作用也是显著的。
附图说明
图1为本发明的通信拓扑的空间结构示意图,图2为本发明的通信拓扑在XY坐标的示意图,图3为两颗临近航天器间的人工势能的示意图,图4为由人工势能得到的两颗临近航天器相对姿态的权系数的示意图。
具体实施方式
具体实施方式一:结合图1和图2说明本实施方式,本实施方式包括具体步骤如下:
选定各航天器姿态运动学及动力学的描述方式:采用修正的罗德里格斯参数来描述各航天器的姿态,得到拉格朗日方程形式的第i颗航天器的姿态动力学模型:
式中:σi即描述第i颗航天器姿态的罗德里格斯参数,τi为第i颗航天器的控制输入力矩; I3是3×3的单位矩阵, 为对称正定矩阵, 为反对称矩阵;
在控制过程中的每个调整姿态周期中,通过下述方法获得每颗航天器的控制力矩实现姿态调整:
步骤一、基于近邻原则设计星间通信拓扑结构;任务需求决定敏感器载荷的参数,进而决定敏感范围δ,位于第i颗航天器敏感范围δ内所有的其它航天器的集合为Ni,该集合Ni内的每一个航天器都视为与所述第i颗航天器发生信息交互、协同,根据所述集合Ni得到通信拓扑结构图,采用Laplacian矩阵描述所述通信拓扑结构图;
现有方法不能保证邻近航天器间的相对姿态始终在预期范围σ内;姿态机动过程中,可能会出现两邻近航天器间姿态偏差过大,星上光学相对测量、通信单元不能正常工作,导致互相无法信息交互,通信链路失效、通信拓扑结构遭到破坏;
步骤二、实现两两邻近航天器间的相对姿态不超出预期的上限σ,采用人工势场法设计比例控制系数pσi,使pσi随着||σi-σj||的变化而变化以达到自行调节的控制作用,||σi-σj||是第i颗航天器与第j颗航天器的相对姿态,pσi为>0的比例控制系数;
根据比例控制的思想,当||σi-σj||变大时,需要相应增大其比例控制系数,施加更大的控制输入以抑制||σi-σj||继续增大乃至使其减小。这与物体在势场中的受力规律是一致的。如普遍存在的万有引力场,在有效距离内,两个物体间的作用力随着间距变化而变化,距离增大时引力呈主导作用,距离减小时斥力呈主导作用;这就是人工势场法借鉴的思想;
步骤三、将步骤二获得的比例控制系数pσi带入公式:
计算获得第i颗航天器的控制输入力矩τi,将该控制输入力矩τi作为第i颗航天器的输入力矩控制其姿态。
应用所述方法解决了大规模编队航天器编队飞行时协同自身姿态的同时保持相对姿态的问题。
本发明能限定邻近航天器间的相对姿态,邻近航天器间的相对姿态被限定在预先设定的范围内,有利于更好地完成某些空间飞行任务,这是基于人工势函数的控制比例系数在起调节作用,当邻近航天器间的相对姿态的模临近预设范围边界值时,该系数将急剧增大以快速抑制相对姿态,同时,能有效防止编队系统中的航天器翻转。
针对大规模航天器编队飞行系统,不受编队航天器数量的限制,能应用于若干颗航天器编队飞行。且便利地利用图Laplaican矩阵刻画了各航天器间的信息交互,能基于此分析系统一致收敛速度、可控可观性。少量航天器受到外界干扰时,编队系统中各航天器的姿态仍会在控制输入的作用下逐渐同步,虽然不能完全消除干扰的影响,但能有效地抑制干扰;有航天器出现故障脱离编队系统、有新的航天器被编入系统,都无需对控制算法进行修改,从网络的角度来讲,少量节点不影响整个网络系统。
具体实施方式二:步骤三中所述的权系数w(σi,σj)的获得方法为:
构造人工势函数V(σi,σj)为:
从人工势函数V(σi,σj)导出权系数w(σi,σj),根据 和 获得:
势场是普遍存在的,任何两个物体间都存在着相互作用,因距离不同而表现为引力或斥力。将该思想引入协调控制中,即可保持两个被控状态的相对值。如图3,既是人工势能。能量是有限的,i和j间不可能越过图3中左边曲线的右顶点,因为右顶点不存在,该曲线是无限接近垂直于x轴的渐近线但永远无法接近的。
实际上,是先构造人工势函数V(σi,σj),再导出w(σi,σj), (σi-σj)与一致控制算法相匹配,这才可以将w(σi,σj)作为与状态相关的权系数,如图4。而这正是构造V(σi,σj)的出发点。
从问题本身来说,我们需要在机动过程中保持两两临近航天器间的相对姿态的模不超过某个上限,因为航天器上的固定相对测量敏感器的视场是有限的,如相机;第一,如果不固定,给相机一个附加的随动系统,随动系统不仅增加航天器的挠性及重量,还需要额外的能源;第二,即使让相机随动,相机无法决策在某个时刻测量哪个临近航天器的信息。那么,人工势场法以系统资源为代价,自主协调临近航天器间的姿态。
该变比例控制系数适用的前提是编队系统分布式通信拓扑结构图必须是无向的,进一步否定了集中式和主从式的编队方案。
Claims (1)
1.一种限定相对姿态的大规模编队航天器姿态一致控制方法,其特征是它包括具体步骤如下:
选定各航天器姿态运动学及动力学的描述方式:采用修正的罗德里格斯参数来描述各航天器的姿态,得到拉格朗日方程形式的第i颗航天器的姿态动力学模型:
式中:σi即描述第i颗航天器姿态的罗德里格斯参数,τi为第i颗航天器的控制输入力矩;
在控制过程中的每个调整姿态周期中,通过下述方法获得每颗航天器的控制力矩实现姿态调整:
步骤一、基于近邻原则设计星间通信拓扑结构;任务需求决定敏感器载荷的参数,进而决定敏感范围δ,位于第i颗航天器敏感范围δ内所有的其它航天器的集合为Ni,该集合Ni内的每一个航天器都视为与所述第i颗航天器发生信息交互、协同,根据所述集合Ni得到通信拓扑结构图,采用Laplacian矩阵描述所述通信拓扑结构图;
步骤二、实现两两邻近航天器间的相对姿态不超出预期的上限σ,采用人工势场法设计比例控制系数pσi,使pσi随着‖σi-σj‖的变化而变化以达到自行调节的控制作用,‖σi-σj‖是第i颗航天器与第j颗航天器的相对姿态,pσi为>0的比例控制系数;
步骤三、将步骤二获得的比例控制系数pσi带入公式:
所述的权系数w(σi,σj)的获得方法为:
构造人工势函数V(σi,σj)为:
从人工势函数V(σi,σj)导出权系数w(σi,σj),根据
计算获得第i颗航天器的控制输入力矩τi,将该控制输入力矩τi作为第i颗航天器的输入力矩控制其姿态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110397491 CN102419597B (zh) | 2011-12-05 | 2011-12-05 | 一种限定相对姿态的大规模编队航天器姿态一致控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110397491 CN102419597B (zh) | 2011-12-05 | 2011-12-05 | 一种限定相对姿态的大规模编队航天器姿态一致控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102419597A CN102419597A (zh) | 2012-04-18 |
CN102419597B true CN102419597B (zh) | 2013-03-13 |
Family
ID=45944039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110397491 Expired - Fee Related CN102419597B (zh) | 2011-12-05 | 2011-12-05 | 一种限定相对姿态的大规模编队航天器姿态一致控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102419597B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104748629B (zh) * | 2015-04-02 | 2016-08-17 | 西北工业大学 | 一种基于人工势场法的水声多目标制导方法 |
CN105512391B (zh) * | 2015-12-04 | 2018-09-25 | 上海新跃仪表厂 | 基于对偶四元数的多星姿轨动力学建模方法及其验证系统 |
CN105955028B (zh) * | 2016-06-02 | 2018-09-07 | 西北工业大学 | 一种航天器在轨规避制导控制一体化算法 |
CN107422741B (zh) * | 2017-04-28 | 2019-12-27 | 西北工业大学 | 基于学习的保预设性能集群飞行分布式姿态追踪控制方法 |
CN107145071B (zh) * | 2017-05-19 | 2019-12-27 | 西北工业大学 | 一种空间系绳系统的拉格朗日动力模型及控制器 |
CN108181916B (zh) * | 2017-12-29 | 2020-04-24 | 清华大学 | 小卫星相对姿态的控制方法及装置 |
CN110110342B (zh) * | 2018-12-11 | 2023-08-25 | 上海航天控制技术研究所 | 一种基于邻近算法的组合体航天器数据驱动控制方法 |
CN110006447B (zh) * | 2019-04-04 | 2020-11-10 | 北京临近空间飞行器系统工程研究所 | 一种无需初始对准的任意姿态mems组合定姿方法 |
CN110598249A (zh) * | 2019-08-02 | 2019-12-20 | 南京航空航天大学 | 大规模航天器集群空间均匀分布方法 |
CN110609572B (zh) * | 2019-08-22 | 2022-09-06 | 合肥工业大学 | 无人机编队信息交互拓扑的启发式优化方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101576750A (zh) * | 2009-04-14 | 2009-11-11 | 上海微小卫星工程中心 | 航天器的姿态跟踪控制系统及方法 |
CN101709973A (zh) * | 2009-11-26 | 2010-05-19 | 哈尔滨工业大学 | 一种大规模编队相对导航方法 |
CN101708780A (zh) * | 2009-11-26 | 2010-05-19 | 哈尔滨工业大学 | 用于目标姿态跟踪的刚性航天器的控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060290311A1 (en) * | 2005-06-28 | 2006-12-28 | Nikhil Chopra | Method and System for Synchronizing Networked Passive Systems |
-
2011
- 2011-12-05 CN CN 201110397491 patent/CN102419597B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101576750A (zh) * | 2009-04-14 | 2009-11-11 | 上海微小卫星工程中心 | 航天器的姿态跟踪控制系统及方法 |
CN101709973A (zh) * | 2009-11-26 | 2010-05-19 | 哈尔滨工业大学 | 一种大规模编队相对导航方法 |
CN101708780A (zh) * | 2009-11-26 | 2010-05-19 | 哈尔滨工业大学 | 用于目标姿态跟踪的刚性航天器的控制方法 |
Non-Patent Citations (2)
Title |
---|
航天器编队飞行多目标姿态快速跟踪鲁棒控制;袁长清等;《应用数学和力学》;20080215;第29卷(第02期);169-180 * |
袁长清等.航天器编队飞行多目标姿态快速跟踪鲁棒控制.《应用数学和力学》.2008,第29卷(第02期),169-180. |
Also Published As
Publication number | Publication date |
---|---|
CN102419597A (zh) | 2012-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102419597B (zh) | 一种限定相对姿态的大规模编队航天器姿态一致控制方法 | |
CN105138010B (zh) | 一种编队卫星分布式有限时间跟踪控制方法 | |
CN109324636B (zh) | 基于二阶一致性和自抗扰的多四旋翼主从式协同编队控制方法 | |
CN110376882A (zh) | 基于有限时间扩张状态观测器的预定性能控制方法 | |
CN105159305B (zh) | 一种基于滑模变结构的四旋翼飞行控制方法 | |
CN105700536A (zh) | 基于绳系拖曳系统的主动星姿态和系绳摆振联合控制方法 | |
CN105353615A (zh) | 一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 | |
CN104898148A (zh) | 一种基于数据压缩和神经网络的低成本ins/gps无缝导航方法 | |
CN108241292A (zh) | 一种基于扩张状态观测器的水下机器人滑模控制方法 | |
Lan et al. | Finite-time control for soft landing on an asteroid based on line-of-sight angle | |
CN109683628A (zh) | 一种基于有限时间分布式速度观测器的航天器相对位置控制方法 | |
CN105759827A (zh) | 一种抑制不期望柔性振动的航天器姿态控制系统 | |
CN106707749A (zh) | 一种针对仿生扑翼飞行机器人的控制方法 | |
CN108508751A (zh) | 一种输入饱和自适应姿态协同跟踪控制方法 | |
CN105116914A (zh) | 一种平流层飞艇解析模型预测路径跟踪控制方法 | |
CN113867143A (zh) | 地外天体安全软着陆解析避障制导方法 | |
CN102866633B (zh) | 微陀螺仪的动态滑模控制系统 | |
Wang et al. | Fixed-time event-triggered sliding mode cooperative path-following control with prescribed performance for USVs based on lumped disturbance observer | |
CN111123702B (zh) | 基于一致性理论的大型柔性航天器分散协调鲁棒控制方法 | |
CN103438892A (zh) | 一种改进的基于ekf的天文自主定轨算法 | |
CN115688337A (zh) | 一种geo航天器多视线融合近场感知队形设计方法 | |
Huang et al. | Optimal satellite formation reconfiguration using co-evolutionary particle swarm optimization in deep space | |
CN111949040A (zh) | 高效利用空间无线资源的卫星编队姿态协同跟踪控制方法 | |
Chen et al. | Adaptive prescribed-time control for coupled 6-DOF spacecraft formation flying | |
Zhang et al. | Intelligent Navigation System based on Big Data Traffic System |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130313 Termination date: 20131205 |