CN108089434A - 一种基于磁强计的皮纳卫星姿态捕获方法 - Google Patents
一种基于磁强计的皮纳卫星姿态捕获方法 Download PDFInfo
- Publication number
- CN108089434A CN108089434A CN201711311392.3A CN201711311392A CN108089434A CN 108089434 A CN108089434 A CN 108089434A CN 201711311392 A CN201711311392 A CN 201711311392A CN 108089434 A CN108089434 A CN 108089434A
- Authority
- CN
- China
- Prior art keywords
- satellite
- attitude
- mrow
- msub
- under
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 16
- 230000005358 geomagnetic field Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 3
- 238000013213 extrapolation Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 abstract description 13
- 230000005484 gravity Effects 0.000 abstract description 2
- 230000036544 posture Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于磁强计的皮纳卫星姿态捕获方法,属于卫星姿态控制技术领域。该方法在卫星失去姿态基准的情况下,以卫星姿态四元数和陀螺常值漂移为状态量、以磁强计测量值为观测量,利用递推的卡尔曼滤波进行姿态确定,首先利用磁力矩器进行速率阻尼,待速率阻尼完成后,该方法已经获取到了较为准确的姿态信息,以此为输入信息利用磁力矩器进行卫星的三轴姿态控制,从而建立卫星对地定向的运行姿态。本发明方法不受太阳信息的限制,也没有姿态条件的约束,具有很强的工程可操作和可实现性。
Description
技术领域
本发明涉及一种卫星姿态捕获的方法,特别是一种基于磁强计的皮纳卫星姿态捕获方法,属于卫星姿态控制技术领域。
背景技术
当卫星在轨失去姿态基准的情况下,卫星往往处于角速度较大、姿态角未知的状态,需要阻尼角速度、获取姿态信息并消除姿态、重新建立对地运行姿态。目前的皮纳卫星在姿态捕获过程中,一般都将速率阻尼与姿态捕获这两大任务的割裂开来,首先进行速率阻尼,在速率阻尼完成后,利用磁强计与太阳敏感器进行双矢量姿态确定从而获取姿态基准进行姿态捕获。这种方法依赖太阳信息,受轨道设计及太阳敏感器安装的限制,在速率阻尼完成后不一定能够确保太阳敏感器见太阳,从而无法实施双矢量姿态确定、无法完成捕获姿态。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种基于磁强计的皮纳卫星姿态捕获方法,不依赖太阳信息,不受轨道设计及太阳敏感器安装的限制,仅依靠磁强计即可在速率阻尼过程中获取卫星姿态,解决了现有技术中受太阳信息约束的问题。
本发明的技术解决方案是:一种基于磁强计的皮纳卫星姿态捕获方法,该方法每个控制周期执行如下步骤:
(1)、以磁强计的测量值为观测量,以惯性坐标系下的卫星姿态四元数和陀螺的常值漂移为状态量,采用递推的卡尔曼滤波方法计算得到惯性坐标系下的卫星姿态四元数估计值和卫星姿态角速度ω信息,进入步骤(2);
(2)、判断惯性坐标系下的卫星的姿态角速度ω幅值|ω|是否小于预设的姿态角速度门限,如果小于,则进入步骤(3);否则,计算控制磁矩M驱动磁力矩器,以消除卫星姿态角速度,在下一个周期到来时,从步骤(1)开始执行;
(3)、将惯性坐标系下的卫星姿态四元数估计值转换为轨道坐标系下的卫星姿态四元数估计qBO,将惯性坐标系下的卫星姿态角速度ω转换成轨道系坐标系下的卫星姿态角速度ωBO,并根据所述轨道系坐标系下的卫星姿态四元数估计值qBO和卫星姿态角速度ωBO计算控制磁矩M驱动磁力矩器,以消除卫星的姿态角,进入步骤(4);
(4)、判断卫星姿态角Φ的幅值|Φ|是否小于预设的姿态角门限,如果小于,则认为捕获成功,结束捕获过程;否则,在下一个周期到来时,从步骤(1)开始执行。
所述步骤(1)计算惯性坐标系下的卫星姿态四元数估计值和卫星的姿态角速度ω信息的具体实现步骤为:
(1.1)、根据控制周期Δt内陀螺输出卫星三轴角度增量Δg和陀螺的常值漂移估计值bk-1/k-1计算卫星的姿态角速度ω:
式中,b的初值为陀螺的常值漂移标称值或者为[0,0,0];
(1.2)、根据上一控制周期的惯性坐标系下的卫星姿态四元数估计值qk-1/k-1外推当前控制周期惯性坐标系下的卫星姿态四元数预估值qk/k-1:
式中,q=[q1,q2,q3,q4],其中q4为标量,q的初始值为q0/0=[0,0,0,1];
(1.3)、计算滤波修正矩阵K:
Pk/k-1=A·Pk-1/k-1·AT+Q
K=Pk/k-1·CT·(C·Pk/k-1·CT+R)-1
式中,P是滤波误差方差阵,Q是状态误差方差阵,与状态方程的随机误差特性相关,R是测量误差方差阵,与测量方程的随机误差特性相关,A是状态转移矩阵,取值为
ωx,ωy,ωz分别为卫星的姿态角速度ω在卫星本体坐标系下的XYZ三轴分量,
ω=[ωx,ωy,ωz],C是观测矩阵,取值为Bx、By、Bz为磁强计测量的地磁场B在卫星本体坐标系下的XYZ三轴的分量,B=[Bx,By,Bz];
(1.4)、根据当前控制周期磁强计测量的地磁场强度Bk与上一个控制周期磁强计测量的地磁场强度Bk-1之差,对惯性坐标系下的卫星姿态四元数预估值qk/k-1进行修正,获取当前控制周期的惯性坐标系下的卫星姿态四元数估值qk/k和陀螺常漂估值bk/k:
δX=K·(Bk-Bk-1);
δq=[δX(0)δX(1)δX(2)]T;
δb=[δX(3)δX(4)δX(5)]T;
qk/k=qk/k-1+E(qk/k-1)·K·δq;
bk/k=bk-1/k-1+K·δb。
所述步骤(2)采用B-dot方法计算控制磁矩M驱动磁力矩器。
所述步骤(3)采用比例+微分方法计算控制磁矩M驱动磁力矩器。
所述预设的姿态角门限取值范围[5°,15°]。
所述预设的姿态角速度门限取值范围[0.1°/s,1°/s]。
本发明与现有技术相比的优点在于:
(1)、本发明以磁强计的测量值确定卫星姿态,不依赖太阳信息,解决了现有技术受轨道设计及太阳敏感器安装的限制而无法获知卫星姿态信息的问题;
(2)、本发明在速率阻尼的同时进行确定姿态运算,当速率阻尼完成的同时,姿态滤波结果已经收敛到真值附近,可直接进行控制消除姿态角的姿态捕获过程,缩短捕获时间;
(3),本发明在速率阻尼阶段采用B-dot方法计算控制磁矩M驱动磁力矩器,简化控制率的设计;
(4)、本发明在姿态捕获阶段采用比例+微分方法计算控制磁矩M驱动磁力矩器,算法稳定可靠、鲁棒性强。
附图说明
图1为本发明方法的流程框图;
图2为本发明实施例中的姿态估计误差示意图;
图3为本发明实施例中的陀螺常值漂移估计误差示意图。
图4为本发明实施例中的对地三轴姿态示意图。
具体实施方式
以下结合附图和具体实施例对本发明进行详细描述。
如图1所示,本发明提供了一种基于磁强计的皮纳卫星姿态捕获方法,该方法每个控制周期执行如下步骤:
(1)、以磁强计的测量值为观测量,以惯性坐标系下的卫星姿态四元数估计值和陀螺的常值漂移估计值为状态量,采用递推的卡尔曼滤波方法计算得到惯性坐标系下的卫星姿态四元数估计值和卫星姿态角速度ω信息,进入步骤(2);
具体实现步骤为:
(1.1)、根据控制周期Δt内陀螺输出卫星三轴角度增量Δg和陀螺的常值漂移估计值bk-1/k-1计算卫星的姿态角速度ω:
式中,b的初值为陀螺的常值漂移标称值或者为[0,0,0];
(1.2)、根据上一控制周期的惯性坐标系下的卫星姿态四元数估计值qk-1/k-1外推当前控制周期惯性坐标系下的卫星姿态四元数预估值qk/k-1:
式中,q=[q1,q2,q3,q4],其中q4为标量,q的
初始值为q0/0=[0,0,0,1];
(1.3)、计算滤波修正矩阵K:
Pk/k-1=A·Pk-1/k-1·AT+Q
K=Pk/k-1·CT·(C·Pk/k-1·CT+R)-1
式中,P是滤波误差方差阵,Q是状态误差方差阵,与状态方程的随机误差特性相关,R是测量误差方差阵,与测量方程的随机误差特性相关,A是状态转移矩阵,取值为ωx,ωy,ωz分别为卫星的姿态角速度ω在卫星本体坐标系下的XYZ三轴分量,ω=[ωx,ωy,ωz],C是观测矩阵,取值为Bx、By、Bz为磁强计测量的地磁场B在卫星本体坐标系下的XYZ三轴的分量,B=[Bx,By,Bz];
(1.4)、根据当前控制周期磁强计测量的地磁场强度Bk与上一个控制周期磁强计测量的地磁场强度Bk-1之差,对惯性坐标系下的卫星姿态四元数预估值qk/k-1进行修正,获取当前控制周期的惯性坐标系下的卫星姿态四元数估值qk/k和陀螺常漂估值bk/k:
δX=K·(Bk-Bk-1);
δq=[δX(0)δX(1)δX(2)]T;
δb=[δX(3)δX(4)δX(5)]T;
qk/k=qk/k-1+E(qk/k-1)·K·δq;
bk/k=bk-1/k-1+K·δb。
上述符号中,不加下标或下标k/k表示本控制周期的值,下标k-1/k-1表示上控制周期的值,下标k/k-1表示上控制周期对本控制周期的预测,下标为x/y/z则表示在卫星x/y/z轴的分量;
本控制周期滤波计算获得的四元数估值qk/k和陀螺常漂估值bk/k,在下一周期计算时,作为上一时刻的qk-1/k-1和bk-1/k-1,再次进行上述计算。
上述计算公式的原理来自于递推卡尔曼滤波算法,状态量为惯性坐标系下的卫星姿态四元数估计值qk/k和陀螺的常值漂移估计值bk/k,观测量为磁强计测量值Bk。
(2)、判断惯性坐标系下的卫星的姿态角速度ω幅值|ω|是否小于预设的姿态角速度门限,如果小于,则进入姿态捕获阶段执行步骤(3);否则,计算控制磁矩M驱动磁力矩器,以消除卫星姿态角速度,在下一个周期到来时,从步骤(1)开始执行,所述预设的姿态角速度门限取值范围[0.1°/s,1°/s]。
姿态控制任务则分为消除角速度的速率阻尼阶段和消除姿态角的姿态捕获阶段,首先进行速率阻尼阶段的控制,待阻尼完成后再转入姿态捕获阶段。这个步骤属于速率阻尼阶段。在速率阻尼阶段,控制目的是消除卫星的角速度。
控制率可以采用B-dot方法,磁力矩器的控制磁矩M计算公式为:
式中Km是磁阻尼系数,与磁力矩器的最大磁矩有关,磁力矩器的最大磁矩越大,磁阻尼系数越大,例如,采用最大磁矩为1Am2的磁力矩器,取值范围为10~50。
上述公式的原理是:在较短的时间间隔Δt内,地磁场的变化非常小,可认为是个常值,则磁强计测量到的地磁场强度B的变化是由卫星转动引起的,因此B的变化信息蕴含了卫星的角速度ω信息,其次,磁力矩器产生的控制磁矩需要与地磁场发生关系才能产生控制力矩,三者为叉乘关系,综合上述即可获取控制磁矩。
(3)、将惯性坐标系下的卫星姿态四元数估计值转换为轨道坐标系下的卫星姿态四元数估计qBO,将惯性坐标系下的卫星姿态角速度ω转换成轨道系坐标系下的卫星姿态角速度ωBO,并根据所述轨道系坐标系下的卫星姿态四元数估计值qBO和卫星姿态角速度ωBO,采用比例+微分方法计算控制磁矩M驱动磁力矩器,以消除卫星的姿态角,进入步骤(4);
在姿态捕获阶段,控制目的是消除卫星的姿态角,控制率可以采用比例+微分方法,比例+微分方法计算磁力矩器的控制磁矩M的公式为:
式中,J是卫星转动惯量阵,Kp是比例控制系数,Kd是微分控制系数,qBO=[qBO1,qBO2,qBO3,qBO4]是卫星在轨道系下的四元数(其中qBO4是标量),由惯性坐标系下四元数估计值q结合轨道系相对于惯性系的转换矩阵COI计算得到,ωBO是卫星相对于轨道系的角速度,由惯性坐标系下的卫星姿态角速度ω结合轨道角速度ω0计算得到。Kp是比例控制系数,Kd是微分控制系数,都为正数,Kp和Kd的值与磁力矩器的最大磁矩有关,磁力矩器的最大磁矩越大,控制系数越大,例如,采用最大磁矩为1Am2的磁力矩器,Kp取值范围为0.2~10,Kd取值范围为200~2000。
该公式的原理是:首先利用四元数的矢量部分作为比例反馈信息,角速度作为微分反馈信息,考虑卫星转动耦合力矩,得到期望的控制力矩,再根据磁矩、地磁场与控制力矩的叉乘关系,最终获取控制磁矩。
(4)、判断卫星姿态角Φ的幅值|Φ|是否小于预设的姿态角门限,如果小于,则认为捕获成功,结束捕获过程;否则,在下一个周期到来时,从步骤(1)开始执行,所述预设的姿态角门限取值范围[5°,15°]。
实施例
以某在轨卫星为例,该卫星运行在轨道高度为500km的太阳同步轨道卫星为例,轨道倾角为97.5°,在星上失去姿态基准的情况下,实际初始姿态为[-45,30,-50]°、初始角速度为[8,-9,7]°/s,陀螺漂移为[-23,21,-22]°/h,利用本发明进行了姿态捕获,其中姿态估计误差如图2所示、陀螺常值漂移估计误差如图3所示,对地三轴姿态如图4所示。
由图2、图3可知,在没有太阳信息的情况下本发明可以持续获取姿态信息,并且姿态估计误差和陀螺常值漂移估计误差随着时间的推移趋近于0,即姿态信息越来越准确;由图4可知,卫星的三轴姿态角趋近于0,在没有太阳信息的情况下实现了对地捕获,表明本发明具有很强的工程可操作和可实现性。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (6)
1.一种基于磁强计的皮纳卫星姿态捕获方法,其特征在于每个控制周期执行如下步骤:
(1)、以磁强计的测量值为观测量,以惯性坐标系下的卫星姿态四元数和陀螺的常值漂移为状态量,采用递推的卡尔曼滤波方法计算得到惯性坐标系下的卫星姿态四元数估计值和卫星姿态角速度ω信息,进入步骤(2);
(2)、判断惯性坐标系下的卫星的姿态角速度ω幅值|ω|是否小于预设的姿态角速度门限,如果小于,则进入步骤(3);否则,计算控制磁矩M驱动磁力矩器,以消除卫星姿态角速度,在下一个周期到来时,从步骤(1)开始执行;
(3)、将惯性坐标系下的卫星姿态四元数估计值转换为轨道坐标系下的卫星姿态四元数估计qBO,将惯性坐标系下的卫星姿态角速度ω转换成轨道系坐标系下的卫星姿态角速度ωBO,并根据所述轨道系坐标系下的卫星姿态四元数估计值qBO和卫星姿态角速度ωBO计算控制磁矩M驱动磁力矩器,以消除卫星的姿态角,进入步骤(4);
(4)、判断卫星姿态角Φ的幅值|Φ|是否小于预设的姿态角门限,如果小于,则认为捕获成功,结束捕获过程;否则,在下一个周期到来时,从步骤(1)开始执行。
2.根据权利要求1所述的一种基于磁强计的皮纳卫星姿态捕获方法,其特征在于所述步骤(1)计算惯性坐标系下的卫星姿态四元数估计值和卫星的姿态角速度ω信息的具体实现步骤为:
(1.1)、根据控制周期Δt内陀螺输出卫星三轴角度增量Δg和陀螺的常值漂移估计值bk-1/k-1计算卫星的姿态角速度ω:
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>g</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
式中,b的初值为陀螺的常值漂移标称值或者为[0,0,0];
(1.2)、根据上一控制周期的惯性坐标系下的卫星姿态四元数估计值qk-1/k-1外推当前控制周期惯性坐标系下的卫星姿态四元数预估值qk/k-1:
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>;</mo>
</mrow>
式中,q=[q1,q2,q3,q4],其中q4为标量,q的初始值为q0/0=[0,0,0,1];
(1.3)、计算滤波修正矩阵K:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mo>+</mo>
<mi>Q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>K</mi>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>+</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,P是滤波误差方差阵,Q是状态误差方差阵,与状态方程的随机误差特性相关,R是测量误差方差阵,与测量方程的随机误差特性相关,A是状态转移矩阵,取值为ωx,ωy,ωz分别为卫星的姿态角速度ω在卫星本体坐标系下的XYZ三轴分量,ω=[ωx,ωy,ωz],C是观测矩阵,取值为Bx、By、Bz为磁强计测量的地磁场B在卫星本体坐标系下的XYZ三轴的分量,B=[Bx,By,Bz];
(1.4)、根据当前控制周期磁强计测量的地磁场强度Bk与上一个控制周期磁强计测量的地磁场强度Bk-1之差,对惯性坐标系下的卫星姿态四元数预估值qk/k-1进行修正,获取当前控制周期的惯性坐标系下的卫星姿态四元数估值qk/k和陀螺常漂估值bk/k:
δX=K·(Bk-Bk-1);
δq=[δX(0) δX(1) δX(2)]T;
δb=[δX(3) δX(4) δX(5)]T;
qk/k=qk/k-1+E(qk/k-1)·K·δq;
bk/k=bk-1/k-1+K·δb。
3.根据权利要求1所述的一种基于磁强计的皮纳卫星姿态捕获方法,其特征在于:所述步骤(2)采用B-dot方法计算控制磁矩M驱动磁力矩器。
4.根据权利要求1所述的一种基于磁强计的皮纳卫星姿态捕获方法,其特征在于:所述步骤(3)采用比例+微分方法计算控制磁矩M驱动磁力矩器。
5.根据权利要求1所述的一种基于磁强计的皮纳卫星姿态捕获方法,其特征在于:所述预设的姿态角门限取值范围[5°,15°]。
6.根据权利要求1所述的一种基于磁强计的皮纳卫星姿态捕获方法,其特征在于:所述预设的姿态角速度门限取值范围[0.1°/s,1°/s]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711311392.3A CN108089434B (zh) | 2017-12-11 | 2017-12-11 | 一种基于磁强计的皮纳卫星姿态捕获方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711311392.3A CN108089434B (zh) | 2017-12-11 | 2017-12-11 | 一种基于磁强计的皮纳卫星姿态捕获方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108089434A true CN108089434A (zh) | 2018-05-29 |
CN108089434B CN108089434B (zh) | 2021-06-11 |
Family
ID=62174666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711311392.3A Active CN108089434B (zh) | 2017-12-11 | 2017-12-11 | 一种基于磁强计的皮纳卫星姿态捕获方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108089434B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110285815A (zh) * | 2019-05-28 | 2019-09-27 | 山东航天电子技术研究所 | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 |
CN116429096A (zh) * | 2023-02-22 | 2023-07-14 | 中国科学院国家空间科学中心 | 一种地磁惯性融合导航方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101381004A (zh) * | 2008-08-20 | 2009-03-11 | 南京航空航天大学 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
JP2011227017A (ja) * | 2010-04-23 | 2011-11-10 | Univ Of Tokyo | 慣性センサ,磁気センサおよび速度計を用いた移動体の姿勢推定装置および姿勢推定方法 |
CN102937450A (zh) * | 2012-10-31 | 2013-02-20 | 北京控制工程研究所 | 一种基于陀螺测量信息的相对姿态确定方法 |
CN104567871A (zh) * | 2015-01-12 | 2015-04-29 | 哈尔滨工程大学 | 一种基于地磁梯度张量的四元数卡尔曼滤波姿态估计方法 |
US20170003751A1 (en) * | 2015-06-30 | 2017-01-05 | Stmicroelectronics S.R.L. | Device and method for determination of angular position in three-dimensional space, and corresponding electronic apparatus |
-
2017
- 2017-12-11 CN CN201711311392.3A patent/CN108089434B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101381004A (zh) * | 2008-08-20 | 2009-03-11 | 南京航空航天大学 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
JP2011227017A (ja) * | 2010-04-23 | 2011-11-10 | Univ Of Tokyo | 慣性センサ,磁気センサおよび速度計を用いた移動体の姿勢推定装置および姿勢推定方法 |
CN102937450A (zh) * | 2012-10-31 | 2013-02-20 | 北京控制工程研究所 | 一种基于陀螺测量信息的相对姿态确定方法 |
CN104567871A (zh) * | 2015-01-12 | 2015-04-29 | 哈尔滨工程大学 | 一种基于地磁梯度张量的四元数卡尔曼滤波姿态估计方法 |
US20170003751A1 (en) * | 2015-06-30 | 2017-01-05 | Stmicroelectronics S.R.L. | Device and method for determination of angular position in three-dimensional space, and corresponding electronic apparatus |
Non-Patent Citations (2)
Title |
---|
张锐 等: ""基于磁强计的微小卫星姿态确定"", 《宇航学报》 * |
郝东: ""微小卫星姿态确定与磁控技术研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110285815A (zh) * | 2019-05-28 | 2019-09-27 | 山东航天电子技术研究所 | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 |
CN110285815B (zh) * | 2019-05-28 | 2023-10-20 | 山东航天电子技术研究所 | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 |
CN116429096A (zh) * | 2023-02-22 | 2023-07-14 | 中国科学院国家空间科学中心 | 一种地磁惯性融合导航方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108089434B (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Searcy et al. | Magnetometer-only attitude determination using novel two-step Kalman filter approach | |
CN106342284B (zh) | 一种飞行载体姿态确定方法 | |
CN104570742B (zh) | 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法 | |
CN108181916B (zh) | 小卫星相对姿态的控制方法及装置 | |
US9643740B2 (en) | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter | |
US20110024571A1 (en) | Gyroless transfer orbit sun acquisition using only wing current measurement feedback | |
CN105571578B (zh) | 一种利用伪观测取代精密转台的原地旋转调制寻北方法 | |
US8321076B2 (en) | On-line inertia estimation for use in controlling an aerospace vehicle | |
CN103090870A (zh) | 一种基于mems传感器的航天器姿态测量方法 | |
CN106597017A (zh) | 一种基于扩展卡尔曼滤波的无人机角加速度估计方法及装置 | |
CN109625335B (zh) | 一种基于角速度估计信息和太阳敏感器的捕获太阳方法 | |
CN103398713A (zh) | 一种星敏感器/光纤惯性设备量测数据同步方法 | |
CN104833375B (zh) | 一种借助星敏感器的imu两位置对准方法 | |
CN102116634A (zh) | 一种着陆深空天体探测器的降维自主导航方法 | |
CN110803304B (zh) | 一种卫星姿态控制系统 | |
CN103019247A (zh) | 一种火星探测器无陀螺自主空间姿态机动控制方法 | |
Li et al. | Low-cost MEMS sensor-based attitude determination system by integration of magnetometers and GPS: A real-data test and performance evaluation | |
CN109870273A (zh) | 基于动量守恒的航天器在轨质心辨识方法 | |
CN101226206A (zh) | 一种对加速度计零位偏差进行在轨标定的方法 | |
CN106767846A (zh) | 三轴稳定卫星不用陀螺的姿态获取方法和系统 | |
CN108089434A (zh) | 一种基于磁强计的皮纳卫星姿态捕获方法 | |
CN108871319B (zh) | 一种基于地球重力场与地磁场序贯修正的姿态解算方法 | |
CN107860382A (zh) | 一种在地磁异常情况下应用ahrs测量姿态的方法 | |
Babcock | CubeSat attitude determination via Kalman filtering of magnetometer and solar cell data | |
CN110667892B (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 |