CN104678408B - 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用系统 - Google Patents

星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用系统 Download PDF

Info

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
based navigation
pulse per
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
Application number
CN201510060813.4A
Other languages
English (en)
Other versions
CN104678408A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201510060813.4A priority Critical patent/CN104678408B/zh
Publication of CN104678408A publication Critical patent/CN104678408A/zh
Application granted granted Critical
Publication of CN104678408B publication Critical patent/CN104678408B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/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中，所述加速度满足以下关系：
acc x , k - 1 = - G M r 3 x k - 1 ( 1 + 3 2 J 2 ( R e r ) 2 ( 1 - 5 ( R e r ) 2 ) ) + ω e 2 x k - 1 + 2 ω e v y , k - 1 acc y , k - 1 = - G M r 3 y k - 1 ( 1 + 3 2 J 2 ( R e r ) 2 ( 1 - 5 ( R e r ) 2 ) ) + ω e 2 y k - 1 - 2 ω e v x , k - 1 acc z , k - 1 = - G M r 3 z k - 1 ( 1 + 3 2 J 2 ( R e r ) 2 ( 3 - 5 ( R e r ) 2 ) ) ;
其中，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所述的授时型星载导航接收机。
CN201510060813.4A 2015-02-05 2015-02-05 星载导航接收机授时方法和授时型星载导航接收机以及星载导航应用系统 Expired - Fee Related CN104678408B (zh)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147924A (zh) * 2020-08-27 2020-12-29 航天东方红卫星有限公司 一种高精度程控任务管理系统

Families Citing this family (22)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN101743453B (zh) 任务后高精确度定位和定向系统
CN101609140B (zh) 一种兼容导航接收机定位系统及其定位方法
US8082099B2 (en) Aircraft navigation using the global positioning system and an attitude and heading reference system
CN102033236B (zh) 一种卫星导航位置速度联合估计方法
CN103487820B (zh) 一种车载捷联/卫星紧组合无缝导航方法
US20130307724A1 (en) Apparatus And Method For Generating GPS Time
CN105607093A (zh) 一种组合导航系统及获取导航坐标的方法
Bock et al. GPS single-frequency orbit determination for low Earth orbiting satellites
CN102679985A (zh) 一种应用星间跟踪的航天器星座分散化自主导航方法
Hwang et al. GPS‐Based Orbit Determination for KOMPSAT‐5 Satellite
WO2020165604A1 (en) Techniques for determining geolocations
CN115356754A (zh) 一种基于gnss和低轨卫星的组合导航定位方法
CN103543454A (zh) 一种嵌入在移动通讯网中的卫星定轨系统
Iiyama et al. Precise Positioning and Timekeeping in a Lunar Orbit via Terrestrial GPS Time-Differenced Carrier-Phase Measurements
CN108205151B (zh) 一种低成本gps单天线姿态测量方法
Capuano et al. GNSS/INS/Star tracker integrated navigation system for Earth-Moon transfer orbit
CN110988932B (zh) 一种提高星载gps接收机实时钟差解算精度的方法
CN115327587A (zh) 基于gnss定位信息的低轨卫星轨道误差修正方法及系统
Liduan et al. Enhanced GPS measurements simulation for space-oriented navigation system design
Tenny et al. Robust Navigation for Urban Air Mobility via Tight Coupling of GNSS with Terrestrial Radionavigation and Inertial Sensing
Liu et al. Measurement-domain cooperative navigation for multi-UAV systems augmented by relative positions
Meng et al. One-way deep space navigation with radiometric and inertial data fusion

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