CN108759798B - 一种高精度航天器精度测量的实现方法 - Google Patents
一种高精度航天器精度测量的实现方法 Download PDFInfo
- Publication number
- CN108759798B CN108759798B CN201810638628.2A CN201810638628A CN108759798B CN 108759798 B CN108759798 B CN 108759798B CN 201810638628 A CN201810638628 A CN 201810638628A CN 108759798 B CN108759798 B CN 108759798B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- payload
- vector
- spacecraft
- prism
- 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.)
- Active
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000013598 vector Substances 0.000 claims abstract description 62
- 230000008859 change Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 31
- 230000009466 transformation Effects 0.000 claims description 27
- 108700041286 delta Proteins 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本发明提供一种高精度航天器精度测量方法,利用光电自准直仪对有效载荷基准棱镜的两个正交法向矢量的角度进行测量,得到棱镜单位矢量在航天器参考坐标系的数学描述,再通过旋转矢量算法获得有效载荷基准棱镜的矢量变化关系,进而解算有效载荷在参考坐标系的姿态角度以及有效载荷之间的相对姿态关系。本发明主要解决了航天器总装、试验等过程中有效载荷角秒级精度测量需求,并且具有测量精度高、适用性强等优点。
Description
技术领域
本发明涉及航天器总装、试验等过程中有效载荷角秒级精度测量,具体涉及一种高精度航天器精度测量方法。
背景技术
为保证有效载荷的正常工作和卫星性能指标的实现,必须对有效载荷、姿控系统测量执行部件以及其它系统相关单机的精度在地面总装、试验等过程中进行测量、调整和控制。一般情况下,在结构上设置正六面体棱镜或多面体棱镜,利用棱镜光轴表征被测载荷或单机的物理特性,安装、装配精度测量主要是通过电子经纬仪建站测量实施,需要两台或多台经纬仪互相瞄准进而计算出棱镜光轴之间的角度关系。
近年来高分辨率对地观测需求呈现不断增长的态势,对光学遥感卫星的空间分辨率、时间分辨率、光谱分辨率及定位精度等提出了更高的要求。航天器有效载荷姿态角度的测量精度从几角分到几十角秒,甚至达到几个角秒。传统电子经纬仪精测方法由于受人为影响因素较大,操作比较复杂、测量环节较多、误差因素复杂,目前测量最高精度可以达到±5″,无法满足角秒级有效载荷姿态角度精测要求。
发明内容
针对航天器有效载荷角秒级精度测量需求,本发明提供了一种高精度航天器精度测量方法。
本发明具有通过以下技术方案实现:
一种高精度航天器精度测量方法,利用光电自准直仪对有效载荷基准棱镜的两个正交法向矢量的角度进行测量,得到棱镜单位矢量在航天器参考坐标系的数学描述,再通过旋转矢量算法获得有效载荷基准棱镜的矢量变化关系,进而解算有效载荷在参考坐标系的姿态角度以及有效载荷之间的相对姿态关系;具体包括如下步骤:
S1、设置有效载荷基准棱镜坐标系与有效载荷局部坐标系一致,采用光电自准直仪建立测角系统,定义测量坐标系{A}和{B},对有效载荷基准棱镜的两个正交法向矢量的角度进行测量;读取自准直仪读数τ1和δ1、τ2和δ2,并统一坐标系得到测量坐标{A}下的基准棱镜矢量描述VA;
S2、设置自准直仪测量坐标系与有效载荷初始局部坐标系重合,对于星体上有效载荷1、2,{XoYoZo}为航天器参考坐标系,{XsiYsiZsi}为第i个有效载荷局部坐标系(i=1,2),[Vix;Viy;Viz]为测得的第i个有效载荷基准棱镜矢量(i=1,2);星体放置在转台上,可以绕轴线ZO旋转。在自准直仪测量坐标系固定不变的基础上,对星体进行一次小角度旋转。
S3、获得有效载荷1局部坐标系下,有效载荷1基准棱镜旋转前初始矢量为V1 S1,第一次旋转后棱镜矢量为V1 S1',以及局部坐标下的转化矩阵T′S1;再统一在航天器参考坐标系下,有效载荷1基准棱镜旋转前初始矢量为V1 O,第一次旋转后棱镜矢量为V1 O',以及航天器参考坐标下的转化矩阵T′O1,其中,有效载荷1坐标系{S1}到卫星参考坐标系{O}的转化矩阵即为需要求解的有效载荷1在参考坐标系的姿态角度。
S5、对于刚性结构而言,在航天器参考坐标系下进行小角度旋转,局部坐标系与公共坐标系的相对关系固定不变,由此建立关系式T′O1=T′O2。
S6、在自准直仪测量坐标系固定不变的基础上,再对星体进行一次小角度旋转,建立第二次旋转的有效载荷矢量变化关系式T″O1=T″O2。
本发明具有以下有益效果:
(1)本发明采用光电自准直仪对有效载荷基准棱镜两个法向矢量的测量可完整描述出航天器参考坐标系内的矢量特征,并通过旋转矢量算法建立有效载荷基准棱镜的矢量变化关系,所有计算均为代数计算,最复杂的求解步骤只是寻找一个三阶实对称矩阵的特征向量,该方法充分利用自准直仪高精度角度测量和代数求解的优势,避免了传统经纬仪建站互相瞄准引入的较大误差,具有测量精度高、计算快捷、操作简单等优点。
(2)光电自准直仪测角精度可达到±0.25″,旋转矢量算法求解精度小于0.2″,考虑环境条件、试验工装、二维转台旋转精度等其他因素的影响,通过合理配置,基于旋转矢量算法仪的航天器精度测量系统误差优于1",该方法有效可行,能够满足航天器地面总装、试验等过程中对有效载荷、姿控系统测量执行部件以及其它系统相关单机的角秒级精度测量要求。
附图说明
图1为有效载荷局部坐标系在航天器参考坐标系下描述示意图;
图2为测量坐标系描述示意图;
图3为旋转矢量算法描述示意图。
图4为本发明实施例的流程图。
具体实施方式
下面对本发明的实施作详细说明,实施具体步骤以本发明技术方案为前提下进行,这里给出详细的实施方式和具体的操作过程。
1.有效载荷局部坐标系定义
有效载荷局部坐标系根据右手定则确定,如图1所示,{XoYoZo}为航天器固定参考坐标系,{XsiYsiZsi}为第i个有效载荷坐标系(i=1,2),使用Xsi、Ysi、Zsi三个坐标轴在航天器参考系中三坐标轴的分量值表示有效载荷局部坐标系与航天器参考坐标系的关系。
例如有效载荷i的Xsi轴,在航天器参考系三个轴方向的值分别为:Xix、Xiy、Xiz,与Xb、Yb、Zb的夹角为α、β、γ,且满足Xix=cosα,Xiy=cosβ,Xiz=cosγ,Xix 2+Xiy 2+Xiz 2=1。则XsiYsiZsi在XbYbZb中的完整表示形式为:
且满足约束:Xix 2+Xiy 2+Xiz 2=1,Yix 2+Yiy 2+Yiz 2=1,Zix 2+Ziy 2+Ziz 2=1。
2.棱镜测量矢量描述
(1)测量坐标系定义
光电自准直仪用于测量有效载荷基准棱镜光轴的角度变化,设置基准棱镜坐标系与有效载荷局部坐标系一致。单台有效载荷基准棱镜指向变化的为三维矢量,需两套二维光电自准直仪进行测量,如图2所示,镜面法向为Vx和Vy,初始状态下测角系统形成两个坐标系{A}和{B},其定义方式如下:
{A}:Vx在水平面的投影方向为X轴,指向镜面;Z轴竖直向上,Y轴由右手法则确定;
{B}:Vy在水平面的投影方向为Y轴,指向镜面;Z轴竖直向上,X轴由右手法则确定;
根据两个坐标系的定义,A和B之间只存在一个绕Z轴的旋转,旋转角度为α(这是个小角度,与棱镜的两个面的夹角有关系,即棱镜的垂直度误差),则坐标系{B}到坐标系{A}的转换矩阵可以定义为:
(2)棱镜矢量描述
初始状态下,自准直仪读数为τ1和δ1、τ2和δ2,在上述坐标系{A}和{B}下,棱镜矢量描述为:
则棱镜三维矢量分别描述为:
(3)航天器参考坐标系下描述
先将测量坐标系{A}统一到有效载荷局部坐标系,如下:
三维棱镜矢量,最后统一在航天器参考坐标下,
3.旋转矢量算法
旋转矢量算法描述如图3所示,数学原理是有效载荷局部坐标系在航天器参考坐标系下的相对关系固定不变。在测量坐标系不变的前提下,旋转航天器参考坐标系,通过自准直仪测量获得有效载荷基准棱镜的矢量变化,进而解算有效载荷在参考坐标系的姿态角度及有效载荷之间的相对姿态角度。
具体算法描述如下:
星体上有效载荷1、2,分别采用两台自准直仪测量有效载荷基准棱镜矢量,{XoYoZo}为航天器参考坐标系,{XsiYsiZsi}为第i个有效载荷局部坐标系(i=1,2),定义自准直仪测量坐标系与有效载荷局部坐标系重合,[Vix;Viy;Viz]为测得的第i个有效载荷基准棱镜矢量(i=1,2)。
星体放置在转台上,可以绕轴线ZO旋转。在自准直仪测量坐标系固定不变的基础上,对星体进行多次的小角度旋转,假设有N次,要求
N≥2
有效载荷1局部坐标系下,有效载荷1基准棱镜旋转前初始矢量为V1 S1,第一次旋转后棱镜矢量为V1 S1',局部坐标下存在转化矩阵T′S1,且有
V1 S1'=V1 S1·T′S1
统一在航天器参考坐标系下,有效载荷1基准棱镜旋转前初始矢量为V1 O,第一次旋转后棱镜矢量为V1 O',航天器参考坐标下存在转化矩阵T′O1,同样有
V1 O'=V1 O·T′O1
联合式子(11),(12),(13),得到转化矩阵T′O1的计算结果
同样地,对于有效载荷2基准棱镜,在航天器参考坐标下存在转化矩阵T′O2
上述转化矩阵均为3×3矩阵,式子中包含9个独立方程。
有效载荷局部坐标系与航天器参考坐标系的转化矩阵即为需要求解的有效载荷在参考坐标系的姿态角度。根据定义,每个有效载荷转化矩阵中含有9个变量,3个约束关系,若要求解中总共18个变量,6个约束关系,至少需要12个独立方程。因此,至少需要两次小角度旋转,联立如下方程组:
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (4)
1.一种高精度航天器精度测量的实现方法,其特征在于:利用光电自准直仪对有效载荷基准棱镜的两个正交法向矢量的角度进行测量,得到棱镜单位矢量在航天器参考坐标系的数学描述,再通过旋转矢量算法获得有效载荷基准棱镜的矢量变化关系,进而解算有效载荷在参考坐标系的姿态角度以及有效载荷之间的相对姿态关系,具体包括如下步骤:
S1、设置有效载荷基准棱镜坐标系与有效载荷局部坐标系一致,采用光电自准直仪建立测角系统,定义测量坐标系{A}和{B},对有效载荷基准棱镜的两个正交法向矢量的角度进行测量;读取自准直仪读数τ1和δ1、τ2和δ2,并统一坐标系得到测量坐标系{A}下的基准棱镜矢量描述VA;
S2、设置自准直仪测量坐标系与有效载荷初始局部坐标系重合,对于星体上有效载荷1、2,{XoYoZo}为航天器参考坐标系,{XsiYsiZsi}为第i个有效载荷局部坐标系(i=1,2),[Vix;Viy;Viz]为测得的第i个有效载荷基准棱镜矢量(i=1,2);星体放置在转台上,可以绕轴线ZO旋转,在自准直仪测量坐标系固定不变的基础上,对星体进行一次小角度旋转;
S3、获得有效载荷1局部坐标系下,有效载荷1基准棱镜旋转前初始矢量为第一次旋转后棱镜矢量为以及局部坐标下的转化矩阵T′S1;再统一在航天器参考坐标系下,有效载荷1基准棱镜旋转前初始矢量为V1 O,第一次旋转后棱镜矢量为V1 O′,以及航天器参考坐标下的转化矩阵T′O1,其中,有效载荷1坐标系{S1}到卫星参考坐标系{O}的转化矩阵即为需要求解的有效载荷1在参考坐标系的姿态角度,
S5、对于刚性结构而言,在航天器参考坐标系下进行小角度旋转,局部坐标系与公共坐标系的相对关系固定不变,由此建立关系式T′O1=T′O2,
S6、在自准直仪测量坐标系固定不变的基础上,再对星体进行一次小角度旋转,建立第二次旋转的有效载荷矢量变化关系式T″O1=T″O2,
2.根据权利要求1所述的一种高精度航天器精度测量的实现方法,其特征在于:所述的光电自准直仪用于测量有效载荷基准棱镜法向矢量的俯仰角和偏转角;所述的有效载荷基准棱镜的坐标系与有效载荷局部坐标系一致。
3.根据权利要求1所述的一种高精度航天器精度测量的实现方法,其特征在于:所述的棱镜单位矢量在航天器参考坐标系的数学描述是先通过数学表达式描述出棱镜的三个法线方向矢量在有效载荷局部坐标系的表示,再利用局部坐标系与参考坐标系的关系,推导出基准棱镜单位矢量在航天器参考坐标系下的表示。
4.根据权利要求1所述的一种高精度航天器精度测量的实现方法,其特征在于:所述的旋转矢量算法是依据有效载荷局部坐标系在航天器参考坐标系下的相对关系固定不变,在测量坐标系不变的前提下,旋转航天器参考坐标系,获得有效载荷基准棱镜的矢量变化关系,解算有效载荷在参考坐标系的姿态角度以及有效载荷之间的相对姿态关系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810638628.2A CN108759798B (zh) | 2018-06-20 | 2018-06-20 | 一种高精度航天器精度测量的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810638628.2A CN108759798B (zh) | 2018-06-20 | 2018-06-20 | 一种高精度航天器精度测量的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108759798A CN108759798A (zh) | 2018-11-06 |
CN108759798B true CN108759798B (zh) | 2021-04-06 |
Family
ID=63979590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810638628.2A Active CN108759798B (zh) | 2018-06-20 | 2018-06-20 | 一种高精度航天器精度测量的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108759798B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109596110A (zh) * | 2018-11-22 | 2019-04-09 | 九江职业技术学院 | 一种机械设备位置检测装置 |
CN109373903A (zh) * | 2018-12-07 | 2019-02-22 | 银河航天(北京)通信技术有限公司 | 航天器部件间姿态关系确定系统和方法 |
CN109724577B (zh) * | 2019-01-30 | 2022-02-11 | 山东省国土测绘院 | 一种单波束与拖曳测量结合的水深数据处理方法 |
CN110514132B (zh) * | 2019-09-26 | 2021-01-05 | 中国科学院长春光学精密机械与物理研究所 | 大尺寸桁架式支撑结构稳定性测试装置及测试方法 |
CN110940296A (zh) * | 2019-10-30 | 2020-03-31 | 中国航天空气动力技术研究院 | 一种高超声速飞行器舵偏角测量方法 |
CN111121734A (zh) * | 2020-01-03 | 2020-05-08 | 中国船舶重工集团公司第七0七研究所 | 惯性设备安装基座变形测量装置及方法 |
CN111623775B (zh) * | 2020-05-15 | 2022-10-04 | 天津时空经纬测控技术有限公司 | 运载体姿态测量系统、方法以及存储介质 |
CN112894490B (zh) * | 2021-01-27 | 2022-07-26 | 天津大学 | 基于旋转l型阵列实现数控机床垂直度误差检测的方法 |
CN113094371B (zh) * | 2021-04-14 | 2023-05-12 | 嘉兴毕格智能科技有限公司 | 一种用户自定义坐标系的实现方法 |
CN113358024B (zh) * | 2021-05-08 | 2022-08-12 | 上海卫星工程研究所 | 一种卫星仪器相对精测数据处理方法及系统 |
CN115112150B (zh) * | 2022-06-29 | 2023-04-18 | 长光卫星技术股份有限公司 | 一种焦平面的指示棱镜标定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104990533A (zh) * | 2015-06-22 | 2015-10-21 | 哈尔滨工业大学 | 卫星地面物理仿真系统超高精度姿态测量方法及装置 |
CN106468544A (zh) * | 2016-08-31 | 2017-03-01 | 上海卫星工程研究所 | 基于光电自准直仪的卫星高精度测角方法 |
CN106767673A (zh) * | 2016-11-29 | 2017-05-31 | 上海卫星工程研究所 | 一种卫星高精度光学敏感载荷的指向测量方法 |
CN107543495A (zh) * | 2017-02-17 | 2018-01-05 | 北京卫星环境工程研究所 | 航天器设备自动准直测量系统、准直方法与测量方法 |
-
2018
- 2018-06-20 CN CN201810638628.2A patent/CN108759798B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104990533A (zh) * | 2015-06-22 | 2015-10-21 | 哈尔滨工业大学 | 卫星地面物理仿真系统超高精度姿态测量方法及装置 |
CN106468544A (zh) * | 2016-08-31 | 2017-03-01 | 上海卫星工程研究所 | 基于光电自准直仪的卫星高精度测角方法 |
CN106767673A (zh) * | 2016-11-29 | 2017-05-31 | 上海卫星工程研究所 | 一种卫星高精度光学敏感载荷的指向测量方法 |
CN107543495A (zh) * | 2017-02-17 | 2018-01-05 | 北京卫星环境工程研究所 | 航天器设备自动准直测量系统、准直方法与测量方法 |
Non-Patent Citations (2)
Title |
---|
航天器姿态系统的自适应鲁棒控制;袁国平;《中国博士学位论文全文数据库 工程科技Ⅱ辑 (月刊 )》;20130115(第1期);第10-178页 * |
角秒级小角度测量技术及其在航天器中的应用研究;王文亮;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑 (月刊 )》;20120315(第3期);第5-89页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108759798A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108759798B (zh) | 一种高精度航天器精度测量的实现方法 | |
CN106525073B (zh) | 一种基于三轴转台的惯性空间陀螺标定试验方法 | |
JP6613236B2 (ja) | 高精度な光ファイバジャイロ慣性測量装置の標定方法 | |
CN106289246B (zh) | 一种基于位置和姿态测量系统的柔性杆臂测量方法 | |
CN103994763B (zh) | 一种火星车的sins/cns深组合导航系统及其实现方法 | |
CN107728182B (zh) | 基于相机辅助的柔性多基线测量方法和装置 | |
CN106767673B (zh) | 一种卫星高精度光学敏感载荷的指向测量方法 | |
CN107042528A (zh) | 一种工业机器人的运动学标定系统及方法 | |
CN108106637B (zh) | 一种分布式pos的精度检校方法和装置 | |
CN103852085B (zh) | 一种基于最小二乘拟合的光纤捷联惯导系统现场标定方法 | |
CN102564455A (zh) | 星敏感器安装误差四位置标定与补偿方法 | |
CN102679999A (zh) | 星敏感器安装误差四位置标定与补偿方法 | |
CN112129322B (zh) | 一种捷联惯组与三轴转台的安装误差检测方法及校正方法 | |
CN105444778A (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
CN108375383A (zh) | 多相机辅助的机载分布式pos柔性基线测量方法和装置 | |
CN210293237U (zh) | 一种小卫星星上设备高精度姿态测量系统 | |
CN105424060B (zh) | 一种飞行器星敏感器与捷联惯组安装误差的测量方法 | |
CN113607188B (zh) | 基于经纬仪叉丝成像的多视场星敏感器的标定系统及方法 | |
CN114659709A (zh) | 一种用于大型带翼航天飞行器的三维质心测量方法 | |
Cui et al. | Novel method of rocket nozzle motion parameters non-contact consistency measurement based on stereo vision | |
CN112097798B (zh) | 高分六号卫星的高分相机的高精度校验方法和装置 | |
CN112484751B (zh) | 一种航天验证器相对大空间试验场坐标系位置姿态测量方法 | |
CN115752391A (zh) | 一种航天器总装过程中的设备安装精度调测方法 | |
CN104864876A (zh) | 一种月球车联合定位方法及系统 | |
CN111780787B (zh) | 一种基于光纤惯性辅助的mems惯组标定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |