CN107861919A - 一种星敏感器定姿星矢量权值实时分配方法 - Google Patents

一种星敏感器定姿星矢量权值实时分配方法 Download PDF

Info

Publication number
CN107861919A
CN107861919A CN201710864352.5A CN201710864352A CN107861919A CN 107861919 A CN107861919 A CN 107861919A CN 201710864352 A CN201710864352 A CN 201710864352A CN 107861919 A CN107861919 A CN 107861919A
Authority
CN
China
Prior art keywords
msub
mrow
msubsup
mtd
star
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
CN201710864352.5A
Other languages
English (en)
Other versions
CN107861919B (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.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology 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 Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN201710864352.5A priority Critical patent/CN107861919B/zh
Publication of CN107861919A publication Critical patent/CN107861919A/zh
Application granted granted Critical
Publication of CN107861919B publication Critical patent/CN107861919B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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
    • G01C21/025Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
    • 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/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Automation & Control Theory (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Astronomy & Astrophysics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Acyclic And Carbocyclic Compounds In Medicinal Compositions (AREA)
  • Navigation (AREA)

Abstract

一种星敏感器定姿星矢量权值实时分配方法,包括以下步骤:S1、根据定姿星观测矢量Wi,参考矢量Vi和星矢量权值其中i=1,2,…,n,n为定姿星颗数,首次解算先验姿态四元数S2、根据首次解算的先验姿态四元数解算先验姿态矩阵A;S3、根据定姿星观测矢量Wi,参考矢量Vi和姿态矩阵A计算星矢量误差;S4、根据矢量误差计算动态权值ai;S5、根据定姿星观测矢量Wi,参考矢量Vi和动态权值ai二次迭代解算姿态四元数q。

Description

一种星敏感器定姿星矢量权值实时分配方法
技术领域
本发明涉及一种高精度星敏感器定姿星矢量权值实时分配方法。
背景技术
星敏感器是基于计算机视觉原理的光电测量仪器,其直接测量的物理量为恒星经光学系统在光电探测器上形成的近圆形光斑的质心,经质心提取、星矢量计算、星图匹配、姿态解算等步骤,最终得到卫星姿轨控系统需要的四元数。姿态四元数用于描述飞行器的姿态信息,准确的姿态矩阵A应满足公式:
式中:n为定姿星颗数,是星敏感器测量坐标系下的观测星矢量;是天球坐标系下的参考星矢量。
从星图成像到输出姿态,引入了各种误差,因此获得星敏感器姿态矩阵的准确值是不可能的,当前的求解思路大多基于1965年Wahba提出的最小损失函数:
式中:ai是恒星矢量的权值。
当前姿态解算算法中,大多以定姿星颗数为依据,采用等权值的方法解算姿态,由于星敏感器装调过程中产生的误差(比如像平面倾斜)及光学系统像差(球差、慧差、象散、场曲、位置色差、倍率色差、畸变),以及恒星星等不一,质心算法提取的恒星矢量即使经过一系列处理,仍存在大小不一的残差。
发明内容
本发明的技术解决问题是:
克服现有技术的缺陷,提出一种计算量小,实现简单的方法,在轨飞行过程中实时估计观测星矢量的误差,以动态分配定姿星矢量的权值,提高姿态解算的精度,适用于高精度星敏感器的定姿星矢量权值的确定。
本发明的技术解决方案是:
一种星敏感器定姿星矢量权值实时分配方法,包括以下步骤:
S1、根据定姿星观测矢量Wi,参考矢量Vi和星矢量权值其中i=1,2,…,n,n为定姿星颗数,首次解算先验姿态四元数
S2、根据首次解算的先验姿态四元数解算先验姿态矩阵A;
S3、根据定姿星观测矢量Wi,参考矢量Vi和姿态矩阵A计算星矢量误差
S4、根据矢量误差计算动态权值ai
S5、根据定姿星观测矢量Wi,参考矢量Vi和动态权值ai二次迭代解算姿态四元数q。
根据本发明的一个实施例,在步骤S1中,根据定姿星观测矢量Wi,参考矢量Vi,和星矢量权值使用QUEST方法解算先验姿态四元数
根据本发明的一个实施例,在步骤S2中,先验姿态矩阵
根据本发明的一个实施例,在步骤S3中,根据等式来计算星矢量误差
根据本发明的一个实施例,在步骤S4中,根据等式来计算动态分配权值ai
根据本发明的一个实施例,在步骤S5中,根据定姿星观测矢量Wi,参考矢量Vi和动态权值ai使用QUEST方法进行二次迭代解算姿态四元数q。
一种存储了指令的计算机可读介质,当所述指令被至少一个处理器执行时,使得所述至少一个处理器执行以下步骤:
S1、根据定姿星观测矢量Wi,参考矢量Vi和星矢量权值其中i=1,2,…,n,n为定姿星颗数,使用QUEST方法首次解算先验姿态四元数
S2、根据首次解算的先验姿态四元数解算先验姿态矩阵A;
S3、根据定姿星观测矢量Wi,参考矢量Vi和姿态矩阵A计算星矢量误差
S4、根据矢量误差计算动态权值ai
S5、根据定姿星观测矢量Wi,参考矢量Vi和动态分配权值ai使用QUEST方法进行二次迭代解算姿态四元数q。
根据本发明的一个实施例,在步骤S2中,先验姿态矩阵
根据本发明的一个实施例,在步骤S3中,根据等式来计算星矢量误差
根据本发明的一个实施例,在步骤S4中,根据等式来计算动态权值ai
本发明与现有技术相比的优点在于:
(1)引入二次迭代,以解决无法直接准确计算每颗星矢量的误差问题;
(2)采用三维的星矢量信息计算误差,相对于转到二维焦平面内计算误差的方法,省去了将恒星矢量从天球坐标系反解到像平面坐标系的过程,节约了计算资源;
(3)采用动态权值的方法,相较于等权值解算姿态,可显著抑制光学系统像差等引起的误差。
附图说明
图1是高精度星敏感器定姿星矢量权值实时分配方法的流程图。
图2是单头星敏感器中的三轴姿态示意图。
图3是多头星敏感器中的三轴姿态示意图。
具体实施方式
下面结合附图,具体描述本发明的实施例。
如图1所示,高精度星敏感器定姿星矢量权值实时分配方法包括以下步骤:
S1、根据定姿星观测矢量Wi,参考矢量Vi和星矢量权值其中i=1,2,…,n,n为定姿星颗数,首次解算先验姿态四元数
在本发明的实施例中,将定姿星观测矢量Wi,参考矢量Vi,和星矢量权值作为QUEST方法的输入,来解算先验姿态四元数
S2、根据首次解算的先验姿态四元数解算先验姿态矩阵A。
在本发明的实施例中,先验姿态矩阵A如以下等式(3)所示:
S3、根据定姿星观测矢量Wi,参考矢量Vi和姿态矩阵A计算星矢量误差
在本发明的实施例中,是根据以下等式(4)来计算的:
S4、根据矢量误差计算动态分配权值ai
在本发明的实施例中,在计算出n个矢量误差后,计算n个矢量误差的倒数的和并根据以下等式(5)来计算动态权值ai
S5、根据定姿星观测矢量Wi,参考矢量Vi和动态权值ai二次迭代解算姿态四元数q。
在本发明的实施例中,将定姿星观测矢量Wi,参考矢量Vi和动态分配权值ai作为QUEST算法的输入,二次迭代解算姿态四元数q。
实施例
如图2所示,在单头星敏感器中,采用本发明的方法,解算姿态四元数q,再将姿态四元数q应用于敏感器,解算得到的三轴姿态相比于现有技术中的使用等权值解算得到的三轴姿态,在X,Y轴上精度均有提高。
如图3所示,在多头星敏感器中,采用本发明的方法,解算姿态四元数q,再将姿态四元数q应用于敏感器,解算得到的三轴姿态相比于现有技术中的使用等权值解算得到的三轴姿态,在X,Y,Z轴上精度均有提高。
综上所述,根据定姿星观测矢量的误差大小,来确定权值ai的动态分配的方法,可以有效提高姿态确定精度。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

Claims (10)

1.一种星敏感器定姿星矢量权值实时分配方法,其特征在于,包括以下步骤:
S1、根据定姿星观测矢量Wi,参考矢量Vi和星矢量权值其中i=1,2,…,n,n为定姿星颗数,首次解算先验姿态四元数
S2、根据首次解算的先验姿态四元数解算先验姿态矩阵A;
S3、根据定姿星观测矢量Wi,参考矢量Vi和姿态矩阵A计算星矢量误差S4、根据矢量误差计算动态权值ai
S5、根据定姿星观测矢量Wi,参考矢量Vi和动态权值ai二次迭代解算姿态四元数q。
2.根据权利要求1所述的星敏感器定姿星矢量权值实时分配方法,其特征在于,在步骤S1中,根据定姿星观测矢量Wi,参考矢量Vi,和星矢量权值使用QUEST方法解算先验姿态四元数
3.根据权利要求2所述的星敏感器定姿星矢量权值实时分配方法,其特征在于,在步骤S2中,先验姿态矩阵
<mrow> <mi>A</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
4.根据权利要求1所述的星敏感器定姿星矢量权值实时分配方法,其特征在于,在步骤S3中,根据等式来计算星矢量误差
5.根据权利要求1所述的星敏感器定姿星矢量权值实时分配方法,其特征在于,在步骤S4中,根据等式来计算动态权值权值ai
6.根据权利要求1所述的星敏感器定姿星矢量权值实时分配方法,其特征在于,在步骤S5中,根据定姿星观测矢量Wi,参考矢量Vi和动态分配权值ai使用QUEST方法进行二次迭代解算姿态四元数q。
7.一种存储了指令的计算机可读介质,其特征在于,当所述指令被至少一个处理器执行时,使得所述至少一个处理器执行以下步骤:
S1、根据定姿星观测矢量Wi,参考矢量Vi和星矢量权值其中i=1,2,…,n,n为定姿星颗数,使用QUEST方法首次解算先验姿态四元数
S2、根据首次解算的先验姿态四元数解算先验姿态矩阵A;
S3、根据定姿星观测矢量Wi,参考矢量Vi和姿态矩阵A计算星矢量误差S4、根据矢量误差计算动态权值ai
S5、根据定姿星观测矢量Wi,参考矢量Vi和动态分配权值ai使用QUEST方法进行二次迭代解算姿态四元数q。
8.根据权利要求7所述的存储了指令的计算机可读介质,其特征在于,在步骤S2中,先验姿态矩阵
<mrow> <mi>A</mi> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <msub> <mi>q</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msubsup> <mi>q</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>q</mi> <mn>2</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>q</mi> <mn>4</mn> <mn>2</mn> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
9.根据权利要求7所述的存储了指令的计算机可读介质,其特征在于,在步骤S3中,根据等式来计算星矢量误差
10.根据权利要求7所述的存储了指令的计算机可读介质,其特征在于,在步骤S4中,根据等式来计算动态权值ai
CN201710864352.5A 2017-09-22 2017-09-22 一种星敏感器定姿星矢量权值实时分配方法 Active CN107861919B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710864352.5A CN107861919B (zh) 2017-09-22 2017-09-22 一种星敏感器定姿星矢量权值实时分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710864352.5A CN107861919B (zh) 2017-09-22 2017-09-22 一种星敏感器定姿星矢量权值实时分配方法

Publications (2)

Publication Number Publication Date
CN107861919A true CN107861919A (zh) 2018-03-30
CN107861919B CN107861919B (zh) 2021-03-26

Family

ID=61698276

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710864352.5A Active CN107861919B (zh) 2017-09-22 2017-09-22 一种星敏感器定姿星矢量权值实时分配方法

Country Status (1)

Country Link
CN (1) CN107861919B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109302256A (zh) * 2018-12-05 2019-02-01 上海航天控制技术研究所 一种星敏感器时间同步信号校时方法和逻辑电路
CN109489689A (zh) * 2018-11-20 2019-03-19 上海航天控制技术研究所 一种基于α-β滤波的星矢量测量误差在轨估计方法
CN109945854A (zh) * 2019-03-26 2019-06-28 上海微小卫星工程中心 一种基于四元数差值融合星敏感器数据的方法
CN111572817A (zh) * 2020-06-08 2020-08-25 北京航天自动控制研究所 一种平台星光修正系数优化计算方法
CN111637879A (zh) * 2020-04-22 2020-09-08 北京控制工程研究所 基于多维差异化误差特性分布的双星敏加权姿态确定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288201A (zh) * 2011-07-06 2011-12-21 清华大学 用于星敏感器的精度测量方法
CN102288200A (zh) * 2011-07-06 2011-12-21 清华大学 用于星敏感器的精度测量系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288201A (zh) * 2011-07-06 2011-12-21 清华大学 用于星敏感器的精度测量方法
CN102288200A (zh) * 2011-07-06 2011-12-21 清华大学 用于星敏感器的精度测量系统

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
凌林本等: "多传感器数据融合时权的最优分配原则", 《中国惯性技术学报》 *
吕沧海等: "《中远程导弹组合导航技术》", 31 March 2014, 国防工业出版社 *
江杰等: "权值动态分配的加权平均法在姿态结算中的应用", 《航空制造技术》 *
王献忠等: "基于地平仪的惯性-天文组合导航", 《导航定位学报》 *
赵剡等: "《高精度卫星导航技术》", 30 April 2016, 北京航空航天大学出版社 *
邢飞等: "《APS CMOS 星敏感器系统原理及实现方法》", 30 April 2017, 国防工业出版社 *
郑循江等: "一种甚高精度星敏感器精度测试方法", 《红外与激光工程》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109489689A (zh) * 2018-11-20 2019-03-19 上海航天控制技术研究所 一种基于α-β滤波的星矢量测量误差在轨估计方法
CN109489689B (zh) * 2018-11-20 2020-11-03 上海航天控制技术研究所 一种基于α-β滤波的星矢量测量误差在轨估计方法
CN109302256A (zh) * 2018-12-05 2019-02-01 上海航天控制技术研究所 一种星敏感器时间同步信号校时方法和逻辑电路
CN109302256B (zh) * 2018-12-05 2020-03-24 上海航天控制技术研究所 一种星敏感器时间同步信号校时方法和逻辑电路
CN109945854A (zh) * 2019-03-26 2019-06-28 上海微小卫星工程中心 一种基于四元数差值融合星敏感器数据的方法
CN111637879A (zh) * 2020-04-22 2020-09-08 北京控制工程研究所 基于多维差异化误差特性分布的双星敏加权姿态确定方法
CN111637879B (zh) * 2020-04-22 2021-10-01 北京控制工程研究所 基于多维差异化误差特性分布的双星敏加权姿态确定方法
CN111572817A (zh) * 2020-06-08 2020-08-25 北京航天自动控制研究所 一种平台星光修正系数优化计算方法

Also Published As

Publication number Publication date
CN107861919B (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
CN107861919A (zh) 一种星敏感器定姿星矢量权值实时分配方法
CN104061929B (zh) 共光路多视场星敏感器及其星体姿态测量方法
Zhou et al. Novel autonomous on-orbit calibration method for star sensors
CN104422425B (zh) 一种不规则外形物体空间姿态动态测量方法
CN106767673B (zh) 一种卫星高精度光学敏感载荷的指向测量方法
CN106248014A (zh) 一种基于单相片的三维坐标测量方法及装置
CN108871373B (zh) 一种基于俯仰滚转台和非线性优化的星敏感器标定方法
CN105526951B (zh) 一种星敏感器原始观测数据预处理方法及系统
CN107869960B (zh) 一种飞机机翼变形测量装置、安装方法及测试方法
CN107588785B (zh) 一种考虑像点误差的星敏感器内外参数简化标定方法
CN109341720A (zh) 一种基于恒星轨迹的遥感相机几何标定方法
CN106679695B (zh) 一种基于星敏感器的焦平面热变形的测试方法
CN113218418A (zh) 一种空间极高精度指向测量仪器热光耦合效应确定系统及方法
CN110686684B (zh) 一种小天体环绕探测器光学协同定轨方法
CN105547286B (zh) 一种复合三视场星敏感器星图模拟方法
CN100348947C (zh) 基于weng模型的星敏感器在轨校准方法
CN112097794B (zh) 遥感卫星载荷平台标定方法及系统
CN110108271B (zh) 一种气动光学效应引起的星光偏折补偿方法
Jiang et al. Optical image generation and high-precision line-of-sight extraction for Mars approach navigation
CN113022898A (zh) 弱引力环境柔性附着系统状态估计方法
Enright et al. Rolling shutter compensation for star trackers
CN109002567A (zh) 基于实际光学系统装调检测光路的计算机辅助装调方法
CN114034470B (zh) 望远镜波前旋转角度的计算方法、装置及望远镜
EP1992907A2 (en) Method for generating pseudo-random pattern designs for optical position determination systems
CN115630254A (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