CN101696885A - 一种提高星敏感器数据处理精度的方法 - Google Patents
一种提高星敏感器数据处理精度的方法 Download PDFInfo
- Publication number
- CN101696885A CN101696885A CN200910044703A CN200910044703A CN101696885A CN 101696885 A CN101696885 A CN 101696885A CN 200910044703 A CN200910044703 A CN 200910044703A CN 200910044703 A CN200910044703 A CN 200910044703A CN 101696885 A CN101696885 A CN 101696885A
- Authority
- CN
- China
- Prior art keywords
- error
- star
- star sensor
- satellite
- attitude
- 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
Links
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种提高星敏感器数据处理精度的方法。将星敏感器测量误差分为三大类:A星敏感器内部变形、星敏感器安装位置改变等因素引起的低频误差;B卫星轨道参数变化、星图变换及卫星在轨运动中的抖动振动等因素引起的中频误差;C星敏感器的高频误差,即服从零均值的高斯白噪声分布的随机误差;经过星敏感器测量误差求解后,姿态估计的偏差已经基本消除,其定姿精度提高了15%以上。
Description
技术领域
本发明涉及一种提高星敏感器数据处理精度的方法。
背景技术
目前,所考虑的星敏感器单位方向矢量测量模型如下
其中r表示在参考坐标系中的某参考矢量,如由星历表给出的某恒星矢量(表示在惯性坐标系中);b表示在星敏感器坐标系中的测量矢量;A为星体坐标系相对参考坐标系的姿态矩阵;Sbs为星敏感器的安装矩阵,即由星敏感器坐标系到星体坐标系的坐标转换矩阵;v为与星敏感器内部设计有关的测量随机误差。然而,星敏感器测量误差不仅包括与星敏感器内部设计有关的误差,还包括与星敏感器安装、在轨环境、星体抖动等因素有关的外部因素引起的外部误差。目前,对于星敏感器测量误差的分析和处理主要集中在与星敏感器内部设计有关的基本误差因素分类及精度估算关系上,很少涉及星敏感器外部综合的误差分析和处理。在星敏感器测量过程中,误差产生机理的不确定性,误差表现形式的多样性,误差关系的复杂性以及要求误差处理的准确性和可靠性都给星敏感器测量误差的特性分析及处理带来了很大的难度。因而造成目前星敏感器测量误差处理效果不理想,进而数据处理精度不高。
发明内容
本发明的目的在于提供一种提高星敏感器数据处理精度的方法,以实现不但能对星敏感内部设计有关的测量误差进行分析与处理,还能对除星敏感器内部误差以外的其他误差进行分析和处理。
本发明的技术方案包括以下步骤:
(1)建立星敏感器姿态测量误差模型:
将星敏感器测量误差分为三大类:一是星敏感器内部变形、星敏感器安装位置改变等因素引起的低频误差;二是由卫星轨道参数变化、星图变换及卫星在轨运动中的抖动振动等因素引起的中频误差;三是星敏感器的高频误差,即服从零均值的高斯白噪声分布的高频误差;三种误差的矫正方法如下:
①低频误差包括二个部分:
(A)安装矩阵误差:其表现形式为:
(B)常值***误差:其表现形式为:
C=[Δs2 Δs2 Δs2]T (3)
②中频误差,其表现形式为:
χ=Bsin(2πft+ζ) (4)
(4)式中,B为中频误差的幅度,f为误差频率,ζ为误差的初相位,且该误差体现在带有误差的安装矩阵上;
③高频误差,其分布函数为:
(2)建立星敏感器姿态测量模型:
把姿态测量数据中所涉及到的低频误差、中频误差和高频误差叠加到传统的星敏感器单位方向矢量测量模型中,得到重构后的姿态测量模型如下表示
(3)结合定姿过程进行误差求解:
利用扩展卡尔曼滤波算法进行测量误差模型参数和卫星姿态的估计,待估参数取为X=(Θ,Δs1,Δs2,χ,σ2),其中Θ为姿态参数。而其中针对步骤1中的姿态测量数据误差模型,设计状态方程,假设低频误差的稳态性和中频误差的时变性,则相应的误差模型参数的状态方程为:
本发明实现了不但能对星敏感内部设计有关的测量误差进行分析与处理,还能对除星敏感器内部误差以外的其他误差进行分析和处理。
具体实施方式
实施例:
(1)建立星敏感器姿态测量误差模型:
将星敏感器测量误差分为三大类:一是星敏感器内部变形、星敏感器安装位置改变等因素引起的低频误差;二是由卫星轨道参数变化、星图变换及卫星在轨运动中的抖动振动等因素引起的中频误差;三是星敏感器的高频误差,即服从零均值的高斯白噪声分布的随机误差;星敏感器测量误差模拟仿真:
假设星敏感器的光轴z_s1和z_s2在本体坐标系中的坐标分别为
(2)根据测量误差源的分析,加入如下各类误差:
①星敏感器高频误差v包括σv=12″的随机误差;
②星敏感器低频误差包括Δs1=3″的安装矩阵误差(乘性噪声)和Δs2=2″的常值***误差(加性噪声);
③中频误差χ=Bsin(2πft+ζ)中,最大振幅B=3″,安装点处的误差频率f=10-2Hz;
(3)星敏感器姿态测量误差的处理:
按照卫星轨道参数为:地球赤道半径Re=6.378×106m,卫星轨道长半轴a=Re+0.6452×106m,轨道倾角i=97.938°,升交点赤经Ω=77°,近心点角距ω=50°,初始时刻平近心点角距M=10°,轨道角速度
利用重构后的星敏感器姿态测量方程和陀螺测量模型,结合误差四元数状态方程和各类模型误差参数的状态方程,利用扩展卡尔曼滤波算法,在卫星姿态确定的过程中,同时求解误差模型参数和卫星姿态参数;
陀螺的精度为:测量噪声均方差σg=0.05°/h,陀螺相关漂移白噪声σd=0.1°/h,相关漂移初始值为d(0)=[0.1 0.1 0.1]T°/h,时间相关常数τx=τy=τz=1h;陀螺常值漂移b=[1 -1 1]T°/h,常值漂移白噪声均方差σb=0.03°/h;滤波时长为3500s,星敏感器采样频率为4Hz。
利用仿真的星敏感器观测数据,滤波过程中不进行误差参数求解,则对应的卫星姿态估计误差见表1:
而利用扩展卡尔曼滤波进行误差参数求解,所得到的定姿结果见表2。
表1不进行误差参数求解下的卫星姿态确定结果
表2星敏感器测量误差求解后的定姿结果
从上表可以看出,经过星敏感器测量误差求解后,姿态估计的偏差已经基本消除,且定姿精度提高了15%以上。
Claims (1)
1.一种提高星敏感器数据处理精度的方法,其特征在于,包括以下步骤:
(一)建立星敏感器姿态测量误差模型:
将星敏感器测量误差分为三大类:A星敏感器内部变形、星敏感器安装位置改变因素
引起的低频误差;B卫星轨道参数变化、星图变换及卫星在轨道运动中的抖动振动因素引起的中频误差;C星敏感器的高频误差,即服从零均值的高斯白噪声分布的随机误差;
(二)根据测量误差源的分析,加入上述低频误差、中频误差、高频误差进行矫正:
①低频误差包括Δs1=3″的安装矩阵误差和Δs2=2″的常值***误差;
②中频误差χ=Bsin(2πft+ζ)中,最大振幅B=3″,安装点处的误差频率f=10-2Hz;
③高频误差v包括σv=12″的随机误差;
(三)星敏感器姿态测量误差的处理:
(1)按照卫星轨道参数为:地球赤道半径Re=6.378×106m,卫星轨道长半轴a=Re+0.6452×106m,轨道倾角i=97.938°,升交点赤经Ω=77°,近心点角距ω=50°,初始时刻平近心点角距M=10°,轨道角速度
(2)利用重构后的星敏感器姿态测量方程和陀螺测量模型,结合现有的误差四元数状态方程和所构建的各类模型误差参数的状态方程,利用扩展卡尔曼滤波算法,在卫星姿态确定的过程中,同时求解误差模型参数和卫星姿态参数;
其中:陀螺的精度为:测量噪声均方差σg=0.05°/h,陀螺相关漂移白噪声σd=0.1°/h,相关漂移初始值为d(0)=[0.10.10.1]T°/h,时间相关常数τx=τy=τz=1h;陀螺常值漂移b=[1 -1 1]T°/h,常值漂移白噪声均方差σb=0.03°/h;滤波时长为3500s,星敏感器采样频率为4Hz;
(3)利用仿真的星敏感器观测数据,滤波过程中利用扩展卡尔曼滤波进行误差参数求解,根据对应的卫星姿态估计误差得出定姿结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100447033A CN101696885B (zh) | 2009-11-05 | 2009-11-05 | 一种提高星敏感器数据处理精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100447033A CN101696885B (zh) | 2009-11-05 | 2009-11-05 | 一种提高星敏感器数据处理精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101696885A true CN101696885A (zh) | 2010-04-21 |
CN101696885B CN101696885B (zh) | 2012-05-23 |
Family
ID=42142001
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100447033A Expired - Fee Related CN101696885B (zh) | 2009-11-05 | 2009-11-05 | 一种提高星敏感器数据处理精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101696885B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176159A (zh) * | 2011-02-28 | 2011-09-07 | 哈尔滨工业大学 | 一种基于状态观测器和等价空间的卫星姿控***故障诊断装置及方法 |
CN102435763A (zh) * | 2011-09-16 | 2012-05-02 | 中国人民解放军国防科学技术大学 | 一种基于星敏感器的航天器姿态角速度测量方法 |
CN102506893A (zh) * | 2011-09-29 | 2012-06-20 | 北京控制工程研究所 | 一种基于地标信息的星敏感器低频误差补偿方法 |
CN102735260A (zh) * | 2012-06-18 | 2012-10-17 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
WO2013004033A1 (zh) * | 2011-07-06 | 2013-01-10 | 清华大学 | 用于星敏感器的精度测量方法及*** |
CN103308074A (zh) * | 2013-06-28 | 2013-09-18 | 上海新跃仪表厂 | 一种基于双星敏在轨数据的精度分析方法 |
CN103309348A (zh) * | 2013-06-28 | 2013-09-18 | 哈尔滨工业大学 | 一种利用二阶卡尔曼滤波算法估计卫星姿态控制***执行机构加性故障大小的方法 |
CN103983278A (zh) * | 2014-05-19 | 2014-08-13 | 中国人民解放军国防科学技术大学 | 一种测量影响卫星姿态确定***精度的方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航***星地联合标定与校正的方法 |
CN106289238A (zh) * | 2016-08-02 | 2017-01-04 | 哈尔滨工业大学 | 一种星敏感器数据优化处理方法 |
CN106568436A (zh) * | 2016-10-26 | 2017-04-19 | 上海航天控制技术研究所 | 一种利用陀螺修正星敏感器轨道周期***误差的方法 |
CN106708072A (zh) * | 2016-11-30 | 2017-05-24 | 上海卫星工程研究所 | 一种天基望远镜高精度姿态确定与控制方法 |
CN107608065A (zh) * | 2017-10-18 | 2018-01-19 | 中国科学院光电技术研究所 | 一种运动平台下基于allan方差分析的望远镜自动指向修正方法 |
CN108195403A (zh) * | 2017-12-28 | 2018-06-22 | 中国人民解放军国防科技大学 | 星敏感器在轨姿态测量数据综合误差模型的构建方法及装置 |
CN112781616A (zh) * | 2020-12-24 | 2021-05-11 | 中国人民解放军国防科技大学 | 星敏感器在轨测量低频误差分析方法、装置和存储介质 |
CN112989637A (zh) * | 2021-05-06 | 2021-06-18 | 中国人民解放军国防科技大学 | 基于分布折算、近似和综合的星敏***可靠性评估方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100575877C (zh) * | 2007-12-12 | 2009-12-30 | 南京航空航天大学 | 基于多信息融合的航天器组合导航方法 |
CN101214860A (zh) * | 2007-12-26 | 2008-07-09 | 北京控制工程研究所 | 轨控过程中自主选取定姿方式的方法 |
-
2009
- 2009-11-05 CN CN2009100447033A patent/CN101696885B/zh not_active Expired - Fee Related
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176159B (zh) * | 2011-02-28 | 2013-05-01 | 哈尔滨工业大学 | 一种基于状态观测器和等价空间的卫星姿控***故障诊断装置及方法 |
CN102176159A (zh) * | 2011-02-28 | 2011-09-07 | 哈尔滨工业大学 | 一种基于状态观测器和等价空间的卫星姿控***故障诊断装置及方法 |
WO2013004033A1 (zh) * | 2011-07-06 | 2013-01-10 | 清华大学 | 用于星敏感器的精度测量方法及*** |
CN102435763B (zh) * | 2011-09-16 | 2013-10-09 | 中国人民解放军国防科学技术大学 | 一种基于星敏感器的航天器姿态角速度测量方法 |
CN102435763A (zh) * | 2011-09-16 | 2012-05-02 | 中国人民解放军国防科学技术大学 | 一种基于星敏感器的航天器姿态角速度测量方法 |
CN102506893A (zh) * | 2011-09-29 | 2012-06-20 | 北京控制工程研究所 | 一种基于地标信息的星敏感器低频误差补偿方法 |
CN102735260A (zh) * | 2012-06-18 | 2012-10-17 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
CN102735260B (zh) * | 2012-06-18 | 2015-02-11 | 航天东方红卫星有限公司 | 一种星敏感器在轨测量误差的确定方法 |
CN103308074B (zh) * | 2013-06-28 | 2015-11-18 | 上海新跃仪表厂 | 一种基于双星敏在轨数据的精度分析方法 |
CN103309348A (zh) * | 2013-06-28 | 2013-09-18 | 哈尔滨工业大学 | 一种利用二阶卡尔曼滤波算法估计卫星姿态控制***执行机构加性故障大小的方法 |
CN103308074A (zh) * | 2013-06-28 | 2013-09-18 | 上海新跃仪表厂 | 一种基于双星敏在轨数据的精度分析方法 |
CN103983278B (zh) * | 2014-05-19 | 2016-07-27 | 中国人民解放军国防科学技术大学 | 一种测量影响卫星姿态确定***精度的方法 |
CN103983278A (zh) * | 2014-05-19 | 2014-08-13 | 中国人民解放军国防科学技术大学 | 一种测量影响卫星姿态确定***精度的方法 |
CN104792340B (zh) * | 2015-05-15 | 2017-08-25 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航***星地联合标定与校正的方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航***星地联合标定与校正的方法 |
CN106289238A (zh) * | 2016-08-02 | 2017-01-04 | 哈尔滨工业大学 | 一种星敏感器数据优化处理方法 |
CN106568436A (zh) * | 2016-10-26 | 2017-04-19 | 上海航天控制技术研究所 | 一种利用陀螺修正星敏感器轨道周期***误差的方法 |
CN106568436B (zh) * | 2016-10-26 | 2019-04-26 | 上海航天控制技术研究所 | 一种利用陀螺修正星敏感器轨道周期***误差的方法 |
CN106708072A (zh) * | 2016-11-30 | 2017-05-24 | 上海卫星工程研究所 | 一种天基望远镜高精度姿态确定与控制方法 |
CN107608065A (zh) * | 2017-10-18 | 2018-01-19 | 中国科学院光电技术研究所 | 一种运动平台下基于allan方差分析的望远镜自动指向修正方法 |
CN107608065B (zh) * | 2017-10-18 | 2020-03-20 | 中国科学院光电技术研究所 | 一种运动平台下基于allan方差分析的望远镜自动指向修正方法 |
CN108195403A (zh) * | 2017-12-28 | 2018-06-22 | 中国人民解放军国防科技大学 | 星敏感器在轨姿态测量数据综合误差模型的构建方法及装置 |
CN108195403B (zh) * | 2017-12-28 | 2020-05-22 | 中国人民解放军国防科技大学 | 星敏感器在轨姿态测量数据综合误差模型的构建方法及装置 |
CN112781616A (zh) * | 2020-12-24 | 2021-05-11 | 中国人民解放军国防科技大学 | 星敏感器在轨测量低频误差分析方法、装置和存储介质 |
CN112781616B (zh) * | 2020-12-24 | 2023-07-25 | 中国人民解放军国防科技大学 | 星敏感器在轨测量低频误差分析方法、装置和存储介质 |
CN112989637A (zh) * | 2021-05-06 | 2021-06-18 | 中国人民解放军国防科技大学 | 基于分布折算、近似和综合的星敏***可靠性评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101696885B (zh) | 2012-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101696885B (zh) | 一种提高星敏感器数据处理精度的方法 | |
CN107525503B (zh) | 基于双天线gps和mimu组合的自适应级联卡尔曼滤波方法 | |
CN103900565B (zh) | 一种基于差分gps的惯导***姿态获取方法 | |
CN102486377B (zh) | 一种光纤陀螺捷联惯导***初始航向的姿态获取方法 | |
Fakharian et al. | Adaptive Kalman filtering based navigation: An IMU/GPS integration approach | |
Wang et al. | Estimation of information sharing error by dynamic deformation between inertial navigation systems | |
CN104374388B (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN103076026B (zh) | 一种捷联惯导***中确定多普勒计程仪测速误差的方法 | |
CN105043415A (zh) | 基于四元数模型的惯性系自对准方法 | |
Xiaojuan et al. | A SINS/CNS deep integrated navigation method based on mathematical horizon reference | |
CN102680996A (zh) | 定位装置以及定位方法 | |
CN104236586B (zh) | 基于量测失准角的动基座传递对准方法 | |
CN102519485B (zh) | 一种引入陀螺信息的二位置捷联惯性导航***初始对准方法 | |
CN101900573B (zh) | 一种实现陆用惯性导航***运动对准的方法 | |
CN102620748A (zh) | 捷联惯导***晃动基座条件下杆臂效应的估计和补偿方法 | |
CN104457748A (zh) | 一种嵌入式瞄准吊舱测姿***及其传递对准方法 | |
CN102981151A (zh) | 相控阵雷达电控波束稳定方法 | |
CN103674064B (zh) | 捷联惯性导航***的初始标定方法 | |
CN103017787A (zh) | 适用于摇摆晃动基座的初始对准方法 | |
CN102116634A (zh) | 一种着陆深空天体探测器的降维自主导航方法 | |
CN106403952A (zh) | 一种动中通低成本组合姿态测量方法 | |
CN114777812B (zh) | 一种水下组合导航***行进间对准与姿态估计方法 | |
CN102393204B (zh) | 一种基于sins/cns的组合导航信息融合方法 | |
CN102768043A (zh) | 一种无外观测量的调制型捷联***组合姿态确定方法 | |
CN108592943A (zh) | 一种基于opreq方法的惯性系粗对准计算方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120523 Termination date: 20141105 |
|
EXPY | Termination of patent right or utility model |