CN101833079A - 全球定位系统原始测量值实时转为精确均匀采样的方法 - Google Patents
全球定位系统原始测量值实时转为精确均匀采样的方法 Download PDFInfo
- Publication number
- CN101833079A CN101833079A CN200910079795A CN200910079795A CN101833079A CN 101833079 A CN101833079 A CN 101833079A CN 200910079795 A CN200910079795 A CN 200910079795A CN 200910079795 A CN200910079795 A CN 200910079795A CN 101833079 A CN101833079 A CN 101833079A
- Authority
- CN
- China
- Prior art keywords
- gps
- time
- 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.)
- Granted
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种全球定位系统原始测量值实时转为精确均匀采样的方法,涉及全球定位技术,利用全球定位系统(GPS,Global Positioning System)接收机实时计算的钟差(接收机时间和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接收机时间系统中名义上按照固定的接收机时间间隔进行采样得到的。
根据实时测量的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实时测量的数据如下表所示:
实际对应的GPS时间和数据如下表所示(测量值没有变化):
实时修正后和GPS时间严格同步且均匀采样的测量值如下表所示:
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 true CN101833079A (zh) | 2010-09-15 |
CN101833079B 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) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102098119A (zh) * | 2011-01-27 | 2011-06-15 | 电子科技大学 | 一种接收机钟差预测方法 |
CN102890279A (zh) * | 2011-07-20 | 2013-01-23 | 中国科学院电子学研究所 | 一种利用全球定位系统进行纳秒级精确时间同步的方法 |
CN107003402A (zh) * | 2014-12-19 | 2017-08-01 | 高通股份有限公司 | 具有用于监测和检测监控区域内的目标的若干收发器的雷达系统 |
CN107003413A (zh) * | 2014-12-04 | 2017-08-01 | 罗伯特·博世有限公司 | 用于确定原始测量值的统计特性的方法和设备 |
CN110568466A (zh) * | 2019-08-19 | 2019-12-13 | 北京自动化控制设备研究所 | 基准站观测量推算方法及载波相位差分定位信息获取方法 |
Family Cites Families (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
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102098119A (zh) * | 2011-01-27 | 2011-06-15 | 电子科技大学 | 一种接收机钟差预测方法 |
CN102098119B (zh) * | 2011-01-27 | 2013-09-25 | 电子科技大学 | 一种接收机钟差预测方法 |
CN102890279A (zh) * | 2011-07-20 | 2013-01-23 | 中国科学院电子学研究所 | 一种利用全球定位系统进行纳秒级精确时间同步的方法 |
CN102890279B (zh) * | 2011-07-20 | 2014-02-12 | 中国科学院电子学研究所 | 一种利用全球定位系统进行纳秒级精确时间同步的方法 |
CN107003413A (zh) * | 2014-12-04 | 2017-08-01 | 罗伯特·博世有限公司 | 用于确定原始测量值的统计特性的方法和设备 |
CN107003402A (zh) * | 2014-12-19 | 2017-08-01 | 高通股份有限公司 | 具有用于监测和检测监控区域内的目标的若干收发器的雷达系统 |
CN107003402B (zh) * | 2014-12-19 | 2020-05-22 | 韦特里西提公司 | 具有用于监测和检测监控区域内的目标的若干收发器的雷达系统 |
CN110568466A (zh) * | 2019-08-19 | 2019-12-13 | 北京自动化控制设备研究所 | 基准站观测量推算方法及载波相位差分定位信息获取方法 |
CN110568466B (zh) * | 2019-08-19 | 2021-09-14 | 北京自动化控制设备研究所 | 基准站观测量推算方法及载波相位差分定位信息获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101833079B (zh) | 2012-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101692163B (zh) | 一种频率标准远程校准方法及其系统 | |
CN102004258B (zh) | 基于多种gnss系统融合的时间频率传递方法和接收机 | |
CN109581856B (zh) | 一种基于高性能晶振频率校准的对时守时方法 | |
CN105607089A (zh) | 一种基于伪距辅助载波相位的北斗共视时间传递的方法 | |
CN101833079B (zh) | 全球定位系统原始测量值实时转为精确均匀采样的方法 | |
CN107679260B (zh) | 一种静止轨道遥感卫星星上时统精度确定方法 | |
CN108958018A (zh) | 一种卫星授时方法和装置、计算机可读存储介质 | |
CN104238352A (zh) | 一种国家基准驾驭的地方时间标准生成系统及方法 | |
CN212278244U (zh) | 一种基于卫星共视的高精度网络授时装置 | |
CN110515109B (zh) | 一种基于多信息源融合自主pnt时间的方法及装置 | |
CN111221007A (zh) | 一种基于北斗共视的新型精密授时装置 | |
CN105578587B (zh) | 一种时间频率比对方法和设备 | |
CN212207982U (zh) | 一种原子钟校频系统 | |
CN105137751A (zh) | 一种计量生产调度平台时间量值的校准系统及其校准方法 | |
CN209486181U (zh) | 一种准确测量频率的系统 | |
KR101725308B1 (ko) | 나노 초 레벨 정밀도를 가진 상대 시간 측정 시스템 | |
CN103698781B (zh) | 一种接收机观测信息的提取方法和装置 | |
CN210742507U (zh) | 一种基于全球导航卫星系统的标准时间频率源装置 | |
CN103901271A (zh) | 一种频率测试方法及系统 | |
CN110149197A (zh) | 一种用于时钟同步系统的高精度同步方法及同步系统 | |
CN202533743U (zh) | 时间分析仪测量基准系统 | |
CN104614981B (zh) | 一种获得原子钟远程校准系统的天频率稳定度的方法 | |
CN212207988U (zh) | 一种时码精度测量装置 | |
CN105589328B (zh) | 时间同步测试方法、测试精度确定方法及装置 | |
CN211180593U (zh) | 时频检测设备 |
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 |