一种基于空间转换的星敏感器精度测试方法
【技术领域】
[0001] 本发明涉及一种星敏感器精度测试技术，具体是一种基于空间转换的星敏感器精 度测试方法。
【背景技术】
[0002] 星敏感器作为空间飞行器的重要姿态敏感器，其测量精度直接决定了在轨姿态确 定的精度，对姿态确定精度要求较高的空间飞行器则尤为重要。而实际在轨飞行的星敏感 器由于平台变形、环境等因素的影响，使得在轨对星敏感器的精度评估变得困难。因此，在 地面对星敏感器精度进行评估则是一项至关重要的工作。现有技术主要通过单星定位精 度、光星模拟器以及外场观星对星敏感器精度进行评估。其中，单星定位精度无法评价星敏 感器自身的星图识别与匹配的精度，只能作为一种预估的手段。光星模拟器由于星点的光 谱范围、星等和位置等与真实星空均存在一定的差别，因此也不能真实地对星敏感器精度 进行测试。而外场观星测试中，普遍使用的是将姿态数据进行曲线拟合，并以拟合值作为真 实值对星敏精度进行测试，拟合效果无法评论，星敏真实精度因此较难获得。
【发明内容】
[0003] 针对现有技术中的缺陷，本发明提供了 了一种基于空间转换的星敏感器精度测试 方法，利用外场观星条件下相对稳定的地速条件，为星敏感器建立相对稳定的工作平台，将 动静空间进行相对变换，并利用外场观星条件下的真实星空，从而真实地对星敏感器精度 进行测试。
[0004] 本发明的目的通过以下技术方案来实现:一种基于空间转换的星敏感器精度测试 方法，包括如下步骤：
[0005] S1、根据星敏感器数据Q(tn)计算星敏感器的运动角速度ω，并利用外场观星的数 据，并对所得角速度ω进行平滑，得到星敏感器转动角速度；
[0006] S2、将星敏感器以所得的角速度ω转动的动态拍图输出的四元数数据Q(tn)转换 为相对to时刻静态拍图的四元数数据Q(tn_t〇);
[0007] S3、根据所得的静态拍图的四元数数据Q(tn-t〇)，通过矢量空间转换，分别获得星 敏感器测量坐标系三轴单位矢量VX= (1，0，0)、VY= (0，1，0)、VZ= (0，0，1)在J2000惯性系下 的单位矢量VXj2_(tn)、VYj2QQ()(tn)、Vzj2QQ()(t n);
[0008] S4、分别计算 J2000 惯性系下单位矢量 Vxj2QQ()(tn)、VYJ2QQ()(tn)、Vzj2QQ()(tn)的平均矢 MVxnooa n Υγποοο：、Vz.i2〇so ;
[0009] S5、以平均矢量作为星敏感器to时刻下三轴矢量在J2000系下的真实矢量，计算真 实矢量与平均矢量之间的矢量夹角％(tn)、θγ(?η)、θζ(?η)，并根据矢量夹角θχ(?η)、θγ(?η)、θζ (tn)进行数据分析，对星敏感器精度进行测试。
[0010]其中，所述星敏感器的运动角度ω通过以下公式计算：
[0011] W = 2((V(U? 你,Mi)-[0 0 0 1]7)/Δ/,;. ?
[0012] 式中：Q(tn)为星每j[感器数据，Δ tn = tn+i-tn，ω =( ωχ, ωγ, ωζ,0)，ig)_为四元数相 乘。
[0013] 其中，所述步骤S2中将星敏感器四元数Q(tn)转换为t0时刻的四元数Q(t n-to)的计 算公式为：
[0014] Q{lii-t,)^QUi!)?\() \t"-〇
[0015] 式中，AQ(tn-to)转换矩阵四元数，计算公式为：
[0016] Δ Q(tn-to) =Q( ω (tn-to))
[0017] = [ ω x/ I I ω I I sin( I I ω I I (tn_to)/2) ωγ/ I I ω I I sin( I I ω I I (tn_to)/2) ω z/ I I ω | |sin(I I ω I I (tn-to)/2)cos(I I ω I I (tn-to)/2)]T。
[0018] 其中，所述步骤S3中星敏感器测量坐标系三轴单位矢量在J2000惯性系下的矢量 四元数计算公式如下：
[0019] Q(V、_(iB)) = _s -i0)(8[l000f 00-心)
[0020] 0(^.,^(/,,)) = 0(^-/,,)0(0 10()]7 '?η-?,Λ)
[0021] Q(VVl：t,)1：!α")) = ρα, 1 0.V ?Q !^"
[0022] 以上矢量四元数的前三列即为每个轴的矢量表达式。
[0023] 其中，所述步骤S4中，J2000惯性系下单位矢量Vxj2_(tn)、VYJ2_(tn)、Vzj2QQ()(tn)的 平均矢里·ν.χΗ〇ο.() s Vyi2〇o.() s Vznoro的计算公式如下： η
[0024] -- - 幻 20()0 n
[0025] ---ξυ) YU' n n
[0026] --- mm n
[0027] 其中，所述步骤S5中星敏感器每个轴向量在J2000系下的向量夹角通过以下公式 计算：
[0028] AcmgX(n) =-- arccos(Vvl：；0!1()(/(i) · )
[0029] AangY(n) = arGC〇s(VYI2〇(!0(/") · V, ?2000)
[0030] AcmgZ(n) ^ arccos(V/J:,)OI)(/") · V/l：0〇(l)
[0031]其中，可以通过对Δ angX、Δ angY、Δ angZ的数据统计分析出星敏感器三轴姿态的 测量精度。
[0032] 与现有技术相比，本发明具有以下优点：
[0033] 利用外场观星条件下相对稳定的地速条件以及真实星空，将星敏感器拍图方式进 行动静空间相对变换，并通过计算矢量夹角的方法测试星敏感器的精度，克服了现有技术 无法模拟真实星空等缺点，为地面对星敏感器精度测试提供了有效的方法。
【附图说明】
[0034]图1为本发明实施例一种基于空间转换的星敏感器精度测试方法的流程图。
【具体实施方式】
[0035]下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术 人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术 人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明 的保护范围。
[0036] 如图1所示，本发明实施例提供了一种基于空间转换的星敏感器精度测试方法，包 括如下步骤：
[0037] S1、根据星敏感器数据Q(tn)计算星敏感器的运动角速度ω，并利用外场观星的数 据，并对所得角速度ω进行平滑，得到星敏感器转动角速度；
[0038] 所述星敏感器的运动角度ω通过以下公式计算：
[0039] ^ = 0 0 i]')/^^
[0040] 式中：Q(tn)为星敏感器数据，Δ tn = tn+i_tn，ω=(ωχ，ωγ，ωζ，〇)，發为四元数相 乘。
[0041] S2、将星敏感器以所得的角速度ω转动的动态拍图输出的四元数数据Q(tn)转换 为相对to时刻静态拍图的四元数数据Q(t n-to);计算公式为：
[0042] %，)
[0043] 式中，AQ(tn-to)转换矩阵四元数，计算公式为：
[0044] Δ Q(tn-to) =Q( ω (tn-to))
[0045] = [ ω