CN113739791A - 基于矢量描述的导航星集确定方法 - Google Patents

基于矢量描述的导航星集确定方法 Download PDF

Info

Publication number
CN113739791A
CN113739791A CN202110891377.0A CN202110891377A CN113739791A CN 113739791 A CN113739791 A CN 113739791A CN 202110891377 A CN202110891377 A CN 202110891377A CN 113739791 A CN113739791 A CN 113739791A
Authority
CN
China
Prior art keywords
vector
edge
view field
declination
circular
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
CN202110891377.0A
Other languages
English (en)
Other versions
CN113739791B (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 Aerospace Automatic Control Research Institute
Original Assignee
Beijing Aerospace Automatic Control Research Institute
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 Beijing Aerospace Automatic Control Research Institute filed Critical Beijing Aerospace Automatic Control Research Institute
Priority to CN202110891377.0A priority Critical patent/CN113739791B/zh
Publication of CN113739791A publication Critical patent/CN113739791A/zh
Application granted granted Critical
Publication of CN113739791B publication Critical patent/CN113739791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Navigation (AREA)

Abstract

本发明涉及基于矢量描述的导航星集确定方法:S1、计算观星设备视场轴向在地心惯性系下的矢量u;S2、计算圆形视场边缘起始点的矢量
Figure DDA0003196227370000012
S3、将以矢量u为中轴线,矢量
Figure DDA0003196227370000011
为母线形成的圆锥体底面边缘,作为圆形视场边缘,根据矢量计算原理获得圆形视场边缘各点矢量Vs的描述模型;S4、计算步骤S3获得的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解;S5、根据圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解,确定圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围;S6、查询恒星星表,得到赤经、赤纬落入步骤S6确定的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围内的恒星集合。本发明计算简单,且适应任意给定星敏感器轴向条件的情况。

Description

基于矢量描述的导航星集确定方法
技术领域
本发明属于天文导航技术领域,涉及基于矢量描述的导航星集确定方法。
背景技术
星光导航具有定姿精度高,自主性强,可靠性高等优点,在航空航天飞行器上得到了广泛的应用。在星光导航地面模拟试验、星图识别方法验证等方面,通常需要计算星敏感器对准某个天区内的星集。可见,传统的基于几何投影的星集计算方法在极区不适应,且计算复杂。
发明内容
本发明解决的技术问题是:克服现有技术的不足,本发明提出了基于矢量描述的导航星集确定方法,视场内导航星集计算简单,且适应任意给定星敏感器轴向条件的情况。
本发明解决技术的方案是:基于矢量描述的导航星集确定方法,该方法包括如下步骤:
S1、根据观星设备视场轴向的赤经和赤纬,计算观星设备视场轴向在地心惯性系下的矢量u;
S2、根据观星设备视场轴向的赤经和赤纬,以观星设备视场中心作为圆心,观星设备视场半径作为圆的半径r,确定圆形视场边缘起始点的赤经αj和赤纬βj,并据此计算圆形视场边缘起始点的矢量
Figure BDA0003196227350000011
S3、将以矢量u为中轴线,矢量
Figure BDA0003196227350000012
为母线形成的圆锥体底面边缘,作为圆形视场边缘,根据矢量计算原理获得圆形视场边缘各点矢量Vs的描述模型;
S4、计算步骤S3获得的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解;
S5、根据圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解,确定圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围;
S6、查询恒星星表,得到赤经、赤纬落入步骤S6确定的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围内的恒星集合。
优选地,所述观星设备包括星敏感器、天文望远镜;所述地心惯性系为J2000地心惯性系。
优选地,所述步骤S1观星设备视场轴向在地心惯性系下的矢量u的计算公式如下:
Figure BDA0003196227350000021
其中,αi为观星设备视场轴向的赤经、βi为观星设备视场轴向的赤纬。
优选地,所述步骤S2中,圆形视场边缘起始点的赤经αj和赤纬βj选择如下:
αj=αi
βj=βi-r。
优选地,所述步骤S2中,圆形视场边缘起始点的矢量
Figure BDA0003196227350000024
为:
Figure BDA0003196227350000022
优选地,所述步骤S3中,圆形视场边缘各点矢量Vs的描述模型如下:
Vs=Vj+u×Vjsinθ+(1-cosθ)u×(u×Vj)
θ=0,Δθ,…nΔθ;
Δθ为间隔,nΔθ<2π。
优选地,所述步骤S4具体如下:
S4.1、计算圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs
Figure BDA0003196227350000023
Figure BDA0003196227350000031
αs=(sinαs)-1,βs=sin-1[Vs(3)]
S4.2、求解圆形视场边缘各点矢量Vs对应的赤经αs唯一解:
Figure BDA0003196227350000032
本发明与现有技术相比的有益效果是:
(1)、本发明利用矢量计算某个视场内的导航星集,计算简单,效率高,除了可以用于地面模拟星图计算外,也可以支撑在线星图识别算法的计算。
(2)、传统的基于几何投影的方法由于存在三角函数除法运算,在极区不适应,本发明则不存在该问题,适应性更强。
附图说明
图1为本发明实施例J2000.0地心惯性坐标系定义;
图2为本发明实施例圆形视场示意图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
实施例1
本发明的技术方案具体实施方式如下所述,其示意图如图1所示。
本实施例提供的基于矢量描述的导航星集确定方法,包括如下步骤:
S1、根据观星设备视场轴向的赤经和赤纬,计算观星设备视场轴向在地心惯性系下的矢量u;所述观星设备包括星敏感器、天文望远镜,本实施例中为星敏感器,所述地心惯性系为J2000地心惯性系。
观星设备视场轴向在地心惯性系下的矢量u的计算公式如下:
Figure BDA0003196227350000033
其中,αi为观星设备视场轴向的赤经,单位为弧度;βi为观星设备视场轴向的赤纬,单位为弧度。
S2、根据观星设备视场轴向的赤经和赤纬,以观星设备视场中心作为圆心,观星设备视场半径作为圆的半径r,确定圆形视场边缘起始点的赤经αj和赤纬βj,并据此计算圆形视场边缘起始点的矢量
Figure BDA0003196227350000041
圆形视场如图2所示,圆形视场边缘起始点的赤经αj和赤纬βj选择如下:
αj=αi
βj=βi-r。
圆形视场边缘起始点的矢量
Figure BDA0003196227350000042
为:
Figure BDA0003196227350000043
S3、将以矢量u为中轴线,矢量
Figure BDA0003196227350000044
为母线形成的圆锥体底面边缘,作为圆形视场边缘,根据矢量计算原理获得圆形视场边缘各点矢量Vs的描述模型;
圆形视场边缘各点矢量Vs的描述模型如下:
Vs=Vj+u×Vjsinθ+(1-cosθ)u×(u×Vj)
θ=0,Δθ,…nΔθ;
Δθ为间隔,nΔθ<2π,Δθ推荐取值范围为0.01°~0.03°。
S4、计算步骤S3获得的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解;
S4.1、计算圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs
Figure BDA0003196227350000045
Figure BDA0003196227350000046
αs=(sinαs)-1,βs=sin-1[Vs(3)]
S4.2、由于上述反三角的值不唯一,根据cosαs和sinαs的数值,需要调整αs的值,因此,求解圆形视场边缘各点矢量Vs对应的赤经αs唯一解:
Figure BDA0003196227350000051
S5、根据圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解,确定圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围;
S6、查询恒星星表,得到赤经、赤纬落入步骤S6确定的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围内的恒星集合。所述恒星表为史密松天体台SAO恒星星表,或其它星表。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。

Claims (7)

1.基于矢量描述的导航星集确定方法,其特征在于包括如下步骤:
S1、根据观星设备视场轴向的赤经和赤纬,计算观星设备视场轴向在地心惯性系下的矢量u;
S2、根据观星设备视场轴向的赤经和赤纬,以观星设备视场中心作为圆心,观星设备视场半径作为圆的半径r,确定圆形视场边缘起始点的赤经αj和赤纬βj,并据此计算圆形视场边缘起始点的矢量
Figure FDA0003196227340000012
S3、将以矢量u为中轴线,矢量
Figure FDA0003196227340000013
为母线形成的圆锥体底面边缘,作为圆形视场边缘,根据矢量计算原理获得圆形视场边缘各点矢量Vs的描述模型;
S4、计算步骤S3获得的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解;
S5、根据圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs唯一解,确定圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围;
S6、查询恒星星表,得到赤经、赤纬落入步骤S6确定的圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs范围内的恒星集合。
2.根据权利要求1所述的基于矢量描述的导航星集确定方法,其特征在于所述观星设备包括星敏感器、天文望远镜;所述地心惯性系为J2000地心惯性系。
3.根据权利要求1所述的基于矢量描述的导航星集确定方法,其特征在于所述步骤S1观星设备视场轴向在地心惯性系下的矢量u的计算公式如下:
Figure FDA0003196227340000011
其中,αi为观星设备视场轴向的赤经、βi为观星设备视场轴向的赤纬。
4.根据权利要求1所述的基于矢量描述的导航星集确定方法,其特征在于所述步骤S2中,圆形视场边缘起始点的赤经αj和赤纬βj选择如下:
αj=αi
βj=βi-r。
5.根据权利要求1所述的基于矢量描述的导航星集确定方法,其特征在于所述步骤S2中,圆形视场边缘起始点的矢量
Figure FDA0003196227340000021
为:
Figure FDA0003196227340000022
6.根据权利要求1所述的基于矢量计算的导航星集确定方法,其特征在于所述步骤S3中,圆形视场边缘各点矢量Vs的描述模型如下:
Vs=Vj+u×Vjsinθ+(1-cosθ)u×(u×Vj)
θ=0,Δθ,…nΔθ;
Δθ为间隔,nΔθ<2π。
7.根据权利要求1所述的基于矢量描述的导航星集确定方法,其特征在于所述步骤S4具体如下:
S4.1、计算圆形视场边缘各点矢量Vs对应的赤经αs、赤纬βs
Figure FDA0003196227340000023
Figure FDA0003196227340000024
αs=(sinαs)-1,βs=sin-1[Vs(3)]
S4.2、求解圆形视场边缘各点矢量Vs对应的赤经αs唯一解:
Figure FDA0003196227340000025
CN202110891377.0A 2021-08-04 2021-08-04 基于矢量描述的导航星集确定方法 Active CN113739791B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110891377.0A CN113739791B (zh) 2021-08-04 2021-08-04 基于矢量描述的导航星集确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110891377.0A CN113739791B (zh) 2021-08-04 2021-08-04 基于矢量描述的导航星集确定方法

Publications (2)

Publication Number Publication Date
CN113739791A true CN113739791A (zh) 2021-12-03
CN113739791B CN113739791B (zh) 2023-07-14

Family

ID=78730054

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110891377.0A Active CN113739791B (zh) 2021-08-04 2021-08-04 基于矢量描述的导航星集确定方法

Country Status (1)

Country Link
CN (1) CN113739791B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001046648A1 (fr) * 1999-12-21 2001-06-28 Hidemi Miura Systeme de navigation dote d'une carte du ciel
US20130013199A1 (en) * 2011-07-06 2013-01-10 Zheng You Method for measuring precision of star sensor and system using the same
CN109579829A (zh) * 2018-11-29 2019-04-05 天津津航技术物理研究所 一种小视场星敏感器短波导航星识别方法
CN109724591A (zh) * 2019-02-27 2019-05-07 北京航空航天大学 一种基于角距筛选的星敏感器假星筛除方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001046648A1 (fr) * 1999-12-21 2001-06-28 Hidemi Miura Systeme de navigation dote d'une carte du ciel
US20130013199A1 (en) * 2011-07-06 2013-01-10 Zheng You Method for measuring precision of star sensor and system using the same
CN109579829A (zh) * 2018-11-29 2019-04-05 天津津航技术物理研究所 一种小视场星敏感器短波导航星识别方法
CN109724591A (zh) * 2019-02-27 2019-05-07 北京航空航天大学 一种基于角距筛选的星敏感器假星筛除方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
BING LI等: "A Neural Network Algorithm For Star Identification" *
MENG NA: "Modified Grid Algorithm for Noisy All-Sky Autonomous Star Identification" *
朱长征;张志龙;沈振康;: "一种新的星模式识别算法" *
李冰;高晓颖;踪华;李林浩;: "可适应多维度信息变化的智能星图识别方法" *
袁晓波;张超;: "大视场星敏感器矢量观测在轨检校模型" *
踪华等: "一种基于模式匹配的自主星图识别算法" *

Also Published As

Publication number Publication date
CN113739791B (zh) 2023-07-14

Similar Documents

Publication Publication Date Title
CN104406583B (zh) 双星敏感器联合确定载体姿态方法
CN106197425B (zh) 基于卫星姿态角的地面目标点位置的计算方法
CN107806874B (zh) 一种视觉辅助的捷联惯导极区初始对准方法
CN108871326B (zh) 一种单轴旋转调制惯性-天文深组合导航方法
CN101750067B (zh) 一种成像式地球敏感器地球扁率修正方法
CN113108783B (zh) 一种无人潜航器惯性/多普勒组合导航方法
CN112880669B (zh) 一种航天器星光折射和单轴旋转调制惯性组合导航方法
CN102901485B (zh) 一种光电经纬仪快速自主定向的方法
CN105910623B (zh) 利用磁强计辅助gnss/mins紧组合系统进行航向校正的方法
CN110986934B (zh) 一体化双轴旋转惯导天文组合导航系统的导航方法及系统
CN112179334A (zh) 基于两步Kalman滤波的星光导航方法及系统
CN103487032A (zh) 低轨空间相机自由指向像移矢量计算方法
CN109407696B (zh) 一种无人机航向角动态校定方法
CN113739791A (zh) 基于矢量描述的导航星集确定方法
CN111879299B (zh) 一种地基望远镜全自动指星方法
CN112284412A (zh) 一种避免欧拉转换奇异导致精度下降的地面静态对准方法
CN109099911B (zh) 一种航空系统导航定位的方法及系统
CN113218390B (zh) 基于姿态和恒星高度角融合的旋转惯性天文组合导航方法
CN113589343B (zh) 一种基于月球成像敏感器的月心矢量和太阳方向提取方法
CN115510378A (zh) 地球静止轨道标称半长轴摄动计算方法
US3439884A (en) Space vehicle guidance system
Zhan et al. Adaptive celestial positioning for the stationary mars rover based on a self-calibration model for the star sensor
CN114353784A (zh) 一种基于运动矢量的制导炮弹空中姿态辨识方法
CN116559917B (zh) 一种定位海面动目标的被动电侦卫星编队构型设计方法
CN114485624A (zh) 一种基于恒星与卫星联合的全天时光学导航方法和装置

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