CN116338754B - 一种固定模糊度的非差非组合gnss时间频率传递方法 - Google Patents
一种固定模糊度的非差非组合gnss时间频率传递方法 Download PDFInfo
- Publication number
- CN116338754B CN116338754B CN202310599251.5A CN202310599251A CN116338754B CN 116338754 B CN116338754 B CN 116338754B CN 202310599251 A CN202310599251 A CN 202310599251A CN 116338754 B CN116338754 B CN 116338754B
- Authority
- CN
- China
- Prior art keywords
- model
- fixed
- time
- ionosphere
- ambiguity
- 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
- 238000000034 method Methods 0.000 title claims description 28
- 238000012546 transfer Methods 0.000 title claims description 23
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 239000005433 ionosphere Substances 0.000 claims description 54
- 230000005540 biological transmission Effects 0.000 claims description 25
- 238000007667 floating Methods 0.000 claims description 25
- 230000001934 delay Effects 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 4
- 239000005436 troposphere Substances 0.000 claims description 3
- 238000002716 delivery method Methods 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract description 3
- 230000006872 improvement Effects 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 101100456255 Arabidopsis thaliana AMC4 gene Proteins 0.000 description 1
- 101100512522 Arabidopsis thaliana AMC7 gene Proteins 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 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/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/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
-
- 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/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
- G01S19/072—Ionosphere corrections
-
- 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/254—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 Doppler shift of satellite signals
-
- 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/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/258—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 the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view
-
- 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/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种固定模糊度的非差非组合GNSS时间频率传递方法,所述固定模糊度的非差非组合GNSS时间频率传递方法包括以下步骤;步骤一、接收机获取相关数据信息,构建卫星钟固定模型;步骤二、将卫星钟固定模型通过S变换进行参数重整,得到满秩卫星钟固定模型;步骤三、对满秩卫星钟固定模型进行模糊度固定,得到精确的接收机钟差。本设计更快的对多频多系统相位模糊度进行固定,模糊度固定后,相位观测可以转换成高精度距离观测值,使解算结果精度更高。
Description
技术领域
本发明涉及一种时间频率传递技术的改进,属于时间频率传递领域,尤其涉及一种固定模糊度的非差非组合GNSS时间频率传递方法。
背景技术
精密的时间频率传递在相对论大地测量学、高分辨率射电天文学以及精密测量物理等研究领域起着至关重要的作用,同时对于金融、空间科学、交通运输和国防等关键基础设施建设也是必不可少的。当前时间传递技术主要有光纤、激光和GNSS等,GNSS因其操作简单、成本低而广泛应用于时间频率传递;目前较常用的GNSS时间频率传递方法为精密单点定位(Precise Point Positioning, PPP),该方法原理是将基线两端接收机分别做PPP处理得到其对应的本地精确时间,这两个时间作差即可知道两端时间的差异,进而实现时间频率传递,其优点是不需要基线两端共视同一组卫星,可在全球范围内实现时间频率传递,但是依赖精密卫星钟差和精密轨道等外部精密产品;传统PPP技术基于消电离层组合(IF),且通常采用双频观测值,但观测值的线性组合不仅会放大了观测噪声,同时在多频场景下存在不灵活的缺点;而且PPP模型中的相位模糊度与卫星相位偏差等其它参数组合在一起并以浮点数形式存在,这导致模糊度丧失了整数特性,使得高精度载波相位观测值的优势未能得到充分利用。即使新出现的非组合PPP解决了消电离层组合PPP放大观测噪声以及在多频场景下不灵活的缺陷,但其模糊度参数也仍为非整数。
申请号为CN202210844843.4,申请日为2022年7月18日的中国专利申请揭示了一种模糊度固定方法以及系统。采用伪距和载波观测值构建无电离层组合观测值并进行滤波得到无电离层组合模糊度;对非组合的伪距和载波观测值进行滤波得到非组合模糊度,并以无电离层组合模糊度作为虚拟模糊度观测值对非组合模糊度进行单差约束;对无电离层组合观测值滤波的浮点解进行单差模糊度固定并从固定的单差模糊度中获取单差电离层估计值;以单差电离层估计值作为虚拟电离层观测值对非组合观测值滤波的浮点解进行单差约束,然后对非组合观测值滤波的浮点解进行单差模糊度固定。本申请可以不需要额外辅助信息,也不需要采用多个观测值频点组合,实现精密单点定位的模糊度快速固定,但是对比文件依旧没有解决多频多系统相位模糊度固定的问题。
公开该背景技术部分的信息仅仅旨在增加对本专利申请的总体背景的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
本发明的目的是克服现有技术中存在多频多系统相位模糊度固定的问题,提供了对多频多系统相位模糊度进行固定的一种固定模糊度的非差非组合GNSS时间频率传递方法。
为实现以上目的,本发明的技术解决方案是:一种固定模糊度的非差非组合GNSS时间频率传递方法,所述固定模糊度的非差非组合GNSS时间频率传递方法包括以下步骤;
步骤一、接收机获取相关数据信息,构建卫星钟固定模型;
步骤二、将卫星钟固定模型通过S变换进行参数重整,得到满秩卫星钟固定模型;
步骤三、对满秩卫星钟固定模型进行模糊度固定,得到精确接收机钟差。
所述接收机获取相关数据信息包括接收机A和接收机B的GNSS码、相位观测值,精密卫星钟差数据、精密卫星轨道数据以及广播星历。
所述接收机A和接收机B具体为:假设现有一条基线AB,基线两端为接收机A和B,共同跟踪同一组卫星s=1,2,…,m,观测频率j=1,2,…,f,则卫星钟固定模型为:
。
所述卫星钟固定模型包括电离层固定模型、电离层加权模型与电离层浮点模型,电离层固定模型用于短基线0-10km的时间传递,电离层加权模型用于中基线10-100km的时间传递,电离层浮点模型用于大于100km的时间传递。
所述电离层固定模型中两个测站的电离层延迟相同,则两个测站只估计同一组电离层延迟;
所述电离层加权模型中两个测站的电离层延迟不同但是有相关性,即两个测站的单差电离层延迟接近0,可引入伪观测值对单差电离层延迟进行约束,提高模型强度;
所述电离层浮点模型认为两个测站的电离层延迟不同且没有相关性,需要独立估计两组电离层延迟;
所述电离层浮点模型如下:
;
基线两端测站的精确位置已知,同时卫星精密轨道由IGS提供,精确计算卫星到接收机的站星距,/>对流层延迟,IF卫星钟差/>由IGS提供;其中/>为天顶对流层湿延迟需要在模型中估计,/>为湿延迟投影函数;/>为倾斜干延迟。
所述电离层浮点模型中观测值如下:
。
所述电离层加权模型如下:
;
为伪观测值,设为0即可,其权重设为/>,其中/>和/>分别表示AB基线长度和卫星高度角;
其中,为站间单差的差分码偏差;
为站间单差电离层延迟。
所述电离层固定模型如下:
;
和/>都是接收机相对卫星参考时间的时间差,在矫正了/>后,最终得到 A和B的时间差异/>,接收机B的时间改正/>后即可与A时间保持时间频率同步。
所述为站间单差的消电离层组合接收机码偏差,其值在恒定环境条件下是不变的,获取方法如下:将接收机A和B时钟接入同一频率源,A和B组成零或短基线,那么有,解算后直接得到/>。
与现有技术相比,本发明的有益效果为:
1、本发明一种固定模糊度的非差非组合GNSS时间频率传递方法中,通过S变换对模糊度参数重整,消除了模糊度其与它参数之间的线性相关,使得整数模糊度参数独立可估,可实现模糊度固定,实现更快的初始化,更快的对多频多系统相位模糊度进行固定,模糊度固定后,相位观测可以转换成高精度距离观测值,使解算结果精度更高。因此,本设计可以对多频多系统相位模糊度进行固定,解算结果精度更高。
2、本发明一种固定模糊度的非差非组合GNSS时间频率传递方法中,通过电离层固定模型、电离层加权模型与电离层浮点模型,可实现短、中、长基线以及全球范围内的时间频率传递,时间频率传递性能更好,可实现更好的时间频率稳定性,同时没有放大观测噪声,且可灵活处理多频数据。因此,本设计时间频率传递性能更好,数据处理灵活。
附图说明
图1是本发明的ROAG-SFER时间差示意图。
图2是本发明中ROAG-SFER的改进Allan方差示意图。
图3是本发明中NIST-AMC4时间差示意图。
图4是本发明中NIST-AMC4的改进Allan方差示意图。
图5是本发明中PTBB-ROAG时间差示意图。
图6是本发明中PTBB-ROAG的改进Allan方差示意图。
具体实施方式
以下结合附图说明和具体实施方式对本发明作进一步详细的说明。
参见图1至图6,一种固定模糊度的非差非组合GNSS时间频率传递方法,所述固定模糊度的非差非组合GNSS时间频率传递方法包括以下步骤;
步骤一、接收机获取相关数据信息,构建卫星钟固定模型;
步骤二、将卫星钟固定模型通过S变换进行参数重整,得到满秩卫星钟固定模型;
步骤三、对满秩卫星钟固定模型进行模糊度固定,得到精确接收机钟差。
所述接收机获取相关数据信息包括接收机A和接收机B的GNSS码、相位观测值,精密卫星钟差数据、精密卫星轨道数据以及广播星历。
所述接收机A和接收机B具体为:假设现有一条基线AB,基线两端为接收机A和B,共同跟踪同一组卫星s=1,2,…,m,观测频率j=1,2,…,f ,则卫星钟固定模型为:
。
所述卫星钟固定模型包括电离层固定模型、电离层加权模型与电离层浮点模型,电离层固定模型用于短基线0-10km的时间传递,电离层加权模型用于中基线10-100km的时间传递,电离层浮点模型用于大于100km的时间传递。
所述电离层固定模型中两个测站的电离层延迟相同,则两个测站只估计同一组电离层延迟;
所述电离层加权模型中两个测站的电离层延迟不同但是有相关性,即两个测站的单差电离层延迟接近0,可引入伪观测值对单差电离层延迟进行约束,提高模型强度;
所述电离层浮点模型认为两个测站的电离层延迟不同且没有相关性,需要独立估计两组电离层延迟。
所述电离层浮点模型如下:
;
基线两端测站的精确位置已知,同时卫星精密轨道由IGS提供,精确计算卫星到接收机的站星距,/>对流层延迟,IF卫星钟差/>由IGS提供;其中/>为天顶对流层湿延迟需要在模型中估计,/>为湿延迟投影函数;/>为倾斜干延迟。
所述电离层浮点模型中观测值如下:
。
所述电离层加权模型如下:
;
为伪观测值,设为0即可,其权重设为/>,其中/>和/>分别表示AB基线长度和卫星高度角;
其中,为站间单差的差分码偏差;
为站间单差电离层延迟。
所述电离层固定模型如下:
。
和/>都是接收机相对卫星参考时间的时间差,在矫正了/>后,最终得到 A和B的时间差异/>,接收机B的时间改正/>后即可与A时间保持时间频率同步。
所述为站间单差的消电离层组合接收机码偏差,其值在恒定环境条件下是不变的,获取方法如下:将接收机A和B时钟接入同一频率源,A和B组成零或短基线,那么有,解算后直接得到/>。
实施例1:
一种固定模糊度的非差非组合GNSS时间频率传递方法,所述固定模糊度的非差非组合GNSS时间频率传递方法包括以下步骤;
步骤一、接收机获取相关数据信息,构建卫星钟固定模型;
步骤二、将卫星钟固定模型通过S变换进行参数重整,得到满秩卫星钟固定模型;
步骤三、对满秩卫星钟固定模型进行模糊度固定,得到精确接收机钟差。
实施例2:
实施例2与实施例1基本相同,其不同之处在于:
所述接收机获取相关数据信息包括接收机A和接收机B的GNSS码、相位观测值,精密卫星钟差数据、精密卫星轨道数据以及广播星历。
所述接收机A和接收机B具体为:假设现有一条基线AB,基线两端为接收机A和B,共同跟踪同一组卫星s=1,2,…,m,观测频率j=1,2,…,f ,则卫星钟固定模型为:
;
卫星钟固定模型中所使用的符号定义如下表所示
所述卫星钟固定模型包括电离层固定模型、电离层加权模型与电离层浮点模型,电离层固定模型用于短基线0-10km的时间传递,电离层加权模型用于中基线10-100km的时间传递,电离层浮点模型用于大于100km的时间传递。
所述电离层固定模型中两个测站的电离层延迟相同,则两个测站只估计同一组电离层延迟;
所述电离层加权模型中两个测站的电离层延迟不同但是有相关性,即两个测站的单差电离层延迟接近0,可引入伪观测值对单差电离层延迟进行约束,提高模型强度;
所述电离层浮点模型认为两个测站的电离层延迟不同且没有相关性,需要独立估计两组电离层延迟。
实施例3:
实施例3与实施例1基本相同,其不同之处在于:
一种固定模糊度的非差非组合GNSS时间频率传递方法,所述电离层浮点模型如下:
;
为了实现高精度的时间和频率传递,基线两端测站的精确位置已知,同时卫星精密轨道由IGS提供,精确计算卫星到接收机的站星距,/>对流层延迟,IF卫星钟差/>由IGS提供,其中/>为天顶对流层湿延迟需要在模型中估计,/>为湿延迟投影函数;/>为倾斜干延迟,可以通过经验模型来提前改正;
所述电离层浮点模型中观测值如下:
所有参数均独立可估,相位模糊度为双差整数形式,可实现模糊度固定,两个测站的电离层独立估计,没有考虑测站A和B电离层延迟的空间相关性。
实施例4:
实施例4与实施例1基本相同,其不同之处在于:
一种固定模糊度的非差非组合GNSS时间频率传递方法,接收机A和接收机B的电离层延迟有较强的空间相关性,可将站间单差电离层延迟作为伪观测值,构建电离层延迟约束模型,提高模型强度,引入单差电离层伪观测值后,所述电离层加权模型如下:
;
为伪观测值,设为0即可,其权重设为/>,其中/>和/>分别表示AB基线长度和卫星高度角;
其中,为站间单差的差分码偏差;
为站间单差电离层延迟。
实施例5:
实施例5与实施例1基本相同,其不同之处在于:
一种固定模糊度的非差非组合GNSS时间频率传递方法,当时间传递的距离(基线长度)不足十公里,可认为站间单差电离层延迟为零,所述电离层固定模型如下:
;
和/>都是接收机相对卫星参考时间的时间差,在矫正了/>后,最终得到 A和B的时间差异/>,接收机B的时间改正/>后即可与A时间保持时间频率同步。
所述为站间单差的消电离层组合接收机码偏差,其值在恒定环境条件下是不变的,获取方法如下:将接收机A和B时钟接入同一频率源,A和B组成零或短基线,那么有,解算后直接得到/>;
得到基线两端的接收机钟差后即可实现时间传递(剩下步骤与前面相同)。但需要注意的是,不管是采用电离层固定模型、电离层加权模型或电离层浮点模型,基线两端需要有共视卫星才能实现模糊度固定,当基线长度到达数千公里或更长时,基线两端已没有足够或没有共视卫星,此时需要在接收机A和B附近分别架设辅助站C和D,A和C、B和D可分别组成零、短或中基线,这时采用电离层加权模型或电离层浮点模型分别解算基线AC和基线BD即可解算出A和B的钟差,进而实现时间频率传递(剩下步骤也与之前相同),但需要注意的是,辅助站C和D只是为了辅助模糊度固定,所以不需要接入原子钟,且这种方法不再受制于基线长度,可实现超长基线或全球范围内的时间频率传递。
实施例6:
实施例6与实施例1基本相同,其不同之处在于:
一种固定模糊度的非差非组合GNSS时间频率传递方法,采用一些实验室公开的短、中、长基线接收机观测数据,这些接收机均接入了原子钟,具体信息如下表所示:
上表中为SCF时间传递实验中涉及的GNSS和原子钟等信息;
短基线:测站ROAG和SFER组成短基线ROAG-SFER,长度124m,采用SCF Ion-Fixed模型和两个站分别采用传统PPP技术进行处理,得到这两个测站的钟差,并计算这两个测站的时间差,结果见图1,对两种方案得到的时间差进行Allan方差估计,结果见图2;图1和2分别给出了ROAG-SFER的时间差以及相应的改进Allan方差;本实验以及后续实验中引入了时间差的历元差作为评估指标,该指标通常用来表征原子钟单位时间内的变化,进而反映时间传递的性能;SCF Ion-Fixe模型相比于传统PPP模型在历元差的标准差上有31.7%的提升(分别为5.65 ps和8.27 ps);此外,SCF模型在频率稳定度方面比传统PPP模型也有明显的提升;以3840 s的频率稳定度为例,两个模型的频率稳定度分别为和/>(图2);此外,SCF模型和传统PPP模型的天稳分别为/>和/>(图2);上述结果表明了SCF模型相比于传统PPP模型的优势;
中长基线实验中:NIST和AMC4组成中基线NIST-AMC4(146.8 km), 分布采用SCFIon-Weighted模型和两个站分别采用传统PPP技术进行处理,得到这两个测站的钟差,并计算这两个测站的时间差,结果见图3。对两种方案得到的时间差进行Allan方差估计,结果见图4;图3为两种方案估计的NIST-AMC4时间差;SCF Ion-Weighted模型和传统PPP历元差的标准差分别为5.61 ps和6.71 ps,有16.4%的提升;图4给出了NIST-AMC4时间差的改进Allan方差,结论为: SCF Ion-Weighted在几十至几千秒的频率稳定度上相比于传统PPP有15%至30%的提升;SCF电离层加权模型比传统PPP在天稳定性上有7.9%的提升;这说明SCF模型在短期稳定度和长期稳定度均比PPP方法有明显的改善;
长基线实验中:PTBB和ROAG组成长基线PTBB-ROAG(2182.3 km), 此时时间链路两端没有足够的共视卫星,无法实现模糊度固定,这时在PTBB和ROAG站附近分别架设辅助站C和D,并采用SCF Ion-Fixed模型解算短基线PTBB-C和ROAG-D,得到PTBB和ROAG的接收机钟差,在做差即可得到时间差,与传统PPP处理方式进行比较,时间差估计结果及Allan方差分析结果见图5-6;图5和图6结果表明,对于几十秒到几千秒的频率稳定度,SCF模型比传统PPP模型有明显的提升。然而,对于4000秒至30000秒的频率稳定度,SCF模型与传统PPP模型的结果相当。此外,SCF模型的天稳定性同样比传统PPP模型的结果好。
以上所述仅为本发明的较佳实施方式,本发明的保护范围并不以上述实施方式为限,但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化,皆应纳入权利要求书中记载的保护范围内。
Claims (7)
1.一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述固定模糊度的非差非组合GNSS时间频率传递方法包括以下步骤:
步骤一、接收机获取相关数据信息,构建卫星钟固定模型;
步骤二、将卫星钟固定模型通过S变换进行参数重整,得到满秩卫星钟固定模型;
步骤三、对满秩卫星钟固定模型进行模糊度固定,得到精确接收机钟差;
所述接收机获取相关数据信息包括接收机A和接收机B的GNSS码、相位观测值,精密卫星钟差数据、精密卫星轨道数据以及广播星历;
所述接收机A和接收机B具体为:假设现有一条基线AB,基线两端为接收机A和B,共同跟踪同一组卫星s=1,2,…,m,观测频率j=1,2,…,f,则卫星钟固定模型为:
所述卫星钟固定模型包括电离层固定模型、电离层加权模型与电离层浮点模型,电离层固定模型用于短基线0-10km的时间传递,电离层加权模型用于中基线10-100km的时间传递,电离层浮点模型用于大于100km的时间传递;
当基线长度大于1000公里时,在接收机A和B附近分别架设辅助站C和D,A和C、B和D可分别组成零、短或中基线,此时采用电离层加权模型或电离层浮点模型分别解算基线AC和基线BD即可解算出A和B的钟差,进而实现时间频率传递。
2.根据权利要求1所述的一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述电离层固定模型中两个测站的电离层延迟相同,则两个测站只估计同一组电离层延迟;
所述电离层加权模型中两个测站的电离层延迟不同但是有相关性,即两个测站的单差电离层延迟接近0,可引入伪观测值对单差电离层延迟进行约束,提高模型强度;
所述电离层浮点模型认为两个测站的电离层延迟不同且没有相关性,需要独立估计两组电离层延迟。
3.根据权利要求1所述的一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述电离层浮点模型如下:
基线两端测站的精确位置已知,同时卫星精密轨道由IGS提供,精确计算卫星到接收机的站星距对流层延迟,IF卫星钟差/>由IGS提供;其中τr为天顶对流层湿延迟需要在模型中估计,/>为湿延迟投影函数;/>为倾斜干延迟。
4.根据权利要求3所述的一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述电离层浮点模型中观测值如下:
5.根据权利要求1所述的一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述电离层加权模型如下:
为伪观测值,设为0即可,其权重设为/>其中lAB和Es分别表示AB基线长度和卫星高度角;
其中,为站间单差的差分码偏差;
为站间单差电离层延迟。
6.根据权利要求1所述的一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述电离层固定模型如下:
dtA和dtB都是接收机相对卫星参考时间的时间差,在矫正了dAB,IF后,最终得到A和B的时间差异接收机B的时间改正dtAB后即可与A时间保持时间频率同步。
7.根据权利要求6所述的一种固定模糊度的非差非组合GNSS时间频率传递方法,其特征在于:所述dAB,IF为站间单差的消电离层组合接收机码偏差,其值在恒定环境条件下是不变的,获取方法如下:将接收机A和B时钟接入同一频率源,A和B组成零或短基线,那么有dtA=dtB,解算后直接得到
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310599251.5A CN116338754B (zh) | 2023-05-25 | 2023-05-25 | 一种固定模糊度的非差非组合gnss时间频率传递方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310599251.5A CN116338754B (zh) | 2023-05-25 | 2023-05-25 | 一种固定模糊度的非差非组合gnss时间频率传递方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116338754A CN116338754A (zh) | 2023-06-27 |
CN116338754B true CN116338754B (zh) | 2023-09-12 |
Family
ID=86880705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310599251.5A Active CN116338754B (zh) | 2023-05-25 | 2023-05-25 | 一种固定模糊度的非差非组合gnss时间频率传递方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116338754B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116893436B (zh) * | 2023-09-11 | 2023-12-05 | 中国科学院精密测量科学与技术创新研究院 | 一种混合频率的rtk定位方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025354A (zh) * | 2019-12-20 | 2020-04-17 | 东南大学 | 基于单差分电离层加权模型的中长基线rtk定位方法 |
CN111190203A (zh) * | 2020-01-14 | 2020-05-22 | 中国科学院测量与地球物理研究所 | 基于非差观测模型和消秩亏的ppp-rtk定位方法 |
CN115267843A (zh) * | 2022-06-14 | 2022-11-01 | 中国科学院精密测量科学与技术创新研究院 | 多频多模gnss高频精密卫星钟差的实时非差估计方法 |
CN115993614A (zh) * | 2022-12-07 | 2023-04-21 | 中国科学院国家授时中心 | 一种削弱北斗系统内偏差的时间传递方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7868820B2 (en) * | 2005-09-09 | 2011-01-11 | Trimble Navigation Limited | Ionosphere modeling apparatus and methods |
DE102012202095A1 (de) * | 2011-02-14 | 2012-08-16 | Trimble Navigation Ltd. | GNSS-Signalverarbeitung mit Ionosphärenmodell für synthetische Referenzdaten |
NL2009695C2 (en) * | 2012-10-25 | 2014-05-06 | Fugro N V | Ppp-rtk method and system for gnss signal based position determination. |
-
2023
- 2023-05-25 CN CN202310599251.5A patent/CN116338754B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025354A (zh) * | 2019-12-20 | 2020-04-17 | 东南大学 | 基于单差分电离层加权模型的中长基线rtk定位方法 |
CN111190203A (zh) * | 2020-01-14 | 2020-05-22 | 中国科学院测量与地球物理研究所 | 基于非差观测模型和消秩亏的ppp-rtk定位方法 |
CN115267843A (zh) * | 2022-06-14 | 2022-11-01 | 中国科学院精密测量科学与技术创新研究院 | 多频多模gnss高频精密卫星钟差的实时非差估计方法 |
CN115993614A (zh) * | 2022-12-07 | 2023-04-21 | 中国科学院国家授时中心 | 一种削弱北斗系统内偏差的时间传递方法 |
Non-Patent Citations (1)
Title |
---|
"Consistency of seven different GNSS global ionosphericmapping techniques during one solar cycle";David Roma-Dollase等;《Journal of Geodesy》;1-16 * |
Also Published As
Publication number | Publication date |
---|---|
CN116338754A (zh) | 2023-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
El-Mowafy et al. | On biases in precise point positioning with multi-constellation and multi-frequency GNSS data | |
CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
Zhang et al. | Performance analysis of triple-frequency ambiguity resolution with BeiDou observations | |
Su et al. | Triple-frequency carrier phase precise time and frequency transfer models for BDS-3 | |
Zhao et al. | Three-carrier ambiguity resolution using the modified TCAR method | |
CN109581455B (zh) | 一种bds和gps融合的三频宽巷紧组合定位方法 | |
CN116338754B (zh) | 一种固定模糊度的非差非组合gnss时间频率传递方法 | |
CN109001781A (zh) | 一种顾及电离层约束的bds三频模糊度解算方法 | |
CN114966760B (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
CN116243591B (zh) | 融合UTC(k)和北斗广播星历的亚纳秒级授时方法 | |
CN113848572A (zh) | 一种基于大气误差增强的多频ppp序贯单历元定位方法 | |
CN115933356B (zh) | 一种虚拟原子钟的高精度时间同步系统和方法 | |
CN111352137A (zh) | 一种顾及广播星历误差的多模gnss异步rtk定位方法 | |
CN114879239B (zh) | 一种增强瞬时ppp固定解的区域三频整数钟差估计方法 | |
CN114935770B (zh) | 一种多历元加快精密单点定位收敛速度的方法及装置 | |
CN115079236A (zh) | 通过低轨增强缩短广域非差非组合ppp-rtk定位收敛时间的方法 | |
CN110988932B (zh) | 一种提高星载gps接收机实时钟差解算精度的方法 | |
CN116466377A (zh) | 一种消除ifb的glonass ppp-rtk方法 | |
CN115308781B (zh) | 基于bdgim辅助的相位平滑伪距高精度时间传递方法 | |
CN116359968A (zh) | 一种联合北斗二号和北斗三号的三频差分定位方法 | |
CN115856945A (zh) | 适用于gnss多频点的精密单点授时方法、装置及介质 | |
Ren et al. | Common-view time transfer with BeiDou new navigation signals | |
CN116893436B (zh) | 一种混合频率的rtk定位方法 | |
US20230077069A1 (en) | System and method for timing synchronization and tracking of satellite signal receivers |
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 |