CN112731493B - 一种精度高的静态伪距单点定位方法及装置 - Google Patents

一种精度高的静态伪距单点定位方法及装置 Download PDF

Info

Publication number
CN112731493B
CN112731493B CN202011514751.7A CN202011514751A CN112731493B CN 112731493 B CN112731493 B CN 112731493B CN 202011514751 A CN202011514751 A CN 202011514751A CN 112731493 B CN112731493 B CN 112731493B
Authority
CN
China
Prior art keywords
epoch
receiver
station
satellite
current epoch
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
CN202011514751.7A
Other languages
English (en)
Other versions
CN112731493A (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.)
Fujian Zhongliang Zhihui Technology Co ltd
Minjiang University
Original Assignee
Fujian Zhongliang Zhihui Technology Co ltd
Minjiang 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 Fujian Zhongliang Zhihui Technology Co ltd, Minjiang University filed Critical Fujian Zhongliang Zhihui Technology Co ltd
Priority to CN202011514751.7A priority Critical patent/CN112731493B/zh
Publication of CN112731493A publication Critical patent/CN112731493A/zh
Application granted granted Critical
Publication of CN112731493B publication Critical patent/CN112731493B/zh
Active 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/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
    • 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
    • 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/28—Satellite selection

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明提供的一种精度高的静态伪距单点定位方法及装置，通过计算每个历元中每颗卫星的高度角和方位角，删除高度角小于预设值的卫星，并根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵，根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置，将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，通过多历元解算并且求平均，有效的削弱伪距误差对单点定位的影响，从而提高定位精度且无需与固定基站配合。

Description

一种精度高的静态伪距单点定位方法及装置
本案是以申请号为201811166743.0，申请日为2018年10月8日，名称为《一种静态伪距单点定位的方法及装置》的专利申请为母案的分案申请。
技术领域
本发明涉及全球卫星导航系统(Global Navigation Satellite System，GNSS)精度高的静态伪距单点定位领域，特别涉及一种精度高的静态伪距单点定位方法及装置。
背景技术
当前，单频接收机因其成本低、尺寸小及重量轻等特点被广泛用在资源调查、徒步旅游和车、船导航等领域。随着全球卫星导航系统相关技术的发展和进步，用户对精度高的静态伪距单点定位的精度要求也越来越高。通常，单频接收机受卫星轨道、大气传播及接收机自身误差的影响，目前精度高的静态伪距单点定位的精度仅在5米左右。
近几年，许多研究人员就如何提高精度高的静态伪距单点定位精度进行了深入研究。其中伪距差分方法和加权最小二乘法最为常见。伪距差分方法可实现厘米级的定位精度，但这种方法需要一个固定基站与之配合，且定位精度受固定基站的位置精度和基线长度的影响，在诸如山区、戈壁及沙漠等环境中应用受到限制；常见的加权最小二乘法，如基于高度角的加权最小二乘法，虽然能在一定程度上提高定位精度，但是在城市峡谷等环境中，低高度角的卫星受多路径效应影响明显，此时，该加权方法对定位精度的提高效果不明显，在精度高的静态伪距单点定位性能的改善上受到限制。
因此需要一种精度高的静态伪距单点定位方法解决上述技术问题。
发明内容
本发明所要解决的技术问题是：提供一种精度高的静态伪距单点定位方法及装置，能够有效削弱伪距误差对单点定位的影响，提高定位精度且无需与固定基站配合。
为了解决上述技术问题，本发明采用的一种技术方案为：
一种精度高的静态伪距单点定位方法，包括步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果。
为了解决上述技术问题，本发明采用的另一种技术方案为：
一种精度高的静态伪距单点定位的装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现以下步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果。
本发明的有益效果在于：通过计算每个历元中每颗卫星的高度角和方位角，删除高度角小于预设值的卫星，并根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵，根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置，将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，通过多历元解算并且求平均，有效的削弱伪距误差对单点定位的影响，从而提高定位精度且无需与固定基站配合。
附图说明
图1为本发明实施例的精度高的静态伪距单点定位方法流程图；
图2为本发明实施例的精度高的静态伪距单点定位的装置的结构示意图；
图3为本发明实施例的在E、N、U方向上的加权策略示意图；
图4为本发明实施例的精度高的静态伪距单点定位方法与高度角加权的最小二乘方法定位的精度对比；
标号说明：
1、精度高的静态伪距单点定位的装置；2、存储器；3、处理器。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果，以下结合实施方式并配合附图予以说明。
本发明最关键的构思在于:通过确定每个历元各卫星在站心坐标系中E、N、U方向上的权重矩阵，根据最小二乘法求得接收机在各历元中的位置后取平均，得到接收机的精度高的静态伪距单点定位结果，有效的削弱伪距误差对单点定位的影响，从而提高定位精度且无需与固定基站配合。
请参照图1，一种精度高的静态伪距单点定位方法，包括步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果。
从上述描述可知，本发明的有益效果在于：通过计算每个历元中每颗卫星的高度角和方位角，删除高度角小于预设值的卫星，并根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵，根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置，将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，通过多历元解算并且求平均，有效的削弱伪距误差对单点定位的影响，从而提高定位精度且无需与固定基站配合。
进一步的，步骤S11包括：
S111、判断当前历元是否为第一历元，若是，则执行步骤S112，否则，执行步骤S113；
S112、根据卫星位置以及接收机位置的初始值将计算卫星的高度角和方位角；
S113、根据卫星位置和上一历元估计的接收机位置计算卫星的高度角和方位角。
由上述描述可知，通过判断当前历元是否为第一历元，并采用不同历元的接收机位置计算卫星的高度角和方位角，保证了后续计算权重矩阵的准确性。
进一步的，步骤S13包括：
S131、根据步骤S12处理后的每颗卫星的高度角和方位角确定当前历元中第i颗卫星在站心坐标系中E、N、U方向上的权值w2 E,i、w2 N,i和w2 U,i，其中i＝(1，2，...，n)，n为当前历元中所述接收机收到的卫星数量；
S132、根据所述E、N、U方向上各颗卫星的权值确定当前历元在E、N、U方向上关于各颗卫星的权重矩阵，所述权重矩阵如下：
WE＝diag{w2 E,1,w2 E,2,...,w2 E,n}
WN＝diag{w2 N,1,w2 N,2,...,w2 N,n}
WU＝diag{w2 U,1,w2 U,2,...,w2 U,n}。
由上述描述可知，通过确定当前历元中每颗卫星在站心坐标系中不同方向上的权值，并根据不同方向上的权值确定当前历元各颗卫星在E、N、U方向上的权重矩阵，不仅提高了定位的准确性，还可满足用户对不同方向上定位精度的需求。
进一步的，步骤S14包括：
S141、分步计算当前历元和上一历元的站心坐标差，计算公式如下：
ΔEE＝(GTWEG)-1GTWEΔρ
ΔEN＝(GTWNG)-1GTWNΔρ
ΔEU＝(GTWUG)-1GTWUΔρ
其中，ΔE为当前历元和上一历元的站心坐标差，ΔEE表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差，ΔEN表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差，ΔEU表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，G为方向余弦矩阵，Δρ为修正后的伪距残差；
S142、根据牛顿迭代法，重复步骤S111至S141，直到小于预设值时，得到所述接收机在当前历元E、N、U方向上的站心坐标的估值如下：
其中，eE,k，nE,k，uE,k分别表示当前历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，eE,k-1，nE,k-1，uE,k-1分别表示上一历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，ΔeE，ΔnE，ΔuE分别表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差；eN,k，nN,k，uN,k分别表示当前历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，eN,k-1，nN,k-1，uN,k-1分别表示上一历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，ΔeN，ΔnN，ΔuN分别表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差；eU,k，nU,k，uU,k分别表示当前历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，eU,k-1，nU,k-1，uU,k-1分别表示上一历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，ΔeU，ΔnU，ΔuU分别表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，k表示迭代次数，且k为大于1的正整数；
S143、取eE,k，nN,k，uU,k作为当前历元接收机在站心坐标系中的坐标，所述接收机的站心坐标如下：
S144、根据所述接收机的站心坐标和坐标变换公式进行坐标变换，得到所述接收机的单点定位结果；
所述坐标变换公式为：其中，x，y，z分别表示接收机在地心地固坐标系X轴、Y轴和Z轴方向上的坐标分量，S为坐标变换矩阵，λ为接收机位置的大地经度，φ为接收机位置的大地纬度。
由上述描述可知，通过分步计算当前历元和上一历元的站心坐标差，并在不断迭代后将不同方向上的最优估值eE,k，nN,k，uU,k作为当前历元接收机坐标，通过坐标变换公式对接收机的站心坐标进行转换，得到接收机的单点定位结果，便于使用者使用。
进一步的，步骤S2包括：
根据取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，所述取平均公式如下：
其中，为各个历元的位置在X轴、Y轴和Z轴方向上坐标分量的平均值，m为历元总数。
由上述描述可知，通过取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，有利于消除定位误差，提高精度高的静态伪距单点定位的精确度。
请参照图2，一种精度高的静态伪距单点定位的装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现以下步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果。
从上述描述可知，本发明的有益效果在于：通过计算每个历元中每颗卫星的高度角和方位角，删除高度角小于预设值的卫星，并根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵，根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置，将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，通过多历元解算并且求平均，有效的削弱伪距误差对单点定位的影响，从而提高定位精度且无需与固定基站配合。
进一步的，步骤S11包括：
S111、判断当前历元是否为第一历元，若是，则执行步骤S112，否则，执行步骤S113；
S112、根据卫星位置以及接收机位置的初始值将计算卫星的高度角和方位角；
S113、根据卫星位置和上一历元估计的接收机位置计算卫星的高度角和方位角。
由上述描述可知，通过判断当前历元是否为第一历元，并采用不同历元的接收机位置计算卫星的高度角和方位角，保证了后续计算权重矩阵的准确性。
进一步的，步骤S13包括：
S131、根据步骤S12处理后的每颗卫星的高度角和方位角确定当前历元中第i颗卫星在站心坐标系中E、N、U方向上的权值w2 E,i、w2 N,i和w2 U,i，其中i＝(1，2，...，n)，n为当前历元中所述接收机收到的卫星数量；
S132、根据所述E、N、U方向上各颗卫星的权值确定当前历元在E、N、U方向上关于各颗卫星的权重矩阵，所述权重矩阵如下：
WE＝diag{w2 E,1,w2 E,2,...,w2 E,n}
WN＝diag{w2 N,1,w2 N,2,...,w2 N,n}
WU＝diag{w2 U,1,w2 U,2,...,w2 U,n}。
由上述描述可知，通过确定当前历元中每颗卫星在站心坐标系中不同方向上的权值，并根据不同方向上的权值确定当前历元各颗卫星在E、N、U方向上的权重矩阵，不仅提高了定位的准确性，还可满足用户对不同方向上定位精度的需求。
进一步的，步骤S14包括：
S141、分步计算当前历元和上一历元的站心坐标差，计算公式如下：
ΔEE＝(GTWEG)-1GTWEΔρ
ΔEN＝(GTWNG)-1GTWNΔρ
ΔEU＝(GTWUG)-1GTWUΔρ
其中，ΔE为当前历元和上一历元的站心坐标差，ΔEE表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差，ΔEN表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差，ΔEU表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，G为方向余弦矩阵，Δρ为修正后的伪距残差；
S142、根据牛顿迭代法，重复步骤S111至S141，直到小于预设值时，得到所述接收机在当前历元E、N、U方向上的站心坐标的估值如下：
其中，eE,k，nE,k，uE,k分别表示当前历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，eE,k-1，nE,k-1，uE,k-1分别表示上一历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，ΔeE，ΔnE，ΔuE分别表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差；eN,k，nN,k，uN,k分别表示当前历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，eN,k-1，nN,k-1，uN,k-1分别表示上一历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，ΔeN，ΔnN，ΔuN分别表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差；eU,k，nU,k，uU,k分别表示当前历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，eU,k-1，nU,k-1，uU,k-1分别表示上一历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，ΔeU，ΔnU，ΔuU分别表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，k表示迭代次数，且k为大于1的正整数；
S143、取eE,k，nN,k，uU,k作为当前历元接收机在站心坐标系中的坐标，所述接收机的站心坐标如下：
S144、根据所述接收机的站心坐标和坐标变换公式进行坐标变换，得到所述接收机的单点定位结果；
所述坐标变换公式为：其中，x，y，z分别表示接收机在地心地固坐标系X轴、Y轴和Z轴方向上的坐标分量，S为坐标变换矩阵，λ为接收机位置的大地经度，φ为接收机位置的大地纬度。
由上述描述可知，通过分步计算当前历元和上一历元的站心坐标差，并在不断迭代后将不同方向上的最优估值eE,k，nN,k，uU,k作为当前历元接收机坐标，通过坐标变换公式对接收机的站心坐标进行转换，得到接收机的单点定位结果，便于使用者使用。
进一步的，步骤S2包括：
根据取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，所述取平均公式如下：
其中，为各个历元的位置在X轴、Y轴和Z轴方向上坐标分量的平均值，m为历元总数。
由上述描述可知，通过取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，有利于消除定位误差，提高精度高的静态伪距单点定位的精确度。
实施例一
请参照图1，一种精度高的静态伪距单点定位方法，包括步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S111、判断当前历元是否为第一历元，若是，则执行步骤S112，否则，执行步骤S113；
S112、根据卫星位置以及接收机位置的初始值将计算卫星的高度角和方位角；
S113、根据卫星位置和上一历元估计的接收机位置计算卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S131、根据步骤S12处理后的每颗卫星的高度角和方位角确定当前历元中第i颗卫星在站心坐标系中E、N、U方向上的权值w2 E,i、w2 N,i和w2 U,i，其中i＝(1，2，...，n)，n为当前历元中所述接收机收到的卫星数量；
S132、根据所述E、N、U方向上各颗卫星的权值确定当前历元在E、N、U方向上关于各颗卫星的权重矩阵，所述权重矩阵如下：
WE＝diag{w2 E,1,w2 E,2,...,w2 E,n}
WN＝diag{w2 N,1,w2 N,2,...,w2 N,n}
WU＝diag{w2 U,1,w2 U,2,...,w2 U,n}；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S141、分步计算当前历元和上一历元的站心坐标差，计算公式如下：
ΔEE＝(GTWEG)-1GTWEΔρ
ΔEN＝(GTWNG)-1GTWNΔρ
ΔEU＝(GTWUG)-1GTWUΔρ
其中，ΔE为当前历元和上一历元的站心坐标差，ΔEE表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差，ΔEN表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差，ΔEU表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，G为方向余弦矩阵，Δρ为修正后的伪距残差；
S142、根据牛顿迭代法，重复步骤S111至S141，直到小于预设值时，得到所述接收机在当前历元E、N、U方向上的站心坐标的估值如下：
其中，eE,k，nE,k，uE,k分别表示当前历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，eE,k-1，nE,k-1，uE,k-1分别表示上一历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，ΔeE，ΔnE，ΔuE分别表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差；eN,k，nN,k，uN,k分别表示当前历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，eN,k-1，nN,k-1，uN,k-1分别表示上一历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，ΔeN，ΔnN，ΔuN分别表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差；eU,k，nU,k，uU,k分别表示当前历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，eU,k-1，nU,k-1，uU,k-1分别表示上一历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，ΔeU，ΔnU，ΔuU分别表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，k表示迭代次数，且k为大于1的正整数；
S143、取eE,k，nN,k，uU,k作为当前历元接收机在站心坐标系中的坐标，所述接收机的站心坐标如下：
S144、根据所述接收机的站心坐标和坐标变换公式进行坐标变换，得到所述接收机的单点定位结果；
所述坐标变换公式为：其中，x，y，z分别表示接收机在地心地固坐标系X轴、Y轴和Z轴方向上的坐标分量，S为坐标变换矩阵，λ为接收机位置的大地经度，φ为接收机位置的大地纬度；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果；
根据取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，所述取平均公式如下：
其中，为各个历元的位置在X轴、Y轴和Z轴方向上坐标分量的平均值，m为历元总数。
实施例二
本实施例将结合具体的应用场景，进一步说明本发明上述精度高的静态伪距单点定位方法是如何实现的：
本发明适用于静态用户，可推广至资源调查和徒步旅游等方面；
1、预设多个历元，并对每个历元分别执行以下步骤：
1.1、基于广播星历计算当前历元中每颗卫星的位置；
1.2、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角θ和方位角α；
1.21、判断当前历元是否为第一历元，若是，则执行步骤1.22，否则，执行步骤1.23；
1.22、根据卫星位置以及接收机位置的初始值将计算卫星的高度角θ和方位角α；
1.23、根据卫星位置和上一历元估计的接收机位置计算卫星的高度角θ和方位角α；
其中，Δx、Δy和Δz为接收机到卫星的观测向量；
1.3、删除高度角小于预设值的卫星，所述预设值优选为15度；
判断当前卫星的高度角是否小于15度，若是，则当前卫星不参与接收机位置解算，否则，跳过当前卫星，继续判断下一颗卫星的高度角是否小于15度；
1.4、根据步骤1.3处理后的卫星的高度角θ和方位角α确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵WE、WN和WU；
一般情况下，卫星的高度角越低，观测量的质量越差，故在确定权值时应给予较小的权重，且本发明中E、N、U方向上的权值为优选加权策略具体参见图3，图3中从左至右依次为E方向上的加权策略、N方向上的加权策略和U方向上的加权策略，其中颜色越深，给予的权重越大，在迭代解算过程中，各个方向上的权值会随着接收机位置的变化而变化；
1.41、根据步骤1.3处理后的每颗卫星的高度角θ和方位角α确定当前历元中第i颗卫星在站心坐标系中E、N、U方向上的权值w2 E,i、w2 N,i和w2 U,i，其中i＝(1，2，...，n)，n为当前历元中所述接收机收到的卫星数量；
优选的，w2 E,i、w2 N,i和w2 U,i的估计策略如下：
1.42、根据所述E、N、U方向上各颗卫星的权值确定当前历元在E、N、U方向上关于各颗卫星的权重矩阵，所述权重矩阵如下：
WE＝diag{w2 E,1,w2 E,2,...,w2 E,n}
WN＝diag{w2 N,1,w2 N,2,...,w2 N,n}
WU＝diag{w2 U,1,w2 U,2,...,w2 U,n}；
1.5、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
1.51、分步计算当前历元和上一历元的站心坐标差，计算公式如下：
ΔEE＝(GTWEG)-1GTWEΔρ
ΔEN＝(GTWNG)-1GTWNΔρ
ΔEU＝(GTWUG)-1GTWUΔρ
其中，ΔE为当前历元和上一历元的站心坐标差，，ΔEE表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差，ΔEN表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差，ΔEU表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，G为方向余弦矩阵，Δρ为修正后的伪距残差；
其中G为雅可比矩阵，(Ix,Iy,Iz)为接收机到卫星的单位观测矢量；
其中(xk,yk,zk)为待求接收机的位置，(xi,yi,zi)为第i颗卫星在地心地固坐标系中的坐标，ρc为经卫星时钟误差、电离层延时误差、对流层延时误差修正后的伪距，公式表述为：
ρc i＝ρi+δti-Ii-Ti
其中，ρi为第i颗卫星到接收机的伪距观测值，δtu为接收机时钟误差引起的等效距离误差，δt i为卫星时钟误差引起的等效距离误差，Ii为电离层延时引起的等效距离误差，Ti为对流层延时引起的等效距离误差，优选的，卫星时钟误差用广播星历中的卫星钟差参数修正，电离层延时误差用Klobuchar模型修正，对流层延误差时用Saastamoinen模型修正，具体的修正模型可根据实际需求进行调整；
1.52、根据牛顿迭代法，重复步骤1.2至1.51，直到小于预设值时，得到所述接收机在当前历元E、N、U方向上的站心坐标的估值如下：
其中，eE,k，nE,k，uE,k分别表示当前历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，eE,k-1，nE,k-1，uE,k-1分别表示上一历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，ΔeE，ΔnE，ΔuE分别表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差；eN,k，nN,k，uN,k分别表示当前历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，eN,k-1，nN,k-1，uN,k-1分别表示上一历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，ΔeN，ΔnN，ΔuN分别表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差；eU,k，nU,k，uU,k分别表示当前历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，eU,k-1，nU,k-1，uU,k-1分别表示上一历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，ΔeU，ΔnU，ΔuU分别表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，k表示迭代次数，且k为大于1的正整数；
优选的，小于预设值是指
1.6根据步骤1.5估计的接收机位置，分别取不同加权策略下E、N、U方向上的最优值经坐标变换后作为所述接收机当前历元的定位结果；
1.61、取eE,k，nN,k，uU,k作为当前历元接收机在站心坐标系中的坐标，所述接收机的站心坐标如下：
1.62、根据所述接收机的站心坐标和坐标变换公式进行坐标变换，得到所述接收机的单点定位结果；
所述坐标变换公式为：其中，x，y，z分别表示接收机在地心地固坐标系X轴、Y轴和Z轴方向上的坐标分量，S为坐标变换矩阵，λ为接收机位置的大地经度，φ为接收机位置的大地纬度。
2、将最小二乘法求得的接收机在所有历元中的定位结果取平均，得到所述接收机的精度高的静态伪距单点定位结果；
根据取平均公式将最小二乘法求得的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，所述取平均公式如下：
其中，为各个历元的定位结果在X轴、Y轴和Z轴方向上坐标分量的平均值，m为历元总数；
若用户需要大地坐标，可将地心地固坐标转换为大地坐标，相应的转换公式为：
其中，λ,h分别为大地纬度、大地经度和大地高，e表示椭球偏心率，N表示基准椭球体的卯酉圆曲率半径；
以上述计算过程为例，利用2016年5月17日在福建省福州市闽江学院广成楼楼顶用一台双频GNSS接收机采集的约26小时数据进行实验，该站点精确坐标用精密单点定位(Precision Point Positioning，PPP)解算得到，数据处理时，卫星系统选用美国的GPS(Global Positioning System)和中国的北斗系统(BeiDou Navigation SatelliteSystem，BDS)，星历使用广播星历，电离层用Klobuchar模型改正，对流层用Saastamoinen模型改正，实验结果表明：相较于高度角加权最小二乘估计方法，本发明所述方法定位精度在水平方向和垂直方向均得到提高；
由图4可知，本发明基于分步加权最小二乘估计的精度高的静态伪距单点定位方法与高度角加权的最小二乘方法定位相比，本发明在水平和垂直分量上定位精度比高度角加权的最小二乘方法定位分别提高了12.1％和24.7％，在三维精度提高了13.7％，具体数值如表1所示：
表1
实施例三
请参照图2，一种精度高的静态伪距单点定位的装置1，包括存储器2、处理器3及存储在存储器2上并可在处理器3上运行的计算机程序，所述处理器3执行所述程序时实现实施例一中的各个步骤。
综上所述，本发明提供的一种精度高的静态伪距单点定位方法及装置，通过计算每个历元中每颗卫星的高度角和方位角，删除高度角小于预设值的卫星，并根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵，根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置，将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，通过多历元解算并且求平均，有效的削弱伪距误差对单点定位的影响，从而提高定位精度且无需与固定基站配合，通过确定当前历元中每颗卫星在站心坐标系中不同方向上的权值，并根据不同方向上的权值确定当前历元各颗卫星在E、N、U方向上的权重矩阵，不仅提高了定位的准确性，还可满足用户对不同方向上定位精度的需求，通过分步计算当前历元和上一历元的站心坐标差，并在不断迭代后将不同方向上的最优估值eE,k，nN,k，uU,k作为当前历元接收机坐标，通过坐标变换公式对接收机的站心坐标进行转换，得到接收机的单点定位结果，便于使用者使用，通过取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，有利于消除定位误差，提高精度高的静态伪距单点定位的精确度。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等同变换，或直接或间接运用在相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (6)

1.一种精度高的静态伪距单点定位方法，其特征在于，包括步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果；
步骤S13包括：
S131、根据步骤S12处理后的每颗卫星的高度角θ和方位角α确定当前历元中第i颗卫星在站心坐标系中E、N、U方向上的权值w2 E,i、w2 N,i和w2 U,i，其中i＝(1，2，...，n)，n为当前历元中所述接收机收到的卫星数量：
S132、根据所述E、N、U方向上各颗卫星的权值确定当前历元在E、N、U方向上关于各颗卫星的权重矩阵，所述权重矩阵如下：
WE＝diag{w2 E,1,w2 E,2,...,w2 E,n}
WN＝diag{w2 N,1,w2 N,2,...,w2 N,n}
WU＝diag{w2 U,1,w2 U,2,...,w2 U,n}；
步骤S14包括：
S141、分步计算当前历元和上一历元的站心坐标差，计算公式如下：
ΔEE＝(GTWEG)-1GTWEΔρ
ΔEN＝(GTWNG)-1GTWNΔρ
ΔEU＝(GTWUG)-1GTWUΔρ
其中，ΔE为当前历元和上一历元的站心坐标差，ΔEE表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差，ΔEN表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差，ΔEU表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，G为方向余弦矩阵，Δρ为修正后的伪距残差；
S142、根据牛顿迭代法，重复步骤S11至S141，直到小于预设值时，得到所述接收机在当前历元E、N、U方向上的站心坐标的估值如下：
其中，eE,k，nE,k，uE,k分别表示当前历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，eE,k-1，nE,k-1，uE,k-1分别表示上一历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，ΔeE，ΔnE，ΔuE分别表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差；eN,k，nN,k，uN,k分别表示当前历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，eN,k-1，nN,k-1，uN,k-1分别表示上一历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，ΔeN，ΔnN，ΔuN分别表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差；eU,k，nU,k，uU,k分别表示当前历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，eU,k-1，nU,k-1，uU,k-1分别表示上一历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，ΔeU，ΔnU，ΔuU分别表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，k表示迭代次数，且k为大于1的正整数；
S143、取eE,k，nN,k，uU,k作为当前历元接收机在站心坐标系中的坐标，所述接收机的站心坐标如下：
S144、根据所述接收机的站心坐标和坐标变换公式进行坐标变换，得到所述接收机的单点定位结果；
所述坐标变换公式为：其中，x，y，z分别表示接收机在地心地固坐标系X轴、Y轴和Z轴方向上的坐标分量，S为坐标变换矩阵，λ为接收机位置的大地经度，φ为接收机位置的大地纬度。
2.根据权利要求1所述的精度高的静态伪距单点定位方法，其特征在于，步骤S11包括：
S111、判断当前历元是否为第一历元，若是，则执行步骤S112，否则，执行步骤S113；
S112、根据卫星位置以及接收机位置的初始值将计算卫星的高度角和方位角；
S113、根据卫星位置和上一历元估计的接收机位置计算卫星的高度角和方位角。
3.根据权利要求1所述的精度高的静态伪距单点定位方法，其特征在于，步骤S2包括：
根据取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，所述取平均公式如下：
其中，为各个历元的位置在X轴、Y轴和Z轴方向上坐标分量的平均值，m为历元总数。
4.一种精度高的静态伪距单点定位的装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
S1、预设多个历元，并对每个历元分别执行以下步骤：
S11、根据当前历元中各颗卫星的位置以及接收机的位置计算每颗卫星的高度角和方位角；
S12、删除高度角小于预设值的卫星；
S13、根据步骤S12处理后的卫星的高度角和方位角确定当前历元各颗卫星在站心坐标系E、N、U方向上的权重矩阵；
S14、根据所述权重矩阵通过最小二乘估计方法分步估计所述接收机在当前历元的位置；
S2、将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果；
步骤S13包括：
S131、根据步骤S12处理后的每颗卫星的高度角θ和方位角α确定当前历元中第i颗卫星在站心坐标系中E、N、U方向上的权值w2 E,i、w2 N,i和w2 U,i，其中i＝(1，2，...，n)，n为当前历元中所述接收机收到的卫星数量：
S132、根据所述E、N、U方向上各颗卫星的权值确定当前历元在E、N、U方向上关于各颗卫星的权重矩阵，所述权重矩阵如下：
WE＝diag{w2 E,1,w2 E,2,...,w2 E,n}
WN＝diag{w2 N,1,w2 N,2,...,w2 N,n}
WU＝diag{w2 U,1,w2 U,2,...,w2 U,n}；
步骤S14包括：
S141、分步计算当前历元和上一历元的站心坐标差，计算公式如下：
ΔEE＝(GTWEG)-1GTWEΔρ
ΔEN＝(GTWNG)-1GTWNΔρ
ΔEU＝(GTWUG)-1GTWUΔρ
其中，ΔE为当前历元和上一历元的站心坐标差，ΔEE表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差，ΔEN表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差，ΔEU表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，G为方向余弦矩阵，Δρ为修正后的伪距残差；
S142、根据牛顿迭代法，重复步骤S11至S141，直到小于预设值时，得到所述接收机在当前历元E、N、U方向上的站心坐标的估值如下：
其中，eE,k，nE,k，uE,k分别表示当前历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，eE,k-1，nE,k-1，uE,k-1分别表示上一历元根据E方向权重矩阵WE所求得的接收机在站心坐标系中的坐标，ΔeE，ΔnE，ΔuE分别表示根据E方向权重矩阵WE所求得的当前历元和上一历元的站心坐标差；eN,k，nN,k，uN,k分别表示当前历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，eN,k-1，nN,k-1，uN,k-1分别表示上一历元根据N方向权重矩阵WN所求得的接收机在站心坐标系中的坐标，ΔeN，ΔnN，ΔuN分别表示根据N方向权重矩阵WN所求得的当前历元和上一历元的站心坐标差；eU,k，nU,k，uU,k分别表示当前历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，eU,k-1，nU,k-1，uU,k-1分别表示上一历元根据U方向权重矩阵WU所求得的接收机在站心坐标系中的坐标，ΔeU，ΔnU，ΔuU分别表示根据U方向权重矩阵WU所求得的当前历元和上一历元的站心坐标差，k表示迭代次数，且k为大于1的正整数；
S143、取eE,k，nN,k，uU,k作为当前历元接收机在站心坐标系中的坐标，所述接收机的站心坐标如下：
S144、根据所述接收机的站心坐标和坐标变换公式进行坐标变换，得到所述接收机的单点定位结果；
所述坐标变换公式为：其中，x，y，z分别表示接收机在地心地固坐标系X轴、Y轴和Z轴方向上的坐标分量，S为坐标变换矩阵，λ为接收机位置的大地经度，φ为接收机位置的大地纬度。
5.根据权利要求4所述的精度高的静态伪距单点定位的装置，其特征在于，步骤S11包括：
S111、判断当前历元是否为第一历元，若是，则执行步骤S112，否则，执行步骤S113；
S112、根据卫星位置以及接收机位置的初始值将计算卫星的高度角和方位角；
S113、根据卫星位置和上一历元估计的接收机位置计算卫星的高度角和方位角。
6.根据权利要求4所述的精度高的静态伪距单点定位的装置，其特征在于，步骤S2包括：
根据取平均公式将计算得到的接收机在各个历元中的位置取平均，得到所述接收机的精度高的静态伪距单点定位结果，所述取平均公式如下：
其中，为各个历元的位置在X轴、Y轴和Z轴方向上坐标分量的平均值，m为历元总数。
CN202011514751.7A 2018-10-08 2018-10-08 一种精度高的静态伪距单点定位方法及装置 Active CN112731493B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011514751.7A CN112731493B (zh) 2018-10-08 2018-10-08 一种精度高的静态伪距单点定位方法及装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202011514751.7A CN112731493B (zh) 2018-10-08 2018-10-08 一种精度高的静态伪距单点定位方法及装置
CN201811166743.0A CN109444931B (zh) 2018-10-08 2018-10-08 一种静态伪距单点定位的方法及装置

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201811166743.0A Division CN109444931B (zh) 2018-10-08 2018-10-08 一种静态伪距单点定位的方法及装置

Publications (2)

Publication Number Publication Date
CN112731493A CN112731493A (zh) 2021-04-30
CN112731493B true CN112731493B (zh) 2023-12-15

Family

ID=65544760

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202011514751.7A Active CN112731493B (zh) 2018-10-08 2018-10-08 一种精度高的静态伪距单点定位方法及装置
CN201811166743.0A Active CN109444931B (zh) 2018-10-08 2018-10-08 一种静态伪距单点定位的方法及装置

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201811166743.0A Active CN109444931B (zh) 2018-10-08 2018-10-08 一种静态伪距单点定位的方法及装置

Country Status (1)

Country Link
CN (2) CN112731493B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110673170A (zh) * 2019-09-17 2020-01-10 闽江学院 一种动态单点定位精度的测试方法及终端
CN111830536A (zh) * 2020-08-10 2020-10-27 南京林业大学 一种结合地形影响的卫星可见性判断方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105589087A (zh) * 2016-01-29 2016-05-18 申研 一种高精度实时卫星定位装置及其方法
CN105629278A (zh) * 2014-11-21 2016-06-01 桂林电子科技大学 一种高精度gnss伪距单点定位的互差中值加权定位方法
CN106501828A (zh) * 2016-09-26 2017-03-15 闽江学院 一种基于模糊逻辑加权的高精度伪距单点定位方法
CN108594275A (zh) * 2018-04-26 2018-09-28 桂林电子科技大学 一种北斗+gps双模单点定位方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4910676B2 (ja) * 2006-12-15 2012-04-04 セイコーエプソン株式会社 演算回路、測位装置、指標算出方法及びプログラム
CN104931981A (zh) * 2015-05-11 2015-09-23 中国科学院光电研究院 一种基于信噪比波动信息的gnss抗多径选星方法
CN105044747B (zh) * 2015-08-27 2017-09-15 交通信息通信技术研究发展中心 一种基于多星共视和滤波的时间同步装置及其方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629278A (zh) * 2014-11-21 2016-06-01 桂林电子科技大学 一种高精度gnss伪距单点定位的互差中值加权定位方法
CN105589087A (zh) * 2016-01-29 2016-05-18 申研 一种高精度实时卫星定位装置及其方法
CN106501828A (zh) * 2016-09-26 2017-03-15 闽江学院 一种基于模糊逻辑加权的高精度伪距单点定位方法
CN108594275A (zh) * 2018-04-26 2018-09-28 桂林电子科技大学 一种北斗+gps双模单点定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
卫星定位中的一种分步加权解算方法;何瑞珠;刘成;黄康;;天文研究与技术(第01期);全文 *

Also Published As

Publication number Publication date
CN112731493A (zh) 2021-04-30
CN109444931B (zh) 2020-12-15
CN109444931A (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
CN110275192B (zh) 一种基于智能手机的高精度单点定位方法与装置
CN107193029B (zh) 北斗三频信号的网络rtk基准站间模糊度快速确定方法
CN110045407B (zh) 一种分布式伪卫星/gnss优化定位方法
CN112285752B (zh) 一种定位精度高的单点定位的方法及装置
CN106842268B (zh) 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法
Li et al. Multi-constellation GNSS PPP instantaneous ambiguity resolution with precise atmospheric corrections augmentation
CN108415049A (zh) 提高网络rtk双差宽巷模糊度固定正确率的方法
CN111751853A (zh) 一种gnss双频载波相位整周模糊度解算方法
US20170269225A1 (en) Navigation Satellite Wide-Lane Bias Determination and Over-Range Adjustment System and Method
CN109085628A (zh) 一种整周模糊度的固定方法及系统
US20050027450A1 (en) Altitude aiding in a satellite positioning system
CN109683182A (zh) 一种引入系统间偏差的gnss多模组合差分定位方法
CN113325453B (zh) 基于参数约束的gnss非差模糊度确定方法及快速定位方法
CN109212562A (zh) 一种生成载波伪距观测量的方法
CN112731493B (zh) 一种精度高的静态伪距单点定位方法及装置
Wu et al. Statistical modeling for the mitigation of GPS multipath delays from day-to-day range measurements
Pan et al. Performance assessment of real-time multi-GNSS integrated PPP with uncombined and ionospheric-free combined observables
CN117826200A (zh) 基于PPP-B2b的海上实时精密定位方法、系统和介质
Du et al. Analyzing the contributions of multi-GNSS and INS to the PPP-AR outage re-fixing
Li et al. A grid-based ionospheric weighted method for PPP-RTK with diverse network scales and ionospheric activity levels
Akpınar et al. Determining the coordinates of control points in hydrographic surveying by the precise point positioning method
Zeng et al. Computationally efficient dual-frequency uncombined precise orbit determination based on IGS clock datum
CN110109163A (zh) 一种附有高程约束的精密单点定位方法
Wielgosz et al. Network approach to precise GPS navigation
Abd Rabbou et al. PPP accuracy enhancement using GPS/GLONASS Observations in Kinematic mode

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