CN110488858A - 一种采用Jordan分解进行编队飞行小推力重构的方法 - Google Patents
一种采用Jordan分解进行编队飞行小推力重构的方法 Download PDFInfo
- Publication number
- CN110488858A CN110488858A CN201910595964.8A CN201910595964A CN110488858A CN 110488858 A CN110488858 A CN 110488858A CN 201910595964 A CN201910595964 A CN 201910595964A CN 110488858 A CN110488858 A CN 110488858A
- Authority
- CN
- China
- Prior art keywords
- indicate
- formation
- value
- configuration
- formula
- 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
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000013598 vector Substances 0.000 claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 239000000446 fuel Substances 0.000 claims abstract description 21
- 238000005457 optimization Methods 0.000 claims abstract description 14
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 230000001133 acceleration Effects 0.000 claims description 22
- 230000017105 transposition Effects 0.000 claims description 6
- 238000003780 insertion Methods 0.000 claims description 4
- 230000037431 insertion Effects 0.000 claims description 4
- 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 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 239000011800 void material Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 230000000737 periodic effect Effects 0.000 abstract description 3
- 238000000926 separation method Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 230000003796 beauty Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
- G05D1/104—Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种采用Jordan分解进行编队飞行小推力重构的方法,包括有:步骤一,建立无摄圆轨道下主星与从星的相对运动模型;步骤二,针对运动系数矩阵进行Jordan分解,得到解耦的相对运动关系的通解;步骤三,将变换后运动状态量的初值视为区分编队构型的不变量,建立由初始轨道构型到目标轨道构型状态量的泛函;步骤四,用多项式描述编队重构的优化系数向量,化简由初始轨道构型到目标轨道构型状态量的泛函;步骤五,以控制燃料消耗总量最小为收敛约束,优化求解编队重构的优化系数向量。本发明应用Jordan分解分离出线性不稳定项和周期项,通过对构型不变量和积分初值的不同定义,以达到最小燃料消耗总量为目标、重构路径初终端为约束,优化求解常值系数向量。
Description
技术领域
本发明涉及卫星编队飞行技术领域,更特别地说,是指一种采用Jordan分解与小推力策略进行编队飞行小推力重构的方法。
背景技术
作为多星系统,航天器编队很大程度上拓展了单一航天器完成任务的局限性,不但能够实现编队中单一航天器的任务功能,而且整个编队可以替代单个大航天器实现更为复杂的任务。在一定条件下,即使在编队中丢失一两个航天器,其他航天器也可协作完成丢失航天器的任务,大大提高编队系统的整体可靠性,具有很强的实际应用价值。由于编队中卫星所携带的燃料有限且无法补给,因此进行编队飞行不同重构策略的研究,具有十分重要的意义。
为了成功地实现编队飞行任务,任务之前和执行任务中的航天器编队轨道设计尤为重要,因此提出了编队构型的相关问题。编队构型机动中主要包含三个主要内容:
构型建立,即将相互独立的航天器群组建为任务所需构型;
构型保持,即通过控制方法,使编队内各航天器在受到轨道摄动等因素的影响下仍可保持稳定的相对位置,维持任务所需构型;
构型重构,即针对编队飞行航天器位置的重分配,根据不同任务的需要而对航天器编队的一颗或多颗航天器进行轨道调整,使航天器编队位置或编队中航天器之间的相对位置发生改变,从而实现编队不同构型之间的转换。
与传统的航天器轨道维持相比,高精度编队飞行的轨道控制对于控制精度具有更高的要求,对于控制频率也有一定的限制。因此,连续小推力轨道控制走进航天研究者的视野,基于连续小推力控制的航天器轨道维持与重构成为众多学者的热门研究内容。作为执行深空探测任务中的常用的一项关键技术,小推力技术具有技术发展成熟、控制精度高、直观的特点。由于新型电推进器以及电磁力和静电力的使用面越来越广,连续小推力方法的应用更为广泛。
此外,在卫星编队重构中,通常都将燃料消耗总量、燃料消耗均衡或者二者的某种加权关系作为编队重构优化指标。
发明内容
鉴于特征值分解在描述相对转道运动模式的不足,本发明采用Jordan分解分离出线性不稳定项和周期项,卫星编队由主星与从星两颗卫星组成,并假设主星运行在无摄圆轨道上。本发明重构方法通过对构型不变量和积分初值的不同定义,优化从星从初始轨道构型到目标轨道构型的重构路径,进而转化满足初终端约束的开环控制规划泛函。基于小推力控制的多项式函数近似求解,以达到最小燃料消耗总量为目标、重构路径初终端为约束,优化求解常值系数向量。
本发明的一种采用Jordan分解进行编队飞行小推力重构的方法,其特征在于包括有下列步骤:
步骤一,建立无摄圆轨道下主星与从星的相对运动模型;
为了维持编队飞行构型采用CW方程来描述主星与从星的相对位置关系,并在所述的CW方程中添加了外部控制加速度u,使得主星与从星满足式(1)的相对运动关系:
A表示Φ的右下角对角元素,且ω表示主星的轨道角速度; -A表示A的相反数;
B表示Φ的左下角对角元素,-B表示B的相反数;
矩阵Φ的特征值是反映式(1)动态行为的主要指标,即式中空白项均为0,U表示Φ的特征向量,U-1表示U的逆矩阵;
步骤二,针对运动系数矩阵进行Jordan分解,得到解耦的相对运动关系的通解;
由步骤一得到矩阵Φ,针对Φ进行Jordan分解V-1ΦV=J,V表示Jordan分解的特征向量,V-1表示V的逆矩阵,得到表征六维运动模态的特征矩阵J,其中
令Z=V-1X,Z表示变换后运动状态量,X表示从星S2相对于主星S1的运动状态量,并将J带入式(1)可得相对运动关系:
当u=0时,由简化的式(2)可得通解式为Z(t)=eJtZ0,Z(t)表示随时间变化的属于X的坐标变换后的状态量,eJt表示状态转移矩阵,J表示六维运动模态的特征矩阵,t表示时间,Z0表示Z的状态初始值;
在u=0时,Z0为常值向量,而状态转移矩阵写作:
所述eJt中的元素分别表征不同的运动模态,其中,cosωt是主星角速度与时间乘积的余弦值,sinωt是主星角速度与时间乘积的正弦值,-sinωt是sinωt的相反数;
步骤三,将变换后运动状态量的初值Z0视为区分编队构型的不变量,建立由初始轨道构型到目标轨道构型状态量的泛函;
在通解Z(t)=eJtZ0中,Z0除了作为积分常数外,还可视为区别编队构型的不变量;状态转移矩阵eJt可以看作为CW方程的基础解系,即从星相对主星的所有运动都视为状态转移矩阵所有列向量的线性组合;
当u≠0时,Z(t)仍可视为状态转移矩阵所有列向量的时变线性组合,即:
Z(t)=eJtZ0(t) (4)
将式(4)求导并代入式(3),可以得到:
联立式(4)和式(5),可得u(t)的最小二乘解为:
将记为M,为同时满足变换后运动状态量的初值微分条件和加速度加载为最小二乘解形式,将P(PTP)-1PT记为N,则:
(N-I)·M=0 (7)
I表示单位矩阵;此时,M表征为N的特征值为1时的特征向量;N具有特征值分布为{0,0,0,1,1,1},其中特征值为1时N的正交特征向量记为E;因此,满足式(7)的随时间变化的任意M可表示为:
M(t)=E·γ(t) (8)
将(8)式带入可得则编队重构问题可表示为由所述γ(t)到编队重构的目标轨道构型状态量Z0(Td)的泛函:
步骤四,用多项式描述编队重构的优化系数向量γ(t),化简由初始轨道构型到目标轨道构型状态量的泛函;
采用多项式近似求解式(9)所建立的编队由初始构型到目标构型的重构路径;则编队重构的优化向量γ(t)可以写作:
将式(10)代入式(9)可以得到:
为了简化,将记为重构路径常数矩阵Gj,可知Gj与重构路径无关,可以离线计算;当j最大值取2时,式(11)可改写为:
Z0(Td)-Z0(0)=[G0E G1E G2E]·[D0 D1 D2]T (12)
步骤五,以控制燃料消耗总量L最小为收敛约束,优化求解编队重构的优化系数向量γ(t);
对于同一重构路径,Dj的选择不唯一,故需根据控制燃料消耗总量选择;将和代入式(6)可得:
则控制燃料消耗总量可以表示为:
对于Dj的选择是将编队重构的路径求解问题转化为具有等式约束的非线性燃料消耗总量优化问题,即求解燃料消耗总量的最小Lmin且Dj满足 Z0(Td)-Z0(0)=[G0E G1EG2E]·[D0 D1 D2]T。
本发明采用Jordan分解进行编队飞行小推力重构的优点在于:
①利用Jordan分解处理相对运动方程,解决传统特征值分解方法处理过程中隐藏长期漂移模态的问题,有利于揭示相对运动的所有运动模态,从而更方便的对相对运动进行分析处理。
②本发明方法采用求解Jordan分解处理后的相对运动学方程可得基础解系,提取提取基础解系的组合系数作为不变量,推导全新的重构模型,从而将编队重构路径问题转换为泛函问题,可直接优化求解,减少计算量。
③采用控制消耗作为优化变量,将重构问题中自变量的选择转化为具有等式约束的非线性优化问题,求解出来的重构路径可以大大减少能量消耗,延长卫星编队整体寿命。
附图说明
图1是两星之间的相对位置示意图。
图2是本发明实施例1中重构前后轨迹以及重构路径图。
图2A是图2的俯视图。
图2B是图2的右视图。
图2C是图2的前视图。
图3是本发明实施例1中重构期间的小推力加速度时间历程图。
具体实施方式
下面将结合附图和实施例对本发明做进一步的详细说明。
当用运动学方程描述编队飞行卫星的相对运动关系时,常采用特征值分解的方法分离出不同的运动模态。但特征值方法并不适用所有的情况,因此在本发明中,采用Jordan分解对相对运动关系进行处理。Jordan分解相关内容参考2013年3月第 2版的《高等代数》,作者李方等,第94-98页。
本发明的一种采用Jordan分解进行编队飞行小推力重构的方法,包括有下列步骤:
步骤一,建立无摄圆轨道下主星与从星的相对运动模型;
此卫星编队由主星与从星两颗卫星组成,并假设主星运行在无摄圆轨道上。参见图1所示的主星轨道坐标系定义如下:主星记为S1,从星记为S2,原点为主星质心, x轴背离地心,z轴垂直于卫星轨道平面,y轴在卫星运动轨道平面内。地心O到主星质心的距离矢量记为r1,地心O到从星质心的距离矢量记为r2,主星质心与从星质心的距离记为ρ。为了维持编队飞行构型采用CW方程来描述主星与从星的相对位置关系,并在所述的CW方程中添加了外部控制加速度u,使得主星与从星满足式 (1)的相对运动关系:
X表示从星S2相对于主星S1的运动状态量。
表示X的变化率。
Φ表示属于X的系数矩阵,即无摄圆轨道下主星与从星的相对运动模型,简称为运动系数矩阵。
u表示小推力发动机提供的对从星施加的控制加速度(简称加速度加载)。
A表示Φ的右下角对角元素,且ω表示主星的轨道角速度; -A表示A的相反数。
B表示Φ的左下角对角元素,-B表示B的相反数。
在本发明中,CW方程参考2015年1月第1版的《航天器编队飞行:动力学、控制、与导航》,作者(美)Kyle T.Alfriend等,张世杰译,第57-58页。
在本发明中,矩阵Φ的特征值是反映式(1)动态行为的主要指标,即式中空白项均为0,U表示Φ的特征向量,U-1表示U的逆矩阵。由Φ的特征值分布可知式(1)只表征了4个维度的周期运动,其中“0”值表明系统临界稳定,没有不稳定项存在。而实际结果表明,相对运动中存在长期漂移项,因此特征值分解不能揭示相对运动全部的运动模态。
在本发明中,也可以存在区间的划分
步骤二,针对运动系数矩阵进行Jordan分解,得到解耦的相对运动关系的通解;
由步骤一得到矩阵Φ,针对Φ进行Jordan分解V-1ΦV=J,V表示Jordan分解的特征向量,V-1表示V的逆矩阵,得到表征六维运动模态的特征矩阵J,其中在本发明中,也可以存在区间的划分
令Z=V-1X,Z表示属于X的坐标变换后的运动状态量(简称为变换后运动状态量),X表示从星S2相对于主星S1的运动状态量,并将J带入式(1)可得相对运动关系:
表示Z的变化率。
P表示u的系数,且P的取值为V的第3列到第6列部分的元素。
在本发明中,当小推力发动机不提供轨道控制加速度(即不考虑u,u=0)时,由简化的式(2)可得通解式为Z(t)=eJtZ0,Z(t)表示随时间变化的属于X的坐标变换后的状态量,eJt表示状态转移矩阵,J表示六维运动模态的特征矩阵,t表示时间, Z0表示Z的状态初始值,在u=0时Z0为常值向量,而状态转移矩阵写作:
所述eJt中的元素分别表征不同的运动模态,其中,cosωt是主星角速度与时间乘积的余弦值,sinωt是主星角速度与时间乘积的正弦值,-sinωt是sinωt的相反数。
步骤三,将变换后运动状态量的初值Z0视为区分编队构型的不变量,建立由初始轨道构型到目标轨道构型状态量的泛函;
在通解Z(t)=eJtZ0中,Z0除了作为积分常数外,还可视为区别编队构型的不变量;状态转移矩阵eJt可以看作为CW方程的基础解系,即从星相对主星的所有运动都视为状态转移矩阵所有列向量的线性组合。当小推力发动机向从星提供轨道控制加速度(即u≠0)时,Z(t)仍可视为状态转移矩阵所有列向量的时变线性组合,即:
Z(t)=eJtZ0(t) (4)
将式(4)求导并代入式(3),可以得到:
表示Z0的变化率,其物理意义是在小推力发动机控制u作用下,卫星编队中从星从初始轨道构型重构到目标轨道构型所走的路径。
联立式(4)和式(5),可得u(t)的最小二乘解为:
u表示小推力发动机提供的对从星施加的控制加速度(简称加速度加载)。
u(t)表示随时间变化的加速度加载。
为了简化说明,将记为M(即),式(6)简化为:
u(t)=(PTP)-1PT·M
为同时满足变换后运动状态量的初值微分条件和加速度加载为最小二乘解形式,将P(PTP)-1PT简化记为N,则:
(N-I)·M=0 (7)
I表示单位矩阵。此时,M表征为N的特征值为1时的特征向量。N具有特征值分布为{0,0,0,1,1,1},其中特征值为1时N的正交特征向量记为E。因此,满足式(7)的随时间变化的任意M可表示为:
M(t)=E·γ(t) (8)
γ(t)表示编队重构的优化系数向量,且为任一3×1维向量。
将(8)式带入可得则编队重构问题可表示为由所述γ(t)到编队重构的目标轨道构型状态量Z0(Td)的泛函:
Td表示编队重构结束时间;
e-Jt为状态转移矩阵eJt的逆;
E为特征值为1时N的正交特征向量;
dt为时间的微分;
Z0(Td)表示编队重构的目标构型状态量;
Z0(0)表示编队重构的初始构型状态量;
步骤四,用多项式描述编队重构的优化系数向量γ(t),化简由初始轨道构型到目标轨道构型状态量的泛函;
本发明是采用多项式近似求解式(9)所建立的编队由初始构型到目标构型的重构路径。则编队重构的优化向量γ(t)可以写作:
j表示编队重构的多项式求解阶数,在本发明中取j最大值取2;
tj表示时间的阶数;
Dj表示当阶数为j时的常值系数向量。
将式(10)代入式(9)可以得到:
为了简化,将记为重构路径常数矩阵Gj,可知Gj与重构路径无关,可以离线计算。当j最大值取2时,式(11)可改写为:
Z0(Td)-Z0(0)=[G0E G1E G2E]·[D0 D1 D2]T (12)
Z0(Td)表示编队重构的目标构型状态量;
Z0(0)表示编队重构的初始构型状态量;
E为特征值为1时N的正交特征向量;
Gj表示重构路径常数矩阵;
G0表示j=0时的重构路径常数矩阵取值;
G1表示j=1时的重构路径常数矩阵取值;
G2表示j=2时的重构路径常数矩阵取值;
D0表示j=0时的常值系数向量取值;
D1表示j=1时的常值系数向量取值;
D2表示j=2时的常值系数向量取值;
上角标T为坐标转置符。
步骤五,以控制燃料消耗总量L最小为收敛约束,优化求解编队重构的优化系数向量γ(t);
对于同一重构路径,Dj的选择不唯一,故需根据控制燃料消耗总量选择。将和代入式(6)可得:
则控制燃料消耗总量可以表示为:
k表示转置后加速度加载的多项式求解阶数,在本发明中最大值取2;
j表示编队重构的多项式求解阶数,在本发明中最大值取2;
表示转置后当阶数为k时的常值系数向量。
因此,在本发明中,对于Dj的选择是将编队重构的路径求解问题转化为具有等式约束的非线性燃料消耗总量优化问题,即求解燃料消耗总量的最小Lmin且Dj满足 Z0(Td)-Z0(0)=[G0E G1E G2E]·[D0 D1 D2]T。在本发明中对于低阶优化问题,通过 Matlab软件“fmincon”和“confuneq”即求解。
实施例1
编队任务中的主星运行于800km圆轨道;重构前后的构型参数选取如下: Z0(0)=[1,0,1,1,0,0]T km(构型中心位于沿迹向7m处、构型半径为3.35km)和 Z0(Td)=[-1,0,0,2,2,0]T km(构型中心位于沿迹向-7m处、构型半径为6.08km);重构周期为1个轨道周期Td=6052.4s。
采用本发明的方法实现由构型Z0(0)到Z0(Td)的重构,优化得到 D0=[-0.0066,0.11,-0.64]×10-4、D1=[-0.1425,0.0003,0.000017]×10-3、 D2=[0.37,0,0.000012]×10-7。则重构前后轨迹以及重构路径如图2、图2A、图2B、图2C所示,重构期间的小推力加速度时间历程如图3所示。累计所需速度增量为ΔV=4.43×10-3m/s,最大控制加速度为1.9×10-3m/s2,且出现在重构完成时刻。分析基于多项式函数拟合的重构路径优化方法的数值仿真结果,可以看出本发明方法可以实现由初始相对轨道到目标相对轨道的重构路径规划,且重构路径的位置坐标与初始轨道目标轨道没有大尺度偏移。这说明在重构过程中从星与主星没有产生较大的相对距离,仍可维持在编队范围内,CW方程仍成立。可以看出,小推力发动机提供的加速度加载是通过多项式进行拟合的,其曲线较为平滑。
本发明采用Jordan分解进行编队飞行小推力重构的方法,通过对构型不变量和积分初值的不同定义,优化从星从初始轨道构型到目标轨道构型的重构路径,进而转化满足初终端约束的开环控制规划泛函。基于小推力控制的多项式函数近似求解,以达到最小燃料消耗总量为目标、重构路径初终端为约束,优化求解常值系数向量。
Claims (3)
1.一种采用Jordan分解进行编队飞行小推力重构的方法,其特征在于包括有下列步骤:
步骤一,建立无摄圆轨道下主星与从星的相对运动模型;
为了维持编队飞行构型采用CW方程来描述主星与从星的相对位置关系,并在所述的CW方程中添加了外部控制加速度u,使得主星与从星满足式(1)的相对运动关系:
X表示从星S2相对于主星S1的运动状态量;
表示X的变化率;
Φ表示运动系数矩阵;
u表示加速度加载;
A表示Φ的右下角对角元素,且ω表示主星的轨道角速度;-A表示A的相反数;
B表示Φ的左下角对角元素,-B表示B的相反数;
矩阵Φ的特征值是反映式(1)动态行为的主要指标,即式中空白项均为0,U表示Φ的特征向量,U-1表示U的逆矩阵;
步骤二,针对运动系数矩阵进行Jordan分解,得到解耦的相对运动关系的通解;
由步骤一得到矩阵Φ,针对Φ进行Jordan分解V-1ΦV=J,V表示Jordan分解的特征向量,V-1表示V的逆矩阵,得到表征六维运动模态的特征矩阵J,其中
令Z=V-1X,Z表示变换后运动状态量,X表示从星S2相对于主星S1的运动状态量,并将J带入式(1)可得相对运动关系:
表示Z的变化率;
P表示u的系数,且P的取值为V的第3列到第6列部分的元素;
当u=0时,由简化的式(2)可得通解式为Z(t)=eJtZ0,Z(t)表示随时间变化的属于X的坐标变换后的状态量,eJt表示状态转移矩阵,J表示六维运动模态的特征矩阵,t表示时间,Z0表示Z的状态初始值;
在u=0时,Z0为常值向量,而状态转移矩阵写作:
所述eJt中的元素分别表征不同的运动模态,其中,cosωt是主星角速度与时间乘积的余弦值,sinωt是主星角速度与时间乘积的正弦值,-sinωt是sinωt的相反数;
步骤三,将变换后运动状态量的初值Z0视为区分编队构型的不变量,建立由初始轨道构型到目标轨道构型状态量的泛函;
在通解Z(t)=eJtZ0中,Z0除了作为积分常数外,还可视为区别编队构型的不变量;状态转移矩阵eJt可以看作为CW方程的基础解系,即从星相对主星的所有运动都视为状态转移矩阵所有列向量的线性组合;
当u≠0时,Z(t)仍可视为状态转移矩阵所有列向量的时变线性组合,即:
Z(t)=eJtZ0(t) (4)
将式(4)求导并代入式(3),可以得到:
表示Z0的变化率,其物理意义是在小推力发动机控制u作用下,卫星编队中从星从初始轨道构型重构到目标轨道构型所走的路径;
联立式(4)和式(5),可得u(t)的最小二乘解为:
u表示加速度加载;
u(t)表示随时间变化的加速度加载;
将记为M,为同时满足变换后运动状态量的初值微分条件和加速度加载为最小二乘解形式,将P(PTP)-1PT记为N,则:
(N-I)·M=0 (7)
I表示单位矩阵;此时,M表征为N的特征值为1时的特征向量;N具有特征值分布为{0,0,0,1,1,1},其中特征值为1时N的正交特征向量记为E;因此,满足式(7)的随时间变化的任意M可表示为:
M(t)=E·γ(t) (8)
γ(t)表示编队重构的优化系数向量,且为任一3×1维向量;
将(8)式带入可得则编队重构问题可表示为由所述γ(t)到编队重构的目标轨道构型状态量Z0(Td)的泛函:
Td表示编队重构结束时间;
e-Jt为状态转移矩阵eJt的逆;
E为特征值为1时N的正交特征向量;
dt为时间的微分;
Z0(Td)表示编队重构的目标构型状态量;
Z0(0)表示编队重构的初始构型状态量;
步骤四,用多项式描述编队重构的优化系数向量γ(t),化简由初始轨道构型到目标轨道构型状态量的泛函;
采用多项式近似求解式(9)所建立的编队由初始构型到目标构型的重构路径;则编队重构的优化向量γ(t)可以写作:
j表示编队重构的多项式求解阶数,在本发明中取j最大值取2;
tj表示时间的阶数;
Dj表示当阶数为j时的常值系数向量;
将式(10)代入式(9)可以得到:
为了简化,将记为重构路径常数矩阵Gj,可知Gj与重构路径无关,可以离线计算;当j最大值取2时,式(11)可改写为:
Z0(Td)-Z0(0)=[G0E G1E G2E]·[D0 D1 D2]T (12)
Z0(Td)表示编队重构的目标构型状态量;
Z0(0)表示编队重构的初始构型状态量;
E为特征值为1时N的正交特征向量;
Gj表示重构路径常数矩阵;
G0表示j=0时的重构路径常数矩阵取值;
G1表示j=1时的重构路径常数矩阵取值;
G2表示j=2时的重构路径常数矩阵取值;
D0表示j=0时的常值系数向量取值;
D1表示j=1时的常值系数向量取值;
D2表示j=2时的常值系数向量取值;
上角标T为坐标转置符;
步骤五,以控制燃料消耗总量L最小为收敛约束,优化求解编队重构的优化系数向量γ(t);
对于同一重构路径,Dj的选择不唯一,故需根据控制燃料消耗总量选择;将和代入式(6)可得:
则,控制燃料消耗总量可以表示为:
k表示转置后加速度加载的多项式求解阶数,在本发明中最大值取2;
j表示编队重构的多项式求解阶数,在本发明中最大值取2;
表示转置后当阶数为k时的常值系数向量;
对于Dj的选择是将编队重构的路径求解问题转化为具有等式约束的非线性燃料消耗总量优化问题,即求解燃料消耗总量的最小Lmin且Dj满足Z0(Td)-Z0(0)=[G0E G1E G2E]·[D0D1 D2]T。
2.根据权利要求1所述的采用Jordan分解进行编队飞行小推力重构的方法,其特征在于:卫星编队由主星与从星两颗卫星组成,并假设主星运行在无摄圆轨道上。
3.根据权利要求1所述的采用Jordan分解进行编队飞行小推力重构的方法,其特征在于:累计所需速度增量为ΔV=4.43×10-3m/s,最大控制加速度为1.9×10-3m/s2,且需要较小的最大加速度。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811129709 | 2018-09-27 | ||
CN2018111297096 | 2018-09-27 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110488858A true CN110488858A (zh) | 2019-11-22 |
CN110488858B CN110488858B (zh) | 2020-05-19 |
Family
ID=68546533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910595964.8A Active CN110488858B (zh) | 2018-09-27 | 2019-07-03 | 一种采用Jordan分解进行编队飞行小推力重构的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110488858B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111189457A (zh) * | 2020-02-25 | 2020-05-22 | 上海航天控制技术研究所 | 一种基于cw方程的解耦变增益自主相对导航方法 |
CN113900442A (zh) * | 2021-10-25 | 2022-01-07 | 北京航空航天大学 | 航天器编队展开重构最优控制求解方法和系统 |
CN115535304A (zh) * | 2022-10-09 | 2022-12-30 | 哈尔滨工业大学 | 一种针对多个编队卫星周期性重访的轨道设计与控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060167596A1 (en) * | 2005-01-24 | 2006-07-27 | Bodin William K | Depicting the flight of a formation of UAVs |
CN103676955A (zh) * | 2013-12-19 | 2014-03-26 | 北京航空航天大学 | 一种实现分布式编队飞行的卫星自主轨道控制系统 |
CN104076819A (zh) * | 2014-07-08 | 2014-10-01 | 清华大学 | 一种圆参考轨道下卫星有界伴飞的边界控制方法 |
CN104317306A (zh) * | 2014-10-28 | 2015-01-28 | 蔡远文 | 卫星编队重构方法 |
CN104898691A (zh) * | 2015-04-29 | 2015-09-09 | 哈尔滨工业大学 | 编队卫星有限时间构型包含控制方法 |
CN105912020A (zh) * | 2016-05-26 | 2016-08-31 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN106125759A (zh) * | 2016-07-21 | 2016-11-16 | 北京理工大学 | 一种地球静止轨道上的绳系‑库仑力混合卫星编队方法 |
-
2019
- 2019-07-03 CN CN201910595964.8A patent/CN110488858B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060167596A1 (en) * | 2005-01-24 | 2006-07-27 | Bodin William K | Depicting the flight of a formation of UAVs |
CN103676955A (zh) * | 2013-12-19 | 2014-03-26 | 北京航空航天大学 | 一种实现分布式编队飞行的卫星自主轨道控制系统 |
CN104076819A (zh) * | 2014-07-08 | 2014-10-01 | 清华大学 | 一种圆参考轨道下卫星有界伴飞的边界控制方法 |
CN104317306A (zh) * | 2014-10-28 | 2015-01-28 | 蔡远文 | 卫星编队重构方法 |
CN104898691A (zh) * | 2015-04-29 | 2015-09-09 | 哈尔滨工业大学 | 编队卫星有限时间构型包含控制方法 |
CN105912020A (zh) * | 2016-05-26 | 2016-08-31 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN106125759A (zh) * | 2016-07-21 | 2016-11-16 | 北京理工大学 | 一种地球静止轨道上的绳系‑库仑力混合卫星编队方法 |
Non-Patent Citations (6)
Title |
---|
YANCHAO HE 等: "Distance-based relative orbital elements determination for formation flying system", 《ACTA ASTRONAUTICA》 * |
张博 等: "一种基于信息一致性的卫星编队协同控制策略", 《航空学报》 * |
徐明 等: "地-月小推力转移轨道研究", 《第三届学术会议论文集》 * |
李明涛 等: "编队构型燃料最优可重构部署轨迹规划方法研究", 《空天资源的可持续发展——第一届中国空天安全会议论文集》 * |
杨雪榕 等: "跟飞编队构型的Lyapunov控制研究", 《第三十届中国控制会议》 * |
黄河 等: "考虑燃料均衡的卫星编队队形重构技术", 《宇航学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111189457A (zh) * | 2020-02-25 | 2020-05-22 | 上海航天控制技术研究所 | 一种基于cw方程的解耦变增益自主相对导航方法 |
CN111189457B (zh) * | 2020-02-25 | 2021-12-07 | 上海航天控制技术研究所 | 一种基于cw方程的解耦变增益自主相对导航方法 |
CN113900442A (zh) * | 2021-10-25 | 2022-01-07 | 北京航空航天大学 | 航天器编队展开重构最优控制求解方法和系统 |
CN115535304A (zh) * | 2022-10-09 | 2022-12-30 | 哈尔滨工业大学 | 一种针对多个编队卫星周期性重访的轨道设计与控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110488858B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Junkins et al. | Optimal spacecraft rotational maneuvers | |
Pontani et al. | Particle swarm optimization applied to space trajectories | |
CN110488858A (zh) | 一种采用Jordan分解进行编队飞行小推力重构的方法 | |
Crassidis et al. | New algorithm for attitude determination using Global Positioning System signals | |
Zhang et al. | Optimal two-impulse rendezvous using constrained multiple-revolution Lambert solutions | |
CN110244767A (zh) | 一种采用有限元方法的编队构型重构优化 | |
CN112051854B (zh) | 一种考虑复杂约束的月球软着陆最优轨迹快速规划方法 | |
Azor et al. | Satellite angular rate estimation from vector measurements | |
Christian et al. | Sequential optimal attitude recursion filter | |
Davenport | Rotations about nonorthogonal axes. | |
Xia et al. | Neuroadaptive saturated control for relative motion based noncooperative spacecraft proximity with prescribed performance | |
Thakur et al. | Spacecraft swarm finite-thrust cooperative control for common orbit convergence | |
CN114253288A (zh) | 多航天器轨道分布式协同跟踪最优控制方法 | |
Wang et al. | Attitude-orbit cooperative control for small-scale spacecraft with chemical propulsion: Adaptive sliding mode control based on neural network | |
CN108628345A (zh) | 一种电磁航天器编队悬停协同控制方法及系统 | |
Wang et al. | Stabilization of coupled orbit–attitude dynamics about an asteroid utilizing Hamiltonian structure | |
CRESHAW et al. | Gravity effects on the rotational motion of a uniaxial artificial satellite. | |
Ren et al. | Optimal low-thrust transfers between libration point orbits | |
Xia et al. | Finite-horizon optimal linear control for autonomous soft landing of small-scale helicopter | |
CN108303874A (zh) | 一种针对绳系空间拖船系统摆振的小推力切换控制方法 | |
Meng et al. | Maneuver guidance and formation maintenance for control of leaderless space-robot teams | |
Breger et al. | J2-modified GVE-based MPC for formation flying spacecraft | |
Wang et al. | Optimal four-impulse rendezvous between coplanar elliptical orbits | |
Koru et al. | Adaptive consensus control of rigid body attitude motion based on rotation matrices | |
CN112009727B (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 |