CN106326576B - 一种任意基准系下的整星偏置角动量的偏航估计方法 - Google Patents
一种任意基准系下的整星偏置角动量的偏航估计方法 Download PDFInfo
- Publication number
- CN106326576B CN106326576B CN201610741694.3A CN201610741694A CN106326576B CN 106326576 B CN106326576 B CN 106326576B CN 201610741694 A CN201610741694 A CN 201610741694A CN 106326576 B CN106326576 B CN 106326576B
- Authority
- CN
- China
- Prior art keywords
- yaw
- angular momentum
- synthesis
- angle
- satellite
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Abstract
一种任意基准系下的整星偏置角动量的偏航估计方法,在滚动角或俯仰角测量由敏感器给出,偏航角没有敏感器测量的情况下,仅通过偏航估计的方法来估算偏航角。在此基础上,整星通过角动量喷气管理的方法来实现基于任意基准系下的卫星三轴姿态稳定控制。本发明方法可根据用户需要或用滚动角测量值或用俯仰角测量作为输入值,具有很强的适应性,即仅通过调整初始的参数设计值,就可以获得需要的偏航估计结果,便于偏航估计的通用规范化方法编制。
Description
技术领域
本发明属于卫星控制技术领域,涉及一种卫星姿态角的预估方法。
背景技术
卫星姿态预估是卫星控制中的常见问题。
现有的卫星偏航估计方法多是基于轨道坐标系的,主要是因为在轨卫星多以轨道坐标系为基准。但是在研型号和背景预研阶段的卫星,则根据有效载荷的任务需求不同,作为本体坐标系参考的基准坐标系也会有所不同。例如,卫星为了适应天线对国土形状的覆盖,要求卫星以东南坐标系为参考基准坐标系,由于东南坐标系与轨道坐标系仅相差一个周期性变化的偏航角,现有的操作方法是在基于轨道坐标系的偏航估计方法中再增加一部分标称偏航角的修正量来实现。然而,未来卫星型号的参考基准坐标系相对轨道坐标系将变得更为复杂,或是要求在不同任务阶段采用不同的参考坐标系,那么现有的基于轨道坐标系的偏航估计的修补方法将不再适用。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种通用性的、可以在任意基准坐标系下预估卫星偏航角的方法。
本发明的技术解决方案是:一种任意基准系下的整星偏置角动量的偏航估计方法,包括如下步骤:
(1)在卫星的滚动轴、俯仰轴和/或偏航轴上分别安装动量轮,同时在卫星的滚动轴和俯仰轴上安装姿态敏感器;
(2)利用所述姿态敏感器,获取卫星的滚动角测量值和俯仰角测量值;
(3)利用滚动角测量值构建如下的偏航估计方程,
观测系数l1,l2满足特征方程的根为负数;
或者利用俯仰角测量值构建如下的偏航估计方程,
观测系数l1,l2满足特征方程的根为负数;
其中e1和e2为偏航角和滚动角的误差状态变量,ωrx,ωry,ωrz为基准系相对惯性系的三个方向的角速度在基准系下的描述;hx,hy,hz分别为由动量轮组在卫星三轴方向产生的角动量,ΔhxΔhy和Δhz分别为三轴角动量偏差值,θb分别为滚动角和俯仰角的偏置值,Tdx,Tdy,Tdz为作用在星体三个方向上的干扰力矩分量,Tdx_com,Tdz_com为滚动方向和偏航方向干扰力矩补偿值,分别为俯仰角估计值、滚动角估计值和偏航角估计值;
(4)采用步骤(3)中建立的偏航估计方程进行偏航角的估计。
所述步骤(4)中,为了进一步提高偏航角估计的精度,在利用建立的偏航估计方程进行偏航角的估计时,还对卫星上滚动轴、俯仰轴和偏航轴上安装的动量轮分别进行角动量卸载,其中:
利用滚动角测量值构建估计方法的卸载方式如下:
俯仰方向卸载槛值δy满足-δy≤h_synthesis_b(2)-hy≤δy;
滚动方向卸载槛值为δx=h_synthesis_i(1)-h0_unload_i(1);
偏航方向卸载槛值为δz=h_synthesis_i(3)-h0_unload_i(3);
利用俯仰角测量值构建估计方法的卸载方式如下:
滚动方向卸载槛值δx满足-δx≤h_synthesis_b(1)-hx≤δx;
俯仰方向卸载槛值为δy=h_synthesis_i(2)-h0_unload_i(2)
偏航方向卸载槛值为δz=h_synthesis_i(3)-h0_unload_i(3)
其中,为星本体角动量设计值在惯性系下的描述,为本体角动量在惯性系下的描述,h_synthesis_i=Cib*h_synthesis_b,Cib为本体坐标系相对于惯性坐标系的姿态矩阵。
本发明与现有技术相比的优点在于:在卫星偏航方向缺少测量的情况下,现有技术是在轨道系下基于整星角动量垂直于轨道面的前提来设计偏航估计器的,本发明方法则扩展了上述前提,一方面整星角动量可以不垂直于轨道面,另一方面,在任意坐标系下都可以进行偏航估计器的设计。另外,本发明还在降低角动量偏差累积效应对偏航估计精度的影响方面,提出了角动量卸载管理方法。现有技术是将星体角动量在轨道系下进行管理,而本发明在设计模型构建时基于的是任意参考坐标系,给出的是惯性系相对任意参考系的相对角速度描述,在设计偏航估计方程中考虑了设计模型中由相对角速度带来的偏差部分,使得估计误差在没有外力矩干扰和残留角动量的影响下经过过渡过程时间后接近零值,即偏航估计值接近真值,因此本发明方法具有更强的通用性。
附图说明
图1为本发明方法的流程框图;
图2为本发明实施例中卫星相对于东南系的偏航估计姿态角;
图3为本发明实施例中卫星相对于东南系的滚动和俯仰测量角;
图4为本发明实施例中卫星相对于轨道系的三轴姿态角;
图5为本发明实施例中惯性系下的星本体角动量设计结果;
图6为本发明实施例中惯性系下的星本体角动量实测结果。
具体实施方式
为了适应未来更广阔更复杂的卫星任务需要,以及从算法的规范化、软件的维护性需求出发,本发明提出了一种适用于任意参考基准坐标系的偏航估计方法,使用者仅需给定当前本体坐标系需要参考的基准坐标系即可。
偏航估计方法主要是基于卫星长周期的运动控制模型,当有三个主要条件成立时才能达到估计效果:其一,基准坐标系相对于惯性坐标系要具备一定的角速度;其二,在此相对角速度方向存在一个可设计的标称角动量;其三,星上喷气等设备要具备角动量管理能力。
本发明方法还能够根据输入的基准坐标系、标称角动量、动量轮配置情况、太阳光压等干扰力矩补偿残余量来判断是否具备偏航估计的条件以及在条件具备下的估计精度。
如图1所示,为本发明方法的流程框图,采用本发明方法,要求卫星上安装有测量滚动方向或者俯仰方向的地球敏感器或太阳敏感器或其它敏感器,动量轮组和喷气装置。本发明方法的主要步骤如下:
(1)整星的动力学方程可以描述为:
其中,H=[Hx Hy Hz]T为整星的三轴角动量(角标x,y,z分别表示俯仰轴、滚动轴和偏航轴),即星体内动量轮组产生的角动量和星体自身产生的角动量之和,可以写成如下形式:
H=[Hx Hy Hz]T=[hx+Ixωx hy+Iyωy hz+Izωz]T (2)
其中,hx,hy,hz分别为由动量轮组在卫星三轴方向产生的角动量,可计算求得,Ixωx,Iyωx,Izωx为星体自身产生的角动量,ω×表示向量ω=[ωx,ωy,ωz]T的反对称矩阵,ω为星体相对惯性坐标系的角速度,可表示为:
其中,为卫星的滚动角、俯仰角和偏航角,为基准坐标系相对惯性坐标系的三个方向的角速度在基准坐标系下的描述。Cbr为卫星本体坐标系相对于基准坐标系的姿态矩阵,对于3-1-2转序下Cbr的数学描述如下:
矩阵A为欧拉角的导数与ω的关系,对于3-1-2转序下A的数学描述如下:
其中,
当卫星本体坐标系与基准坐标系之间的夹角为小角度时,A为单位矩阵。式(3)可写成如下形式:
为的反对称矩阵;
将式(2)、(3)、(4)代入式(1)中,且忽略二阶小量后,经整理得
其中,h=[hx hy hz]T,Td=[Tdx,Tdy,Tdz]T为外界干扰力矩,如太阳光压、气动力矩等,Tdx,Tdy,Tdz为作用在星体三轴方向上的干扰力矩分量。
本发明估计方法的输入可以选择采用滚动角测量,也可以选择采用俯仰角测量。
在不考虑角动量累积的情况下,若采用滚动角测量,则设置
h=[0;hy;0]T (6a)
若采用俯仰角测量,则设置
h=[hx;0;0]T (6b)
其中,hx和hy分别为动量轮组产生的俯仰方向角动量偏置设计值和滚动方向角动量偏置设计值。
但在实际工作过程中,受到外界干扰,星体和动量轮之间通过角动量交换,使得动量轮吸收了外界的干扰角动量,三个方向的角动量就会偏离设计值(6a)或(6b),则实际三个方向的角动量可以写成如下形式:
对应于式(6a),实际状况下的三轴角动量为
hreal=[Δhx hy+ΔhyΔhz]T (7a)
对应于式(6b),实际状况下的三轴角动量为
hreal=[hx+ΔhxΔhyΔhz]T (7b)
ΔhxΔhy和Δhz分别为卫星三轴角动量偏差值。
(2)若采用滚动角测量,采用理论设计值(6a),则方程(5)可以简化为
其中,为滚动角测量值。方程(8)中仅考虑偏航和滚动的耦合设计方程部分,俯仰方向有独立的测量值和控制量,因此不予在方程中给出。
若采用(7a),则实际对象表述为
进入步骤(3);
若采用俯仰角测量,采用理论设计值(6b),则待设计的对象为
其中,θm为俯仰角测量值。方程(10)中仅考虑偏航和俯仰的耦合设计方程部分,滚动方向有独立的测量值和控制量,因此不予在方程中给出。
若采用(7b),则实际对象表述为
进入步骤(4);
(3)根据式(8)设计的偏航估计方程为
其中,分别为滚动角估计值和偏航角估计值。
利用式(9)和(12),则有
其中,
待设计的观测系数l1,l2满足以下特征方程根为负数:
θb为俯仰角偏置值,该偏置值为控制任务需要预期设置值,Tdx_com,Tdz_com为滚动方向和偏航方向干扰力矩补偿值。偏航估计的稳态精度为
进入步骤(5);
(4)根据式(10)设计的偏航估计方程为
利用式(11)和(16),则有
其中,为俯仰角估计值,
待设计的观测系数l1,l2满足以下特征方程根为负数:
为滚动角偏置值,该偏置值为控制任务需要预期设置值,Tdy_com,Tdz_com为俯仰方向和偏航方向干扰力矩补偿值。偏航估计的稳态精度为
进入步骤(5);
(5)角动量卸载。由式子(15)和式子(18)可知,当星体姿态接近目标姿态且干扰力矩被补偿后,影响偏航估计误差的因素主要为ΔhxΔhyΔhz,即
由式(15)得到的偏航误差为:
由式(18)得到的偏航误差为:
若能消除ΔhxΔhyΔhz,则偏航估计值将趋近于真值。因此,需要进行角动量卸载,使得ΔhxΔhyΔhz维持在给定的范围内,达到偏航估计精度要求。
动量轮转速合成的三轴角动量在卫星本体坐标系下的描述记为h_synthesis_b,则其在惯性坐标系下的描述为
h_synthesis_i=Cib*h_synthesis_b
其中,Cib为本体坐标系相对于惯性坐标系的姿态矩阵,h_synthesis_i为卫星本体角动量在惯性坐标系下的描述。h_synthesis_b由动量轮转速合成方法如下:
h_synthesis_b=KmCMW*[rMW1,rMW2,…,rMWn]T+Iω (21)
其中,CMW为动量轮安装矩阵,rMWi,i=1,2,…,n为动量轮i的转速,Km为转速与角动量常值转换系数,根据不同类型的动量轮而不同。I=diag[Ix,Iy,Iz]为星体三轴转动惯量,ω=[ωx,ωy,ωz]T为陀螺测量的星体三轴角速度。
若采用滚动角测量,角动量设计值在卫星本体坐标系下的描述为h=[0;hy;0]T,则有
h0_unload_i=Cib*[0,h_synthesis_b(2),0]T (22)
其中,h0_unload_i为星本体角动量设计值在惯性坐标系下的描述,由于干扰角动量对卫星影响分析一般是以一个轨道周期为基础的,因此动量轮卸载基线也应该按照一个轨道周期来计算。这里,-δy≤h_synthesis_b(2)-hy≤δy,δy为本体俯仰轴方向的角动量卸载槛值。滚动轴和偏航轴的角动量卸载槛值分别为:
δx=h_synthesis_i(1)-h0_unload_i(1)
δz=h_synthesis_i(3)-h0_unload_i(3)
其中,h_synthesis_i(k),k=1,2,3表示向量的第k个元素,以下类同。
若采用俯仰角测量,角动量设计值在本体系下描述为h=[hx;0;0]T,则有
h0_unload_i=Cib*[h_synthesis_b(1),0,0]T (23)
其中,-δx≤h_synthesis_b(1)-hx≤δx,δx为本体滚动轴方向的角动量卸载槛值。俯仰轴和偏航轴的角动量卸载槛值分别为:
δy=h_synthesis_i(2)-h0_unload_i(2)
δz=h_synthesis_i(3)-h0_unload_i(3)
角动量卸载本质就是将由干扰力矩引起的偏离标称值的那部分角动量通过其他执行机构给消除掉,例如,喷气推力器、磁力矩器等。
实施例
此处选用卫星的基准坐标系为东南坐标系,偏航估计基于东南坐标系设计。
东南坐标系OXdYdZd定义:O为卫星质心,Zd轴指向地心,Xd轴指向卫星运行方向,且平行于赤道平面,Yn轴与Xn、Zn轴右手正交。
卫星运行在地球同步轨道上,轨道倾角为5.5度。
初始时间(北京时):2015年8月6日16时40分33秒
卫星位置(km):-16504.192625 -38628.441000 -3752.763227
卫星速度(km/s):2.821086 -1.215998 0.110196
东南系相对于惯性系的角速度描述为
其中,ω0为轨道角速度,约为-0.00416度/秒;ψl为轨道系相对于东南系的偏航角,为日周期正余弦变化,最大幅值为5.5度。
采用滚动角测量输入时,其角动量设计值为
将(24)和(25)代入下式中
经整理后得,
偏航估计器为
Tdx_com和Tdz_com为滚动轴和偏航轴方向的太阳光压力矩补偿分量。
仿真时间长度为10天,仿真结果见图2~图6。
图2为卫星相对于东南系的偏航估计姿态角。图中,横坐标为时间,纵坐标为偏航估计姿态角,该估计值为本方法估计所得。图3为卫星相对于东南系的滚动和俯仰测量角。图中,横坐标为时间,纵坐标为分别为滚动角测量值和俯仰角测量值。图4为卫星相对于轨道系的三轴姿态角,图中,横坐标为时间,纵坐标分别为滚动角、俯仰角和偏航角,这里给出该图的目的是为了能够清楚的表明偏航估计效果,即偏航估计在基于东南系下设计,在东南系的偏航角近零,而东南系相对于轨道系主要差别在于偏航角的日周期性变化,其幅值等于轨道倾角,从图中可以看出在轨道系下,偏航角呈现周期性变化,幅值为5.5。图5为在惯性系下的星本体角动量设计结果。图中,横坐标为时间,纵坐标分别为惯性系下的星本体滚动轴角动量、俯仰轴角动量和偏航轴角动量,星体俯仰方向角动量受到外界干扰力矩的作用呈现周期性变化,那么该角动量在惯性系下的描述为图5所示。图6在惯性系下的星本体角动量实测结果,图中,横坐标为时间,纵坐标分别为惯性系下星体实际角动量结果。从图中可以看出,当超过给定界限7Nms后,通过喷气方式对星体角动量进行卸载。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (2)
1.一种任意基准系下的整星偏置角动量的偏航估计方法,其特征在于包括如下步骤:
(1)在卫星的滚动轴、俯仰轴和/或偏航轴上分别安装动量轮,同时在卫星的滚动轴和俯仰轴上安装姿态敏感器;
(2)利用所述姿态敏感器,获取卫星的滚动角测量值和俯仰角测量值θm;
(3)利用滚动角测量值构建如下的偏航估计方程,
观测系数l1,l2满足特征方程的根为负数;
或者利用俯仰角测量值构建如下的偏航估计方程,
观测系数l1,l2满足特征方程的根为负数;
其中ωrx,ωry,ωrz为基准系相对惯性系的三个方向的角速度在基准系下的描述;hx,hy,hz分别为由动量轮组在卫星三轴方向产生的角动量,θb分别为滚动角和俯仰角的偏置值,Tdx,Tdy,Tdz为作用在星体三个方向上的干扰力矩分量,Tdx_com,Tdz_com为滚动方向和偏航方向干扰力矩补偿值,分别为俯仰角估计值、滚动角估计值和偏航角估计值;
(4)采用步骤(3)中建立的偏航估计方程进行偏航角的估计。
2.根据权利要求1所述的一种任意基准系下的整星偏置角动量的偏航估计方法,其特征在于:所述步骤(4)中,在利用建立的偏航估计方程进行偏航角的估计时,还对卫星上滚动轴、俯仰轴和偏航轴上安装的动量轮分别进行角动量卸载,其中:
利用滚动角测量值构建估计方法的卸载方式如下:
俯仰方向卸载槛值δy满足-δy≤h_synthesis_b(2)-hy≤δy;
滚动方向卸载槛值为δx=h_synthesis_i(1)-h0_unload_i(1);
偏航方向卸载槛值为δz=h_synthesis_i(3)-h0_unload_i(3);
利用俯仰角测量值构建估计方法的卸载方式如下:
滚动方向卸载槛值δx满足-δx≤h_synthesis_b(1)-hx≤δx;
俯仰方向卸载槛值为δy=h_synthesis_i(2)-h0_unload_i(2)
偏航方向卸载槛值为δz=h_synthesis_i(3)-h0_unload_i(3)
其中,为星本体角动量设计值在惯性系下的描述,为本体角动量在惯性系下的描述,h_synthesis_i=Cib*h_synthesis_b,Cib为本体坐标系相对于惯性坐标系的姿态矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610741694.3A CN106326576B (zh) | 2016-08-26 | 2016-08-26 | 一种任意基准系下的整星偏置角动量的偏航估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610741694.3A CN106326576B (zh) | 2016-08-26 | 2016-08-26 | 一种任意基准系下的整星偏置角动量的偏航估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106326576A CN106326576A (zh) | 2017-01-11 |
CN106326576B true CN106326576B (zh) | 2019-07-12 |
Family
ID=57791198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610741694.3A Active CN106326576B (zh) | 2016-08-26 | 2016-08-26 | 一种任意基准系下的整星偏置角动量的偏航估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106326576B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107499536B (zh) * | 2017-08-17 | 2020-04-07 | 上海航天控制技术研究所 | 一种偏置动量卫星高精度跟踪控制偏流角的方法 |
CN108959734B (zh) * | 2018-06-15 | 2023-08-29 | 北京控制工程研究所 | 一种基于实时递推太阳光压力矩辨识方法及系统 |
CN110688731B (zh) * | 2019-08-26 | 2020-11-20 | 北京控制工程研究所 | 一种并联式指向平台的扰动建模与抑制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1174982A (zh) * | 1996-01-12 | 1998-03-04 | 环球星有限合伙人公司 | 用于轨道偏航控制的动态偏置 |
CN101852605A (zh) * | 2010-06-10 | 2010-10-06 | 南京航空航天大学 | 基于简化自适应滤波的磁测微小卫星姿态确定方法 |
CN103153790A (zh) * | 2010-10-01 | 2013-06-12 | 希尔克瑞斯特实验室公司 | 使用运动传感器和附接至装置的磁力计的测量数据估计该装置在重力参照系中的偏航角的设备和方法 |
CN104020778A (zh) * | 2014-06-18 | 2014-09-03 | 哈尔滨工业大学 | 基于跟踪时间-能耗最优轨线的挠性卫星姿态机动控制方法 |
CN105116910A (zh) * | 2015-09-21 | 2015-12-02 | 中国人民解放军国防科学技术大学 | 一种对地面点凝视成像的卫星姿态控制方法 |
-
2016
- 2016-08-26 CN CN201610741694.3A patent/CN106326576B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1174982A (zh) * | 1996-01-12 | 1998-03-04 | 环球星有限合伙人公司 | 用于轨道偏航控制的动态偏置 |
CN101852605A (zh) * | 2010-06-10 | 2010-10-06 | 南京航空航天大学 | 基于简化自适应滤波的磁测微小卫星姿态确定方法 |
CN103153790A (zh) * | 2010-10-01 | 2013-06-12 | 希尔克瑞斯特实验室公司 | 使用运动传感器和附接至装置的磁力计的测量数据估计该装置在重力参照系中的偏航角的设备和方法 |
CN104020778A (zh) * | 2014-06-18 | 2014-09-03 | 哈尔滨工业大学 | 基于跟踪时间-能耗最优轨线的挠性卫星姿态机动控制方法 |
CN105116910A (zh) * | 2015-09-21 | 2015-12-02 | 中国人民解放军国防科学技术大学 | 一种对地面点凝视成像的卫星姿态控制方法 |
Non-Patent Citations (1)
Title |
---|
三轴稳定卫星姿态控制系统研究;褚庆军;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》;20061215(第12期);第4章 |
Also Published As
Publication number | Publication date |
---|---|
CN106326576A (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106873611B (zh) | 一种多通道线性自抗扰控制器的设计方法 | |
CN110579740B (zh) | 一种基于自适应联邦卡尔曼滤波的无人船组合导航方法 | |
CN105700536B (zh) | 基于绳系拖曳系统的主动星姿态和系绳摆振联合控制方法 | |
CN106915477B (zh) | 一种姿态控制方法 | |
CN104792340B (zh) | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 | |
CN102608596B (zh) | 一种用于机载惯性/多普勒雷达组合导航系统的信息融合方法 | |
CN107450588B (zh) | 一种挠性航天器姿态控制系统的强抗扰控制方法 | |
CN102692225B (zh) | 一种用于低成本小型无人机的姿态航向参考系统 | |
CN101726295B (zh) | 考虑加速度补偿和基于无迹卡尔曼滤波的惯性位姿跟踪方法 | |
CN103822633B (zh) | 一种基于二阶量测更新的低成本姿态估计方法 | |
CN108227485B (zh) | 一种带有执行器饱和的空间机器人自抗扰控制方法 | |
CN103175530B (zh) | 一种航空遥感惯性稳定平台耦合力矩估计与补偿方法 | |
CN103092209B (zh) | 一种基于动量轮控制的全姿态捕获方法 | |
CN102519470B (zh) | 多级嵌入式组合导航系统及导航方法 | |
CN105629732B (zh) | 一种考虑控制受限的航天器姿态输出反馈跟踪控制方法 | |
CN103112603B (zh) | 欠驱动高速自旋卫星建立正常姿态的方法 | |
CN108181916B (zh) | 小卫星相对姿态的控制方法及装置 | |
CN104483973A (zh) | 基于滑模观测器的低轨挠性卫星姿态跟踪控制方法 | |
CN110162855B (zh) | 遥感卫星星上旋转载荷动态精度分析及误差分配方法 | |
CN105929836B (zh) | 用于四旋翼飞行器的控制方法 | |
CN103488081B (zh) | 惯性稳定平台控制方法 | |
CN106326576B (zh) | 一种任意基准系下的整星偏置角动量的偏航估计方法 | |
CN103955224B (zh) | 一种用于相对运动视线跟踪的姿态控制方法 | |
CN108959734B (zh) | 一种基于实时递推太阳光压力矩辨识方法及系统 | |
CN106052682A (zh) | 一种混合式惯性导航系统及导航方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |