CN107085225B - 一种基于gnss信息的冷启动航天器在轨导航初值建立方法 - Google Patents
一种基于gnss信息的冷启动航天器在轨导航初值建立方法 Download PDFInfo
- Publication number
- CN107085225B CN107085225B CN201710501331.7A CN201710501331A CN107085225B CN 107085225 B CN107085225 B CN 107085225B CN 201710501331 A CN201710501331 A CN 201710501331A CN 107085225 B CN107085225 B CN 107085225B
- Authority
- CN
- China
- Prior art keywords
- gnss
- spacecraft
- initial value
- coordinate system
- orbit
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开一种基于GNSS信息的冷启动航天器在轨导航初值建立方法,所述方法包括:S1:获取连续多帧的航天器GNSS输出的时间、位置、速度;S2:判断时间有效性;S3:判断数据有效性;S4:轨道外推,建立在轨导航初值,本发明能够使航天器在冷启动状态下利用GNSS信息获取在轨导航初值,可以有效降低由于GNSS数据的不准确带来的风险,且本发明可以在准确获取在轨导航初值的前提下有效减少航天器与运载的接口,降低成本。
Description
技术领域
本发明涉及飞行器控制领域。更具体地,涉及一种基于GNSS信息的冷启动航天器在轨导航初值建立方法。
背景技术
航天器的在轨导航初值建立是航天器实施在轨导航的必要条件。传统的航天器在轨导航初值建立主要有两种方式:一是在航天器起飞之前,利用发射惯性系初始位置作为导航初值,并利用精确的发射场地理位置信息以及对时信息将发射惯性系航天器位置速度转换到赤道惯性系位置速度,在起飞后的整个运载段航天器需要全程加电并持续导航,在星箭分离时刻记录航天器位置和速度,将其作为航天器在轨导航的轨道初值;另一种是通过基准传递的方法,即在航天器星箭分离前,通过运载主惯导系统向航天器子惯导系统传递轨道初值,航天器子惯导系统接受到轨道初值后,启动在轨导航,在星箭分离时刻记录航天器位置和速度,将其作为航天器在轨导航的轨道初值。以上的方法中,第一种方法需要航天器起飞前处于加电状态,并保持运载段持续加电导航;第二种方法需要航天器与运载进行数据交互,需要一些接口和电缆支持,增加了整个系统的重量和成本,降低了系统的可靠性。对于与运载只有机械接口的航天器,在与运载分离后才加电,加电时航天器无任何先验信息,属于冷启动范畴。在这样的情况下,现有技术方案无法实现轨道初值建立。可以通过利用航天器GNSS设备提供的信息进行轨道初值建立,然而,GNSS接收设备容易受到各种短时间电磁因素、环境因素干扰,影响定时、定位功能;同时航天器是在嵌入式环境下实时进行导航初值建立,缺乏外部其他比较基准,无法进行必要的验证,其定时、定位误差可能带入整个导航系统中,导致导航错误,进而导致整个航天器的飞行任务失败。
因此,需要提供一种抗干扰能力强、误差小的基于GNSS信息的冷启动航天器在轨导航初值建立方法。
发明内容
本发明的一个目的在于提供一种基于GNSS信息的冷启动航天器在轨导航初值建立方法,提高基于GNSS信息的冷启动航天器在轨导航初值建立的抗干扰能力,克服干扰因素,并进行验证,减小误差,提高导航准确度。
为达到上述目的,本发明采用下述技术方案:
本发明公开了一种基于GNSS信息的冷启动航天器在轨导航初值建立方法,所述方法包括:
S1:获取航天器GNSS输出数据;
S2:判断时间有效性;
S3:判断数据有效性;
S4:轨道外推,建立在轨导航初值。
优选地,所述S1获取连续多帧的航天器GNSS输出的时间、位置、速度数据。
优选地,S2包括:
tB=Week×7×24×3600+WeekSec
令
其中,Weekj、WeekSecj分别为GNSS输出的连续多帧的周、周秒,j为帧数,取值范围为1~N,N为采集的帧数;
若满足则通过时间有效性判断,否则,舍弃该帧数据,不再进行后续的数据有效性判断。
优选地,所述S3包括:
S31:将WGS84坐标系的航天器GNSS输出的位置、速度转换到J2000坐标系下的位置、速度为
其中,Ciw为从WGS84坐标系到J2000坐标系的转换矩阵,ωD为地球转速,为WGS84坐标系的GNSS输出的位置坐标,为WGS84坐标系的GNSS输出的速度坐标,上角标j表示第j帧数,取值范围为1~N,,N为采集的帧数,下角标w表示为WGS84坐标系下的数值,为J2000坐标系的GNSS输出的位置坐标,为J2000坐标系的GNSS输出的速度坐标;
S32:计算轨道根数,判断轨道根数的有效性,若满足
(aj-2≥Re+100000)&&(aj-1≥Re+100000)&&(aj≥Re+100000)
(|aj-1-aj-2|≤Δa)&&(|aj-aj-1|≤Δa)
(|ej-1-ej-2|≤Δe)&&(|ej-ej-1|≤Δe)
(|Ωj-1-Ωj-2|≤ΔΩ)&&(|Ωj-Ωj-1|≤ΔΩ)
(|ij-1-ij-2|≤Δi)&&(|ij-ij-1|≤Δi)
(|ωj-1-ωj-2|≤Δω)&&(|ωj-ωj-1|≤Δω)
其中,Re为地球半径,aj,ej,Ωj,ij,ωj为轨道根数,Δa、Δe、ΔΩ、Δi、Δω分别为轨道根数稳定性判断阈值,根据实际情况进行设定;
则数据有效性判断通过,否则舍弃该帧数据。
优选地,假设第K帧时通过数据有效性判断,则航天器的在轨导航初值为
tB=Weekk×7×24×3600+WeekSeck
其中,Weekk、WeekSeck分别为GNSS输出的K帧时的周、周秒,ΔtBG为当前时刻与GNSS数据天秒时刻的差值,ΔT为采样时间间隔,为J2000坐标系的GNSS输出的K帧时的位置坐标,为J2000坐标系的GNSS输出的K帧时的速度坐标,x0,y0,z0航天器的在轨导航初值的位置坐标,vx0,vy0,vz0为航天器的在轨导航初值的速度坐标。
本发明的有益效果如下:
本发明所述技术方案能够使航天器在冷启动状态下利用GNSS信息获取在轨导航初值,无需地面和运载段加电,不依赖运载基准传递数据,可以有效降低由于GNSS数据的不准确带来的风险,减少误差,且本发明的计算方法抗干扰性强,定位准确。本发明可以在准确获取在轨导航初值的前提下有效减少航天器与运载的接口,降低成本。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出本发明一种基于GNSS信息的冷启动航天器在轨导航初值建立方法具体实施例的流程图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
如图1所示,本发明公开了一种基于GNSS信息的冷启动航天器在轨导航初值建立方法,该方法包括:
S1:获取航天器GNSS输出数据。获取连续多帧的航天器GNSS输出的时间、位置、速度。采集GNSS输出的连续N帧的周、周秒、位置、速度,将其分别记录为Weekj,WeekSecj,(j=1:N)。其中,上角标j表示第j帧数,下角标w表示为WGS84坐标系下的数值。
S2:判断时间有效性。从第二帧开始,进行时间有效性判断。具体的,
tB=Week×7×24×3600+WeekSec
令
其中,Weekj、WeekSecj分别为GNSS输出的连续多帧的周、周秒,tB为其转成成的天秒,j为帧数,取值范围为1~N,N为采集的帧数。Week和WeekSec为航天器经过对时后根据晶振计算得到的精确的星上周和周秒。
若满足则通过时间有效性判断,否则,舍弃该帧数据,不再进行后续的数据有效性判断。
S3:判断数据有效性。所述S3包括:
S31:将WGS84坐标系的航天器GNSS输出的位置、速度转换到J2000坐标系下的位置、速度为
其中,Ciw为从WGS84坐标系到J2000坐标系的转换矩阵,ωD为地球转速,为WGS84坐标系的GNSS输出的位置坐标,为WGS84坐标系的GNSS输出的速度坐标,上角标j表示第j帧数,取值范围为1~N,,N为采集的帧数,下角标w表示为WGS84坐标系下的数值,为J2000坐标系的GNSS输出的位置坐标,为J2000坐标系的GNSS输出的速度坐标;
S32:计算轨道根数,判断轨道根数的有效性。轨道根数aj,ej,Ωj,ij,ωj的计算方法为成熟算法,不再赘述。若满足
(aj-2≥Re+100000)&&(aj-1≥Re+100000)&&(aj≥Re+100000)
(|aj-1-aj-2|≤Δa)&&(|aj-aj-1|≤Δa)
(|ej-1-ej-2|≤Δe)&&(|ej-ej-1|≤Δe)
(|Ωj-1-Ωj-2|≤ΔΩ)&&(|Ωj-Ωj-1|≤ΔΩ)
(|ij-1-ij-2|≤Δi)&&(|ij-ij-1|≤Δi)
(|ωj-1-ωj-2|≤Δω)&&(|ωj-ωj-1|≤Δω)
其中,Re为地球半径,aj,ej,Ωj,ij,ωj为轨道根数,Δa、Δe、ΔΩ、Δi、Δω分别为轨道根数稳定性判断阈值,根据实际情况进行设定;
则数据有效性判断通过,否则舍弃该帧数据。
S4:轨道外推,建立在轨导航初值。假设第K帧时通过数据有效性判断,则航天器的在轨导航初值为
tB=Weekk×7×24×3600+WeekSeck
其中,Weekk、WeekSeck分别为GNSS输出的K帧时的周、周秒,ΔtBG为当前时刻与GNSS数据天秒时刻的差值,ΔT为采样时间间隔,为J2000坐标系的GNSS输出的K帧时的位置坐标,为J2000坐标系的GNSS输出的K帧时的速度坐标,x0,y0,z0航天器的在轨导航初值的位置坐标,vx0,vy0,vz0为航天器的在轨导航初值的速度坐标。
下面通过一个实例来对本发明作进一步地说明,续采集5帧GNSS数据,如表1所示。
表1 GNSS原始输出
表1的GNSS原始输出数据经过时间有效性判断后的结果如表2所示。
表2 时间有效性计算结果
注:-表示无需计算。
表1的GNSS输出的WGS84坐标系的数据经过转换后得到的J2000系的结果为表3所示。
表3 J2000坐标系下的结果
表3的GNSS原始数据经过数据有效性判断后,可以得到轨道六根数,下面仅以半长轴a的结果作为示例,对于其他轨道根数的判断类似。结果为表4所示。
表4 数据有效性判断结果
由表4可知,数据有效性性判断通过。
经过外推,得到
x0=-1788880.7829806;
y0=-5472827.5559055;
z0=-3561727.5989835;
vx0=5686.1912620268;
vy0=-3977.2912748817;
vz0=3273.3158471800。
本发明利用GNSS接收机在定时解算过程是个具有固定周期、每个周期之间的数据具有连续性的特点,通过GNSS输出绝对时间连续多帧数据间隔关系进行逻辑判断,确保GNSS时间的有效性;同时,利用航天器轨道根数中的aj,ej,Ωj,ij,ωj在短时间内基本保持不变的特性,将GNSS输出的WGS84坐标系下的位置、速度转换为J2000系下的位置速度,并计算其轨道根数,根据轨道根数的稳定性,确保GNSS输出位置速度的有效性。在时间和位置速度均有效的前提下,对位置速度进行轨道外推,建立在轨导航初值。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (4)
1.一种基于GNSS信息的冷启动航天器在轨导航初值建立方法,其特征在于,所述方法包括:
S1:获取航天器GNSS输出数据;
S2:判断时间有效性,包括:
tB=Week×7×24×3600+WeekSec
令
其中,Weekj、WeekSecj分别为GNSS输出的连续多帧的周、周秒,j为帧数,取值范围为1~N,N为采集的帧数;
若满足则通过时间有效性判断,否则,舍弃该帧数据,不再进行后续的数据有效性判断;
S3:判断数据有效性;
S4:轨道外推,建立在轨导航初值。
2.根据权利要求1所述的冷启动航天器在轨导航初值建立方法,其特征在于,所述S1获取连续多帧的航天器GNSS输出的时间、位置、速度数据。
3.根据权利要求1所述的冷启动航天器在轨导航初值建立方法,其特征在于,所述S3包括:
S31:将WGS84坐标系的航天器GNSS输出的位置、速度转换到J2000坐标系下的位置、速度为
其中,Ciw为从WGS84坐标系到J2000坐标系的转换矩阵,ωD为地球转速,为WGS84坐标系的GNSS输出的位置坐标,为WGS84坐标系的GNSS输出的速度坐标,上角标j表示第j帧数,取值范围为1~N, N为采集的帧数,下角标w表示为WGS84坐标系下的数值,为J2000坐标系的GNSS输出的位置坐标,为J2000坐标系的GNSS输出的速度坐标;
S32:计算轨道根数,判断轨道根数的有效性,若满足
(aj-2≥Re+100000)&&(aj-1≥Re+100000)&&(aj≥Re+100000)
(|aj-1-aj-2|≤△a)&&(|aj-aj-1|≤△a)
(|ej-1-ej-2|≤△e)&&(|ej-ej-1|≤△e)
(|Ωj-1-Ωj-2|≤△Ω)&&(|Ωj-Ωj-1|≤△Ω)
(|ij-1-ij-2|≤△i)&&(|ij-ij-1|≤△i)
(|ωj-1-ωj-2|≤△ω)&&(|ωj-ωj-1|≤△ω)
其中,Re为地球半径,aj,ej,Ωj,ij,ωj为轨道根数,△a、△e、△Ω、△i、△ω分别为轨道根数稳定性判断阈值,根据实际情况进行设定;
则数据有效性判断通过,否则舍弃该帧数据。
4.根据权利要求1所述的冷启动航天器在轨导航初值建立方法,其特征在于,假设第K帧时通过数据有效性判断,则航天器的在轨导航初值为
tB=Weekk×7×24×3600+WeekSeck
其中,Weekk、WeekSeck分别为GNSS输出的K帧时的周、周秒,△tBG为当前时刻与GNSS数据天秒时刻的差值,△T为采样时间间隔,为J2000坐标系的GNSS输出的K帧时的位置坐标,为J2000坐标系的GNSS输出的K帧时的速度坐标,x0,y0,z0航天器的在轨导航初值的位置坐标,vx0,vy0,vz0为航天器的在轨导航初值的速度坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710501331.7A CN107085225B (zh) | 2017-06-27 | 2017-06-27 | 一种基于gnss信息的冷启动航天器在轨导航初值建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710501331.7A CN107085225B (zh) | 2017-06-27 | 2017-06-27 | 一种基于gnss信息的冷启动航天器在轨导航初值建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107085225A CN107085225A (zh) | 2017-08-22 |
CN107085225B true CN107085225B (zh) | 2019-10-25 |
Family
ID=59606675
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710501331.7A Active CN107085225B (zh) | 2017-06-27 | 2017-06-27 | 一种基于gnss信息的冷启动航天器在轨导航初值建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107085225B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110011719B (zh) * | 2019-03-05 | 2020-05-29 | 上海卫星工程研究所 | 一种利用gps定轨数据的星上轨道获取控制方法 |
CN113188571B (zh) * | 2021-04-30 | 2024-04-30 | 无锡卡尔曼导航技术有限公司 | 一种用于组合导航动基座初始对准的载体前进/倒车判定方法 |
CN113536547B (zh) * | 2021-06-25 | 2023-07-14 | 北京控制工程研究所 | 一种可靠的多源轨道外推自主选用方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140070722A (ko) * | 2012-11-26 | 2014-06-11 | 한국전자통신연구원 | 다중속도시스템 결합 장치 및 운용 방법 |
CN104777490A (zh) * | 2015-03-18 | 2015-07-15 | 广东工业大学 | 一种导航卫星信号接收机及其冷启动方法 |
CN105929410A (zh) * | 2016-04-14 | 2016-09-07 | 上海交通大学 | 一种gnss接收机状态调度方法和gnss接收机 |
CN105954776A (zh) * | 2016-05-14 | 2016-09-21 | 四川中卫北斗科技有限公司 | 一种导航信号接收方法和接收机 |
CN106199661A (zh) * | 2016-06-29 | 2016-12-07 | 西北工业大学 | 基于阵列天线确定微小卫星相对位置和相对姿态的方法 |
CN106680843A (zh) * | 2016-12-12 | 2017-05-17 | 大唐半导体设计有限公司 | 一种接收机及实现接收机启动的方法 |
-
2017
- 2017-06-27 CN CN201710501331.7A patent/CN107085225B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140070722A (ko) * | 2012-11-26 | 2014-06-11 | 한국전자통신연구원 | 다중속도시스템 결합 장치 및 운용 방법 |
CN104777490A (zh) * | 2015-03-18 | 2015-07-15 | 广东工业大学 | 一种导航卫星信号接收机及其冷启动方法 |
CN105929410A (zh) * | 2016-04-14 | 2016-09-07 | 上海交通大学 | 一种gnss接收机状态调度方法和gnss接收机 |
CN105954776A (zh) * | 2016-05-14 | 2016-09-21 | 四川中卫北斗科技有限公司 | 一种导航信号接收方法和接收机 |
CN106199661A (zh) * | 2016-06-29 | 2016-12-07 | 西北工业大学 | 基于阵列天线确定微小卫星相对位置和相对姿态的方法 |
CN106680843A (zh) * | 2016-12-12 | 2017-05-17 | 大唐半导体设计有限公司 | 一种接收机及实现接收机启动的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107085225A (zh) | 2017-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111045068B (zh) | 一种基于非导航卫星信号的低轨卫星自主定轨定姿方法 | |
CN110487301A (zh) | 一种雷达辅助机载捷联惯性导航系统初始对准方法 | |
CN108871336B (zh) | 一种车辆位置估算系统及方法 | |
CN107085225B (zh) | 一种基于gnss信息的冷启动航天器在轨导航初值建立方法 | |
CN101916915B (zh) | 动中通卫星通信系统天线波束跟踪装置及跟踪方法 | |
CN102679985B (zh) | 一种应用星间跟踪的航天器星座分散化自主导航方法 | |
CN106643709B (zh) | 一种海上运载体的组合导航方法及装置 | |
CN101733746A (zh) | 空间机器人非合作目标自主识别与捕获方法 | |
CN109823571A (zh) | 一种遥感微纳卫星的多阶段姿态控制方法 | |
CN108519110B (zh) | 基于图像的空间非合作目标自主导航在轨验证系统 | |
CN103017772B (zh) | 一种基于可观性分析的光学和脉冲星融合自主导航方法 | |
CN101344391A (zh) | 基于全功能太阳罗盘的月球车位姿自主确定方法 | |
CN105160125A (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN103592533A (zh) | 一种基于小卫星信息系统的数传天线整星测试方法 | |
CN103674021A (zh) | 基于捷联惯导与星敏感器的组合导航系统及方法 | |
CN107449444A (zh) | 一种多星图姿态关联的星敏感器内参数标定方法及其装置 | |
CN108512590A (zh) | 一种卫星姿轨控分系统与gnss子系统的联试系统及方法 | |
CN106410410A (zh) | 一种具有物理水平平台的vsat天线系统卫星捕获跟踪方法 | |
CN104048664A (zh) | 一种导航卫星星座自主定轨的方法 | |
CN108931791A (zh) | 卫惯紧组合钟差修正系统和方法 | |
CN102944238B (zh) | 一种行星探测器接近目标过程中相对位置确定方法 | |
CN112926208A (zh) | 一种低轨卫星被动对月定标时机计算方法和系统 | |
CN109269504A (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
CN102914306A (zh) | 一种双探头星敏感器及其设计方法 | |
Park et al. | Effects of initial attitude estimation errors on loosely coupled smartphone GPS/IMU integration system |
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 |