CN112115545B - 一种基于二项式近似模型的挠性卫星姿态机动控制方法 - Google Patents

一种基于二项式近似模型的挠性卫星姿态机动控制方法 Download PDF

Info

Publication number
CN112115545B
CN112115545B CN202010798931.6A CN202010798931A CN112115545B CN 112115545 B CN112115545 B CN 112115545B CN 202010798931 A CN202010798931 A CN 202010798931A CN 112115545 B CN112115545 B CN 112115545B
Authority
CN
China
Prior art keywords
satellite
attitude
equation
flexible
sampling
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
Application number
CN202010798931.6A
Other languages
English (en)
Other versions
CN112115545A (zh
Inventor
吕金虎
陈磊
刘克新
谷海波
朱国梁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202010798931.6A priority Critical patent/CN112115545B/zh
Publication of CN112115545A publication Critical patent/CN112115545A/zh
Application granted granted Critical
Publication of CN112115545B publication Critical patent/CN112115545B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • B64G1/245Attitude control algorithms for spacecraft attitude control
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Remote Sensing (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Combustion & Propulsion (AREA)
  • Chemical & Material Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提供一种基于二项式近似模型的挠性卫星姿态机动控制方法,可以在部分系统参数的条件下利用对输出的采样和输入的零阶保持实现对卫星姿态的快速转动控制。此方法可应用于卫星姿态系统上,既可用来实现姿态稳定也可用来跟踪既定目标,且因控制算法复杂度低计算效率高,不同于传统控制算法可用于系统参数未知系统,对各类外界干扰具有较强的鲁棒性;具有应用灵活、功能较强的特点,具有广阔的应用前景。

Description

一种基于二项式近似模型的挠性卫星姿态机动控制方法
技术领域
本发明涉及导卫星姿态控制领域,尤其基于简化模型的卫星姿态机动控 制方法。
背景技术
随着现代卫星的快速发展,卫星的姿态机动能力得到了越来越多的关注。 对卫星姿态控制律的可靠性提出了更高的要求。
目前,采用传统的控制器对简单卫星施加控制,实现姿态机动。卫星的 任务需求增加造成卫星的结构日趋复杂,形成多结构的挠性体。由于传统的 控制方法在面对复杂挠性体时有局限性,在现代卫星应用上有潜在风险。
因此,如何在低成本、易移植的情况下,实现挠性卫星姿态机动成为亟 待解决的问题。
发明内容
本发明提供一种基于二项式近似模型的挠性卫星姿态机动控制方法,用 以解决如何在低成本、易移植的情况下,完成挠性卫星姿态机动这一问题。
一种基于二项式近似模型的挠性卫星姿态机动控制方法,其特征在于, 包括如下步骤:
步骤一、设卫星飞轮力矩遵从采样控制律,基于卫星三轴姿态的采样值、 期望路径规划及卫星本体转动惯量、挠性附件耦合系数,构建卫星姿态机动 系统的近似模型,具体为:
挠性卫星姿态系统的精确模型为:
Figure BDA0002626660460000011
其中,
Figure BDA0002626660460000021
代表卫星三轴姿态角速度;φ,θ,
Figure BDA0002626660460000022
分别 为卫星滚动、俯仰、偏航轴;t为连续时间变量;
Figure BDA0002626660460000023
为三轴姿态角速度一 阶导数;I为卫星刚体转动惯量;F为卫星刚体与挠性附件的耦合矩阵系数; η(t)为挠性附件的振动模态坐标向量;
Figure BDA0002626660460000024
Figure BDA0002626660460000025
分别为挠性附件的振动模 态坐标向量的一阶与二阶导数;
Figure BDA00026266604600000214
为挠性附件的振动模态频率矩阵;ξ为挠 性附件的振动阻尼比;u(t)为卫星飞轮力矩;并且有
Figure BDA0002626660460000026
将式(1)中的第二个模态方程改写为:
Figure BDA0002626660460000027
其中
Figure BDA0002626660460000028
为挠性附件振动的模态频率矩阵的平方;将式(3)代入式(1)中的第 一个方程得到
Figure BDA0002626660460000029
利用采样时间为T的采样器对三轴姿态角进行如下采样:
Figure BDA00026266604600000210
其中k为自然数,代表离散时间采样点,
Figure BDA00026266604600000211
表示向下取整;与此同时, 控制输入根据离散控制律进行如下零阶保持:
u(t)=u(kT)=u[k] (6)
如此式(4)得到离散为:
Figure BDA00026266604600000212
其中,ω[k],u[k],
Figure BDA00026266604600000213
η[k]分别代表三轴姿态角速度、飞轮力矩、挠 性附件的振动模态坐标向量一阶导数及自身在时间k处的采样值,ω[k+1]代 表下一刻k处的三轴姿态角速度采样值,ω×[k]为式(2)的k处的采样值;
给定移位算子函数Δ:ω[k]→ω[k]-ω[k-1],其中:表示将函数映射到三 轴姿态角速度采样值上,ω[k-1]代表三轴姿态角速度前一步的采样值,将移 位算子在式(7)等号两端映射n-1次得到:
Figure BDA0002626660460000031
其中Θ[k]为卫星姿态机动系统的近似模型的冗余项,
Figure BDA0002626660460000032
为二项式系数, A=(I-FFT)-1为反馈矩阵系数,n≥2为采样时间及卫星姿态机动系统的近 似模型阶次,一般应为自然数,并且是可调参数,具体为:
Figure BDA0002626660460000033
其中Ω[k]为卫星自身参数项,具体为:
Figure BDA0002626660460000034
考虑到挠性附件的振动模态与卫星姿态机动系统路径的连续性,式(10) 进一步表达为:
Θ[k]=TnΩ(n-1)[k] (11)
由于Θ[k]正比于采样时间T的n次方,当采样时间T<1时,忽略此项;由 此得到卫星姿态机动系统的近似模型为:
Figure BDA0002626660460000035
步骤二、将步骤一中得到的卫星姿态机动系统的近似模型
转化为期望路径误差模型,基于此模型设计控制器,使得卫星飞轮根据 指定信号旋转完成卫星三轴跟踪期望路径并稳定,具体为:
设三轴期望路径为v(t);依据式(6)对期望路径进行采样后,从式(1)中的 左右两端减去v[k+1]得到期望路径的误差模型:
Figure BDA0002626660460000036
其中
Figure BDA0002626660460000041
为期望路径误差的估计值,
e[k]=ω[k]-v[k]为期望路径误差的实际值;
式(12)中的
Figure RE-GDA0002767696560000042
表示期望路径计算式的部分 已经忽略;忽略原因为卫星姿态期望路径为连续可导信号,此项等于Tnv(n)[k], 认为足够小或等于0;
完成卫星姿态机动系统任务即使卫星三轴角速度跟踪预期信号,也就是 使式(12)中的期望路径误差
Figure BDA0002626660460000043
设计的控制律为:
Figure BDA0002626660460000044
本发明具有如下有益效果
本发明中的方法可以通过部分卫星系统参数实现对挠性卫星姿态机动的 控制,大大提高了控制器设计的简易性与有效性。此方法可应用于多种类型 的卫星姿态控制上,既可用来进行姿态稳定也可用来进行姿态机动,且因算 法计算复杂度低计算效率高,不同于传统控制算法可用于在线机动优化,对 于具有不确定性的卫星姿态系统有较强的鲁棒性。因此,本发明中的方法具 有应用灵活、功能较强的特点,具有广阔的工程应用前景。
附图说明
图1为本发明的方法流程图。
图2为本发明一实施例的控制效果图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
当前,在挠性卫星姿态控制研究中,如何设计简单易行的控制律完成卫 星的姿态机动是增强卫星应用的重要问题。在本发明中:
如图1所示,首先引入卫星三轴的机动期望路径,并利用采样器对卫星 的姿态陀螺仪的输出进行采样,获得三轴姿态的采样值;结合正负交错的二 项式系数与三轴姿态采样值,构建卫星姿态的离散近似模型,在此基础上与 期望路径的采样值进行移位做差,获得期望误差的近似模型,进一步的根据 此模型设计采样控制律;将此控制律通过零阶保持的方法转化为连续信号后 施加于卫星的飞轮力矩,卫星飞轮将根据指定命令完成卫星的姿态机动。
对符号的注释:本发明中凡是出现符号上面带“·”的表达式,其采用课 本上表达式,即为该符号的求导,例如下文中
Figure BDA0002626660460000051
对ω的求导。
本发明为控制卫星完成姿态机动的方法,以太阳同步卫星为例,如图1 所示,具体的过程为:
步骤一、假设太阳同步卫星飞轮力矩遵从采样控制律,基于卫星三轴姿 态的采样值、期望路径规划及卫星本体转动惯量、挠性附件太阳帆板耦合系 数构建卫星姿态机动系统的近似离散模型,具体为:
太阳同步卫星的连续模型为:
Figure BDA0002626660460000052
其中:
Figure BDA0002626660460000053
代表卫星三轴姿态角速度;φ,θ,
Figure BDA0002626660460000054
分别为卫星滚动、俯仰、偏航轴;I为太阳同步卫星刚体转动惯量;F为卫星刚体与 挠性附件太阳帆板的耦合矩阵系数;η(t)为太阳帆板的振动模态坐标向量;
Figure BDA0002626660460000055
为太阳帆板振动的模态频率矩阵;ξ为太阳帆板的振动阻尼比;u(t)为卫星飞 轮力矩;并且有
Figure BDA0002626660460000056
将(16)式中的第二个模态方程改写为:
Figure BDA0002626660460000057
代入(16)式中的第一个方程可得
Figure BDA0002626660460000058
利用采样时间为T的采样器对三轴姿态角进行如下采样:
Figure BDA0002626660460000061
与此同时,控制输入离散控制律进行如下零阶保持:
u(t)=u(kT):=u[k](21)
如此(19)式可离散为:
Figure BDA0002626660460000062
给定移位算子Δ:ω[k]→ω[k]-ω[k-1],将移位算子在式(22)等号两端映射n-1次可得:
Figure BDA0002626660460000063
其中:
Figure BDA0002626660460000064
Figure BDA0002626660460000065
考虑到模态与卫星姿态机动路径的连续性,(24)式可进一步表达为:
Θ[k]=TnΩ(n-1)[k] (26)
由于Θ[k]正比于采样时间T的n次方,当T<1时,便可忽略此项。由此可得 太阳同步卫星姿态的近似离散模型为:
Figure BDA0002626660460000066
步骤二、将步骤一中得到的近似离散模型转化为期望路径误差模型,基 于此模型设计控制器,使得太阳同步卫星飞轮根据指定信号旋转完成卫星三 轴跟踪期望路径并稳定,具体为:
假设三轴期望路径为v。从式(27)中的左右两端减去v[k+1]可得期望路径 误差模型:
Figure BDA0002626660460000071
其中
Figure BDA0002626660460000072
部分已经忽略。忽略原因为一般卫星姿态期望路径为连续可导信号,此项等于Tnv(n)[k],一般可认为足够小或等于0。
完成太阳同步卫星姿态机动任务即使卫星三轴角速度跟踪预期信号,也 就是使式(28)中的期望路径误差
Figure BDA0002626660460000073
设计的控制律为:
Figure BDA0002626660460000074
以下为太阳同步卫星姿态机动控制的验证:
假设卫星三轴机动的期望路径为:
Figure BDA0002626660460000075
这里的角速度单位为rad/s。卫星参数如下表(仿真中太阳帆板振动模态取前 五阶):
表1:太阳同步卫星系统参数
Figure BDA0002626660460000076
Figure BDA0002626660460000081
仿真中采用采样时间T=0.01(s),模型阶次n=2。卫星三轴姿态如图二所 示,图中“roll”,“pitch”,“yaw”分别代表滚动、俯仰、偏航轴。可以 看到,控制器完成了卫星姿态的机动与稳定。

Claims (1)

1.一种基于二项式近似模型的挠性卫星姿态机动控制方法,其特征在于,包括如下步骤:
步骤一、设卫星飞轮力矩遵从采样控制律,基于卫星三轴姿态的采样值、期望路径规划及卫星本体转动惯量、挠性附件耦合系数,构建卫星姿态机动系统的近似模型,具体为:
挠性卫星姿态系统的精确模型为:
Figure FDA0003057953880000011
其中,
Figure FDA0003057953880000012
代表卫星三轴姿态角速度;φ,θ,
Figure FDA0003057953880000013
分别为卫星滚动、俯仰、偏航轴;t为连续时间变量;
Figure FDA0003057953880000014
为三轴姿态角速度一阶导数;I为卫星刚体转动惯量;F为卫星刚体与挠性附件的耦合矩阵系数;η(t)为挠性附件的振动模态坐标向量;
Figure FDA0003057953880000015
Figure FDA0003057953880000016
分别为挠性附件的振动模态坐标向量的一阶与二阶导数;
Figure FDA0003057953880000017
为挠性附件的振动模态频率矩阵;ξ为挠性附件的振动阻尼比;u(t)为卫星飞轮力矩;并且有
Figure FDA0003057953880000018
将式(1)中的第二个模态方程改写为:
Figure FDA0003057953880000019
其中
Figure FDA00030579538800000110
为挠性附件振动的模态频率矩阵的平方;将式(3)代入式(1)中的第一个方程得到
Figure FDA00030579538800000111
利用采样时间为T的采样器对三轴姿态角进行如下采样:
Figure FDA0003057953880000021
其中k为自然数,代表离散时间采样点,
Figure FDA0003057953880000022
表示向下取整;与此同时,控制输入根据离散控制律进行如下零阶保持:
u(t)=u(kT)=u[k] (6)
如此式(4)得到离散为:
Figure FDA0003057953880000023
其中,ω[k],u[k],
Figure FDA0003057953880000024
η[k]分别代表三轴姿态角速度、飞轮力矩、挠性附件的振动模态坐标向量一阶导数及自身在时间k处的采样值,ω[k+1]代表下一刻k处的三轴姿态角速度采样值,ω×[k]为式(2)的k处的采样值;
给定移位算子函数Δ:ω[k]→ω[k]-ω[k-1],其中:表示将函数映射到三轴姿态角速度采样值上,ω[k-1]代表三轴姿态角速度前一步的采样值,将移位算子在式(7)等号两端映射n-1次得到:
Figure FDA0003057953880000025
其中Θ[k]为卫星姿态机动系统的近似模型的冗余项,
Figure FDA0003057953880000026
为二项式系数,A=(I-FFT)-1为反馈矩阵系数,n≥2为采样时间及卫星姿态机动系统的近似模型阶次,应为自然数,并且是可调参数,具体为:
Figure FDA0003057953880000027
其中Ω[k]为卫星自身参数项,具体为:
Figure FDA0003057953880000028
考虑到挠性附件的振动模态与卫星姿态机动系统路径的连续性,式(9)进一步表达为:
Θ[k]=TnΩ(n-1)[k] (11)
由于Θ[k]正比于采样时间T的n次方,当采样时间T<1时,忽略此项;由此得到卫星姿态机动系统的近似模型为:
Figure FDA0003057953880000031
步骤二、将步骤一中得到的卫星姿态机动系统的近似模型转化为期望路径误差模型,基于此模型设计控制器,使得卫星飞轮根据指定信号旋转完成卫星三轴跟踪期望路径并稳定,具体为:
设三轴期望路径为v(t);依据式(6)对期望路径进行采样后,从式(1)中的左右两端减去v[k+1]得到期望路径的误差模型:
Figure FDA0003057953880000032
其中
Figure FDA0003057953880000033
为期望路径误差的估计值,
e[k]=ω[k]-v[k]为期望路径误差的实际值;
式(12)中的
Figure FDA0003057953880000034
表示期望路径计算式的部分已经忽略;忽略原因为卫星姿态期望路径为连续可导信号,式(12)中的v[k+1]=Tnv(n)[k],认为足够小或等于0;
完成卫星姿态机动系统任务即使卫星三轴角速度跟踪预期信号,也就是使式(12)中的期望路径误差
Figure FDA0003057953880000035
设计的控制律为:
Figure FDA0003057953880000036
CN202010798931.6A 2020-08-11 2020-08-11 一种基于二项式近似模型的挠性卫星姿态机动控制方法 Active CN112115545B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010798931.6A CN112115545B (zh) 2020-08-11 2020-08-11 一种基于二项式近似模型的挠性卫星姿态机动控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010798931.6A CN112115545B (zh) 2020-08-11 2020-08-11 一种基于二项式近似模型的挠性卫星姿态机动控制方法

Publications (2)

Publication Number Publication Date
CN112115545A CN112115545A (zh) 2020-12-22
CN112115545B true CN112115545B (zh) 2021-08-24

Family

ID=73804789

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010798931.6A Active CN112115545B (zh) 2020-08-11 2020-08-11 一种基于二项式近似模型的挠性卫星姿态机动控制方法

Country Status (1)

Country Link
CN (1) CN112115545B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113568421B (zh) * 2021-06-17 2023-06-30 北京控制工程研究所 兼顾机动控制与稳定控制的挠性卫星姿态控制器及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777469A (zh) * 2016-11-15 2017-05-31 上海航天控制技术研究所 用于挠性卫星挠性振动评估的地面物理仿真验证系统和验证方法
CN109815916A (zh) * 2019-01-28 2019-05-28 成都蝉远科技有限公司 一种基于卷积神经网络算法的植被种植区域识别方法及系统
CN110500995A (zh) * 2019-07-12 2019-11-26 武汉大学 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法
CN111319793A (zh) * 2020-02-25 2020-06-23 上海航天控制技术研究所 敏捷卫星多种工作模式姿态基准通用设计方法
CN111332501A (zh) * 2020-02-26 2020-06-26 上海航天控制技术研究所 一种卫星姿态控制系统的多控制多基准的设计方法
CN111483615A (zh) * 2020-04-17 2020-08-04 中国科学院微小卫星创新研究院 一种量子卫星对地面光学站指向姿态实现方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8472735B2 (en) * 2010-09-30 2013-06-25 The Charles Stark Draper Laboratory, Inc. Attitude estimation with compressive sampling of starfield data
CN107479566A (zh) * 2017-08-25 2017-12-15 中国科学院长春光学精密机械与物理研究所 基于三段式路径规划的挠性卫星姿态控制方法
US10902165B2 (en) * 2018-01-09 2021-01-26 The Charles Stark Draper Laboratory, Inc. Deployable development platform for autonomous vehicle (DDPAV)

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777469A (zh) * 2016-11-15 2017-05-31 上海航天控制技术研究所 用于挠性卫星挠性振动评估的地面物理仿真验证系统和验证方法
CN109815916A (zh) * 2019-01-28 2019-05-28 成都蝉远科技有限公司 一种基于卷积神经网络算法的植被种植区域识别方法及系统
CN110500995A (zh) * 2019-07-12 2019-11-26 武汉大学 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法
CN111319793A (zh) * 2020-02-25 2020-06-23 上海航天控制技术研究所 敏捷卫星多种工作模式姿态基准通用设计方法
CN111332501A (zh) * 2020-02-26 2020-06-26 上海航天控制技术研究所 一种卫星姿态控制系统的多控制多基准的设计方法
CN111483615A (zh) * 2020-04-17 2020-08-04 中国科学院微小卫星创新研究院 一种量子卫星对地面光学站指向姿态实现方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A distributed cooperative guidance law for salvo attack of multiple anti-ship missiles;Zhang Youan,;《Chinese Journal of Aeronautics》;20150829;第28卷(第5期);第1438–1450页 *
挠性卫星的姿态机动滚动优化控制;范国伟;《红外与激光工程》;20141231;第43卷(第s期);第108-115页 *

Also Published As

Publication number Publication date
CN112115545A (zh) 2020-12-22

Similar Documents

Publication Publication Date Title
CN110794863B (zh) 一种控制性能指标可定制的重型运载火箭姿态控制方法
CN103869704B (zh) 基于扩展雅克比矩阵的空间机器人星臂协调控制方法
CN103412491B (zh) 一种挠性航天器特征轴姿态机动指数时变滑模控制方法
CN100363851C (zh) 航天器大角度机动控制的单框架力矩陀螺群的最优控制法
CN109283934B (zh) 基于旋转路径质量的航天器多约束姿态机动优化方法
CN110329546B (zh) 一种考虑引力姿轨耦合效应的小天体着陆轨迹优化方法
CN109343550B (zh) 一种基于滚动时域估计的航天器角速度的估计方法
CN105116910A (zh) 一种对地面点凝视成像的卫星姿态控制方法
CN104898642A (zh) 一种用于航天器姿态控制算法的集成测试仿真系统
CN111498147B (zh) 挠性航天器的有限时间分段滑模姿态跟踪控制算法
CN108181806B (zh) 基于采样输出的空间机器人位置与姿态自抗扰控制方法
CN108959734B (zh) 一种基于实时递推太阳光压力矩辨识方法及系统
CN105865459A (zh) 一种考虑视线角约束的小天体接近段制导方法
CN112115545B (zh) 一种基于二项式近似模型的挠性卫星姿态机动控制方法
CN111459188A (zh) 一种基于四元数的多旋翼非线性飞行控制方法
CN113268059A (zh) 一种基于有限时间扩张状态观测器的多无人艇编队控制方法
CN109901402B (zh) 一种基于航向平滑技术的自主水下机器人路径跟踪方法
CN114138010A (zh) 一种基于加权偏差的多智能体高阶滑模编队控制方法
CN110488854B (zh) 一种基于神经网络估计的刚性飞行器固定时间姿态跟踪控制方法
CN103412571A (zh) 一种基于反馈线性化技术的航天器相对姿态控制方法
Lasemi et al. Spacecraft attitude control: Application of fine trajectory linearization control
Ma et al. Deep learning-based nonlinear model predictive control of the attitude manoeuvre of a barbell electric sail through voltage regulation
McCann et al. Rigid body pose estimation on TSE (3) for spacecraft with unknown moments of inertia
CN103438892B (zh) 一种改进的基于ekf的天文自主定轨算法
Aucoin Real-time optimal slew maneuver planning for small satellites

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