本发明公开了一种用星敏感器星点位置变化估算星体角速度的方法，包括步骤：计算星图转角α；计算星图平移量(Δx,Δy)；根据星图转角α和星图平移量(Δx,Δy)计算惯性角速度。本发明还提供相应的系统。本发明可以在敏感器完成星图识别前提供角速度信息，解决了以往算法在没有连续有效姿态四元数的条件下不能完成角速度输出功能，大幅度放宽了限制条件。角速度的算法比较简单，总体上计算量很小，对于星敏感器处理器速度，各项功能的正常完成影响很小。在控制系统实际使用方面，可以在飞行器发射入轨阶段，能够快速、并可靠地输出角速度信息，用于速率阻尼、姿态快速机动等各种控制模式，从提高系统可靠性等工程应用方面，有着很大实际意义。
1.一种用星敏感器星点位置变化估算星体角速度的方法，其特征在于，包括如下步骤：步骤1：计算星图转角α；
步骤3：根据所述星图转角α和星图平移量(Δx,Δy)计算惯性角速度；
第i和第j个星点间相对位置矢量 的计算公式为：式中：(xi(k),yi(k))为k时刻第i个星点的标定位置，i＝1,2,...,n，(xj(k),yj(k))为k时刻第j个星点的标定位置，j＝1,2,...,n；
k时刻的星点相对位置矢量 和k+1时刻的星点相对位置矢量 之间的夹角αij(k+1)的计算公式如下：得到 个角度值，对这 个角度值进行平均处理，得到第k时刻到k+1时刻星图的转动角度α(k+1)的计算公式如下：所述步骤2，具体为，
第k时刻各星点在空间中的位置(x′i(k)，y′i(k))计算公式如下：x′z(k)＝xi(k)·cosα(k+1)-yi(k)·sinα(k+1)y′i(k)＝xi(k)·sinα(k+1)+yi(k)·cosα(k+1)星点位置平移量计算公式如下：
Δxi(k+1)＝xi(k+1)-x′i(k)Δyi(k+1)＝yi(k+1)-y′i(k)对n个平移值求平均，得到k+1时刻星图的位置平移为：所述步骤3，具体为，
惯性角速度wi(k+1)计算公式如下：其中，Tb为第k到k+1之间的时间间隔。
2.一种用星敏感器星点位置变化估算星体角速度的系统，其特征在于，包括如下装置：第一计算装置，用于计算星图转角α；
第三计算装置，用于根据所述星图转角α和星图平移量(Δx,Δy)计算惯性角速度；
第i和第j个星点间相对位置矢量 的计算公式为：式中：(xi(k),yi(k))为k时刻第i个星点的标定位置，i＝1,2,...,n，为k时刻第j个星点的标定位置，j＝1,2,...,n；
k时刻的星点相对位置矢量 和k+1时刻的星点相对位置矢量 之间的夹角αij(k+1)的计算公式如下：得到 个角度值，对这 个角度值进行平均处理，得到第k时刻到k+1时刻星图的转动角度α(k+1)的计算公式如下：所述第二计算装置，具体为，
第k时刻各星点在空间中的位置(x′i(k)，y′i(k))计算公式如下：x′i(k)＝xi(k)·cosα(k+1)-yi(k)·sinα(k+1)y′i(k)＝xi(k)·sinα(k+1)+yi(k)·cosα(k+1)星点位置平移量计算公式如下：
Δxi(k+1)＝xi(k+1)-x′i(k)Δyi(k+1)＝yi(k+1)-y′i(k)对n个平移值求平均，得到k+1时刻星图的位置平移为：所述第三计算装置，具体为，
惯性角速度wi(k+1)计算公式如下：其中，Tb为第k到k+1之间的时间间隔。
用星敏感器星点位置变化估算星体角速度的方法及系统
技术领域
[0001] 本发明涉及一种利用星敏感器星点位置变化估算星体角速度的技术。
背景技术
[0002] 星敏感器姿态测量精度高、产品性能可靠，在空间飞行器上的应用广泛。部分星敏感器产品根据前后拍的姿态的变化估算并输出惯性角速度，完成这种方法的条件是星敏感器需要运行在正常跟踪、并姿态确定数据有连续有效。
发明内容
[0003] 针对现有技术存在的不足，本发明解决的技术问题是提出一种采用星敏感器中星点位置变化来估算角速度的方法，通过星点位置变化与角度变化之间的关系估算角速度，可消除星敏感器须有连续姿态确定数据的比较苛刻的条件，以改善产品使用适应性。
[0004] 为解决上述技术问题，本发明是通过以下的技术方案实现的。
[0005] 根据本发明的一个方面，提供一种用星敏感器星点位置变化估算星体角速度的方法，包括如下步骤：
[0006] 步骤1：计算星图转角α；
[0007] 步骤2：计算星图平移量(Δx,Δy)；
[0008] 步骤3：根据所述星图转角α和星图平移量(Δx,Δy)计算惯性角速度。
[0009] 优选地，所述步骤1，具体为，
[0010] 第i和第j个星点间相对位置矢量 的计算公式为：
[0011]
[0012] 式中：(xi(k),yi(k))为k时 刻第i个星点的标定位 置，i＝1,2,...,n，yi(k)-yj(k))为k时刻第j个星点的标定位置，j＝1,2,...,n；
[0013] k时刻的星点相对位置矢量 和k+1时刻的星点相对位置矢量 之间的夹角αij(k+1)的计算公式如下：
[0014]
[0015] 得到 个角度值，对这 个角度值进行平均处理，得到第k时刻到k+1时刻星图的转动角度α(k+1)的计算公式如下：
[0016]
[0017] 优选地，所述步骤2，具体为，
[0018] 第k时刻各星点在空间中的位置(xi′(k)，yi′(k))计算公式如下：
[0019] xi′(k)＝xi(k)·cosα(k+1)-yi(k)·sinα(k+1)
[0020] yi′(k)＝xi(k)·sinα(k+1)+yi(k)·cosα(k+1)
[0021] 星点位置平移量计算公式如下：
[0022] Δxi(k+1)＝xi(k+1)-xi′(k)
[0023] Δyi(k+1)＝yi(k+1)-yi′(k)
[0024] 对n个平移值求平均，得到k+1时刻星图的位置平移为：
[0025]
[0026]
[0027] 优选地，所述步骤3，具体为，
[0028] 惯性角速度wi(k+1)计算公式如下：
[0029]
[0030] 其中，Tb为第k到k+1之间的时间间隔。
[0031] 根据本发明的另一个方面，还提供一种用星敏感器星点位置变化估算星体角速度的系统，包括如下装置：
[0032] 第一计算装置，用于计算星图转角α；
[0033] 第二计算装置，用于计算星图平移量(Δx,Δy)；
[0034] 第三计算装置，用于根据所述星图转角α和星图平移量(Δx,Δy)计算惯性角速度。
[0035] 优选地，所述第一计算装置，具体为，
[0036] 第i和第j个星点间相对位置矢量 的计算公式为：
[0037]
[0038] 式中：(xi(k),yi(k))为k时 刻第i个星点的标定位 置，i＝1,2,...,n，yi(k)-yj(k))为k时刻第j个星点的标定位置，j＝1,2,...,n；
[0039] k时刻的星点相对位置矢量 和k+1时刻的星点相对位置矢量 之间的夹角αij(k+1)的计算公式如下：
[0040]
[0041] 得到 个角度值，对这 个角度值进行平均处理，得到第k时刻到k+1时刻星图的转动角度α(k+1)的计算公式如下：
[0042]
[0043] 优选地，所述第二计算装置，具体为，
[0044] 第k时刻各星点在空间中的位置(xi′(k)，yi′(k))计算公式如下：
[0045] xi′(k)＝xi(k)·cosα(k+1)-yi(k)·sinα(k+1)
[0046] yi′(k)＝xi(k)·sinα(k+1)+yi(k)·cosα(k+1)
[0047] 星点位置平移量计算公式如下：
[0048] Δxi(k+1)＝xi(k+1)-xi′(k)
[0049] Δyi(k+1)＝yi(k+1)-yi′(k)
[0050] 对n个平移值求平均，得到k+1时刻星图的位置平移为：
[0051]
[0052]
[0053] 优选地，所述第三计算装置，具体为，
[0054] 惯性角速度wi(k+1)计算公式如下：
[0055]
[0056] 其中，Tb为第k到k+1之间的时间间隔。
[0057] 星敏感器采用星点位移算法计算惯性角速度的方法，可以在敏感器完成星图识别前提供角速度信息，解决了以往算法在没有连续有效姿态四元数的条件下不能完成角速度输出功能，大幅度放宽了限制条件。
[0058] 角速度的算法比较简单，比通常的角速度算法计算量略为增加，但总体上计算量还是很小，对于星敏感器处理器速度，各项功能的正常完成影响很小。
[0059] 在控制系统实际使用方面，可以在飞行器发射入轨阶段，能够快速、并可靠地输出角速度信息，用于速率阻尼、姿态快速机动等各种控制模式，对提高系统可靠性等工程应用方面，有着很大实际意义。
附图说明
[0060] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
[0061] 图1为星点位置变化示意图。
具体实施方式
[0062] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
[0063] 1、星图转角α的计算
[0064] 星敏感器CCD上像点位置变化，是由于敏感器坐标系绕某个空间轴转动引起。根据星点相对位置关系保持固定，敏感器坐标系绕某空间轴旋转一个角度，敏感器CCD测得的星点图形形状不变。星点标定位置变化可以分解为绕Zs轴转角α，与平面内的位置平动(Δx,Δy)二部分，如图1所示，设某星点开始位置(x1,y1)，绕Os转动质心位置到(x′,y′)，再通过平移到(x2,y2)。
[0065] 根据敏感器坐标系运动关系，星图平移量Δx、Δy和转动角度α，是以3、1、2转序下得到的欧拉转角，分别对应于俯仰角、滚动角和偏航角。
[0066] 图1中可以看出，线段A1B1绕Os转动α得到A′B′，因此只要计算线段A1B1与A′B′间的夹角就可得到转角α。点A1和B1为移动前的位置，敏感器可以测量得到，而点A′和B′只是星点位置移动前后理论分解的坐标，不能直接测量得到。
[0067] 由于A′B′通过平移得到A2B2，A′B′与A2B2平行，因此A1B1与A′B′间的夹角，等于A1B1与A2B2之间的夹角，而点A2、B2位置可以测量得到，因此通过计算星点相对位置矢量A1B1和A2B2的夹角得到转动角度。
[0068] 通常情况下敏感器中有多个恒星点，当经过处理判断敏感器视场中有效的星点数目为n，k时刻测量得到的标定位置分别为(xi(k),yi(k))，i＝1,2,...,n。通过矢量计算，第i和第j个星点间相对位置矢量为
[0069]
[0070] 其中，i＝1,2,...,n-1，j＝i+1,2,...,n。通过(1)式计算k和k+1时刻的星点相对位置矢量 和 两个矢量之间的夹角可以用余弦定理计算得到：
[0071]
[0072] 可以得到 个角度值，对这些值的进行平均处理以提高估算精度，得到第k时刻到k+1时刻星图的转动角度α(k+1)：
[0073]
[0074] 2、星图平移量(Δx，Δy)的计算
[0075] 第k时刻的n个星点标定位置(xi(k)，yi(k))，i＝1,2,...,n，因星图绕光轴Os转动角度为α(k+1)，由此各星点的标定位置同样绕Os转动α(k+1)角度，得到各星点在空间中的位置(xi′(k)，yi′(k))：
[0076] xi′(k)＝xi(k)·cosα(k+1)-yi(k)·sinα(k+1) (4)
[0077] yi′(k)＝xi(k)·sinα(k+1)+yi(k)·cosα(k+1) (5)
[0078] 敏感器在k+1时刻测得的星点标定位置为(xi(k+1)，yi(k+1))，与第k时刻转动α(k+1)得到各星点在空间中的位置比较，可以得到n个星点位置平移量：
[0079] Δxi(k+1)＝xi(k+1)-xi′(k) (6)
[0080] Δyi(k+1)＝yi(k+1)-yi′(k) (7)
[0081] 对n个平移值求平均，得到k+1时刻星图的位置平移为：
[0082]
[0083]
[0084] 3、惯性角速度的估算
[0085] 根据转动运动，敏感器绕OsXs轴正向转动，像点的纵坐标将变大，Δy为正值；敏感器绕OsYs正向转动，像点的横坐标将变小，Δx为负值；敏感器绕OsZs正向转动，星图将反相转动，α为负值。
[0086] 若星敏感器曝光时刻第k到k+1之间的时间间隔为Tb，敏感器中的星图位置变化分解为星图位置绕Zs轴转角α(k+1)和平移(Δx(k+1),Δy(k+1))。
[0087] 因此相对于敏感器坐标系，飞行器的惯性角速度wi(k+1)可以用下式进行测算：
[0088]
[0089] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。
