CN105005312A - 一种基于最大角加速度和最大角速度卫星规划轨迹方法 - Google Patents
一种基于最大角加速度和最大角速度卫星规划轨迹方法 Download PDFInfo
- Publication number
- CN105005312A CN105005312A CN201510367515.XA CN201510367515A CN105005312A CN 105005312 A CN105005312 A CN 105005312A CN 201510367515 A CN201510367515 A CN 201510367515A CN 105005312 A CN105005312 A CN 105005312A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- phi
- max
- maximum angular
- planned trajectory
- 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
Landscapes
- External Artificial Organs (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种基于最大角加速度和最大角速度卫星规划轨迹方法,属于卫星机动轨迹规划领域。现有的规划轨迹确定方法不能充分利用执行机构的机动能力,且不能保证机动时间最短的问题。一种基于最大角加速度和最大角速度卫星规划轨迹方法,设定与目标姿态对应的目标坐标系,计算卫星由初始姿态机动至目标姿态的欧拉轴em和转角Φm;获得规划轨迹的最大角加速度和最大角速度的约束方程;由表示获得受飞轮最大角动量限制的计算使机动时间tm取最小值时规划轨迹的最大角速度并通过规划轨迹的最大角速度求出规划轨迹的最大角加速度本发明能够保证规划轨迹充分利用飞轮的能力,以使机动时间最短。
Description
技术领域
本发明涉及一种基于最大角加速度和最大角速度卫星规划轨迹方法。
背景技术
对地观测卫星特别是具有快速响应能力的成像卫星,在轨运行时需要迅速获取地面不同目标信息或者地面同一目标连续/间断信息,因此,这类卫星都必须具备快速大角度姿态机动的能力。
卫星实现大角度姿态机动的前提是机动轨迹规划。在以往的研究中,研究人员往往人为设定规划轨迹的最大角加速度和最大角速度,通过试凑的方法,使设计出的规划轨迹在执行机构的能力范围之内,以实现大角度机动。这类方法给出的规划轨迹往往不能充分利用执行机构的机动能力,而且不能保证机动时间最短。因此,有必要给出规划轨迹的最大角加速度和最大角速度的设计方法,使执行机构的能力得到充分利用,从而最大程度地缩短机动时间。
发明内容
本发明的目的是为了解决现有的规划轨迹确定方法不能充分利用执行机构的机动能力,而且不能保证机动时间最短的问题,而提出一种基于最大角加速度和最大角速度卫星规划轨迹方法。
一种基于最大角加速度和最大角速度卫星规划轨迹方法,所述方法通过以下步骤实现:
步骤一、设定与目标姿态对应的目标坐标系obxtytzt,根据:
Φm=2arccosqm0 (2),
计算卫星由初始姿态机动至目标姿态的欧拉轴em和转角Φm;式中,qm0表示卫星由初始姿态机动至目标姿态的机动四元数Qm的标部,qm表示卫星由初始姿态机动至目标姿态的机动四元数Qm的矢部,Qm为四维矢量Qm=[qm0 qm]T,且
式中,Qi表示初始姿态四元数,Qt为表示目标姿态四元数;
步骤二、获得规划轨迹的最大角加速度和最大角速度的约束方程:
其中,
式中,表示规划轨迹的最大角加速度,表示规划轨迹的最大角速度,表示角加速度力矩项的系数,表示角速度力矩项的系数,I是3×3的矩阵,表示卫星相对其质心的转动惯量,em表示欧拉轴,ωoI表示轨道角速度,H0表示初始时刻卫星系统的角动量,Twmax是标量,为飞轮力矩包络半径;
步骤三、将步骤二中涉及的角速度力矩项的系数分解为和两部分,并获得和的模值:
其中,是与平行的角速度力矩项系数,是与垂直的角速度力矩项系数;
步骤四、由规划轨迹的最大角速度表示规划轨迹的最大角加速度
步骤五、获得受飞轮最大角动量限制的最大角速度的表达式:
式中,Hwmax表示飞轮角动量包络半径,下角标H表示最大角速度受飞轮角动量上限约束;
步骤六、计算最短机动时间tm:
并联合式(9)得到以规划轨迹的最大角速度为自变量的机动时间tm的函数:
结合如下约束:
求出使机动时间tm取最小值时规划轨迹的最大角速度
步骤七、将步骤六获得的规划轨迹的最大角速度代入式(9),求出规划轨迹的最大角加速度
本发明的有益效果为:
本发明实现了一种设计卫星规划轨迹的最大角加速度和最大角速度的方法,通过限定项的取值范围,使得飞轮的力矩得以充分利用,通过获得的规划轨迹的最大角速度保证所设计的规划轨迹充分利用飞轮的能力,以使得机动时间最短。
本发明方法在分析飞轮力矩空间的利用过程中设计各明确物理含义的物理量,易于理解;使得本发明方法简单易行,简化计算过程,使得卫星在轨运行时快速获取地面目标的信息和相应。
附图说明
图1为本发明的流程图;
图2为具体实施方式步骤三、四所述在飞轮力矩包络内的力矩示意图;
图3为具体实施方式步骤六所述的规划轨迹的角加速度曲线和角速度曲线;且图中上半部分图示为规划轨迹的角加速度曲线示意图,下半部分图示为规划轨迹的角加速度曲线示意图;
图4为算例的规划轨迹角加速度曲线;
图5为算例的规划轨迹角速度曲线;
图6为算例的卫星姿态角速度曲线;
图7为算例的卫星姿态角曲线;
图8为算例的飞轮力矩曲线;
图9为算例的飞轮力矩模值曲线;
图10为算例的飞轮角动量曲线;
图11为算例的飞轮角动量模值曲线;
具体实施方式
具体实施方式一:
本实施方式的一种基于最大角加速度和最大角速度卫星规划轨迹方法,结合图1至图3所示内容,所述方法通过以下步骤实现:
步骤一、设定与目标姿态对应的目标坐标系obxtytzt,根据:
Φm=2arccosqm0 (2),
计算卫星由初始姿态机动至目标姿态的欧拉轴em和转角Φm;式中,qm0表示卫星由初始姿态机动至目标姿态的机动四元数Qm的标部,qm表示卫星由初始姿态机动至目标姿态的机动四元数Qm的矢部,Qm为四维矢量Qm=[qm0 qm]T,且
式中,Qi表示初始姿态四元数,Qt为表示目标姿态四元数;
步骤二、获得规划轨迹的最大角加速度和最大角速度的约束方程:
其中,
式中,表示规划轨迹的最大角加速度,表示规划轨迹的最大角速度,表示角加速度力矩项的系数,表示角速度力矩项的系数,I是3×3的矩阵,表示卫星相对其质心的转动惯量,em表示欧拉轴,ωoI表示轨道角速度,H0表示初始时刻卫星系统的角动量,Twmax是标量,为飞轮力矩包络半径;
步骤三、将步骤二中涉及的角速度力矩项的系数分解为和两部分,并获得和的模值:
其中,是与平行的角速度力矩项系数,是与垂直的角速度力矩项系数;
步骤四、由规划轨迹的最大角速度表示规划轨迹的最大角加速度
步骤五、获得受飞轮最大角动量限制的最大角速度的表达式:
式中,Hwmax表示飞轮角动量包络半径,下角标H表示最大角速度受飞轮角动量上限约束;
步骤六、计算最短机动时间tm:
并联合式(9)得到以规划轨迹的最大角速度为自变量的机动时间tm的函数:
结合如下约束:
求出使机动时间tm取最小值时规划轨迹的最大角速度
步骤七、将步骤六获得的规划轨迹的最大角速度代入式(9),求出规划轨迹的最大角加速度
具体实施方式二:
与具体实施方式一不同的是,本实施方式的一种基于最大角加速度和最大角速度卫星规划轨迹方法,步骤二所述获得规划轨迹的最大角加速度和最大角速度的约束方程过程为,忽略外界干扰力矩的影响,卫星的动力学方程为:通常,卫星的一次机动时间远小于其轨道周期,那么,在机动过程中,卫星本体系相对轨道系的角速度ωbo的模值远大于轨道角速度ωo的模值,再结合角动量守恒定理:IωbI+Hw=H0,动力学方程简化为:
假设卫星能够对规划轨迹进行很好地跟踪,则:
将式(16),式(17)代入式(15),得到直观体现角加速度和角速度与飞轮力矩Tw关系的动力学方程: 令最大角加速度和最大角速度的取值使得式(18)右端的飞轮力矩Tw的模值达到上限,以充分利用飞轮的力矩空间,即:
(-Iem)和(Iem ×ωoI-em ×H0)为已知量,令(-Iem)由表示,(Iem ×ωoI-em ×H0)由表示,将式(19)简化为:
具体实施方式三:
与具体实施方式一或二不同的是,本实施方式的一种基于最大角加速度和最大角速度卫星规划轨迹方法,步骤三所述获得和的模值的过程为,如图2所示,由于由平行于的和垂直于的两部分组成,从而有:
具体实施方式四:
与具体实施方式三不同的是,本实施方式的一种基于最大角加速度和最大角速度卫星规划轨迹方法,步骤四所述由规划轨迹的最大角速度表示规划轨迹的最大角加速度的过程为,由图2和式(19)可知,最大角速度的存在使飞轮力矩Tw无法被充分利用,导致规划轨迹的最大角加速度的可用力矩包络从Twmax缩减为Twmax-△T,则减少量△T表示为:
令sinα≈α,以充分利用飞轮力矩Tw的空间,加快机动速度,使角加速度力矩项尽量大,则式(20)近似的表示为:
根据图2可知,规划轨迹的最大角加速度可表示为:
具体实施方式五:
与具体实施方式一、二或四不同的是,本实施方式的一种基于最大角加速度和最大角速度卫星规划轨迹方法,步骤五所述获得受飞轮最大角动量限制的最大角速度的表达式的过程为,
由于卫星在机动过程中,卫星本体角动量的变化量是由飞轮吸收获得,所以,规划轨迹的最大角速度不能超过飞轮转速的上限;根据动量守恒定理:IωbI+Hw=H0,并假设在机动过程中ωbo的模值远大于ωo的模值,得到飞轮角动量上限对最大角速度的限制方程的表达式:
具体实施方式六:
与具体实施方式五不同的是,本实施方式的一种基于最大角加速度和最大角速度卫星规划轨迹方法,步骤六所述求出使机动时间tm取最小值的规划轨迹的最大角速度的过程为,
第一,根据图3得机动时间tm表示为:
将式(9)代入式(11),得到以规划轨迹的最大角速度为自变量的机动时间tm的函数:
控制最大角速度满足
第二,由图3所示的角速度和角加速度的变化规律知,规划轨迹的最大角速度和规划轨迹的最大角加速度的取值满足即: 避免在跟踪过程中角速度未增加到最大值前开始减速;
第三,利用Matlab工具箱中的非线性优化方法,得到带有约束的机动时间tm的最小值及其对应的规划轨迹的最大角速度
实施例:
结合图4~11说明本具体实施例。
进行卫星姿态轨迹跟踪的仿真,说明本发明设计的规划轨迹的最大加速度和最大速度可以充分利用执行机构,使机动时间最短。卫星初始姿态为[0° 10° 10°]T,目标姿态为[0° 30° 40°]T。飞轮采用四斜装构型,单个飞轮可产生的力矩和角动量分别为0.2Nm,15Nm·s,则飞轮的力矩包络半径和角动量包络半径分别为0.267Nm,20Nm·s。卫星转动惯量为I=diag(2000 1000 3000)kg·m2,轨道角速度为ωoI=[0 -0.001 0]Trad/s,系统初始角动量为H0=[0 0 0]TNm·s。姿态角速度测量误差为三轴10-4°/s,姿态角测量误差为三轴5×10-4°。干扰力矩为 PD控制器参数为Kp=0.04I3,Kd=0.36I3,I3为3×3的单位矩阵。
根据步骤一,得到em=[-0.287 0.544 0.789]T,Φm=0.627rad。
根据步骤二、三,得到
根据步骤五,得到
根据步骤六,得到tm=134.702s,
根据步骤七,得到
由图5~7可以得出,卫星的姿态角速度和姿态角能够很好的跟踪规划轨迹,而且,卫星在加速过程结束后,经过短暂的匀速过程,立马进入减速过程,显然能使机动时间尽量缩短。由图8~9可以看出,飞轮力矩模值为0.237Nm,而力矩包络为0.267Nm,力矩空间能够较为充分的利用。由图10~11可以看出,飞轮角动量模值为13.031Nm·s,而角动量包络为20Nm·s,角动量空间的利用不如力矩空间利用得充分,原因在于,根据角动量约束求得的 而根据 约束求得的 的值减少了15%,所以使得角动量空间不能得到有效利用。
Claims (6)
1.一种基于最大角加速度和最大角速度卫星规划轨迹方法,其特征在于:所述方法通过以下步骤实现:
步骤一、设定与目标姿态对应的目标坐标系obxtytzt,根据:
Φm=2arccosqm0 (2),
计算卫星由初始姿态机动至目标姿态的欧拉轴em和转角Φm;式中,qm0表示卫星由初始姿态机动至目标姿态的机动四元数Qm的标部,qm表示卫星由初始姿态机动至目标姿态的机动四元数Qm的矢部,Qm为四维矢量Qm=[qm0 qm]T,且
式中,Qi表示初始姿态四元数,Qt为表示目标姿态四元数;
步骤二、获得规划轨迹的最大角加速度和最大角速度的约束方程:
其中,
NΦ=Iem ×ωoI-em ×H0 (6),
式中,表示规划轨迹的最大角加速度,表示规划轨迹的最大角速度,表示角加速度力矩项的系数,表示角速度力矩项的系数,I是3×3的矩阵,表示卫星相对其质心的转动惯量,em表示欧拉轴,ωoI表示轨道角速度,H0表示初始时刻卫星系统的角动量,Twmax是标量,为飞轮力矩包络半径;
步骤三、将步骤二中涉及的角速度力矩项的系数分解为和两部分,并获得和的模值:
其中,是与平行的角速度力矩项系数,是与垂直的角速度力矩项系数;
步骤四、由规划轨迹的最大角速度表示规划轨迹的最大角加速度
步骤五、获得受飞轮最大角动量限制的最大角速度的表达式:
式中,Hwmax表示飞轮角动量包络半径,下角标H表示最大角速度受飞轮角动量上限约束;
步骤六、计算最短机动时间tm:
并联合式(9)得到以规划轨迹的最大角速度为自变量的机动时间tm的函数:
求出使机动时间tm取最小值时规划轨迹的最大角速度
步骤七、将步骤六获得的规划轨迹的最大角速度代入式(9),求出规划轨迹的最大角加速度
2.根据权利要求1所述一种基于最大角加速度和最大角速度卫星规划轨迹方法,其特征在于:步骤二所述获得规划轨迹的最大角加速度和最大角速度的约束方程过程为,忽略外界干扰力矩的影响,卫星的动力学方程为:在机动过程中,卫星本体系相对轨道系的角速度ωbo的模值大于轨道角速度ωo的模值,结合角动量守恒定理:IωbI+Hw=H0,卫星的动力学方程式(14)简化为:
设卫星对规划轨迹进行跟踪,则:
将式(16),式(17)代入式(15),得到直观体现角加速度和角速度与飞轮力矩Tw关系的动力学方程:
令最大角加速度和最大角速度的取值使得式(18)右端的飞轮力矩Tw的模值达到上限,以充分利用飞轮的力矩空间,即: (-Iem)和(Iem ×ωoI-em ×H0)为已知量,令(-Iem)由表示,(Iem ×ωoI-em ×H0)由表示,将式(19)简化为:
3.根据权利要求1或2所述一种基于最大角加速度和最大角速度卫星规划轨迹方法,其特征在于:步骤三所述获得和的模值的过程为,由于由平行于的和垂直于的两部分组成,从而有:
4.根据权利要求3所述一种基于最大角加速度和最大角速度卫星规划轨迹方法,其特征在于:步骤四所述由规划轨迹的最大角速度表示规划轨迹的最大角加速度的过程为,由式(19)可知,规划轨迹的最大角加速度的可用力矩包络从Twmax缩减为Twmax-△T,则减少量△T表示为
令sinα≈α,则式(20)近似的表示为:
规划轨迹的最大角加速度可表示为:
5.根据权利要求1、2或4所述一种基于最大角加速度和最大角速度卫星规划轨迹方法,其特征在于:步骤五所述获得受飞轮最大角动量限制的最大角速度的表达式的过程为,
卫星在机动过程中,规划轨迹的最大角速度不能超过飞轮转速的上限;根据动量守恒定理:IωbI+Hw=H0,并假设在机动过程中ωbo的模值远大于ωo的模值,得到飞轮角动量上限对最大角速度的限制方程的表达式:
6.根据权利要求5所述一种基于最大角加速度和最大角速度卫星规划轨迹方法,其特征在于:步骤六所述求出使机动时间tm取最小值的规划轨迹的最大角速度的过程为,
第一,机动时间tm表示为:
将式(9)代入式(11),得到以规划轨迹的最大角速度为自变量的机动时间tm的函数:
控制最大角速度满足:
第二,由角速度和角加速度的变化规律知,规划轨迹的最大角速度和规划轨迹的最大角加速度的取值满足即:
第三,利用Matlab工具箱中的非线性优化方法,得到带有约束的机动时间tm的最小值及其对应的规划轨迹的最大角速度
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510367515.XA CN105005312B (zh) | 2015-06-29 | 2015-06-29 | 一种基于最大角加速度和最大角速度卫星规划轨迹方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510367515.XA CN105005312B (zh) | 2015-06-29 | 2015-06-29 | 一种基于最大角加速度和最大角速度卫星规划轨迹方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105005312A true CN105005312A (zh) | 2015-10-28 |
CN105005312B CN105005312B (zh) | 2017-11-03 |
Family
ID=54378019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510367515.XA Active CN105005312B (zh) | 2015-06-29 | 2015-06-29 | 一种基于最大角加速度和最大角速度卫星规划轨迹方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105005312B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106200660A (zh) * | 2016-07-29 | 2016-12-07 | 北京控制工程研究所 | 一种用于同步轨道sar卫星的姿态机动轨迹计算方法 |
CN108319143A (zh) * | 2018-02-11 | 2018-07-24 | 北京控制工程研究所 | 一种航天器对动坐标系机动目标实时规划方法 |
CN109039222A (zh) * | 2018-07-19 | 2018-12-18 | 中国科学院长春光学精密机械与物理研究所 | 一种加速力矩受限的转台运动曲线在线规划方法及系统 |
CN109407691A (zh) * | 2019-01-08 | 2019-03-01 | 深圳航天东方红海特卫星有限公司 | 一种微小卫星精密闭环轨道控制设备 |
CN109634294A (zh) * | 2018-12-10 | 2019-04-16 | 上海航天控制技术研究所 | 一种基于机动能力辨识的抗干扰四元数姿态机动路径规划方法 |
CN109774977A (zh) * | 2019-03-28 | 2019-05-21 | 上海微小卫星工程中心 | 一种基于四元数的时间最优的卫星姿态快速机动方法 |
CN113761664A (zh) * | 2021-09-15 | 2021-12-07 | 长光卫星技术有限公司 | 一种遥感卫星成像期间飞轮转速优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101934863A (zh) * | 2010-09-29 | 2011-01-05 | 哈尔滨工业大学 | 基于磁力矩器和飞轮的卫星姿态全方位控制方法 |
CN103235597A (zh) * | 2013-04-09 | 2013-08-07 | 北京理工大学 | 一种航天器的快速姿态机动快速稳定联合控制方法 |
CN103808323A (zh) * | 2012-11-07 | 2014-05-21 | 上海航天控制工程研究所 | 一种卫星姿态跟踪机动的余弦过渡加速度路径方法 |
CN103940433A (zh) * | 2014-05-12 | 2014-07-23 | 哈尔滨工业大学 | 一种基于改进的自适应平方根ukf算法的卫星姿态确定方法 |
-
2015
- 2015-06-29 CN CN201510367515.XA patent/CN105005312B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101934863A (zh) * | 2010-09-29 | 2011-01-05 | 哈尔滨工业大学 | 基于磁力矩器和飞轮的卫星姿态全方位控制方法 |
CN103808323A (zh) * | 2012-11-07 | 2014-05-21 | 上海航天控制工程研究所 | 一种卫星姿态跟踪机动的余弦过渡加速度路径方法 |
CN103235597A (zh) * | 2013-04-09 | 2013-08-07 | 北京理工大学 | 一种航天器的快速姿态机动快速稳定联合控制方法 |
CN103940433A (zh) * | 2014-05-12 | 2014-07-23 | 哈尔滨工业大学 | 一种基于改进的自适应平方根ukf算法的卫星姿态确定方法 |
Non-Patent Citations (3)
Title |
---|
CHENFENG ZHOU, 等: "Attitude Control of Flexible Spacecraft Considering Rigid-Flex Coupling Parameters", 《PROCEEDINGS OF 2012 INTERNATIONAL CONFERENCE ON MODELLING, IDENTIFICATION AND CONTROL》 * |
HYUNGJOO YOON,等: "Adaptive Spacecraft Attitude Tracking Control with Actuator Uncertainties", 《THE JOURNAL OF THE ASTRONAUTICAL SCIENCES》 * |
耿云海,等: "使用混合执行机构的快速机动卫星力矩分配算法", 《宇航学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106200660A (zh) * | 2016-07-29 | 2016-12-07 | 北京控制工程研究所 | 一种用于同步轨道sar卫星的姿态机动轨迹计算方法 |
CN108319143A (zh) * | 2018-02-11 | 2018-07-24 | 北京控制工程研究所 | 一种航天器对动坐标系机动目标实时规划方法 |
CN109039222A (zh) * | 2018-07-19 | 2018-12-18 | 中国科学院长春光学精密机械与物理研究所 | 一种加速力矩受限的转台运动曲线在线规划方法及系统 |
CN109634294A (zh) * | 2018-12-10 | 2019-04-16 | 上海航天控制技术研究所 | 一种基于机动能力辨识的抗干扰四元数姿态机动路径规划方法 |
CN109634294B (zh) * | 2018-12-10 | 2022-01-04 | 上海航天控制技术研究所 | 一种基于机动能力辨识的抗干扰四元数姿态机动路径规划方法 |
CN109407691A (zh) * | 2019-01-08 | 2019-03-01 | 深圳航天东方红海特卫星有限公司 | 一种微小卫星精密闭环轨道控制设备 |
CN109407691B (zh) * | 2019-01-08 | 2022-03-18 | 深圳航天东方红海特卫星有限公司 | 一种微小卫星精密闭环轨道控制设备 |
CN109774977A (zh) * | 2019-03-28 | 2019-05-21 | 上海微小卫星工程中心 | 一种基于四元数的时间最优的卫星姿态快速机动方法 |
CN113761664A (zh) * | 2021-09-15 | 2021-12-07 | 长光卫星技术有限公司 | 一种遥感卫星成像期间飞轮转速优化方法 |
CN113761664B (zh) * | 2021-09-15 | 2024-01-30 | 长光卫星技术股份有限公司 | 一种遥感卫星成像期间飞轮转速优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105005312B (zh) | 2017-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105005312A (zh) | 一种基于最大角加速度和最大角速度卫星规划轨迹方法 | |
CN103092208B (zh) | 基于sgcmg和rw的航天器高精度快速姿态机动方法 | |
CN104990549B (zh) | 一种采煤机‑液压支架相对定位的方法及装置 | |
CN103759742B (zh) | 基于模糊自适应控制技术的捷联惯导非线性对准方法 | |
CN102163057B (zh) | 一种自治飞艇平面路径跟踪控制方法 | |
CN103955223B (zh) | 一种深空探测软着陆过程的姿轨耦合控制方法 | |
CN102878872B (zh) | 一种针对导引头失锁情况的制导信息处理方法 | |
CN105929842A (zh) | 一种基于动态速度调节的欠驱动uuv平面轨迹跟踪控制方法 | |
CN103293957A (zh) | 一种相对动坐标系进行路径规划的卫星姿态机动方法 | |
CN104281150A (zh) | 一种姿态机动的轨迹规划方法 | |
CN104932512A (zh) | 一种基于mimo非线性不确定反步法的四旋翼位姿控制方法 | |
CN103231810B (zh) | 一种利用卫星俯仰轴姿态机动卸载俯仰轴角动量的方法 | |
CN105629732B (zh) | 一种考虑控制受限的航天器姿态输出反馈跟踪控制方法 | |
CN103955218A (zh) | 一种基于非线性控制理论的无人艇轨迹跟踪控制装置及方法 | |
CN104443432B (zh) | 一种卫星有限推力共面圆轨道自主轨道转移制导方法 | |
CN109211230B (zh) | 一种基于牛顿迭代法的炮弹姿态和加速度计常值误差估计方法 | |
CN102923317B (zh) | 一种适用于卫星姿态角速度阻尼的欠驱动控制方法 | |
CN104730506B (zh) | 一种合成孔径雷达卫星的全零多普勒姿态导引方法 | |
CN104749955B (zh) | 一种独轮机器人的有限时间自平衡控制方法 | |
CN108319143A (zh) | 一种航天器对动坐标系机动目标实时规划方法 | |
CN104176275A (zh) | 一种使用动量轮与磁力矩器联合的速率阻尼方法 | |
CN107063244A (zh) | 一种飞行器飞行过程模拟方法 | |
CN104085539A (zh) | 成像定标的姿态控制方法 | |
CN102435193B (zh) | 一种捷联惯性导航系统的高精度初始对准方法 | |
CN102323759A (zh) | 翼伞自主归航半实物仿真系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |