CN101833079B - 全球定位***原始测量值实时转为精确均匀采样的方法 - Google Patents
全球定位***原始测量值实时转为精确均匀采样的方法 Download PDFInfo
- Publication number
- CN101833079B CN101833079B CN2009100797959A CN200910079795A CN101833079B CN 101833079 B CN101833079 B CN 101833079B CN 2009100797959 A CN2009100797959 A CN 2009100797959A CN 200910079795 A CN200910079795 A CN 200910079795A CN 101833079 B CN101833079 B CN 101833079B
- Authority
- CN
- China
- Prior art keywords
- time
- gps
- measurement value
- receiver
- real time
- 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
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种全球定位***原始测量值实时转为精确均匀采样的方法,涉及全球定位技术,利用全球定位***(GPS,Global PositioningSystem)接收机实时计算的钟差(接收机时间和GPS时间之间的时间差),以及实时测量的多普勒值对接收机实时采集的伪距测量和载波相位测量值分别进行修正,从而可以得到和GPS时间精确同步、均匀采样的伪距测量和载波相位测量值,使全球定位***接收机实时采集的原始测量值转换成和GPS时间***精确同步且均匀采样的测量值。本发明巧妙应用了接收机钟差和多普勒值与GPS原始测量值之间的内在关系,实现了实时获取精确、均匀采样的伪距和载波相位等原始测量值,大大方便了GPS数据的实时定位、定向等处理工作。
Description
技术领域
本发明涉及全球定位技术,是一种全球定位***原始测量值实时转为精确均匀采样的方法。
背景技术
随着全球定位***(GPS,Global Positioning System)的技术和应用的深入,利用GPS进行实时定位和定向成为军事、科研、生活中的重要手段。由于单台GPS接收机定位的典型精度只能达到5~10米(中误差),因此实际高精度(亚米级、厘米级)应用中都需要两台或多台GPS接收机进行同步观测并通过差分处理来实现。
在应用两台或多台GPS接收机进行同步观测和数据处理时,都会要求两台或多台GPS接收机的原始测量值是精确同步且是均匀采样的测量值,这样才能应用差分方法进行处理。常规的方法是在差分处理前,将其中的一台GPS接收机的测量值首先转化成和GPS时间***严格同步的测量值,并以该接收机对应的采样时间作为参考,按照它的采样时间和采样间隔对其余GPS接收机所采集的原始测量值进行重新采样处理,重新采样处理的方法一般是数值计算中用到的数学内插方法。它的缺点是处理时比较复杂、计算量大,尤其是同步观测的GPS接收机数量多时,重新采样的测量值也不均匀。此外如果是对GPS进行实时处理,则在对GPS测量值重新采样时,需要采用当前采样时刻前后一段时间内的多个测量值,实时性差。
发明内容
本发明的目的是提供一种全球定位***原始测量值实时转为精确均匀采样的方法,以解决GPS接收机采集的原始测量值和GPS时间***不精确同步且采样不均匀的问题,从而得到采样均匀且和GPS时间同步的测量值。
为达到上述目的,本发明的技术解决方案是:
一种全球定位***原始测量值实时转为精确均匀采样的方法,利用至少四颗全球定位***卫星;其包括:
步骤1:确定全球定位***接收机钟差;
步骤2:从全球定位***接收机中实时提取接收机所接收的多普勒测量值;
步骤3:根据全球定位***接收机的钟差和多普勒测量值对接收机实时采集的原始测量值进行实时修正,得到和全球定位***时间***严格同步且采样间隔均匀的实时测量值。
所述的精确均匀采样的方法,其所述步骤1,是从全球定位***接收机中实时提取接收机时间***和全球定位***时间***的时间差:接收机钟差,或者根据全球定位***接收机所接收的伪距测量值实时计算出接收机钟差。
所述的精确均匀采样的方法,其特征在于,所述步骤3的原始测量值为伪距测量值和载波相位测量值,分别按不同的修正方法修正:
A)实时测量的任一颗卫星的伪距测量值修正公式为:
P′i=Pi-dT×Di×λ
其中,dT为钟差;Pi为伪距测量值,单位为m;多普勒值Di,单位Hz;λ为对应的全球定位***无线电波波长;
B)实时测量的任一颗卫星的载波相位测量值修正公式为:
C′i=Ci-dT×Di
其中,dT为钟差;Ci为载波相位测量值,单位为周,多普勒值Di,单位Hz。
本发明方法充分考虑了GPS接收机采样时间和GPS时间***之间的内在关系,应用了多普勒测量值和伪距、载波相位测量值之间的物理关系,根据接收机钟差和多普勒测量值对实时采集的伪距和载波相位测量值进行修正,修正后的测量值和GPS时间***精确同步,而且是均匀的采样。
本发明方法大大方便多台GPS测量值的实时差分数据,保证各GPS接收机的测量值是和GPS时间***精确同步、均匀采样的特点,省去了内插和重新采样。
附图说明
图1是用本发明方法对GPS接收机原始测量值进行实时修正的示意框图。
具体实施方式
一种全球定位***(GPS)原始测量值实时转为精确均匀采样的方法,包括:
步骤1:从GPS接收机中实时提取接收机时间***和GPS时间***的时间差:接收机钟差,或者根据接收机所接收的伪距测量值实时计算出接收机的钟差;
步骤2:从GPS接收机中实时提取接收机所接收的多普勒测量值;
步骤3:根据接收机的钟差和多普勒测量值对GPS接收机实时采集的伪距测量值和载波相位测量值按照不同的修正方法进行实时修正,就得到了和GPS时间***严格同步且采样间隔均匀的伪距测量值和载波相位测测量值。
本发明全球定位***(GPS)原始测量值实时转为精确均匀采样的方法,是对不同GPS接收机分别进行各自的修正,这样修正的各台GPS接收机的伪距和载波相位测量值都和GPS时间***严格同步,而且彼此之间的采样间隔也完全一致。
下面将结合附图1对本发明方法加以详细说明,应指出的是,所描述的实施例仅旨在便于对本发明的理解,而对本发明不起任何限定作用。
GPS***采用的GPS时间是一种高精度的时间***,一般精度可达到几十纳秒,而GPS接收机内部一般采用石英钟晶振来维持自己的时间***,其精度一般只能达到几百微妙甚至亚毫秒级的精度。
GPS接收机采集的原始数据是按照接收机自身的时间***进行数据的实时采样。因此GPS接收机所提供的实时测量值是在GPS接收机时间***中名义上按照固定的接收机时间间隔进行采样得到的。
假设,Ti r表示GPS接收机采集伪距和载波相位对应的接收机时间,此时刻接收机采集的原始数据分别为:伪距测量值Pi(单位为m)、载波相位测量值Ci(单位为周)。
Ti G表示接收机时间对应真正的GPS时间。
则 为接收机时间和GPS时间之间的时钟差,简称钟差,钟差dT是随着不同的采样时刻而不断变化的参数。
根据实时测量的4颗以上(含4颗)GPS卫星和接收机之间的伪距值,就可以实时计算接收机的钟差dT,而且一般接收机可以实时给出接收机的钟差参数。
因此GPS接收机实时采集的原始测量值(伪距和载波相位)所对应的GPS时刻和修正后实际固定采样间隔的GPS时刻存在的时间差就等于接收机的钟差dT。为了得到固定采样间隔的GPS时刻所对应的原始测量值,需要对钟差dT时间内引起的伪距测量值和载波相位测量值进行实时修正补偿。
GPS接收机在实时采集伪距测量值和载波相位测量值的同时,还可以采集到每一颗卫星和接收机之间的多普勒值Di(单位Hz),根据物理上的多普勒定义可以知道,多普勒显示了GPS接收机和卫星之间的距离变化率,因此可以根据多普勒数值来修正钟差dT时间内伪距测量值和载波相位测量值的变化量,设修正后的伪距为P′i、载波相位为C′i,且对应的GPS无线电波波长为λ,则有:
修正后每颗卫星的伪距测量值:P′i=Pi-dT×Di×λ
修正后每颗卫星的载波相位测量值:C′i=Ci-dT×Di
对每个采样时刻的伪距和载波相位测量值都按照提取的钟差和多普勒值Di进行实时修正,这样修正的GPS测量值都是对应于实际的GPS时间,而且是按照固定的采样间隔进行修正的,因此它们是均匀采样的测量值。
举例说明如下:
GPS实时测量的数据如下表所示:
接收机采样时刻对应的接收机时间(秒) | 伪距测量值(米) | 载波相位测量值(周) | 多普勒值(Hz) | 接收机钟差(秒) |
51.0000 | 24023042.804 | 126241967.688 | 1028.502 | 0.00015 |
52.0000 | 24022847.338 | 126240939.648 | 1027.506 | 0.00023 |
53.0000 | 24022651.364 | 126239912.581 | 1026.671 | 0.00031 |
实际对应的GPS时间和数据如下表所示(测量值没有变化):
接收机采样时刻对应的GPS时间(秒) | 伪距测量值(米) | 载波相位测量值(周) | 多普勒值(Hz) | 接收机钟差(秒) |
50.99985000 | 24023042.804 | 126241967.688 | 1028.502 | 0.00015 |
51.99977000 | 24022847.338 | 126240939.648 | 1027.506 | 0.00023 |
52.99969000 | 24022651.364 | 126239912.581 | 1026.671 | 0.00031 |
实时修正后和GPS时间严格同步且均匀采样的测量值如下表所示:
接收机采样时刻对应的GPS时间(秒) | 伪距测量值(米) | 载波相位测量值(周) |
51.000000000 | 24023042.775 | 126241967.534 |
52.000000000 | 24022847.102 | 126240939.412 |
53.000000000 | 24022651.303 | 126239912.263 |
Claims (3)
1.一种全球定位***原始测量值实时转为精确均匀采样的方法,利用至少四颗全球定位***卫星;其特征在于,包括:
步骤1:确定全球定位***接收机钟差;
步骤2:从全球定位***接收机中实时提取接收机所接收的多普勒测量值;
步骤3:根据全球定位***接收机的钟差和多普勒测量值对接收机实时采集的原始测量值进行实时修正,得到和全球定位***时间***严格同步且采样间隔均匀的实时测量值。
2.如权利要求1所述的全球定位***原始测量值实时转为精确均匀采样的方法,其特征在于,步骤1中所述全球定位***接收机钟差,是从全球定位***接收机中实时提取的接收机时间***与全球定位***时间***之间的时间差,或者是根据全球定位***接收机所接收的伪距测量值实时计算出的接收机钟差。
3.如权利要求1所述的全球定位***原始测量值实时转为精确均匀采样的方法,其特征在于,所述步骤3的原始测量值为伪距测量值和载波相位测量值,分别按不同的修正方法修正:
A)实时测量的任一颗卫星的伪距测量值修正公式为:
P′i=Pi-dT×Di×λ
其中,dT为钟差;Pi为伪距测量值,单位为m;多普勒值Di,单位Hz;λ为对应的全球定位***无线电波波长;
B)实时测量的任一颗卫星的载波相位测量值修正公式为:
C′i=Ci-dT×Di
其中,dT为钟差;Ci为载波相位测量值,单位为周,多普勒值Di,单位Hz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100797959A CN101833079B (zh) | 2009-03-11 | 2009-03-11 | 全球定位***原始测量值实时转为精确均匀采样的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100797959A CN101833079B (zh) | 2009-03-11 | 2009-03-11 | 全球定位***原始测量值实时转为精确均匀采样的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101833079A CN101833079A (zh) | 2010-09-15 |
CN101833079B true CN101833079B (zh) | 2012-08-22 |
Family
ID=42717219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100797959A Expired - Fee Related CN101833079B (zh) | 2009-03-11 | 2009-03-11 | 全球定位***原始测量值实时转为精确均匀采样的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101833079B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102098119B (zh) * | 2011-01-27 | 2013-09-25 | 电子科技大学 | 一种接收机钟差预测方法 |
CN102890279B (zh) * | 2011-07-20 | 2014-02-12 | 中国科学院电子学研究所 | 一种利用全球定位***进行纳秒级精确时间同步的方法 |
DE102014224911A1 (de) * | 2014-12-04 | 2016-06-09 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Bestimmen statistischer Eigenschaften von Rohmesswerten |
US10012725B2 (en) * | 2014-12-19 | 2018-07-03 | Qualcomm Incorporated | Systems, methods, and apparatus for living object protection having extended functionality in wireless power transfer applications |
CN110568466B (zh) * | 2019-08-19 | 2021-09-14 | 北京自动化控制设备研究所 | 基准站观测量推算方法及载波相位差分定位信息获取方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7133480B2 (en) * | 2001-03-09 | 2006-11-07 | Leica Geosystems Inc. | Method and apparatus for processing digitally sampled signals at a resolution finer than that of a sampling clock |
-
2009
- 2009-03-11 CN CN2009100797959A patent/CN101833079B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7133480B2 (en) * | 2001-03-09 | 2006-11-07 | Leica Geosystems Inc. | Method and apparatus for processing digitally sampled signals at a resolution finer than that of a sampling clock |
Also Published As
Publication number | Publication date |
---|---|
CN101833079A (zh) | 2010-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061696B (zh) | 一种确定导航卫星轨道和钟差的方法 | |
CN102004258B (zh) | 基于多种gnss***融合的时间频率传递方法和接收机 | |
CN109581856B (zh) | 一种基于高性能晶振频率校准的对时守时方法 | |
CN101692163A (zh) | 一种频率标准远程校准方法及其*** | |
CN101833079B (zh) | 全球定位***原始测量值实时转为精确均匀采样的方法 | |
CN108363084A (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
Su et al. | Triple-frequency carrier phase precise time and frequency transfer models for BDS-3 | |
CN114994727B (zh) | 一种实现高精度时间校准及卫星定位的设备 | |
CN108958018A (zh) | 一种卫星授时方法和装置、计算机可读存储介质 | |
CN104238352A (zh) | 一种国家基准驾驭的地方时间标准生成***及方法 | |
CN110515109B (zh) | 一种基于多信息源融合自主pnt时间的方法及装置 | |
CN113325446B (zh) | 一种多模共频gnss载波相位时间传递方法及*** | |
CN107679260B (zh) | 一种静止轨道遥感卫星星上时统精度确定方法 | |
CN105578587B (zh) | 一种时间频率比对方法和设备 | |
CN212207982U (zh) | 一种原子钟校频*** | |
CN103148845A (zh) | 一种卫星组合罗盘及方位角和俯仰角测量计算方法 | |
Li et al. | Calibration errors in determining slant Total Electron Content (TEC) from multi-GNSS data | |
CN105137751A (zh) | 一种计量生产调度平台时间量值的校准***及其校准方法 | |
CN103543454A (zh) | 一种嵌入在移动通讯网中的卫星定轨*** | |
CN210742507U (zh) | 一种基于全球导航卫星***的标准时间频率源装置 | |
KR101725308B1 (ko) | 나노 초 레벨 정밀도를 가진 상대 시간 측정 시스템 | |
CN101907715B (zh) | 一种卫星导航数字基带信号的精密延迟实现方法 | |
CN110244332A (zh) | 一种天线阵大气相位扰动修正方法及*** | |
CN103901271A (zh) | 一种频率测试方法及*** | |
Byram et al. | High-precision GNSS orbit, clock and EOP estimation at the United States Naval Observatory |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent for invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: You Hongjian Inventor before: You Hongjian |
|
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: 20120822 Termination date: 20150311 |
|
EXPY | Termination of patent right or utility model |