一种对地三轴稳定观测卫星的磁强计在轨标定系统及标定方法技术方案

本发明专利技术提出了一种对地三轴稳定观测卫星的磁强计在轨标定系统及标定方法，首先根据有效的星敏数据，计算卫星本体坐标系下的计算磁场；接着根据有效的磁强计测量数据，得到卫星本体坐标系下的测量磁场；最后，根据计算磁场和测量磁场信息，结合设计的标定方法对磁强计进行标定；在卫星工作于对地三轴稳定模式时，准确的估计出测量标度系数和测量偏差系数，对磁强计测量值进行校准；使用校准后的磁强计信息，可以保证高效率的卸载及精确的双矢量姿态确定。确定。确定。

【技术实现步骤摘要】
一种对地三轴稳定观测卫星的磁强计在轨标定系统及标定方法


[0001]本专利技术属于航空航天
，具体地，涉及一种对地三轴稳定观测卫星的磁强计在轨标定系统及标定方法。

技术介绍

[0002]随着微电子技术的发展和卫星应用技术的成熟，微小卫星不仅能够实现传统大型航天器的功能，而且在设计灵活性、研制速度、研制成本和风险承受能力等方面也具有独特优势。目前，微小卫星技术被广泛应用于遥感、通信、导航、气象观测和科学试验等各个领域。
[0003]磁强计在实际工作中，地磁信号和各类误差信号耦合在一起，组成传感器输出，这将严重影响磁强计的精度。根据三轴磁强计的输出模型，可将外界干扰等效成三轴磁场数据的零偏误差和标度因数；三轴磁强计是卫星不可缺失的重要传感器，用于测量卫星所处环境的磁场强度大小和方向，在卫星长期稳定模式下可以完成姿态确定和动量卸载的任务，其性能关系着微小卫星在轨运行的稳定性和安全性。在传统卫星设计过程中，大多采用测量精度较高，稳定性较好的高精度磁强计。为了控制研制成本，微小卫星的姿态控制系统多采用质量轻、成本低的工业级磁强计。低成本的工业级磁强计由于自身性能的限制和相对较差的抗干扰能力，其测量的磁场信息存在一定的偏差，这导致卫星的动量卸载效率下降，以及双矢量姿态确定存在定姿偏差。

技术实现思路

[0004]针对长期工作于对地三轴稳定模式的微小卫星，本专利技术提出了一种对地三轴稳定观测卫星的磁强计在轨标定系统及标定方法，在卫星长期对地三轴稳定工作模式下，保证高效率的角动量卸载功能，并提供高精度的双矢量姿态确定信息，实现高指向精度和稳定度的控制。
[0005]本专利技术通过以下技术方案实现：
[0006]一种对地三轴稳定观测卫星的磁强计在轨标定方法：所述方法具体包括以下步骤：
[0007]步骤1，根据惯性坐标和卫星本体坐标，建立星载磁强计的测量模型；获取磁强计测量磁场，并判断磁强计是否有效，无效则重新获取测量磁场，有效则进行步骤2；
[0008]步骤2，建立星载计算磁场模型；获取星敏感器姿态数据及标志，判断星敏感器是否有效，无效则重新获取姿态数据及标志，有效则进行步骤3；
[0009]步骤3，将计算磁场变换至卫星本体坐标系；
[0010]步骤4，建立对地三轴稳定模式下的标定模型，将磁强计测量磁场变换至卫星本体坐标系；
[0011]步骤5，进行磁强计标定，最终得到测量标度系数K和测量偏差系数b
e
的所有变量。
[0012]进一步地，在步骤1中，
[0013]定义F
I
表示惯性坐标系，F
B
表示卫星本体坐标系；构建工业级磁强计的测量模型为：
[0014][0015]式中：B
m
为3
×
1的向量，表示磁强计坐标系下的测量磁场强度；B
b
为3
×
1的向量，表示卫星本体坐标系下校准后的测量磁场强度；
[0016]η
m
为3
×
1的向量，表示磁强计的测量随机噪声；K＝diag(k1,k2,k3)为3
×
3矩阵，表示磁强计的测量标度系数；
[0017]b为3
×
1的向量，表示磁强计的测量偏差系数；为3
×
3矩阵，表示卫星惯性坐标系至本体坐标系的旋转变换矩阵；
[0018]将公式简化可以得到：
[0019][0020]令表示卫星本体坐标系下校准前的测量磁场强度；表示卫星本体坐标系下总的测量偏差系数；表示卫星本体坐标系下总的测量随机噪声；
[0021]根据上述定义，三轴磁强计的测量模型最终简化为：
[0022]B
b
＝KB
bm
+b
e
+η
bm (3)。
[0023][0024]进一步地，在步骤2中，
[0025]星载的计算磁场模型采用国际地磁场参考模型(International Geomagnetic Reference Field，IGRF)，在地面或上空的任意一点，由IGRF模型将主磁场的标量位函数用球谐级数表示为：
[0026][0027]式中，a＝6371.2km为地球平均半径，r为球体的参考中心到球体表面的计算点的距离，r＝a+h，h为地程高；θ为地心余纬度；
[0028]λ为自格林尼治向东算起的经度；为经过标准化的施密特球谐系数，
[0029]为施密特准归一化n次m阶勒让德函数；
[0030]通过位场转换原理，可以将地磁场三个轴上的分量表示为以下形式：
[0031][0032]式中，及其一阶导数表示为：
[0033][0034][0035]假设通过IGRF模型得到的卫星所处位置的惯性坐标系下的磁场强度为B
Ir
，则
[0036][0037]式中表示地心坐标系至惯性坐标系的转换矩阵，且存在随机测量噪声为η
Ir
。
[0038]进一步地，在步骤3中，
[0039]卫星本体坐标系F
B
相对惯性坐标系F
I
的姿态表示为单位四元数并满足约束条件
[0040]惯性坐标系至卫星本体坐标系的旋转矩阵R(Q)有如下关系：
[0041][0042]式中，
[0043]则计算磁场在卫星本体坐标系下的分量B
bc
为：
[0044]B
bc
＝R(Q)(B
Ir
+η
Ir
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0045]定义B
br
＝R(Q)B
Ir
，η
br
＝R(Q)η
Ir
，则计算磁场的最终简化为：
[0046]B
bc
＝B
br
+η
br
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)。
[0047][0048]进一步地，在步骤4中，
[0049]采用递推最小二乘法对磁强计的参数进行标定，根据公式和(11)可以得到：
[0050][0051]式中：η＝η
bm
‑
η
br
；
[0052]向量B
br
,B
bm
,b
e
,η和矩阵K具体展开形式为：
[0053][0054]根据公式(12)和(13)可以得到：
[0055][0056]进一步地，在步骤5中，
[0057]根据公式(14)，定义如下变量：
[0058][0059]使用最下最小二乘法对变量x进行估计：
[0060][0061]分别令x＝[k
21 k
22 k
23 b
e本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种对地三轴稳定观测卫星的磁强计在轨标定方法，其特征在于：所述方法具体包括以下步骤：步骤1，根据惯性坐标和卫星本体坐标，建立星载磁强计的测量模型；获取磁强计测量磁场，并判断磁强计是否有效，无效则重新获取测量磁场，有效则进行步骤2；步骤2，建立星载计算磁场模型；获取星敏感器姿态数据及标志，判断星敏感器是否有效，无效则重新获取姿态数据及标志，有效则进行步骤3；步骤3，将计算磁场变换至卫星本体坐标系；步骤4，建立对地三轴稳定模式下的标定模型，将磁强计测量磁场变换至卫星本体坐标系；步骤5，进行磁强计标定，最终得到测量标度系数K和测量偏差系数b
e
的所有变量。2.根据权利要求1所述标定方法，其特征在于：在步骤1中，定义F
I
表示惯性坐标系，F
B
表示卫星本体坐标系；构建工业级磁强计的测量模型为：式中：B
m
为3
×
1的向量，表示磁强计坐标系下的测量磁场强度；B
b
为3
×
1的向量，表示卫星本体坐标系下校准后的测量磁场强度；η
m
为3
×
1的向量，表示磁强计的测量随机噪声；K＝diag(k1,k2,k3)为3
×
3矩阵，表示磁强计的测量标度系数；b为3
×
1的向量，表示磁强计的测量偏差系数；为3
×
3矩阵，表示卫星惯性坐标系至本体坐标系的旋转变换矩阵；将公式简化可以得到：令表示卫星本体坐标系下校准前的测量磁场强度；表示卫星本体坐标系下总的测量偏差系数；表示卫星本体坐标系下总的测量随机噪声；根据上述定义，三轴磁强计的测量模型最终简化为：B
b
＝KB
bm
+b
e
+η
bm (3)。3.根据权利要求2所述标定方法，其特征在于：在步骤2中，星载的计算磁场模型采用国际地磁场参考模型(International GeomagneticReference Field，IGRF)，在地面或上空的任意一点，由IGRF模型将主磁场的标量位函数用球谐级数表示为：式中，a＝6371.2km为地球平均半径，r为球体的参考中心到球体表面的计算点的距离，r＝a+h，h为地程高；θ为地心余纬度；λ为自格林尼治向东算起的经度；为经过标准化的施密特球谐系数，为施密特准归一化n次m阶勒让德函数；通过位场转换原理，可以将地磁场三个轴上的分量表示为以下形式：
式中，及其一阶导数表示为：表示为：假设通过IGRF模型得到的卫星所处位置的惯性坐标系下的磁场强度为B
Ir
，则式中表示地心坐标系至惯性坐标系的转换矩阵，且存在随机测量噪声为η
Ir
。4.根据权利要求3所述标定方法，其特征在于：在步骤3中，卫星本体坐标系F
B
相对惯性坐标系F
I
的姿态表示为单...

【专利技术属性】
技术研发人员：沈莹，孔令波，戴路，
申请(专利权)人：长光卫星技术股份有限公司，
类型：发明
国别省市：