一种星敏感器低频误差在轨实时补偿方法
【技术领域】
[0001] 本发明属于航天器姿态确定领域，涉及一种卫星平台的高精度姿态确定。
【背景技术】
[0002] 高精度指向卫星广泛采用星敏感器测量卫星姿态。对于卫星姿态确定系统，在过 去的几十年中，人们通常将星敏感器的测量误差近似为高斯白噪声。在对姿态确定精度要 求不高的场合，这种近似是合理的。近年来，高分辨率对地观测卫星的发展，对卫星姿态确 定精度提出了越来越高的要求。在这种情况下，将星敏感器测量误差近似为高斯白噪声进 行处理的做法不能满足卫星姿态确定的精度要求。
[0003] 事实上，星敏感器的测量误差可分为几个部分，包括随机噪声、常值偏差和低频误 差等。其中，测量噪声包括散粒噪声、探测器噪声和量化噪声等，这类误差属于随机误差。对 于高精度星敏感器而言，测量噪声的量级通常为几角秒。星敏感器常值偏差可视为等效安 装偏差，是星敏感器自身及其在卫星上的安装面结构变形造成的。常值偏差的存在使得星 敏感器的实际安装矩阵偏离其标称值，最终表现为姿态确定的系统误差。系统误差的量级 通常为数十角秒。星敏感器低频误差指的是安装在航天器上的星敏感器受冷热交变的空间 热环境以及视场变化影响而产生的依轨道周期变化的误差，是影响高精度姿态确定的主要 因素之一。在三轴稳定卫星对地定向运行模式下，卫星运动的角速度等于轨道角速度。卫 星在轨运行过程中受太阳照射角度呈现周期性变化趋势，自身结构受热不均匀，光学系统 受其影响，造成周期性的测量误差。此外，星敏感器观测恒星过程中，其光轴所指向的天区 随卫星轨道运动发生周期性变化，导致光学系统误差、标定误差和星表误差随恒星进出视 场而发生变化，也会造成周期性误差。对于高精度星敏感器而言，低频误差的量级可达数角 秒。
[0004] 在现有的卫星星敏感器测量数据使用中，针对随机噪声，通过星敏/陀螺联合滤 波的方法消弱。对安装误差，在卫星发射前，利用专门的检测设备对星敏感器安装位置进行 标定，可以部分补偿安装偏差的影响，受测量仪器精度限制，安装偏差地面标定精度通常在 15〃左右。但卫星发射过程中的振动变形和器件老化等因素还会造成星敏感器安装矩阵的 变化。星敏感器的相对安装偏差可通过在轨标定的方法进行补偿。低频误差不能直接用传 统的卡尔曼滤波进行处理，也很难通过一般的标定方法完全消除。因此，一般对低频误差不 加考虑。不考虑低频误差存在以下不足：
[0005] 1、实际的姿态确定精度不高
[0006] 在目前的控制系统仿真和在轨型号的姿态确定中，不考虑低频误差，只将消除随 机噪声的精度当作卫星的姿态确定精度。这种简化模型不能准确反映卫星的姿态确定精 度，无法准确定量评估姿态确定精度对有效载荷工作性能的影响。
[0007] 2、影响卫星指向精度
[0008] 由于星敏感器低频误差的存在，在指向误差分配中，需要分配一定欲度。对于要求 越来越高的卫星，是提高指向精度的一个重要限制因素。
[0009] 3、影响遥感图像定位精度
[0010] 在地面图像处理时，由于低频误差的存在，遥感图像定位精度受到影响。低轨卫 星，轨道高度600km，5〃的低频误差影响的定位精度约15m。国外已实现优于IOm的定位精 度，其必然减小了星敏感器低频误差影响。
【发明内容】
[0011] 本发明解决的技术问题是：克服现有技术的不足，提供了一种星敏感器低频误差 在轨实时补偿方法，为实现卫星高精度姿态确定提供了保障。
[0012] 本发明的技术解决方案是：
[0013] 一种星敏感器低频误差在轨实时补偿方法，包括如下步骤：
[0014] (1)将星敏感器低频误差表示为傅立叶级数的形式，再通过该傅立叶级数确定相 对惯性坐标系的误差四元数何= # Sq ;误差四元数柯装订值为[0 0 0 1]；
[0015] (2)设置低频误差补偿标志，且该低频误差补偿标志的值由地面遥控控制，低频误 差补偿标志的值为1或0 ;
[0016] 若低频误差补偿标志值为1且星敏感器数据有效，则使用步骤（1)中所述傅立叶 级数计算误差四元数&7，修正星敏感器测量得到的惯性姿态四元数，从而标定星敏感器低 频误差；
[0017] 若低频误差补偿标志值为1且星敏感器数据无效，则使用上一拍星敏感器有效时 计算的经过误差修正的惯性姿态四元数；
[0018] 若低频误差补偿标志为0,则不进行误差四元数计算，直接使用星敏感器测量得到 的惯性姿态四元数；
[0019] (3)根据卫星在轨运行获得的星敏感器测量信息和有效载荷成像信息，获得星敏 感器相对有效载荷姿态基准的误差；
[0020] (4)根据步骤（3)中得到的星敏感器相对有效载荷姿态基准的误差，注入修改星 上软件装订的系数hFEx、Θ Θ 使用修改后的系数，在轨实时计算星敏感器误差四 元玫&7,按照步骤（2)的方式进行实时补偿。
[0021] 所述步骤（1)中星敏感器低频误差表示为傅立叶级数的形式，具体为：
【主权项】
1. 一种星敏感器低频误差在轨实时补偿方法，其特征在于包括如下步骤： (1) 将星敏感器低频误差表示为傅立叶级数的形式，再通过该傅立叶级数确定相对 惯性坐标系的误差四元数= 误差四元数邱装订值为[0 O O 1]； (2) 设置低频误差补偿标志，且该低频误差补偿标志的值由地面遥控控制，低频误差补 偿标志的值为1或0 ; 若低频误差补偿标志值为1且星敏感器数据有效，则使用步骤（1)中所述傅立叶级数 计算误差四元数&7，修正星敏感器测量得到的惯性姿态四元数，从而标定星敏感器低频 误差； 若低频误差补偿标志值为1且星敏感器数据无效，则使用上一拍星敏感器有效时计算 的经过误差修正的惯性姿态四元数； 若低频误差补偿标志为0,则不进行误差四元数计算，直接使用星敏感器测量得到的惯 性姿态四元数； (3) 根据卫星在轨运行获得的星敏感器测量信息和有效载荷成像信息，获得星敏感器 相对有效载荷姿态基准的误差； (4) 根据步骤（3)中得到的星敏感器相对有效载荷姿态基准的误差，注入修改星上软 件装订的系数hFEx、和Θ 使用修改后的系数，在轨实时计算星敏感器误差四元数 &7,按照步骤（2)的方式进行实时补偿。
2. 根据权利要求1所述的一种星敏感器低频误差在轨实时补偿方法，其特 征在于：所述步骤（1)中星敏感器低频误差表示为傅立叶级数的形式，具体为： ?<,+i>mcos(wr) + A,sin(wr)]，其中，α和β分别表示傅立叶级数中余弦和正弦信号的振 m=\ 2ττ 幅，M傅立叶级数项的个数，取正整数，τ = Coije为星敏感器低频误差 的傅立叶级数基频，T为卫星的轨道周期，tAbs为绝对星时； 所述通过该傅立叶级数确定相对惯性坐标系的误差四元数 Sq = [(>>ν Λ/, 具体为： 5 Qx= Φ 9 LFEx； δ qy= Φ θ LFEy; 5 Qz= Φ θ LFEz； - φ ~Sq~ ； 其中，φ = [I cos(T) sin(O cos(2t) sin(2O ... cos(Mt) sin(MO]， ^ LFEx 一 [ α LFExO Q LFExl ^ LFExl Q LFEx2 ^ LFEx2 · · · Q LFExM ^ LFExM^ ? 9 LFEy 一 [ α LFEyO Q LFEyl ^ LFEyl Q LFEy2 ^ LFEy2 · · · Q LFEyM ^ LFEylJ， G LFEz 一 [ α LFEzO Q LFEzl ^ LFEzl Q LFEz2 ^ LFEz2 · · · Q LFEzM ^ LFEzM^ ? Q LFExO'' Q LFEyO ^ Q LFEzO 为误差的常值分量， α LFExm、α LFEym、α LFEzm、^ LFExm、^ LFEym、 β_ζπ分别表示误差的余弦和正弦 信号振幅，初始装订值均为0。
3. 根据权利要求1所述的一种星敏感器低频误差在轨实时补偿方法，其特征在于：所 述的步骤（3)中星敏感器相对有效载荷姿态基准的误差§为：
，其中，Z是多个轨道周期观测获得的参数辨识的样 本序列，
uB，n为有效载荷观测的地标矢量或恒星星光矢量，η = 1，2,…，N
%"为根据卫星与地标矢量的位置关系计算得到 或者根据卫星与恒星之间位置关系计算得到
：第η次的星敏感器测量值；
cos (2 τ n) sin (2 τ η) ... cos (Μ τ η) sin (Μ τ η)], τη= ω LFEtAbs n，tAbs n为星敏感器第η次 测量所对应的绝对星时。
4. 根据权利要求3所述的一种星敏感器低频误差在轨实时补偿方法，其特征在于：u 为根据卫星与地标矢量的位置关系计算得到或者根据卫星与恒星之间位置关系计算得到， 具体为： 若有效载荷观测的为地标，则u,
CIF，n为观测 地标时刻对应的地固系到惯性系的方向余弦阵，λη、αη为所观测地标的地理经度、地理炜 度，
为观测地标时刻对应的卫星位置在惯性系下的表示，Re为地球半径；
若有效载荷观测的为恒星，则 RAn、DECnS所观测恒星的 赤经、赤炜。
【专利摘要】一种星敏感器低频误差在轨实时补偿方法，针对星敏感器受热变形等因素产生的低频误差，以成像有效载荷为姿态基准，通过观测得到星敏感器相对有效载荷的姿态变化，使用傅立叶级数进行拟合，在轨实时补偿星敏感器低频误差。实时补偿低频误差后，可提高卫星姿态确定精度，消除低频误差影响，有利于卫星高精度成像和图像定位。
【IPC分类】G01C21-02, G01C25-00
【公开号】CN104729537
【申请号】CN201510122691
【发明人】汤亮, 陈守磊, 熊凯, 刘一武, 王楷, 胡少春, 郭廷荣
【申请人】北京控制工程研究所
【公开日】2015年6月24日
【申请日】2015年3月19日