CN109358352A - 一种星载实时掩星预报方法 - Google Patents
一种星载实时掩星预报方法 Download PDFInfo
- Publication number
- CN109358352A CN109358352A CN201811166157.6A CN201811166157A CN109358352A CN 109358352 A CN109358352 A CN 109358352A CN 201811166157 A CN201811166157 A CN 201811166157A CN 109358352 A CN109358352 A CN 109358352A
- Authority
- CN
- China
- Prior art keywords
- occultation
- spacecraft
- elv
- speed
- antenna
- 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
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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/08—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing integrity information, e.g. health of satellites or quality of ephemeris data
-
- 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
- G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种星载实时掩星预报方法,包括先根据掩星探测仪的定位设备获取航天器的位置速度,并计算航天器轨道参数;然后计算航天器和所有导航星在10s后的位置速度,并计算所有导航星的俯仰角、方位角、切点高度及航天器的天线方位角,判断得到的俯仰角Elv、切点高度Hight及天线方位角是否满足掩星事件的门限,若满足,则将时间外推30s和1min,分别重复上述计算和判断过程。若全部满足,则判定该导航星为掩星。本发明通过惯导设备采用连续弧段判断掩星事件,增强了预报结果的准确性,方法简单易行,对于实时预报具有重要应用价值。
Description
技术领域
本发明涉及大气掩星探测技术领域,具体涉及一种星载实时掩星预报方法,利用星载大气掩星探测设备和初始化硬件设备,提前预报掩星事件发生,为快速捕获掩星事件提供初始条件的方法。
背景技术
GPS/MET实验的成功,验证了GNSS无线电掩星技术对于地球大气探测的具有重大价值。随着GNSS无线电掩星技术的发展,掩星技术被越来越多国家重视。掩星事件是指GNSS卫星的导航信号经过地球表面的电离层和大气层而产生弯折的物理变化现象。掩星技术是指掩星事件发生时,掩星信号被掩星接收机捕获并跟踪,从而通过分析采集到的观测量来反演信号传播路径中的一些参量的技术。
星载掩星接收机和掩星相对速度较大,因此掩星现象通常持续时间短,因此对于星载掩星探测仪捕获掩星事件,就需快速捕获,利用掩星预报算法,可以提高掩星捕获速度,另外对于不同的轨道掩星事件又具有随机性,因此结合天线指向,利用惯导设备获取当前姿态,通过精准预报,也能够合理设置通道,节省硬件资源。
现有掩星预报算法主要是运用仿真软件计算,然后提前置入设备中,未考虑实际应用中的可见性和实时性,因此准确率不足,从而导致浪费硬件资源。
发明内容
本发明要解决的技术问题是提供一种星载实时掩星预报方法,克服现有技术的不足,通过轨道动力学方程以及坐标系转换,在保证GNSS导航星可视的情况下,计算连续弧段内GNSS导航星是否满足掩星事件要求,解决了星载掩星探测仪实时精确预报的需求,有效的提高计算效率和预报准确度。
为解决上述技术问题,本发明提供了一种星载实时掩星预报方法,步骤如下:
(1)获取航天器在WGS-84坐标系下的轨道参数
通过掩星探测仪的定位设备获取航天器的位置速度P=[x,y,z],V=[vx,vy,vz]和当前时间T,并根据该位置速度计算航天器的轨道六根数S;
所述轨道六根数S由以下公式:
构建中间变量,
O=[cos S[1] sin S[1] 0]
Od=[-cos S[0]*O[1] cos S[0]*O[0] sin S[0]]
构建中间变量,
S[4]=Phi-S[2]*sin Phi-S[3]*cos Phi
给出,其中r,v为位置速度的模;
RR为法向矢量,并将RR单位化,计算公式如下:
RR=P×V
其中P=[x,y,z],V=[vx,vy,vz],为航天器的位置速度。
(2)获取航天器和所有GNSS导航星在10s后的位置速度
航天器10s后的位置速度获取方法:利用步骤(1)中获取的轨道六根数S,根据轨道动力学建立航天器运动方程,以当前时间T+10s带入,计算10s后的轨道六根数S10,再利用轨道根数转直角坐标系公式,可以获得10s后的航天器位置速度P10,V10;
各GNSS导航星10s后的位置速度获取方法:通过掩星探测仪的定位设备获取的各GNSS导航星的广播电文,根据各星座公布ICD接口文件获得所有导航星在当前时间T之后10S的位置速度PGNSS10,VGNSS10。
(3)获取所有GNSS导航星的俯仰角Elv、方位角及切点高度h
方位角和俯仰角Ely,计算公式如下:
式中xRELnew,yRELnew,zRELnew是10s后GNSS导航星相对于航天器的相对矢量RREL10转化为东北天坐标下的相对矢量RRELnew的各轴分量,由以下转化公式计算而来:
RRELnew=RREL10*C
式中C为转换矩阵,用来将WGS-84坐标系下航天器的位置速度计算转换到东北天坐标系;RREL10为GNSS导航星相对于航天器的相对矢量,由如下公式计算得到:
RREL10=PGNSS10-P10
式中P10为步骤(2)中得到的航天器10s后的位置速度;PGNSS10为步骤(2)中得到的各GNSS导航星在当前时间T之后10S的位置速度;
切点高度h的计算,通过如下的面积公式来计算:
h*||RREL||=||P10||*||PGNSS10||*sin(θ)
式中θ为LEO-地心-GNSS之间的夹角,通过三角形余弦定理计算获得;
再计算切点-地心-LEO之间的夹角运用的公式如下:
以角度旋转GNSS位置坐标到切点位置,计算切点坐标RTanget,并转化为经纬高[Lat,Lon,Hight]。
(4)判断是否符合门限
利用步骤(3)中得到的俯仰角Elv和切点高度Hight判断是否满足掩星事件的门限;如果满足掩星事件门限,则再利用惯导姿态信息计算航天器天线相对于东北天的方位角,判断导航星是否满足天线方位角约束。
(5)如果步骤(4)中的天线方位角满足约束,则外推多个固定时点,分别继续计算步骤(1)~(3)在选定的固定时点时是否满足步骤(4);若全部满足,则判定该导航星为掩星。
所述步骤(4)中的掩星事件的门限为:Elvmin≤Elv≤Elvmax,Hightmin≤Hight≤Hightmax;其中:Elvmax为掩星的俯仰角门限上限,Elvmin为掩星的俯仰角门限下限:所述Elvmax及Elvmin由公式Elvmax=w/2+ae和Elvmin=Elvtanget-2计算,式中w为天线的幅角,ae为天线安装角俯仰角,az为方位角,Elvtanget为掩星事件切地球的最低俯仰角。所述最低俯仰角Elvtanget由公式计算,式中切点最低高度a为地球半径,H为航天器的高度。
Hightmin和Hightmax则根据所需要的观测的大气高度选取。如果关注底层大气,可设置为[1,10]Km,如果关注中层大气,可设置为[10,80]Km。
所述步骤(4)中的天线方位角约束具体为:天线方位角az应该在天线幅角w范围内,即掩星实际的方位角和天线的安装角方位角的差dangle的绝对值小于w/2,天线的安装方位角由天线的安装角和惯导提供的姿态角相加获得。
所述步骤(5)中的固定时点为30S和1min。
本发明与现有技术相比的有益效果:
(1)本发明采用连续弧段判断掩星事件,避免了单次判断造成的误判或者漏判;
(2)本发明通过惯导设备来判断导航星的可见性,增强了预报结果的准确性,方法简单易行,对于实时预报具有重要应用价值;
(3)本发明综合考虑了卫星实际的姿态和掩星天线的安装方向,增强了掩星事件预报的准确性。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细地描述,其中:
图1为本发明在WGS-84坐标系下的几何示意图。
图2为本发明星载实时掩星预报方法在水平面的投影图。
图3为本发明星载实时掩星预报方法的流程图。
图4为本发明的实施场景设计。
具体实施方式
下面结合实施例附图,对本发明的技术方案进行清楚、完整的描述,显然,描述的实施例仅仅是本发明的一个具体的实施例,不是全部的实施例。下述实施例是说明性的,不是限制性的,不能以下述实施例来限定本发明的保护范围。
以下结合图1到图3阐述本发明星载实时掩星预报方法,详细步骤如下:
(1)获取航天器在WGS-84坐标系下的轨道参数
通过掩星探测仪的定位设备获取航天器的位置速度P=[x,y,z],V=[vx,vy,vz]和当前时间T,并根据该位置速度计算航天器的轨道六根数S;
所述轨道六根数S由以下公式:
构建中间变量,
O=[cos S[1] sin S[1] 0]
Od=[-cos S[0]*O[1] cos S[0]*O[0] sin S[0]]
构建中间变量,
S[4]=Phi-S[2]*sin Phi-S[3]*cos Phi
给出,其中r,v为位置速度的模;
RR为法向矢量,并将RR单位化,计算公式如下:
RR=P×V
其中P=[x,y,z],V=[vx,vy,vz],为航天器的位置速度;
(2)获取航天器和所有GNSS导航星在10s后的位置速度
航天器10s后的位置速度获取方法:利用步骤(1)中获取的轨道六根数S,根据轨道动力学建立航天器运动方程,以当前时间T+10s带入,计算10s后的轨道六根数S10,再利用轨道根数转直角坐标系公式,可以获得10s后的航天器位置速度P10,V10;
各GNSS导航星10s后的位置速度获取方法:通过掩星探测仪的定位设备获取的各GNSS导航星的广播电文,根据各星座公布ICD接口文件获得所有导航星在当前时间T之后10S的位置速度PGNSS10,VGNSS10;
(3)获取所有GNSS导航星的俯仰角Elv、方位角Azm及切点高度h
方位角和俯仰角Elv,计算公式如下:
式中xRELnew,yRELnew,zRELnew是10s后GNSS导航星相对于航天器的相对矢量RREL10转化为东北天坐标下的相对矢量RRELnew的各轴分量,由以下转化公式计算而来:
RRELnew=RREL10*C
式中C为转换矩阵,用来将WGS-84坐标系下航天器的位置速度计算转换到东北天坐标系;RREL10为GNSS导航星相对于航天器的相对矢量,由如下公式计算得到:
RREL10=PGNSS10-P10
式中P10为步骤(2)中得到的航天器10s后的位置速度;PGNSS10为步骤(2)中得到的各GNSS导航星在当前时间T之后10S的位置速度;
切点高度h的计算,通过如下的面积公式来计算:
h*||RREL||=||P10||*||PGNSS10||*sin(θ)
式中θ为LEO-地心-GNSS之间的夹角,通过三角形余弦定理计算获得;
再计算切点-地心-LEO之间的夹角
以角度旋转GNSS位置坐标到切点位置,计算切点坐标RTanget,并转化为经纬高[Lat,Lon,Hight];
(4)判断是否符合门限
利用步骤(3)中得到的俯仰角Elv和切点高度Hight判断是否满足掩星事件的门限;如果满足掩星事件门限,则再利用惯导姿态信息计算航天器天线相对于东北天的方位角,判断导航星是否满足天线方位角Azm约束;
(5)如果步骤(4)中的天线方位角满足约束,则外推30S和1min,分别继续计算步骤(1)~(3)在30S和1min时是否满足步骤(4);若全部满足,则判定该导航星为掩星。
所述步骤(4)中的掩星事件的判断门限为:Elvmin≤Elv≤Elvmax,Hightmin≤Hight≤Hightmax;
其中:Elvmax为掩星的俯仰角门限上限,Elvmin为掩星的俯仰角门限下限:所述Elvmax及Elvmin由如下公式计算:
Elvmax=w/2+ae
E1vmin=Elvtanget-2
其中:w为天线的幅角,ae为天线安装角俯仰角,az为方位角,Elvtanget为掩星事件切地球的最低俯仰角;所述最低俯仰角Elvtanget由如下公式计算:
其中:切点最低高度a为地球半径,H为航天器的高度;
其中Hightmin和Hightmax根据所需要的观测的大气高度选取,对于底层大气关注的话,可设置为[1,10]Km,对于中层大气关注的话,可设置为[10,80]Km。
该方法中,步骤(5)中的天线方位角约束具体为:方位角az应该在天线幅角w范围内,即掩星实际的方位角和天线的安装角方位角的差dangle的绝对值小于w/2,天线的安装方位角由天线的安装角和惯导提供的姿态角相加获得。
实施例
仿真场景如图4所示设计,卫星设计轨道高度400km,轨道倾角设计为55度,姿态设计为对地定向,三轴稳定,UTC时间设为2017年8月5号4:00:00至2017年8月6号4:00:00,导航星座以GPS星座为例,星数选区Prn1,航天器接收天线幅角w为40度,天线安装俯仰角为-20,掩星切角为-30,则天线俯仰角门限为[-22,0],由于场景设计为三轴稳定,所以天线和运动方向保持重合,因此天线方位角门限为[-20,20],如果天线有姿态机动,需引入惯导设备,加上姿态机动导致的天线方位变化。以时间60s为步长,通过本发明计算全天Prn1号GPS卫星掩星事件如下表所示:
通过半实物仿真,得出以上数据,可以证实掩星探测仪可以根据预报结果快速、准确跟踪掩星事件。
本发明提供的星载实时掩星预报方法,增强了预报结果的准确性,方法简单易行,对于掩星实时预报具有重要应用价值。
上面结合附图对本发明的实施方式作了详细说明,本发明说明书中未作详细阐述的内容属于本领域专业技术人员的公知技术。
Claims (4)
1.一种星载实时掩星预报方法,其特征在于,包括以下步骤:
(1)获取航天器在WGS-84坐标系下的轨道参数
通过掩星探测仪的定位设备获取航天器的位置速度P=[x,y,z],V=[vx,vy,vz]和当前时间T,并根据该位置速度计算航天器的轨道六根数S;
所述轨道六根数S由以下公式:
构建中间变量,
O=[cos S[1] sin S[1] 0]
Od=[-cos S[0]*O[1] cos S[0]*O[0] sin S[0]]
构建中间变量,
S[4]=Phi-S[2]*sin Phi-S[3]*cos Phi
给出,其中r,v为位置速度的模;
RR为法向矢量,并将RR单位化,计算公式如下:
RR=P×V
其中P=[x,y,z],V=[vx,vy,vz],为航天器的位置速度;
(2)获取航天器和所有GNSS导航星在10s后的位置速度
航天器10s后的位置速度获取方法:利用步骤(1)中获取的轨道六根数S,根据轨道动力学建立航天器运动方程,以当前时间T+i0s带入,计算10s后的轨道六根数S10,再利用轨道根数转直角坐标系公式,可以获得10s后的航天器位置速度P10,V10;
各GNSS导航星10s后的位置速度获取方法:通过掩星探测仪的定位设备获取的各GNSS导航星的广播电文,根据各星座公布ICD接口文件获得所有导航星在当前时间T之后10S的位置速度PGNSS10,VGNSS10;
(3)获取所有GNSS导航星的俯仰角Elv、方位角Azm及切点高度h
方位角和俯仰角Elv,计算公式如下:
式中xRELnew,是10s后GNSS导航星相对于航天器的相对矢量RREL10转化为东北天坐标下的相对矢量RRELnew的各轴分量,由以下转化公式计算而来:
RRELnew=RREL10*C
式中C为转换矩阵,用来将WGS-84坐标系下航天器的位置速度计算转换到东北天坐标系;RREL10为GNSS导航星相对于航天器的相对矢量,由如下公式计算得到:
RREL10=PGNSS10-P10
式中P10为步骤(2)中得到的航天器10s后的位置速度;PGNSS10为步骤(2)中得到的各GNSS导航星在当前时间T之后10S的位置速度;
切点高度h的计算,通过如下的面积公式来计算:
h*||RREL||=||P10||*||PGNSS10||*sin(θ)
式中θ为LEO-地心-GNSS之间的夹角,通过三角形余弦定理计算获得;
再计算切点-地心-LEO之间的夹角运用的公式如下:
以角度旋转GNSS位置坐标到切点位置,计算切点坐标RTanget,并转化为经纬高[Lat,Lon,Hight];
(4)判断是否符合门限
利用步骤(3)中得到的俯仰角Elv和切点高度Hight判断是否满足掩星事件的门限;如果满足掩星事件门限,则再利用惯导姿态信息计算航天器天线相对于东北天的方位角,判断导航星是否满足天线方位角约束;
(5)如果步骤(4)中的天线方位角满足约束,则外推多个固定时点,分别继续计算步骤(1)~(3)在选定的固定时点时是否满足步骤(4);若全部满足,则判定该导航星为掩星。
2.按照权利要求1所述的星载实时掩星预报方法,其特征在于,所述步骤(4)中的掩星事件的门限为:Elvmin≤Elv≤Elvmax,HightminHight≤Hightmax;
其中:Elvmax为掩星的俯仰角门限上限,Elvmin为掩星的俯仰角门限下限:所述Elvmax及Elvmin由如下公式计算:
Elvmax=w/2-ae
Elvmin=Elvtanget-2
其中:w为天线的幅角,ae为天线安装角俯仰角,az为方位角,Elvtanget为掩星事件切地球的最低俯仰角;所述最低俯仰角Elvtanget由如下公式计算:
其中:切点最低高度a为地球半径,H为航天器的高度;
其中Hightmin和Hightmax根据所需要的观测的大气高度选取,对于底层大气关注的话,可设置为[1,10]Km,对于中层大气关注的话,可设置为[10,80]Km。
3.按照权利要求1所述的星载实时掩星预报方法,其特征在于,所述步骤(4)中的天线方位角约束具体为:天线方位角az应该在天线幅角w范围内,即掩星实际的方位角和天线的安装角方位角的差dangle的绝对值小于w/2,天线的安装方位角由天线的安装角和惯导提供的姿态角相加获得。
4.按照权利要求1所述的星载实时掩星预报方法,其特征在于,所述步骤(5)中的固定时间段为30S和1min。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811166157.6A CN109358352B (zh) | 2018-09-30 | 2018-09-30 | 一种星载实时掩星预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811166157.6A CN109358352B (zh) | 2018-09-30 | 2018-09-30 | 一种星载实时掩星预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358352A true CN109358352A (zh) | 2019-02-19 |
CN109358352B CN109358352B (zh) | 2021-05-18 |
Family
ID=65348341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811166157.6A Active CN109358352B (zh) | 2018-09-30 | 2018-09-30 | 一种星载实时掩星预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358352B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109946728A (zh) * | 2019-03-13 | 2019-06-28 | 北京遥感设备研究所 | 一种适用于卫星用户站数字跟踪接收机的程序跟踪方法 |
CN111366953A (zh) * | 2020-03-20 | 2020-07-03 | 中国西安卫星测控中心 | 一种北斗Ka复杂相控阵天线跟踪可见报生成算法 |
CN111949616A (zh) * | 2020-09-08 | 2020-11-17 | 天津云遥宇航科技有限公司 | 一种gnss掩星数据地面实时反演演示系统 |
CN112346086A (zh) * | 2020-10-23 | 2021-02-09 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
CN114413911A (zh) * | 2022-02-09 | 2022-04-29 | 中国科学院微小卫星创新研究院 | 掩星探测自主任务姿态导引方法及系统 |
CN115826007A (zh) * | 2023-02-17 | 2023-03-21 | 天津云遥宇航科技有限公司 | 一种星基全球导航卫星系统中性大气掩星事件预报方法 |
CN116299574A (zh) * | 2023-05-11 | 2023-06-23 | 天津云遥宇航科技有限公司 | 一种基于高度角的glonass掩星对应参考星prn校正方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288974A (zh) * | 2011-05-10 | 2011-12-21 | 航天恒星科技有限公司 | 一种卫星导航掩星信号生成系统 |
CN103542854A (zh) * | 2013-11-02 | 2014-01-29 | 中国人民解放军国防科学技术大学 | 基于星载处理器的自主定轨方法 |
CN103605140A (zh) * | 2013-11-06 | 2014-02-26 | 航天恒星科技有限公司 | 一种掩星信号快速恢复方法 |
CN103852812A (zh) * | 2012-12-07 | 2014-06-11 | 深圳航天东方红海特卫星有限公司 | 一种前向模拟掩星方法 |
WO2014121197A2 (en) * | 2013-02-01 | 2014-08-07 | NanoSatisfi Inc. | System and method for widespread low cost orbital satellite access |
CN104216031A (zh) * | 2014-09-15 | 2014-12-17 | 航天东方红卫星有限公司 | 一种可变步长掩星预报方法 |
CN104332707A (zh) * | 2014-10-27 | 2015-02-04 | 西安空间无线电技术研究所 | 一种用于低轨星载天线跟踪地面站的方法 |
CN104933316A (zh) * | 2015-07-02 | 2015-09-23 | 华东交通大学 | 一种基于双参数混合正则化的电离层掩星反演方法 |
CN105891856A (zh) * | 2016-04-07 | 2016-08-24 | 中国科学院国家空间科学中心 | 一种gnss无线电掩星探测仪 |
CN106291560A (zh) * | 2016-10-27 | 2017-01-04 | 中国科学院国家空间科学中心 | 一种兼容gnss‑r探测技术的gnss无线电掩星探测仪 |
CN106772445A (zh) * | 2016-11-25 | 2017-05-31 | 北京天工科仪空间技术有限公司 | 一种电离层活动空间观测系统 |
CN107238846A (zh) * | 2017-04-25 | 2017-10-10 | 清华大学 | 一种基于glonass历书参数的卫星位置与速度预报方法 |
CN107561557A (zh) * | 2017-08-18 | 2018-01-09 | 中国科学院遥感与数字地球研究所 | 一种掩星探测仪大气成分反演方法 |
US20180038994A1 (en) * | 2016-08-02 | 2018-02-08 | International Business Machines Corporation | Techniques to Improve Global Weather Forecasting Using Model Blending and Historical GPS-RO Dataset |
CN108519110A (zh) * | 2018-04-26 | 2018-09-11 | 北京空间飞行器总体设计部 | 基于图像信息的空间非合作目标自主相对导航在轨验证系统 |
-
2018
- 2018-09-30 CN CN201811166157.6A patent/CN109358352B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288974A (zh) * | 2011-05-10 | 2011-12-21 | 航天恒星科技有限公司 | 一种卫星导航掩星信号生成系统 |
CN103852812A (zh) * | 2012-12-07 | 2014-06-11 | 深圳航天东方红海特卫星有限公司 | 一种前向模拟掩星方法 |
WO2014121197A2 (en) * | 2013-02-01 | 2014-08-07 | NanoSatisfi Inc. | System and method for widespread low cost orbital satellite access |
CN103542854A (zh) * | 2013-11-02 | 2014-01-29 | 中国人民解放军国防科学技术大学 | 基于星载处理器的自主定轨方法 |
CN103605140A (zh) * | 2013-11-06 | 2014-02-26 | 航天恒星科技有限公司 | 一种掩星信号快速恢复方法 |
CN104216031A (zh) * | 2014-09-15 | 2014-12-17 | 航天东方红卫星有限公司 | 一种可变步长掩星预报方法 |
CN104332707A (zh) * | 2014-10-27 | 2015-02-04 | 西安空间无线电技术研究所 | 一种用于低轨星载天线跟踪地面站的方法 |
CN104933316A (zh) * | 2015-07-02 | 2015-09-23 | 华东交通大学 | 一种基于双参数混合正则化的电离层掩星反演方法 |
CN105891856A (zh) * | 2016-04-07 | 2016-08-24 | 中国科学院国家空间科学中心 | 一种gnss无线电掩星探测仪 |
US20180038994A1 (en) * | 2016-08-02 | 2018-02-08 | International Business Machines Corporation | Techniques to Improve Global Weather Forecasting Using Model Blending and Historical GPS-RO Dataset |
CN106291560A (zh) * | 2016-10-27 | 2017-01-04 | 中国科学院国家空间科学中心 | 一种兼容gnss‑r探测技术的gnss无线电掩星探测仪 |
CN106772445A (zh) * | 2016-11-25 | 2017-05-31 | 北京天工科仪空间技术有限公司 | 一种电离层活动空间观测系统 |
CN107238846A (zh) * | 2017-04-25 | 2017-10-10 | 清华大学 | 一种基于glonass历书参数的卫星位置与速度预报方法 |
CN107561557A (zh) * | 2017-08-18 | 2018-01-09 | 中国科学院遥感与数字地球研究所 | 一种掩星探测仪大气成分反演方法 |
CN108519110A (zh) * | 2018-04-26 | 2018-09-11 | 北京空间飞行器总体设计部 | 基于图像信息的空间非合作目标自主相对导航在轨验证系统 |
Non-Patent Citations (6)
Title |
---|
S. B. HEALYA: "Assimilating GPS radio occultation measurements with two-dimensional bending angle observation operators", 《QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY》 * |
Y-H. KUO: "Assimilation of GPS Radio Occultation Data for Numerical Weather Prediction", 《SPECIAL ISSUE OF TERRESTRIAL, ATMOSPHERIC AND OCEANIC SCIENCE》 * |
冉隆燧: "《航天工程设计实践》", 31 August 2012 * |
吕秋杰: "掩星探测的发生条件及预报算法研究", 《2010年第二十三届全国空间探测学术交流会论文》 * |
周文: "利用COSMIC掩星资料研究青藏高原地区", 《气象学报》 * |
张训械: "山基无线电掩星模拟", 《电波科学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109946728A (zh) * | 2019-03-13 | 2019-06-28 | 北京遥感设备研究所 | 一种适用于卫星用户站数字跟踪接收机的程序跟踪方法 |
CN109946728B (zh) * | 2019-03-13 | 2021-10-26 | 北京遥感设备研究所 | 一种适用于卫星用户站数字跟踪接收机的程序跟踪方法 |
CN111366953B (zh) * | 2020-03-20 | 2022-11-18 | 中国西安卫星测控中心 | 一种北斗Ka复杂相控阵天线跟踪可见报生成算法 |
CN111366953A (zh) * | 2020-03-20 | 2020-07-03 | 中国西安卫星测控中心 | 一种北斗Ka复杂相控阵天线跟踪可见报生成算法 |
CN111949616A (zh) * | 2020-09-08 | 2020-11-17 | 天津云遥宇航科技有限公司 | 一种gnss掩星数据地面实时反演演示系统 |
CN111949616B (zh) * | 2020-09-08 | 2023-05-26 | 天津云遥宇航科技有限公司 | 一种gnss掩星数据地面实时反演演示系统 |
CN112346086A (zh) * | 2020-10-23 | 2021-02-09 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
CN112346086B (zh) * | 2020-10-23 | 2024-03-15 | 西安空间无线电技术研究所 | 一种基于临近空间浮空平台的高效快速掩星预报方法 |
CN114413911A (zh) * | 2022-02-09 | 2022-04-29 | 中国科学院微小卫星创新研究院 | 掩星探测自主任务姿态导引方法及系统 |
CN114413911B (zh) * | 2022-02-09 | 2023-09-22 | 中国科学院微小卫星创新研究院 | 掩星探测自主任务姿态导引方法及系统 |
CN115826007A (zh) * | 2023-02-17 | 2023-03-21 | 天津云遥宇航科技有限公司 | 一种星基全球导航卫星系统中性大气掩星事件预报方法 |
CN116299574A (zh) * | 2023-05-11 | 2023-06-23 | 天津云遥宇航科技有限公司 | 一种基于高度角的glonass掩星对应参考星prn校正方法 |
CN116299574B (zh) * | 2023-05-11 | 2023-08-15 | 天津云遥宇航科技有限公司 | 一种基于高度角的glonass掩星对应参考星prn校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109358352B (zh) | 2021-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109358352A (zh) | 一种星载实时掩星预报方法 | |
CN109059916A (zh) | 一种基于惯导的浮空器掩星预报方法 | |
KR102049371B1 (ko) | 선박 보조 도킹 방법과 시스템 | |
Abulude et al. | Global Positioning System And It's Wide Applications | |
CN101893440A (zh) | 基于星敏感器的天文自主导航方法 | |
JP2009041932A (ja) | 移動体用測位装置 | |
CN202471085U (zh) | 一种导航装置 | |
CN103412198A (zh) | 船舶防护电场的三维空间分布特性测量装置及测量方法 | |
CN104216031B (zh) | 一种可变步长掩星预报方法 | |
CN108896957A (zh) | 一种无人机操控信号源的定位系统及方法 | |
CN101140322A (zh) | 用于布局卫星导航地基区域增强系统的装置及方法 | |
CN112762935B (zh) | 一种基于船体姿态监测的定位方法及系统 | |
CN111060111A (zh) | 一种低轨卫星入轨初期定轨方法 | |
CN112945242B (zh) | 一种卫星在轨自主规划任务最佳时间与姿态的方法 | |
CN116996115B (zh) | 一种低轨卫星接收时间窗口计算方法、装置及设备 | |
CN205374736U (zh) | 一种基于动态监测的海上溢油预测系统 | |
CN104199072A (zh) | 空间飞行物发现与轨迹标记系统及其实现轨迹标记的方法 | |
CN114659496B (zh) | 一种用于船载北斗一体机倾斜监测的方法 | |
CN106324600A (zh) | 风场测量方法及装置 | |
Trzuskowsky et al. | ANCHOR: navigation, routing and collision warning during operations in harbors | |
CN115291264A (zh) | 一种适用于公路隧道高速行驶车辆的北斗定位信号模拟方法 | |
González-Rocha et al. | Multirotor-assisted measurements of wind-induced drift of irregularly shaped objects in aquatic environments | |
CN113687398A (zh) | 一种基于北斗差分定位的无人机跟瞄通信系统及方法 | |
Kim et al. | 3D Grid-based Global Positioning System Satellite Signal Shadowing Range Modeling in Urban Area. | |
Strandjord et al. | Framework and techniques for cooperative group situational awareness in urban environments |
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 |