CN104210677A - 一种用于磁控小卫星的磁强计补偿因子的确定方法 - Google Patents
一种用于磁控小卫星的磁强计补偿因子的确定方法 Download PDFInfo
- Publication number
- CN104210677A CN104210677A CN201410418758.7A CN201410418758A CN104210677A CN 104210677 A CN104210677 A CN 104210677A CN 201410418758 A CN201410418758 A CN 201410418758A CN 104210677 A CN104210677 A CN 104210677A
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- magnetometer
- centerdot
- moonlet
- whole star
- 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.)
- Pending
Links
Landscapes
- Measuring Magnetic Variables (AREA)
Abstract
本发明一种用于磁控小卫星的磁强计补偿因子的确定方法,包括步骤如下:整星放置在零磁试验室,调整卫星本体坐标系与试验室磁场坐标系三轴指向一致,利用磁强计测量得到在剩磁场环境下的整星剩磁场值;在零磁试验室中模拟真实的地磁场,通过线圈通电流的方法产生可控磁场,在卫星的三轴方向别加载模拟的地磁场强度值,并利用磁强计测量得到该模拟地磁场环境下的整星的磁感应强度;确定磁强计补偿因子。本发明最大的创新点在于解决了磁控小卫星的磁场测量精度问题,该问题的解决是该类需要磁强计的测量值以定姿的小卫星能够圆满完成航天任务的必要前提,具有非常关键的工程实际意义。
Description
技术领域
本发明涉及一种用于磁控小卫星的磁强计补偿因子的确定方法,属于空间环境试验测试技术领域。
背景技术
利用空间磁场环境控制卫星的姿态的技术,由于节省能源燃料、减轻卫星质量,降低发射成本而受到微小卫星平台的青睐。地磁场采用国际参考地磁场模型(IGRF)来描述,地磁场的强度和矢量是位置的函数,当确定了卫星的轨道,就可以求得地磁场矢量在轨道坐标系的分量(可记为BO).而磁强计的测量值是地磁场矢量在卫星本体的坐标系的分量(可记为Bb),理论上两者的关系应为:
Bb=CboBO
式中Cbo是卫星姿态方向余弦矩阵,这样就可以在磁强计测量值与卫星姿态运动学方程之间建立数学关系,再利用滤波算法即可得到卫星姿态角。
同时,卫星由于结构和性能的需要总要使用一些永磁材料和感磁材料,因此总会有一定的剩磁矩,卫星上仪器内和仪器之间连线中的电流也会产生磁矩。剩磁特性,可近似认为不随外界磁场变化而变,相当于有一个固定不变的磁矩,该剩磁可通过添加相反的硬磁性材料方式对剩磁进行补偿。感磁特性,是随外界磁场变化而变化的磁场,磁性的方向、大小都与外界磁场有关。
对于只能使用磁强计的数据进行定姿的一类磁控小卫星,磁强计的测量值是地磁场矢量在卫星本体的坐标系的分量,在磁强计探测器位置处的航天器固有磁场必须很低和稳定,以便估计仪器的探测精度,否则被测量的弱小值以及在时间和空间上的变化,会被航天器本身的磁场所淹没。因此,在该类卫星的磁试验中,必须要解决在轨运行中整星剩磁矩、动态磁矩对磁强计测量数据的磁噪声干扰问题;以及磁感应器件的测量精度问题。所以现在急需一种在卫星轨道运行中对磁强计测量结果进行补偿从而获得能够对真实磁场情况进行描述的方法。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种用于磁控小卫星的磁强计补偿因子的确定方法,通过地面模拟地磁场试验得到在卫星轨道运行中磁强计测量结果的补偿,从而能够实现空间环境地磁场环境的真实描述。
本发明的技术解决方案是:
一种用于磁控小卫星的磁强计补偿因子的确定方法包括步骤如下:
(1)整星放置在零磁试验室,调整卫星本体坐标系与试验室磁场坐标系三轴指向一致,利用磁强计测量得到在剩磁场环境下的整星剩磁场值Bx0、By0、Bz0;
(2)在零磁试验室中模拟真实的地磁场,通过线圈通电流的方法产生可控磁场,在卫星的三轴方向X、Y、Z分别加载模拟的地磁场强度值Bmx、Bmy、Bmz,并利用磁强计测量得到该模拟地磁场环境下的整星的磁感应强度Bxc、Byc、Bzc;
(3)确定磁强计补偿因子;
磁强计补偿因子K的计算公式如下:
其中,
所述步骤(1)中测量得到整星剩磁场值Bx0、By0、Bz0的具体方式如下:
(1a)将整星置于零磁场环境;
(1b)将卫星本体坐标系的三个坐标轴与零磁场环境的三个坐标轴对应;
(1c)对整星加电,利用磁强计测量出卫星不同的在轨工作模式下的整星剩磁场值Bx0、By0、Bz0,同时将测量得到多组剩磁场值Bx0、By0、Bz0进行求平均计算,得到一组整星剩磁场平均值,测量得到的多组剩磁场值Bx0、By0、Bz0与平均值误差在一定范围内的剩磁场值Bx0、By0、Bz0即可作为最终的整星剩磁场值Bx0、By0、Bz0。
本发明与现有技术相比的有益效果是:
(1)本发明最大的创新点在于解决了磁控小卫星的磁场测量精度问题,该问题的解决是该类需要磁强计的测量值以定姿的小卫星能够圆满完成航天任务的必要前提,具有非常关键的工程实际意义。
(2)本发明通过模拟实验计算得到了磁强计测量结果的补偿因子,解决了在轨运行中整星剩磁矩、动态磁矩对磁强计测量数据的磁噪声干扰问题;以及磁感应器件的测量精度问题,本发明简单可行,易于实现,尤其对于空间复杂环境来说,越简单的方法,运行的可靠性越高。
(3)在目前现有的地磁场环境的模拟中,是通过加载一个三轴叠加的电磁场,而本发明在零磁环境中模拟真实的地磁场,在磁强计的三轴分别加载模拟的地磁场强度值,并利用磁强计分别测量得到该模拟地磁场环境下的整星的三轴磁感应强度,这是本发明的一大亮点,能够准确分析三轴电磁场的强度值,进而能够分析三轴补偿因子,分析更加准确,大大提高了工作效率。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
如图1所示,本发明一种用于磁控小卫星的磁强计补偿因子的确定方法,包括步骤如下:
(1)整星放置在零磁试验室,调整卫星本体坐标系与试验室磁场坐标系三轴指向一致,利用磁强计测量得到在剩磁场环境下的整星剩磁场值Bx0、By0、Bz0;
测量得到整星剩磁场值Bx0、By0、Bz0的具体方式如下:
(1a)将整星置于零磁场环境;
(1b)将卫星本体坐标系的三个坐标轴与零磁场环境的三个坐标轴对应;
(1c)对整星加电,利用磁强计测量出卫星不同的在轨工作模式下的整星剩磁场值Bx0、By0、Bz0,同时将测量得到多组剩磁场值Bx0、By0、Bz0进行求平均计算,得到一组整星剩磁场平均值,测量得到的多组剩磁场值Bx0、By0、Bz0与平均值误差在一定范围内的剩磁场值Bx0、By0、Bz0即可作为最终的整星剩磁场值Bx0、By0、Bz0。
(2)在零磁试验室中模拟真实的地磁场,通过线圈通电流的方法产生可控磁场,在卫星的三轴方向X、Y、Z分别加载模拟的地磁场强度值Bmx、Bmy、Bmz,并利用磁强计测量得到该模拟地磁场环境下的整星的磁感应强度Bxc、Byc、Bzc;
例如:模拟地磁场环境的具体输出如下:
(2a)模拟地磁场,使得磁场Y、Z轴输出磁场为零,在磁强计的X轴分别加0、-10000、-20000、-30000、-40000nT;在星上各部件处于正常工作状态时利用磁强计进行测量,得到Bxc、Byc、Bzc如表1所示:
表1 X轴磁场数据
从表1中可以看出
(2b)模拟地磁场,使得磁场X、Z轴输出磁场为零,在磁强计的Y轴分别加0、-20000、-10000、10000、20000nT;在星上各部件处于正常工作状态时利用磁强计进行测量,得到Bxc、Byc、Bzc如表2所示:
表2 Y轴磁场数据
从表2中可以看出
(2c)模拟地磁场,使得磁场X、Y轴输出磁场为零,在Z轴分别加0、10000、20000、30000、40000nT;在星上各部件处于正常工作状态时利用磁强计进行测量,得到Bxc、Byc、Bzc如表3所示:
表3 Z轴磁场数据
从表1中可以看出
(3)确定磁强计补偿因子;
磁场是一个矢量叠加场,要完全描述它需要3个完全独立的坐标,同时,对于测量元件磁强计的测量结果也表示为直角坐标系中X、Y、Z三个方向的场的相互叠加量。因此,在有剩磁矩的环境下,磁强计的零位会发生变化,但是由于剩磁矩不随环境磁场变化而变化,因此确定了该磁强计补偿因子算法。
另外,在有感磁的环境下,感应磁矩是随环境磁场变化而变化的,并且其变化与环境磁场大小或正负、方向有关,感应磁场对磁强计测量的关系表示为下式中矩阵K和模拟地磁场强度值Bmx、Bmy、Bmz。其中,矩阵K的物理意义在于Kxx为整星X轴的感磁引起的磁强计X方向的自感应系数,Kyx为整星Y轴的感磁引起的磁强计X方向的自感应系数、Kzx整星Z轴的感磁引起的磁强计X方向的自感应系数(包含了磁强计产品在整星的安装误差、产品自身的不平行度影响);Kxy、Kyy、Kzy、Kxz、Kyz、Kzz的物理意义以此类推。
磁强计补偿因子K的计算公式如下:
其中,
例如,本发明已经成功应用在相应的型号上,以下面一组具体数例进行说明:表4为磁强计补偿因子K的值:
表4 K阵系数
利用磁强计补偿因子K补偿磁强计后对于模拟磁场的测量结果:
表5 模拟磁场的测量结果
可以明显看出,三方向加任意磁场,修正前磁场测量最大误差为712nT,修正后三轴最大误差不超过83.1nT,补偿效果明显,达到了工程预期试验目的。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (2)
1.一种用于磁控小卫星的磁强计补偿因子的确定方法,其特征在于步骤如下:
(1)整星放置在零磁试验室,调整卫星本体坐标系与试验室磁场坐标系三轴指向一致,利用磁强计测量得到在剩磁场环境下的整星剩磁场值Bx0、By0、Bz0;
(2)在零磁试验室中模拟真实的地磁场,通过线圈通电流的方法产生可控磁场,在卫星的三轴方向X、Y、Z分别加载模拟的地磁场强度值Bmx、Bmy、Bmz,并利用磁强计测量得到该模拟地磁场环境下的整星的磁感应强度Bxc、Byc、Bzc;
(3)确定磁强计补偿因子;
磁强计补偿因子K的计算公式如下:
其中,
2.根据权利要求1所述的一种用于磁控小卫星的磁强计补偿因子的确定方法,其特征在于在:所述步骤(1)中测量得到整星剩磁场值Bx0、By0、Bz0的具体方式如下:
(1a)将整星置于零磁场环境;
(1b)将卫星本体坐标系的三个坐标轴与零磁场环境的三个坐标轴对应;
(1c)对整星加电,利用磁强计测量出卫星不同的在轨工作模式下的整星剩磁场值Bx0、By0、Bz0,同时将测量得到多组剩磁场值Bx0、By0、Bz0进行求平均计算,得到一组整星剩磁场平均值,测量得到的多组剩磁场值Bx0、By0、Bz0与平均值误差在一定范围内的剩磁场值Bx0、By0、Bz0即可作为最终的整星剩磁场值Bx0、By0、Bz0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410418758.7A CN104210677A (zh) | 2014-08-22 | 2014-08-22 | 一种用于磁控小卫星的磁强计补偿因子的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410418758.7A CN104210677A (zh) | 2014-08-22 | 2014-08-22 | 一种用于磁控小卫星的磁强计补偿因子的确定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104210677A true CN104210677A (zh) | 2014-12-17 |
Family
ID=52092677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410418758.7A Pending CN104210677A (zh) | 2014-08-22 | 2014-08-22 | 一种用于磁控小卫星的磁强计补偿因子的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104210677A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104748762A (zh) * | 2015-03-13 | 2015-07-01 | 西北工业大学 | 一种高性能地磁场模拟装置设计及制作方法 |
CN105182254A (zh) * | 2015-07-31 | 2015-12-23 | 江汉大学 | 静态弱磁场检测装置 |
CN106556852A (zh) * | 2016-11-21 | 2017-04-05 | 重庆卓观科技有限公司 | 一种基于定位姿态信息模拟磁罗盘的方法与磁强计模拟器 |
CN108802653A (zh) * | 2018-06-15 | 2018-11-13 | 上海卫星工程研究所 | 基于地磁场特性的卫星磁强计磁环境分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1355421A (zh) * | 2001-12-28 | 2002-06-26 | 清华大学 | 应力分布的磁检测方法 |
US20050283988A1 (en) * | 2004-06-11 | 2005-12-29 | Yamaha Corporation | Method and apparatus for measuring magnetic offset of geomagnetic sensor and portable electronic apparatus |
CN101887068A (zh) * | 2010-06-01 | 2010-11-17 | 中北大学 | 三轴矢量传感器及两轴矢量传感器的标定补偿方法 |
CN102589537A (zh) * | 2012-03-05 | 2012-07-18 | 无锡汉和航空技术有限公司 | 一种有磁环境下无人机的电子罗盘校正方法 |
-
2014
- 2014-08-22 CN CN201410418758.7A patent/CN104210677A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1355421A (zh) * | 2001-12-28 | 2002-06-26 | 清华大学 | 应力分布的磁检测方法 |
US20050283988A1 (en) * | 2004-06-11 | 2005-12-29 | Yamaha Corporation | Method and apparatus for measuring magnetic offset of geomagnetic sensor and portable electronic apparatus |
CN101887068A (zh) * | 2010-06-01 | 2010-11-17 | 中北大学 | 三轴矢量传感器及两轴矢量传感器的标定补偿方法 |
CN102589537A (zh) * | 2012-03-05 | 2012-07-18 | 无锡汉和航空技术有限公司 | 一种有磁环境下无人机的电子罗盘校正方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104748762A (zh) * | 2015-03-13 | 2015-07-01 | 西北工业大学 | 一种高性能地磁场模拟装置设计及制作方法 |
CN104748762B (zh) * | 2015-03-13 | 2018-06-01 | 西北工业大学 | 一种高性能地磁场模拟装置设计及制作方法 |
CN105182254A (zh) * | 2015-07-31 | 2015-12-23 | 江汉大学 | 静态弱磁场检测装置 |
CN105182254B (zh) * | 2015-07-31 | 2018-01-05 | 江汉大学 | 静态弱磁场检测装置 |
CN106556852A (zh) * | 2016-11-21 | 2017-04-05 | 重庆卓观科技有限公司 | 一种基于定位姿态信息模拟磁罗盘的方法与磁强计模拟器 |
CN108802653A (zh) * | 2018-06-15 | 2018-11-13 | 上海卫星工程研究所 | 基于地磁场特性的卫星磁强计磁环境分析方法 |
CN108802653B (zh) * | 2018-06-15 | 2020-11-13 | 上海卫星工程研究所 | 基于地磁场特性的卫星磁强计磁环境分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107024674B (zh) | 一种基于递推最小二乘法的磁强计现场快速标定方法 | |
Fang et al. | A novel calibration method of magnetic compass based on ellipsoid fitting | |
CN105425764B (zh) | 一种基于动态地磁场模拟的三轴磁强计闭环测试***及方法 | |
CN106353824A (zh) | 航空磁通门磁梯度张量仪的***校正及磁干扰补偿融合方法 | |
Pang et al. | Calibration of a fluxgate magnetometer array and its application in magnetic object localization | |
CN105891755A (zh) | 航空器挂载式磁通门磁梯度张量仪的校正方法 | |
CN104210677A (zh) | 一种用于磁控小卫星的磁强计补偿因子的确定方法 | |
CN106997035A (zh) | 一种基于磁梯度不变量的磁梯度计校正方法 | |
CN104199115B (zh) | 一种基于线性模型的地磁矢量测量误差综合补偿方法 | |
Pang et al. | Calibration of three-axis magnetometers with differential evolution algorithm | |
Wu et al. | Constrained total least-squares calibration of three-axis magnetometer for vehicular applications | |
CN107544042A (zh) | 一种磁力计阵列校正方法 | |
CN104406610A (zh) | 一种磁力计实时校正装置及方法 | |
CN109444774B (zh) | 基于小信号法测量水下航行器磁场噪声系数的方法 | |
Lassakeur et al. | Magnetic cleanliness program on CubeSats and nanosatellites for improved attitude stability | |
Chi et al. | Calibration of triaxial magnetometer with ellipsoid fitting method | |
CN113156355B (zh) | 一种超导全张量磁梯度测量装置的磁干扰补偿方法 | |
Liu et al. | Compensation of geomagnetic vector measurement system with differential magnetic field method | |
Yin et al. | Recalibration of the CHAMP satellite magnetic field measurements | |
Ke et al. | Magnetometer compensation scheme and experimental results on ZDPS-1A pico-satellite | |
Wang et al. | Compensation for mobile carrier magnetic interference in a SQUID-based full-tensor magnetic gradiometer using the flower pollination algorithm | |
Lassakeur et al. | Enhanced attitude stability and control for CubeSats by real-time on-orbit determination of their dynamic magnetic moment | |
Gao et al. | A calibration method for the misalignment error between inertial navigation system and tri-axial magnetometer in three-component magnetic measurement system | |
Dorveaux et al. | On-the-field calibration of an array of sensors | |
Heilig | Determining the orthogonality error of coil systems by means of a scalar magnetometer: application to delta inclination–delta declination (dIdD) magnetometers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141217 |