CN104678408B - 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用*** - Google Patents
星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用*** Download PDFInfo
- Publication number
- CN104678408B CN104678408B CN201510060813.4A CN201510060813A CN104678408B CN 104678408 B CN104678408 B CN 104678408B CN 201510060813 A CN201510060813 A CN 201510060813A CN 104678408 B CN104678408 B CN 104678408B
- Authority
- CN
- China
- Prior art keywords
- satellite
- time
- navigation receiver
- borne
- epoch
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000006872 improvement Effects 0.000 claims abstract description 31
- 230000001360 synchronised effect Effects 0.000 claims abstract description 7
- 230000001174 ascending effect Effects 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims description 31
- 230000000630 rising effect Effects 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000036962 time dependent Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 238000010295 mobile communication Methods 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
-
- 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
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
- G04R20/04—Tuning or receiving; Circuits therefor
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
- G04R20/06—Decoding time data; Circuits therefor
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种星载导航接收机授时方法,包括:产生秒脉冲;在所述秒脉冲上升沿采集导航卫星信号中的卫星发射时刻和观测载波相位;根据所述卫星发射时刻和观测载波相位计算伪距作为观测伪距,构建解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差;对所述频差和钟差进行质量改善处理;根据改善后的频差和钟差调整所述秒脉冲的相位和频率,使所述秒脉冲上升沿与标准时间整秒起始时刻同步;将所述调整后的秒脉冲广播给卫星上的其它电子***;在下一个所述调整后的秒脉冲到来之前,将所述位置、速度、时间广播给卫星上的所述其它电子***;以及生成和更新与时间相关的遥测量,并在导航卫星请求对应遥测量时给出应答。
Description
技术领域
本发明涉及一种导航接收机授时方法与***,尤其是一种星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用***。
背景技术
在各种无线电定位技术中,全球导航卫星***(Global Navigation SatelliteSystem, GNSS)是最基本的手段,它是所有导航卫星***的全称,目前主要包括美国的全球定位***(Global Positioning System, GPS),俄罗斯的全球导航卫星***(GlobalNavigation Satellite System, GLONASS),欧洲的伽利略***(Galileo),中国的北斗(Compass)。全球导航卫星***接收机基本工作原理是:接收到导航卫星发送无线电信号并提取伪矩、载波相位等观测量,并根据来自4个以上卫星的观测量进行解算,解算的结果是接收机的位置、速度和时间。
卫星导航接收机通常只是某一个设备的一部分,这一设备称作宿主设备。低轨道卫星上通常安装星载导航接收机,这时低轨道卫星就是宿主设备。星载导航接收机是安装在低轨道卫星上并接收全球导航卫星***的卫星信号并为宿主低轨道卫星提供位置、速度和时间的卫星导航接收机。低轨道卫星通常轨道高度不超过1000公里。近年来,星载导航接收机也越来越多的应用在更高轨道的接收机上。
卫星导航接收机为宿主设备提供时间的功能称作授时。授时需求广泛存在,比如在轨道卫星移动通信中,为了实现轨道卫星与地面设备的频率和时间一致性,往往需要星载导航接收机实现地面静态授时型导航接收机一样的高授时精度。随着移动互联网的迅速发展,轨道卫星移动通信***的核心应用向互联网变迁。地面移动通信***往往被等价的称作移动互联网,轨道卫星移动通信***也常常被称作空间互联网。随着空间互联网的发展,轨道卫星移动通信***正在迎来新一轮更大的发展,因此星载授时型导航接收机将会得到更多的使用。但在轨道卫星移动通信***中,由于宿主卫星以接近第一宇宙速度的速度运行,导致导航接收机环路带宽更宽、噪声更大,这为授时型星载导航接收机带来了比传统导航接收机的更多挑战。
发明内容
有鉴于此,确有必要提供一种授时精度高且适用于有精密授时需求的低轨卫星的星载导航接收机授时方法以及授时型星载导航接收机。
一种星载导航接收机授时方法,包括以下步骤:
S1,星载导航接收机产生秒脉冲;
S2,星载导航接收机在所述秒脉冲上升沿采集导航卫星信号中的观测信息,所述观测信息包括卫星发射时刻和观测载波相位;
S3,根据所述卫星发射时刻和观测载波相位计算伪距作为观测伪距,构建解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差;
S4,对所述频差进行质量改善处理,具体包括以下步骤:
S41,设定初始化历元数;
S42,比较当前历元数和所述初始化历元数;
S43,当时,设定经过质量改善处理的当前历元的频差等于当前历元的频差,即:;以及
S44,当时,设定经过质量改善处理的当前历元的频差,其中,为滤波系数,且;
S5,对所述钟差进行质量改善处理,包括以下步骤:
S51,设定初始化历元数以及钟差限定量;
S52,比较当前历元数和初始化历元数;
S53,当时,设定经过质量改善处理的当前历元的钟差等于当前历元的钟差,即:;以及
S54,当时,判断当前历元的与差值是否超出最大范围;
S55,若超过了所述最大范围,且成立,则令;
S56,若超过了所述最大范围,且成立,则令,以及;
S57,若未超出所述最大范围,则令;
S6,根据改善后的频差和钟差调整所述秒脉冲的相位和频率,使所述秒脉冲上升沿与标准时间整秒起始时刻同步;
S7,星载导航接收机将所述调整后的秒脉冲广播给卫星上的其它电子***;
S8,在下一个所述调整后的秒脉冲到来之前,所述星载导航接收机将所述位置、速度、时间广播给卫星上的所述其它电子***;以及
S9,生成和更新与时间相关的遥测量,并在导航卫星请求对应遥测量时给出应答。
一种用于实现上述星载导航接收机授时方法的授时型星载导航接收机,包括:
天线,用于接收卫星导航信号;
射频前端模块,与所述天线连接,接收并处理所述卫星导航信号,所述处理包括滤波、变频以及模数采样;
捕获跟踪模块,与所述射频前端模块连接,接收经所述射频前端处理后的卫星导航信号,并捕获和跟踪目标卫星的导航信号;
观测量获取模块,接收所述跟踪后的目标卫星的导航信号以及所述秒脉冲,并在所述秒脉冲的上升沿采集该跟踪后的目标卫星的导航信号,获得所述观测信息;
解算模块,接收所述观测信息,并根据所述卫星发射时刻和观测载波相位计算所述伪距,构建所述解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差;
钟差和频差质量改善模块,用对所述解算模块解算获得的频差和钟差进行质量改善处理;
遥测模块,用于生成和更新与时间相关的所述遥测量,并在星上计算机请求对应遥测量时给出应答;以及
秒脉冲模块,用于产生所述秒脉冲,并根据所述质量改善处理后的频差和钟差调节所述秒脉冲的频率和相位,使所述秒脉冲的上升沿与标准时间整秒起始时刻同步。
一种星载导航应用***,包括:星上计算机、测控通信模块、被授时载荷以及地面通信***,该星载导航应用***进一步包括所述授时型星载导航接收机。
与现有技术相比较,本发明实施例提供的星载导航接收机授时方法与授时型星载导航接收机通过充分利用轨道动力学获得高精度的定位结果和原始钟差;并采用所述质量改善处理方法改善钟差和频差,有效地提高了授时精度以及星载导航接收机的***可靠性,尤其适用于有精密授时需求的低轨卫星。此外,本发明实施例提供的星载导航接收机授时方法与授时型星载导航接收机经过在轨飞行验证,72小时授时精度方均根误差小于12ns。
附图说明
图1为本发明实施例提供的星载导航接收机授时方法的流程图。
图2为本发明实施例提供的星载导航接收机授时方法中解算方法的流程图。
图3为本发明实施例提供的星载导航接收机授时方法中频差质量改善处理的流程图。
图4为本发明实施例提供的星载导航接收机授时方法中钟差质量改善处理的流程图。
图5为本发明实施例提供的授时型星载导航接收机以及星载导航应用***的功能连接框图。
主要元件符号说明
授时型星载导航接收机 | 100 |
天线 | 102 |
射频前端模块 | 104 |
捕获和跟踪模块 | 106 |
观测量获取模块 | 108 |
解算模块 | 110 |
钟差和频差质量改善模块 | 112 |
遥测模块 | 114 |
秒脉冲模块 | 116 |
星载导航应用*** | 10 |
星上计算机 | 200 |
测控通信模块 | 300 |
被授时载荷 | 400 |
地面通信*** | 500 |
星上数据总线 | 600 |
秒脉冲数据总线 | 700 |
如下具体实施方式将结合上述附图进一步说明本发明。
具体实施方式
以下将结合附图详细说明本发明实施例提供的星载导航接收机授时方法和授时型星载导航接收机。
请参阅图1,本发明实施例提供一种星载导航接收机授时方法,包括以下步骤:
S1,星载导航接收机产生秒脉冲;
S2,星载导航接收机在所述秒脉冲上升沿采集导航卫星信号中的观测信息,所述观测信息包括卫星发射时刻和观测载波相位;
S3,根据所述卫星发射时刻和观测载波相位计算伪距作为观测伪距,构建解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差;
S4,对所述频差进行质量改善处理;
S5,对所述钟差进行质量改善处理;
S6,根据改善后的频差和钟差调整所述秒脉冲的相位和频率,使所述秒脉冲上升沿与标准时间整秒起始时刻同步;
S7,星载导航接收机将所述调整后的秒脉冲广播给卫星上的其它电子***;
S8,在下一个所述调整后的秒脉冲到来之前,所述星载导航接收机将所述位置、速度、时间广播给卫星上的所述其它电子***;以及
S9,生成和更新与时间相关的遥测量,并在导航卫星请求对应遥测量时给出应答。
在上述步骤S3中,所述星载导航接收机可利用所述卫星发射时刻和观测载波相位计算伪距,然后构建解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差。
请参阅图2,所述步骤S3具体包括:
S31,设定初始化历元数;
S32,比较当前历元数和所述初始化历元数;
S33,当时,进行初始化,具体包括:使用迭代最小二乘法获得当前历元状态量,所述当前历元状态量包括位置、速度、钟差以及频差;
S34,当时,进行异常观测量剔除以获得保留伪距和保留载波相位;
S35,根据所述保留伪距和保留载波相位构建解算方程,并进行解算,得到当前历元的状态量;
S36,对所述当前历元的状态量进行自主完整性校验,若校验通过,则采纳作为当前历元的状态量,否则,采纳作为当前历元的状态量;以及
S37,合理性校验:若发生连续自主完整性校验不通过或定位结果明显不合理,则回到步骤S31,重新开始初始化。
首先需要说明的是,所述步骤S34-S37均是在时才进行,当时始终初始化。
在上述步骤S31中,所述初始化历元数的设定没有特殊限定,可根据需要自行设定。在上述步骤S33中,当时,初始化完成。
上述步骤S34进一步包括:
S341,根据位置、速度和加速度的运动学关系和前一历元状态量估算卫星当前历元的状态量;
S342,根据所述估算的当前历元的状态量计算预测伪距和预测载波相位,以及
S343,比较所述观测伪距和预测伪距以及比较所述观测载波相位和预测载波相位来确定所述保留伪距和保留载波相位。
在上述步骤S341中,加速度满足以下关系:
;
以上加速度表达式,在开普勒轨道定律基础上,考虑了地球非球形因素摄动中影响最大的J2项摄动,且为地心地固(ECEF)坐标系下的加速度表达式。其中地心地固坐标系的定义为:以地心作为坐标系原点,Z轴指向协议地球北极,X轴指向格林尼治参考子午面与地球赤道的一个交点,Y轴与X轴和Z轴构成右手直角坐标系。
其中,,,,为轨道动力学中与地球相关的参量,通常为万有引力常量与地球质量乘积,为J2项摄动参量,为地球半径,为地球自转角速度,为卫星到地心的距离。
此外,加速度积分即可以得到速度,速度积分就可以得到前一时刻到当前时刻移动的距离,利用该移动的距离再加上前一时刻的位置就可以得到当前位置。此外,可以用前一时刻的钟差和频差作为当前时刻的钟差和频差。从而获得所述预测的当前历元的状态量。此外,所述星载导航接收机的时间(标准时间)等于所述星载导航接收机本地时间减去钟差后的值。
在上述步骤S342中,所述预测伪距和预测载波相位可通过如下方法进行计算:
所述当前历元状态量包括用户位置、用户速度,已知卫星位置,卫星速度,是卫星编号,则预测伪距为(为前一历元的频差,t为历元间隔),预测载波相位为。
在上述步骤S343中,通过比较预测值(预测伪距和预测载波相位)与观测值(观测伪距和观测载波相位),将与预测值相差较大的观测值剔除。具体地,剔除过程包括以下步骤:
S3431,设定伪距门限范围以及载波相位门限范围;
S3432,如果观测伪距和预测伪距的差值超出所述伪距门限范围,则选取所述预测伪距作为所述保留伪距,否则选取所述观测伪距作为所述保留伪距;同时,如果所述观测的载波相位和预测载波相位的差值超出所述载波相位门限值,则选取所述预测载波相位作为所述保留载波相位,否则,选取所述观测的载波相位作为所述保留载波相位。
所述伪距门限范围和所述载波相位门限范围可根据需要来限定,优选地,所述伪距门限范围可以为3米至30米,更为优选地,所述伪距门限范围为10米。所述载波相位门限范围可以为0.1米至10米。优选地,所述载波相位门限范围为1米。
在上述步骤S35中,根据伪距定位原理构建所述解算方程,所述解算方程为:
其中,为保留伪距,为接收机与卫星的距离,为接收机钟差,为卫星钟差,为电离层延时,为对流层延时,为伪距测量误差。其中,时间参量、、、均以该时间参量乘以光速后的距离量表示。进一步地,可将所述解算方程简化为校正后的伪距观测方程:
,其中,ρc为校正后的伪距。
在上述步骤S35中,可采用非线性滤波法进行解算,得到当前历元的状态量。优选地,所述非线性滤波法可为Cubature卡尔曼滤波法(CKF)。
具体地,以表示以为均值、为方差的高斯分布,定义为当前历元观测量(伪距、卫星位置、卫星速度),为当前历元随机***噪声,服从,为当前历元随机观测噪声,服从。则离散非线性***可表示为:
其中,为***状态方程,为***观测方程。则CKF状态更新与观测更新方法如下:
对于状态更新,假设历元后验概率密度已知,对误差协方差进行Cholesky分解,有。则Cubature采样点可计算如下:
其中,为状态向量维数。通过状态方程将Cubature采样点传播到历元:
估计历元时的状态预测值和状态误差协方差预测值:
对于观测更新,对误差协方差进行Cholesky分解,有。则Cubature采样点可计算如下:
通过观测方程传播Cubature采样点:
估计历元时的观测量预测值及其自协方差和互协方差矩阵:
估计卡尔曼增益:
并最终得到历元时的状态估计值和状态误差协方差估计值:
。
在上述步骤S36中,可采用最小二乘残差自主完整性校验方法,即最小残差平方和法。其中残差向量的每一分量定义为星载导航接收机对相应卫星的校正后伪距值减去几何距离预测值与星载导航接收机钟差之和:
而残差平方和为标量,定义为残差向量与其转置的乘积:
的大小体现着各个测量值之间的一致性。由于定位解算中有四个独立有效的控制方程,因此在理论上服从自由度为的分布,其中为有效卫星的个数。根据分布概率密度函数,可以由设定的误警率确定的门限值,当超过门限值时则认为各测量值之间不一致,也就是其中某些测量值出现了错误,即表示校验不通过。
可见卫星数越多且卫星几何分布越好时,自主完整性校验算法的可靠性就越高。优选地,当同一时刻星载导航接收机可见卫星数大于等于5时,所述自主完整性校验才能够有效进行。
在上述步骤S37中,可设定自主完整性校验连续不通过最大历元数,当连续个历元发生自主完整性校验不通过时,或定位结果明显与常识或预期不符时,说明合理性校验不通过,即定位解算部分功能出现错误,此时重新返回步骤S31,开始初始化过程,以保证定位功能可以正常有效运行。
此外,若未出现合理性校验不通过的现象,可令,循环进入步骤S34,进行下一历元的解算方程构建和接收机状态解算工作。
请参阅图3,在上述步骤S4中,对各个历元的频差进行质量改善处理,具体包括以下步骤:
S41,设定初始化历元数;
S42,比较当前历元数和所述初始化历元数;
S43,当时,设定经过质量改善处理的当前历元的频差等于当前历元的频差,即:;以及
S44,当时,设定经过质量改善处理的当前历元的频差,其中,为滤波系数,且。
在上述步骤S41中,所述初始化历元数可根据需要自行设定,如可为10个历元、15个历元。
在上述步骤S43中,当前历元数小于所述初始化历元数时,即当已经历的历元数未达到给定的初始化历元数时,不对当前历元的频差进行质量改善处理。
在上述步骤S44中,当当前历元数大于等于所述初始化历元数时,才对当前历元的频差进行质量改善处理。通过加权平均的方式可使所述星载导航接收机频差的变化得到平滑,从而可减小随机噪声误差的影响以及提高所述星载导航接收机频差的精度。
请参阅图4,在上述步骤S5中,对所述钟差进行质量改善处理,具体包括以下步骤:
S51,设定初始化历元数以及钟差限定量;
S52,比较当前历元数和初始化历元数;
S53,当时,设定经过质量改善处理的当前历元的钟差等于当前历元的钟差,即:;以及
S54,当时,判断当前历元的与差值是否超出最大范围;
S55,若超过了所述最大范围,且成立,则令;
S56,若超过了所述最大范围,且成立,则令;以及
S57,若未超出所述最大范围,则令。
在上述步骤S51,所述钟差限定量用于在时,将当前历元的钟差与经过改善处理的钟差的差值限定在所述最大范围内。所述钟差限定量的取值范围为100纳秒至1微秒。
在上述步骤S53中,当已经历的历元数未达到所述设定的初始化历元数时,不对当前历元的频差进行质量改善处理。
在上述步骤S55-S57中,通过将当前历元的钟差与当前历元质量改善处理后的频差的差值限定在给定的最大范围内,对超出范围的值,以范围边界取代之。该方法保证钟差与频差的差值在给定范围内,从而提高钟差精度。
在上述步骤S7中,所述星载导航接收机可将所述调整后的秒脉冲通过秒脉冲广播总线广播给星上其它电子***。所述秒脉冲广播总线优选为多点低电压差分总线,该多点低电压差分总线可以使所述秒脉冲以低电压差分信号传输,从而可实现几百Mbps的传输速率,并且低电压差分信号的低压幅和低电流驱动输出可进一步降低噪声和减少功耗。此外,采用多点低电压差分总线传输还可减小电磁干扰。
在上述步骤S8中,所述星载导航接收机将所述位置、速度、时间转换为满足卫星总线设计要求的信息格式,并通过卫星数据总线广播给卫星上其它电子***。其中,所述信息格式包括精确到秒的时间、速度、位置、精度几何因子以及校验和。所述校验和是对精确到秒的时间、速度、位置、精度几何因子的校验结果。
在上述步骤S9中,,所述星载导航接收机生成和更新与时间相关的遥测量,并在星上计算机请求对应遥测量时给出应答的方法,其中星载导航接收机生成和更新的与时间相关的遥测量包括:升交点地方时,降交点地方时,以及所述{时间,位置,速度}组成的三元组。
升交点为星载导航接收机沿轨道从南向北运动时与赤道面的焦点,该点地方时即为升交点地方时。降交点为星载导航接收机沿轨道从北向南运动时与赤道面的焦点,该点地方时即为降交点地方时。所述时间、位置、速度三元组表示星载导航接收机经过解算和自主完整性、合理性校验后得到的接收机状态解算结果。
在上述步骤S9中,获得所述升交点地方时和降交点地方时的方法为:用地心地固(ECEF)坐标系表示解算结果中的位置,当当前历元最终解算结果的位置的Z轴解算结果与上一历元相比符号发生变化时,进行如下判断:如果Z轴速度为正值,则判断此时的时间为升交点地方时;否则,认为此时的时间为降交点地方时。
当当前历元最终解算结果的位置的Z轴解算结果与上一历元相比符号发生变化时,说明星载导航接收机在该历元通过了轨道与赤道面交点,此时若Z轴速度为正值,根据地心地固坐标系定义,可知星载导航接收机在该历元的运动方向为由南向北运动,从而此时的时间为升交点地方时;若Z轴速度为负值,根据地心地固坐标系定义,可知星载导航接收机在该历元的运动方向为由北向南运动,从而此时的时间为降交点地方时。
在上述步骤S9中,更新所述遥测量中所述{时间、位置、速度}三元组的方法为:
对之前缓存的{时间,位置,速度}三元组的遥测状态进行检测,只有当之前缓存的{时间,位置,速度}三元组已经被星上计算机遥测完成后,才用当前历元的{时间,位置,速度}三元组替代之前缓存的{时间,位置,速度}三元组,否则不进行更新。
请参阅图5,本发明实施例进一步提供一种授时型星载导航接收机100以实现所述星载导航接收机授时方法。
所述授时型星载导航接收机100包括天线102、射频前端模块104、捕获跟踪模块106、观测量获取模块108、解算模块110、钟差和频差质量改善模块112、遥测模块114以及秒脉冲模块116。
所述天线102用于接收卫星导航信号。所述射频前端模块104与所述天线102连接,接收并处理所述卫星导航信号,所述处理包括滤波、变频以及模数采样。所述变频可以为下变频。所述捕获跟踪模块106与所述射频前端模块104连接,接收经所述射频前端104处理后的卫星导航信号,并捕获和跟踪目标卫星的导航信号。所述观测量获取模块108接收所述跟踪后的目标卫星的导航信号以及秒脉冲,并在所述秒脉冲的上升沿采集该跟踪后的目标卫星的导航信号,获得所述观测信息。所述解算模块110接收所述观测信息,并根据所述卫星发射时刻和观测载波相位计算伪距,构建所述解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差。所述钟差和频差质量改善模块112用于所述解算模块110解算获得频差和钟差进行质量改善处理。所述遥测模块114用于生成和更新与所述时间相关的遥测量,并在星上计算机请求对应遥测量时给出应答。所述秒脉冲模块116用于产生所述秒脉冲,并根据所述质量改善处理后的频差和钟差调节所述秒脉冲的频率和相位,使所述秒脉冲的上升沿与标准时间整秒起始时刻同步。
本发明进一步提供一种星载导航应用***10,包括所述授时型星载导航接收机100、星上计算机200、测控通信模块300、被授时载荷400以及地面通信***500。
所述遥测模块114通过星上数据总线600与所述星上计算机200连接以使所述星上计算机200请求对应遥测量数据时给出应答。所述秒脉冲模块116通过秒脉冲数据总线700将调整后的所述秒脉冲传输到所述被授时载荷400。所述测控通信模块300用于在所述星上计算机200与地面通信***500之间实现指令与数据传输。
所述星上计算机200用于执行卫星上日常管理和记录功能,包括对所述授时型星载导航接收机100的日常管理和记录。所述被授时载荷400可为卫星上需要授时信息的专门应用载荷,如所述被授时载荷400可以为移动通信载荷。所述地面通信***500为地面人员对卫星进行远程控制和跟踪的输入输出界面。该地面通信***500可对所述授时型星载导航接收机100进行远程控制和跟踪。
本发明实施例提供的星载导航接收机授时方法与授时型星载导航接收机通过充分利用轨道动力学获得高精度的定位结果和原始钟差;并采用所述质量改善处理方法改善钟差和频差,有效地提高了授时精度以及星载导航接收机的***可靠性,并且可在轨计算卫星上需要的秒脉冲、升交点地方时、降交点地方时和标准时间等多种时间。尤其适用于有精密授时需求的低轨卫星。此外,本发明实施例提供的星载导航接收机授时方法与授时型星载导航接收机经过在轨飞行验证,采用第三方模拟环境和国家标准测试,结果显示:72小时授时精度方均根误差小于12ns。
另外,本领域技术人员还可以在本发明精神内做其它变化,当然,这些依据本发明精神所做的变化,都应包含在本发明所要求保护的范围之内。
Claims (10)
1.一种星载导航接收机授时方法,包括以下步骤:
S1,星载导航接收机产生秒脉冲;
S2,星载导航接收机在所述秒脉冲上升沿采集导航卫星信号中的观测信息,所述观测信息包括卫星发射时刻和观测载波相位;
S3,根据所述卫星发射时刻和观测载波相位计算伪距作为观测伪距,构建解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差;
S4,对所述频差进行质量改善处理,具体包括以下步骤:
S41,设定初始化历元数Nf;
S42,比较当前历元数k和所述初始化历元数Nf;
S43,当k<Nf时,设定经过质量改善处理的当前历元的频差等于当前历元的频差dfk,即:以及
S44,当k≥Nf时,设定经过质量改善处理的当前历元的频差其中,α为滤波系数,且0<α≤1;
S5,对所述钟差进行质量改善处理,包括以下步骤:
S51,设定初始化历元数Nt以及钟差限定量Δmax;
S52,比较当前历元数k和初始化历元数Nt;
S53,当k<Nt时,设定经过质量改善处理的当前历元的钟差等于当前历元的钟差dtk,即:以及
S54,当k≥Nt时,判断当前历元的dtk与差值是否超出最大范围[-Δmax,Δmax];
S55,若超过了所述最大范围,且成立,则令
S56,若超过了所述最大范围,且成立,则令以及;
S57,若未超出所述最大范围,则令
S6,根据改善后的频差和钟差调整所述秒脉冲的相位和频率,使所述秒脉冲上升沿与标准时间整秒起始时刻同步;
S7,星载导航接收机将所述调整后的秒脉冲广播给卫星上的其它电子***;
S8,在下一个所述调整后的秒脉冲到来之前,所述星载导航接收机将所述位置、速度、时间广播给卫星上的所述其它电子***;以及
S9,生成和更新与时间相关的遥测量,并在导航卫星请求对应遥测量时给出应答。
2.如权利要求1所述的星载导航接收机授时方法,其特征在于,所述步骤S3进一步包括:
S31,设定初始化历元数Np;
S32,比较当前历元数k和所述初始化历元数Np;
S33,当k<Np时,进行初始化,具体包括:使用迭代最小二乘法获得当前历元状态量Xk,所述当前历元状态量Xk包括位置、速度、钟差以及频差;
S34,当k≥Np时,进行异常观测量剔除以获得保留伪距和保留载波相位;
S35,根据所述保留伪距和保留载波相位构建解算方程,并进行解算,得到当前历元的状态量Xk;
S36,对所述当前历元的状态量Xk进行自主完整性校验,若校验通过,则采纳Xk作为当前历元的状态量,否则,采纳X'k作为当前历元的状态量;以及
S37,合理性校验:若发生连续自主完整性校验不通过或定位结果明显不合理,则回到步骤S31,重新开始初始化;若未出现合理性校验不通过的现象,可令k=k+1,循环进入步骤S34,进行下一历元的解算方程构建和接收机状态解算工作;
其中,所述步骤S34进一步包括以下步骤:
S341,根据位置、速度和加速度的运动学关系和前一历元状态量Xk-1估算卫星当前历元的状态量X'k;
S342,根据所述估算的当前历元的状态量X'k计算预测伪距和预测载波相位;以及
S343,比较所述观测伪距和预测伪距以及比较所述观测载波相位和预测载波相位,如果观测伪距和预测伪距的差值超出所述伪距门限范围,则选取所述预测伪距作为所述保留伪距,否则选取所述观测伪距作为所述保留伪距;同时,如果所述观测的载波相位和预测载波相位的差值超出所述载波相位门限值,则选取所述预测载波相位作为所述保留载波相位,否则,选取所述观测的载波相位作为所述保留载波相位。
3.如权利要求2所述的星载导航接收机授时方法,其特征在于,在上述步骤S341中,所述加速度满足以下关系:
其中,GM=3.986005×1014m3/s2,J2=1.082628×10-3,Re=6378137m,ωe=7.292115e-5rad/s,r为卫星到地心的距离,所述accx,k-1代表上一时刻地心地固坐标系中星载导航接收机X轴的加速度,所述accy,k-1代表上一时刻地心地固坐标系中星载导航接收机Y轴的加速度,所述accz,k-1代表上一时刻地心地固坐标系中星载导航接收机Z轴的加速度;所述xk-1为解算获得所述星载导航接收机的上一时刻位置的X轴分量,所述yk-1为解算获得所述星载导航接收机的上一时刻位置的Y轴分量,所述zk-1为为解算获得所述星载导航接收机的上一时刻位置的Z轴分量,所述vy,k-1为解算获得所述星载导航接收机的上一时刻速度的Y轴分量,所述vx,k-1为解算获得所述星载导航接收机的上一时刻速度的Z轴分量。
4.如权利要求2所述的星载导航接收机授时方法,其特征在于,采用Cubature卡尔曼滤波法进行解算。
5.如权利要求2所述的星载导航接收机授时方法,其特征在于,在上述步骤S36中,采用最小二乘残差平方和法实现自主完整性校验。
6.如权利要求1所述的星载导航接收机授时方法,其特征在于,所述星载导航接收机生成和更新的与时间相关的遥测量包括:升交点地方时、降交点地方时、以及时间、位置、速度组成的三元组。
7.如权利要求6所述的星载导航接收机授时方法,其特征在于,获得所述升交点地方时和降交点地方时的方法为:用地心地固坐标系表示解算结果中的位置,当当前历元最终解算结果的位置的Z轴解算结果与上一历元相比符号发生变化时,进行如下判断:如果Z轴速度为正值,则判断此时的时间为所述升交点地方时;否则,认为此时的时间为所述降交点地方时。
8.如权利要求6所述的星载导航接收机授时方法,其特征在于,更新所述遥测量中所述时间、位置、速度三元组的方法为:对之前缓存的时间,位置,速度三元组的遥测状态进行检测,只有当之前缓存的时间,位置,速度三元组已经被星上计算机遥测完成后,才用当前历元的时间,位置,速度三元组替代之前缓存的时间,位置,速度三元组,否则不进行更新。
9.一种用于实现如权利要求1-8项中任意一项所述的星载导航接收机授时方法的授时型星载导航接收机,其特征在于,包括:
天线,用于接收卫星导航信号;
射频前端模块,与所述天线连接,接收并处理所述卫星导航信号,所述处理包括滤波、变频以及模数采样;
捕获跟踪模块,与所述射频前端模块连接,接收经所述射频前端处理后的卫星导航信号,并捕获和跟踪目标卫星的导航信号;
观测量获取模块,接收所述跟踪后的目标卫星的导航信号以及所述秒脉冲,并在所述秒脉冲的上升沿采集该跟踪后的目标卫星的导航信号,获得所述观测信息;
解算模块,接收所述观测信息,并根据所述卫星发射时刻和观测载波相位计算所述观测伪距,构建所述解算方程,并通过解算获得所述星载导航接收机的位置、速度、时间、频差和钟差;
钟差和频差质量改善模块,用对所述解算模块解算获得的频差和钟差进行质量改善处理;
遥测模块,用于生成和更新与时间相关的所述遥测量,并在星上计算机请求对应遥测量时给出应答;以及
秒脉冲模块,用于产生所述秒脉冲,并根据所述质量改善处理后的频差和钟差调节所述秒脉冲的频率和相位,使所述秒脉冲的上升沿与标准时间整秒起始时刻同步。
10.一种星载导航应用***,包括:星上计算机、测控通信模块、被授时载荷以及地面通信***,其特征在于,进一步包括如权利要求9所述的授时型星载导航接收机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510060813.4A CN104678408B (zh) | 2015-02-05 | 2015-02-05 | 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510060813.4A CN104678408B (zh) | 2015-02-05 | 2015-02-05 | 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104678408A CN104678408A (zh) | 2015-06-03 |
CN104678408B true CN104678408B (zh) | 2017-02-22 |
Family
ID=53313732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510060813.4A Expired - Fee Related CN104678408B (zh) | 2015-02-05 | 2015-02-05 | 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104678408B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112147924A (zh) * | 2020-08-27 | 2020-12-29 | 航天东方红卫星有限公司 | 一种高精度程控任务管理*** |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106713820B (zh) * | 2016-11-29 | 2020-02-07 | 苏州科达特种视讯有限公司 | 卫星视频会议的广播频率切换方法及装置 |
CN109752735A (zh) * | 2017-11-08 | 2019-05-14 | 泰斗微电子科技有限公司 | 基于实时差分技术的时间同步方法和授时接收机*** |
CN108089431B (zh) * | 2018-01-24 | 2019-11-08 | 桂林电子科技大学 | 一种自适应卫星导航授时方法及装置 |
CN108345013B (zh) * | 2018-04-16 | 2023-09-01 | 南京天际易达通信技术有限公司 | 一种提高卫星导航信号接收灵敏度的方法 |
CN108828630A (zh) * | 2018-04-18 | 2018-11-16 | 交通运输部水运科学研究所 | 基带芯片性能评估方法及装置 |
CN108667507B (zh) * | 2018-05-31 | 2020-11-20 | 中国空间技术研究院 | 一种卫星智能化开放式综合电子*** |
CN109085616B (zh) * | 2018-08-30 | 2022-10-21 | 桂林电子科技大学 | 一种卫星授时方法、装置及存储介质 |
CN109633700B (zh) * | 2018-12-03 | 2022-11-22 | 天津津航计算技术研究所 | 一种多gps接收机授时精度测试方法 |
CN109521448B (zh) * | 2018-12-18 | 2020-01-17 | 清华大学 | 基于轨道根数预测的星载导航接收机定位授时方法和装置 |
CN109634092B (zh) * | 2019-01-16 | 2021-04-09 | 和芯星通科技(北京)有限公司 | 一种基于gnss接收机的授时方法及gnss接收机 |
CN110289903B (zh) * | 2019-06-26 | 2020-04-14 | 清华大学 | 一种基于链路测量的卫星网络自主时间同步方法和*** |
CN112147644B (zh) * | 2019-06-28 | 2024-02-23 | 清华大学 | 星地协同确定时空基准的方法、装置、设备及存储介质 |
CN110907954B (zh) * | 2019-10-22 | 2022-07-15 | 深圳航天东方红海特卫星有限公司 | 一种标校数据广播装置 |
CN111464397B (zh) * | 2020-04-08 | 2021-04-30 | 清华大学 | 一种双向距离和钟差测量的方法和*** |
CN111596327B (zh) * | 2020-04-16 | 2021-09-07 | 中国人民解放军军事科学院国防科技创新研究院 | 提高低轨卫星时钟性能的装置及方法 |
CN111398994B (zh) * | 2020-04-26 | 2020-10-30 | 清华大学 | 中轨通信卫星定位授时方法和装置 |
CN112817022B (zh) * | 2020-12-23 | 2022-04-12 | 浙江吉利控股集团有限公司 | 一种低轨卫星时频同步方法、***、电子设备和存储介质 |
US12007487B1 (en) * | 2020-12-31 | 2024-06-11 | Waymo Llc | Global positioning system time verification for autonomous vehicles |
CN113093250B (zh) * | 2021-06-04 | 2021-08-27 | 腾讯科技(深圳)有限公司 | 伪距观测数据修复方法、定位信息确定方法及装置 |
CN113325447B (zh) * | 2021-06-04 | 2023-07-28 | 华北电力大学 | 通过北斗gps信号测算时标的时间同一准对方法和装置 |
CN113254107B (zh) * | 2021-06-23 | 2021-10-22 | 星河动力(北京)空间科技有限公司 | 遥测数据的解算方法、装置、电子设备及计算机存储介质 |
CN113341686B (zh) * | 2021-07-02 | 2022-05-13 | 长沙学院 | 一种单星多历元校时方法、装置、计算机设备和存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176190A (zh) * | 2013-03-06 | 2013-06-26 | 西北工业大学 | 基于卫星导航和Kalman滤波的高精度授时***及方法 |
CN104049528A (zh) * | 2014-06-30 | 2014-09-17 | 北京北斗星通导航技术股份有限公司 | 一种北斗授时方法及卫星导航接收机 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013112428A1 (en) * | 2012-01-23 | 2013-08-01 | Fast Location.Net, Llc | Method and system for processing positioning signals based on predetermined message data segment |
-
2015
- 2015-02-05 CN CN201510060813.4A patent/CN104678408B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103176190A (zh) * | 2013-03-06 | 2013-06-26 | 西北工业大学 | 基于卫星导航和Kalman滤波的高精度授时***及方法 |
CN104049528A (zh) * | 2014-06-30 | 2014-09-17 | 北京北斗星通导航技术股份有限公司 | 一种北斗授时方法及卫星导航接收机 |
Non-Patent Citations (2)
Title |
---|
An Uplink Timing Synchronization Method for GEO Mobile SAT-LTE System;Jiajia Wang et al;《International Conference on Advanced Communications Technology,ICACT2014》;20140219;1045-1049 * |
卫星授时校频***研究;党晓圆;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20100115;C031-12 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112147924A (zh) * | 2020-08-27 | 2020-12-29 | 航天东方红卫星有限公司 | 一种高精度程控任务管理*** |
Also Published As
Publication number | Publication date |
---|---|
CN104678408A (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104678408B (zh) | 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用*** | |
CN109521448B (zh) | 基于轨道根数预测的星载导航接收机定位授时方法和装置 | |
US6424914B1 (en) | Fully-coupled vehicle positioning method and system thereof | |
EP3003861B1 (en) | Satellite orbital determination (od) using doppler and kepler orbital elements | |
US7817087B2 (en) | Method and apparatus for relative navigation using reflected GPS signals | |
Yang et al. | GPS-based onboard real-time orbit determination for LEO satellites using consider Kalman filter | |
US20190041527A1 (en) | Gps-based navigation system using a nonlinear discrete-time tracking filter | |
CN111965685B (zh) | 一种基于多普勒信息的低轨卫星/惯性组合导航定位方法 | |
US10908296B2 (en) | Method for calculating a speed of an aircraft, method for calculating a protection radius, positioning system and associated aircraft | |
Hwang et al. | GPS‐Based Orbit Determination for KOMPSAT‐5 Satellite | |
Chiaradia et al. | Onboard and Real‐Time Artificial Satellite Orbit Determination Using GPS | |
Wang et al. | GNSS-based orbit determination method and flight performance for geostationary satellites | |
Iiyama et al. | Precise Positioning and Timekeeping in a Lunar Orbit via Terrestrial GPS Time-Differenced Carrier-Phase Measurements | |
Iiyama et al. | Terrestrial GPS time-differenced carrier-phase positioning of lunar surface users | |
El-Mowafy et al. | The impact of orbital and clock errors on positioning from LEO constellations and proposed orbital solutions | |
Racelis et al. | Impact of cascading faults on mega-constellation-augmented GNSS PPP integrity | |
Bhamidipati et al. | Time transfer from GPS for designing a SmallSat-based lunar navigation satellite system | |
Liduan et al. | Enhanced GPS measurements simulation for space-oriented navigation system design | |
Defeng et al. | Reduced dynamic orbit determination using differenced phase in adjacent epochs for spaceborne dual-frequency GPS | |
CN115327587A (zh) | 基于gnss定位信息的低轨卫星轨道误差修正方法及*** | |
Willis et al. | DORIS satellite phase center determination and consequences on the derived scale of the Terrestrial Reference Frame | |
Xia et al. | Estimation and evaluation of COSMIC radio occultation excess phase using undifferenced measurements | |
Bhardwaj et al. | An assessment of ionospheric delay correction at L5 and S1 frequencies for NavIC Satellite System | |
Öz et al. | Performance assessment of a turn around ranging in communication satellite orbit determination | |
Bae | Near real-time precise orbit determination of low earth orbit satellites using an optimal GPS triple-differencing technique |
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: 20170222 Termination date: 20220205 |
|
CF01 | Termination of patent right due to non-payment of annual fee |