CN106092092A - 面向脉冲星导航系统的分数阶可观测性分析方法 - Google Patents

面向脉冲星导航系统的分数阶可观测性分析方法 Download PDF

Info

Publication number
CN106092092A
CN106092092A CN201610387902.4A CN201610387902A CN106092092A CN 106092092 A CN106092092 A CN 106092092A CN 201610387902 A CN201610387902 A CN 201610387902A CN 106092092 A CN106092092 A CN 106092092A
Authority
CN
China
Prior art keywords
observability
navigation
analysis
fractional order
spacecraft
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
CN201610387902.4A
Other languages
English (en)
Other versions
CN106092092B (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.)
Wuhan University of Science and Engineering WUSE
Wuhan University of Science and Technology WHUST
Original Assignee
Wuhan University of Science and Engineering WUSE
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 Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN201610387902.4A priority Critical patent/CN106092092B/zh
Publication of CN106092092A publication Critical patent/CN106092092A/zh
Application granted granted Critical
Publication of CN106092092B publication Critical patent/CN106092092B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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。
CN201610387902.4A 2016-06-02 2016-06-02 面向脉冲星导航系统的分数阶可观测性分析方法 Expired - Fee Related CN106092092B (zh)

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 true CN106092092A (zh) 2016-11-09
CN106092092B 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)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106679675A (zh) * 2016-12-29 2017-05-17 北京理工大学 一种基于相对测量信息的火星最终接近段自主导航方法
CN113222095A (zh) * 2021-04-08 2021-08-06 浙江大学 基于进化计算的电力系统分数阶预测辅助状态估计方法
CN114577221A (zh) * 2022-01-27 2022-06-03 北京空间飞行器总体设计部 一种只用观测信息的自主导航系统可观测性解析判定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101762272A (zh) * 2010-01-18 2010-06-30 哈尔滨工业大学 一种基于可观测度分析的深空自主导航方法
US20130179128A1 (en) * 2012-01-05 2013-07-11 General Electric Company System And Method For Selecting A Representative Sensor Set Of A Power Plant
CN103591956A (zh) * 2013-11-27 2014-02-19 北京理工大学 一种基于可观测性分析的深空探测器自主导航方法
CN103616024A (zh) * 2013-11-27 2014-03-05 北京理工大学 一种行星探测进入段自主导航系统可观测度确定方法
CN103870714A (zh) * 2014-04-03 2014-06-18 武汉科技大学 一种基于高阶可观阵的导航系统可观测分析方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101762272A (zh) * 2010-01-18 2010-06-30 哈尔滨工业大学 一种基于可观测度分析的深空自主导航方法
US20130179128A1 (en) * 2012-01-05 2013-07-11 General Electric Company System And Method For Selecting A Representative Sensor Set Of A Power Plant
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)

* Cited by examiner, † Cited by third party
Title
刘劲等: "Fast non-linearly constrained least square joint estimation of position and velocity for X-ray pulsar-based navigation", 《IET RADAR, SONAR AND NAVIGATION》 *
张邦楚: "基于分数阶微积分的飞航式导弹控制系统设计方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *
王奕迪等: "基于单探测器的X射线脉冲星深空导航算法", 《力学学报》 *
郑广楼等: "单脉冲星自主导航系统可观测性分析", 《应用科学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106679675A (zh) * 2016-12-29 2017-05-17 北京理工大学 一种基于相对测量信息的火星最终接近段自主导航方法
CN106679675B (zh) * 2016-12-29 2019-10-29 北京理工大学 一种基于相对测量信息的火星最终接近段自主导航方法
CN113222095A (zh) * 2021-04-08 2021-08-06 浙江大学 基于进化计算的电力系统分数阶预测辅助状态估计方法
CN113222095B (zh) * 2021-04-08 2022-01-18 浙江大学 基于进化计算的电力系统分数阶预测辅助状态估计方法
CN114577221A (zh) * 2022-01-27 2022-06-03 北京空间飞行器总体设计部 一种只用观测信息的自主导航系统可观测性解析判定方法
CN114577221B (zh) * 2022-01-27 2024-03-26 北京空间飞行器总体设计部 一种只用观测信息的自主导航系统可观测性解析判定方法

Also Published As

Publication number Publication date
CN106092092B (zh) 2018-08-21

Similar Documents

Publication Publication Date Title
Nebeker Calculating the weather: Meteorology in the 20th century
Klee et al. Simulation of dynamic systems with MATLAB® and Simulink®
CN103299308B (zh) 利用实体建模计算任意形状密闭容器中的液位
Bennett Inverse methods in physical oceanography
Rodgers Estimable functions of age, period, and cohort effects
Salas et al. Linking expertise and naturalistic decision making
CN106092092B (zh) 面向脉冲星导航系统的分数阶可观测性分析方法
Li et al. An extended Kalman filter-based attitude tracking algorithm for star sensors
Mourtos et al. Open-ended problem-solving skills in thermal-fluids engineering
Possolo et al. Evaluating the uncertainty of input quantities in measurement models
CN108680187A (zh) 基于可见光源的x射线脉冲星导航地面验证系统
Mayorova et al. “Vivid mathematics” as a general vector of multidisciplinary STEM education for future aerospace engineers
Badiru et al. Handbook of measurements: benchmarks for systems accuracy and precision
CN105865432B (zh) 一种针对陀螺仪多源噪声的混合滤波方法与测试平台
CN106651727A (zh) 航班信息的处理及展示方法、系统
Simon A Student's Introduction to Engineering Design: Pergamon Unified Engineering Series
CN104408326B (zh) 一种对深空探测自主导航滤波算法的评估方法
Castro-Fernández et al. Three-dimensional unsteady aerodynamic analysis of a rigid-framed delta kite applied to airborne wind energy
WO2020213547A1 (ja) 気象予測データ作成プログラム、気象予測データ作成方法、及び、移動体
CN103870714B (zh) 一种基于高阶可观阵的导航系统可观测分析方法
Winegrad Celebrating the birth of modern computing: the fiftieth anniversary of a discovery at the Moore School of Engineering of the University of Pennsylvania
McClelland Towards Detecting Atmospheric Coherent Structures using Small Fixed-Wing Unmanned Aircraft
Gil Low Cost Physics Experiments using ITC
Hodge et al. Monte Carlo analysis of inaccuracies in estimated aircraft parameters caused by unmodeled flight instrumentation errors
Indicators Part A

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