CN113252074B - 一种基于双波束激光实时测量的星敏感器标定方法 - Google Patents
一种基于双波束激光实时测量的星敏感器标定方法 Download PDFInfo
- Publication number
- CN113252074B CN113252074B CN202110515754.0A CN202110515754A CN113252074B CN 113252074 B CN113252074 B CN 113252074B CN 202110515754 A CN202110515754 A CN 202110515754A CN 113252074 B CN113252074 B CN 113252074B
- Authority
- CN
- China
- Prior art keywords
- laser
- star sensor
- vector information
- coordinate system
- emission source
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Manufacturing & Machinery (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Navigation (AREA)
Abstract
本发明涉及航天器姿态控制技术领域,提出了一种基于双波束激光实时测量的星敏感器标定方法,包括:由激光发射和折返系统将激光发射源发射的两束激光折返进入星敏感器的视场;计算所述激光在激光发射源坐标系和折返镜坐标系中的矢量信息;计算所述激光在星敏感器坐标系中的矢量信息;计算星敏感器与激光发射源的相对姿态矩阵;以及计算星敏感器与光学载荷系统的相对姿态矩阵。本发明可以避免力学变形误差和热弹变形误差,并且可以使姿态确定系统能够在任何条件下精确测量,为航天器姿态精度的提高提供了一个新的思路。
Description
技术领域
本发明总的来说涉及航天器姿态控制技术领域,具体而言涉及一种基于双波束激光实时测量的星敏感器标定方法。
背景技术
航天器的姿态确定精度是航天器姿态控制的重要基础,也是航天器执行各类任务成功与否的关键。各类空间任务都对航天器的姿态确定精度提出了很高的要求,而通常情况下星敏感器是航天器各种姿态测量部件中精度最高的。
星敏感器在轨实时测量星敏光学测量坐标系相对于惯性空间的姿态信息,经由星敏感器和光学载荷系统间的相对姿态矩阵转换,得到满足任务需求的光学载荷系统相对于惯性空间的姿态指向信息o现有星敏感器的姿态确定精度可达亚角秒级,若相对姿态信息测量精度较差则会严重影响姿态确定系统的最终测量精度。
因此,需要测量星敏感器和光学载荷系统间的相对姿态矩阵信息,以完成星敏感器标定。现有技术中主要基于地面标定和在轨标定两种方法测量星敏感器和光学载荷系统间的相对姿态信息,下面对它们作介绍。
地面标定方法:采用地面激光测量系统测量星敏感器安装棱镜和光学载荷系统安装棱镜之间的相对姿态信息并一次性写入姿态确定软件,测量精度为角秒级。但是在航天器入轨发射过程中会发生结构力学变形导致相对安装信息发生较大变化,地面标定相对姿态矩阵信息精度因此变差。发射阶段的力学变形达百角秒量级,因此地面标定方法无法保证航天器在轨高精度的姿态确定。
在轨标定方法:由星敏感器和光学载荷系统分别同时获得相对惯性系的姿态信息,以两者各自的姿态信息之差表示相对姿态矩阵信息。该方法通常在航天器发射入轨后进行一次或周期性在轨标定完成相对姿态矩阵信息测量,可以避免发射过程中的力学变形带来的误差,但标定周期内航天器的热弹变形引入的误差无法消除,也限制了星敏感器高精度定姿能力。热弹变形导致的定姿误差达十角秒量级,严重影响了高精度姿态测量。
发明内容
针对现有技术中进行星敏感器姿态测量时,航天器发射阶段的力学变形和在轨阶段的热弹变形会引入误差项的问题,本发明提出一种基于双波束激光实时测量的星敏感器标定方法包括下列步骤:
由激光发射和折返系统将激光发射源发射的第一激光以及第二激光折返进入星敏感器的视场,其中所述激光发射和折返系统包括激光发射源和折返镜;
根据激光发射和折返系统的安装信息计算所述激光在激光发射源坐标系和折返镜坐标系中的矢量信息;
根据所述激光在星敏感器的敏感面上的激光成像位置信息计算所述激光在星敏感器坐标系中的矢量信息;
根据所述激光在激光发射源坐标系、折返镜坐标系和星敏感器坐标系中的矢量信息计算星敏感器与激光发射源的相对姿态矩阵;以及
根据星敏感器与激光发射源的相对姿态矩阵和光学载荷系统与激光发射源的相对姿态矩阵计算星敏感器与光学载荷系统的相对姿态矩阵。
在本发明一个实施例中规定,根据激光发射和折返系统的安装信息计算所述激光在激光发射源坐标系和折返镜坐标系中的矢量信息包括下列步骤:
以O-X0Y0Z0表示激光发射源坐标系,O-XiYiZi表示折返光路上的第i个折返镜的折返镜坐标系,将所述激光在O-X0Y0Z0中表示为下式:
根据激光发射和折返系统的安装信息得到O-XiYiZi相对于O-Xi-1Yi-1Zi-1的转置矩阵信息,表示为:Ai,i-1;
以M表示折返前后激光矢量的变换矩阵,表示为:M=diag(1,1,-1);
计算所述激光在第1个折返镜反射后在O-X1Y1Z1中的矢量信息,表示为下式:
计算所述激光在第i个折返镜反射后在O-XiYiZi中的矢量信息,表示为下式:
以O-XnYnZn表示所述激光在折返进入星敏感器的视场前的最后一个折返镜的折返镜坐标系,计算所述激光在O-XnYnZn中的矢量信息,表示为下式:
在本发明一个实施例中规定,根据所述激光在星敏感器的敏感面上的激光成像位置信息计算所述激光在星敏感器坐标系中的矢量信息包括下列步骤:
以O-XsYsZs表示星敏感器坐标系,以(x1,y1)表示第一激光在星敏感器的敏感面上的激光成像位置信息,以(x2,y2)表示第二激光在星敏感器的敏感面上的激光成像位置信息,将所述激光在O-XsYsZs中的矢量信息表示为下式:
在本发明一个实施例中规定,根据所述激光在激光发射源坐标系、折返镜坐标系和星敏感器坐标系中的矢量信息计算星敏感器与激光发射源的相对姿态矩阵包括下列步骤:
计算O-XsYsZs相对于O-XnYnZn的转置矩阵信息As,n,表示为下式:
AS,n=Ms·(Mn)T;以及
计算星敏感器与激光发射源的相对姿态矩阵,表示为下式:
AS,0=AS,n·An,n-1…A2,1·A1,0。
在本发明一个实施例中规定,根据星敏感器与激光发射源的相对姿态矩阵和光学载荷系统与激光发射源的相对姿态矩阵计算星敏感器与光学载荷系统的相对姿态矩阵包括下列步骤:
以AZ,0表示光学载荷系统与激光发射源的相对姿态矩阵,计算星敏感器与光学载荷系统的相对姿态矩阵,表示为下式:
AZ,S=AZ,0·(AS,0)T。
本发明至少具有如下有益效果:针对现有技术中进行星敏感器标定时,航天器发射阶段的力学变形和在轨阶段的热弹变形会引入误差项的问题,利用激光发射和转折系统精确测量了星敏感器与光学载荷系统之间的相对姿态矩阵,避免了力学变形误差和热弹变形误差,大大提高了所述相对姿态矩阵的测量精度;实时测量星敏感器和光学载荷系统之间的安装矩阵,使姿态确定系统能够在任何条件下精确测量光学载荷系统的姿态指向信息;为航天器姿态精度的提高提供了一个新的思路。本发明能够适用于大多数具有高精度姿态确定需求的航天器上,有较好的工程应用前景和推广价值。
附图说明
为进一步阐明本发明的各实施例的以上和其它优点和特征,将参考附图来呈现本发明的各实施例的更具体的描述。可以理解,这些附图只描绘本发明的典型实施例,因此将不被认为是对其范围的限制。在附图中,为了清楚明了,相同或相应的部件将用相同或类似的标记表示。
图1示出了本发明一个实施例中的整体流程图。
图2示出了本发明一个实施例中由激光发射和折返系统将激光发射源发射的两束激光折返进入星敏感器的视场的光路图。
具体实施方式
应当指出,各附图中的各组件可能为了图解说明而被夸大地示出,而不一定是比例正确的。在各附图中,给相同或功能相同的组件配备了相同的附图标记。
在本发明中,除非特别指出,“布置在…上”、“布置在…上方”以及“布置在…之上”并未排除二者之间存在中间物的情况。此外,“布置在…上或上方”仅仅表示两个部件之间的相对位置关系,而在一定情况下、如在颠倒产品方向后,也可以转换为“布置在…下或下方”,反之亦然。
在本发明中,各实施例仅仅旨在说明本发明的方案,而不应被理解为限制性的。
在本发明中,除非特别指出,量词“一个”、“一”并未排除多个元素的场景。
在此还应当指出,在本发明的实施例中,为清楚、简单起见,可能示出了仅仅一部分部件或组件,但是本领域的普通技术人员能够理解,在本发明的教导下,可根据具体场景需要添加所需的部件或组件。另外,除非另行说明,本发明的不同实施例中的特征可以相互组合。例如,可以用第二实施例中的某特征替换第一实施例中相对应或功能相同或相似的特征,所得到的实施例同样落入本申请的公开范围或记载范围。
在此还应当指出,在本发明的范围内,“相同”、“相等”、“等于”等措辞并不意味着二者数值绝对相等,而是允许一定的合理误差,也就是说,所述措辞也涵盖了“基本上相同”、“基本上相等”、“基本上等于”。以此类推,在本发明中,表方向的术语“垂直于”、“平行于”等等同样涵盖了“基本上垂直于”、“基本上平行于”的含义。
另外,本发明的各方法的步骤的编号并未限定所述方法步骤的执行顺序。除非特别指出,各方法步骤可以以不同顺序执行。
下面结合具体实施方式参考附图进一步阐述本发明。
图1示出了本实施例的整体流程。
本实施例采用高精度星敏感器以及激光发射和折返系统,包括如下步骤:
以O-X0Y0Z0表示激光发射源坐标系,O-XiYiZi表示折返光路上的第i个折返镜的折返镜坐标系,将所述激光在O-X0Y0Z0中表示为下式:
根据激光发射和折返系统的安装信息得到O-XiYiZi相对于O-Xi-1Yi-1Zi-1的转置矩阵信息,表示为:Ai,i-1;
以M表示折返前后矢量的变换矩阵,表示为:M=diag(1,1,-1);
计算所述激光在第1个折返镜反射后在O-X1Y1Z1中的矢量信息,表示为下式:
计算所述激光在第i个折返镜反射后在O-XiYiZi中的矢量信息,表示为下式:
以O-XnYnZn表示所述激光在折返进入星敏感器的视场前的最后一个折返镜的折返镜坐标系,计算所述激光在O-XnYnZn中的矢量信息,表示为下式:
根据所述激光在星敏感器的敏感面上的激光成像位置信息计算所述激光在星敏感器坐标系中的矢量信息包括下列步骤:
以O-XsYsZs表示星敏感器坐标系,以(x1,y1)表示第一激光在星敏感器的敏感面上的激光成像位置信息,以(x2,y2)表示第二激光在星敏感器的敏感面上的激光成像位置信息,将所述激光在O-XsYsZs中的矢量信息表示为下式:
图2示出了本实施例中由激光发射和折返系统将激光发射源发射的两束激光折返进入星敏感器的视场的光路图。
根据所述激光在激光发射源坐标系、折返镜坐标系和星敏感器坐标系中的矢量信息计算星敏感器与激光发射源的相对姿态矩阵包括下列步骤:
计算O-XsYsZs相对于O-XnYnZn的转置矩阵信息As,n,表示为下式:
AS,n=MS·(Mn)T;以及
计算星敏感器与激光发射源的相对姿态矩阵,表示为下式:
AS,0=AS,n·An,n-1…A2,1·A1,0。
根据星敏感器与激光发射源的相对姿态矩阵和光学载荷系统与激光发射源的相对姿态矩阵计算星敏感器与光学载荷系统的相对姿态矩阵,包括下列步骤:
以AZ,0表示光学载荷系统与激光发射源的相对姿态矩阵,计算星敏感器与光学载荷系统的相对姿态矩阵,表示为下式:
AZ,S=AZ,0·(AS,0)T。
尽管上文描述了本发明的各实施例,但是,应该理解,它们只是作为示例来呈现的,而不作为限制。对于相关领域的技术人员显而易见的是,可以对其做出各种组合、变型和改变而不背离本发明的精神和范围。因此,此处所公开的本发明的宽度和范围不应被上述所公开的示例性实施例所限制,而应当仅根据所附权利要求书及其等同替换来定义。
Claims (5)
1.一种基于双波束激光实时测量的星敏感器标定方法,其特征在于包括下列步骤:
由激光发射和折返系统将激光发射源发射的第一激光以及第二激光折返进入星敏感器的视场,其中所述激光发射和折返系统包括激光发射源和折返镜;
根据激光发射和折返系统的安装信息计算所述激光在激光发射源坐标系和折返镜坐标系中的矢量信息;
根据所述激光在星敏感器的敏感面上的激光成像位置信息计算所述激光在星敏感器坐标系中的矢量信息;
根据所述激光在激光发射源坐标系、折返镜坐标系和星敏感器坐标系中的矢量信息计算星敏感器与激光发射源的相对姿态矩阵;以及
根据星敏感器与激光发射源的相对姿态矩阵和光学载荷系统与激光发射源的相对姿态矩阵计算星敏感器与光学载荷系统的相对姿态矩阵。
2.根据权利要求1所述的基于双波束激光实时测量的星敏感器标定方法,其中根据激光发射和折返系统的安装信息计算所述激光在激光发射源坐标系和折返镜坐标系中的矢量信息包括下列步骤:
以O-X0Y0Z0表示激光发射源坐标系,O-XiYiZi表示折返光路上的第i个折返镜的折返镜坐标系,将所述激光在O-X0Y0Z0中表示为下式:
根据激光发射和折返系统的安装信息得到O-XiYiZi相对于O-Xi-1Yi-1Zi-1的转置矩阵信息,表示为:Ai,i-1;
以M表示折返前后矢量的变换矩阵,表示为:M=diag(1,1,-1);
计算所述激光在第1个折返镜反射后在O-X1Y1Z1中的矢量信息,表示为下式:
计算所述激光在第i个折返镜反射后在O-XiYiZi中的矢量信息,表示为下式:
以O-XnYnZn表示所述激光在折返进入星敏感器的视场前的最后一个折返镜的折返镜坐标系,计算所述激光在O-XnYnZn中的矢量信息,表示为下式:
5.根据权利要求4所述的基于双波束激光实时测量的星敏感器标定方法,其中根据星敏感器与激光发射源的相对姿态矩阵和光学载荷系统与激光发射源的相对姿态矩阵计算星敏感器与光学载荷系统的相对姿态矩阵包括下列步骤:
以AZ,0表示光学载荷系统与激光发射源的相对姿态矩阵,计算星敏感器与光学载荷系统的相对姿态矩阵,表示为下式:
AZ,S=AZ,0·(AS,0)T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110515754.0A CN113252074B (zh) | 2021-05-12 | 2021-05-12 | 一种基于双波束激光实时测量的星敏感器标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110515754.0A CN113252074B (zh) | 2021-05-12 | 2021-05-12 | 一种基于双波束激光实时测量的星敏感器标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113252074A CN113252074A (zh) | 2021-08-13 |
CN113252074B true CN113252074B (zh) | 2023-04-07 |
Family
ID=77222908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110515754.0A Active CN113252074B (zh) | 2021-05-12 | 2021-05-12 | 一种基于双波束激光实时测量的星敏感器标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113252074B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105157700A (zh) * | 2015-07-07 | 2015-12-16 | 中国空间技术研究院 | 用于空间光学相机的无地面控制点定位方法和装置 |
CN106918349A (zh) * | 2016-11-22 | 2017-07-04 | 上海航天控制技术研究所 | 星敏感器在轨相对安装误差标定系统及标定补偿方法 |
CN107270939A (zh) * | 2017-06-20 | 2017-10-20 | 上海航天控制技术研究所 | 一种星敏感器相对安装标定方法及系统 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7142981B2 (en) * | 2003-08-05 | 2006-11-28 | The Boeing Company | Laser range finder closed-loop pointing technology of relative navigation, attitude determination, pointing and tracking for spacecraft rendezvous |
CN101699222B (zh) * | 2009-11-02 | 2011-04-20 | 中国人民解放军国防科学技术大学 | 一种星敏感器校准装置及应用于高精度星敏感器校准的方法 |
CN103674023B (zh) * | 2013-12-26 | 2014-10-15 | 中国人民解放军国防科学技术大学 | 一种基于陀螺精确角度关联的星敏感器动态测姿方法 |
CN104792340B (zh) * | 2015-05-15 | 2017-08-25 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
CN106767673B (zh) * | 2016-11-29 | 2019-06-21 | 上海卫星工程研究所 | 一种卫星高精度光学敏感载荷的指向测量方法 |
CN106896819A (zh) * | 2017-02-23 | 2017-06-27 | 上海航天控制技术研究所 | 基于三星敏感器的卫星姿态确定方法和系统 |
CN112082574B (zh) * | 2020-09-04 | 2023-05-12 | 中国科学院微小卫星创新研究院 | 星敏感器的校正方法及系统 |
-
2021
- 2021-05-12 CN CN202110515754.0A patent/CN113252074B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105157700A (zh) * | 2015-07-07 | 2015-12-16 | 中国空间技术研究院 | 用于空间光学相机的无地面控制点定位方法和装置 |
CN106918349A (zh) * | 2016-11-22 | 2017-07-04 | 上海航天控制技术研究所 | 星敏感器在轨相对安装误差标定系统及标定补偿方法 |
CN107270939A (zh) * | 2017-06-20 | 2017-10-20 | 上海航天控制技术研究所 | 一种星敏感器相对安装标定方法及系统 |
Non-Patent Citations (2)
Title |
---|
Smith NH 等.Laser Reference Sensor Alignment Estimation Using Reference Signal Observations.《JOURNAL OF SPACECRAFT AND ROCKETS》.2014,第第51卷卷(第第1期期),第48-56页. * |
高迪.光学导航敏感器标定系统装调方法研究.《《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》》.2013,第2013年卷卷(第第03期期),C031-93. * |
Also Published As
Publication number | Publication date |
---|---|
CN113252074A (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chung et al. | ARGOS testbed: study of multidisciplinary challenges of future spaceborne interferometric arrays | |
US20110157600A1 (en) | Optical wave-front recovery for active and adaptive imaging control | |
CN104318047B (zh) | 一种激光标刻系统几何畸变的校正方法 | |
US20100053634A1 (en) | Direct solve image based wave-front sensing | |
CN108519671B (zh) | 一种拼接望远镜系统相位平移误差的闭环校正控制方法 | |
CN113252074B (zh) | 一种基于双波束激光实时测量的星敏感器标定方法 | |
CN113687521A (zh) | 基于波前校正的低像差高精度二维光电自准直方法与装置 | |
US20010037185A1 (en) | Method for determining the state variables of a moving rigid body in space | |
US20100053635A1 (en) | Direct solve image based wave-front sensing | |
CN109855652B (zh) | 星载激光测高仪指向角误差为非常数时的在轨标定方法 | |
CN113091732B (zh) | 用于卫星性能指标在轨评估的恒星捕获方法及系统 | |
US5285074A (en) | Dynamic compensation of non-linear electron beam landing angle in variable axis lenses | |
CN113447043A (zh) | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 | |
CN112097794B (zh) | 遥感卫星载荷平台标定方法及系统 | |
Burge et al. | Development of surface metrology for the Giant Magellan Telescope primary mirror | |
CN115630254A (zh) | 高精度星敏感器光学系统内参数微变化量的在轨标定方法 | |
Delabie et al. | Testing and calibrating an advanced cubesat attitude determination and control system | |
CN109506645B (zh) | 一种星敏感器安装矩阵地面精测方法 | |
CN110427686B (zh) | 一种考虑约束条件的大刚体位移参量计算方法 | |
CN112254682A (zh) | 一种大量程光电自准直仪的三维测角方法 | |
Howard et al. | Optical modeling activities for the James Webb Space Telescope (JWST) project: II. Determining image motion and wavefront error over an extended field of view with a segmented optical system | |
LeCroy et al. | Effects of optical artifacts in a laser-based spacecraft navigation sensor | |
CN115060227B (zh) | 卫星在轨视轴指向变化高精度测量的组合夹角测量方法 | |
CN113297745B (zh) | 一种基于短弧拟合位置的双弧段轨道改进方法 | |
Sun et al. | Analysis and calibration of star sensor's image plane displacement |
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 |