CN106092092B - 面向脉冲星导航系统的分数阶可观测性分析方法 - Google Patents
面向脉冲星导航系统的分数阶可观测性分析方法 Download PDFInfo
- Publication number
- CN106092092B CN106092092B CN201610387902.4A CN201610387902A CN106092092B CN 106092092 B CN106092092 B CN 106092092B CN 201610387902 A CN201610387902 A CN 201610387902A CN 106092092 B CN106092092 B CN 106092092B
- Authority
- CN
- China
- Prior art keywords
- observability
- navigation
- spacecraft
- pulsar
- matrix
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 238000005259 measurement Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 230000007704 transition Effects 0.000 claims description 12
- 238000010276 construction Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 230000004069 differentiation Effects 0.000 claims description 5
- 238000002360 preparation method Methods 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 description 8
- 238000007792 addition Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- 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
-
- 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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明提供一种面向脉冲星导航系统的分数阶可观测性分析方法,包括预备阶段和系统可观测性分析阶段,所述预备阶段包括建立导航滤波所需的航天器的轨道动力学模型和脉冲星导航模型,所述系统可观测性分析阶段包括测量模型对时间求分数阶微分,获取分数阶测量矩阵,构造可观测性分析矩阵,对条件数重新排序,计算指数加权条件数,获得可观测性分析结果。根据本发明技术方案获得的可观测性分析结果更为精确,充分利用了历史测量信息,能够体现轨道要素对导航精度的影响,并且计算简单,实现方便,节约系统资源,能够高效获取分析结果。
Description
技术领域
本发明属于航天器自主导航领域,特别涉及一种面向脉冲星导航系统的分数阶可观测性分析方法。
背景技术
在航天器自主导航领域,系统可观测性分析是至关重要的。可观测性分析结果可反映航天器导航性能。这体现在以下方面:(1)若可观测性矩阵是满秩的,自主导航系统可正常工作;否则,系统无法正常工作。(2)可观测性矩阵的条件数能反映系统的好坏,即条件数较大时导航系统性能差,反之性能好。
可观测性分析方法可分为两个步骤:可观测性矩阵的构造,以及条件数的计算。已有学者提出了两种可观测性矩阵构造方法。在PWCS(piece-wise constant systems)中,分别对测量模型和状态转移模型求导,可获得测量矩阵和状态转移矩阵。利用这两个矩阵可构造可观测性矩阵。值得注意的是,测量矩阵和状态转移矩阵与当前的测量值和状态相关。因此,该可观测性模型与历史数据关系不大。在2014年,崔平远教授提出了基于李导数的可观测性矩阵构造方法。由于李导数注重局部信息,因此李导数方法不能反映整个导航系统的性能。总之,以上两个可观测性方法不能充分利用历史测量信息。传统的条件数用于反映导航系统性能。在2009年,房建成院士提出了混合条件数方法。
X射线脉冲星导航是一种新兴的航天器自主导航方法。在此系统中,测量模型可近似为线性模型,状态转移模型为非线性。在这种情况下,由PWCS和李导数方法构造的可观测性矩阵相等,并且是常值。他们不能反映轨道对导航性能的影响。即对于不同轨道,可观测矩阵是相同的。因此,对于不同的轨道,传统条件数和混合条件数都无法给出不同的分析结果。从以上的分析可以看出,可观测性矩阵的构造至关重要。
发明内容
本发明提出了一种面向脉冲星导航系统的分数阶可观测性分析方法,旨在为脉冲星导航系统提供精确的评估。
本发明技术方案提供一种面向脉冲星导航系统的分数阶可观测性分析方法,包括预备阶段和系统可观测性分析阶段,
所述预备阶段,包括建立导航滤波所需的各种模型,包括以下步骤,
步骤A1,建立航天器的轨道动力学模型,实现如下,
设航天器的状态矢量X为,
其中,r=[x,y,z]T和v=[vx,vy,vz]T分别为航天器的位置和速度矢量,x,y,z分别为航天器的位置在三轴上的分量,vx,vy,vz分别为航天器的速度在三轴上的分量;
则航天器的轨道动力学模型为,
其中,分别为r,v的导数,|r|表示矢量r的模,ΔF表示其他摄动力造成的加速度;
表示为,
其中,是状态矢量X的导数,为时刻t的f(X,t)为航天器的状态转移模型,μ是火星的引力常数,ω(t)为时刻t航天器的导航系统噪声;
步骤A2,建立脉冲星导航模型,
其中,tj和分别是脉冲到达航天器和太阳系质心的时间,nj是第j颗脉冲星方位矢量,D0是脉冲星到太阳系质心的距离,b是太阳系质心相对于太阳质心的位置矢量,c是光速,σj是TOA测量噪声,|·|表示矢量的模,j表示脉冲星编号,j=1,2,…J,J表示导航脉冲星数量;
所述系统可观测性分析阶段包括构造可观测性分析矩阵和计算条件数,获得可观测性分析结果,包括以下步骤,
步骤B1,分数阶可观测性矩阵构造包括以下步骤,
步骤B11,测量模型对时间t求分数阶微分,
其中,v为分数阶微分的阶数,0<v<1;T为脉冲星导航滤波周期,Z为脉冲星导航测量值,n表示导航脉冲星方位矢量矩阵;r(t)为t时刻航天器的位置矢量,Γ(.)为Gamma函数,表示Z的v阶Grunwald-Letnikov分数阶导数;m!表示截止数变量m的阶乘;
根据脉冲到达航天器和太阳系质心的时间tj和定义如下,
步骤B12,对状态矢量X求导,得到分数阶测量矩阵如下式,
其中,03×3为0矩阵,I6×6为单位矩阵,F是状态转移模型,
其中,矩阵S的表达式为,
步骤B13,构造分数阶可观测性矩阵O,
步骤B2,指数加权条件数计算包括以下步骤,
步骤B21,计算分数阶可观测性矩阵的条件数,对条件数重新排序,得到序列c(i),i=1,2,…N,c(1)≤c(2)≤c(3)≤…≤c(N),N为条件数的总数;
步骤B22,计算指数加权条件数C,指数加权条件数C越大表明导航性能越好,反之导航性能越差,
其中,a和b是预设的系数,e是数学常数。
而且,a和b的取值分别为20和150。
本发明充分利用了分数阶微分的记忆性,构造了分数阶可观测性矩阵,其分析结果可体现轨道要素对脉冲星导航的影响。因此,本发明对航天器自主导航具有重要的实际意义。
本发明与现有技术相比的优点在于:
(1)可观测性分析结果更为精确。条件数是用于判断导航性能好坏的一个重要指标。传统PWCS和李导数方法都不能体现轨道要素对脉冲星导航性能的影响。而本发明根据分数阶微分的记忆性,从全局出发,充分利用了历史测量信息,构造分数阶可观测性矩阵。分数阶可观测性分析矩阵能体现轨道要素对导航精度的影响,且分析结果与仿真结果一致。
(2)分数阶可观测性分析矩阵及指数加权条件数计算量小。基于李导数的方法涉及多次进行复杂的李导数运算,计算量很大。本发明技术方案仅涉及几次矩阵乘法和加法,计算简单,实现方便,节约系统资源,能够高效获取分析结果。
附图说明
图1为本发明实施例的分数阶可观测性分析流程图。
具体实施方式
本发明技术方案可采用计算机软件方式支持自动运行流程。以下结合附图和实施例详细说明本发明技术方案。
本发明首先,构造分数阶测量矩阵,这样可以充分利用历史测量信息。然后,结合状态转移矩阵,构造分数阶可观测性矩阵。最后,提出指数加权条件数,将一个周期内变化的条件数合成为一个值,利用该值来评价导航系统精度。
本发明实施例采用火星探测器,导航脉冲星为PSR B0531+21,B1821-24,和B1937+21。涉及的火星卫星轨道,如表1所示。
表1金星快车初始轨道参数
本发明提供一种面向脉冲星导航系统的分数阶可观测性分析方法,包括预备阶段和系统可观测性分析阶段。
所述预备阶段,包括建立导航滤波所需的各种模型,包括以下步骤,
步骤A1,建立航天器的轨道动力学模型,实现如下,
设航天器的状态矢量X为,
其中,r=[x,y,z]T和v=[vx,vy,vz]T分别为航天器的位置和速度矢量,x,y,z分别为航天器的位置在三轴上的分量,vx,vy,vz分别为航天器的速度在三轴上的分量;
则航天器的轨道动力学模型为,
其中,分别为r,v的导数,|r|表示矢量r的模,ΔF表示其他行星、太阳光压等其他摄动力造成的加速度。
式(二)表示为,
其中,是状态矢量X的导数,为时刻t的f(X,t)为航天器的状态转移模型,μ是火星的引力常数,ω(t)为时刻t航天器的导航系统噪声;
步骤A2,建立脉冲星导航模型;
其中,tj和分别是脉冲到达航天器和太阳系质心的时间;nj是第j颗脉冲星方位矢量;D0是脉冲星到太阳系质心的距离,b是太阳系质心相对于太阳质心的位置矢量,c是光速,σj是TOA(到达时间,time-of-arrival)测量噪声,|·|表示矢量的模。j表示脉冲星编号,j=1,2,…J,J表示导航脉冲星数量。
所述系统可观测性分析阶段包括构造可观测性分析矩阵和计算条件数,获得可观测性分析结果。
步骤B1,分数阶可观测性矩阵构造包括以下步骤:
步骤B11,测量模型对时间t求分数阶微分。
其中,v为分数阶微分的阶数,0<v<1。本发明实施例中,分数阶阶数取0.2。T为脉冲星导航滤波周期,本发明实施例中,取值为300s。Z为脉冲星导航测量值,n=[n1,n2,…nJ]T表示导航脉冲星方位矢量矩阵。r(t)为t时刻航天器的位置矢量,Γ(.)为Gamma函数,表示Z的v阶Grunwald-Letnikov分数阶导数。m!表示截止数变量m的阶乘。
根据脉冲到达航天器和太阳系质心的时间tj和定义如下,
步骤B12,上式对状态矢量X求导,得到分数阶测量矩阵如下式。
其中,03×3为0矩阵,I6×6为单位矩阵,F是状态转移模型,即状态转移矩阵:
其中,矩阵S的表达式为:
步骤B13,利用式(七,八)构造分数阶可观测性矩阵O。
步骤B2,指数加权条件数计算包括以下步骤:
步骤B21,计算分数阶可观测性矩阵的条件数,具体计算为常见矩阵运算方式,本发明不予赘述;对条件数重新排序,得到序列c(i),i=1,2,…N,c(1)≤c(2)≤c(3)≤…≤c(N),N为条件数的总数。
步骤B22,计算指数加权条件数C。
其中,a和b是预设的系数,e是数学常数。具体实施时,本领域技术人员可自行预设a和b的取值,可采用经验值。本发明实施例中采用优选的取值,分别为20和150。
本发明的可观测性分析结果为指数加权条件数,该值与导航性能息息相关,具体表现为:指数加权条件数越大,表明导航性能越好;反之,导航性能越差。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (2)
1.一种面向脉冲星导航系统的分数阶可观测性分析方法,其特征在于:包括预备阶段和系统可观测性分析阶段,
所述预备阶段,包括建立导航滤波所需的各种模型,包括以下步骤,
步骤A1,建立航天器的轨道动力学模型,实现如下,
设航天器的状态矢量X为,
其中,r=[x,y,z]T和v=[vx,vy,vz]T分别为航天器的位置和速度矢量,x,y,z分别为航天器的位置在三轴上的分量,vx,vy,vz分别为航天器的速度在三轴上的分量;
则航天器的轨道动力学模型为,
其中,分别为r,v的导数,|r|表示矢量r的模,ΔF表示其他摄动力造成的加速度;
表示为,
其中,是状态矢量X的导数,为时刻t的f(X,t)为航天器的状态转移模型,μ是火星的引力常数,ω(t)为时刻t航天器的导航系统噪声;
步骤A2,建立脉冲星导航模型,
其中,tj和分别是脉冲到达航天器和太阳系质心的时间,nj是第j颗脉冲星方位矢量,D0是脉冲星到太阳系质心的距离,b是太阳系质心相对于太阳质心的位置矢量,c是光速,σj是TOA测量噪声,|·|表示矢量的模,j表示脉冲星编号,j=1,2,…J,J表示导航脉冲星数量;
所述系统可观测性分析阶段包括构造可观测性分析矩阵和计算条件数,获得可观测性分析结果,包括以下步骤,
步骤B1,分数阶可观测性矩阵构造包括以下步骤,
步骤B11,测量模型对时间t求分数阶微分,
其中,v为分数阶微分的阶数,0<v<1;T为脉冲星导航滤波周期,Z为脉冲星导航测量值,n表示导航脉冲星方位矢量矩阵;r(t)为t时刻航天器的位置矢量,Γ(.)为Gamma函数,表示Z的v阶Grunwald-Letnikov分数阶导数;m!表示截止数变量m的阶乘;
根据脉冲到达航天器和太阳系质心的时间tj和定义如下,
步骤B12,对状态矢量X求导,得到分数阶测量矩阵如下式,
其中,03×3为0矩阵,I6×6为单位矩阵,F是状态转移模型,
其中,矩阵S的表达式为,
步骤B13,构造分数阶可观测性矩阵O,
步骤B2,指数加权条件数计算包括以下步骤,
步骤B21,计算分数阶可观测性矩阵的条件数,对条件数重新排序,得到序列c(i),i=1,2,…N,c(1)≤c(2)≤c(3)≤…≤c(N),N为条件数的总数;
步骤B22,计算指数加权条件数C,指数加权条件数C越大表明导航性能越好,反之导航性能越差,
其中,a和b是预设的系数,e是数学常数。
2.根据权利要求1所述面向脉冲星导航系统的分数阶可观测性分析方法,其特征在于:a和b的取值分别为20和150。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610387902.4A CN106092092B (zh) | 2016-06-02 | 2016-06-02 | 面向脉冲星导航系统的分数阶可观测性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610387902.4A CN106092092B (zh) | 2016-06-02 | 2016-06-02 | 面向脉冲星导航系统的分数阶可观测性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106092092A CN106092092A (zh) | 2016-11-09 |
CN106092092B true CN106092092B (zh) | 2018-08-21 |
Family
ID=57448155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610387902.4A Expired - Fee Related CN106092092B (zh) | 2016-06-02 | 2016-06-02 | 面向脉冲星导航系统的分数阶可观测性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106092092B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106679675B (zh) * | 2016-12-29 | 2019-10-29 | 北京理工大学 | 一种基于相对测量信息的火星最终接近段自主导航方法 |
CN113222095B (zh) * | 2021-04-08 | 2022-01-18 | 浙江大学 | 基于进化计算的电力系统分数阶预测辅助状态估计方法 |
CN114577221B (zh) * | 2022-01-27 | 2024-03-26 | 北京空间飞行器总体设计部 | 一种只用观测信息的自主导航系统可观测性解析判定方法 |
CN114623832B (zh) * | 2022-01-28 | 2024-07-23 | 北京空间飞行器总体设计部 | 自主导航系统可观测能力降维表征及解析判定方法和系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101762272A (zh) * | 2010-01-18 | 2010-06-30 | 哈尔滨工业大学 | 一种基于可观测度分析的深空自主导航方法 |
CN103591956A (zh) * | 2013-11-27 | 2014-02-19 | 北京理工大学 | 一种基于可观测性分析的深空探测器自主导航方法 |
CN103616024A (zh) * | 2013-11-27 | 2014-03-05 | 北京理工大学 | 一种行星探测进入段自主导航系统可观测度确定方法 |
CN103870714A (zh) * | 2014-04-03 | 2014-06-18 | 武汉科技大学 | 一种基于高阶可观阵的导航系统可观测分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130179128A1 (en) * | 2012-01-05 | 2013-07-11 | General Electric Company | System And Method For Selecting A Representative Sensor Set Of A Power Plant |
-
2016
- 2016-06-02 CN CN201610387902.4A patent/CN106092092B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101762272A (zh) * | 2010-01-18 | 2010-06-30 | 哈尔滨工业大学 | 一种基于可观测度分析的深空自主导航方法 |
CN103591956A (zh) * | 2013-11-27 | 2014-02-19 | 北京理工大学 | 一种基于可观测性分析的深空探测器自主导航方法 |
CN103616024A (zh) * | 2013-11-27 | 2014-03-05 | 北京理工大学 | 一种行星探测进入段自主导航系统可观测度确定方法 |
CN103870714A (zh) * | 2014-04-03 | 2014-06-18 | 武汉科技大学 | 一种基于高阶可观阵的导航系统可观测分析方法 |
Non-Patent Citations (4)
Title |
---|
Fast non-linearly constrained least square joint estimation of position and velocity for X-ray pulsar-based navigation;刘劲等;《IET Radar, Sonar and Navigation》;20141231;第8卷(第9期);第1154-1163页 * |
单脉冲星自主导航系统可观测性分析;郑广楼等;《应用科学学报》;20080930;第26卷(第5期);第506-510页 * |
基于分数阶微积分的飞航式导弹控制系统设计方法研究;张邦楚;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20070115;第C032-33页 * |
基于单探测器的X射线脉冲星深空导航算法;王奕迪等;《力学学报》;20120930;第44卷(第5期);第912-918页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106092092A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106092092B (zh) | 面向脉冲星导航系统的分数阶可观测性分析方法 | |
Irrgang et al. | Milky Way mass models for orbit calculations | |
Pancoast et al. | The lick AGN monitoring project 2011: dynamical modeling of the broad-line region in Mrk 50 | |
Bovy et al. | A direct dynamical measurement of the Milky Way's disk surface density profile, disk scale length, and dark matter profile at 4 kpc≲ R≲ 9 kpc | |
CN104236546B (zh) | 一种卫星星光折射导航误差确定与补偿方法 | |
Reese et al. | Regular patterns in the acoustic spectrum of rapidly rotating stars | |
EP3293484B1 (en) | Astronomical velocity measurement and ground radio based mars capture phase navigation method | |
CN109255096A (zh) | 一种基于微分代数的地球同步卫星轨道不确定演化方法 | |
CN103900577B (zh) | 一种面向编队飞行的相对导航测速及组合导航方法 | |
Takahashi et al. | Morphology driven density distribution estimation for small bodies | |
CN111191368B (zh) | 一种连续小推力行星际转移轨道优化方法和装置 | |
CN103218482B (zh) | 一种动力学系统中不确定参数的估计方法 | |
CN102735260A (zh) | 一种星敏感器在轨测量误差的确定方法 | |
CN108917772A (zh) | 基于序列图像的非合作目标相对导航运动估计方法 | |
RU2432592C1 (ru) | Моделирующий комплекс для проверки системы управления беспилотного летательного аппарата | |
CN117933095B (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
Lima et al. | A 3% determination of H0 at intermediate redshifts | |
CN110006462B (zh) | 基于奇异值分解的星敏感器在轨标定方法 | |
CN103296995B (zh) | 任意维高阶(≥4阶)无味变换与无味卡尔曼滤波方法 | |
Terrée et al. | Addressing the gas kinetics Boltzmann equation with branching-path statistics | |
CN113204909B (zh) | 基于地基观测光度信号的卫星几何特征与姿态估计方法 | |
Pecora et al. | Multipoint turbulence analysis with helioswarm | |
Luo et al. | Emergent Gravity Fails to Explain Color-dependent Galaxy–Galaxy Lensing Signal from SDSS DR7 | |
CN111428912A (zh) | 一种基于支持向量机的火星探测器轨道预测方法及系统 | |
CN103870714B (zh) | 一种基于高阶可观阵的导航系统可观测分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180821 Termination date: 20190602 |