专利名称：：一种提高星敏感器数据处理精度的方法
技术领域：
：本发明涉及一种提高星敏感器数据处理精度的方法。
背景技术：
：目前，所考虑的星敏感器单位方向矢量测量模型如下611]TV(1)其中r表示在参考坐标系中的某参考矢量，如由星历表给出的某恒星矢量(表示在惯性坐标系中)；b表示在星敏感器坐标系中的测量矢量；A为星体坐标系相对参考坐标系的姿态矩阵；Sbs为星敏感器的安装矩阵，即由星敏感器坐标系到星体坐标系的坐标转换矩阵；v为与星敏感器内部设计有关的测量随机误差。然而，星敏感器测量误差不仅包括与星敏感器内部设计有关的误差，还包括与星敏感器安装、在轨环境、星体抖动等因素有关的外部因素引起的外部误差。目前，对于星敏感器测量误差的分析和处理主要集中在与星敏感器内部设计有关的基本误差因素分类及精度估算关系上，很少涉及星敏感器外部综合的误差分析和处理。在星敏感器测量过程中，误差产生机理的不确定性，误差表现形式的多样性，误差关系的复杂性以及要求误差处理的准确性和可靠性都给星敏感器测量误差的特性分析及处理带来了很大的难度。因而造成目前星敏感器测量误差处理效果不理想，进而数据处理精度不高。
发明内容本发明的目的在于提供一种提高星敏感器数据处理精度的方法，以实现不但能对星敏感内部设计有关的测量误差进行分析与处理，还能对除星敏感器内部误差以外的其他误差进行分析和处理。本发明的技术方案包括以下步骤(1)建立星敏感器姿态测量误差模型将星敏感器测量误差分为三大类一是星敏感器内部变形、星敏感器安装位置改变等因素引起的低频误差；二是由卫星轨道参数变化、星图变换及卫星在轨运动中的抖动振动等因素引起的中频误差；三是星敏感器的高频误差，即服从零均值的高斯白噪声分布的高频误差；三种误差的矫正方法如下CM=①低频误差包括二个部分(A)安装矩阵误差其表现形式为cos(A5》sin(As)0、—sin(As^)cos(As^)0\001,(B)常值系统误差其表现形式为C=[As2As2As2]T②中频误差，其表现形式为(2)(3)x=Bsin(2Jift+"(4)(4)式中，B为中频误差的幅度，f为误差频率，4为误差的初相位，且该误差体现在带有误差的安装矩阵上；③高频误差，其分布函数为1^p(v)-exp{-(5)6:■。fc(6)2;rcr'、2cr2(2)建立星敏感器姿态领把姿态测量数据中所涉及到的低频误差、中频误差和高频误差叠加到传统的星敏感器单位方向矢量测量模型中，得到重构后的姿态测量模型如下表示cos(+义)sin(As^+0、一sin(+力cos(As1!+力0\00ly(3)结合定姿过程进行误差求解利用扩展卡尔曼滤波算法进行测量误差模型参数和卫星姿态的估计，待估参数取为X二(，ASl，As2，x，02)，其中为姿态参数。而其中针对步骤1中的姿态测量数据误差模型，设计状态方程，假设低频误差的稳态性和中频误差的时变性，则相应的误差模型参数的状态方程为A^00=Ai2(f)=0(7)i=2;z/Bcos(2;z/+《)本发明实现了不但能对星敏感内部设计有关的测量误差进行分析与处理，还能对除星敏感器内部误差以外的其他误差进行分析和处理。具体实施方式实施例(1)建立星敏感器姿态测量误差模型将星敏感器测量误差分为三大类一是星敏感器内部变形、星敏感器安装位置改变等因素引起的低频误差；二是由卫星轨道参数变化、星图变换及卫星在轨运动中的抖动振动等因素引起的中频误差；三是星敏感器的高频误差，即服从零均值的高斯白噪声分布的随机误差；星敏感器测量误差模拟仿真假设星敏感器的光轴z_sl和z_s2在本体坐标系中的坐标分别为丁丄—fcl=[cos320,0,sin320]]x一fe2=[-cos32。，0,-sin32。](2)根据测量误差源的分析，加入如下各类误差①星敏感器高频误差v包括ov=12〃的随机误差；②星敏感器低频误差包括ASl=3〃的安装矩阵误差(乘性噪声)禾PAs2=2〃的常值系统误差(加性噪声)；③中频误差x=Bsin(+4)中，最大振幅B二3〃，安装点处的误差频率f=10—2Hz;(3)星敏感器姿态测量误差的处理按照卫星轨道参数为地球赤道半径&=6.378Xl(fm，卫星轨道长半轴a=Re+0.6452Xl()6m，轨道倾角i=97.938°，升交点赤经Q=77°，近心点角距w=50°，初始时刻平近心点角距M=10°，轨道角速度w。=;利用重构后的星敏感器姿态测量方程和陀螺测量模型，结合误差四元数状态方程和各类模型误差参数的状态方程，利用扩展卡尔曼滤波算法，在卫星姿态确定的过程中，同时求解误差模型参数和卫星姿态参数；陀螺的精度为测量噪声均方差Og=0.05°/h，陀螺相关漂移白噪声od=0.1°/h，相关漂移初始值为d(0)=T°/h，时间相关常数、=、=、二lh;陀螺常值漂移b二[1-11]T°/h，常值漂移白噪声均方差ob=0.03°/h;滤波时长为3500s，星敏感器采样频率为4Hz。利用仿真的星敏感器观测数据，滤波过程中不进行误差参数求解，则对应的卫星姿态估计误差见表1:而利用扩展卡尔曼滤波进行误差参数求解，所得到的定姿结果见表2。表1不进行误差参数求解下的卫星姿态确定结果<table>tableseeoriginaldocumentpage5</column></row><table>表2星敏感器测量误差求解后的定姿结果<table>tableseeoriginaldocumentpage5</column></row><table>从上表可以看出，经过星敏感器测量误差求解后，姿态估计的偏差已经基本消除，且定姿精度提高了15%以上。权利要求一种提高星敏感器数据处理精度的方法，其特征在于，包括以下步骤(一)建立星敏感器姿态测量误差模型将星敏感器测量误差分为三大类A星敏感器内部变形、星敏感器安装位置改变因素引起的低频误差；B卫星轨道参数变化、星图变换及卫星在轨道运动中的抖动振动因素引起的中频误差；C星敏感器的高频误差，即服从零均值的高斯白噪声分布的随机误差；(二)根据测量误差源的分析，加入上述低频误差、中频误差、高频误差进行矫正①低频误差包括Δs1＝3″的安装矩阵误差和Δs2＝2″的常值系统误差；②中频误差χ＝Bsin(2πft+ζ)中，最大振幅B＝3″，安装点处的误差频率f＝10-2Hz；③高频误差v包括σv＝12″的随机误差；(三)星敏感器姿态测量误差的处理(1)按照卫星轨道参数为地球赤道半径Re＝6.378×106m，卫星轨道长半轴a＝Re+0.6452×106m，轨道倾角i＝97.938°，升交点赤经Ω＝77°，近心点角距ω＝50°，初始时刻平近心点角距M＝10°，轨道角速度(2)利用重构后的星敏感器姿态测量方程和陀螺测量模型，结合现有的误差四元数状态方程和所构建的各类模型误差参数的状态方程，利用扩展卡尔曼滤波算法，在卫星姿态确定的过程中，同时求解误差模型参数和卫星姿态参数；其中陀螺的精度为测量噪声均方差σg＝0.05°/h，陀螺相关漂移白噪声σd＝0.1°/h，相关漂移初始值为d(0)＝T°/h，时间相关常数τx＝τy＝τz＝1h；陀螺常值漂移b＝[1-11]T°/h，常值漂移白噪声均方差σb＝0.03°/h；滤波时长为3500s，星敏感器采样频率为4Hz；(3)利用仿真的星敏感器观测数据，滤波过程中利用扩展卡尔曼滤波进行误差参数求解，根据对应的卫星姿态估计误差得出定姿结果。全文摘要本发明公开了一种提高星敏感器数据处理精度的方法。将星敏感器测量误差分为三大类A星敏感器内部变形、星敏感器安装位置改变等因素引起的低频误差；B卫星轨道参数变化、星图变换及卫星在轨运动中的抖动振动等因素引起的中频误差；C星敏感器的高频误差，即服从零均值的高斯白噪声分布的随机误差；经过星敏感器测量误差求解后，姿态估计的偏差已经基本消除，其定姿精度提高了15％以上。文档编号G01C21/18GK101696885SQ20091004470公开日2010年4月21日申请日期2009年11月5日优先权日2009年11月5日发明者刘洋,周海银,杨建坤,段晓君,潘晓刚,王炯琦,矫媛媛,郭雪姣申请人:中国人民解放军国防科学技术大学;