一种基于查表法确定空间磁场强度的方法技术

本发明专利技术公开了一种基于查表法确定空间磁场强度的方法，该方法包含：(1)建立高精度磁场表；(2)修正卫星飞行轨迹降交点地理经度间隔Gdjg以及即将飞过的降交点地理经度lamda；(3)在星箭分离后，星载软件首次调用磁场表子程序时，已经飞过的降交点地理经度ccl由lamda赋值；其它情况下，当卫星飞到降交点时更新已经飞过的降交点地理经度ccl；(4)当ccl发生变化时，计算卫星当前位置在磁场表中所对应磁场轨道的编号cca及其降交点地理经度；(5)根据卫星当前位置在磁场表中的cca查找出卫星当前位置四周的4个网格点；(6)基于线性插值的方法计算卫星当前位置在轨道坐标系下的磁场强度。本发明专利技术的方法不受阶数的限制，简单易实现，运算快且计算精度较高。

【技术实现步骤摘要】
一种基于查表法确定空间磁场强度的方法
本专利技术涉及一种确定空间磁场强度的方法，具体涉及一种基于查表法确定空间磁场强度的方法。
技术介绍
地磁场是重要的地球物理场，在低轨卫星上地磁场应用前景十分广泛。将地磁矢量与地心矢量或者太阳矢量通过双矢量定姿的方式可粗略的估算卫星姿态。在常规大卫星中，该姿态可作为安全模式时的备份姿态基准，在对于定姿精度要求不高的小型化卫星上，更是可作为稳态时的主份姿态基准。对于三轴对地定向卫星，地磁矢量在卫星本体系下的投影由三轴磁强计测量得到，而该矢量在轨道坐标系下的投影通常根据国际地磁场模型(IGRF)理论计算得到。理论计算得到的磁场强度，除了可配合磁强计用于双矢量定姿，还可单独作为磁强计的备份，在卫星姿态稳定时用于磁卸载。IGRF是描述全球地磁场及其长期变化的，而且是国际上通用的全球地磁标准模型，在该模型中，地磁场的标量位用球谐级数表示，目前最新的球谐级数的阶数为13阶。通过该模型计算磁场强度时所取的阶数越高，则得到的磁场强度精度越高，但球谐级数由多项多阶的三角函数组成，随着球谐级数阶数的增加，计算量以二次方的方式激增。根据国内外已发表的相关文献可知，目前卫星上通常直接在星载软件中按照IGRF模型公式计算磁场强度，受限于星载计算机的运算能力，球谐级数的阶数一般取1阶或者3阶。但是，上述方法均存在一定的局限性。阶数取1阶，虽然公式较为简单，计算量小，但计算得到的磁场强度最大误差高达28％；阶数取3阶，虽然最大误差下降到了8％，但误差仍然偏大且占用了星载计算机相对较多地运行时间。
技术实现思路
本专利技术的目的是提供一种基于查表法确定空间磁场强度的方法，该方法解决了现有方法的误差偏大的问题，不受阶数的限制，算法简单，星载软件容易实现，运算时间较短且计算出的磁场强度精度较高。为了达到上述目的，本专利技术提供了一种基于查表法确定空间磁场强度的方法，该方法包含：(1)建立高精度磁场表：在卫星飞行轨迹经过的球面上按照一定的时间间隔进行网格划分，根据国际地磁场模型最高阶公式计算出网格点上的三轴磁场强度，形成地磁场数据表，将其装载在卫星星载软件中；(2)修正卫星飞行轨迹降交点地理经度间隔Gdjg以及即将飞过的降交点地理经度lamda；(3)在星箭分离后，星载软件首次调用磁场表子程序时，已经飞过的降交点地理经度ccl由即将飞过的降交点地理经度lamda赋值；其它情况下，当卫星飞到降交点时更新已经飞过的降交点地理经度ccl；(4)当已经飞过的降交点地理经度ccl发生变化时，计算卫星当前位置在磁场表中所对应磁场轨道的编号cca及其降交点地理经度，否则保持cca不变；其中，降交点地理经度根据步骤(3)获得，磁场轨道的编号cca与网格划分相对应；(5)根据卫星当前位置在磁场表中所对应磁场轨道的编号cca查找出卫星当前位置四周的4个网格点；(6)基于线性插值的方法计算卫星当前位置在轨道坐标系下的磁场强度。优选地，在步骤(1)中，所述网格划分的方式为：横向以降交点地理经度等间隔分成n轨，n为卫星每天飞行圈数向下取整；纵向以卫星飞行的一轨时间等间隔取m个点，作为轨道钟。优选地，横向每两轨之间的地理经度间隔ΔL为360°/n，地理经度与cca的关系为：地理经度＝(cca-1)*360°/n；纵向中，一轨时间等间隔取m为70～80。优选地，所述网格划分的数据点从东经0°的降交点时刻开始按顺序排列。优选地，在步骤(2)中，根据卫星实际运行的轨道高度修正降交点地理经度间隔Gdjg；根据卫星精测的轨道参数适时修正即将飞过的降交点地理经度lamda。优选地，在步骤(3)中，更新已经飞过的降交点地理经度ccl，其通过如下公式(1)获得：ccl＝ccl-Gdjg(1)。优选地，在步骤(6)中，在4个网格点中，通过处于同一纵向位置的两个网格点根据线性插值方法计算两个中间点的磁场强度，两个中间点分别中间点一和中间点二，中间点一和中间点二与卫星当前位置处于同一横向位置，根据计算得到的两个中间点的磁场强度通过线性插值方法获得卫星当前位置的磁场强度。优选地，在步骤(6)中，所述线性插值方法的公式为：在对4个网格点计算处于同一纵向位置的两个网格点的中间点时，式(2)中，Y1为横向上处于卫星当前位置前一网格点在磁场表中对应的磁场强度，Y2为横向上处于卫星当前位置后一网格点在磁场表中对应的磁场强度；X为卫星已经飞过的降交点地理经度ccl，Y为待求的中间点的磁场强度；X1为卫星当前位置对应磁场表中下一条磁场轨道的降交点地理经度；X2为卫星当前位置对应磁场表中上一条磁场轨道的降交点地理经度。在通过两个中间点计算卫星当前位置的磁场强度时，中间点一为卫星当前位置下一轨道钟对应的网格点，中间点二为卫星当前位置上一轨道钟对应的网格点，式(2)中，Y1为中间点一的磁场强度，Y2为中间点二的磁场强度；X为卫星当前时刻对应的轨道钟；Y为卫星当前位置的磁场强度；X1为中间点一在磁场表中对应的轨道钟；X2为中间点二在磁场表中对应的轨道钟。优选地，所述的轨道钟在卫星经过降交点时从0开始累加计时。本专利技术的基于查表法确定空间磁场强度的方法，解决了现有方法的误差偏大的问题，具有以下优点：本专利技术的方法通过建立高精度磁场表，将其预置在星载软件中，卫星在轨运行时通过查表的方式获取当前位置磁场强度的方式，解决了星上在计算磁场强度理论值时受限于星载计算机能力而只能采用低阶公式从而导致精度较低的问题。本专利技术的方法算法简单，星载软件容易实现，运算时间较短且计算出的磁场强度精度较高。附图说明图1为本专利技术的地磁场查表线性插值计算示意图。具体实施方式下面将对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。一种基于查表法确定空间磁场强度的方法，该方法包含：(1)建立高精度磁场表：在卫星飞行轨迹经过的球面上按照一定的时间间隔进行网格划分，根据国际地磁场模型(IGRF)最高阶公式计算出网格点上的三轴磁场强度，形成地磁场数据表，将其装载在卫星星载软件中；(2)修正卫星飞行轨迹降交点地理经度间隔Gdig以及即将飞过的降交点地理经度lamda；(3)在星箭分离后，星载软件首次调用磁场表子程序时，已经飞过的降交点地理经度ccl由即将飞过的降交点地理经度lamda赋值；其它情况下，当卫星飞到降交点时更新已经飞过的降交点地理经度ccl；(4)当已经飞过的降交点地理经度ccl发生变化时，计算卫星当前位置在磁场表中所对应磁场轨道的编号cca及其降交点地理经度，否则保持cca不变；其中，降交点地理经度根据步骤(3)获得，其中，磁场轨道的编号cca与网格划分相对应；(5)根据卫星当前位置在磁场表中所对应磁场轨道的编号cca查找出卫星当前位置四周的4个网格点；如图1所示，为本专利技术的地磁场查表线性插值计算示意图，卫星当前位置四周的4个网格点为A、B、C、D；(6)基于线性插值的方法计算卫星当前位置在轨道坐标系下的磁场强度。进一步地，在步骤(1)中，网格划分的方式为：横向以降交点(即卫星由北向南运行时，其轨道面与赤道面的交点)地理经度等间隔分成n本文档来自技高网...

【技术保护点】
1.一种基于查表法确定空间磁场强度的方法，其特征在于，该方法包含：(1)建立高精度磁场表：在卫星飞行轨迹经过的球面上按照一定的时间间隔进行网格划分，根据国际地磁场模型最高阶公式计算出网格点上的三轴磁场强度，形成地磁场数据表，将其装载在卫星星载软件中；(2)修正卫星飞行轨迹降交点地理经度间隔Gdjg以及即将飞过的降交点地理经度lamda；(3)在星箭分离后，星载软件首次调用磁场表子程序时，已经飞过的降交点地理经度ccl由即将飞过的降交点地理经度lamda赋值；其它情况下，当卫星飞到降交点时更新已经飞过的降交点地理经度ccl；(4)当已经飞过的降交点地理经度ccl发生变化时，计算卫星当前位置在磁场表中所对应磁场轨道的编号cca及其降交点地理经度，否则保持cca不变；其中，降交点地理经度根据步骤(3)获得，磁场轨道的编号cca与网格划分相对应；(5)根据卫星当前位置在磁场表中所对应磁场轨道的编号cca查找出卫星当前位置四周的4个网格点；(6)基于线性插值的方法计算卫星当前位置在轨道坐标系下的磁场强度。

【技术特征摘要】
1.一种基于查表法确定空间磁场强度的方法，其特征在于，该方法包含：(1)建立高精度磁场表：在卫星飞行轨迹经过的球面上按照一定的时间间隔进行网格划分，根据国际地磁场模型最高阶公式计算出网格点上的三轴磁场强度，形成地磁场数据表，将其装载在卫星星载软件中；(2)修正卫星飞行轨迹降交点地理经度间隔Gdjg以及即将飞过的降交点地理经度lamda；(3)在星箭分离后，星载软件首次调用磁场表子程序时，已经飞过的降交点地理经度ccl由即将飞过的降交点地理经度lamda赋值；其它情况下，当卫星飞到降交点时更新已经飞过的降交点地理经度ccl；(4)当已经飞过的降交点地理经度ccl发生变化时，计算卫星当前位置在磁场表中所对应磁场轨道的编号cca及其降交点地理经度，否则保持cca不变；其中，降交点地理经度根据步骤(3)获得，磁场轨道的编号cca与网格划分相对应；(5)根据卫星当前位置在磁场表中所对应磁场轨道的编号cca查找出卫星当前位置四周的4个网格点；(6)基于线性插值的方法计算卫星当前位置在轨道坐标系下的磁场强度。2.根据权利要求1所述的基于查表法确定空间磁场强度的方法，其特征在于，在步骤(1)中，所述网格划分的方式为：横向以降交点地理经度等间隔分成n轨，n为卫星每天飞行圈数向下取整；纵向以卫星飞行的一轨时间等间隔取m个点，作为轨道钟。3.根据权利要求2所述的基于查表法确定空间磁场强度的方法，其特征在于，横向每两轨之间的地理经度间隔ΔL为360°/n，地理经度与cca的关系为：地理经度＝(cca-1)*360°/n；纵向中，一轨时间等间隔取m为70～80。4.根据权利要求2所述的基于查表法确定空间磁场强度的方法，其特征在于，所述网格划分的数据点从东经0°的降交点时刻开始按顺序排列。5.根据权利要求1所述的基于查表法确定空间磁场强度的方法，其特征在于，在步骤(2)中，根据卫星实际运行的轨...

【专利技术属性】
技术研发人员：李利亮，朱文山，李圣文，牛睿，刘川，
申请(专利权)人：上海航天控制技术研究所，
类型：发明
国别省市：上海,31