CN115390587B - 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 - Google Patents
一种非预设参数的欠迹向驱动多航天器编队重构控制方法 Download PDFInfo
- Publication number
- CN115390587B CN115390587B CN202211012265.4A CN202211012265A CN115390587B CN 115390587 B CN115390587 B CN 115390587B CN 202211012265 A CN202211012265 A CN 202211012265A CN 115390587 B CN115390587 B CN 115390587B
- Authority
- CN
- China
- Prior art keywords
- sliding mode
- error
- linear transformation
- relative
- parameters
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 19
- 230000009466 transformation Effects 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000014509 gene expression Effects 0.000 claims abstract description 19
- 238000013461 design Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012546 transfer Methods 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
本发明公开了一种非预设参数的欠迹向驱动多航天器编队重构控制方法,所方法包括将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的非预设参数线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的非预设参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
Description
技术领域
本发明涉及欠驱动航天器控制技术领域,具体涉及一种非预设参数的欠迹向驱动多航天器编队重构控制方法。
背景技术
多颗航天器被发射到太空中,首先需要进行构型初始化,接着进行队形的保持,然后根据任务的调整再进行队形重新构造,这些步骤涉及航天器在运行过程中的变轨技术,整体形态发生变化(数量、形状等)。无论是编队飞行、星座构建、集群任务、在轨服务、甚至航天器群博弈,都要求每颗航天器从当前位置转移到指定目标位置,因此重构是多航天器任务执行层面。航天器在运行过程中,由于太空垃圾、宇宙粒子、机械故障等原因可能导致个别推机器失效,此时全驱动控制算法无法应对欠驱动情况,因此研究欠驱动航天器控制算法具有实际意义。
欠驱动情形可分为两种:欠径向驱动和欠迹向驱动,线性化的动力学方程在没有沿轨迹输入的情况下存在一个不可控的特征值,根据当前欠迹向驱动的研究结论指出欠迹向情形不可控状态变量等于圆轨道编队自然构型的保持条件且恒为零,所以不影响系统的能空性。欠迹向驱动系统仅包含5维变量,控制输入只有径向和法向2个自由度,因此设计稳定可控的欠迹向驱动控制方法是一大难点。
进一步的说,当前的面向输入饱和的迹向欠驱动航天器编队重构的控制方法忽视了一个问题,即为了简化控制器的推导,预先设定了欠驱动线性变换矩阵中的某些参数的值。这种操作不仅缺乏理论依据,还可能会限制系统状态之间的耦合关系,预先决定系统的稳定重构精度,从而影响所提出的控制器的普适性。本发明提供的是一种非预设参数的控制器,系统状态之间的耦合关系不像在以前的方法中那样受预设参数的限制,系统收敛性的推导过程及相关表达式更为简洁,同时可以预估制控制器在三个通道位置和速度的控制精度,对实际工程更具借鉴意义。
发明内容
针对现有技术的不足,本发明旨在提供一种非预设参数的欠迹向驱动多航天器编队重构控制方法,通过本发明的方法可以在仅包含径向和法向两个自由度推力时完成卫星编队重构任务。
为了实现上述目的,本发明采用如下技术方案:
一种非预设参数的欠迹向驱动多航天器编队重构控制方法,所方法包括:
S1将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的非预设参数线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;
S2利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;
S3分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
需要说明的是,所述步骤S1包括:
S1.1假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:
其中,eiρ=[ex ey ez]T为x,y,z三个方向的位置误差,为相应的速度误差;将欠迹向驱动动力学方程分解为一个3维和一个2维变量构成的状态表达式:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T;
S1.2由于而ei2a∈R2,为了使状态变量的自由度与控制输入的自由度相匹配,所以设计一个2×3阶非预设参数线性变换矩阵:
其中,a3,b3,f3均为常数;
使用P21对进行线性变换得到一个2阶状态变量:
需要说明的是,利用经线性化变换后的变量设计滑模面,其中:
S2.1根据滑模面设计规则选择如下滑模面:
其中,c2>0为常数。
S2.2等效控制律ui1可通过对滑模面求导得到参数形式的等效控制律:
其中,
Θ1,Θ2分别表示Θ的第一项和第二项;
选择开关控制为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2;
需要说明的是,所述步骤S3包括:
S3.1当k2-δim||si||-γ>0或k1-δim||si||-1>0成立时,闭环系统将在有限时间收敛到期望轨迹,滑模面的收敛误差精度可以被描述为:
而线性化后的状态变量和/>的收敛误差表达式为:
S3.2当线性变换矩阵P21中的非预设参数满足
a21=b3/a3>0;
可以的得到相对位置跟踪误差在x,y,z三个通道的收敛误差为:
相对位置和相对运动速度跟踪误差在迹向的收敛误差为:
相对位置和相对运动速度跟踪误差在法向的收敛误差为:
设置滑模控制参数及满足线性变换矩阵P21中的参数值,根据三个通道的相对位置和相对速度收敛误差的表达式,预估控制器(11)式在x,y,z三个通道的相对位置和相对速度的控制精度。
与现有技术相比,本发明提供的是一种非预设参数的控制器,系统状态之间的耦合关系不像以前的方法中那样受预设参数的限制,系统收敛性的推导过程及相关表达式更为简洁,同时可以预估制控制器在三个通道位置和速度的控制精度,对实际工程更具借鉴意义。
附图说明
图1为本发明的基于线性变换的欠迹向驱动多航天器编队重构控制方法的流程图;
图2为本发明仿真试验中跟随者航天器的控制输入示意图,其中,图2a为跟随者航天器1的控制输入、图2b为跟随者航天器2的控制输入,图2c跟随者航天器3的控制输入;
图3为本发明仿真试验中跟随者航天器的相对于领航者航天器的速度示意图,其中,图3a为跟随者航天器1的相对于领航者航天器的速度、图3b为跟随者航天器1的相对于领航者航天器的速度、图3c为跟随者航天器1的相对于领航者航天器的速度;
图4为本发明仿真试验中3颗跟随者航天器在xy平面、xz平面、yz平面和三维空间的轨迹示意图,其中,图4a为3颗跟随者航天器在xy平面的轨迹、图4b为3颗跟随者航天器在xz平面的轨迹、图4c为3颗跟随者航天器在yz平面的轨迹、图4d为3颗跟随者航天器在三维空间轨迹。
具体实施方式
以下将对本发明作进一步的描述,需要说明的是,以下实施例以本技术方案为前提,给出了详细的实施方式和具体的操作过程,但本发明的保护范围并不限于本实施例。
如图1所示,本发明为一种非预设参数的欠迹向驱动多航天器编队重构控制方法,所方法包括:
S1将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;
S2利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;
S3分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
进一步的,本发明的所述步骤S1包括:
S1.1假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:
其中,为为x,y,z三个方向的位置误差,/>为相应的速度误差;将欠迹向驱动动力学方程分解为一个3维和一个2维变量构成的状态表达式:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T;
S1.2由于而ei2a∈R2,所以设计一个2×3阶非预设参数线性变换矩阵:
其中,a3,b3,f3均为常数;
使用P21对进行线性变换得到一个2阶状态变量:
进一步的,本发明利用经线性化变换后的变量设计滑模面,其中:
S2.1根据滑模面设计规则选择如下滑模面:
其中,c2>0为常数。
S2.2等效控制律ui1可通过对滑模面求导得到参数形式的等效控制律:
其中,
Θ1,Θ2分别表示Θ的第一项和第二项;
选择开关控制为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2;
进一步的,本发明的所述步骤S3包括:
S3.1当k2-δim||si||-γ>0或k1-δim||si||-1>0成立时,闭环系统将在有限时间收敛到期望轨迹,滑模面的收敛误差精度可以被描述为:
而线性化后的状态变量和/>的收敛误差表达式为:
S3.2当线性变换矩阵P21中的参数满足
a21=b3/a3>0;
可以的得到相对位置跟踪误差在x,y,z三个通道的收敛误差为:
相对位置和相对运动速度跟踪误差在迹向的收敛误差为:
|ey|≤Δey=ζin /a 21;
相对位置和相对运动速度跟踪误差在法向的收敛误差为:
设置滑模控制参数及满足线性变换矩阵P21中的参数值,根据三个通道的相对位置和相对速度收敛误差的表达式,预估控制器(11)式在x,y,z三个通道的相对位置和相对速度的控制精度。
实施例1
步骤1.1:假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:其中eiρ=[ex ey ez]T为位置误差,/>为速度误差。已知欠迹向驱动动力学方程为其中,
将其改写为由一个3维变量和一个2维变量的形式:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T (4)
步骤1.2:由于为三阶变量而ei2a为二阶变量,于是设计一个2×3阶非预设参数线性变换矩阵:
其中,a3,b3,f3为常数。
对进行线性变换得到一个2阶状态变量:
步骤2.1:通过线性变量设计滑模面:
其中,c2>0为常数,并且可知:
步骤2.2:对滑模面求导:
则可以得到参数形式的效控制律ui1:
其中,
Θ1,Θ2分别表示Θ的第一项和第二项。
开关控制选择为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2,(11)
步骤3.1:选取Lyapunovp函数V对时间求导后得到:
其中,δi=(c2P21+P21A21)di2u+di2a,||δi||≤δi2m=(c2||P12||+||P12A21||+1)dim,而(12)式可以转化为:
上式可以进一步改写为如下两种形式:
所以当k2-δim||si||-γ>0或k1-δim||si||-1>0成立时,闭环系统将在有限时间收敛到期望轨迹,滑模面的收敛误差重构精度可以被描述为:
对于式(7),当滑模面收敛至时,系统的动力学方程可以构建为:
其中,ζip=-P21dj2u,且有ζip≤-||P21||dim。进一步,将上式改写为如下两种形式:
当和/>成立时,线性状态变量/>和/>的重构精度的表达式为:
步骤3.2:接下来根据线性变换矩阵P21求取闭环系统在x,y,z三个通道的收敛误差。根据式(6)可以得到所以法向的相对位置收敛误差可以表示为:
对(6)式第二项求导,可得法向的相对速度收敛误差:
根据(6)式第一项可得一个一阶系统:
其中,根据一阶系统稳定的充要条件,可知当
a21=b3/a3>0,(25)
此时式(24)可以写为如下两种表达式:
当或/>成立时,迹向的相对位置误差和相对速度收敛误差为:
对欠迹向驱动系统在迹向的动力学表达式积分:
令上式可表示为:
对于自然构型编队,整个重构过程一直满足关系因此故(30)式化简为:
于是在径向的相对位置收敛误差为:
径向的相对速度收敛误差为:
步骤3.3:线性变换矩阵P21和滑模控制法的控制器参数设置为:a3=-454.5,b3=-0.4,f3=1,c3=0.003γ=0.5,k1=0.003,k2=0.000001,c2=0.003。假设扰动的边界值为δi2m=1×10-6m/s2,根据式(16)、(20)-(23)、(27)-(28)和(32)-(33)可以预估滑模面、线性化后的相对位置误差、线性化后相对速度误差及径向、迹向和法向的相对位置和相对速度的收敛误差上届依次为: 即在实际的仿真中,重构后期的稳定重构精度小于上述数值。
仿真试验
本发明的实验结果如下:
选择一般圆轨道为当前队形轨道,投影圆轨道为期望构型轨道:
航天器期望的相对运动位置矢量为当前相对位置矢量一般圆轨道和投影圆轨道的半径分别设置为rig=500m,rip=1000m,3颗跟随者航天器的初始队形相角/>分别为/>而期望构型的相角/>为/>领航者航天器的轨道要素选择为:半长轴6878.137km、偏心率为0、轨道倾角42deg、升交点赤经-60deg、纬度幅角30deg。
图2(a)-(c)为3颗跟随者航天器在径向和法向的控制输入,数值大小为10-3量级。图3(a)-(c)表示3颗跟随者航天器相对于领航者的速度变化情况,相对速度的大小影响着控制器参数的设置,在本实施案例中,在10-3量级控制输入对应于100量级相对运动速度。图4(a)-(c)中三颗跟随者航天器在xy、xz、yz平面上的轨迹;图4(d)为3D重构轨迹,每颗跟随者在图2所示径向和法向两个方向的作用下转移至目标点,最终在投影圆轨道上形成了期望的编队构型,如果控制过大则重构轨迹将变成直线。
对于本领域的技术人员来说,可以根据以上的技术方案和构思,给出各种相应的改变和变形,而所有的这些改变和变形,都应该包括在本发明权利要求的保护范围之内。
Claims (1)
1.一种非预设参数的欠迹向驱动多航天器编队重构控制方法,其特征在于,所述方法包括:
S1将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的非预设参数线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;
所述S1包括:
S1.1假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:
其中,eiρ=[ex ey ez]T为x,y,z三个方向的位置误差,为相应的速度误差;将欠迹向驱动动力学方程分解为一个3维和一个2维变量构成的状态表达式:
其中,
A24=[02×2];
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T;
S1.2由于而ei2a∈R2,所以设计一个2×3阶的非预设参数线性变换矩阵:
其中,a3,b3,f3均为常数;
为了使状态变量的自由度与控制输入的自由度相匹配,使用P21对进行线性变换得到一个2阶状态变量:
S2利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;
利用经线性化变换后的变量设计滑模面,其中:
S2.1根据滑模面设计规则选择如下滑模面s:
其中,c2>0为常数;
S2.2等效控制律ui1可通过对滑模面求导得到参数形式的等效控制律:
其中,
Θ1,Θ2分别表示Θ的第一项和第二项;
选择开关控制为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2;
S3分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的非预设未知参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度;
所述S3包括:
S3.1当k2-δim||si||-γ>0或k1-δim||si||-1>0成立时,闭环系统将在有限时间收敛到期望轨迹,滑模面的收敛误差精度可以被描述为:
而线性化后的状态变量和/>的收敛误差表达式为:
S3.2当线性变换矩阵P21中的参数满足
a21=b3/a3>0;
可以的得到相对位置跟踪误差在x,y,z三个通道的收敛误差为:
相对位置和相对运动速度跟踪误差在迹向的收敛误差为:
相对位置和相对运动速度跟踪误差在法向的收敛误差为:
设置滑模控制参数及满足线性变换矩阵P21中的参数值,根据三个通道的相对位置和相对速度收敛误差的表达式,预估控制器在x,y,z三个通道的相对位置和相对速度的控制精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211012265.4A CN115390587B (zh) | 2022-08-23 | 2022-08-23 | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211012265.4A CN115390587B (zh) | 2022-08-23 | 2022-08-23 | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115390587A CN115390587A (zh) | 2022-11-25 |
CN115390587B true CN115390587B (zh) | 2023-12-15 |
Family
ID=84120418
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211012265.4A Active CN115390587B (zh) | 2022-08-23 | 2022-08-23 | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115390587B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105739511A (zh) * | 2016-01-31 | 2016-07-06 | 中国人民解放军国防科学技术大学 | 一种缺失迹向控制的欠驱动航天器悬停渐近控制方法 |
CN105786013A (zh) * | 2016-03-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 面向输入饱和的迹向欠驱动航天器编队重构控制方法 |
CN106773713A (zh) * | 2017-01-17 | 2017-05-31 | 北京航空航天大学 | 针对欠驱动海洋航行器的高精度非线性路径跟踪控制方法 |
CN114237272A (zh) * | 2021-12-15 | 2022-03-25 | 中国空间技术研究院 | 一种多卫星同步编队重构控制方法 |
-
2022
- 2022-08-23 CN CN202211012265.4A patent/CN115390587B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105739511A (zh) * | 2016-01-31 | 2016-07-06 | 中国人民解放军国防科学技术大学 | 一种缺失迹向控制的欠驱动航天器悬停渐近控制方法 |
CN105786013A (zh) * | 2016-03-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 面向输入饱和的迹向欠驱动航天器编队重构控制方法 |
CN106773713A (zh) * | 2017-01-17 | 2017-05-31 | 北京航空航天大学 | 针对欠驱动海洋航行器的高精度非线性路径跟踪控制方法 |
CN114237272A (zh) * | 2021-12-15 | 2022-03-25 | 中国空间技术研究院 | 一种多卫星同步编队重构控制方法 |
Non-Patent Citations (2)
Title |
---|
"Formation Reconfiguration Hovering Technique for Underactuated Satellite";Jiang Shao 等;Proceedings of the 2021 5th Chinese Conference on Swarm Intelligence and Cooperative Control;第439-448页 * |
"欠驱动航天器编队重构的切换神经网络控制";于金龙 等;宇航学报;第43卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115390587A (zh) | 2022-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Horri et al. | Practical implementation of attitude-control algorithms for an underactuated satellite | |
JP4511390B2 (ja) | 人工衛星の姿勢制御装置 | |
CN109782759B (zh) | 一种轮式移动机器人的近似解耦、快速轨迹跟踪控制方法 | |
Dinh et al. | Adaptive tracking control of a quadrotor unmanned vehicle | |
CN107263466A (zh) | 空间机器人基于二次规划问题的基座无扰控制方法 | |
CN103869704A (zh) | 基于扩展雅克比矩阵的空间机器人星臂协调控制方法 | |
Zha et al. | Quaternion-based nonlinear trajectory tracking control of a quadrotor unmanned aerial vehicle | |
CN111506095B (zh) | 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法 | |
CN112550770A (zh) | 一种基于凸优化的火箭软着陆轨迹规划方法 | |
CN116277036B (zh) | 一种柔性基、柔性臂空间机器人快速容错抑振控制方法 | |
Sun | Saturated adaptive output-constrained control of cooperative spacecraft rendezvous and docking | |
Pollini et al. | Simulation and robust backstepping control of a quadrotor aircraft | |
CN115390587B (zh) | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 | |
CN113485396A (zh) | 一种航天器交会对接最后逼近段相对轨道与姿态跟踪控制方法 | |
CN114237272A (zh) | 一种多卫星同步编队重构控制方法 | |
CN109062044B (zh) | 一种终端迭代学习对接控制方法 | |
CN115390576B (zh) | 一种非预设参数的欠径向推力多卫星编队重构控制方法 | |
Benosman et al. | Output trajectory tracking for a switched nonlinear non-minimum phase system: The vstol aircraft | |
Sah et al. | Optimal finite thrust rendezvous maneuvers of chaser spacecraft towards capturing target object | |
Yoshimura et al. | Global trajectory design for position and attitude control of an underactuated satellite | |
Бурнашев et al. | Micro-satellite orientation control system with rotation wheels | |
CN111897352B (zh) | 基于控制力矩输出约束的姿态机动动态优化方法及介质 | |
Kang et al. | Robust nonsingular fixed time terminal sliding mode control for atmospheric pollution detection lidar scanning mechanism | |
CN112346472B (zh) | 航天器交会系统的保成本鲁棒增益调度控制器设计方法 | |
Yin et al. | Attitude stability control of capturing spacecraft without angular velocity |
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 |