CN102288201A - 用于星敏感器的精度测量方法 - Google Patents

用于星敏感器的精度测量方法 Download PDF

Info

Publication number
CN102288201A
CN102288201A CN 201110189264 CN201110189264A CN102288201A CN 102288201 A CN102288201 A CN 102288201A CN 201110189264 CN201110189264 CN 201110189264 CN 201110189264 A CN201110189264 A CN 201110189264A CN 102288201 A CN102288201 A CN 102288201A
Authority
CN
China
Prior art keywords
coordinate system
under
star
direction vector
star sensor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 201110189264
Other languages
English (en)
Other versions
CN102288201B (zh
Inventor
尤政
邢飞
孙婷
张高飞
李滨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Tianyin star technology limited liability company
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN2011101892642A priority Critical patent/CN102288201B/zh
Priority to US13/386,979 priority patent/US8433515B2/en
Priority to PCT/CN2011/077977 priority patent/WO2013004033A1/zh
Publication of CN102288201A publication Critical patent/CN102288201A/zh
Application granted granted Critical
Publication of CN102288201B publication Critical patent/CN102288201B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Navigation (AREA)

Abstract

本发明公开了一种用于星敏感器的精度测量方法，包括：1)将星敏感器固定在地球上；2)向所述星敏感器输入测试开始时间相对于J2000.0的时间T；3)根据导航星在J2000.0坐标系下的赤纬和赤经以及视运动参数(α′，δ′)确定在J2000.0直角坐标系下的方向矢量；4)将导航星在J2000.0直角坐标系下的方向矢量转换为历元黄道坐标系下的方向矢量；5)将历元黄道坐标系下的方向矢量转变成天球坐标系下的方向矢量(vCRFT)；6)将导航星从天球坐标系下的方向矢量变到地固坐标系下的方向矢量(vTRF)，并基于所述地固坐标系下的方向矢量(vTRF)，获得所述星敏感器的精度。根据本发明的方法，很容易实现该星敏感器的精度测量方法。

Description

用于星敏感器的精度测量方法
技术领域
本发明属于姿态传感器技术领域，尤其涉及一种用于星敏感器的精度测量方法。
背景技术
星敏感器以精度高、功耗低、体积小等优点成为目前航天器最具竞争力的姿态敏感器件。目前，星敏感器的定姿精度已经可以达到10″，某些型号的星敏感器精度甚至可以达到1″水平，高精度是星敏感器得以迅速发展和广泛应用的关键因素。随着星敏感器精度越来越高，对精度测量方法也提出了更高的要求。传统的测试方法主要基于星模拟器及精密转台，需要转台的位置精度比星敏感器的测量精度再高一个数量级，即达到亚角秒的量级水平，这种设备价格昂贵，操作过程复杂。同时，实验室通过转台标定时，以星模拟器作为测量基准，但实现光谱范围、星等和位置精度皆满足要求的全天球星模拟器难度很大，星模拟器与真实星空的导航星还有较大差距，还不能完全模拟真实星空情况，使实验室测试的真实性和准确性难以得到人们的信服。
因此，找到一个易实现的、能满足精度要求的星敏感器精度测量方法就显得十分重要和迫切。
发明内容
本发明旨在至少解决上述技术问题之一。
为此，本发明需要提供一种用于星敏感器的精度测量方法，所述精度测量方法能够很容易实现、解决传统的测试方法操作复杂，需要价格昂贵的精密转台和星模拟器的困扰，同时测量结果较转台式测量方法更具有准确性和真实性，且测试精度能满足星敏感器的要求。
根据本发明的一方面，提供了一种用于星敏感器的精度测量方法，包括如下步骤：1)将星敏感器固定在地球上，且使得星敏感器的主轴指向天顶，所述星敏感器可输入时间参数且存储有导航星表和导航星的视运动参数；2)向所述星敏感器输入测试开始时间相对于J2000.0时刻的当前时刻T；3)根据星敏感器中的导航星在J2000.0坐标系下的赤纬和赤经(α，δ)以及在两个方向上的视运动参数
Figure BDA0000074120090000021
来确定导航星在当前时刻T在J2000.0直角坐标系下的方向矢量；4)将导航星在当前时刻T在J2000.0直角坐标系下的方向矢量转换为历元黄道坐标系下的方向矢量；5)将历元黄道坐标系下的方向矢量转变成当前时刻T下的天球坐标系下的方向矢量(vCRFT)；以及6)根据实际拍摄时刻(T+Δt)将导航星在当前时刻T从天球坐标系下的方向矢量(vCRFT)变到实际拍摄时刻(T+Δt)在地固坐标系下的方向矢量(vTRF)，并基于所述地固坐标系下的方向矢量(vTRF)，获得所述星敏感器的精度。
由此，在本发明的上述精度测量方法中，通过利用地球本身自转的精密性，将星敏感器固连于地球，使星敏感器的主轴正对天顶进行观测，星敏感器随着地球的一起运动(Ω＝7.292115×10-5rad/s)，星敏感器测量值的角度变化与之相对应，而存储在星敏感器星表内的导航星是在J2000.0坐标系(CRFJ2000)下的坐标，由于星敏感器的三轴精度不一致性，其指向精度比滚转精度高一个量级，为保证测量指向精度的准确性和高精度，将星敏感器中导航星的坐标转换到当前测量时刻地固坐标系(TRF)下的坐标，这样就消除了地球滚转轴对指向精度的影响，此时测量星敏感器的输出结果理论上为恒定值，即为星敏感器坐标系相对于地固坐标系的安装矩阵，以此矩阵为基础可以测量出星敏感器主轴在地固坐标系中的变化，进而测量出星敏感器的指向轴精度。
根据本发明的一个实施例，在所述步骤3)中，在所述当前时刻T下，导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)为：
Figure BDA0000074120090000022
根据本发明的一个实施例，在所述步骤4)中，历元黄道坐标系下的方向矢量(vERF)基于所述导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)且将所述J2000.0坐标系绕着所述J2000.0坐标系的X轴逆时针旋转23°26′21″的方向变换之后获得：
Figure BDA0000074120090000023
根据本发明的一个实施例，将导航星在历元黄道坐标系下的方向矢量(vERF)转变成当前时刻T下的天球坐标系下的方向矢量通过下述获得：
将历元黄道坐标下的方向矢量(vERF)绕其Z轴顺时针方向转动
Figure BDA0000074120090000024
接着绕第一次转动后的坐标系的X轴顺时针方向转动
Figure BDA0000074120090000025
接着绕第二次旋转后的坐标系的X轴逆时针方向旋转εA；
接着绕第三次旋转后的坐标系的Z轴顺时针方向旋转
Figure BDA0000074120090000031
以及
接着绕第四次旋转后的坐标系的X轴顺时针方向旋转以获得含有章动项的当前时刻(T)的天球坐标系下的方向矢量(vCRFT)，其中
Figure BDA0000074120090000033
Figure BDA0000074120090000034
分别表示黄经章动和斜章动。
根据本发明的一个实施例，所述导航星在天球坐标系下的方向矢量(vCRFT)通过下述公式获得：
Figure BDA0000074120090000035
Figure BDA0000074120090000036
其中Rx、Rz为绕X轴和Z轴旋转的坐标变换基。
根据本发明的一个实施例，根据IAU2000B章动模型，εA与黄经章动和斜章动
Figure BDA0000074120090000038
分别为：
Figure BDA0000074120090000039
Figure BDA00000741200900000310
Figure BDA00000741200900000311
其中，
Figure BDA00000741200900000313
Figure BDA00000741200900000314
t为从J2000.0开始的儒略世纪数并基于时刻T获得；
幅角αi为幅角的线性组合：
Figure BDA00000741200900000316
式中，nik为整数，Fk为与太阳月亮位置有关的Delaunay幅角。
根据本发明的一个实施例，所述步骤(6)进一步包括：
(61)根据实际拍摄时刻(T+Δt)将导航星矢量从T坐标系转到实际拍摄时刻(T+Δt)地固坐标系下的方向矢量(vTRF)；
(62)根据所述地固坐标系下的方向矢量(vTRF)通过QUEST方法求解星敏感器的最优姿态矩阵
Figure BDA00000741200900000317
以及
(63)计算实际拍摄时刻
Figure BDA00000741200900000318
的星敏感器主轴指向矢量以及
(64)计算实际拍摄时刻
Figure BDA00000741200900000320
的星敏感器主轴指向矢量的夹角(αij)，以获得所述星敏感器的指向精度。
根据本发明的一个实施例，导航星在地固坐标系下的方向矢量(vTRF)通过将所述导航星在天球坐标系下的方向矢量(vCRFT)绕天球坐标系的Z轴以Ω＝7.292115×10-5rad/s逆时针旋转获得：
Figure BDA0000074120090000041
Figure BDA0000074120090000042
根据本发明的一个实施例，所述最优姿态矩阵
Figure BDA0000074120090000043
通过使得下面的目标函数
Figure BDA0000074120090000044
达到最小值而获得：
Figure BDA0000074120090000045
其中，wi，vi分别表示导航星在星敏感器感器坐标系下的方向矢量和在地固坐标系下的方向矢量，αi表示加权系数，满足；
根据本发明的一个实施例，所述星敏感器主轴指向矢量
Figure BDA0000074120090000047
为：
根据本发明的一个实施例，所述星敏感器主轴指向矢量的夹角(αij)为：
其中，iγj。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是恒星在天球球面坐标系和直角坐标系中的坐标矢量示意图；
图2是根据本发明的星敏感器的成像原理图；
图3为地球在天球系统中运动的主要坐标系参数示意图；
图4显示了根据本发明的用于星敏感器的精度测量方法的天球赤道坐标系、历元天球黄道坐标系、地固坐标系和星敏感器坐标系的示意图；
图5显示了根据本发明的用于星敏感器的精度测量方法的流程图；
图6显示了根据本发明的用于星敏感器的精度测量系统的示意图；
图7显示了根据本发明的星敏感器精度测量单元的结构框图；以及
图8显示了根据本发明的用于测量星敏感器的指向精度的示意图。。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
需要说明的是，此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。进一步地，在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
为了详细地阐述本发明的用于星敏感器精度测试的方法和系统，下面将首先介绍根据本发明的一个实施例的星敏感器的工作原理。
星敏感器测量原理
星敏感器姿态通常指的是相对某一指定坐标系的指向，最常用的是采用相对于天球惯性坐标系的指向。星敏感器依靠测量航天器坐标系中导航星的指向来确定星敏感器所在的航天器相对于惯性空间的姿态。在工作状态下，首先测量到导航星在星敏感器坐标系中的矢量，然后通过已经获得的星图来进行识别得到该导航星在惯性坐标系下相对应的矢量。通过比较两个坐标系统中相应导航星的矢量关系，就可以得到从惯性坐标系到航天器坐标系的变换矩阵，即航天器在惯性坐标系中的姿态。
恒星是星敏感器进行工作的参考基准。经过多年大量的天文观测，每颗恒星都在天球1’中具有各自相对固定的位置。图1为恒星在天球球面坐标系和直角坐标系中的坐标矢量示意图。如图1中所示，以天球球面坐标的赤经和赤纬来表示，该恒星的在天球球面坐标系中的坐标可以记作(α，δ)。根据直角坐标与球面坐标的关系，可以得到每颗恒星在天球直角坐标系下的方向矢量为：
Figure BDA0000074120090000061
从星库中选出满足星敏感器成像条件的恒星组成导航星，并由此构成导航星表。根据本发明的一个实施例，该导航星表可以在制造的过程中一次性地固化到星敏感器的存储器中。
当星敏感器1处于天球坐标系中的某一姿态矩阵为A时，利用星敏感器的小孔成像原理，可以通过星敏感器1的镜头2测量得到导航星si(其对应天球坐标系下的方向矢量为vi)在星敏感器坐标系内的方向矢量为wi，如图2中所示。
如图2中所示，星敏感器1的主轴中心在探测器上的位置(x0，y0)，导航星si在星敏感器1的探测器3上的位置坐标为(xi，yi)，星敏感器的焦距为f，则可以得到wi矢量的表达式如下：
Figure BDA0000074120090000062
在理想情况下具有如下关系：
Figure BDA0000074120090000063
其中：A为星敏感器姿态矩阵。
当观测量多于两颗星时，可以直接通过例如QUEST的方法对星敏感器的姿态矩阵A进行求解，即使得下面的目标函数J(Aq)达到最小值来求出最优姿态矩阵Aq：
Figure BDA0000074120090000064
其中，αi表示加权系数，满足；
Figure BDA0000074120090000065
这样，可以计算获得星敏感器在惯性空间中的最优姿态矩阵Aq。
由此可以看出，在真实的星敏感器测量系统中需要高精度的导航星，同时为了保证星敏感器视场的覆盖性，需要转动系统来实现导航星出现在视场的不同位置上，为此传统的标定和测试方法通过单星模拟器以及高精度的转台来实现星点在不同视场下的成像，进而实现系统的标定和测试。为了更加真实和全面的覆盖整个系统，根据本发明的一个实施例，发明人利用了真实星空(导航星表)和地球自转的模式相结合，从而使得用于星敏感器的精度测量更加真实和准确。
下面将详细描述为了地球的运动，以用于根据本发明的星敏感器的高精度测量和分析。
地球的运动规律
本发明的测量方法是将地球的精密运动作为星敏感器的精度测量基准，对于地球在惯性空间的运动需要严格的分析和计算。图3为地球在天球坐标系统中运动的主要坐标系参数。
如图3，以地球为中心作任意半径的一假想大球面称“天球”，地球赤道平面与天球相交的圆称为“天赤道”，地球绕日公转的轨道平面与天球相交的圆称为“黄道”。天赤道与黄道相交于两点，太阳视行从天赤道以南进入天赤道以北与天赤道的交点叫春分点。太阳视行从天赤道以北进入天赤道以南与天赤道的交点叫秋分点。太阳从春分点出发，沿黄道运行一周回到春分点称为一个“回归年”。
如果地轴不改变方向，二分点不动，回归年与恒星年相等。但地轴绕黄极缓慢进动，赤道面与黄道面的交线也以同一周期在黄道面上旋转，如图3所示，天北极以23°26′21″为半径按顺时针方向绕黄北极转动。由于地球的自转方向与地轴的进动方向相反，使春分点每年产生一个微小的西移，天文学上称之为岁差。现代天文学的测量和计算结果表明，地球每年的岁差为50.29″，这样大约25765年北天极绕黄北极旋转一周。
与陀螺的运动模型相似，地球自转轴在进行进动的同时，也在进行着章动，其形成原因较为复杂，笼统的认为是地球附近的其他行星和月亮等对于地球的引力造成的，现代天文学测量结果显示，章动的周期为18.6年(6798天)，在黄道上的黄经章动分量是17.24″，垂直于黄道的斜章动是9.21″。因而使得天体的坐标如赤经、赤纬等都发生变化。
地球的自转轴还存在着极移等现象，但是其周期性的变化都在0.1″以下，因此相对于星敏感器的精度测试可以忽略不计。
地球在惯性空间的运动包括本身围绕地轴的自转外，还主要包括地轴围绕黄北极的进动，地轴的章动和极移。地球围绕太阳的公转不产生地轴在惯性空间的变化，对星敏感器的测试不会产生影响。
系统坐标系的建立
下面将对本发明中所使用的天球赤道坐标系、历元天球黄道坐标系、地固坐标系和星敏感器坐标系这四个坐标系系统进行详细说明。
1)天球赤道坐标系：使用CRF(Celestial Reference Frame)表示，考虑到岁差和章动的影响，天球赤道坐标系是与时间相关的。为系统分析方便，国际上建立了J2000.0天球赤道坐标系，简称J2000.0坐标系，使用符号CRFJ2000表示，如图4中的CRFJ2000坐标系所示。J2000.0坐标系是在公元2000年1月1日地球力学时12时建立的天球赤道坐标系，Z轴指向地球的北极，X轴指向建立时刻的春分点，Y轴与X轴、Z轴满足右手定则。星敏感器有关导航星的信息都是基于此而建立。在星敏感器中的导航星位置都用此坐标系表示。由于岁差和章动等影响，不同时刻的天球坐标系会发生相应的旋转。某一时刻的天球坐标系需要在J2000.0的基础上消除岁差和章动的影响才可获得，使用符号CRFT表示。
2)历元天球黄道坐标系：用ERF(Ecliptic Reference Frame)来表示，如图4中的XERF、YERF和ZERF所标示。定义建立在公元2000年1月1日地球力学时12时，并保持固定不变。地球绕太阳的公转轨道被称之为黄道，以地心为中心，以指向建立时刻的春分点为X轴，以垂直于黄道平面为Z轴，Y轴与X轴、Z轴满足右手定则，J2000坐标系的X轴与黄道坐标系的X轴一致，历元天球黄道坐标系的Z轴与J2000坐标系的Z轴夹角为23°26′21″，天球赤道坐标系绕着历元天球黄道坐标系的Z轴以每年50.29″的速度旋转，称之为岁差。
3)地固坐标系：地固坐标系的坐标轴定义和天球坐标系一致，但区别是，随着地球运动，地固坐标系围绕着地球的Z轴(即天球坐标系的Z轴)作近似匀速转动，角速度为Ω＝7.292115×10-5rad/s。地固坐标系使用如图4中所示的TRF(Terrestrial Reference Frame)来表示。
4)星敏感器坐标系：星敏感器坐标系固连于星敏感器上，并与之一同运动。其中心为星敏感器的探测器中心。X轴和Y轴分别平行于探测器的行和列，Z轴与另外两轴满足右手定则，用SCF表示(Star tracker Coordinate Frame)，如图4中的XSCF、YSCF和ZSCF所示。在使用时，将星敏感器与地球固定在一起，随着地固坐标系一起运动。
星敏感器所测量的导航星都是恒星，距离非常遥远，因此上述的4个坐标系统的坐标原点都可以认为是在同一点，坐标系之间的变换就只有旋转变换了。旋转变换的基本方法如下：
设x，y，z为原坐标系下的坐标，
Figure BDA0000074120090000081
为坐标系发生旋转之后的坐标，则
其中坐标系分别绕X轴、Y轴、Z轴旋转的坐标变换基为：
Figure BDA0000074120090000091
Figure BDA0000074120090000092
Figure BDA0000074120090000093
本发明的发明人在长期的研究中发现，通过利用地球本身自转的精密性，将星敏感器固连于地球，使星敏感器的主轴正对天顶进行观测，星敏感器随着地球的一起运动(Ω＝7.292115×10-5rad/s)，星敏感器测量值的角度变化与之相对应，而存储在星敏感器星表内的导航星是在J2000.0坐标系(CRFJ2000)下的坐标，由于星敏感器的三轴精度不一致性，其指向精度较滚转精度高一个量级，为保证测量指向精度的准确性和高精度，将星敏感器中导航星的坐标转换到当前测量时刻地固坐标系(TRF)下的坐标，这样就消除了地球滚转轴对指向精度的影响，此时测量星敏感器的输出结果理论上为恒定值，即星敏感器坐标系相对于地固坐标系的安装矩阵。以此矩阵为基础，可以测量出星敏感器主轴在地固坐标系中的变化，进而测量出星敏感器的指向轴精度。
下面将参照附图来详细描述本发明的星敏感器、用于星敏感器的精度测量方法和系统。
根据本发明的星敏感器1，所述星敏感器1可接收时间。具体而言，该星敏感器1可以包括：存储器(未示出)。所述存储器中存储有由导航星所构成的导航星表，且该星敏感器1中存储有与导航星相关联的导航星视运动参数。
根据本发明的星敏感器1，由于该星敏感器1可以具有星表转换功能以及输入时间参数，以方便在使用星敏感器1的过程中、利用本发明的方法和系统来对所述星敏感器1的精度进行测量。为方便实施本发明，所述导航星表可以基于J2000.0坐标系所形成。该星敏感器用于将基于J2000.0坐标系的导航星表转换成基于地固坐标系的导航星表。
根据本发明的一个实施例，所述导航星表包括各导航星的视运动参数。在制造的过程中，出于后续方便的考虑，所述导航星表可以一次固化在所述存储器4中。
下面将参照图5来说明用于星敏感器的精度测量方法。如图5中所示，该精度测量方法可以包括如下步骤：
1)将星敏感器固定在地球上，且使得星敏感器的主轴指向天顶，所述星敏感器可输入时间参数(步骤S1)。在该步骤S1中，通过将星敏感器固定在地球上，为尽量减小大气等影响，将星敏感器正对天顶，这样星敏感器就可以随着地球的运动输出相应的姿态和图像信息。星敏感器的精度测试问题就转换为星敏感器的测量结果与地球的转动进行精确比对的问题。
2)向所述星敏感器输入测试开始时间相对于J2000.0时刻的时间T(步骤S2)；
3)根据星敏感器中的导航星在J2000.0坐标系下的赤纬和赤经(α，δ)以及在两个方向上的视运动参数来确定导航星在当前时刻在J2000.0直角坐标系下的方向矢量(步骤S3)；
4)将导航星在当前时刻在J2000.0直角坐标系下的方向矢量转换为历元黄道坐标系下的方向矢量(步骤S4)；
5)将历元黄道坐标系下的方向矢量转变成T时刻下的天球坐标系下的方向矢量(vCRFT)(步骤S5)；
6)根据实际拍摄时刻(T+Δt)将导航星在T时刻从天球坐标系下的方向矢量(vCRFT)变到实际拍摄时刻(T+Δt)在地固坐标系下的方向矢量(vTRF)，并基于所述地固坐标系下的方向矢量(vTRF)，获得所述星敏感器的精度(步骤S6)。
由此，在本发明的上述精度测量方法中，通过利用地球本身自转的精密性，将星敏感器固连于地球，使星敏感器的主轴正对天顶进行观测，星敏感器随着地球的一起运动(Ω＝7.292115×10-5rad/s)，星敏感器测量值的角度变化与之相对应，而存储在星敏感器星表内的导航星是在J2000.0坐标系(CRFJ2000)下的坐标，由于星敏感器的三轴精度不一致性，其指向精度较滚转精度高一个量级，为保证测量指向精度的准确性和高精度，将星敏感器中导航星的坐标转换到当前测量时刻地固坐标系(TRF)下的坐标，这样就消除了地球滚转轴对指向精度的影响，此时测量星敏感器的输出结果理论上为恒定值，即星敏感器坐标系相对于地固坐标系的安装矩阵，以此矩阵为基础可以测量出星敏感器主轴在地固坐标系中的变化，进而测量出星敏感器的指向轴精度。
下面将详细描述上述精度测量方法中的各步骤。
在步骤S3中，在所述时间T下，导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)：
在所述步骤S4中，历元黄道坐标系下的方向矢量(vERF)基于所述导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)和将所述J2000.0坐标系绕X轴逆时针方向转动23°26′21″的方向变换之后获得：
Figure BDA0000074120090000111
根据本发明的一个实施例，将导航星在历元黄道坐标系下的方向矢量(vERF)转变成T时刻下的天球坐标系下的方向矢量通过下述获得：
将历元黄道坐标下的方向矢量(vERF)绕Z轴顺时针方向转动此时，岁差的影响已经消除，接着绕X轴顺时针方向转动23°26′21″；将坐标系绕X轴逆时针方向旋转εA，将坐标系绕Z轴顺时针方向旋转
Figure BDA0000074120090000113
绕X轴顺时针方向旋转
Figure BDA0000074120090000114
此时可以获得含有章动项的时刻T的天球坐标系下的方向矢量(vCRFT)，其中
Figure BDA0000074120090000115
Figure BDA0000074120090000116
分别表示黄经章动和斜章动。
具体而言，在该步骤中，所述导航星在天球坐标系下的方向矢量(vCRFT)通过下述公式获得：
Figure BDA0000074120090000117
其中Rx、Rz为绕X轴和Z轴旋转的坐标变换基，如前所述。
根据本发明的一个实施例，根据IAU2000B章动模型，εA与黄经章动
Figure BDA0000074120090000119
和斜章动
Figure BDA00000741200900001110
分别为：
Figure BDA00000741200900001111
Figure BDA00000741200900001112
Figure BDA00000741200900001113
其中，
Figure BDA00000741200900001114
Figure BDA00000741200900001115
Figure BDA00000741200900001116
t为从J2000.0开始的儒略世纪数并基于时刻T获得，式中的求和符号表示77个正弦余弦项的和，每一项均为一个正弦项和一个余弦项相加。此外，在上式中，幅角αi为幅角的线性组合：
Figure BDA00000741200900001117
Figure BDA00000741200900001118
式中，nik为整数，Fk为与太阳月亮位置有关的Delaunay幅角，具体而言，在上式中：
Figure BDA0000074120090000123
Figure BDA0000074120090000124
Figure BDA0000074120090000125
进一步地，章动表达式中的nik及的前10项在下述的表1、2中列出。其余的参数取值可以在国际地球自转和参考系服务(International Earth Rotation and ReferenceSystems Service)的网站：http://www.iers.org中查到。
章动表达式中的系数可以从《天球参考系变换及其应用》中查到(出版社：科学出版社；作者：李广宇；ISBN：9787030285102；出版年月：2010.08)。最终得到的系数的前10项如下表1和表2所示。
表1：章动量级数前10项幅角的系数
Figure BDA0000074120090000127
表2：章动量级数前10项的系数
Figure BDA0000074120090000128
根据本发明的一个实施例，所述步骤S6可以进一步包括：
(61)根据实际拍摄时刻T+Δt将导航星矢量从T时刻天球坐标系转到T+Δt时刻地固坐标系下的方向矢量(vTRF)；
(62)根据所述地固坐标系下的方向矢量(vTRF)通过QUEST方法求解星敏感器的最优姿态矩阵以及
(63)计算实际拍摄时刻
Figure BDA0000074120090000132
的星敏感器主轴指向矢量
Figure BDA0000074120090000133
以及
(64)计算实际拍摄时刻
Figure BDA0000074120090000134
的星敏感器主轴指向矢量的夹角(αij)，以获得所述星敏感器的指向精度。
导航星在地固坐标系下的方向矢量(vTRF)通过将所述导航星在天球坐标系下的方向矢量(vCRFT)绕天球坐标系的Z轴以Ω＝7.292115×10-5rad/s逆时针旋转获得：
Figure BDA0000074120090000135
Figure BDA0000074120090000136
根据本发明的一个实施例，所述最优姿态矩阵
Figure BDA0000074120090000137
通过使得下面的目标函数达到最小值而获得：
Figure BDA0000074120090000139
其中，wi，vi分别表示导航星在星敏感器感器坐标系下的方向矢量和在地固坐标系下的方向矢量，αi表示加权系数，满足；
所述星敏感器主轴指向矢量
Figure BDA00000741200900001311
为：
Figure BDA00000741200900001312
根据本发明的一个实施例，所述星敏感器主轴指向矢量的夹角(αij)为：
Figure BDA00000741200900001313
其中，iγj，统计αij即可以表示星敏感器的精度的评价标准。
在上述精度测量方法中，其中步骤S1-S5只需进行一次，步骤S6需要时刻转换，方可得到随着实际拍摄时刻
Figure BDA00000741200900001314
而变化的任意时刻的导航星相对于地固坐标系下的坐标数据，通过求解星敏感器的最优姿态矩阵计算不同时刻的星敏感器主轴指向
Figure BDA00000741200900001316
计算不同时刻星敏感器主轴指向矢量的夹角αij，统计αij即可以表示星敏感器指向轴的精度，如图8中所示。其中在图8中，星敏感器的指向轴11发生在星敏感器1随着地球4的自转而测量星空的过程中会发生角度的变化，并且这个角度变化之间的夹角(即星敏感器1的主轴指向矢量之间的夹角)可以用作表示该星敏感器1的指向精度。
下面将参照图6来详细描述根据本发明的一个实施例的用于测量星敏感器的精度测量系统。如图6中所示，该精度测量系统100可以包括：星敏感器1、固定器102和星敏感器精度测量单元103。星敏感器1可以包括导航星表和用于接收输入测试开始时间的时间输入接口101，且所述星敏感器1的主轴指向天顶，所述导航星表包括导航星视运动参数。固定器102用于固定所述星敏感器，其可以例如为三脚架。如前所述，通过将星敏感器1固定在地球上，为尽量减小大气等影响，将星敏感器正对天顶，这样星敏感器就可以随着地球的运动输出相应的姿态和图像信息。星敏感器的精度测试问题就转换为星敏感器的测量结果与地球的转动进行精确比对的问题。
在本发明的精度测量系统中，星敏感器精度测量单元103用于测量所述导航星的精度，其中通过所述时间输入接口向所述星敏感器输入测试开始时间相对于J2000.0时刻的时间T，根据星敏感器中的导航星在J2000.0坐标系下的赤纬和赤经(α，δ)以及在两个方向上的视运动参数
Figure BDA0000074120090000141
来确定导航星在当前时刻在J2000.0直角坐标系下的方向矢量，将导航星在当前时刻在J2000.0直角坐标系下的方向矢量转换为历元黄道坐标系下的方向矢量，将历元黄道坐标系下的方向矢量转变成T时刻下的天球坐标系下的方向矢量(vCRFT)，根据实际拍摄时刻(T+Δt)将导航星在T时刻从天球坐标系下的方向矢量(vCRFT)变到实际拍摄时刻(T+Δt)在地固坐标系下的方向矢量(vTRF)，并基于所述地固坐标系下的方向矢量(vTRF)获得所述星敏感器的精度。
根据本发明的上述精度测量系统，通过利用地球本身自转的精密性，将星敏感器1固连于地球，使星敏感器的主轴正对天顶进行观测，星敏感器随着地球的一起运动(Ω＝7.292115×10-5rad/s)，星敏感器测量值的角度变化与之相对应，而存储在星敏感器星表内的导航星是在J2000.0坐标系(CRFJ2000)下的坐标，由于星敏感器的三轴精度不一致性，其指向精度较滚转精度高一个量级，为保证测量指向精度的准确性和高精度，将星敏感器中导航星的坐标转换到当前测量时刻地固坐标系(TRF)下的坐标，这样就消除了地球滚转轴对指向精度的影响，此时测量星敏感器的输出结果理论上为恒定值，即星敏感器坐标系相对于地固坐标系的安装矩阵，以此矩阵为基础可以测量出星敏感器主轴在地固坐标系中的变化，进而测量出星敏感器的指向轴精度。
如图6中所示，该精度测量系统可以进一步包括：遮光罩104，所述遮光罩104套设在星敏感器1上，用于去除环境杂光的干扰。
根据本发明的一个实施例，如图7中所示，所述星敏感器精度测量单元103进一步包括：直角坐标方向矢量获取模块105，所述直角坐标方向矢量获取模块1031在所述时间T下通过下述公式获得所述导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)：
如图7中所示，所述星敏感器精度测量单元103进一步包括：历元黄道坐标系方向矢量(vERF)获取模块1032，所述历元黄道坐标系方向矢量获取模块1032基于所述导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)和将所述J2000.0坐标系绕X轴逆时针方向转动23°26′21″的方向变换之后获得：
Figure BDA0000074120090000152
进一步地，所述星敏感器精度测量单元103可以进一步包括：天球坐标系方向矢量获取模块1033，所述天球坐标系方向矢量获取模块1033通过下述将导航星在历元黄道坐标系下的方向矢量(vERF)转变成T时刻下的天球坐标系下的方向矢量：
将历元黄道坐标下的方向矢量(vERF)绕其Z轴顺时针方向转动
Figure BDA0000074120090000153
接着绕第一次转动后的坐标系的X轴顺时针方向转动
Figure BDA0000074120090000154
接着绕第二次旋转后的坐标系的X轴逆时针方向旋转εA；
接着绕第三次旋转后的坐标系的Z轴顺时针方向旋转以及
接着绕第四次旋转后的坐标系的X轴顺时针方向旋转
Figure BDA0000074120090000156
以获得含有章动项的当前时刻(T)的天球坐标系下的方向矢量(vCRFT)，其中
Figure BDA0000074120090000157
分别表示黄经章动和斜章动。
具体而言，所述天球坐标系方向矢量获取模块1033通过下述公式获得所述导航星在天球坐标系下的方向矢量(vCRFT)：
Figure BDA0000074120090000159
Figure BDA00000741200900001510
其中Rx、Rz为绕X轴和Z轴旋转的坐标变换基，如前所述。
根据本发明的一个实施例，根据IAU2000B章动模型，εA与黄经章动
Figure BDA00000741200900001511
和斜章动
Figure BDA00000741200900001512
分别为：
Figure BDA00000741200900001513
Figure BDA0000074120090000161
Figure BDA0000074120090000162
其中，
Figure BDA0000074120090000164
Figure BDA0000074120090000165
t为从J2000.0开始的儒略世纪数并基于时刻T获得，式中的求和符号表示77个正弦余弦项的和，每一项均为一个正弦项和一个余弦项相加。此外，在上式中，幅角αi为幅角的线性组合：
Figure BDA0000074120090000166
Figure BDA0000074120090000167
式中，nik为整数，Fk为与太阳月亮位置有关的Delaunay幅角。上述参数的各取值可以参见前述的精度测量方法中的详细说明，此处为简洁起见，不再赘述。
根据本发明的一个实施例，所述星敏感器精度测量单元103根据实际拍摄时刻
Figure BDA0000074120090000168
将导航星矢量从T时刻天球坐标系转到T+Δt时刻地固坐标系下的方向矢量(vTRF)；根据所述地固坐标系下的方向矢量(vTRF)通过QUEST方法求解星敏感器的最优姿态矩阵
Figure BDA0000074120090000169
计算实际拍摄时刻的星敏感器主轴指向矢量
Figure BDA00000741200900001611
以及计算实际拍摄时刻
Figure BDA00000741200900001612
的星敏感器主轴指向矢量的夹角(αij)，以获得所述星敏感器的指向精度。
根据本发明的一个实施例，所述星敏感器精度测量单元进一步包括：地固坐标系方向矢量获取模块1034，所述地固坐标系方向矢量获取模块1034通过将所述导航星在天球坐标系下的方向矢量(vCRFT)绕天球坐标系的Z轴以Ω＝7.292115×10-5rad/s逆时针旋转获得导航星在地固坐标系下的方向矢量(vTRF)：
Figure BDA00000741200900001613
Figure BDA00000741200900001614
根据本发明的一个实施例，所述最优姿态矩阵
Figure BDA00000741200900001615
通过使得下面的目标函数
Figure BDA00000741200900001616
达到最小值而获得：
Figure BDA00000741200900001617
其中，wi，vi分别表示导航星在星敏感器感器坐标系下的方向矢量和在地固坐标系下的方向矢量，αi表示加权系数，满足；
Figure BDA0000074120090000171
根据本发明的一个实施例，所述星敏感器主轴指向矢量为：
Figure BDA0000074120090000173
根据本发明的一个实施例，所述星敏感器主轴指向矢量的夹角(αij)为：
Figure BDA0000074120090000174
其中，iγj，统计αij即可以表示星敏感器的精度的评价标准。
通过求解星敏感器的最优姿态矩阵
Figure BDA0000074120090000175
计算不同时刻的星敏感器主轴指向
Figure BDA0000074120090000176
计算不同时刻星敏感器主轴指向矢量的夹角αij，统计αij即可以表示星敏感器指向轴的精度。
在本发明的该精度测量系统100中，还包括星敏感器精度输出单元105，该星敏感器精度输出单元105可以用于输出星敏感器精度测量单元103所测量的星敏感器主轴指向精度。如图6中所示，该系统100在操作中通过对实际星空的连续测量，利用星敏感器精度测量单元103即可以获得该星敏感器1的主轴指向精度。
在本发明的精度测量方法和系统中，通过利用地球本身自转的精密性，将星敏感器固连于地球，使星敏感器的主轴正对天顶进行观测。通过利用坐标变化并利用实时检测的结果，解决了传统的测试方法和系统中操作复杂、需要价格昂贵的精密转台和星模拟器的困扰，同时测量结果较转台式测量方法和系统更具有准确性，且更具有真实性，测试精度满足要求、过程简便、易于实现。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (11)

1.一种用于星敏感器的精度测量方法，包括如下步骤：
1)将星敏感器固定在地球上，且使得星敏感器的主轴指向天顶，所述星敏感器可输入时间参数且存储有导航星表和导航星的视运动参数；
2)向所述星敏感器输入测试开始时间相对于J2000.0时刻的当前时刻(T)；
3)根据星敏感器中的导航星在J2000.0坐标系下的赤纬和赤经(α，δ)以及在两个方向上的视运动参数
Figure FDA0000074120080000011
来确定导航星在当前时刻(T)在J2000.0直角坐标系下的方向矢量；
4)将导航星在当前时刻(T)在J2000.0直角坐标系下的方向矢量转换为历元黄道坐标系下的方向矢量；
5)将历元黄道坐标系下的方向矢量转变成当前时刻(T)下的天球坐标系下的方向矢量(vCRFT)；以及
6)根据实际拍摄时刻(T+Δt)将导航星在当前时刻(T)从天球坐标系下的方向矢量(vCRFT)变到实际拍摄时刻(T+Δt)在地固坐标系下的方向矢量(vTRF)，并基于所述地固坐标系下的方向矢量(vTRF)，获得所述星敏感器的精度。
2.根据权利要求1所述的精度测量方法，其特征在于，在所述步骤3)中，在所述当前时刻(T)下，导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)为：
3.根据权利要求2所述的精度测量方法，其特征在于，在所述步骤(4)中，历元黄道坐标系下的方向矢量(vERF)基于所述导航星在J2000.0直角坐标系下的方向矢量(vCRFJ2000)且将所述J2000.0坐标系绕绕着所述J2000.0坐标系的X轴逆时针方向转动23°26′21″的变换之后获得：
4.根据权利要求3所述的精度测量方法，其特征在于，将导航星在历元黄道坐标系下的方向矢量(vERF)转变成当前时刻(T)下的天球坐标系下的方向矢量(vCRFT)通过下述获得：
将历元黄道坐标下的方向矢量(vERF)绕其Z轴顺时针方向转动
Figure FDA0000074120080000021
接着绕第一次转动后的坐标系的X轴顺时针方向转动
Figure FDA0000074120080000022
接着绕第二次旋转后的坐标系的X轴逆时针方向旋转εA；
接着绕第三次旋转后的坐标系的Z轴顺时针方向旋转
Figure FDA0000074120080000023
以及
接着绕第四次旋转后的坐标系的X轴顺时针方向旋转
Figure FDA0000074120080000024
以获得含有章动项的当前时刻(T)的天球坐标系下的方向矢量(vCRFT)，其中
Figure FDA0000074120080000025
分别表示黄经章动和斜章动。
5.根据权利要求4所述的精度测量方法，其特征在于，所述导航星在天球坐标系下的方向矢量(vCRFT)通过下述公式获得：
Figure FDA0000074120080000027
其中Rx、Rz为坐标变换基。
6.根据权利要求5所述的精度测量方法，其特征在于，根据IAU2000B章动模型，εA与黄经章动
Figure FDA0000074120080000029
和斜章动
Figure FDA00000741200800000210
分别为：
Figure FDA00000741200800000211
Figure FDA00000741200800000213
其中，
Figure FDA00000741200800000214
Figure FDA00000741200800000216
t为从J2000.0开始的儒略世纪数并基于当前时刻(T)获得；
幅角αi为幅角的线性组合：
Figure FDA00000741200800000218
式中，nik为整数，Fk为与太阳月亮位置有关的Delaunay幅角。
7.根据权利要求1所述的精度测量方法，其特征在于，所述步骤(6)进一步包括：
(61)根据实际拍摄时刻(T+Δt)将导航星矢量从当前时刻(T)天球坐标系转到实际拍摄时刻(T+Δt)的地固坐标系下的方向矢量(vTRF)；
(62)根据所述地固坐标系下的方向矢量(vTRF)通过QUEST方法求解星敏感器的最优姿态矩阵
Figure FDA0000074120080000031
以及
(63)计算实际拍摄时刻
Figure FDA0000074120080000032
的星敏感器主轴指向矢量
Figure FDA0000074120080000033
以及
(64)计算实际拍摄时刻的星敏感器主轴指向矢量的夹角(αij)，以获得所述星敏感器的指向精度。
8.根据权利要求7所述的精度测量方法，其特征在于，导航星在地固坐标系下的方向矢量(vTRF)通过将所述导航星在天球坐标系下的方向矢量(vCRFT)绕天球坐标系的Z轴以Ω＝7.292115×10-5rad/s逆时针旋转获得：
Figure FDA0000074120080000035
Figure FDA0000074120080000036
9.根据权利要求8所述的精度测量方法，其特征在于，所述最优姿态矩阵
Figure FDA0000074120080000037
通过使得下面的目标函数
Figure FDA0000074120080000038
达到最小值而获得：
Figure FDA0000074120080000039
其中，wi，vi分别表示导航星在星敏感器感器坐标系下的方向矢量和在地固坐标系下的方向矢量，αi表示加权系数，满足；
Figure FDA00000741200800000310
10.根据权利要求9所述的精度测量方法，其特征在于，所述星敏感器的主轴指向矢量
Figure FDA00000741200800000311
为：
Figure FDA00000741200800000312
11.根据权利要求10所述的精度测量方法，其特征在于，所述星敏感器的主轴指向矢量的夹角(αij)为：
Figure FDA00000741200800000313
其中，iγj。
CN2011101892642A 2011-07-06 2011-07-06 用于星敏感器的精度测量方法 Active CN102288201B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN2011101892642A CN102288201B (zh) 2011-07-06 2011-07-06 用于星敏感器的精度测量方法
US13/386,979 US8433515B2 (en) 2011-07-06 2011-08-03 Method for measuring precision of star sensor and system using the same
PCT/CN2011/077977 WO2013004033A1 (zh) 2011-07-06 2011-08-03 用于星敏感器的精度测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101892642A CN102288201B (zh) 2011-07-06 2011-07-06 用于星敏感器的精度测量方法

Publications (2)

Publication Number Publication Date
CN102288201A true CN102288201A (zh) 2011-12-21
CN102288201B CN102288201B (zh) 2012-08-22

Family

ID=45334822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101892642A Active CN102288201B (zh) 2011-07-06 2011-07-06 用于星敏感器的精度测量方法

Country Status (1)

Country Link
CN (1) CN102288201B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607597A (zh) * 2012-03-08 2012-07-25 清华大学 星敏感器的三轴精度表述与测量方法
CN102706363A (zh) * 2012-05-25 2012-10-03 清华大学 一种高精度星敏感器的精度测量方法
CN104280049A (zh) * 2014-10-20 2015-01-14 北京控制工程研究所 一种高精度星敏感器外场精度测试方法
CN105371870A (zh) * 2015-11-27 2016-03-02 北京空间飞行器总体设计部 一种基于星图数据的星敏感器在轨精度测量方法
CN107085634A (zh) * 2017-04-12 2017-08-22 上海航天控制技术研究所 快速计算太阳光与太阳同步卫星星敏感器最小夹角的方法
CN107861919A (zh) * 2017-09-22 2018-03-30 上海航天控制技术研究所 一种星敏感器定姿星矢量权值实时分配方法
CN111879299A (zh) * 2020-08-17 2020-11-03 中国科学院上海天文台 一种地基望远镜全自动指星方法
CN112019745A (zh) * 2020-08-31 2020-12-01 苏州振旺光电有限公司 一种获取天区目标图像的方法和天文摄影设备
CN115236961A (zh) * 2022-04-02 2022-10-25 北京航天飞行控制中心 火星时间的确定方法、装置、存储介质以及电子装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02147813A (ja) * 1988-11-29 1990-06-06 Toshiba Corp ２軸スターセンサ
JPH03164398A (ja) * 1989-11-21 1991-07-16 Toshiba Corp ２軸スターセンサ
US6108593A (en) * 1997-07-09 2000-08-22 Hughes Electronics Corporation Method and apparatus for estimating attitude sensor bias in a satellite
KR20010011662A (ko) * 1999-07-29 2001-02-15 이우백 셔터를 이용한 3축자세 제어용 별 감지 장치
US6691033B1 (en) * 2000-07-26 2004-02-10 Hughes Electronics Corporation System and method for calibrating inter-star-tracker misalignments in a stellar inertial attitude determination system
CN100348460C (zh) * 2005-10-12 2007-11-14 北京航空航天大学 一种基于星场的星敏感器校准方法
CN101846510A (zh) * 2010-05-28 2010-09-29 北京航空航天大学 一种基于星敏感器和陀螺的高精度卫星姿态确定方法
CN101893440A (zh) * 2010-05-19 2010-11-24 哈尔滨工业大学 基于星敏感器的天文自主导航方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02147813A (ja) * 1988-11-29 1990-06-06 Toshiba Corp ２軸スターセンサ
JPH03164398A (ja) * 1989-11-21 1991-07-16 Toshiba Corp ２軸スターセンサ
US6108593A (en) * 1997-07-09 2000-08-22 Hughes Electronics Corporation Method and apparatus for estimating attitude sensor bias in a satellite
KR20010011662A (ko) * 1999-07-29 2001-02-15 이우백 셔터를 이용한 3축자세 제어용 별 감지 장치
US6691033B1 (en) * 2000-07-26 2004-02-10 Hughes Electronics Corporation System and method for calibrating inter-star-tracker misalignments in a stellar inertial attitude determination system
CN100348460C (zh) * 2005-10-12 2007-11-14 北京航空航天大学 一种基于星场的星敏感器校准方法
CN101893440A (zh) * 2010-05-19 2010-11-24 哈尔滨工业大学 基于星敏感器的天文自主导航方法
CN101846510A (zh) * 2010-05-28 2010-09-29 北京航空航天大学 一种基于星敏感器和陀螺的高精度卫星姿态确定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《光学精密工程》 20090731 李葆华等 天文导航中的星敏感器技术 全文 1-11 第17卷, 第7期 *
《宇航学报》 20020531 杨博 航天器星敏感器自主定位方法及精度分析 全文 1-11 第23卷, 第3期 *
《清华大学学报（自然科学版）》 20051231 邢飞等 星敏感器参数分析与自主校正 全文 1-11 第45卷, 第11期 *
《航天控制》 20011231 杨博等 一种用星敏感器自主定位方法的精度分析 全文 1-11 , 第1期 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607597B (zh) * 2012-03-08 2014-05-21 清华大学 星敏感器的三轴精度表述与测量方法
CN102607597A (zh) * 2012-03-08 2012-07-25 清华大学 星敏感器的三轴精度表述与测量方法
CN102706363A (zh) * 2012-05-25 2012-10-03 清华大学 一种高精度星敏感器的精度测量方法
CN102706363B (zh) * 2012-05-25 2015-04-15 清华大学 一种高精度星敏感器的精度测量方法
CN104280049A (zh) * 2014-10-20 2015-01-14 北京控制工程研究所 一种高精度星敏感器外场精度测试方法
CN104280049B (zh) * 2014-10-20 2017-04-05 北京控制工程研究所 一种高精度星敏感器外场精度测试方法
CN105371870B (zh) * 2015-11-27 2018-11-23 北京空间飞行器总体设计部 一种基于星图数据的星敏感器在轨精度测量方法
CN105371870A (zh) * 2015-11-27 2016-03-02 北京空间飞行器总体设计部 一种基于星图数据的星敏感器在轨精度测量方法
CN107085634A (zh) * 2017-04-12 2017-08-22 上海航天控制技术研究所 快速计算太阳光与太阳同步卫星星敏感器最小夹角的方法
CN107861919A (zh) * 2017-09-22 2018-03-30 上海航天控制技术研究所 一种星敏感器定姿星矢量权值实时分配方法
CN107861919B (zh) * 2017-09-22 2021-03-26 上海航天控制技术研究所 一种星敏感器定姿星矢量权值实时分配方法
CN111879299A (zh) * 2020-08-17 2020-11-03 中国科学院上海天文台 一种地基望远镜全自动指星方法
CN112019745A (zh) * 2020-08-31 2020-12-01 苏州振旺光电有限公司 一种获取天区目标图像的方法和天文摄影设备
CN112019745B (zh) * 2020-08-31 2021-08-24 苏州振旺光电有限公司 一种获取天区目标图像的方法和天文摄影设备
CN115236961A (zh) * 2022-04-02 2022-10-25 北京航天飞行控制中心 火星时间的确定方法、装置、存储介质以及电子装置
CN115236961B (zh) * 2022-04-02 2023-09-12 北京航天飞行控制中心 火星时间的确定方法、装置、存储介质以及电子装置

Also Published As

Publication number Publication date
CN102288201B (zh) 2012-08-22

Similar Documents

Publication Publication Date Title
CN102288201B (zh) 用于星敏感器的精度测量方法
CN102288200B (zh) 用于星敏感器的精度测量系统
CN102252673B (zh) 一种星敏感器在轨光行差的修正方法
US8433515B2 (en) Method for measuring precision of star sensor and system using the same
CN101344391B (zh) 基于全功能太阳罗盘的月球车位姿自主确定方法
CN102706363B (zh) 一种高精度星敏感器的精度测量方法
Zhang Star identification
CN103323026B (zh) 星敏感器和有效载荷的姿态基准偏差估计与修正方法
CN109738919A (zh) 一种用于gps接收机自主预测星历的方法
CN103727937A (zh) 一种基于星敏感器的舰船姿态确定方法
Haiyong et al. Image smearing modeling and verification for strapdown star sensor
CN103591966A (zh) 一种星模拟器测试平台及测试标定方法
CN105446346A (zh) 遥感卫星对月相对定标姿态调整方法
CN102607597B (zh) 星敏感器的三轴精度表述与测量方法
Ye et al. Looking vector direction analysis for the Moon-based Earth observation optical sensor
CN201629094U (zh) 辅助观星系统
CN107421504B (zh) 月基对地观测光电影像的拍摄时间计算方法
CN108896038A (zh) 月基光学传感器成像方法
CN104458653A (zh) 一种测量大天顶距处的大气折射值的方法和系统
Yang et al. Laboratory Test System Design for Star Sensor Performance Evaluation.
Kaplan New technology for celestial navigation
Zhang et al. A new method of single celestial-body sun positioning based on theory of mechanisms
Zhang et al. Random Star Recognition Algorithm Based on Image Total Station and Its Application to Astronomical Positioning
Yu et al. Differential measurement of atmospheric refraction using a telescope with double fields of view
Zhan et al. Absolute positioning based on the sun for mars rover

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20111221

Assignee: Beijing Tianyin star technology limited liability company

Assignor: Tsinghua University

Contract record no.: 2017990000122

Denomination of invention: Precision measurement method for star sensor

Granted publication date: 20120822

License type: Exclusive License

Record date: 20170407

EE01 Entry into force of recordation of patent licensing contract
TR01 Transfer of patent right

Effective date of registration: 20171206

Address after: 100192 3 floor, No. 2, No. 27, Longgang Road, Haidian District, Beijing

Patentee after: Beijing Tianyin star technology limited liability company

Address before: 100084 Haidian District 100084-82 mailbox Beijing

Patentee before: Tsinghua University