本发明涉及基于角速度跟踪的变轨姿态保持控制方法，具体涉及一种大推力变轨期间利用陀螺信息的姿态跟踪控制方法。
背景技术：
由于运载能力受限，无法直接将卫星送入预定轨道。尤其是高轨卫星一般需要依靠自身变轨进入目标轨道。转移段变轨控制是一般高轨卫星一个过程，变轨期间的姿态控制是实现成功变轨的基础。
目前国内常用的方法是采用地球敏感器和太阳敏感器分别确定地球和太阳的指向，联合确定在变轨期间卫星的姿态指向。该方法突出的矛盾是其一对太阳、地球、卫星三者的关系约束强，导致卫星的发射窗口大大受限；其二，部分卫星变轨姿态下使用的地球敏感器在进入目标轨道后的正常姿态下地球是不可见的，引起的很大的资源浪费。
技术实现要素：
本发明的目的是提供一种基于角速度跟踪的变轨姿态保持控制方法，通过充分利用星上资源，解决高轨卫星变轨期间的姿态控制问题，实现高精度高可靠控制。
为了实现以上目的，本发明通过以下技术方案实现：
一种基于角速度跟踪的变轨姿态保持控制方法，包含以下过程：根据变轨策略地面仿真模拟变轨过程，计算整个变轨过程中卫星x、y与z三轴姿态角速度。通过二次曲线拟合理论计算变轨期间的三轴惯性姿态角速度，并上注二次曲线系数。在实施过程中，根据星敏和陀螺信息建立初始变轨指向姿态。在变轨期间利用陀螺信息连续跟踪理论变轨期间理论惯性姿态。
优选地，所述计算整个变轨过程中卫星x、y与z三轴姿态惯性角速度进一步为地面根据精度需求建立推力器模型和轨道动力学模型，并根据变轨策略确定的变轨开始时间、时间长度仿真模拟变轨过程计算得到。
优选地，根据上述整个变轨过程中卫星x、y与z三轴姿态惯性角速度，
采用二次曲线拟合依次计算变轨期间的理论三轴惯性姿态角速度，拟合算法如下：需要拟合的参数y，需要拟合的时间为[t0,tf]，拟合的原始数据记为yi(ti∈[t0,tf])。
二次多项式拟合的系数阵
其中，
得到三轴惯性姿态角速度系数
优选地，根据上述中拟合的二次曲线系数，在变轨期间利用陀螺信息连续跟踪理论变轨期间理论惯性姿态进一步包含以下过程：
采集陀螺数据，进行数据有效性诊断；在数据有效的情况下，将陀螺测量到角速度扣除跟踪角速度作为角速度偏差，利用该角速度偏差的积分作为姿态偏差，进行卫星姿态控制；其中跟踪角速度计算如下：
当前时间记为t，则三轴的跟踪角速度分别为x轴姿态角速度为
本发明与现有技术相比具有以下优点：
本发明利用陀螺信息，实现变轨姿态的连续跟踪，无需对使用条件和发射窗口进行约束，同时不需要变轨期间轨道预报进行预报，有助于提升卫星平台控制能力。对太阳光照、敏感器视场均无约束，具有精度高，适应性强的特点。
附图说明
图1为本发明基于角速度跟踪的变轨姿态保持控制方法的流程图。
具体实施方式
以下结合附图，通过详细说明一个较佳的具体实施例，对本发明做进一步阐述。
如图1所示，本发明一种基于角速度跟踪的变轨姿态保持控制方法，包含以下过程：
s1、根据变轨策略地面仿真模拟变轨过程，计算整个变轨过程中卫星三轴姿态角速度。
s2、二次曲线拟合理论变轨期间的三轴惯性姿态角速度，并上注二次曲线系数。
s3、在实施过程中根据星敏和陀螺信息建立初始变轨指向姿态。
s4、在变轨期间利用陀螺信息连续跟踪理论变轨期间理论惯性姿态。
关于步骤s1进一步包含以下过程，地面根据精度需求建立推力器模型和轨道动力学模型，并根据变轨策略确定的变轨开始时间、时间长度仿真模拟变轨过程，计算整个变轨过程中卫星三轴姿态惯性角速度。
关于步骤s2进一步包含以下过程，利用二次曲线拟合理论变轨期间的三轴惯性角速度，并上注二次曲线系数；其中，
二次多项式拟合的算法：
需要拟合的参数y，需要拟合的时间为[t0,tf]，拟合的原始数据记为yi(ti∈[t0,tf])。
对二次多项式拟合的系数阵
式中，
依次对变轨开始时刻到变轨结束时刻的卫星x、y、z三轴惯性姿态角速度拟合，得到二次曲线系数
关于步骤s3进一步包含以下过程，在实施过程中变轨开始时刻前，星上姿轨控分系统完成根据星敏和陀螺信息建立初始变轨指向姿态。
关于步骤s4进一步包含以下过程，在变轨期间利用陀螺信息连续跟踪理论变轨期间理论惯性姿态。采集陀螺数据，进行数据有效性诊断。在数据有效的情况下，将陀螺测量到角速度扣除跟踪角速度作为角速度偏差，利用该角速度偏差的积分作为姿态偏差，进行卫星姿态控制。其中跟踪角速度计算如下：
当前时间记为t，则卫星x、y、z三轴的跟踪角速度分别为x轴姿态角速度为
综上所述，本发明利用全陀螺信息实现变轨，无需对使用条件和发射窗口进行约束，同时不需要变轨期间轨道预报进行预报，有助于提升卫星平台控制能力。
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。