CN112034402A - 微纳卫星的剩磁和剩磁矩联合标定方法 - Google Patents
微纳卫星的剩磁和剩磁矩联合标定方法 Download PDFInfo
- Publication number
- CN112034402A CN112034402A CN202010685634.0A CN202010685634A CN112034402A CN 112034402 A CN112034402 A CN 112034402A CN 202010685634 A CN202010685634 A CN 202010685634A CN 112034402 A CN112034402 A CN 112034402A
- Authority
- CN
- China
- Prior art keywords
- remanence
- micro
- moment
- nano satellite
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种微纳卫星的剩磁和剩磁矩联合标定方法,该方法包括:建立微纳卫星动力学模型;根据动力学模型建立UKF滤波模型;将陀螺仪输出的速度增量作为测量基准,利用UKF算法确定微纳卫星的剩磁和剩磁矩。本发明的微纳卫星剩磁和剩磁矩联合标定方法,通过在轨算法标定,可以去除整星剩磁干扰和采集引起的磁强计测量误差,有效提高磁强计的测量精度,满足在轨应用磁强计测量精度的要求。在轨标定过程降低了微纳卫星的研制成本,标定效率较高。通过陀螺仪的输出作为测量基准,操作简单,准确性高。
Description
技术领域
本发明涉及航天器控制技术领域,尤其涉及一种微纳卫星的剩磁和剩磁矩联合标定方法。
背景技术
磁强计是微纳卫星最常用的敏感器,因具有价格低廉、质量轻、体积小、采样电路简单、性能稳定等优点,在卫星中得到广泛应用。为了降低卫星本体扰动磁场对磁强计测量值的影响,大卫星通常将磁强计通过长杆伸出到卫星外部,由于磁场强度与磁体距离的三次方呈反比,因此增加距离可以降低磁场强度,从而提高磁强计的测量精度。但是微纳卫星由于收到体积、质量的限制,通常只能将磁强计安装在微纳卫星的星体内部,与其他部组件之间的距离较近。因此,当星体内的其他组件工作产生磁场时,会严重影响磁强计的测量精度。当磁强计测量的地磁场与真实磁场存在较大偏差,基于该测量值计算的卫星姿态精度较低。此外,微纳卫星通常采用体装太阳电池阵,磁强计安装在太阳电池阵形成的空腔内,其产生的时变磁场对磁强计的影响更不可忽略。
发明内容
为解决上述现有技术中存在的技术问题,本发明提供了一种微纳卫星的剩磁和剩磁矩联合标定方法。具体技术方案如下:
一种微纳卫星的剩磁和剩磁矩联合标定方法,所述方法包括:
建立微纳卫星动力学模型;
根据动力学模型建立UKF滤波模型;
将陀螺仪输出的速度增量作为测量基准,利用UKF算法确定微纳卫星的剩磁和剩磁矩。
可选地,动力学模型包括:
式中,τb表示时间常数,ηb为系统噪声。
可选地,动力学模型还包括:
其中,τm表示时间常数,ηm为系统噪声。
可选地,动力学模型还包括:
式中,J表示卫星转动惯量矩阵,m0表示卫星磁矩,Bm表示磁强计的测量输出值,B0表示磁强计的测量偏值,Tc表示主动控制力矩,ω表示卫星姿态角速度。
可选地,根据动力学模型建立UKF滤波模型,进一步包括:
选取系统的状态变量为
建立系统动力学方程为
其中,
其中,03表示3阶零矩阵;
在[tk,tk+1]时间区间内,对系统动力学方程进行离散化可得:
Xk+1|k=Xk+f(Xk)Δt
其中,Δt表示采样间隔时间。
可选地,将陀螺仪输出的速度增量作为观测方程输出,则有Z=Δω;
系统观测方程为:
Zk=HkXk
其中,观测矩阵Hk=[I3 03],I3表示3阶单位阵;
根据系统观测方程和离散化后的系统动力学方程,根据标准UKF滤波流程就确定状态参数B0、m0的值。
本发明技术方案的主要优点如下:
本发明的微纳卫星的剩磁和剩磁矩联合标定方法,通过在轨算法标定,可以去除整星剩磁干扰和采集引起的磁强计测量误差,有效提高磁强计的测量精度,满足在轨应用磁强计测量精度的要求。且在轨标定过程降低了微纳卫星的研制成本,标定效率较高。通过陀螺仪的输出作为测量基准,操作简单,准确性高。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明一实施例提供的微纳卫星的剩磁和剩磁矩联合标定方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下结合附图,详细说明本发明实施例提供的技术方案。
本发明实施例提供了一种微纳卫星的剩磁和剩磁矩联合标定方法,如附图1所示,该方法包括:
建立微纳卫星动力学模型;
根据动力学模型建立UKF滤波模型;
将陀螺仪输出的速度增量作为测量基准,利用UKF算法确定微纳卫星的剩磁和剩磁矩。
本发明实施例提供的微纳卫星的剩磁和剩磁矩联合标定方法,通过在轨算法标定,可以去除整星剩磁干扰和采集引起的磁强计测量误差,有效提高磁强计的测量精度,满足在轨应用磁强计测量精度的要求。且在轨标定过程降低了微纳卫星的研制成本,标定效率较高。通过陀螺仪的输出作为测量基准,操作简单,准确性高。
具体地,以下对本发明实施例提供的微纳卫星的剩磁和剩磁矩联合标定方法的具体原理进行说明:
由于磁强计作为敏感器件经过严格的标定,所以磁强计本身的零偏远小于由剩磁引起的偏差。因此,磁强计的输出可以表示为:
Bm=Bb+B0 (1)
其中,Bm为磁强计的测量输出值,Bb为卫星当前位置地磁场强度真值,B0为磁强计的测量偏值。
由于微纳卫星上的剩磁是慢时变的,则磁强计的测量偏差值B0可以建模为:
其中,τb表示时间常数,ηb为系统噪声。
此外,微纳卫星的磁矩m0也是慢时变的,同理可以建模为:
其中,τm表示时间常数,ηm为系统噪声。
剩磁和剩磁矩相互作用,会对微纳卫星的姿态产生干扰力矩Td,卫星受到的剩磁干扰力矩Td可以表示为:
Td=m0×Bm-m0×B0 (4)
由于微纳卫星一般均为刚性的,因此微纳卫星的姿态动力学模型可以表示为:
其中,ω表示卫星姿态角速度,J表示卫星转动惯量矩阵,Tc表示主动控制力矩。
根据动力学模型式(2)、式(3)和式(7)建立UKF滤波模型。
选取系统的状态变量为:
建立系统动力学方程为:
其中,
03表示3阶零矩阵。
在[tk,tk+1]时间区间内,对上式进行离散化可得:
Xk+1|k=Xk+f(Xk)Δt (10)
其中,Δt表示采样间隔时间。
将陀螺仪输出的速度增量作为观测方程输出,则有
Z=Δω (11)
则系统观测方程可以表示为:
Zk=HkXk (12)
其中,观测矩阵Hk=diag([I3 03 03]),I3表示3阶单位阵。
根据式(10)和式(12),根据标准UKF滤波流程就能确定状态参数B0、m0的值,从而在轨估计出微纳卫星上的剩磁和剩磁矩。
需要说明的是,在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。此外,本文中“前”、“后”、“左”、“右”、“上”、“下”均以附图中表示的放置状态为参照。
最后应说明的是:以上实施例仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种微纳卫星的剩磁和剩磁矩联合标定方法,其特征在于,所述方法包括:
建立微纳卫星动力学模型;
根据动力学模型建立UKF滤波模型;
将陀螺仪输出的速度增量作为测量基准,利用UKF算法确定微纳卫星的剩磁和剩磁矩。
6.根据权利要求5所述的微纳卫星的剩磁和剩磁矩联合标定方法,其特征在于,将陀螺仪输出的速度增量作为观测方程输出,则有Z=Δω;
系统观测方程为:
Zk=HkXk
其中,观测矩阵Hk=[I3 03],I3表示3阶单位阵;
根据系统观测方程和离散化后的系统动力学方程,根据标准UKF滤波流程就确定状态参数B0、m0的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010685634.0A CN112034402B (zh) | 2020-07-16 | 2020-07-16 | 微纳卫星的剩磁和剩磁矩联合标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010685634.0A CN112034402B (zh) | 2020-07-16 | 2020-07-16 | 微纳卫星的剩磁和剩磁矩联合标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112034402A true CN112034402A (zh) | 2020-12-04 |
CN112034402B CN112034402B (zh) | 2021-12-28 |
Family
ID=73579377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010685634.0A Active CN112034402B (zh) | 2020-07-16 | 2020-07-16 | 微纳卫星的剩磁和剩磁矩联合标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112034402B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113471319A (zh) * | 2021-06-29 | 2021-10-01 | 中国电子科技集团公司第十八研究所 | 一种空间用太阳电池阵剩磁矩的消除结构及方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226561A (zh) * | 2007-12-28 | 2008-07-23 | 南京航空航天大学 | 用于微型航天器姿态轨道控制系统的微型仿真支持系统及工作方法 |
JP2014025715A (ja) * | 2012-07-24 | 2014-02-06 | Fujikura Ltd | 磁気測定データ校正装置及び方位角計測装置 |
CN103954286A (zh) * | 2014-04-24 | 2014-07-30 | 南京航空航天大学 | 微小卫星磁传感器多误差模型在轨迭代标定方法 |
CN107270940A (zh) * | 2017-07-11 | 2017-10-20 | 深圳清华大学研究院 | 三轴磁强计的输出偏差估计方法 |
CN107422281A (zh) * | 2017-04-28 | 2017-12-01 | 中国人民解放军国防科学技术大学 | 微纳卫星磁场测量偏差修正方法及装置 |
CN109655070A (zh) * | 2018-12-28 | 2019-04-19 | 清华大学 | 一种遥感微纳卫星的多模式姿态确定方法 |
CN110736484A (zh) * | 2019-11-18 | 2020-01-31 | 北京理工大学 | 基于陀螺仪及磁传感器融合的背景磁场标定方法 |
-
2020
- 2020-07-16 CN CN202010685634.0A patent/CN112034402B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226561A (zh) * | 2007-12-28 | 2008-07-23 | 南京航空航天大学 | 用于微型航天器姿态轨道控制系统的微型仿真支持系统及工作方法 |
JP2014025715A (ja) * | 2012-07-24 | 2014-02-06 | Fujikura Ltd | 磁気測定データ校正装置及び方位角計測装置 |
CN103954286A (zh) * | 2014-04-24 | 2014-07-30 | 南京航空航天大学 | 微小卫星磁传感器多误差模型在轨迭代标定方法 |
CN107422281A (zh) * | 2017-04-28 | 2017-12-01 | 中国人民解放军国防科学技术大学 | 微纳卫星磁场测量偏差修正方法及装置 |
CN107270940A (zh) * | 2017-07-11 | 2017-10-20 | 深圳清华大学研究院 | 三轴磁强计的输出偏差估计方法 |
CN109655070A (zh) * | 2018-12-28 | 2019-04-19 | 清华大学 | 一种遥感微纳卫星的多模式姿态确定方法 |
CN110736484A (zh) * | 2019-11-18 | 2020-01-31 | 北京理工大学 | 基于陀螺仪及磁传感器融合的背景磁场标定方法 |
Non-Patent Citations (2)
Title |
---|
华冰 等: "基于光照区太阳敏感器和陀螺辅助修正的微小卫星磁测技术", 《中国惯性技术学报》 * |
郁丰 等: "微小卫星剩磁在轨标定技术研究", 《航空学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113471319A (zh) * | 2021-06-29 | 2021-10-01 | 中国电子科技集团公司第十八研究所 | 一种空间用太阳电池阵剩磁矩的消除结构及方法 |
CN113471319B (zh) * | 2021-06-29 | 2022-06-07 | 中国电子科技集团公司第十八研究所 | 一种空间用太阳电池阵剩磁矩的消除结构及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112034402B (zh) | 2021-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | Fast complementary filter for attitude estimation using low-cost MARG sensors | |
CN106873611B (zh) | 一种多通道线性自抗扰控制器的设计方法 | |
US8326533B2 (en) | Apparatus and methodology for calibration of a gyroscope and a compass included in a handheld device | |
CN106990426B (zh) | 一种导航方法和导航装置 | |
CN108195404B (zh) | 一种陀螺仪零点漂移的校准方法 | |
CN107024674B (zh) | 一种基于递推最小二乘法的磁强计现场快速标定方法 | |
CN100405014C (zh) | 一种载体姿态测量方法 | |
CN101290326B (zh) | 石英挠性加速度计测量组件的参数辨识标定方法 | |
CN108036785A (zh) | 一种基于直接法与惯导融合的飞行器位姿估计方法 | |
CN102914972B (zh) | 基于模型整体逼近的微陀螺仪rbf网络自适应控制方法 | |
CN100559123C (zh) | 一种mems陀螺仪的差分测量方法 | |
CN108592952A (zh) | 基于杆臂补偿与正反倍速率同时标定多mimu误差的方法 | |
CN101726295A (zh) | 考虑加速度补偿和基于无迹卡尔曼滤波的惯性位姿跟踪方法 | |
CN110160554A (zh) | 一种基于寻优法的单轴旋转捷联惯导系统标定方法 | |
CN109655070A (zh) | 一种遥感微纳卫星的多模式姿态确定方法 | |
CN109709628B (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
CN105043348A (zh) | 基于卡尔曼滤波的加速度计陀螺仪水平角度测量方法 | |
CN104483973A (zh) | 基于滑模观测器的低轨挠性卫星姿态跟踪控制方法 | |
CN113203429B (zh) | 一种陀螺仪温度漂移误差的在线估计及补偿方法 | |
CN110361031A (zh) | 一种基于回溯理论的imu全参数误差快速标定方法 | |
Li et al. | Low-cost MEMS sensor-based attitude determination system by integration of magnetometers and GPS: A real-data test and performance evaluation | |
CN111121820B (zh) | 基于卡尔曼滤波的mems惯性传感器阵列融合方法 | |
CN112034402B (zh) | 微纳卫星的剩磁和剩磁矩联合标定方法 | |
Chou et al. | Two‐step optimal filter design for the low‐cost attitude and heading reference systems | |
Wang et al. | An efficient calibration method for triaxial gyroscope |
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 |