基于微分修正的太阳同步冻结轨道双脉冲推力标定方法技术

本发明专利技术公开的基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，首先将两脉冲速度增量反解为解析解，然后采用微分修正的方法以半长轴实际改变量为目标函数，对两脉冲控制量进行迭代寻优，从而得到实际速度增量，再以理论速度增量和反解的实际速度增量求比的形式，给出两脉冲每次变轨时推力器效率的评估结果，同时给出第一脉冲控后

【技术实现步骤摘要】
基于微分修正的太阳同步冻结轨道双脉冲推力标定方法


[0001]本专利技术属于航天测量与控制方法
，具体涉及基于微分修正的太阳同步冻结轨道双脉冲推力标定方法
。

技术介绍

[0002]太阳同步冻结轨道，是一种通过约束半场轴
a、
偏心率
e
和近地点幅角
ω
，实现航天器经过同一纬度圈时具有相同高度的轨道，这种轨道特性可以确保航天器上的载荷具有稳定的对地条件，因此在航天器系统中得到了广泛的应用
。
该类型的轨道捕获控制一般需要进行多批次的同轨双脉冲控制，即在一个轨道周期的时间内进行两次变轨
。
根据高斯轨道动力学方程可知，按照偏心率和半长轴改变量的关系，会采用同向或反向双脉冲控制两种形式
。
对于同向双脉冲变轨，第1脉冲位置可以指定，位置选定后速度增量的大小以及第二个速度增量的大小和位置就完全确定了；对于反向双脉冲变轨，两次变轨位置由偏心率矢量改变量决定，两次速度增量也是确定的，且两脉冲执行位置差
180
°
。
这种典型的变轨形式，由于测控条件的绝对约束，不能保证对两次变轨均进行实时监视或者对单次机动分别进行定轨，这就给轨控效果和同类型推力器的效率标定工作带来了一定的困难
。
[0003]对于肼燃料推进系统的航天器，其主要工作特点是落压式
、
推力渐变，所产生的推力是以贮箱压力为自变量的多项式函数，而多项式函数的系数为装订常数，在地面模拟真空状态下通过实验测试得到，而且在航天器的整个寿命期内保持不变
。
考虑航天器姿态控制误差
、
点火时长控制误差
、
测控系统定轨精度误差等因素的影响，随着轨控次数的增加，贮箱压力会不断降低，而仍以该多项式函数表示的推力来计算轨道机动点火时长，产生的实际执行速度增量与理论所需速度增量便会出现一定的误差，从而导致轨道控制精度降低，因此需要考虑综合误差因素对推力器效率进行在轨标定
、
对其迭代修正，从而完成高精度的航天器轨道控制任务
。
[0004]对于使用同类型推力器的同轨双脉冲轨道控制任务，由于很难满足两脉冲中间的测定轨条件，因此不能采用传统的基于速度增量或者半长轴变化的方法对推力效率进行标定，只能对两次机动进行联合标定，但这样也会遇到无法标定的情况，如第1脉冲的速度增量为
‑
0.092m/s
，第2脉冲的速度增量为
0.097m/s
，两次合成的速度增量为
0.005
米
/
秒，两次机动控后偏心率
e
和近地点幅角
ω
会产生一定的改变，但半长轴的改变量仅为米级，两次机动的误差可能会相互抵消也可能相互叠加，因此无法评估推力效率
。
当然，也可以通过对两次脉冲控制中间增加测控资源的方法以进行定轨，但这样会增长轨控完成总时间，同轨一个轨道周期内的两次控制会被人为的增加两到三个轨道周期，造成了为提高轨控精度而多占用测控资源的矛盾
。
[0005]目前，根据公开的文献和专利资料，针对同类型推力器的推力器标定
、
轨控效果评估方法的研究主要包括：高永等在论文
《
卫星单元肼推进系统轨控过程不确定性研究
》
中采用数值方法对肼推进系统轨控过程进行了仿真研究，将贮箱温度
、
压力对轨控效果的敏感性进行了定量分析；唐歌实等在论文
《
嫦娥一号卫星轨控标定方法研究与实现
》
中针对“嫦
娥一号”卫星提出了一种综合利用控前控后精密定轨
、
轨控过程遥测姿态数据
、
遥测加速度计测量数据对沉底发动机
、
轨控发动机进行标定的方法；
Joan
和
Farheen
等在论文
《Thruster
‑
Specific Force Estimation and Trending of Cassini Hydrazine Thrusters at Saturn》、《Cassini Thruster Calibration Algorithm Using Reaction Wheel Biasing Data》
中研究了土星探测器
Cassini
号的在轨数据，提出了利用反作用轮偏置数据标定推力器的方法；张莹等在论文
《
基于夏氏最小二乘的轨道控制力系数辨识
》
中提出了一种基于夏氏最小二乘法对航天器轨道机动过程中推力系数进行辨识的方法，可根据在轨工作实际对推力系数进行优化预测
。
上述推力标定方法，均需要基于轨控过程中的遥测姿态数据
、
反作用轮偏置数据和控前控后精密定轨结果进行推力估计，且只能用于单次轨控时的标定，不能解决同轨两次控制的推力分别标定问题
。
杨维廉等在论文
《
冻结轨道及其应用
》、《
月球卫星的冻结轨道
》、《
火星卫星的冻结轨道研究
》
中对地月火星卫星的冻结轨道特性
、
轨道控制方法进行了深入的研究；杨嘉墀等在图书
《
航天器轨道动力学与控制
》
中介绍了太阳同步冻结轨道的控制方法，脉冲式变轨只是瞬间改变其速度矢量，位置矢量不发生改变；完备等在专利
《
编队任务多脉冲控制条件下的推力器在轨自主标定方法
》
中研究了基于星间相对导航数据和牛顿迭代方法实现编队任务两脉冲控制条件下的推力器在轨自主标定技术，在
InSAR
绕飞卫星任务中实现了星上推力器效率的自主标定，提高了编队任务自主控制的精度；孙守明等在专利
《
一种偏心率冻结同轨双脉冲控制的中间轨道确定方法
》
中研究了基于双脉冲控前
、
控后轨道反演双脉冲速度增量的解析模型，同时给出了中间轨道的确定方法，但由于是解析解，未考虑两次轨控之间存在三体摄动
、
日月引力对轨道的影响，其速度增量的解析解精度较低
。
[0006]上述文献表明，太阳同步回归冻结轨道的轨道捕获控制为公开的
、
已知技术，而推力器性能的在轨标定问题是提高轨道控制精度的决定性问题，具有重要的研究意义和价值
。

技术实现思路

[0007]本专利技术的目的在于提供基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，能够提高航天器轨道控制精度
。
[0008]本专利技术所采用的技术方案是：基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，首先将双脉冲速度增量反解为解析解，然后采用微分修正的方法以半长轴实际改变量为目标函数，对双脉冲控制量进行本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，其特征在于，首先将双脉冲速度增量反解为解析解，然后采用微分修正的方法以半长轴实际改变量为目标函数，对双脉冲控制量进行迭代寻优，得到实际速度增量，再以理论速度增量和反解的实际速度增量求比的形式，给出双脉冲每次变轨时推力器效率的评估结果，同时给出第1脉冲控后
、
第2脉冲控前的轨道参数
。2.
如权利要求1所述的基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，其特征在于，包括以下步骤：步骤
1、
将第1脉冲点火点
J2000
地心惯性系位置速度矢量第2脉冲熄火点
J2000
地心惯性系位置速度矢量分别转换为
J2000
惯性系瞬根和步骤
2、
根据步骤1所得第1脉冲点火点
J2000
地心惯性系瞬时轨道参数计算第1脉冲点火点相位角步骤
3、
根据步骤1所得第2脉冲熄火点
J2000
地心惯性系瞬时轨道参数计算第2脉冲熄火点相位角步骤
4、
根据第1脉冲点火点
、
第2脉冲熄火点
J2000
地心惯性系轨道参数平根计算偏心率矢量改变量
Δ
e
x
、
Δ
e
y
和轨道平半长轴变化量
Δ
a
；步骤
5、
根据第1脉冲点火点
J2000
地心惯性系速度矢量步骤1所得第1脉冲点火点
J2000
地心惯性系瞬时轨道参数步骤2所得第1脉冲点火点相位角及步骤4所得偏心率矢量改变量
Δ
e
x
、
Δ
e
y
和轨道平半长轴变化量
Δ
a
，反解两脉冲轨道坐标系下实际速度增量的参考解析解步骤
6、
在步骤2所得第1脉冲点火点相位角位置施加步骤5所得轨道坐标系下的参考速度增量得到第1脉冲熄火点瞬时参考轨道步骤
7、
将步骤6所得第1脉冲熄火点的瞬时参考轨道进行轨道外推至第2脉冲点火点相位角位置，由于脉冲式变轨是瞬时完成，即第2脉冲点火点相位角与步骤3所得第2脉冲熄火点相位角相同，得到第2脉冲点火点瞬时参考轨道步骤
8、
基于步骤7所得第2脉冲点火点瞬时参考轨道在第2脉冲点火点相位角位置施加步骤5所得轨道坐标系下的参考速度增量得到第2脉冲熄火点瞬时参考轨道步骤
9、
计算步骤8所得第2脉冲熄火点瞬时参考轨道与步骤1所得实际熄火点瞬时轨道的半长轴差
Δ
a
′
；步骤
10、
采用微分迭代的方法，修正以步骤9计算的
Δ
a
′
最小为目标函数，在两脉冲点火位置施加修正后的速度增量，重复步骤6～
10
，得到修正后的双脉冲速度增量
步骤
11、
根据步骤
10
所得双脉冲速度增量计算本次同轨双脉冲控制的推力器标定系数和两次控制推力效率
k1、k2，更新下一次双脉冲控制即将使用的推力器标定系数；步骤
12、
根据步骤
11
所得推力器标定系数，对下一次同轨双脉冲轨道控制的点火时长
T1、T2进行修正
。3.
如权利要求2所述的基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，其特征在于，所述步骤1中，中，式中，
f
rv2orb
(
·
,
·
)
是将
J2000
惯性系位置
、
速度矢量转换为
Kepler
瞬时根数的函数；分别为第1脉冲点火点的瞬时半长轴
、
偏心率
、
倾角
、
升交点赤经
、
近地点幅角
、
平近点角，分别为第2脉冲熄火点的半长轴
、
偏心率
、
倾角
、
升交点赤经
、
近地点幅角
、
平近点角
。4.
如权利要求3所述的基于微分修正的太阳同步冻结轨道双脉冲推力标定方法，其特征在于，所述步骤2中，中，且式中，为第1脉冲点火点真近点角；
f
TruA
(
·
,
·
)
为计算真近点角的函数；
fmod(
·
,
·
)
为求余的函数；为第1脉冲熄火点相位角；所述步骤3中，中...

【专利技术属性】
技术研发人员：马宏，伍升钢，孙守明，曹静，任登高，邢楠，李超，叶修松，杨彪，高云鹏，
申请(专利权)人：中国西安卫星测控中心，
类型：发明
国别省市：