CN108072887B - 单基站远距离海上实时动态定位方法 - Google Patents

单基站远距离海上实时动态定位方法 Download PDF

Info

Publication number
CN108072887B
CN108072887B CN201711083653.0A CN201711083653A CN108072887B CN 108072887 B CN108072887 B CN 108072887B CN 201711083653 A CN201711083653 A CN 201711083653A CN 108072887 B CN108072887 B CN 108072887B
Authority
CN
China
Prior art keywords
observation
satellite
simplified
time
station
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
Application number
CN201711083653.0A
Other languages
English (en)
Other versions
CN108072887A (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.)
Tongji University
Original Assignee
Tongji 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 Tongji University filed Critical Tongji University
Priority to CN201711083653.0A priority Critical patent/CN108072887B/zh
Publication of CN108072887A publication Critical patent/CN108072887A/zh
Application granted granted Critical
Publication of CN108072887B publication Critical patent/CN108072887B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/425—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between signals derived from different satellite radio beacon positioning systems
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明提供了一种单基站远距离海上实时动态定位的方法，利用北斗导航***的短报文功能传递简化后的差分改正数，解决海上通讯不便甚至没有通讯的难题，再利用北斗与GPS两个卫星导航***的多个频率多种观测值，针对海上的特殊环境，开发了稳定，快速，精准的流动端实时定位算法，最后实现了相应的基准站改正数播发软件以及流动站实时定位软件，搭建了一套基于单基站的长距离海上实时动态定位***。

Description

单基站远距离海上实时动态定位方法
技术领域
本发明涉及GNSS精密定位领域，尤其是海上的远距离高精度实时动态定位。
背景技术
自从GPS***问世以来，全球卫星导航***(GNSS)已广泛应用于各种户外定位与导航应用当中。实时高精度定位是GNSS的重要研究领域。目前，实时动态(RTK)技术，精密单点定位(PPP)技术以及他们的扩展网络RTK技术与PPP-RTK技术是高精度GNSS定位技术的代表。但是这些技术在海上应用中都有着各自的不足。
PPP技术需要精密星历、精密钟等信息，而这些信息都是在特定网站上延迟一周左右发布。至于一些超快速精密星历与IGS(International GPS Service)提供的实时精密星历、精密钟服务，在远海环境中也很难实时获取到。事实上，目前大部分所谓的实时PPP技术只是实现了算法上的实时而并没有投入实用中，仅在事后处理中模拟了实时的情况而已。此外，已有的一些可在全球范围内提供实时精密定位服务的商业软件，主要是通过卫星通讯来传递一些改正数信息，收费高昂且只能维持分米级的定位精度。不难看出，尽管实时PPP技术在理论上与实验中能够实现远海甚至全球的实时精密定位，但在实际应用中，还有一些问题需要解决。
目前RTK技术理论上的单站有效服务范围可达100公里甚至以上。在不考虑固定模糊度的前提下，适当放宽对精度的要求，RTK的单站服务范围还可以进一步扩展至几百公里。这意味着利用RTK技术在理论山可以为沿海几百公里范围的海域提供高精度实时动态定位。可事实上，受到通讯手段的限制，在远海没有网络的环境中，改正数通过电台发送无线电波进行传输，有效的传输距离只有三十公里左右，这使得RTK的单站服务范围远达不到它的极限距离。
与PPP相比，RTK似乎更不适合用于海上定位。然而，PPP技术不仅受到通讯条件的限制，还需要昂贵的硬件网基础设施支撑。此外RTK技术的单站服务范围若能到达几百公里便能满足许多远海的海上应用。因此，若通讯问题能得到妥善解决，由于开支较小且建设简便，RTK技术更适合用于海上的高精度导航。但是使用RTK技术实现单站几百公里服务范围的海上高精度实时动态定位还有许多问题需要解决。
发明内容
申请人研究发现，对于海上的通讯问题，卫星通讯是一种可行的解决方案。但是卫星通讯的费用较为昂贵，相比之下，北斗短报文通讯的服务费用要便宜很多，但是北斗短报文通讯的带宽很小，一次能传输的信息量十分有限。如果不对需要传递的改正数信息进行简化压缩处理，是无法使用北斗短报文传输改正数信息的。但是对改正数信息进行简化处理时需要尽可能避免对流动站的定位效果产生不良影响。因此，如何使用北斗短报文来传递基准站的改正数信息，并且不会影响到流动站的定位效果是一个重要问题。
申请人进一步研究发现，在海上，观测环境不同于陆地，多路径效应较为明显，并且设备需要安置在船只上，而船只上经常会有一些其它设备干扰到卫星信号，因此，在这种环境下得到的流动站观测数据往往会含有若干个粗差。对于单个粗差，可以在单点定位中利用相关分析法较好地探测出。但对于两个甚至更多粗差，由于卫星定位中的多余观测数较少，粗差探测的效果很差，经常会出现误判与漏判的情况。因此，如何尽可能地减少多个粗差对定位结果的影响是流动站定位算法需要解决的问题。
为了解决上述问题，本发明的目的在于提供一种单基站远距离海上实时动态定位的方法与***，利用北斗短报文传递改正数信息并考虑海上的特殊观测环境，以满足各种海上工程项目对定位的需求。
为了达到上述目的，本发明提供了一种单基站远距离海上实时动态定位方法，包括以下步骤：
在服务端实时接收基准站RTCM格式数据流，并对其进行解码，获得星历数据与当前历元实时的观测数据。对数据进行预处理，剔除不合格的观测数据。然后利用简化改正数算法对观测数据进行简化处理，从而获得简化改正数。
判断当前历元算得的简化改正数是否符合要求，若符合要求，则将改正数进行ASCII编码后进入下一历元的观测数据处理流程；否则，将弃用当前历元的简化改正数并进入下一历元的观测数据处理流程。
将编码后的简化改正数打包拆分成若干个短报文，交予服务端的北斗短报文天线进行发送，再由用户端的北斗短报文天线依次接收服务端天线发出的若干个北斗短报文，并将这些短报文重新拼接成为编码后的简化改正数，交予用户端定位软件使用。
用户端定位软件需同时接收流动站的RTCM格式数据流以及短报文天线发来的编码后简化改正数。从数据流中读取星历与当前历元流动站的观测数据，同样对数据进行预处理，剔除不合格的观测数据。对编码后简化改正数进行解码，判断简化改正数是否合格，若合格，则进行海上实时动态(Ocean-RTK，ORTK)定位；否则，只进行单点定位。最后，吐出定位结果后进入下一历元的观测数据处理流程。
优选的，对基准站观测数据与用户端数据进行预处理的步骤包括单点定位、伪距观测值检核及相关分析法粗差探测。伪距观测值检核是指利用已知的基准站坐标或者预测的用户端流动站坐标与星历中解算的卫星坐标反推对应的卫地距，从而可以与对应的伪距观测值进行比较，如果伪距观测值与反推的卫地距之间相差较大，则认为观测值含有粗差，应予以抛弃。相关分析法粗差探测是指在通过单点定位得到观测值残差序列后，利用基于相关分析的粗差理论分析观测值中可能包含的粗差，若探测到某观测值含有粗差则予以抛弃。
优选的，利用简化改正数算法对观测数据进行简化处理的步骤包括，在原始观测值中扣除算得的卫地距，消除观测值中由于卫星运动而产生的几何距离变化项；再在观测值中扣除对流层延迟的模型改正量，尽量减小原始观测值的数值变化，同时尽可能减小了其中的伪距观测值的绝对值大小，计算公式如下所示：
Figure BDA0001459565830000041
Figure BDA0001459565830000042
其中，
Figure BDA0001459565830000043
与
Figure BDA0001459565830000044
是改正后的相位与伪距观测值，
Figure BDA0001459565830000045
与
Figure BDA0001459565830000046
是原始相位与伪距观测值，下标r代表基准站，上标s代表某一卫星，
Figure BDA0001459565830000047
与
Figure BDA0001459565830000048
为单点定位算得的卫地距、接收机钟差与卫星钟差，
Figure BDA0001459565830000049
为对流层延迟的模型改正量。
最后再在其中的相位观测值中扣除其四舍五入整数部分而只保留小数，从而使相位观测值只在-0.5到0.5之间取值，计算公式如下所示：
Figure BDA00014595658300000410
Figure BDA00014595658300000411
其中，
Figure BDA00014595658300000412
是相位简化观测值，[*]表示对*进行四舍五入运算，λ是相位观测值
Figure BDA00014595658300000413
对应的波长。实际应用中，整数
Figure BDA00014595658300000414
在相位观测值
Figure BDA00014595658300000415
没有发生周跳时可以一直沿用下去，而不必每次都重新计算。
得到的简化改正数中的观测值包括，改正后伪距观测值
Figure BDA00014595658300000416
与相位简化观测值
Figure BDA00014595658300000417
优选的，判断当前历元算得的简化改正数是否符合要求的准则在于判断伪距或者相位观测值的绝对值大小是否超过1000，若不超过1000，则符合要求，否则不符合要求。
优选的，对改正数进行ASCII编码是指将所有的简化改正数信息，包括改正数观测时间、***标志、卫星数与简化后观测值等，转换成二进制的形式，然后将二进制形式的简化改正数信息再转换成ASCII字符串的形式进行传递。例如当用3个二进制位表示***标志为1的形式为001，用5个二进制位来表示卫星数为8的形式为01000，合并***标志与卫星数信息为00101000，转换成ASCII字符即为‘40’。具体的编码方法包括，二进制改正数的第0位到第9位用于存储简化改正数观测时间的后三位，第10位到第12位存储***标志，第13到第17位存储卫星数；从第18位开始存储每一颗观测到的卫星的信息，6位用于存储卫星号，5位用于存储卫星星历时间，1位用于存储
Figure BDA0001459565830000051
是否改变的标志，21位用于存储一个相位或者伪距的观测值。得到二进制改正数后，每8位可以转换成一个ASCII字符，最后不足8位则用0补齐8位后转换成ASCII字符。
优选的，对编码后的简化改正数打包拆分的步骤包括，将简化改正数信息分成若干个短报文，每个短报文都需要添加校验符以确保他们能够在用户端完成的恢复出来，同时，每个短报文的大小不能超过服务端北斗短报文天线一次可以发送的最大数据量。
优选的，由用户端将接收的若干个北斗短报文重新拼接成为编码后的简化改正数的步骤包括，利用每个短报文的校验符检核接收到的短报文是否属于同一个简化改正数信息或者同一简化改正数信息的短报文是否有缺漏，检核后，若不存在上述问题，则根据校验符，按顺序拼接若干个属于同一简化改正数信息的短报文，获得一条完整的编码后的简化改正数信息，否则，等待直至一条简化改正数信息的所有短报文被用户端接收到。
优选的，在用户端判断简化改正数是否合格的方法是，计算简化改正数的观测时间与流动站观测数据的观测时间之差，若时间差小于1分钟则简化改正数合格，否则，简化改正数不合格。
优选的，流动站的ORTK定位基于GPS与北斗两个导航***的双频相位与伪距观测值的无电离层(Ionosphere-Free,IF)组合进行，其步骤包括，利用流动站观测数据与接收到的简化改正数构造站间差分观测值，公式如下所示：
Figure BDA0001459565830000052
Figure BDA0001459565830000053
其中，
Figure BDA0001459565830000054
与
Figure BDA0001459565830000055
为简化后的站间差分相位与伪距观测值，
Figure BDA0001459565830000056
与
Figure BDA0001459565830000057
是基准站的相位简化观测值与改正后伪距观测值，
Figure BDA0001459565830000058
与
Figure BDA0001459565830000059
为流动站的原始相位和伪距观测值,
Figure BDA00014595658300000510
与
Figure BDA00014595658300000511
为单点定位算得的流动站近似卫地距、接收机钟差与卫星钟差，
Figure BDA00014595658300000512
为流动站对流层延迟的模型改正量。
得到站间差分观测值后，选取一颗合适的参考卫星，进行星间差分，得到双差观测值，公式如下所示：
Figure BDA0001459565830000061
Figure BDA0001459565830000062
其中，上标s代表选定的参考卫星，上标v代表其它任一卫星
Figure BDA0001459565830000063
与
Figure BDA0001459565830000064
代表简化后双差相位与伪距观测值，
Figure BDA0001459565830000065
与
Figure BDA0001459565830000066
代表卫星v的简化后站间差分相位与伪距观测值。
之后，考虑观测值有两个频率，可以对两个频率的双差观测值进行一个线性组合，得到IF组合的双差观测值，公式如下所示：
Figure BDA0001459565830000067
Figure BDA0001459565830000068
其中，
Figure BDA0001459565830000069
与
Figure BDA00014595658300000610
为IF组合的相位与伪距双差观测值，f1为第一个频率观测值的频率，f2为第二个频率观测值的频率，
Figure BDA00014595658300000611
与
Figure BDA00014595658300000612
为第一个频率的相位与伪距双差观测值，
Figure BDA00014595658300000613
与
Figure BDA00014595658300000614
为第二个频率的相位与伪距双差观测值。
由于不考虑残留的对流层延迟，基于IF组合双差观测值需要解算的参数只有流动站的位置参数与IF组合双差模糊度参数，可以得到如下所示的误差方程：
Figure BDA00014595658300000615
Figure BDA00014595658300000616
其中，
Figure BDA00014595658300000617
是双差卫地距，里面包含了三个未知的流动站位置参数，
Figure BDA00014595658300000618
为窄巷波长，c代表光速，
Figure BDA00014595658300000619
为IF组合双差模糊度参数，
Figure BDA00014595658300000620
与
Figure BDA00014595658300000621
是残留的相位和伪距观测值误差。
将误差方程写成矩阵形式：
Figure BDA00014595658300000622
其中ΦIF与PIF表示多个卫星相位与伪距双差观测值组成的观测向量，A是流动站位置参数的系数矩阵，Im是m阶单位矩阵，m是模糊度参数的个数，x代表流动站的位置参数向量，a代表模糊度参数向量，
Figure BDA00014595658300000623
与
Figure BDA00014595658300000624
是残留误差构成的误差向量。
流动站位置参数每一刻都由于流动站的运动而不相同，但模糊度参数是不随时间变化的，因此可以利用这一特性，联合多个时刻的观测值信息进行序贯解算。序贯解算中，流动站位置参数与模糊度参数的计算公式如下所示：
Figure BDA0001459565830000071
Figure BDA0001459565830000072
Figure BDA0001459565830000073
Figure BDA0001459565830000074
Figure BDA0001459565830000075
Figure BDA0001459565830000076
Figure BDA0001459565830000077
Figure BDA0001459565830000078
其中，
Figure BDA0001459565830000079
与
Figure BDA00014595658300000710
分别是两个频率原始相位与伪距观测值的先验中误差；Qdd是双差观测值的协因数矩阵，与频率无关；
Figure BDA00014595658300000711
与
Figure BDA00014595658300000712
分别是IF组合的伪距与相位双差观测值的协方差矩阵；下标k代表时刻，
Figure BDA00014595658300000713
与Ak分别代表k时刻的模糊度参数解，流动站位置参数解与流动站位置参数的系数矩阵；
Figure BDA00014595658300000714
是k时刻模糊度参数解的协方差矩阵；
Figure BDA00014595658300000715
是k时刻流动站位置参数解的协方差矩阵；当k＝0时，
Figure BDA00014595658300000716
和
Figure BDA00014595658300000717
为零向量或零矩阵。
优选的，流动站的ORTK定位有别于常规的RTK定位，它需要额外考虑流动站观测数据与简化改正数的观测时刻不一致以及海上多路径效应严重、观测环境特殊两个问题。观测数据与简化改正数的观测时刻不一致会导致各自的卫星钟差无法通过站间差分消除，因此必须在单点定位的时候，利用卫星的星历算得准确的各自观测时刻的卫星钟差，再如公式(1a)与(1b)所示，在原始的观测值中扣除。并且，必须保证流动站与基准站在计算卫星钟差时选用相同时刻的卫星星历。在海上，观测环境不同于陆地，多路径效应较为明显，并且设备需要安置在船只上，而船只上经常会有一些其它设备干扰到卫星信号，因此，在这种环境下得到的流动站观测数据往往会含有若干个粗差。对于单个粗差，可以在单点定位中利用相关分析法较好地探测出。但对于两个甚至更多粗差，由于卫星定位中的多余观测数较少，粗差探测的效果很差，经常会出现误判与漏判的情况，此时，由于ORTK定位基于序贯解算，因此某一时刻的粗差探测失败会影响到后续所有时刻的参数解算质量。为了避免出现这种情况，ORTK不仅在单点定位中加入相关分析法粗差探测，在之后的ORTK定位中也会利用IF组合的伪距双差观测值进行定位与粗差探测，并且会检查最后ORTK定位解算得到的参数的验后中误差。若粗差探测失败，参数的验后中误差会明显偏大甚至超过正常值几个数量级，此时将跳过该时刻而不把结果更新到序贯解算的参数中去，并且该时刻吐出单点定位的结果。
在本发明提供的单基站远距离海上实时动态定位方法中，首先在服务端实时接收基准站RTCM格式数据流，并对其进行解码，获得星历数据与当前历元实时的观测数据。再对数据进行预处理，剔除不合格的观测数据。然后利用简化改正数算法对观测数据进行简化处理，从而获得简化改正数。之后判断当前历元算得的简化改正数是否符合要求，若符合要求，则将改正数进行ASCII编码后进入下一历元的观测数据处理流程；否则，将弃用当前历元的简化改正数并进入下一历元的观测数据处理流程。接着将编码后的简化改正数打包拆分成若干个短报文，交予服务端的北斗短报文天线进行发送，再由用户端的北斗短报文天线依次接收服务端天线发出的若干个北斗短报文，并将这些短报文重新拼接成为编码后的简化改正数，交予用户端定位软件使用。而用户端定位软件需同时接收流动站的RTCM格式数据流以及短报文天线发来的编码后简化改正数。从数据流中读取星历与当前历元流动站的观测数据，同样对数据进行预处理，剔除不合格的观测数据。对编码后简化改正数进行解码，判断简化改正数是否合格，若合格，则进行海上实时动态(Ocean-RTK，ORTK)定位；否则，只进行单点定位。最后，吐出定位结果后进入下一历元的观测数据处理流程。本发明利用北斗短报文解决了海上通讯不便的问题；利用简化改正数算法以及相应的编码方法解决了北斗短报文通讯带宽过小，不适用于传递正常改正数信息的问题；利用多次粗差探测与验后中误差检验的方法，在多个粗差条件下，尽可能地避免了定位效果的持续性恶化，从而使用户端流动站定位算法能够适应海上的特殊环境。
附图说明
图1为本发明一实施例中单基站远距离海上实时动态定位方法的流程图。
具体实施方式
下面将结合示意图对本发明的具体实施方式进行更详细的描述。根据下列描述和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
本发明一实施例提供了一种单基站远距离海上实时动态定位方法，具体的，如图1所示，图1为本发明一实施例中单基站远距离海上实时动态定位方法的流程图，所述方法包括以下步骤：首先，在服务端实时接收基准站RTCM格式数据流，并对其进行解码，获得星历数据与当前历元实时的观测数据。如图1中的步骤S1所示。
然后对数据进行预处理，剔除不合格的观测数据，并利用简化改正数算法对观测数据进行简化处理，从而获得简化改正数。如图1中的步骤S2所示。
具体的，获得基准站观测数据后的数据预处理的步骤包括单点定位、伪距观测值检核及相关分析法粗差探测。伪距观测值检核是指利用已知的基准站坐标与星历中解算的卫星坐标反推对应的卫地距，从而可以与对应的伪距观测值进行比较，如果伪距观测值与反推的卫地距之间相差较大，则认为观测值含有粗差，应予以抛弃。相关分析法粗差探测是指在通过单点定位得到观测值残差序列后，利用基于相关分析的粗差理论分析观测值中可能包含的粗差，若探测到某观测值含有粗差则予以抛弃。基于相关分析的粗差理论是武汉大学的施闯等人总结提出的一套粗差探测理论，因此该部分已属于现有技术。
而利用简化改正数算法对观测数据进行简化处理的步骤包括，在原始观测值中扣除算得的卫地距，消除观测值中由于卫星运动而产生的几何距离变化项；再在观测值中扣除对流层延迟的模型改正量，尽量减小原始观测值的数值变化，同时尽可能减小了其中的伪距观测值的绝对值大小，计算公式如下所示：
Figure BDA0001459565830000101
Figure BDA0001459565830000102
其中，
Figure BDA0001459565830000103
与
Figure BDA0001459565830000104
是改正后的相位与伪距观测值，
Figure BDA0001459565830000105
与
Figure BDA0001459565830000106
是原始相位与伪距观测值，下标r代表基准站，上标s代表某一卫星，
Figure BDA0001459565830000107
与
Figure BDA0001459565830000108
为单点定位算得的卫地距、接收机种差与卫星钟差，
Figure BDA0001459565830000109
为对流层延迟的模型改正量。
最后再在其中的相位观测值中扣除其四舍五入整数部分而只保留小数，从而使相位观测值只在-0.5到0.5之间取值，计算公式如下所示：
Figure BDA00014595658300001010
Figure BDA00014595658300001011
其中，
Figure BDA00014595658300001012
是相位简化观测值，[*]表示对*进行四舍五入运算，λ是相位观测值
Figure BDA00014595658300001013
对应的波长。实际应用中，整数
Figure BDA00014595658300001014
在相位观测值
Figure BDA00014595658300001015
没有发生周跳时可以一直沿用下去，而不必每次都重新计算。
得到的简化改正数中的观测值包括，改正后伪距观测值
Figure BDA00014595658300001016
与相位简化观测值
Figure BDA00014595658300001017
步骤S3：判断当前历元算得的简化改正数是否符合要求。若符合要求，则将改正数进行ASCII编码后进入下一历元的观测数据处理流程，如步骤S4所示；否则，将弃用当前历元的简化改正数并进入下一历元的观测数据处理流程。
具体的，判断当前历元算得的简化改正数是否符合要求的准则在于判断伪距或者相位观测值的绝对值大小是否超过1000，若不超过1000，则符合要求，否则不符合要求。
对改正数进行ASCII编码是指将所有的简化改正数信息，包括改正数观测时间、***标志、卫星数与简化后观测值等，转换成二进制的形式，然后将二进制形式的简化改正数信息再转换成ASCII字符串的形式进行传递。例如当用3个二进制位表示***标志为1的形式为001，用5个二进制位来表示卫星数为8的形式为01000，合并***标志与卫星数信息为00101000，转换成ASCII字符即为‘40’。具体的编码方法包括，二进制改正数的第0位到第9位用于存储简化改正数观测时间的后三位，第10位到第12位存储***标志，第13到第17位存储卫星数；从第18位开始存储每一颗观测到的卫星的信息，6位用于存储卫星号，5位用于存储卫星星历时间，1位用于存储
Figure BDA0001459565830000111
是否改变的标志，21位用于存储一个相位或者伪距的观测值。得到二进制改正数后，每8位可以转换成一个ASCII字符，最后不足8位则用0补齐8位后转换成ASCII字符。
步骤S5：将编码后的简化改正数打包拆分成若干个短报文。
具体的，对编码后的简化改正数打包拆分的步骤包括，将简化改正数信息分成若干个短报文，每个短报文都需要添加校验符以确保他们能够在用户端完成的恢复出来，同时，每个短报文的大小不能超过服务端北斗短报文天线一次可以发送的最大数据量。
步骤S6：将打包后的数据交予服务端的北斗短报文天线进行发送。
具体的，使用北斗短报文传输改正数信息时需要注意北斗短报文的三个特点。第一，一张北斗卡一次只允许用户播发一条数据量不大于78.5字节的短报文。第二，一张北斗卡一分钟只能使用一次，一旦发出一条北斗短报文便将休眠一分钟，发出北斗短报文的时间很短，可以不计。第三，可以融合多张北斗卡来连续播发改正数信息，但最快只能一秒钟发出一条短报文。基于这三个特点，结合实际需求，选定需要融合的北斗卡数并确定相应地播发规则。本发明提供的实施例中，一个历元的简化改正数信息的数据量大概在300字节左右，考虑到改正数信息的时间延迟小于半分钟为宜，因此选择融合15张北斗卡并且每隔15秒播发一次改正数信息，从而使延迟在15-30秒左右。时间延迟的计算公式如下所示：
m<tj<60+2m-k-(j-1)d(6a)
m<ti<2m+d(6a)
Figure BDA0001459565830000121
Figure BDA0001459565830000122
Figure BDA0001459565830000123
其中，tj与ti代表第j与第i条简化改正数信息的时间延迟,且0<i<j；n是简化改正数信息经编码后的数据量，单位是字节；k是融合的北斗卡数。floor是向下取整运算。以上公式给出的是一分钟以内播发出的改正数信息时间延迟计算公式，一分钟以后播发的改正数信息的时间延迟与之近似。改正数的播发规则是播发出一条改正数信息后间隔d秒再播发下一条改正数信息，直到第j条改正数信息发出。
步骤S7：由用户端的北斗短报文天线依次接收服务端天线发出的若干个北斗短报文，并将这些短报文重新拼接成为编码后的简化改正数，交予用户端定位软件解码后使用。
具体的，由用户端将接收的若干个北斗短报文重新拼接成为编码后的简化改正数的步骤包括，利用每个短报文的校验符检核接收到的短报文是否属于同一个简化改正数信息或者同一简化改正数信息的短报文是否有缺漏，检核后，若不存在上述问题，则根据校验符，按顺序拼接若干个属于同一简化改正数信息的短报文，获得一条完整的编码后的简化改正数信息，否则，等待直至一条简化改正数信息的所有短报文被用户端接收到。
步骤S8：用户端定位软件需同时接收流动站的RTCM格式数据流以及短报文天线发来的编码后简化改正数。从数据流中读取星历与当前历元流动站的观测数据，同样对数据进行预处理，剔除不合格的观测数据并对编码后简化改正数进行解码。
具体的，对用户端数据进行预处理的步骤包括，单点定位、伪距观测值检核及相关分析法粗差探测。伪距观测值检核是指利用预测的流动站坐标与星历中解算的卫星坐标反推对应的近似卫地距，从而可以与对应的伪距观测值进行比较，如果伪距观测值与反推的近似卫地距之间相差较大，则认为观测值含有粗差，应予以抛弃。相关分析法粗差探测是指在通过单点定位得到观测值残差序列后，使用相关分析法分析观测值中可能包含的粗差，若探测到某观测值含有粗差则予以抛弃。
步骤S9：判断简化改正数是否合格。若合格，则进行海上实时动态(Ocean-RTK，ORTK)定位，如步骤S10所示；否则，只进行单点定位，如步骤S11所示。最后，吐出定位结果后进入下一历元的观测数据处理流程。
具体的，在用户端判断简化改正数是否合格的方法是，计算简化改正数的观测时间与流动站观测数据的观测时间之差，若时间差小于1分钟则简化改正数合格，否则，简化改正数不合格。
而流动站的ORTK定位基于GPS与北斗两个导航***的双频相位与伪距观测值的无电离层(Ionosphere-Free,IF)组合进行，其步骤包括，利用流动站观测数据与接收到的简化改正数构造站间差分观测值，公式如下所示：
Figure BDA0001459565830000131
Figure BDA0001459565830000132
其中，
Figure BDA0001459565830000133
与
Figure BDA0001459565830000134
为简化后的站间差分相位与伪距观测值，
Figure BDA0001459565830000135
与
Figure BDA0001459565830000136
是基准站的相位简化观测值与改正后伪距观测值，
Figure BDA0001459565830000137
与
Figure BDA0001459565830000138
为流动站的原始相位和伪距观测值,
Figure BDA0001459565830000139
与
Figure BDA00014595658300001310
为单点定位算得的流动站近似卫地距、接收机钟差与卫星钟差，
Figure BDA00014595658300001311
为流动站对流层延迟的模型改正量。
得到站间差分观测值后，选取一颗合适的参考卫星，进行星间差分，得到双差观测值，公式如下所示：
Figure BDA00014595658300001312
Figure BDA00014595658300001313
其中，上标s代表选定的参考卫星，上标v代表其它任一卫星
Figure BDA00014595658300001314
与
Figure BDA00014595658300001315
代表简化后双差相位与伪距观测值，
Figure BDA00014595658300001316
与
Figure BDA00014595658300001317
代表卫星v的简化后站间差分相位与伪距观测值。
之后，考虑观测值有两个频率，可以对两个频率的双差观测值进行一个线性组合，得到IF组合的双差观测值，公式如下所示：
Figure BDA00014595658300001318
Figure BDA0001459565830000141
其中，
Figure BDA0001459565830000142
与
Figure BDA0001459565830000143
为IF组合的相位与伪距双差观测值，f1为第一个频率观测值的频率，f2为第二个频率观测值的频率，
Figure BDA0001459565830000144
与
Figure BDA0001459565830000145
为第一个频率的相位与伪距双差观测值，
Figure BDA0001459565830000146
与
Figure BDA0001459565830000147
为第二个频率的相位与伪距双差观测值。
由于不考虑残留的对流层延迟，基于IF组合双差观测值需要解算的参数只有流动站的位置参数与IF组合双差模糊度参数，可以得到如下所示的误差方程：
Figure BDA0001459565830000148
Figure BDA0001459565830000149
其中，
Figure BDA00014595658300001410
是双差卫地距，里面包含了三个未知的流动站位置参数，
Figure BDA00014595658300001411
为窄巷波长，c代表光速，
Figure BDA00014595658300001412
为IF组合双差模糊度参数，
Figure BDA00014595658300001413
与
Figure BDA00014595658300001414
是残留的相位和伪距观测值误差。
将误差方程写成矩阵形式：
Figure BDA00014595658300001415
其中ΦIF与PIF表示多个卫星相位与伪距双差观测值组成的观测向量，A是流动站位置参数的系数矩阵，Im是m阶单位矩阵，m是模糊度参数的个数，x代表流动站的位置参数向量，a代表模糊度参数向量，
Figure BDA00014595658300001416
与
Figure BDA00014595658300001417
是残留误差构成的误差向量。
流动站位置参数每一刻都由于流动站的运动而不相同，但模糊度参数是不随时间变化的，因此可以利用这一特性，联合多个时刻的观测值信息进行序贯解算。序贯解算中，流动站位置参数与模糊度参数的计算公式如下所示：
Figure BDA00014595658300001418
Figure BDA00014595658300001419
Figure BDA00014595658300001420
Figure BDA00014595658300001421
Figure BDA00014595658300001422
Figure BDA00014595658300001423
Figure BDA00014595658300001424
Figure BDA0001459565830000151
其中，
Figure BDA0001459565830000152
与
Figure BDA0001459565830000153
分别是两个频率原始相位与伪距观测值的先验中误差；Qdd是双差观测值的协因数矩阵，与频率无关；
Figure BDA0001459565830000154
与
Figure BDA0001459565830000155
分别是IF组合的伪距与相位双差观测值的协方差矩阵；下标k代表时刻，
Figure BDA0001459565830000156
与Ak分别代表k时刻的模糊度参数解，流动站位置参数解与流动站位置参数的系数矩阵；
Figure BDA0001459565830000157
是k时刻模糊度参数解的协方差矩阵；
Figure BDA0001459565830000158
是k时刻流动站位置参数解的协方差矩阵；当k＝0时，
Figure BDA0001459565830000159
和
Figure BDA00014595658300001510
为零向量或零矩阵。
优选的，流动站的ORTK定位有别于常规的RTK定位，它需要额外考虑流动站观测数据与简化改正数的观测时刻不一致以及海上多路径效应严重、观测环境特殊两个问题。观测数据与简化改正数的观测时刻不一致会导致各自的卫星钟差无法通过站间差分消除，因此必须在单点定位的时候，利用卫星的星历算得准确的各自观测时刻的卫星钟差，再如公式(1a)与(1b)所示，在原始的观测值中扣除。并且，必须保证流动站与基准站在计算卫星钟差时选用相同时刻的卫星星历。在海上，观测环境不同于陆地，多路径效应较为明显，并且设备需要安置在船只上，而船只上经常会有一些其它设备干扰到卫星信号，因此，在这种环境下得到的流动站观测数据往往会含有若干个粗差。对于单个粗差，可以在单点定位中利用相关分析法较好地探测出。但对于两个甚至更多粗差，由于卫星定位中的多余观测数较少，粗差探测的效果很差，经常会出现误判与漏判的情况，此时，由于ORTK定位基于序贯解算，因此某一时刻的粗差探测失败会影响到后续所有时刻的参数解算质量。为了避免出现这种情况，ORTK不仅在单点定位中加入相关分析法粗差探测，在之后的ORTK定位中也会利用IF组合的伪距双差观测值进行定位与粗差探测，并且会检查最后ORTK定位解算得到的参数的验后中误差。若粗差探测失败，参数的验后中误差会明显偏大甚至超过正常值几个数量级，此时将跳过该时刻而不把结果更新到序贯解算的参数中去，并且该时刻吐出单点定位的结果。
与现有技术相比，本发明实施例所提供的单基站远距离海上实时动态定位方法至少具有以下有益效果：
第一，本发明实施例所提供的单基站远距离海上实时动态定位方法与IF组合观测值进行，再除中高纬度区域不需要考虑电离层对于定位的影响，理论上，定位精度只受残留的对流层误差影响，因而定位结果更为稳定。
第二，本发明实例所提供的单基站远距离海上实时动态定位方法在理论上并不固定窄巷模糊度，在距基站几十公里范围内可提供与固定模糊度时等效精度的定位服务，在距基站更远甚至几百公里的范围，可以提供精度略低的亚分米级定位服务。
第三，本发明实例所提供的单基站远距离海上实时动态定位方法在理论上没有忽略了天顶方向残留对流层延迟参数，从而在牺牲了一定高程方向的定位精度，但是提高了定位结果的可靠性，缩短了定位结果的收敛时间。
第四，本发明实例所提供的单基站远距离海上实时动态定位方法可以处理改正数信息与流动站观测数据观测时刻不一致的问题，从而使得改正数信息在用户端的有效使用时间更长，减少了用户端对改正数信息的需求频率。
第五，本发明利用北斗短报文解决了海上通讯不便，无法传输改正数信息的问题。相较于常规的卫星通讯，北斗短报文的成本更低。此外，本发明还提供了一种简化改正数算法以及相应的编码方法，从而大幅降低了改正数信息传递对通讯带宽的要求，使得通过北斗短报文传递改正数信息成为可能。
第六，本发明针对海上的特殊观测环境，开发了特定的用户端流动站定位软件，可以尽可能避免观测值中含有多个粗差的情形对于定位效果产生持续性的影响，从而使流动站的定位结果为稳定、可靠。
综上，在本发明实施例提供的单基站远距离海上实时动态定位方法中，首先在服务端实时接收基准站RTCM格式数据流，并对其进行解码，获得星历数据与当前历元实时的观测数据。再对数据进行预处理，剔除不合格的观测数据。然后利用简化改正数算法对观测数据进行简化处理，从而获得简化改正数。之后判断当前历元算得的简化改正数是否符合要求，若符合要求，则将改正数进行ASCII编码后进入下一历元的观测数据处理流程；否则，将弃用当前历元的简化改正数并进入下一历元的观测数据处理流程。接着将编码后的简化改正数打包拆分成若干个短报文，交予服务端的北斗短报文天线进行发送，再由用户端的北斗短报文天线依次接收服务端天线发出的若干个北斗短报文，并将这些短报文重新拼接成为编码后的简化改正数，交予用户端定位软件使用。而用户端定位软件需同时接收流动站的RTCM格式数据流以及短报文天线发来的编码后简化改正数。从数据流中读取星历与当前历元流动站的观测数据，同样对数据进行预处理，剔除不合格的观测数据。对编码后简化改正数进行解码，判断简化改正数是否合格，若合格，则进行海上实时动态(Ocean-RTK，ORTK)定位；否则，只进行单点定位。最后，吐出定位结果后进入下一历元的观测数据处理流程。本发明利用北斗短报文解决了海上通讯不便的问题；利用简化改正数算法以及相应的编码方法解决了北斗短报文通讯带宽过小，不适用于传递正常改正数信息的问题；利用多次粗差探测与验后中误差检验的方法，在多个粗差条件下，尽可能地避免了定位效果的持续性恶化，从而使用户端流动站定位算法能够适应海上的特殊环境。
上述仅为本发明的优选实施例而已，并不对本发明起到任何限制作用。任何所属技术领域的技术人员，在不脱离本发明的技术方案的范围内，对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动，均属未脱离本发明的技术方案的内容，仍属于本发明的保护范围之内。

Claims (8)

1.一种单基站远距离海上实时动态定位方法，其特征在于，包括以下步骤：
在服务端实时接收基准站RTCM格式数据流，并对其进行解码，获得星历数据与当前历元实时的观测数据；对数据进行预处理，剔除不合格的观测数据；然后利用简化改正数算法对观测数据进行简化处理，从而获得简化改正数；
判断当前历元算得的简化改正数是否符合要求，若符合要求，则将改正数进行ASCII编码后进入下一历元的观测数据处理流程；否则，将弃用当前历元的简化改正数并进入下一历元的观测数据处理流程；
将编码后的简化改正数打包拆分成若干个短报文，交予服务端的北斗短报文天线进行发送，再由用户端的北斗短报文天线依次接收服务端天线发出的若干个北斗短报文，并将这些短报文重新拼接成为编码后的简化改正数，交予用户端定位软件使用；
用户端定位软件需同时接收流动站的RTCM格式数据流以及短报文天线发来的编码后简化改正数；从数据流中读取星历与当前历元流动站的观测数据，同样对数据进行预处理，剔除不合格的观测数据；对编码后简化改正数进行解码，判断简化改正数是否合格，若合格，则进行海上实时动态(Ocean-RTK，ORTK)定位；否则，只进行单点定位；最后，吐出定位结果后进入下一历元的观测数据处理流程；
利用简化改正数算法对观测数据进行简化处理的步骤包括，在原始观测值中扣除算得的卫地距，消除观测值中由于卫星运动而产生的几何距离变化项；再在观测值中扣除对流层延迟的模型改正量，尽量减小原始观测值的数值变化，同时尽可能减小了其中的伪距观测值的绝对值大小，计算公式如下所示：
Figure FDA0002946719660000011
Figure FDA0002946719660000012
其中，
Figure FDA0002946719660000013
与
Figure FDA0002946719660000014
是改正后的相位与伪距观测值，
Figure FDA0002946719660000015
与
Figure FDA0002946719660000016
是原始相位与伪距观测值，下标r代表基准站，上标s代表某一卫星，
Figure FDA0002946719660000017
与
Figure FDA0002946719660000018
为单点定位算得的卫地距、接收机种差与卫星钟差，
Figure FDA0002946719660000021
为对流层延迟的模型改正量；
最后再在其中的相位观测值中扣除其四舍五入整数部分而只保留小数，从而使相位观测值只在-0.5到0.5之间取值，计算公式如下所示：
Figure FDA0002946719660000022
Figure FDA0002946719660000023
其中，
Figure FDA0002946719660000024
是相位简化观测值，[*]表示对*进行四舍五入运算，λ是相位观测值
Figure FDA0002946719660000025
对应的波长；实际应用中，整数
Figure FDA0002946719660000026
在相位观测值
Figure FDA0002946719660000027
没有发生周跳时可以一直沿用下去，而不必每次都重新计算；
得到的简化改正数中的观测值包括，改正后伪距观测值
Figure FDA0002946719660000028
与相位简化观测值
Figure FDA0002946719660000029
流动站的ORTK定位有别于常规的RTK定位，它需要额外考虑流动站观测数据与简化改正数的观测时刻不一致以及海上多路径效应严重、观测环境特殊两个问题；观测数据与简化改正数的观测时刻不一致会导致各自的卫星钟差无法通过站间差分消除，因此必须在单点定位的时候，利用卫星的星历算得准确的各自观测时刻的卫星钟差，再如公式(1a)与(1b)所示，在原始的观测值中扣除；并且，必须保证流动站与基准站在计算卫星钟差时选用相同时刻的卫星星历；在海上，观测环境不同于陆地，多路径效应较为明显，并且设备需要安置在船只上，而船只上经常会有一些其它设备干扰到卫星信号，因此，在这种环境下得到的流动站观测数据往往会含有若干个粗差；对于单个粗差，可以在单点定位中利用相关分析法较好地探测出；但对于两个甚至更多粗差，由于卫星定位中的多余观测数较少，粗差探测的效果很差，经常会出现误判与漏判的情况，此时，由于ORTK定位基于序贯解算，因此某一时刻的粗差探测失败会影响到后续所有时刻的参数解算质量；为了避免出现这种情况，ORTK不仅在单点定位中加入相关分析法粗差探测，在之后的ORTK定位中也会利用IF组合的伪距双差观测值进行定位与粗差探测，并且会检查最后ORTK定位解算得到的参数的验后中误差；若粗差探测失败，参数的验后中误差会明显偏大甚至超过正常值几个数量级，此时将跳过该时刻而不把结果更新到序贯解算的参数中去，并且该时刻吐出单点定位的结果。
2.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，对基准站观测数据与用户端数据进行预处理的步骤包括单点定位、伪距观测值检核及相关分析法粗差探测；伪距观测值检核是指利用已知的基准站坐标或者预测的用户端流动站坐标与星历中解算的卫星坐标反推对应的卫地距，从而可以与对应的伪距观测值进行比较，如果伪距观测值与反推的卫地距之间相差较大，则认为观测值含有粗差，应予以抛弃；相关分析法粗差探测是指在通过单点定位得到观测值残差序列后，利用基于相关分析的粗差理论分析观测值中可能包含的粗差，若探测到某观测值含有粗差则予以抛弃。
3.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，判断当前历元算得的简化改正数是否符合要求的准则在于判断伪距或者相位观测值的绝对值大小是否超过1000，若不超过1000，则符合要求，否则不符合要求。
4.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，对改正数进行ASCII编码是指将所有的简化改正数信息，包括改正数观测时间、***标志、卫星数与简化后观测值等，转换成二进制的形式，然后将二进制形式的简化改正数信息再转换成ASCII字符串的形式进行传递；当用3个二进制位表示***标志为1的形式为001，用5个二进制位来表示卫星数为8的形式为01000，合并***标志与卫星数信息为00101000，转换成ASCII字符即为‘40’；具体的编码方法包括，二进制改正数的第0位到第9位用于存储简化改正数观测时间的后三位，第10位到第12位存储***标志，第13到第17位存储卫星数；从第18位开始存储每一颗观测到的卫星的信息，6位用于存储卫星号，5位用于存储卫星星历时间，1位用于存储
Figure FDA0002946719660000031
是否改变的标志，21位用于存储一个相位或者伪距的观测值；得到二进制改正数后，每8位可以转换成一个ASCII字符，最后不足8位则用0补齐8位后转换成ASCII字符。
5.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，对编码后的简化改正数打包拆分的步骤包括，将简化改正数信息分成若干个短报文，每个短报文都需要添加校验符以确保他们能够在用户端完成的恢复出来，同时，每个短报文的大小不能超过服务端北斗短报文天线一次可以发送的最大数据量；使用北斗短报文传输改正数信息时需要注意北斗短报文的三个特点；第一，一张北斗卡一次只允许用户播发一条数据量不大于78.5字节的短报文；第二，一张北斗卡一分钟只能使用一次，一旦发出一条北斗短报文便将休眠一分钟，发出北斗短报文的时间很短，可以不计；第三，可以融合多张北斗卡来连续播发改正数信息，但最快只能一秒钟发出一条短报文；基于这三个特点，结合实际需求，选定需要融合的北斗卡数并确定相应地播发规则，从而使得改正数信息的时间延迟不会影响最终的定位效果；改正数信息的时间延迟是指简化改正数信息的观测时刻相对于流动站观测数据的观测时刻的时间延迟。
6.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，由用户端将接收的若干个北斗短报文重新拼接成为编码后的简化改正数的步骤包括，利用每个短报文的校验符检核接收到的短报文是否属于同一个简化改正数信息或者同一简化改正数信息的短报文是否有缺漏，检核后，若不存在上述问题，则根据校验符，按顺序拼接若干个属于同一简化改正数信息的短报文，获得一条完整的编码后的简化改正数信息，否则，等待直至一条简化改正数信息的所有短报文被用户端接收到。
7.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，在用户端判断简化改正数是否合格的方法是，计算简化改正数的观测时间与流动站观测数据的观测时间之差，若时间差小于1分钟则简化改正数合格，否则，简化改正数不合格。
8.如权利要求1所述的单基站远距离海上实时动态定位方法，其特征在于，流动站的ORTK定位基于GPS与北斗两个导航***的双频相位与伪距观测值的无电离层(Ionosphere-Free,IF)组合进行，其步骤包括，利用流动站观测数据与接收到的简化改正数构造站间差分观测值，公式如下所示：
Figure FDA0002946719660000041
Figure FDA0002946719660000042
其中，
Figure FDA0002946719660000051
与
Figure FDA0002946719660000052
为简化后的站间差分相位与伪距观测值，
Figure FDA0002946719660000053
与
Figure FDA0002946719660000054
是基准站的相位简化观测值与改正后伪距观测值，
Figure FDA0002946719660000055
与
Figure FDA0002946719660000056
为流动站的原始相位和伪距观测值,
Figure FDA0002946719660000057
与
Figure FDA0002946719660000058
为单点定位算得的流动站近似卫地距、接收机钟差与卫星钟差，
Figure FDA0002946719660000059
为流动站对流层延迟的模型改正量；
得到站间差分观测值后，选取一颗合适的参考卫星，进行星间差分，得到双差观测值，公式如下所示：
Figure FDA00029467196600000510
Figure FDA00029467196600000511
其中，上标s代表选定的参考卫星，上标v代表其它任一卫星
Figure FDA00029467196600000512
与
Figure FDA00029467196600000513
代表简化后双差相位与伪距观测值，
Figure FDA00029467196600000514
与
Figure FDA00029467196600000515
代表卫星v的简化后站间差分相位与伪距观测值；
之后，考虑观测值有两个频率，可以对两个频率的双差观测值进行一个线性组合，得到IF组合的双差观测值，公式如下所示：
Figure FDA00029467196600000516
Figure FDA00029467196600000517
其中，
Figure FDA00029467196600000518
与
Figure FDA00029467196600000519
为IF组合的相位与伪距双差观测值，f1为第一个频率观测值的频率，f2为第二个频率观测值的频率，
Figure FDA00029467196600000520
与
Figure FDA00029467196600000521
为第一个频率的相位与伪距双差观测值，
Figure FDA00029467196600000522
与
Figure FDA00029467196600000523
为第二个频率的相位与伪距双差观测值；
由于不考虑残留的对流层延迟，基于IF组合双差观测值需要解算的参数只有流动站的位置参数与IF组合双差模糊度参数，可以得到如下所示的误差方程：
Figure FDA00029467196600000524
Figure FDA00029467196600000525
其中，
Figure FDA00029467196600000526
是双差卫地距，里面包含了三个未知的流动站位置参数，
Figure FDA00029467196600000527
为窄巷波长，c代表光速，
Figure FDA00029467196600000528
为IF组合双差模糊度参数，
Figure FDA00029467196600000529
与
Figure FDA00029467196600000530
是残留的相位和伪距观测值误差；
将误差方程写成矩阵形式：
Figure FDA00029467196600000531
其中ΦIF与PIF表示多个卫星相位与伪距双差观测值组成的观测向量，A是流动站位置参数的系数矩阵，Im是m阶单位矩阵，m是模糊度参数的个数，x代表流动站的位置参数向量，a代表模糊度参数向量，
Figure FDA0002946719660000061
与
Figure FDA0002946719660000062
是残留误差构成的误差向量；
流动站位置参数每一刻都由于流动站的运动而不相同，但模糊度参数是不随时间变化的，因此可以利用这一特性，联合多个时刻的观测值信息进行序贯解算；序贯解算中，流动站位置参数与模糊度参数的计算公式如下所示：
Figure FDA0002946719660000063
Figure FDA0002946719660000064
Figure FDA0002946719660000065
Figure FDA0002946719660000066
Figure FDA0002946719660000067
Figure FDA0002946719660000068
Figure FDA0002946719660000069
Figure FDA00029467196600000610
其中，
Figure FDA00029467196600000611
与
Figure FDA00029467196600000612
分别是两个频率原始相位与伪距观测值的先验中误差；Qdd是双差观测值的协因数矩阵，与频率无关；
Figure FDA00029467196600000613
与
Figure FDA00029467196600000614
分别是IF组合的伪距与相位双差观测值的协方差矩阵；下标k代表时刻，
Figure FDA00029467196600000615
与Ak分别代表k时刻的模糊度参数解，流动站位置参数解与流动站位置参数的系数矩阵；
Figure FDA00029467196600000616
是k时刻模糊度参数解的协方差矩阵；
Figure FDA00029467196600000617
是k时刻流动站位置参数解的协方差矩阵；当k＝0时，
Figure FDA00029467196600000618
和
Figure FDA00029467196600000619
为零向量或零矩阵。
CN201711083653.0A 2017-11-07 2017-11-07 单基站远距离海上实时动态定位方法 Active CN108072887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711083653.0A CN108072887B (zh) 2017-11-07 2017-11-07 单基站远距离海上实时动态定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711083653.0A CN108072887B (zh) 2017-11-07 2017-11-07 单基站远距离海上实时动态定位方法

Publications (2)

Publication Number Publication Date
CN108072887A CN108072887A (zh) 2018-05-25
CN108072887B true CN108072887B (zh) 2021-04-30

Family

ID=62159739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711083653.0A Active CN108072887B (zh) 2017-11-07 2017-11-07 单基站远距离海上实时动态定位方法

Country Status (1)

Country Link
CN (1) CN108072887B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108616977B (zh) * 2018-06-21 2020-06-16 同济大学 一种基于最小二乘平差的uwb多点时延校正方法
CN108983267B (zh) * 2018-07-23 2022-04-12 中国石油大学(华东) 利用单台套北斗短报文实现远海精密定位的方法
CN109709591B (zh) * 2018-12-07 2021-04-20 中国科学院光电研究院 一种面向智能终端的gnss高精度定位方法
CN109521443B (zh) * 2018-12-29 2021-04-23 广东电网有限责任公司 一种探测星历异常的方法
CN109856659B (zh) * 2019-01-21 2021-02-12 同济大学 海床基定位授时与数据回收***及方法
CN109901196B (zh) * 2019-03-22 2023-04-21 中国石油大学(华东) 利用北斗短报文传输多频gnss数据和rts数据的方法
CN110988948B (zh) * 2019-11-07 2021-11-02 北京航空航天大学 一种基于动对动相对定位场景中完好性分析方法
CN111045042B (zh) * 2019-12-20 2022-03-04 西安空间无线电技术研究所 一种基于“云-端”架构的ppp-rtk增强方法与***
CN111385743B (zh) * 2020-03-11 2021-03-19 北京邮电大学 一种位置估计方法、装置及电子设备
CN111896988B (zh) * 2020-07-28 2022-11-08 中国石油大学(华东) 远海实时定位方法、***、介质、卫星导航、定位***
CN112051597A (zh) * 2020-09-27 2020-12-08 西安领航遥感技术有限公司 一种精密单点定位方法与装置
CN113075710B (zh) * 2020-10-21 2021-11-19 中移(上海)信息通信科技有限公司 定位差分信息获取方法、装置、设备及计算机存储介质
CN114935768B (zh) * 2022-07-13 2022-11-04 武汉大学 一种基于单基站的虚拟参考站的构建方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8259008B2 (en) * 2008-11-17 2012-09-04 Qualcomm Incorporated DGNSS correction for positioning
US8174437B2 (en) * 2009-07-29 2012-05-08 Hemisphere Gps Llc System and method for augmenting DGNSS with internally-generated differential correction
CN104133235B (zh) * 2014-07-28 2016-08-24 辽宁工程技术大学 一种gps/bds非差网络rtk误差改正数传输方法
EP3361287B1 (en) * 2015-10-29 2021-10-13 Huawei Technologies Co., Ltd. Positioning method, base station and server of a mobile network
CN106371115B (zh) * 2016-06-27 2019-01-25 上海达华测绘有限公司 基于北斗短报文技术多模多频海上精密定位方法
CN106249264B (zh) * 2016-07-19 2019-03-22 深圳市安煋信息技术有限公司 一种定位方法、***及导航装置
CN107153209B (zh) * 2017-07-06 2019-07-30 武汉大学 一种短弧段低轨导航卫星实时精密定轨方法

Also Published As

Publication number Publication date
CN108072887A (zh) 2018-05-25

Similar Documents

Publication Publication Date Title
CN108072887B (zh) 单基站远距离海上实时动态定位方法
CN110187364B (zh) 一种低轨导航增强精密改正数据生成、上注***及方法
CN106371115B (zh) 基于北斗短报文技术多模多频海上精密定位方法
CN111290005B (zh) 载波相位的差分定位方法、装置、电子设备及存储介质
CN108919305B (zh) 交通运输中北斗地基增强带状服务方法及***
CN106873009B (zh) 利用非组合ppp辅助的长基线模糊度解算方法
US9354319B2 (en) Ambiguity windowing in communications among global navigation system satellite receivers
CN104459737A (zh) 一种基于信标实时差分的卫星定位方法
CN108076662A (zh) 具有使用未组合公式来解算模糊度的能力的gnss接收机
CN108387169B (zh) 一种基于实时大气产品的gnss形变监测***
CN104459745A (zh) 一种多星座长基线网络rtk部分模糊度快速解算方法
Tomasz GNSS-warp software for real-time precise point positioning
KR101419339B1 (ko) 위성항법 보정시스템의 보정정보 처리 장치 및 그 방법
CN111290004A (zh) 伪距差分定位方法、装置、电子设备及存储介质
CN103728643A (zh) 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法
CN105204048A (zh) 一种北斗兼容gps/glonass的网络rtk参考站间模糊度快速解算方法
CN101435861A (zh) 弱信号搜星环境下的gps信号处理方法
CN104898145A (zh) 一种基于半周模糊度的模糊度固定方法和***
AU2009322084A1 (en) A system and method of reference position determination
CN113325446A (zh) 一种多模共频gnss载波相位时间传递方法及***
KR101067416B1 (ko) 위성을 이용한 상대 측위 방법 및 상대 측위 시스템
CN114460608A (zh) 一种gnss差分数据的压缩编码与使用策略
CN104898144A (zh) 一种北斗***参考站整周模糊度单历元确定方法
CN104749587A (zh) 接收机伪距故障监测方法和接收机
CN104991265A (zh) 一种北斗卫星导航***用户统一性定位方法

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